Deepred-Mt: Deep representation learning for predicting C-to-U RNA editing in plant mitochondria

Analysis of the complete mitochondrial genome of Panax quinquefolius reveals shifts from cis-splicing to trans-splicing of intron cox2i373

Panax quinquefolius is a perennial plant with medicinal values. In this study, we assembled the complete mitochondrial genome (mitogenome) of P. quinquefolius using PMAT assembler. The total length of P. quinquefolius mitogenome is 573,154 bp. We annotated a total of 34 protein-coding genes (PCGs), 35 tRNA genes, and 6 rRNA genes in this mitogenome. The analysis of repetitive elements shows that there are 153 SSRs, 24 tandem repeats and 242 pairs of dispersed repeats this mitogenome. Also, we found 24 homologous sequences with a total length of 64,070 bp among its mitogenome and plastome, accounting for 41.05 % of the plastome, and 11.18 % of the mitogenome, showing a remarkable frequent sequence dialogue between plastome and mitogenomes. Besides, a total of 583 C to U RNA editing sites on 34 PCGs of high confidence were predicted by using Deepred-mt. We also inferred the phylogenetic relationships of P. quinquefolius and other angiosperms based on mitochondrial PCGs. Finally, we observed a shift from cis- to trans-splicing in P. quinquefolius for two mitochondrial introns, namely cox2i373 and nad1i728, and a pair of 48 bp short repetitive sequences may be associated with the breaking and rearrangement of the cox2i373 intron. The fragmentation of the cox2i373 intron was further confirmed by our PCR amplification experiments. In summary, our report on the P. quinquefolius mitogenome provides a new perspective on the intron evolution of the mitogenome.

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.

Assembly and comparative analysis of the first complete mitochondrial genome of Salix psammophila, a good windbreak and sand fixation shrub

Salix psammophila, commonly known as the sandlive willow, is a vital shrub species within the Salicaceae family, particularly significant for its ecological role in regions susceptible to desertification and sandy soils. In this study, we assembled the complete S. psammophila mitochondrial genome using Pacbio HiFi third-generation sequencing data. The genome was found to be a typical single circular structure, with a total length of 715,555 bp and a GC content of 44.89%. We annotated 33 unique protein-coding genes (PCGs), which included 24 core mitochondrial genes and 9 variable genes, as well as 18 tRNA genes (5 of which were multicopy genes) and 3 rRNA genes. Comparative analysis of the PCGs from the mitochondrial genomes of S. psammophila, Populus deltoides, Populus simonii, Salix wilsonii, and Salix suchowensis revealed that these genes are relatively conserved within the Salicaceae family, with variability primarily occurring in the ribosomal protein genes. The absence of the rps14, which encodes a ribosomal protein, may have played a role in the evolution of stress tolerance in Salicaceae plants. Additionally, we identified 232 SSRs, 19 tandem repeat sequences, and 236 dispersed repeat sequences in the S. psammophila mitochondrial genome, with palindromic and forward repeats being the most abundant. The longest palindromic repeat measured 260 bp, while the longest forward repeat was 86,068 bp. Furthermore, 324 potential RNA editing sites were discovered, all involving C-to-U edits, with the nad4 having the highest number of edits. These findings provide valuable insights into the phylogenetic and genetic research of Salicaceae plants.

Characterization of the complete mitogenome of Tiarella polyphylla, commonly known as Asian foamflower: insights into the multi-chromosomes structure and DNA transfers

BACKGROUND

Tiarella polyphylla D. Don has been traditionally used to cure asthma and skin eruptions. However, the sequence and the structure of the mitogenome of T. polyphylla remained elusive, limiting the genomic and evolution analysis based on the mitogenome.

RESULTS

Using a combination of Illumina and Nanopore sequencing reads, we de novo assembled the complete mitogenome of T. polyphylla. In addition to unveiling the major configuration of the T. polyphylla mitogenome was three circular chromosomes with lengths of 430,435 bp, 126,943 bp, and 55,269 bp, we revealed five (R01-R05) and one (R06) repetitive sequence could mediate the intra- and inter-chromosomal recombination, respectively. Furthermore, we identified 208 short and 25 long tandem segments, seven cp-derived mtDNAs, 106 segments of mtDNAs transferred to the nuclear genome, and 653 predicted RNA editing sites. Based on the sequence of the mitogenomes, we obtained the resolved phylogeny of the seven Saxifragales species.

CONCLUSIONS

These results presented the mitogenome features and expanded its potential applications in phylogenetics, species identification, and cytoplasmic male sterility (CMS) in the future.

The assembly and comparative analysis of the first complete mitogenome of Lindera aggregata

Lindera aggregata, a member belongs to the genus Lindera of Lauraceae family. Its roots and leaves have been used as traditional Chinese medicine or functional food for thousands of years. However, its mitochondrial genome has not been explored. Our aim is to sequence and assemble the mitogenome of L. aggregata to elucidate the genetic mechanism and evolutionary pathway. The results had shown that the mitogenome was extremely complex and had a unique multi-branched conformation with total size of 912,473 bp. Comprehensive analysis of protein coding genes of 7 related species showed that there were 40 common genes in their mitogenome. Interestingly, positive selection had become an important factor in the evolution of ccmB, ccmFC, rps10, rps11 and rps7 genes. Furthermore, our data highlighted the repeated trend of homologous fragment migrations between chloroplast and mitochondrial organelles, and 38 homologous fragments were identified. Phylogenetic analysis identified a tree that was basically consistent with the phylogeny of Laurales species described in the APG IV system. To sum up, this study will be helpful to the study of population genetics and evolution of Lindera species.

Assembly and comparative analysis of the complete mitochondrial genome of Vaccinium carlesii Dunn

Vaccinium L. is an important fruit tree with nutritional, medicinal, and ornamental values. However, the mitochondrial (mt) genome of Vaccinium L. remains largely unexplored. Vaccinium carlesii Dunn is an endemic wild resource in China, which is crucial for blueberry breeding. The V. carlesii mt genomes were sequenced using Illumina and Nanopore, which total length was 636,904 bp with 37 protein coding genes, 20 tRNA genes, and three rRNA genes. We found four pairs of long repeat fragments homologous recombination mediated the generation of substructures in the V. carlesii mt genome. We predicted 383 RNA editing sites, all converting cytosine (C) to uracil (U). According to the phylogenetic analysis, V. carlesii and V. macrocarpon of the Ericaceae exhibited the closest genetic relationship. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of Vaccinium germplasm resources.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

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

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

Complete mitochondrial genome of Melia azedarach L., reveals two conformations generated by the repeat sequence mediated recombination

Melia azedarach is a species of enormous value of pharmaceutical industries. Although the chloroplast genome of M. azedarach has been explored, the information of mitochondrial genome (Mt genome) remains surprisingly limited. In this study, we used a hybrid assembly strategy of BGI short-reads and Nanopore long-reads to assemble the Mt genome of M. azedarach. The Mt genome of M. azedarach is characterized by two circular chromosomes with 350,142 bp and 290,387 bp in length, respectively, which encodes 35 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes. A pair of direct repeats (R1 and R2) were associated with genome recombination, resulting in two conformations based on the Sanger sequencing and Oxford Nanopore sequencing. Comparative analysis identified 19 homologous fragments between Mt and chloroplast genome, with the longest fragment of 12,142 bp. The phylogenetic analysis based on PCGs were consist with the latest classification of the Angiosperm Phylogeny Group. Notably, a total of 356 potential RNA editing sites were predicted based on 35 PCGs, and the editing events lead to the formation of the stop codon in the rps10 gene and the start codons in the nad4L and atp9 genes, which were verified by PCR amplification and Sanger sequencing. Taken together, the exploration of M. azedarach gap-free Mt genome provides a new insight into the evolution research and complex mitogenome architecture.

De Novo Hybrid Assembly Unveils Multi-Chromosomal Mitochondrial Genomes in Ludwigia Species, Highlighting Genomic Recombination, Gene Transfer, and RNA Editing Events

Biological invasions have been identified as the fifth cause of biodiversity loss, and their subsequent dispersal represents a major ecological challenge. The aquatic invasive species Ludwigia grandiflora subsp. hexapetala (Lgh) and Ludwigia peploides subsp. montevidensis (Lpm) are largely distributed in aquatic environments in North America and in Europe. However, they also present worrying terrestrial forms that are able to colonize wet meadows. To comprehend the mechanisms of the terrestrial adaptation of Lgh and Lpm, it is necessary to develop their genomic resources, which are currently poorly documented. We performed de novo assembly of the mitogenomes of Lgh and Lpm through hybrid assemblies, combining short reads (SR) and/or long reads (LR) before annotating both mitogenomes. We successfully assembled the mitogenomes of Lgh and Lpm into two circular molecules each, resulting in a combined total length of 711,578 bp and 722,518 bp, respectively. Notably, both the Lgh and Lpm molecules contained plastome-origin sequences, comprising 7.8% of the mitochondrial genome length. Additionally, we identified recombinations that were mediated by large repeats, suggesting the presence of multiple alternative conformations. In conclusion, our study presents the first high-quality mitogenomes of Lpm and Lgh, which are the only ones in the Myrtales order found as two circular molecules.

Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome

KEY MESSAGE

We reported the mitochondrial genome of Cinnamomum camphora for the first time, revealing frequent rearrangement events in the non-coding regions of Magnoliids mitochondrial genomes. As one of the representative species in the Lauraceae family of Magnoliids, Cinnamomum camphora holds significant economic and ecological value. In this study, the mitochondrial genome (mitogenome) of C. camphora was complete assembled and annotated using PacBio HiFi sequencing. The C. camphora mitogenome is characterized by a branch structure, spans 900,894 bp, and contains 43 protein-coding genes (PCGs), 24 tRNAs, and 3 rRNAs. Most of these PCGs are under purifying selection, with only two (ccmFc and rps7) exhibiting signs of positive selection. The C. camphora mitogenome contains numerous repetitive sequences and intracellular gene transfers, with a total of 36 mitochondrial plastid DNAs, amounting to a combined length of 23,816 bp. Comparative analysis revealed that the non-coding regions of Magnoliids mitogenomes have undergone frequent rearrangements during evolution, but the coding sequences remain highly conserved (more than 98% similarity for protein-coding sequences). Furthermore, a maximum-likelihood phylogenetic tree was reconstructed based on 25 PCGs from 23 plant mitogenomes. The analysis supports the closest relationship between C. camphora and C. chekiangense, consistent with the APG IV classification system. This study elucidates the unique evolutionary features of the C. camphora mitogenome, which will provide valuable insights into the study of genetics and evolution of the family Lauraceae.

Mitochondrial genome complexity in Stemona sessilifolia: nanopore sequencing reveals chloroplast gene transfer and DNA rearrangements

Mitochondria are semi-autonomous organelles in eukaryotic cells with their own genome. Plant mitogenomes differ from animal mitogenomes in size, structure, and repetitive DNA sequences. Despite larger sizes, plant mitogenomes do not have significantly more genes. They exhibit diverse structures due to variations in size, repetitive DNA, recombination frequencies, low gene densities, and reduced nucleotide substitution rates. In this study, we analyzed the mitochondrial genome of Stemona sessilifolia using Nanopore and Illumina sequencing. De-novo assembly and annotation were conducted using Unicycler, Geseq, tRNAscan-SE and BLASTN, followed by codon usage, repeat sequence, RNA-editing, synteny, and phylogenetic analyses. S. sessilifolia's mitogenome consisted of one linear contig and six circular contigs totaling 724,751 bp. It had 39 protein-coding genes, 27 tRNA genes, and 3 rRNA genes. Transfer of chloroplast sequences accounted for 13.14% of the mitogenome. Various analyses provided insights into genetic characteristics, evolutionary dynamics, and phylogenetic placement. Further investigations can explore transferred genes' functions and RNA-editing's role in mitochondrial gene expression in S. sessilifolia.

Analysis of the Rhodomyrtus tomentosa mitochondrial genome: Insights into repeat-mediated recombination and intra-cellular DNA transfer

Plant mitochondrial genomes participate in encoding proteins crucial to the major producers of ATP in the cell and replication and heredity of their own DNA. The sequences and structure of the plant mitochondrial genomes profoundly impact these fundamental processes, and studies of plant mitochondrial genomes are needed. We reported the complete sequences of the Rhodomyrtus tomentosa mitochondrial genome here, totaling 400,482 bp. Nanopore ONT reads and PCR amplification provided evidence for recombination mediated by the eight repeat pairs for the R. tomentosa mitochondrial genome. Thirty-eight genes were identified in the R. tomentosa mitochondrial genome. Comparative analyses of the mitochondrial genome and plastome and PCR amplification suggest that five fragments of mitochondrial plastid DNA were unfunctional sequences resulting from intracellular gene transfer. Phylogenetic analysis based on each and all of the 27 mitochondrial protein-coding genes of nine Myrtales species revealed that R. tomentosa always clustered with other species of Myrtaceae. This study uncovered the enormous complexity of the R. tomentosa mitochondrial genome, the active repeat-mediated recombinations, the presence of mitochondrial plastid DNAs, and the topological incongruence of Myrtales among the single-gene trees.

De Novo Assembly and Comparative Analysis of Mitochondrial Genomes of Two Pueraria montana Varieties

Pueraria montana is a species with important medicinal value and a complex genetic background. In this study, we sequenced and assembled the mitochondrial (mt) genomes of two varieties of P. montana. The mt genome lengths of P. montana var. thomsonii and P. montana var. montana were 457,390 bp and 456,731 bp, respectively. Both P. montana mitogenomes showed a multi-branched structure consisting of two circular molecules, with 56 genes annotated, comprising 33 protein-coding genes, 18 tRNA genes (trnC-GCA and trnM-CAU are multi-copy genes), and 3 rRNA genes. Then, 207 pairs of long repeats and 96 simple sequence repeats (SSRs) were detected in the mt genomes of P. montana, and 484 potential RNA-editing sites were found across the 33 mitochondrial protein-coding genes of each variety. Additionally, a syntenic sequence analysis showed a high collinearity between the two mt genomes. This work is the first to analyze the mt genomes of P. montana. It can provide information that can be used to analyze the structure of mt genomes of higher plants and provide a foundation for future comparative genomic studies and evolutionary biology research in related species.

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

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

The first complete chloroplast genome of Thalictrum fargesii: insights into phylogeny and species identification

Introduction

Thalictrum fargesii is a medicinal plant belonging to the genus Thalictrum of the Ranunculaceae family and has been used in herbal medicine in the Himalayan regions of China and India. This species is taxonomically challenging because of its morphological similarities to other species within the genus. Thus, herbal drugs from this species are frequently adulterated, substituted, or mixed with other species, thereby endangering consumer safety.

Methods

The present study aimed to sequence and assemble the entire chloroplast (cp) genome of T. fargesii using the Illumina HiSeq 2500 platform to better understand the genomic architecture, gene composition, and phylogenetic relationships within the Thalictrum.

Results and discussion

The cp genome was 155,929 bp long and contained large single-copy (85,395 bp) and small single-copy (17,576 bp) regions that were segregated by a pair of inverted repeat regions (26,479 bp) to form a quadripartite structure. The cp genome contains 133 genes, including 88 protein-coding genes (PCGs), 37 tRNA genes, and 8 rRNA genes. Additionally, this genome contains 64 codons that encode 20 amino acids, the most preferred of which are alanine and leucine. We identified 68 SSRs, 27 long repeats, and 242 high-confidence C-to-U RNA-editing sites in the cp genome. Moreover, we discovered seven divergent hotspot regions in the cp genome of T. fargesii, among which ndhD-psaC and rpl16-rps3 may be useful for developing molecular markers for identifying ethnodrug species and their contaminants. A comparative study with eight other species in the genus revealed that pafI and rps19 had highly variable sites in the cp genome of T. fargesii. Additionally, two special features, (i) the shortest length of the ycf1 gene at the IRA-SSC boundary and (ii) the distance between the rps19 fragment and trnH at the IRA-LSC junction, distinguish the cp genome of T. fargesii from those of other species within the genus. Furthermore, phylogenetic analysis revealed that T. fargesii was closely related to T. tenue and T. petaloidium.

Conclusion

Considering all these lines of evidence, our findings offer crucial molecular and evolutionary information that could play a significant role in further species identification, evolution, and phylogenetic studies on T. fargesii.

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

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.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant

BACKGROUND

Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported.

RESULTS

In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted.

CONCLUSIONS

This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

Genomic insights into the clonal reproductive Opuntia cochenillifera: mitochondrial and chloroplast genomes of the cochineal cactus for enhanced understanding of structural dynamics and evolutionary implications

Background

The cochineal cactus (Opuntia cochenillifera), notable for its substantial agricultural and industrial applications, predominantly undergoes clonal reproduction, which presents significant challenges in breeding and germplasm innovation. Recent developments in mitochondrial genome engineering offer promising avenues for introducing heritable mutations, potentially facilitating selective sexual reproduction through the creation of cytoplasmic male sterile genotypes. However, the lack of comprehensive mitochondrial genome information for Opuntia species hinders these efforts. Here, we intended to sequence and characterize its mitochondrial genome to maximize the potential of its genomes for evolutionary studies, molecular breeding, and molecular marker developments.

Results

We sequenced the total DNA of the O. cochenillifera using DNBSEQ and Nanopore platforms. The mitochondrial genome was then assembled using a hybrid assembly strategy using Unicycler software. We found that the mitochondrial genome of O. cochenillifera has a length of 1,156,235 bp, a GC content of 43.06%, and contains 54 unique protein-coding genes and 346 simple repeats. Comparative genomic analysis revealed 48 homologous fragments shared between mitochondrial and chloroplast genomes, with a total length of 47,935 bp. Additionally, the comparison of mitochondrial genomes from four Cactaceae species highlighted their dynamic nature and frequent mitogenomic reorganizations.

Conclusion

Our study provides a new perspective on the evolution of the organelle genome and its potential application in genetic breeding. These findings offer valuable insights into the mitochondrial genetics of Cactaceae, potentially facilitating future research and breeding programs aimed at enhancing the genetic diversity and adaptability of O. cochenillifera by leveraging its unique mitochondrial genome characteristics.

The first mitochondrial genome of Calophyllum soulattri Burm.f

Calophyllum soulattri Burm.f. is traditionally used to treat skin infections and reduce rheumatic pain, yet genetic and genomic studies are still limited. Here, we present the first complete mitochondrial genome of C. soulattri. It is 378,262 bp long with 43.97% GC content, containing 55 genes (30 protein-coding, 5 rRNA, and 20 tRNA). Repeat analysis of the mitochondrial genome revealed 194 SSRs, mostly mononucleotides, and 266 pairs of dispersed repeats ( ≥ 30 bp) that were predominantly palindromic. There were 23 homologous fragments found between the mitochondrial and plastome genomes. We also predicted 345 C-to-U RNA editing sites from 30 protein-coding genes (PCGs) of the C. soulatrii mitochondrial genome. These RNA editing events created the start codon of nad1 and the stop codon of ccmFc. Most PCGs of the C. soulattri mitochondrial genome underwent negative selection, but atp4 and ccmB experienced positive selection. Phylogenetic analyses showed C. soulattri is a sister taxon of Garcinia mangostana. This study has shed light on C. soulattri's evolution and Malpighiales' phylogeny. As the first complete mitochondrial genome in Calophyllaceae, it can be used as a reference genome for other medicinal plant species within the family for future genetic studies.

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.

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

RNA editing events and expression profiles of mitochondrial protein-coding genes in the endemic and endangered medicinal plant, Corydalis saxicola

Corydalis saxicola, an endangered medicinal plant endemic to karst habitats, is widely used in Traditional Chinese Medicine to treat hepatitis, abdominal pain, bleeding hemorrhoids and other conditions. However, to date, the mitochondrial (mt) genome of C. saxicola has not been reported, which limits our understanding of the genetic and biological mechanisms of C. saxicola. Here, the mt genome of C. saxicola was assembled by combining the Nanopore and Illumina reads. The mt genome of C. saxicola is represented by a circular chromosome which is 587,939 bp in length, with an overall GC content of 46.50%. 40 unique protein-coding genes (PCGs), 22 tRNA genes and three rRNA genes were identified. Codon usage of the PCGs was investigated and 167 simple sequence repeats were identified. Twelve homologous fragments were identified between the mt and ct genomes of C. saxicola, accounting for 1.04% of the entire mt genome. Phylogenetic examination of the mt genomes of C. saxicola and 30 other taxa provided an understanding of their evolutionary relationships. We also predicted 779 RNA editing sites in 40 C. saxicola mt PCGs and successfully validated 506 (65%) of these using PCR amplification and Sanger sequencing. In addition, we transcriptionally profiled 24 core mt PCGs in C. saxicola roots treated with different concentrations of CaCl2, as well as in other organs. These investigations will be useful for effective utilization and molecular breeding, and will also provide a reference for further studies of the genus Corydalis.

Decoding the complete organelle genomic architecture of Stewartia gemmata: an early-diverging species in Theaceae

BACKGROUND

Theaceae, comprising 300 + species, holds significance in biodiversity, economics, and culture, notably including the globally consumed tea plant. Stewartia gemmata, a species of the earliest diverging tribe Stewartieae, is critical to offer insights into Theaceae's origin and evolutionary history.

RESULT

We sequenced the complete organelle genomes of Stewartia gemmata using short/long reads sequencing technologies. The chloroplast genome (158,406 bp) exhibited a quadripartite structure including the large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs); 114 genes encoded 80 proteins, 30 tRNAs, and four rRNAs. The mitochondrial genome (681,203 bp) exhibited alternative conformations alongside a monocyclic structure: 61 genes encoding 38 proteins, 20 tRNAs, three rRNAs, and RNA editing-impacting genes, including ATP6, RPL16, COX2, NAD4L, NAD5, NAD7, and RPS1. Comparative analyses revealed frequent recombination events and apparent rRNA gene gains and losses in the mitochondrial genome of Theaceae. In organelle genomes, the protein-coding genes exhibited a strong A/U bias at codon endings; ENC-GC3 analysis implies selection-driven codon bias. Transposable elements might facilitate interorganelle sequence transfer. Phylogenetic analysis confirmed Stewartieae's early divergence within Theaceae, shedding light on organelle genome characteristics and evolution in Theaceae.

CONCLUSIONS

We studied the detailed characterization of organelle genomes, including genome structure, composition, and repeated sequences, along with the identification of lateral gene transfer (LGT) events and complexities. The discovery of a large number of repetitive sequences and simple sequence repeats (SSRs) has led to new insights into molecular phylogenetic markers. Decoding the Stewartia gemmata organellar genome provides valuable genomic resources for further studies in tea plant phylogenomics and evolutionary biology.

Multi-Chromosomal mitochondrial genome of medicinal plant Acorus tatarinowii (Acoraceae): Firstly reported from Acorales Order

Acorus tatarinowii Schott (A. tatarinowii), a well-known traditional Chinese medicinal plant renowned for its high medicinal value, but its mitochondrial genome (mitogenome) is still unexplored. In this study, we meticulously assembled the complete mitochondrial genome of A. tatarinowii using a combination of Illumina short reads and Oxford Nanopore long reads. Our findings revealed that A. tatarinowii possesses a complex chromosomal structural mitogenome, comprising two linear chromosomes and seven circular chromosomes. This mitogenome spans 1.81 Mb in length with a GC content of 38.29 %. Notably, it contained 24 unique mitochondrial core genes, seven unique variable genes, 17 tRNA genes, and three rRNA genes. Analyses of codon usage, most protein-coding genes (PCGs) exhibited a common codon usage preference, with RSCU values greater than 1, and the codon with the highest RSCU value was UAA(End, 1.90). We conducted a thorough analysis of repeat sequences, the distribution of repetitive sequences in nine mitochondrial chromosomes showed distinct patterns. Moreover, we identified 82 and 12 homologous fragments by comparing the sequences of chloroplast and nuclear genomes to the A. tatarinowii mitogenome, respectively. Lastly, We predicted a total of 234 potential RNA editing sites in 28 unique PCGs and discovered that the nad4 gene has been edited the most often, at 26 times. Our results contribute to the enrichment of mitochondrial genome resources for Acoraceae, and the mitogenome also can be used as a reference for other species.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae)

KEY MESSAGE

Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plants feature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604 bp and 46.83%, 288,122 bp and 44.79%, and 287,054 bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editing sites. The Mli mitogenome features four conformations likely mediated by the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

The mitochondrial genomes of Panax notoginseng reveal recombination mediated by repeats associated with DNA replication

Panax notoginseng is one of the most valuable medicinal species. However, its mitochondrial genome has not been reported yet. We aimed to determine the mitogenome sequence of P. notoginseng. We de novo assembled the mitogenome with Illumina short reads and Nanopore long reads. The mitochondrial genome of P. notoginseng has a multipartite structure consisting of interconversion between a "master circle" and numerous "subgenomic circles" through recombinations mediated by 64 pairs of repetitive sequences. Among the multipartite structure, seven subgenomic circles were best supported. Six of the seven subgenomic circles shared an 852 bp conserved fragment. The complete mitogenome of P. notoginseng was 662,479 bp long including 34 mitochondrial protein-coding genes (PCGs), three rRNA, and 19 tRNA genes. We identified 166 microsatellite repeats and 26 long-tandem repeats. Phylogenetic analysis resolved a tree that was mostly congruent with the phylogeny of Apiales species described in the APG IV system and the tree built with the chloroplast genome sequences. A total of 12 mitochondrial plastid DNA fragments were identified. Lastly, we predicted 591C-to-U RNA editing sites in the coding regions of mitochondrial PCGs. The mitochondrial genome will lay the foundation for understanding the evolution of Panax species.

A task-specific encoding algorithm for RNAs and RNA-associated interactions based on convolutional autoencoder

RNAs play essential roles in diverse physiological and pathological processes by interacting with other molecules (RNA/protein/compound), and various computational methods are available for identifying these interactions. However, the encoding features provided by existing methods are limited and the existing tools does not offer an effective way to integrate the interacting partners. In this study, a task-specific encoding algorithm for RNAs and RNA-associated interactions was therefore developed. This new algorithm was unique in (a) realizing comprehensive RNA feature encoding by introducing a great many of novel features and (b) enabling task-specific integration of interacting partners using convolutional autoencoder-directed feature embedding. Compared with existing methods/tools, this novel algorithm demonstrated superior performances in diverse benchmark testing studies. This algorithm together with its source code could be readily accessed by all user at: https://idrblab.org/corain/ and https://github.com/idrblab/corain/.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

Background

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

Results

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

Conclusion

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

The complete mitochondrial genome of Amorphophallus albus and development of molecular markers for five Amorphophallus species based on mitochondrial DNA

Introduction

Amorphophallus albus is an herbaceous, cormous, perennial plant used as a food source and traditional medicine in Asia.

Methods

In this study, we assembled and annotated the complete mitochondrial genome (mitogenome) of A. albus. Then we analyzed the repeated elements and mitochondrial plastid sequences (MTPTs), predicted RNA editing sites in mitochondrial protein-coding genes (PCGs). Lastly, we inferred the phylogenetic relationships of A. albus and other angiosperms based on mitochondrial PCGs, and designed two molecular markers based on mitochondrial DNA.

Results and discussion

The complete mitogenome of A. albus consists of 19 circular chromosomes. And the total length of A. albus mitogenome is 537,044 bp, with the longest chromosome measuring 56,458 bp and the shortest measuring 12,040 bp. We identified and annotated a total of 36 protein-coding genes (PCGs), 21 tRNA genes, and 3 rRNA genes in the mitogenome. Additionally, we analyzed mitochondrial plastid DNAs (MTPTs) and identified 20 MTPTs between the two organelle genomes, with a combined length of 22,421 bp, accounting for 12.76% of the plastome. Besides, we predicted a total of 676 C to U RNA editing sites on 36 protein-coding genes of high confidence using Deepred-mt. Furthermore, extensive genomic rearrangement was observed between A. albus and the related mitogenomes. We conducted phylogenetic analyses based on mitochondrial PCGs to determine the evolutionary relationships between A. albus and other angiosperms. Finally, we developed and validated two molecular markers, Ai156 and Ai976, based on two intron regions (nad2i156 and nad4i976) respectively. The discrimination success rate was 100 % in validation experiments for five widely grown konjac species. Our results reveal the multi-chromosome mitogenome of A. albus, and the developed markers will facilitate molecular identification of this genus.

Assembly of the Complete Mitochondrial Genome of Pereskia aculeata Revealed That Two Pairs of Repetitive Elements Mediated the Recombination of the Genome

Pereskia aculeata is a potential new crop species that has both food and medicinal (antinociceptive activity) properties. However, comprehensive genomic research on P. aculeata is still lacking, particularly concerning its organelle genome. In this study, P. aculeata was studied to sequence the mitochondrial genome (mitogenome) and to ascertain the assembly, informational content, and developmental expression of the mitogenome. The findings revealed that the mitogenome of P. aculeata is circular and measures 515,187 bp in length with a GC content of 44.05%. It contains 52 unique genes, including 33 protein-coding genes, 19 tRNA genes, and three rRNA genes. Additionally, the mitogenome analysis identified 165 SSRs, primarily consisting of tetra-nucleotides, and 421 pairs of dispersed repeats with lengths greater than or equal to 30, which were mainly forward repeats. Based on long reads and PCR experiments, we confirmed that two pairs of long-fragment repetitive elements were highly involved with the mitogenome recombination process. Furthermore, there were 38 homologous fragments detected between the mitogenome and chloroplast genome, and the longest fragment was 3962 bp. This is the first report on the mitogenome in the family Cactaceae. The decoding of the mitogenome of P. aculeata will provide important genetic materials for phylogenetic studies of Cactaceae and promote the utilization of species germplasm resources.

iPReditor-CMG: Improving a predictive RNA editor for crop mitochondrial genomes using genomic sequence features and an optimal support vector machine

In crops, RNA editing is one of the most important post-transcriptional processes in which specific cytidines (C) in virtually all mitochondrial protein-coding genes are converted to uridines (U). Despite extensive recent research in RNA editing, exploring all of the C-to-U editing events efficiently on the genomic scale remains challengeable. Developing accurate prediction methods for the detection of RNA editing sites would dramatically reduce experimental determination. Therefore, we propose a novel method, iPReditor-CMG (improved predictive RNA editor for crop mitochondrial genomes), to predict crop mitochondrial editing sites using genome sequence and an optimised support vector machine (SVM). We first selected three mitochondrial genomes with known RNA editing sites from Arabidopsis thaliana, Brassica napus and Oryza sativa, released by NCBI, as the training and test sets. The genes and their transcripts from self-sequenced tobacco mitochondrial ATPase were selected as the validation set. The iPReditor-CMG first coded the genome sequences as numerical vectors and then performed an efficient feature selection on the high-dimensional feature space, where the SVM was employed in feature selection and following modelling. The average independent prediction accuracy of intraspecific editing sites across three species was 0.85, and up to 0.91 in A. thaliana, which outperformed the reference models. For the interspecific independent prediction, the prediction accuracy between dicotyledons was 0.78 and the accuracy between dicotyledons and monocotyledons was 0.56, which implies that there might be similarity in the C-to-U editing mechanism in close relatives. Finally, the best model was identified with an independent test accuracy of 0.91 and an AUC of 0.88, which suggested that five unreported feature sequences, i.e. TGACA, ACAAC, GTAGA, CCGTT and TAACA, are closely associated with the editing phenomenon. Multiple tests supported that the iPReditor-CMG could be effectively applied to predict editing sites in crop mitochondria, which may further contribute to understanding the mechanisms of site editing and post-transcriptional events in crop mitochondria.

Molecular Diversity and Phylogeny Reconstruction of Genus Colobanthus (Caryophyllaceae) Based on Mitochondrial Gene Sequences

Mitochondrial genomes have become an interesting object of evolutionary and systematic study both for animals and plants, including angiosperms. Although the framework of the angiosperm phylogeny was built on the information derived from chloroplast and nuclear genes, mitochondrial sequences also revealed their usefulness in solving the phylogenetic issues at different levels of plant systematics. Here, we report for the first time the complete sequences of 26 protein-coding genes of eight Colobanthus species (Caryophyllaceae). Of these, 23 of them represented core mitochondrial genes, which are directly associated with the primary function of that organelle, and the remaining three genes represented a facultative set of mitochondrial genes. Comparative analysis of the identified genes revealed a generally high degree of sequence conservation. The Ka/Ks ratio was <1 for most of the genes, which indicated purifying selection. Only for rps12 was Ka/Ks > 1 in all studied species, suggesting positive selection. We identified 146−165 potential RNA editing sites in genes of the studied species, which is lower than in most angiosperms. The reconstructed phylogeny based on mitochondrial genes was consistent with the taxonomic position of the studied species, showing the separate character of the family Caryophyllaceae and close relationships between all studied Colobanthus species, with C. lycopodioides sharing less similarity.

DeepMC-iNABP: Deep learning for multiclass identification and classification of nucleic acid-binding proteins

Nucleic acid-binding proteins (NABPs), including DNA-binding proteins (DBPs) and RNA-binding proteins (RBPs), play vital roles in gene expression. Accurate identification of these proteins is crucial. However, there are two existing challenges: one is the problem of ignoring DNA- and RNA-binding proteins (DRBPs), and the other is a cross-predicting problem referring to DBP predictors predicting DBPs as RBPs, and vice versa. In this study, we proposed a computational predictor, called DeepMC-iNABP, with the goal of solving these difficulties by utilizing a multiclass classification strategy and deep learning approaches. DBPs, RBPs, DRBPs and non-NABPs as separate classes of data were used for training the DeepMC-iNABP model. The results on test data collected in this study and two independent test datasets showed that DeepMC-iNABP has a strong advantage in identifying the DRBPs and has the ability to alleviate the cross-prediction problem to a certain extent. The web-server of DeepMC-iNABP is freely available at http://www.deepmc-inabp.net/. The datasets used in this research can also be downloaded from the website.