Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle

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.

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.

Mitonuclear interactions modulate nutritional preference

In nature, organisms are faced with constant nutritional options which fuel key life-history traits. Studies have shown that species can actively make nutritional decisions based on internal and external cues. Metabolism itself is underpinned by complex genomic interactions involving components from both nuclear and mitochondrial genomes. Products from these two genomes must coordinate how nutrients are extracted, used and recycled. Given the complicated nature of metabolism, it is not well understood how nutritional choices are affected by mitonuclear interactions. This is under the rationale that changes in genomic interactions will affect metabolic flux and change physiological requirements. To this end we used a large Drosophila mitonuclear genetic panel, comprising nine isogenic nuclear genomes coupled to nine mitochondrial haplotypes, giving a total of 81 different mitonuclear genotypes. We use a capillary-based feeding assay to screen this panel for dietary preference between carbohydrate and protein. We find significant mitonuclear interactions modulating nutritional choices, with these epistatic interactions also being dependent on sex. Our findings support the notion that complex genomic interactions can place a constraint on metabolic flux. This work gives us deeper insights into how key metabolic interactions can have broad implications on behaviour.

Transcriptomic evidence indicates that montane leaf beetles prioritize digestion and reproduction in a sex-specific manner during emergence from dormancy

During winter, many organisms conserve resources by entering dormancy, suppressing metabolism and biosynthesis. The transition out of winter dormancy to summer activity requires a quick reversal of this suppression, in order to exploit now-favorable environmental conditions. To date, mechanisms by which winter climate variation affects this transition remains unelucidated. Here we experimentally manipulated snow cover for naturally overwintering montane leaf beetles (Chrysomela aeneicollis), and profiled changes in gene expression during the transition out of dormancy in spring. Upon emergence, beetles up-regulate transcripts associated with digestion and nutrient acquisition and down regulate those associated with lipid metabolism, suggesting a shift away from utilizing stored lipid and towards digestion of carbohydrate-rich host plant tissue. Development of digestive capacity is followed by up-regulation of transcripts associated with reproduction; a transition that occurs earlier in females than males. Snow manipulation strongly affected the ground thermal regime and correspondingly gene expression profiles, with beetles showing a delayed up-regulation of reproduction in the dry compared to snowy plots. This suggests that winter conditions can alter the timing and prioritization of processes during emergence from dormancy, potentially magnifying the effects of declining snow cover in the Sierra's and other snowy mountains.

Multi-locus genomic signatures of local adaptation to snow across the landscape in California populations of a willow leaf beetle

Organisms living in mountains contend with extreme climatic conditions, including short growing seasons and long winters with extensive snow cover. Anthropogenic climate change is driving unprecedented, rapid warming of montane regions across the globe, resulting in reduced winter snowpack. Loss of snow as a thermal buffer may have serious consequences for animals overwintering in soil, yet little is known about how variability in snowpack acts as a selective agent in montane ecosystems. Here, we examine genomic variation in California populations of the leaf beetle Chrysomela aeneicollis, an emerging natural model system for understanding how organisms respond to climate change. We used a genotype-environment association approach to identify genomic signatures of local adaptation to microclimate in populations from three montane regions with variable snowpack and a coastal region with no snow. We found that both winter-associated environmental variation and geographical distance contribute to overall genomic variation across the landscape. We identified non-synonymous variation in novel candidate loci associated with cytoskeletal function, ion transport and membrane stability, cellular processes associated with cold tolerance in other insects. These findings provide intriguing evidence that variation in snowpack imposes selective gradients in montane ecosystems.

A chromosome-scale genome assembly and evaluation of mtDNA variation in the willow leaf beetle Chrysomela aeneicollis

The leaf beetle Chrysomela aeneicollis has a broad geographic range across Western North America but is restricted to cool habitats at high elevations along the west coast. Central California populations occur only at high altitudes (2,700-3,500 m) where they are limited by reduced oxygen supply and recent drought conditions that are associated with climate change. Here, we report a chromosome-scale genome assembly alongside a complete mitochondrial genome and characterize differences among mitochondrial genomes along a latitudinal gradient over which beetles show substantial population structure and adaptation to fluctuating temperatures. Our scaffolded genome assembly consists of 21 linkage groups; one of which we identified as the X chromosome based on female/male whole genome sequencing coverage and orthology with Tribolium castaneum. We identified repetitive sequences in the genome and found them to be broadly distributed across all linkage groups. Using a reference transcriptome, we annotated a total of 12,586 protein-coding genes. We also describe differences in putative secondary structures of mitochondrial RNA molecules, which may generate functional differences important in adaptation to harsh abiotic conditions. We document substitutions at mitochondrial tRNA molecules and substitutions and insertions in the 16S rRNA region that could affect intermolecular interactions with products from the nuclear genome. This first chromosome-level reference genome will enable genomic research in this important model organism for understanding the biological impacts of climate change on montane insects.

Mapping mitonuclear epistasis using a novel recombinant yeast population

Genetic variation in mitochondrial and nuclear genomes can perturb mitonuclear interactions and lead to phenotypic differences between individuals and populations. Despite their importance to most complex traits, it has been difficult to identify the interacting mitonuclear loci. Here, we present a novel advanced intercrossed population of Saccharomyces cerevisiae yeasts, called the Mitonuclear Recombinant Collection (MNRC), designed explicitly for detecting mitonuclear loci contributing to complex traits. For validation, we focused on mapping genes that contribute to the spontaneous loss of mitochondrial DNA (mtDNA) that leads to the petite phenotype in yeast. We found that rates of petite formation in natural populations are variable and influenced by genetic variation in nuclear DNA, mtDNA and mitonuclear interactions. We mapped nuclear and mitonuclear alleles contributing to mtDNA stability using the MNRC by integrating a term for mitonuclear epistasis into a genome-wide association model. We found that the associated mitonuclear loci play roles in mitotic growth most likely responding to retrograde signals from mitochondria, while the associated nuclear loci with main effects are involved in genome replication. We observed a positive correlation between growth rates and petite frequencies, suggesting a fitness tradeoff between mitotic growth and mtDNA stability. We also found that mtDNA stability was correlated with a mobile mitochondrial GC-cluster that is present in certain populations of yeast and that selection for nuclear alleles that stabilize mtDNA may be rapidly occurring. The MNRC provides a powerful tool for identifying mitonuclear interacting loci that will help us to better understand genotype-phenotype relationships and coevolutionary trajectories.

G × G × E effect on phenotype expression in a non-conventional model organism, the unicellular slime mould Physarum polycephalum

In metazoans, the expression of key phenotypic traits is sensitive to two- and three-way interactions between variation in mitochondrial DNA, nuclear DNA and the external environment. Whether gene-by-environment interactions affect phenotypes in single-celled eukaryotes is poorly studied, except in a few species of yeast and fungi. We developed a genetic panel of the unicellular slime mould, Physarum polycephalum containing strains differing in mitochondrial and nuclear DNA haplotypes. The panel also included two strains harbouring a selfishly replicating mitochondrial-fusion (mF) plasmid that could affect phenotype expression. We assayed movement and growth rate differences among the strains across two temperature regimes: 24° and 28°C. We found that the slime mould's growth rate, but not movement, is affected by G × G × E interactions. Predictably, mtDNA × nDNA interactions significantly affected both traits. The inter-trait correlation across the strains in each temperature regime was positive. Surprisingly, the mF plasmid had no negative effects on our chosen traits. Our study is the first to demonstrate genetic regulation of phenotype expression in a unicellular slime mould. The genetic effect on phenotypes manifests via epistatic interactions with the thermal environment, thus shedding new light on the role of G × G × E interactions in trait evolution in protists.

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.

Snow modulates winter energy use and cold exposure across an elevation gradient in a montane ectotherm

Snow insulates the soil from air temperature, decreasing winter cold stress and altering energy use for organisms that overwinter in the soil. As climate change alters snowpack and air temperatures, it is critical to account for the role of snow in modulating vulnerability to winter climate change. Along elevational gradients in snowy mountains, snow cover increases but air temperature decreases, and it is unknown how these opposing gradients impact performance and fitness of organisms overwintering in the soil. We developed experimentally validated ecophysiological models of cold and energy stress over the past decade for the montane leaf beetle Chrysomela aeneicollis, along five replicated elevational transects in the Sierra Nevada mountains in California. Cold stress peaks at mid-elevations, while high elevations are buffered by persistent snow cover, even in dry years. While protective against cold, snow increases energy stress for overwintering beetles, particularly at low elevations, potentially leading to mortality or energetic tradeoffs. Declining snowpack will predominantly impact mid-elevation populations by increasing cold exposure, while high elevation habitats may provide refugia as drier winters become more common.

Positive selection on mitochondria may eliminate heritable microbes from arthropod populations

Diverse eukaryotic taxa carry facultative heritable symbionts, microbes that are passed from mother to offspring. These symbionts are coinherited with mitochondria, and selection favouring either new symbionts, or new symbiont variants, is known to drive loss of mitochondrial diversity as a correlated response. More recently, evidence has accumulated of episodic directional selection on mitochondria, but with currently unknown consequences for symbiont evolution. We therefore employed a population genetic mean field framework to model the impact of selection on mitochondrial DNA (mtDNA) upon symbiont frequency for three generic scenarios of host–symbiont interaction. Our models predict that direct selection on mtDNA can drive symbionts out of the population where a positively selected mtDNA mutation occurs initially in an individual that is uninfected with the symbiont, and the symbiont is initially at low frequency. When, by contrast, the positively selected mtDNA mutation occurs in a symbiont-infected individual, the mutation becomes fixed and in doing so removes symbiont variation from the population. We conclude that the molecular evolution of symbionts and mitochondria, which has previously been viewed from a perspective of selection on symbionts driving the evolution of a neutral mtDNA marker, should be reappraised in the light of positive selection on mtDNA.

Mitonuclear mismatch alters nuclear gene expression in naturally introgressed Rhinolophus bats

Background

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

Results

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

Conclusion

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

Supplementary Information

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

Extensive alternative splicing triggered by mitonuclear mismatch in naturally introgressed Rhinolophus bats

Mitochondrial function needs strong interactions of mitochondrial and nuclear (mitonuclear) genomes, which can be disrupted by mitonuclear mismatch due to mitochondrial DNA (mtDNA) introgression between two formerly isolated populations or taxa. This mitonuclear disruption may cause severe cellular stress in mismatched individuals. Gene expression changes and alternative splicing (AS) are two important transcriptional regulations to respond to environmental or cellular stresses. We previously identified a naturally introgressed population in the intermediate horseshoe bat (Rhinolophus affinis). Individuals from this population belong to R. a. himalayanus and share almost identical nuclear genetic background; however, some of them had mtDNA from another subspecies (R. a. macrurus). With this unique natural system, we examined gene expression changes in six tissues between five mitonuclear mismatched and five matched individuals. A small number of differentially expressed genes (DEGs) were identified, and functional enrichment analysis revealed that most DEGs were related to immune response although some may be involved in response to oxidative stress. In contrast, we identified extensive AS events and alternatively spliced genes (ASGs) between mismatched and matched individuals. Functional enrichment analysis revealed that multiple ASGs were directly or indirectly associated with energy production in mitochondria which is vital for survival of organism. To our knowledge, this is the first study to examine the role of AS in responding to cellular stress caused by mitonuclear mismatch in natural populations. Our results suggest that AS may play a more important role than gene expression regulation in responding to severe environmental or cellular stresses.

A wAlbB Wolbachia transinfection displays stable phenotypic effects across divergent Aedes aegypti mosquito backgrounds

Aedes mosquitoes harboring intracellular Wolbachia bacteria are being released in arbovirus and mosquito control programs. With releases taking place around the world, understanding the contribution of host variation to Wolbachia phenotype is crucial. We generated a Wolbachia transinfection (wAlbB^Q) in Aedes aegypti and performed backcrossing to introduce the infection into Australian or Malaysian nuclear backgrounds. Whole Wolbachia genome sequencing shows that the wAlbB^Q transinfection is near-identical to the reference wAlbB genome, suggesting few changes since the infection was first introduced to Ae. aegypti over 15 years ago. However, these sequences were distinct from other available wAlbB genome sequences, highlighting the potential diversity of wAlbB in natural Ae. albopictus populations. Phenotypic comparisons demonstrate effects of wAlbB infection on egg hatch and nuclear background on fecundity and body size, but no interactions between wAlbB infection and nuclear background for any trait. The wAlbB infection was stable at high temperatures and showed perfect maternal transmission and cytoplasmic incompatibility regardless of host background. Our results demonstrate the stability of wAlbB across host backgrounds and point to its long-term effectiveness for controlling arbovirus transmission and mosquito populations. Importance Wolbachia bacteria are being used to control the transmission of dengue and other arboviruses by mosquitoes. For Wolbachia release programs to be effective globally, Wolbachia infections must be stable across mosquito populations from different locations. In this study, we transferred Wolbachia (strain wAlbB) to Aedes aegypti mosquitoes with an Australian genotype and introduced the infection to Malaysian mosquitoes through backcrossing. We found that the phenotypic effects of Wolbachia are stable across both mosquito backgrounds. We sequenced the genome of wAlbB and found very few genetic changes despite spending over 15 years in a novel mosquito host. Our results suggest that the effects of Wolbachia infections are likely to remain stable across time and host genotype.

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

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

Maternal reproductive output and F1 hybrid fitness may influence contact zone dynamics

The outcome of secondary contact between divergent lineages or species may be influenced by both the reproductive traits of parental species and the fitness of offspring; however, their relative contributions have rarely been evaluated, particularly in longer-lived vertebrate species. We performed pure and reciprocal laboratory crosses between Ctenophorus decresii (tawny dragon) and C. modestus (swift dragon) to examine how parental reproductive traits and ecologically relevant offspring fitness traits may explain contact zone dynamics in the wild. The two species meet in a contact zone of post-F1 hybrids with asymmetric backcrossing and predominantly C. modestus mtDNA haplotypes. We found no evidence for reduced parental fecundity or offspring fitness for F1 hybrid crosses. However, maternal reproductive strategy differed between species, irrespective of the species of their mate. Ctenophorus modestus females had higher fecundity and produced more and larger clutches with lower embryonic mortality. Parental species also influenced sex ratios and offspring traits, with C. modestus ♀ × C. decresii ♂ hybrids exhibiting higher trait values for more fitness measures (growth rate, sprint speed, bite force) than offspring from all other pairings. Together, these patterns are consistent with the prevalence of C. modestus mtDNA in the contact zone, and asymmetric backcrossing likely reflects fitness effects that manifest in the F2 generation. Our results highlight how parental species can influence multiple offspring traits in different ways, which together may combine to influence offspring fitness and shape contact zone dynamics.

Mitochondria and the thermal limits of ectotherms

Temperature is a critical abiotic factor shaping the distribution and abundance of species, but the mechanisms that underpin organismal thermal limits remain poorly understood. One possible mechanism underlying these limits is the failure of mitochondrial processes, as mitochondria play a crucial role in animals as the primary site of ATP production. Conventional measures of mitochondrial performance suggest that these organelles can function at temperatures much higher than those that limit whole-organism function, suggesting that they are unlikely to set organismal thermal limits. However, this conclusion is challenged by recent data connecting sequence variation in mitochondrial genes to whole-organism thermal tolerance. Here, we review the current state of knowledge of mitochondrial responses to thermal extremes and ask whether they are consistent with a role for mitochondrial function in shaping whole-organism thermal limits. The available data are fragmentary, but it is possible to draw some conclusions. There is little evidence that failure of maximal mitochondrial oxidative capacity as assessed in vitro sets thermal limits, but there is some evidence to suggest that temperature effects on ATP synthetic capacity may be important. Several studies suggest that loss of mitochondrial coupling is associated with the thermal limits for organismal growth, although this needs to be rigorously tested. Most studies have utilized isolated mitochondrial preparations to assess the effects of temperature on these organelles, and there remain many untapped opportunities to address these questions using preparations that retain more of their biological context to better connect these subcellular processes with whole-organism thermal limits.

Digest: Mitonuclear interactions modulate life history phenotypes in the wild

Do mitonuclear interactions impact life history traits? Rank et al. found that these genomic interactions are of great importance in wild populations of the leaf beetle Chrysomela aeneicollis and may explain why populations are highly differentiated.