The biogeography of mitochondrial and nuclear discordance in animals

Population genomics reveals that an anthropophilic population of Aedes aegypti mosquitoes in West Africa recently gave rise to American and Asian populations of this major disease vector

BACKGROUND

The mosquito Aedes aegypti is the main vector of dengue, Zika, chikungunya and yellow fever viruses. This major disease vector is thought to have arisen when the African subspecies Ae. aegypti formosus evolved from being zoophilic and living in forest habitats into a form that specialises on humans and resides near human population centres. The resulting domestic subspecies, Ae. aegypti aegypti, is found throughout the tropics and largely blood-feeds on humans.

RESULTS

To understand this transition, we have sequenced the exomes of mosquitoes collected from five populations from around the world. We found that Ae. aegypti specimens from an urban population in Senegal in West Africa were more closely related to populations in Mexico and Sri Lanka than they were to a nearby forest population. We estimate that the populations in Senegal and Mexico split just a few hundred years ago, and we found no evidence of Ae. aegypti aegypti mosquitoes migrating back to Africa from elsewhere in the tropics. The out-of-Africa migration was accompanied by a dramatic reduction in effective population size, resulting in a loss of genetic diversity and rare genetic variants.

CONCLUSIONS

We conclude that a domestic population of Ae. aegypti in Senegal and domestic populations on other continents are more closely related to each other than to other African populations. This suggests that an ancestral population of Ae. aegypti evolved to become a human specialist in Africa, giving rise to the subspecies Ae. aegypti aegypti. The descendants of this population are still found in West Africa today, and the rest of the world was colonised when mosquitoes from this population migrated out of Africa. This is the first report of an African population of Ae. aegypti aegypti mosquitoes that is closely related to Asian and American populations. As the two subspecies differ in their ability to vector disease, their existence side by side in West Africa may have important implications for disease transmission.

Genomewide SNP data reveal cryptic phylogeographic structure and microallopatric divergence in a rapids-adapted clade of cichlids from the Congo River

The lower Congo River is a freshwater biodiversity hot spot in Africa characterized by some of the world's largest rapids. However, little is known about the evolutionary forces shaping this diversity, which include numerous endemic fishes. We investigated phylogeographic relationships in Teleogramma, a small clade of rheophilic cichlids, in the context of regional geography and hydrology. Previous studies have been unable to resolve phylogenetic relationships within Teleogramma due to lack of variation in nuclear genes and discrete morphological characters among putative species. To sample more broadly across the genome, we analysed double-digest restriction-associated sequencing (ddRAD) data from 53 individuals across all described species in the genus. We also assessed body shape and mitochondrial variation within and between taxa. Phylogenetic analyses reveal previously unrecognized lineages and instances of microallopatric divergence across as little as ~1.5 km. Species ranges appear to correspond to geographic regions broadly separated by major hydrological and topographic barriers, indicating these features are likely important drivers of diversification. Mitonuclear discordance indicates one or more introgressive hybridization events, but no clear evidence of admixture is present in nuclear genomes, suggesting these events were likely ancient. A survey of female fin patterns hints that previously undetected lineage-specific patterning may be acting to reinforce species cohesion. These analyses highlight the importance of hydrological complexity in generating diversity in certain freshwater systems, as well as the utility of ddRAD-Seq data in understanding diversification processes operating both below and above the species level.

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.

Ancient and modern genome shuffling: Reticulate mito-nuclear phylogeny of four related allopatric species of Gyrodactylus von Nordmann, 1832 (Monogenea: Gyrodactylidae), ectoparasites on the Eurasian minnow Phoxinus phoxinus (L.) (Cyprinidae)

Phylogenetic analyses including four allopatric species of Gyrodactylus von Nordmann, 1832 on the Eurasian minnow Phoxinus phoxinus (L.) (Cyprinidae) revealed incongruence between the nuclear ITS1-5.8S-ITS2 and mitochondrial cox1 phylogenies due to ancient hybridisation. Gyrodactylus pannonicus Molnár, 1968 was sampled close to its type-locality, the upper reaches of River Tisza, tributary of Danube in the Black Sea Basin. Faunistic search detected three new related species with maximum composite likelihood distances in cox1 between 16.8-23.2% (tentatively 1.3 to 1.8 My of divergence). Gyrodactylus albolacustris n. sp. recorded in the White Sea Basin, eastern Baltic Basin and Mongolia was close to G. pannonicus in the nuclear ITS (divergence of 0.9%), but diverged in cox1 by 19.8%. The Mongolian isolate of G. albolacustris n. sp. diverged from the European isolates in cox1 by 8.9%, suggesting 0.7 My of isolation. The two other new species differed from G. pannonicus by >4% in ITS and some large indels in ITS1, and by >20% in cox1. Gyrodactylus danastriae n. sp. was found in River Strwiąż, a tributary of the River Dniester (Black Sea Basin) and was characterised by smaller size of anchors and by 29-41 bp dimorphic insertion in ITS1. Gyrodactylus botnicus n. sp. is considered endemic in the Baltic Basin, but was also found in the White Sea Basin as a postglacial immigrant, where it had hybridised with G. albolacustris n. sp. in spite of the high divergence in ITS (3.9%) and cox1 (22%). The discordant nuclear and mitochondrial phylogenies revealed an ancient mitochondrial introgression: G. albolacustris n. sp. was derived from a hybridisation combining proto-pannonicus ITS with proto-danastriae mitochondria, perhaps 1.3 My ago. The postglacial hybridisation of G. albolacustris n. sp. (as the donor of mtDNA^alb and ITS^alb) and G. botnicus n. sp. (donor of the ITS^bot) offered a model of shuffling of the genomic components: the process of the homogenisation and stabilisation of nuclear ITS (concerted evolution) and the lineage sorting has hardly begun.

Range expansion underlies historical introgressive hybridization in the Iberian hare

Introgressive hybridization is an important and widespread evolutionary process, but the relative roles of neutral demography and natural selection in promoting massive introgression are difficult to assess and an important matter of debate. Hares from the Iberian Peninsula provide an appropriate system to study this question. In its northern range, the Iberian hare, Lepus granatensis, shows a northwards gradient of increasing mitochondrial DNA (mtDNA) introgression from the arctic/boreal L. timidus, which it presumably replaced after the last glacial maximum. Here, we asked whether a south-north expansion wave of L. granatensis into L. timidus territory could underlie mtDNA introgression, and whether nuclear genes interacting with mitochondria (“mitonuc” genes) were affected. We extended previous RNA-sequencing and produced a comprehensive annotated transcriptome assembly for L. granatensis. We then genotyped 100 discovered nuclear SNPs in 317 specimens spanning the species range. The distribution of allele frequencies across populations suggests a northwards range expansion, particularly in the region of mtDNA introgression. We found no correlation between variants at 39 mitonuc genes and mtDNA introgression frequency. Whether the nuclear and mitochondrial genomes coevolved will need a thorough investigation of the hundreds of mitonuc genes, but range expansion and species replacement likely promoted massive mtDNA introgression.

Widespread introgression in deep-sea hydrothermal vent mussels

Background

The analysis of hybrid zones is crucial for gaining a mechanistic understanding of the process of speciation and the maintenance of species boundaries. Hybrid zones have been studied intensively in terrestrial and shallow-water ecosystems, but very little is known about their occurrence in deep-sea environments. Here we used diagnostic, single nucleotide polymorphisms in combination with one mitochondrial gene to re-examine prior hypotheses about a contact zone involving deep-sea hydrothermal vent mussels, Bathymodiolus azoricus and B. puteoserpentis, living along the Mid-Atlantic Ridge.

Results

Admixture was found to be asymmetric with respect to the parental species, while introgression was more widespread geographically than previously recognized. Admixed individuals with a majority of alleles from one of the parental species were most frequent in habitats corresponding to that species. Mussels found at a geographically intermediate vent field constituted a genetically mixed population that showed no evidence for hybrid incompatibilities, a finding that does not support a previously inferred tension zone model.

Conclusions

Our analyses indicate that B. azoricus and B. puteoserpentis hybridize introgressively across a large geographic area without evidence for general hybrid incompatibilities. While these findings shed new light onto the genetic structure of this hybrid zone, many aspects about its nature still remain obscure. Our study sets a baseline for further research that should primarily focus on the acquisition of additional mussel samples and environmental data, a detailed exploration of vent areas and hidden populations as well as genomic analyses in both mussel hosts and their bacterial symbionts.

Electronic supplementary material

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

When Viruses Don't Go Viral: The Importance of Host Phylogeographic Structure in the Spatial Spread of Arenaviruses

Many emerging infections are RNA virus spillovers from animal reservoirs. Reservoir identification is necessary for predicting the geographic extent of infection risk, but rarely are taxonomic levels below the animal species considered as reservoir, and only key circumstances in nature and methodology allow intrinsic virus-host associations to be distinguished from simple geographic (co-)isolation. We sampled and genetically characterized in detail a contact zone of two subtaxa of the rodent Mastomys natalensis in Tanzania. We find two distinct arenaviruses, Gairo and Morogoro virus, each spatially confined to a single M. natalensis subtaxon, only co-occurring at the contact zone’s centre. Inter-subtaxon hybridization at this centre and a continuum of quality habitat for M. natalensis show that both viruses have the ecological opportunity to spread into the other substaxon’s range, but do not, strongly suggesting host-intrinsic barriers. Such barriers could explain why human cases of another M. natalensis-borne arenavirus, Lassa virus, are limited to West Africa.

Reservoirs of zoonotic viruses are usually equated with a particular wildlife species. It is rarely assessed whether genetic groups below the species level may instead represent the actual reservoir, though this would have major implications on estimations of the zoonosis’ spatial distribution. Here we investigate whether geographically and genetically distinct subtaxa of the widespread African rodent Mastomys natalensis carry distinct arenaviruses, by sampling in detail across a contact zone of two of these subtaxa. Ongoing hybridization shows that individuals of the subtaxa are in direct physical contact, in principle allowing viral exchange, yet neither of the two arenaviruses -Gairo and Morogoro virus- were found to have crossed the zone. Such intraspecific genetic barriers to arenavirus spatial spread have important implications for our understanding of the related Lassa arenavirus, a pathogen potentially lethal to humans of which Mastomys natalensis is also the main reservoir. Although Lassa virus appears to infect several secondary hosts, its distribution is restricted to West Africa and matches that of another M. natalensis subtaxon. Our data thus indicates that it is because of M. natalensis intraspecific distinctions that the human Lassa fever endemic area has not expanded to the rest of sub-Saharan Africa.

A combination of long term fragmentation and glacial persistence drove the evolutionary history of the Italian wall lizard Podarcis siculus

BACKGROUND

The current distribution of genetic diversity is the result of a vast array of microevolutionary processes, including short-term demographic and ecological mechanisms and long-term allopatric isolation in response to Quaternary climatic fluctuations. We investigated past processes that drove the population differentiation and spatial genetic distribution of the Italian wall lizard Podarcis siculus by means of sequences of mitochondrial cytb (n = 277 from 115 localities) and nuclear mc1r and β-fibint7genes (n = 262 and n = 91, respectively) from all its distribution range. The pattern emerging from the genetic data was compared with current and past (last glacial maximum) species distribution modeling (SDM).

RESULTS

We identified seven deeply divergent parapatric clades which presumably remained isolated in different refugia scattered mainly throughout the Tyrrhenian coast. Conversely, the Adriatic coast showed only two haplogroups with low genetic variability. These results appear to agree with the SDM prediction at the last glacial maximum (LGM) indicating a narrow area of habitat suitability along the Tyrrhenian coast and much lower suitability along the Adriatic one. However, the considerable land exposure of the Adriatic coastline favored a glacial colonization of the Balkan Peninsula.

CONCLUSIONS

Our population-level historical demography showed a common trend consistent with glacial expansions and regional persistence during the last glacial maximum. This complex genetic signature appears to be inconsistent with the expectation of the expansion-contraction model and post-LGM (re)colonizations from southern refugia. Hence it is one of an increasing number of cases in which these assumptions are not met, indicating that long-term fragmentation and pre-LGM events such as glacial persistence were more prominent in shaping genetic variation in this temperate species.

Population genetics of the Manila clam (Ruditapes philippinarum) introduced in North America and Europe

Globally, the Manila clam (Ruditapes philippinarum) stands as the second most important bivalve species in fisheries and aquaculture. Native to the Pacific coast of Asia, it is now well-established in North America and Europe, where its on-going management reflects local economic interests. The historic record of transfers spans the 20th century and suggests sequential movement from Japan to North America, as a hitch-hiker on oysters, and then intentional introduction in Europe, but global genetic data are missing. We have studied mitochondrial DNA and microsatellite markers in nine populations from Asia, North America and Europe. The results from the two types of markers indicated a good concordance of present-day genetic structure with the reported history of clam transfers across continents, and no evidence of relevant concealed introductions from continental Asia in Europe and North America. However, European populations showed a loss of genetic variability and significant genetic differentiation as compared to their American counterparts. Our study shows that in spite of the increasing ease for species to spread out of their native range, in the case of the Manila clam this has not resulted in new invasion waves in the two studied continents.

Genomic Signatures of Historical Allopatry and Ecological Divergence in an Island Lizard

Geographical variation among contiguous populations is frequently attributed to ecological divergence or historical isolation followed by secondary contact. Distinguishing between these effects is key to studies of incipient speciation and could be revealed by different genomic signatures. We used RAD-seq analyses to examine morphologically divergent populations of the endemic lizard (Gallotia galloti) from the volcanic island of Tenerife. Previous analyses have suggested ecological and historical causes to explain the morphological diversity. Analyses of 276,483 single nucleotide polymorphisms (SNPs) from >20 Mbp of the genome revealed one genetically divergent population from Anaga, a region associated with divergent mtDNA lineages in other Tenerife endemics. This population also has a high number of private alleles, and its divergence can be explained by historical isolation. Bayesian outlier analyses identified a small proportion of SNPs as candidates for selection (0.04%) which were strongly differentiated between xeric and mesic habitat types. Individual testing for specific xeric–mesic selection using an alternative approach also supported ecological divergence in a similarly small proportion of SNPs. The study indicates the roles of both historical isolation and ecological divergence in shaping genomic diversity in G. galloti. However, north–south morphological divergence appears solely associated with the latter and likely involves a relatively small proportion of the genome.

Introgression and selection shaped the evolutionary history of sympatric sister-species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister-species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site-associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister-species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174-0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.

Matching loci surveyed to questions asked in phylogeography

Although mitochondrial DNA (mtDNA) has long been used for assessing genetic variation within and between populations, its workhorse role in phylogeography has been criticized owing to its single-locus nature. The only choice for testing mtDNA results is to survey nuclear loci, which brings into contrast the difference in locus effective size and coalescence times. Thus, it remains unclear how erroneous mtDNA-based estimates of species history might be, especially for evolutionary events in the recent past. To test the robustness of mtDNA and nuclear sequences in phylogeography, we provide one of the largest paired comparisons of summary statistics and demographic parameters estimated from mitochondrial, five Z-linked and 10 autosomal genes of 30 avian species co-distributed in the Caucasus and Europe. The results suggest that mtDNA is robust in estimating inter-population divergence but not in intra-population diversity, which is sensitive to population size change. Here, we provide empirical evidence showing that mtDNA was more likely to detect population divergence than any other single locus owing to its smaller Ne and thus faster coalescent time. Therefore, at least in birds, numerous studies that have based their inferences of phylogeographic patterns solely on mtDNA should not be readily dismissed.

Exon-intron structure and sequence variation of the calreticulin gene among Rhipicephalus sanguineus group ticks

Background

Calreticulin proteins (CRTs) are important components of tick saliva, which is involved in the blood meal success, pathogen transmission and host allergic responses. The characterization of the genes encoding for salivary proteins, such as CRTs, is pivotal to understand the mechanisms of tick-host interaction during blood meal and to develop tick control strategies based on their inhibition. In hard ticks, crt genes were shown to have only one intron with conserved position among species. In this study we investigated the exon-intron structure and variation of the crt gene in Rhipicephalus spp. ticks in order to assess the crt exon-intron structure and the potential utility of crt gene as a molecular marker.

Methods

We sequenced the exon-intron region of crt gene in ticks belonging to so-called tropical and temperate lineages of Rhipicephalus sanguineus (sensu lato), Rhipicephalus sp. I, Rhipicephalus sp. III, Rhipicephalus sp. IV, R. guilhoni, R. muhsamae and R. turanicus. Genetic divergence and phylogenetic relationships between the sequences obtained were estimated.

Results

All individuals belonging to the tropical lineage of R. sanguineus (s.l.), R. guilhoni, R. muhsamae, R. turanicus, Rhipicephalus sp. III and Rhipicephalus sp. IV analysed showed crt intron-present alleles. However, both crt intron-present and intron-absent alleles were found in Rhipicephalus sp. I and the temperate lineage of R. sanguineus (s.l.), showing the occurrence of an intraspecific intron presence-absence polymorphism. Phylogenetic relationships among the crt intron-present sequences showed distinct lineages for all taxa, with the tropical and temperate lineages of R. sanguineus (s.l.) being more closely related to each other.

Conclusions

We expanded previous studies about the characterization of crt gene in hard ticks. Our results highlighted a previously overlooked variation in the crt structure among Rhipicephalus spp., and among hard ticks in general. Notably, the intron presence/absence polymorphism observed herein can be a candidate study-system to investigate the early stages of intron gain/loss before fixation at species level and some debated questions about intron evolution. Finally, the sequence variation observed supports the suitability of the crt gene for molecular recognition of Rhipicephalus spp. and for phylogenetic studies in association with other markers.

Electronic supplementary material

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

DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland

We present a DNA barcoding study on the insect order Orthoptera that was generated in collaboration between four barcoding projects in three countries, viz. Barcoding Fauna Bavarica (Germany), German Barcode of Life, Austrian Barcode of Life and Swiss Barcode of Life. Our data set includes 748 COI sequences from 127 of the 162 taxa (78.4%) recorded in the three countries involved. Ninety-three of these 122 species (76.2%, including all Ensifera) can be reliably identified using DNA barcodes. The remaining 26 caeliferan species (families Acrididae and Tetrigidae) form ten clusters that share barcodes among up to five species, in three cases even across different genera, and in six cases even sharing individual barcodes. We discuss incomplete lineage sorting and hybridization as most likely causes of this phenomenon, as the species concerned are phylogenetically young and hybridization has been previously observed. We also highlight the problem of nuclear mitochondrial pseudogenes (numts), a known problem in the barcoding of orthopteran species, and the possibility of Wolbachia infections. Finally, we discuss the possible taxonomic implications of our barcoding results and point out future research directions.

Ancient Mitochondrial Capture as Factor Promoting Mitonuclear Discordance in Freshwater Fishes: A Case Study in the Genus Squalius (Actinopterygii, Cyprinidae) in Greece

Hybridization and incomplete lineage sorting are common confounding factors in phylogeny and speciation resulting in mitonuclear disparity. Mitochondrial introgression, a particular case of hybridization, may, in extreme cases, lead to replacement of the mitochondrial genome of one species with that of another (mitochondrial capture). We investigated mitochondrial introgression involving two species of the cyprinid genus Squalius in the western Peloponnese region of Greece using molecular and morphological data. We found evidence of complete mitochondrial introgression of Squalius keadicus into two populations recognized as Squalius peloponensis from the Miras and Pamissos River basins and a divergence of mitochondrial genomes of S. keadicus from the Evrotas basin from that of the introgressed populations dating from the Pleistocene. Secondary contact among basins is a possible factor in connection of the species and the introgression event. Morphological analyses support the hypothesis of mitochondrial introgression, as S. keadicus was different from the other three populations recognized as S. peloponensis, although significant differences were found among the four populations. Isolation by geographical barriers arose during Pleistocene in the western Peloponnese were the source of the evolution of the two reciprocally monophyletic subclades found in the S. keadicus mitochondrial clade, and the morphological differences found among the four populations. Along with the lack of structure in the nuclear genome in the three populations ascribed to S. peloponensis, this suggests an incipient speciation process occurring in these Squalius species in the western Peloponnese.

Discrete phenotypes are not underpinned by genome-wide genetic differentiation in the squat lobster Munida gregaria (Crustacea: Decapoda: Munididae): a multi-marker study covering the Patagonian shelf

BACKGROUND

DNA barcoding has demonstrated that many discrete phenotypes are in fact genetically distinct (pseudo)cryptic species. Genetically identical, isogenic individuals, however, can also express similarly different phenotypes in response to a trigger condition, e.g. in the environment. This alternative explanation to cryptic speciation often remains untested because it requires considerable effort to reject the hypothesis that the observed underlying genetic homogeneity of the different phenotypes may be trivially caused by too slowly evolving molecular markers. The widespread squat lobster Munida gregaria comprises two discrete ecotypes, gregaria s. str. and subrugosa, which were long regarded as different species due to marked differences in morphological, ecological and behavioral traits. We studied the morphometry and genetics of M. gregaria s. l. and tested (1) whether the phenotypic differences remain stable after continental-scale sampling and inclusion of different life stages, (2) and whether each phenotype is underpinned by a specific genotype.

RESULTS

A total number of 219 gregaria s. str. and subrugosa individuals from 25 stations encompassing almost entire range in South America were included in morphological and genetic analyses using nine unlinked hypervariable microsatellites and new COI sequences. Results from the PCA and using discriminant functions demonstrated that the morphology of the two forms remains discrete. The mitochondrial data showed a shallow, star-like haplotype network and complete overlap of genetic distances within and among ecotypes. Coalescent-based species delimitation methods, PTP and GMYC, coherently suggested that haplotypes of both ecotypes forms a single species. Although all microsatellite markers possess sufficient genetic variation, AMOVA, PCoA and Bayesian clustering approaches revealed no genetic clusters corresponding to ecotypes or geographic units across the entire South-American distribution. No evidence of isolation-by-distance could be detected for this species in South America.

CONCLUSIONS

Despite their pronounced bimodal morphologies and different lifestyles, the gregaria s. str. and subrugosa ecotypes form a single, dimorphic species M. gregaria s. l.. Based on adequate geographic coverage and multiple independent polymorphic loci, there is no indication that each phenotype may have a unique genetic basis, leaving phenotypic plasticity or localized genomic islands of speciation as possible explanations.

Genomic diversity in Onchocerca volvulus and its Wolbachia endosymbiont

Ongoing elimination efforts have altered the global distribution of Onchocerca volvulus, the agent of river blindness, and further population restructuring is expected as efforts continue. Therefore, a better understanding of population genetic processes and their effect on biogeography is needed to support elimination goals. We describe O. volvulus genome variation in 27 isolates from the early 1990s (before widespread mass treatment) from four distinct locales: Ecuador, Uganda, the West African forest and the West African savanna. We observed genetic substructuring between Ecuador and West Africa and between the West African forest and savanna bioclimes, with evidence of unidirectional gene flow from savanna to forest strains. We identified forest:savanna-discriminatory genomic regions and report a set of ancestry informative loci that can be used to differentiate between forest, savanna and admixed isolates, which has not previously been possible. We observed mito-nuclear discordance possibly stemming from incomplete lineage sorting. The catalogue of the nuclear, mitochondrial and endosymbiont DNA variants generated in this study will support future basic and translational onchocerciasis research, with particular relevance for ongoing control programmes, and boost efforts to characterize drug, vaccine and diagnostic targets.

Molecular phylogeny and patterns of diversification in syngnathid fishes

The family Syngnathidae is a large and diverse clade of morphologically unique bony fishes, with 57 genera and 300 described species of seahorses, pipefishes, pipehorses, and seadragons. They primarily inhabit shallow coastal waters in temperate and tropical oceans, and are characterized by a fused jaw, male brooding, and extraordinary crypsis. Phylogenetic relationships within the Syngnathidae remain poorly resolved due to lack of generic taxon sampling, few diagnostic morphological characters, and limited molecular data. The phylogenetic placement of the threatened, commercially exploited seahorses remains a topic of intense interest, with conflicting topologies based on morphology and predominantly mitochondrial genetic data. In this study, we integrate eight nuclear and mitochondrial markers and 17 morphological characters to investigate the phylogenetic structure of the family Syngnathidae at the generic level. We include 91 syngnathid species representing 48 of the 57 recognized genera, all major ocean basins, and a broad array of temperate and tropical habitats including rocky and coral reefs, sand and silt, mangroves, seagrass beds, estuaries, and rivers. Maximum likelihood and Bayesian analyses of 5160bp from eight loci produced high congruence among alternate topologies, defining well-supported and sometimes novel clades. We present a hypothesis that confirms a deep phylogenetic split between lineages with trunk- or tail-brood pouch placement, and provides significant new insights into the morphological evolution and biogeography of this highly derived fish clade. Based on the fundamental division between lineages - the tail brooding "Urophori" and the trunk brooding "Gastrophori" - we propose a revision of Syngnathidae classification into only two subfamilies: the Nerophinae and the Syngnathinae. We find support for distinct principal clades within the trunk-brooders and tail-brooders, the latter of which include seahorses, seadragons, independent lineages of pipehorses, and clades that originated in southern Australia and the Western Atlantic. We suggest the seahorse genus Hippocampus is of Indo-Pacific origin and its sister clade is an unexpected grouping of several morphologically disparate Indo-Pacific genera, including the Pacific pygmy pipehorses. Taxonomic revision is required for multiple genera, particularly to reflect deep evolutionary splits in nominal lineages from the Atlantic versus the Indo-Pacific.

Population subdivision of hydrothermal vent polychaete Alvinella pompejana across equatorial and Easter Microplate boundaries

Background

The Equator and Easter Microplate regions of the eastern Pacific Ocean exhibit geomorphological and hydrological features that create barriers to dispersal for a number of animals associated with deep-sea hydrothermal vent habitats. This study examined effects of these boundaries on geographical subdivision of the vent polychaete Alvinella pompejana. DNA sequences from one mitochondrial and eleven nuclear genes were examined in samples collected from ten vent localities that comprise the species’ known range from 23°N latitude on the East Pacific Rise to 38°S latitude on the Pacific Antarctic Ridge.

Results

Multi-locus genotypes inferred from these sequences clustered the individual worms into three metapopulation segments — the northern East Pacific Rise (NEPR), southern East Pacific Rise (SEPR), and northeastern Pacific Antarctic Ridge (PAR) — separated by the Equator and Easter Microplate boundaries. Genetic diversity estimators were negatively correlated with tectonic spreading rates. Application of the isolation-with-migration (IMa2) model provided information about divergence times and demographic parameters. The PAR and NEPR metapopulation segments were estimated to have split roughly 4.20 million years ago (Mya) (2.42–33.42 Mya, 95 % highest posterior density, (HPD)), followed by splitting of the SEPR and NEPR segments about 0.79 Mya (0.07–6.67 Mya, 95 % HPD). Estimates of gene flow between the neighboring regions were mostly low (2 Nm < 1). Estimates of effective population size decreased with southern latitudes: NEPR > SEPR > PAR.

Conclusions

Highly effective dispersal capabilities allow A. pompejana to overcome the temporal instability and intermittent distribution of active hydrothermal vents in the eastern Pacific Ocean. Consequently, the species exhibits very high levels of genetic diversity compared with many co-distributed vent annelids and mollusks. Nonetheless, its levels of genetic diversity in partially isolated populations are inversely correlated with tectonic spreading rates. As for many other vent taxa, this pioneering colonizer is similarly affected by local rates of habitat turnover and by major dispersal filters associated with the Equator and the Easter Microplate region.

Electronic supplementary material

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

Phylogeographic structure in long-tailed voles (Rodentia: Arvicolinae) belies the complex Pleistocene history of isolation, divergence, and recolonization of Northwest North America's fauna

Quaternary climate fluctuations restructured biodiversity across North American high latitudes through repeated episodes of range contraction, population isolation and divergence, and subsequent expansion. Identifying how species responded to changing environmental conditions not only allows us to explore the mode and tempo of evolution in northern taxa, but also provides a basis for forecasting future biotic response across the highly variable topography of western North America. Using a multilocus approach under a Bayesian coalescent framework, we investigated the phylogeography of a wide‐ranging mammal, the long‐tailed vole, Microtus longicaudus. We focused on populations along the North Pacific Coast to refine our understanding of diversification by exploring the potentially compounding roles of multiple glacial refugia and more recent fragmentation of an extensive coastal archipelago. Through a combination of genetic data and species distribution models (SDMs), we found that historical climate variability influenced contemporary genetic structure, with multiple isolated locations of persistence (refugia) producing multiple divergent lineages (Beringian or northern, southeast Alaska or coastal, and southern or continental) during glacial advances. These vole lineages all occur along the North Pacific Coast where the confluence of numerous independent lineages in other species has produced overlapping zones of secondary contact, collectively a suture zone. Finally, we detected high levels of neoendemism due to complex island geography that developed in the last 10,000 years with the rising sea levels of the Holocene.

Diversification in wild populations of the model organism Anolis carolinensis: A genome-wide phylogeographic investigation

The green anole (Anolis carolinensis) is a lizard widespread throughout the southeastern United States and is a model organism for the study of reproductive behavior, physiology, neural biology, and genomics. Previous phylogeographic studies of A. carolinensis using mitochondrial DNA and small numbers of nuclear loci identified conflicting and poorly supported relationships among geographically structured clades; these inconsistencies preclude confident use of A. carolinensis evolutionary history in association with morphological, physiological, or reproductive biology studies among sampling localities and necessitate increased effort to resolve evolutionary relationships among natural populations. Here, we used anchored hybrid enrichment of hundreds of genetic markers across the genome of A. carolinensis and identified five strongly supported phylogeographic groups. Using multiple analyses, we produced a fully resolved species tree, investigated relative support for each lineage across all gene trees, and identified mito‐nuclear discordance when comparing our results to previous studies. We found fixed differences in only one clade—southern Florida restricted to the Everglades region—while most polymorphisms were shared between lineages. The southern Florida group likely diverged from other populations during the Pliocene, with all other diversification during the Pleistocene. Multiple lines of support, including phylogenetic relationships, a latitudinal gradient in genetic diversity, and relatively more stable long‐term population sizes in southern phylogeographic groups, indicate that diversification in A. carolinensis occurred northward from southern Florida.

The complex evolutionary history of big-eared horseshoe bats (Rhinolophus macrotis complex): insights from genetic, morphological and acoustic data

Palaeoclimatic oscillations and different landscapes frequently result in complex population-level structure or the evolution of cryptic species. Elucidating the potential mechanisms is vital to understanding speciation events. However, such complex evolutionary patterns have rarely been reported in bats. In China, the Rhinolophus macrotis complex contains a large form and a small form, suggesting the existence of a cryptic bat species. Our field surveys found these two sibling species have a continuous and widespread distribution with partial sympatry. However, their evolutionary history has received little attention. Here, we used extensive sampling, morphological and acoustic data, as well as different genetic markers to investigate their evolutionary history. Genetic analyses revealed discordance between the mitochondrial and nuclear data. Mitochondrial data identified three reciprocally monophyletic lineages: one representing all small forms from Southwest China, and the other two containing all large forms from Central and Southeast China, respectively. The large form showed paraphyly with respect to the small form. However, clustering analyses of microsatellite and Chd1 gene sequences support two divergent clusters separating the large form and the small form. Moreover, morphological and acoustic analyses were consistent with nuclear data. This unusual pattern in the R. macrotis complex might be accounted for by palaeoclimatic oscillations, shared ancestral polymorphism and/or interspecific hybridization.

Expanding anchored hybrid enrichment to resolve both deep and shallow relationships within the spider tree of life

BACKGROUND

Despite considerable effort, progress in spider molecular systematics has lagged behind many other comparable arthropod groups, thereby hindering family-level resolution, classification, and testing of important macroevolutionary hypotheses. Recently, alternative targeted sequence capture techniques have provided molecular systematics a powerful tool for resolving relationships across the Tree of Life. One of these approaches, Anchored Hybrid Enrichment (AHE), is designed to recover hundreds of unique orthologous loci from across the genome, for resolving both shallow and deep-scale evolutionary relationships within non-model systems. Herein we present a modification of the AHE approach that expands its use for application in spiders, with a particular emphasis on the infraorder Mygalomorphae.

RESULTS

Our aim was to design a set of probes that effectively capture loci informative at a diversity of phylogenetic timescales. Following identification of putative arthropod-wide loci, we utilized homologous transcriptome sequences from 17 species across all spiders to identify exon boundaries. Conserved regions with variable flanking regions were then sought across the tick genome, three published araneomorph spider genomes, and raw genomic reads of two mygalomorph taxa. Following development of the 585 target loci in the Spider Probe Kit, we applied AHE across three taxonomic depths to evaluate performance: deep-level spider family relationships (33 taxa, 327 loci); family and generic relationships within the mygalomorph family Euctenizidae (25 taxa, 403 loci); and species relationships in the North American tarantula genus Aphonopelma (83 taxa, 581 loci). At the deepest level, all three major spider lineages (the Mesothelae, Mygalomorphae, and Araneomorphae) were supported with high bootstrap support. Strong support was also found throughout the Euctenizidae, including generic relationships within the family and species relationships within the genus Aptostichus. As in the Euctenizidae, virtually identical topologies were inferred with high support throughout Aphonopelma.

CONCLUSIONS

The Spider Probe Kit, the first implementation of AHE methodology in Class Arachnida, holds great promise for gathering the types and quantities of molecular data needed to accelerate an understanding of the spider Tree of Life by providing a mechanism whereby different researchers can confidently and effectively use the same loci for independent projects, yet allowing synthesis of data across independent research groups.

Genetic diversity and phylogeography of Siberian roe deer, Caproulus pygargus, in central and peripheral populations

Current understanding of phylogeographical structure and genetic diversity of Siberian roe deer remains limited mainly due to small sample size and/or low geographical coverage in previous studies. Published data suggest at least two phylogroups: western (Ural Mountains and Western Siberia) and eastern (east from lake Baikal, including the Korean peninsula), but their phylogenetic relationship remains unclear. Combined sequences of cytochrome b (1140 bp) and the mtDNA control region (963 bp) were analyzed from 219 Siberian roe deer from 12 locations in Russia, Mongolia, and South Korea, which cover a large part of its range, to assess genetic diversity and phylogeographical status. Special emphasis was placed on the demographic history and genetic features of central, peripheral, and isolated populations. Results of median-joining network and phylogenetic tree analyses indicate that Siberian roe deer from the Urals to the Pacific Ocean are genetically diverse and that geographical distribution and composition of haplogroups coincide with previously described ranges of the subspecies Capreolus pygargus pygargus and Capreolus pygargus tianschanicus. We found that peripheral populations in the northwestern parts of the species range (Urals), as well as the isolated population from Jeju Island, are genetically distinct from those in the core part of the range, both in terms of genetic diversity and quantitative composition of haplogroups. We also found that northwestern (Urals) and northern (Yakutia) peripheral populations share the same haplogroup and fall into the same phylogenetic clade with the isolated population from Jeju Island. This finding sheds light on the taxonomic status of the Jeju Island population and leads to hypotheses about the discordance of morphological and genetic evolution in isolated populations and specific genetic features of peripheral populations.

Microsatellite markers and cytoplasmic sequences reveal contrasting pattern of spatial genetic structure in the red algae species complex Mazzaella laminarioides

Mazzaella laminarioides is a common haploid-diploid red alga that forms dense beds. This alga has a wide distributional range, covering 3,500 km of the Chilean coast, but is restricted to high rocky intertidal zones. Recently, the existence of three highly divergent genetic lineages was demonstrated for this taxon, and two cytoplasmic markers were used to determine that these lineages are distributed in strict parapatry. Here, using 454 next-generation sequencing, we developed polymorphic microsatellite loci that cross amplify in all three cytoplasmic lineages. Six sites (i.e., two sites within each lineage) were analyzed using nine microsatellite loci. Our work shows that, although substantial cytoplasmic differentiation occurs within M. laminarioides, the microsatellite loci did not retrieve three nuclear genetic clusters as expected. Indeed, while the northernmost and southernmost cytoplasmic lineages form two strongly divergent nuclear groups characterized by diagnostic alleles, the third cytoplasmic lineage did not form a third nuclear independent group. It is possible that inter-lineage gene exchange has occurred, particularly at sites along the contact zone between the different cytoplasmic lineages. This nuclear-cytoplasmic incongruence in M. laminarioides could be explained by incomplete lineage sorting of the nuclear genes or asymmetric introgressive hybridization between the lineages. Finally, highly significant heterozygote deficiencies (suggesting occurrence of intergametophytic selfing) were observed in the three small northernmost sites while the large southernmost sites generally approached panmixia.

Importance of incomplete lineage sorting and introgression in the origin of shared genetic variation between two closely related pines with overlapping distributions

Genetic variation shared between closely related species may be due to retention of ancestral polymorphisms because of incomplete lineage sorting (ILS) and/or introgression following secondary contact. It is challenging to distinguish ILS and introgression because they generate similar patterns of shared genetic diversity, but this is nonetheless essential for inferring accurately the history of species with overlapping distributions. To address this issue, we sequenced 33 independent intron loci across the genome of two closely related pine species (Pinus massoniana Lamb. and Pinus hwangshanensis Hisa) from Southeast China. Population structure analyses revealed that the species showed slightly more admixture in parapatric populations than in allopatric populations. Levels of interspecific differentiation were lower in parapatry than in allopatry. Approximate Bayesian computation suggested that the most likely speciation scenario explaining this pattern was a long period of isolation followed by a secondary contact. Ecological niche modeling suggested that a gradual range expansion of P. hwangshanensis during the Pleistocene climatic oscillations could have been the cause of the overlap. Our study therefore suggests that secondary introgression, rather than ILS, explains most of the shared nuclear genomic variation between these two species and demonstrates the complementarity of population genetics and ecological niche modeling in understanding gene flow history. Finally, we discuss the importance of contrasting results from markers with different dynamics of migration, namely nuclear, chloroplast and mitochondrial DNA.

Cyto-nuclear discordance suggests complex evolutionary history in the cave-dwelling salamander, Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.

Combined hybridization and mitochondrial capture shape complex phylogeographic patterns in hybridogenetic Cataglyphis desert ants

Some species of Cataglyphis desert ants have evolved a hybridogenetic mode of reproduction at the social scale. In hybridogenetic populations, two distinct genetic lineages coexist. Non-reproductive offspring (workers) are hybrids of the two lineages, whereas sexual offspring (males and new queens) are produced by parthenogenesis and belong to the mother queen lineage. How this unusual reproductive system affects phylogeographic patterns and speciation processes remains completely unknown to date. Using one mitochondrial and four nuclear genes, we examined the phylogenetic relationships between three species of Cataglyphis (C. hispanica, C. humeya and C. velox) where complex DNA inheritance through social hybridogenesis may challenge phylogenetic inference. Our results bring two important insights. First, our data confirm a hybridogenetic mode of reproduction across the whole distribution range of the species C. hispanica. In contrast, they do not provide support for hybridogenesis in the populations sampled of C. humeya and C. velox. This suggests that these populations are not hybridogenetic, or that hybridogenesis is too recent to result in reciprocally monophyletic lineages on nuclear genes. Second, due to mitochondrial introgression between lineages (Darras and Aron, 2015), the faster-evolving COI marker is not lineage specific, hence, unsuitable to further investigate the segregation of lineages in the species studied. Different mitochondrial haplotypes occur in each locality sampled, resulting in strongly structured populations. This micro-allopatric structure leads to over-splitting species delimitation on mitochondrial gene, as every locality could potentially be considered a putative species; haploweb analyses of nuclear markers, however, yield species delimitations that are consistent with morphology. Overall, this study highlights how social hybridogenesis varies across species and shapes complex phylogeographic patterns.

An examination of species limits in the Aulacorhynchus "prasinus" toucanet complex (Aves: Ramphastidae)

The number of species recognized in Aulacorhynchus toucanets has varied tremendously over the past century. Revisors seem to disagree on whether head and bill coloration are useful indicators of species limits, especially in the A. "prasinus" complex. Using morphometrics, I tested the hypothesis that the major color-based subspecific groups of A. "prasinus" sensu lato are simply "cookie-cutter" (i.e., morphologically nearly identical) toucanets with different head and bill colorations. Univariate and multivariate analyses show that they are not simply morphological replicates of different colors: a complex array of morphometric similarities and dissimilarities occur between the major subspecific groups, and these variations differ between the sexes. Latitude and longitude had a small but significant association with female (but not male) PC1 and PC2. Hybridization and intergradation were also considered using plumage and bill characters as a surrogate to infer gene flow. Hybridization as indicated by phenotype appears to be substantial between A. "p." cyanolaemus and A. "p." atrogularis and nonexistent between other major groups, although from genetic evidence it is likely rare between A. "p." albivitta and A. "p." cyanolaemus. The congruence and complexities of the morphological and color changes occurring among these groups suggest that ecological adaptation (through natural selection) and social selection have co-occurred among these groups and that species limits are involved. Further, hybridization is not evident at key places, despite in many cases (hypothetical) opportunity for gene flow. Consequently, I recommend that this complex be recognized as comprising five biological species: A. wagleri, prasinus, caeruleogularis, albivitta, and atrogularis. Four of these also have valid subspecies within them, and additional work may eventually support elevation of some of these subspecies to full species. Species limits in South America especially need more study.

Phylogeographic Patterns in Africa and High Resolution Delineation of Genetic Clades in the Lion (Panthera leo)

Comparative phylogeography of African savannah mammals shows a congruent pattern in which populations in West/Central Africa are distinct from populations in East/Southern Africa. However, for the lion, all African populations are currently classified as a single subspecies (Panthera leo leo), while the only remaining population in Asia is considered to be distinct (Panthera leo persica). This distinction is disputed both by morphological and genetic data. In this study we introduce the lion as a model for African phylogeography. Analyses of mtDNA sequences reveal six supported clades and a strongly supported ancestral dichotomy with northern populations (West Africa, Central Africa, North Africa/Asia) on one branch, and southern populations (North East Africa, East/Southern Africa and South West Africa) on the other. We review taxonomies and phylogenies of other large savannah mammals, illustrating that similar clades are found in other species. The described phylogeographic pattern is considered in relation to large scale environmental changes in Africa over the past 300,000 years, attributable to climate. Refugial areas, predicted by climate envelope models, further confirm the observed pattern. We support the revision of current lion taxonomy, as recognition of a northern and a southern subspecies is more parsimonious with the evolutionary history of the lion.

High regional genetic diversity and lack of host-specificity in Ostrinia nubilalis (Lepidoptera: Crambidae) as revealed by mtDNA variation

The European corn borer (Ostrinia nubilalis) infests a wide array of host plants and is considered one of the most serious pests of maize in Europe. Recent studies suggest that individuals feeding on maize in Europe should be referred to O. nubilalis (sensu nov.), while those infesting dicots as Ostrinia scapulalis (sensu nov.). We test if the clear genetic distinctiveness among individuals of O. nubilalis living on maize vs. dicots is tracked by mitochondrial DNA (mtDNA). We used fragments of COI and COII genes of 32 individuals traditionally recognized as O. nubilalis collected on three host plants, maize, mugwort and hop, growing in different parts of Poland. In addition, we reconstructed the mtDNA phylogeny of Ostrinia species based on our data and sequences retrieved from GenBank to assess host and/or biogeographic patterns. We also compared haplotype variation found in Poland (east-central Europe) with other regions (Anatolia, Eastern Europe, Balkans, Far East, North America). Our study showed high mtDNA diversity of O. nubilalis in Poland in comparison with other regions and revealed rare haplotypes likely of Asian origin. We did not find distinct mtDNA haplotypes in larvae feeding on maize vs. dicotyledonous plants. Phylogenetic analyses showed an apparent lack of mtDNA divergence among putatively distinct lineages belonging to the O. nubilalis group as identical haplotypes are shared by Asian and European individuals. We argue that human-mediated dispersal, hybridization and sporadic host jumps are likely responsible for the lack of a geographic pattern in mtDNA variation.

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population-specific mitonuclear coadaptations create between-population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.