The complete nucleotide sequence and multipartite organization of the tobacco mitochondrial genome: comparative analysis of mitochondrial genomes in higher plants

H2O2 negatively regulates aluminum resistance via oxidation and degradation of the transcription factor STOP1

Aluminum (Al) stress triggers the accumulation of hydrogen peroxide (H2O2) in roots. However, whether H2O2 plays a regulatory role in aluminum resistance remains unclear. In this study, we show that H2O2 plays a crucial role in regulation of Al resistance, which is modulated by the mitochondrion-localized pentatricopeptide repeat protein REGULATION OF ALMT1 EXPRESSION 6 (RAE6). Mutation in RAE6 impairs the activity of complex I of the mitochondrial electron transport chain, resulting in the accumulation of H2O2 and increased sensitivity to Al. Our results suggest that higher H2O2 concentrations promote the oxidation of SENSITIVE TO PROTON RHIZOTOXICITY 1 (STOP1), an essential transcription factor that promotes Al resistance, thereby promoting its degradation by enhancing the interaction between STOP1 and the F-box protein RAE1. Conversely, decreasing H2O2 levels or blocking the oxidation of STOP1 leads to greater STOP1 stability and increased Al resistance. Moreover, we show that the thioredoxin TRX1 interacts with STOP1 to catalyze its chemical reduction. Thus, our results highlight the importance of H2O2 in Al resistance and regulation of STOP1 stability in Arabidopsis (Arabidopsis thaliana).

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 comparative analysis of the complete mitochondrial genome of Pinellia ternata

Pinellia ternata is an important natural medicinal herb in China. However, it is susceptible to withering when exposed to high temperatures during growth, which limits its tuber production. Mitochondria usually function in stress response. The P . ternata mitochondrial (mt) genome has yet to be explored. Therefore, we integrated PacBio and Illumina sequencing reads to assemble and annotate the mt genome of P . ternata . The circular mt genome of P . ternata is 876 608bp in length and contains 38 protein-coding genes (PCGs), 20 tRNA genes and three rRNA genes. Codon usage, sequence repeats, RNA editing and gene migration from chloroplast (cp) to mt were also examined. Phylogenetic analysis based on the mt genomes of P . ternata and 36 other taxa revealed the taxonomic and evolutionary status of P . ternata . Furthermore, we investigated the mt genome size and GC content by comparing P . ternata with the other 35 species. An evaluation of non-synonymous substitutions and synonymous substitutions indicated that most PCGs in the mt genome underwent negative selection. Our results provide comprehensive information on the P . ternata mt genome, which may facilitate future research on the high-temperature response of P . ternata and provide new molecular insights on the Araceae family.

Targeted knockout of a conserved plant mitochondrial gene by genome editing

Fusion proteins derived from transcription activator-like effectors (TALEs) have emerged as genome editing tools for mitochondria. TALE nucleases (TALENs) have been applied to delete chimaeric reading frames and duplicated (redundant) genes but produced complex genomic rearrangements due to the absence of non-homologous end-joining. Here we report the targeted deletion of a conserved mitochondrial gene, nad9, encoding a subunit of respiratory complex I. By generating a large number of TALEN-mediated mitochondrial deletion lines, we isolated, in addition to mutants with rearranged genomes, homochondriomic mutants harbouring clean nad9 deletions. Characterization of the knockout plants revealed impaired complex I biogenesis, male sterility and defects in leaf and flower development. We show that these defects can be restored by expressing a functional Nad9 protein from the nuclear genome, thus creating a synthetic cytoplasmic male sterility system. Our data (1) demonstrate the feasibility of using genome editing to study mitochondrial gene functions by reverse genetics, (2) highlight the role of complex I in plant development and (3) provide proof-of-concept for the construction of synthetic cytoplasmic male sterility systems for hybrid breeding by genome editing.

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.

Cytoplasmic male sterility and abortive seed traits generated through mitochondrial genome editing coupled with allotopic expression of atp1 in tobacco

Allotopic expression is the term given for the deliberate relocation of gene function from an organellar genome to the nuclear genome. We hypothesized that the allotopic expression of an essential mitochondrial gene using a promoter that expressed efficiently in all cell types except those responsible for male reproduction would yield a cytoplasmic male sterility (CMS) phenotype once the endogenous mitochondrial gene was inactivated via genome editing. To test this, we repurposed the mitochondrially encoded atp1 gene of tobacco to function in the nucleus under the transcriptional control of a CaMV 35S promoter (construct 35S:nATP1), a promoter that has been shown to be minimally expressed in early stages of anther development. The endogenous atp1 gene was eliminated (Δatp1) from 35S:nATP1 tobacco plants using custom-designed meganucleases directed to the mitochondria. Vegetative growth of most 35S:nATP1/Δatp1 plants appeared normal, but upon flowering produced malformed anthers that failed to shed pollen. When 35S:nATP1/Δatp1 plants were cross-pollinated, ovary/capsule development appeared normal, but the vast majority of the resultant seeds were small, largely hollow and failed to germinate, a phenotype akin to the seedless trait known as stenospermocarpy. Characterization of the mitochondrial genomes from three independent Δatp1 events suggested that spontaneous recombination over regions of microhomology and substoichiometric shifting were the mechanisms responsible for atp1 elimination and genome rearrangement in response to exposure to the atp1-targeting meganucleases. Should the results reported here in tobacco prove to be translatable to other crop species, then multiple applications of allotopic expression of an essential mitochondrial gene followed by its elimination through genome editing can be envisaged. Depending on the promoter(s) used to drive the allotopic gene, this technology may have potential application in the areas of: (1) CMS trait development for use in hybrid seed production; (2) seedless fruit production; and (3) transgene containment.

The complete mitochondrial genome of Allium fistulosum L. (Amaryllidaceae)

Allium fistulosum L. (Linnaeus, Carolus, 1753) is an aromatic vegetable with health benefits and medicinal value. In this study, the complete mitochondrial genome of A. fistulosum was determined. Circular mitochondrial DNA (mtDNA) was 382,053 bp in size, encoded 44 genes, and contained 26 protein-coding genes (PCGs), 14 tRNAs, and four rRNAs. Phylogenetic analysis of amino acid sequences of the 26 PCGs revealed that the closest relationship was between A. fistulosum and Allium cepa. The high-quality mitochondrial genomic sequences of A. fistulosum presented in this study will serve as a useful resource for a range of genetic, functional, evolutionary, and comparative genomic studies on this species of the Amaryllidaceae family.

Development of InDels markers for the identification of cytoplasmic male sterility in Sorghum by complete chloroplast genome sequences analysis

Cytoplasmic male sterility (CMS) is predominantly used for F1 hybrid breeding and seed production in Sorghum. DNA markers to distinguish between normal fertile (CMS-N) and sterile (CMS-S) male cytoplasm can facilitate F1 hybrid cultivar development in Sorghum breeding programs. In this study, the complete chloroplast (cp) genome sequences of CMS-S and Korean Sorghum cultivars were obtained using next-generation sequencing. The de novo assembled genome size of ATx623, the CMS-S line of the chloroplast, was 140,644bp. When compared to the CMS-S and CMS-N cp genomes, 19 single nucleotide polymorphisms (SNPs) and 142 insertions and deletions (InDels) were identified, which can be used for marker development for breeding, population genetics, and evolution studies. Two InDel markers with sizes greater than 20 bp were developed to distinguish cytotypes based on the copy number variation of lengths as 28 and 22 bp tandem repeats, respectively. Using the newly developed InDel markers with five pairs of CMS-S and their near isogenic maintainer line, we were able to easily identify their respective cytotypes. The InDel markers were further examined and applied to 1,104 plants from six Korean Sorghum cultivars to identify variant cytotypes. Additionally, the phylogenetic analysis of seven Sorghum species with complete cp genome sequences, including wild species, indicated that CMS-S and CMS-N contained Milo and Kafir cytotypes that might be hybridized from S. propinquum and S. sudanese, respectively. This study can facilitate F1 hybrid cultivar development by providing breeders with reliable tools for marker-assisted selection to breed desirable Sorghum varieties.

Plant mitochondrial introns as genetic markers - conservation and variation

Plant genomes are comprised of nuclear, plastid and mitochondrial components characterized by different patterns of inheritance and evolution. Genetic markers from the three genomes provide complementary tools for investigations of inheritance, genetic relationships and phenotypic contributions. Plant mitochondrial genomes are challenging for universal marker development because they are highly variable in terms of size, gene order and intergenic sequences and highly conserved with respect to protein-coding sequences. PCR amplification of introns with primers that anneal to conserved, flanking exons is effective for the development of polymorphic nuclear genome markers. The potential for plant mitochondrial intron polymorphisms to distinguish between congeneric species or intraspecific varieties has not been systematically investigated and is possibly constrained by requirements for intron secondary structure and interactions with co-evolved organelle intron splicing factors. To explore the potential for broadly applicable plant mitochondrial intron markers, PCR primer sets based upon conserved sequences flanking 11 introns common to seven angiosperm species were tested across a range of plant orders. PCR-amplified introns were screened for indel polymorphisms among a group of cross-compatible Citrus species and relatives; two Raphanus sativus mitotypes; representatives of the two Phaseolus vulgaris gene pools; and congeneric pairs of Cynodon, Cenchrus, Solanum, and Vaccinium species. All introns were successfully amplified from each plant entry. Length polymorphisms distinguishable by gel electrophoresis were common among genera but infrequent within genera. Sequencing of three introns amplified from 16 entries identified additional short indel polymorphisms and nucleotide substitutions that separated Citrus, Cynodon, Cenchrus and Vaccinium congeners, but failed to distinguish Solanum congeners or representatives of the Phaseolus vulgaris major gene pools. The ability of primer sets to amplify a wider range of plant species' introns and the presence of intron polymorphisms that distinguish congeners was confirmed by in silico analysis. While mitochondrial intron variation is limited in comparison to nuclear introns, these exon-based primer sets provide robust tools for the amplification of mitochondrial introns across a wide range of plant species wherein useful polymorphisms can be identified.

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

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

Polymeric structure of the Cannabis sativa L. mitochondrial genome identified with an assembly graph model

Cannabis sativa L. belongs to the family Cannabaceae in Rosales. It has been widely used as medicines, building materials, and textiles. Elucidating its genome is critical for molecular breeding and synthetic biology study. Many studies have shown that the mitochondrial genomes (mitogenomes) and even chloroplast genomes (plastomes) had complex polymeric structures. Using the Nanopore sequencing platform, we sequenced, assembled, and analyzed its mitogenome and plastome. The resulting unitig graph suggested that the mitogenome had a complex polymeric structure. However, a gap-free, circular sequence was further assembled from the unitig graph. In contrast, a circular sequence representing the plastome was obtained. The mitogenome major conformation was 415,837 bp long, and the plastome was 153,927 bp long. To test if the repeat sequences promote recombination, which corresponds to the branch points in the structure, we tested the sequences around repeats by long-read mapping. Among 208 pairs of predicted repeats, the mapping results supported the presence of cross-over around 25 pairs of repeats. Subsequent PCR amplification confirmed the presence of cross-over around 15 of the 25 repeats. By comparing the mitogenome and plastome sequences, we identified 19 mitochondria plastid DNAs, including seven complete genes (trnW-CCA, trnP-UGG, psbJ, trnN-GUU, trnD-GUC, trnH-GUG, trnM-CAU) and nine gene fragments. Furthermore, the selective pressure analysis results showed that five genes (atp1, ccmB, ccmC, cox1, nad7) had 19 positively selected sites. Lastly, we predicted 28 RNA editing sites. A total of 8 RNA editing sites located in the coding regions were successfully validated by PCR amplification and Sanger sequencing, of which four were synonymous, and four were nonsynonymous. In particular, the RNA editing events appeared to be tissue-specific in C. sativa mitogenome. In summary, we have confirmed the major confirmation of C. sativa mitogenome and characterized its structural features in detail. These results provide critical information for future variety breeding and resource development for C. sativa.

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.

Establishing Physalis as a Solanaceae model system enables genetic reevaluation of the inflated calyx syndrome

The highly diverse Solanaceae family contains several widely studied models and crop species. Fully exploring, appreciating, and exploiting this diversity requires additional model systems. Particularly promising are orphan fruit crops in the genus Physalis, which occupy a key evolutionary position in the Solanaceae and capture understudied variation in traits such as inflorescence complexity, fruit ripening and metabolites, disease and insect resistance, self-compatibility, and most notable, the striking inflated calyx syndrome (ICS), an evolutionary novelty found across angiosperms where sepals grow exceptionally large to encapsulate fruits in a protective husk. We recently developed transformation and genome editing in Physalis grisea (groundcherry). However, to systematically explore and unlock the potential of this and related Physalis as genetic systems, high-quality genome assemblies are needed. Here, we present chromosome-scale references for P. grisea and its close relative Physalis pruinosa and use these resources to study natural and engineered variations in floral traits. We first rapidly identified a natural structural variant in a bHLH gene that causes petal color variation. Further, and against expectations, we found that CRISPR-Cas9-targeted mutagenesis of 11 MADS-box genes, including purported essential regulators of ICS, had no effect on inflation. In a forward genetics screen, we identified huskless, which lacks ICS due to mutation of an AP2-like gene that causes sepals and petals to merge into a single whorl of mixed identity. These resources and findings elevate Physalis to a new Solanaceae model system and establish a paradigm in the search for factors driving ICS.

Comparative analysis of mitochondrial genomes of two alpine medicinal plants of Gentiana (Gentianaceae)

Gentiana crassicaulis and G. straminea are alpine plants of Gentiana with important medicinal value and complex genetic backgrounds. In this study, the mitochondrial genomes (mtDNAs) of these two species were sequenced. The mtDNAs of G. crassicaulis and G. straminea are 368,808 and 410,086 bp long, respectively, 52 and 49 unique genes are annotated in the two species, and the gene arrangement varies widely. Compared to G. crassicaulis, G. straminea loses three effective genes, namely atp6, trnG-GCC and trnV-GAC. As a pseudogene, the atp6 gene of G. straminea is incomplete, which is rare in higher plants. We detected 1696 and 1858 pairs of long repeats and 213 SSRs and 250 SSs in the mtDNAs of G. crassicaulis and G. straminea, respectively. There are 392 SNPs and 18 InDels between the two genomes, and syntenic sequence and structural variation analysis show low collinearity between the two genomes. Chloroplast DNA transferring to mtDNA is observed in both species, and 46,511 and 55,043 bp transferred segments containing three tRNA genes are identified, respectively. Comparative analysis of mtDNAs of G. crassicaulis, G. straminea and four species of Gentianales determined 18 core genes, and there is no specific gene in G. crassicaulis and G. straminea. The phylogenetic tree based on mtDNAs places Gentianaceae in a branch of Gentianales. This study is the first to analyze the mtDNAs of Gentianaceae, which could provide information for analysis of the structure of mtDNAs of higher plants and phylogenetic research of Gentianaceae and Gentianales.

Assembly and comparative analysis of the complete mitochondrial genome of Salix wilsonii using PacBio HiFi sequencing

Salix L. (willows) is one of the most taxonomically complex genera of flowering plants, including shrubs, tall trees, bushes, and prostrate plants. Despite the high species diversity, only five mitochondrial genomes (mitogenomes) have been released in this genus. Salix wilsonii is an important ornamental and economic willow tree in section Wilsonia of the genus Salix. In this study, the S. wilsonii mitogenome was assembled into a typical circular structure with a size of 711,456 bp using PacBio HiFi sequencing. A total of 58 genes were annotated in the S. wilsonii mitogenome, including 33 protein-coding genes (PCGs), 22 tRNAs, and 3 rRNAs. In the S. wilsonii mitogenome, four genes (mttB, nad3, nad4, and sdh4) were found to play important roles in its evolution through selection pressure analysis. Collinearity analysis of six Salix mitogenomes revealed high structural variability. To determine the evolutionary position of S. wilsonii, we conducted a phylogenetic analysis of the mitogenomes of S. wilsonii and 12 other species in the order Malpighiales. Results strongly supported the segregation of S. wilsonii and other five Salix species with 100% bootstrap support. The comparative analysis of the S. wilsonii mitogenome not only sheds light on the functional and structural features of S. wilsonii but also provides essential information for genetic studies of the genus Salix.

A High-Quality Genome Assembly of the Mitochondrial Genome of the Oil-Tea Tree Camellia gigantocarpa (Theaceae)

Camellia gigantocarpa is one of the oil-tea trees whose seeds can be used to extract high-quality vegetable oil. To date, there are no data on the mitochondrial genome of the oil-tea tree, in contrast to the tea-tree C. sinensis, which belongs to the same genus. In this paper, we present the first complete mitochondrial genomes of C. gigantocarpa obtained using PacBio Hi-Fi (high-fidelity) and Hi-C sequencing technologies to anchor the 970,410 bp genome assembly into a single sequence. A set of 44 protein-coding genes, 22 non-coding genes, 746 simple sequence repeats (SSRs), and more than 201 kb of repetitive sequences were annotated in the genome assembly. The high percentage of repetitive sequences in the mitochondrial genome of C. gigantocarpa (20.81%) and C.sinensis (22.15%, tea tree) compared to Arabidopsis thaliana (4.96%) significantly increased the mitogenome size in the genus Camellia. The comparison of the mitochondrial genomes between C. gigantocarpa and C. sinensis revealed genes exhibit high variance in gene order and low substitution rate within the genus Camellia. Information on the mitochondrial genome provides a better understanding of the structure and evolution of the genome in Camellia and may contribute to further study of the after-ripening process of oil-tea trees.

Evaluation of Intracellular Gene Transfers from Plastome to Nuclear Genome across Progressively Improved Assemblies for Arabidopsis thaliana and Oryza sativa

DNA originating from organellar genomes are regularly discovered in nuclear sequences during genome assembly. Nevertheless, such insertions are sometimes omitted during the process of nuclear genome assembly because the inserted DNA is assigned to organellar genomes, leading to a systematic underestimation of their frequency. With the rapid development of high-throughput sequencing technology, more inserted fragments from organelle genomes can now be detected. Therefore, it is necessary to be aware of the insertion events from organellar genomes during nuclear genome assembly to properly attribute the impact and rate of such insertions in the evolution of nuclear genomes. Here, we investigated the impact of intracellular gene transfer (IGT) from the plastome to the nuclear genome using genome assemblies that were refined through time with technological improvements from two model species, Arabidopsis thaliana and Oryza sativa. We found that IGT from the plastome to the nuclear genome is a dynamic and ongoing process in both A. thaliana and O. sativa, and mostly occurred recently, as the majority of transferred sequences showed over 95% sequence similarity with plastome sequences of origin. Differences in the plastome-to-nuclear genome IGT between A. thaliana and O. sativa varied among the different assembly versions and were associated with the quality of the nuclear genome assembly. IGTs from the plastome to nuclear genome occurred more frequently in intergenic regions, which were often associated with transposable elements (TEs). This study provides new insights into intracellular genome evolution and nuclear genome assembly by characterizing and comparing IGT from the plastome into the nuclear genome for two model plant species.

Revealing the Complete Chloroplast Genome of an Andean Horticultural Crop, Sweet Cucumber (Solanum muricatum), and Its Comparison with Other Solanaceae Species

Sweet cucumber (Solanum muricatum) sect. Basarthrum is a neglected horticultural crop native to the Andean region. It is naturally distributed very close to other two Solanum crops of high importance, potatoes, and tomatoes. To date, molecular tools for this crop remain undetermined. In this study, the complete sweet cucumber chloroplast (cp) genome was obtained and compared with seven Solanaceae species. The cp genome of S. muricatum was 155,681 bp in length and included a large single copy (LSC) region of 86,182 bp and a small single-copy (SSC) region of 18,360 bp, separated by a pair of inverted repeats (IR) regions of 25,568 bp. The cp genome possessed 87 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and one pseudogene. Furthermore, 48 perfect microsatellites were identified. These repeats were mainly located in the noncoding regions. Whole cp genome comparative analysis revealed that the SSC and LSC regions showed more divergence than IR regions. Similar to previous studies, our phylogenetic analysis showed that S. muricatum is a sister species to members of sections Petota + Lycopersicum + Etuberosum. We expect that this first sweet cucumber chloroplast genome will provide potential molecular markers and genomic resources to shed light on the genetic diversity and population studies of S. muricatum, which will allow us to identify varieties and ecotypes. Finally, the features and the structural differentiation will provide us with information about the genes of interest, generating tools for the most precise selection of the best individuals of sweet cucumber, in less time and with fewer resources.

Rearrangement and domestication as drivers of Rosaceae mitogenome plasticity

Background

The mitochondrion is an important cellular component in plants and that functions in producing vital energy for the cell. However, the evolution and structure of mitochondrial genomes (mitogenomes) remain unclear in the Rosaceae family. In this study, we assembled 34 Rosaceae mitogenomes and characterized genome variation, rearrangement rate, and selection signal variation within these mitogenomes.

Results

Comparative analysis of six genera from the Amygdaloideae and five from the Rosoideae subfamilies of Rosaceae revealed that three protein-coding genes were absent from the mitogenomes of five Rosoideae genera. Positive correlations between genome size and repeat content were identified in 38 Rosaceae mitogenomes. Twenty repeats with high recombination frequency (> 50%) provided evidence for predominant substoichiometric conformation of the mitogenomes. Variations in rearrangement rates were identified between eleven genera, and within the Pyrus, Malus, Prunus, and Fragaria genera. Based on population data, phylogenetic inferences from Pyrus mitogenomes supported two distinct maternal lineages of Asian cultivated pears. A Pyrus-specific deletion (DEL-D) in selective sweeps was identified based on the assembled genomes and population data. After the DEL-D sequence fragments originally arose, they may have experienced a subsequent doubling event via homologous recombination and sequence transfer in the Amygdaloideae; afterwards, this variant sequence may have significantly expanded to cultivated groups, thereby improving adaptation during the domestication process.

Conclusions

This study characterizes the variations in gene content, genome size, rearrangement rate, and the impact of domestication in Rosaceae mitogenomes and provides insights into their structural variation patterns and phylogenetic relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01383-3.

Comparative analysis of the plastid and mitochondrial genomes of Artemisia giraldii Pamp

Artemisia giraldii Pamp. is an herbaceous plant distributed only in some areas in China. To understand the evolutionary relationship between plastid and mitochondria in A. giraldii, we sequenced and analysed the plastome and mitogenome of A. giraldii on the basis of Illumina and Nanopore DNA sequencing data. The mitogenome was 194,298 bp long, and the plastome was 151,072 bp long. The mitogenome encoded 56 genes, and the overall GC content was 45.66%. Phylogenetic analysis of the two organelle genomes revealed that A. giraldii is located in the same branching position. We found 13 pairs of homologous sequences between the plastome and mitogenome, and only one of them might have transferred from the plastid to the mitochondria. Gene selection pressure analysis in the mitogenome showed that ccmFc, nad1, nad6, atp9, atp1 and rps12 may undergo positive selection. According to the 18 available plastome sequences, we found 17 variant sites in two hypervariable regions that can be used in completely distinguishing 18 Artemisia species. The most interesting discovery was that the mitogenome of A. giraldii was only 43,226 bp larger than the plastome. To the best of our knowledge, this study represented one of the smallest differences between all sequenced mitogenomes and plastomes from vascular plants. The above results can provide a reference for future taxonomic and molecular evolution studies of Asteraceae species.

Characterization and phylogenetic analysis of the complete mitochondrial genome sequence of Diospyros oleifera, the first representative from the family Ebenaceae

Plant mitochondrial genomes are a valuable source of genetic information for a better understanding of phylogenetic relationships. However, no mitochondrial genome of any species in Ebenaceae has been reported. In this study, we reported the first mitochondrial genome of an Ebenaceae model plant Diospyros oleifera. The mitogenome was 493,958 bp in length, contained 39 protein-coding genes, 27 transfer RNA genes, and 3 ribosomal RNA genes. The rps2 and rps11 genes were missing in the D. oleifera mt genome, while the rps10 gene was identified. The length of the repetitive sequence in the D. oleifera mt genome was 31 kb, accounting for 6.33%. A clear bias in RNA-editing sites were found in the D. oleifera mt genome. We also detected 28 chloroplast-derived fragments significantly associated with D. oleifera mt genes, indicating intracellular tRNA genes transferred frequently from chloroplasts to mitochondria in D. oleifera. Phylogenetic analysis based on the mt genomes of D. oleifera and 27 other taxa reflected the exact evolutionary and taxonomic status of D. oleifera. Ka/Ks analysis revealed that 95.16% of the protein-coding genes in the D. oleifera mt genome had undergone negative selections. But, the rearrangement of mitochondrial genes has been widely occur among D. oleifera and these observed species. These results will lay the foundation for identifying further evolutionary relationships within Ebenaceae.

Graph-based models of the Oenothera mitochondrial genome capture the enormous complexity of higher plant mitochondrial DNA organization

Plant mitochondrial genomes display an enormous structural complexity, as recombining repeat-pairs lead to the generation of various sub-genomic molecules, rendering these genomes extremely challenging to assemble. We present a novel bioinformatic data-processing pipeline called SAGBAC (Semi-Automated Graph-Based Assembly Curator) that identifies recombinogenic repeat-pairs and reconstructs plant mitochondrial genomes. SAGBAC processes assembly outputs and applies our novel ISEIS (Iterative Sequence Ends Identity Search) algorithm to obtain a graph-based visualization. We applied this approach to three mitochondrial genomes of evening primrose (Oenothera), a plant genus used for cytoplasmic genetics studies. All identified repeat pairs were found to be flanked by two alternative and unique sequence-contigs defining so-called 'double forks', resulting in four possible contig-repeat-contig combinations for each repeat pair. Based on the inferred structural models, the stoichiometry of the different contig-repeat-contig combinations was analyzed using Illumina mate-pair and PacBio RSII data. This uncovered a remarkable structural diversity of the three closely related mitochondrial genomes, as well as substantial phylogenetic variation of the underlying repeats. Our model allows predicting all recombination events and, thus, all possible sub-genomes. In future work, the proposed methodology may prove useful for the investigation of the sub-genome organization and dynamics in different tissues and at various developmental stages.

CFM6 is an Essential CRM Protein Required for the Splicing of nad5 Transcript in Arabidopsis Mitochondria

Plant chloroplast RNA splicing and ribosome maturation (CRM)-domain-containing proteins are capable of binding RNA to facilitate the splicing of group I or II introns in chloroplasts, but their functions in mitochondria are less clear. In the present study, Arabidopsis thaliana CFM6, a protein with a single CRM domain, was expressed in most plant tissues, particularly in flower tissues, and restricted to mitochondria. Mutation of CFM6 causes severe growth defects, including stunted growth, curled leaves, delayed embryogenesis and pollen development. CFM6 functions specifically in the splicing of group II intron 4 of nad5, which encodes a subunit of mitochondrial complex I, as evidenced by the loss of nad5 intron 4 splicing and high accumulation of its pretranscripts in cfm6 mutants. The phenotypic and splicing defects of cfm6 were rescued in transgenic plants overexpressing 35S::CFM6-YFP. Splicing failure in cfm6 also led to the loss of complex I activity and to its improper assembly. Moreover, dysfunction of complex I induced the expression of proteins or genes involved in alternative respiratory pathways in cfm6. Collectively, CFM6, a previously uncharacterized CRM domain-containing protein, is specifically involved in the cis-splicing of nad5 intron 4 and plays a pivotal role in mitochondrial complex I biogenesis and normal plant growth.

A complete sequence of mitochondrial genome of Neolamarckia cadamba and its use for systematic analysis

Neolamarckia cadamba is an important tropical and subtropical tree for timber industry in southern China and is also a medicinal plant because of the secondary product cadambine. N. cadamba belongs to Rubiaceae family and its taxonomic relationships with other species are not fully evaluated based on genome sequences. Here, we report the complete sequences of mitochondrial genome of N. cadamba, which is 414,980 bp in length and successfully assembled in two genome circles (109,836 bp and 305,144 bp). The mtDNA harbors 83 genes in total, including 40 protein-coding genes (PCGs), 31 transfer RNA genes, 6 ribosomal RNA genes, and 6 other genes. The base composition of the whole genome is estimated as 27.26% for base A, 22.63% for C, 22.53% for G, and 27.56% for T, with the A + T content of 54.82% (54.45% in the small circle and 54.79% in the large circle). Repetitive sequences account for ~ 0.14% of the whole genome. A maximum likelihood (ML) tree based on DNA sequences of 24 PCGs supports that N. cadamba belongs to order Gentianales. A ML tree based on rps3 gene of 60 species in family Rubiaceae shows that N. cadamba is more related to Cephalanthus accidentalis and Hymenodictyon parvifolium and belongs to the Cinchonoideae subfamily. The result indicates that N. cadamba is genetically distant from the species and genera of Rubiaceae in systematic position. As the first sequence of mitochondrial genome of N. cadamba, it will provide a useful resource to investigate genetic variation and develop molecular markers for genetic breeding in the future.

Shifts from cis-to trans-splicing of five mitochondrial introns in Tolypanthus maclurei

Shifts from cis-to trans-splicing of mitochondrial introns tend to correlate with relative genome rearrangement rates during vascular plant evolution, as is particularly apparent in some lineages of gymnosperms. However, although many angiosperms have also relatively high mitogenomic rearrangement rates, very few cis-to trans-splicing shifts except for five trans-spliced introns shared in seed plants have been reported. In this study, we sequenced and characterized the mitogenome of Tolypanthus maclurei, a hemiparasitic plant from the family Loranthaceae (Santalales). The mitogenome was assembled into a circular chromosome of 256,961 bp long, relatively small compared with its relatives from Santalales. It possessed a gene content of typical angiosperm mitogenomes, including 33 protein-coding genes, three rRNA genes and ten tRNA genes. Plastid-derived DNA fragments took up 9.1% of the mitogenome. The mitogenome contained one group I intron (cox1i729) and 23 group II introns. We found shifts from cis-to trans-splicing of five additional introns in its mitogenome, of which two are specific in T. maclurei. Moreover, atp1 is a chimeric gene and phylogenetic analysis indicated that a 356 bp region near the 3′ end of atp1 of T. maclurei was acquired from Lamiales via horizontal gene transfer. Our results suggest that shifts to trans-splicing of mitochondrial introns may not be uncommon among angiosperms.

Mitochondrial genes are involved in the fertility transformation of the thermosensitive male-sterile line YS3038 in wheat

Heterosis can improve the stress resistance, quality, and yield of crops, and the male sterility of wheat can be utilized to accelerate the breeding process of hybrid. To determine whether mitochondrial genes are involved in the fertility of K-type cytoplasmic male-sterile (CMS) line and the YS-type thermosensitive male-sterile (TMS) line in wheat, we sequenced and assembled the mitochondrial genomes of K519A, 519B, and YS3038 by next-generation sequencing (NGS). The non-synonymous mutations were analyzed, and the first-generation sequencing was conducted to verify the non-synonymous mutation sites. Furthermore, the expression patterns of genes with non-synonymous mutations were analyzed. Finally, the candidate genes were silenced by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) to test the functions of the candidate genes. The results revealed that the mitochondrial genomes of K519A, 519B, and YS3038 were 420,543, 433,560, and 452,567 bp in length, respectively. Besides, 33, 31, and 37 protein-coding genes were identified in K519A, 519B, and YS3038, respectively. There were 14 protein-coding genes and 83 open reading frame (ORF) sequences that differed between K519A and 519B and 10 protein-coding genes and 122 ORF sequences that differed between K519A and YS3038. At the binucleate stage, seven genes (nad6, ORF256, ORF216, ORF138, atp6, nad3, and cox1) were downregulated in K519A compared with 519B, and 10 genes (nad6, atp6, cox3, atp8, nad3, cox1, rps3, ORF216, ORF138, and ORF224) were downregulated in YS3038 compared with K519A. Besides, six genes (nad6, ORF138, cox3, cox1, rps3, and ORF224) were downregulated under fertile conditions relative to sterile conditions in YS3038. Gene silencing analysis showed that the silencing of cox1 significantly reduced the seed setting rate of YS3038, indicating that the cox1 gene may be involved in the fertility transformation of YS3038.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01252-x.

Cotton Fiber Development Requires the Pentatricopeptide Repeat Protein GhIm for Splicing of Mitochondrial nad7 mRNA

Pentatricopeptide repeat (PPR) proteins encoded by nuclear genomes can bind to organellar RNA and are involved in the regulation of RNA metabolism. However, the functions of many PPR proteins remain unknown in plants, especially in polyploidy crops. Here, through a map-based cloning strategy and Clustered regularly interspaced short palindromic repeats/cas9 (CRISPR/cas9) gene editing technology, we cloned and verified an allotetraploid cotton immature fiber (im) mutant gene (GhImA) encoding a PPR protein in chromosome A03, that is associated with the non-fluffy fiber phenotype. GhImA protein targeted mitochondrion and could bind to mitochondrial nad7 mRNA, which encodes the NAD7 subunit of Complex I. GhImA and its homolog GhImD had the same function and were dosage-dependent. GhImA in the im mutant was a null allele with a 22 bp deletion in the coding region. Null GhImA resulted in the insufficient GhIm dosage, affected mitochondrial nad7 pre-mRNA splicing, produced less mature nad7 transcripts, and eventually reduced Complex I activities, up-regulated alternative oxidase metabolism, caused reactive oxygen species (ROS) burst and activation of stress or hormone response processes. This study indicates that the GhIm protein participates in mitochondrial nad7 splicing, affects respiratory metabolism, and further regulates cotton fiber development via ATP supply and ROS balance.

Transcriptional Landscape and Splicing Efficiency in Arabidopsis Mitochondria

Plant mitochondrial transcription is initiated from multiple promoters without an apparent motif, which precludes their identification in other species based on sequence comparisons. Even though coding regions take up only a small fraction of plant mitochondrial genomes, deep RNAseq studies uncovered that these genomes are fully or nearly fully transcribed with significantly different RNA read depth across the genome. Transcriptomic analysis can be a powerful tool to understand the transcription process in diverse angiosperms, including the identification of potential promoters and co-transcribed genes or to study the efficiency of intron splicing. In this work, we analyzed the transcriptional landscape of the Arabidopsis mitochondrial genome (mtDNA) based on large-scale RNA sequencing data to evaluate the use of RNAseq to study those aspects of the transcription process. We found that about 98% of the Arabidopsis mtDNA is transcribed with highly different RNA read depth, which was elevated in known genes. The location of a sharp increase in RNA read depth upstream of genes matched the experimentally identified promoters. The continuously high RNA read depth across two adjacent genes agreed with the known co-transcribed units in Arabidopsis mitochondria. Most intron-containing genes showed a high splicing efficiency with no differences between cis and trans-spliced introns or between genes with distinct splicing mechanisms. Deep RNAseq analyses of diverse plant species will be valuable to recognize general and lineage-specific characteristics related to the mitochondrial transcription process.

Characterization of the Mitochondrial Genome of a Wheat AL-Type Male Sterility Line and the Candidate CMS Gene

Heterosis utilization is very important in hybrid seed production. An AL-type cytoplasmic male sterile (CMS) line has been used in wheat-hybrid seed production, but its sterility mechanism has not been explored. In the present study, we sequenced and verified the candidate CMS gene in the AL-type sterile line (AL18A) and its maintainer line (AL18B). In the late uni-nucleate stage, the tapetum cells of AL18A showed delayed programmed cell death (PCD) and termination of microspore at the bi-nucleate stage. As compared to AL18B, the AL18A line produced 100% aborted pollens. The mitochondrial genomes of AL18A and AL18B were sequenced using the next generation sequencing such as Hiseq and PacBio. It was found that the mitochondrial genome of AL18A had 99% similarity with that of Triticum timopheevii, AL18B was identical to that of Triticum aestivum cv. Chinese Yumai. Based on transmembrane structure prediction, 12 orfs were selected as candidate CMS genes, including a previously suggested orf256. Only the lines harboring orf279 showed sterility in the transgenic Arabidopsis system, indicating that orf279 is the CMS gene in the AL-type wheat CMS lines. These results provide a theoretical basis and data support to further analyze the mechanism of AL-type cytoplasmic male sterility in wheat.

Computational Detection of Plant RNA Editing Events

Computers are able to systematically exploit RNA-seq data allowing us to efficiently detect RNA editing sites in a genome-wide scale. This chapter introduces a very flexible computational framework for detecting RNA editing sites in plant organelles. This framework comprises three major steps: RNA-seq data processing, RNA read alignment, and RNA editing site detection. Each step is discussed in sufficient detail to be implemented by the reader. As a study case, the framework will be used with publicly available sequencing data to detect C-to-U RNA editing sites in the coding sequences of the mitochondrial genome of Nicotiana tabacum.

Plastid genomics of Nicotiana (Solanaceae): insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica)

Species of the genus Nicotiana (Solanaceae), commonly referred to as tobacco plants, are often cultivated as non-food crops and garden ornamentals. In addition to the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species: N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes had the typical quadripartite structure, consisting of a pair of inverted repeat (IR) regions (25,323-25,369 bp each) separated by a large single-copy (LSC) region (86,510-86,716 bp) and a small single-copy (SSC) region (18,441-18,555 bp). Comparative analyses of Nicotiana plastid genomes with currently available Solanaceae genome sequences showed similar GC and gene content, codon usage, simple sequence and oligonucleotide repeats, RNA editing sites, and substitutions. We identified 20 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which could be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. Relaxed molecular clock analyses estimated the speciation event between N. rustica and N. knightiana appeared 0.56 Ma (HPD 0.65-0.46). Biogeographical analysis supported a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicate that four genes involved in different plastid functions, including DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations in these genes might have contributed to better survival and superior adaptations during the evolutionary history of tobacco species.

Characterization and Analysis of the Mitochondrial Genome of Common Bean (Phaseolus vulgaris) by Comparative Genomic Approaches

The common bean (Phaseolus vulgaris) is a major source of protein and essential nutrients for humans. To explore the genetic diversity and phylogenetic relationships of P. vulgaris, its complete mitochondrial genome (mitogenome) was sequenced and assembled. The mitogenome is 395,516 bp in length, including 31 unique protein-coding genes (PCGs), 15 transfer RNA (tRNA) genes, and 3 ribosomal RNA (rRNA) genes. Among the 31 PCGs, four genes (mttB, nad1, nad4L, and rps10) use ACG as initiation codons, which are altered to standard initiation codons by RNA editing. In addition, the termination codon CGA in the ccmFC gene is converted to UGA. Selective pressure analysis indicates that the ccmB, ccmFC, rps1, rps10, and rps14 genes were under evolutionary positive selection. The proportions of five amino acids (Phe, Leu, Pro, Arg, and Ser) in the whole amino acid profile of the proteins in each mitogenome can be used to distinguish angiosperms from gymnosperms. Phylogenetic analyses show that P. vulgaris is evolutionarily closer to the Glycininae than other leguminous plants. The results of the present study not only provide an important opportunity to conduct further genomic breeding studies in the common bean, they also provide valuable information for future evolutionary and molecular studies of leguminous plants.

Comparative Transcriptome Analysis Reveals the Potential Mechanism of Abortion in Tobacco sua-Cytoplasmic Male Sterility

sua-CMS (cytoplasmic male sterility) is the only male sterile system in tobacco breeding, but the mechanism of abortion is unclear. Cytological characteristics show that abortion in the sua-CMS line msZY occurs before the differentiation of sporogenous cells. In this study, a comparative transcriptomic analysis was conducted on flower buds at the abortion stage of msZY and its male fertile control ZY. A total of 462 differentially expressed genes were identified in msZY and ZY, which were enriched via protein processing in the endoplasmic reticulum (ER), oxidative phosphorylation, photosynthesis, and circadian rhythm-plant by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. Most genes were downregulated in the ER stress pathway, heat-shock protein family, F1F0-ATPase encoding by the mitochondrial genome, and differentiation of stamens. Genes in the programmed cell death (PCD) pathway were upregulated in msZY. The transcriptome results were consistent with those of qRT-PCR. Ultrastructural and physiological analyses indicted active vacuole PCD and low ATP content in msZY young flower buds. We speculated that PCD and a deficiency in ATP synthesis are essential for the abortion of sua-CMS. This study reveals the potential mechanism of abortion of tobacco sua-CMS.

Comparative Analysis of the Mitochondrial Genomes of Nicotiana tabacum: Hints Toward the Key Factors Closely Related to the Cytoplasmic Male Sterility Mechanism

Background

Cytoplasmic male sterility (CMS) is a complex phenomenon of plant sterility that can produce non-functional pollen. It is caused by mutation, rearrangement or recombination in the mitochondrial genome. So far, the systematic structural characteristics of the changes in the mitochondrial genome from the maintainer lines to the CMS lines have not been reported in tobacco.

Results

The mitochondrial genomes of the flower buds from both CMS lines and maintainer lines of two Nicotiana tabacum cultivars (YY85, sYY85, ZY90, and sZY90) were sequenced using the PacBio and Illumina Hiseq technology, and several findings were produced by comparative analysis based on the de novo sequencing. (1) The genomes of the CMS lines were larger, and the different areas were mostly non-coding regions. (2) A large number of rearrangement regions were detected in the CMS lines, with many translocation regions. (3) Thirteen gene clusters were shared by the four mitochondrial genomes, among which two of the gene clusters, nad2-sdh3 and nad6-rps4, were far from each other in the CMS lines. (4) Thirty-three protein-coding genes were conserved in four mitochondrial genomes. However, nad3 was detected one additional copy in the maintainer lines, and sequence differences were revealed in the four candidate genes (atp6, cox2, nad2, and sdh3). Importantly, the evolutionary tree based on the four genes could be used to distinguish the CMS lines and the maintainer lines well for the sequenced mitochondrial genomes of the tobacco. (5) Sixteen CMS-specific open reading frames (ORFs) were found, three of which (orf91, orf115b, and orf100) were previously reported. (6) The differences in intensity of the protein–protein (PPI) interaction in ATP6 were further verified using the yeast two-hybrid analysis.

Conclusion

Although the majority of the sequences, genes and gene clusters were shared by the mitochondrial genomes of the maintainer and the CMS lines in tobacco, extensive structural variations identified with comprehensive analysis based on the mitochondrial genomes, including rearrangement, gene order, the mitochondrial genome expansion and shrinkage events, might be related to CMS. Additionally, the candidate protein-coding genes and CMS-specific ORFs were closely associated with the CMS mechanism. Verification experiments of one of the candidate genes were performed, and the validity of our research results was supported.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BACKGROUND

Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.

RESULTS

Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.

CONCLUSIONS

We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Whole Plastome Sequencing Within Silene Section Psammophilae Reveals Mainland Hybridization and Divergence With the Balearic Island Populations

Reconstructing the phylogenetic relationships within Caryophyllaceae tribe Sileneae has been obscured by hybridization and incomplete lineage sorting. Silene is the largest genus in the Caryophyllaceae, and unraveling its evolutionary history has been particularly challenging. In order to infer the phylogenetic relationships among the five species in Silene section Psammophilae, we have performed a genome skimming approach to acquire the complete plastid genome (cpDNA), nuclear ribosomal cistron (nrDNA), and partial mitochondrial genome (mtDNA). We have included 26 populations, representing the range of each species' distribution. This section includes five morphologically similar species endemic to the Iberian Peninsula and Balearic Islands (Ibiza and Formentera), yet some of them occupy distinct edaphic habitats (e.g. maritime sands, calcareous sandstones). In addition to phylogeographic analyses, genetic structuring using the chloroplast data set was inferred with Discriminant Analysis of Principal Components (DAPC), analyses of molecular variance (AMOVA), and a partial Mantel test. Reference-guided assembly of 50 bp single-end and 250 bp paired-end Illumina reads produced the nearly complete cpDNA genome (154 kbp), partial mtDNA genome (from 81 to 114 kbp), and the nrDNA cistron (6.4 kbp). Selected variable regions of the cpDNA and mtDNA assemblies were confirmed by Sanger sequencing. Phylogenetic analyses of the mainland populations reveal incongruence among the three genomes. None of the three data sets produced relationships consistent with taxonomy or geography. In contrast, Silene cambessedesii, present in the Balearic Islands, is the only species that forms a strongly supported monophyletic clade in the cpDNA genome and is strongly differentiated with respect to the remaining taxa of the Iberian Peninsula. These results contrast with those obtained for mainland populations. Across the entire analysis, only one well-supported mainland clade of Silene littorea and Silene stockenii emerges from the southern region of the Iberian Peninsula. DAPC and AMOVA results suggest the absence of genetic structure among mainland populations of Silene section Psammophilae, whereas partial Mantel test discarded spatial correlation of genetic differentiation. The widespread incongruence between morphology-based taxonomic boundaries and phylogeography suggests a history of interspecific hybridization, in which only a substantial geographic barrier, like isolation by the Mediterranean Sea, was sufficient to create and maintain species boundaries in Silene section Psammophilae.

Intergenomic gene transfer in diploid and allopolyploid Gossypium

BACKGROUND

Intergenomic gene transfer (IGT) between nuclear and organellar genomes is a common phenomenon during plant evolution. Gossypium is a useful model to evaluate the genomic consequences of IGT for both diploid and polyploid species. Here, we explore IGT among nuclear, mitochondrial, and plastid genomes of four cotton species, including two allopolyploids and their model diploid progenitors (genome donors, G. arboreum: A₂ and G. raimondii: D₅).

RESULTS

Extensive IGT events exist for both diploid and allotetraploid cotton (Gossypium) species, with the nuclear genome being the predominant recipient of transferred DNA followed by the mitochondrial genome. The nuclear genome has integrated 100 times more foreign sequences than the mitochondrial genome has in total length. In the nucleus, the integrated length of chloroplast DNA (cpDNA) was between 1.87 times (in diploids) to nearly four times (in allopolyploids) greater than that of mitochondrial DNA (mtDNA). In the mitochondrion, the length of nuclear DNA (nuDNA) was typically three times than that of cpDNA. Gossypium mitochondrial genomes integrated three nuclear retrotransposons and eight chloroplast tRNA genes, and incorporated chloroplast DNA prior to divergence between the diploids and allopolyploid formation. For mitochondrial chloroplast-tRNA genes, there were 2-6 bp conserved microhomologies flanking their insertion sites across distantly related genera, which increased to 10 bp microhomologies for the four cotton species studied. For organellar DNA sequences, there are source hotspots, e.g., the atp6-trnW intergenic region in the mitochondrion and the inverted repeat region in the chloroplast. Organellar DNAs in the nucleus were rarely expressed, and at low levels. Surprisingly, there was asymmetry in the survivorship of ancestral insertions following allopolyploidy, with most numts (nuclear mitochondrial insertions) decaying or being lost whereas most nupts (nuclear plastidial insertions) were retained.

CONCLUSIONS

This study characterized and compared intracellular transfer among nuclear and organellar genomes within two cultivated allopolyploids and their ancestral diploid cotton species. A striking asymmetry in the fate of IGTs in allopolyploid cotton was discovered, with numts being preferentially lost relative to nupts. Our results connect intergenomic gene transfer with allotetraploidy and provide new insight into intracellular genome evolution.

CFM9, a Mitochondrial CRM Protein, Is Crucial for Mitochondrial Intron Splicing, Mitochondria Function and Arabidopsis Growth and Stress Responses

Although the importance of chloroplast RNA splicing and ribosome maturation (CRM) domain-containing proteins has been established for chloroplast RNA metabolism and plant development, the functional role of CRM proteins in mitochondria remains largely unknown. Here, we investigated the role of a mitochondria-targeted CRM protein (At3g27550), named CFM9, in Arabidopsis thaliana. Confocal analysis revealed that CFM9 is localized in mitochondria. The cfm9 mutant exhibited delayed seed germination, retarded growth and shorter height compared with the wild type under normal conditions. The growth-defect phenotypes were more manifested upon high salinity, dehydration or ABA application. Complementation lines expressing CFM9 in the mutant background fully recovered the wild-type phenotypes. Notably, the mutant had abnormal mitochondria, increased hydrogen peroxide and reduced respiration activity, implying that CFM9 is indispensable for normal mitochondrial function. More important, the splicing of many intron-containing genes in mitochondria was defective in the mutant, suggesting that CFM9 plays a crucial role in the splicing of mitochondrial introns. Collectively, our results provide clear evidence emphasizing that CFM9 is an essential factor in the splicing of mitochondrial introns, which is crucial for mitochondrial biogenesis and function and the growth and development of Arabidopsis.

Elucidating genomic patterns and recombination events in plant cybrid mitochondria

KEY MESSAGE

Cybrid plant mitochondria undergo homologous recombination, mainly BIR, keep a single allele for each gene, and maintain exclusive sequences of each parent and a single copy of the homologous regions. The maintenance of a dynamic equilibrium between the mitochondrial and nuclear genomes requires continuous communication and a high level of compatibility between them, so that alterations in one genetic compartment need adjustments in the other. The co-evolution of nuclear and mitochondrial genomes has been poorly studied, even though the consequences and effects of this interaction are highly relevant for human health, as well as for crop improvement programs and for genetic engineering. The mitochondria of plants represent an excellent system to understand the mechanisms of genomic rearrangements, chimeric gene formation, incompatibility between nucleus and cytoplasm, and horizontal gene transfer. We carried out detailed analyses of the mtDNA of a repeated cybrid between the solanaceae Nicotiana tabacum and Hyoscyamus niger. The mtDNA of the cybrid was intermediate between the size of the parental mtDNAs and the sum of them. Noticeably, most of the homologous sequences inherited from both parents were lost. In contrast, the majority of the sequences exclusive of a single parent were maintained. The mitochondrial gene content included a majority of N. tabacum derived genes, but also chimeric, two-parent derived, and H. niger-derived genes in a tobacco nuclear background. Any of these alterations in the gene content could be the cause of CMS in the cybrid. The parental mtDNAs interacted through 28 homologous recombination events and a single case of illegitimate recombination. Three main homologous recombination mechanisms were recognized in the cybrid mitochondria. Break induced replication (BIR) pathway was the most frequent. We propose that BIR could be one of the mechanisms responsible for the loss of the majority of the repeated regions derived from H. niger.

The complete organelle genomes of Physochlaina orientalis: Insights into short sequence repeats across seed plant mitochondrial genomes

Short repeats (SR) play an important role in shaping seed plant mitochondrial genomes (mtDNAs). However, their origin, distribution, and relationships across the different plant lineages remain unresolved. We focus on the angiosperm family Solanaceae that shows great variation in repeat content and extend the study to a wide diversity of seed plants. We determined the complete nucleotide sequences of the organellar genomes of the medicinal plant Physochlaina orientalis (Solanaceae), member of the tribe Hyoscyameae. To understand the evolution of the P. orientalis mtDNA we made comparisons with those of five other Solanaceae. P. orientalis mtDNA presents the largest mitogenome (∼685 kb in size) among the Solanaceae and has an unprecedented 8-copy repeat family of ∼8.2 kb in length and a great number of SR arranged in tandem-like structures. We found that the SR in the Solanaceae share a common origin, but these only expanded in members of the tribe Hyoscyameae. We discuss a mechanism that could explain SR formation and expansion in P. orientalis and Hyoscyamus niger. Finally, the great increase in plant mitochondrial data allowed us to systematically extend our repeat analysis to a total of 136 seed plants to characterize and analyze for the first time families of SR among seed plant mtDNAs.

Repeats of Unusual Size in Plant Mitochondrial Genomes: Identification, Incidence and Evolution

Plant mitochondrial genomes have excessive size relative to coding capacity, a low mutation rate in genes and a high rearrangement rate. They also have abundant non-tandem repeats often including pairs of large repeats which cause isomerization of the genome by recombination, and numerous repeats of up to several hundred base pairs that recombine only when the genome is stressed by DNA damaging agents or mutations in DNA repair pathway genes. Early work on mitochondrial genomes led to the suggestion that repeats in the size range from several hundred to a few thousand base pair are underrepresented. The repeats themselves are not well-conserved between species, and are not always annotated in mitochondrial sequence assemblies. We systematically identified and compared these repeats, which are important clues to mechanisms of DNA maintenance in mitochondria. We developed a tool to find and curate non-tandem repeats larger than 50bp and analyzed the complete mitochondrial sequences from 157 plant species. We observed an interesting difference between taxa: the repeats are larger and more frequent in the vascular plants. Analysis of closely related species also shows that plant mitochondrial genomes evolve in dramatic bursts of breakage and rejoining, complete with DNA sequence gain and loss. We suggest an adaptive explanation for the existence of the repeats and their evolution.

Comparative study on mitogenomes of green tide algae

Since 2007, the annual green tide disaster in the Yellow Sea has brought serious economic losses to China. There is no research on the genetic similarities of four constituent species of green tide algae at the genomic level. We previously determined the mitochondrial genomes of Ulva prolifera, Ulva linza and Ulva flexuosa. In the present work, the mitochondrial genome of another green tide (Ulva compressa) was sequenced and analyzed. With the length of 62,311 bp, it contained 29 encoding genes, 26 tRNAs and 10 open reading frames. By comparing these four mitochondrial genomes, we found that U. compressa was quite different from the other three types of Ulva species. However, there were similarities between U. prolifera and U. linza in the number, distribution and homology of open reading frames, evolutionary and codon variation of tRNA, evolutionary relationship and selection pressure of coding genes. Repetitive sequence analysis of simple sequence repeats, tandem repeat and forward repeats further supposed that they have evolved from the same origin. In addition, we directly analyzed gene homologies and translocation of four green tide algae by Mauve alignment. There were gene order rearrangements among them. With fast-evolving genomes, these four green algal mitochondria have both conservatism and variation, thus opening another window for the understanding of origin and evolution of Ulva.

Characterization of genes specific to sua-CMS in Nicotiana tabacum

KEY MESSAGE

Six unique ORFs were characterized in tobacco plants with sua-CMS sterile cytoplasm, identifying the mtDNA basis for pollen sterility. sua-CMS (cytoplasmic male sterility), the most widely used sterile system in tobacco hybrids, is the only CMS type identified as having no negative effects on agronomic or quality traits in tobacco (Nicotiana tabacum) and as being fully male sterile. CMS is often associated with alterations of mitochondrial DNA (mtDNA), including novel chimeric open reading frames (ORFs), which result from rearrangement and recombination. Here, we obtained 34 mitochondrial ORFs in the sua-CMS line msZhongyan100 (sZY) by BLAST analysis. When we amplified these mitochondrial ORFs in seven tobacco CMS lines including sua-, glu-, rep-, rus-, tab1-, tab2-, and tab3-CMS types and in fertile tobacco, we found that six ORFs-orf82, orf103, orf115a, orf91, orf115b, and orf100-were located in three small regions (m-sr) of the mitochondrial genome of sZY and were unique to the sua-CMS line. We further amplified the m-sr fragments in three different backcross populations of the seven types of CMS, three F1 hybrids with sua-CMS sterile cytoplasm, two sua-CMS lines, and 284 fertile tobacco accessions. The ORFs were specific to plants with the sua-CMS background. All six unique ORFs were chimeric and had no homology with the mitochondrial genomes of fertile tobacco. Transcript analysis revealed that the ORFs were highly expressed in the anthers and floral buds of sZY. These six ORFs were specific to sua-CMS and could be used as molecular markers to identify sua-CMS lines, which is useful for improving breeding for heterosis in tobacco.

Complete sequence of kenaf (Hibiscus cannabinus) mitochondrial genome and comparative analysis with the mitochondrial genomes of other plants

Plant mitochondrial (mt) genomes are species specific due to the vast of foreign DNA migration and frequent recombination of repeated sequences. Sequencing of the mt genome of kenaf (Hibiscus cannabinus) is essential for elucidating its evolutionary characteristics. In the present study, single-molecule real-time sequencing technology (SMRT) was used to sequence the complete mt genome of kenaf. Results showed that the complete kenaf mt genome was 569,915 bp long and consisted of 62 genes, including 36 protein-coding, 3 rRNA and 23 tRNA genes. Twenty-five introns were found among nine of the 36 protein-coding genes, and five introns were trans-spliced. A comparative analysis with other plant mt genomes showed that four syntenic gene clusters were conserved in all plant mtDNAs. Fifteen chloroplast-derived fragments were strongly associated with mt genes, including the intact sequences of the chloroplast genes psaA, ndhB and rps7. According to the plant mt genome evolution analysis, some ribosomal protein genes and succinate dehydrogenase genes were frequently lost during the evolution of angiosperms. Our data suggest that the kenaf mt genome retained evolutionarily conserved characteristics. Overall, the complete sequencing of the kenaf mt genome provides additional information and enhances our better understanding of mt genomic evolution across angiosperms.

Assembly of the Mitochondrial Genome in the Campanulaceae Family Using Illumina Low-Coverage Sequencing

Platycodongrandiflorus (balloon flower) and Codonopsislanceolata (bonnet bellflower) are important herbs used in Asian traditional medicine, and both belong to the botanical family Campanulaceae. In this study, we designed and implemented a de novo DNA sequencing and assembly strategy to map the complete mitochondrial genomes of the first two members of the Campanulaceae using low-coverage Illumina DNA sequencing data. We produced a total of 28.9 Gb of paired-end sequencing data from the genomic DNA of P.grandiflorus (20.9 Gb) and C.lanceolata (8.0 Gb). The assembled mitochondrial genome of P.grandiflorus was found to consist of two circular chromosomes; the master circle contains 56 genes, and the minor circle contains 42 genes. The C.lanceolata mitochondrial genome consists of a single circle harboring 54 genes. Using a comparative genome structure and a pattern of repeated sequences, we show that the P.grandiflorus minor circle resulted from a recombination event involving the direct repeats of the master circle. Our dataset will be useful for comparative genomics and for evolutionary studies, and will facilitate further biological and phylogenetic characterization of species in the Campanulaceae.

Next generation sequencing reveals packaging of host RNAs by brome mosaic virus

Although RNA viruses evolved the mechanisms of specific encapsidation, miss-packaging of cellular RNAs has been reported in such RNA virus systems as flock house virus or cucumber necrosis virus. To find out if brome mosaic virus (BMV), a tripartite RNA virus, can package cellular RNAs, BMV was propagated in barley and in Nicotiana benthamiana hosts, purified by cesium chloride (CsCl) gradient ultracentrifugation followed by nuclease treatment to remove any contaminating cellular (host) RNAs. The extracted virion RNA was then sequenced by using next-generation sequencing (NGS RNA-Seq) with the Illumina protocol. Bioinformatic analysis revealed the content of host RNAs ranging from 0.07% for BMV extracted from barley to 0.10% for the virus extracted from N. benthamiana. The viruses from two sources appeared to co-encapsidate different patterns of host-RNAs, including ribosomal RNAs (rRNAs), messenger RNAs (mRNAs) but also mitochondrial and plastid RNAs and, interestingly, transposable elements, both transposons and retrotransposons. Our data reveal that BMV virions can carry host RNAs, having a potential to mediate horizontal gene transfer (HGT) in plants.

Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria

Our understanding of the dynamic and evolution of RNA editing in angiosperms is in part limited by the few editing sites identified to date. This study identified 10,217 editing sites from 17 diverse angiosperms. Our analyses confirmed the universality of certain features of RNA editing, and offer new evidence behind the loss of editing sites in angiosperms. RNA editing is a post-transcriptional process that substitutes cytidines (C) for uridines (U) in organellar transcripts of angiosperms. These substitutions mostly take place in mitochondrial messenger RNAs at specific positions called editing sites. By means of publicly available RNA-seq data, this study identified 10,217 editing sites in mitochondrial protein-coding genes of 17 diverse angiosperms. Even though other types of mismatches were also identified, we did not find evidence of non-canonical editing processes. The results showed an uneven distribution of editing sites among species, genes, and codon positions. The analyses revealed that editing sites were conserved across angiosperms but there were some species-specific sites. Non-synonymous editing sites were particularly highly conserved (~ 80%) across the plant species and were efficiently edited (80% editing extent). In contrast, editing sites at third codon positions were poorly conserved (~ 30%) and only partially edited (~ 40% editing extent). We found that the loss of editing sites along angiosperm evolution is mainly occurring by replacing editing sites with thymidines, instead of a degradation of the editing recognition motif around editing sites. Consecutive and highly conserved editing sites had been replaced by thymidines as result of retroprocessing, by which edited transcripts are reverse transcribed to cDNA and then integrated into the genome by homologous recombination. This phenomenon was more pronounced in eudicots, and in the gene cox1. These results suggest that retroprocessing is a widespread driving force underlying the loss of editing sites in angiosperm mitochondria.

Maize Dek37 Encodes a P-type PPR Protein That Affects cis-Splicing of Mitochondrial nad2 Intron 1 and Seed Development

Mitochondrial group II introns require the participation of numerous nucleus-encoded general and specific factors to achieve efficient splicing in vivo Pentatricopeptide repeat (PPR) proteins have been implicated in assisting group II intron splicing. Here, we identified and characterized a new maize seed mutant, defective kernel 37 (dek37), which has significantly delayed endosperm and embryo development. Dek37 encodes a classic P-type PPR protein that targets mitochondria. The dek37 mutation causes no detectable DEK37 protein in mutant seeds. Mitochondrial transcripts analysis indicated that dek37 mutation decreases splicing efficiency of mitochondrial nad2 intron 1, leading to reduced assembly and NADH dehydrogenase activity of complex I. Transmission Electron Microscopy (TEM) revealed severe morphological defects of mitochondria in dek37 Transcriptome analysis of dek37 endosperm indicated enhanced expression in the alternative respiratory pathway and extensive differentially expressed genes related to mitochondrial function. These results indicated that Dek37 is involved in cis-splicing of mitochondrial nad2 intron 1 and is required for complex I assembly, mitochondrial function, and seed development in maize.

Comparison of mitochondrial gene expression and polysome loading in different tobacco tissues

BACKGROUND

To investigate translational regulation of gene expression in plant mitochondria, a mitochondrial polysome isolation protocol was established for tobacco to investigate polysomal mRNA loading as a proxy for translational activity. Furthermore, we developed an oligonucleotide based microarray platform to determine the level of Nicotiana tabacum and Arabidopsis thaliana mitochondrial mRNA.

RESULTS

Microarray analysis of free and polysomal mRNAs was used to characterize differences in the levels of free transcripts and ribosome-bound mRNAs in various organs of tobacco plants. We have observed higher mitochondrial transcript levels in young leaves, flowers and floral buds as compared to fully expanded leaves and roots. A similar pattern of abundance was observed for ribosome-bound mitochondrial mRNAs in these tissues. However, the accumulation of the mitochondrial protein COX2 was found to be inversely related to that of its ribosome-bound mRNA.

CONCLUSIONS

Our results indicate that the association of mitochondrial mRNAs to ribosomes is largely determined by the total transcript level of a gene. However, at least for Cox2, we demonstrated that the level of ribosome-bound mRNA is not reflected by the amount of COX2 protein.

The complete mitochondrial genome sequences of potato (Solanum tuberosum L., Solanaceae)

Potato (Solanum tuberosum) from the Solanaceae is the fourth most important food crop worldwide. In this study, five complete mitochondrial genome sequences of S. tuberosum were characterized through de novo assembly of whole genome sequencing data. The resulting circular mitochondrial DNA molecules ranged from 49,171 bp to 297,014 bp in size and contained a total of 80 non-redundant genes, comprising 34 protein-coding genes, 24 hypothetical open reading frames, 19 tRNA genes, and 3 rRNA genes. Phylogenetic analysis using common protein-coding sequences confirmed that S. tuberosum belongs to the Solanoideae subfamily in the Solanaceae family.