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.     

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.     

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.     

Using DNA taxonomy to investigate the ecological determinants of plankton diversity: explaining the occurrence of Synchaeta spp. (Rotifera,Monogononta) in mountain lakes

1. The occurrence of unresolved complexes of cryptic species may hinder the identification of the main ecological drivers of biodiversity when different cryptic taxa have different ecological requirements. 2. We assessed factors influencing the occurrence of Synchaeta species (monogonont rotifers) in 17 waterbodies of the Trentino‐South Tyrol region in the Eastern Alps. To do so, we compared the results of using unresolved complexes of cryptic species, as is common practice in limnological studies based on morphological taxonomy, and having resolved cryptic complexes, made possible by DNA taxonomy. 3. To identify cryptic species, we used the generalised mixed Yule coalescent (GMYC) model. We investigated the relationship between the environment and the occurrence of Synchaeta spp. by multivariate ordination using two definitions of the units of diversity, namely (i) unresolved species complexes (morphospecies) and (ii) putative cryptic species (GMYC entities). Our expectation was that resolving complexes of cryptic species could provide more information than using morphospecies. 4. As expected, DNA taxonomy provided greater taxonomic resolution than morphological taxonomy. Further, environmental‐based multivariate ordination on cryptic species explained a significantly higher proportion of variance than that based on morphospecies. Occurrence of GMYC entities was related to total phosphorus (TP), whereas no relationship could be found between morphospecies and the environment. Moreover, different cryptic species within the same morphospecies showed different, and even opposite, preferences for TP. In addition, the wide geographical distribution of haplotypes and cryptic species indicated the absence of barriers to dispersal in Synchaeta.     

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double‐digest restriction‐site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex‐biased asymmetries, NUMTs, natural selection, introgression or Wolbachia‐mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.     

Barcoding rotifer biodiversity in Mediterranean ponds using diapausing egg banks

The biodiversity of Mediterranean freshwater bodies is among the most threatened worldwide; therefore, its accurate estimation is an urgent issue. However, traditional methods are likely to underestimate freshwater zooplankton biodiversity due to its high species seasonality and cryptic diversity. We test the value of applying DNA barcoding to diapausing egg banks, in combination with the creation of a reference collection of DNA barcodes using adult individual samples, to characterize rotifer communities. We use monogonont rotifers from two lakes in Doñana National Park and one from Ruidera Natural Park in Spain as models to create a reference collection of DNA barcodes for taxonomically diagnosed adult individuals sampled from the water column, to compare with the sequences obtained from individual eggs from the diapausing egg banks. We apply two different approaches to carry out DNA taxonomy analyses, the generalized mixed Yule coalescent method (GMYC) and the Automatic Barcode Gap Discovery (ABGD), to the obtained sequences and to publicly available rotifer sequences. We obtained a total of 210 new rotifer COI sequences from all three locations (151 diapausing eggs and 59 adults). Both GMYC and ABGD generated the same 35 operational taxonomic units (OTUs), revealing four potential cryptic species. Most sequences obtained from diapausing eggs (85%) clustered with sequences obtained from morphologically diagnosed adults. Our approach, based on a single sediment sample, retrieved estimates of rotifer biodiversity higher than or similar to those of previous studies based on a number of seasonal samples. This study shows that DNA barcoding of diapausing egg banks is an effective aid to characterize rotifer diversity in Mediterranean freshwater bodies.     

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

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

Discordance between nuclear and mitochondrial genomes in sexual and asexual lineages of the freshwater snail Potamopyrgus antipodarum

The presence and extent of mitonuclear discordance in coexisting sexual and asexual lineages provides insight into 1) how and when asexual lineages emerged, and 2) the spatial and temporal scales at which the ecological and evolutionary processes influencing the evolution of sexual and asexual reproduction occur. Here, we used nuclear single‐nucleotide polymorphism (SNP) markers and a mitochondrial gene to characterize phylogeographic structure and the extent of mitonuclear discordance in Potamopyrgus antipodarum. This New Zealand freshwater snail is often used to study the evolution and maintenance of sex because obligately sexual and obligately asexual individuals often coexist. While our data indicate that sexual and asexual P. antipodarum sampled from the same lake population are often genetically similar, suggesting recent origin of these asexuals from sympatric sexual P. antipodarum, we also found significantly more population structure in sexuals vs. asexuals. This latter result suggests that some asexual lineages originated in other lakes and/or in the relatively distant past. When comparing mitochondrial and nuclear population genetic structure, we discovered that one mitochondrial haplotype (‘1A’) was rare in sexuals, but common and widespread in asexuals. Haplotype 1A frequency and nuclear genetic diversity were not associated, suggesting that the commonness of this haplotype cannot be attributed entirely to genetic drift and pointing instead to a role for selection.     

Species delimitation of the white‐tailed rubythroat Calliope pectoralis complex (Aves,Muscicapidae) using an integrative taxonomic approach

Our knowledge of the systematics and taxonomy of Asian birds has improved much in the last two decades, and the number of recognised species has increased significantly as a result of in‐depth studies using an integrative taxonomic approach. The Sino‐Himalayan mountains harbor a high level of passerine diversity. Several allopatric or parapatric taxa that are currently treated as subspecies of polytypic species within that region are likely to deserve full species status, and thus their taxonomic status needs to be revisited. Based on analyses of multilocus data, vocalizations and morphology, we propose that the white‐tailed rubythroat Calliope pectoralis should be treated as two species, the Himalayan rubythroat C. pectoralis sensu stricto in the Tian Shan and Himalayan mountains, and the Chinese rubythroat C. tschebaiewi in the mountains of southwestern and north‐central China. According to our dating analyses based on mitochondrial loci, these two species diverged approximately 2.2 million yr ago. We further found that C. tschebaiewi was paraphyletic to C. pectoralis sensu stricto in nuclear data, which demonstrates a state of mitonuclear discordance that warrants further work. Our results suggest that geographic changes and glacial cycles in the Pleistocene may have caused allopatric divergence in the C. pectoralis complex. Our study stresses the importance of applying an integrative taxonomy approach to fully unravel the true avian diversity in the Sino‐Himalayan Mountains.     

The toad’s warts: Discordance creates bumpy expectations of mitochondrial‐nuclear evolution between species

Discordance between the mitochondrial and nuclear genomes is a prevalent phenomenon in nature, in which the underlying processes responsible are considered to be important in shaping genetic variation in natural populations. Among the evolutionary processes that best explain such genomic mismatches incomplete lineage sorting and introgression are commonly identified, however, many studies are unable to distinguish between these hypotheses, which has become a major challenge in the field. In this issue of Molecular Ecology, Firneno et al. (2020) present an elegant exploration of mitochondrial‐nuclear discordance in Mesoamerican toads. Integrating genome‐scale and spatial data to test between these hypotheses within an empirical model testing framework, they find strong support that incomplete lineage sorting explains the observed discordance. Their work, along with many previous articles in Molecular Ecology, highlights the commonality of mito‐nuclear discordance among species despite the expectations of tightly concerted mitochondrial and nuclear genome evolution. It is increasingly clear that the nuclear genomes of many species are (at least for short periods of evolutionary time) functionally compatible with multiple, divergent mitochondrial haplotypes. As such, we suggest future research not only seeks to understand the processes causing spatial mito‐nuclear discordance (e.g. incomplete lineage sorting, introgression), but also explores those that maintain discordance through time and space (e.g. relaxed selection on mito‐nuclear interactions, heterozygosity, population demographics). We also discuss the vital role that taxonomy plays in interpreting patterns of mito‐nuclear discordance when data‐consistent yet differing taxonomies are used, such as treating allopatrically distributed taxa as multiple isolated populations versus multiple micro‐endemic species.     

Positive and purifying selection in mitochondrial genomes of a bird with mitonuclear discordance

Diversifying selection on metabolic pathways can reduce intraspecific gene flow and promote population divergence. An opportunity to explore this arises from mitonuclear discordance observed in an Australian bird Eopsaltria australis. Across >1500 km, nuclear differentiation is low and latitudinally structured by isolation by distance, whereas two highly divergent, parapatric mitochondrial lineages (>6.6% in ND2) show a discordant longitudinal geographic pattern and experience different climates. Vicariance, incomplete lineage sorting and sex‐biased dispersal were shown earlier to be unlikely drivers of the mitonuclear discordance; instead, natural selection on a female‐linked trait was the preferred hypothesis. Accordingly, here we tested for signals of positive, divergent selection on mitochondrial genes in E. australis. We used codon models and physicochemical profiles of amino acid replacements to analyse complete mitochondrial genomes of the two mitochondrial lineages in E. australis, its sister species Eopsaltria griseogularis, and outgroups. We found evidence of positive selection on at least five amino acids, encoded by genes of two oxidative phosphorylation pathway complexes NADH dehydrogenase (ND4 and ND4L) and cytochrome bc₁ (cyt‐b) against a background of widespread purifying selection on all mitochondrial genes. Three of these amino acid replacements were fixed in ND4 of the geographically most widespread E. australis lineage. The other two replacements were fixed in ND4L and cyt‐b of the geographically more restricted E. australis lineage. We discuss whether this selection may reflect local environmental adaptation, a by‐product of other selective processes, or genetic incompatibilities, and propose how these hypotheses can be tested in future.     

Widespread decoupling of mtDNA variation and species integrity in Grammia tiger moths (Lepidoptera: Noctuidae)

Abstract. We investigate the diversity of the North American tiger moth genus Grammia Rambur (Lepidoptera: Noctuidae) by comparing mitochondrial DNA (mtDNA) ‘barcode’ fragments of cytochrome oxidase I with non‐molecular characters such as morphology, ecology, behaviour and distribution. Mitochondrial DNA genealogy is strikingly at odds with morpho‐species taxonomy for most of the 28 sampled species, as haplotypic polyphyly not only is taxonomically widespread, but involves multiple shared haplotypes among two to four species. Morpho‐ecological traits show that those species sharing haplotypes are often not closely related. Furthermore, high mtDNA divergences occur within species. Haplotypic variation is highly discordant with species taxonomy, but variation at a continental scale reveals significant geographic structuring of haplogroups, transcending morpho‐species boundaries. A nested clade analysis and comparison of non‐molecular with mtDNA data indicate that most discordance between mtDNA and taxonomy in Grammia is explained best by taxonomically and geographically widespread ongoing hybridization events resulting in mtDNA introgression. We hypothesize that broad areas of sympatry, interspecifically compatible genitalic structure, and species overlap in pheromone components facilitate hybridization, with disparate interspecies abundances promoting mitochondrial introgression. The molecular evolution of Grammia challenges the view that interspecific gene exchange occurs rarely and is restricted to recently diverged species. These results show the value of mtDNA in detecting cryptic hybridization, while highlighting the inherent dangers of drawing taxonomic conclusions based solely on mtDNA.     

Molecular phylogeny and DNA barcoding confirm cryptic species in the African freshwater oyster Etheria elliptica Lamarck, 1807 (Bivalvia: Etheriidae)

Recent molecular approaches to taxonomy have led to a steady increase in the identification of cryptic species. Within the Etheriidae, the species Etheria elliptica (freshwater oyster) is widespread and common and exists in most of the major African drainages. Within the African freshwater ecosystems, there are major threats to biodiversity and cryptic species complicate conservation strategies; unknown species exist and no conservation status has been assigned. Our objective here was to determine if E. elliptica from several locations in the Congo drainage are correctly classified as representing a single species. We analysed the genetic diversity at two mitochondrial loci (COI and 16S) and two nuclear loci (H3 and 28S), and estimated evolutionary relationships using phylogenetic and DNA barcoding techniques. Bayesian inference yielded three cryptic species of Etheria, and mismatch analysis revealed discrete differences between the cryptic species. We identified three cryptic species within these collections, and evidence indicates that the third species may resolve further with more sampling. In conclusion, the taxonomic history of E. elliptica makes finding cryptic species unsurprising. However, molecular studies such as this may finally help to resolve the number of species within this genus.     

Extreme levels of hidden diversity in microscopic animals (Rotifera) revealed by DNA taxonomy

Knowledge and estimates of species richness at all scales are biased both by our understanding of the evolutionary processes shaping diversity and by the methods used to delineate the basic units of diversity. DNA taxonomy shows that diversity may be underestimated by traditional taxonomy, especially for microscopic animals. The effects of such hidden diversity are usually overlooked in ecological studies. Here, we estimate hidden diversity in bdelloid rotifers, a group of microscopic animals. We analyse cryptic diversity using a coalescent approach to infer taxonomical units from phylogenetic trees. Cryptic diversity was measured for eight traditional species of bdelloid rotifers and the results compared to that of the monogonont rotifer Brachionus plicatilis species complex, which is well studied and for which cross-mating experiments have been performed to explicitly define some of the species boundaries. A taxonomic inflation of up to 34 potential cryptic taxa was found in bdelloids. Cryptic taxa within each traditional species may be spatially isolated, but do not have narrower ecological niches. The species deemed as generalists exhibit the highest cryptic diversity. Cryptic diversity based on molecular characterization is commonly found in animals; nevertheless, the amount of cryptic diversity in bdelloids is much higher than in other groups analysed so far, maybe because of their peculiar parthenogenetic reproduction, other than microscopic size. We discuss this hypothesis in the light of the available empirical evidence from other groups of microscopic invertebrates, such as tardigrades and mites, which share size, habitat heterogeneity, potential for dispersal, and/or parthenogenetic reproduction.     

Cryptic diversity with wide salinity tolerance in the putative euryhaline Testudinella clypeata (Rotifera,Monogononta)

Aquatic faunas in fresh, brackish, and salt waters are usually well defined and differ amongst these three habitats. Nonetheless, some animals are known to be euryhaline, namely present across wide salinity ranges. The wide tolerance of putative euryhaline species has, however been refuted in some cases by DNA taxonomy, which has uncovered cryptic diversity with narrow ecological niches. We aim to improve knowledge on the putative euryhalinism of microinvertebrates and test whether it might actually be a real phenomenon or if euryhaline species are mostly a consequence of our previous inability to identify cryptic species with narrow salinity ranges, as discovered in Brachionus plicatilis. Using morphological analyses and DNA taxonomy, we investigated the species reality and distribution of a putative euryhaline rotifer species, Testudinella clypeata, and evaluated whether cryptic species are ecologically and/or geographically segregated. Different DNA taxonomy approaches concurred in revealing the presence of seven cryptic species within the T. clypeata morphospecies, which, in contrast to what has been previously detected, are actually euryhaline. Moreover, differences in analysed morphological traits were not significantly different amongst cryptic species. This suggests that DNA taxonomy improves our estimates of the actual diversity of microscopic species, in contrast to the morphological approach. © 2013 The Linnean Society of London     

Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.)

Distinct genetic markers should show similar patterns of differentiation between species reflecting their common evolutionary histories, yet there are increasing examples of differences in the biogeographic distribution of species‐specific nuclear (nuDNA) and mitochondrial DNA (mtDNA) variants within and between species. Identifying the evolutionary processes that underlie these anomalous patterns of genetic differentiation is an important goal. Here, we analyse the putative mitonuclear discordance observed between sister species of mole salamanders (Ambystoma barbouri and A. texanum) in which A. barbouri‐specific mtDNA is found in animals located within the range of A. texanum. We test three hypotheses for this discordance (undetected range expansion, mtDNA introgression, and hybridization) using nuDNA and mtDNA data analysed with methods that varied in the parameters estimated and the timescales measured. Results from a Bayesian clustering technique (structure ), bidirectional estimates of gene flow (migrate ‐n and IMa2) and phylogeny‐based methods (*beast , buck y) all support the conclusion that the discordance is due to geographically restricted mtDNA introgression from A. barbouri into A. texanum. Limited data on species‐specific tooth morphology match this conclusion. Significant differences in environmental conditions exist between sites where A. texanum with and without A. barbouri‐like mtDNA occur, suggesting a possible role for selection in the process of introgression. Overall, our study provides a general example of the value of using complimentary analyses to make inferences of the directionality, timescale, and source of mtDNA introgression in animals.     

Integrative insect taxonomy based on morphology,mitochondrial DNA,and hyperspectral reflectance profiling

Integrative taxonomy is considered a reliable taxonomic approach of closely related and cryptic species by integrating different sources of taxonomic data (genetic, ecological, and morphological characters). In order to infer the boundaries of seven species of the evacanthine leafhopper genus Bundera Distant, 1908 (Hemiptera: Cicadellidae), an integrated analysis based on morphology, mitochondrial DNA, and hyperspectral reflectance profiling (37 spectral bands from 411–870 nm) was conducted. Despite their morphological similarities, the genetic distances of the cytochrome c oxidase subunit I (COI) gene among the tested species are relatively large (5.8–17.3%). The species‐specific divergence of five morphologically similar species (Bundera pellucida and Bundera spp. 1–4) was revealed in mitochondrial DNA data and reflectance profiling. A key to identifying males is provided, and their morphological characters are described. Average reflectance profiles from the dorsal side of specimens were classified based on linear discriminant analysis. Cross‐validation of reflectance‐based classification revealed that the seven species could be distinguished with 91.3% classification accuracy. This study verified the feasibility of using hyperspectral imaging data in insect classification, and our work provides a good example of using integrative taxonomy in studies of closely related and cryptic species. © 2015 The Linnean Society of London     

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

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

Conservative and compensatory evolution in oxidative phosphorylation complexes of angiosperms with highly divergent rates of mitochondrial genome evolution

Interactions between nuclear and mitochondrial gene products are critical for eukaryotic cell function. Nuclear genes encoding mitochondrial‐targeted proteins (N‐mt genes) experience elevated rates of evolution, which has often been interpreted as evidence of nuclear compensation in response to elevated mitochondrial mutation rates. However, N‐mt genes may be under relaxed functional constraints, which could also explain observed increases in their evolutionary rate. To disentangle these hypotheses, we examined patterns of sequence and structural evolution in nuclear‐ and mitochondrial‐encoded oxidative phosphorylation proteins from species in the angiosperm genus Silene with vastly different mitochondrial mutation rates. We found correlated increases in N‐mt gene evolution in species with fast‐evolving mitochondrial DNA. Structural modeling revealed an overrepresentation of N‐mt substitutions at positions that directly contact mutated residues in mitochondrial‐encoded proteins, despite overall patterns of conservative structural evolution. These findings support the hypothesis that selection for compensatory changes in response to mitochondrial mutations contributes to the elevated rate of evolution in N‐mt genes. We discuss these results in light of theories implicating mitochondrial mutation rates and mitonuclear coevolution as drivers of speciation and suggest comparative and experimental approaches that could take advantage of heterogeneity in rates of mtDNA evolution across eukaryotes to evaluate such theories.     

Sex‐specific mitonuclear epistasis and the evolution of mitochondrial bioenergetics,ageing, and life history in seed beetles

The role of mitochondrial DNA for the evolution of life‐history traits remains debated. We examined mitonuclear effects on the activity of the multisubunit complex of the electron transport chain (ETC) involved in oxidative phosphorylation (OXPHOS) across lines of the seed beetle Acanthoscelides obtectus selected for a short (E) or a long (L) life for more than >160 generations. We constructed and phenotyped mitonuclear introgression lines, which allowed us to assess the independent effects of the evolutionary history of the nuclear and the mitochondrial genome. The nuclear genome was responsible for the largest share of divergence seen in ageing. However, the mitochondrial genome also had sizeable effects, which were sex‐specific and expressed primarily as epistatic interactions with the nuclear genome. The effects of mitonuclear disruption were largely consistent with mitonuclear coadaptation. Variation in ETC activity explained a large proportion of variance in ageing and life‐history traits and this multivariate relationship differed somewhat between the sexes. In conclusion, mitonuclear epistasis has played an important role in the laboratory evolution of ETC complex activity, ageing, and life histories and these are closely associated. The mitonuclear architecture of evolved differences in life‐history traits and mitochondrial bioenergetics was sex‐specific.