MSH1 is required for maintenance of the low mutation rates in plant mitochondrial and plastid genomes

Plant Organellar MSH1 Is a Displacement Loop-Specific Endonuclease

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.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BACKGROUND

Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.

RESULTS

We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.

CONCLUSIONS

The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Plasticity of repetitive sequences demonstrated by the complete mitochondrial genome of Eucalyptus camaldulensis

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

Collective mitochondrial dynamics resolve conflicting cellular tensions: From plants to general principles

Mitochondria play diverse and essential roles in eukaryotic cells, and plants are no exception. Plant mitochondria have several differences from their metazoan and fungal cousins: they often exist in a fragmented state, move rapidly on actin rather than microtubules, have many plant-specific metabolic features and roles, and usually contain only a subset of the complete mtDNA genome, which itself undergoes frequent recombination. This arrangement means that exchange and complementation is essential for plant mitochondria, and recent work has begun to reveal how their collective dynamics and resultant "social networks" of encounters support this exchange, connecting plant mitochondria in time rather than in space. This review will argue that this social network perspective can be extended to a "societal network", where mitochondrial dynamics are an essential part of the interacting cellular society of organelles and biomolecules. Evidence is emerging that mitochondrial dynamics allow optimal resolutions to competing cellular priorities; we will survey this evidence and review potential future research directions, highlighting that plant mitochondria can help reveal and test principles that apply across other kingdoms of life. In parallel with this fundamental cell biology, we also highlight the translational "One Health" importance of plant mitochondrial behaviour - which is exploited in the production of a vast amount of crops consumed worldwide - and the potential for multi-objective optimisation to understand and rationally re-engineer the evolved resolutions to these tensions.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA

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.

Investigating low frequency somatic mutations in Arabidopsis with Duplex Sequencing

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.

Stochastic organelle genome segregation through Arabidopsis development and reproduction

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.

UV damage induces production of mitochondrial DNA fragments with specific length profiles

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.

Base Excision DNA Repair in Plants: Arabidopsis and Beyond

Base excision DNA repair (BER) is a key pathway safeguarding the genome of all living organisms from damage caused by both intrinsic and environmental factors. Most present knowledge about BER comes from studies of human cells, E. coli, and yeast. Plants may be under an even heavier DNA damage threat from abiotic stress, reactive oxygen species leaking from the photosynthetic system, and reactive secondary metabolites. In general, BER in plant species is similar to that in humans and model organisms, but several important details are specific to plants. Here, we review the current state of knowledge about BER in plants, with special attention paid to its unique features, such as the existence of active epigenetic demethylation based on the BER machinery, the unexplained diversity of alkylation damage repair enzymes, and the differences in the processing of abasic sites that appear either spontaneously or are generated as BER intermediates. Understanding the biochemistry of plant DNA repair, especially in species other than the Arabidopsis model, is important for future efforts to develop new crop varieties.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Characterization and development of a plastid genome base editor, ptpTALECD

The modification of photosynthesis-related genes in plastid genomes may improve crop yields. Recently, we reported that a plastid-targeting base editor named ptpTALECD, in which a cytidine deaminase DddA functions as the catalytic domain, can homoplasmically substitute a targeted C to T in plastid genomes of Arabidopsis thaliana. However, some target Cs were not substituted. In addition, although ptpTALECD could substitute Cs on the 3' side of T and A, it was unclear whether it could also substitute Cs on the 3' side of G and C. In this study, we identified the preferential positions of the substituted Cs in ptpTALECD-targeting sequences in the Arabidopsis plastid genome. We also found that ptpTALECD could substitute Cs on the 3' side of all four bases in plastid genomes of Arabidopsis. More recently, a base editor containing an improved version of DddA (DddA11) was reported to substitute Cs more efficiently, and to substitute Cs on the 3' side of more varieties of bases in human mitochondrial genomes than a base editor containing DddA. Here, we also show that ptpTALECD_v2, in which a modified version of DddA11 functions as the catalytic domain, more frequently substituted Cs than ptpTALECD in the Arabidopsis plastid genome. We also found that ptpTALECD_v2 tended to substitute Cs at more positions than ptpTALECD. Our results reveal that ptpTALECD can cause a greater variety of codon changes and amino acid substitutions than previously thought, and that ptpTALECD and ptpTALECD_v2 are useful tools for the targeted base editing of plastid genomes.

Transcriptional Profiling Reveals the Wheat Defences against Fusarium Head Blight Disease Regulated by a NAC Transcription Factor

The wheat NAC transcription factor TaNACL-D1 enhances resistance to the economically devastating Fusarium head blight (FHB) disease. The objective of this study was to decipher the alterations in gene expression, pathways and biological processes that led to enhanced resistance as a result of the constitutive expression of TaNACL-D1 in wheat. Transcriptomic analysis was used to determine the genes and processes enhanced in wheat due to TaNACL-D1 overexpression, both in the presence and absence of the causal agent of FHB, Fusarium graminearum (0- and 1-day post-treatment). The overexpression of TaNACL-D1 resulted in more pronounced transcriptional reprogramming as a response to fungal infection, leading to the enhanced expression of genes involved in detoxification, immune responses, secondary metabolism, hormone biosynthesis, and signalling. The regulation and response to JA and ABA were differentially regulated between the OE and the WT. Furthermore, the results suggest that the OE may more efficiently: (i) regulate the oxidative burst; (ii) modulate cell death; and (iii) induce both the phenylpropanoid pathway and lignin synthesis. Thus, this study provides insights into the mode of action and downstream target pathways for this novel NAC transcription factor, further validating its potential as a gene to enhance FHB resistance in wheat.

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.

Origin of minicircular mitochondrial genomes in red algae

Eukaryotic organelle genomes are generally of conserved size and gene content within phylogenetic groups. However, significant variation in genome structure may occur. Here, we report that the Stylonematophyceae red algae contain multipartite circular mitochondrial genomes (i.e., minicircles) which encode one or two genes bounded by a specific cassette and a conserved constant region. These minicircles are visualized using fluorescence microscope and scanning electron microscope, proving the circularity. Mitochondrial gene sets are reduced in these highly divergent mitogenomes. Newly generated chromosome-level nuclear genome assembly of Rhodosorus marinus reveals that most mitochondrial ribosomal subunit genes are transferred to the nuclear genome. Hetero-concatemers that resulted from recombination between minicircles and unique gene inventory that is responsible for mitochondrial genome stability may explain how the transition from typical mitochondrial genome to minicircles occurs. Our results offer inspiration on minicircular organelle genome formation and highlight an extreme case of mitochondrial gene inventory reduction.

When and why are mitochondria paternally inherited?

In contrast with nuclear genes that are passed on through both parents, mitochondrial genes are maternally inherited in most species, most of the time. The genetic conflict stemming from this transmission asymmetry is well-documented, and there is an abundance of population-genetic theory associated with it. While occasional or aberrant paternal inheritance occurs, there are only a few cases where exclusive paternal inheritance of mitochondrial genomes is the evolved state. Why this is remains poorly understood. By examining commonalities between species with exclusive paternal inheritance, we discuss what they may tell us about the evolutionary forces influencing mitochondrial inheritance patterns. We end by discussing recent technological advances that make exploring the causes and consequences of paternal inheritance feasible.

Rate accelerations in plastid and mitochondrial genomes of Cyperaceae occur in the same clades

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes.

WHIRLY protein functions in plants

Environmental stresses pose a significant threat to food security. Understanding the function of proteins that regulate plant responses to biotic and abiotic stresses is therefore pivotal in developing strategies for crop improvement. The WHIRLY (WHY) family of DNA-binding proteins are important in this regard because they fulfil a portfolio of important functions in organelles and nuclei. The WHY1 and WHY2 proteins function as transcription factors in the nucleus regulating phytohormone synthesis and associated growth and stress responses, as well as fulfilling crucial roles in DNA and RNA metabolism in plastids and mitochondria. WHY1, WHY2 (and WHY3 proteins in Arabidopsis) maintain organelle genome stability and serve as auxiliary factors for homologous recombination and double-strand break repair. Our understanding of WHY protein functions has greatly increased in recent years, as has our knowledge of the flexibility of their localization and overlap of functions but there is no review of the topic in the literature. Our aim in this review was therefore to provide a comprehensive overview of the topic, discussing WHY protein functions in nuclei and organelles and highlighting roles in plant development and stress responses. In particular, we consider areas of uncertainty such as the flexible localization of WHY proteins in terms of retrograde signalling connecting mitochondria, plastids, and the nucleus. Moreover, we identify WHY proteins as important targets in plant breeding programmes designed to increase stress tolerance and the sustainability of crop yield in a changing climate.

Ancient endosymbiont-mediated transmission of a selfish gene provides a model for overcoming barriers to gene transfer into animal mitochondrial genomes

In contrast to bilaterian animals, non-bilaterian mitochondrial genomes contain atypical genes, often attributed to horizontal gene transfer (HGT) as an ad hoc explanation. Although prevalent in plants, HGT into animal mitochondrial genomes is rare, lacking suitable explanatory models for their occurrence. HGT of the mismatch DNA repair gene (mtMutS) from giant viruses to octocoral (soft corals and their kin) mitochondrial genomes provides a model for how barriers to HGT to animal mitochondria may be overcome. A review of the available literature suggests that this HGT was mediated by an alveolate endosymbiont infected with a lysogenic phycodnavirus that enabled insertion of the homing endonuclease containing mtMutS into octocoral mitochondrial genomes. We posit that homing endonuclease domains and similar selfish elements play a crucial role in such inter-domain gene transfers. Understanding the role of selfish genetic elements in HGT has the potential to aid development of tools for manipulating animal mitochondrial DNA.

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.

Coordinated regulation of the mitochondrial retrograde response by circadian clock regulators and ANAC017

Mitochondrial retrograde signaling (MRS) supports photosynthetic function under a variety of conditions. Induction of mitochondrial dysfunction with myxothiazol (a specific inhibitor of the mitochondrial bc₁ complex) or antimycin A (an inhibitor of the mitochondrial bc₁ complex and cyclic electron transport in the chloroplast under light conditions) in the light and dark revealed diurnal control of MRS. This was evidenced by (1) significantly enhanced binding of ANAC017 to promoters in the light compared with the dark in Arabidopsis plants treated with myxothiazol (but not antimycin A), (2) overlap in the experimentally determined binding sites for ANAC017 and circadian clock regulators in the promoters of ANAC013 and AOX1a, (3) a diurnal expression pattern for ANAC017 and transcription factors it regulates, (4) altered expression of ANAC017-regulated genes in circadian clock mutants with and without myxothiazol treatment, and (5) a decrease in the magnitude of LHY and CCA1 expression in an ANAC017-overexpressing line and protein-protein interaction between ANAC017 and PIF4. This study also shows a large difference in transcriptome responses to antimycin A and myxothiazol in the dark: these responses are ANAC017 independent, observed in shoots and roots, similar to biotic challenge and salicylic acid responses, and involve ERF and ZAT transcription factors. This suggests that antimycin A treatment stimulates a second MRS pathway that is mediated or converges with salicylic acid signaling and provides a merging point with chloroplast retrograde signaling.

Control of plastid inheritance by environmental and genetic factors

The genomes of cytoplasmic organelles (mitochondria and plastids) are maternally inherited in most eukaryotes, thus excluding organellar genomes from the benefits of sexual reproduction and recombination. The mechanisms underlying maternal inheritance are largely unknown. Here we demonstrate that two independently acting mechanisms ensure maternal inheritance of the plastid (chloroplast) genome. Conducting large-scale genetic screens for paternal plastid transmission, we discovered that mild chilling stress during male gametogenesis leads to increased entry of paternal plastids into sperm cells and strongly increased paternal plastid transmission. We further show that the inheritance of paternal plastid genomes is controlled by the activity of a genome-degrading exonuclease during pollen maturation. Our data reveal that (1) maternal inheritance breaks down under specific environmental conditions, (2) an organelle exclusion mechanism and a genome degradation mechanism act in concert to prevent paternal transmission of plastid genes and (3) plastid inheritance is determined by complex gene-environment interactions.

Long-read sequencing characterizes mitochondrial and plastid genome variants in Arabidopsis msh1 mutants

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.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

Altered collective mitochondrial dynamics in the Arabidopsis msh1 mutant compromising organelle DNA maintenance

Mitochondria form highly dynamic populations in the cells of plants (and almost all eukaryotes). The characteristics and benefits of this collective behaviour, and how it is influenced by nuclear features, remain to be fully elucidated. Here, we use a recently developed quantitative approach to reveal and analyse the physical and collective 'social' dynamics of mitochondria in an Arabidopsis msh1 mutant where the organelle DNA maintenance machinery is compromised. We use a newly created line combining the msh1 mutant with mitochondrially targeted green fluorescent protein (GFP), and characterize mitochondrial dynamics with a combination of single-cell time-lapse microscopy, computational tracking, and network analysis. The collective physical behaviour of msh1 mitochondria is altered from that of the wild type in several ways: mitochondria become less evenly spread, and networks of inter-mitochondrial encounters become more connected, with greater potential efficiency for inter-organelle exchange-reflecting a potential compensatory mechanism for the genetic challenge to the mitochondrial DNA population, supporting more inter-organelle exchange. We find that these changes are similar to those observed in friendly, where mitochondrial dynamics are altered by a physical perturbation, suggesting that this shift to higher connectivity may reflect a general response to mitochondrial challenges, where physical dynamics of mitochondria may be altered to control the genetic structure of the mtDNA population.

Social networks in the single cell

This article comments on: Chustecki JM, Etherington RD, Gibbs DJ, Johnston IG. 2022. Altered collective mitochondrial dynamics in the Arabidopsis msh1 mutant compromising organelle DNA maintenance. Journal of Experimental Botany 73, 5428-5439. Plant mitochondrial DNA (mtDNA) can become damaged in many ways. A major repair mechanism is homologous recombination, which requires an undamaged DNA template. Presumably, this template comes from a different mitochondrion in the same cell. Plant mitochondria undergo fission and fusion to form transient networks which could allow the exchange of genetic information. To test this hypothesis, Chustecki et al. (2022) used msh1 mutants with defective DNA repair, and showed that mitochondrial interactions increased, revealing a link between the physical and genetic behavior of mitochondria.

Sorting of mitochondrial and plastid heteroplasmy in Arabidopsis is extremely rapid and depends on MSH1 activity

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.

Cytonuclear coevolution in a holoparasitic plant with highly disparate organellar genomes

Contrasting substitution rates in the organellar genomes of Lophophytum agree with the DNA repair, replication, and recombination gene content. Plastid and nuclear genes whose products form multisubunit complexes co-evolve. The organellar genomes of the holoparasitic plant Lophophytum (Balanophoraceae) show disparate evolution. In the plastid, the genome has been severely reduced and presents a > 85% AT content, while in the mitochondria most protein-coding genes have been replaced by homologs acquired by horizontal gene transfer (HGT) from their hosts (Fabaceae). Both genomes carry genes whose products form multisubunit complexes with those of nuclear genes, creating a possible hotspot of cytonuclear coevolution. In this study, we assessed the evolutionary rates of plastid, mitochondrial and nuclear genes, and their impact on cytonuclear evolution of genes involved in multisubunit complexes related to lipid biosynthesis and proteolysis in the plastid and those in charge of the oxidative phosphorylation in the mitochondria. Genes from the plastid and the mitochondria (both native and foreign) of Lophophytum showed extremely high and ordinary substitution rates, respectively. These results agree with the biased loss of plastid-targeted proteins involved in angiosperm organellar repair, replication, and recombination machinery. Consistent with the high rate of evolution of plastid genes, nuclear-encoded subunits of plastid complexes showed disproportionate increases in non-synonymous substitution rates, while those of the mitochondrial complexes did not show different rates than the control (i.e. non-organellar nuclear genes). Moreover, the increases in the nuclear-encoded subunits of plastid complexes were positively correlated with the level of physical interaction they possess with the plastid-encoded ones. Overall, these results suggest that a structurally-mediated compensatory factor may be driving plastid-nuclear coevolution in Lophophytum, and that mito-nuclear coevolution was not altered by HGT.

Temporal dynamics of mildly deleterious nonsynonymous substitutions in mitochondrial gene sequences in rodents and moles

We have previously estimated the evolutionary rate (number of substitutions/site/million years) of mitochondrial cytochrome b gene (Cytb) sequences in rodents and moles to be about 0.11 at more recent divergence times of tens of thousands of years, and to decrease rapidly to about 0.03 at more distant divergence times. Because this time dependency is thought to be caused by the removal of mildly deleterious substitutions in later generations, we focused in this study on the abundance of nonsynonymous substitutions. We collected 23 haplogroups of Cytb with signals of late Quaternary population expansion events from rodents and moles and categorized them into three groups for comparison based on predicted expansion start time: 5,000-15,000 years ago (Group I), ca. 53,000 years ago (Group II) and 130,000-230,000 years ago (Group III). We counted the numbers of nonsynonymous and synonymous substitutions in all haplogroups. The rates of nonsynonymous substitutions were lowest in Groups II and III (0.08-0.22), whereas those in Group I varied markedly. We further classified Group I into two subgroups based on high (0.29-0.43) and low (0.09-0.20) nonsynonymous substitution rates, which were likely to be associated with the start of the expansion within 10,000 years and at around 15,000 years ago, respectively. The Group II and III networks had two- or three-step star-shaped structures and tended to exhibit frequent and less frequent nonsynonymous substitutions on exterior and interior branches, respectively. Based on temporal dynamics, nonsynonymous mitochondrial DNA (mtDNA) substitutions in small mammals accounted for at most 40% of all substitutions during the early evolutionary stage and then rapidly declined, dropping to approximately 15%. The results of this study provide a good explanation of the time-dependent trend in the mtDNA evolution rate predicted in previous work.

Combining Proteomics and Metabolomics to Analyze the Effects of Spaceflight on Rice Progeny

Spaceflight is a special abiotic stress, the biological effect mechanism of which on contemporary rice has been clarified, However, its effect on offspring rice was still unclear. In order to understand the response mechanism of F2 generation plants to space flight, this study used SJ-10 recoverable satellite to carry DN423 rice seeds for 12.5 days in orbit flight. After returning to the ground, the plants were then planted to F2 generation to explore the biological effect mechanism. Our research showed that in the F2 generation of TLS, the rice plant height of the space flight group increased by 33.8%, the ear length and thousand-grain weight decreased by 9.7 and 4.6%, respectively, and the grain number per panicle increased by 6.5%. Moreover, related proteins that control changes in agronomic traits have been identified. The changes of MDA, H2O2, soluble sugar, electron leakage and antioxidant enzyme activity confirmed the stress response in F2 generation plants. ITRAQ and LC-MS technology were used to reveal the change pattern of protein levels and metabolite levels in F2 generation plants, 389 and 405 proteins were identified as differentially abundant proteins in TLS and TS, respectively. In addition, there were 124 and 125 metabolites that changed during these two periods. The proteome and metabolome result further confirmed that the F2 generation plants still retained the memory of space flight stress, and retained the memory of space flight stress through genome instability. Oxidative stress signals activated sugar signals to rebuild metabolic networks to adapt to space flight stress. The reconstruction of energy metabolism, amino acid metabolism, phenylalanine metabolism, and flavonoid metabolism played an important role in the process of adapting to space flight stress. The results of this study broaden the perspective of space biological effects and provide a basis for studying the effects of abiotic stress on plant progeny.

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.

Complete sequence of a 641-kb insertion of mitochondrial DNA in the Arabidopsis thaliana nuclear genome

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 Complicated Nature of Somatic mtDNA Mutations in Aging

Mitochondria are the main source of energy used to maintain cellular homeostasis. This aspect of mitochondrial biology underlies their putative role in age-associated tissue dysfunction. Proper functioning of the electron transport chain (ETC), which is partially encoded by the extra-nuclear mitochondrial genome (mtDNA), is key to maintaining this energy production. The acquisition of de novo somatic mutations that interrupt the function of the ETC have long been associated with aging and common diseases of the elderly. Yet, despite over 30 years of study, the exact role(s) mtDNA mutations play in driving aging and its associated pathologies remains under considerable debate. Furthermore, even fundamental aspects of age-related mtDNA mutagenesis, such as when mutations arise during aging, where and how often they occur across tissues, and the specific mechanisms that give rise to them, remain poorly understood. In this review, we address the current understanding of the somatic mtDNA mutations, with an emphasis of when, where, and how these mutations arise during aging. Additionally, we highlight current limitations in our knowledge and critically evaluate the controversies stemming from these limitations. Lastly, we highlight new and emerging technologies that offer potential ways forward in increasing our understanding of somatic mtDNA mutagenesis in the aging process.

The Barley Chloroplast Mutator (cpm) Mutant: All Roads Lead to the Msh1 Gene

The barley chloroplast mutator (cpm) is a nuclear gene mutant that induces a wide spectrum of cytoplasmically inherited chlorophyll deficiencies. Plastome instability of cpm seedlings was determined by identification of a particular landscape of polymorphisms that suggests failures in a plastome mismatch repair (MMR) protein. In Arabidopsis, MSH genes encode proteins that are in charge of mismatch repair and have anti-recombination activity. In this work, barley homologs of these genes were identified, and their sequences were analyzed in control and cpm mutant seedlings. A substitution, leading to a premature stop codon and a truncated MSH1 protein, was identified in the Msh1 gene of cpm plants. The relationship between this mutation and the presence of chlorophyll deficiencies was established in progenies from crosses and backcrosses. These results strongly suggest that the mutation identified in the Msh1 gene of the cpm mutant is responsible for the observed plastome instabilities. Interestingly, comparison of mutant phenotypes and molecular changes induced by the barley cpm mutant with those of Arabidopsis MSH1 mutants revealed marked differences.

The Development of Mitochondrial Gene Editing Tools and Their Possible Roles in Crop Improvement for Future Agriculture

We are living in the era of genome editing. Nowadays, targeted editing of the plant nuclear DNA is prevalent in basic biological research and crop improvement since its first establishment a decade ago. However, achieving the same accomplishment for the plant mitochondrial genome has long been deemed impossible. Recently, the pioneer studies on editing plant mitogenome have been done using the mitochondria-targeted transcription activator-like effector nucleases (mitoTALENs) in rice, rapeseed, and Arabidopsis. It is well documented that mitochondria play essential roles in plant development and stress tolerance, particularly, in cytoplasmic male sterility widely used in production of hybrids. The success of mitochondrial genome editing enables studying the fundamentals of mitochondrial genome. Furthermore, mitochondrial RNA editing (mostly by nuclear-encoded pentatricopeptide repeat (PPR) proteins) in a sequence-specific manner can simultaneously change the production of translatable mitochondrial mRNA. Moreover, direct editing of the nuclear-encoding mitochondria-targeted factors required for plant mitochondrial genome dynamics and recombination may facilitate genetic manipulation of plant mitochondria. Here, the present state of knowledge on editing the plant mitochondrial genome is reviewed.

The minicircular and extremely heteroplasmic mitogenome of the holoparasitic plant Rhopalocnemis phalloides

The plastid and nuclear genomes of parasitic plants exhibit deeply altered architectures,^1-13 whereas the few examined mitogenomes range from deeply altered to conventional.^14-20 To provide further insight on mitogenome evolution in parasitic plants, we report the highly modified mitogenome of Rhopalocnemis phalloides, a holoparasite in Balanophoraceae. Its mitogenome is uniquely arranged in 21 minicircular chromosomes that vary in size from 4,949 to 7,861 bp, with a total length of only 130,713 bp. All chromosomes share an identical 896 bp conserved region, with a large stem-loop that acts as the origin of replication, flanked on each side by hypervariable and semi-conserved regions. Similar minicircular structures with shared and unique regions have been observed in parasitic animals and free-living protists,^21-24 suggesting convergent structural evolution. Southern blots confirm both the minicircular structure and the replication origin of the mitochondrial chromosomes. PacBio reads provide evidence for chromosome recombination and rolling-circle replication for the R. phalloides mitogenome. Despite its small size, the mitogenome harbors a typical set of genes and introns within the unique regions of each chromosome, yet introns are the smallest among seed plants and ferns. The mitogenome also exhibits extreme heteroplasmy, predominantly involving short indels and more complex variants, many of which cause potential loss-of-function mutations for some gene copies. All heteroplasmic variants are transcribed, and functional and nonfunctional protein-coding variants are spliced and RNA edited. Our findings offer a unique perspective into how mitogenomes of parasitic plants can be deeply altered and shed light on plant mitogenome replication.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

Plant mitochondria - past, present and future

The study of plant mitochondria started in earnest around 1950 with the first isolations of mitochondria from animal and plant tissues. The first 35 years were spent establishing the basic properties of plant mitochondria and plant respiration using biochemical and physiological approaches. A number of unique properties (compared to mammalian mitochondria) were observed: (i) the ability to oxidize malate, glycine and cytosolic NAD(P)H at high rates; (ii) the partial insensitivity to rotenone, which turned out to be due to the presence of a second NADH dehydrogenase on the inner surface of the inner mitochondrial membrane in addition to the classical Complex I NADH dehydrogenase; and (iii) the partial insensitivity to cyanide, which turned out to be due to an alternative oxidase, which is also located on the inner surface of the inner mitochondrial membrane, in addition to the classical Complex IV, cytochrome oxidase. With the appearance of molecular biology methods around 1985, followed by genomics, further unique properties were discovered: (iv) plant mitochondrial DNA (mtDNA) is 10-600 times larger than the mammalian mtDNA, yet it only contains approximately 50% more genes; (v) plant mtDNA has kept the standard genetic code, and it has a low divergence rate with respect to point mutations, but a high recombinatorial activity; (vi) mitochondrial mRNA maturation includes a uniquely complex set of activities for processing, splicing and editing (at hundreds of sites); (vii) recombination in mtDNA creates novel reading frames that can produce male sterility; and (viii) plant mitochondria have a large proteome with 2000-3000 different proteins containing many unique proteins such as 200-300 pentatricopeptide repeat proteins. We describe the present and fairly detailed picture of the structure and function of plant mitochondria and how the unique properties make their metabolism more flexible allowing them to be involved in many diverse processes in the plant cell, such as photosynthesis, photorespiration, CAM and C4 metabolism, heat production, temperature control, stress resistance mechanisms, programmed cell death and genomic evolution. However, it is still a challenge to understand how the regulation of metabolism and mtDNA expression works at the cellular level and how retrograde signaling from the mitochondria coordinates all those processes.

Mitochondrial mutations in Caenorhabditis elegans show signatures of oxidative damage and an AT-bias

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.

Assessment of Mitochondrial DNA Copy Number, Stability, and Repair in Arabidopsis

Mitochondrial functions depend on the proper maintenance and expression of the mitochondrial genome (mtDNA). Therefore, understanding mtDNA replication and repair requires methods to assess its integrity. Mutations or chemical treatments that affect processes involved in the maintenance or stability of the mtDNA can affect its global copy number, but also the relative abundance of different genomic regions or the frequency of illegitimate recombination across repeated sequences. These can be conveniently tested by quantitative PCR (qPCR). Arabidopsis thaliana offers several advantages for studying these processes, because of the extensive collections of mutants, natural accessions and other genetic resources available from stock centers. Here we describe protocols we routinely use to explore changes in mtDNA copy number and relative stoichiometry in Arabidopsis mutants of genes involved in the replication, repair and recombination of the mtDNA.

Multiple Losses of MSH1, Gain of mtMutS, and Other Changes in the MutS Family of DNA Repair Proteins in Animals

MutS is a key component of the mismatch repair (MMR) pathway. Members of the MutS protein family are present in prokaryotes, eukaryotes, and viruses. Six MutS homologs (MSH1-6) have been identified in yeast, of which three function in nuclear MMR, while MSH1 functions in mitochondrial DNA repair. MSH proteins are believed to be well conserved in animals, except for MSH1-which is thought to be lost. Two intriguing exceptions to this general picture have been found, both in the class Anthozoa within the phylum Cnidaria. First, an ortholog of the yeast-MSH1 was reported in one hexacoral species. Second, a MutS homolog (mtMutS) has been found in the mitochondrial genome of all octocorals. To understand the origin and potential functional implications of these exceptions, we investigated the evolution of the MutS family both in Cnidaria and in animals in general. Our study confirmed the acquisition of octocoral mtMutS by horizontal gene transfer from a giant virus. Surprisingly, we identified MSH1 in all hexacorals and several sponges and placozoans. By contrast, MSH1 orthologs were lacking in other cnidarians, ctenophores, and bilaterian animals. Furthermore, while we identified MSH2 and MSH6 in nearly all animals, MSH4, MSH5, and, especially, MSH3 were missing in multiple species. Overall, our analysis revealed a dynamic evolution of the MutS family in animals, with multiple losses of MSH1, MSH3, some losses of MSH4 and MSH5, and a gain of the octocoral mtMutS. We propose that octocoral mtMutS functionally replaced MSH1 that was present in the common ancestor of Anthozoa.

An Insight Into the Mechanism of Plant Organelle Genome Maintenance and Implications of Organelle Genome in Crop Improvement: An Update

Besides the nuclear genome, plants possess two small extra chromosomal genomes in mitochondria and chloroplast, respectively, which contribute a small fraction of the organelles’ proteome. Both mitochondrial and chloroplast DNA have originated endosymbiotically and most of their prokaryotic genes were either lost or transferred to the nuclear genome through endosymbiotic gene transfer during the course of evolution. Due to their immobile nature, plant nuclear and organellar genomes face continuous threat from diverse exogenous agents as well as some reactive by-products or intermediates released from various endogenous metabolic pathways. These factors eventually affect the overall plant growth and development and finally productivity. The detailed mechanism of DNA damage response and repair following accumulation of various forms of DNA lesions, including single and double-strand breaks (SSBs and DSBs) have been well documented for the nuclear genome and now it has been extended to the organelles also. Recently, it has been shown that both mitochondria and chloroplast possess a counterpart of most of the nuclear DNA damage repair pathways and share remarkable similarities with different damage repair proteins present in the nucleus. Among various repair pathways, homologous recombination (HR) is crucial for the repair as well as the evolution of organellar genomes. Along with the repair pathways, various other factors, such as the MSH1 and WHIRLY family proteins, WHY1, WHY2, and WHY3 are also known to be involved in maintaining low mutation rates and structural integrity of mitochondrial and chloroplast genome. SOG1, the central regulator in DNA damage response in plants, has also been found to mediate endoreduplication and cell-cycle progression through chloroplast to nucleus retrograde signaling in response to chloroplast genome instability. Various proteins associated with the maintenance of genome stability are targeted to both nuclear and organellar compartments, establishing communication between organelles as well as organelles and nucleus. Therefore, understanding the mechanism of DNA damage repair and inter compartmental crosstalk mechanism in various sub-cellular organelles following induction of DNA damage and identification of key components of such signaling cascades may eventually be translated into strategies for crop improvement under abiotic and genotoxic stress conditions. This review mainly highlights the current understanding as well as the importance of different aspects of organelle genome maintenance mechanisms in higher plants.

Development and Evaluation of a Droplet Digital PCR Assay for 8p23 β-Defensin Cluster Copy Number Determination

Introduction

It has been proposed that the copy number (CN) variation (CNV) in β-defensin genes (DEFB) on human chromosome 8p23 determines phenotypic differences in inflammatory diseases. However, no method for accurate and easy DEFB CN quantification is yet available.

Objective

Droplet digital polymerase chain reaction (ddPCR) is a novel method for CNV quantification and has been used for genes such as CCL4L, CCL3L1, AMY1, and HER2. However, to date, no ddPCR assay has been available for DEFB CN determination. In the present study, we aimed to develop and evaluate such a ddPCR assay.

Methods

The assay was designed using DEFB4 and RPP30 as the target and the reference gene, respectively. To evaluate the assay, 283 DNA samples with known CNs previously determined using the multiple ligation-dependent probe amplification (MLPA) method, the current gold standard, were used as standards. To discover the optimal DNA template amount, we tested 80 to 2.5 ng DNA by a serial of one to two dilutions of eight samples. To evaluate the reproducibility of the assay, 31 samples were repeated to calculate the intra- and inter-assay variations. To further validate the reliability of the assay, the CNs of all 283 samples were determined using ddPCR. To compare results with those using quantitative PCR (qPCR), DEFB CNs for 48 samples were determined using qPCR with the same primers and probes.

Results

In a one-dimensional plot, the positive and negative droplets were clearly separated in both DEFB4 and RPP30 detection channels. In a two-dimensional plot, four populations of droplets were observed. The 20 ng template DNA proved optimal, with either high (80 ng) or low (10, 5, 2.5 ng) template input leading to ambiguous or inaccurate results. For the 31 standard samples, DEFB CNs were accurately determined with small intra- and inter-assay variations (coefficient of variation < 0.04 for both). In the validation cohort, ddPCR provided the correct CN for all 283 samples with high confidence. qPCR measurements for the 48 samples produced noisy data with high uncertainty and low accuracy.

Conclusions

ddPCR is an accurate, reproducible, easy-to-use, cheap, high-throughput method for DEFB CN determination. ddPCR could be applied to DEFB CN quantification in large-scale case–control studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40291-021-00546-2.

Detecting de novo mitochondrial mutations in angiosperms with highly divergent evolutionary rates

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.

Mitochondrial Short-Term Plastic Responses and Long-Term Evolutionary Dynamics in Animal Species

How do species respond or adapt to environmental changes? The answer to this depends partly on mitochondrial epigenetics and genetics, new players in promoting adaptation to both short- and long-term environmental changes. In this review, we explore how mitochondrial epigenetics and genetics mechanisms, such as mtDNA methylation, mtDNA-derived noncoding RNAs, micropeptides, mtDNA mutations, and adaptations, can contribute to animal plasticity and adaptation. We also briefly discuss the challenges in assessing mtDNA adaptive evolution. In sum, this review covers new advances in the field of mitochondrial genomics, many of which are still controversial, and discusses processes still somewhat obscure, and some of which are still quite speculative and require further robust experimentation.

Plant Mitochondria are a Riddle Wrapped in a Mystery Inside an Enigma

A fundamental paradox motivates the study of plant mitochondrial genomics: the mutation rate is very low (lower than in the nucleus) but the rearrangement rate is high. A landmark paper published in Journal of Molecular Evolution in 1988 established these facts and revealed the paradox. Jeffrey Palmer and Laura Herbon did a prodigious amount of work in the pre-genome sequencing era to identify both the high frequency of rearrangements between closely related species, and the low frequency of mutations, observations that have now been confirmed many times by sequencing. This paper was also the first to use molecular data on rearrangements as a phylogenetic trait to build a parsimonious tree. The work was a technical tour-de-force, its findings are still at the heart of plant mitochondrial genomics, and the underlying molecular mechanisms that produce this paradox are still not completely understood.

Mutational Pressure Drives Differential Genome Conservation in Two Bacterial Endosymbionts of Sap-Feeding Insects

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.