The limits of mtDNA phylogeography: complex patterns of population history in a highly structured Iberian lizard are only revealed by the use of nuclear markers

Integrative analysis reveals the divergence and speciation between sister Sooty Copper butterflies Lycaena bleusei and L. tityrus

The comparison of closely related taxa is cornerstone in biology, as understanding mechanisms leading up to differentiation in relation to extant shared characters are powerful tools in interpreting the evolutionary process. Hotspots of biodiversity such as the west-Mediterranean, where many lineages meet are ideal grounds to study these processes. We set to explore the interesting example of Sooty Copper butterflies: widespread Eurasian Lycaena tityrus (Poda, 1761) comes into contact in Iberia with closely related and local endemic, L. bleusei (Oberthür, 1884), which hasn't always been considered a distinct species. An integrative analysis was designed, combining the use of extensive molecular data (five genes), geometric morphometrics analyses, verified and up-to-date distribution data, and environmental niche modelling, aimed at deciphering their true relationship, their placement within European Lycaena and trace their evolutionary history. We revealed several levels of differentiation: L. bleusei and L. tityrus appear to be reciprocally monophyletic independent gene-pools, distinct in all genes analysed, having mutually diverged 4.8 Ma ago. L. tityrus but not L. bleusei, further displays a genetic structure compatible with several glacial refugia, where populations assignable to infraspecific taxa surface. Conversely, L. bleusei shows a loss in mtDNA diversity in relation to nuDNA. Morphological analyses differentiate both species according to size and shape but also discriminate strong seasonal and sexual traits and a geographical phenotype segregation in L. tityrus. Finally, updated distribution and its modelling for current and glacial timeframes reveal both species respond differently to environmental variables, defining a mostly parapatric distribution and an overlapping belt where sympatry was recovered. During the last glacial maximum, a wider expansion in L. bleusei distribution explains current isolated populations. Our study highlights the importance of gathering several lines of evidence when deciphering the relationships between closely related populations in the fringe of cryptic species realm.

Genetic structuring among colonies of a pantropical seabird: Implication for subspecies validation and conservation

Investigations of the genetic structure of populations over the entire range of a species yield valuable information about connectivity among populations. Seabirds are an intriguing taxon in this regard because they move extensively when not breeding, facilitating intermixing of populations, but breed consistently on the same isolated islands, restricting gene flow among populations. The degree of genetic structuring of populations varies extensively among seabird species but they have been understudied in their tropical ranges. Here, we address this across a broad spatial scale by using microsatellite and mitochondrial data to explore the population connectivity of 13 breeding populations representing the six subspecies of the white-tailed tropicbird (Phaethon lepturus) in the Atlantic, Indian, and Pacific Oceans. Our primary aim was to identify appropriate conservation units for this little known species. Three morphometric characters were also examined in the subspecies. We found a clear pattern of population structuring with four genetic groups. The most ancient and the most isolated group was in the northwestern Atlantic Ocean. The South Atlantic populations and South Mozambique Channel population on Europa were genetically isolated and may have had a common ancestor. Birds from the Indo-Pacific region showed unclear and weak genetic differentiation. This structuring was most well defined from nuclear and mtDNA markers but was less well resolved by morphological data. The validity of classifying white-tailed tropicbirds into six distinct subspecies is discussed in light of our new findings. From a conservation standpoint our results highlight that the three most threatened conservation units for this species are the two subspecies of the tropical North and South Atlantic Oceans and that of Europa Island in the Indian Ocean.

Population genetics of wild Macaca fascicularis with low-coverage shotgun sequencing of museum specimens

OBJECTIVES

Long-tailed macaques (Macaca fascicularis) are widely distributed throughout the mainland and islands of Southeast Asia, making them a useful model for understanding the complex biogeographical history resulting from drastic changes in sea levels throughout the Pleistocene. Past studies based on mitochondrial genomes (mitogenomes) of long-tailed macaque museum specimens have traced their colonization patterns throughout the archipelago, but mitogenomes trace only the maternal history. Here, our objectives were to trace phylogeographic patterns of long-tailed macaques using low-coverage nuclear DNA (nDNA) data from museum specimens.

METHODS

We performed population genetic analyses and phylogenetic reconstruction on nuclear single nucleotide polymorphisms (SNPs) from shotgun sequencing of 75 long-tailed macaque museum specimens from localities throughout Southeast Asia.

RESULTS

We show that shotgun sequencing of museum specimens yields sufficient genome coverage (average ~1.7%) for reconstructing population relationships using SNP data. Contrary to expectations of divergent results between nuclear and mitochondrial genomes for a female philopatric species, phylogeographical patterns based on nuclear SNPs proved to be closely similar to those found using mitogenomes. In particular, population genetic analyses and phylogenetic reconstruction from the nDNA identify two major clades within M. fascicularis: Clade A includes all individuals from the mainland along with individuals from northern Sumatra, while Clade B consists of the remaining island-living individuals, including those from southern Sumatra.

CONCLUSIONS

Overall, we demonstrate that low-coverage sequencing of nDNA from museum specimens provides enough data for examining broad phylogeographic patterns, although greater genome coverage and sequencing depth would be needed to distinguish between very closely related populations, such as those throughout the Philippines.

Phylogeny of caucasian rock lizards (Darevskia) and other true lizards based on mitogenome analysis: Optimisation of the algorithms and gene selection

We generated a phylogeny for Caucasian rock lizards (Darevskia), and included six other families of true lizards (Lacertini), based on complete mitochondrial genome analysis. Next-generation sequencing (NGS) of genomic DNA was used to obtain 16 new mitogenomes of Darevskia. These, along with 35 sequences downloaded from GenBank: genera Darevskia, Zootoca, Podarcis, Phoenicolacerta, Takydromus, Lacerta, and Eremias-were used in the analysis. All four analytical methods (Bayesian Inference, BI; Maximum Likelihood, ML; Maximum Parsimony, MP; and Neighbor-Joining, NJ) showed almost congruent intra-generic topologies for Darevskia and other lizard genera. However, ML and NJ methods on one side, and BI and MP methods on the other harvested conflicting phylogenies. The ML/NJ topology supports earlier published separation of Darevskia into three mitochondrial clades (Murphy, Fu, Macculloch, Darevsky, and Kupinova, 2000), but BI and MP topologies support that the basal branching occurred between D. parvula from the western Lesser Caucasus and the rest of Darevskia. All topologies altered the phylogenetic position of some individual species, including D. daghestanica, D. derjugini, and D. chlorogaster. Reanalysis after excluding four saturated genes from the data set, and excluding genus Eremias gives fully convergent topologies. The most basal branching for true lizards was between Far Eastern Takydromus and the Western Eurasian genera (BI). Comparing phylogenetic performance of individual genes relative to whole mitogenome data, concatenated 16S RNA (the least saturated gene in our analyses) and Cytochrome b genes generate a robust phylogeny that is fully congruent with that based on the complete mitogenome.

Genetic diversity and inferred dispersal history of the Schlegel's Japanese Gecko (Gekko japonicus) in Northeast Asia based on population genetic analyses and paleo-species distribution modelling

To understand the genetic diversity and dispersal history of Schlegel's Japanese gecko (Gekko japonicus), we performed genetic analyses and paleo-species distributional modelling. For the genetic analysis, we analyzed mitochondrial DNA (mtDNA) (cytochrome b [Cytb] and NADH dehydrogenase 2 [ND2]) and seven microsatellite loci of 353 individuals from 11 populations (2 east coast China, 4 west and central coast Japan and 5 Korea). For the paleo-species distribution modelling, we used 432 occurrence data points (125 China, 291 Japan and 16 Korea) over the Pleistocene and Holocene. China is inferred to be the source population, which had higher genetic diversity (mtDNA) and more private alleles (mtDNA) compared to Japanese and Korean populations. Differences between the three counties were very small in the mtDNA haplotype network despite some genetic structure among the three countries. Microsatellite analysis inferred that genetic exchange has actively occurred among the Chinese, Japanese and Korean populations. Suitable habitats in Japan should have been plentiful by the mid-Holocene, but have only recently become available in Korea. These results suggest that dispersal of G. japonicus occurred after the Holocene warming from the east coast of China to the west and central coasts of Japan and Korea, and gene flow is actively occurring among the three countries.

Population Genetic Analyses Using 10 New Polymorphic Microsatellite Loci Confirms Genetic Subdivision within the Olm, Proteus anguinus

We provide a comparative population genetic study of the elusive amphibian, Proteus anguinus, by comparing the genetic diversity and divergence among 4 cave populations (96 individuals) sampled in the Dinaric Karst of Croatia. We developed 10 variable microsatellite markers using pyrosequencing and applied them to the 4 selected populations belonging to 4 different cave systems. The results showed strong genetic differentiation between the 4 caves corroborating with previous findings suggesting that Proteus might comprise several unrecognized taxa. Our results confirmed that gene flow should be high within the caves, whereas it is low between hydrographic systems since geological periods. Finally, we conclude that the high genetic subdivision suggests the necessity of treating the 4 studied Proteus populations as evolutionary significant units.

Evolutionary history of the endemic water shrew Neomys anomalus: Recurrent phylogeographic patterns in semi-aquatic mammals of the Iberian Peninsula

The Cabrera's water shrew (Neomys anomalus) is a small semi-aquatic mammal whose taxonomic status was recently elevated from subspecies to species; as a consequence of this change, this species is now endemic to the Iberian Peninsula. In this study, we looked at its evolutionary history by combining phylogeography, the spatial distribution of genetic diversity, and species distribution modeling. To perform these analyses, we used noninvasive samples collected across the species distribution range and sequenced partial mitochondrial cytochrome b and D-loop genes. Maximum-likelihood and Bayesian phylogenetic trees derived from these sequences indicated that N. anomalus is divided into two main phylogroups that correlate strongly with geography, with two contact zones between the groups that showed limited spatial mixing between them. River basins were responsible for only a small percentage of the structure of the genetic diversity of this species despite its riparian habitat. The nucleotide diversity variation map showed the highest genetic diversity to be in the north of the Iberian Peninsula. Finally, species distribution modeling allowed the inference of an optimal area during the Last Interglacial in the north of the Iberian Peninsula, and multiple glacial refugia during the Last Glacial Maximum. The phylogeographic pattern of N. anomalus is strikingly similar to that of another semi-aquatic Iberian mammal, the Pyrenean desman (Galemys pyrenaicus), revealing how Pleistocene glaciations could have had equivalent effects on species of similar ecology and distribution. This phylogeographic structure is consistent with N. anomalus having been isolated for long periods in multiple glacial refugia within the Iberian Peninsula, in agreement with the "refugia-within-refugia" hypothesis, and further supporting its status as a distinct species.

Phylogeographic analysis reveals high genetic structure with uniform phenotypes in the paper wasp Protonectarina sylveirae (Hymenoptera: Vespidae)

Swarm-founding wasps are endemic and common representatives of neotropical fauna and compose an interesting social tribe of vespids, presenting both complex social characteristics and uncommon traits for a eusocial group, such as the absence of castes with distinct morphology. The paper wasp Protonectarina sylveirae (Saussure) presents a broad distribution from Brazil, Argentina and Paraguay, occurring widespread in the Atlantic rainforest and arboreal Caatinga, being absent in the Amazon region. Given the peculiar distribution among swarm-founding wasps, an integrative approach to reconstruct the evolutionary history of P. sylveirae in a spatial-temporal framework was performed to investigate: the presence of genetic structure and its relationship with the geography, the evolution of distinct morphologic lineages and the possible historical event(s) in Neotropical region, which could explain the observed phylogeographic pattern. Individuals of P. sylveirae were obtained from populations of 16 areas throughout its distribution for DNA extraction and amplification of mitochondrial genes 12S, 16S and COI. Analysis of genetic diversity, construction of haplotype net, analysis of population structure and dating analysis of divergence time were performed. A morphometric analysis was also performed using 8 measures of the body of the adult (workers) to test if there are morphological distinction among populations. Thirty-five haplotypes were identified, most of them exclusively of a group and a high population structure was found. The possibility of genetic divergence because of isolation by distance was rejected. Morphological analysis pointed to a great uniformity in phenotypes, with only a small degree of differentiation between populations of south and the remaining. Divergence time analysis showed a Middle/Late Miocene origin, a period where an extensive marine ingression occurred in South America. Divergence of haplogroups began from the Plio/Pleistocene boundary and the last glacial maximum most likely modeled the current distribution of species, even though it was not the cause of genetic breaks.

Ecotypes and evolutionary significant units in endangered North African gazelles

Conservation planning of threatened taxa relies upon accurate data on systematics, ecological traits and suitable habitats. The genus Gazella includes taxa with distinct morphologies and ecological traits, but close phylogenetic relationships. The North African Gazella cuvieri and Gazella leptoceros loderi share morphological and physiological characters but the former is darker and found in mountain areas, while the latter is lighter and associated with sand dunes. Here we aim to assess the genetic distinctiveness of these taxa, to characterize their ecological niches and to identify potential occurrence areas, by analysing 327 samples across North-West Africa. Phylogenetic analyses based on mitochondrial (CYTB) and five nuclear gene fragments (KCAS, LAC, SPTBN1, PRKCI and THYR) show that both taxa comprise a single monophyletic group. However, ecological niche-based modelling suggests that populations of these taxa occupy distinct geographic areas and specific environments. Predicted areas of sympatry were restricted, as a consequence of local sharp transitions in climatic traits. The lack of genetic differentiation between these taxa suggests they should be lumped into G. cuvieri, while ecological and morphological differences indicate they correspond to distinct ecotypes. Conservation planning of G. cuvieri should consider the preservation of both mountain and lowland ecotypes to maintain the overall adaptive potential of the species. This integrative approach provides valuable insights in identifying evolutionary units and should be extended to other gazelles.

Postglacial dispersal patterns and mitochondrial genetic structure of the Pyrenean desman (Galemys pyrenaicus) in the northwestern region of the Iberian Peninsula

The genetic structure of small semiaquatic animals may be influenced by dispersal across both rivers and land. The relative importance of these two modes of dispersal may vary across different species and with ecological conditions and evolutionary periods. The Pyrenean desman (Galemys pyrenaicus) is an endemic mammal of the Iberian Peninsula with a strong phylogeographic structure and semiaquatic habits, thus making it an ideal model to study the effects of river and overland dispersal on its genetic structure. Thanks to different types of noninvasive samples, we obtained an extensive sampling of the Pyrenean desman from the northwestern region of the Iberian Peninsula and sequenced two mitochondrial DNA fragments. We then analyzed, using an isolation‐by‐distance approach, the correlation between phylogenetic distances and geographical distances measured along both river networks and land to infer the relative importance of river and overland dispersal. We found that the correlations in the whole area and in a large basin were consistent with an effect of overland dispersal, which may be due to the postglacial colonization of new territories using terrestrial corridors and, possibly, a more extensive fluvial network that may have been present during the Holocene. However, in a small basin, likely to be less influenced by the impact of ancient postglacial dispersal, the correlations suggested significant overall effects of both overland and river dispersal, as expected for a semiaquatic mammal. Therefore, different scales and geographical regions reflect different aspects of the evolutionary history and ecology of this semiaquatic species using this isolation‐by‐distance method. The results we obtained may have crucial implications for the conservation of the Pyrenean desman because they reinforce the importance of interbasin dispersal for this species in the studied area and the need to protect the whole riverine ecosystem, including rivers, upland streams and terrestrial corridors between basins.

Integrative inference of population history in the Ibero-Maghrebian endemic Pleurodeles waltl (Salamandridae)

Inference of population histories from the molecular signatures of past demographic processes is challenging, but recent methodological advances in species distribution models and their integration in time-calibrated phylogeographic studies allow detailed reconstruction of complex biogeographic scenarios. We apply an integrative approach to infer the evolutionary history of the Iberian ribbed newt (Pleurodeles waltl), an Ibero-Maghrebian endemic with populations north and south of the Strait of Gibraltar. We analyzed an extensive multilocus dataset (mitochondrial and nuclear DNA sequences and ten polymorphic microsatellite loci) and found a deep east-west phylogeographic break in Iberian populations dating back to the Plio-Pleistocene. This break is inferred to result from vicariance associated with the formation of the Guadalquivir river basin. In contrast with previous studies, North African populations showed exclusive mtDNA haplotypes, and formed a monophyletic clade within the Eastern Iberian lineage in the mtDNA genealogy. On the other hand, microsatellites failed to recover Moroccan populations as a differentiated genetic cluster. This is interpreted to result from post-divergence gene flow based on the results of IMA2 and Migrate analyses. Thus, Moroccan populations would have originated after overseas dispersal from the Iberian Peninsula in the Pleistocene, with subsequent gene flow in more recent times, implying at least two trans-marine dispersal events. We modeled the distribution of the species and of each lineage, and projected these models back in time to infer climatically favourable areas during the mid-Holocene, the last glacial maximum (LGM) and the last interglacial (LIG), to reconstruct more recent population dynamics. We found minor differences in climatic favourability across lineages, suggesting intraspecific niche conservatism. Genetic diversity was significantly correlated with the intersection of environmental favourability in the LIG and LGM, indicating that populations of P. waltl are genetically more diverse in regions that have remained environmentally favourable through the last glacial cycle, particularly southern Iberia and northern Morocco. This study provides novel insights into the relative roles of geology and climate on the biogeography of a biodiversity hotspot.

Islands within an island: Population genetic structure of the endemic Sardinian newt, Euproctus platycephalus

The identification of historic and contemporary barriers to dispersal is central to the conservation of endangered amphibians, but may be hindered by their complex life history and elusive nature. The complementary information generated by mitochondrial (mtDNA) and microsatellite markers generates a valuable tool in elucidating population structure and the impact of habitat fragmentation. We applied this approach to the study of an endangered montane newt, Euproctus platycephalus. Endemic to the Mediterranean island of Sardinia, it is threatened by anthropogenic activity, disease, and climate change. We have demonstrated a clear hierarchy of structure across genetically divergent and spatially distinct subpopulations. Divergence between three main mountain regions dominated genetic partitioning with both markers. Mitochondrial phylogeography revealed a deep division dating to ca. 1 million years ago (Mya), isolating the northern region, and further differentiation between the central and southern regions ca. 0.5 Mya, suggesting an association with Pleistocene severe glacial oscillations. Our findings are consistent with a model of southward range expansion during glacial periods, with postglacial range retraction to montane habitat and subsequent genetic isolation. Microsatellite markers revealed further strong population structure, demonstrating significant divergence within the central region, and partial differentiation within the south. The northern population showed reduced genetic diversity. Discordance between mitochondrial and microsatellite markers at this scale indicated a further complexity of population structure, in keeping with male-biased dispersal and female philopatry. Our study underscores the need to elucidate cryptic population structure in the ecology and conservation strategies for endangered island-restricted amphibians, especially in the context of disease and climate change.

Population structure and dispersal routes of an invasive parasite, Fascioloides magna, in North America and Europe

Background

Fascioloides magna (Trematoda: Fasciolidae) is an important liver parasite of a wide range of free-living and domestic ruminants; it represents a remarkable species due to its large spatial distribution, invasive character, and potential to colonize new territories. The present study provides patterns of population genetic structure and admixture in F. magna across all enzootic regions in North America and natural foci in Europe, and infers migratory routes of the parasite on both continents.

Methods

In total, 432 individuals from five North American enzootic regions and three European foci were analysed by 11 microsatellite loci. Genetic data were evaluated by several statistical approaches: (i) the population genetic structure of F. magna was inferred using program STRUCTURE; (ii) the genetic interrelationships between populations were analysed by PRINCIPAL COORDINATES ANALYSIS; and (iii) historical dispersal routes in North America and recent invasion routes in Europe were explored using MIGRATE.

Results

The analysis of dispersal routes of the parasite in North America revealed west-east and south-north lineages that partially overlapped in the central part of the continent, where different host populations historically met. The exact origin of European populations of F. magna and their potential translocation routes were determined. Flukes from the first European focus, Italy, were related to F. magna from northern Pacific coast, while parasites from the Czech focus originated from south-eastern USA, particularly South Carolina. The Danube floodplain forests (third and still expanding focus) did not display relationship with any North American population; instead the Czech origin of the Danube population was indicated. A serial dilution of genetic diversity along the dispersion route across central and eastern Europe was observed. The results of microsatellite analyses were compared to previously acquired outputs from mitochondrial haplotype data and correlated with past human-directed translocations and natural migration of the final cervid hosts of F. magna.

Conclusions

The present study revealed a complex picture of the population genetic structure and interrelationships of North American and European populations, global distribution and migratory routes of F. magna and an origin of European foci.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1811-z) contains supplementary material, which is available to authorized users.

Phylogeny and population genetic structure of the ant genus Acropyga (Hymenoptera : Formicidae) in Papua New Guinea

Spatial isolation and geological history are important factors in the diversification and population differentiation of species. Here we describe distributional patterns of ants in the genus Acropyga across Papua New Guinea (PNG), a highly biodiverse but little-studied region. We estimate phylogenetic relationships among currently recognised species of Acropyga and assess population genetic structure of the widespread species, A. acutiventris, across lowland areas of the island. We find that species of Acropyga present in PNG diversified during the Pliocene, between six and two million years ago. Most species now exhibit a patchy distribution that does not show a strong signal of geological history. However, the population genetic structure of the widespread species A. acutiventris has been influenced by geography, habitat association and, possibly, historical habitat fragmentation. There is a significant effect of isolation-by-distance within continuous lowland forest, and proximity to Australia has had a larger impact in structuring populations of A. acutiventris in PNG than has the Central Papuan Cordillera. This study is the first to describe population genetic patterns of an ant species in Papua New Guinea.

Revisiting the Iberian honey bee (Apis mellifera iberiensis) contact zone: maternal and genome-wide nuclear variations provide support for secondary contact from historical refugia

Dissecting diversity patterns of organisms endemic to Iberia has been truly challenging for a variety of taxa, and the Iberian honey bee is no exception. Surveys of genetic variation in the Iberian honey bee are among the most extensive for any honey bee subspecies. From these, differential and complex patterns of diversity have emerged, which have yet to be fully resolved. Here, we used a genome-wide data set of 309 neutrally tested single nucleotide polymorphisms (SNPs), scattered across the 16 honey bee chromosomes, which were genotyped in 711 haploid males. These SNPs were analysed along with an intergenic locus of the mtDNA, to reveal historical patterns of population structure across the entire range of the Iberian honey bee. Overall, patterns of population structure inferred from nuclear loci by multiple clustering approaches and geographic cline analysis were consistent with two major clusters forming a well-defined cline that bisects Iberia along a northeastern-southwestern axis, a pattern that remarkably parallels that of the mtDNA. While a mechanism of primary intergradation or isolation by distance could explain the observed clinal variation, our results are more consistent with an alternative model of secondary contact between divergent populations previously isolated in glacial refugia, as proposed for a growing list of other Iberian taxa. Despite current intense honey bee management, human-mediated processes have seemingly played a minor role in shaping Iberian honey bee genetic structure. This study highlights the complexity of the Iberian honey bee patterns and reinforces the importance of Iberia as a reservoir of Apis mellifera diversity.

Population genetic structure, linkage disequilibrium and effective population size of conserved and extensively raised village chicken populations of Southern Africa

Extensively raised village chickens are considered a valuable source of biodiversity, with genetic variability developed over thousands of years that ought to be characterized and utilized. Surveys that can reveal a population's genetic structure and provide an insight into its demographic history will give valuable information that can be used to manage and conserve important indigenous animal genetic resources. This study reports population diversity and structure, linkage disequilibrium and effective population sizes of Southern African village chickens and conservation flocks from South Africa. DNA samples from 312 chickens from South African village and conservation flocks (n = 146), Malawi (n = 30) and Zimbabwe (n = 136) were genotyped using the Illumina iSelect chicken SNP60K BeadChip. Population genetic structure analysis distinguished the four conservation flocks from the village chicken populations. Of the four flocks, the Ovambo clustered closer to the village chickens particularly those sampled from South Africa. Clustering of the village chickens followed a geographic gradient whereby South African chickens were closer to those from Zimbabwe than to chickens from Malawi. Different conservation flocks seemed to have maintained different components of the ancestral genomes with a higher proportion of village chicken diversity found in the Ovambo population. Overall population LD averaged over chromosomes ranged from 0.03 ± 0.07 to 0.58 ± 0.41 and averaged 0.15 ± 0.16. Higher LD, ranging from 0.29 to 0.36, was observed between SNP markers that were less than 10 kb apart in the conservation flocks. LD in the conservation flocks steadily decreased to 0.15 (PK) and 0.24 (VD) at SNP marker interval of 500 kb. Genomewide LD decay in the village chickens from Malawi, Zimbabwe and South Africa followed a similar trend as the conservation flocks although the mean LD values for the investigated SNP intervals were lower. The results suggest low effective population sizes particularly in the conservation flocks. The utility and limitations of the iselect chicken SNP60K in village chicken populations is discussed.

Assessing the phylogeographic history of the montane caddisfly Thremma gallicum using mitochondrial and restriction-site-associated DNA (RAD) markers

Repeated Quaternary glaciations have significantly shaped the present distribution and diversity of several European species in aquatic and terrestrial habitats. To study the phylogeography of freshwater invertebrates, patterns of intraspecific variation have been examined primarily using mitochondrial DNA markers that may yield results unrepresentative of the true species history. Here, population genetic parameters were inferred for a montane aquatic caddisfly, Thremma gallicum, by sequencing a 658-bp fragment of the mitochondrial CO1 gene, and 12,514 nuclear RAD loci. T. gallicum has a highly disjunct distribution in southern and central Europe, with known populations in the Cantabrian Mountains, Pyrenees, Massif Central, and Black Forest. Both datasets represented rangewide sampling of T. gallicum. For the CO1 dataset, this included 352 specimens from 26 populations, and for the RAD dataset, 17 specimens from eight populations. We tested 20 competing phylogeographic scenarios using approximate Bayesian computation (ABC) and estimated genetic diversity patterns. Support for phylogeographic scenarios and diversity estimates differed between datasets with the RAD data favouring a southern origin of extant populations and indicating the Cantabrian Mountains and Massif Central populations to represent highly diverse populations as compared with the Pyrenees and Black Forest populations. The CO1 data supported a vicariance scenario (north-south) and yielded inconsistent diversity estimates. Permutation tests suggest that a few hundred polymorphic RAD SNPs are necessary for reliable parameter estimates. Our results highlight the potential of RAD and ABC-based hypothesis testing to complement phylogeographic studies on non-model species.

Mitochondrial DNA variation and phylogenetic relationships among five tuna species based on sequencing of D-loop region

In order to assess the DNA sequence variation and phylogenetic relationship among five tuna species (Auxis thazard, Euthynnus affinis, Katsuwonus pelamis, Thunnus tonggol, and T. albacares) out of all four tuna genera, partial sequences of the mitochondrial DNA (mtDNA) D-loop region were analyzed. The estimate of intra-specific sequence variation in studied species was low, ranging from 0.027 to 0.080 [Kimura's two parameter distance (K2P)], whereas values of inter-specific variation ranged from 0.049 to 0.491. The longtail tuna (T. tonggol) and yellowfin tuna (T. albacares) were found to share a close relationship (K2P = 0.049) while skipjack tuna (K. pelamis) was most divergent studied species. Phylogenetic analysis using Maximum-Likelihood (ML) and Neighbor-Joining (NJ) methods supported the monophyletic origin of Thunnus species. Similarly, phylogeny of Auxis and Euthynnus species substantiate the monophyly. However, results showed a distinct origin of K. pelamis from genus Thunnus as well as Auxis and Euthynnus. Thus, the mtDNA D-loop region sequence data supports the polyphyletic origin of tuna species.

Genetic structure of the endangered Greater Short-horned Lizard (Phrynosoma hernandesi) in Canada: evidence from mitochondrial and nuclear genes

The northern edge of the range of the Greater Short-horned Lizard (Phrynosoma hernandesi Girard, 1858) occurs in western Canada, where the species has “endangered” status and exhibits a patchy distribution. Phylogenetic inference and genetic analyses were employed to investigate the genetic structure of P. hernandesi throughout its Canadian range. One nuclear and two mitochondrial DNA genes were sequenced from 94 lizard tail tips. Overall, sequences from lizards from both Alberta and Saskatchewan displayed very little variability, and the consistent clustering of all the P. hernandesi mitochondrial and nuclear DNA sequences from Canada in both phylogenetic and population genetic analyses is consistent with the lizards from all sampled localities having originated from a single glacial refugium, and with being, until recently (or currently) interconnected genetically. The genetic data obtained so far furnish no information useful for interpreting the species’ present-day patchy distribution patterns or for formulating conservation strategies.

Real-time assessment of hybridization between wolves and dogs: combining noninvasive samples with ancestry informative markers

Wolves and dogs provide a paradigmatic example of the ecological and conservation implications of hybridization events between wild and domesticated forms. However, our understanding of such implications has been traditionally hampered by both high genetic similarity and the difficulties in obtaining tissue samples (TS), which limit our ability to assess ongoing hybridization events. To assess the occurrence and extension of hybridization in a pack of wolf-dog hybrids in northwestern Iberia, we compared the power of 52 nuclear markers implemented on TS with a subset of 13 ancestry informative markers (AIMs) typed in noninvasive samples (NIS). We demonstrate that the 13 AIMs are as accurate as the 52 markers that were chosen without regard to the power to differentiate between wolves and dogs, also having the advantage of being rapidly screened on NIS. The efficiency of AIMs significantly outperformed ten random sets of similar size and an additional commercial set of 18 markers. Bayesian clustering analysis implemented on AIMs and NIS identified nine hybrids, two wolves and two dogs. Four hybrids were unambiguously assigned to F1xWolf backcrosses. Our approach (AIMs + NIS) overcomes previous difficulties related to sample availability and informative power of markers, allowing a quick identification of wolf-dog hybrids in the first phases of hybridization episodes. This provides managers with a reliable tool to evaluate hybridization and estimate the success of their actions. This approach may be easily adapted for other pairs of wild/domesticated species, thus improving our understanding of the introgression of domestication genes into natural populations.

Phylogeography and postglacial recolonization of Europe by Rhinolophus hipposideros: evidence from multiple genetic markers

The demographic history of Rhinolophus hipposideros (lesser horseshoe bat) was reconstructed across its European, North African and Middle-Eastern distribution prior to, during and following the most recent glaciations by generating and analysing a multimarker data set. This data set consisted of an X-linked nuclear intron (Bgn; 543 bp), mitochondrial DNA (cytb-tRNA-control region; 1630 bp) and eight variable microsatellite loci for up to 373 individuals from 86 localities. Using this data set of diverse markers, it was possible to determine the species' demography at three temporal stages. Nuclear intron data revealed early colonization into Europe from the east, which pre-dates the Quaternary glaciations. The mtDNA data supported multiple glacial refugia across the Mediterranean, the largest of which were found in the Ibero-Maghreb region and an eastern location (Anatolia/Middle East)-that were used by R. hipposideros during the most recent glacial cycles. Finally, microsatellites provided the most recent information on these species' movements since the Last Glacial Maximum and suggested that lineages that had diverged into glacial refugia, such as in the Ibero-Maghreb region, have remained isolated. These findings should be used to inform future conservation management strategies for R. hipposideros and show the power of using a multimarker data set for phylogeographic studies.

Molecular phylogeny of the Thyropygus allevatus group of giant millipedes and some closely related groups

Giant cylindrical millipedes of the family Harpagophoridae, especially species of the genus Thyropygus, are broadly distributed in Thailand and nearby countries. They show a great deal of variation in body size, color patterns and gonopodal characters. Phylogenetic analyses of 26 nominate species from six genera in the subfamilies Harpagophorinae and Rhynchoproctinae, as well as nine new morphotypes (regarded as new species), were performed with the DNA sequences from two mitochondrial gene fragments (16S rRNA and COI). The genus Thyropygus (Harpagophorinae) was recovered as monophyletic under all analyses, whilst the representatives of Rhynchoproctinae also formed a monophyletic group. However, the analyses suggested that the tribe Gonoplectini should be elevated to a separate subfamily, Gonoplectinae. The molecular analyses were largely (but not totally) congruent with, and so supported the usefulness of, gonopodal characters for the classification and identification of harpagophorid millipedes, and additionally supported previous studies on the delimitation of species and subgroups. This is the first molecular study inside the family Harpagophoridae and provides the basis for further studies of the evolutionary processes and biogeographic patterns of millipedes in Southeast Asia.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression

BACKGROUND

The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure.

RESULTS

The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France.

CONCLUSIONS

The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

Causes and evolutionary consequences of population subdivision of an Iberian mountain lizard, Iberolacerta monticola

Aim

The study of the factors that influence population connectivity and spatial distribution of genetic variation is crucial for understanding speciation and for predicting the effects of landscape modification and habitat fragmentation, which are considered severe threats to global biodiversity. This dual perspective is obtained from analyses of subalpine mountain species, whose present distribution may have been shaped both by cyclical climate changes over ice ages and anthropogenic perturbations of their habitats. Here, we examine the phylogeography, population structure and genetic diversity of the lacertid lizard Iberolacerta monticola, an endemism considered to be facing a high risk of extinction in several populations.

Location

Northwestern quadrant of the Iberian Peninsula.

Methods

We analyzed the mtDNA variation at the control region (454 bp) and the cytochrome b (598 bp) loci, as well as at 10 nuclear microsatellite loci from 17 populations throughout the distribution range of the species.

Results

According to nuclear markers, most sampling sites are defined as distinct, genetically differentiated populations, and many of them show traces of recent bottlenecks. Mitochondrial data identify a relatively old, geographically restricted lineage, and four to six younger geographically vicariant sister clades, whose origin may be traced back to the mid-Pleistocene revolution, with several subclades possibly associated to the mid-Bruhnes transition. Geographic range fragmentation of one of these clades, which includes lowland sites, is very recent, and most likely due to the accelerated loss of Atlantic forests by human intervention.

Main Conclusions

Altogether, the data fit a “refugia within refugia” model, some lack of pattern uniformity notwithstanding, and suggest that these mountains might be the cradles of new species of Iberolacerta. However, the changes operated during the Holocene severely compromise the long-term survival of those genetic lineages more exposed to the anthropogenic perturbations of their habitats.

Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus

Background

Species with strict ecological requirements may provide new insights into the forces that shaped the geographic variation of genetic diversity. The Pyrenean desman, Galemys pyrenaicus, is a small semi-aquatic mammal that inhabits clean streams of the northern half of the Iberian Peninsula and is endangered in most of its geographic range, but its genetic structure is currently unknown. While the stringent ecological demands derived from its aquatic habitat might have caused a partition of the genetic diversity among river basins, Pleistocene glaciations would have generated a genetic pattern related to glacial refugia.

Results

To study the relative importance of historical and ecological factors in the genetic structure of G. pyrenaicus, we used mitochondrial and intronic sequences of specimens covering most of the species range. We show, first, that the Pyrenean desman has very low levels of genetic diversity compared to other mammals. In addition, phylogenetic and dating analyses of the mitochondrial sequences reveal a strong phylogeographic structure of a Middle Pleistocene origin, suggesting that the main lineages arose during periods of glacial isolation. Furthermore, both the spatial distribution of nuclear and mitochondrial diversity and the results of species distribution modeling suggest the existence of a major glacial refugium in the northwestern part of the Iberian Peninsula. Finally, the main mitochondrial lineages show a striking parapatric distribution without any apparent exchange of mitochondrial haplotypes between the lineages that came into secondary contact (although with certain permeability to nuclear genes), indicating incomplete mixing after the post-glacial recolonization. On the other hand, when we analyzed the partition of the genetic diversity among river basins, the Pyrenean desman showed a lower than expected genetic differentiation among main rivers.

Conclusions

The analysis of mitochondrial and intronic markers in G. pyrenaicus showed the predominant effects of Pleistocene glaciations on the genetic structure of this species, while the distribution of the genetic diversity was not greatly influenced by the main river systems. These results and, particularly, the discovery of a marked phylogeographic structure, may have important implications for the conservation of the Pyrenean desman.

Analyses of the population structure in a global collection of Phytophthora nicotianae isolates inferred from mitochondrial and nuclear DNA sequences

Genetic variation within the heterothallic cosmopolitan plant pathogen Phytophthora nicotianae was determined in 96 isolates from a wide range of hosts and geographic locations by characterizing four mitochondrial (10% of the genome) and three nuclear loci. In all, 52 single-nucleotide polymorphisms (SNPs) (an average of 1 every 58 bp) and 313 sites with gaps representing 5,450 bases enabled the identification of 50 different multilocus mitochondrial haplotypes. Similarly, 24 SNPs (an average of 1 every 69 bp), with heterozygosity observed at each locus, were observed in three nuclear regions (hyp, scp, and β-tub) differentiating 40 multilocus nuclear genotypes. Both mitochondrial and nuclear markers revealed a high level of dispersal of isolates and an inconsistent geographic structuring of populations. However, a specific association was observed for host of origin and genetic grouping with both nuclear and mitochondrial sequences. In particular, the majority of citrus isolates from Italy, California, Florida, Syria, Albania, and the Philippines clustered in the same mitochondrial group and shared at least one nuclear allele. A similar association was also observed for isolates recovered from Nicotiana and Solanum spp. The present study suggests an important role of nursery populations in increasing genetic recombination within the species and the existence of extensive phenomena of migration of isolates that have been likely spread worldwide with infected plant material.

Continuing fragmentation of a widespread species by geographical barriers as initial step in a land snail radiation on crete

The phylogeographic structure of the land snail Xerocrassa mesostena on Crete inferred from AFLP markers and mitochondrial cox1 sequences can be explained by three mechanisms: gene flow restriction, population expansion and leptokurtic dispersal. Gene flow restriction by geographic barriers caused subdivision of the gene pool into distinct clusters. Population expansion was probably facilitated by deforestation of Crete in the postglacial. Newly available areas were colonized by leptokurtic dispersal, i.e. slow active expansion resulting in isolation by distance within the clusters and occasional long distance dispersal events that resulted in departures from the isolation by distance model. Less than one percent of the AFLP markers show correlations with environmental variables. Random phylogeographic breaks in the distribution of the mitochondrial haplotype groups indicate that single locus markers, especially mitochondrial DNA, might result in a misleading picture of the phylogeographic structure of a species. Restriction of gene flow between metapopulations caused by geographical barriers can interact with sexual selection resulting in the differentiation of these metapopulations into separate species without noticeable ecological differentiation. Evidence for gene flow between parapatrically distributed evolutionary units representing different stages of the speciation process suggests that the ongoing process of fragmentation of the X. mesostena complex might be an example for parapatric speciation. The lack of ecological differentiation between these units confirms theoretical predictions that divergent selection for local adaptation is not required for rapid speciation.

Fine-scale genetic breaks driven by historical range dynamics and ongoing density-barrier effects in the estuarine seaweed Fucus ceranoides L

BACKGROUND

Factors promoting the emergence of sharp phylogeographic breaks include restricted dispersal, habitat discontinuity, physical barriers, disruptive selection, mating incompatibility, genetic surfing and secondary contact. Disentangling the role of each in any particular system can be difficult, especially when species are evenly distributed across transition zones and dispersal barriers are not evident. The estuarine seaweed Fucus ceranoides provides a good example of highly differentiated populations along its most persistent distributional range at the present rear edge of the species distribution, in NW Iberia. Intrinsic dispersal restrictions are obvious in this species, but have not prevented F. ceranoides from vastly expanding its range northwards following the last glaciation, implying that additional factors are responsible for the lack of connectivity between neighbouring southern populations. In this study we analyze 22 consecutive populations of F. ceranoides along NW Iberia to investigate the processes generating and maintaining the observed high levels of regional genetic divergence.

RESULTS

Variation at seven microsatellite loci and at mtDNA spacer sequences was concordant in revealing that Iberian F. ceranoides is composed of three divergent genetic clusters displaying nearly disjunct geographical distributions. Structure and AFC analyses detected two populations with an admixed nuclear background. Haplotypic diversity was high in the W sector and very low in the N sector. Within each genetic cluster, population structure was also pervasive, although shallower.

CONCLUSIONS

The deep divergence between sectors coupled with the lack of support for a role of oceanographic barriers in defining the location of breaks suggested 1) that the parapatric genetic sectors result from the regional reassembly of formerly vicariant sub-populations, and 2) that the genetic discontinuities at secondary contact zones (and elsewhere) are maintained despite normal migration rates. We conclude that colonization and immigration, as sources of gene-flow, have very different genetic effects. Migration between established populations is effectively too low to prevent their differentiation by drift or to smooth historical differences inherited from the colonization process. F. ceranoides, but possibly low-dispersal species in general, appear to be unified to a large extent by historical, non-equilibrium processes of extinction and colonization, rather than by contemporary patterns of gene flow.

Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula: multiple refugia, range expansions and secondary contact zones

BACKGROUND

The Iberian Peninsula is recognized as an important refugial area for species survival and diversification during the climatic cycles of the Quaternary. Recent phylogeographic studies have revealed Iberia as a complex of multiple refugia. However, most of these studies have focused either on species with narrow distributions within the region or species groups that, although widely distributed, generally have a genetic structure that relates to pre-Quaternary cladogenetic events. In this study we undertake a detailed phylogeographic analysis of the lizard species, Lacerta lepida, whose distribution encompasses the entire Iberian Peninsula. We attempt to identify refugial areas, recolonization routes, zones of secondary contact and date demographic events within this species.

RESULTS

Results support the existence of 6 evolutionary lineages (phylogroups) with a strong association between genetic variation and geography, suggesting a history of allopatric divergence in different refugia. Diversification within phylogroups is concordant with the onset of the Pleistocene climatic oscillations. The southern regions of several phylogroups show a high incidence of ancestral alleles in contrast with high incidence of recently derived alleles in northern regions. All phylogroups show signs of recent demographic and spatial expansions. We have further identified several zones of secondary contact, with divergent mitochondrial haplotypes occurring in narrow zones of sympatry.

CONCLUSIONS

The concordant patterns of spatial and demographic expansions detected within phylogroups, together with the high incidence of ancestral haplotypes in southern regions of several phylogroups, suggests a pattern of contraction of populations into southern refugia during adverse climatic conditions from which subsequent northern expansions occurred. This study supports the emergent pattern of multiple refugia within Iberia but adds to it by identifying a pattern of refugia coincident with the southern distribution limits of individual evolutionary lineages. These areas are important in terms of long-term species persistence and therefore important areas for conservation.

Interspecific X-chromosome and mitochondrial DNA introgression in the Iberian hare: selection or allele surfing?

Introgression from a resident species into an invading one is predicted to occur through the demographic process of "allele surfing," and to particularly affect genomic regions transmitted by the lower migrating sex, such as mtDNA. This could explain that northern Iberian populations of Lepus granatensis harbor high frequencies of mtDNA from L. timidus, an arctic hare it replaced there after deglaciation. We report that variation of introgressed timidus-like mtDNA reflects several predicted effects of this process: increasing frequency and diversity in the direction of expansion, strong perpendicular phylogeographic structure and signs of postglacial demographic growth. However, demographic inferences for the granatensis and timidus-like mtDNA lineages suggest the latter may have outcompeted the former in northern Iberia. Autosomal introgression occurs at low frequencies and species-wide rather than only in the north. If this difference with mtDNA resulted from sex-biased migration, an intermediate pattern should prevail for the X-chromosome, but we report species-wide and high-frequency introgression of an X-fragment. Either selection favored this ubiquitous X-introgression, or more complex postglacial expansion patterns prevailed, with different consequences depending on the genomic and geographic region. This illustrates the difficulty of distinguishing demographic and selective effects and the need for genome and species-wide based demographic models.

Population genetic patterns revealed by microsatellite data challenge the mitochondrial DNA based taxonomy of Astyanax in Mexico (Characidae, Teleostei)

Astyanax has become an important model system for evolutionary studies of cave animals. We investigated correlations of population genetic patterns revealed by microsatellite data and phylogeographic patterns shown by mitochondrial DNA sequences in Mexican cave and surface fish of the genus Astyanax (Characidae, Teleostei) to improve the understanding of the colonization history of this neotropical fish in Central and North America and to assess a recent taxonomic classification. The distribution of nuclear genotypes is not congruent with that of the mitochondrial clades. Admixture analyses suggest there has been nuclear gene flow between populations defined by different mitochondrial clades. The microsatellite data indicate that there was mitochondrial capture of a cave population from adjacent populations. Furthermore, gene flow also occurred between populations belonging to different nuclear genotypic clusters. This indicates that neither the nuclear genotypic clusters nor the mitochondrial clades represent independent evolutionary units, although the mitochondrial divergences are high and in a range usually characteristic for different fish species. This conclusion is supported by the presence of morphologically intermediate forms. Our analyses show that the Trans-Mexican Volcanic Belt limited gene flow, but has been crossed by Astyanax several times. In Yucatán, where obvious geographic barriers are missing, the incongruence between the distribution of nuclear and mitochondrial markers reflects random colonization events caused by inundations or marine transgressions resulting in random phylogeographic breaks. Thus, conclusions about the phylogeographic history and even more about the delimitation of species should not be based on single genetic markers.