Do relaxed selection and habitat temperature facilitate biased mitogenomic introgression in a narrowly endemic fish?

Hybridization in the Anthropocene - how pollution and climate change disrupt mate selection in freshwater fish

Chemical pollutants and/or climate change have the potential to break down reproductive barriers between species and facilitate hybridization. Hybrid zones may arise in response to environmental gradients and secondary contact between formerly allopatric populations, or due to the introduction of non-native species. In freshwater ecosystems, field observations indicate that changes in water quality and chemistry, due to pollution and climate change, are correlated with an increased frequency of hybridization. Physical and chemical disturbances of water quality can alter the sensory environment, thereby affecting chemical and visual communication among fish. Moreover, multiple chemical compounds (e.g. pharmaceuticals, metals, pesticides, and industrial contaminants) may impair fish physiology, potentially affecting phenotypic traits relevant for mate selection (e.g. pheromone production, courtship, and coloration). Although warming waters have led to documented range shifts, and chemical pollution is ubiquitous in freshwater ecosystems, few studies have tested hypotheses about how these stressors may facilitate hybridization and what this means for biodiversity and species conservation. Through a systematic literature review across disciplines (i.e. ecotoxicology and evolutionary biology), we evaluate the biological interactions, toxic mechanisms, and roles of physical and chemical environmental stressors (i.e. chemical pollution and climate change) in disrupting mate preferences and inducing interspecific hybridization in freshwater fish. Our study indicates that climate change-driven changes in water quality and chemical pollution may impact visual and chemical communication crucial for mate choice and thus could facilitate hybridization among fishes in freshwater ecosystems. To inform future studies and conservation management, we emphasize the importance of further research to identify the chemical and physical stressors affecting mate choice, understand the mechanisms behind these interactions, determine the concentrations at which they occur, and assess their impact on individuals, populations, species, and biological diversity in the Anthropocene.

Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

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.

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.

Interspecific Gene Flow and Mitochondrial Genome Capture During the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae)

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, A. opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting. Here we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene-tree, species-tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene-tree recovers the same relationships previously inferred for this group, which is strikingly different from the species-tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model, so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group, however coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent adaptive radiations.

Mitogenomics of the endemic Ethiopian rats: looking for footprints of adaptive evolution in sky islands

Organisms living in high altitude must adapt to environmental conditions with hypoxia and low temperature, e.g. by changes in the structure and function of proteins associated with oxidative phosphorylation in mitochondria. Here we analysed the signs of adaptive evolution in 27 mitogenomes of endemic Ethiopian rats (Stenocephalemys), where individual species adapted to different elevation. Significant signals of positive selection were detected in 10 of the 13 mitochondrial protein-coding genes, with a majority of functional substitutions in the NADH dehydrogenase complex. Higher frequency of positively selected sites was found in phylogenetic lineages corresponding to Afroalpine specialists.

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747-Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA-sequencing, RAD-sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA-sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO-lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage-specific selection on mitochondria-related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.

Population Genomic Evidence Reveals Subtle Patterns of Differentiation in the Trophically Polymorphic Cuatro Ciénegas Cichlid, Herichthys minckleyi

In recent decades, an increased understanding of molecular ecology has led to a reinterpretation of the role of gene flow during the evolution of reproductive isolation and biological novelty. For example, even in the face of ongoing gene flow strong selection may maintain divergent polymorphisms, or gene flow may introduce novel biological diversity via hybridization and introgression from a divergent species. Herein, we elucidate the evolutionary history and genomic basis of a trophically polymorphic trait in a species of cichlid fish, Herichthys minckleyi. We explored genetic variation at 3 hierarchical levels; between H. minckleyi (n = 69) and a closely related species Herichthys cyanoguttatus (n = 10), between H. minckleyi individuals from 2 geographic locations, and finally between individuals with alternate morphotypes at both a genome-wide and locus-specific scale. We found limited support for the hypothesis that the H. minckleyi polymorphism is the result of ongoing hybridization between the 2 species. Within H. minckleyi we found evidence of geographic genetic structure, and using traditional population genetic analyses found that individuals of alternate morphotypes within a pool appear to be panmictic. However, when we used a locus-specific approach to examine the relationship between multi-locus genotype, tooth size, and geographic sampling, we found the first evidence for molecular genetic differences between the H. minckleyi morphotypes.

Reconciling the Mitonuclear Compatibility Species Concept with Rampant Mitochondrial Introgression

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

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

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