Supplementary Material Table S1 Pairwise comparisons of inter-island genetic differentiation in the louse Degeeriella regalis. Data are based on 1099 bp of mtDNA sequences (partial COI and 12S) collected from 110 individual lice from 8 Galápagos hawk island populations. Pairwise Fst values estimated from Kimura 2-parameter genetic distances are below the diagonal and the P values of those comparisons are above the diagonal.
| Fernandina
| Pinta
| Marchena
| Isabela
| Santiago
| Pinzón
| Española
| Santa Fe
| Fernandina | -
| <0.00001
| <0.00001
| 0.264
| 0.409
| <0.00001
| <0.00001
| <0.00001
| Pinta
| 0.7718
| -
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| Marchena
| 0.8889
| 0.8396
| -
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| Isabela
| 0.0200
| 0.7946
| 0.8887
| -
| 0.072
| <0.00001
| <0.00001
| <0.00001
| Santiago
| 0.0001
| 0.6860
| 0.7641
| 0.0272
| -
| <0.001
| <0.00001
| <0.00001
| Pinzón
| 0.8282
| 0.8778
| 0.9669
| 0.8074
| 0.5354
| -
| <0.00001
| <0.001
| Española
| 0.9670
| 0.9550
| 0.9874
| 0.9622
| 0.9064
| 1.0
| -
| <0.001
| Santa Fe
| 0.9192
| 0.8808
| 0.8573
| 0.9180
| 0.8265
| 0.9724
| 0.9876
| -
| Table S2 Pairwise comparisons of inter-island genetic differentiation in the louse Colpocephalum turbinatum. Data are based on 950 bp of mtDNA sequences (partial COI and 12S) collected from 127 individual lice from 8 Galápagos hawk island populations. Pairwise Fst values estimated from Kimura 2-parameter genetic distances are below the diagonal and the P values of those comparisons are above the diagonal.
| Fernandina
| Pinta
| Marchena
| Isabela
| Santiago
| Pinzón
| Española
| Santa Fe
| Fernandina
| -
| <0.00001
| <0.00001
| 0.279
| <0.00001
| <0.01
| <0.05
| <0.00001
| Pinta
| 0.3895
| -
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| Marchena
| 0.7953
| 0.9417
| -
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| <0.00001
| Isabela
| 0.0286
| 0.5143
| 0.8128
| -
| <0.00001
| <0.01
| <0.05
| <0.00001
| Santiago
| 0.3785
| 0.8200
| 0.8734
| 0.4295
| -
| <0.01
| 0.108
| <0.00001
| Pinzón
| 0.3387
| 0.8702
| 0.8061
| 0.3336
| 0.1946
| -
| 0.0750
| <0.001
| Española
| 0.2606
| 0.9
| 0.8038
| 0.2541
| 0.1013
| 0.1447
| -
| <0.001
| Santa Fe
| 0.6162
| 0.9
| 0.8955
| 0.6646
| 0.8462
| 0.7874
| 0.7948
| -
|
Table S3. Pairwise comparisons of inter-island genetic differentiation in the lousefly Icosta nigra. Data are based on 937 bp of mtDNA sequences (partial COI and 12S) collected from 117 individual flies from 7 Galápagos hawk island populations. Pairwise Fst values estimated from Kimura 2-parameter genetic distances are below the diagonal and the P values of those comparisons are above the diagonal.
| Fernandina
| Pinta
| Isabela
| Santiago
| Pinzón
| Española
| Santa Fe
| Fernandina | -
| <0.00001
| 1
| 0.173
| <0.00001
| <0.00001
| <0.00001
| Pinta
| 1.0
| -
| <0.00001
| <0.00001
| 0.092
| 0.216
| 0.260
| Isabela
| -0.0171
| 0.8830
| -
| 0.160
| <0.00001
| <0.00001
| <0.00001
| Santiago
| 0.0953
| 0.7402
| 0.0353
| -
| <0.00001
| <0.00001
| <0.00001
| Pinzón
| 0.7692
| 0.1624
| 0.6561
| 0.5089
| -
| 0.1051
| 0.2770
| Española
| 0.8571
| 0.0051
| 0.7585
| 0.6358
| 0.0072
| -
| 0.4711
| Santa Fe
| 0.9161
| 0.241
| 0.7878
| 0.6712
| 0.0874
| 0.0244
| -
|
Table S4 Mdiv estimates of the non-equilibrium migration rates (Nefm below diagonals) and estimates of the divergence times (above diagonal) between population pairs of Buteo galapagoensis, using the variable 911 bp mtDNA dataset (divergence times were standardized by taking the product of * t, where = 2Nef).
| Fernandina
| Marchena
| Santiago
| Isabela
| Pinzón
| Pinta
| Santa Fe
| Española
| Fernandina
| -
| 1.08
| 1.15
| 0.46
| 1.34
| 1.09
| 1.48
| 1.56
| Marchena
| 0.02
| -
| 0
| 1.23
| 1.2
| 0
| 0
| 2.68
| Santiago
| 0.02
| 0.14
| -
| 1.03
| 1.48
| 0
| 0
| 2.77
| Isabela
| 0.40
| 0.10
| 0.06
| -
| 1.41
| 1.12
| 0.85
| 2.85
| Pinzón
| 0.02
| 0.02
| 0.02
| 0.62
| -
| 2.10
| 1.63
| 3.38
| Pinta
| 0.02
| 0.02
| 0.02
| 0.10
| 0.02
| -
| 0
| 3.17
| Santa Fe
| 0.02
| 0.02
| 1.76
| 0.10
| 0.02
| 2.74
| -
| 2.97
| Española
| 0.02
| 0.02
| 0.02
| 0.02
| 0.02
| 0.02
| 0.02
| -
|
Table S5 Mdiv estimates of the non-equilibrium migration rates (Nefm below diagonals) and estimates of the divergence times (above diagonal) between population pairs of Degeeriella regalis (divergence times were standardized by taking the product of * t, where = 2Nef).
| Fernandina
| Marchena
| Santiago
| Isabela
| Pinzón
| Pinta
| Santa Fe
| Española
| Fernandina
| -
| 2.03
| 0.19
| 0.48
| 0.88
| 1.83
| 4.28
| 7.84
| Marchena
| 0.01
| -
| 2.56
| 2.58
| 3.84
| 1.60
| 0.93
| 7.71
| Santiago
| 2.76
| 0.12
| -
| 0.39
| 2.04
| 1.33
| 3.17
| 8.06
| Isabela
| 1.19
| 0.01
| 0.9
| -
| 1.21
| 0.96
| 4.11
| 8.93
| Pinzón
| 0.01
| 0.01
| 0.43
| 0.01
| -
| 2.24
| 4.66
| 6.55
| Pinta
| 0.01
| 0.012
| 0.16
| 0.02
| 0.01
| -
| 3.35
| 6.80
| Santa Fe
| 0.01
| 0.09
| 0.03
| 0.01
| 0.01
| 0.01
| -
| 9.92
| Española
| 0.01
| 0.01
| 0.01
| 0.01
| 0.01
| 0.01
| 0.01
| -
|
Table S6 Mdiv estimates of the non-equilibrium migration rates (Nefm below diagonals) and estimates of the divergence times (above diagonal) between population pairs of Colpocephalum turbinatum (divergence times were standardized by taking the product of * t, where = 2Nef).
| Fernandina
| Marchena
| Santiago
| Isabela
| Pinzón
| Pinta
| Santa Fe
| Española
| Fernandina
| -
| 3.24
| 0.69
| 0.11
| 0.71
| 0.23
| 1.78
| 0.77
| Marchena
| 0.01
| -
| 3.56
| 2.74
| 3.08
| 4.01
| 4.93
| 2.76
| Santiago
| 0.17
| 0.01
| -
| 0.64
| 0.05
| 0.84
| 2.39
| 0.39
| Isabela
| 1.1
| 0.01
| 0.2
| -
| 0.36
| 0.30
| 1.25
| 0.62
| Pinzón
| 0.22
| 0.01
| 3.85
| 0.27
| -
| 1.05
| 2.33
| 0.02
| Pinta
| 0.02
| 0.01
| 0.01
| 0.08
| 0.01
| -
| 0.91
| 0.73
| Santa Fe
| 0.02
| 0.01
| 0.03
| 0.02
| 0.01
| 0.01
| -
| 1.97
| Española
| 0.17
| 0.01
| 1.09
| 0.42
| 4.99
| 0.02
| 0.02
| -
| Table S7 Mdiv estimates of the non-equilibrium migration rates (Nefm below diagonals) and estimates of the divergence times (above diagonal) between population pairs of Icosta nigra (divergence times were standardized by taking the product of * t, where = 2Nef).
| Fernandina
| Santiago
| Isabela
| Pinzón
| Pinta
| Santa Fe
| Española
| Fernandina
| -
| 0.09
| 0.01
| 0.27
| 0.50
| 0.48
| 1.13
| Santiago
| 0.3
| -
| 0.01
| 1.05
| 0.38
| 0.42
| 1.84
| Isabela
| 6.54
| 7.72
| -
| 0.46
| 0.81
| 0.72
| 0.63
| Pinzón
| 0.08
| 0.8
| 0.46
| -
| 0.08
| 0.18
| 0.06
| Pinta
| 0.01
| 0.09
| 0.11
| 0.4
| -
| 0.13
| 0.02
| Santa Fe
| 0.01
| 0.21
| 0.16
| 0.28
| 0.68
| -
| 0.42
| Española
| 0.09
| 0.34
| 0.28
| 4.75
| 3.67
| 9.8
| -
| |