DNA Repair and the Stability of the Plant Mitochondrial Genome

De novo assembly of the mitochondrial genome of Glycyrrhiza glabra and identification of two types of homologous recombination configurations caused by repeat sequences

BACKGROUND

Glycyrrhiza glabra, which is widely used in medicine and therapy, is known as the 'king of traditional Chinese medicine'. In this study, we successfully assembled and annotated the mitochondrial and chloroplast genomes of G. glabra via high-throughput sequencing technology, combining the advantages of short-read (Illumina) and long-read (Oxford Nanopore) sequencing.

RESULTS

We revealed the ring structure of the mitochondrial genome, which spans 421,293 bp with 45.1% GC content and 56 annotated genes. Notably, we identified 514 repetitive sequences, including 123 Simple sequence repeats (SRs), 3 Tndem sequence repeats (TSRs), and 388 Dispersed sequence repeats (DSRs). We identified 79 out of the 388 DSRs as potentially involved in homologous recombination. We identified five forward repeats and four palindromic repeats that facilitate homologous recombination and induce alterations in the mitochondrial genome structure. We corroborated this finding via polymerase chain reaction (PCR). Furthermore, we identified chloroplast-derived sequence fragments within the mitochondrial genome, offering novel insights into the evolutionary history of plant mitochondrial genomes. We predicted 460 potential RNA editing sites, primarily involving cytosine-to-uracil transitions. This study reveals the complexity of repetitive sequence-mediated homologous recombination in the mitochondrial genome of G. glabra and provides new insights into its structure, function, and evolution.

CONCLUSIONS

These findings have important implications for conservation biology, population genetics, and evolutionary studies, underscoring the role of repetitive sequences in genome dynamics and highlighting the need for further research on mitochondrial genome evolution and function in plants.

The evolutionarily diverged single-stranded DNA-binding proteins SSB1/SSB2 differentially affect the replication, recombination and mutation of organellar genomes in Arabidopsis thaliana

Single-stranded DNA-binding proteins (SSBs) play essential roles in the replication, recombination and repair processes of organellar DNA molecules. In Arabidopsis thaliana, SSBs are encoded by a small family of two genes (SSB1 and SSB2). However, the functional divergence of these two SSB copies in plants remains largely unknown, and detailed studies regarding their roles in the replication and recombination of organellar genomes are still incomplete. In this study, phylogenetic, gene structure and protein motif analyses all suggested that SSB1 and SSB2 probably diverged during the early evolution of seed plants. Based on accurate long-read sequencing results, ssb1 and ssb2 mutants had decreased copy numbers for both mitochondrial DNA (mtDNA) and plastid DNA (ptDNA), accompanied by a slight increase in structural rearrangements mediated by intermediate-sized repeats in mt genome and small-scale variants in both genomes. Our findings provide an important foundation for further investigating the effects of DNA dosage in the regulation of mutation frequencies in plant organellar genomes.

Assembly and analysis of the complete mitochondrial genome of Carya illinoinensis to provide insights into the conserved sequences of tRNA genes

Carya illinoinensis is an economically important nut tree, and its chloroplast (cp.) genome has been reported; however, its mitochondrial (mt) genome remains unknown. In the present study, we assembled the first mt genome of C. illinoinensis. The circular mt genome of C. illinoinensis is 495,205 bp long, with 37 protein-coding genes(PCGs), 24 tRNA genes, and 3 rRNA genes. All the tRNAs could be folded into typical cloverleaf secondary structures, with lengths of 58-88 bp. A conserved U-U-C-x-A-x2 consensus nucleotide sequence was discovered in the Ψ-loops of tRNA sequences. In addition, 447 dispersed repeats were detected, as well as found 482 RNA editing sites and 9,960 codons in the mt genome. Furthermore, a total of 27 DNA sequences with a length of 43,277 bp were transferred from the cp. to the mt genome, and eight integrated cp-derived genes (trnL-CAA, trnV-GAC, trnD-GUC, trnW-CCA, trnN-GUU, trnH-GUG, trnM-CAU, and rps7) were identified. We also obtained 1,086 hits, including 364.023 kp of nuclear genome sequences, that were transferred to the mt genome. To determine the evolutionary position of C. illinoinensis, we conducted a phylogenetic analysis of the mitogenomes of C. illinoinensis and 14 other taxa. The results strongly suggested that C. illinoinensis and Fagus sylvatica formed a single clade with 100% bootstrap support. This study sequenced comprehensive data on the C. illinoinensis mitochondrial genome and provided insights into the conserved sequences of tRNA genes, which could facilitate evolutionary research in other Carya trees in the future.

The complete mitochondrial genome assembly of Capsicum pubescens reveals key evolutionary characteristics of mitochondrial genes of two Capsicum subspecies

BACKGROUND

Pepper (Capsicum pubescens), one of five domesticated pepper species, has unique characteristics, such as numerous hairs on the epidermis of its leaves and stems, black seeds, and vibrant purple flowers. To date, no studies have reported on the complete assembly of the mitochondrial genome (mitogenome) of C. pubescens. Understanding the mitogenome is crucial for further research on C. pubescens.

RESULTS

In our study, we successfully assembled the first mitogenome of C. pubescens, which was assigned the GenBank accession number OP957066. This mitogenome has a length of 454,165 bp and exhibits the typical circular structure observed in most mitogenomes. We annotated a total of 70 genes, including 35 protein-coding genes (PCGs), 30 tRNA genes, 3 rRNA genes, and 2 pseudogenes. Compared to the other three pepper mitogenomes (KJ865409, KJ865410, and MN196478), C. pubescens OP957066 exhibited four unique PCGs (atp4, atp8, mttB, and rps1), while two PCGs (rpl10 and rps3) were absent. Notably, each of the three pepper mitogenomes from C. annuum (KJ865409, KJ865410, and MN196478) experienced the loss of four PCGs (atp4, atp8, mttB, and rps1). To further explore the evolutionary relationships, we reconstructed a phylogenetic tree using the mitogenomes of C. pubescens and fourteen other species. Structural comparison and synteny analysis of the above four pepper mitogenomes revealed that C. pubescens shares high sequence similarity with KJ865409 and that C. pubescens has rearranged with the other three pepper mitogenomes. Interestingly, we observed 72 similar sequences between the mitochondrial and chloroplast genomes, which accounted for 12.60% of the mitogenome, with a total length of 57,207 bp. These sequences encompassed 12 tRNA genes and the rRNA gene (rrn18). Remarkably, selective pressure analysis suggested that the nad5 gene underwent obvious positive selection. Furthermore, a single-base mutation in three genes (nad1, nad2, and nad4) resulted in an amino acid change.

CONCLUSION

This study provides a high-quality mitogenome of pepper, providing valuable molecular data for future investigations into the exchange of genetic information between pepper organelle genomes.

Mitochondrial genome insights into the spatio-temporal distribution and genetic diversity of Dendrobium hancockii Rolfe (Orchidaceae)

Introduction

With its distinctive evolutionary rate and inheritance patterns separate from the nuclear genome, mitochondrial genome analysis has become a prominent focus of current research. Dendrobium hancockii Rolfe, a species of orchid with both medicinal and horticultural value, will benefit from the application of the fully assembled and annotated mitochondrial genome. This will aid in elucidating its phylogenetic relationships, comparative genomics, and population genetic diversity.

Methods

Based on sequencing results from Illumina combined with PacBio and Nanopore, the mitochondrial genome map of D. hancockii was constructed. Comparative analysis was conducted from the perspectives of phylogeny across multiple species, selection pressure on protein-coding genes, and homologous segments. The population diversity of D. hancockii was analyzed using single nucleotide polymorphism (SNP) data from the mitochondrial genome and single-copy nuclear genes.

Results and discussion

This research constructed a circular mitochondrial map for D. hancockii, spanning 523,952 bp, containing 40 unique protein-coding genes, 37 transfer RNA genes, and 4 ribosomal RNA genes. Comparative analysis of mitochondrial genes from 26 land plants revealed a conserved gene cluster, "rpl16-ccmFn-rps3-rps19," particularly within the Dendrobium genus. The mitochondrial genome of D. hancockii exhibits a lower point mutation rate but significant structural variation. Analysis of 103 resequencing samples identified 19,101 SNP sites, dividing D. hancockii into two major groups with limited gene flow between them, as supported by population diversity, genetic structure analysis, principal component analysis, and phylogenetic trees. The geographical distribution and genetic differentiation of D. hancockii into two major groups suggest a clear phytogeographical division, likely driven by ancient geological or climatic events. The close alignment of mitochondrial data with nuclear gene data highlights the potential of the mitochondrial genome for future studies on genetic evolution in this species.

Repeat-mediated recombination results in Complex DNA structure of the mitochondrial genome of Trachelospermum jasminoides

BACKGROUND

Trachelospermum jasminoides has medicinal and ornamental value and is widely distributed in China. Although the chloroplast genome has been documented, the mitochondrial genome has not yet been studied.

RESULTS

The mitochondrial genome of T. jasminoides was assembled and functionally annotated using Illumina and nanopore reads. The mitochondrial genome comprises a master circular molecular structure of 605,764 bp and encodes 65 genes: 39 protein-coding genes, 23 transfer RNA (tRNA) genes and 3 ribosomal RNA genes. In addition to the single circular conformation, we found many alternative conformations of the T. jasminoides mitochondrial genome mediated by 42 repetitive sequences. Six repetitive sequences (DRS01-DRS06) were supported by nanopore long reads, polymerase chain reaction (PCR) amplifications, and Sanger sequencing of the PCR products. Eleven homologous fragments were identified by comparing the mitochondrial and chloroplast genome sequences, including three complete tRNA genes. Moreover, 531 edited RNA sites were identified in the protein-coding sequences based on RNA sequencing data, with nad4 having the highest number of sites (54).

CONCLUSION

To our knowledge, this is the first description of the mitochondrial genome of T. jasminoides. Our results demonstrate the existence of multiple conformations. These findings lay a foundation for understanding the genetics and evolutionary dynamics of Apocynaceae.

Methods for assembling complex mitochondrial genomes in land plants

The large size and complex structural rearrangements inherent in the mitochondrial genomes of land plants pose challenges for their sequencing. Originally, the assembly of these genomes required the cloning of mitochondrial DNA fragments followed by Sanger sequencing. Subsequently, the advent of next-generation sequencing significantly expedited the process. This review highlights examples of plant mitochondrial genome assembly employing various technologies, including 454 sequencing, Illumina short sequencing reads, and Pacific Biosciences or Oxford Nanopore Technology long sequencing reads. The combination of short and long reads in hybrid assembly has proven to be the most efficient approach for achieving reliable assemblies of land plant mitochondrial genomes.

Horizontally transferred mitochondrial DNA tracts become circular by microhomology-mediated repair pathways

The holoparasitic plant Lophophytum mirabile exhibits remarkable levels of mitochondrial horizontal gene transfer (HGT). Gathering comparative data from other individuals and host plants can provide insights into the HGT process. We sequenced the mitochondrial genome (mtDNA) from individuals of two species of Lophophytum and from mimosoid hosts. We applied a stringent phylogenomic approach to elucidate the origin of the whole mtDNAs, estimate the timing of the transfers, and understand the molecular mechanisms involved. Ancestral and recent HGT events replaced and enlarged the multichromosomal mtDNA of Lophophytum spp., with the foreign DNA ascending to 74%. A total of 14 foreign mitochondrial chromosomes originated from continuous regions in the host mtDNA flanked by short direct repeats. These foreign tracts are circularized by microhomology-mediated repair pathways and replicate independently until they are lost or they eventually recombine with other chromosomes. The foreign noncoding chromosomes are variably present in the population and likely evolve by genetic drift. We present the 'circle-mediated HGT' model in which foreign mitochondrial DNA tracts become circular and are maintained as plasmid-like molecules. This model challenges the conventional belief that foreign DNA must be integrated into the recipient genome for successful HGT.

Screening clusters of charged residues in plants' mitochondrial proteins and biological significance

Protein function is dependent on charge interactions and charge biased regions, which are involved in a wide range of cellular and biochemical processes. We report the development of a new algorithm implemented in Python and its use to identify charge clusters CC (NegativeCC: NCC, PositiveCC: PCC and MixedCC: MCC) and compare their presence in mitochondrial proteins of plant groups. To characterize the resulting CC, statistical, structural and functional analyses were conducted. The screening of 105 399 protein sequences showed that 2.6 %, 0.48 % and 0.03 % of the proteins contain NCC, PCC and MCC, respectively. Mitochondrial proteins encoded by the nuclear genome of green algae have the biggest proportion of both PCC (1.6 %) and MCC (0.4 %) and mitochondrial proteins coded by the nuclear genome of other plants group have the highest portion of NCC (7.5 %). The mapping of the identified CC showed that they are mainly located in the terminal regions of the protein. Annotation showed that proteins with CC are classified as binding proteins, are included in the transmembrane transport processes, and are mainly located in the membrane. The CC scanning revealed the presence of 2373 and 784 sites and 192 and 149 motif profiles within NCC and PCC, respectively. The investigation of CC within pentatricopeptide repeat-containing proteins revealed that they are involved in correct and specific RNA editing. CC were proven to play a key role in providing insightful structural and functional information of complex protein assemblies which could be useful in biotechnological applications.

Structural analysis of the mitochondrial genome of Santalum album reveals a complex branched configuration

BACKGROUND

Santalum album L. is an evergreen tree which is mainly distributes throughout tropical and temperate regions. And it has a great medicinal and economic value.

RESULTS

In this study, the complete mitochondrial genome of S. album were assembled and annotated, which could be descried by a complex branched structure consisting of three contigs. The lengths of these three contigs are 165,122 bp, 93,430 bp and 92,491 bp. We annotated 34 genes coding for proteins (PCGs), 26 tRNA genes, and 4 rRNA genes. The analysis of repeated elements shows that there are 89 SSRs and 242 pairs of dispersed repeats in S. album mitochondrial genome. Also we found 20 MTPTs among the chloroplast and mitochondria. The 20 MTPTs sequences span a combined length of 22,353 bp, making up 15.52 % of the plastome, 6.37 % of the mitochondrial genome. Additionally, by using the Deepred-mt tool, we found 628 RNA editing sites in 34 PCGs. Moreover, significant genomic rearrangement is observed between S. album and its associated mitochondrial genomes. Finally, based on mitochondrial genome PCGs, we deduced the phylogenetic ties between S. album and other angiosperms.

CONCLUSIONS

We reported the mitochondrial genome from Santalales for the first time, which provides a crucial genetic resource for our study of the evolution of mitochondrial genome.

Evolution and maintenance of mtDNA gene content across eukaryotes

Across eukaryotes, most genes required for mitochondrial function have been transferred to, or otherwise acquired by, the nucleus. Encoding genes in the nucleus has many advantages. So why do mitochondria retain any genes at all? Why does the set of mtDNA genes vary so much across different species? And how do species maintain functionality in the mtDNA genes they do retain? In this review, we will discuss some possible answers to these questions, attempting a broad perspective across eukaryotes. We hope to cover some interesting features which may be less familiar from the perspective of particular species, including the ubiquity of recombination outside bilaterian animals, encrypted chainmail-like mtDNA, single genes split over multiple mtDNA chromosomes, triparental inheritance, gene transfer by grafting, gain of mtDNA recombination factors, social networks of mitochondria, and the role of mtDNA dysfunction in feeding the world. We will discuss a unifying picture where organismal ecology and gene-specific features together influence whether organism X retains mtDNA gene Y, and where ecology and development together determine which strategies, importantly including recombination, are used to maintain the mtDNA genes that are retained.

Accumulation of Large Lineage-Specific Repeats Coincides with Sequence Acceleration and Structural Rearrangement in Plantago Plastomes

Repeats can mediate rearrangements and recombination in plant mitochondrial genomes and plastid genomes. While repeat accumulations are linked to heightened evolutionary rates and complex structures in specific lineages, debates persist regarding the extent of their influence on sequence and structural evolution. In this study, 75 Plantago plastomes were analyzed to investigate the relationships between repeats, nucleotide substitution rates, and structural variations. Extensive repeat accumulations were associated with significant rearrangements and inversions in the large inverted repeats (IRs), suggesting that repeats contribute to rearrangement hotspots. Repeats caused infrequent recombination that potentially led to substoichiometric shifting, supported by long-read sequencing. Repeats were implicated in elevating evolutionary rates by facilitating localized hypermutation, likely through DNA damage and repair processes. This study also observed a decrease in nucleotide substitution rates for loci translocating into IRs, supporting the role of biased gene conversion in maintaining lower substitution rates. Combined with known parallel changes in mitogenomes, it is proposed that potential dysfunction in nuclear-encoded genes associated with DNA replication, recombination, and repair may drive the evolution of Plantago organellar genomes. These findings contribute to understanding how repeats impact organellar evolution and stability, particularly in rapidly evolving plant lineages.

Assembly and evolutionary analysis of the complete mitochondrial genome of Trichosanthes kirilowii, a traditional Chinese medicinal plant

Trichosanthes kirilowii (T. kirilowii) is a valuable plant used for both medicinal and edible purposes. It belongs to the Cucurbitaceae family. However, its phylogenetic position and relatives have been difficult to accurately determine due to the lack of mitochondrial genomic information. This limitation has been an obstacle to the potential applications of T. kirilowii in various fields. To address this issue, Illumina and Nanopore HiFi sequencing were used to assemble the mitogenome of T. kirilowii into two circular molecules with sizes of 245,700 bp and 107,049 bp, forming a unique multi-branched structure. The mitogenome contains 61 genes, including 38 protein-coding genes (PCGs), 20 tRNAs, and three rRNAs. Within the 38 PCGs of the T. kirilowii mitochondrial genome, 518 potential RNA editing sites were identified. The study also revealed the presence of 15 homologous fragments that span both the chloroplast and mitochondrial genomes. The phylogenetic analysis strongly supports that T. kirilowii belongs to the Cucurbitaceae family and is closely related to Luffa. Collinearity analysis of five Cucurbitaceae mitogenomes shows a high degree of structural variability. Interestingly, four genes, namely atp1, ccmFC, ccmFN, and matR, played significant roles in the evolution of T. kirilowii through selection pressure analysis. The comparative analysis of the T. kirilowii mitogenome not only sheds light on its functional and structural features but also provides essential information for genetic studies of the genus of Cucurbitaceae.

The Complete Chloroplast Genome Sequence of the Medicinal Moss Rhodobryum giganteum (Bryaceae, Bryophyta): Comparative Genomics and Phylogenetic Analyses

Rhodobryum giganteum (Bryaceae, Bryophyta), a rare medicinal bryophyte, is valued for its cardiovascular therapeutic properties in traditional Chinese medicine. This study presents the first complete chloroplast genome sequence of R. giganteum, including its assembly and annotation. The circular chloroplast genome of R. giganteum is 124,315 bp in length, displaying a typical quadripartite structure with 128 genes: 83 protein-coding genes, 37 tRNAs, and 8 rRNAs. Analyses of codon usage bias, repetitive sequences, and simple sequence repeats (SSRs) revealed an A/U-ending codon preference, 96 repetitive sequences, and 385 SSRs in the R. giganteum chloroplast genome. Nucleotide diversity analysis identified 10 high mutational hotspots. Ka/Ks ratio analysis suggested potential positive selection in rpl20, rps18, petG, and psbM genes. Phylogenetic analysis of whole chloroplast genomes from 38 moss species positioned R. giganteum within Bryales, closely related to Rhodobryum laxelimbatum. This study augments the chloroplast genomic data for Bryales and provides a foundation for molecular marker development and genetic diversity analyses in medicinal bryophytes.

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). 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.

The complete mitochondrial genome of Pontederia crassipes: using HiFi reads to investigate genome recombination and gene transfer from chloroplast genome

Water hyacinth (Pontederia crassipes Mart.) is a monocotyledonous aquatic plant renowned for its rapid growth, extensive proliferation, biological invasiveness, and ecological resilience to variations in pH, nutrients, and temperature. The International Union for Conservation of Nature (IUCN) has listed P. crassipes among the top 100 invasive species. However, comprehensive genomic information, particularly concerning its mitochondrial genome (mitogenome), remains surprisingly limited. In this study, the complete mitogenome of P. crassipes was analyzed using bioinformatics approaches. The mitogenome is 399,263 bp long and contains 38 protein-coding genes (PCGs), 24 tRNA genes, and 3 rRNA genes. Sequence analysis revealed that the complete mitogenome of the species contains 3,289 dispersed repeats, and 765 RNA editing sites in protein-coding genes. The P. crassipes mitogenome possessed un-conserved structures, including extensive sequence transfer between its chloroplasts and mitochondria. Our study on the mitogenome of P. crassipes offers critical insights into its evolutionary patterns and phylogenetic relationships with related taxa. This research enhances our understanding of this invasive species, known for its significant biomass and rapid overgrowth in aquatic environments.

Bacterial cell volume regulation and the importance of cyclic di-AMP

SUMMARYNucleotide-derived second messengers are present in all domains of life. In prokaryotes, most of their functionality is associated with general lifestyle and metabolic adaptations, often in response to environmental fluctuations of physical parameters. In the last two decades, cyclic di-AMP has emerged as an important signaling nucleotide in many prokaryotic lineages, including Firmicutes, Actinobacteria, and Cyanobacteria. Its importance is highlighted by the fact that both the lack and overproduction of cyclic di-AMP affect viability of prokaryotes that utilize cyclic di-AMP, and that it generates a strong innate immune response in eukaryotes. In bacteria that produce the second messenger, most molecular targets of cyclic di-AMP are associated with cell volume control. Besides, other evidence links the second messenger to cell wall remodeling, DNA damage repair, sporulation, central metabolism, and the regulation of glycogen turnover. In this review, we take a biochemical, quantitative approach to address the main cellular processes that are directly regulated by cyclic di-AMP and show that these processes are very connected and require regulation of a similar set of proteins to which cyclic di-AMP binds. Altogether, we argue that cyclic di-AMP is a master regulator of cell volume and that other cellular processes can be connected with cyclic di-AMP through this core function. We further highlight important directions in which the cyclic di-AMP field has to develop to gain a full understanding of the cyclic di-AMP signaling network and why some processes are regulated, while others are not.

The Whole Mitochondrial Genome Sequence of Dendrobium loddigesii Rolfe, an Endangered Orchid Species in China, Reveals a Complex Multi-Chromosome Structure

Dendrobium loddigesii is a precious traditional Chinese medicine with high medicinal and ornamental value. However, the characterization of its mitochondrial genome is still pending. Here, we assembled the complete mitochondrial genome of D. loddigesii and discovered that its genome possessed a complex multi-chromosome structure. The mitogenome of D. loddigesii consisted of 17 circular subgenomes, ranging in size from 16,323 bp to 56,781 bp. The total length of the mitogenome was 513,356 bp, with a GC content of 43.41%. The mitogenome contained 70 genes, comprising 36 protein-coding genes (PCGs), 31 tRNA genes, and 3 rRNA genes. Furthermore, we detected 403 repeat sequences as well as identified 482 RNA-editing sites and 8154 codons across all PCGs. Following the sequence similarity analysis, 27 fragments exhibiting homology to both the mitogenome and chloroplast genome were discovered, accounting for 9.86% mitogenome of D. loddigesii. Synteny analysis revealed numerous sequence rearrangements in D. loddigesii and the mitogenomes of related species. Phylogenetic analysis strongly supported that D. loddigesii and D. Amplum formed a single clade with 100% bootstrap support. The outcomes will significantly augment the orchid mitochondrial genome database, offering profound insights into Dendrobium's intricate mitochondrial genome architecture.

Insights into structure, codon usage, repeats, and RNA editing of the complete mitochondrial genome of Perilla frutescens (Lamiaceae)

Perilla frutescens (L.) Britton, a member of the Lamiaceae family, stands out as a versatile plant highly valued for its unique aroma and medicinal properties. Additionally, P. frutescens seeds are rich in Îś-linolenic acid, holding substantial economic importance. While the nuclear and chloroplast genomes of P. frutescens have already been documented, the complete mitochondrial genome sequence remains unreported. To this end, the sequencing, annotation, and assembly of the entire Mitochondrial genome of P. frutescens were hereby conducted using a combination of Illumina and PacBio data. The assembled P. frutescens mitochondrial genome spanned 299,551 bp and exhibited a typical circular structure, involving a GC content of 45.23%. Within the genome, a total of 59 unique genes were identified, encompassing 37 protein-coding genes, 20 tRNA genes, and 2 rRNA genes. Additionally, 18 introns were observed in 8 protein-coding genes. Notably, the codons of the P. frutescens mitochondrial genome displayed a notable A/T bias. The analysis also revealed 293 dispersed repeat sequences, 77 simple sequence repeats (SSRs), and 6 tandem repeat sequences. Moreover, RNA editing sites preferentially produced leucine at amino acid editing sites. Furthermore, 70 sequence fragments (12,680 bp) having been transferred from the chloroplast to the mitochondrial genome were identified, accounting for 4.23% of the entire mitochondrial genome. Phylogenetic analysis indicated that among Lamiaceae plants, P. frutescens is most closely related to Salvia miltiorrhiza and Platostoma chinense. Meanwhile, inter-species Ka/Ks results suggested that Ka/Ks < 1 for 28 PCGs, indicating that these genes were evolving under purifying selection. Overall, this study enriches the mitochondrial genome data for P. frutescens and forges a theoretical foundation for future molecular breeding research.

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

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.

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.

Comprehensive analysis of the mitochondrial genome of Rehmannia glutinosa: insights into repeat-mediated recombinations and RNA editing-induced stop codon acquisition

Rehmannia glutinosa is an economically significant medicinal plant. Yet, the structure and sequence of its mitochondrial genome has not been published, which plays a crucial role in evolutionary analysis and regulating respiratory-related macromolecule synthesis. In this study, the R. glutinosa mitogenome was sequenced employing a combination of Illumina short reads and Nanopore long reads, with subsequent assembly using a hybrid strategy. We found that the predominant configuration of the R. glutinosa mitogenome comprises two circular chromosomes. The primary structure of the mitogenome encompasses two mitochondrial chromosomes corresponding to the two major configurations, Mac1-1 and Mac1-2. The R. glutinosa mitogenome encoded an angiosperm-typical set of 24 core genes, nine variable genes, three rRNA genes, and 15 tRNA genes. A phylogenetic analysis using the 16 shared protein-coding genes (PCG) yielded a tree consistent with the phylogeny of Lamiales species and two outgroup taxa. Mapping RNA-seq data to the coding sequences (CDS) of the PCGs revealed 507 C-to-U RNA editing sites across 31 PCGs of the R. glutinosa mitogenome. Furthermore, one start codon (nad4L) and two stop codons (rpl10 and atp6) were identified as products of RNA editing events in the R. glutinosa mitogenome.

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing

BACKGROUND

Solanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date.

RESULTS

PacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum.

CONCLUSIONS

We provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

Elucidating the multichromosomal structure within the Brasenia schreberi mitochondrial genome through assembly and analysis

Brasenia schreberi, a plant species traditionally utilized in Chinese medicine and cuisine, represents an early evolutionary stage among flowering plants (angiosperms). While the plastid genome of this species has been published, its mitochondrial genome (mitogenome) has not been extensively explored, with a notable absence of thorough comparative analyses of its organellar genomes. In our study, we had assembled the entire mitogenome of B. schreberi utilizing the sequencing data derived from both Illumina platform and Oxford Nanopore. The B. schreberi mitogenome mostly exists as six circular DNA molecules, with the largest being 628,257 base pairs (bp) and the smallest 110,220 bp, amounting to 1.49 megabases (Mb). Then we annotated the mitogenome of B. schreberi. The mitogenome encompasses a total of 71 genes: 40 of these are coding proteins genes (PCGs), 28 are genes for transfer RNA (tRNA), and the remaining 3 are genes for ribosomal RNA (rRNA). In the analysis of codon usage, we noted a unique codon preference specific to each amino acid. The most commonly used codons exhibited an average RSCU of 1.36, indicating a noticeable bias in codon selection. In the repeat sequence analysis, a total of 553 simple sequence repeats (SSRs) were identified, 1,822 dispersed repeats (comprising 1,015 forward and 807 palindromic repeats), and 608 long terminal repeats (LTRs). Additionally, in the analysis of homologous sequences between organelle genomes, we detected 38 homologous sequences derived from the plastid genome, each exceeding 500 bp, within the B. schreberi mitochondrial genome. Notably, ten tRNA genes (trnC-GCA, trnM-CAU, trnI-CAU, trnQ-UUG, trnN-GUU, trnT-GGU, trnW-CCA, trnA-UGC, trnI-GAU, and trnV-GAC) appear to have been completely transferred from the chloroplast to the mitogenome. Utilizing the Deepred-mt to predict the RNA editing sites in the mitogenome, we have identified 675 high-quality RNA editing sites in the 40 mitochondrial PCGs. In the final stage of our study, we performed an analysis of colinearity and inferred the phylogenetic relationship of B. schreberi with other angiosperms, utilizing the mitochondrial PCGs as a basis. The results showed that the non-coding regions of the B. schreberi mitogenome are characterized by an abundance of repetitive sequences and exogenous sequences, and B. schreberi is more closely related with Euryale ferox.

The complete mitochondrial genome of Aglaia odorata, insights into its genomic structure and RNA editing sites

Aglaia odorata, native to Guangdong, Guangxi, and Hainan provinces in China, has long been utilized as an herbal remedy in ancient China. In this study, we assembled and annotated the complete mitochondrial genome (mitogenome) of A. odorata, which spans a total length of 537,321 bp. Conformation of the A. odorata recombination was verified through PCR experiments and Sanger sequencing. We identified and annotated 35 protein-coding genes (PCGs), 22 tRNA genes, and 3 rRNA genes within the mitogenome. Analysis of repeated elements revealed the presence of 192 SSRs, 29 pairs of tandem repeats, and 333 pairs of dispersed repeats in the A. odorata mitogenome. Additionally, we analyzed codon usage and mitochondrial plastid DNAs (MTPTs). Twelve MTPTs between the plastome and mitogenome of A. odorata were identified, with a combined length of 2,501 bp, accounting for 0.47% of the mitogenome. Furthermore, 359 high-confidence C to U RNA editing sites were predicted on PCGs, and four selected RNA editing sites were specially examined to verify the creation of start and/or stop codons. Extensive genomic rearrangement was observed between A. odorata and related mitogenomes. Phylogenetic analysis based on mitochondrial PCGs were conducted to elucidate the evolutionary relationships between A. odorata and other angiosperms.

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.

Comparative analysis of chloroplast and mitochondrial genomes of sweet potato provides evidence of gene transfer

The increasing number of plant mitochondrial DNA genomes (mtDNA) sequenced reveals the extent of transfer from both chloroplast DNA genomes (cpDNA) and nuclear DNA genomes (nDNA). This study created a library and assembled the chloroplast and mitochondrial genomes of the leafy sweet potato better to understand the extent of mitochondrial and chloroplast gene transfer. The full-length chloroplast genome of the leafy sweet potato (OM808940) is 161,387 bp, with 132 genes annotated, including 87 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. The mitochondrial genome (OM808941) was 269,578 bp in length and contained 69 functional genes, including 39 protein-coding genes, 6 rRNA genes, and 24 tRNA genes. 68 SSR loci were found in the leafy sweet potato organelle genome, including 54 in the chloroplast genome and 14 in the mitochondria genome. In the sweet potato mitochondrial genome, most genes have RNA editing sites, and the conversion ratio from hydrophilic amino acids to hydrophobic amino acids is the highest, reaching 47.12%. Horizontal transfer occurs in the sweet potato organelle genome and nuclear genome. 40 mitochondrial genome segments share high homology with 14 chloroplast genome segments, 33 of which may be derived from chloroplast genome horizontal transfer. 171 mitochondrial genome sequences come from the horizontal transfer of nuclear genome. The phylogenetic analysis of organelle genes revealed that the leafy sweet potato was closely related to the tetraploid wild species Ipomoea tabascana and the wild diploid species Ipomoea trifida.

MtDNA deletions and aging

Aging is the major risk factor in most of the leading causes of mortality worldwide, yet its fundamental causes mostly remain unclear. One of the clear hallmarks of aging is mitochondrial dysfunction. Mitochondria are best known for their roles in cellular energy generation, but they are also critical biosynthetic and signaling organelles. They also undergo multiple changes with organismal age, including increased genetic errors in their independent, circular genome. A key group of studies looking at mice with increased mtDNA mutations showed that premature aging phenotypes correlated with increased deletions but not point mutations. This generated an interest in mitochondrial deletions as a potential fundamental cause of aging. However, subsequent studies in different models have yielded diverse results. This review summarizes the research on mitochondrial deletions in various organisms to understand their possible roles in causing aging while identifying the key complications in quantifying deletions across all models.

Charge cluster occurrence in land plants' mitochondrial proteomes with functional and structural insights

The Charge Clusters (CCs) are involved in key functions and are distributed according to the organism, the protein's type, and the charge of amino acids. In the present study, we have explored the occurrence, position, and annotation as a first large-scale study of the CCs in land plants mitochondrial proteomes. A new python script was used for data curation. The Finding Clusters Charge in Protein Sequences Program was performed after adjusting the reading window size. A 44316 protein sequences belonging to 52 species of land plants were analysed. The occurrence of Negative Charge Clusters (NCCs) (1.2%) is two times more frequent than the Positive Charge Clusters (PCCs) (0.64%). Moreover, 39 and 30 NCCs were conserved in 88 and 41 proteins in intra and in inter proteomes respectively, while 14 and 21 PCCs were conserved in 53 and 85 protein sequences in intra and inter proteomes consecutively. Sequences carrying mixed CCs are rare (0.12%). Despite this low abundance, CCs play a crucial role in protein function. The CCs tend to be located mainly in the terminal regions of proteins which guarantees specific protein targeting and import into the mitochondria. In addition, the functional annotation of CCs according to Gene Ontology shows that CCs are involved in binding functions of either proteins or macromolecules which are deployed in different metabolic and cellular processes such as RNA editing and transcription. This study may provide valuable information while considering the CCs in understanding the environmental adaptation of plants.Communicated by Ramaswamy H. Sarma.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

Primase promotes the competition between transcription and replication on the same template strand resulting in DNA damage

Transcription-replication conflicts (TRCs), especially Head-On TRCs (HO-TRCs) can introduce R-loops and DNA damage, however, the underlying mechanisms are still largely unclear. We previously identified a chloroplast-localized RNase H1 protein AtRNH1C that can remove R-loops and relax HO-TRCs for genome integrity. Through the mutagenesis screen, we identify a mutation in chloroplast-localized primase ATH that weakens the binding affinity of DNA template and reduces the activities of RNA primer synthesis and delivery. This slows down DNA replication, and reduces competition of transcription-replication, thus rescuing the developmental defects of atrnh1c. Strand-specific DNA damage sequencing reveals that HO-TRCs cause DNA damage at the end of the transcription unit in the lagging strand and overexpression of ATH can boost HO-TRCs and exacerbates DNA damage. Furthermore, mutation of plastid DNA polymerase Pol1A can similarly rescue the defects in atrnh1c mutants. Taken together these results illustrate a potentially conserved mechanism among organisms, of which the primase activity can promote the occurrence of transcription-replication conflicts leading to HO-TRCs and genome instability.

Comparative mitochondrial genome analysis reveals a candidate ORF for cytoplasmic male sterility in tropical onion

Cytoplasmic male sterility (CMS) has been widely exploited for hybrid seed production in onions (Allium cepa L.). In contrast to long-day onion cultivars, short-day onion has not yet been investigated for mitochondrial genome structure and DNA rearrangements associated with CMS activity. Here, we report the 3,16,321 bp complete circular mitochondrial genome of tropical onion CMS line (97A). Due to the substantial number of repetitive regions, the assembled mitochondrial genome of maintainer line (97B) remained linear with 15 scaffolds. Additionally, 13 and 20 chloroplast-derived fragments with a size ranging from 143 to 13,984 bp and 153-17,725 bp were identified in the 97A and 97B genomes, respectively. Genome annotation revealed 24 core protein-coding genes along with 24 and 28 tRNA genes in the mitochondrial genomes of 97A and 97B, respectively. Furthermore, comparative genome analysis of the 97A and 97B mitochondrial genomes showed that gene content was almost similar except for the chimeric ORF725 gene which is the extended form of the COX1 gene.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03850-2.

Phylogeny and genetic variations of the three genome compartments in haptophytes shed light on the rapid evolution of coccolithophores

Haptophyte algae, including coccolithophores, play key roles in global carbon cycling and ecosystem. They exhibit exceptional morphological and functional diversity. However, their phylogeny is mostly based on short markers and genome researches are always limited to few species, hindering a better understanding about their evolution and diversification. In this study, by assembling 69 new plastid genomes, 65 new mitochondrial genomes, and 55 nuclear drafts, we systematically analyzed their genome variations and built the most comprehensive phylogenies in haptophytes and Noelaerhabdaceae, with the latter is the family of the model coccolithophore Emiliania huxleyi. The haptophyte genomes vary significantly in size, gene content, and structure. We detected phylogenetic incongruence of Prymnesiales between genome compartments. In Noelaerhabdaceae, by including Reticulofenestra sessilis and a proper outgroup, we found R. sessilis was not the basal taxon of this family. Noelaerhabdaceae strains have very similar genomic features and conserved sequences, but different gene content and dynamic structure. We speculate that was caused by DNA double-strand break repairs. Our results provide valuable genetic resources and new insights into the evolution of haptophytes, especially coccolithophores.

Assembly and analysis of the complete mitochondrial genome of Forsythia suspensa (Thunb.) Vahl

BACKGROUND

Forsythia suspensa (Thunb.) Vahl is a valuable ornamental and medicinal plant. Although the nuclear and chloroplast genomes of F. suspensa have been published, its complete mitochondrial genome sequence has yet to be reported. In this study, the genomic DNA of F. suspensa yellowish leaf material was extracted, sequenced by using a mixture of Illumina Novaseq6000 short reads and Oxford Nanopore PromethION long reads, and the sequencing data were assembled and annotated.

RESULT

The F. suspensa mitochondrial genome was obtained in the length of 535,692 bp with a circular structure, and the GC content was 44.90%. The genome contains 60 genes, including 36 protein-coding genes, 21 tRNA genes, and three rRNA genes. We further analyzed RNA editing of the protein-coding genes, relative synonymous codon usage, and sequence repeats based on the genomic data. There were 25 homologous sequences between F. suspensa mitochondria and chloroplast genome, which involved the transfer of 8 mitochondrial genes, and 9473 homologous sequences between mitochondrial and nuclear genomes. Analysis of the nucleic acid substitution rate, nucleic acid diversity, and collinearity of protein-coding genes of the F. suspensa mitochondrial genome revealed that the majority of genes may have undergone purifying selection, exhibiting a slower rate of evolution and a relatively conserved structure. Analysis of the phylogenetic relationships among different species revealed that F. suspensa was most closely related to Olea europaea subsp. Europaea.

CONCLUSION

In this study, we sequenced, assembled, and annotated a high-quality F. suspensa mitochondrial genome. The results of this study will enrich the mitochondrial genome data of Forsythia, lay a foundation for the phylogenetic development of Forsythia, and promote the evolutionary analysis of Oleaceae species.

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.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications

BACKGROUND

Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported.

RESULTS

This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3).

CONCLUSIONS

This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.

Complete Mitochondrial Genome Assembly of an Upland Wild Rice Species, Oryza granulata and Comparative Mitochondrial Genomic Analyses of the Genus Oryza

Wild upland rice species, including Oryza granulata, possess unique characteristics that distinguish them from other Oryza species. For instance, O. granulata characteristically has a GG genome and is accordingly classified as a basal lineage of the genus Oryza. Here, we deployed a versatile hybrid approach by integrating Illumina and PacBio sequencing data to generate a high-quality mitochondrial genome (mitogenome) assembly for O. granulata. The mitogenome of O. granulata was 509,311 base pairs (bp) with sixty-seven genes comprising two circular chromosomes, five ribosomal RNA (rRNA) coding genes, twenty-five transfer RNA (tRNA) coding genes, and thirty-seven genes coding for proteins. We identified a total of 378 simple sequence repeats (SSRs). The genome also contained 643 pairs of dispersed repeats comprising 340 palindromic and 303 forward. In the O. granulata mitogenome, the length of 57 homologous fragments in the chloroplast genome occupied 5.96% of the mitogenome length. Collinearity analysis of three Oryza mitogenomes revealed high structural variability and frequent rearrangements. Phylogenetic analysis showed that, compared to other related genera, O. granulata had the closest genetic relationship with mitogenomes reported for all members of Oryza, and occupies a position at the base of the Oryza phylogeny. Comparative analysis of complete mitochondrial genome assemblies for Oryza species revealed high levels of mitogenomic diversity, providing a foundation for future conservation and utilization of wild rice biodiversity.

Assembly and comparative analysis of the complete mitochondrial genome of Viburnum chinshanense

BACKGROUND

Viburnum chinshanense is an endemic species found exclusively in the North-Central and South-Central regions of China. This species is a lush garden ornamental tree and is extensively utilized for vegetation restoration in rocky desertification areas.

RESULTS

In this study, we obtained 13.96 Gb of Oxford Nanopore data for the whole genome, and subsequently, by combining Illumina short-reads, we successfully assembled the complete mitochondrial genome (mitogenome) of the V. chinshanense using a hybrid assembly strategy. The assembled genome can be described as a circular genome. The total length of the V. chinshanense mitogenome measures 643,971 bp, with a GC content of 46.18%. Our annotation efforts have revealed a total of 39 protein-coding genes (PCGs), 28 tRNA genes, and 3 rRNA genes within the V. chinshanense mitogenome. The analysis of repeated elements has identified 212 SSRs, 19 long tandem repeat elements, and 325 pairs of dispersed repeats in the V. chinshanense mitogenome. Additionally, we have investigated mitochondrial plastid DNAs (MTPTs) and identified 21 MTPTs within the mitogenome and plastidial genome. These MTPTs collectively span a length of 9,902 bp, accounting for 1.54% of the mitogenome. Moreover, employing Deepred-mt, we have confidently predicted 623 C to U RNA editing sites across the 39 protein-coding genes. Furthermore, extensive genomic rearrangements have been observed between V. chinshanense and the mitogenomes of related species. Interestingly, we have also identified a bacterial-derived tRNA gene (trnC-GCA) in the V. chinshanense mitogenome. Lastly, we have inferred the phylogenetic relationships of V. chinshanense with other angiosperms based on mitochondrial PCGs.

CONCLUSIONS

This study marks the first report of a mitogenome from the Viburnum genus, offering a valuable genomic resource for exploring the evolution of mitogenomes within the Dipsacales order.

Genetic diversity of Coffea arabica L. mitochondrial genomes caused by repeat- mediated recombination and RNA editing

Background

Coffea arabica L. is one of the most important crops widely cultivated in 70 countries across Asia, Africa, and Latin America. Mitochondria are essential organelles that play critical roles in cellular respiration, metabolism, and differentiation. C. arabica's nuclear and chloroplast genomes have been reported. However, its mitochondrial genome remained unreported. Here, we intended to sequence and characterize its mitochondrial genome to maximize the potential of its genomes for evolutionary studies, molecular breeding, and molecular marker developments.

Results

We sequenced the total DNA of C. arabica using Illumina and Nanopore platforms. We then assembled the mitochondrial genome with a hybrid strategy using Unicycler software. We found that the mitochondrial genome comprised two circular chromosomes with lengths of 867,678 bp and 153,529 bp, encoding 40 protein-coding genes, 26 tRNA genes, and three rRNA genes. We also detected 270 Simple Sequence Repeats and 34 tandem repeats in the mitochondrial genome. We found 515 high-scoring sequence pairs (HSPs) for a self-to-self similarity comparison using BLASTn. Three HSPs were found to mediate recombination by the mapping of long reads. Furthermore, we predicted 472 using deep-mt with the convolutional neural network model. Then we randomly validated 90 RNA editing events by PCR amplification and Sanger sequencing, with the majority being non-synonymous substitutions and only three being synonymous substitutions. These findings provide valuable insights into the genetic characteristics of the C. arabica mitochondrial genome, which can be helpful for future study on coffee breeding and mitochondrial genome evolution.

Conclusion

Our study sheds new light on the evolution of C. arabica organelle genomes and their potential use in genetic breeding, providing valuable data for developing molecular markers that can improve crop productivity and quality. Furthermore, the discovery of RNA editing events in the mitochondrial genome of C. arabica offers insights into the regulation of gene expression in this species, contributing to a better understanding of coffee genetics and evolution.

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.

Structural and gene composition variation of the complete mitochondrial genome of Mammillaria huitzilopochtli (Cactaceae, Caryophyllales), revealed by de novo assembly

BACKGROUND

Structural descriptions of complete genomes have elucidated evolutionary processes in angiosperms. In Cactaceae (Caryophyllales), a high structural diversity of the chloroplast genome has been identified within and among genera. In this study, we assembled the first mitochondrial genome (mtDNA) for the short-globose cactus Mammillaria huitzilopochtli. For comparative purposes, we used the published genomes of 19 different angiosperms and the gymnosperm Cycas taitungensis as an external group for phylogenetic issues.

RESULTS

The mtDNA of M. huitzilopochtli was assembled into one linear chromosome of 2,052,004 bp, in which 65 genes were annotated. These genes account for 57,606 bp including 34 protein-coding genes (PCGs), 27 tRNAs, and three rRNAs. In the non-coding sequences, repeats were abundant, with a total of 4,550 (179,215 bp). In addition, five complete genes (psaC and four tRNAs) of chloroplast origin were documented. Negative selection was estimated for most (23) of the PCGs. The phylogenetic tree showed a topology consistent with previous analyses based on the chloroplast genome.

CONCLUSIONS

The number and type of genes contained in the mtDNA of M. huitzilopochtli were similar to those reported in 19 other angiosperm species, regardless of their phylogenetic relationships. Although other Caryophyllids exhibit strong differences in structural arrangement and total size of mtDNA, these differences do not result in an increase in the typical number and types of genes found in M. huitzilopochtli. We concluded that the total size of mtDNA in angiosperms increases by the lengthening of the non-coding sequences rather than a significant gain of coding genes.

Genome Survey and Chromosome-Level Draft Genome Assembly of Glycine max var. Dongfudou 3: Insights into Genome Characteristics and Protein Deficiencies

Dongfudou 3 is a highly sought-after soybean variety due to its lack of beany flavor. To support molecular breeding efforts, we conducted a genomic survey using next-generation sequencing. We determined the genome size, complexity, and characteristics of Dongfudou 3. Furthermore, we constructed a chromosome-level draft genome and speculated on the molecular basis of protein deficiency in GmLOX1, GmLOX2, and GmLOX3. These findings set the stage for high-quality genome analysis using third-generation sequencing. The estimated genome size is approximately 1.07 Gb, with repetitive sequences accounting for 72.50%. The genome is homozygous and devoid of microbial contamination. The draft genome consists of 916.00 Mb anchored onto 20 chromosomes, with annotations of 46,446 genes and 77,391 transcripts, achieving Benchmarking Single-Copy Orthologue (BUSCO) completeness of 99.5% for genome completeness and 99.1% for annotation. Deletions and substitutions were identified in the three GmLox genes, and they also lack corresponding active proteins. Our proposed approach, involving k-mer analysis after filtering out organellar DNA sequences, is applicable to genome surveys of all plant species, allowing for accurate assessments of size and complexity. Moreover, the process of constructing chromosome-level draft genomes using closely related reference genomes offers cost-effective access to valuable information, maximizing data utilization.

Complete mitochondrial genome of Agrostis stolonifera: insights into structure, Codon usage, repeats, and RNA editing

BACKGROUND

Plants possess mitochondrial genomes that are large and complex compared to animals. Despite their size, plant mitochondrial genomes do not contain significantly more genes than their animal counterparts. Studies into the sequence and structure of plant mitochondrial genomes heavily imply that the main mechanism driving replication of plant mtDNA, and offer valuable insights into plant evolution, energy production, and environmental adaptation.

RESULTS

This study presents the first comprehensive analysis of Agrostis stolonifera's mitochondrial genome, characterized by a branched structure comprising three contiguous chromosomes, totaling 560,800 bp with a GC content of 44.07%. Annotations reveal 33 unique protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The predominant codons for alanine and glutamine are GCU and CAA, respectively, while cysteine and phenylalanine exhibit weaker codon usage biases. The mitogenome contains 73, 34, and 23 simple sequence repeats (SSRs) on chromosomes 1, 2, and 3, respectively. Chromosome 1 exhibits the most frequent A-repeat monomeric SSR, whereas chromosome 2 displays the most common U-repeat monomeric SSR. DNA transformation analysis identifies 48 homologous fragments between the mitogenome and chloroplast genome, representing 3.41% of the mitogenome's total length. The PREP suite detects 460 C-U RNA editing events across 33 mitochondrial PCGs, with the highest count in the ccmFn gene and the lowest in the rps7 gene. Phylogenetic analysis confirms A. stolonifera's placement within the Pooideae subfamily, showing a close relationship to Lolium perenne, consistent with the APG IV classification system. Numerous homologous co-linear blocks are observed in A. stolonifera's mitogenomes and those of related species, while certain regions lack homology.

CONCLUSIONS

The unique features and complexities of the A. stolonifera mitochondrial genome, along with its similarities and differences to related species, provide valuable insights into plant evolution, energy production, and environmental adaptation. The findings from this study significantly contribute to the growing body of knowledge on plant mitochondrial genomes and their role in plant biology.

Comparative analysis of the complete mitogenome of Geoffroea decorticans: a native tree surviving in the Atacama Desert

Chañar (Geoffroea decorticans (Gill., ex Hook. & Arn.) Burkart) has been highly significant for indigenous people in the Atacama Desert for over 3,000 years. Through evolutionary processes, the G. decorticans mitogenome likely underwent changes facilitating its adaptation to the extreme conditions of the Atacama Desert. Here, we compare the mitochondrial genome of G. decorticans with those of other Papilionoideae family species. The complete mitogenome of G. decorticans was sequenced and assembled, making it the first in the genus Geoffroea. The mitogenome contained 383,963 base pairs, consisting of 33 protein coding genes, 21 transfer RNA genes, and 3 ribosomal RNA genes. The Chañar mitogenome is relatively compact, and has two intact genes (sdh4 and nad1) which were not observed in most other species. Additionally, Chañar possessed the highest amount of mitochondrial DNA of plastid origin among angiosperm species. The phylogenetic analysis of the mitogenomes of Chañar and 12 other taxa displayed a high level of consistency in taxonomic classification, when compared to those of the plastid genome. Atp8 was subjected to positive selection, while the ccmFc and rps1 were subjected to neutral selection. This study provides valuable information regarding its ability to survive the extreme environmental conditions of the Atacama Desert.

Assembly and analysis of the first complete mitochondrial genome of Punica granatum and the gene transfer from chloroplast genome

Pomegranate (Punica granatum L.) is one of the oldest fruits with edible, medicinal and ornamental values. However, there is no report on the mitochondrial genome of pomegranate. In this study, the mitochondrial genome of P. granatum was sequenced, assembled and analyzed in detail, while the chloroplast genome was assembled using the same set of data. The results showed that the P. granatum mitogenome had a multi branched structure, using BGI + Nanopore mixed assembly strategy. The total genome length was 404,807 bp, with the GC content of 46.09%, and there were 37 protein coding genes, 20 tRNA genes and three rRNA genes. In the whole genome, 146 SSRs were identified. Besides, 400 pairs of dispersed repeats were detected, including 179 palindromic, 220 forward and one reverse. In the P. granatum mitochondrial genome, 14 homologous fragments of chloroplast genome were found, accounting for 0.54% of the total length. Phylogenetic analysis showed that among the published mitochondrial genomes of related genera, P. granatum had the closest genetic relationship with Lagerstroemia indica of Lythraceae. The 580 and 432 RNA editing sites were predicted on 37 protein coding genes of mitochondrial genome using BEDTools software and online website PREPACT respectively, but all were from C to U, of which ccmB and nad4 gene were most frequently edited, with 47 sites. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of pomegranate germplasm resources.

Oxidative and Glycation Damage to Mitochondrial DNA and Plastid DNA during Plant Development

Oxidative damage to plant proteins, lipids, and DNA caused by reactive oxygen species (ROS) has long been studied. The damaging effects of reactive carbonyl groups (glycation damage) to plant proteins and lipids have also been extensively studied, but only recently has glycation damage to the DNA in plant mitochondria and plastids been reported. Here, we review data on organellar DNA maintenance after damage from ROS and glycation. Our focus is maize, where tissues representing the entire range of leaf development are readily obtained, from slow-growing cells in the basal meristem, containing immature organelles with pristine DNA, to fast-growing leaf cells, containing mature organelles with highly-fragmented DNA. The relative contributions to DNA damage from oxidation and glycation are not known. However, the changing patterns of damage and damage-defense during leaf development indicate tight coordination of responses to oxidation and glycation events. Future efforts should be directed at the mechanism by which this coordination is achieved.

Rickettsial DNA and a trans-splicing rRNA group I intron in the unorthodox mitogenome of the fern Haplopteris ensiformis

Plant mitochondrial genomes can be complex owing to highly recombinant structures, lack of gene syntenies, heavy RNA editing and invasion of chloroplast, nuclear or even foreign DNA by horizontal gene transfer (HGT). Leptosporangiate ferns remained the last major plant clade without an assembled mitogenome, likely owing to a demanding combination of the above. We here present both organelle genomes now for Haplopteris ensiformis. More than 1,400 events of C-to-U RNA editing and over 500 events of reverse U-to-C edits affect its organelle transcriptomes. The Haplopteris mtDNA is gene-rich, lacking only the ccm gene suite present in ancestral land plant mitogenomes, but is highly unorthodox, indicating extraordinary recombinogenic activity. Although eleven group II introns known in disrupted trans-splicing states in seed plants exist in conventional cis-arrangements, a particularly complex structure is found for the mitochondrial rrnL gene, which is split into two parts needing reassembly on RNA level by a trans-splicing group I intron. Aside from ca. 80 chloroplast DNA inserts that complicated the mitogenome assembly, the Haplopteris mtDNA features as an idiosyncrasy 30 variably degenerated protein coding regions from Rickettiales bacteria indicative of heavy bacterial HGT on top of tRNA genes of chlamydial origin.

De Novo Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Chaenomeles speciosa (Sweet) Nakai Revealed the Existence of Two Structural Isomers

As a valuable Chinese traditional medicinal species, Chaenomeles speciosa (Sweet) Nakai (C. speciosa) is a natural resource with significant economic and ornamental value. However, its genetic information is not well understood. In this study, the complete mitochondrial genome of C. speciosa was assembled and characterized to explore the repeat sequences, recombination events, rearrangements, and IGT, to predict RNA editing sites, and to clarify the phylogenetic and evolutionary relationship. The C. speciosa mitochondrial genome was found to have two circular chromosomes as its major conformation, with a total length of 436,464 bp and 45.2% GC content. The mitochondrial genome contained 54 genes, including 33 unique protein-coding genes, 18 tRNAs, and 3 rRNA genes. Seven pairs of repeat sequences involving recombination events were analyzed. Both the repeat pairs, R1 and R2, played significant roles in mediating the major and minor conformations. In total, 18 MTPTs were identified, 6 of which were complete tRNA genes. There were 454 RNA editing sites in the 33 protein-coding sequences predicted by the PREPACT3 program. A phylogenetic analysis based on 22 species of mitochondrial genomes was constructed and indicated highly conserved PCG sequences. Synteny analyses showed extensive genomic rearrangements in the mitochondrial genome of C. speciosa and closely related species. This work is the first to report the C. speciosa mitochondrial genome, which is of great significance for conducting additional genetic studies on this organism.

Two-Step Identification of N-, S-, R- and T-Cytoplasm Types in Onion Breeding Lines Using High-Resolution Melting (HRM)-Based Markers

High-resolution melting (HRM) analysis is a powerful detection method for fast, high-throughput post-PCR analysis. A two-step HRM marker system was developed for identification of the N-, S-, R- and T-cytoplasms of onion. In the first step for the identification of N-, S- and R-cytoplasms, one forward primer was designed to the identical sequences of both cox1 and orf725 genes, and two reverse primers specific to the polymorphic sequences of cox1 and orf725 genes were used. For the second step, breeding lines with N-cytoplasm were evaluated with primers developed from the orfA501 sequence to distinguish between N- and T-cytoplasms. An amplicon with primers to the mitocondrial atp9 gene was used as an internal control. The two-step HRM marker system was tested using 246 onion plants. HRM analysis showed that the most common source of CMS, often used by Russian breeders, was S-cytoplasm; the rarest type of CMS was R-cytoplasm; and the proportion of T-cytoplasm among the analyzed breeding lines was 20.5%. The identification of the cytoplasm of a single plant by phenotype takes from 4 to 8 years. The HRM-based system enables quick and easy distinguishing of the four types of onion cytoplasm.