Evolutionary genetics of the mitochondrial genome: insights from Drosophila

The Mitochondrial Genome of Melipona fasciculata (Apidae, Meliponini): Genome Organization and Comparative Analyses, Phylogenetic Implications and Divergence Time Estimations

The native stingless bee Melipona fasciculata is economically and ecologically important to the Brazilian Northeast, providing a sustainable source of income to family farmers and being considered an effective pollinator in most ecosystems and crops. This study describes, for the first time, the mitogenome of the species and its phylogenetic position. The mitochondrial genome was sequenced using a MiSeq Sequencer (Illumina Inc.) and compared with other GenBank bee mitogenomes. The length of the mitochondrial DNA, excluding most of the control region, is 14,753 bp, and contains 13 protein-coding genes (PCGs), 21 transfer RNAs, 2 ribosomal RNAs (16S and 12S), and 1 AT-rich region. The GC-content of the M. fasciculata mitogenome was 13.4%. Of the 36 coding regions, 12 tRNAs and 9 PCGs were encoded in the heavy strand, and 9 tRNAs, 4 PCGs and 2 rRNAs were encoded in the light strand. The relative orientation and gene order was the same as other stingless bee mitogenomes. Phylogenetic inference produced well-resolved relationships with high statistical support for concordant branch topologies, under different optimization schemes and model parameters, within and among Melipona, Bombus, Apis, and related clades of Hymenoptera. In general, our divergence time estimates, which were based on the concatenated gene sequences (PCGs + rRNAs) from various groups, overlapped estimations captured by Bayesian analysis from different studies. The divergence time among Melipona species was estimated to occur during the Oligocene, approximately 24 Mya (95% HPD 14-36 Mya). Our results represent a valuable addition to help understanding not only the taxonomy and evolution of Brazilian stingless bee species, but also to uncover historical dispersal and isolation patterns in Meliponinae.

Convergence of autism proteins at the cilium

Hundreds of high confidence autism genes have been identified, yet the relevant etiological mechanisms remain unclear. Gene ontology analyses have repeatedly identified enrichment of proteins with annotated functions in gene expression regulation and neuronal communication. However, proteins are often pleiotropic and these annotations are inherently incomplete. Our recent autism functional genetics work has suggested that these genes may share a common mechanism at the cilium, a membrane-bound organelle critical for neurogenesis, brain patterning, and neuronal activity-all processes strongly implicated in autism. Moreover, autism commonly co-occurs with conditions that are known to involve ciliary-related pathologies, including congenital heart disease, hydrocephalus, and blindness. However, the role of autism genes at the cilium has not been systematically investigated. Here we demonstrate that autism proteins spanning disparate functional annotations converge in expression, localization, and function at cilia, and that patients with pathogenic variants in these genes have cilia-related co-occurring conditions and biomarkers of disrupted ciliary function. This degree of convergence among genes spanning diverse functional annotations strongly suggests that cilia are relevant to autism, as well as to commonly co-occurring conditions, and that this organelle should be explored further for therapeutic potential.

Rapamycin's lifespan effect is modulated by mito-nuclear epistasis in Drosophila

The macrolide drug rapamycin is a benchmark anti-ageing drug, which robustly extends lifespan of diverse organisms. For any health intervention, it is paramount to establish whether benefits are distributed equitably among individuals and populations, and ideally to match intervention to recipients' needs. However, how responses to rapamycin vary is surprisingly understudied. Here we investigate how among-population variation in both mitochondrial and nuclear genetics shapes rapamycin's effects on lifespan. We show that epistatic "mito-nuclear" interactions, between mitochondria and nuclei, modulate the response to rapamycin treatment. Differences manifest as differential demographic effects of rapamycin, with altered age-specific mortality rate. However, a fitness cost of rapamycin early in life does not show a correlated response, suggesting that mito-nuclear epistasis can decouple costs and benefits of treatment. These findings suggest that a deeper understanding of how variation in mitochondrial and nuclear genomes shapes physiology may facilitate tailoring of anti-ageing therapy to individual need.

Palmitoylated Importin α Regulates Mitotic Spindle Orientation Through Interaction with NuMA

Regulation of cell division orientation is a fundamental process critical to differentiation and tissue homeostasis. Microtubules emanating from the mitotic spindle pole bind a conserved complex of proteins at the cell cortex which orients the spindle and ultimately the cell division plane. Control of spindle orientation is of particular importance in developing tissues, such as the developing brain. Misorientation of the mitotic spindle and thus subsequent division plane misalignment can contribute to improper segregation of cell fate determinants in developing neuroblasts, leading to a rare neurological disorder known as microcephaly. We demonstrate that the nuclear transport protein importin α, when palmitoylated, plays a critical role in mitotic spindle orientation through localizing factors, such as NuMA, to the cell cortex. We also observe craniofacial developmental defects in Xenopus laevis when importin α palmitoylation is abrogated, including smaller head and brains, a hallmark of spindle misorientation and microcephaly. These findings characterize not only a role for importin α in spindle orientation, but also a broader role for importin α palmitoylation which has significance for many cellular processes.

Molecular Phylogenetics and Mitochondrial Genomic Evolution in the Endemic Genus Pielomastax (Orthoptera: Eumastacoidea) in China

BACKGROUND/OBJECTIVES

The genus Pielomastax Chang (Orthoptera: Eumastacoidea, 1937) is endemic to China, which is mainly distributed in low- and medium-altitude areas in central and eastern China. However, there are relatively few molecular data studies on the genus Pielomastax.

METHODS

In this study, three species of the genus Pielomastax were collected from Hubei and Henan, China, namely Pielomastax sp., Pielomastax shennongjiaensis Wang (1995) and Pielomastax tenuicerca Hsia and Liu (1989). Both Pielomastax sp. and Pielomastax shennongjiaensis were collected from the Shennongjia area of Hubei, but they exhibit some differences in morphological characteristics.

RESULTS

We obtained the mitochondrial genome structures of the three species, which were similar to those of the published mitochondrial genome structures of species within Eumastacoidea with 37 typical mitochondrial genes, including 13 PCGs, 22 tRNAs, and 2 ribosomal RNAs. The results of the maximum likelihood (ML) tree and the Bayesian inference (BI) tree showed that the families Eumastacidae, Chorotypidae and Episactinae in Eumastacoidea are a monophyletic group, and Thericleinae and Episactinae are sister clades. The time-calibrated phylogeny results indicated that the divergence time between Thericleinae and Episactinae was 95.58 Ma (56.71-128.02 Ma).

CONCLUSIONS

These phylogenetic tree results indicate that Pielomastax sp. and Pielomastax shennongjiaensis are the same species. And the time-calibrated phylogeny tree and the species distribution map of the genus Pielomastax indicate that the species of the genus Pielomastax spread from eastern to central China and diversified. These studies fill the gap in molecular data for the genus Pielomastax and the taxonomic status of Episactidae.

Targeted single cell expression profiling identifies integrators of sleep and metabolic state

Animals modulate sleep in accordance with their internal and external environments. Metabolic cues are particularly potent regulators of sleep, allowing animals to alter their sleep timing and amount depending on food availability and foraging duration. The fruit fly, Drosophila melanogaster, suppresses sleep in response to acute food deprivation, presumably to forage for food. This process is dependent on a single pair of Lateral Horn Leucokinin (LHLK) neurons, that secrete the neuropeptide Leucokinin. These neurons signal to insulin producing cells and suppress sleep under periods of starvation. The identification of individual neurons that modulate sleep-metabolism interactions provides the opportunity to examine the cellular changes associated with sleep modulation. Here, we use single-cell sequencing of LHLK neurons to examine the transcriptional responses to starvation. We validate that a Patch-seq approach selectively isolates RNA from individual LHLK neurons. Single-cell CEL-Seq comparisons of LHLK neurons between fed and 24-hr starved flies identified 24 genes that are differentially expressed in accordance with starvation state. In total, 12 upregulated genes and 12 downregulated genes were identified. Gene-ontology analysis showed an enrichment for Attacins, a family of anti-microbial peptides, along with several transcripts with diverse roles in regulating cellular function. Targeted knockdown of differentially expressed genes identified multiple genes that function within LHLK neurons to regulate sleep-metabolism interactions. Functionally validated genes include an essential role for the E3 ubiquitin Ligase insomniac, the sorbitol dehydrogenase Sodh1, as well as AttacinC and AttacinB in starvation-induced sleep suppression. Taken together, these findings provide a pipeline for identifying novel regulators of sleep-metabolism interactions within individual neurons.

Endemic Radiation of African Moonfish, Selene dorsalis (Gill 1863), in the Eastern Atlantic: Mitogenomic Characterization and Phylogenetic Implications of Carangids (Teleostei: Carangiformes)

This study offers an in-depth analysis of the mitochondrial genome of Selene dorsalis (Gill 1863), a species native to the Eastern Atlantic Ocean. The circular mitochondrial DNA molecule measures 16,541 base pairs and comprises 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and a control region (CR). The nucleotide composition exhibits a notable adenine-thymine (AT) bias, accounting for 53.13%, which aligns with other species in the Carangidae family. Most PCGs initiate with the ATG codon, with the exception of Cytochrome C oxidase subunit I, which starts with GTG. Analysis of relative synonymous codon usage reveals that leucine and serine are the most prevalent amino acids in the mitochondrial genome of S. dorsalis and its congeners (S. vomer and S. setapinnis). All tRNAs display the typical cloverleaf structure, though tRNA Serine (S1) lacks a dihydrouracil arm. Pairwise comparisons of synonymous and nonsynonymous substitutions for all PCGs yielded values below '1', indicating strong purifying selection. The CR spans 847 bp, representing 5.12% of the mitochondrial genome, and is characterized by high AT content (62.81%). It is situated between tRNA-Pro (TGG) and tRNA-Phe (GAA). The CR contains conserved sequence blocks, with CSB-1 being the longest at 22 bp and CSB-D the shortest at 18 bp. Phylogenetic analysis, using Bayesian and Maximum-likelihood trees constructed from concatenated PCGs across 72 species, successfully differentiates S. dorsalis from other carangids. This study also explores how ocean currents and gyres might influence lineage diversification and parapatric speciation of Selene species between the Atlantic and Pacific Oceans. These results highlight the importance of the mitochondrial genome in elucidating the structural organization and evolutionary dynamics of S. dorsalis and its relatives within marine ecosystems.

Mitochondrial Variation in Anopheles gambiae and Anopheles coluzzii: Phylogeographic Legacy and Mitonuclear Associations With Metabolic Resistance to Pathogens and Insecticides

Mitochondrial DNA has been a popular marker in phylogeography, phylogeny, and molecular ecology, but its complex evolution is increasingly recognized. Here, we investigated mitochondrial DNA variation in Anopheles gambiae and Anopheles coluzzii, in relation to other species in the Anopheles gambiae complex, by assembling the mitogenomes of 1,219 mosquitoes across Africa. The mitochondrial DNA phylogeny of the Anopheles gambiae complex was consistent with previously reported highly reticulated evolutionary history, revealing important discordances with the species tree. The three most widespread species (An. gambiae, An. coluzzii, and Anopheles arabiensis), known for extensive historical introgression, could not be discriminated based on mitogenomes. Furthermore, a monophyletic clustering of the three saltwater-tolerant species (Anopheles merus, Anopheles melas, and Anopheles bwambae) in the Anopheles gambiae complex also suggested that introgression and possibly selection shaped mitochondrial DNA evolution. Mitochondrial DNA variation in An. gambiae and An. coluzzii across Africa revealed significant partitioning among populations and species. A peculiar mitochondrial DNA lineage found predominantly in An. coluzzii and in the hybrid taxon of the African "far-west" exhibited divergence comparable to the interspecies divergence in the Anopheles gambiae complex, with a geographic distribution matching closely An. coluzzii's geographic range. This phylogeographic relict of the An. coluzzii and An. gambiae split was associated with population and species structure, but not with the rare Wolbachia occurrence. The lineage was significantly associated with single nucleotide polymorphisms in the nuclear genome, particularly in genes associated with pathogen and insecticide resistance. These findings underline potential mitonuclear coevolution history and the role played by mitochondria in shaping metabolic responses to pathogens and insecticides in Anopheles.

Hard Ticks as Vectors: The Emerging Threat of Tick-Borne Diseases in India

Hard ticks (Ixodidae) play a critical role in transmitting various tick-borne diseases (TBDs), posing significant global threats to human and animal health. Climatic factors influence the abundance, diversity, and vectorial capacity of tick vectors. It is imperative to have a comprehensive understanding of hard ticks, pathogens, eco-epidemiology, and the impact of climatic changes on the transmission dynamics of TBDs. The distribution and life cycle patterns of hard ticks are influenced by diverse ecological factors that, in turn, can be impacted by changes in climate, leading to the expansion of the tick vector's range and geographical distribution. Vector competence, a pivotal aspect of vectorial capacity, involves the tick's ability to acquire, maintain, and transmit pathogens. Hard ticks, by efficiently feeding on diverse hosts and manipulating their immunity through their saliva, emerge as competent vectors for various pathogens, such as viruses, parasites and bacteria. This ability significantly influences the success of pathogen transmission. Further exploration of genetic diversity, population structure, and hybrid tick vectors is crucial, as they play a substantial role in influencing vector competence and complicating the dynamics of TBDs. This comprehensive review deals with important TBDs in India and delves into a profound understanding of hard ticks as vectors, their biology, and the factors influencing their vector competence. Given that TBDs continue to pose a substantial threat to global health, the review emphasizes the urgency of investigating tick control strategies and advancing vaccine development. Special attention is given to the pivotal role of population genetics in comprehending the genetic diversity of tick populations and providing essential insights into their adaptability to environmental changes.

Cellular function of the GndA small open reading frame-encoded polypeptide during heat shock

Over the past 15 years, hundreds of previously undiscovered bacterial small open reading frame (sORF)-encoded polypeptides (SEPs) of fewer than fifty amino acids have been identified, and biological functions have been ascribed to an increasing number of SEPs from intergenic regions and small RNAs. However, despite numbering in the dozens in Escherichia coli, and hundreds to thousands in humans, same-strand nested sORFs that overlap protein coding genes in alternative reading frames remain understudied. In order to provide insight into this enigmatic class of unannotated genes, we characterized GndA, a 36-amino acid, heat shock-regulated SEP encoded within the +2 reading frame of the gnd gene in E. coli K-12 MG1655. We show that GndA pulls down components of respiratory complex I (RCI) and is required for proper localization of a RCI subunit during heat shock. At high temperature GndA deletion (ΔGndA) cells exhibit perturbations in cell growth, NADH+/NAD ratio, and expression of a number of genes including several associated with oxidative stress. These findings suggest that GndA may function in maintenance of homeostasis during heat shock. Characterization of GndA therefore supports the nascent but growing consensus that functional, overlapping genes occur in genomes from viruses to humans.

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BACKGROUND

Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.

RESULTS

Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.

CONCLUSIONS

Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

The role of mitochondria in sex- and age-specific gene expression in a species without sex chromosomes

Mitochondria perform an array of functions, many of which involve interactions with gene products encoded by the nucleus. These mitochondrial functions, particularly those involving energy production, can be expected to differ between sexes and across ages. Here, we measured mitochondrial effects on sex- and age-specific gene expression in parental and reciprocal F1 hybrids between allopatric populations of Tigriopus californicus with over 20% mitochondrial DNA divergence. Because the species lacks sex chromosomes, sex-biased mitochondrial effects are not confounded by the effects of sex chromosomes. Results revealed pervasive sex differences in mitochondrial effects, including effects on energetics and aging involving nuclear interactions throughout the genome. Using single-individual RNA sequencing, sex differences were found to explain more than 80% of the variance in gene expression. Males had higher expression of mitochondrial genes and mitochondrially targeted proteins (MTPs) involved in oxidative phosphorylation (OXPHOS), while females had elevated expression of non-OXPHOS MTPs, indicating strongly sex-dimorphic energy metabolism at the whole organism level. Comparison of reciprocal F1 hybrids allowed insights into the nature of mito-nuclear interactions, showing both mitochondrial effects on nuclear expression, and nuclear effects on mitochondrial expression. While based on a small set of crosses, sex-specific increases in mitochondrial expression with age were associated with longer life. Network analyses identified nuclear components of strong mito-nuclear interactions and found them to be sexually dimorphic. These results highlight the profound impact of mitochondria and mito-nuclear interactions on sex- and age-specific gene expression.

The first complete mitochondrional genome of Anopheles gibbinsi using a skimming sequencing approach

Mosquitoes belonging to the genus Anopheles are the only vectors of human malaria. Anopheles gibbinsi has been linked to malaria transmission in Kenya, with recent collections in Zambia reporting the mosquito species exhibiting zoophilic and exophilic behavioral patterns with occasional contact with humans. Given the paucity of genetic data, and challenges to identification and molecular taxonomy of the mosquitoes belonging to the Anopheles genus; we report the first complete mitochondrial genome of An. gibbinsi using a genome skimming approach. An Illumina Novaseq 6000 platform was used for sequencing, the length of the mitochondrial genome was 15401 bp, with 78.5% AT content comprised of 37 genes. Phylogenetic analysis by maximum likelihood using concatenation of the 13 protein coding genes demonstrated that An. marshallii was the closest relative based on existing sequence data. This study demonstrates that the skimming approach is an inexpensive and efficient approach for mosquito species identification and concurrent taxonomic rectification, which may be a useful alternative for generating reference sequence data for evolutionary studies among the Culicidae.

Complete mitochondrial genome of critically endangered catfish Hemibagrus punctatus (Jerdon, 1849) and comparative analysis for insights into the phylogeny of hemibagrids through mitogenomic approach

BACKGROUND

Hemibagrus punctatus (Jerdon, 1849) is a critically endangered bagrid catfish endemic to the Western Ghats of India, whose population is declining due to anthropogenic activities. The current study aims to compare the mitogenome of H. punctatus with that of other Bagrid catfishes and provide insights into their evolutionary relationships.

METHODS AND RESULTS

Samples were collected from Hemmige Karnataka, India. In the present study, the mitogenome of H. punctatus was successfully assembled, and its phylogenetic relationships with other Bagridae species were studied. The total genomic DNA of samples was extracted following the phenol-chloroform isoamyl alcohol method. Samples were sequenced, and the Illumina paired-end reads were assembled to a contig length of 16,517 bp. The mitochondrial genome was annotated using MitoFish and MitoAnnotator (Iwasaki et al., 2013). A robust phylogenetic analysis employing NJ (Maximum composite likelihood) and ASAP methods supports the classification of H. punctatus within the Bagridae family, which validates the taxonomic status of this species. In conclusion, this research enriches our understanding of H. punctatus mitogenome, shedding light on its evolutionary dynamics within the Bagridae family and contributing to the broader knowledge of mitochondrial genes in the context of evolutionary biology.

CONCLUSIONS

The study's findings contribute to a better understanding of the mitogenome of H. punctatus and provide insights into the evolutionary relationships within other Hemibagrids.

Mother's Curse effects on lifespan and aging

The Mother's Curse hypothesis posits that mothers curse their sons with harmful mitochondria, because maternal mitochondrial inheritance makes selection blind to mitochondrial mutations that harm only males. As a result, mitochondrial function may be evolutionarily optimized for females. This is an attractive explanation for ubiquitous sex differences in lifespan and aging, given the prevalence of maternal mitochondrial inheritance and the established relationship between mitochondria and aging. This review outlines patterns expected under the hypothesis, and traits most likely to be affected, chiefly those that are sexually dimorphic and energy intensive. A survey of the literature shows that evidence for Mother's Curse is limited to a few taxonomic groups, with the strongest support coming from experimental crosses in Drosophila. Much of the evidence comes from studies of fertility, which is expected to be particularly vulnerable to male-harming mitochondrial mutations, but studies of lifespan and aging also show evidence of Mother's Curse effects. Despite some very compelling studies supporting the hypothesis, the evidence is quite patchy overall, with contradictory results even found for the same traits in the same taxa. Reasons for this scarcity of evidence are discussed, including nuclear compensation, factors opposing male-specific mutation load, effects of interspecific hybridization, context dependency and demographic effects. Mother's Curse effects may indeed contribute to sex differences, but the complexity of other contributing factors make Mother's Curse a poor general predictor of sex-specific lifespan and aging.

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.

The role of mitochondria in sex- and age-specific gene expression in a species without sex chromosomes

Mitochondria perform an array of functions, many of which involve interactions with gene products encoded by the nucleus. These mitochondrial functions, particularly those involving energy production, can be expected to differ between sexes and across ages. Here we measured mitochondrial effects on sex- and age-specific gene expression in parental and reciprocal F1 hybrids between allopatric populations of Tigriopus californicus with over 20% mitochondrial DNA divergence. Because the species lacks sex chromosomes, sex-biased mitochondrial effects are not confounded by the effects of sex chromosomes. Using single-individual RNA sequencing, sex differences were found to explain more than 80% of the variance in gene expression. Males had higher expression of mitochondrial genes and mitochondrially targeted proteins (MTPs) involved in oxidative phosphorylation (OXPHOS), while females had elevated expression of non-OXPHOS MTPs, indicating strongly sex-dimorphic energy metabolism at the whole organism level. Comparison of reciprocal F1 hybrids allowed insights into the nature of mito-nuclear interactions, showing both mitochondrial effects on nuclear expression, as well as nuclear effects on mitochondrial expression. Across both sexes, increases in mitochondrial expression with age were associated with longer life. Network analyses identified nuclear components of strong mito-nuclear interactions, and found them to be sexually dimorphic. These results highlight the profound impact of mitochondria and mito-nuclear interactions on sex- and age-specific gene expression.

A Comparative Analysis of Selection Pressures Suffered by Mitochondrial Genomes in Two Planthopper Species with Divergent Climate Distributions

Mitochondrial DNA (mtDNA) has been widely used as a valuable tool in studies related to evolution and population genetics, under the implicit assumption of neutral evolution. However, recent studies suggest that natural selection also plays a significant role in shaping mitochondrial genome evolution, although the specific driving forces remain elusive. In this study, we aimed to investigate whether and how climate influences mitochondrial genome evolution by comparing the selection pressures acting on mitochondrial genomes between two rice planthoppers, Sogatella furcifera (Horváth) and Laodelphax striatellus (Fallén), which have different climate distributions. We employed the dN/dS method, MK test and Tajima's D tests for our analysis. Our results showed that the mitochondrial genomes of the two species appear to undergo predominantly purifying selection, consistent with the nearly neutral evolution model. However, we observed varied degrees of purifying selection among the 13 protein-coding genes. Notably, ND1, ND2, ND6, COIII, and ATP8 exhibited significantly stronger purifying selection and greater divergence between the two species compared to the other genes. Additionally, we observed relatively stronger purifying selection in the mitochondrial genomes of S. furcifera compared to L. striatellus. This difference could be attributed to varying metabolic requirements arising from distinct habitats or other factors that are unclear here. Furthermore, we speculate that mito-nuclear epistatic interactions may play a role in maintaining nonsynonymous polymorphisms, particularly for COI and COII. Overall, our results shed some light on the influence of climate on mitochondrial genome evolution.

Mitochondrial genome diversity across the subphylum Saccharomycotina

Introduction

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available.

Methods

By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades.

Results

Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements.

Discussion

As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.

Mitochondrial Genome Diversity across the Subphylum Saccharomycotina

Eukaryotic life depends on the functional elements encoded by both the nuclear genome and organellar genomes, such as those contained within the mitochondria. The content, size, and structure of the mitochondrial genome varies across organisms with potentially large implications for phenotypic variance and resulting evolutionary trajectories. Among yeasts in the subphylum Saccharomycotina, extensive differences have been observed in various species relative to the model yeast Saccharomyces cerevisiae, but mitochondrial genome sampling across many groups has been scarce, even as hundreds of nuclear genomes have become available. By extracting mitochondrial assemblies from existing short-read genome sequence datasets, we have greatly expanded both the number of available genomes and the coverage across sparsely sampled clades. Comparison of 353 yeast mitochondrial genomes revealed that, while size and GC content were fairly consistent across species, those in the genera Metschnikowia and Saccharomyces trended larger, while several species in the order Saccharomycetales, which includes S. cerevisiae, exhibited lower GC content. Extreme examples for both size and GC content were scattered throughout the subphylum. All mitochondrial genomes shared a core set of protein-coding genes for Complexes III, IV, and V, but they varied in the presence or absence of mitochondrially-encoded canonical Complex I genes. We traced the loss of Complex I genes to a major event in the ancestor of the orders Saccharomycetales and Saccharomycodales, but we also observed several independent losses in the orders Phaffomycetales, Pichiales, and Dipodascales. In contrast to prior hypotheses based on smaller-scale datasets, comparison of evolutionary rates in protein-coding genes showed no bias towards elevated rates among aerobically fermenting (Crabtree/Warburg-positive) yeasts. Mitochondrial introns were widely distributed, but they were highly enriched in some groups. The majority of mitochondrial introns were poorly conserved within groups, but several were shared within groups, between groups, and even across taxonomic orders, which is consistent with horizontal gene transfer, likely involving homing endonucleases acting as selfish elements. As the number of available fungal nuclear genomes continues to expand, the methods described here to retrieve mitochondrial genome sequences from these datasets will prove invaluable to ensuring that studies of fungal mitochondrial genomes keep pace with their nuclear counterparts.