The biogeography of mitochondrial and nuclear discordance in animals

Are members of the Anopheles fluviatilis complex conspecific?

Anopheles fluviatilis sensu lato, a primary malaria vector in India, has been identified to be comprised of four cryptic species, provisionally designated as species S, T, U and V. However, Kumar et al. (Mol Ecol Resour, 2013;13:354-61) considered all of the then known three members of this species complex (S, T and U) conspecific. The specific status of species S and T was refuted based on the lack of sufficient barcode gap in mitochondrial-CO1 and the perceived presence of heterozygotes in populations as detected through one of the two species-specific PCR assays employed for the cryptic species identification. The existence of species U was refuted claiming that earlier investigations have already refuted their existence. Here we discuss problems associated with the CO1-based barcode approach for delimitation of cryptic species, the perceived heterozygosity between species S and T based on a species-specific PCR assay, and interpretation of published reports. We demonstrated that fixed differences do exist in the ITS2-rDNA sequence of species S and T with no evidence of heterozygotes in sympatric populations and, that the observed heterozygosity by Kumar et al. in the ITS2-based species diagnostic PCR is due to the high mispriming tendency of the T-specific primer with species S. We infer that mitochondrial DNA-based 'barcoding gap', an arbitrary threshold recommended for species delimitation, alone, is inadequate to delimit the members of An. fluviatilis complex.

Different rates of pollen and seed gene flow cause branch-length and geographic cytonuclear discordance within Asian butternuts

Topological cytonuclear discordance is commonly observed in plant phylogenetic and phylogeographic studies, yet few studies have attempted to detect two other forms of cytonuclear discordance (branch length and geographical) and to uncover the causes of the discordance. We used the whole nuclear and chloroplast genome data from 80 individual Asian butternuts to reveal the pattern and processes of cytonuclear discordance. Our findings indicate that the chloroplast genome had substantially deeper divergence (branch-length discordance) and a steeper cline in the contact zone (geographic discordance) compared with the nuclear genome. After various hypothesis have been tested, the results suggest that incomplete lineage sorting, positive selection and cytonuclear incompatibility are probably insufficient to explain this pattern. However, isolation-by-distance analysis and gene flow estimation point to a much higher level of gene flow by pollen compared with by seeds, which may have slowed down lineage divergence and mediated wider contact for nuclear genome compared with the chloroplast genome. Altogether, this study highlights a critical role of sex-biased dispersal in causing discordance between the nuclear and plastid genome of Asian butternuts. Given its ubiquity among plants, asymmetric gene flow should be given a high priority in future studies of cytonuclear discordance.

Northward range extension for Durvillaea poha bull kelp: Response to tectonic disturbance?

Understanding the forces that shape species distributions is increasingly important in a fast-changing world. Although major disturbance events can adversely affect natural populations, they can also present new opportunities, for example by opening up habitat for colonization by other lineages. Following extensive geographic sampling, we use genomic data to infer a range extension following disturbance for an ecologically important intertidal macroalgal species. Specifically, we genotyped 288 southern bull kelp (Durvillaea) plants from 28 localities across central New Zealand. All specimens from the North Island were expected to be D. antarctica, but unexpectedly 10 samples from four sites were identified as D. poha. Extensive sampling from the northern South Island (105 samples at five locations) confirmed the absence of D. poha north of the Kaikōura Peninsula. The North Island specimens of D. poha therefore reveal a biogeographic disjunction, some 150 km northeast of the nearest (South Island) population of this species. Based on strong geographic correspondence between these North Island samples and historic disturbance, we infer that tectonic upheaval, particularly earthquake-generated landslides, likely extirpated local D. antarctica and created an opportunity for a northward range expansion event by D. poha. Close phylogenomic relationships between this new North Island population and South Island samples support a geologically recent northward expansion, rather than a deeper evolutionary origin. These findings indicate the potential of large-scale disturbances to facilitate sudden biogeographic range expansions, and they emphasize the ability of genomic analyses with fine-scale sampling to reveal long-lasting signatures of past disturbance, dispersal, and colonization.

Genomic Evidence for Rare Hybridization and Large Demographic Changes in the Evolutionary Histories of Four North American Dove Species

Introductions and invasions provide opportunities for interaction and hybridization between colonists and closely related native species. We investigate this phenomenon using the mitochondrial DNA COI and 81,416 base-pairs of overlapping nuclear variation to examine the evolutionary histories and signatures of hybridization among introduced feral Rock Pigeon and Eurasian Collared-Dove and native White-winged and Mourning doves in southwestern North America. First, we report all four species to be highly divergent across loci (overall pair-wise species Φ_ST range = 0.17-0.70) and provide little evidence for gene flow at evolutionary timescales. Despite this, evidence from multiple population genetics analyses supports the presence of six putative contemporary late-stage hybrids among the 182 sampled individuals. These putative hybrids contain various ancestry combinations, but all involve the most populous species, the Mourning Dove. Next, we use a novel method to reconstruct demographic changes through time using partial genome sequence data. We identify recent, species-specific fluctuations in population size that are likely associated with changing environments since the Miocene and suggest that these fluctuations have influenced the genetic diversity of each dove species in ways that may impact their future persistence. Finally, we discuss the importance of using multiple marker types when attempting to infer complex evolutionary histories and propose important considerations when analyzing populations that were recently established or of domestic origins.

Molecular phylogeny of the Anopheles hyrcanus group (Diptera: Culicidae) based on rDNA-ITS2 and mtDNA-COII

BACKGROUND

The Anopheles hyrcanus group, which includes 25 species, is widely distributed in the Oriental and Palaearctic regions. Given the difficulty in identifying cryptic or sibling species based on their morphological characteristics, molecular identification is regarded as an important complementary approach to traditional morphological taxonomy. The aim of this study was to reconstruct the phylogeny of the Hyrcanus group using DNA barcoding markers in order to determine the phylogenetic correlations of closely related taxa and to compare these markers in terms of identification efficiency and genetic divergence among species.

METHODS

Based on data extracted from the GenBank database and data from the present study, we used 399 rDNA-ITS2 sequences of 19 species and 392 mtDNA-COII sequences of 14 species to reconstruct the molecular phylogeny of the Hyrcanus group across its worldwide range. We also compared the performance of rDNA-ITS2 against that of mtDNA-COII to assess the genetic divergence of closely related species within the Hyrcanus group.

RESULTS

Average interspecific divergence for the rDNA-ITS2 sequence (0.376) was 125-fold higher than the average intraspecies divergence (0.003), and average interspecific divergence for the mtDNA-COII sequence (0.055) was eightfold higher than the average intraspecies divergence (0.007). The barcoding gap ranged from 0.015 to 0.073 for rDNA-ITS2, and from 0.017 to 0.025 for mtDNA-COII. Two sets of closely related species, namely, Anophels lesteri and An. paraliae, and An. sinensis, An. belenrae and An. kleini, were resolved by rDNA-ITS2. In contrast, the relationship of An. sinensis/An. belenrae/An. kleini was poorly defined in the COII tree. The neutrality test and mismatch distribution revealed that An. peditaeniatus, An. hyrcanus, An. sinensis and An. lesteri were likely to undergo hitchhiking or population expansion in accordance with both markers. In addition, the population of an important vivax malaria vector, An. sinensis, has experienced an expansion after a bottleneck in northern and southern Laos.

CONCLUSIONS

The topology of the Hyrcanus group rDNA-ITS2 and mtDNA-COII trees conformed to the morphology-based taxonomy for species classification rather than for that for subgroup division. rDNA-ITS2 is considered to be a more reliable diagnostic tool than mtDNA-COII in terms of investigating the phylogenetic correlation between closely related mosquito species in the Hyrcanus group. Moreover, the population expansion of an important vivax malaria vector, An. sinensis, has underlined a potential risk of malaria transmission in northern and southern Laos. This study contributes to the molecular identification of the Anopheles hyrcanus group in vector surveillance.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

Mitochondrial function involves the interplay between mitochondrial and nuclear genomes. Such mitonuclear interactions can be disrupted by the introgression of mitochondrial DNA between taxa or divergent populations. Previous studies of several model systems (e.g. Drosophila) indicate that the disruption of mitonuclear interactions, termed mitonuclear mismatch, can alter nuclear gene expression, yet few studies have focused on natural populations.

Results

Here we study a naturally introgressed population in the secondary contact zone of two subspecies of the intermediate horseshoe bat (Rhinolophus affinis), in which individuals possess either mitonuclear matched or mismatched genotypes. We generated transcriptome data for six tissue types from five mitonuclear matched and five mismatched individuals. Our results revealed strong tissue-specific effects of mitonuclear mismatch on nuclear gene expression with the largest effect seen in pectoral muscle. Moreover, consistent with the hypothesis that genes associated with the response to oxidative stress may be upregulated in mitonuclear mismatched individuals, we identified several such gene candidates, including DNASE1L3, GPx3 and HSPB6 in muscle, and ISG15 and IFI6 in heart.

Conclusion

Our study reveals how mitonuclear mismatch arising from introgression in natural populations is likely to have fitness consequences. Underlying the processes that maintain mitonuclear discordance is a step forward to understand the role of mitonuclear interactions in population divergence and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00424-x.

Delimitation despite discordance: Evaluating the species limits of a confounding species complex in the face of mitonuclear discordance

The delimitation of species is an essential pursuit of biology, and proper taxonomies are crucial for the assessment and conservation management of organismal diversity. However, delimiting species can be hindered by a number of factors including highly conserved morphologies (e.g., cryptic species), differences in criteria of species concepts, lineages being in the early stages of the speciation or divergence process, and discordance between gene topologies (e.g., mitonuclear discordance). Here we use a taxonomically confounded species complex of toads in Central America that exhibits extensive mitonuclear discordance to test delimitation hypotheses. Our investigation integrates mitochondrial sequences, nuclear SNPs, morphology, and macroecological data to determine which taxonomy best explains the divergence and evolutionary relationships among these toads. We found that a three species taxonomy following the distributions of the nuclear SNP haplotypes offers the best explanation of the species in this complex based off of the integrated data types. Due to the taxonomic instability of this group, we also discuss conservation concerns in the face of improper taxonomic delimitation. Our study provides an empirical and integrative hypothesis testing framework to assess species delimitation hypotheses in the face of cryptic morphology and mitonuclear discordance and highlights the importance that a stable taxonomy has over conservation-related actions.

Opportunities and challenges of macrogenetic studies

The rapidly emerging field of macrogenetics focuses on analysing publicly accessible genetic datasets from thousands of species to explore large-scale patterns and predictors of intraspecific genetic variation. Facilitated by advances in evolutionary biology, technology, data infrastructure, statistics and open science, macrogenetics addresses core evolutionary hypotheses (such as disentangling environmental and life-history effects on genetic variation) with a global focus. Yet, there are important, often overlooked, limitations to this approach and best practices need to be considered and adopted if macrogenetics is to continue its exciting trajectory and reach its full potential in fields such as biodiversity monitoring and conservation. Here, we review the history of this rapidly growing field, highlight knowledge gaps and future directions, and provide guidelines for further research.

Evidence for multiple colonisations and Wolbachia infections shaping the genetic structure of the widespread butterfly Polyommatus icarus in the British Isles

The paradigm of isolation in southern refugia during glacial periods followed by expansions during interglacials, producing limited genetic differentiation in northern areas, dominates European phylogeography. However, the existence of complex structured populations in formerly glaciated areas, and islands connected to mainland areas during glacial maxima, call for alternative explanations. We reconstructed the mtDNA phylogeography of the widespread Polyommatus Icarus butterfly with an emphasis on the formerly glaciated and connected British Isles. We found distinct geographical structuring of CO1 haplogroups, with an ancient lineage restricted to the marginal European areas, including Northern Scotland and Outer Hebrides. Population genomic analyses, using ddRADSeq genomic markers, also reveal substantial genetic structuring within Britain. However, there is negligble mito-nuclear concordance consistent with independent demographic histories of mitochondrial versus nuclear DNA. While mtDNA-Wolbachia associations in northern Britain could account for the geographic structuring of mtDNA across most of the British Isles, for nuclear DNA markers (derived from ddRADseq data) butterflies from France cluster between northern and southern British populations - an observation consistent with a scenario of multiple recolonisation. Taken together our results suggest that contemporary mtDNA structuring in the British Isles (and potentially elsewhere in Europe) largely results from Wolbachia infections, however, nuclear genomic structuring suggests a history of at least two distinct colonisations. This two-stage colonisation scenario has previously been put forth to explain genetic diversity and structuring in other British flora and fauna. Additionally, we also present preliminary evidence for potential Wolbachia-induced feminization in the Outer Hebrides.

Ancient divergence of Indian and Tibetan wolves revealed by recombination-aware phylogenomics

The grey wolf (Canis lupus) expanded its range across Holarctic regions during the late Pleistocene. Consequently, most grey wolves share recent (<100,000 years ago) maternal origins corresponding to a widespread Holarctic clade. However, two deeply divergent (200,000-700,000 years ago) mitochondrial clades are restricted, respectively, to the Indian subcontinent and the Tibetan Plateau, where remaining wolves are endangered. No genome-wide analysis had previously included wolves corresponding to the mitochondrial Indian clade or attempted to parse gene flow and phylogeny. We sequenced four Indian and two Tibetan wolves and included 31 additional canid genomes to resolve the phylogenomic history of grey wolves. Genomic analyses revealed Indian and Tibetan wolves to be distinct from each other and from broadly distributed wolf populations corresponding to the mitochondrial Holarctic clade. Despite gene flow, which was reflected disproportionately in high-recombination regions of the genome, analyses revealed Indian and Tibetan wolves to be basal to Holarctic grey wolves, in agreement with the mitochondrial phylogeny. In contrast to mitochondrial DNA, however, genomic findings suggest the possibility that the Indian wolf could be basal to the Tibetan wolf, a discordance potentially reflecting selection on the mitochondrial genome. Together, these findings imply that southern regions of Asia have been important centers for grey wolf evolution and that Indian and Tibetan wolves represent evolutionary significant units (ESUs). Further study is needed to assess whether these ESUs warrant recognition as distinct species. This question is especially urgent regarding the Indian wolf, which represents one of the world's most endangered wolf populations.

Concordance of the spectral properties of dorsal wing scales with the phylogeographic structure of European male Polyommatus icarus butterflies

The males of more than 80% of the Lycaenidae species belonging to the tribe Polyommatini exhibit structural coloration on their dorsal wing surfaces. These colors have a role in reinforcement in prezygotic reproductive isolation. The species-specific colors are produced by the cellular self-assembly of chitin/air nanocomposites. The spectral position of the reflectance maximum of such photonic nanoarchitectures depends on the nanoscale geometric dimensions of the elements building up the nanostructure. Previous work showed that the coloration of male Polyommatus icarus butterflies in the Western and Eastern Palearctic exhibits a characteristic spectral difference (20 nm). We investigated the coloration and the de novo developed DNA microsatellites of 80 P. icarus specimens from Europe from four sampling locations, spanning a distance of 1621 km. Remarkably good concordance was found between the spectral properties of the blue sexual signaling color (coincident within 5 nm) and the population genetic structure as revealed by 10 microsatellites for the P. icarus species.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Phylogenomic Assessment of Biodiversity Using a Reference-Based Taxonomy: An Example With Horned Lizards (Phrynosoma)

Phylogenomic investigations of biodiversity facilitate the detection of fine-scale population genetic structure and the demographic histories of species and populations. However, determining whether or not the genetic divergence measured among populations reflects species-level differentiation remains a central challenge in species delimitation. One potential solution is to compare genetic divergence between putative new species with other closely related species, sometimes referred to as a reference-based taxonomy. To be described as a new species, a population should be at least as divergent as other species. Here, we develop a reference-based taxonomy for Horned Lizards (Phrynosoma; 17 species) using phylogenomic data (ddRADseq data) to provide a framework for delimiting species in the Greater Short-horned Lizard species complex (P. hernandesi). Previous species delimitation studies of this species complex have produced conflicting results, with morphological data suggesting that P. hernandesi consists of five species, whereas mitochondrial DNA support anywhere from 1 to 10 + species. To help address this conflict, we first estimated a time-calibrated species tree for P. hernandesi and close relatives using SNP data. These results support the paraphyly of P. hernandesi; we recommend the recognition of two species to promote a taxonomy that is consistent with species monophyly. There is strong evidence for three populations within P. hernandesi, and demographic modeling and admixture analyses suggest that these populations are not reproductively isolated, which is consistent with previous morphological analyses that suggest hybridization could be common. Finally, we characterize the population-species boundary by quantifying levels of genetic divergence for all 18 Phrynosoma species. Genetic divergence measures for western and southern populations of P. hernandesi failed to exceed those of other Phrynosoma species, but the relatively small population size estimated for the northern population causes it to appear as a relatively divergent species. These comparisons underscore the difficulties associated with putting a reference-based approach to species delimitation into practice. Nevertheless, the reference-based approach offers a promising framework for the consistent assessment of biodiversity within clades of organisms with similar life histories and ecological traits.

Phylogenomic systematics of the spotted skunks (Carnivora, Mephitidae, Spilogale): Additional species diversity and Pleistocene climate change as a major driver of diversification

Four species of spotted skunks (Carnivora, Mephitidae, Spilogale) are currently recognized: Spilogale angustifrons, S. gracilis, S. putorius, and S. pygmaea. Understanding species boundaries within this group is critical for effective conservation given that regional populations or subspecies (e.g., S. p. interrupta) have experienced significant population declines. Further, there may be currently unrecognized diversity within this genus as some taxa (e.g., S. angustifrons) and geographic regions (e.g., Central America) never have been assessed using DNA sequence data. We analyzed species limits and diversification patterns in spotted skunks using multilocus nuclear (ultraconserved elements) and mitochondrial (whole mitogenomes and single gene analysis) data sets from broad geographic sampling representing all currently recognized species and subspecies. We found a high degree of genetic divergence among Spilogale that reflects seven distinct species and eight unique mitochondrial lineages. Initial divergence between S. pygmaea and all other Spilogale occurred in the Early Pliocene (∼ 5.0 million years ago). Subsequent diversification of the remaining Spilogale into an "eastern" and a "western" lineage occurred during the Early Pleistocene (∼1.5 million years ago). These two lineages experienced temporally coincident patterns of diversification at ∼0.66 and ∼0.35 million years ago into two and ultimately three distinct evolutionary units, respectively. Diversification was confined almost entirely within the Pleistocene during a timeframe characterized by alternating glacial-interglacial cycles, with the origin of this diversity occurring in northeastern Mexico and the southwestern United States of America. Mitochondrial-nuclear discordance was recovered across three lineages in geographic regions consistent with secondary contact, including a distinct mitochondrial lineage confined to the Sonoran Desert. Our results have direct consequences for conservation of threatened populations, or species, as well as for our understanding of the evolution of delayed implantation in this enigmatic group of small carnivores.

Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.

Massive gene rearrangement in mitogenomes of phytoseiid mites

Mitochondrial (mt) gene sequences have been widely used to infer phylogeny in animals. The relative order of mt genes in the mitogenome can also be a useful marker for evolution, but the propensity of mt gene rearrangements vary tremendously among taxa. Ticks and mites in Acari exemplify this trend as some families retain the ancestral arthropod gene order, while others show highly divergent gene orders. Mites in Phytoseiidae, many of which are effective biological control agents, show some of the most divergent gene orders. However, the diversity of mitogenome order within this family is little known. We thus sequenced three mt genomes of phytoseiid mites from two of the most speciose genera: Amblyseius swirskii (Athias-Henriot), Amblyseius tsugawai (Ehara) and Neoseiulus womersleyi (Schicha). We find differences in mt GC skew and nucleotide composition, especially between N. womersleyi and the two Amblyseius species. Each species within Phytoseiidae (including three previously available sequences) present a unique gene order. Phytoseiid mitogenomes show some of the highest numbers of breakpoints when compared to the ancestral arthropod order (up to 33), as well as high numbers of breakpoints within the family (14-30). This suggests a history of massive, ongoing mitogenome rearrangements in the family. Phylogenetic analyses of mt sequences confirm that the degree of gene rearrangements follows phylogenetic relatedness. We discuss possible causes for the high degree of mt gene rearrangement within phytoseiid mites as well as selection in the mt and nuclear genome tied to the independent evolution of many diverse feeding strategies in the family. Finally, we suggest N. womersleyi should be used instead of the synonym Amblyseius pseudolongispinosus.

Sympatric occurrence of deeply diverged mitochondrial DNA lineages in Siberian geometrid moths (Lepidoptera: Geometridae): cryptic speciation, mitochondrial introgression, secondary admixture or effect of Wolbachia?

The divergent sympatric mitochondrial lineages within traditionally recognized species present a challenge regularly faced by taxonomists and evolutionary biologists. We encountered this problem when studying the Siberian geometrid moths, Alcis deversata and Thalera chlorosaria. Within each of these species we found two deeply diverged mitochondrial lineages that demonstrated a level of genetic differentiation exceeding the standard interspecific DNA barcode threshold. Using analyses of nuclear genes, morphology, ecological preferences and Wolbachia endosymbionts, we tested five hypotheses that might explain the mitochondrial pattern observed: cryptic speciation, ancestral polymorphism, interspecific mitochondrial introgression, secondary admixture of allopatrically evolved populations and an effect of intracellular Wolbachia endosymbionts. We demonstrate that in A. deversata and Th. chlorosaria the mitochondrial differences are not correlated with differences in nuclear genes, morphology, ecology and Wolbachia infection status, thus not supporting the hypothesis of cryptic species and an effect of Wolbachia. Mitochondrial introgression can lead to a situation in which one species has both its own mitochondrial lineage and the lineage obtained from another species. We found this situation in the species pair Alcis repandata and Alcis extinctaria. We conclude that the mitochondrial heterogeneity in A. deversata and Th. chlorosaria is most likely to be attributable to the secondary admixture of allopatrically evolved populations.

Whole-genome sequencing reveals asymmetric introgression between two sister species of cold-resistant leaf beetles

A large number of genetic variation studies have identified cases of mitochondrial genome introgression in animals, indicating that reproductive barriers among closely related species are often permeable. Because of its sheer size, the impact of hybridization on the evolution of the nuclear genome is more difficult to apprehend. Only a few studies have explored it recently thanks to recent progress in DNA sequencing and genome assembly. Here, we analysed whole-genome sequence variation among multiple individuals of two sister species of leaf beetles inside their hybrid zone, in which asymmetric mitochondrial genome introgression had previously been established. We used a machine learning approach based on computer simulations for training to identify regions of the nuclear genome that were introgressed. We inferred asymmetric introgression of ≈2% of the genome, in the same direction that was observed for the mitochondrial genome. Because a previous study based on a reduced-representation sequencing approach was not able to detect this introgression, we conclude that whole-genome sequencing is necessary when the fraction of the introgressed genome is small. We also analysed the whole-genome sequence of a hybrid individual, demonstrating that hybrids have the capacity to backcross with the species for which virtually no introgression was observed. Our data suggest that one species has recently invaded the range of the other and/or some alleles that where transferred from the invaded into the invading species could be under positive selection and may have favoured the adaptation of the invading species to the Alpine environment.

Genomic data support management of anadromous Arctic Char fisheries in Nunavik by highlighting neutral and putatively adaptive genetic variation

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice-free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation-by-distance reflecting contemporary gene flow. Results from three gene-environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

Endosymbiont Capture, a Repeated Process of Endosymbiont Transfer with Replacement in Trypanosomatids Angomonas spp

Trypanosomatids of the subfamily Strigomonadinae bear permanent intracellular bacterial symbionts acquired by the common ancestor of these flagellates. However, the cospeciation pattern inherent to such relationships was revealed to be broken upon the description of Angomonas ambiguus, which is sister to A. desouzai, but bears an endosymbiont genetically close to that of A. deanei. Based on phylogenetic inferences, it was proposed that the bacterium from A. deanei had been horizontally transferred to A. ambiguus. Here, we sequenced the bacterial genomes from two A. ambiguus isolates, including a new one from Papua New Guinea, and compared them with the published genome of the A. deanei endosymbiont, revealing differences below the interspecific level. Our phylogenetic analyses confirmed that the endosymbionts of A. ambiguus were obtained from A. deanei and, in addition, demonstrated that this occurred more than once. We propose that coinfection of the same blowfly host and the phylogenetic relatedness of the trypanosomatids facilitate such transitions, whereas the drastic difference in the occurrence of the two trypanosomatid species determines the observed direction of this process. This phenomenon is analogous to organelle (mitochondrion/plastid) capture described in multicellular organisms and, thereafter, we name it endosymbiont capture.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.

Cytonuclear discordance in the crowned-sparrows, Zonotrichia atricapilla and Zonotrichia leucophrys

The golden-crowned (Zonotrichia atricapilla) and white-crowned (Z. leucophrys) sparrows have been presented as a compelling case for rapid speciation. They display divergence in song and plumage with overlap in their breeding ranges implying reproductive isolation, but have almost identical mitochondrial genomes. Previous research proposed hybridization and subsequent mitochondrial introgression as an alternate explanation, but lacked robust nuclear gene trees to distinguish between introgression and incomplete lineage sorting. We test for signatures of these processes between Z. atricapilla and Z. leucophrys, and investigate the relationships among Z. leucophrys subspecies, using mitochondrial sequencing and a reduced representation nuclear genomic dataset. Contrary to the paraphyly evident in mitochondrial gene trees, we confirmed the reciprocal monophyly of Z. atricapilla and Z. leucophrys using large panels of single nucleotide polymorphisms (SNPs). The pattern of cytonuclear discordance is consistent with limited, historical hybridization and mitochondrial introgression, rather than a recent origin and incomplete lineage sorting between recent sister species. We found evidence of nuclear phylogeographic structure within Z. leucophrys with two distinct clades. Altogether, our results indicate deeper divergences between Z. atricapilla and Z. leucophrys than inferred using mitochondrial markers. Our results demonstrate the limitations of relying solely on mitochondrial DNA for taxonomy, and raise questions about the possibility of selection on the mitochondrial genome during temperature oscillations (e.g. during the Pleistocene). Historical mitochondrial introgression facilitated by past environmental changes could cause erroneous dating of lineage splitting in other taxa when based on mitochondrial DNA alone.

Molecular Phylogeny of Cryptonanus (Didelphidae: Thylamyini): Evidence for a recent and complex diversification in South American open biomes

Systematic revisions of South American marsupials have contributed to our knowledge about genus and species diversity in the last decades, including studies of the most recently described genus Cryptonanus (Didelphidae), currently comprising four recognized species. Herein we provide the first phylogeny for these mouse opossums based on comprehensive sampling, including representatives from all nominal taxa, encompassing most of the geographic distribution of the genus while also extending its known range. The taxonomic status of Cryptonanus species was explored by analyses of multiple mitochondrial and nuclear DNA markers to assess phylogenetic relationships and to provide divergence time estimates, species delimitations and biogeographical hypotheses. Cryptonanus monophyly remained highly supported despite the inclusion of abundant new data from more than a hundred specimens, comprising 10 independent evolutionary lineages. Species-complexes within valid nominal taxa reveal higher species richness in the genus. Based on divergence estimates from a dated phylogeny, we suggest that Cryptonanus diversified along the Quaternary, with speciation events occurring well into the Pleistocene. The best supported biogeographical hypothesis endorses speciation by vicariance and subset speciation across open formations in shaping the evolutionary history of this didelphid genus, strongly associated with dry tropical landscapes of South America.

Evaluating evidence of mitonuclear incompatibilities with the sex chromosomes in an avian hybrid zone

The exploration of hybrid zones and the intergenomic conflicts exposed through hybridization provide windows into the processes of divergence and speciation. Sex chromosomes and mitonuclear incompatibilities have strong associations with the genetics of hybrid dysfunction. In ZW sex-determining systems, maternal co-inheritance of the mitochondrial and W chromosomes immediately exposes incompatibilities between these maternal contributions of one species and the Z chromosome of another. We analyze mitochondrial and Z chromosome admixture in the long-tailed finch (Poephila acuticauda) of Australia, where hybridizing subspecies differ prominently in Z chromosome genotype and in bill color, yet the respective centers of geographic admixture for these two traits are offset by 350 km. We report two well-defined mitochondrial clades that diverged ∼0.5 million years ago. Mitochondrial contact is geographically co-located within a hybrid zone of Z chromosome admixture and is displaced from bill color admixture by nearly 400 km. Consistent with Haldane's rule expectations, hybrid zone females are significantly less likely than males to carry an admixed Z chromosome or have mismatched Z-mitochondrial genotypes. Furthermore, there are significantly fewer than expected mitonuclear mismatches in hybrid zone females and paternal backcross males. Results suggest a potential for mitonuclear/sex chromosome incompatibilities in the emergence of reproductive isolation in this system.

Evolutionary history of the oriental fire-bellied toad (Bombina orientalis) in Northeast China

The evolutionary history of a species is generally affected by the combination of geological events and climate fluctuations. By analyzing the population features, genetic structure and the effective population historical dynamics of existing species, the population evolutionary history can be reestablished. In recent years, geological evidence shows that the Yilan-Yitong fault zone located in Northeast Asia experienced strong and frequent geological changes in the late Quaternary period. Species population history has been shaped by the combination of the complex climatic conditions of the Quaternary and Pleistocene glacial interglacial cycles and palaeogeological events in Northeast Asia and it has become a research focus for evolutionary biology researchers. In this study, mitochondrial and microsatellite molecular markers were used to reveal the population features, genetic structure, and the effective population historical dynamics of the Oriental fire-bellied toad (Bombina orientalis). The results showed that the strong seismic activity of the Yilan-Yitong fault zone in the late Quaternary period was the main reason for the population differentiation of Oriental fire-bellied toad in northeast China. The Quaternary Pleistocene glacial interglacial cycles led to the significant bottleneck effect of the western population located in the Maoer mountain area. As a result, the western population has low genetic diversity. Recent gene flow between eastern and western populations and historical evidence of population expansion proved that the dispersal behavior of the western populations was the main cause of the low genetic diversity and mitochondrial and nuclear discordance. Human economic activity may be the mainly driving factor. These evidences showed that the comprehensive influence of geology, climate, human activities and other factors should be considered in the process of exploring the evolutionary history of species.

Phylogeography and genetic diversity of the widespread katydid Ducetia japonica (Thunberg, 1815) across China

Habitat fragmentation can lower migration rates and genetic connectivity among remaining populations of native species. Ducetia japonica is one of the most widespread katydids in China, but little is known about its genetic structure and phylogeographic distribution. We combined the five-prime region of cytochrome c oxidase subunit I (COI-5P), 11 newly developed microsatellite loci coupled with an ecological niche model (ENM) to examine the genetic diversity and population structure of D. japonica in China and beyond to Laos and Singapore. Both Bayesian inference (BI) and haplotype network methods revealed six mitochondrial COI-5P lineages. The distribution of COI-5P haplotypes may not demonstrate significant phylogeographic structure (N ST > G ST, p > .05). The STRUCTURE analysis based on microsatellite data also revealed six genetic clusters, but discordant with those obtained from COI-5P haplotypes. For both COI-5P and microsatellite data, Mantel tests revealed a significant positive correlation between geographic and genetic distances in mainland China. Bayesian skyline plot (BSP) analyses indicated that the population size of D. japonica's three major mitochondrial COI-5P lineages were seemingly not affected by last glacial maximum (LGM, 0.015-0.025 Mya). The ecological niche models showed that the current distribution of D. japonica was similar to the species' distribution during the LGM period and only slightly extended in northern China. Further phylogeographic studies based on more extensive sampling are needed to identify specific locations of glacial refugia in northern China.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

Landscape resistance constrains hybridization across contact zones in a reproductively and morphologically polymorphic salamander

Explicitly accounting for phenotypic differentiation together with environmental heterogeneity is crucial to understand the evolutionary dynamics in hybrid zones. Species showing intra-specific variation in phenotypic traits that meet across environmentally heterogeneous regions constitute excellent natural settings to study the role of phenotypic differentiation and environmental factors in shaping the spatial extent and patterns of admixture in hybrid zones. We studied three environmentally distinct contact zones where morphologically and reproductively divergent subspecies of Salamandra salamandra co-occur: the pueriparous S. s. bernardezi that is mostly parapatric to its three larviparous subspecies neighbours. We used a landscape genetics framework to: (i) characterise the spatial location and extent of each contact zone; (ii) assess patterns of introgression and hybridization between subspecies pairs; and (iii) examine the role of environmental heterogeneity in the evolutionary dynamics of hybrid zones. We found high levels of introgression between parity modes, and between distinct phenotypes, thus demonstrating the evolution to pueriparity alone or morphological differentiation do not lead to reproductive isolation between these highly divergent S. salamandra morphotypes. However, we detected substantial variation in patterns of hybridization across contact zones, being lower in the contact zone located on a topographically complex area. We highlight the importance of accounting for spatial environmental heterogeneity when studying evolutionary dynamics of hybrid zones.

The Asian plethodontid salamander preserves historical genetic imprints of recent northern expansion

The Korean Peninsula, located at the southern tip of Northeast Asia, has never been covered by ice sheets and was a temperate refugium during the Pleistocene. Karsenia koreana, the sole Asian plethodontid salamander species, occurs only on the southern half of the Korean Peninsula and is thought to have found various climatic refugia. Despite its phylogenetic and biogeographic importance, no population-level genetic analysis has been performed on this species. Here we study the population genetic structure of K. koreana using mitochondrial and microsatellite loci to understand the recent historical dispersion process that shaped its current distribution. Overall, the genetic distance between populations correlated well with the spatial distance, and the genetic structure among populations showed signs of a unilateral northward expansion from a southernmost refugium population. Given the distinct genetic structure formed among the populations, the level of historical gene flow among populations appears to have been very low. As the estimated effective population size of K. koreana was also small, these results suggest that the small, restricted populations of K. koreana are extremely vulnerable to environmental changes that may require high levels of genetic diversity to cope with. Thus, special management strategies are needed to preserve these remnant populations.

Ring distribution patterns-diversification or speciation? Comparative phylogeography of two small mammals in the mountains surrounding the Sichuan Basin

Studying the genetic differentiation in a unique geographical area contributes to understanding the process of speciation. Here, we explore the spatial genetic structure and underlying formation mechanism of two congeneric small mammal species (Apodemus draco and A. chevrieri), which are mainly distributed in the mountains surrounding the lowland Sichuan Basin, southwest China. We applied a set of comparative phylogeographical analyses to determine their genetic diversification patterns, combining mitochondrial (Cytb and COI) and nuclear (microsatellite loci) markers, with dense sampling throughout the range (411 A. draco from 21 sites and 191 A. chevrieri from 22 sites). Moreover, we performed three complementary statistical methods to investigate the correlation between genotype and geographical and environmental components, and predicted the potential suitable distributional range under the present and historical climate conditions. Our results suggest that both species have experienced allopatric differentiation and admixture in historical periods, resulting in a ring-shape diversification, under the barrier effect of the Sichuan Basin. We infer that the tectonic events of the Qinghai-Tibetan Plateau and climatic oscillations during the Quaternary played an important role on the genetic divergence of the two species by providing environmental heterogeneity and geographical variation. Our study reveals a case of two sympatric small mammals following a ring-shaped diversification pattern and provides insight into the process of differentiation.

Reticulate evolution in Conidae: Evidence of nuclear and mitochondrial introgression

Conidae is a hyperdiverse family of marine snails that has many hallmarks of adaptive radiation. Hybridization and introgression may contribute to such instances of rapid diversification by generating novel gene combinations that facilitate exploitation of distinct niches. Here we evaluated whether or not these mechanisms may have contributed to the evolutionary history of a subgenus of Conidae (Virroconus). Several observations hint at evidence of past introgression for members of this group, including incongruence between phylogenetic relationships inferred from mitochondrial gene sequences and morphology and widespread sympatry of many Virroconus species in the Indo-West Pacific. We generated and analyzed transcriptome data of Virroconus species to (i) infer a robust nuclear phylogeny, (ii) assess mitochondrial and nuclear gene tree discordance, and (iii) formally test for introgression of nuclear loci. We identified introgression of mitochondrial genomes and nuclear gene regions between ancestors of one pair of Virroconus species, and mitochondrial introgression between another pair. We also found evidence of adaptive introgression of conotoxin venom loci between a third pair of species. Together, our results demonstrate that hybridization and introgression impacted the evolutionary history of Virroconus and hence may have contributed to the adaptive radiation of Conidae.

Redescription of Emplectonema viride - a ubiquitous intertidal hoplonemertean found along the West Coast of North America

Emplectonema viride Stimpson, 1857, a barnacle predator, is one of the most common and conspicuous intertidal nemerteans found along the West Coast of North America from Alaska to California, but it is currently referred to by the wrong name. Briefly described without designation of type material or illustrations, the species was synonymized with the Atlantic look-alike, Emplectonema gracile (Johnston, 1837) by Coe. Here we present morphological and molecular evidence that E. viride is distinct from E. gracile. The two species exhibit differences in color of live specimens and egg size and are clearly differentiated with species delimitation analyses based on sequences of the partial regions of the 16S rRNA and cytochrome c oxidase subunit I genes. In order to improve nomenclatural stability, we re-describe E. viride based on specimens from the southern coast of Oregon and discuss which species should be the type species of the genus. Emplectonema viride was one of the two species originally included in the genus Emplectonema Stimpson, 1857, but subsequent synonymization of E. viride with E. gracile resulted in acceptance of the Atlantic species, E. gracile, as the type species of the genus. We resurrect E. viride Stimpson, 1857 and following Corrêa's designation, this should be the type species of the genus Emplectonema.

Complex reticulate evolution of speckled brush-furred rats (Lophuromys) in the Ethiopian centre of endemism

The Ethiopian highlands represent a remarkable biodiversity 'hot spot' with a very high number of endemic species, even among vertebrates. Ethiopian representatives of a species complex of speckled brush-furred rats (Lophuromys flavopunctatus sensu lato) inhabit highland habitats ranging from low-elevation forests to Afroalpine grasslands. These may serve as a suitable model for understanding evolutionary processes leading to high genetic and ecological diversity in montane biodiversity hot spots. Here, we analyse the most comprehensive genetic data set of this group, comprising 315 specimens (all nine putative Ethiopian Lophuromys taxa sampled across most of their distribution ranges) genotyped at one mitochondrial and four nuclear markers, and thousands of SNPs from ddRAD sequencing. We performed phylogenetic analyses, delimited species and mapped their distribution and estimated divergence time between species (under the species-tree framework) and mitochondrial lineages. We found significant incongruence between mitochondrial and nuclear phylogenies, most probably caused by multiple interspecific introgression events. We discuss alternative scenarios of Ethiopian Lophuromys evolution, from retention of ancestral polymorphism to hybridization upon secondary contact of partially reproductively isolated lineages leading to reticulate evolution. Finally, we use the diversity of the speckled brush-furred rats for the description of the main biogeographic patterns in the fauna of the Ethiopian highlands.

Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds

Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.

Estimating hybridization as a consequence of climatic fluctuations: A case study of two geometridae species

In some cases, the phylogenetic analysis based on the mitochondrial DNA (mtDNA) and the nuclear DNA (ncDNA) are discordant. There are three major causes of the discordance within insects, including hybridization, incomplete lineage sorting (ILS) and infection by Wolbachia. In this study, we used a combination of multilocus and coalescent analyses to explore these processes occurred during the evolutionary history of Limbatochlamys rosthorni Rothschild, 1894 and Limbatochlamys pararosthorni Han and Xue, 2005. The ncDNA phylogenetic tree supported two reciprocally monophyletic species, whereas the mtDNA results failed to reveal such a structure and revealed an extensive level of admixture between two species. Because of very low Wolbachia infection rates (<20%), we firstly excluded this reason for the discordance. The fixed nucleotide differences and large genetic distances (1.5-2.5%) at the ncDNA genes suggested that the lineage sorting process between these two species is nearly complete and two species have experienced a prolonged period of independent evolution. Thus, we secondly excluded ILS. Sharing haplotypes, mtDNA gene flow occurring and the transitional samples with morphological features supported hybridization. The distribution contraction during glaciations and postglacial distribution expansion might have facilitated hybridization. Taken together, our study indicates that the current genetic structure of L. rosthorni and L. pararosthorni is the results of contraction and fragmentation into separated refugia during glaciations, followed by postglacial expansion and admixture.

Male-biased dispersal in a fungus-gardening ant symbiosis

For nearly all organisms, dispersal is a fundamental life-history trait that can shape their ecology and evolution. Variation in dispersal capabilities within a species exists and can influence population genetic structure and ecological interactions. In fungus-gardening (attine) ants, co-dispersal of ants and mutualistic fungi is crucial to the success of this obligate symbiosis. Female-biased dispersal (and gene flow) may be favored in attines because virgin queens carry the responsibility of dispersing the fungi, but a paucity of research has made this conclusion difficult. Here, we investigate dispersal of the fungus-gardening ant Trachymyrmex septentrionalis using a combination of maternally (mitochondrial DNA) and biparentally inherited (microsatellites) markers. We found three distinct, spatially isolated mitochondrial DNA haplotypes; two were found in the Florida panhandle and the other in the Florida peninsula. In contrast, biparental markers illustrated significant gene flow across this region and minimal spatial structure. The differential patterns uncovered from mitochondrial DNA and microsatellite markers suggest that most long-distance ant dispersal is male-biased and that females (and concomitantly the fungus) have more limited dispersal capabilities. Consequently, the limited female dispersal is likely an important bottleneck for the fungal symbiont. This bottleneck could slow fungal genetic diversification, which has significant implications for both ant hosts and fungal symbionts regarding population genetics, species distributions, adaptive responses to environmental change, and coevolutionary patterns.

Prevalence of Genetically Complex Leishmania Strains With Hybrid and Mito-Nuclear Discordance

Approximately 20 Leishmania species are known to cause cutaneous, mucocutaneous, and visceral disorders in humans. Identification of the causative species in infected individuals is important for appropriate treatment and a favorable prognosis because infecting species are known to be the major determinant of clinical manifestations and may affect treatments for leishmaniasis. Although Leishmania species have been conventionally identified by multilocus enzyme electrophoresis, genetic analysis targeting kinetoplast and nuclear DNA (kDNA and nDNA, respectively) is now widely used for this purpose. Recently, we conducted countrywide epidemiological studies of leishmaniasis in Ecuador and Peru to reveal prevalent species using PCR-RFLP targeting nDNA, and identified unknown hybrid parasites in these countries together with species reported previously. Furthermore, comparative analyses of kDNA and nDNA revealed the distribution of parasites with mismatches between these genes, representing the first report of mito-nuclear discordance in protozoa. The prevalence of an unexpectedly high rate (~10%) of genetically complex strains including hybrid strains, in conjunction with the observation of mito-nuclear discordance, suggests that genetic exchange may occur more frequently than previously thought in natural Leishmania populations. Hybrid Leishmania strains resulting from genetic exchanges are suggested to cause more severe clinical symptoms when compared with parental strains, and to have increased transmissibility by vectors of the parental parasite species. Therefore, it is important to clarify how such genetic exchange influences disease progression and transmissibility by sand flies in nature. In addition, our aim was to identify where and how the genetic exchange resulting in the formation of hybrid and mito-nuclear discordance occurs.