Multiple losses and transfers to the nucleus of two mitochondrial succinate dehydrogenase genes during angiosperm evolution

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica

The family Schisandraceae is a basal angiosperm plant group distributed in East and Southeast Asia and includes many medicinal plant species such as Schisandra chinensis. In this study, mitochondrial genomes (mitogenomes) of two species, Schisandra repanda and Kadsura japonica, in the family were characterized through de novo assembly using sequencing data obtained with Oxford Nanopore and Illumina sequencing technologies. The mitogenomes of S. repanda were assembled into one circular contig (571,107 bp) and four linear contigs (10,898-607,430 bp), with a total of 60 genes: 38 protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The mitogenomes of K. japonica were assembled into five circular contigs (211,474-973,503 bp) and three linear contigs (8,010-72,712 bp), with a total of 66 genes: 44 PCGs, 19 tRNA genes, and 3 rRNA genes. The mitogenomes of the two species had complex structural features with high repeat numbers and chloroplast-derived sequences, as observed in other plant mitogenomes. Phylogenetic analysis based on PCGs revealed the taxonomical relationships of S. repanda and K. japonica with other species from Schisandraceae. Finally, molecular markers were developed to distinguish between S. repanda, K. japonica, and S. chinensis on the basis of InDel polymorphisms present in the mitogenomes. The mitogenomes of S. repanda and K. japonica will be valuable resources for molecular and taxonomic studies of plant species that belong to the family Schisandraceae.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

Native and foreign mitochondrial and nuclear encoded proteins conform the OXPHOS complexes of a holoparasitic plant

The intimate contact between the holoparasitic plant Lophophytum mirabile (Balanophoraceae) and its host plant (Fabaceae) facilitates the exchange of genetic information, increasing the frequency of horizontal gene transfer (HGT). Lophophytum stands out because it acquired a large number of mitochondrial genes (greater than 20) from its legume host that replaced the majority of the native homologs. These foreign genes code for proteins that form multisubunit enzyme complexes, such as those in the oxidative phosphorylation system (OXPHOS) and cytochrome c maturation (ccm) system, together with dozens of nuclear-encoded subunits. However, the existence and the origin of the nuclear subunits that form the major part of the OXPHOS and ccm system in Lophophytum remain unknown. It was proposed that nuclear-encoding genes whose products interact with foreign mitochondrial proteins are also foreign, minimizing the incompatibilities that could arise in the assembly and functioning of these multiprotein complexes. We identified a nearly complete set of OXPHOS and ccm system subunits evolving under selective constraints in the transcriptome of Lophophytum, indicating that OXPHOS is functional and resembles that of free-living angiosperms. Maximum Likelihood phylogenetic analyses revealed a single case of HGT in the nuclear genes, which results in mosaic OXPHOS and ccm system in Lophophytum. These observations raise new questions about the evolution and physiology of this parasitic plant. A putative case of cooperation between two foreign (one mitochondrial and one nuclear) genes is presented.

Discovery of succinate dehydrogenase candidate fungicides via lead optimization for effective resistance management of Fusarium oxysporum f. sp. capsici

Fusarium wilt of chili caused by the fungus Fusarium oxysporum f. sp. capsici (FCO) severely reduces the production of chili worldwide. There is growing evidence of resistance to commercial fungicides targeting succinate dehydrogenase (Sdh) of FCO soliciting the development of new Sdh inhibitors (SdhIs). In the current work, optimized docking and virtual screening were used to mine twelve SdhIs from the ZINC database, followed by in vitro antifungal evaluation on spore and radial mycelium development. Four new promising SdhIs exhibiting a mean mycelium inhibition rate greater than 85.6% (F = 155.8, P = 0.001, P < 0.05) were observed on ten strains of virulent and resistant FCO. Importantly, three of the discovered molecules exhibited potent spore germination inhibition (≥ 80%, P = 0.01, P < 0.05) compared to the commonly used fungicide penthiopyrad. A significant positive correlation (r* ≥ 0.67, P < 0.05) between the activities of the newly discovered SdhIs compared to penthiopyrad against all tested FCO strains indicated a broad-spectrum fungicidal activity. The current findings indicate that the four SdhI's discovered could judiciously replace certain commercial SdhIs that some FCO displays resistance to.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03157-8.

Complete mitogenomes of four Trichiurus species: A taxonomic review of the T.lepturus species complex

Four Trichiurus species, T.japonicus, T.lepturus, T.nanhaiensis, and T.brevis, from the coasts of the China Seas, have been identified and their entire mitochondrial genomes (mitogenomes) have been sequenced by next-generation sequencing technology. A comparative analysis of five mitogenomes was conducted, including the mitogenome of T.gangeticus. The mitogenomes contained 16.568-16.840 bp and encoded 36 typical mitochondrial genes (13 protein-coding, 2 ribosomal RNA-coding, and 21 transfer RNA-coding genes) and two typical noncoding control regions. Although tRNAPro is absent from Trichiurus mitogenomes, when compared with the 22 tRNAs reported in other vertebrates, the gene arrangements in the mitogenomes of the studied species are consistent with those in most teleost mitogenomes. The full-length sequences and protein-coding genes (PCGs) in the mitogenomes of the five species had obvious AT biases and negative GC skew values. Our study indicate that the specimens in the Indian Ocean are neither T.lepturus nor T.nanhaiensis but they are T.gangeticus; the Trichiurus species composition in the Indian Ocean is totally different from that in Pacific and Atlantic oceans; there are at least two Trichiurus species in Indian Ocean; and the worldwide systematics and diversity of the genus Trichiurus need to be reviewed.

Both Conifer II and Gnetales are characterized by a high frequency of ancient mitochondrial gene transfer to the nuclear genome

BACKGROUND

Mitochondrial gene transfer/loss is common in land plants, and therefore the fate of missing mitochondrial genes has attracted more and more attention. The gene content of gymnosperm mitochondria varies greatly, supplying a system for studying the evolutionary fate of missing mitochondrial genes.

RESULTS

Here, we studied the tempo and pattern of mitochondrial gene transfer/loss in gymnosperms represented by all 13 families, using high-throughput sequencing of both DNA and cDNA. All 41 mitochondrial protein-coding genes were found in cycads, Ginkgo and Pinaceae, whereas multiple mitochondrial genes were absent in Conifer II and Gnetales. In Conifer II, gene transfer from mitochondria to the nucleus followed by loss of the mitochondrial copy was common, but complete loss of a gene in both mitochondrial and nuclear genomes was rare. In contrast, both gene transfer and loss were commonly found in Gnetales. Notably, in Conifer II and Gnetales, the same five mitochondrial genes were transferred to the nuclear genome, and these gene transfer events occurred, respectively, in ancestors of the two lineages. A two-step transfer mechanism (retroprocessing and subsequent DNA-mediated gene transfer) may be responsible for mitochondrial gene transfer in Conifer II and Gnetales. Moreover, the mitochondrial gene content variation is correlated with gene length, GC content, hydrophobicity, and nucleotide substitution rates in land plants.

CONCLUSIONS

This study reveals a complete evolutionary scenario for variations of mitochondrial gene transferring in gymnosperms, and the factors responsible for mitochondrial gene content variation in land plants.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans

Background

The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered.

Results

We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome.

Conclusions

Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01035-y.

Extensive variation in nucleotide substitution rate and gene/intron loss in mitochondrial genomes of Pelargonium

Geraniaceae organelle genomes have been shown to exhibit several highly unusual features compared to most other photosynthetic angiosperms. This includes massively rearranged plastomes with considerable size variation, extensive gene and intron loss, accelerated rates of nucleotide substitutions in both mitogenomes and plastomes, and biparental inheritance and cytonuclear incompatibility of the plastome. Most previous studies have focused on plastome evolution with mitogenome comparisons limited to only a few taxa or genes. In this study, mitogenomes and transcriptomes were examined for 27 species of Geraniales, including 13 species of Pelargonium. Extensive gene and intron losses were detected across the Geraniales with Pelargonium representing the most gene depauperate lineage in the family. Plotting these events on the Geraniaceae phylogenetic tree showed that gene losses occurred multiple times, whereas intron losses more closely reflected the relationships among taxa. In addition, P. australe acquired an intron by horizontal transfer. Comparisons of nucleotide substitution rates in Pelargonium showed that synonymous changes in nuclear genes were much lower than in mitochondrial genes. This is in contrast to the previously published studies that indicated that nuclear genes have 16 fold higher rates than mitochondrial genes across angiosperms. Elevated synonymous substitutions occurred for each mitochondrial gene in Pelargonium with the highest values 783 and 324 times higher than outgroups and other Geraniaceae, respectively. Pelargonium is one of four unrelated genera of angiosperms (Ajuga, Plantago and Silene) that have experienced highly accelerated nucleotide substitutions in mitogenomes. It is distinct from most angiosperms in also having elevated substitution rates in plastid genes but the cause of rate accelerations in Pelargonium plastomes and mitogenomes may be different.

Mitochondrial Genome of Fagus sylvatica L. as a Source for Taxonomic Marker Development in the Fagales

European beech, Fagus sylvatica L., is one of the most important and widespread deciduous tree species in Central Europe and is widely managed for its hard wood. The complete DNA sequence of the mitochondrial genome of Fagus sylvatica L. was assembled and annotated based on Illumina MiSeq reads and validated using long reads from nanopore MinION sequencing. The genome assembled into a single DNA sequence of 504,715 bp in length containing 58 genes with predicted function, including 35 protein-coding, 20 tRNA and three rRNA genes. Additionally, 23 putative protein-coding genes were predicted supported by RNA-Seq data. Aiming at the development of taxon-specific mitochondrial genetic markers, the tool SNPtax was developed and applied to select genic SNPs potentially specific for different taxa within the Fagales. Further validation of a small SNP set resulted in the development of four CAPS markers specific for Fagus, Fagaceae, or Fagales, respectively, when considering over 100 individuals from a total of 69 species of deciduous trees and conifers from up to 15 families included in the marker validation. The CAPS marker set is suitable to identify the genus Fagus in DNA samples from tree tissues or wood products, including wood composite products.

The Mitochondrial Genome of Eleusine indica and Characterization of Gene Content within Poaceae

Plant mitochondrial (mt) genome assembly provides baseline data on size, structure, and gene content, but resolving the sequence of these large and complex organelle genomes remains challenging due to fragmentation, frequent recombination, and transfers of DNA from neighboring plastids. The mt genome for Eleusine indica (Poaceae: goosegrass) is comprehensibly analyzed here, providing key reference data for an economically significant invasive species that is also the maternal parent of the allotetraploid crop Finger millet (Eleusine coracana). The assembled E. indica genome contains 33 protein coding genes, 6 rRNA subunits, 24 tRNA, 8 large repetitive regions 15 kb of transposable elements across a total of 520,691 bp. Evidence of RNA editing and loss of rpl2, rpl5, rps14, rps11, sdh4, and sdh3 genes is evaluated in the context of an updated survey of mt genomic gene content across the grasses through an analysis of publicly available data. Hypothesized patterns of Poaceae mt gene loss are examined in a phylogenetic context to clarify timing, showing that rpl2 was transferred to the nucleus from the mitochondrion prior to the origin of the PACMAD clade.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

Fluctuations in Fabaceae mitochondrial genome size and content are both ancient and recent

BACKGROUND

Organelle genome studies of Fabaceae, an economically and ecologically important plant family, have been biased towards the plastid genome (plastome). Thus far, less than 15 mitochondrial genome (mitogenome) sequences of Fabaceae have been published, all but four of which belong to the subfamily Papilionoideae, limiting the understanding of size variation and content across the family. To address this, four mitogenomes were sequenced and assembled from three different subfamilies (Cercidoideae, Detarioideae and Caesalpinioideae).

RESULTS

Phylogenetic analysis based on shared mitochondrial protein coding regions produced a fully resolved and well-supported phylogeny that was completely congruent with the plastome tree. Comparative analyses suggest that two kinds of mitogenome expansions have occurred in Fabaceae. Size expansion of four genera (Tamarindus, Libidibia, Haematoxylum, and Leucaena) in two subfamilies (Detarioideae and Caesalpinioideae) occurred in relatively deep nodes, and was mainly caused by intercellular gene transfer and/or interspecific horizontal gene transfer (HGT). The second, more recent expansion occurred in the Papilionoideae as a result of duplication of native mitochondrial sequences. Family-wide gene content analysis revealed 11 gene losses, four (rps2, 7, 11 and 13) of which occurred in the ancestor of Fabaceae. Losses of the remaining seven genes (cox2, rpl2, rpl10, rps1, rps19, sdh3, sdh4) were restricted to specific lineages or occurred independently in different clades. Introns of three genes (cox2, ccmFc and rps10) showed extensive lineage-specific length variation due to large sequence insertions and deletions. Shared DNA analysis among Fabaceae mitogenomes demonstrated a substantial decay of intergenic spacers and provided further insight into HGT between the mimosoid clade of Caesalpinioideae and the holoparasitic Lophophytum (Balanophoraceae).

CONCLUSION

This study represents the most exhaustive analysis of Fabaceae mitogenomes so far, and extends the understanding the dynamic variation in size and gene/intron content. The four newly sequenced mitogenomes reported here expands the phylogenetic coverage to four subfamilies. The family has experienced multiple mitogenome size fluctuations in both ancient and recent times. The causes of these size variations are distinct in different lineages. Fabaceae mitogenomes experienced extensive size fluctuation by recruitment of exogenous DNA and duplication of native mitochondrial DNA.

Mitochondrial complex II of plants: subunit composition, assembly, and function in respiration and signaling

Complex II [succinate dehydrogenase (succinate-ubiquinone oxidoreductase); EC 1.3.5.1; SDH] is the only enzyme shared by both the electron transport chain and the tricarboxylic acid (TCA) cycle in mitochondria. Complex II in plants is considered unusual because of its accessory subunits (SDH5-SDH8), in addition to the catalytic subunits of SDH found in all eukaryotes (SDH1-SDH4). Here, we review compositional and phylogenetic analysis and biochemical dissection studies to both clarify the presence and propose a role for these subunits. We also consider the wider functional and phylogenetic evidence for SDH assembly factors and the reports from plants on the control of SDH1 flavination and SDH1-SDH2 interaction. Plant complex II has been shown to influence stomatal opening, the plant defense response and reactive oxygen species-dependent stress responses. Signaling molecules such as salicyclic acid (SA) and nitric oxide (NO) are also reported to interact with the ubiquinone (UQ) binding site of SDH, influencing signaling transduction in plants. Future directions for SDH research in plants and the specific roles of its different subunits and assembly factors are suggested, including the potential for reverse electron transport to explain the succinate-dependent production of reactive oxygen species in plants and new avenues to explore the evolution of plant mitochondrial complex II and its utility.

PacBio-Based Mitochondrial Genome Assembly of Leucaena trichandra (Leguminosae) and an Intrageneric Assessment of Mitochondrial RNA Editing

Reconstructions of vascular plant mitochondrial genomes (mt-genomes) are notoriously complicated by rampant recombination that has resulted in comparatively few plant mt-genomes being available. The dearth of plant mitochondrial resources has limited our understanding of mt-genome structural diversity, complex patterns of RNA editing, and the origins of novel mt-genome elements. Here, we use an efficient long read (PacBio) iterative assembly pipeline to generate mt-genome assemblies for Leucaena trichandra (Leguminosae: Caesalpinioideae: mimosoid clade), providing the first assessment of non-papilionoid legume mt-genome content and structure to date. The efficiency of the assembly approach facilitated the exploration of alternative structures that are common place among plant mitochondrial genomes. A compact version (729 kbp) of the recovered assemblies was used to investigate sources of mt-genome size variation among legumes and mt-genome sequence similarity to the legume associated root holoparasite Lophophytum. The genome and an associated suite of transcriptome data from select species of Leucaena permitted an in-depth exploration of RNA editing in a diverse clade of closely related species that includes hybrid lineages. RNA editing in the allotetraploid, Leucaena leucocephala, is consistent with co-option of nearly equal maternal and paternal C-to-U edit components, generating novel combinations of RNA edited sites. A preliminary investigation of L. leucocephala C-to-U edit frequencies identified the potential for a hybrid to generate unique pools of alleles from parental variation through edit frequencies shared with one parental lineage, those intermediate between parents, and transgressive patterns.

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.

The RNA Editing Factor SlORRM4 Is Required for Normal Fruit Ripening in Tomato

RNA editing plays a key posttranscriptional role in gene expression. Existing studies on cytidine-to-uridine RNA editing in plants have focused on maize (Zea mays), rice (Oryza sativa), and Arabidopsis (Arabidopsis thaliana). However, the importance and regulation of RNA editing in several critical agronomic processes are not well understood, a notable example of which is fruit ripening. Here, we analyzed the expression profile of 33 RNA editing factors and identified 11 putative tomato (Solanum lycopersicum) fruit ripening-related factors. A rapid virus-induced gene silencing assay indicated that the organelle RNA recognition motif-containing protein SlORRM4 affected tomato fruit ripening. Knocking out SlORRM4 expression using a clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 genome editing strategy delayed tomato fruit ripening by lowering respiratory rate and ethylene production. Additionally, the expression of numerous genes associated with fruit ripening and mitochondrial functions changed significantly when SlORRM4 was knocked out. Moreover, the loss of SlORRM4 function significantly reduced RNA editing of many mitochondrial transcripts, leading to low-level expression of some core subunits that are critical for mitochondrial complex assembly (i.e. Nad3, Cytc1, and COX II). Taken together, these results indicate that SlORRM4 is involved in RNA editing of transcripts in ripening fruit that influence mitochondrial function and key aspects of fruit ripening.

Mitochondrial genome evolution in Alismatales: Size reduction and extensive loss of ribosomal protein genes

The order Alismatales is a hotspot for evolution of plant mitochondrial genomes characterized by remarkable differences in genome size, substitution rates, RNA editing, retrotranscription, gene loss and intron loss. Here we have sequenced the complete mitogenomes of Zostera marina and Stratiotes aloides, which together with previously sequenced mitogenomes from Butomus and Spirodela, provide new evolutionary evidence of genome size reduction, gene loss and transfer to the nucleus. The Zostera mitogenome includes a large portion of DNA transferred from the plastome, yet it is the smallest known mitogenome from a non-parasitic plant. Using a broad sample of the Alismatales, the evolutionary history of ribosomal protein gene loss is analyzed. In Zostera almost all ribosomal protein genes are lost from the mitogenome, but only some can be found in the nucleus.

Assembly and comparative analysis of complete mitochondrial genome sequence of an economic plant Salix suchowensis

Willow is a widely used dioecious woody plant of Salicaceae family in China. Due to their high biomass yields, willows are promising sources for bioenergy crops. In this study, we assembled the complete mitochondrial (mt) genome sequence of S. suchowensis with the length of 644,437 bp using Roche-454 GS FLX Titanium sequencing technologies. Base composition of the S. suchowensis mt genome is A (27.43%), T (27.59%), C (22.34%), and G (22.64%), which shows a prevalent GC content with that of other angiosperms. This long circular mt genome encodes 58 unique genes (32 protein-coding genes, 23 tRNA genes and 3 rRNA genes), and 9 of the 32 protein-coding genes contain 17 introns. Through the phylogenetic analysis of 35 species based on 23 protein-coding genes, it is supported that Salix as a sister to Populus. With the detailed phylogenetic information and the identification of phylogenetic position, some ribosomal protein genes and succinate dehydrogenase genes are found usually lost during evolution. As a native shrub willow species, this worthwhile research of S. suchowensis mt genome will provide more desirable information for better understanding the genomic breeding and missing pieces of sex determination evolution in the future.

SDH6 and SDH7 Contribute to Anchoring Succinate Dehydrogenase to the Inner Mitochondrial Membrane in Arabidopsis thaliana

The succinate dehydrogenase complex (complex II) is a highly conserved protein complex composed of the SDH1 to SDH4 subunits in bacteria and in the mitochondria of animals and fungi. The reason for the occurrence of up to four additional subunits in complex II of plants, termed SDH5 to SDH8, so far is a mystery. Here, we present a biochemical approach to investigate the internal subunit arrangement of Arabidopsis (Arabidopsis thaliana) complex II. Using low-concentration detergent treatments, the holo complex is dissected into subcomplexes that are analyzed by a three-dimensional gel electrophoresis system. Protein identifications by mass spectrometry revealed that the largest subcomplex (IIa) represents the succinate dehydrogenase domain composed of SDH1 and SDH2. Another subcomplex (IIb) is composed of the SDH3, SDH4, SDH6, and SDH7 subunits. All four proteins include transmembrane helices and together form the membrane anchor of complex II. Sequence analysis revealed that SDH3 and SDH4 lack helices conserved in other organisms. Using homology modeling and phylogenetic analyses, we present evidence that SDH6 and SDH7 substitute missing sequence stretches of SDH3 and SDH4 in plants. Together with SDH5, which is liberated upon dissection of complex II into subcomplexes, SDH6 and SDH7 also add some hydrophilic mass to plant complex II, which possibly inserts further functions into this smallest protein complex of the oxidative phosphorylation system (which is not so small in plants).

Most Compositae (Asteraceae) are descendants of a paleohexaploid and all share a paleotetraploid ancestor with the Calyceraceae

PREMISE OF THE STUDY

Like many other flowering plants, members of the Compositae (Asteraceae) have a polyploid ancestry. Previous analyses found evidence for an ancient duplication or possibly triplication in the early evolutionary history of the family. We sought to better place this paleopolyploidy in the phylogeny and assess its nature.

METHODS

We sequenced new transcriptomes for Barnadesia, the lineage sister to all other Compositae, and four representatives of closely related families. Using a recently developed algorithm, MAPS, we analyzed nuclear gene family phylogenies for evidence of paleopolyploidy.

KEY RESULTS

We found that the previously recognized Compositae paleopolyploidy is also in the ancestry of the Calyceraceae. Our phylogenomic analyses uncovered evidence for a successive second round of genome duplication among all sampled Compositae except Barnadesia.

CONCLUSIONS

Our analyses of new samples with new tools provide a revised view of paleopolyploidy in the Compositae. Together with results from a high density Lactuca linkage map, our results suggest that the Compositae and Calyceraceae have a common paleotetraploid ancestor and that most Compositae are descendants of a paleohexaploid. Although paleohexaploids have been previously identified, this is the first example where the paleotetraploid and paleohexaploid lineages have survived over tens of millions of years. The complex polyploidy in the ancestry of the Compositae and Calyceraceae represents a unique opportunity to study the long-term evolutionary fates and consequences of different ploidal levels.

The Whole Genome Assembly and Comparative Genomic Research of Thellungiella parvula (Extremophile Crucifer) Mitochondrion

The complete nucleotide sequences of the mitochondrial (mt) genome of an extremophile species Thellungiella parvula (T. parvula) have been determined with the lengths of 255,773 bp. T. parvula mt genome is a circular sequence and contains 32 protein-coding genes, 19 tRNA genes, and three ribosomal RNA genes with a 11.5% coding sequence. The base composition of 27.5% A, 27.5% T, 22.7% C, and 22.3% G in descending order shows a slight bias of 55% AT. Fifty-three repeats were identified in the mitochondrial genome of T. parvula, including 24 direct repeats, 28 tandem repeats (TRs), and one palindromic repeat. Furthermore, a total of 199 perfect microsatellites have been mined with a high A/T content (83.1%) through simple sequence repeat (SSR) analysis and they were distributed unevenly within this mitochondrial genome. We also analyzed other plant mitochondrial genomes' evolution in general, providing clues for the understanding of the evolution of organelles genomes in plants. Comparing with other Brassicaceae species, T. parvula is related to Arabidopsis thaliana whose characters of low temperature resistance have been well documented. This study will provide important genetic tools for other Brassicaceae species research and improve yields of economically important plants.

Genome Sequences of Populus tremula Chloroplast and Mitochondrion: Implications for Holistic Poplar Breeding

Complete Populus genome sequences are available for the nucleus (P. trichocarpa; section Tacamahaca) and for chloroplasts (seven species), but not for mitochondria. Here, we provide the complete genome sequences of the chloroplast and the mitochondrion for the clones P. tremula W52 and P. tremula x P. alba 717-1B4 (section Populus). The organization of the chloroplast genomes of both Populus clones is described. A phylogenetic tree constructed from all available complete chloroplast DNA sequences of Populus was not congruent with the assignment of the related species to different Populus sections. In total, 3,024 variable nucleotide positions were identified among all compared Populus chloroplast DNA sequences. The 5-prime part of the LSC from trnH to atpA showed the highest frequency of variations. The variable positions included 163 positions with SNPs allowing for differentiating the two clones with P. tremula chloroplast genomes (W52, 717-1B4) from the other seven Populus individuals. These potential P. tremula-specific SNPs were displayed as a whole-plastome barcode on the P. tremula W52 chloroplast DNA sequence. Three of these SNPs and one InDel in the trnH-psbA linker were successfully validated by Sanger sequencing in an extended set of Populus individuals. The complete mitochondrial genome sequence of P. tremula is the first in the family of Salicaceae. The mitochondrial genomes of the two clones are 783,442 bp (W52) and 783,513 bp (717-1B4) in size, structurally very similar and organized as single circles. DNA sequence regions with high similarity to the W52 chloroplast sequence account for about 2% of the W52 mitochondrial genome. The mean SNP frequency was found to be nearly six fold higher in the chloroplast than in the mitochondrial genome when comparing 717-1B4 with W52. The availability of the genomic information of all three DNA-containing cell organelles will allow a holistic approach in poplar molecular breeding in the future.

Rapid evolutionary divergence of Gossypium barbadense and G. hirsutum mitochondrial genomes

Background

The mitochondrial genome from upland cotton, G. hirsutum, was previously sequenced. To elucidate the evolution of mitochondrial genomic diversity within a single genus, we sequenced the mitochondrial genome from Sea Island cotton (Gossypium barbadense L.).

Methods

Mitochondrial DNA from week-old etiolated seedlings was extracted from isolated organelles using discontinuous sucrose density gradient method. Mitochondrial genome was sequenced with Solexa using paired-end, 90 bp read. The clean reads were assembled into contigs using ABySS and finished via additional fosmid and BAC sequencing. Finally, the genome was annotated and analyzed using different softwares.

Results

The G. barbadense (Sea Island cotton) mitochondrial genome was fully sequenced (677,434-bp) and compared to the mitogenome of upland cotton. The G. barbadense mitochondrial DNA contains seven more genes than that of upland cotton, with a total of 40 protein coding genes (excluding possible pseudogenes), 6 rRNA genes, and 29 tRNA genes. Of these 75 genes, atp1, mttB, nad4, nad9, rrn5, rrn18, and trnD(GTC)-cp were each represented by two identical copies. A single 64 kb repeat was largely responsible for the 9 % difference in genome size between the two mtDNAs. Comparison of genome structures between the two mitochondrial genomes revealed 8 rearranged syntenic regions and several large repeats. The largest repeat was missing from the master chromosome in G. hirsutum. Both mitochondrial genomes contain a duplicated copy of rps3 (rps3-2) in conjunction with a duplication of repeated sequences. Phylogenetic and divergence considerations suggest that a 544-bp fragment of rps3 was transferred to the nuclear genome shortly after divergence of the A- and D- genome diploid cottons.

Conclusion

These results highlight the insights to the evolution of structural variation between Sea Island and upland cotton mitochondrial genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1988-0) contains supplementary material, which is available to authorized users.

Parallel Structural Evolution of Mitochondrial Ribosomes and OXPHOS Complexes

The five macromolecular complexes that jointly mediate oxidative phosphorylation (OXPHOS) in mitochondria consist of many more subunits than those of bacteria, yet, it remains unclear by which evolutionary mechanism(s) these novel subunits were recruited. Even less well understood is the structural evolution of mitochondrial ribosomes (mitoribosomes): while it was long thought that their exceptionally high protein content would physically compensate for their uniquely low amount of ribosomal RNA (rRNA), this hypothesis has been refuted by structural studies. Here, we present a cryo-electron microscopy structure of the 73S mitoribosome from Neurospora crassa, together with genomic and proteomic analyses of mitoribosome composition across the eukaryotic domain. Surprisingly, our findings reveal that both structurally and compositionally, mitoribosomes have evolved very similarly to mitochondrial OXPHOS complexes via two distinct phases: A constructive phase that mainly acted early in eukaryote evolution, resulting in the recruitment of altogether approximately 75 novel subunits, and a reductive phase that acted during metazoan evolution, resulting in gradual length-reduction of mitochondrially encoded rRNAs and OXPHOS proteins. Both phases can be well explained by the accumulation of (slightly) deleterious mutations and deletions, respectively, in mitochondrially encoded rRNAs and OXPHOS proteins. We argue that the main role of the newly recruited (nuclear encoded) ribosomal- and OXPHOS proteins is to provide structural compensation to the mutationally destabilized mitochondrially encoded components. While the newly recruited proteins probably provide a selective advantage owing to their compensatory nature, and while their presence may have opened evolutionary pathways toward novel mitochondrion-specific functions, we emphasize that the initial events that resulted in their recruitment was nonadaptive in nature. Our framework is supported by population genetic studies, and it can explain the complete structural evolution of mitochondrial ribosomes and OXPHOS complexes, as well as many observed functions of individual proteins.

Complete sequences of organelle genomes from the medicinal plant Rhazya stricta (Apocynaceae) and contrasting patterns of mitochondrial genome evolution across asterids

Background

Rhazya stricta is native to arid regions in South Asia and the Middle East and is used extensively in folk medicine to treat a wide range of diseases. In addition to generating genomic resources for this medicinally important plant, analyses of the complete plastid and mitochondrial genomes and a nuclear transcriptome from Rhazya provide insights into inter-compartmental transfers between genomes and the patterns of evolution among eight asterid mitochondrial genomes.

Results

The 154,841 bp plastid genome is highly conserved with gene content and order identical to the ancestral organization of angiosperms. The 548,608 bp mitochondrial genome exhibits a number of phenomena including the presence of recombinogenic repeats that generate a multipartite organization, transferred DNA from the plastid and nuclear genomes, and bidirectional DNA transfers between the mitochondrion and the nucleus. The mitochondrial genes sdh3 and rps14 have been transferred to the nucleus and have acquired targeting presequences. In the case of rps14, two copies are present in the nucleus; only one has a mitochondrial targeting presequence and may be functional. Phylogenetic analyses of both nuclear and mitochondrial copies of rps14 across angiosperms suggests Rhazya has experienced a single transfer of this gene to the nucleus, followed by a duplication event. Furthermore, the phylogenetic distribution of gene losses and the high level of sequence divergence in targeting presequences suggest multiple, independent transfers of both sdh3 and rps14 across asterids. Comparative analyses of mitochondrial genomes of eight sequenced asterids indicates a complicated evolutionary history in this large angiosperm clade with considerable diversity in genome organization and size, repeat, gene and intron content, and amount of foreign DNA from the plastid and nuclear genomes.

Conclusions

Organelle genomes of Rhazya stricta provide valuable information for improving the understanding of mitochondrial genome evolution among angiosperms. The genomic data have enabled a rigorous examination of the gene transfer events. Rhazya is unique among the eight sequenced asterids in the types of events that have shaped the evolution of its mitochondrial genome. Furthermore, the organelle genomes of R. stricta provide valuable genomic resources for utilizing this important medicinal plant in biotechnology applications.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-405) contains supplementary material, which is available to authorized users.

In vitro import and assembly of the nucleus-encoded mitochondrial subunit III of cytochrome c oxidase (Cox3)

The cox3 gene, encoding subunit III of cytochrome c oxidase (Cox3) is in mitochondrial genomes except in chlorophycean algae, where it is localized in the nucleus. Therefore, algae like Chlamydomonas reinhardtii, Polytomella sp. and Volvox carteri, synthesize the Cox3 polypeptide in the cytosol, import it into mitochondria, and integrate it into the cytochrome c oxidase complex. In this work, we followed the in vitro internalization of the Cox3 precursor by isolated, import-competent mitochondria of Polytomella sp. In this colorless alga, the precursor Cox3 protein is synthesized with a long, cleavable, N-terminal mitochondrial targeting sequence (MTS) of 98 residues. In an import time course, a transient Cox3 intermediate was identified, suggesting that the long MTS is processed more than once. The first processing step is sensitive to the metalo-protease inhibitor 1,10-ortophenantroline, suggesting that it is probably carried out by the matrix-located Mitochondrial Processing Protease. Cox3 is readily imported through an energy-dependent import pathway and integrated into the inner mitochondrial membrane, becoming resistant to carbonate extraction. Furthermore, the imported Cox3 protein was assembled into cytochrome c oxidase, as judged by the presence of a labeled band co-migrating with complex IV in Blue Native Electrophoresis. A model for the biogenesis of Cox3 in chlorophycean algae is proposed. This is the first time that the in vitro mitochondrial import of a cytosol-synthesized Cox3 subunit is described.

Succinate dehydrogenase: the complex roles of a simple enzyme

Succinate dehydrogenase (SDH) oxidises succinate to fumarate as a component of the tricarboxylic acid cycle and ubiquinone to ubiquinol in the mitochondrial electron transport chain. Studies of SDH mutants have revealed far-reaching effects of altering succinate oxidation in plant cells. The plant SDH complex composition, structure and assembly are all beginning to be understood but the implications of the divergence across eukaryotes is still unclear. We propose an integration of the reported physiological roles of SDH in plants which influence photosynthesis, the function of stomata, root elongation and fungal defence. Future SDH research needed in plants should involve tissue-specific studies of mutants, analysis of the pathways induced by succinate-dependent reactive oxygen species generation and assessment of the impact of succinate accumulation on metabolism.

The complete sequence of the mitochondrial genome of Butomus umbellatus--a member of an early branching lineage of monocotyledons

In order to study the evolution of mitochondrial genomes in the early branching lineages of the monocotyledons, i.e., the Acorales and Alismatales, we are sequencing complete genomes from a suite of key taxa. As a starting point the present paper describes the mitochondrial genome of Butomus umbellatus (Butomaceae) based on next-generation sequencing data. The genome was assembled into a circular molecule, 450,826 bp in length. Coding sequences cover only 8.2% of the genome and include 28 protein coding genes, four rRNA genes, and 12 tRNA genes. Some of the tRNA genes and a 16S rRNA gene are transferred from the plastid genome. However, the total amount of recognized plastid sequences in the mitochondrial genome is only 1.5% and the amount of DNA transferred from the nucleus is also low. RNA editing is abundant and a total of 557 edited sites are predicted in the protein coding genes. Compared to the 40 angiosperm mitochondrial genomes sequenced to date, the GC content of the Butomus genome is uniquely high (49.1%). The overall similarity between the mitochondrial genomes of Butomus and Spirodela (Araceae), the closest relative yet sequenced, is low (less than 20%), and the two genomes differ in size by a factor 2. Gene order is also largely unconserved. However, based on its phylogenetic position within the core alismatids Butomus will serve as a good reference point for subsequent studies in the early branching lineages of the monocotyledons.

The mitochondrial genome of soybean reveals complex genome structures and gene evolution at intercellular and phylogenetic levels

Determining mitochondrial genomes is important for elucidating vital activities of seed plants. Mitochondrial genomes are specific to each plant species because of their variable size, complex structures and patterns of gene losses and gains during evolution. This complexity has made research on the soybean mitochondrial genome difficult compared with its nuclear and chloroplast genomes. The present study helps to solve a 30-year mystery regarding the most complex mitochondrial genome structure, showing that pairwise rearrangements among the many large repeats may produce an enriched molecular pool of 760 circles in seed plants. The soybean mitochondrial genome harbors 58 genes of known function in addition to 52 predicted open reading frames of unknown function. The genome contains sequences of multiple identifiable origins, including 6.8 kb and 7.1 kb DNA fragments that have been transferred from the nuclear and chloroplast genomes, respectively, and some horizontal DNA transfers. The soybean mitochondrial genome has lost 16 genes, including nine protein-coding genes and seven tRNA genes; however, it has acquired five chloroplast-derived genes during evolution. Four tRNA genes, common among the three genomes, are derived from the chloroplast. Sizeable DNA transfers to the nucleus, with pericentromeric regions as hotspots, are observed, including DNA transfers of 125.0 kb and 151.6 kb identified unambiguously from the soybean mitochondrial and chloroplast genomes, respectively. The soybean nuclear genome has acquired five genes from its mitochondrial genome. These results provide biological insights into the mitochondrial genome of seed plants, and are especially helpful for deciphering vital activities in soybean.

Succinate dehydrogenase assembly factor 2 is needed for assembly and activity of mitochondrial complex II and for normal root elongation in Arabidopsis

Mitochondria complex II (succinate dehydrogenase, SDH) plays a central role in respiratory metabolism as a component of both the electron transport chain and the tricarboxylic acid cycle. We report the identification of an SDH assembly factor by analysis of T-DNA insertions in At5g51040, a protein with unknown function that was identified by mass spectrometry analysis as a low abundance mitochondrial protein. This gene is co-expressed with a number of genes encoding mitochondrial proteins, including SDH1-1, and has low partial sequence similarity to human SDHAF2, a protein required for flavin-adenine dinucleotide (FAD) insertion into SDH. In contrast to observations of other SDH deficient lines in Arabidopsis, the sdhaf2 line did not affect photosynthetic rate or stomatal conductance, but instead showed inhibition of primary root elongation with early lateral root emergence, presumably due to the low SDH activity caused by the reduced abundance of SDHAF2. Both roots and leaves showed succinate accumulation but different responses in the abundance of other organic acids and amino acids assayed. Isolated mitochondria showed lowered SDH1 protein abundance, lowered maximal SDH activity and less protein-bound flavin-adenine dinucleotide (FAD) at the molecular mass of SDH1 in the gel separation. The short root phenotype and SDH function of sdhaf2 was fully complemented by transformation with SDHAF2. Application of the SDH inhibitor, malonate, phenocopied the sdhaf2 root architecture in WT. Whole root respiratory assays showed no difference between WT and sdhaf2, but micro-respirometry of the tips of roots clearly showed low oxygen consumption in sdhaf2 which could explain a metabolic deficit responsible for root tip growth.

The mitochondrial genome of an aquatic plant, Spirodela polyrhiza

BACKGROUND

Spirodela polyrhiza is a species of the order Alismatales, which represent the basal lineage of monocots with more ancestral features than the Poales. Its complete sequence of the mitochondrial (mt) genome could provide clues for the understanding of the evolution of mt genomes in plant.

METHODS

Spirodela polyrhiza mt genome was sequenced from total genomic DNA without physical separation of chloroplast and nuclear DNA using the SOLiD platform. Using a genome copy number sensitive assembly algorithm, the mt genome was successfully assembled. Gap closure and accuracy was determined with PCR products sequenced with the dideoxy method.

CONCLUSIONS

This is the most compact monocot mitochondrial genome with 228,493 bp. A total of 57 genes encode 35 known proteins, 3 ribosomal RNAs, and 19 tRNAs that recognize 15 amino acids. There are about 600 RNA editing sites predicted and three lineage specific protein-coding-gene losses. The mitochondrial genes, pseudogenes, and other hypothetical genes (ORFs) cover 71,783 bp (31.0%) of the genome. Imported plastid DNA accounts for an additional 9,295 bp (4.1%) of the mitochondrial DNA. Absence of transposable element sequences suggests that very few nuclear sequences have migrated into Spirodela mtDNA. Phylogenetic analysis of conserved protein-coding genes suggests that Spirodela shares the common ancestor with other monocots, but there is no obvious synteny between Spirodela and rice mtDNAs. After eliminating genes, introns, ORFs, and plastid-derived DNA, nearly four-fifths of the Spirodela mitochondrial genome is of unknown origin and function. Although it contains a similar chloroplast DNA content and range of RNA editing as other monocots, it is void of nuclear insertions, active gene loss, and comprises large regions of sequences of unknown origin in non-coding regions. Moreover, the lack of synteny with known mitochondrial genomic sequences shed new light on the early evolution of monocot mitochondrial genomes.

The cytosol-synthesized subunit II (Cox2) precursor with the point mutation W56R is correctly processed in yeast mitochondria to rescue cytochrome oxidase

Deletion of the yeast mitochondrial gene COX2 encoding subunit 2 (Cox2) of cytochrome c oxidase (CcO) results in loss of respiration (Δcox2 strain). Supekova et al. (2010) [1] transformed a Δcox2 strain with a vector expressing Cox2 with a mitochondrial targeting sequence (MTS) and the point mutation W56R (Cox2(W56R)), restoring respiratory growth. Here, the CcO carrying the allotopically-expressed Cox2(W56R) was characterized. Yeast mitochondria from the wild-type (WT) and the Δcox2+Cox2(W56R) strains were subjected to Blue Native electrophoresis. In-gel activity of CcO and spectroscopic quantitation of cytochromes revealed that only 60% of CcO is present in the complemented strain, and that less CcO is found associated in supercomplexes as compared to WT. CcOs from the WT and the mutant exhibited similar subunit composition, although activity was 20-25% lower in the enzyme containing Cox2(W56R) than in the one with Cox2(WT). Tandem mass spectrometry confirmed that W(56) was substituted by R(56) in Cox2(W56R). In addition, Cox2(W56R) exhibited the same N-terminus than Cox2(WT), indicating that the MTS of Oxa1 and the leader sequence of 15 residues were removed from Cox2(W56R) during maturation. Thus, Cox2(W56R) is identical to Cox2(WT) except for the point mutation W56R. Mitochondrial Cox1 synthesis is strongly reduced in Δcox2 mutants, but the Cox2(W56R) complemented strain led to full restoration of Cox1 synthesis. We conclude that the cytosol-synthesized Cox2(W56R) follows a rate-limiting process of import, maturation or assembly that yields lower steady-state levels of CcO. Still, the allotopically-expressed Cox2(W56R) restores CcO activity and allows mitochondrial Cox1 synthesis to advance at WT levels.

The mitochondrial genome of Malus domestica and the import-driven hypothesis of mitochondrial genome expansion in seed plants

Mitochondrial genomes of spermatophytes are the largest of all organellar genomes. Their large size has been attributed to various factors; however, the relative contribution of these factors to mitochondrial DNA (mtDNA) expansion remains undetermined. We estimated their relative contribution in Malus domestica (apple). The mitochondrial genome of apple has a size of 396 947 bp and a one to nine ratio of coding to non-coding DNA, close to the corresponding average values for angiosperms. We determined that 71.5% of the apple mtDNA sequence was highly similar to sequences of its nuclear DNA. Using nuclear gene exons, nuclear transposable elements and chloroplast DNA as markers of promiscuous DNA content in mtDNA, we estimated that approximately 20% of the apple mtDNA consisted of DNA sequences imported from other cell compartments, mostly from the nucleus. Similar marker-based estimates of promiscuous DNA content in the mitochondrial genomes of other species ranged between 21.2 and 25.3% of the total mtDNA length for grape, between 23.1 and 38.6% for rice, and between 47.1 and 78.4% for maize. All these estimates are conservative, because they underestimate the import of non-functional DNA. We propose that the import of promiscuous DNA is a core mechanism for mtDNA size expansion in seed plants. In apple, maize and grape this mechanism contributed far more to genome expansion than did homologous recombination. In rice the estimated contribution of both mechanisms was found to be similar.

De novo assembly of the carrot mitochondrial genome using next generation sequencing of whole genomic DNA provides first evidence of DNA transfer into an angiosperm plastid genome

BACKGROUND

Sequence analysis of organelle genomes has revealed important aspects of plant cell evolution. The scope of this study was to develop an approach for de novo assembly of the carrot mitochondrial genome using next generation sequence data from total genomic DNA.

RESULTS

Sequencing data from a carrot 454 whole genome library were used to develop a de novo assembly of the mitochondrial genome. Development of a new bioinformatic tool allowed visualizing contig connections and elucidation of the de novo assembly. Southern hybridization demonstrated recombination across two large repeats. Genome annotation allowed identification of 44 protein coding genes, three rRNA and 17 tRNA. Identification of the plastid genome sequence allowed organelle genome comparison. Mitochondrial intergenic sequence analysis allowed detection of a fragment of DNA specific to the carrot plastid genome. PCR amplification and sequence analysis across different Apiaceae species revealed consistent conservation of this fragment in the mitochondrial genomes and an insertion in Daucus plastid genomes, giving evidence of a mitochondrial to plastid transfer of DNA. Sequence similarity with a retrotransposon element suggests a possibility that a transposon-like event transferred this sequence into the plastid genome.

CONCLUSIONS

This study confirmed that whole genome sequencing is a practical approach for de novo assembly of higher plant mitochondrial genomes. In addition, a new aspect of intercompartmental genome interaction was reported providing the first evidence for DNA transfer into an angiosperm plastid genome. The approach used here could be used more broadly to sequence and assemble mitochondrial genomes of diverse species. This information will allow us to better understand intercompartmental interactions and cell evolution.

Cloning and characterization of the gene for the iron-sulfur subunit of succinate dehydrogenase from the violet root rot fungus, Helicobasidium mompa

The sdhB gene, encoding the iron-sulfur protein (Ip) subunit of succinate dehydrogenase (Sdh, EC 1.3.99.1), has been cloned from the violet root rot fungus, Helicobasidium mompa, and characterized. The promoter region contains a CCAAT box, TATA-like box, and CT-rich region. The gene is interrupted by eight introns and is predicted to encode a polypeptide of 291 amino acid residues. The putative amino acid sequence of the encoded product of sdhB gene from H. mompa shows high homology to the other known sdhB genes and is 79% identical to the Ip subunit of SdhB of Uromyces fabae. Three cysteine-rich clusters associated with the iron-sulfur centers involved in electron transport were particularly well conserved. One of these clusters contains a critical histidine residue implicated in carboxin sensitivity in the basidiomycetes. Only one copy of the gene was present in the genome of H. mompa, and reverse transcription (RT)-PCR analysis of mRNA expression showed that the sdhB gene was transcribed in potato dextrose broth.

Recent insertion of a 52-kb mitochondrial DNA segment in the wheat lineage

The assembly of a 1.3-Mb size region of the wheat genome has provided the opportunity to study a recent nuclear mitochondrial DNA insertion (NUMT). In the present study, we have studied two bacterial artificial chromosomes (BACs) and characterized a 52-kb NUMT segment from the tetraploid and hexaploid wheat BAC libraries. The conserved orthologous NUMT regions from tetraploid and hexaploid wheat Langdon and Chinese Spring shared identical gene haplotypes even though mutations (insertions, deletions, and substitutions) had occurred. The 52-kb NUMT was present in hexaploid variety Chinese Spring, but absent in variety Hope, by sequence comparison of their corresponding region. Amplifying the NUMT junctions using a set of the wheat materials including diploid, tetraploid, and hexaploid lines showed that none of the diploid wheat carried the region and only some tetraploid and hexaploid wheat were positive for the NUMT. Age estimation of the NUMT displayed the mean ages of Langdon NUMT and Chinese Spring NUMT to be 378,000 and 416,000 years ago, respectively. Reverse transcription PCR and sequencing of the nad7 gene showed 28 C → U RNA editing sites and four partial editing sites, as expected for mitochondrial DNA expression. Specific SNPs discriminated between cDNA from the nucleus and the mitochondria and suggested that the nuclear copy was not expressed. The mitochondrial DNA studied was inserted into the genome quite recently within the wheat lineage and gave rise to the non-coding nuclear nad7 gene. The NUMT segment could be lost and acquired frequently during the wheat evolution.

The mitochondrial genome of the legume Vigna radiata and the analysis of recombination across short mitochondrial repeats

The mitochondrial genomes of seed plants are exceptionally fluid in size, structure, and sequence content, with the accumulation and activity of repetitive sequences underlying much of this variation. We report the first fully sequenced mitochondrial genome of a legume, Vigna radiata (mung bean), and show that despite its unexceptional size (401,262 nt), the genome is unusually depauperate in repetitive DNA and "promiscuous" sequences from the chloroplast and nuclear genomes. Although Vigna lacks the large, recombinationally active repeats typical of most other seed plants, a PCR survey of its modest repertoire of short (38–297 nt) repeats nevertheless revealed evidence for recombination across all of them. A set of novel control assays showed, however, that these results could instead reflect, in part or entirely, artifacts of PCR-mediated recombination. Consequently, we recommend that other methods, especially high-depth genome sequencing, be used instead of PCR to infer patterns of plant mitochondrial recombination. The average-sized but repeat- and feature-poor mitochondrial genome of Vigna makes it ever more difficult to generalize about the factors shaping the size and sequence content of plant mitochondrial genomes.

Functional and composition differences between mitochondrial complex II in Arabidopsis and rice are correlated with the complex genetic history of the enzyme

Complex II plays a central role in mitochondrial metabolism as a component of both the electron transport chain and the tricarboxylic acid cycle. However, the composition and function of the plant enzyme has been elusive and differs from the well-characterised enzymes in mammals and bacteria. Herewith, we demonstrate that mitochondrial Complex II from Arabidopsis and rice differ significantly in several aspects: (1) Stability-Rice complex II in contrast to Arabidopsis is not stable when resolved by native electrophoresis and activity staining. (2) Composition-Arabidopsis complex II contains 8 subunits, only 7 of which have homologs in the rice genome. SDH 1 and 2 subunits display high levels of amino acid identity between two species, while the remainder of the subunits are not well conserved at a sequence level, indicating significant divergence. (3) Gene expression-the pairs of orthologous SDH1 and SDH2 subunits were universally expressed in both Arabidopsis and rice. The very divergent genes for SDH3 and SDH4 were co-expressed in both species, consistent with their functional co-ordination to form the membrane anchor. The plant-specific SDH5, 6 and 7 subunits with unknown functions appeared to be differentially expressed in both species. (4) Biochemical regulation -succinate-dependent O(2) consumption and SDH activity of isolated Arabidopsis mitochondria were substantially stimulated by ATP, but a much more minor effect of ATP was observed for the rice enzyme. The ATP activation of succinate-dependent reduction of DCPIP in frozen-thawed and digitonin-solubilised mitochondrial samples, and with or without the uncoupler CCCP, indicate that the differential ATP effect on SDH is not via the protonmotive force but likely due to an allosteric effect on the plant SDH enzyme itself, in contrast to the enzyme in other organisms.

The complete mitochondrial genome sequence of the hornwort Megaceros aenigmaticus shows a mixed mode of conservative yet dynamic evolution in early land plant mitochondrial genomes

Land plants possess some of the most unusual mitochondrial genomes among eukaryotes. However, in early land plants these genomes resemble those of green and red algae or early eukaryotes. The question of when during land plant evolution the dramatic change in mtDNAs occurred remains unanswered. Here we report the first completely sequenced mitochondrial genome of the hornwort, Megaceros aenigmaticus, a member of the sister group of vascular plants. It is a circular molecule of 184,908 base pairs, with 32 protein genes, 3 rRNA genes, 17 tRNA genes, and 30 group II introns. The genome contains many genes arranged in the same order as in those of a liverwort, a moss, several green and red algae, and Reclinomonas americana, an early-branching eukaryote with the most ancestral form of mtDNA. In particular, the gene order between mtDNAs of the hornwort and Physcomitrella patens (moss) differs by only 8 inversions and translocations. However, the hornwort mtDNA possesses 4 derived features relative to green alga mtDNAs--increased genome size, RNA editing, intron gains, and gene losses--which were all likely acquired during the origin and early evolution of land plants. Overall, this genome and those of other 2 bryophytes show that mitochondrial genomes in early land plants, unlike their seed plant counterparts, exhibit a mixed mode of conservative yet dynamic evolution.

Novel mitochondrial complex II isolated from Trypanosoma cruzi is composed of 12 peptides including a heterodimeric Ip subunit

Mitochondrial respiratory enzymes play a central role in energy production in aerobic organisms. They differentiated from the alpha-proteobacteria-derived ancestors by adding noncatalytic subunits. An exception is Complex II (succinate: ubiquinone reductase), which is composed of four alpha-proteobacteria-derived catalytic subunits (SDH1-SDH4). Complex II often plays a pivotal role in adaptation of parasites in host organisms and would be a potential target for new drugs. We purified Complex II from the parasitic protist Trypanosoma cruzi and obtained the unexpected result that it consists of six hydrophilic (SDH1, SDH2N, SDH2C, and SDH5-SDH7) and six hydrophobic (SDH3, SDH4, and SDH8-SDH11) nucleus-encoded subunits. Orthologous genes for each subunit were identified in Trypanosoma brucei and Leishmania major. Notably, the iron-sulfur subunit was heterodimeric; SDH2N and SDH2C contain the plant-type ferredoxin domain in the N-terminal half and the bacterial ferredoxin domain in the C-terminal half, respectively. Catalytic subunits (SDH1, SDH2N plus SDH2C, SDH3, and SDH4) contain all key residues for binding of dicarboxylates and quinones, but the enzyme showed the lower affinity for both substrates and inhibitors than mammalian enzymes. In addition, the enzyme binds protoheme IX, but SDH3 lacks a ligand histidine. These unusual features are unique in the Trypanosomatida and make their Complex II a target for new chemotherapeutic agents.

Comparative analysis of structural diversity and sequence evolution in plant mitochondrial genes transferred to the nucleus

The transfer of functional mitochondrial genes to the nucleus is an ongoing process during plant evolution that has made a major impact on cytonuclear interactions and mitochondrial genome evolution. Analysis of evolutionarily recent transfers in plants provides insights into the evolutionary dynamics of the process and how transferred genes become functional in the nucleus. Here, we report 42 new transferred genes in various angiosperms, including 9 separate transfers of the succinate dehydrogenase gene sdh3. We performed comparative analyses of gene structures and sequence evolution of 77 genes transferred to the nucleus in various angiosperms, including multiple transfers of 10 genes in different lineages. Many genes contain mitochondrial targeting presequences, and potentially 5' cis-regulatory elements, that were acquired from pre-existing nuclear genes for mitochondrial proteins to create chimeric gene structures. In eight separate cases, the presequence was acquired from either the hsp70 chaperonin gene or the hsp22 chaperonin gene. The most common location of introns is in the presequence, and the least common is in the region transferred from the mitochondrion. Several genes have an intron between the presequence and the core region, or an intron in the 5'UTR (untranslated region) or 3'UTR, suggesting presequence and/or regulatory element acquisition by exon shuffling. Both synonymous and nonsynonymous substitution rates have increased considerably in the transferred genes compared with their mitochondrial counterparts, and the degree of rate acceleration varies by gene, species, and evolutionary timing of transfer. Pairwise and branchwise K(a)/K(s) analysis identified four genes with evidence for positive selection, but positive selection is generally uncommon in transferred genes. This study provides a detailed portrayal of structural and sequence evolution in mitochondrial genes transferred to the nucleus, revealing the frequency of different mechanisms for how presequences and introns are acquired and showing how the sequences of transferred genes evolve after movement between cellular genomes.

Transfer of rice mitochondrial ribosomal protein L6 gene to the nucleus: acquisition of the 5'-untranslated region via a transposable element

Background

The mitochondria of contemporary organisms contain fewer genes than the ancestral bacteria are predicted to have contained. Because most of the mitochondrial proteins are encoded in the nucleus, the genes would have been transferred from the mitochondrion to the nucleus at some stage of evolution and they must have acquired cis-regulatory elements compatible with eukaryotic gene expression. However, most of such processes remain unknown.

Results

The ribosomal protein L6 gene (rpl6) has been lost in presently-known angiosperm mitochondrial genomes. We found that each of the two rice rpl6 genes (OsRpl6-1 and OsRpl6-2) has an intron in an identical position within the 5'-untranslated region (UTR), which suggests a duplication of the rpl6 gene after its transfer to the nucleus. Each of the predicted RPL6 proteins lacks an N-terminal extension as a mitochondrial targeting signal. Transient assays using green fluorescent protein indicated that their mature N-terminal coding regions contain the mitochondrial targeting information. Reverse transcription-PCR analysis showed that OsRpl6-2 expresses considerably fewer transcripts than OsRpl6-1. This might be the result of differences in promoter regions because the 5'-noncoding regions of the two rpl6 genes differ at a point close to the center of the intron. There are several sequences homologous to the region around the 5'-UTR of OsRpl6-1 in the rice genome. These sequences have characteristics similar to those of the transposable elements (TE) belonging to the PIF/Harbinger superfamily.

Conclusion

The above evidences suggest a novel mechanism in which the 5'-UTR of the transferred mitochondrial gene was acquired via a TE. Since the 5'-UTRs and introns within the 5'-UTRs often contain transcriptional and posttranscriptional cis-elements, the transferred rice mitochondrial rpl6 gene may have acquired its cis-element from a TE.

Heat stress enhances the accumulation of polyadenylated mitochondrial transcripts in Arabidopsis thaliana

Background

Polyadenylation of RNA has a decisive influence on RNA stability. Depending on the organisms or subcellular compartment, it either enhances transcript stability or targets RNAs for degradation. In plant mitochondria, polyadenylation promotes RNA degradation, and polyadenylated mitochondrial transcripts are therefore widely considered to be rare and unstable. We followed up a surprising observation that a large number of mitochondrial transcripts are detectable in microarray experiments that used poly(A)-specific RNA probes, and that these transcript levels are significantly enhanced after heat treatment.

Methodology/Principal Findings

As the Columbia genome contains a complete set of mitochondrial genes, we had to identify polymorphisms to differentiate between nuclear and mitochondrial copies of a mitochondrial transcript. We found that the affected transcripts were uncapped transcripts of mitochondrial origin, which were polyadenylated at multiple sites within their 3′region. Heat-induced enhancement of these transcripts was quickly restored during a short recovery period.

Conclusions/Significance

Our results show that polyadenylated transcripts of mitochondrial origin are more stable than previously suggested, and that their steady-state levels can even be significantly enhanced under certain conditions. As many microarrays contain mitochondrial probes, due to the frequent transfer of mitochondrial genes into the genome, these effects need to be considered when interpreting microarray data.

Do recent findings in plant mitochondrial molecular and population genetics have implications for the study of gynodioecy and cytonuclear conflict?

The coexistence of females and hermaphrodites in plant populations, or gynodioecy, is a puzzle recognized by Darwin. Correns identified cytoplasmic inheritance of one component of sex expression, now known as cytoplasmic male sterility (CMS). Lewis established cytonuclear inheritance of gynodioecy as an example of genetic conflict. Although biologists have since developed an understanding of the mechanisms allowing the joint maintenance of CMS and nuclear male fertility restorer genes, puzzles remain concerning the inheritance of sex expression and mechanisms governing the origination of CMS. Much of the theory of gynodioecy rests on the assumption of maternal inheritance of the mitochondrial genome. Here we review recent studies of the genetics of plant mitochondria, and their implications for the evolution and transmission of CMS. New studies of intragenomic recombination provide a plausible origin for the chimeric ORFs that characterize CMS. Moreover, evidence suggests that nonmaternal inheritance of mitochondria may be more common than once believed. These findings may have consequences for the maintenance of cytonuclear polymorphism, mitochondrial recombination, generation of gynomonoecious phenotypes, and interpretation of experimental crosses. Finally we point out that CMS can alter the nature of the cytonuclear conflict that may have originally selected for uniparental inheritance.

Extensive variation in nuclear mitochondrial DNA content between the genomes of Phytophthora sojae and Phytophthora ramorum

Fragments of mitochondrial DNA (mtDNA) transferred to the nuclear genome are called nuclear mitochondrial DNAs (NUMTs). We report here a comparison of NUMT content between genomes from two species of the same genus. Analysis of the genomes of Phytophthora sojae and P. ramorum revealed large differences in the NUMT content of the two genomes: 16.27 x 10(-3) and 2.28 x 10(-3)% of each genome, respectively. Substantial differences also exist between the two species in the sizes of the NUMTs found in each genome, with ranges of 20 to 405 bp for P. sojae and 19 to 137 bp for P. ramorum. Furthermore, in P. sojae, fragments from the mitochondrial genes rns, rnl, coxl, and nad (various subunits) are found most frequently, whereas P. ramorum NUMTs most often originate from the cox3, rpsl4, nad4, and nad5 genes. The large differences in the presumptive mtDNA insertions suggest that the insertions occurred subsequent to the divergence of the two species, and this is supported by sequence comparisons among the NUMTs and the mtDNA sequences of the two species. P. sojae mtDNA sequences inserted in the nuclear genome appear to have been altered as a result of insertions, deletions, inversions, and translocations and provide insights into active mechanisms of sequence divergence in this plant pathogen. No clear examples were found of NUMTs forming functional nuclear genes or of NUMTs inserted into exons or introns of any nuclear gene.

Analysis of Euglena gracilis plastid-targeted proteins reveals different classes of transit sequences

The plastid of Euglena gracilis was acquired secondarily through an endosymbiotic event with a eukaryotic green alga, and as a result, it is surrounded by a third membrane. This membrane complexity raises the question of how the plastid proteins are targeted to and imported into the organelle. To further explore plastid protein targeting in Euglena, we screened a total of 9,461 expressed sequence tag (EST) clusters (derived from 19,013 individual ESTs) for full-length proteins that are plastid localized to characterize their targeting sequences and to infer potential modes of translocation. Of the 117 proteins identified as being potentially plastid localized whose N-terminal targeting sequences could be inferred, 83 were unique and could be classified into two major groups. Class I proteins have tripartite targeting sequences, comprising (in order) an N-terminal signal sequence, a plastid transit peptide domain, and a predicted stop-transfer sequence. Within this class of proteins are the lumen-targeted proteins (class IB), which have an additional hydrophobic domain similar to a signal sequence and required for further targeting across the thylakoid membrane. Class II proteins lack the putative stop-transfer sequence and possess only a signal sequence at the N terminus, followed by what, in amino acid composition, resembles a plastid transit peptide. Unexpectedly, a few unrelated plastid-targeted proteins exhibit highly similar transit sequences, implying either a recent swapping of these domains or a conserved function. This work represents the most comprehensive description to date of transit peptides in Euglena and hints at the complex routes of plastid targeting that must exist in this organism.

Patterns of partial RNA editing in mitochondrial genes of Beta vulgaris

RNA editing is a process that modifies the information in transcripts of almost all angiosperm mitochondrial protein-coding genes. In order to determine the frequency and distribution of mitochondrial RNA editing in Beta vulgaris, cDNAs were sequenced and compared to the published genome sequence. 357 C to U conversions were identified across the 31 known protein genes and pseudogenes in Beta, the fewest so far for a plant mitochondrial genome. Editing patterns in the putative gene orf518 indicate that it is most likely a functional ccmC homolog, indicating that patterns of editing can be a useful determinant of gene functionality. orf518 also contains a triplicated repeat region whose members are nearly identical yet differentially edited, most likely due to differences in the sequence context of the editing sites. In addition, we show that partial editing in Beta is common at silent editing sites but rare at nonsilent editing sites, extending previous observations to a complete plant mitochondrial genome. Finally, the degree of partial editing observed for certain genes was dependent on the choice of primers used, demonstrating that care must be taken when designing primers for use in editing studies.

Pervasive survival of expressed mitochondrial rps14 pseudogenes in grasses and their relatives for 80 million years following three functional transfers to the nucleus

Background

Many mitochondrial genes, especially ribosomal protein genes, have been frequently transferred as functional entities to the nucleus during plant evolution, often by an RNA-mediated process. A notable case of transfer involves the rps14 gene of three grasses (rice, maize, and wheat), which has been relocated to the intron of the nuclear sdh2 gene and which is expressed and targeted to the mitochondrion via alternative splicing and usage of the sdh2 targeting peptide. Although this transfer occurred at least 50 million years ago, i.e., in a common ancestor of these three grasses, it is striking that expressed, nearly intact pseudogenes of rps14 are retained in the mitochondrial genomes of both rice and wheat. To determine how ancient this transfer is, the extent to which mitochondrial rps14 has been retained and is expressed in grasses, and whether other transfers of rps14 have occurred in grasses and their relatives, we investigated the structure, expression, and phylogeny of mitochondrial and nuclear rps14 genes from 32 additional genera of grasses and from 9 other members of the Poales.

Results

Filter hybridization experiments showed that rps14 sequences are present in the mitochondrial genomes of all examined Poales except for members of the grass subfamily Panicoideae (to which maize belongs). However, PCR amplification and sequencing revealed that the mitochondrial rps14 genes of all examined grasses (Poaceae), Cyperaceae, and Joinvilleaceae are pseudogenes, with all those from the Poaceae sharing two 4-NT frameshift deletions and all those from the Cyperaceae sharing a 5-NT insertion (only one member of the Joinvilleaceae was examined). cDNA analysis showed that all mitochondrial pseudogenes examined (from all three families) are transcribed, that most are RNA edited, and that surprisingly many of the edits are reverse (U→C) edits. Putatively nuclear copies of rps14 were isolated from one to several members of each of these three Poales families. Multiple lines of evidence indicate that the nuclear genes are probably the products of three independent transfers.

Conclusion

The rps14 gene has, most likely, been functionally transferred from the mitochondrion to the nucleus at least three times during the evolution of the Poales. The transfers in Cyperaceae and Poaceae are relatively ancient, occurring in the common ancestor of each family, roughly 80 million years ago, whereas the putative Joinvilleaceae transfer may be the most recent case of functional organelle-to-nucleus transfer yet described in any organism. Remarkably, nearly intact and expressed pseudogenes of rps14 have persisted in the mitochondrial genomes of most lineages of Poaceae and Cyperaceae despite the antiquity of the transfers and of the frameshift and RNA editing mutations that mark the mitochondrial genes as pseudogenes. Such long-term, nearly pervasive survival of expressed, apparent pseudogenes is to our knowledge unparalleled in any genome. Such survival probably reflects a combination of factors, including the short length of rps14, its location immediately downstream of rpl5 in most plants, and low rates of nucleotide substitutions and indels in plant mitochondrial DNAs. Their survival also raises the possibility that these rps14 sequences may not actually be pseudogenes despite their appearance as such. Overall, these findings indicate that intracellular gene transfer may occur even more frequently in angiosperms than already recognized and that pseudogenes in plant mitochondrial genomes can be surprisingly resistant to forces that lead to gene loss and inactivation.