Cytonuclear integration and co-evolution

Wheat Cybrid Plants, OryzaWheat, Regenerated from Wheat-Rice Hybrid Zygotes via in Vitro Fertilization System Possess Wheat-Rice Hybrid Mitochondria

Hybridization generates biodiversity, and wide hybridization plays a pivotal role in enhancing and broadening the useful attributes of crops. The hybridization barrier between wheat and rice, the two most important cereals, was recently overcome by in vitro production of allopolyploid wheat-rice hybrid zygotes, which can develop and grow into mature plants. In the study, genomic sequences and compositions of the possible hybrid plants were investigated through short- and long-read sequencing analyses and fluorescence in situ hybridization (FISH)-based visualization. The possible hybrid possessed whole wheat nuclear and cytoplasmic DNAs and rice mitochondrial (mt) DNA, along with variable retention rates of rice mtDNA ranging from 11% to 47%. The rice mtDNA retained in the wheat cybrid, termed Oryzawheat, can be transmitted across generations. In addition to mitochondrial hybridization, translocation of rice chromosome 1 into wheat chromosome 6A was detected in a F1 hybrid individual. OryzaWheat can provide a new horizon for utilizing inter-subfamily genetic resources among wheat and rice belonging to different subfamilies, Pooideae and Ehrhartoideae, respectively.

Comprehensive analysis of the Lycopodium japonicum mitogenome reveals abundant tRNA genes and cis-spliced introns in Lycopodiaceae species

Lycophytes and ferns represent one of the earliest-diverging lineages of vascular plants, with the Lycopodiaceae family constituting the basal clade among lycophytes. In this research, we successfully assembled and annotated the complete Lycopodium japonicum Thunb. (L. japonicum) mitochondrial genome (mitogenome) utilizing PacBio HiFi sequencing data, resulting in a single circular molecule with a size of 454,458 bp. 64 unique genes were annotated altogether, including 34 protein-coding genes, 27 tRNAs and 3 rRNAs. It also contains 32 group II introns, all of which undergo cis-splicing. We identified 195 simple sequence repeats, 1,948 dispersed repeats, and 92 tandem repeats in the L. japonicum mitogenome. Collinear analysis indicated that the mitogenomes of Lycopodiaceae are remarkably conserved compared to those of other vascular plants. We totally identified 326 RNA editing sites in 31 unique protein-coding genes with 299 sites converting cytosine to uracil and 27 sites the reverse. Notably, the L. japonicum mitogenome has small amounts foreign DNA from plastid or nuclear origin, accounting for only 2.81% of the mitogenome. The maximum likelihood phylogenetic analysis based on 23 diverse land plant mitogenomes and plastid genomes supports the basal position of lycophytes within vascular plants and they form a sister clade to all other vascular lineages, which is consistent with the PPG I classification system. As the first reported mitogenome of Lycopodioideae subfamily, this study enriches our understanding of Lycopodium mitogenomes, and sets the stage for future research on mitochondrial diversity and evolution within the lycophytes and ferns.

Assembly and comparative analysis of the complete mitochondrial genome of Fritillaria ussuriensis Maxim. (Liliales: Liliaceae), an endangered medicinal plant

BACKGROUND

Fritillaria ussuriensis is an endangered medicinal plant known for its notable therapeutic properties. Unfortunately, its population has drastically declined due to the destruction of forest habitats. Thus, effectively protecting F. ussuriensis from extinction poses a significant challenge. A profound understanding of its genetic foundation is crucial. To date, research on the complete mitochondrial genome of F. ussuriensis has not yet been reported.

RESULTS

The complete mitochondrial genome of F. ussuriensis was sequenced and assembled by integrating PacBio and Illumina sequencing technologies, revealing 13 circular chromosomes totaling 737,569 bp with an average GC content of 45.41%. A total of 55 genes were annotated in this mitogenome, including 2 rRNA genes, 12 tRNA genes, and 41 PCGs. The mitochondrial genome of F. ussuriensis contained 192 SSRs and 4,027 dispersed repeats. In the PCGs of F. ussuriensis mitogenome, 90.00% of the RSCU values exceeding 1 exhibited a preference for A-ended or U-ended codons. In addition, 505 RNA editing sites were predicted across these PCGs. Selective pressure analysis suggested negative selection on most PCGs to preserve mitochondrial functionality, as the notable exception of the gene nad3 showed positive selection. Comparison between the mitochondrial and chloroplast genomes of F. ussuriensis revealed 20 homologous fragments totaling 8,954 bp. Nucleotide diversity analysis revealed the variation among genes, and gene atp9 was the most notable. Despite the conservation of GC content, mitogenome sizes varied significantly among six closely related species, and colinear analysis confirmed the lack of conservation in their genomic structures. Phylogenetic analysis indicated a close relationship between F. ussuriensis and Lilium tsingtauense.

CONCLUSIONS

In this study, we sequenced and annotated the mitogenome of F. ussuriensis and compared it with the mitogenomes of other closely related species. In addition to genomic features and evolutionary position, this study also provides valuable genomic resources to further understand and utilize this medicinal plant.

Identification of a longevity gene through evolutionary rate covariation of insect mito-nuclear genomes

Oxidative phosphorylation, essential for energy metabolism and linked to the regulation of longevity, involves mitochondrial and nuclear genes. The functions of these genes and their evolutionary rate covariation (ERC) have been extensively studied, but little is known about whether other nuclear genes not targeted to mitochondria evolutionarily and functionally interact with mitochondrial genes. Here we systematically examined the ERC of mitochondrial and nuclear benchmarking universal single-copy ortholog (BUSCO) genes from 472 insects, identifying 75 non-mitochondria-targeted nuclear genes. We found that the uncharacterized gene CG11837-a putative ortholog of human DIMT1-regulates insect lifespan, as its knockdown reduces median lifespan in five diverse insect species and Caenorhabditis elegans, whereas its overexpression extends median lifespans in fruit flies and C. elegans and enhances oxidative phosphorylation gene activity. Additionally, DIMT1 overexpression protects human cells from cellular senescence. Together, these data provide insights into the ERC of mito-nuclear genes and suggest that CG11837 may regulate longevity across animals.

A Systematic Review and Developmental Perspective on Origin of CMS Genes in Crops

Cytoplasmic male sterility (CMS) arises from the incompatibility between the nucleus and cytoplasm as typical representatives of the chimeric structures in the mitochondrial genome (mitogenome), which has been extensively applied for hybrid seed production in various crops. The frequent occurrence of chimeric mitochondrial genes leading to CMS is consistent with the mitochondrial DNA (mtDNA) evolution. The sequence conservation resulting from faithfully maternal inheritance and the chimeric structure caused by frequent sequence recombination have been defined as two major features of the mitogenome. However, when and how these chimeric mitochondrial genes appear in the context of the highly conserved reproduction of mitochondria is an enigma. This review, therefore, presents the critical view of the research on CMS in plants to elucidate the mechanisms of this phenomenon. Generally, distant hybridization is the main mechanism to generate an original CMS source in natural populations and in breeding. Mitochondria and mitogenomes show pleomorphic and dynamic changes at key stages of the life cycle. The promitochondria in dry seeds develop into fully functioning mitochondria during seed imbibition, followed by massive mitochondria or mitogenome fusion and fission in the germination stage along with changes in the mtDNA structure and quantity. The mitogenome stability is controlled by nuclear loci, such as the nuclear gene Msh1. Its suppression leads to the rearrangement of mtDNA and the production of heritable CMS genes. An abundant recombination of mtDNA is also often found in distant hybrids and somatic/cybrid hybrids. Since mtDNA recombination is ubiquitous in distant hybridization, we put forward a hypothesis that the original CMS genes originated from mtDNA recombination during the germination of the hybrid seeds produced from distant hybridizations to solve the nucleo-cytoplasmic incompatibility resulting from the allogenic nuclear genome during seed germination.

De novo assembly and comprehensive analysis of the mitochondrial genome of Taxus wallichiana reveals different repeats mediate recombination to generate multiple conformations

Taxus plants are the exclusive source of paclitaxel, an anticancer drug with significant medicinal and economic value. Interspecies hybridization and gene introgression during evolution have obscured distinctions among Taxus species, complicating their phylogenetic classification. While the chloroplast genome of Taxus wallichiana, a widely distributed species in China, has been sequenced, its mitochondrial genome (mitogenome) remains uncharacterized.We sequenced and assembled the T. wallichiana mitogenome using BGI short reads and Nanopore long reads, facilitating comparisons with other gymnosperm mitogenomes. The T. wallichiana mitogenome spanning 469,949 bp, predominantly forms a circular configuration with a GC content of 50.51%, supplemented by 3 minor configurations mediated by one pair of LRs and two pairs of IntRs. It includes 32 protein-coding genes, 7 tRNA genes, and 3 rRNA genes, several of which exist in multiple copies.We detailed the mitogenome's structure, codon usage, RNA editing, and sequence migration between organelles, constructing a phylogenetic tree to elucidate evolutionary relationships. Unlike typical gymnosperm mitochondria, T. wallichiana shows no evidence of mitochondrial-plastid DNA transfer (MTPT), highlighting its unique genomic architecture. Synteny analysis indicated extensive genomic rearrangements in T. wallichiana, likely driven by recombination among abundant repetitive sequences. This study offers a high-quality T. wallichiana mitogenome, enhancing our understanding of gymnosperm mitochondrial evolution and supporting further cultivation and utilization of Taxus species.

Highly active repeat-mediated recombination in the mitogenome of the aquatic grass Hygroryza aristata

BACKGROUND

Floating bamboo (Hygroryza aristata) is an endangered species with a narrow native distribution and is renowned for its unique aesthetic qualities, which holds significant ecological and ornamental value. However, the lack of genetic information research, with only one complete plastome available, significantly hampers conservation efforts and further research for this species.

RESULTS

In this research, we sequenced and assembled the organelle genomes of floating bamboo, including the mitogenome (587,847 bp) and plastome (135,675 bp). The mitogenome can recombine into various configurations, which are mediated by 25 repeat pairs (13 SRs, 6 MRs, 1 LR, and 5 CRs). LR1 and SR5 are particularly notable as they have the ability to combine with other contigs, forming complex repeat units that facilitate further homologous recombination. The rate of homologous recombination varies significantly among species, yet there is still a pronounced positive correlation observed between the length of these repeat pairs and the rate of recombination they mediate. The mitogenome integrates seven intact protein-coding genes from the chloroplast. The codon usage patterns in both organelles are similar, with a noticeable bias towards C and T on the third codon. The gene map of Poales shows the entire loss of rpl6, succinate dehydrogenase subunits (sdh3 and sdh4). Additionally, the BOP clade retained more variable genes compared to the PACMAD clade.

CONCLUSIONS

We provided a high-quality and well-annotated mitogenome for floating bamboo and demonstrated the presence of diverse configurations. Our study has revealed the correlation between repeat length and their corresponding recombination rate despite variations among species. Although the mitogenome can potentially exist in the form of a unicircular in vivo, this occurrence is rare and may not be stable.

Graph-based mitochondrial genomes of three foundation species in the Saccharum genus

KEY MESSAGE

We reported the graph-based mitochondrial genomes of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) for the first time. The results revealed pan-structural variation and evolutionary processes in the mitochondrial genomes within Saccharum. Saccharum belongs to the Andropogoneae, and cultivars species in Saccharum contribute nearly 80% of sugar production in the world. To explore the genomic studies in Saccharum, we assembled 15 complete mitochondrial genomes (mitogenome) of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) using Illumina and Oxford Nanopore Technologies sequencing data. The mitogenomes of the three species were divided into a total of eight types based on contig numbers and linkages. All mitogenomes in the three species encoded 51 unique genes, including 32 protein-coding, 3 ribosomal RNA (rRNA) and 16 transfer RNA (tRNA) genes. The existence of long and short-repeat-mediated recombinations in the mitogenome of S. officinarum and S. robustum was revealed and confirmed through PCR validation. Furthermore, employing comparative genomics and phylogenetic analyses of the organelle genomes, we unveiled the evolutionary relationships and history of the major interspecific lineages in Saccharum genus. Phylogenetic analyses of homologous fragments between S. officinarum and S. robustum showed that S. officinarum and S. robustum are phylogenetically distinct and that they were likely parallel rather than domesticated. The variations between ancient (S. sinense and S. barberi) and modern cultivated species (S. hybrid) possibly resulted from hybridization involving different S. officinarum accessions. Lastly, this project reported the first graph-based mitogenomes of three Saccharum species, and a systematic comparison of the structural organization, evolutionary processes, and pan-structural variation of the Saccharum mitogenomes revealed the differential features of the Saccharum mitogenomes.

Plant organellar genomes: much done, much more to do

Plastids and mitochondria are the only organelles that possess genomes of endosymbiotic origin. In recent decades, advances in sequencing technologies have contributed to a meteoric rise in the number of published organellar genomes, and have revealed greatly divergent evolutionary trajectories. In this review, we quantify the abundance and distribution of sequenced plant organellar genomes across the plant tree of life. We compare numerous genomic features between the two organellar genomes, with an emphasis on evolutionary trajectories, transfers, the current state of organellar genome editing by transcriptional activator-like effector nucleases (TALENs), transcription activator-like effector (TALE)-mediated deaminase, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas), as well as genetic transformation. Finally, we propose future research to understand these different evolutionary trajectories, and genome-editing strategies to promote functional studies and eventually improve organellar genomes.

Genome-wide patterns of homoeologous gene flow in allotetraploid coffee

Premise

Allopolyploidy-a hybridization-induced whole-genome duplication event-has been a major driver of plant diversification. The extent to which chromosomes pair with their proper homolog vs. with their homoeolog in allopolyploids varies across taxa, and methods to detect homoeologous gene flow (HGF) are needed to understand how HGF has shaped polyploid lineages.

Methods

The ABBA-BABA test represents a classic method for detecting introgression between closely related species, but here we developed a modified use of the ABBA-BABA test to characterize the extent and direction of HGF in allotetraploid Coffea arabica.

Results

We found that HGF is abundant in the C. arabica genome, with both subgenomes serving as donors and recipients of variation. We also found that HGF is highly maternally biased in plastid-targeted-but not mitochondrial-targeted-genes, as would be expected if plastid-nuclear incompatibilities exist between the two parent species.

Discussion

Together, our analyses provide a simple framework for detecting HGF and new evidence consistent with selection favoring overwriting of paternally derived alleles by maternally derived alleles to ameliorate plastid-nuclear incompatibilities. Natural selection therefore appears to shape the direction and intensity of HGF in allopolyploid coffee, indicating that cytoplasmic inheritance has long-term consequences for polyploid lineages.

Ancestry affects the transcription of small mitochondrial RNAs in human lymphocytes

Mitochondrial mutations have been linked to changes in phenotypes such as fertility or longevity, however, these changes have been often inconsistent across populations for unknown reasons. A hypothesis that could explain this inconsistency is that some still uncharacterized mitochondrial products are mediating the phenotypic changes across populations. It has been hypothesized that one such product could be the small RNAs encoded in the mitochondrial genome, thus this work will provide new evidence for their existence and function. By using data from the 1000 genome project and knowledge from previously characterized nuclear small RNAs, this study found that 10 small RNAs encoded in tRNA fragments are consistently expressed in 450 individuals from five different populations. Furthermore, this study demonstrated that the expression of some small mitochondrial RNAs is different in individuals of African ancestry, similar to what was observed before in nuclear and mitochondria mRNAs. Lastly, we investigate the causes behind these differences in expression, showing that at least one of the mt-tRFs might be regulated by TRMT10B. The analyses presented in this work further support the small mitochondrial RNAs as functional molecules, and their population-specific expression supports the hypothesis that they act as a mediator between the nucleus and mitochondria differently across populations.

Comprehensive analysis of the complete mitochondrial genome of Lilium tsingtauense reveals a novel multichromosome structure

KEY MESSAGE

Lilium tsingtauense mitogenome comprises 27 independent chromosome molecules, it undergoes frequent genomic recombination, and the rate of recombination and mutation between different repetitive sequences affects the formation of multichromosomal structures. Given the extremely large genome of Lily, which likely harbors additional genetic resources, it serves as an ideal material for studying the phylogenetic evolution of organisms. Although the Lilium chloroplast genome has been documented, the sequence of its mitochondrial genome (mitogenome) remains uncharted. Using BGI short reads and Nanopore long reads, we sequenced, assembled, and annotated the mitogenome of Lilium tsingtauense. This effort culminated in the characterization of Lilium's first complete mitogenome. Comparative analysis with other angiosperms revealed the unique multichromosomal structure of the L. tsingtauense mitogenome, spanning 1,125,108 bp and comprising 27 independent circular chromosomes. It contains 36 protein-coding genes, 12 tRNA genes, and 3 rRNA genes, with a GC content of 44.90%. Notably, three chromosomes in the L. tsingtauense mitogenome lack identifiable genes, hinting at the potential existence of novel genes and noncoding elements. The high degree of observed genome fragmentation implies frequent reorganization, with recombination and mutation rates among diverse repetitive sequences likely driving the formation of multichromosomal structures. Our comprehensive analysis, covering genome size, coding genes, structure, RNA editing, repetitive sequences, and sequence migration, sheds light on the evolutionary and molecular biology of multichromosomal mitochondria in Lilium. This high-quality mitogenome of L. tsingtauense not only enriches our understanding of multichromosomal mitogenomes but also establishes a solid foundation for future genome breeding and germplasm innovation in Lilium.

Paternal leakage of plastids rescues inter-lineage hybrids in Silene nutans

BACKGROUND AND AIMS

Organelle genomes are usually maternally inherited in angiosperms. However, biparental inheritance has been observed, especially in hybrids resulting from crosses between divergent genetic lineages. When it concerns the plastid genome, this exceptional mode of inheritance might rescue inter-lineage hybrids suffering from plastid-nuclear incompatibilities. Genetically differentiated lineages of Silene nutans exhibit strong postzygotic isolation owing to plastid-nuclear incompatibilities, highlighted by inter-lineage hybrid chlorosis and mortality. Surviving hybrids can exhibit variegated leaves, which might indicate paternal leakage of the plastid genome. We tested whether the surviving hybrids inherited the paternal plastid genome and survived thanks to paternal leakage.

METHODS

We characterized the leaf phenotype (fully green, variegated or white) of 504 surviving inter-lineage hybrids obtained from a reciprocal cross experiment among populations of four genetic lineages (W1, W2, W3 and E1) of S. nutans from Western Europe and genotyped 560 leaf samples (both green and white leaves for variegated hybrids) using six lineage-specific plastid single nucleotide polymorphisms.

KEY RESULTS

A high proportion of the surviving hybrids (≤98 %) inherited the paternal plastid genome, indicating paternal leakage. The level of paternal leakage depended on cross type and cross direction. The E1 and W2 lineages as maternal lineages led to the highest hybrid mortality and to the highest paternal leakage from W1 and W3 lineages in the few surviving hybrids. This was consistent with E1 and W2 lineages, which contained the most divergent plastid genomes. When W3 was the mother, more hybrids survived, and no paternal leakage was detected.

CONCLUSIONS

By providing a plastid genome potentially more compatible with the hybrid nuclear background, paternal leakage has the potential to rescue inter-lineage hybrids from plastid-nuclear incompatibilities. This phenomenon might slow down the speciation process, provided hybrid survival and reproduction can occur in the wild.

Polyploid plants take cytonuclear perturbations in stride

Hybridization in plants is often accompanied by nuclear genome doubling (allopolyploidy), which has been hypothesized to perturb interactions between nuclear and organellar (mitochondrial and plastid) genomes by creating imbalances in the relative copy number of these genomes and producing genetic incompatibilities between maternally derived organellar genomes and the half of the allopolyploid nuclear genome from the paternal progenitor. Several evolutionary responses have been predicted to ameliorate these effects, including selection for changes in protein sequences that restore cytonuclear interactions; biased gene retention/expression/conversion favoring maternal nuclear gene copies; and fine-tuning of relative cytonuclear genome copy numbers and expression levels. Numerous recent studies, however, have found that evolutionary responses are inconsistent and rarely scale to genome-wide generalities. The apparent robustness of plant cytonuclear interactions to allopolyploidy may reflect features that are general to allopolyploids such as the lack of F2 hybrid breakdown under disomic inheritance, and others that are more plant-specific, including slow sequence divergence in organellar genomes and preexisting regulatory responses to changes in cell size and endopolyploidy during development. Thus, cytonuclear interactions may only rarely act as the main barrier to establishment of allopolyploid lineages, perhaps helping to explain why allopolyploidy is so pervasive in plant evolution.

The first two whole mitochondrial genomes for the genus Dactylis species: assembly and comparative genomics analysis

BACKGROUND

Orchardgrass (Dactylis glomerata L.), a perennial forage, has the advantages of rich leaves, high yield, and good quality and is one of the most significant forage for grassland animal husbandry and ecological management in southwest China. Mitochondrial (mt) genome is one of the major genetic systems in plants. Studying the mt genome of the genus Dactylis could provide more genetic information in addition to the nuclear genome project of the genus.

RESULTS

In this study, we sequenced and assembled two mitochondrial genomes of Dactylis species of D. glomerata (597, 281 bp) and D. aschersoniana (613, 769 bp), based on a combination of PacBio and Illumina. The gene content in the mitochondrial genome of D. aschersoniana is almost identical to the mitochondrial genome of D. glomerata, which contains 22-23 protein-coding genes (PCGs), 8 ribosomal RNAs (rRNAs) and 30 transfer RNAs (tRNAs), while D. glomerata lacks the gene encoding the Ribosomal protein (rps1) and D. aschersoniana contains one pseudo gene (atp8). Twenty-three introns were found among eight of the 30 protein-coding genes, and introns of three genes (nad 1, nad2, and nad5) were trans-spliced in Dactylis aschersoniana. Further, our mitochondrial genome characteristics investigation of the genus Dactylis included codon usage, sequences repeats, RNA editing and selective pressure. The results showed that a large number of short repetitive sequences existed in the mitochondrial genome of D. aschersoniana, the size variation of two mitochondrial genomes is due largely to the presence of a large number of short repetitive sequences. We also identified 52-53 large fragments that were transferred from the chloroplast genome to the mitochondrial genome, and found that the similarity was more than 70%. ML and BI methods used in phylogenetic analysis revealed that the evolutionary status of the genus Dactylis.

CONCLUSIONS

Thus, this study reveals the significant rearrangements in the mt genomes of Pooideae species. The sequenced Dactylis mt genome can provide more genetic information and improve our evolutionary understanding of the mt genomes of gramineous plants.

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.

Innovative computational tools provide new insights into the polyploid wheat genome

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

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.

Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Stronger evidence for relaxed selection than adaptive evolution in high-elevation animal mtDNA

Mitochondrial (mt) genes are the subject of many adaptive hypotheses due to the key role of mitochondria in energy production and metabolism. One widespread adaptive hypothesis is that selection imposed by life at high elevation leads to the rapid fixation of beneficial alleles in mtDNA, reflected in the increased rates of mtDNA evolution documented in many high-elevation species. However, the assumption that fast mtDNA evolution is caused by positive, rather than relaxed purifying selection has rarely been tested. Here, we calculated the dN/dS ratio, a metric of nonsynonymous substitution bias, and explicitly tested for relaxed selection in the mtDNA of over 700 species of terrestrial vertebrates, freshwater fishes, and arthropods, with information on elevation and latitudinal range limits, range sizes, and body sizes. We confirmed that mitochondrial genomes of high-elevation taxa have slightly higher dN/dS ratios compared to low-elevation relatives. High-elevation species tend to have smaller ranges, which predict higher dN/dS ratios and more relaxed selection across species and clades, while absolute elevation and latitude do not predict higher dN/dS. We also find a positive relationship between body mass and dN/dS, supporting a role for small effective population size leading to relaxed selection. We conclude that higher mt dN/dS among high-elevation species is more likely to reflect relaxed selection due to smaller ranges and reduced effective population size than adaptation to the environment. Our results highlight the importance of rigorously testing adaptive stories against non-adaptive alternative hypotheses, especially in mt genomes.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Complete mitochondrial genome of the endangered Prunus pedunculata (Prunoideae, Rosaceae) in China: characterization and phylogenetic analysis

Introduction

Prunus pedunculata (Prunoideae: Rosaceae), a relic shrub with strong resistance and multiple application values, is endangered in China. Extensive research had been devoted to gene expression, molecular markers, plastid genome analysis, and genetic background investigations of P. pedunculata. However, the mitochondrial genome of this species has not been systematically described, owing to the complexity of the plant mitogenome.

Methods

In the present research, the complete mitochondrial genome of P. pedunculata was assembled, annotated, and characterized. The genomic features, gene content and repetitive sequences were analyzed. The genomic variation and phylogenetic analysis have been extensively enumerated.

Results and discussion

The P. pedunculata mitogenome is a circular molecule with a total length of 405,855 bp and a GC content of 45.63%, which are the smallest size and highest GC content among the known Prunus mitochondrial genomes. The mitogenome of P. pedunculata encodes 62 genes, including 34 unique protein-coding genes (PCGs, excluding three possible pseudogenes), three ribosomal RNA genes, and 19 transfer RNA genes. The mitogenome is rich in repetitive sequences, counting 112 simple sequence repeats, 15 tandem repeats, and 50 interspersed repetitive sequences, with a total repeat length of 11,793 bp, accounting for 2.91% of the complete genome. Leucine (Leu) was a predominant amino acid in PCGs, with a frequency of 10.67%, whereas cysteine (Cys) and tryptophan (Trp) were the least adopted. The most frequently used codon was UUU (Phe), with a relative synonymous codon usage (RSCU) value of 1.12. Selective pressure was calculated based on 20 shared PCGs in the mitogenomes of the 32 species, most of which were subjected to purifying selection (Ka/Ks < 1), whereas ccmC and ccmFn underwent positive selection. A total of 262 potential RNA editing sites in 26 PCGs were identified. Furthermore, 56 chloroplast-derived fragments were ascertained in the mitogenome, ranging from 30 to 858 bp, and were mainly located across IGS (intergenic spacer) regions or rRNA genes. These findings verify the occurrence of intracellular gene transfer events from the chloroplast to the mitochondria. Furthermore, the phylogenetic relationship of P. pedunculata was supported by the mitogenome data of 30 other taxa of the Rosaceae family. Understanding the mitochondrial genome characteristics of P. pedunculata is of great importance to promote comprehension of its genetic background and this study provides a basis for the genetic breeding of Prunus.

Mitochondrial alternative oxidase across the tree of life: Presence, absence, and putative cases of lateral gene transfer

The alternative oxidase (AOX) is a terminal oxidase in the electron transport system that plays a role in mitochondrial bioenergetics. The past 20 years of research shows AOX has a wide yet patchy distribution across the tree of life. AOX has been suggested to have a role in stress tolerance, growth, and development in plants, but less is known about its function in other groups, including animals. In this study, we analyzed the taxonomic distribution of AOX across >2800 species representatives from prokaryotes and eukaryotes and developed a standardized workflow for finding and verifying the authenticity of AOX sequences. We found that AOX is limited to proteobacteria among prokaryotes, but is widely distributed in eukaryotes, with the highest prevalence in plants, fungi, and protists. AOX is present in many invertebrates, but is absent in others including most arthropods, and is absent from vertebrates. We found aberrant AOX sequences associated with some animal groups. Some of these aberrant AOXs were contaminants, but we also found putative cases of lateral gene transfer of AOX from fungi and protists to nematodes, springtails, fungus gnats, and rotifers. Our findings provide a robust and detailed analysis of the distribution of AOX and a method for identifying and verifying putative AOX sequences, which will be useful as more sequence data becomes available on public repositories.

Diversity and Evolution of Mitochondrial Translation Apparatus

The evolution of mitochondria has proceeded independently in different eukaryotic lines, which is reflected in the diversity of mitochondrial genomes and mechanisms of their expression in eukaryotic species. Mitochondria have lost most of bacterial ancestor genes by transferring them to the nucleus or eliminating them. However, mitochondria of almost all eukaryotic cells still retain relatively small genomes, as well as their replication, transcription, and translation apparatuses. The dependence on the nuclear genome, specific features of mitochondrial transcripts, and synthesis of highly hydrophobic membrane proteins in the mitochondria have led to significant changes in the translation apparatus inherited from the bacterial ancestor, which retained the basic structure necessary for protein synthesis but became more specialized and labile. In this review, we discuss specific properties of translation initiation in the mitochondria and how the evolution of mitochondria affected the functions of main factors initiating protein biosynthesis in these organelles.

Assembly and comparative analysis of the complete mitochondrial genome of Isopyrum anemonoides (Ranunculaceae)

Ranunculaceae is a large family of angiosperms comprising 2500 known species-a few with medicinal and ornamental values. Despite this, only two mitochondrial genomes (mitogenomes) of the family have been released in GenBank. Isopyrum anemonoides is a medicinal plant belonging to the family Ranunculaceae, and its chloroplast genome has recently been reported; however, its mitogenome remains unexplored. In this study, we assembled and analyzed the complete mitochondrial genome of I. anemonoides and performed a comparative analysis against different Ranunculaceae species, reconstructing the phylogenetic framework of Isopyrum. The circular mitogenome of I. anemonoides has a length of 206,722 bp, with a nucleotide composition of A (26.4%), T (26.4%), C (23.6%), and G (23.6%), and contains 62 genes, comprising 37 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and three ribosomal RNA (rRNA) genes. Abundantly interspersed repetitive and simple sequence repeat (SSR) loci were detected in the I. anemonoides mitogenome, with tetranucleotide repeats accounting for the highest proportion of SSRs. By detecting gene migration, we observed gene exchange between the chloroplast and mitogenome in I. anemonoides, including six intact tRNA genes, six PCG fragments, and fragments from two rRNA genes. Comparative mitogenome analysis of three Ranunculaceae species indicated that the PCG contents were conserved and the GC contents were similar. Selective pressure analysis revealed that only two genes (nad1 and rpl5) were under positive selection during their evolution in Ranunculales, and two specific RNA editing sites (atp6 and mttB) were detected in the I. anemonoides mitogenome. Moreover, a phylogenetic analysis based on the mitogenomes of I. anemonoides and the other 15 taxa accurately reflected the evolutionary and taxonomic status of I. anemonoides. Overall, this study provides new insights into the genetics, systematics, and evolution of mitochondrial evolution in Ranunculaceae, particularly I. anemonoides.

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

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.

Mutation pressure, drift, and the pace of molecular coevolution

Most aspects of the molecular biology of cells involve tightly coordinated intermolecular interactions requiring specific recognition at the nucleotide and/or amino acid levels. This has led to long-standing interest in the degree to which constraints on interacting molecules result in conserved vs. accelerated rates of sequence evolution, with arguments commonly being made that molecular coevolution can proceed at rates exceeding the neutral expectation. Here, a fairly general model is introduced to evaluate the degree to which the rate of evolution at functionally interacting sites is influenced by effective population sizes (Ne), mutation rates, strength of selection, and the magnitude of recombination between sites. This theory is of particular relevance to matters associated with interactions between organelle- and nuclear-encoded proteins, as the two genomic environments often exhibit dramatic differences in the power of mutation and drift. Although genes within low Ne environments can drive the rate of evolution of partner genes experiencing higher Ne, rates exceeding the neutral expectation require that the former also have an elevated mutation rate. Testable predictions, some counterintuitive, are presented on how patterns of coevolutionary rates should depend on the relative intensities of drift, selection, and mutation.

What is the potential impact of genetic divergence of plastid ribosomal genes between Silene nutans lineages in hybrids? An in silico approach using the 3D structure of the plastid ribosome

Introduction

Following the integration of cyanobacteria into the eukaryotic cells, many genes were transferred from the plastid to the nucleus. As a result, plastid complexes are encoded both by plastid and nuclear genes. Tight co-adaptation is required between these genes as plastid and nuclear genomes differ in several characteristics, such as mutation rate and inheritance patterns. Among these are complexes from the plastid ribosome, composed of two main subunits: a large and a small one, both composed of nuclear and plastid gene products. This complex has been identified as a potential candidate for sheltering plastid-nuclear incompatibilities in a Caryophyllaceae species, Silene nutans. This species is composed of four genetically differentiated lineages, which exhibit hybrid breakdown when interlineage crosses are conducted. As this complex is composed of numerous interacting plastid-nuclear gene pairs, in the present study, the goal was to reduce the number of gene pairs that could induce such incompatibilities.

Method

We used the previously published 3D structure of the spinach ribosome to further elucidate which of the potential gene pairs might disrupt plastid-nuclear interactions within this complex. After modeling the impact of the identified mutations on the 3D structure, we further focused on one strongly mutated plastid-nuclear gene pair: rps11-rps21. We used the centrality measure of the mutated residues to further understand if the modified interactions and associated modified centralities might be correlated with hybrid breakdown.

Results and discussion

This study highlights that lineage-specific mutations in essential plastid and nuclear genes might disrupt plastid-nuclear protein interactions of the plastid ribosome and that reproductive isolation correlates with changes in residue centrality values. Because of this, the plastid ribosome might be involved in hybrid breakdown in this system.

CRISPR-Cas9-based repeat depletion for high-throughput genotyping of complex plant genomes

High-throughput genotyping enables the large-scale analysis of genetic diversity in population genomics and genome-wide association studies that combine the genotypic and phenotypic characterization of large collections of accessions. Sequencing-based approaches for genotyping are progressively replacing traditional genotyping methods because of the lower ascertainment bias. However, genome-wide genotyping based on sequencing becomes expensive in species with large genomes and a high proportion of repetitive DNA. Here we describe the use of CRISPR-Cas9 technology to deplete repetitive elements in the 3.76-Gb genome of lentil (Lens culinaris), 84% consisting of repeats, thus concentrating the sequencing data on coding and regulatory regions (single-copy regions). We designed a custom set of 566,766 gRNAs targeting 2.9 Gbp of repeats and excluding repetitive regions overlapping annotated genes and putative regulatory elements based on ATAC-seq data. The novel depletion method removed ∼40% of reads mapping to repeats, increasing those mapping to single-copy regions by ∼2.6-fold. When analyzing 25 million fragments, this repeat-to-single-copy shift in the sequencing data increased the number of genotyped bases of ∼10-fold compared to nondepleted libraries. In the same condition, we were also able to identify ∼12-fold more genetic variants in the single-copy regions and increased the genotyping accuracy by rescuing thousands of heterozygous variants that otherwise would be missed because of low coverage. The method performed similarly regardless of the multiplexing level, type of library or genotypes, including different cultivars and a closely related species (L. orientalis). Our results showed that CRISPR-Cas9-driven repeat depletion focuses sequencing data on single-copy regions, thus improving high-density and genome-wide genotyping in large and repetitive genomes.

Identification of hidden associations among eukaryotic genes through statistical analysis of coevolutionary transitions

Coevolution at the gene level, as reflected by correlated events of gene loss or gain, can be revealed by phylogenetic profile analysis. The optimal method and metric for comparing phylogenetic profiles, especially in eukaryotic genomes, are not yet established. Here, we describe a procedure suitable for large-scale analysis, which can reveal coevolution based on the assessment of the statistical significance of correlated presence/absence transitions between gene pairs. This metric can identify coevolution in profiles with low overall similarities and is not affected by similarities lacking coevolutionary information. We applied the procedure to a large collection of 60,912 orthologous gene groups (orthogroups) in 1,264 eukaryotic genomes extracted from OrthoDB. We found significant cotransition scores for 7,825 orthogroups associated in 2,401 coevolving modules linking known and unknown genes in protein complexes and biological pathways. To demonstrate the ability of the method to predict hidden gene associations, we validated through experiments the involvement of vertebrate malate synthase-like genes in the conversion of (S)-ureidoglycolate into glyoxylate and urea, the last step of purine catabolism. This identification explains the presence of glyoxylate cycle genes in metazoa and suggests an anaplerotic role of purine degradation in early eukaryotes.

Sex differences in adult lifespan and aging rate across mammals: a test of the 'Mother Curse hypothesis'

In many animal species, including humans, males have shorter lifespan and show faster survival aging than females. This differential increase in mortality between sexes could result from the accumulation of deleterious mutations in the mitochondrial genome of males due to the maternal mode of mtDNA inheritance. To date, empirical evidence supporting the existence of this mechanism - called the Mother Curse hypothesis - remains largely limited to a few study cases in humans and Drosophila. In this study, we tested whether the Mother Curse hypothesis accounts for sex differences in lifespan and aging rate across 128 populations of mammals (60 and 68 populations studied in wild and captive conditions, respectively) encompassing 104 species. We found that adult lifespan decreases with increasing mtDNA neutral substitution rate in both sexes in a similar way in the wild - but not in captivity. Moreover, the aging rate marginally increased with neutral substitution rate in males and females in the wild. Overall, these results indicate that the Mother Curse hypothesis is not supported across mammals. We further discuss the implication of these findings for our understanding of the evolution of sex differences in mortality and aging.

Cellular and environmental dynamics influence species-specific extents of organelle gene retention

Mitochondria and plastids rely on many nuclear-encoded genes, but retain small subsets of the genes they need to function in their own organelle DNA (oDNA). Different species retain different numbers of oDNA genes, and the reasons for these differences are not completely understood. Here, we use a mathematical model to explore the hypothesis that the energetic demands imposed by an organism's changing environment influence how many oDNA genes it retains. The model couples the physical biology of cell processes of gene expression and transport to a supply-and-demand model for the environmental dynamics to which an organism is exposed. The trade-off between fulfilling metabolic and bioenergetic environmental demands, and retaining genetic integrity, is quantified for a generic gene encoded either in oDNA or in nuclear DNA. Species in environments with high-amplitude, intermediate-frequency oscillations are predicted to retain the most organelle genes, whereas those in less dynamic or noisy environments the fewest. We discuss support for, and insight from, these predictions with oDNA data across eukaryotic taxa, including high oDNA gene counts in sessile organisms exposed to day-night and intertidal oscillations (including plants and algae) and low counts in parasites and fungi.

Variation in mitogenome structural conformation in wild and cultivated lineages of sorghum corresponds with domestication history and plastome evolution

BACKGROUND

Mitochondria are organelles within eukaryotic cells that are central to the metabolic processes of cellular respiration and ATP production. However, the evolution of mitochondrial genomes (mitogenomes) in plants is virtually unknown compared to animal mitogenomes or plant plastids, due to complex structural variation and long stretches of repetitive DNA making accurate genome assembly more challenging. Comparing the structural and sequence differences of organellar genomes within and between sorghum species is an essential step in understanding evolutionary processes such as organellar sequence transfer to the nuclear genome as well as improving agronomic traits in sorghum related to cellular metabolism.

RESULTS

Here, we assembled seven sorghum mitochondrial and plastid genomes and resolved reticulated mitogenome structures with multilinked relationships that could be grouped into three structural conformations that differ in the content of repeats and genes by contig. The grouping of these mitogenome structural types reflects the two domestication events for sorghum in east and west Africa.

CONCLUSIONS

We report seven mitogenomes of sorghum from different cultivars and wild sources. The assembly method used here will be helpful in resolving complex genomic structures in other plant species. Our findings give new insights into the structure of sorghum mitogenomes that provides an important foundation for future research into the improvement of sorghum traits related to cellular respiration, cytonuclear incompatibly, and disease resistance.

Single-molecule Sequencing of an Animal Mitochondrial Genome Reveals Chloroplast-like Architecture and Repeat-mediated Recombination

Recent advances in long-read sequencing technology have allowed for single-molecule sequencing of entire mitochondrial genomes, opening the door for direct investigation of the mitochondrial genome architecture and recombination. We used PacBio sequencing to reassemble mitochondrial genomes from two species of New Zealand freshwater snails, Potamopyrgus antipodarum and Potamopyrgus estuarinus. These assemblies revealed a ∼1.7 kb structure within the mitochondrial genomes of both species that was previously undetected by an assembly of short reads and likely corresponding to a large noncoding region commonly present in the mitochondrial genomes. The overall architecture of these Potamopyrgus mitochondrial genomes is reminiscent of the chloroplast genomes of land plants, harboring a large single-copy (LSC) region and a small single-copy (SSC) region separated by a pair of inverted repeats (IRa and IRb). Individual sequencing reads that spanned across the Potamopyrgus IRa-SSC-IRb structure revealed the occurrence of a "flip-flop" recombination. We also detected evidence for two distinct IR haplotypes and recombination between them in wild-caught P. estuarinus, as well as extensive intermolecular recombination between single-nucleotide polymorphisms in the LSC region. The chloroplast-like architecture and repeat-mediated mitochondrial recombination we describe here raise fundamental questions regarding the origins and commonness of inverted repeats in cytoplasmic genomes and their role in mitochondrial genome evolution.

Host-symbiont interactions in Angomonas deanei include the evolution of a host-derived dynamin ring around the endosymbiont division site

The trypanosomatid Angomonas deanei is a model to study endosymbiosis. Each cell contains a single β-proteobacterial endosymbiont that divides at a defined point in the host cell cycle and contributes essential metabolites to the host metabolism. Additionally, one endosymbiont gene, encoding an ornithine cyclodeaminase (OCD), was transferred by endosymbiotic gene transfer (EGT) to the nucleus. However, the molecular mechanisms mediating the intricate host/symbiont interactions are largely unexplored. Here, we used protein mass spectrometry to identify nucleus-encoded proteins that co-purify with the endosymbiont. Expression of fluorescent fusion constructs of these proteins in A. deanei confirmed seven host proteins to be recruited to specific sites within the endosymbiont. These endosymbiont-targeted proteins (ETPs) include two proteins annotated as dynamin-like protein and peptidoglycan hydrolase that form a ring-shaped structure around the endosymbiont division site that remarkably resembles organellar division machineries. The EGT-derived OCD was not among the ETPs, but instead localizes to the glycosome, likely enabling proline production in the glycosome. We hypothesize that recalibration of the metabolic capacity of the glycosomes that are closely associated with the endosymbiont helps to supply the endosymbiont with metabolites it is auxotrophic for and thus supports the integration of host and endosymbiont metabolic networks. Hence, scrutiny of endosymbiosis-induced protein re-localization patterns in A. deanei yielded profound insights into how an endosymbiotic relationship can stabilize and deepen over time far beyond the level of metabolite exchange.

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.

J-Domain Proteins Orchestrate the Multifunctionality of Hsp70s in Mitochondria: Insights from Mechanistic and Evolutionary Analyses

Mitochondrial J-domain protein (JDP) co-chaperones orchestrate the function of their Hsp70 chaperone partner(s) in critical organellar processes that are essential for cell function. These include folding, refolding, and import of mitochondrial proteins, maintenance of mitochondrial DNA, and biogenesis of iron-sulfur cluster(s) (FeS), prosthetic groups needed for function of mitochondrial and cytosolic proteins. Consistent with the organelle's endosymbiotic origin, mitochondrial Hsp70 and the JDPs' functioning in protein folding and FeS biogenesis clearly descended from bacteria, while the origin of the JDP involved in protein import is less evident. Regardless of their origin, all mitochondrial JDP/Hsp70 systems evolved unique features that allowed them to perform mitochondria-specific functions. Their modes of functional diversification and specialization illustrate the versatility of JDP/Hsp70 systems and inform our understanding of system functioning in other cellular compartments.

Translation in Mitochondrial Ribosomes

Mitochondrial protein synthesis is essential for the life of aerobic eukaryotes. Without it, oxidative phosphorylation cannot be coupled. Evolution has shaped a battery of factors and machinery that are key to production of just a handful of critical proteins. In this general concept chapter, we attempt to briefly summarize our current knowledge of the overall process in mitochondria from a variety of species, breaking this down to the four parts of translation: initiation, elongation, termination, and recycling. Where appropriate, we highlight differences between species and emphasize gaps in our understanding. Excitingly, with the current revolution in cryoelectron microscopy and mitochondrial genome editing, it is highly likely that many of these gaps will be resolved in the near future. However, the absence of a faithful in vitro reconstituted system to study mitochondrial translation is still problematic.

Three Parts of the Plant Genome: On the Way to Success in the Production of Recombinant Proteins

Recombinant proteins are the most important product of current industrial biotechnology. They are indispensable in medicine (for diagnostics and treatment), food and chemical industries, and research. Plant cells combine advantages of the eukaryotic protein production system with simplicity and efficacy of the bacterial one. The use of plants for the production of recombinant proteins is an economically important and promising area that has emerged as an alternative to traditional approaches. This review discusses advantages of plant systems for the expression of recombinant proteins using nuclear, plastid, and mitochondrial genomes. Possibilities, problems, and prospects of modifications of the three parts of the genome in light of obtaining producer plants are examined. Examples of successful use of the nuclear expression platform for production of various biopharmaceuticals, veterinary drugs, and technologically important proteins are described, as are examples of a high yield of recombinant proteins upon modification of the chloroplast genome. Potential utility of plant mitochondria as an expression system for the production of recombinant proteins and its advantages over the nucleus and chloroplasts are substantiated. Although these opportunities have not yet been exploited, potential utility of plant mitochondria as an expression system for the production of recombinant proteins and its advantages over the nucleus and chloroplasts are substantiated.

Evolutionary inference across eukaryotes identifies universal features shaping organelle gene retention

Mitochondria and plastids power complex life. Why some genes and not others are retained in their organelle DNA (oDNA) genomes remains a debated question. Here, we attempt to identify the properties of genes and associated underlying mechanisms that determine oDNA retention. We harness over 15k oDNA sequences and over 300 whole genome sequences across eukaryotes with tools from structural biology, bioinformatics, machine learning, and Bayesian model selection. Previously hypothesized features, including the hydrophobicity of a protein product, and less well-known features, including binding energy centrality within a protein complex, predict oDNA retention across eukaryotes, with additional influences of nucleic acid and amino acid biochemistry. Notably, the same features predict retention in both organelles, and retention models learned from one organelle type quantitatively predict retention in the other, supporting the universality of these features-which also distinguish gene profiles in more recent, independent endosymbiotic relationships. A record of this paper's transparent peer review process is included in the supplemental information.

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.

Structure of a mitochondrial ribosome with fragmented rRNA in complex with membrane-targeting elements

Mitoribosomes of green algae display a great structural divergence from their tracheophyte relatives, with fragmentation of both rRNA and proteins as a defining feature. Here, we report a 2.9 Å resolution structure of the mitoribosome from the alga Polytomella magna harbouring a reduced rRNA split into 13 fragments. We found that the rRNA contains a non-canonical reduced form of the 5S, as well as a permutation of the LSU domain I. The mt-5S rRNA is stabilised by mL40 that is also found in mitoribosomes lacking the 5S, which suggests an evolutionary pathway. Through comparison to other ribosomes with fragmented rRNAs, we observe that the pattern is shared across large evolutionary distances, and between cellular compartments, indicating an evolutionary convergence and supporting the concept of a primordial fragmented ribosome. On the protein level, eleven peripherally associated HEAT-repeat proteins are involved in the binding of 3' rRNA termini, and the structure features a prominent pseudo-trimer of one of them (mL116). Finally, in the exit tunnel, mL128 constricts the tunnel width of the vestibular area, and mL105, a homolog of a membrane targeting component mediates contacts with an inner membrane bound insertase. Together, the structural analysis provides insight into the evolution of the ribosomal machinery in mitochondria.

Interaction of client—the scaffold on which FeS clusters are build—with J-domain protein Hsc20 and its evolving Hsp70 partners

In cells molecular chaperone systems consisting of Hsp70 and its obligatory J-domain protein (JDP) co-chaperones transiently interact with a myriad of client proteins—with JDPs typically recruiting their partner Hsp70 to interact with particular clients. The fundamentals of this cyclical interactions between JDP/Hsp70 systems and clients are well established. Much less is known about other aspects of JDP/Hsp70 system function, including how such systems evolved over time. Here we discuss the JDP/Hsp70 system involved in the biogenesis of iron-sulfur (FeS) clusters. Interaction between the client protein, the scaffold on which clusters are built, and its specialized JDP Hsc20 has stayed constant. However, the system’s Hsp70 has changed at least twice. In some species Hsc20’s Hsp70 partner interacts only with the scaffold, in others it has many JDP partners in addition to Hsc20 and interacts with many client proteins. Analysis of this switching of Hsp70 partners has provided insight into the insulation of JDP/Hsp70 systems from one another that can occur when more than one Hsp70 is present in a cellular compartment, as well as how competition among JDPs is balanced when an Hsp70 partner is shared amongst a number of JDPs. Of particularly broad relevance, even though the scaffold’s interactions with Hsc20 and Hsp70 are functionally critical for the biogenesis of FeS cluster-containing proteins, it is the modulation of the Hsc20-Hsp70 interaction per se that allows Hsc20 to function with such different Hsp70 partners.

Expression of nuclear-encoded, haptophyte-derived ftsH genes support extremely rapid PSII repair and high-light photoacclimation in Karenia brevis (Dinophyceae)

Karenia brevis, a neurotoxic dinoflagellate that produces brevetoxins, is endemic to the Gulf of Mexico and can grow at high irradiances typical of surface waters found there. To build upon a growing number of studies addressing high-light tolerance in K. brevis, specific photobiology and molecular mechanisms underlying this capacity were evaluated in culture. Since photosystem II (PSII) repair cycle activity can be crucial to high light tolerance in plants and algae, the present study assessed this capacity in K. brevis and characterized the ftsH-like genes which are fundamental to this process. Compared with cultures grown in low-light, cultures grown in high-light showed a 65-fold increase in PSII photoinactivation, a ∼50-fold increase in PSII repair, enhanced nonphotochemical quenching (NPQ), and depressed F_v/F_m. Repair rates were among the fastest reported in phytoplankton. Publicly available K. brevis transcriptomes (MMETSP) were queried for ftsH-like sequences and refined with additional sequencing from two K. brevis strains. The genes were phylogenetically related to haptophyte orthologs, implicating acquisition during tertiary endosymbiosis. RT-qPCR of three of the four ftsH-like homologs revealed that poly-A tails predominated in all homologs, and that the most highly expressed homolog had a 5' splice leader and amino-acid motifs characteristic of chloroplast targeting, indicating nuclear encoding for this plastid-targeted gene. High-light cultures showed a ∼1.5-fold upregulation in mRNA expression of the thylakoid-associated genes. Overall, in conjunction with NPQ mechanisms, rapid PSII repair mediated by a haptophyte-derived ftsH prevents chronic photoinhibition in K. brevis. Our findings continue to build the case that high-light photobiology-supported by the acquisition and maintenance of tertiary endosymbiotic genes-is critical to the success of K. brevis in the Gulf of Mexico.

UPRmt and coordinated UPRER in type 2 diabetes

The mitochondrial unfolded protein response (UPR^mt) is a molecular mechanism that maintains mitochondrial proteostasis under stress and is closely related to various metabolic diseases, such as type 2 diabetes (T2D). Similarly, the unfolded protein response of the endoplasmic reticulum (UPR^ER) is responsible for maintaining proteomic stability in the endoplasmic reticulum (ER). Since the mitochondria and endoplasmic reticulum are the primary centers of energy metabolism and protein synthesis in cells, respectively, a synergistic mechanism must exist between UPR^mt and UPR^ER to cooperatively resist stresses such as hyperglycemia in T2D. Increasing evidence suggests that the protein kinase RNA (PKR)-like endoplasmic reticulum kinase (PERK) signaling pathway is likely an important node for coordinating UPR^mt and UPR^ER. The PERK pathway is activated in both UPR^mt and UPR^ER, and its downstream molecules perform important functions. In this review, we discuss the mechanisms of UPR^mt, UPR^ER and their crosstalk in T2D.

Organellar transcripts dominate the cellular mRNA pool across plants of varying ploidy levels

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling.

Inheritance through the cytoplasm

Most heritable information in eukaryotic cells is encoded in the nuclear genome, with inheritance patterns following classic Mendelian segregation. Genomes residing in the cytoplasm, however, prove to be a peculiar exception to this rule. Cytoplasmic genetic elements are generally maternally inherited, although there are several exceptions where these are paternally, biparentally or doubly-uniparentally inherited. In this review, we examine the diversity and peculiarities of cytoplasmically inherited genomes, and the broad evolutionary consequences that non-Mendelian inheritance brings. We first explore the origins of vertical transmission and uniparental inheritance, before detailing the vast diversity of cytoplasmic inheritance systems across Eukaryota. We then describe the evolution of genomic organisation across lineages, how this process has been shaped by interactions with the nuclear genome and population genetics dynamics. Finally, we discuss how both nuclear and cytoplasmic genomes have evolved to co-inhabit the same host cell via one of the longest symbiotic processes, and all the opportunities for intergenomic conflict that arise due to divergence in inheritance patterns. In sum, we cannot understand the evolution of eukaryotes without understanding hereditary symbiosis.

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.