Correction of Persistent Errors in Arabidopsis Reference Mitochondrial Genomes

Comparative analysis of the complete mitochondrial genomes of Firmiana danxiaensis and F. kwangsiensis (Malvaceae), two endangered Firmiana species in China

MAIN CONCLUSION

We reported the mitogenomes of F. danxiaensis and F. kwangsiensis for the first time. Mitogenome size and structure differ significantly between them. Firmiana danxiaensis and F. kwangsiensis belong to the Firmiana genus and are distributed in the Danxia and Karst regions of southern China. Both species have been designated as endangered. Currently, the chloroplast genomes of F. danxiaensis and F. kwangsiensis have been sequenced, but the mitochondrial genome (mitogenome) of these two species has not been reported. To further understand the mitogenome characteristics, evolution, and phylogeny of F. danxiaensis and F. kwangsiensis, we assembled the mitogenomes of these two species based on a combination of Illumina and Nanopore sequencing methods. The mitogenome of F. danxiaensis exhibits a branching structure consisting of nine circular molecules with a total length of 938,890 bp, while the F. kwangsiensis has a circular structure with a length of 736,334 bp. Compared to F. kwangsiensis, F. danxiaensis has more tRNA genes, SSRs, tandem repeats, and dispersed repeats, while the codon use patterns are similar in these two species. There were 24 and 23 homologous sequences between mitogenome and chloroplast genome of F. danxiaensis and F. kwangsiensis, accounting for 0.37% and 0.49% of the mitogenome, respectively. In addition, the Ka/Ks ratio and the nucleic acid diversity analysis revealed that most of the mitochondria protein-coding genes in F. danxiaensis and F. kwangsiensis are highly conserved and may have undergone purifying selection. Furthermore, the collinear and comparative analysis showed that extensive genomic rearrangement events existed among the Malvaceae species. Lastly, a phylogenetic tree based on shared mitochondrial PCGs of 29 species revealed that F. danxiaensis and F. kwangsiensis form a sister group with high support values. Overall, the current study reports two mitogenomes (F. danxiaensis and F. kwangsiensis) in the Firmiana genus for the first time, which will help enhance comprehension of the mitogenome evolutionary patterns within Firmiana and promote the evolutionary and comparative genomic analyses within Malvaceae species.

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (∼15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

The molecular basis and evolution of the organellar RNA editosome by complementary DYW deaminases in seed plants

The DYW deaminase domain catalyzes the conversion of cytidines (C) to uridines (U) in RNA editing of plant organelles. While the DYW subgroup contains a complete DYW deaminase domain at the C-terminus, the E2 and E+ subgroups rely on complementary deaminases, in which catalytic activity depends on interactions with short DYW proteins, such as DYW1, DYW2, and MITOCHONDRIAL EDITING FACTOR 8 (MEF8)/MITOCHONDRIAL EDITING FACTOR 8 SIMILAR (MEF8S). Although orthogonal RNA editing in bacteria by a DYW subgroup pentatricopeptide repeat (PPR) has been reported, attempts to activate the DYW deaminase through molecular complementation in bacteria have been unsuccessful, leaving its molecular basis unresolved. In this study, we reconstituted the simplest editosome in Escherichia coli, composed of PPR56PPRE1E2-CRR4PG and DYW1 alone. Systematical mutational analysis of the PG-box of CHLORORESPIRATORY REDUCTION 4 (CRR4) in bacteria and in planta revealed the critical role of serine, isoleucine, and phenylalanine residues in DYW deaminase complementation and catalysis. CRR4-like PPR proteins, termed the "PG-type" characterized by the PG-box with these 3 key amino acid residues at the C-terminus, are minor in angiosperms but constitute one of the major subgroups in gymnosperms. Putative orthologs of Arabidopsis thaliana DYW1 are present in limited angiosperm species, suggesting that in other species, other short DYW proteins serve as the interaction partners for PG-type PPR proteins. Our findings reveal a minimal functional editosome module, shedding light on the conserved and diverse mechanisms of RNA editing in plant organelles.

The high-quality chromosome-level genome assembly of Dracocephalum rupestre Hance

Dracocephalum rupestre Hance is China's traditional herbal medicine in the family Labiatae with numerous health benefits, including anti-inflammatory, antiviral and anti-tumor activities. However, the genus Dracocephalum has no reference genome currently, which restricts the research on the breeding, cultivation and exploration of medicinal properties in D. rupestre. Thus, we present the high-quality chromosome-level genome assembly of D. rupestre using a combination of Pacbio HiFi sequencing and Hi-C scaffolding technologies. The final genome was 435.45 Mb with a contig N50 of 49.83 Mb and a scaffold N50 of 59.06 Mb. The assembled sequences were anchored to 7 chromosomes with an integration efficiency of 96.96%. Furthermore, we predicted 25,865 protein-coding genes, 98.23% of which were functionally annotated. These results offer valuable resources for understanding the genetic basis of the unique phenotypes of D. rupestre and will facilitate further study of the functional genomics of this species.

The maize mTERF18 regulates transcriptional termination of the mitochondrial nad6 gene and is essential for kernel development

Mitochondria are semi-autonomous organelles present in eukaryotic cells, containing their own genome and transcriptional machinery. However, their functions are intricately linked to proteins encoded by the nuclear genome. Mitochondrial transcription termination factors (mTERFs) are nucleic acid-binding proteins involved in RNA splicing and transcription termination within plant mitochondria and chloroplasts. Despite their recognized importance, the specific roles of mTERF proteins in maize remain largely unexplored. Here, we clone and functionally characterize the maize mTERF18 gene. Our findings reveal that mTERF18 mutations lead to severely undifferentiated embryos, resulting in abortive phenotypes. Early kernel exhibits abnormal basal endosperm transfer layer and a significant reduction in both starch and protein accumulation in mterf18. We identify the mTERF18 gene through mapping-based cloning and validate this gene through allelic tests. mTERF18 is widely expressed across various maize tissues and encodes a highly conserved mitochondrial protein. Transcriptome data reveal that mTERF18 mutations disrupt transcriptional termination of the nad6 gene, leading to undetectable levels of Nad6 protein and reduced complex I assembly and activity. Furthermore, transmission electron microscopy observation of mterf18 endosperm uncover severe mitochondrial defects. Collectively, these findings highlight the critical role of mTERF18 in mitochondrial gene transcription termination and its pivotal impact on maize kernel development.

Near-complete assembly and comprehensive annotation of the wheat Chinese Spring genome

A complete reference genome assembly is crucial for biological research and genetic improvement. Owing to its large size and highly repetitive nature, there are numerous gaps in the globally used wheat Chinese Spring (CS) genome assembly. In this study, we generated a 14.46 Gb near-complete assembly of the CS genome, with a contig N50 of over 266 Mb and an overall base accuracy of 99.9963%. Among the 290 gaps that remained (26, 257, and 7 gaps from the A, B, and D subgenomes, respectively), 278 were extremely high-copy tandem repeats, whereas the remaining 12 were transposable-element-associated gaps. Four chromosome assemblies were completely gap-free, including chr1D, chr3D, chr4D, and chr5D. Extensive annotation of the near-complete genome revealed 151 405 high-confidence genes, of which 59 180 were newly annotated, including 7602 newly assembled genes. Except for the centromere of chr1B, which has a gap associated with superlong GAA repeat arrays, the centromeric sequences of all of the remaining 20 chromosomes were completely assembled. Our near-complete assembly revealed that the extent of tandem repeats, such as simple-sequence repeats, was highly uneven among different subgenomes. Similarly, the repeat compositions of the centromeres also varied among the three subgenomes. With the genome sequences of all six types of seed storage proteins (SSPs) fully assembled, the expression of ω-gliadin was found to be contributed entirely by the B subgenome, whereas the expression of the other five types of SSPs was most abundant from the D subgenome. The near-complete CS genome will serve as a valuable resource for genomic and functional genomic research and breeding of wheat as well as its related species.

Flipping the switch on some of the slowest mutating genomes: Direct measurements of plant mitochondrial and plastid mutation rates in msh1 mutants

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

Comparative analysis of mitochondrial genomes of Stemona tuberosa lour. reveals heterogeneity in structure, synteny, intercellular gene transfer, and RNA editing

BACKGROUND

Stemona tuberosa, a vital species in traditional Chinese medicine, has been extensively cultivated and utilized within its natural distribution over the past decades. While the chloroplast genome of S. tuberosa has been characterized, its mitochondrial genome (mitogenome) remains unexplored.

RESULTS

This paper details the assembly of the complete S. tuberosa mitogenome, achieved through the integration of Illumina and Nanopore sequencing technologies. The assembled mitogenome is 605,873 bp in size with a GC content of 45.63%. It comprises 66 genes, including 38 protein-coding genes, 25 tRNA genes, and 3 rRNA genes. Our analysis delved into codon usage, sequence repeats, and RNA editing within the mitogenome. Additionally, we conducted a phylogenetic analysis involving S. tuberosa and 17 other taxa to clarify its evolutionary and taxonomic status. This study provides a crucial genetic resource for evolutionary research within the genus Stemona and other related genera in the Stemonaceae family.

CONCLUSION

Our study provides the inaugural comprehensive analysis of the mitochondrial genome of S. tuberosa, revealing its unique multi-branched structure. Through our investigation of codon usage, sequence repeats, and RNA editing within the mitogenome, coupled with a phylogenetic analysis involving S. tuberosa and 17 other taxa, we have elucidated its evolutionary and taxonomic status. These investigations provide a crucial genetic resource for evolutionary research within the genus Stemona and other related genera in the Stemonaceae family.

The mitochondrial genome of Carex pseudochinensis H. Lév. & Vaniot, an endemic sedge in Korea

Carex pseudochinensis H. Lév. & Vaniot is an endemic species in Korea and is included in the clade of section Paludosae in the recent classification system. We present the complete mitochondrial genome sequence of C. pseudochinensis based on the POLAP pipeline with both long- and short-read sequences. The mitochondrial genome is 997,628 bp in length, containing two large regions of 536.94 and 419.04 kbp, respectively, and a pair of direct repeat regions of about 20.25 kbp. The genome contains 57 genes, including 31 protein-coding genes, 20 tRNAs, and 6 rRNAs. Phylogenetic analysis based on mitochondrial proteomes, including those from ten species of related taxa, confirmed a close phylogenetic relationship between C. breviculmis and C. pseudochinensis.

Characterization, comparative phylogenetic, and gene transfer studies of Rosa chinensis var. spontanea organelle genomes

Rosa chinensis var. spontanea, an important horticultural species, holds high ornamental value. To investigate the phylogenetic relationships and organelle genome data of R. chinensis var. spontanea and its related species, the organelle genomes were constructed and sequenced in this study. The complete mitochondrial genome consists of a linear DNA molecule with a length of 270,143 base pairs and contains 55 genes. The whole chloroplast genome harbors 131 annotated genes. There are 39 and 54 simple sequence repeats (SSRs) in the mitochondrial and chloroplast genomes of R. chinensis var. spontanea, respectively. The phylogenetic tree of the mitochondrial genome indicates that R. chinensis var. spontanea is related to R. chinensis. The codon usage biases of one mitochondrial gene and nine chloroplast genes in Rosaceae were affected by mutation, while those of other genes have undergone selection. Overall, our findings suggest that the information within the organelle genomes can be utilized for constructing molecular markers, conducting genetic evolution research, and interpreting the phylogenetics of plants. The results of this study will contribute to a more in - depth investigation of the population genetics and evolution of Rosaceae.

Chromosome-level genome assembly of the tetraploid medicinal and natural dye plant Persicaria tinctoria

Persicaria tinctoria (2n = 40) is an important traditional medicinal plant and natural dye source within the genus Persicaria. P. tinctoria has been utilized for its antibacterial, antiviral, anti-inflammatory, and tumor treatment properties. Additionally, it has served as a natural blue dye for thousands of years worldwide, and continues to be used in countries such as China and Japan. Here, we assembled a tetraploid chromosome-scale genome of P. tinctoria, organized into two subgenomes: subgenome A, which contains 10 pseudochromosomes with a genome size of 888.67 Mb and a scaffold N50 of 90.56 Mb, and subgenome B, which also comprises 10 pseudochromosomes with a genome size of 771.58 Mb and a scaffold N50 of 76.84 Mb. Repeat sequences constitute 77.9% of the genome. A total of 76,742 high-confidence protein-coding genes were annotated, with 94.28% of these genes assigned functional annotations. This high-quality genome assembly of P. tinctoria will provide valuable genomic resources for studying the biosynthesis and evolution of indigoids in indigo plants, as well as for further research on the Polygonaceae family.

Species-wide inventory of Arabidopsis thaliana organellar variation reveals ample phenotypic variation for photosynthetic performance

Efforts to improve photosynthetic performance are increasingly employing natural genetic variation. However, genetic variation in the organellar genomes (plasmotypes) is often disregarded due to the difficulty of studying the plasmotypes and the lack of evidence that this is a worthwhile investment. Here, we systematically phenotyped plasmotype diversity using Arabidopsis thaliana as a model species. A reanalysis of whole-genome resequencing data of 1,541 representative accessions shows that the genetic diversity among the mitochondrial genomes is eight times lower than among the chloroplast genomes. Plasmotype diversity of the accessions divides the species into two major phylogenetic clusters, within which highly divergent subclusters are distinguished. We combined plasmotypes from 60 A. thaliana accessions with the nuclear genomes (nucleotypes) of four A. thaliana accessions to create a panel of 232 cytonuclear genotypes (cybrids). The cybrid plants were grown in a range of different light and temperature conditions and phenotyped using high-throughput phenotyping platforms. Analysis of the phenotypes showed that several plasmotypes alone or in interaction with the nucleotypes have significant effects on photosynthesis and that the effects are highly dependent on the environment. Moreover, we introduce Plasmotype Association Studies (PAS) as a method to reveal plasmotypic effects. Within A. thaliana, several organellar variants can influence photosynthetic phenotypes, which emphasizes the valuable role this variation has on improving photosynthetic performance. The increasing feasibility of producing cybrids in various species calls for further research into how these phenotypes may support breeding goals in crop species.

The landscape of Arabidopsis tRNA aminoacylation

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

Comparative analysis of mitochondrial genomes of invasive weed Mikania micrantha and its indigenous congener Mikania cordata

Mikania micrantha and Mikania cordata are two distinct species in China. The former is notorious as one of the top 100 worst invasive species, whereas the latter is an indigenous species harmless to native plants or the environment. They form an ideal congener pair for comparative studies aimed at deeply understanding the invasion mechanisms of the exotic weed. In this study, we have assembled and annotated the mitogenomes of both species using Illumina and PacBio sequencing data and compared their characteristic differences. The complete mitogenome of M. micrantha is a double-stranded DNA with a length of 336,564 bp, while the mitogenome of M. cordata exhibits a branching structure, consisting of two small circular molecules and six linear molecules, with a combined length totaling 335,444 bp. Compared to M. cordata, M. micrantha has less SSRs, tandem repeats, dispersed repeats, mitochondrial protein coding genes (PCGs). The two plants show similar codon usage patterns. This comparative study has revealed the structure and function of the mitogenomes of the two species and laid a solid foundation for investigating the effects of gene loss and duplication on the development of invasive traits in M. micrantha.

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.

PCIS1, Encoded by a Pentatricopeptide Protein Co-expressed Gene, Is Required for Splicing of Three Mitochondrial nad Transcripts in Angiosperms

Group II introns are large catalytic RNAs, which reside mainly within genes encoding respiratory complex I (CI) subunits in angiosperms' mitochondria. Genetic and biochemical analyses led to the identification of many nuclear-encoded factors that facilitate the splicing of the degenerated organellar introns in plants. Here, we describe the analysis of the pentatricopeptide repeat (PPR) co-expressed intron splicing-1 (PCIS1) factor, which was identified in silico by its co-expression pattern with many PPR proteins. PCIS1 is well conserved in land plants but has no sequence similarity with any known protein motifs. PCIS1 mutant lines are arrested in embryogenesis and can be maintained by the temporal expression of the gene under the embryo-specific ABI3 promoter. The pABI3::PCIS1 mutant plants display low germination and stunted growth phenotypes. RNA-sequencing and quantitative RT-PCR analyses of wild-type and mutant plants indicated that PCIS1 is a novel splicing cofactor that is pivotal for the maturation of several nad transcripts in Arabidopsis mitochondria. These phenotypes are tightly associated with respiratory CI defects and altered plant growth. Our data further emphasize the key roles of nuclear-encoded cofactors that regulate the maturation and expression of mitochondrial transcripts for the biogenesis of the oxidative phosphorylation system, and hence for plant physiology. The discovery of novel splicing factors other than typical RNA-binding proteins suggests further complexity of splicing mechanisms in plant mitochondria.

Systematic analysis of bZIP gene family in Suaeda australis reveal their roles under salt stress

BACKGROUND

Suaeda australis is one of typical halophyte owing to high levels of salt tolerance. In addition, the bZIP gene family assumes pivotal functions in response to salt stress. However, there are little reports available regarding the bZIP gene family in S. australis.

RESULTS

In this study, we successfully screened 44 bZIP genes within S. australis genome. Subsequently, we conducted an extensive analysis, encompassing investigations into chromosome location, gene structure, phylogenetic relationship, promoter region, conserved motif, and gene expression profile. The 44 bZIP genes were categorized into 12 distinct groups, exhibiting an uneven distribution among the 9 chromosomes of S. australis chromosomes, but one member (Sau23745) was mapped on unanchored scaffolds. Examination of cis-regulatory elements revealed that bZIP promoters were closely related to anaerobic induction, transcription start, and light responsiveness. Comparative transcriptome analysis between ST1 and ST2 samples identified 2,434 DEGs, which were significantly enriched in some primary biological pathways related to salt response-regulating signaling based on GO and KEGG enrichment analysis. Expression patterns analyses clearly discovered the role of several differently expressed SabZIPs, including Sau08107, Sau08911, Sau11415, Sau16575, and Sau19276, which showed higher expression levels in higher salt concentration than low concentration and a response to salt stress. These expression patterns were corroborated through RT-qPCR analysis. The six differentially expressed SabZIP genes, all localized in the nucleus, exhibited positive regulation involved in the salt stress response. SabZIP14, SabZIP26, and SabZIP36 proteins could bind to the promoter region of downstream salt stress-related genes and activate their expressions.

CONCLUSIONS

Our findings offer valuable insights into the evolutionary trajectory of the bZIP gene family in S. australis and shed light on their roles in responding to salt stress. In addition to fundamental genomic information, these results would serve as a foundational framework for future investigations into the regulation of salt stress responses in S. australis.

Plant mitochondrial genome map (PMGmap): A software tool for the comprehensive visualization of coding, noncoding and genome features of plant mitochondrial genomes

Tools for visualizing genomes are essential for investigating genomic features and their interactions. Currently, tools designed originally for animal mitogenomes and plant plastomes are used to visualize the mitogens of plants but cannot accurately display features specific to plant mitogenomes, such as nonlinear exon arrangement for genes, the prevalence of functional noncoding features and complex chromosomal architecture. To address these problems, a software package, plant mitochondrial genome map (PMGmap), was developed using the Python programming language. PMGmap can draw genes at exon levels; draw cis- and trans-splicing gene maps, noncoding features and repetitive sequences; and scale genic regions by using the scaling of the genic regions on the mitogenome (SAGM) algorithm. It can also draw multiple chromosomes simultaneously. Compared with other state-of-the-art tools, PMGmap showed better performance in visualizing 405 plant mitogenomes, showing potential as an invaluable tool for plant mitogenome research. The web and container versions and the source code of PMGmap can be accessed through the following link: http://www.1kmpg.cn/pmgmap.

Systematic characterization of CsbZIP transcription factors in Camelina sativa and functional analysis of CsbZIP-A12 mediating regulation of unsaturated fatty acid-enriched oil biosynthesis

The basic leucine zipper (bZIP) transcription factors (TFs) function importantly in numerous life processes in plants. However, bZIP members and their biological roles remain unknown in Camelina sativa, a worldwide promising oil crop. Here, 220 CsbZIP proteins were identified in camelina and classified into thirteen groups. Two and 347 pairs of tandem and segmental duplication genes were detected to be underwent purification selection, with segmental duplication as the main driven-force of CsbZIP gene family expansion. Most CsbZIP genes displayed a tissue-specific expression pattern. Particularly, CsbZIP-A12 significantly positively correlated with many FA/oil biosynthesis-related genes, indicating CsbZIP-A12 may regulate lipid biosynthesis. Notably, yeast one-hybrid (Y1H), β-Glucuronidase (GUS), dual-luciferase (LUC) and EMSA assays evidenced that CsbZIP-A12 located in nucleus interacted with the promoters of CsSAD2-3 and CsFAD3-3 genes responsible for unsaturated fatty acid (UFA) synthesis, thus activating their transcriptions. Overexpression of CsbZIP-A12 led to an increase of total lipid by 3.275 % compared to the control, followed with oleic and α-linolenic acid levels enhanced by 3.4 % and 5.195 %, and up-regulated the expressions of CsSAD2-3, CsFAD3-3 and CsPDAT2-3 in camelina seeds. Furthermore, heterogeneous expression of CsbZIP-A12 significantly up-regulated the expressions of NtSAD2, NtFAD3 and NtPDAT genes in tobacco plants, thereby improving the levels of total lipids and UFAs in both leaves and seeds without negative effects on other agronomic traits. Together, our findings suggest that CsbZIP-A12 upregulates FA/oil biosynthesis by activating CsSAD2-3 and CsFAD3-3 as well as possible other related genes. These data lay a foundation for further functional analyses of CsbZIPs, providing new insights into the TF-based lipid metabolic engineering to increase vegetable oil yield and health-beneficial quality in oilseeds.

Assembly and comparative analysis of the complete mitochondrial genome of Brassica rapa var. Purpuraria

BACKGROUND

Purple flowering stalk (Brassica rapa var. purpuraria) is a widely cultivated plant with high nutritional and medicinal value and exhibiting strong adaptability during growing. Mitochondrial (mt) play important role in plant cells for energy production, developing with an independent genetic system. Therefore, it is meaningful to assemble and annotate the functions for the mt genome of plants independently. Though there have been several reports referring the mt genome of in Brassica species, the genome of mt in B. rapa var. purpuraria and its functional gene variations when compared to its closely related species has not yet been addressed.

RESULTS

The mt genome of B. rapa var. purpuraria was assembled through the Illumina and Nanopore sequencing platforms, which revealed a length of 219,775 bp with a typical circular structure. The base composition of the whole B. rapa var. purpuraria mt genome revealed A (27.45%), T (27.31%), C (22.91%), and G (22.32%). 59 functional genes, composing of 33 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes, were annotated. The sequence repeats, codon usage, RNA editing, nucleotide diversity and gene transfer between the cp genome and mt genome were examined in the B. rapa var. purpuraria mt genome. Phylogenetic analysis show that B. rapa var. Purpuraria was closely related to B. rapa subsp. Oleifera and B. juncea. Ka/Ks analysis reflected that most of the PCGs in the B. rapa var. Purpuraria were negatively selected, illustrating that those mt genes were conserved during evolution.

CONCLUSIONS

The results of our findings provide valuable information on the B.rapa var. Purpuraria genome, which might facilitate molecular breeding, genetic variation and evolutionary researches for Brassica species in the future.

Root Primordium Defective 1 Encodes an Essential PORR Protein Required for the Splicing of Mitochondria-Encoded Group II Introns and for Respiratory Complex I Biogenesis

Cellular respiration involves complex organellar metabolic activities that are pivotal for plant growth and development. Mitochondria contain their own genetic system (mitogenome, mtDNA), which encodes key elements of the respiratory machinery. Plant mtDNAs are notably larger than their counterparts in Animalia, with complex genome organization and gene expression characteristics. The maturation of the plant mitochondrial transcripts involves extensive RNA editing, trimming and splicing events. These essential processing steps rely on the activities of numerous nuclear-encoded cofactors, which may also play key regulatory roles in mitochondrial biogenesis and function and hence in plant physiology. Proteins that harbor the plant organelle RNA recognition (PORR) domain are represented in a small gene family in plants. Several PORR members, including WTF1, WTF9 and LEFKOTHEA, are known to act in the splicing of organellar group II introns in angiosperms. The AT4G33495 gene locus encodes an essential PORR protein in Arabidopsis, termed ROOT PRIMORDIUM DEFECTIVE 1 (RPD1). A null mutation of At.RPD1 causes arrest in early embryogenesis, while the missense mutant lines, rpd1.1 and rpd1.2, exhibit a strong impairment in root development and retarded growth phenotypes, especially under high-temperature conditions. Here, we further show that RPD1 functions in the splicing of introns that reside in the coding regions of various complex I (CI) subunits (i.e. nad2, nad4, nad5 and nad7), as well as in the maturation of the ribosomal rps3 pre-RNA in Arabidopsis mitochondria. The altered growth and developmental phenotypes and modified respiration activities are tightly correlated with respiratory chain CI defects in rpd1 mutants.

Polly: An R package for genotyping microsatellites and detecting highly polymorphic DNA markers from short-read data

Highly polymorphic markers, such as microsatellites, are invaluable for the study of natural populations. However, contemporary methods for genotyping highly polymorphic variants have serious drawbacks that impede their efficiency. We created Polly, an R package with C++ source code that uses Illumina short-read data to genotype microsatellites, detect highly polymorphic variants and identify clusters of highly polymorphic SNPs, indels and microsatellites. We tested Polly on short-read data from Xiphophorus birchmanni (Teleostei: Poeciliidae) and Arabidopsis thaliana, finding it to be efficient and accurate both for microsatellite genotyping and polymorphic marker detection. This program can be applied to any diploid population for which there exists short-read data and at least one scaffolded reference genome.

Mitochondrial genome variation and intergenomic sequence transfers in Hevea species

Among the Hevea species, rubber tree (Hevea brasiliensis) is the most important source of natural rubber. In previous studies, we sequenced the complete nuclear and chloroplast genomes of Hevea species, providing an invaluable resource for studying their phylogeny, disease resistance, and breeding. However, given that plant mitochondrial genomes are more complex and more difficult to assemble than that of the other organelles, little is known about their mitochondrial genome, which limits the comprehensive understanding of Hevea genomic evolution. In this study, we sequenced and assembled the mitochondrial genomes of four Hevea species. The four mitochondrial genomes had consistent GC contents, codon usages and AT skews. However, there were significant differences in the genome lengths and sequence repeats. Specifically, the circular mitochondrial genomes of the four Hevea species ranged from 935,732 to 1,402,206 bp, with 34-35 unique protein-coding genes, 35-38 tRNA genes, and 6-13 rRNA genes. In addition, there were 17,294-46,552 bp intergenomic transfer fragments between the chloroplast and mitochondrial genomes, consisting of eight intact genes (psaA, rrn16S, tRNA-Val, rrn5S, rrn4.5S, tRNA-Arg, tRNA-Asp, and tRNA-Asn), intergenic spacer regions and partial gene sequences. The evolutionary position of Hevea species, crucial for understanding its adaptive strategies and relation to other species, was verified by phylogenetic analysis based on the protein-coding genes in the mitochondrial genomes of 21 Malpighiales species. The findings from this study not only provide valuable insights into the structure and evolution of the Hevea mitochondrial genome but also lay the foundation for further molecular, evolutionary studies, and genomic breeding studies on rubber tree and other Hevea species, thereby potentially informing conservation and utilization strategies.

Mitochondrial atp1 mRNA knockdown by a custom-designed pentatricopeptide repeat protein alters ATP synthase

Spontaneous mutations are rare in mitochondria and the lack of mitochondrial transformation methods has hindered genetic analyses. We show that a custom-designed RNA-binding pentatricopeptide repeat (PPR) protein binds and specifically induces cleavage of ATP synthase subunit1 (atp1) mRNA in mitochondria, significantly decreasing the abundance of the Atp1 protein and the assembled F1Fo ATP synthase in Arabidopsis (Arabidopsis thaliana). The transformed plants are characterized by delayed vegetative growth and reduced fertility. Five-fold depletion of Atp1 level was accompanied by a decrease in abundance of other ATP synthase subunits and lowered ATP synthesis rate of isolated mitochondria, but no change to mitochondrial electron transport chain complexes, adenylates, or energy charge in planta. Transcripts for amino acid transport and a variety of stress response processes were differentially expressed in lines containing the PPR protein, indicating changes to achieve cellular homeostasis when ATP synthase was highly depleted. Leaves of ATP synthase-depleted lines showed higher respiratory rates and elevated steady-state levels of numerous amino acids, most notably of the serine family. The results show the value of using custom-designed PPR proteins to influence the expression of specific mitochondrial transcripts to carry out reverse genetic studies on mitochondrial gene functions and the consequences of ATP synthase depletion on cellular functions in Arabidopsis.

PMAT: an efficient plant mitogenome assembly toolkit using low-coverage HiFi sequencing data

Complete mitochondrial genomes (mitogenomes) of plants are valuable resources for nucleocytoplasmic interactions, plant evolution, and plant cytoplasmic male sterile line breeding. However, the complete assembly of plant mitogenomes is challenging due to frequent recombination events and horizontal gene transfers. Previous studies have adopted Illumina, PacBio, and Nanopore sequencing data to assemble plant mitogenomes, but the poor assembly completeness, low sequencing accuracy, and high cost limit the sampling capacity. Here, we present an efficient assembly toolkit (PMAT) for de novo assembly of plant mitogenomes using low-coverage HiFi sequencing data. PMAT has been applied to the de novo assembly of 13 broadly representative plant mitogenomes, outperforming existing organelle genome assemblers in terms of assembly accuracy and completeness. By evaluating the assembly of plant mitogenomes from different sequencing data, it was confirmed that PMAT only requires 1× HiFi sequencing data to obtain a complete plant mitogenome. The source code for PMAT is available at https://github.com/bichangwei/PMAT. The developed PMAT toolkit will indeed accelerate the understanding of evolutionary variation and breeding application of plant mitogenomes.

Detection and editing of the updated Arabidopsis plastid- and mitochondrial-encoded proteomes through PeptideAtlas

Arabidopsis (Arabidopsis thaliana) ecotype Col-0 has plastid and mitochondrial genomes encoding over 100 proteins. Public databases (e.g. Araport11) have redundancy and discrepancies in gene identifiers for these organelle-encoded proteins. RNA editing results in changes to specific amino acid residues or creation of start and stop codons for many of these proteins, but the impact of RNA editing at the protein level is largely unexplored due to the complexities of detection. Here, we assembled the nonredundant set of identifiers, their correct protein sequences, and 452 predicted nonsynonymous editing sites of which 56 are edited at lower frequency. We then determined accumulation of edited and/or unedited proteoforms by searching ∼259 million raw tandem MS spectra from ProteomeXchange, which is part of PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/). We identified all mitochondrial proteins and all except 3 plastid-encoded proteins (NdhG/Ndh6, PsbM, and Rps16), but no proteins predicted from the 4 ORFs were identified. We suggest that Rps16 and 3 of the ORFs are pseudogenes. Detection frequencies for each edit site and type of edit (e.g. S to L/F) were determined at the protein level, cross-referenced against the metadata (e.g. tissue), and evaluated for technical detection challenges. We detected 167 predicted edit sites at the proteome level. Minor frequency sites were edited at low frequency at the protein level except for cytochrome C biogenesis 382 at residue 124 (Ccb382-124). Major frequency sites (>50% editing of RNA) only accumulated in edited form (>98% to 100% edited) at the protein level, with the exception of Rpl5-22. We conclude that RNA editing for major editing sites is required for stable protein accumulation.

Assembly and characterization of the complete mitochondrial genome of Ventilago leiocarpa

KEY MESSAGE

We reported the mitochondrial genome of Ventilago leiocarpa for the first time. Two and one sites lead to the generation of stop and stat codon through editing were verified. Ventilago leiocarpa, a member of the Rhamnaceae family, is frequently utilized in traditional medicine due to the medicinal properties of its roots. In this study, we successfully assembled the mitogenome of V. leiocarpa using both BGI short reads and Nanopore long reads. This mitogenome has a total length of 331,839 bp. The annotated results showed 36 unique protein-coding, 16 tRNA and 3 rRNA genes in this mitogenome. Furthermore, we confirmed the presence of a branched structure through the utilization of long reads mapping, PCR amplification, and Sanger sequencing. Specifically, the ctg1 can form a single circular molecule or combine with ctg4 to form a linear molecule. Likewise, ctg2 can form a single circular molecule or can be connected to ctg4 to form a linear molecule. Subsequently, through a comparative analysis of the mitogenome and cpgenome sequences, we identified ten mitochondrial plastid sequences (MTPTs), including two complete protein-coding genes and five complete tRNA genes. The existence of MTPTs was verified by long reads. Colinear analysis showed that the mitogenomes of Rosales were highly divergent in structure. Finally, we identified 545 RNA editing sites involving 36 protein-coding genes by Deepred-mt. To validate our findings, we conducted PCR amplification and Sanger sequencing, which confirmed the generation of stop codons in atp9-223 and rps10-391, as well as the generation of a start codon in nad4L-2. This project reported the complex structure and RNA editing event of the V. Leiocarpa mitogenome, which will provide valuable information for the study of mitochondrial gene expression.

Assembly and analysis of the complete mitochondrial genome of the Chinese wild dwarf almond (Prunus tenella)

Background: The wild dwarf almond (Prunus tenella) is one of the national key grade II-protected wild plants in China. It is a relic deciduous forest species from the middle Eocene of the ancient Mediterranean Sea and is also known as a "living fossil of plants." It is distributed in Southeast Europe, West Asia, Central Asia, Siberia, and Xinjiang (Tacheng) and other areas of China. The plant grows on arid slopes, steppes, depressions, and valleys at an altitude of 1,200 m. The seeds of wild dwarf almonds are frost resistant and contain oil and bitter lentil glycosides, which possess medicinal value. Additionally, the seeds of wild dwarf almonds can be used as the original material for breeding new varieties of almonds and obtain ornamental flowers and trees. Results: The complete mitochondrial genome of P. tenella was sequenced and assembled using two sequencing platforms, namely, Illumina Novaseq6000 and Oxford Nanopore PromethION. The assembled genome was 452,158-bp long with a typical loop structure. The total number of A, T, C, and G bases in the genome was 122,066 (26.99%), 124,114 (27.45%), 103,285 (22.84%), and 102,693 (22.71%), respectively, with a GC content of 45.55%. A total of 63 unique genes, including 36 protein-coding genes, 24 tRNA genes, and 3 rRNA genes, were identified in the genome. Furthermore, codon usage, sequence duplication, RNA editing, and mitochondrial and chloroplast DNA fragment transfer events in the genome were analyzed. A phylogenetic tree was also constructed using 30 protein-coding genes that are common to the mitochondrial genomes of 24 species, which indicated that the genome of wild lentils is highly conserved with those of apples and pears belonging to Rosaceae. Conclusion: Assembly and annotation of the P. tenella mitochondrial genome provided comprehensive information about the mitochondrial genome of wild dwarf almonds, This study provides information on the mitochondrial genome of Prunus species and serves as a reference for further evolutionary studies on wild dwarf almonds.

Detection of the Arabidopsis Proteome and Its Post-translational Modifications and the Nature of the Unobserved (Dark) Proteome in PeptideAtlas

This study describes a new release of the Arabidopsis thaliana PeptideAtlas proteomics resource (build 2023-10) providing protein sequence coverage, matched mass spectrometry (MS) spectra, selected post-translational modifications (PTMs), and metadata. 70 million MS/MS spectra were matched to the Araport11 annotation, identifying ∼0.6 million unique peptides and 18,267 proteins at the highest confidence level and 3396 lower confidence proteins, together representing 78.6% of the predicted proteome. Additional identified proteins not predicted in Araport11 should be considered for the next Arabidopsis genome annotation. This release identified 5198 phosphorylated proteins, 668 ubiquitinated proteins, 3050 N-terminally acetylated proteins, and 864 lysine-acetylated proteins and mapped their PTM sites. MS support was lacking for 21.4% (5896 proteins) of the predicted Araport11 proteome: the "dark" proteome. This dark proteome is highly enriched for E3 ligases, transcription factors, and for certain (e.g., CLE, IDA, PSY) but not other (e.g., THIONIN, CAP) signaling peptides families. A machine learning model trained on RNA expression data and protein properties predicts the probability that proteins will be detected. The model aids in discovery of proteins with short half-life (e.g., SIG1,3 and ERF-VII TFs) and for developing strategies to identify the missing proteins. PeptideAtlas is linked to TAIR, tracks in JBrowse, and several other community proteomics resources.

The complete mitochondrial genome of the Kandelia obovata Sheue, H.Y.Liu & J.W.H.Yong (Rhizophoraceae)

Kandelia obovata Sheue, H.Y.Liu & J.W.H.Yong is one of the most cold-resistant true mangrove species, and it is widely distributed from the South China Sea to southern Japan. In the current study, the complete mitochondrial genome sequence of K. obovata was assembled using Illumina reads. It is the first mitochondrial genome of the Kandelia genus within the family Rhizophoraceae to be sequenced. The mitochondrial genome size is 312,146 bp with a total of 49 predicted genes, including 29 protein-coding genes, 17 transfer RNA genes, and 3 ribosomal RNA genes. The overall GC content of the genome is 41.87%. A phylogenetic tree constructed using nine complete mitochondrial genomes revealed that K. obovata is more closely related to Bruguiera species. This study enriches the plastid genome of Kandelia, furnishing valuable genetic insights for the investigation of evolutionary and population genetics in Kandelia and other mangrove species.

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

Introduction

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

Methods

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

Results and discussion

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

MOMS: A pipeline for scaffolding using multi-optical maps

Here, we report a new multi-optical maps scaffolder (MOMS) aiming at utilizing complementary information among optical maps labelled by distinct enzymes. This pipeline was designed for data structure organization, scaffolding by path traversal, gap-filling and molecule reuse of optical maps. Our testing showed that this pipeline has uncapped enzyme tolerance in scaffolding. This means that there are no inbuilt limits as to the number of maps generated by different enzymes that can be utilized by MOMS. For the genome assembly of the human GM12878 cell line, MOMS significantly improved the contiguity and completeness with an up to 144-fold increase of scaffold N50 compared with initial assemblies. Benchmarking on the genomes of human and O. sativa showed that MOMS is more effective and robust compared with other optical-map-based scaffolders. We believe this pipeline will contribute to high-fidelity chromosome assembly and chromosome-level evolutionary analysis.

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.

Multiple configurations of the plastid and mitochondrial genomes of Caragana spinosa

MAIN CONCLUSION

In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes. Caragana spinosa belongs to the Papilionoidea subfamily and has significant pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled and analyzed the plastome and mitogenome of C. spinosa using the Illumina and Nanopore DNA sequencing data. The plastome of C. spinosa was 129,995 bp belonging to the inverted repeat lacking clade (IRLC), which contained 77 protein-coding genes, 29 tRNA genes, and four rRNA genes. The mitogenome was 378,373 bp long and encoded 54 unique genes, including 33 protein-coding, three ribosomal RNA (rRNA), and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, alternative conformations mediated by one and four repetitive sequences in the plastome and mitogenome were identified and validated, respectively. The inverted repeat (PDR12, the 12th dispersed repeat sequence in C. spinosa plastome) of plastome mediating recombinant was conserved in the genus Caragana. Furthermore, we identified 14 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. A phylogenetic analysis of protein-coding genes extracted from the plastid and mitochondrial genomes revealed congruent topologies. Analyses of sequence divergence found one intergenic region, trnN-GUU-ycf1, exhibiting a high degree of variation, which can be used to develop novel molecular markers to distinguish the nine Caragana species accurately. This plastome and mitogenome of C. spinosa could provide critical information for the molecular breeding of C. spinosa and be used as a reference genome for other species of Caragana. In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes.

The complete mitochondrial genome of the biodiesel plant Jatropha curcas L

Jatropha curcas (Linnaeus, 1753) is a plant species in the order Malpighiales and the family Euphorbiaceae and is native to the tropical regions of America, such as Mexico and Argentina. Currently, this plant species inhabits tropical and subtropical regions of the world. Jatropha has been widely used as a biofuel plant to produce high-quality diesel engine fuel. In this study, the complete mitochondrial genome sequence of J. curcas was assembled into 561,839 bp circular nucleotides with a GC content of 44.6%. The mitochondrial genome of J. curcas comprises 33 known protein-coding genes, 22 tRNA genes, three rRNA genes, one ncRNA gene, and 85 open reading frame genes. Phylogenetic analysis showed this species is closely related to the castor bean (Ricinus communis).

An mTRAN-mRNA interaction mediates mitochondrial translation initiation in plants

Plant mitochondria represent the largest group of respiring organelles on the planet. Plant mitochondrial messenger RNAs (mRNAs) lack Shine-Dalgarno-like ribosome-binding sites, so it is unknown how plant mitoribosomes recognize mRNA. We show that "mitochondrial translation factors" mTRAN1 and mTRAN2 are land plant-specific proteins, required for normal mitochondrial respiration chain biogenesis. Our studies suggest that mTRANs are noncanonical pentatricopeptide repeat (PPR)-like RNA binding proteins of the mitoribosomal "small" subunit. We identified conserved Adenosine (A)/Uridine (U)-rich motifs in the 5' regions of plant mitochondrial mRNAs. mTRAN1 binds this motif, suggesting that it is a mitoribosome homing factor to identify mRNAs. We demonstrate that mTRANs are likely required for translation of all plant mitochondrial mRNAs. Plant mitochondrial translation initiation thus appears to use a protein-mRNA interaction that is divergent from bacteria or mammalian mitochondria.

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.

Nuclear DNA segments homologous to mitochondrial DNA are obstacles for detecting heteroplasmy in sugar beet (Beta vulgaris L.)

Heteroplasmy, the coexistence of multiple mitochondrial DNA (mtDNA) sequences in a cell, is well documented in plants. Next-generation sequencing technology (NGS) has made it feasible to sequence entire genomes. Thus, NGS has the potential to detect heteroplasmy; however, the methods and pitfalls in heteroplasmy detection have not been fully investigated and identified. One obstacle for heteroplasmy detection is the sequence homology between mitochondrial-, plastid-, and nuclear DNA, of which the influence of nuclear DNA segments homologous to mtDNA (numt) need to be minimized. To detect heteroplasmy, we first excluded nuclear DNA sequences of sugar beet (Beta vulgaris) line EL10 from the sugar beet mtDNA sequence. NGS reads were obtained from single plants of sugar beet lines NK-195BRmm-O and NK-291BRmm-O and mapped to the unexcluded mtDNA regions. More than 1000 sites exhibited intra-individual polymorphism as detected by genome browsing analysis. We focused on a 309-bp region where 12 intra-individual polymorphic sites were closely linked to each other. Although the existence of DNA molecules having variant alleles at the 12 sites was confirmed by PCR amplification from NK-195BRmm-O and NK-291BRmm-O, these variants were not always called by six variant-calling programs, suggesting that these programs are inappropriate for intra-individual polymorphism detection. When we changed the nuclear DNA reference, a numt absent from EL10 was found to include the 309-bp region. Genetic segregation of an F2 population from NK-195BRmm-O x NK-291BRmm-O supported the numt origin of the variant alleles. Using four references, we found that numt detection exhibited reference dependency, and extreme polymorphism of numts exists among sugar beet lines. One of the identified numts absent from EL10 is also associated with another intra-individual polymorphic site in NK-195mm-O. Our data suggest that polymorphism among numts is unexpectedly high within sugar beets, leading to confusion about the true degree of heteroplasmy.

A chromosome-scale and haplotype-resolved genome assembly of carnation (Dianthus caryophyllus) based on high-fidelity sequencing

Dianthus caryophyllus is an economic species often considered excellent cut flowers and is suitable for bouquets and gardens. Here, we assembled the haplotype-resolved genome of D. caryophyllus 'Aili' at the chromosome level for the first time. The total lengths of the two assembled haplotypes of carnation were 584.88 Mb for haplotype genome 1 (hap1) and 578.78 Mb for haplotype genome 2 (hap2), respectively. We predicted a total of 44,098 and 42,425 protein-coding genes, respectively. The remarkable structure variation was identified between two haplotypes. Moreover, we identified 403.80 Mb of transposable elements (TEs) in hap1, which accounted for 69.34% of the genome. In contrast, hap2 had 402.70 Mb of TEs, representing 69.61% of the genome. Long terminal repeats were the predominant transposable elements. Phylogenetic analysis showed that the species differentiation time between carnation and gypsophila was estimated to be ~54.43 MYA. The unique gene families of carnation genomes were identified in 'Aili' and previously published 'Francesco' and 'Scarlet Queen'. The assembled and annotated haplotype-resolved D. caryophyllus genome not only promises to facilitate molecular biology studies but also contributes to genome-level evolutionary studies.

Assembly and analysis of the mitochondrial genome of Prunella vulgaris

Prunella vulgaris (Lamiaceae) is widely distributed in Eurasia. Former studies have demonstrated that P. vulgaris has a wide range of pharmacological effects. Nevertheless, no complete P. vulgaris mitochondrial genome has been reported, which limits further understanding of the biology of P. vulgaris. Here, we assembled the first complete mitochondrial genome of P. vulgaris using a hybrid assembly strategy based on sequencing data from both Nanopore and Illumina platforms. Then, the mitochondrial genome of P. vulgaris was analyzed comprehensively in terms of gene content, codon preference, intercellular gene transfer, phylogeny, and RNA editing. The mitochondrial genome of P. vulgaris has two circular structures. It has a total length of 297, 777 bp, a GC content of 43.92%, and 29 unique protein-coding genes (PCGs). There are 76 simple sequence repeats (SSRs) in the mitochondrial genome, of which tetrameric accounts for a large percentage (43.4%). A comparative analysis between the mitochondrial and chloroplast genomes revealed that 36 homologous fragments exist in them, with a total length of 28, 895 bp. The phylogenetic analysis showed that P. vulgaris belongs to the Lamiales family Lamiaceae and P. vulgaris is closely related to Salvia miltiorrhiza. In addition, the mitochondrial genome sequences of seven species of Lamiaceae are unconservative in their alignments and undergo frequent genome reorganization. This work reports for the first time the complete mitochondrial genome of P. vulgaris, which provides useful genetic information for further Prunella studies.

Mapping the Arabidopsis thaliana proteome in PeptideAtlas and the nature of the unobserved (dark) proteome; strategies towards a complete proteome

This study describes a new release of the Arabidopsis thaliana PeptideAtlas proteomics resource providing protein sequence coverage, matched mass spectrometry (MS) spectra, selected PTMs, and metadata. 70 million MS/MS spectra were matched to the Araport11 annotation, identifying ∼0.6 million unique peptides and 18267 proteins at the highest confidence level and 3396 lower confidence proteins, together representing 78.6% of the predicted proteome. Additional identified proteins not predicted in Araport11 should be considered for building the next Arabidopsis genome annotation. This release identified 5198 phosphorylated proteins, 668 ubiquitinated proteins, 3050 N-terminally acetylated proteins and 864 lysine-acetylated proteins and mapped their PTM sites. MS support was lacking for 21.4% (5896 proteins) of the predicted Araport11 proteome - the 'dark' proteome. This dark proteome is highly enriched for certain ( e.g. CLE, CEP, IDA, PSY) but not other ( e.g. THIONIN, CAP,) signaling peptides families, E3 ligases, TFs, and other proteins with unfavorable physicochemical properties. A machine learning model trained on RNA expression data and protein properties predicts the probability for proteins to be detected. The model aids in discovery of proteins with short-half life ( e.g. SIG1,3 and ERF-VII TFs) and completing the proteome. PeptideAtlas is linked to TAIR, JBrowse, PPDB, SUBA, UniProtKB and Plant PTM Viewer.

MSP1 encodes an essential RNA-binding pentatricopeptide repeat factor required for nad1 maturation and complex I biogenesis in Arabidopsis mitochondria

Mitochondrial biogenesis relies on nuclearly encoded factors, which regulate the expression of the organellar-encoded genes. Pentatricopeptide repeat (PPR) proteins constitute a major gene family in angiosperms that are pivotal in many aspects of mitochondrial (mt)RNA metabolism (e.g. trimming, splicing, or stability). Here, we report the analysis of MITOCHONDRIA STABILITY/PROCESSING PPR FACTOR1 (MSP1, At4g20090), a canonical PPR protein that is necessary for mitochondrial functions and embryo development. Loss-of-function allele of MSP1 leads to seed abortion. Here, we employed an embryo-rescue method for the molecular characterization of msp1 mutants. Our analyses reveal that msp1 embryogenesis fails to proceed beyond the heart/torpedo stage as a consequence of a nad1 pre-RNA processing defect, resulting in the loss of respiratory complex I activity. Functional complementation confirmed that msp1 phenotypes result from a disruption of the MSP1 gene. In Arabidopsis, the maturation of nad1 involves the processing of three RNA fragments, nad1.1, nad1.2, and nad1.3. Based on biochemical analyses and mtRNA profiles of wild-type and msp1 plants, we concluded that MSP1 facilitates the generation of the 3' terminus of nad1.1 transcript, a prerequisite for nad1 exons a-b splicing. Our data substantiate the importance of mtRNA metabolism for the biogenesis of the respiratory system during early plant life.

Plant organellar RNA maturation

Plant organellar RNA metabolism is run by a multitude of nucleus-encoded RNA-binding proteins (RBPs) that control RNA stability, processing, and degradation. In chloroplasts and mitochondria, these post-transcriptional processes are vital for the production of a small number of essential components of the photosynthetic and respiratory machinery-and consequently for organellar biogenesis and plant survival. Many organellar RBPs have been functionally assigned to individual steps in RNA maturation, often specific to selected transcripts. While the catalog of factors identified is ever-growing, our knowledge of how they achieve their functions mechanistically is far from complete. This review summarizes the current knowledge of plant organellar RNA metabolism taking an RBP-centric approach and focusing on mechanistic aspects of RBP functions and the kinetics of the processes they are involved in.

Mitochondrial gene defects in Arabidopsis can broadly affect mitochondrial gene expression through copy number

How mitochondria regulate the expression of their genes is poorly understood, partly because methods have not been developed for stably transforming mitochondrial genomes. In recent years, the disruption of mitochondrial genes has been achieved in several plant species using mitochondria-localized TALEN (mitoTALEN). In this study, we attempted to disrupt the NADH dehydrogenase subunit7 (NAD7) gene, a subunit of respiratory chain complex I, in Arabidopsis (Arabidopsis thaliana) using the mitoTALEN method. In some of the transformants, disruption of NAD7 was accompanied by severe growth inhibition and lethality, suggesting that NAD7 has an essential function in Arabidopsis. In addition, the mitochondrial genome copy number and overall expression of genes encoding mitochondrial proteins were generally increased by nad7 knockout. Similar increases were also observed in mutants with decreased NAD7 transcripts and with dysfunctions of other mitochondrial respiratory complexes. In these mutants, the expression of nuclear genes involved in mitochondrial translation or protein transport was induced in sync with mitochondrial genes. Mitochondrial genome copy number was also partly regulated by the nuclear stress-responsive factors NAC domain containing protein 17 and Radical cell death 1. These findings suggest the existence of overall gene-expression control through mitochondrial genome copy number in Arabidopsis and that disruption of single mitochondrial genes can have additional broad consequences in both the nuclear and mitochondrial genomes.

Master graph: an essential integrated assembly model for the plant mitogenome based on a graph-based framework

Unlike the typical single circular structure of most animal mitochondrial genomes (mitogenome), the drastic structural variation of plant mitogenomes is a result of a mixture of molecules of various sizes and structures. Obtaining the full panoramic plant mitogenome is still considered a roadblock in evolutionary biology. In this study, we developed a graph-based sequence assembly toolkit (GSAT) to construct the pan-structural landscape of plant mitogenome with high-quality mitochondrial master graphs (MMGs) for model species including rice (Oryza sativa) and thale cress (Arabidopsis thaliana). The rice and thale cress MMGs have total lengths of 346 562 and 358 041 bp, including 9 and 6 contigs and 12 and 8 links, respectively, and could be further divided into 6 and 3 minimum master circles and 4 and 2 minimum secondary circles separately. The nuclear mitochondrial DNA segments (NUMTs) in thale cress strongly affected the frequency evaluation of the homologous structures in the mitogenome, while the effects of NUMTs in rice were relatively weak. The mitochondrial plastid DNA segments (MTPTs) in both species had no effects on the assessment of the MMGs. All potential recombinant structures were evaluated, and the findings revealed that all, except for nuclear-homologous structures, MMG structures are present at a much higher frequency than non-MMG structures are. Investigations of potential circular and linear molecules further supported multiple dominant structures in the mitogenomes and could be completely summarized in the MMG. Our study provided an efficient and accurate model for assembling and applying graph-based plant mitogenomes to assess their pan-structural variations.

Characterization and phylogenetic analysis of the complete mitochondrial genome sequence of Photinia serratifolia

Plant mitochondrial genomes (mitogenomes) are a valuable source of genetic information for a better understanding of phylogenetic relationships. However, no mitogenome of any species in the genus of Photinia has been reported. In this study, using NGS sequencing, we reported the mitogenome assembly and annotation of Photinia serratifolia, which is 473,579 bp in length, contains 38 protein-coding genes, 23 tRNAs, and 6 rRNAs, with 61 genes have no introns. The rps2 and rps11 genes are missing in the P. serratifolia mitogenome. Although there are more editing sites (488) in the P. serratifolia mitogenome than in most angiosperms, fewer editing types were found in the P. serratifolia mitogenome, showing a clear bias in RNA-editing. Phylogenetic analysis based on the mitogenomes of P. serratifolia and 8 other taxa of the Rosaceae family reflected the exact evolutionary and taxonomic status of P. serratifolia. However, Ka/Ks analysis revealed that 72.69% of the protein-coding genes in the P. serratifolia mitogenome had undergone negative selections, reflecting the importance of those genes in the P. serratifolia mitogenome. Collectively, these results will provide valuable information for the evolution of P. serratifolia and provide insight into the evolutionary relationships within Photinia and the Rosaceae family.

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

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

DEMETHYLATION REGULATOR 1 regulates DNA demethylation of the nuclear and mitochondrial genomes

Active DNA demethylation effectively modulates gene expression during plant development and in response to stress. However, little is known about the upstream regulatory factors that regulate DNA demethylation. We determined that the demethylation regulator 1 (demr1) mutant exhibits a distinct DNA methylation profile at selected loci queried by methylation-sensitive polymerase chain reaction and globally based on whole-genome bisulfite sequencing. Notably, the transcript levels of the DNA demethylase gene REPRESSOR OF SILENCING 1 (ROS1) were lower in the demr1 mutant. We established that DEMR1 directly binds to the ROS1 promoter in vivo and in vitro, and the methylation level in the DNA methylation monitoring sequence of ROS1 promoter decreased by 60% in the demr1 mutant. About 40% of the hyper-differentially methylated regions (DMRs) in the demr1 mutant were shared with the ros1-4 mutant. Genetic analysis indicated that DEMR1 acts upstream of ROS1 to positively regulate abscisic acid (ABA) signaling during seed germination and seedling establishment stages. Surprisingly, the loss of DEMR1 function also caused a rise in methylation levels of the mitochondrial genome, impaired mitochondrial structure and an early flowering phenotype. Together, our results show that DEMR1 is a novel regulator of DNA demethylation of both the nuclear and mitochondrial genomes in response to ABA and plant development in Arabidopsis.

The mitochondrial genome sequence of Abies alba Mill. reveals a high structural and combinatorial variation

BACKGROUND

Plant mitogenomes vary widely in size and genomic architecture. Although hundreds of plant mitogenomes of angiosperm species have already been sequence-characterized, only a few mitogenomes are available from gymnosperms. Silver fir (Abies alba) is an economically important gymnosperm species that is widely distributed in Europe and occupies a large range of environmental conditions. Reference sequences of the nuclear and chloroplast genome of A. alba are available, however, the mitogenome has not yet been assembled and studied.

RESULTS

Here, we used paired-end Illumina short reads generated from a single haploid megagametophyte in combination with PacBio long reads from high molecular weight DNA of needles to assemble the first mitogenome sequence of A. alba. Assembly and scaffolding resulted in 11 mitogenome scaffolds, with the largest scaffold being 0.25 Mbp long. Two of the scaffolds displayed a potential circular structure supported by PCR. The total size of the A. alba mitogenome was estimated at 1.43 Mbp, similar to the size (1.33 Mbp) of a draft assembly of the Abies firma mitogenome. In total, 53 distinct genes of known function were annotated in the A. alba mitogenome, comprising 41 protein-coding genes, nine tRNA, and three rRNA genes. The proportion of highly repetitive elements (REs) was 0.168. The mitogenome seems to have a complex and dynamic structure featured by high combinatorial variation, which was specifically confirmed by PCR for the contig with the highest mapping coverage. Comparative analysis of all sequenced mitogenomes of gymnosperms revealed a moderate, but significant positive correlation between mitogenome size and proportion of REs.

CONCLUSIONS

The A. alba mitogenome provides a basis for new comparative studies and will allow to answer important structural, phylogenetic and other evolutionary questions. Future long-read sequencing with higher coverage of the A. alba mitogenome will be the key to further resolve its physical structure. The observed positive correlation between mitogenome size and proportion of REs will be further validated once available mitogenomes of gymnosperms would become more numerous. To test whether a higher proportion of REs in a mitogenome leads to an increased recombination and higher structural complexity and variability is a prospective avenue for future research.

De Novo Hybrid Assembly of the Salvia miltiorrhiza Mitochondrial Genome Provides the First Evidence of the Multi-Chromosomal Mitochondrial DNA Structure of Salvia Species

Salvia miltiorrhiza has been an economically important medicinal plant. Previously, an S. miltiorrhiza mitochondrial genome (mitogenome) assembled from Illumina short reads, appearing to be a single circular molecule, has been published. Based on the recent reports on the plant mitogenome structure, we suspected that this conformation does not accurately represent the complexity of the S. miltiorrhiza mitogenome. In the current study, we assembled the mitogenome of S. miltiorrhiza using the PacBio and Illumina sequencing technologies. The primary structure of the mitogenome contained two mitochondrial chromosomes (MC1 and MC2), which corresponded to two major conformations, namely, Mac1 and Mac2, respectively. Using two approaches, including (1) long reads mapping and (2) polymerase chain reaction amplification followed by Sanger sequencing, we observed nine repeats that can mediate recombination. We predicted 55 genes, including 33 mitochondrial protein-coding genes (PCGs), 3 rRNA genes, and 19 tRNA genes. Repeat analysis identified 112 microsatellite repeats and 3 long-tandem repeats. Phylogenetic analysis using the 26 shared PCGs resulted in a tree that was congruent with the phylogeny of Lamiales species in the APG IV system. The analysis of mitochondrial plastid DNA (MTPT) identified 16 MTPTs in the mitogenome. Moreover, the analysis of nucleotide substitution rates in Lamiales showed that the genes atp4, ccmB, ccmFc, and mttB might have been positively selected. The results lay the foundation for future studies on the evolution of the Salvia mitogenome and the molecular breeding of S. miltiorrhiza.