Ultra-deep sequencing reveals dramatic alteration of organellar genomes in Physcomitrella patens due to biased asymmetric recombination

Mitochondrial genome structural variants and candidate cytoplasmic male sterility-related gene in sugarcane

BACKGROUND

Sugarcane is a crucial crop for both sugar and bioethanol production. The nobilization breeding and utilization of wild germplasm have significantly enhanced its productivity. However, the pollen sterility in Saccharum officinarum restricts its role to being a female parent in crosses with Saccharum spontaneum during nobilization breeding, resulting in a narrow genetic basis for modern sugarcane cultivars. Mitochondria, often referred to as the intracellular "energy factories", provide energy for plant life activities, and are also implicated in cytoplasmic male sterility (CMS).

RESULTS

We performed mitochondrial genome assembly and structural analysis of two Saccharum founding species. We discovered that the proportions of repeat sequences are the primary factor contributing to the variations in mitochondrial genome structure and size between the two Saccharum species. Heterologous expression of the mitochondrial chimeric gene ORF113, which is highly expressed in male-sterile S. officinarum flowers, significantly inhibits growth and ATP synthesis in yeast cells, making it a key candidate CMS-related gene in sugarcane. Furthermore, we developed two co-dominant simple sequence repeat (SSR) markers based on the mitochondrial genome, which can effectively distinguish the cytoplasmic types of the two Saccharum species.

CONCLUSION

In this study, we identified structural variants and developed SSR molecular markers in the mitochondrial genomes of both S. officinarum and S. spontaneum. We also identified a novel mitochondrial chimeric ORF as a key candidate CMS-related gene. These findings offer valuable insights into variety identification, genetic resource development, and cross-breeding strategies in sugarcane.

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.

The complete mitochondrial genome of Castanopsis carlesii and Castanea henryi reveals the rearrangement and size differences of mitochondrial DNA molecules

BACKGROUND

Castanopsis carlesii is a dominant tree species in subtropical evergreen broad-leaved forests and holds significant ecological value. It serves as an excellent timber tree species and raw material for cultivating edible fungi. Henry Chinquapin (Castanea henryi) wood is known for its hardness and resistance to water and moisture, making it an exceptional timber species. Additionally, its fruit has a sweet and fruity taste, making it a valuable food source. However, the mitogenomes of these species have not been previously reported. To gain a better understanding of them, this study successfully assembled high-quality mitogenomes of C. carlesii and Ca. henryi for the first time.

RESULTS

Our research reveals that the mitochondrial DNA (mtDNA) of C. carlesii exhibits a unique multi-branched conformation, while Ca. henryi primarily exists in the form of two independent molecules that can be further divided into three independent molecules through one pair of long repetitive sequences. The size of the mitogenomes of C. carlesii and Ca. henryi are 592,702 bp and 379,929 bp respectively, which are currently the largest and smallest Fagaceae mitogenomes recorded thus far. The primary factor influencing mitogenome size is dispersed repeats. Comparison with published mitogenomes from closely related species highlights differences in size, gene loss patterns, codon usage preferences, repetitive sequences, as well as mitochondrial plastid DNA segments (MTPTs).

CONCLUSIONS

Our study enhances the understanding of mitogenome structure and evolution in Fagaceae, laying a crucial foundation for future research on cell respiration, disease resistance, and other traits in this family.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

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.

The first complete mitochondrial genome of Carex (C. breviculmis): a significantly expanded genome with highly structural variations

MAIN CONCLUSION

The first complete mitochondrial genome of Carex (C. breviculmis) was sequenced and assembled, and its genomic signature was analyzed and the possible conformations of its mitochondrial genome were validated. Carex breviculmis is a very adaptable grass that is highly resistant to environmental stresses such as drought and low light. It is also admired as a landscape plant with high development prospects and scientific research value. In this study, the mitochondrial genome of C. breviculmis was assembled using Pacbio and Illumina sequencing data. We detected 267 pairs of repeats and found that three pairs of repeats could mediate the recombination of its mitochondrial genome and formed four possible conformations, of which we verified the two conformations mediated by the shortest pair of repeats using PCR amplification and Sanger sequencing. The major conformation of the C. breviculmis mitochondrial genome is a 1,414,795 bp long circular molecule with 33 annotated protein-coding genes, 15 tRNA genes, and three rRNA genes. We detected a total of 25 homologous sequences between the chloroplast and mitochondrial genomes, corresponding to 0.40% of the mitochondrial genome. Combined with the low GC content (41.24%), we conclude that the reduction in RNA editing sites in the C. breviculmis mitochondrial genome may be due to an accumulation of point mutations in C-to-T or retroprocessing events within the genome. The relatively low number of RNA editing sites in its mitochondrial genome could serve as important material for subsequent studies on the selection pressure of RNA editing in angiosperms. A maximum likelihood analysis based on 23 conserved mitochondrial genes from 28 species reflects an accurate evolutionary and taxonomic position of C. breviculmis. This research provided us with a comprehensive understanding of the mitochondrial genome of Carex and also provided important information for its molecular breeding.

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.

The evolution of extremely diverged plastomes in Selaginellaceae (lycophyte) is driven by repeat patterns and the underlying DNA maintenance machinery

Two factors are proposed to account for the unusual features of organellar genomes: the disruptions of organelle-targeted DNA replication, repair, and recombination (DNA-RRR) systems in the nuclear genome and repetitive elements in organellar genomes. Little is known about how these factors affect organellar genome evolution. The deep-branching vascular plant family Selaginellaceae is known to have a deficient DNA-RRR system and convergently evolved organellar genomes. However, we found that the plastid genome (plastome) of Selaginella sinensis has extremely accelerated substitution rates, a low GC content, pervasive repeat elements, a dynamic network structure, and it lacks direct or inverted repeats. Unexpectedly, its organelle DNA-RRR system is short of a plastid-targeted Recombinase A1 (RecA1) and a mitochondrion-targeted RecA3, in line with other explored Selaginella species. The plastome contains a large collection of short- and medium-sized repeats. Given the absence of RecA1 surveillance, we propose that these repeats trigger illegitimate recombination, accelerated mutation rates, and structural instability. The correlations between repeat quantity and architectural complexity in the Selaginella plastomes support these conclusions. We, therefore, hypothesize that the interplay of the deficient DNA-RRR system and the high repeat content has led to the extraordinary divergence of the S. sinensis plastome. Our study not only sheds new light on the mechanism of plastome divergence by emphasizing the power of cytonuclear integration, but it also reconciles the longstanding contradiction on the effects of DNA-RRR system disruption on genome structure evolution.

Complex Physical Structure of Complete Mitochondrial Genome of Quercus acutissima (Fagaceae): A Significant Energy Plant

Quercus acutissima Carruth. is a Chinese important energy plant with high ecological and economic values. While the species chloroplast genome has been reported, its mitochondrial genome (mitogenome) is still unexplored. Here, we assembled and annotated the Q. acutissima mitogenome, and we compared its characteristic differences with several closely related species. The Q. acutissima mitogenome’s main structure is branched with three distinguished contigs (linear molecule 1, circular molecule 2, and circular molecule 3) with 448,982 bp total length and 45.72% GC content. The mitogenome contained 51 genes, including 32 protein-coding, 16 tRNA and 3 rRNA genes. We examined codon usage, repeated sequences, genome recombination, chloroplast to mitochondrion DNA transformation, RNA editing, and synteny in the Q. acutissima mitogenome. Phylogenetic trees based on 29 species mitogenomes clarified the species classification. Our results provided comprehensive information of Q. acutissima mitogenome, and they are expected to provide valuable information for Fagaceae evolutionary biology and to promote the species germplasm utilization.

The complete mitochondrial genome of okra (Abelmoschus esculentus): using nanopore long reads to investigate gene transfer from chloroplast genomes and rearrangements of mitochondrial DNA molecules

BACKGROUND

Okra (Abelmoschus esculentus L. Moench) is an economically important crop and is known for its slimy juice, which has significant scientific research value. The A. esculentus chloroplast genome has been reported; however, the sequence of its mitochondrial genome is still lacking.

RESULTS

We sequenced the plastid and mitochondrial genomes of okra based on Illumina short reads and Nanopore long reads and conducted a comparative study between the two organelle genomes. The plastid genome of okra is highly structurally conserved, but the mitochondrial genome of okra has been confirmed to have abundant subgenomic configurations. The assembly results showed that okra's mitochondrial genome existed mainly in the form of two independent molecules, which could be divided into four independent molecules through two pairs of long repeats. In addition, we found that four pairs of short repeats could mediate the integration of the two independent molecules into one complete molecule at a low frequency. Subsequently, we also found extensive sequence transfer between the two organelles of okra, where three plastid-derived genes (psaA, rps7 and psbJ) remained intact in the mitochondrial genome. Furthermore, psbJ, psbF, psbE and psbL were integrated into the mitochondrial genome as a conserved gene cluster and underwent pseudogenization as nonfunctional genes. Only psbJ retained a relatively complete sequence, but its expression was not detected in the transcriptome data, and we speculate that it is still nonfunctional. Finally, we characterized the RNA editing events of protein-coding genes located in the organelle genomes of okra.

CONCLUSIONS

In the current study, our results not only provide high-quality organelle genomes for okra but also advance our understanding of the gene dialogue between organelle genomes and provide information to breed okra cultivars efficiently.

Genome skimming approach reveals the gene arrangements in the chloroplast genomes of the highly endangered Crocus L. species: Crocus istanbulensis (B.Mathew) Rukšāns

Crocus istanbulensis (B.Mathew) Rukšāns is one of the most endangered Crocus species in the world and has an extremely limited distribution range in Istanbul. Our recent field work indicates that no more than one hundred individuals remain in the wild. In the present study, we used genome skimming to determine the complete chloroplast (cp) genome sequences of six C. istanbulensis individuals collected from the locus classicus. The cp genome of C. istanbulensis has 151,199 base pairs (bp), with a large single-copy (LSC) (81,197 bp), small single copy (SSC) (17,524 bp) and two inverted repeat (IR) regions of 26,236 bp each. The cp genome contains 132 genes, of which 86 are protein-coding (PCGs), 8 are rRNA and 38 are tRNA genes. Most of the repeats are found in intergenic spacers of Crocus species. Mononucleotide repeats were most abundant, accounting for over 80% of total repeats. The cp genome contained four palindrome repeats and one forward repeat. Comparative analyses among other Iridaceae species identified one inversion in the terminal positions of LSC region and three different gene (psbA, rps3 and rpl22) arrangements in C. istanbulensis that were not reported previously. To measure selective pressure in the exons of chloroplast coding sequences, we performed a sequence analysis of plastome-encoded genes. A total of seven genes (accD, rpoC2, psbK, rps12, ccsA, clpP and ycf2) were detected under positive selection in the cp genome. Alignment-free sequence comparison showed an extremely low sequence diversity across naturally occurring C. istanbulensis specimens. All six sequenced individuals shared the same cp haplotype. In summary, this study will aid further research on the molecular evolution and development of ex situ conservation strategies of C. istanbulensis.

OUP accepted manuscript

Despite recent methodology and reference database improvements for taxonomic profiling tools, metagenomic assembly and genomic binning remain important pillars of metagenomic analysis workflows. In case reference information is lacking, genomic binning is considered to be a state-of-the-art method in mixed culture metagenomic data analysis. In this light, our previously published tool BusyBee Web implements a composition-based binning method efficient enough to function as a rapid online utility. Handling assembled contigs and long nanopore generated reads alike, the webserver provides a wide range of supplementary annotations and visualizations. Half a decade after the initial publication, we revisited existing functionality, added comprehensive visualizations, and increased the number of data analysis customization options for further experimentation. The webserver now allows for visualization-supported differential analysis of samples, which is computationally expensive and typically only performed in coverage-based binning methods. Further, users may now optionally check their uploaded samples for plasmid sequences using PLSDB as a reference database. Lastly, a new application programming interface with a supporting python package was implemented, to allow power users fully automated access to the resource and integration into existing workflows. The webserver is freely available under: https://www.ccb.uni-saarland.de/busybee.