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

Evolution and maintenance of mtDNA gene content across eukaryotes

Across eukaryotes, most genes required for mitochondrial function have been transferred to, or otherwise acquired by, the nucleus. Encoding genes in the nucleus has many advantages. So why do mitochondria retain any genes at all? Why does the set of mtDNA genes vary so much across different species? And how do species maintain functionality in the mtDNA genes they do retain? In this review, we will discuss some possible answers to these questions, attempting a broad perspective across eukaryotes. We hope to cover some interesting features which may be less familiar from the perspective of particular species, including the ubiquity of recombination outside bilaterian animals, encrypted chainmail-like mtDNA, single genes split over multiple mtDNA chromosomes, triparental inheritance, gene transfer by grafting, gain of mtDNA recombination factors, social networks of mitochondria, and the role of mtDNA dysfunction in feeding the world. We will discuss a unifying picture where organismal ecology and gene-specific features together influence whether organism X retains mtDNA gene Y, and where ecology and development together determine which strategies, importantly including recombination, are used to maintain the mtDNA genes that are retained.

Accelerated mitochondrial evolution and asymmetric fitness of hybrids contribute to the persistence of Helix thessalica in the Helix pomatia range

Interbreeding and introgression between recently diverged species is common. However, the processes that prevent these species from merging where they co-occur are not well understood. We studied the mechanisms that allowed an isolated group of populations of the snail Helix thessalica to persist within the range of the related Helix pomatia despite high gene flow. Using genomic cline analysis, we found that the nuclear gene flow between the two taxa across the mosaic hybrid zone was not different from that expected under neutral admixture, but that the exchange of mtDNA was asymmetric. Tests showed that there is relaxed selection in the mitochondrial genome of H. thessalica and that the substitution rate is elevated compared to that of H. pomatia. A lack of hybrids that combine the mtDNA of H. thessalica with a mainly (>46%) H. pomatia genomic background indicates that the nuclear-encoded mitochondrial proteins of H. pomatia are not well adapted to the more rapidly evolving proteins and RNAs encoded by the mitochondrion of H. thessalica. The presumed reduction of fitness of hybrids with the fast-evolving mtDNA of H. thessalica and a high H. pomatia ancestry, similar to 'Darwin's Corollary to Haldane's rule', resulted in a relative loss of H. pomatia nuclear ancestry compared to H. thessalica ancestry in the hybrid zone. This probably prevents the H. thessalica populations from merging quickly with the surrounding H. pomatia populations and supports the hypothesis that incompatibilities between rapidly evolving mitochondrial genes and nuclear genes contribute to speciation.

Phylogeography and paleoclimatic range dynamics explain variable outcomes to contact across a species' range

Replicability of divergence after contact is a poorly characterized process, particularly in the contexts of phylogeography and postglacial range dynamics within species. Using contact zones located at the leading-, mid- and rear-edges of a species' range, we examined variation in outcomes to contact between divergent lineages of Campanula americana. We investigated whether contact zones vary in quantity and directionality of gene flow, how phylogeographic structure differs between contact zones, and how historic range dynamics may affect outcomes to contact. We found that all contact zones formed at similar times via primary contact yet detected significant admixture in only the rear-edge (RE) contact zone. In the northern leading-edge contact zone and the mid-range Virginia contact zone, gene flow was minimal and asymmetric. In the southern RE contact zone, gene flow was strong and symmetric. Asymmetric admixture in the leading-edge and Virginia contact zones matches the directionality of a known cosmopolitan cytonuclear incompatibility between lineages of C. americana. Our results emphasize the dependence of speciation processes on phylogeographic structure, evolutionary history and range dynamics.

GENETIC AND TRAIT VARIABILITY OF GYRINICOLA REVEALS THE EXISTENCE OF AT LEAST FOUR SPECIES WITHIN THE UNITED STATES

The tadpole-dwelling pinworm, Gyrinicola batrachiensis (Walton, 1929) Adamson, 1981 was recognized as the sole representative of the genus across Canada and the United States. However, evaluation of the morphology of these parasites across their range revealed considerable morphological variability that suggested diagnosable morphotypes. These morphotypes were associated with different species of anurans, several of which occurred in sympatry. Herein we use an extensive geographic sampling across the United States to obtain the morphotypes, screen their genetic diversity, and analyze this information using an integrative approach. We reconstructed their phylogeny using nuclear ribosomal partial genes 18S and 28S, ITS1, 5.8S, and ITS2, as well as 5 mitochondrial genes generated with Next-Generation sequencing technology. This phylogeny reveals 3 well-resolved lineages, which upon the use of a statistical approach (bPTP [Bayesian implementation of the Poisson tree processes]) supports the delimitation of 4 distinct groups equivalent to species. These putative species groups were tested using morphological characteristics paired with a MANOVA and canonical variate analysis. Results suggest that at least 4 species of Gyrinicola are present within North America, resulting in the resurrection of G. armatus (Walton, 1933) and the description of 2 new species.

Mitonuclear interactions impact aerobic metabolism in hybrids and may explain mitonuclear discordance in young, naturally hybridizing bird lineages

Understanding genetic incompatibilities and genetic introgression between incipient species are major goals in evolutionary biology. Mitochondrial genes evolve rapidly and exist in dense gene networks with coevolved nuclear genes, suggesting that mitochondrial respiration may be particularly susceptible to disruption in hybrid organisms. Mitonuclear interactions have been demonstrated to contribute to hybrid dysfunction between deeply divergent taxa crossed in the laboratory, but there are few empirical examples of mitonuclear interactions between younger lineages that naturally hybridize. Here, we use controlled hybrid crosses and high-resolution respirometry to provide the first experimental evidence in a bird that inter-lineage mitonuclear interactions impact mitochondrial aerobic metabolism. Specifically, respiration capacity of the two mitodiscordant backcrosses (with mismatched mitonuclear combinations) differs from one another, although they do not differ significantly from the parental groups or mitoconcordant backcrosses as we would expect of mitonuclear disruptions. In the wild hybrid zone between these subspecies, the mitochondrial cline centre is shifted west of the nuclear cline centre, which is consistent with the direction of our experimental results. Our results therefore demonstrate asymmetric mitonuclear interactions that impact the capacity of cellular mitochondrial respiration and may help to explain the geographic discordance between mitochondrial and nuclear genomes observed in the wild.

Genome Evolution and Introgression in the New Zealand mud Snails Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa

We have sequenced, assembled, and analyzed the nuclear and mitochondrial genomes and transcriptomes of Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa, two prosobranch snail species native to New Zealand that together span the continuum from estuary to freshwater. These two species are the closest known relatives of the freshwater species Potamopyrgus antipodarum-a model for studying the evolution of sex, host-parasite coevolution, and biological invasiveness-and thus provide key evolutionary context for understanding its unusual biology. The P. estuarinus and P. kaitunuparaoa genomes are very similar in size and overall gene content. Comparative analyses of genome content indicate that these two species harbor a near-identical set of genes involved in meiosis and sperm functions, including seven genes with meiosis-specific functions. These results are consistent with obligate sexual reproduction in these two species and provide a framework for future analyses of P. antipodarum-a species comprising both obligately sexual and obligately asexual lineages, each separately derived from a sexual ancestor. Genome-wide multigene phylogenetic analyses indicate that P. kaitunuparaoa is likely the closest relative to P. antipodarum. We nevertheless show that there has been considerable introgression between P. estuarinus and P. kaitunuparaoa. That introgression does not extend to the mitochondrial genome, which appears to serve as a barrier to hybridization between P. estuarinus and P. kaitunuparaoa. Nuclear-encoded genes whose products function in joint mitochondrial-nuclear enzyme complexes exhibit similar patterns of nonintrogression, indicating that incompatibilities between the mitochondrial and the nuclear genome may have prevented more extensive gene flow between these two species.

Mitogenomes do not substantially improve phylogenetic resolution in a young non-model adaptive radiation of freshwater gastropods

BACKGROUND

Species flocks in ancient lakes, and particularly those arising from adaptive radiation, make up the bulk of overall taxonomic and morphological diversity in these insular ecosystems. For these mostly young species assemblages, classical mitochondrial barcoding markers have so far been key to disentangle interspecific relationships. However, with the rise and further development of next-generation sequencing (NGS) methods and mapping tools, genome-wide data have become an increasingly important source of information even for non-model groups.

RESULTS

Here, we provide, for the first time, a comprehensive mitogenome dataset of freshwater gastropods endemic to Sulawesi and thus of an ancient lake invertebrate species flock in general. We applied low-coverage whole-genome sequencing for a total of 78 individuals including 27 out of the 28 Tylomelania morphospecies from the Malili lake system as well as selected representatives from Lake Poso and adjacent catchments. Our aim was to assess whether mitogenomes considerably contribute to the phylogenetic resolution within this young species flock. Interestingly, we identified a high number of variable and parsimony-informative sites across the other 'non-traditional' mitochondrial loci. However, although the overall support was very high, the topology obtained was largely congruent with previously published single-locus phylogenies. Several clades remained unresolved and a large number of species was recovered polyphyletic, indicative of both rapid diversification and mitochondrial introgression.

CONCLUSIONS

This once again illustrates that, despite the higher number of characters available, mitogenomes behave like a single locus and thus can only make a limited contribution to resolving species boundaries, particularly when introgression events are involved.

Friend turned foe: selfish behavior of a spontaneously arising mitochondrial deletion in an experimentally evolved Caenorhabditis elegans population

Selfish mitochondrial DNA (mtDNA) mutations are variants that can proliferate within cells and enjoy a replication or transmission bias without fitness benefits for the host. mtDNA deletions in Caenorhabditis elegans can reach high heteroplasmic frequencies despite significantly reducing fitness, illustrating how new mtDNA variants can give rise to genetic conflict between different levels of selection and between the nuclear and mitochondrial genomes. During a mutation accumulation experiment in C. elegans, a 1,034-bp deletion originated spontaneously and reached an 81.7% frequency within an experimental evolution line. This heteroplasmic mtDNA deletion, designated as meuDf1, eliminated portions of 2 protein-coding genes (coxIII and nd4) and tRNA-thr in entirety. mtDNA copy number in meuDf1 heteroplasmic individuals was 35% higher than in individuals with wild-type mitochondria. After backcrossing into a common genetic background, the meuDf1 mitotype was associated with reduction in several fitness traits and independent competition experiments found a 40% reduction in composite fitness. Experiments that relaxed individual selection by single individual bottlenecks demonstrated that the deletion-bearing mtDNA possessed a strong transmission bias, thereby qualifying it as a novel selfish mitotype.

Conservation Mitonuclear Replacement: Facilitated mitochondrial adaptation for a changing world

Most species will not be able to migrate fast enough to cope with climate change, nor evolve quickly enough with current levels of genetic variation. Exacerbating the problem are anthropogenic influences on adaptive potential, including the prevention of gene flow through habitat fragmentation and the erosion of genetic diversity in small, bottlenecked populations. Facilitated adaptation, or assisted evolution, offers a way to augment adaptive genetic variation via artificial selection, induced hybridization, or genetic engineering. One key source of genetic variation, particularly for climatic adaptation, are the core metabolic genes encoded by the mitochondrial genome. These genes influence environmental tolerance to heat, drought, and hypoxia, but must interact intimately and co-evolve with a suite of important nuclear genes. These coadapted mitonuclear genes form some of the important reproductive barriers between species. Mitochondrial genomes can and do introgress between species in an adaptive manner, and they may co-introgress with nuclear genes important for maintaining mitonuclear compatibility. Managers should consider the relevance of mitonuclear genetic variability in conservation decision-making, including as a tool for facilitating adaptation. I propose a novel technique dubbed Conservation Mitonuclear Replacement (CmNR), which entails replacing the core metabolic machinery of a threatened species-the mitochondrial genome and key nuclear loci-with those from a closely related species or a divergent population, which may be better-adapted to climatic changes or carry a lower genetic load. The most feasible route to CmNR is to combine CRISPR-based nuclear genetic editing with mitochondrial replacement and assisted reproductive technologies. This method preserves much of an organism's phenotype and could allow populations to persist in the wild when no other suitable conservation options exist. The technique could be particularly important on mountaintops, where rising temperatures threaten an alarming number of species with almost certain extinction in the next century.

Cytonuclear Interactions and Subgenome Dominance Shape the Evolution of Organelle-Targeted Genes in the Brassica Triangle of U

The interaction and coevolution between nuclear and cytoplasmic genomes are one of the fundamental hallmarks of eukaryotic genome evolution and, 2 billion yr later, are still major contributors to the formation of new species. Although many studies have investigated the role of cytonuclear interactions following allopolyploidization, the relative magnitude of the effect of subgenome dominance versus cytonuclear interaction on genome evolution remains unclear. The Brassica triangle of U features 3 diploid species that together have formed 3 separate allotetraploid species on similar evolutionary timescales, providing an ideal system for understanding the contribution of the cytoplasmic donor to hybrid polyploid. Here, we investigated the evolutionary pattern of organelle-targeted genes in Brassica carinata (BBCC) and 2 varieties of Brassica juncea (AABB) at the whole-genome level, with particular focus on cytonuclear enzyme complexes. We found partial evidence that plastid-targeted genes experience selection to match plastid genomes, but no obvious corresponding signal in mitochondria-targeted genes from these 2 separately formed allopolyploids. Interestingly, selection acting on plastid genomes always reduced the retention rate of plastid-targeted genes encoded by the B subgenome, regardless of whether the Brassica nigra (BB) subgenome was contributed by the paternal or maternal progenitor. More broadly, this study illustrates the distinct selective pressures experienced by plastid- and mitochondria-targeted genes, despite a shared pattern of inheritance and natural history. Our study also highlights an important role for subgenome dominance in allopolyploid genome evolution, even in genes whose function depends on separately inherited molecules.

Mitochondrial genetic variation as a potential mediator of intraspecific behavioural diversity

Mitochondrial genes play an essential role in energy metabolism. Variation in the mitochondrial DNA (mtDNA) sequence often exists within species, and this variation can have consequences for energy production and organismal life history. Yet, despite potential links between energy metabolism and the expression of animal behaviour, mtDNA variation has been largely neglected to date in studies investigating intraspecific behavioural diversity. We outline how mtDNA variation and interactions between mitochondrial and nuclear genotypes may contribute to the expression of individual-to-individual behavioural differences within populations, and why such effects may lead to sex differences in behaviour. We contend that integration of the mitochondrial genome into behavioural ecology research may be key to fully understanding the evolutionary genetics of animal behaviour.

A lethal mitonuclear incompatibility in complex I of natural hybrids

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1-3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4-6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Association Analysis Provides Insights into Plant Mitonuclear Interactions

Cytonuclear interaction refers to the complex and ongoing process of coevolution between nuclear and organelle genomes, which are responsible for cellular respiration, photosynthesis, lipid metabolism, etc. and play a significant role in adaptation and speciation. There have been a large number of studies to detect signatures of cytonuclear interactions. However, identification of the specific nuclear and organelle genetic polymorphisms that are involved in these interactions within a species remains relatively rare. The recent surge in whole genome sequencing has provided us an opportunity to explore cytonuclear interaction from a population perspective. In this study, we analyzed a total of 3,439 genomes from 7 species to identify signals of cytonuclear interactions by association (linkage disequilibrium) analysis of variants in both the mitochondrial and nuclear genomes across flowering plants. We also investigated examples of nuclear loci identified based on these association signals using subcellular localization assays, gene editing, and transcriptome sequencing. Our study provides a novel perspective on the investigation of cytonuclear coevolution, thereby enriching our understanding of plant fitness and offspring sterility.

Mitochondrial function and sexual selection: can physiology resolve the 'lek paradox'?

Across many taxa, males use elaborate ornaments or complex displays to attract potential mates. Such sexually selected traits are thought to signal important aspects of male 'quality'. Female mating preferences based on sexual traits are thought to have evolved because choosy females gain direct benefits that enhance their lifetime reproductive success (e.g. greater access to food) and/or indirect benefits because high-quality males contribute genes that increase offspring fitness. However, it is difficult to explain the persistence of female preferences when males only provide genetic benefits, because female preferences should erode the heritable genetic variation in fitness that sexually selected traits signal. This 'paradox of the lek' has puzzled evolutionary biologists for decades, and inspired many hypotheses to explain how heritable variation in sexually selected traits is maintained. Here, we discuss how factors that affect mitochondrial function can maintain variation in sexually selected traits despite strong female preferences. We discuss how mitochondrial function can influence the expression of sexually selected traits, and we describe empirical studies that link the expression of sexually selected traits to mitochondrial function. We explain how mothers can affect mitochondrial function in their offspring by (a) influencing their developmental environment through maternal effects and (b) choosing a mate to increase the compatibility of mitochondrial and nuclear genes (i.e. the 'mitonuclear compatibility model of sexual selection'). Finally, we discuss how incorporating mitochondrial function into models of sexual selection might help to resolve the paradox of the lek, and we suggest avenues for future research.

Multiple contact zones and karyotypic evolution in a neotropical frog species complex

Previous studies of DNA sequence and karyotypic data have revealed high genetic diversity in the Physalaemus cuvieri - Physalaemus ephippifer species complex-a group of small leptodactylid frogs in South America. To date, seven major genetic lineages have been recognized in this group, with species delimitation tests supporting four to seven of them as valid species. Among these, only P. ephippifer shows heteromorphic sex chromosomes, but the implications of cytogenetic divergence for the evolution of this group are unknown. We analyzed karyotypic, mitochondrial DNA, and 3RAD genomic data to characterize a putative contact zone between P. ephippifer and P. cuvieri Lineage 1, finding evidence for admixture and karyotypic evolution. We also describe preliminary evidence for admixture between two other members of this species complex-Lineage 1 and Lineage 3 of P. cuvieri. Our study sheds new light on evolutionary relationships in the P. cuvieri - P. ephippifer species complex, suggesting an important role of karyotypic divergence in its evolutionary history and underscoring the importance of hybridization as a mechanism of sex chromosome evolution in amphibians.

Genomic landscape of introgression from the ghost lineage in a gobiid fish uncovers the generality of forces shaping hybrid genomes

Extinct lineages can leave legacies in the genomes of extant lineages through ancient introgressive hybridization. The patterns of genomic survival of these extinct lineages provide insight into the role of extinct lineages in current biodiversity. However, our understanding on the genomic landscape of introgression from extinct lineages remains limited due to challenges associated with locating the traces of unsampled 'ghost' extinct lineages without ancient genomes. Herein, we conducted population genomic analyses on the East China Sea (ECS) lineage of Chaenogobius annularis, which was suspected to have originated from ghost introgression, with the aim of elucidating its genomic origins and characterizing its landscape of introgression. By combining phylogeographic analysis and demographic modelling, we demonstrated that the ECS lineage originated from ancient hybridization with an extinct ghost lineage. Forward simulations based on the estimated demography indicated that the statistic γ of the HyDe analysis can be used to distinguish the differences in local introgression rates in our data. Consistent with introgression between extant organisms, we found reduced introgression from extinct lineage in regions with low recombination rates and with functional importance, thereby suggesting a role of linked selection that has eliminated the extinct lineage in shaping the hybrid genome. Moreover, we identified enrichment of repetitive elements in regions associated with ghost introgression, which was hitherto little known but was also observed in the re-analysis of published data on introgression between extant organisms. Overall, our findings underscore the unexpected similarities in the characteristics of introgression landscapes across different taxa, even in cases of ghost introgression.

Assessing the phylogenetic relationship among varieties of Toona ciliata (Meliaceae) in sympatry with chloroplast genomes

Toona ciliata is an endangered species due to over-cutting and low natural regeneration in China. Its genetic conservation is of an increasing concern. However, several varieties are recognized according to the leaf and flower traits, which complicates genetic conservation of T. ciliata. Here, we sequenced the whole chloroplast genome sequences of three samples for each of four varieties (T. ciliata var. ciliata, T. ciliata var. yunnanensis, T. ciliata var. pubescens, and T. ciliata var. henryi) in sympatry and assessed their phylogenetic relationship at a fine spatial scale. The four varieties had genome sizes ranged from 159,546 to 159,617 bp and had small variations in genome structure. Phylogenomic analysis indicated that the four varieties were genetically well-mixed in branch groups. Genetic diversity from the whole chloroplast genome sequences of 12 samples was low among varieties (average π = 0.0003). Besides, we investigated genetic variation of 58 samples of the four varieties in sympatry using two markers (psaA and trnL-trnF) and showed that genetic differentiation was generally insignificant among varieties (Ф st = 0%-5%). Purifying selection occurred in all protein-coding genes except for the ycf2 gene that was under weak positive selection. Most amino acid sites in all protein-coding genes were under purifying selection except for a few sites that were under positive selection. The chloroplast genome-based phylogeny did not support the morphology-based classification. The overall results implicated that a conservation strategy based on the T. ciliata complex rather than on intraspecific taxon was more appropriate.

Signals of differential introgression in the genome of natural hybrids of Caribbean anoles

Hybridization facilitates recombination between divergent genetic lineages and can be shaped by both neutral and selective processes. Upon hybridization, loci with no net fitness effects introgress randomly from parental species into the genomes of hybrid individuals. Conversely, alleles from one parental species at some loci may provide a selective advantage to hybrids, resulting in patterns of introgression that do not conform to random expectations. We investigated genomic patterns of differential introgression in natural hybrids of two species of Caribbean anoles, Anolis pulchellus and A. krugi in Puerto Rico. Hybrids exhibit A. pulchellus phenotypes but possess A. krugi mitochondrial DNA, originated from multiple, independent hybridization events, and appear to have replaced pure A. pulchellus across a large area in western Puerto Rico. Combining genome-wide SNP datasets with bioinformatic methods to identify signals of differential introgression in hybrids, we demonstrate that the genomes of hybrids are dominated by pulchellus-derived alleles and show only 10%-20% A. krugi ancestry. The majority of A. krugi loci in hybrids exhibit a signal of non-random differential introgression and include loci linked to genes involved in development and immune function. Three of these genes (delta like canonical notch ligand 1, jagged1 and notch receptor 1) affect cell differentiation and growth and interact with mitochondrial function. Our results suggest that differential non-random introgression for a subset of loci may be driven by selection favouring the inheritance of compatible mitochondrial and nuclear-encoded genes in hybrids.

Biogeography and the evolution of acoustic communication in the polyploid North American grey treefrog complex

After polyploid species are formed, interactions between diploid and polyploid lineages may generate additional diversity in novel cytotypes and phenotypes. In anurans, mate choice by acoustic communication is the primary method by which individuals identify their own species and assess suitable mates. As such, the evolution of acoustic signals is an important mechanism for contributing to reproductive isolation and diversification in this group. Here, we estimate the biogeographical history of the North American grey treefrog complex, consisting of the diploid Hyla chrysoscelis and the tetraploid Hyla versicolor, focusing specifically on the geographical origin of whole genome duplication and the expansion of lineages out of glacial refugia. We then test for lineage-specific differences in mating signals by applying comparative methods to a large acoustic data set collected over 52 years that includes >1500 individual frogs. Along with describing the overall biogeographical history and call diversity, we found evidence that the geographical origin of H. versicolor and the formation of the midwestern polyploid lineage are both associated with glacial limits, and that the southwestern polyploid lineage is associated with a shift in acoustic phenotype relative to the diploid lineage with which they share a mitochondrial lineage. In H. chrysoscelis, we see that acoustic signals are largely split by Eastern and Western lineages, but that northward expansion along either side of the Appalachian Mountains is associated with further acoustic diversification. Overall, results of this study provide substantial clarity on the evolution of grey treefrogs as it relates to their biogeography and acoustic communication.

Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

Variation in Rice Plastid Genomes in Wide Crossing Reveals Dynamic Nucleo-Cytoplasmic Interaction

Plastid genomes (plastomes) of angiosperms are well known for their relative stability in size, structure, and gene content. However, little is known about their heredity and variations in wide crossing. To such an end, the plastomes of five representative rice backcross inbred lines (BILs) developed from crosses of O. glaberrima/O. sativa were analyzed. We found that the size of all plastomes was about 134,580 bp, with a quadripartite structure that included a pair of inverted repeat (IR) regions, a small single-copy (SSC) region and a large single-copy (LSC) region. They contained 76 protein genes, 4 rRNA genes, and 30 tRNA genes. Although their size, structure, and gene content were stable, repeat-mediated recombination, gene expression, and RNA editing were extensively changed between the maternal line and the BILs. These novel discoveries demonstrate that wide crossing causes not only nuclear genomic recombination, but also plastome variation in plants, and that the plastome plays a critical role in coordinating the nuclear-cytoplasmic interaction.

Evolutionary genetics of the mitochondrial genome: insights from Drosophila

Mitochondria are key to energy conversion in virtually all eukaryotes. Intriguingly, despite billions of years of evolution inside the eukaryote, mitochondria have retained their own small set of genes involved in the regulation of oxidative phosphorylation (OXPHOS) and protein translation. Although there was a long-standing assumption that the genetic variation found within the mitochondria would be selectively neutral, research over the past 3 decades has challenged this assumption. This research has provided novel insight into the genetic and evolutionary forces that shape mitochondrial evolution and broader implications for evolutionary ecological processes. Many of the seminal studies in this field, from the inception of the research field to current studies, have been conducted using Drosophila flies, thus establishing the species as a model system for studies in mitochondrial evolutionary biology. In this review, we comprehensively review these studies, from those focusing on genetic processes shaping evolution within the mitochondrial genome, to those examining the evolutionary implications of interactions between genes spanning mitochondrial and nuclear genomes, and to those investigating the dynamics of mitochondrial heteroplasmy. We synthesize the contribution of these studies to shaping our understanding of the evolutionary and ecological implications of mitochondrial genetic variation.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

A Test of Haldane's Rule in Neodiprion Sawflies and Implications for the Evolution of Postzygotic Isolation in Haplodiploids

AbstractHaldane's rule-a pattern in which hybrid sterility or inviability is observed in the heterogametic sex of an interspecific cross-is one of the most widely obeyed rules in nature. Because inheritance patterns are similar for sex chromosomes and haplodiploid genomes, Haldane's rule may apply to haplodiploid taxa, predicting that haploid male hybrids will evolve sterility or inviability before diploid female hybrids. However, there are several genetic and evolutionary mechanisms that may reduce the tendency of haplodiploids to obey Haldane's rule. Currently, there are insufficient data from haplodiploids to determine how frequently they adhere to Haldane's rule. To help fill this gap, we crossed a pair of haplodiploid hymenopteran species (Neodiprion lecontei and Neodiprion pinetum) and evaluated the viability and fertility of female and male hybrids. Despite considerable divergence, we found no evidence of reduced fertility in hybrids of either sex, consistent with the hypothesis that hybrid sterility evolves slowly in haplodiploids. For viability, we found a pattern opposite to that of Haldane's rule: hybrid females, but not males, had reduced viability. This reduction was most pronounced in one direction of the cross, possibly due to a cytoplasmic-nuclear incompatibility. We also found evidence of extrinsic postzygotic isolation in hybrids of both sexes, raising the possibility that this form or reproductive isolation tends to emerge early in speciation in host-specialized insects. Our work emphasizes the need for more studies on reproductive isolation in haplodiploids, which are abundant in nature but underrepresented in the speciation literature.

Mitonuclear Interactions and the Origin of Macaque Societies

In most eukaryotes, aerobic respiration requires interactions between autosomally encoded genes (Ninteract genes) and mitochondrial DNA, RNA, and protein. In species where females are philopatric, contrasting distributions of genetic variation in mitochondrial and nuclear genomes create variation in mitonuclear interactions that may be subject to natural selection. To test this expectation, we turned to a group with extreme female philopatry: the macaque monkeys. We examined four genomic data sets from (1) wild caught and (2) captive populations of rhesus macaque, which is the most widely distributed nonhuman primate, and (3) the stump-tailed macaque and (4) a subspecies of longtail macaque, both of whose mitochondrial DNA is introgressed from a highly diverged ancestor. We identified atypically long runs of homozygosity, low polymorphism, high differentiation, and/or rapid protein evolution associated with Ninteract genes compared with non-Ninteract genes. These metrics suggest a subset of Ninteract genes were independently subject to atypically pervasive natural selection in multiple species. These findings suggest that natural selection on mitonuclear interactions could have influenced several aspects of macaque societies including species diversity, ecological breadth, female-biased adult sex ratio and demography, sexual dimorphism, and mitonuclear phylogenomics.

Species delimitation using genomic data to resolve taxonomic uncertainties in a speciation continuum of pelagic seabirds

Speciation is a continuous and complex process shaped by the interaction of numerous evolutionary forces. Despite the continuous nature of the speciation process, the implementation of conservation policies relies on the delimitation of species and evolutionary significant units (ESUs). Puffinus shearwaters are globally distributed and threatened pelagic seabirds. Due to remarkable morphological status the group has been under intense taxonomic debate for the past three decades. Here, we use double digest Restriction-Site Associated DNA sequencing (ddRAD-Seq) to genotype species and subspecies of North Atlantic and Mediterranean Puffinus shearwaters across their entire geographical range. We assess the phylogenetic relationships and population structure among and within the group, evaluate species boundaries, and characterise the genomic landscape of divergence. We find that current taxonomies are not supported by genomic data and propose a more accurate taxonomy by integrating genomic information with other sources of evidence. Our results show that several taxon pairs are at different stages of a speciation continuum. Our study emphasises the potential of genomic data to resolve taxonomic uncertainties, which can help to focus management actions on relevant taxa, even if they do not necessarily coincide with the taxonomic rank of species.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

High heterogeneity in genomic differentiation between phenotypically divergent songbirds: a test of mitonuclear co-introgression

Comparisons of genomic variation among closely related species often show more differentiation in mitochondrial DNA (mtDNA) and sex chromosomes than in autosomes, a pattern expected due to the differing effective population sizes and evolutionary dynamics of these genomic components. Yet, introgression can cause species pairs to deviate dramatically from general differentiation trends. The yellowhammer (Emberiza citrinella) and pine bunting (E. leucocephalos) are hybridizing avian sister species that differ greatly in appearance and moderately in nuclear DNA, but that show no mtDNA differentiation. This discordance is best explained by adaptive mtDNA introgression-a process that can select for co-introgression at nuclear genes with mitochondrial functions (mitonuclear genes). To better understand these discordant differentiation patterns and characterize nuclear differentiation in this system, we investigated genome-wide differentiation between allopatric yellowhammers and pine buntings and compared it to what was seen previously in mtDNA. We found significant nuclear differentiation that was highly heterogeneous across the genome, with a particularly wide differentiation peak on the sex chromosome Z. We further investigated mitonuclear gene co-introgression between yellowhammers and pine buntings and found support for this process in the direction of pine buntings into yellowhammers. Genomic signals indicative of co-introgression were common in mitonuclear genes coding for subunits of the mitoribosome and electron transport chain complexes. Such introgression of mitochondrial DNA and mitonuclear genes provides a possible explanation for the patterns of high genomic heterogeneity in genomic differentiation seen among some species groups.

Unexpected genetic integrity boosts hope for the conservation of the red-legged partridge (Alectoris rufa, Galliformes) in Italy

The red-legged partridge (Alectoris rufa) is a medium-sized galliform endemic to southwestern Europe. In the easternmost part of the species' range, the population inhabiting Elba Island (Tuscan Archipelago National Park, Italy) is of undisputed conservation value. While we found nuclear and maternal DNA introgression with the exotic chukar partridge (A. chukar) in previous studies based on microsatellite DNA (n = 25) and two mitochondrial markers (n = 103), respectively, we disclosed a limited or null admixture in a few Elban partridges (n = 4) in a recent genomic investigation relying on 168,675 Single Nucleotide Polymorphisms (SNPs). We herein carried out an extended microsatellite DNA survey including additional 65 samples (total, 90) and six loci (total, 11) to determine both spatial structure and genetic integrity of local A. rufa. A sharp divergence between the subpopulations inhabiting the two sides of the island was disclosed, and the microsatellites indicated that all Elban partridges were not admixed with the chukar, thus fully reflecting the picture inferred using SNPs. We hypothesized that the spreading of chukar genes was constrained by negative selection, with the persistence of only the maternal lineage being indicative of thermal adaptation. The two subpopulations should be treated as distinct Management Units, and an envisaged plan to secure a stock onto nearby Pianosa Island could not only warrant endurance of the Elban population but also establish a source of valuable founders for the ex-situ management of the species in Italy. Our study exemplifies how a suitable samples/loci combination is the key to solve wildlife issues dealing with introgression.

Characterising Mitochondrial Capture in an Iberian Shrew

Mitochondrial introgression raises questions of biogeography and of the extent of reproductive isolation and natural selection. Previous phylogenetic work on the Sorex araneus complex revealed apparent mitonuclear discordance in Iberian shrews, indicating past hybridisation of Sorex granarius and the Carlit chromosomal race of S. araneus, enabling introgression of the S. araneus mitochondrial genome into S. granarius. To further study this, we genetically typed 61 Sorex araneus/coronatus/granarius from localities in Portugal, Spain, France, and Andorra at mitochondrial, autosomal, and sex-linked loci and combined our data with the previously published sequences. Our data are consistent with earlier data indicating that S. coronatus and S. granarius are the most closely related of the three species, confirming that S. granarius from the Central System mountain range in Spain captured the mitochondrial genome from a population of S. araneus. This mitochondrial capture event can be explained by invoking a biogeographical scenario whereby S. araneus was in contact with S. granarius during the Younger Dryas in central Iberia, despite the two species currently having disjunct distributions. We discuss whether selection favoured S. granarius with an introgressed mitochondrial genome. Our data also suggest recent hybridisation and introgression between S. coronatus and S. granarius, as well as between S. araneus and S. coronatus.

Genomic Signatures of Mitonuclear Coevolution in Mammals

Mitochondrial (mt) and nuclear-encoded proteins are integrated in aerobic respiration, requiring co-functionality among gene products from fundamentally different genomes. Different evolutionary rates, inheritance mechanisms, and selection pressures set the stage for incompatibilities between interacting products of the two genomes. The mitonuclear coevolution hypothesis posits that incompatibilities may be avoided if evolution in one genome selects for complementary changes in interacting genes encoded by the other genome. Nuclear compensation, in which deleterious mtDNA changes are offset by compensatory nuclear changes, is often invoked as the primary mechanism for mitonuclear coevolution. Yet, direct evidence supporting nuclear compensation is rare. Here, we used data from 58 mammalian species representing eight orders to show strong correlations between evolutionary rates of mt and nuclear-encoded mt-targeted (N-mt) proteins, but not between mt and non-mt-targeted nuclear proteins, providing strong support for mitonuclear coevolution across mammals. N-mt genes with direct mt interactions also showed the strongest correlations. Although most N-mt genes had elevated dN/dS ratios compared to mt genes (as predicted under nuclear compensation), N-mt sites in close contact with mt proteins were not overrepresented for signs of positive selection compared to noncontact N-mt sites (contrary to predictions of nuclear compensation). Furthermore, temporal patterns of N-mt and mt amino acid substitutions did not support predictions of nuclear compensation, even in positively selected, functionally important residues with direct mitonuclear contacts. Overall, our results strongly support mitonuclear coevolution across ∼170 million years of mammalian evolution but fail to support nuclear compensation as the major mode of mitonuclear coevolution.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

PREMISE

Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives.

METHODS

To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution.

RESULTS

Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species.

CONCLUSIONS

The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

PREMISE

Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives.

METHODS

To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution.

RESULTS

Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species.

CONCLUSIONS

The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

Highly active repeat-mediated recombination in the mitogenome of the holoparasitic plant Aeginetia indica

Mitogenomes of most flowering plants evolve slowly in sequence, but rapidly in structure. The rearrangements in structure are mainly caused by repeat-mediated recombination. However, patterns of repeat-mediated recombination vary substantially among plants, and to provide a comprehensive picture, characterization of repeat-mediated recombination should extend to more plant species, including parasitic plants with a distinct heterotrophic lifestyle. Here we assembled the mitogenome of the holoparasitic plant Aeginetia indica (Orobanchaceae) using Illumina sequencing reads. The mitogenome was assembled into a circular chromosome of 420,362 bp, 18,734 bp longer than that of another individual of A. indica which was assembled before as a linear molecule. Synteny analysis between the two mitogenomes revealed numerous rearrangements, unique regions of each individual and 0.2% sequence divergence in their syntenic regions. The A. indica mitogenome contains a gene content typical of flowering plants (33 protein-coding, 3 rRNA, and 17 tRNA genes). Repetitive sequences >30 bp in size totals 57,060 bp, representing 13.6% of the mitogenome. We examined recombination mediated by repeats >100 bp in size and found highly active recombination for all the repeats, including a very large repeat of ~16 kb. Recombination between these repeats can form much smaller subgenomic circular chromosomes, which may lead to rapid replication of mitochondrial DNA and thus be advantageous for A. indica with a parasitic lifestyle. In addition, unlike some other parasitic plants, A. indica shows no evidence for horizontal gene transfer of protein-coding genes in its mitogenome.

What is reproductive isolation?

Reproductive isolation (RI) is a core concept in evolutionary biology. It has been the central focus of speciation research since the modern synthesis and is the basis by which biological species are defined. Despite this, the term is used in seemingly different ways, and attempts to quantify RI have used very different approaches. After showing that the field lacks a clear definition of the term, we attempt to clarify key issues, including what RI is, how it can be quantified in principle, and how it can be measured in practice. Following other definitions with a genetic focus, we propose that RI is a quantitative measure of the effect that genetic differences between populations have on gene flow. Specifically, RI compares the flow of neutral alleles in the presence of these genetic differences to the flow without any such differences. RI is thus greater than zero when genetic differences between populations reduce the flow of neutral alleles between populations. We show how RI can be quantified in a range of scenarios. A key conclusion is that RI depends strongly on circumstances-including the spatial, temporal and genomic context-making it difficult to compare across systems. After reviewing methods for estimating RI from data, we conclude that it is difficult to measure in practice. We discuss our findings in light of the goals of speciation research and encourage the use of methods for estimating RI that integrate organismal and genetic approaches.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

Variation in cytonuclear expression accommodation among allopolyploid plants

Cytonuclear coevolution is a common feature among plants, which coordinates gene expression and protein products between the nucleus and organelles. Consequently, lineage-specific differences may result in incompatibilities between the nucleus and cytoplasm in hybrid taxa. Allopolyploidy is also a common phenomenon in plant evolution. The hybrid nature of allopolyploids may result in cytonuclear incompatibilities, but the massive nuclear redundancy created during polyploidy affords additional avenues for resolving cytonuclear conflict (i.e., cytonuclear accommodation). Here we evaluate expression changes in organelle-targeted nuclear genes for six allopolyploid lineages that represent four genera (i.e., Arabidopsis, Arachis, Chenopodium, and Gossypium) and encompass a range in polyploid ages. Because incompatibilities between the nucleus and cytoplasm could potentially result in biases toward the maternal homoeolog and/or maternal expression level, we evaluate patterns of homoeolog usage, expression bias, and expression-level dominance in cytonuclear genes relative to the background of non-cytonuclear expression changes and to the diploid parents. Although we find subsets of cytonuclear genes in most lineages that match our expectations of maternal preference, these observations are not consistent among either allopolyploids or categories of organelle-targeted genes. Our results indicate that cytonuclear expression evolution may be subtle and variable among genera and genes, likely reflecting a diversity of mechanisms to resolve nuclear-cytoplasmic incompatibilities in allopolyploid species.

Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations

Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

Increased resolution in the face of conflict: phylogenomics of the Neotropical bellflowers (Campanulaceae: Lobelioideae), a rapid plant radiation

BACKGROUND AND AIMS

The centropogonid clade (Lobelioideae: Campanulaceae) is an Andean-centred rapid radiation characterized by repeated convergent evolution of morphological traits, including fruit type and pollination syndromes. While previous studies have resolved relationships of lineages with fleshy fruits into subclades, relationships among capsular species remain unresolved. This lack of resolution has impeded reclassification of non-monophyletic genera, whose current taxonomy relies heavily on traits that have undergone convergent evolution.

METHODS

Targeted sequence capture using a probe-set recently developed for the centropogonid clade was used to obtain phylogenomic data from DNA extracted from both silica-dried and herbarium leaf tissue. These data were used to infer relationships among species using concatenated and partitioned species tree methods, and to quantify gene tree discordance.

KEY RESULTS

While silica-dried leaf tissue resulted in longer assembled sequence data, the inclusion of herbarium samples improved taxonomic representation. Relationships among baccate lineages are similar to those inferred in previous studies, although they differ for lineages within and among capsular clades. We improve the phylogenetic resolution of Siphocampylus, which forms ten groups of closely related species which we informally name. Two subclades of Siphocampylus and two individual species are rogue taxa whose placement differs widely across analyses. Gene tree discordance (including cytonuclear discordance) is rampant.

CONCLUSIONS

This first phylogenomic study of the centropogonid clade considerably improves our understanding of relationships in this rapid radiation. Differences across analyses and the possibility of additional lineage discoveries still hamper a solid and stable reclassification. Rapid morphological innovation corresponds with a high degree of phylogenomic complexity, including cytonuclear discordance, nuclear gene tree conflict and well-supported differences between analyses based on different nuclear loci. Together, these results point to a potential role of hemiplasy underlying repeated convergent evolution. This hallmark of rapid radiations is probably present in many other species-rich Andean plant radiations.

Global Patterns of subgenome evolution in organelle-targeted genes of six allotetraploid angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary “mismatches” between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Mitochondria as environments for the nuclear genome in Drosophila: mitonuclear G×G×E

Mitochondria evolved from a union of microbial cells belonging to distinct lineages that were likely anaerobic. The evolution of eukaryotes required a massive reorganization of the 2 genomes and eventual adaptation to aerobic environments. The nutrients and oxygen that sustain eukaryotic metabolism today are processed in mitochondria through coordinated expression of 37 mitochondrial genes and over 1000 nuclear genes. This puts mitochondria at the nexus of gene-by-gene (G×G) and gene-by-environment (G×E) interactions that sustain life. Here we use a Drosophila model of mitonuclear genetic interactions to explore the notion that mitochondria are environments for the nuclear genome, and vice versa. We construct factorial combinations of mtDNA and nuclear chromosomes to test for epistatic interactions (G×G), and expose these mitonuclear genotypes to altered dietary environments to examine G×E interactions. We use development time and genome-wide RNAseq analyses to assess the relative contributions of mtDNA, nuclear chromosomes, and environmental effects on these traits (mitonuclear G×G×E). We show that the nuclear transcriptional response to alternative mitochondrial "environments" (G×G) has significant overlap with the transcriptional response of mitonuclear genotypes to altered dietary environments. These analyses point to specific transcription factors (e.g., giant) that mediated these interactions, and identified coexpressed modules of genes that may account for the overlap in differentially expressed genes. Roughly 20% of the transcriptome includes G×G genes that are concordant with G×E genes, suggesting that mitonuclear interactions are part of an organism's environment.

Historical mitochondrial genome introgression confounds species delimitation: Evidence from phylogenetic inference in the odorrana grahami species complex

Species delimitation is essential to informing conservation policy and understanding ecological and evolutionary processes. Most of our recent gains in knowledge on animal diversity rely on morphological characteristics and mitochondrial (mt) DNA variation. Concordant results based on both have led to an unprecedented acceleration in the identification of new species and enriched the field of taxonomy. However, discordances are also found commonly between morphological and mtDNA evidence. This confounds species delimitation, especially when gene flow or mitochondrial genome introgression has occurred. Here we illustrate how mitochondrial genome introgression among species of the Odorrana grahami complex confounds species delimitation using the combined evidence of morphological characters, mitochondrial variation, and thousands of nuclear single nucleotide polymorphisms (SNPs) from genotyping-by-sequencing (GBS). Fifty-eight samples across the distribution of the O. grahami complex were included. The mtDNA matrilineal genealogy indicated two clades, with O. grahami and O. junlianensis clustered together. In contrast, all nuclear evidence including gene trees, species trees, and genetic structure analyses based on GBS data support three species with distinct genetic clusters. These three distinct genetic clusters also correspond to distinct morphological characters. They affirm the distinct taxonomic entities of both O. grahami and O. junlianensis, as well as a third clade distinct from either. Which species the third clade belongs to remains unclear and will require further testing. The nuclear genomic loci contradict the COI evidence, with indications of rampant historical mitochondrial genome introgression among the species of the O. grahami complex. These discordant signals previously confused species delimitation efforts in this group. Based on these findings, we recommend the integration of independent data, especially nuclear genomic evidence, in species delimitation so as to be robust against the pitfalls of mitochondrial introgression.

Rapid sequence evolution is associated with genetic incompatibilities in the plastid Clp complex

KEY MESSAGE

Replacing the native clpP1 gene in the Nicotiana plastid genome with homologs from different donor species showed that the extent of genetic incompatibilities depended on the rate of sequence evolution. The plastid caseinolytic protease (Clp) complex plays essential roles in maintaining protein homeostasis and comprises both plastid-encoded and nuclear-encoded subunits. Despite the Clp complex being retained across green plants with highly conserved protein sequences in most species, examples of extremely accelerated amino acid substitution rates have been identified in numerous angiosperms. The causes of these accelerations have been the subject of extensive speculation but still remain unclear. To distinguish among prevailing hypotheses and begin to understand the functional consequences of rapid sequence divergence in Clp subunits, we used plastome transformation to replace the native clpP1 gene in tobacco (Nicotiana tabacum) with counterparts from another angiosperm genus (Silene) that exhibits a wide range in rates of Clp protein sequence evolution. We found that antibiotic-mediated selection could drive a transgenic clpP1 replacement from a slowly evolving donor species (S. latifolia) to homoplasmy but that clpP1 copies from Silene species with accelerated evolutionary rates remained heteroplasmic, meaning that they could not functionally replace the essential tobacco clpP1 gene. These results suggest that observed cases of rapid Clp sequence evolution are a source of epistatic incompatibilities that must be ameliorated by coevolutionary responses between plastid and nuclear subunits.

The role of mitonuclear incompatibilities in allopatric speciation

Aerobic metabolism in eukaryotic cells requires extensive interactions between products of the nuclear and mitochondrial genomes. Rapid evolution of the mitochondrial genome, including fixation of both adaptive and deleterious mutations, creates intrinsic selection pressures favoring nuclear gene mutations that maintain mitochondrial function. As this process occurs independently in allopatry, the resulting divergence between conspecific populations can subsequently be manifest in mitonuclear incompatibilities in inter-population hybrids. Such incompatibilities, mitonuclear versions of Bateson-Dobzhansky-Muller incompatibilities that form the standard model for allopatric speciation, can potentially restrict gene flow between populations, ultimately resulting in varying degrees of reproductive isolation. The potential role of mitonuclear incompatibilities in speciation is further enhanced where mtDNA substitution rates are elevated compared to the nuclear genome and where population structure maintains allopatry for adequate time to evolve multiple mitonuclear incompatibilities. However, the fact that mitochondrial introgression occurs across species boundaries has raised questions regarding the efficacy of mitonuclear incompatibilities in reducing gene flow. Several scenarios now appear to satisfactorily explain this phenomenon, including cases where differences in mtDNA genetic load may drive introgression or where co-introgression of coadapted nuclear genes may support the function of introgressed mtDNA. Although asymmetries in reproductive isolation between taxa are consistent with mitonuclear incompatibilities, interactions between autosomes and sex chromosomes yield similar predictions that are difficult to disentangle. With regard to establishing reproductive isolation while in allopatry, existing studies clearly suggest that mitonuclear incompatibilities can contribute to the evolution of barriers to gene flow. However, there is to date relatively little definitive evidence supporting a primary role for mitonuclear incompatibilities in the speciation process.

Phylogenomics for Chagas Disease Vectors of the Rhodnius Genus (Hemiptera, Triatominae): What We Learn From Mito-Nuclear Conflicts and Recommendations

We provide in this study a very large DNA dataset on Rhodnius species including 36 samples representing 16 valid species of the three Rhodnius groups, pictipes, prolixus and pallescens. Samples were sequenced at low-depth with whole-genome shotgun sequencing (Illumina technology). Using phylogenomics including 15 mitochondrial genes (13.3 kb), partial nuclear rDNA (5.2 kb) and 51 nuclear protein-coding genes (36.3 kb), we resolve sticking points in the Rhodnius phylogeny. At the species level, we confirmed the species-specific status of R. montenegrensis and R. marabaensis and we agree with the synonymy of R. taquarussuensis with R. neglectus. We also invite to revisit the species-specific status of R. milesi that is more likely R. nasutus. We proposed to define a robustus species complex that comprises the four close relative species: R. marabaensis, R. montenegrensis, R. prolixus and R. robustus. As Psammolestes tertius was included in the Rhodnius clade, we strongly recommend reclassifying this species as R. tertius. At the Rhodnius group level, molecular data consistently supports the clustering of the pictipes and pallescens groups, more related to each other than they are to the prolixus group. Moreover, comparing mitochondrial and nuclear tree topologies, our results demonstrated that various introgression events occurred in all the three Rhodnius groups, in laboratory strains but also in wild specimens. We demonstrated that introgressions occurred frequently in the prolixus group, involving the related species of the robustus complex but also the pairwise R. nasutus and R. neglectus. A genome wide analysis highlighted an introgression event in the pictipes group between R. stali and R. brethesi and suggested a complex gene flow between the three species of the pallescens group, R. colombiensis, R. pallescens and R. ecuadoriensis. The molecular data supports also a sylvatic distribution of R. prolixus in Brazil (Pará state) and the monophyly of R. robustus. As we detected extensive introgression events and selective pressure on mitochondrial genes, we strongly recommend performing separate mitochondrial and nuclear phylogenies and to take advantages of mito-nuclear conflicts in order to have a comprehensive evolutionary vision of this genus.

Selection for biparental inheritance of mitochondria under hybridization and mitonuclear fitness interactions

Uniparental inheritance (UPI) of mitochondria predominates over biparental inheritance (BPI) in most eukaryotes. However, examples of BPI of mitochondria, or paternal leakage, are becoming increasingly prevalent. Most reported cases of BPI occur in hybrids of distantly related sub-populations. It is thought that BPI in these cases is maladaptive; caused by a failure of female or zygotic autophagy machinery to recognize divergent male-mitochondrial DNA ‘tags’. Yet recent theory has put forward examples in which BPI can evolve under adaptive selection, and empirical studies across numerous metazoan taxa have demonstrated outbreeding depression in hybrids attributable to disruption of population-specific mitochondrial and nuclear genotypes (mitonuclear mismatch). Based on these developments, we hypothesize that BPI may be favoured by selection in hybridizing populations when fitness is shaped by mitonuclear interactions. We test this idea using a deterministic, simulation-based population genetic model and demonstrate that BPI is favoured over strict UPI under moderate levels of gene flow typical of hybridizing populations. Our model suggests that BPI may be stable, rather than a transient phenomenon, in hybridizing populations.

Evolutionary innovations in Antarctic brittle stars linked to glacial refugia

The drivers behind evolutionary innovations such as contrasting life histories and morphological change are central questions of evolutionary biology. However, the environmental and ecological contexts linked to evolutionary innovations are generally unclear. During the Pleistocene glacial cycles, grounded ice sheets expanded across the Southern Ocean continental shelf. Limited ice-free areas remained, and fauna were isolated from other refugial populations. Survival in Southern Ocean refugia could present opportunities for ecological adaptation and evolutionary innovation. Here, we reconstructed the phylogeographic patterns of circum-Antarctic brittle stars Ophionotus victoriae and O. hexactis with contrasting life histories (broadcasting vs brooding) and morphology (5 vs 6 arms). We examined the evolutionary relationship between the two species using cytochrome c oxidase subunit I (COI) data. COI data suggested that O. victoriae is a single species (rather than a species complex) and is closely related to O. hexactis (a separate species). Since their recent divergence in the mid-Pleistocene, O. victoriae and O. hexactis likely persisted differently throughout glacial maxima, in deep-sea and Antarctic island refugia, respectively. Genetic connectivity, within and between the Antarctic continental shelf and islands, was also observed and could be linked to the Antarctic Circumpolar Current and local oceanographic regimes. Signatures of a probable seascape corridor linking connectivity between the Scotia Sea and Prydz Bay are also highlighted. We suggest that survival in Antarctic island refugia was associated with increase in arm number and a switch from broadcast spawning to brooding in O. hexactis, and propose that it could be linked to environmental changes (such as salinity) associated with intensified interglacial-glacial cycles.