A reassessment of explanations for discordant introgressions of mitochondrial and nuclear genomes

Phylogenomic test of mitochondrial clues to archaic ancestors in a group of hybridizing swallowtail butterflies

Genomics has revolutionized our understanding of hybridization and introgression, but most of the early evidence for these processes came from studies of mitochondrial introgression. To expand these evolutionary insights from mitochondrial patterns, we evaluate phylogenetic discordance across the nuclear genomes of a hybridizing system, the Papilio machaon group of swallowtail butterflies. This species group contains three hybrid lineages (P. brevicauda, P. joanae, and P. m. kahli) that are geographically disjunct across North America and have complete fixation of a mitochondrial lineage that is otherwise primarily found in P. m. hudsonianus, a boreal subspecies of the Holarctic P. machaon. Genome-wide nuclear markers place the three hybrid lineages as a monophyletic group that is sister to P. polyxenes/P. zelicaon rather than P. machaon, although ancient hybridization between a subspecies of P. machaon and the ancestor of these three lineages is also shown by their greater nuclear affinity to P. m. hudsonianus than to other subspecies of P. machaon. Individuals from contemporary hybrid swarms in Alberta, where mitochondrial DNA fixation has not occurred, were more intermediate between their respective parent species, demonstrating diversity in mito-nuclear discordance following hybrid interactions. Our new phylogenetic findings for the P. machaon species group also include: subspecific paraphyly within P. machaon itself across its Holarctic distribution; paraphyly of P. zelicaon relative to P. polyxenes; and more divergent placement of a Mediterranean species, P. hospiton. These results provide the first comprehensive genomic evaluation of relationships within this species group and provide insight into the evolutionary dynamics of hybridization and mitochondrial introgression.

Successful mating and hybridisation in two closely related flatworm species despite significant differences in reproductive morphology and behaviour

Reproductive traits are some of the fastest diverging characters and can serve as reproductive barriers. The free-living flatworm Macrostomum lignano, and its congener M. janickei are closely related, but differ substantially in their male intromittent organ (stylet) morphology. Here, we examine whether these morphological differences are accompanied by differences in behavioural traits, and whether these could represent barriers to successful mating and hybridization between the two species. Our data shows that the two species differ in many aspects of their mating behaviour. Despite these differences, the species mate readily with each other in heterospecific pairings. Although both species have similar fecundity in conspecific pairings, the heterospecific pairings revealed clear postmating barriers, as few heterospecific pairings produced F1 hybrids. These hybrids had a stylet morphology that was intermediate between that of the parental species, and they were fertile. Finally, using a mate choice experiment, we show that the nearly two-fold higher mating rate of M. lignano caused it to mate more with conspecifics, leading to assortative mating, while M. janickei ended up mating more with heterospecifics. Thus, while the two species can hybridize, the mating rate differences could possibly lead to higher fitness costs for M. janickei compared to M. lignano.

Finding complexity in complexes: Assessing the causes of mitonuclear discordance in a problematic species complex of Mesoamerican toads

Mitonuclear discordance is a frequently encountered pattern in phylogeographic studies and occurs when mitochondrial and nuclear DNA display conflicting signals. Discordance among these genetic markers can be caused by several factors including confounded taxonomies, gene flow, and incomplete lineage sorting. In this study, we present a strong case of mitonuclear discordance in a species complex of toads (Bufonidae: Incilius coccifer complex) found in the Chortís Block of Central America. To determine the cause of mitonuclear discordance in this complex, we used spatially explicit genetic data to test species limits and relationships, characterize demographic history, and quantify gene flow. We found extensive mitonuclear discordance among the three recognized species within this group, especially in populations within the Chortís Highlands of Honduras. Our data reveal nuclear introgression within the Chortís Highlands populations that was most probably driven by cyclical range expansions due to climatic fluctuations. Though we determined introgression occurred within the nuclear genome, our data suggest that it is not the key factor in driving mitonuclear discordance in the entire species complex. Rather, due to a lack of discernible geographic pattern between mitochondrial and nuclear DNA, as well as a relatively recent divergence time of this complex, we concluded that mitonuclear discordance has been caused by incomplete lineage sorting. Our study provides a framework to test sources of mitonuclear discordance and highlights the importance of using multiple marker types to test species boundaries in cryptic species.

Robbery in progress: Historical museum collections bring to light a mitochondrial capture within a bird species widespread across southern Australia, the Copperback Quail-thrush Cinclosoma clarum

We surveyed mitochondrial, autosomal, and Z chromosome diversity within and between the Copperback Quail-thrush Cinclosoma clarum and Chestnut Quail-thrush C. castanotum, which together span the arid and semi-arid zones of southern Australia, and primarily from specimens held in museum collections. We affirm the recent taxonomic separation of the two species and then focus on diversity within the more widespread of the two species, C. clarum. To guide further study of the system and what it offers to understanding the genomics of the differentiation and speciation processes, we develop and present a hypothesis to explain mitonuclear discordance that emerged in ourdata. Following a period of historical allopatry, secondary contact has resulted in an eastern mitochondrial genome replacing the western mitochondrial genome in western populations. This is predicted under a population-level invasion in the opposite direction, that of the western population invading the range of the eastern one. Mitochondrial captures can be driven by neutral, demographic processes, or adaptive mechanisms, and we favor the hypothesized capture being driven by neutral means. We cannot fully reject the adaptive process but suggest how these alternatives may be further tested. We acknowledge an alternative hypothesis, which finds some support in phenotypic data published elsewhere, namely that outcomes of secondary contact have been more complex than our current genomic data suggest. Discriminating and reconciling these two alternative hypotheses, which may not be mutually exclusive, could be tested with closer sampling at levels of population, individual, and nucleotide than has so far been possible. This would be further aided by knowledge of the genetic basis to phenotypic variation described elsewhere.

Genome-wide data reveal discordant mitonuclear introgression in the intermediate horseshoe bat (Rhinolophus affinis)

Closely related taxa often exhibit mitonuclear discordance attributed to introgression of mitochondrial DNA (mtDNA), yet few studies have considered the underlying causes of mtDNA introgression. Here we test for demographic versus adaptive processes as explanations for mtDNA introgression in three subspecies of the intermediate horseshoe bat (Rhinolophus affinis). We generated sequences of 1692 nuclear genes and 13 mitochondrial protein-coding genes for 48 individuals. Phylogenetic reconstructions based on 320 exon sequences and 2217 single nucleotide polymorphisms (SNPs) both revealed conflicts between the species tree and mtDNA tree. These results, together with geographic patterns of mitonuclear discordance, and shared identical or near-identical mtDNA sequences, suggest extensive introgression of mtDNA between the two parapatric mainland subspecies. Under demographic hypotheses, we would also expect to uncover traces of ncDNA introgression, however, population structure and gene flow analyses revealed little nuclear admixture. Furthermore, we found inconsistent estimates of the timing of population expansion and that of the most recent common ancestor for the clade containing introgressed haplotypes. Without a clear demographic explanation, we also examined whether introgression likely arises from adaptation. We found that five mtDNA genes contained fixed amino acid differences between introgressed and non-introgressed individuals, including putative positive selection found in one codon, although this did not show introgression. While our evidence for rejecting demographic hypotheses is arguably stronger than that for rejecting adaptation, we find no definitive support for either explanation. Future efforts will focus on larger-scale resequencing to decipher the underlying causes of discordant mitonuclear introgression in this system.

MitoFinder: Efficient automated large-scale extraction of mitogenomic data in target enrichment phylogenomics

Thanks to the development of high-throughput sequencing technologies, target enrichment sequencing of nuclear ultraconserved DNA elements (UCEs) now allows routine inference of phylogenetic relationships from thousands of genomic markers. Recently, it has been shown that mitochondrial DNA (mtDNA) is frequently sequenced alongside the targeted loci in such capture experiments. Despite its broad evolutionary interest, mtDNA is rarely assembled and used in conjunction with nuclear markers in capture-based studies. Here, we developed MitoFinder, a user-friendly bioinformatic pipeline, to efficiently assemble and annotate mitogenomic data from hundreds of UCE libraries. As a case study, we used ants (Formicidae) for which 501 UCE libraries have been sequenced whereas only 29 mitogenomes are available. We compared the efficiency of four different assemblers (IDBA-UD, MEGAHIT, MetaSPAdes, and Trinity) for assembling both UCE and mtDNA loci. Using MitoFinder, we show that metagenomic assemblers, in particular MetaSPAdes, are well suited to assemble both UCEs and mtDNA. Mitogenomic signal was successfully extracted from all 501 UCE libraries, allowing us to confirm species identification using CO1 barcoding. Moreover, our automated procedure retrieved 296 cases in which the mitochondrial genome was assembled in a single contig, thus increasing the number of available ant mitogenomes by an order of magnitude. By utilizing the power of metagenomic assemblers, MitoFinder provides an efficient tool to extract complementary mitogenomic data from UCE libraries, allowing testing for potential mitonuclear discordance. Our approach is potentially applicable to other sequence capture methods, transcriptomic data and whole genome shotgun sequencing in diverse taxa. The MitoFinder software is available from GitHub (https://github.com/RemiAllio/MitoFinder).

Variable outcomes of hybridization between declining Alosa alosa and Alosa fallax

Hybridization dynamics between co-occurring species in environments where human-mediated changes take place are important to quantify for furthering our understanding of human impacts on species evolution and for informing management. The allis shad Alosa alosa (Linnaeus, 1758) and twaite shad Alosa fallax (Lacépède, 1803), two clupeids sister species, have been severely impacted by human activities across Europe. The shrinkage of A. alosa distribution range along with the decline of the remaining populations' abundance threatens its persistence. The main objective was to evaluate the extent of hybridization and introgression between those interacting species. We developed a set of 77 species-specific SNP loci that allowed a better resolution than morphological traits as they enabled the detection of hybrids up to the third generation. Variable rates of contemporary hybridization and introgression patterns were detected in 12 studied sites across the French Atlantic coast. Mitochondrial markers revealed a cyto-nuclear discordance almost invariably involving A. alosa individuals with an A. fallax mitochondrial DNA and provided evidence of historical asymmetric introgression. Overall, contemporary and historical introgression revealed by nuclear and mitochondrial markers strongly suggests that a transfer of genes occurs from A. fallax toward A. alosa genome since at least four generations. Moreover, the outcomes of introgression greatly depend on the catchments where local processes are thought to occur. Undoubtedly, interspecific interaction and gene flow should not be overlooked when considering the management of those species.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

'Mainland-island' population structure of a terrestrial salamander in a forest-bocage landscape with little evidence for in situ ecological speciation

Adaptation to different ecological environments can, through divergent selection, generate phenotypic and genetic differences between populations, and eventually give rise to new species. The fire salamander (Salamandra salamandra) has been proposed to represent an early stage of ecological speciation, driven by differential habitat adaptation through the deposition and development of larvae in streams versus ponds in the Kottenforst near Bonn (Germany). We set out to test this hypothesis of ecological speciation in an area different from the one where it was raised and we took the opportunity to explore for drivers of genetic differentiation at a landscape scale. A survey over 640 localities demonstrated the species' presence in ponds and streams across forests, hilly terrain and areas with hedgerows ('bocage'). Genetic variation at 14 microsatellite loci across 41 localities in and around two small deciduous forests showed that salamander effective population sizes were higher in forests than in the bocage, with panmixia in the forests (Fst < 0.010) versus genetic drift or founder effects in several of the small and more or less isolated bocage populations (Fst > 0.025). The system fits the 'mainland-island' metapopulation model rather than indicating adaptive genetic divergence in pond versus stream larval habitats. A reanalysis of the Kottenforst data indicated that microsatellite genetic variation fitted a geographical rather than an environmental axis, with a sharp transition from a western pond-breeding to an eastern, more frequently stream-breeding group of populations. A parallel changeover in mitochondrial DNA exists but remains to be well documented. The data support the existence of a hybrid zone following secondary contact of differentiated lineages, more so than speciation in situ.

Genome-scale data reveal the role of hybridization in lichen-forming fungi

Advancements in molecular genetics have revealed that hybridization may be common among plants, animals, and fungi, playing a role in evolutionary dynamics and speciation. While hybridization has been well-documented in pathogenic fungi, the effects of these processes on speciation in fungal lineages with different life histories and ecological niches are largely unexplored. Here we investigated the potential influence of hybridization on the emergence of morphologically and reproductively distinct asexual lichens. We focused on vagrant forms (growing obligately unattached to substrates) within a clade of rock-dwelling, sexually reproducing species in the Rhizoplaca melanophthalma (Lecanoraceae, Ascomycota) species complex. We used phylogenomic data from both mitochondrial and nuclear genomes to infer evolutionary relationships and potential patterns of introgression. We observed multiple instances of discordance between the mitochondrial and nuclear trees, including the clade comprising the asexual vagrant species R. arbuscula, R. haydenii, R. idahoensis, and a closely related rock-dwelling lineage. Despite well-supported phylogenies, we recovered strong evidence of a reticulated evolutionary history using a network approach that incorporates both incomplete lineage sorting and hybridization. These data suggest that the rock-dwelling western North American subalpine endemic R. shushanii is potentially the result of a hybrid speciation event, and introgression may have also played a role in other taxa, including vagrant species R. arbuscula, R. haydenii and R. idahoensis. We discuss the potential roles of hybridization in terms of generating asexuality and novel morphological traits in lichens. Furthermore, our results highlight the need for additional study of reticulate phylogenies when investigating species boundaries and evolutionary history, even in cases with well-supported topologies inferred from genome-scale data.

Integrating hybrid zone analyses in species delimitation: lessons from two anuran radiations of the Western Mediterranean

Molecular ecologists often rely on phylogenetic evidence for assessing the species-level systematics of newly discovered lineages. Alternatively, the extent of introgression at phylogeographic transitions can provide a more direct test to assign candidate taxa into subspecies or species categories. Here, we compared phylogenetic versus hybrid zone approaches of species delimitation in two groups of frogs from the Western Mediterranean region (Discoglossus and Pelodytes), by using genomic data (ddRAD). In both genera, coalescent analyses recovered almost all nominal taxa as "species". However, the least-diverged pairs D. g. galganoi/jeanneae and P. punctatus/hespericus admix over hundreds of kilometers, suggesting that they have not yet developed strong reproductive isolation and should be treated as conspecifics. In contrast, the comparatively older D. scovazzi/pictus and P. atlanticus/ibericus form narrow contact zones, consistent with species distinctiveness. Due to their complementarity, we recommend taxonomists to combine phylogenomics with hybrid zone analyses to scale the gray zone of speciation, i.e., the evolutionary window separating widely admixing lineages versus nascent reproductively isolated species. The radically different transitions documented here conform to the view that genetic incompatibilities accumulating with divergence generate a weak barrier to gene flow for long periods of time, until their effects multiply and the speciation process then advances rapidly. Given the variability of the gray zone among taxonomic groups, at least from our current abilities to measure it, we recommend to customize divergence thresholds within radiations to categorize lineages for which no direct test of speciation is possible.

The roles of vicariance and dispersal in the differentiation of two species of the Rhinella marina species complex

The high levels of Neotropical biodiversity are commonly associated with the intense Neogene-Quaternary geological events and climate dynamics. Here, we investigate the evolutionary history of two species of Neotropical closely related amphibians (R. horribilis and R. marina). We combine published data with new mitochondrial DNA sequences and multiple nuclear markers, including 12 microsatellites. The phylogenetic analyses showed support for grouping the samples in two main clades; R. horribilis (Central America and Mexico) and R. marina (South America east of the Andes). However, the phylogenetic inferences also show an evident mito-nuclear discordance. We use Approximate Bayesian Computation (ABC) to test the role of different events in the diversification between the two groups recovered. We found that both species were affected primarily by a recent Pleistocene divergence, which was similar to the divergence estimate revealed by the Isolation-with-Migration model, under persistent bidirectional gene flow through time. We provide the first evidence that R. horribilis is differentiated from the South American R. marina at the nuclear level supporting the taxonomic status of R. horribilis, which has been controversial for more than a century.

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.

Phylotranscriptomic Insights into the Diversification of Endothermic Thunnus Tunas

Birds, mammals, and certain fishes, including tunas, opahs and lamnid sharks, are endothermic, conserving internally generated, metabolic heat to maintain body or tissue temperatures above that of the environment. Bluefin tunas are commercially important fishes worldwide, and some populations are threatened. They are renowned for their endothermy, maintaining elevated temperatures of the oxidative locomotor muscle, viscera, brain and eyes, and occupying cold, productive high-latitude waters. Less cold-tolerant tunas, such as yellowfin tuna, by contrast, remain in warm-temperate to tropical waters year-round, reproducing more rapidly than most temperate bluefin tuna populations, providing resiliency in the face of large-scale industrial fisheries. Despite the importance of these traits to not only fisheries but also habitat utilization and responses to climate change, little is known of the genetic processes underlying the diversification of tunas. In collecting and analyzing sequence data across 29,556 genes, we found that parallel selection on standing genetic variation is associated with the evolution of endothermy in bluefin tunas. This includes two shared substitutions in genes encoding glycerol-3 phosphate dehydrogenase, an enzyme that contributes to thermogenesis in bumblebees and mammals, as well as four genes involved in the Krebs cycle, oxidative phosphorylation, β-oxidation, and superoxide removal. Using phylogenetic techniques, we further illustrate that the eight Thunnus species are genetically distinct, but found evidence of mitochondrial genome introgression across two species. Phylogeny-based metrics highlight conservation needs for some of these species.

Reconciling the Mitonuclear Compatibility Species Concept with Rampant Mitochondrial Introgression

The mitonuclear compatibility species concept defines a species as a population that is genetically isolated from other populations by uniquely coadapted mitochondrial (mt) and nuclear genes. A key prediction of this hypothesis is that the mt genotype of each species will be functionally distinct and that introgression of mt genomes will be prevented by mitonuclear incompatibilities that arise when heterospecific mt and nuclear genes attempt to cofunction to enable aerobic respiration. It has been proposed, therefore, that the observation of rampant introgression of mt genotypes from one species to another constitutes a strong refutation of the mitonuclear speciation. The displacement of a mt genotype from a nuclear background with which it co-evolved to a foreign nuclear background will necessarily lead to fitness loss due to mitonuclear incompatibilities. Here I consider two potential benefits of mt introgression between species that may, in some cases, overcome fitness losses arising from mitonuclear incompatibilities. First, the introgressed mt genotype may be better adapted to the local environment than the native mt genotype such that higher fitness is achieved through improved adaptation via introgression. Second, if the mitochondria of the recipient taxa carry a high mutational load, then introgression of a foreign, less corrupt mt genome may enable the recipient taxa to escape its mutational load and gain a fitness advantage. Under both scenarios, fitness gains from novel mt genotypes could theoretically compensate for the fitness that is lost via mitonuclear incompatibility. I also consider the role of endosymbionts in non-adaptive rampant introgression of mt genomes. I conclude that rampant introgression is not necessarily evidence against the idea of tight mitonuclear coadaptation or the mitonuclear compatibility species concept. Rampant mt introgression will typically lead to erasure of species but in some cases could lead to hybrid speciation.

Phylogeographic diversification and postglacial range dynamics shed light on the conservation of the kelp Saccharina japonica

Studies of postglacial range shifts could enhance our understanding of seaweed species' responses to climate change and hence facilitate the conservation of natural resources. However, the distribution dynamics and phylogeographic diversification of the commercially and ecologically important kelp Saccharina japonica in the Northwest Pacific (NWP) are still poorly surveyed. In this study, we analyzed the evolutionary history of S. japonica using two mitochondrial markers and 24 nuclear microsatellites. A STRUCTURE analysis revealed two partially isolated lineages: lineage H, which is scattered along the coast of Japan; and lineage P, which occurs along the west coast of the Japan Sea. Ecological niche modeling projections to the Last Glacial Maximum (LGM) revealed that the southern coasts of the Japan Sea and the Pacific side of the Oshima and Honshu Peninsulas provided the most suitable habitats for S. japonica, implying that these regions served as ancient refugia during the LGM. Ancient isolation in different refugia may explain the observed divergence between lineages P and H. An approximate Bayesian computation analysis indicated that the two lineages experienced post-LGM range expansion and that postglacial secondary contact occurred in Sakhalin. Model projections into the year 2,100 predicted that S. japonica will shift northwards and lose its genetic diversity center on the Oshima Peninsula in Hokkaido and Shimokita Peninsula in Honshu. The range shifts and evolutionary history of S. japonica improve our understanding of how climate change impacted the distribution range and diversity of this species and provide useful information for the conservation of natural resources under ongoing environmental change in the NWP.

Identifying hybrids & the genomics of hybridization: Mallards & American black ducks of Eastern North America

Resolving evolutionary relationships and establishing population structure depends on molecular diagnosability that is often limited for closely related taxa. Here, we use 3,200 ddRAD-seq loci across 290 mallards, American black ducks, and putative hybrids to establish population structure and estimate hybridization rates. We test between traditional assignment probability and accumulated recombination events based analyses to assign hybrids to generational classes. For hybrid identification, we report the distribution of recombination events complements ADMIXTURE simulation by extending resolution past F4 hybrid status; however, caution against hybrid assignment based on accumulated recombination events due to an inability to resolve F1 hybrids. Nevertheless, both analyses suggest that there are relatively few backcrossed stages before a lineage's hybrid ancestry is lost and the offspring are effectively parental again. We conclude that despite high rates of observed interspecific hybridization between mallards and black ducks in the middle part of the 20th century, our results do not support the predicted hybrid swarm. Conversely, we report that mallard samples genetically assigned to western and non-western clusters. We indicate that these non-western mallards likely originated from game-farm stock, suggesting landscape level gene flow between domestic and wild conspecifics.

Diversification and speciation in tree frogs from the Maghreb (Hyla meridionalis sensu lato), with description of a new African endemic

Comparative molecular studies emphasized a new biogeographic paradigm for the terrestrial fauna of North Africa, one of the last uncharted ecoregions of the Western Palearctic: two independent east-west divisions across the Maghreb. Through a comprehensive phylogeography, we assessed how this model suits the genetic diversification documented for the tree frog Hyla meridionalis sensu lato. Analyses of mtDNA variation and thousands of nuclear loci confirmed the old split (low-Pliocene) between Tunisian and Algerian populations. These lineages meet but barely admix in the eastern Maghreb (Algerian-Tunisian border), a sign of putatively advanced reproductive isolation. In the western Maghreb, we report a Pleistocene divergence between Moroccan and Algerian populations. Tree frogs thus follow both predictions: a double east-west break that gave rise to two suture zones characteristic of North-African phylogeography. Moreover, some intraspecific mtDNA variation is not mirrored by the nuclear data, emphasizing that evolutionary units should always be designated by multilocus approaches. Last but not least, we describe the Tunisian lineage as a new species endemic to Africa.

Population Connectivity and Traces of Mitochondrial Introgression in New Zealand Black-Billed Gulls (Larus bulleri)

Black-billed gulls (Larus bulleri) are endemic to New Zealand and are suspected to be undergoing substantial population declines. They primarily breed on open gravel beds in braided rivers of the South Island—a habitat that is diminishing and becoming increasingly modified. Although management of this species is increasing, little has been published on their movements and demographics. In this study, both mitochondrial DNA (mtDNA) control region domain I and nuclear single nucleotide polymorphisms (SNPs) were examined to help understand the connectivity and population structure of black-billed gulls across the country and to help inform management decisions. Mitochondrial DNA showed no population structure, with high haplotype and low nucleotide diversity, and analyses highlighted mitochondrial introgression with the closely related red-billed gulls (Larus novaehollandiae scopulinus). Nuclear DNA analyses, however, identified two groups, with Rotorua birds in the North Island being distinct from the rest of New Zealand, and isolation-by-distance evident across the South Island populations. Gene flow primarily occurs between nearby colonies with a stepwise movement across the landscape. The importance from a genetic perspective of the more isolated North Island birds (1.6% of total population) needs to be further evaluated. From our results, we infer that the South Island black-billed gull management should focus on maintaining several populations within each region rather than focusing on single specific colonies or river catchments. Future study is needed to investigate the genetic structure of populations at the northern limit of the species’ range, and identify the mechanisms behind, and extent of, the hybridisation between red-billed and black-billed gulls.

Genome-wide markers reveal a complex evolutionary history involving divergence and introgression in the Abert's squirrel (Sciurus aberti) species group

Background

Genetic introgression between divergent lineages is now considered more common than previously appreciated, with potentially important consequences for adaptation and speciation. Introgression is often asymmetric between populations and patterns can vary for different types of loci (nuclear vs. organellar), complicating phylogeographic reconstruction. The taxonomy of the ecologically specialized Abert’s squirrel species group has been controversial, and previous studies based on mitochondrial data have not fully resolved the evolutionary relationships among populations. Moreover, while these studies identified potential areas of secondary contact between divergent lineages, the possibility for introgression has not been tested.

Results

We used RAD-seq to unravel the complex evolutionary history of the Abert’s squirrel species group. Although some of our findings reinforce inferences based on mitochondrial data, we also find significant areas of discordance. Discordant signals generally arise from previously undetected introgression between divergent populations that differentially affected variation at mitochondrial and nuclear loci. Most notably, our results support earlier claims (disputed by mitochondrial data) that S. aberti kaibabensis, found only on the north rim of the Grand Canyon, is highly divergent from other populations. However, we also detected introgression of S. aberti kaibabensis DNA into other S. aberti populations, which likely accounts for the previously inferred close genetic relationship between this population and those south of the Grand Canyon.

Conclusions

Overall, the evolutionary history of Abert’s squirrels appears to be shaped largely by divergence during periods of habitat isolation. However, we also found evidence for interbreeding during periods of secondary contact resulting in introgression, with variable effects on mitochondrial and nuclear markers. Our results support the emerging view that populations often diversify under scenarios involving both divergence in isolation and gene flow during secondary contact, and highlight the value of genome-wide datasets for resolving such complex evolutionary histories.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1248-4) contains supplementary material, which is available to authorized users.

Hybridization and extensive mitochondrial introgression among fire salamanders in peninsular Italy

Discordance between mitochondrial and nuclear patterns of population genetic structure is providing key insights into the eco-evolutionary dynamics between and within species, and their assessment is highly relevant to biodiversity monitoring practices based on DNA barcoding approaches. Here, we investigate the population genetic structure of the fire salamander Salamandra salamandra in peninsular Italy. Both mitochondrial and nuclear markers clearly identified two main population groups. However, nuclear and mitochondrial zones of geographic transition between groups were located 600 km from one another. Recent population declines in central Italy partially erased the genetic imprints of past hybridization dynamics. However, the overall pattern of genetic variation, together with morphological and fossil data, suggest that a rampant mitochondrial introgression triggered the observed mitonuclear discordance, following a post-glacial secondary contact between lineages. Our results clearly show the major role played by reticulate evolution in shaping the structure of Salamandra salamandra populations and, together with similar findings in other regions of the species' range, contribute to identify the fire salamander as a particularly intriguing case to investigate the complexity of mechanisms triggering patterns of mitonuclear discordance in animals.

The genomic impact of historical hybridization with massive mitochondrial DNA introgression

BACKGROUND

The extent to which selection determines interspecific patterns of genetic exchange enlightens the role of adaptation in evolution and speciation. Often reported extensive interspecific introgression could be selection-driven, but also result from demographic processes, especially in cases of invasive species replacements, which can promote introgression at their invasion front. Because invasion and selective sweeps similarly mold variation, population genetics evidence for selection can only be gathered in an explicit demographic framework. The Iberian hare, Lepus granatensis, displays in its northern range extensive mitochondrial DNA introgression from L. timidus, an arctic/boreal species that it replaced locally after the last glacial maximum. We use whole-genome sequencing to infer geographic and genomic patterns of nuclear introgression and fit a neutral model of species replacement with hybridization, allowing us to evaluate how selection influenced introgression genome-wide, including for mtDNA.

RESULTS

Although the average nuclear and mtDNA introgression patterns contrast strongly, they fit a single demographic model of post-glacial invasive replacement of timidus by granatensis. Outliers of elevated introgression include several genes related to immunity, spermatogenesis, and mitochondrial metabolism. Introgression is reduced on the X chromosome and in low recombining regions.

CONCLUSIONS

General nuclear and mtDNA patterns of introgression can be explained by purely demographic processes. Hybrid incompatibilities and interplay between selection and recombination locally modulate levels of nuclear introgression. Selection promoted introgression of some genes involved in conflicts, either interspecific (parasites) or possibly cytonuclear. In the latter case, nuclear introgression could mitigate the potential negative effects of alien mtDNA on mitochondrial metabolism and male-specific traits.

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build-up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double-digest restriction site-associated DNA sequencing (ddRADseq) to obtain and analyse genomic-scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.

Multi-locus phylogenetics, lineage sorting, and reticulation in Pinus subsection Australes

PREMISE OF THE STUDY

Both incomplete lineage sorting and reticulation have been proposed as causes of phylogenetic incongruence. Disentangling these factors may be most difficult in long-lived, wind-pollinated plants with large population sizes and weak reproductive barriers.

METHODS

We used solution hybridization for targeted enrichment and massive parallel sequencing to characterize low-copy-number nuclear genes and high-copy-number plastomes (Hyb-Seq) in 74 individuals of Pinus subsection Australes, a group of ~30 New World pine species of exceptional ecological and economic importance. We inferred relationships using methods that account for both incomplete lineage sorting and reticulation.

KEY RESULTS

Concatenation- and coalescent-based trees inferred from nuclear genes mainly agreed with one another, but they contradicted the plastid DNA tree in recovering the Attenuatae (the California closed-cone pines) and Oocarpae (the egg-cone pines of Mexico and Central America) as monophyletic and the Australes sensu stricto (the southern yellow pines) as paraphyletic to the Oocarpae. The plastid tree featured some relationships that were discordant with morphological and geographic evidence and species limits. Incorporating gene flow into the coalescent analyses better fit the data, but evidence supporting the hypothesis that hybridization explains the non-monophyly of the Attenuatae in the plastid tree was equivocal.

CONCLUSIONS

Our analyses document cytonuclear discordance in Pinus subsection Australes. We attribute this discordance to ancient and recent introgression and present a phylogenetic hypothesis in which mostly hierarchical relationships are overlain by gene flow.