Search

Plant mitochondrial genomes (mitogenomes) experience remarkable levels of horizontal gene transfer (HGT), including the recent discovery that orchids anciently acquired DNA from fungal mitogenomes. Thus far, however, there is no evidence that any of the genes from this interkingdom HGT are functional in orchid mitogenomes. Here, we applied a specialized sequencing approach to the orchid Corallorhiza maculata and found that some fungal-derived tRNA genes in the transferred region are transcribed, post-transcriptionally modified, and aminoacylated. In contrast, all the transferred protein-coding sequences appear to be pseudogenes. These findings show that fungal HGT has altered the composition of the orchid mitochondrial tRNA pool and suggest that these foreign tRNAs function in translation. The exceptional capacity of tRNAs for HGT and functional replacement is further illustrated by the diversity of tRNA genes in the C. maculata mitogenome, which also include genes of plastid and bacterial origin in addition to their native mitochondrial counterparts., (© The Author(s) 2025. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

"Big data" generated from outsourced or centralized facilities often lacks methodological information. Here, we outline how and why researchers, service providers, and other parties should report on methodology and sample metadata to improve scientific reproducibility., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2025 Sloan, Stenglein. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Plant mitochondrial and plastid genomes have exceptionally slow rates of sequence evolution, and recent work has identified an unusual member of the MutS gene family ("plant MSH1 ") as being instrumental in preventing point mutations in these genomes. However, the effects of disrupting MSH1 -mediated DNA repair on "germline" mutation rates have not been quantified. Here, we used Arabidopsis thaliana mutation accumulation (MA) lines to measure mutation rates in msh1 mutants and matched wild type (WT) controls. We detected 124 single nucleotide variants (SNVs: 49 mitochondrial and 75 plastid) and 668 small insertions and deletions (indels: 258 mitochondrial and 410 plastid) in msh1 MA lines. In striking contrast, we did not find any organelle mutations in the WT MA lines, and reanalysis of data from a much larger WT MA experiment also failed to detect any variants. The observed number of SNVs in the msh1 MA lines corresponds to estimated mutation rates of 6.1×10 -7 and 3.2 ×10 -6 per bp per generation in mitochondrial and plastid genomes, respectively. These rates exceed those of species known to have very high mitochondrial mutation rates (e.g., nematodes and fruit flies) by an order of magnitude or more and are on par with estimated rates in humans despite the generation times of A. thaliana being nearly 100-fold shorter. Therefore, disruption of a single plant-specific genetic factor in A. thaliana is sufficient to erase or even reverse the enormous difference in organelle mutation rates between plants and animals.

Eukaryotic nuclear genomes often encode distinct sets of translation machinery for function in the cytosol vs. organelles (mitochondria and plastids). This raises questions about why multiple translation systems are maintained even though they are capable of comparable functions and whether they evolve differently depending on the compartment where they operate. These questions are particularly interesting in plants because translation machinery, including aminoacyl-transfer RNA (tRNA) synthetases (aaRS), is often dual-targeted to the plastids and mitochondria. These organelles have different functions, with much higher rates of translation in plastids to supply the abundant, rapid-turnover proteins required for photosynthesis. Previous studies have indicated that plant organellar aaRS evolve more slowly compared to mitochondrial aaRS in eukaryotes that lack plastids. Thus, we investigated the evolution of nuclear-encoded organellar and cytosolic aaRS and tRNA maturation enzymes across a broad sampling of angiosperms, including nonphotosynthetic (heterotrophic) plant species with reduced plastid gene expression, to test the hypothesis that translational demands associated with photosynthesis constrain the evolution of enzymes involved in organellar tRNA metabolism. Remarkably, heterotrophic plants exhibited wholesale loss of many organelle-targeted aaRS and other enzymes, even though translation still occurs in their mitochondria and plastids. These losses were often accompanied by apparent retargeting of cytosolic enzymes and tRNAs to the organelles, sometimes preserving aaRS-tRNA charging relationships but other times creating surprising mismatches between cytosolic aaRS and mitochondrial tRNA substrates. Our findings indicate that the presence of a photosynthetic plastid drives the retention of specialized systems for organellar tRNA metabolism., Competing Interests: Competing interests statement:The authors declare no competing interest.

The widely distributed MutS gene family functions in recombination, DNA repair, and protein translation. Multiple evolutionary processes have expanded this gene family in plants relative to other eukaryotes. Here, we investigate the origins and functions of these plant-specific genes. Cyanobacterial-like MutS1 and MutS2 genes were ancestrally gained via plastid endosymbiotic gene transfer. MutS1 was subsequently lost in seed plants, whereas MutS2 was duplicated in Viridiplantae (i.e., land plants and green algae). Viridiplantae also have two anciently duplicated copies of the eukaryotic MSH6 gene and acquired MSH1 via horizontal gene transfer - potentially from a nucleocytovirus. Despite sharing a name, "plant MSH1" is not directly related to the MSH1 gene in some fungi and animals, which may be an ancestral eukaryotic gene acquired via mitochondrial endosymbiosis and subsequently lost in most eukaryotes. There has been substantial progress in understanding the functions of plant MSH1 and MSH6 genes, but the cyanobacterial-like MutS1 and MutS2 genes remain uncharacterized. Known functions of bacterial homologs and predicted protein structures, including fusions to diverse nuclease domains, provide hypotheses about potential molecular mechanisms. Because most plant-specific MutS proteins are mitochondrial and/or plastid-targeted, the expansion of this family has played a large role in shaping plant organelle genetics., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

The function of transfer RNAs (tRNAs) depends on enzymes that cleave primary transcript ends, add a 3' CCA tail, introduce post-transcriptional base modifications, and charge (aminoacylate) mature tRNAs with the correct amino acid. Maintaining an available pool of the resulting aminoacylated tRNAs is essential for protein synthesis. High-throughput sequencing techniques have recently been developed to provide a comprehensive view of aminoacylation state in a tRNA-specific fashion. However, these methods have never been applied to plants. Here, we treated Arabidopsis thaliana RNA samples with periodate and then performed tRNA-seq to distinguish between aminoacylated and uncharged tRNAs. This approach successfully captured every tRNA isodecoder family and detected expression of additional tRNA-like transcripts. We found that estimated aminoacylation rates and CCA tail integrity were significantly higher on average for organellar (mitochondrial and plastid) tRNAs than for nuclear/cytosolic tRNAs. Reanalysis of previously published human cell line data showed a similar pattern. Base modifications result in nucleotide misincorporations and truncations during reverse transcription, which we quantified and used to test for relationships with aminoacylation levels. We also determined that the Arabidopsis tRNA-like sequences (t-elements) that are cleaved from the ends of some mitochondrial messenger RNAs have post-transcriptionally modified bases and CCA-tail addition. However, these t-elements are not aminoacylated, indicating that they are only recognized by a subset of tRNA-interacting enzymes and do not play a role in translation. Overall, this work provides a characterization of the baseline landscape of plant tRNA aminoacylation rates and demonstrates an approach for investigating environmental and genetic perturbations to plant translation machinery., (© 2024 The Author(s). The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Intragenomic mutation rates can vary dramatically due to transcription-associated mutagenesis or transcription-coupled repair, which vary based on local epigenomic modifications that are nonuniformly distributed across genomes. One feature associated with decreased mutation is higher expression level, which depends on environmental cues. To understand the magnitude of expression-dependent mutation rate variation, we perturbed expression through a heat treatment in Arabidopsis thaliana. We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using duplex sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty. Somatic mutations in plants may be useful for understanding drivers of DNA damage and repair in the germline since plants experience late germline segregation and both somatic and germline cells share common repair machinery. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild-type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8 × 10-8 and 2.6 × 10-8, respectively), while MMR mutant frequencies were significantly elevated (1.13 × 10-6). Interestingly, in the MMR mutant lines, there was no difference in the somatic mutation frequencies between temperature treatments or between highly versus lowly expressed genes. The extremely low somatic variant frequencies in WT plants indicate that larger datasets will be needed to address fundamental evolutionary questions about whether environmental change leads to gene-specific changes in mutation rate., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Premise: A complicating factor in analyzing allopolyploid genomes is the possibility of physical interactions between homoeologous chromosomes during meiosis, resulting in either crossover (homoeologous exchanges) or non-crossover products (homoeologous gene conversion). Homoeologous gene conversion was first described in cotton by comparing SNP patterns in sequences from two diploid progenitors with those from the allopolyploid subgenomes. These analyses, however, did not explicitly consider other evolutionary scenarios that may give rise to similar SNP patterns as homoeologous gene conversion, creating uncertainties about the reality of the inferred gene conversion events., Methods: Here, we use an expanded phylogenetic sampling of high-quality genome assemblies from seven allopolyploid Gossypium species (all derived from the same polyploidy event), four diploid species (two closely related to each subgenome), and a diploid outgroup to derive a robust method for identifying potential genomic regions of gene conversion and homoeologous exchange., Results: We found little evidence for homoeologous gene conversion in allopolyploid cottons, and that only two of the 40 best-supported events were shared by more than one species. We did, however, reveal a single, shared homoeologous exchange event at one end of chromosome 1, which occurred shortly after allopolyploidization but prior to divergence of the descendant species., Conclusions: Overall, our analyses demonstrated that homoeologous gene conversion and homoeologous exchanges are uncommon in Gossypium, affecting between zero and 24 genes per subgenome (0.0-0.065%) across the seven species. More generally, we highlighted the potential problems of using simple four-taxon tests to investigate patterns of homoeologous gene conversion in established allopolyploids., (© 2024 The Author(s). American Journal of Botany published by Wiley Periodicals LLC on behalf of Botanical Society of America.)

UV light is a potent mutagen that induces bulky DNA damage in the form of cyclobutane pyrimidine dimers (CPDs). Photodamage and other bulky lesions occurring in nuclear genomes can be repaired through nucleotide excision repair (NER), where incisions on both sides of a damaged site precede the removal of a single-stranded oligonucleotide containing the damage. Mitochondrial genomes (mtDNAs) are also susceptible to damage from UV light, but current evidence suggests that the only way to eliminate bulky mtDNA damage is through mtDNA degradation. Damage-containing oligonucleotides excised during NER can be captured with antidamage antibodies and sequenced (XR-seq) to produce high-resolution maps of active repair locations following UV exposure. We analyzed previously published datasets from Arabidopsis thaliana, Saccharomyces cerevisiae, and Drosophila melanogaster to identify reads originating from the mtDNA (and plastid genome in A. thaliana). In A. thaliana and S. cerevisiae, the mtDNA-mapping reads have unique length distributions compared to the nuclear-mapping reads. The dominant fragment size was 26 nt in S. cerevisiae and 28 nt in A. thaliana with distinct secondary peaks occurring in regular intervals. These reads also show a nonrandom distribution of di-pyrimidines (the substrate for CPD formation) with TT enrichment at positions 7-8 of the reads. Therefore, UV damage to mtDNA appears to result in production of DNA fragments of characteristic lengths and positions relative to the damaged location. The mechanisms producing these fragments are unclear, but we hypothesize that they result from a previously uncharacterized DNA degradation pathway or repair mechanism in mitochondria., Competing Interests: Conflicts of interest The author(s) declare no conflicts of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made towards understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana . We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA , recA1 , or recA3 . In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

MutS HOMOLOG 1 (MSH1) is an organellar-targeted protein that obstructs ectopic recombination and the accumulation of mutations in plant organellar genomes. MSH1 also modulates the epigenetic status of nuclear DNA, and its absence induces a variety of phenotypic responses. MSH1 is a member of the MutS family of DNA mismatch repair proteins but harbors an additional GIY-YIG nuclease domain that distinguishes it from the rest of this family. How MSH1 hampers recombination and promotes fidelity in organellar DNA inheritance is unknown. Here, we elucidate its enzymatic activities by recombinantly expressing and purifying full-length MSH1 from Arabidopsis thaliana (AtMSH1). AtMSH1 is a metalloenzyme that shows a strong binding affinity for displacement loops (D-loops). The DNA-binding abilities of AtMSH1 reside in its MutS domain and not in its GIY-YIG domain, which is the ancillary nickase of AtMSH1. In the presence of divalent metal ions, AtMSH1 selectively executes multiple incisions at D-loops, but not other DNA structures including Holliday junctions or dsDNA, regardless of the presence or absence of mismatches. The selectivity of AtMSH1 to dismantle D-loops supports the role of this enzyme in preventing recombination between short repeats. Our results suggest that plant organelles have evolved novel DNA repair routes centered around the anti-recombinogenic activity of MSH1., (© The Author(s) 2023. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis., Competing Interests: Competing interests statement:The authors declare no competing interest.

Mutations are the source of novel genetic diversity but can also lead to disease and maladaptation. The conventional view is that mutations occur randomly with respect to their environment-specific fitness consequences. However, intragenomic mutation rates can vary dramatically due to transcription coupled repair and based on local epigenomic modifications, which are non-uniformly distributed across genomes. One sequence feature associated with decreased mutation is higher expression level, which can vary depending on environmental cues. To understand whether the association between expression level and mutation rate creates a systematic relationship with environment-specific fitness effects, we perturbed expression through a heat treatment in Arabidopsis thaliana . We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using Duplex Sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty in plant mutation research. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8×10 -8 and 2.6×10 -8 , respectively), while MMR mutant frequencies were significantly elevated (1.13×10 -6 ). These results show that somatic variant frequencies are extremely low in WT plants, indicating that larger datasets will be needed to address the fundamental evolutionary question as to whether environmental change leads to gene-specific changes in mutation rate.

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Organelle DNA (oDNA) in mitochondria and plastids is vital for plant (and eukaryotic) life. Selection against damaged oDNA is mediated in part by segregation - sorting different oDNA types into different cells in the germline. Plants segregate oDNA very rapidly, with oDNA recombination protein MSH1 a key driver of this segregation, but we have limited knowledge of the dynamics of this segregation within plants and between generations. Here, we reveal how oDNA evolves through Arabidopsis thaliana development and reproduction. We combine stochastic modelling, Bayesian inference, and model selection with new and existing tissue-specific oDNA measurements from heteroplasmic Arabidopsis plant lines through development and between generations. Segregation proceeds gradually but continually during plant development, with a more rapid increase between inflorescence formation and the next generation. When MSH1 is compromised, the majority of observed segregation can be achieved through partitioning at cell divisions. When MSH1 is functional, mtDNA segregation is far more rapid; we show that increased oDNA gene conversion is a plausible mechanism quantitatively explaining this acceleration. These findings reveal the quantitative, time-dependent details of oDNA segregation in Arabidopsis. We also discuss the support for different models of the plant germline provided by these observations., (© 2023 The Authors New Phytologist © 2023 New Phytologist Foundation.)

UV light is a potent mutagen that induces bulky DNA damage in the form of cyclobutane pyrimidine dimers (CPDs). In eukaryotic cells, photodamage and other bulky lesions occurring in nuclear genomes (nucDNAs) can be repaired through nucleotide excision repair (NER), where dual incisions on both sides of a damaged site precede the removal of a single-stranded oligonucleotide containing the damage. Mitochondrial genomes (mtDNAs) are also susceptible to damage from UV light, but current views hold that the only way to eliminate bulky DNA damage in mtDNAs is through mtDNA degradation. Damage-containing oligonucleotides excised during NER can be captured with anti-damage antibodies and sequenced (XR-seq) to produce high resolution maps of active repair locations following UV exposure. We analyzed previously published datasets from Arabidopsis thaliana, Saccharomyces cerevisiae , and Drosophila melanogaster to identify reads originating from the mtDNA (and plastid genome in A. thaliana ). In A. thaliana and S. cerevisiae , the mtDNA-mapping reads have unique length distributions compared to the nuclear-mapping reads. The dominant fragment size was 26 nt in S. cerevisiae and 28 nt in A. thaliana with distinct secondary peaks occurring in 2-nt ( S. cerevisiae ) or 4-nt ( A. thaliana ) intervals. These reads also show a nonrandom distribution of di-pyrimidines (the substrate for CPD formation) with TT enrichment at positions 7-8 of the reads. Therefore, UV damage to mtDNA appears to result in production of DNA fragments of characteristic lengths and positions relative to the damaged location. We hypothesize that these fragments may reflect the outcome of a previously uncharacterized mechanism of NER-like repair in mitochondria or a programmed mtDNA degradation pathway.

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on Pacific Bioscience HiFi, Hi-C, and Bionano technologies. The assembly produced 10 scaffolds (1 per chromosome) with a total length of 862 Mb and only ∼1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The number of tRNAs encoded in plant mitochondrial genomes varies considerably. Ongoing loss of bacterial-like mitochondrial tRNA genes in many lineages necessitates the import of nuclear-encoded counterparts that share little sequence similarity. Because tRNAs are involved in highly specific molecular interactions, this replacement process raises questions about the identity and trafficking of enzymes necessary for the maturation and function of newly imported tRNAs. In particular, the aminoacyl-tRNA synthetases (aaRSs) that charge tRNAs are usually divided into distinct classes that specialize on either organellar (mitochondrial and plastid) or nuclear-encoded (cytosolic) tRNAs. Here, we investigate the evolution of aaRS subcellular localization in a plant lineage (Sileneae) that has experienced extensive and rapid mitochondrial tRNA loss. By analyzing full-length mRNA transcripts (PacBio Iso-Seq), we found predicted retargeting of many ancestrally cytosolic aaRSs to the mitochondrion and confirmed these results with colocalization microscopy assays. However, we also found cases where aaRS localization does not appear to change despite functional tRNA replacement, suggesting evolution of novel interactions and charging relationships. Therefore, the history of repeated tRNA replacement in Sileneae mitochondria reveals that differing constraints on tRNA/aaRS interactions may determine which of these alternative coevolutionary paths is used to maintain organellar translation in plant cells., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

There is growing evidence that cytonuclear incompatibilities (i.e. disruption of cytonuclear coadaptation) might contribute to the speciation process. In a former study, we described the possible involvement of plastid-nuclear incompatibilities in the reproductive isolation between four lineages of Silene nutans (Caryophyllaceae). Because organellar genomes are usually cotransmitted, we assessed whether the mitochondrial genome could also be involved in the speciation process, knowing that the gynodioecious breeding system of S. nutans is expected to impact the evolutionary dynamics of this genome. Using hybrid capture and high-throughput DNA sequencing, we analyzed diversity patterns in the genic content of the organellar genomes in the four S. nutans lineages. Contrary to the plastid genome, which exhibited a large number of fixed substitutions between lineages, extensive sharing of polymorphisms between lineages was found in the mitochondrial genome. In addition, numerous recombination-like events were detected in the mitochondrial genome, loosening the linkage disequilibrium between the organellar genomes and leading to decoupled evolution. These results suggest that gynodioecy shaped mitochondrial diversity through balancing selection, maintaining ancestral polymorphism and, thus, limiting the involvement of the mitochondrial genome in evolution of hybrid inviability between S. nutans lineages., (© 2023 The Authors New Phytologist © 2023 New Phytologist Foundation.)

Mitochondrial DNA heteroplasmy samples can shed light on vital developmental and genetic processes shaping mitochondrial DNA populations. The sample means and sample variance of a set of heteroplasmy observations are typically used both to estimate bottleneck sizes and to perform fits to the theoretical "Kimura" distribution in seeking evidence for mitochondrial DNA selection. However, each of these applications raises problems. Sample statistics do not generally provide optimal fits to the Kimura distribution and so can give misleading results in hypothesis testing, including false positive signals of selection. Using sample variance can give misleading results for bottleneck size estimates, particularly for small samples. These issues can and do lead to false positive results for mitochondrial DNA mechanisms-all published experimental datasets we re-analyzed, reported as displaying departures from the Kimura model, do not in fact give evidence for such departures. Here we outline a maximum likelihood approach that is simple to implement computationally and addresses all of these issues. We advocate the use of maximum likelihood fits and explicit hypothesis tests, not fits and Kolmogorov-Smirnov tests via summary statistics, for ongoing work with mitochondrial DNA heteroplasmy., Competing Interests: Conflicts of interest The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Eukaryotes maintain separate protein translation systems for nuclear and organellar genes, including distinct sets of tRNAs and aminoacyl-tRNA synthetases (aaRSs). In animals, mitochondrial-targeted aaRSs are expressed at lower levels and are less conserved in sequence than cytosolic aaRSs involved in translation of nuclear mRNAs, likely reflecting lower translational demands in mitochondria. In plants, translation is further complicated by the presence of plastids, which share most aaRSs with mitochondria. In addition, plant mitochondrial tRNA pools have a dynamic history of gene loss and functional replacement by tRNAs from other compartments. To investigate the consequences of these distinctive features of translation in plants, we analyzed sequence evolution in angiosperm aaRSs. In contrast to previously studied eukaryotic systems, we found that plant organellar and cytosolic aaRSs exhibit only a small difference in expression levels, and organellar aaRSs are slightly more conserved than cytosolic aaRSs. We hypothesize that these patterns result from high translational demands associated with photosynthesis in mature chloroplasts. We also investigated aaRS evolution in Sileneae, an angiosperm lineage with extensive mitochondrial tRNA replacement and aaRS retargeting. We predicted positive selection for changes in aaRS sequence resulting from these recent changes in subcellular localization and tRNA substrates but found little evidence for accelerated sequence divergence. Overall, the complex tripartite translation system in plant cells appears to have imposed more constraints on the long-term evolutionary rates of organellar aaRSs compared with other eukaryotic lineages, and plant aaRS protein sequences appear largely robust to more recent perturbations in subcellular localization and tRNA interactions., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Society for Molecular Biology and 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., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The abundant repeats in plant mitochondrial genomes can cause rapid genome rearrangements and are also a major obstacle in short-read sequencing studies. Nuclear-encoded proteins such as MSH1 are known to suppress the generation of repeat-associated mitochondrial genome variants, but our understanding of these mechanisms has been constrained by the limitations of short-read technologies. Here, we used highly accurate long-read sequencing (PacBio HiFi) to characterize mitochondrial and plastid genome variants in Arabidopsis thaliana msh1 mutant individuals. The HiFi reads provided a global view of recombination dynamics with detailed quantification of parental and crossover recombination products for both large and small repeats. We found that recombination breakpoints were distributed relatively evenly across the length of repeated sequences and detected widespread internal exchanges of sequence variants between pairs of imperfect repeats in the mitochondrial genome of msh1 mutants. Long-read assemblies of mitochondrial genomes from seven other A. thaliana wild-type accessions differed by repeat-mediated structural rearrangements similar to those observed in msh1 mutants, but they were all in a simple low-heteroplasmy state. The Arabidopsis plastid genome generally lacks small repeats and exhibited a very different pattern of variant accumulation in msh1 mutants compared with the mitochondrial genome. Our data illustrate the power of HiFi technology in studying repeat-mediated recombination in plant organellar genomes and improved the sequence resolution for recombinational processes suppressed by MSH1. Plant organellar genomes can undergo rapid rearrangements. Long-read sequencing provides a detailed and quantitative view of mitochondrial and plastid genome variants normally suppressed by MSH1, advancing our understanding of plant organellar genome dynamics., (© 2022 Society for Experimental Biology and John Wiley & Sons Ltd.)

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

The fate of new mitochondrial and plastid mutations depends on their ability to persist and spread among the numerous organellar genome copies within a cell (heteroplasmy). The extent to which heteroplasmies are transmitted across generations or eliminated through genetic bottlenecks is not well understood in plants, in part because their low mutation rates make these variants so infrequent. Disruption of MutS Homolog 1 ( MSH1 ), a gene involved in plant organellar DNA repair, results in numerous de novo point mutations, which we used to quantitatively track the inheritance of single nucleotide variants in mitochondrial and plastid genomes in Arabidopsis . We found that heteroplasmic sorting (the fixation or loss of a variant) was rapid for both organelles, greatly exceeding rates observed in animals. In msh1 mutants, plastid variants sorted faster than those in mitochondria and were typically fixed or lost within a single generation. Effective transmission bottleneck sizes ( N ) for plastids and mitochondria were N ∼ 1 and 4, respectively. Restoring MSH1 function further increased the rate of heteroplasmic sorting in mitochondria ( N ∼ 1.3), potentially because of its hypothesized role in promoting gene conversion as a mechanism of DNA repair, which is expected to homogenize genome copies within a cell. Heteroplasmic sorting also favored GC base pairs. Therefore, recombinational repair and gene conversion in plant organellar genomes can potentially accelerate the elimination of heteroplasmies and bias the outcome of this sorting process.

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.

Intracellular transfers of mitochondrial DNA continue to shape nuclear genomes. Chromosome 2 of the model plant Arabidopsis thaliana contains one of the largest known nuclear insertions of mitochondrial DNA (numts). Estimated at over 600 kb in size, this numt is larger than the entire Arabidopsis mitochondrial genome. The primary Arabidopsis nuclear reference genome contains less than half of the numt because of its structural complexity and repetitiveness. Recent data sets generated with improved long-read sequencing technologies (PacBio HiFi) provide an opportunity to finally determine the accurate sequence and structure of this numt. We performed a de novo assembly using sequencing data from recent initiatives to span the Arabidopsis centromeres, producing a gap-free sequence of the Chromosome 2 numt, which is 641 kb in length and has 99.933% nucleotide sequence identity with the actual mitochondrial genome. The numt assembly is consistent with the repetitive structure previously predicted from fiber-based fluorescent in situ hybridization. Nanopore sequencing data indicate that the numt has high levels of cytosine methylation, helping to explain its biased spectrum of nucleotide sequence divergence and supporting previous inferences that it is transcriptionally inactive. The original numt insertion appears to have involved multiple mitochondrial DNA copies with alternative structures that subsequently underwent an additional duplication event within the nuclear genome. This work provides insights into numt evolution, addresses one of the last unresolved regions of the Arabidopsis reference genome, and represents a resource for distinguishing between highly similar numt and mitochondrial sequences in studies of transcription, epigenetic modifications, and de novo mutations., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

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

The chloroplast chaperone CLPC1 unfolds and delivers substrates to the stromal CLPPRT protease complex for degradation. We previously used an in vivo trapping approach to identify interactors with CLPC1 in Arabidopsis thaliana by expressing a STREPII-tagged copy of CLPC1 mutated in its Walker B domains (CLPC1-TRAP) followed by affinity purification and mass spectrometry. To create a larger pool of candidate substrates, adaptors, or regulators, we carried out a far more sensitive and comprehensive in vivo protein trapping analysis. We identified 59 highly enriched CLPC1 protein interactors, in particular proteins belonging to families of unknown functions (DUF760, DUF179, DUF3143, UVR-DUF151, HugZ/DUF2470), as well as the UVR domain proteins EXE1 and EXE2 implicated in singlet oxygen damage and signaling. Phylogenetic and functional domain analyses identified other members of these families that appear to localize (nearly) exclusively to plastids. In addition, several of these DUF proteins are of very low abundance as determined through the Arabidopsis PeptideAtlas http://www.peptideatlas.org/builds/arabidopsis/ showing that enrichment in the CLPC1-TRAP was extremely selective. Evolutionary rate covariation indicated that the HugZ/DUF2470 family coevolved with the plastid CLP machinery suggesting functional and/or physical interactions. Finally, mRNA-based coexpression networks showed that all 12 CLP protease subunits tightly coexpressed as a single cluster with deep connections to DUF760-3. Coexpression modules for other trapped proteins suggested specific functions in biological processes, e.g., UVR2 and UVR3 were associated with extraplastidic degradation, whereas DUF760-6 is likely involved in senescence. This study provides a strong foundation for discovery of substrate selection by the chloroplast CLP protease system., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

While the chloroplast (plastid) is known for its role in photosynthesis, it is also involved in many other metabolic pathways essential for plant survival. As such, plastids contain an extensive suite of enzymes required for non-photosynthetic processes. The evolution of the associated genes has been especially dynamic in flowering plants (angiosperms), including examples of gene duplication and extensive rate variation. We examined the role of ongoing gene duplication in two key plastid enzymes, the acetyl-CoA carboxylase (ACCase) and the caseinolytic protease (Clp), responsible for fatty acid biosynthesis and protein turnover, respectively. In plants, there are two ACCase complexes-a homomeric version present in the cytosol and a heteromeric version present in the plastid. Duplications of the nuclear-encoded homomeric ACCase gene and retargeting of one resultant protein to the plastid have been previously reported in multiple species. We find that these retargeted homomeric ACCase proteins exhibit elevated rates of sequence evolution, consistent with neofunctionalization and/or relaxation of selection. The plastid Clp complex catalytic core is composed of nine paralogous proteins that arose via ancient gene duplication in the cyanobacterial/plastid lineage. We show that further gene duplication occurred more recently in the nuclear-encoded core subunits of this complex, yielding additional paralogs in many species of angiosperms. Moreover, in six of eight cases, subunits that have undergone recent duplication display increased rates of sequence evolution relative to those that have remained single copy. We also compared substitution patterns between pairs of Clp core paralogs to gain insight into post-duplication evolutionary routes. These results show that gene duplication and rate variation continue to shape the plastid proteome., (Copyright © 2022 Elsevier Inc. All rights reserved.)

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

In most eukaryotes, transfer RNAs (tRNAs) are one of the very few classes of genes remaining in the mitochondrial genome, but some mitochondria have lost these vestiges of their prokaryotic ancestry. Sequencing of mitogenomes from the flowering plant genus Silene previously revealed a large range in tRNA gene content, suggesting rapid and ongoing gene loss/replacement. Here, we use this system to test longstanding hypotheses about how mitochondrial tRNA genes are replaced by importing nuclear-encoded tRNAs. We traced the evolutionary history of these gene loss events by sequencing mitochondrial genomes from key outgroups (Agrostemma githago and Silene [=Lychnis] chalcedonica). We then performed the first global sequencing of purified plant mitochondrial tRNA populations to characterize the expression of mitochondrial-encoded tRNAs and the identity of imported nuclear-encoded tRNAs. We also confirmed the utility of high-throughput sequencing methods for the detection of tRNA import by sequencing mitochondrial tRNA populations in a species (Solanum tuberosum) with known tRNA trafficking patterns. Mitochondrial tRNA sequencing in Silene revealed substantial shifts in the abundance of some nuclear-encoded tRNAs in conjunction with their recent history of mt-tRNA gene loss and surprising cases where tRNAs with anticodons still encoded in the mitochondrial genome also appeared to be imported. These data suggest that nuclear-encoded counterparts are likely replacing mitochondrial tRNAs even in systems with recent mitochondrial tRNA gene loss, and the redundant import of a nuclear-encoded tRNA may provide a mechanism for functional replacement between translation systems separated by billions of years of evolutionary divergence., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e., the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CG→AT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Furthermore, we found an excess of G→T and C→T changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

The plant genome is partitioned across three distinct subcellular compartments: the nucleus, mitochondria, and plastids. Successful coordination of gene expression among these organellar genomes and the nuclear genome is critical for plant function and fitness. Whole genome duplication (WGD) events in the nucleus have played a major role in the diversification of land plants and are expected to perturb the relative copy number (stoichiometry) of nuclear, mitochondrial, and plastid genomes. Thus, elucidating the mechanisms whereby plant cells respond to the cytonuclear stoichiometric imbalance that follows WGDs represents an important yet underexplored question in understanding the evolutionary consequences of genome doubling. We used droplet digital PCR to investigate the relationship between nuclear and organellar genome copy numbers in allopolyploids and their diploid progenitors in both wheat and Arabidopsis. Polyploids exhibit elevated organellar genome copy numbers per cell, largely preserving the cytonuclear stoichiometry observed in diploids despite the change in nuclear genome copy number. To investigate the timescale over which cytonuclear stoichiometry may respond to WGD, we also estimated the organellar genome copy number in Arabidopsis synthetic autopolyploids and in a haploid-induced diploid line. We observed corresponding changes in organellar genome copy number in these laboratory-generated lines, indicating that at least some of the cellular response to cytonuclear stoichiometric imbalance is immediate following WGD. We conclude that increases in organellar genome copy numbers represent a common response to polyploidization, suggesting that maintenance of cytonuclear stoichiometry is an important component in establishing polyploid lineages., (© 2021 Society for Experimental Biology and John Wiley & Sons Ltd.)

The angiosperm genus Silene is a model system for several traits of ecological and evolutionary significance in plants, including breeding system and sex chromosome evolution, host-pathogen interactions, invasive species biology, heavy metal tolerance, and cytonuclear interactions. Despite its importance, genomic resources for this large genus of approximately 850 species are scarce, with only one published whole-genome sequence (from the dioecious species Silene latifolia). Here, we provide genomic and transcriptomic resources for a hermaphroditic representative of this genus (S. noctiflora), including a PacBio Iso-Seq transcriptome, which uses long-read, single-molecule sequencing technology to analyze full-length mRNA transcripts. Using these data, we have assembled and annotated high-quality full-length cDNA sequences for approximately 14,126 S. noctiflora genes and 25,317 isoforms. We demonstrated the utility of these data to distinguish between recent and highly similar gene duplicates by identifying novel paralogous genes in an essential protease complex. Furthermore, we provide a draft assembly for the approximately 2.7-Gb genome of this species, which is near the upper range of genome-size values reported for diploids in this genus and threefold larger than the 0.9-Gb genome of Silene conica, another species in the same subgenus. Karyotyping confirmed that S. noctiflora is a diploid, indicating that its large genome size is not due to polyploidization. These resources should facilitate further study and development of this genus as a model in plant ecology and evolution., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Nuclear and plastid (chloroplast) genomes experience different mutation rates, levels of selection, and transmission modes, yet key cellular functions depend on their coordinated interactions. Functionally related proteins often show correlated changes in rates of sequence evolution across a phylogeny [evolutionary rate covariation (ERC)], offering a means to detect previously unidentified suites of coevolving and cofunctional genes. We performed phylogenomic analyses across angiosperm diversity, scanning the nuclear genome for genes that exhibit ERC with plastid genes. As expected, the strongest hits were highly enriched for genes encoding plastid-targeted proteins, providing evidence that cytonuclear interactions affect rates of molecular evolution at genome-wide scales. Many identified nuclear genes functioned in post-transcriptional regulation and the maintenance of protein homeostasis (proteostasis), including protein translation (in both the plastid and cytosol), import, quality control, and turnover. We also identified nuclear genes that exhibit strong signatures of coevolution with the plastid genome, but their encoded proteins lack organellar-targeting annotations, making them candidates for having previously undescribed roles in plastids. In sum, our genome-wide analyses reveal that plastid-nuclear coevolution extends beyond the intimate molecular interactions within chloroplast enzyme complexes and may be driven by frequent rewiring of the machinery responsible for maintenance of plastid proteostasis in angiosperms., (� American Society of Plant Biologists 2021. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Although plant mitochondrial genomes typically show low rates of sequence evolution, levels of divergence in certain angiosperm lineages suggest anomalously high mitochondrial mutation rates. However, de novo mutations have never been directly analyzed in such lineages. Recent advances in high-fidelity DNA sequencing technologies have enabled detection of mitochondrial mutations when still present at low heteroplasmic frequencies. To date, these approaches have only been performed on a single plant species (Arabidopsis thaliana). Here, we apply a high-fidelity technique (Duplex Sequencing) to multiple angiosperms from the genus Silene, which exhibits extreme heterogeneity in rates of mitochondrial sequence evolution among close relatives. Consistent with phylogenetic evidence, we found that Silene latifolia maintains low mitochondrial variant frequencies that are comparable with previous measurements in Arabidopsis. Silene noctiflora also exhibited low variant frequencies despite high levels of historical sequence divergence, which supports other lines of evidence that this species has reverted to lower mitochondrial mutation rates after a past episode of acceleration. In contrast, S. conica showed much higher variant frequencies in mitochondrial (but not in plastid) DNA, consistent with an ongoing bout of elevated mitochondrial mutation rates. Moreover, we found an altered mutational spectrum in S. conica heavily biased towards AT→GC transitions. We also observed an unusually low number of mitochondrial genome copies per cell in S. conica, potentially pointing to reduced opportunities for homologous recombination to accurately repair mismatches in this species. Overall, these results suggest that historical fluctuations in mutation rates are driving extreme variation in rates of plant mitochondrial sequence evolution., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Compared with free-living bacteria, endosymbionts of sap-feeding insects have tiny and rapidly evolving genomes. Increased genetic drift, high mutation rates, and relaxed selection associated with host control of key cellular functions all likely contribute to genome decay. Phylogenetic comparisons have revealed massive variation in endosymbiont evolutionary rate, but such methods make it difficult to partition the effects of mutation versus selection. For example, the ancestor of Auchenorrhynchan insects contained two obligate endosymbionts, Sulcia and a betaproteobacterium (BetaSymb; called Nasuia in leafhoppers) that exhibit divergent rates of sequence evolution and different propensities for loss and replacement in the ensuing ∼300 Ma. Here, we use the auchenorrhynchan leafhopper Macrosteles sp. nr. severini, which retains both of the ancestral endosymbionts, to test the hypothesis that differences in evolutionary rate are driven by differential mutagenesis. We used a high-fidelity technique known as duplex sequencing to measure and compare low-frequency variants in each endosymbiont. Our direct detection of de novode novo mutations reveals that the rapidly evolving endosymbiont (Nasuia) has a much higher frequency of single-nucleotide variants than the more stable endosymbiont (Sulcia) and a mutation spectrum that is potentially even more AT-biased than implied by the 83.1% AT content of its genome. We show that indels are common in both endosymbionts but differ substantially in length and distribution around repetitive regions. Our results suggest that differences in long-term rates of sequence evolution in Sulcia versus BetaSymb, and perhaps the contrasting degrees of stability of their relationships with the host, are driven by differences in mutagenesis., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Although tRNA structure is one of the most conserved and recognizable shapes in molecular biology, aberrant tRNAs are frequently found in the mitochondrial genomes of metazoans. The extremely degenerate structures of several mitochondrial tRNAs (mt-tRNAs) have led to doubts about their expression and function. Mites from the arachnid superorder Acariformes are predicted to have some of the shortest mt-tRNAs, with a complete loss of cloverleaf-like shape. While performing mitochondrial isolations and recently developed tRNA-seq methods in plant tissue, we inadvertently sequenced the mt-tRNAs from a common plant pest, the acariform mite Tetranychus urticae , to a high enough coverage to detect all previously annotated T. urticae tRNA regions. The results not only confirm expression, CCA-tailing and post-transcriptional base modification of these highly divergent tRNAs, but also revealed paired sense and antisense expression of multiple T. urticae mt-tRNAs. Mirrored expression of mt-tRNA genes has been hypothesized but not previously demonstrated to be common in any system. We discuss the functional roles that these divergent tRNAs could have as both decoding molecules in translation and processing signals in transcript maturation pathways, as well as how sense-antisense pairs add another dimension to the bizarre tRNA biology of mitochondrial genomes., (© The Author(s) 2021. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana , which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair ( POLIA , POLIB , MSH1 , RECA3 , UNG , FPG , and OGG1 ) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisen mitochondrial and plastid mutations even at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineage within the mutS mismatch repair family that we find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of double-stranded breaks and repair pathways based on faithful homologous recombination., Competing Interests: The authors declare no competing interest.

Transfer RNAs (tRNAs) remain one of the very few classes of genes still encoded in the mitochondrial genome. These key components of the protein translation system must interact with a large enzymatic network of nuclear-encoded gene products to maintain mitochondrial function. Plants have an evolutionarily dynamic mitochondrial tRNA population, including ongoing tRNA gene loss and replacement by both horizontal gene transfer from diverse sources and import of nuclear-expressed tRNAs from the cytosol. Thus, plant mitochondria represent an excellent model for understanding how anciently divergent genes can act as "interchangeable parts" during the evolution of complex molecular systems. In particular, understanding the integration of the mitochondrial translation system with elements of the corresponding machinery used in cytosolic protein synthesis is a key area for eukaryotic cellular evolution. Here, we review the increasingly detailed phylogenetic data about the evolutionary history of mitochondrial tRNA gene loss, transfer, and functional replacement that has created extreme variation in mitochondrial tRNA populations across plant species. We describe emerging tRNA-seq methods with promise for refining our understanding of the expression and subcellular localization of tRNAs. Finally, we summarize current evidence and identify open questions related to coevolutionary changes in nuclear-encoded enzymes that have accompanied turnover in mitochondrial tRNA populations., (Copyright © 2020 Elsevier B.V. and Mitochondria Research Society. All rights reserved.)

DNA sequence alignments have provided the majority of data for inferring phylogenetic relationships with both concatenation and coalescent methods. However, DNA sequences are susceptible to extensive homoplasy, especially for deep divergences in the Tree of Life. Retroelement insertions have emerged as a powerful alternative to sequences for deciphering evolutionary relationships because these data are nearly homoplasy-free. In addition, retroelement insertions satisfy the "no intralocus-recombination" assumption of summary coalescent methods because they are singular events and better approximate neutrality relative to DNA loci commonly sampled in phylogenomic studies. Retroelements have traditionally been analyzed with parsimony, distance, and network methods. Here, we analyze retroelement data sets for vertebrate clades (Placentalia, Laurasiatheria, Balaenopteroidea, Palaeognathae) with 2 ILS-aware methods that operate by extracting, weighting, and then assembling unrooted quartets into a species tree. The first approach constructs a species tree from retroelement bipartitions with ASTRAL, and the second method is based on split-decomposition with parsimony. We also develop a Quartet-Asymmetry test to detect hybridization using retroelements. Both ILS-aware methods recovered the same species-tree topology for each data set. The ASTRAL species trees for Laurasiatheria have consecutive short branch lengths in the anomaly zone whereas Palaeognathae is outside of this zone. For the Balaenopteroidea data set, which includes rorquals (Balaenopteridae) and gray whale (Eschrichtiidae), both ILS-aware methods resolved balaeonopterids as paraphyletic. Application of the Quartet-Asymmetry test to this data set detected 19 different quartets of species for which historical introgression may be inferred. Evidence for introgression was not detected in the other data sets., (© The American Genetic Association 2019. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Introgressive hybridization results in the transfer of genetic material between species, often with fitness implications for the recipient species. The development of statistical methods for detecting the signatures of historical introgression in whole-genome data has been a major area of focus. Although existing techniques are able to identify the taxa that exchanged genes during introgression using a four-taxon system, most methods do not explicitly distinguish which taxon served as donor and which as recipient during introgression (i.e., polarization of introgression directionality). Existing methods that do polarize introgression are often only able to do so when there is a fifth taxon available and that taxon is sister to one of the taxa involved in introgression. Here, we present divergence-based introgression polarization (DIP), a method for polarizing introgression using patterns of sequence divergence across whole genomes, which operates in a four-taxon context. Thus, DIP can be applied to infer the directionality of introgression when additional taxa are not available. We use simulations to show that DIP can polarize introgression and identify potential sources of bias in the assignment of directionality, and we apply DIP to a well-described hominin introgression event., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The mechanisms of sequence divergence in angiosperm mitochondrial genomes have long been enigmatic. In particular, it is difficult to reconcile the rapid divergence of intergenic regions that can make non-coding sequences almost unrecognizable even among close relatives with the unusually high levels of sequence conservation found in genic regions. It has been hypothesized that different mutation and repair mechanisms act on genic and intergenic sequences or alternatively that mutational input is relatively constant but that selection has strikingly different effects on these respective regions. To test these alternative possibilities, we analyzed mtDNA divergence within Arabidopsis thaliana , including variants from the 1001 Genomes Project and changes accrued in published mutation accumulation (MA) lines. We found that base-substitution frequencies are relatively similar for intergenic regions and synonymous sites in coding regions, whereas indel and nonsynonymous substitutions rates are greatly depressed in coding regions, supporting a conventional model in which mutation/repair mechanisms are consistent throughout the genome but differentially filtered by selection. Most types of sequence and structural changes were undetectable in 10-generation MA lines, but we found significant shifts in relative copy number across mtDNA regions for lines grown under stressed vs. benign conditions. We confirmed quantitative variation in copy number across the A. thaliana mitogenome using both whole-genome sequencing and droplet digital PCR, further undermining the classic but oversimplified model of a circular angiosperm mtDNA structure. Our results suggest that copy number variation is one of the most fluid features of angiosperm mitochondrial genomes., (Copyright © 2020 Wu et al.)

Background: Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods., Results: Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception., Conclusions: These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Angiosperm mitochondrial (mt) genes are generally slow-evolving, but multiple lineages have undergone dramatic accelerations in rates of nucleotide substitution and extreme changes in mt genome structure. While molecular evolution in these lineages has been investigated, very little is known about their mt function. Some studies have suggested altered respiration in individual taxa, although there are several reasons why mt variation might be neutral in others. Here, we develop a new protocol to characterize respiration in isolated plant mitochondria and apply it to species of Silene with mt genomes that are rapidly evolving, highly fragmented, and exceptionally large (~11 Mbp). This protocol, complemented with traditional measures of plant fitness, cytochrome c oxidase activity assays, and fluorescence microscopy, was also used to characterize inter- and intraspecific variation in mt function. Contributions of the individual "classic" OXPHOS complexes, the alternative oxidase, and external NADH dehydrogenases to overall mt respiratory flux were found to be similar to previously studied angiosperms with more typical mt genomes. Some differences in mt function could be explained by inter- and intraspecific variation. This study suggests that Silene species with peculiar mt genomes still show relatively normal mt respiration. This may be due to strong purifying selection on mt variants, coevolutionary responses in the nucleus, or a combination of both. Future experiments should explore such questions using a comparative framework and investigating other lineages with unusual mitogenomes., (Copyright © 2019 Elsevier B.V. and Mitochondria Research Society. All rights reserved.)

Genomic datasets sometimes support conflicting phylogenetic relationships when different tree-building methods are applied. Coherent interpretations of such results are enabled by partitioning support for controversial relationships among the constituent genes of a phylogenomic dataset. For the supermatrix (=concatenation) approach, several methods that measure the distribution of support and conflict among loci were introduced over 15 years ago. More recently, partitioned coalescence support (PCS) was developed for phylogenetic coalescence methods that account for incomplete lineage sorting and use the summed fits of gene trees to estimate the species tree. Here, we automate computation of PCS to permit application of this index to genome-scale matrices that include hundreds of loci. Reanalyses of four phylogenomic datasets for amniotes, land plants, skinks, and angiosperms demonstrate how PCS scores can be used to: (1) compare conflicting results favored by alternative coalescence methods, (2) identify outlier gene trees that have a disproportionate influence on the resolution of contentious relationships, (3) assess the effects of missing data in species-tree analysis, and (4) clarify biases in commonly-implemented coalescence methods and support indices. We show that key phylogenomic conclusions from these analyses often hinge on just a few gene trees and that results can be driven by specific biases of a particular coalescence method and/or the differential weight placed on gene trees with high versus low taxon sampling. The attribution of exceptionally high weight to some gene trees and very low weight to other gene trees counters the basic logic of phylogenomic coalescence analysis; even clades in species trees with high support according to commonly used indices (likelihood-ratio test, bootstrap, Bayesian local posterior probability) can be unstable to the removal of only one or two gene trees with high PCS. Computer simulations cannot adequately describe all of the contingencies and complexities of empirical genetic data. PCS scores complement simulation work by providing specific insights into a particular dataset given the assumptions of the phylogenetic coalescence method that is applied. In combination with standard measures of nodal support, PCS provides a more complete understanding of the overall genomic evidence for contested evolutionary relationships in species trees., (Copyright © 2019 Elsevier Inc. All rights reserved.)