The biogeography of mitochondrial and nuclear discordance in animals

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male-biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal-biased X-chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.

Signatures of mito-nuclear discordance in Schistosoma turkestanicum indicate a complex evolutionary history of emergence in Europe

High levels of molecular diversity were identified in mitochondrial cytochrome c oxidase (cox1) gene sequences of Schistosoma turkestanicum from Hungary. These cox1 sequences were all specific to Hungary which contrasted with the low levels of diversity seen in the nuclear internal transcribed spacer region (ITS) sequences, the majority of which were shared between China and Iran isolates. Measures of within and between host molecular variation within S. turkestanicum showed there to be substantial differences in molecular diversity, with cox1 being significantly more diverse than the ITS. Measures of haplotype frequencies revealed that each host contained its own subpopulation of genetically unique parasites with significant levels of differentiation. Pairwise mismatch analysis of cox1 sequences indicated S. turkestanicum populations to have a bimodal pairwise difference distribution and to be stable unlike the ITS sequences, which appeared to have undergone a recent population expansion event. Positive selection was also detected in the cox1 sequences, and biochemical modelling of the resulting protein illustrated significant mutational events causing an alteration to the isoelectric point of the cox1 protein, potentially altering metabolism. The evolutionary signature from the cox1 indicates local adaptation and long establishment of S. turkestanicum in Hungary with continual introgression of nuclear genes from Asian isolates. These processes have led to the occurrence of mito-nuclear discordance in a schistosome population.

Snails in the desert: Species diversification of Theba (Gastropoda: Helicidae) along the Atlantic coast of NW Africa

The spatial subdivision of species often plays a pivotal role in speciation. Across their entire range, species are rarely panmictic and crucial consequences of spatial subdivision are (1) random genetic drift including historical factors, (2) uniform selection, and (3) divergent selection. Each of these consequences may result in geographic variation and eventually reproductive isolation, but their relative importance in speciation is still unclear. In this study, we used a combination of genetic, morphological, and climatic data to obtain a comprehensive picture of differentiation among three closely related, parapatrically distributed taxa of the land snail genus Theba occurring along the Atlantic coasts of South Morocco and Western Sahara. We conducted Mantel and partial Mantel tests to relate phenotypic and genotypic variation of these species to geography and/or climate. As null hypothesis for an evolutionary scenario, we assumed nonadaptive speciation and expected a pattern of isolation by distance among taxa. Rejection of the null hypothesis would indicate isolation by environment due to adaptation. Generally, genetic drift plays an important role but is rarely considered as sole driver of speciation. It is the combination of drift and selection that predominantly drives speciation. This study, however, provides a potential example, in which nonadaptive speciation, that is, genetic drift, is apparently the main driver of shaping the diversity of Theba in NW Africa. Restriction of gene flow between populations caused by geographic isolation probably has played an important role. Climate oscillations during the Plio- and Pleistocene may have led to repeated ecological changes in NW Africa and disruptions of habitats promoting differentiation by geographic isolation. The inferred evolutionary scenario, however, did not fully explain the incongruence between the AFLP- and mtDNA-tree topologies. This incongruence might indicate past hybridization among the studied Theba forms.

Is Galba schirazensis (Mollusca, Gastropoda) an intermediate host of Fasciola hepatica (Trematoda, Digenea) in Ecuador?

Fasciolosis is a widely distributed disease in livestock in South America but knowledge about the epidemiology and the intermediate hosts is relatively scarce in Ecuador. For three months, lymnaeid snails were sampled (n = 1482) in Pichincha Province at two sites located in a highly endemic area. Snails were identified (based on morphology and ITS-2 sequences) and the infection status was established through microscopic dissection and a multiplex polymerase chain reaction (PCR)-based technique. Techniques based on morphology were not useful to accurately name the collected snail species. Comparison with available DNA sequences showed that a single snail species was collected, Galba schirazensis. Live rediae were observed in 1.75% (26/1482) and Fasciola sp. DNA was detected in 6% (89/1482) of collected snails. The COX-1 region permitted identification of the parasite as Fasciola hepatica. The relative sensitivity and specificity of the microscope study, compared to PCR results, were 25.84% and 99.78%, respectively. The mean size of the snails recorded positive for F. hepatica through crushing and microscopy was significantly higher than the mean size of negative snails, but there was no such difference in PCR-positive snails. The role of G. schirazensis as an intermediate host of F. hepatica in Ecuador is discussed and the hypothesis of an adaptation of the parasite to this invasive snail is proposed. For the first time, an epidemiological survey based on molecular biology-based techniques assessed the possible role of lymnaeid snails in the epidemiology of fasciolosis in Ecuador.

Human-mediated secondary contact of two tortoise lineages results in sex-biased introgression

Human-mediated secondary contact of recently diverged taxa offers valuable opportunities for studying the evolutionary mechanisms involved in the establishment and maintenance of genetic boundaries between taxa. We used mitochondrial and microsatellite markers to examine a recently introduced population of the spur-thighed tortoise (Testudo graeca) of mixed origin in the Doñana National Park (SW Spain). The earliest records of tortoises in Doñana trace back to the 18th century, but several population reinforcements in the 20th century with animals from Morocco are well-documented. Consequently, different genetic lineages, which represent distinct subspecies, are thought to co-exist there. Our results confirmed the presence of distinct lineages by revealing that tortoises of the subspecies T. g. marokkensis were introduced into a local allochthonous T. g. graeca population. Unexpectedly, T. g. marokkensis haplotypes exclusively appeared in males, and admixture levels were statistically sex-biased toward males. The sex ratio of the population deviated from parity, with males being 2.36-fold more abundant than females. Our results indicated that population reinforcements had a strong effect on the genetic composition of this population and aggravated its sex ratio deviation. We predict that this sex-biased pattern of introgression is ephemeral and advocated to the near loss of T. g. marokkensis haplotypes.

Geographically distinct patterns of reproductive isolation and hybridization in two sympatric species of the Jaera albifrons complex (marine isopods)

Sympatric species that in some populations hybridize and in other populations remain reproductively isolated open interesting research possibilities for the study of hybridization and speciation. Here, we test for such a situation in two littoral isopods (Jaera albifrons and J. praehirsuta) that occur in mixed populations and where past morphological descriptions suggested that the two species are generally reproductively isolated except in rare populations where hybridization may be happening. Using field surveys and microsatellite genetic structure analyses in two regions from France (Normandy and Brittany), we confirmed that introgressive hybridization occurs in a subset of mixed J. albifrons/J. praehirsuta populations (region Normandy) where the two species are found in the same habitat (pebbles on the shore). Moreover, we found that introgression in these populations is differential, 21 of 23 microsatellite markers showing little genetic divergence between species (hierarchical analysis of molecular variance F_CT = 0.017) while the remaining two loci were strongly differentiated (F_CT = 0.428). By contrast, J. albifrons and J. praehirsuta in mixed populations from region Brittany occupied distinct habitats (pebbles and seaweeds, respectively) with little overlap and showed stronger genetic divergence (F_CT = 0.132). In hybridizing populations, the majority of individuals show morphological traits that are characteristic of one or the other species. This raises the question of the forces that act to maintain this polymorphism, noting that hybridizing populations seem to be geographically isolated from potential source parental populations and show no detectable habitat divergence between species.

Quaternary history, population genetic structure and diversity of the cold-adapted Alpine newt Ichthyosaura alpestris in peninsular Italy

Mediterranean peninsulas are major biodiversity hotspots, and cold-adapted species are an important component of this biodiversity. However, cold-adapted species contributed surprisingly little to our knowledge of the intimate links between Quaternary environmental changes, species' responses to these changes, and current patterns of intraspecific biodiversity. Here, we investigated the genetic structure and evolutionary history of a cold-adapted amphibian, the Alpine newt Ichthyosaura alpestris, within the Italian peninsula. Nuclear and mitochondrial markers consistently identified three distinct genetic lineages, whose divergence dates to the Early Pleistocene (1.9 and 0.8 million years ago). Our results show that the Italian peninsula provided multiple Pleistocene refugia to this cold-adapted species, and suggest that allopatric fragmentation followed by secondary admixture have been key events in the formation of its current pattern of genetic diversity. Indeed, estimates of population genetic diversity clearly identified contact populations as those achieving the highest levels of diversity. Such concordance among cold-adapted and temperate species in terms of processes triggering the formation of regional patterns of genetic diversity provides strong support for the hypothesis that gene exchange between divergent lineages, rather than long-term stability of refugial populations, has been the main step toward the formation of hotspots of intraspecific biodiversity.

New insight into hybridization and unidirectional introgression between Ammodytes japonicus and Ammodytes heian (Trachiniformes, Ammodytidae)

Based on northern (NOL) and southern (SOL) mitochondrial lineages, recently, it proposed the new species Ammodytes heian and revived the species name Ammodytes japonicus to describe sand lances from the northwestern Pacific Ocean. This study used molecular methods to investigate genetic relationships between the two sand lance species in Korea and Japan. In total, 154 specimens were collected from four locations in Korea (Baengnyeongdo in the Yellow Sea, Tongyeong in the Korean Strait, and Jumunjin and Gijang in the East Sea), and 50 specimens were collected from a single location in Japan (Wakkanai in the Okhotsk Sea). Mitochondrial DNA analysis demonstrated that the individuals from Baengnyeongdo and Tongyeong all belonged to the SOL, whereas those from Gijang, Jumunjin, and Wakkanai included individuals from both the NOL and SOL (over 75% NOL). Population structure analyses were performed on the same individuals using seven microsatellite DNA markers. The population structure analysis based on 201 specimens identified two clusters (named as northern group and southern group), with the admixture proportion (q) of < 0.1 for the northern group in the Backyeongdo and Tongyeong sand lances and < 0.1 for the southern group in the Wakkanai sand lances. The high heterogeneity indicated that the former was probably A. japonicus and the latter probably A. heian. However, the admixture proportion in the Jumunjin and Gijang sand lances was 0.71–0.75 for the southern group, indicating that hybridization and unidirectional introgression from SOL to NOL occurs in southwestern margin of the East Sea. Our findings illustrate the speciation process based on different patterns of gene flow between Korean and Japanese sand lance, which is strongly influenced by both the paleo-climatic change and the contemporary local oceanic current pattern.

Squamate Conserved Loci (SqCL): A unified set of conserved loci for phylogenomics and population genetics of squamate reptiles

The identification of conserved loci across genomes, along with advances in target capture methods and high-throughput sequencing, has helped spur a phylogenomics revolution by enabling researchers to gather large numbers of homologous loci across clades of interest with minimal upfront investment in locus design. Target capture for vertebrate animals is currently dominated by two approaches-anchored hybrid enrichment (AHE) and ultraconserved elements (UCE)-and both approaches have proven useful for addressing questions in phylogenomics, phylogeography and population genomics. However, these two sets of loci have minimal overlap with each other; moreover, they do not include many traditional loci that that have been used for phylogenetics. Here, we combine across UCE, AHE and traditional phylogenetic gene locus sets to generate the Squamate Conserved Loci set, a single integrated probe set that can generate high-quality and highly complete data across all three loci types. We use these probes to generate data for 44 phylogenetically disparate taxa that collectively span approximately 33% of terrestrial vertebrate diversity. Our results generated an average of 4.29 Mb across 4709 loci per individual, of which an average of 2.99 Mb was sequenced to high enough coverage (≥10×) to use for population genetic analyses. We validate the utility of these loci for both phylogenomic and population genomic questions, provide a comparison among these locus sets of their relative usefulness and suggest areas for future improvement.

A genomic footprint of hybrid zone movement in crested newts

Speciation typically involves a stage in which species can still exchange genetic material. Interspecific gene flow is facilitated by the hybrid zones that such species establish upon secondary contact. If one member of a hybridizing species pair displaces the other, their hybrid zone would move across the landscape. Although theory predicts that moving hybrid zones quickly stagnate, hybrid zones tracked over one or a few decades do not always follow such a limitation. This suggests that hybrid zones have the potential to traverse considerable distances over extended periods of time. When hybrid zones move, introgression is predicted to result in biased gene flow of selectively neutral alleles, from the receding species into the advancing species. We test for such a genomic footprint of hybrid zone movement in a pair of crested newt species (genus Triturus) for which we have a priori support for westward hybrid zone movement. We perform a multilocus phylogeographical survey and conduct Bayesian clustering analysis, estimation of ancestry and heterozygosity, and geographical cline analysis. In a 600 km wide area east of the present day hybrid zone a genomic footprint constitutes empirical evidence consistent with westward hybrid zone movement. The crested newt case suggests that hybrid zone movement can occur over an extensive span of time and space. Inferring hybrid zone movement provides fundamental insight into historical biogeography and the speciation process, and we anticipate that hybrid zones will prove to be far more mobile than currently appreciated.

Historical mitochondrial diversity in African leopards (Panthera pardus) revealed by archival museum specimens

Once found throughout Africa and Eurasia, the leopard (Panthera pardus) was recently uplisted from Near Threatened to Vulnerable by the International Union for the Conservation of Nature (IUCN). Historically, more than 50% of the leopard's global range occurred in continental Africa, yet sampling from this part of the species' distribution is only sparsely represented in prior studies examining patterns of genetic variation at the continental or global level. Broad sampling to determine baseline patterns of genetic variation throughout the leopard's historical distribution is important, as these measures are currently used by the IUCN to direct conservation priorities and management plans. By including data from 182 historical museum specimens, faecal samples from ongoing field surveys, and published sequences representing sub-Saharan Africa, we identify previously unrecognized genetic diversity in African leopards. Our mtDNA data indicates high levels of divergence among regional populations and strongly differentiated lineages in West Africa on par with recent studies of other large vertebrates. We provide a reference benchmark of genetic diversity in African leopards against which future monitoring can be compared. These findings emphasize the utility of historical museum collections in understanding the processes that shape present biodiversity. Additionally, we suggest future research to clarify African leopard taxonomy and to differentiate between delineated units requiring monitoring or conservation action.

The on-again, off-again relationship between mitochondrial genomes and species boundaries

The study of reproductive isolation and species barriers frequently focuses on mitochondrial genomes and has produced two alternative and almost diametrically opposed narratives. On one hand, mtDNA may be at the forefront of speciation events, with co-evolved mitonuclear interactions responsible for some of the earliest genetic incompatibilities arising among isolated populations. On the other hand, there are numerous cases of introgression of mtDNA across species boundaries even when nuclear gene flow is restricted. We argue that these seemingly contradictory patterns can result from a single underlying cause. Specifically, the accumulation of deleterious mutations in mtDNA creates a problem with two alternative evolutionary solutions. In some cases, compensatory or epistatic changes in the nuclear genome may ameliorate the effects of mitochondrial mutations, thereby establishing coadapted mitonuclear genotypes within populations and forming the basis of reproductive incompatibilities between populations. Alternatively, populations with high mitochondrial mutation loads may be rescued by replacement with a more fit, foreign mitochondrial haplotype. Coupled with many nonadaptive mechanisms of introgression that can preferentially affect cytoplasmic genomes, this form of adaptive introgression may contribute to the widespread discordance between mitochondrial and nuclear genealogies. Here, we review recent advances related to mitochondrial introgression and mitonuclear incompatibilities, including the potential for cointrogression of mtDNA and interacting nuclear genes. We also address an emerging controversy over the classic assumption that selection on mitochondrial genomes is inefficient and discuss the mechanisms that lead lineages down alternative evolutionary paths in response to mitochondrial mutation accumulation.

Correlated patterns of genetic diversity and differentiation across an avian family

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole-genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.

Genetic structure among greater white-fronted goose populations of the Pacific Flyway

An understanding of the genetic structure of populations in the wild is essential for long‐term conservation and stewardship in the face of environmental change. Knowledge of the present‐day distribution of genetic lineages (phylogeography) of a species is especially important for organisms that are exploited or utilize habitats that may be jeopardized by human intervention, including climate change. Here, we describe mitochondrial (mtDNA) and nuclear genetic (microsatellite) diversity among three populations of a migratory bird, the greater white‐fronted goose (Anser albifrons), which breeds discontinuously in western and southwestern Alaska and winters in the Pacific Flyway of North America. Significant genetic structure was evident at both marker types. All three populations were differentiated for mtDNA, whereas microsatellite analysis only differentiated geese from the Cook Inlet Basin. In sexual reproducing species, nonrandom mate selection, when occurring in concert with fine‐scale resource partitioning, can lead to phenotypic and genetic divergence as we observed in our study. If mate selection does not occur at the time of reproduction, which is not uncommon in long‐lived organisms, then mechanisms influencing the true availability of potential mates may be obscured, and the degree of genetic and phenotypic diversity may appear incongruous with presumed patterns of gene flow. Previous investigations revealed population‐specific behavioral, temporal, and spatial mechanisms that likely influence the amount of gene flow measured among greater white‐fronted goose populations. The degree of observed genetic structuring aligns well with our current understanding of population differences pertaining to seasonal movements, social structure, pairing behavior, and resource partitioning.

Frequent discordance between morphology and mitochondrial DNA in a species group of European water beetles (Coleoptera: Dytiscidae)

The Hydroporus memnonius species group includes both widespread and range restricted diving beetle taxa in the western Palaearctic, some of which have been divided into a number of geographical subspecies. Of these, Hydroporus necopinatus is distributed in the far west of Europe, from central Spain to southern Britain, and has been split into three subspecies, occurring in Iberia (necopinatus sst.), France (robertorum) and England (roni) respectively—the last of these being a rare example of an insect taxon apparently endemic to northern Europe. Here we explore inter-relationships between populations and subspecies of H. necopinatus and related members of the Hydroporus melanarius subgroup, using mitochondrial COI sequence data. We reveal widespread discordance between mitochondrial DNA sequence variation and morphology in areas where H. necopinatus and H. melanarius come into contact, consistent with historical introgressive hybridization between these taxa. In light of this discordance, the lack of clear genetic divergence between H. necopinatus subspecies, and the fact that both robertorum and roni are morphologically intermediate between H. necopinatus sstr. and H. melanarius, we suggest that these taxa may be of hybridogenic origin, rather than representing discrete evolutionary lineages.

The biogeographic and evolutionary history of an endemic clade of Middle American sparrows: Melozone and Aimophila (Aves: Passerellidae)

The large number of endemic species in Middle America is frequently attributed to the interplay of geographical barriers and historical climatic changes in the region. This process promotes genetic divergence between populations, and given enough time, may yield new species. Animals that inhabit mid-elevation or highland habitats may be disproportionately affected in this way. Genetic analyses of animals in this region allow us to better understand how historical patterns of isolation have influenced the generation of new species in this biodiversity hotspot. We studied the biogeography and systematics of two closely related genera of sparrows (Passerellidae): Melozone and Aimophila. Collectively, this group is distributed from the southwestern United States and southward as far as central Costa Rica. We sampled 81 individuals of 8 Melozone and 2 Aimophila species, from 19 localities distributed throughout their ranges. We reconstructed phylogenetic relationships and time-calibrated species trees using multilocus sequence data comprised of one mitochondrial gene and five nuclear genes. We conducted an ancestral area reconstruction analysis to determine the probability of ancestral range at each divergent event. Despite analyzing six loci, we were unable to obtain a fully resolved phylogenetic tree. We recovered four main lineages: lineage 1 includes four Melozone species distributed north of Isthmus of Tehuantepec (M. albicollis, M. crissalis, M. aberti, M. fusca); lineage 2 includes three Melozone species distributed south of the Isthmus of Tehuantepec (M. biarcuata, M. cabanisi, M. leucotis); lineage 3 lineage consists of a single species endemic to the Pacific coast of Mexico (M. kieneri); and lineage 4 includes the more widely distributed sparrows in the genus Aimophila. Our analyses suggest that these genera probably originated during the late Miocene in the Madrean Highlands of southern Mexico. We identified dispersal as the prevalent cause of speciation in this clade with most lineages dispersing to their current distributions from southern Mexico either to the north following a developing and expanding Madro-Tertiary flora, or to the south across the Isthmus of Tehuantepec. A similar pattern of dispersal from this biogeographic region has been reported in other taxa including fishes, reptiles, and birds. Our results reveal that the four lineages identified represent geographically coherent and ecologically similar assemblages of taxa. Finally, when our genetic results are considered, along with apparent differences in morphology and song, the allopatric forms M. b. cabanisi and M. l. occipitalis warrant recognition as biological species.

Integrative Approaches for Studying Mitochondrial and Nuclear Genome Co-evolution in Oxidative Phosphorylation

In animals, interactions among gene products of mitochondrial and nuclear genomes (mitonuclear interactions) are of profound fitness, evolutionary, and ecological significance. Most fundamentally, the oxidative phosphorylation (OXPHOS) complexes responsible for cellular bioenergetics are formed by the direct interactions of 13 mitochondrial-encoded and ∼80 nuclear-encoded protein subunits in most animals. It is expected that organisms will develop genomic architecture that facilitates co-adaptation of these mitonuclear interactions and enhances biochemical efficiency of OXPHOS complexes. In this perspective, we present principles and approaches to understanding the co-evolution of these interactions, with a novel focus on how genomic architecture might facilitate it. We advocate that recent interdisciplinary advances assist in the consolidation of links between genotype and phenotype. For example, advances in genomics allow us to unravel signatures of selection in mitochondrial and nuclear OXPHOS genes at population-relevant scales, while newly published complete atomic-resolution structures of the OXPHOS machinery enable more robust predictions of how these genes interact epistatically and co-evolutionarily. We use three case studies to show how integrative approaches have improved the understanding of mitonuclear interactions in OXPHOS, namely those driving high-altitude adaptation in bar-headed geese, allopatric population divergence in Tigriopus californicus copepods, and the genome architecture of nuclear genes coding for mitochondrial functions in the eastern yellow robin.

Digging up the roots of an insular hotspot of genetic diversity: decoupled mito-nuclear histories in the evolution of the Corsican-Sardinian endemic lizard Podarcis tiliguerta

BACKGROUND

Mediterranean islands host a disproportionately high level of biodiversity and endemisms. Growing phylogeographic evidence on island endemics has unveiled unexpectedly complex patterns of intra-island diversification, which originated at diverse spatial and temporal scales. We investigated multilocus genetic variation of the Corsican-Sardinian endemic lizard Podarcis tiliguerta with the aim of shedding more light on the evolutionary processes underlying the origin of Mediterranean island biodiversity.

RESULTS

We analysed DNA sequences of mitochondrial (12S and nd4) and nuclear (acm4 and mc1r) gene fragments in 174 individuals of P. tiliguerta from 81 localities across the full range of the species in a geographic and genealogical framework. We found surprisingly high genetic diversity both at mitochondrial and nuclear loci. Seventeen reciprocally monophyletic allopatric mitochondrial haplogroups were sharply divided into four main mitochondrial lineages (two in Corsica and two in Sardinia) of Miocene origin. In contrast, shallow divergence and shared diversity within and between islands was observed at the nuclear loci. We evaluated alternative biogeographic and evolutionary scenarios to explain such profound discordance in spatial and phylogenetic patterning between mitochondrial and nuclear genomes. While neutral models provided unparsimonious explanations for the observed pattern, the hypothesis of environmental selection driving mitochondrial divergence in the presence of nuclear gene flow is favoured.

CONCLUSIONS

Our study on the genetic variation of P. tiliguerta reveals surprising levels of diversity underlining a complex phylogeographic pattern with a striking example of mito-nuclear discordance. These findings have profound implications, not only for the taxonomy and conservation of P. tiliguerta. Growing evidence on deep mitochondrial breaks in absence of geographic barriers and of climatic factors associated to genetic variation of Corsican-Sardinian endemics warrants additional investigation on the potential role of environmental selection driving the evolution of diversity hotspots within Mediterranean islands.

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.