Genomic recombination of the mitochondrial atp6 gene in Arabidopsis thaliana at the protein processing site creates two different presequences

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.

Targeted gene disruption of ATP synthases 6-1 and 6-2 in the mitochondrial genome of Arabidopsis thaliana by mitoTALENs

We recently achieved targeted disruptions of cytoplasmic male sterility (CMS)-associated genes in the mitochondrial genomes of rice and rapeseed by using mitochondria-targeted transcription activator-like effector nucleases (mitoTALENs). It was the first report of stable and heritable targeted gene modification of plant mitochondrial genomes. Here, we attempted to use mitoTALENs to disrupt two mitochondrial genes in the model plant Arabidopsis thaliana(Arabidopsis) using three different promoters and two types of TALENs. The targets were the two isoforms of the ATP synthase subunit 6 gene, atp6-1 and atp6-2. Each of these genes was successfully deleted and the mitochondrial genomes were recovered in a homoplasmic state. The nuclear genome also has a copy of atp6-1, and we were able to confirm that it was the mitochondrial gene and not the nuclear pseudogene that was knocked out. Among the three mitoTALEN promoters tried, the RPS5A promoter was the most effective. Conventional mitoTALENs were more effective than single-molecule mito-compactTALENs. Targeted mitochondrial gene deletion was achieved by crossing as well as by floral-dip transformation to introduce the mitoTALEN constructs into the nucleus. The gene disruptions were caused by large (kb-size) deletions. The ends of the remaining sequences were connected to distant loci, mostly by illegitimate homologous recombinations between repeats.

The Mitochondrial Genome Assembly of Fennel (Foeniculum vulgare) Reveals Two Different atp6 Gene Sequences in Cytoplasmic Male Sterile Accessions

Cytoplasmic male sterility (CMS) has always aroused interest among researchers and breeders, being a valuable resource widely exploited not only to breed F1 hybrid varieties but also to investigate genes that control stamen and pollen development. With the aim of identifying candidate genes for CMS in fennel, we adopted an effective strategy relying on the comparison between mitochondrial genomes (mtDNA) of both fertile and sterile genotypes. mtDNA raw reads derived from a CMS genotype were assembled in a single molecule (296,483 bp), while a draft mtDNA assembly (166,124 nucleotides, 94 contigs) was performed using male fertile sample (MF) sequences. From their annotation and alignment, two atp6-like sequences were identified. atp6-, the putative mutant copy with a 300 bp truncation at the 5'-end, was found only in the mtDNA of CMS samples, while the wild type copy (atp6+) was detected only in the MF mtDNA. Further analyses (i.e., reads mapping and Sanger sequencing), revealed an atp6+ copy also in CMS samples, probably in the nuclear DNA. However, qPCRs performed on different tissues proved that, despite its availability, atp6+ is expressed only in MF samples, while apt6- mRNA was always detected in CMS individuals. In the light of these findings, the energy deficiency model could explain the pollen deficiency observed in male sterile flower. atp6- could represent a gene whose mRNA is translated into a not-fully functional protein leading to suboptimal ATP production that guarantees essential cellular processes but not a high energy demand process such as pollen development. Our study provides novel insights into the fennel mtDNA genome and its atp6 genes, and paves the way for further studies aimed at understanding their functional roles in the determination of male sterility.

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.

Genome-Wide Identification and Transcriptional Expression Profiles of the F-box Gene Family in Common Walnut (Juglans regia L.)

The common walnut (or Persian walnut), Juglans regia L., is an economically important temperate tree species valued for both its edible nut and high-quality wood. F-box gene family members are involved in plant development, which includes regulating plant development, reproduction, cellular protein degradation, response to biotic and abiotic stresses, and flowering. However, in common walnut (J. regia), there are no reports about the F-box gene family. Here, we report a genome-wide identification of J. regia F-box genes and analyze their phylogeny, duplication, microRNA, pathway, and transcriptional expression profile. In this study, 74 F-box genes were identified and clustered into three groups based on phylogenetic analysis and eight subfamilies based on special domains in common walnut. These common walnut F-box genes are distributed on 31 different pseudo-chromosomes. The gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and microRNA profiles showed that the F-box gene family might play a critical role in the flowering of common walnut. The expressions were significantly higher in female flowers and male flowers compared with leaf and hull tissues at a transcriptome level. The results revealed that the expressions of the F-box gene in female flowers were positively correlated with male flowers, but there was no correlation between any other tissue combinations in common walnut. Our results provided insight into the general characteristics of the F-box genes in common walnut.

Different lengths, copies and expression levels of the mitochondrial atp6 gene in male sterile and fertile lines of carrot (Daucus carota L.)

The male-sterile carrot is an effective material for carrot breeding. The atp6 gene is involved in carrot fertility. However, the differences in lengths, copies, and expression profiles of the atp6 gene in fertile and male-sterile lines of carrot are unclear. In this study, one copy atp6 gene was found in the mtDNAs of 'Kuroda' (fertility, 954 bp) and 'Wuye-BY' (male sterility, 819 bp) carrot lines, while two copies atp6 genes (Wuye-L and Wuye-D, 954 bp and 819 bp, respectively) were found in the mtDNA of 'Wuye' (fertility). Two putative conserved domains have been detected in the carrot atp6 protein. Evolutionary analysis showed that the atp6 protein sequences of Wuye-L and Kuroda were clustered in the same branch, while Wuye-D and Wuye-BY were clustered in the same branch. The atp6 gene was higher expressed in the flowers of 'Kuroda' and 'Wuye' (Wuye-L), while lower expressed in 'Wuye-BY' and 'Wuye' (Wuye-D).

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 silico analysis of SSRs in mitochondrial genomes of plants

Simple sequence repeats (SSRs) or microsatellites constitute a countable portion of genomes. However, the significance of SSRs in organelle genomes has not been completely understood. The availability of organelle genome sequences allows us to understand the organization of SSRs in their genic and intergenic regions. In the current study we surveyed the patterns of SSRs in mitochondrial genomes of different taxa of plants. A total of 16 mitochondrial genomes, from algae to angiosperms, have been considered to analyze the pattern of simple sequence repeats present in them. Based on study, the mononucleotide repeats of A/T were found to be more prevalent in mitochondrial genomes over other repeat types. The dinucleotides repeats, TA/AT, were the second most numerous, whereas tri-, tetra-, and pentanucleotide repeats were in less number and present in intronic or intergenic portions only. Mononucleotide repeats prevailed in protein-coding exonic portions of all organisms. These results indicates that microsatellite pattern in mitochondrial genomes is different from nuclear genomes and also focuses on organization and diversity at SSR locuses in mitochondrial genomes. This is the novel report of microsatellite polymorphism in plant mitochondrion on whole genome level.

Phylogenetic analysis of mitochondrial substitution rate variation in the angiosperm tribe Sileneae

BACKGROUND

Recent phylogenetic studies have revealed that the mitochondrial genome of the angiosperm Silene noctiflora (Caryophyllaceae) has experienced a massive mutation-driven acceleration in substitution rate, placing it among the fastest evolving eukaryotic genomes ever identified. To date, it appears that other species within Silene have maintained more typical substitution rates, suggesting that the acceleration in S. noctiflora is a recent and isolated evolutionary event. This assessment, however, is based on a very limited sampling of taxa within this diverse genus.

RESULTS

We analyzed the substitution rates in 4 mitochondrial genes (atp1, atp9, cox3 and nad9) across a broad sample of 74 species within Silene and related genera in the tribe Sileneae. We found that S. noctiflora shares its history of elevated mitochondrial substitution rate with the closely related species S. turkestanica. Another section of the genus (Conoimorpha) has experienced an acceleration of comparable magnitude. The phylogenetic data remain ambiguous as to whether the accelerations in these two clades represent independent evolutionary events or a single ancestral change. Rate variation among genes was equally dramatic. Most of the genus exhibited elevated rates for atp9 such that the average tree-wide substitution rate for this gene approached the values for the fastest evolving branches in the other three genes. In addition, some species exhibited major accelerations in atp1 and/or cox3 with no correlated change in other genes. Rates of non-synonymous substitution did not increase proportionally with synonymous rates but instead remained low and relatively invariant.

CONCLUSION

The patterns of phylogenetic divergence within Sileneae suggest enormous variability in plant mitochondrial mutation rates and reveal a complex interaction of gene and species effects. The variation in rates across genomic and phylogenetic scales raises questions about the mechanisms responsible for the evolution of mutation rates in plant mitochondrial genomes.

A mitochondrial gene involved in cytochrome c maturation (ccmC) is expressed as a precursor with a long NH2-terminal extension in sugar beet

Extensive genome rearrangement is one of the major mechanisms of angiosperm mitochondrial evolution. As a by-product, some angiosperm mitochondrial genes exhibit divergent organization, but not all of these genes have been fully characterized. Sugar beet ccmC, which plays an important role in cytochrome c maturation, harbors a unique extended NH(2) terminal region of 277 amino acid residues (N-extension) instead of a conserved translational initiation codon. The 5' termini of two major RNA species were determined by primer extension analysis, which revealed that the larger transcript covered the entire N-extension. Nucleotide sequencing of the cDNA revealed that a total of 31 C-to-U RNA editing events occurred in the N-extension and the ccmC-homologous region (ccmC-core region), resulting in improvement of amino acid sequence conservation. Antiserum was raised against a synthetic peptide corresponding to the ccmC-core region and was used for protein gel blot analysis of sugar beet and radish mitochondrial proteins. The detected 29.5-kDa signal band is shared by sugar beet and radish. Two additional larger signal bands are exclusively detected from sugar beet. The largest signal band is also detected by anti-N-extension antiserum. Our results indicate that sugar beet ccmC is translated as a long precursor with N-extension.

Resolving and identifying protein components of plant mitochondrial respiratory complexes using three dimensions of gel electrophoresis

Analyzing highly hydrophobic proteins is a challenge for identification protocols based on gel separation and mass spectrometry. We combined Blue Native and 2D tricine gel electrophoresis to allow separation and identification of respiratory complex subunits from Arabidopsis mitochondria. We identified many of the highly hydrophobic mitochondrion-encoded subunits (GRAVY scores between +0.6 to +1.4) and also found a number of nucleus-encoded proteins associated with complex I for the first time in plants.

Comparisons among two fertile and three male-sterile mitochondrial genomes of maize

We have sequenced five distinct mitochondrial genomes in maize: two fertile cytotypes (NA and the previously reported NB) and three cytoplasmic-male-sterile cytotypes (CMS-C, CMS-S, and CMS-T). Their genome sizes range from 535,825 bp in CMS-T to 739,719 bp in CMS-C. Large duplications (0.5-120 kb) account for most of the size increases. Plastid DNA accounts for 2.3-4.6% of each mitochondrial genome. The genomes share a minimum set of 51 genes for 33 conserved proteins, three ribosomal RNAs, and 15 transfer RNAs. Numbers of duplicate genes and plastid-derived tRNAs vary among cytotypes. A high level of sequence conservation exists both within and outside of genes (1.65-7.04 substitutions/10 kb in pairwise comparisons). However, sequence losses and gains are common: integrated plastid and plasmid sequences, as well as noncoding "native" mitochondrial sequences, can be lost with no phenotypic consequence. The organization of the different maize mitochondrial genomes varies dramatically; even between the two fertile cytotypes, there are 16 rearrangements. Comparing the finished shotgun sequences of multiple mitochondrial genomes from the same species suggests which genes and open reading frames are potentially functional, including which chimeric ORFs are candidate genes for cytoplasmic male sterility. This method identified the known CMS-associated ORFs in CMS-S and CMS-T, but not in CMS-C.

The organization of mitochondrial atp6 gene region in male fertile and CMS lines of pepper (Capsicum annuum L.)

The mitochondrial atp6 gene in male fertile (N) and CMS (S) pepper has previously been compared and was found to be present in two copies (Kim et al. in J Kor Soc Hort Sci 42:121-127 2001). In the current study, these atp6 copies were amplified by an inverse PCR technique, and the coding region as well as the 5' and 3' flanking regions were sequenced. The atp6 copies in CMS pepper were detected as one intact gene and one pseudogene, truncated at the 3' coding region. When the atp6 genes in pepper were compared to other plant species, pepper, potato, and petunia all possessed a sequence of 12 identical amino acids at the 3' extended region, which was considered a hallmark of the Solanaceae family. Northern blot analysis showed differences in mRNA band patterns between CMS and restorer lines, indicating that atp6 gene is one of the candidates for CMS in pepper.

The 5'-leader sequence of sugar beet mitochondrial atp6 encodes a novel polypeptide that is characteristic of Owen cytoplasmic male sterility

Cytoplasmic male sterility (CMS) is a mitochondrially encoded trait, which is characterized by a failure of plants to produce viable pollen. We have investigated the protein profile of mitochondria from sugar beet plants with normal (fertile) or CMS cytoplasm, and observed that a 35-kDa polypeptide is expressed in Owen CMS plants but not in normal plants. The variant 35-kDa polypeptide was found in CMS mitochondria placed in five different nuclear backgrounds. Interestingly, this polypeptide proved to be antigenically related to a 387-codon ORF (preSatp6) that is fused in-frame with the downstream atp6. The presequence extension of the atp6 ORF is commonly found in higher plants, but whether or not it is normally expressed has hitherto remained unclear. Our study is thus the first to demonstrate that the atp6 presequence is actually translated in mitochondria. We also observed that preSATP6 is a mitochondrial membrane protein that assembles into a homogeneous 200-kDa protein complex. In organello translation experiments in the presence of protease inhibitors showed a reduction in the abundance of mature preSATP6 with time, suggesting that the mature preSATP6 may be derived by proteolytic processing of a translation product of the preSatp6/Satp6 ORF.

The pea mitochondrial atp6: RNA editing and similarity of presequences in the Vicieae tribe

The atp6 gene has been identified as a single-copy sequence in the mitochondrial genome of the pea. An unexpected finding concerns the atp6 5' extension which is known to be poorly conserved at the sequence level, even between closely related plant species. We have shown that the presequences of ATP6 from the pea and other species belonging to the Vicieae tribe of Fabaceae (broad bean, hairy vetch) share a sequence similarity which extends to long 5' untranslated transcript termini. The reason for the observed conservation is not clear but may simply reflect the close phylogenetic relationship of species from the Vicieae tribe. The result of editing analysis indicates the occurrence of fully and partially edited transcripts of atp6 in the pea mitochondria. The majority of the editing sites are targeted to the last transmembrane domain of the pea ATP6, important in proton translocation and interactions with other subunits of ATP synthase.

An Arabidopsis homologue of bacterial RecA that complements an E. coli recA deletion is targeted to plant mitochondria

Homologous recombination results in the exchange and rearrangement of DNA, and thus generates genetic variation in living organisms. RecA is known to function in all bacteria as the central enzyme catalyzing strand transfer and has functional homologues in eukaryotes. Most of our knowledge of homologous recombination in eukaryotes is limited to processes in the nucleus. The mitochondrial genomes of higher plants contain repeated sequences that are known to undergo frequent rearrangements and recombination events. However, very little is known about the proteins involved or the biochemical mechanisms of DNA recombination in plant mitochondria. We provide here the first report of an Arabidopsis thaliana homologue of Escherichia coli RecA that is targeted to mitochondria. The mt recA gene has a putative mitochondrial presequence identified from the A. thaliana genome database. This nuclear gene encodes a predicted product that shows highest sequence homology to chloroplast RecA and RecA proteins from proteobacteria. When fused to the GFP coding sequence, the predicted presequence was able to target the fusion protein to isolated mitochondria but not to chloroplasts. The mitochondrion-specific localization of the mt recA gene product was confirmed by Western analysis using polyclonal antibodies raised against a synthetic peptide from a unique region of the mature mtRecA. The Arabidopsis mt recA gene partially complemented a recA deletion in E. coli, enhancing survival after exposure to DNA-damaging agents. These results suggest a possible role for mt recA in homologous recombination and/or repair in Arabidopsis mitochondria.

New insights into the role of CcmC, CcmD and CcmE in the haem delivery pathway during cytochrome c maturation by a complete mutational analysis of the conserved tryptophan-rich motif of CcmC

Maturation of c-type cytochromes in Escherichia coli is a complex process requiring eight membrane proteins encoded by the ccmABCDEFGH operon. CcmE is a mediator of haem delivery. It binds haem transiently at a conserved histidine residue and releases it for directed transfer to apocytochrome c. CcmC, an integral membrane protein with six transmembrane helices, is necessary and sufficient to incorporate haem covalently into CcmE. CcmC contains a highly conserved tryptophan-rich motif, WGXXWXWD, in its second periplasmic loop. Here, we present the results of a systematic mutational analysis of this motif. Changes of the non-conserved T121 and W122 to A resulted in wild-type CcmC activity. Changes of the single amino acids W119A, G120A, W123A, W125I and D126A or of the spacing within the motif by deleting V124 (DeltaV124) inhibited the covalent haem incorporation into CcmE. Enhanced expression of ccmD suppressed this mutant phenotype by increasing the amounts of CcmC and CcmE polypeptides in the membrane. The DeltaV124 mutant showed the strongest defect of all single mutants. Mutants in which six residues of the tryptophan-rich motif were changed showed no residual CcmC activity. This phenotype was independent of the level of ccmD expression. Our results demonstrate the functional importance of the tryptophan-rich motif for haem transfer to CcmE. We propose that the three membrane proteins CcmC, CcmD and CcmE interact directly with each other, establishing a cytoplasm to periplasm haem delivery pathway for cytochrome c maturation.

The mitochondrial genome of Arabidopsis is composed of both native and immigrant information

Plants contain large mitochondrial genomes, which are several times as complex as those in animals, fungi or algae. However, genome size is not correlated with information content. The mitochondrial genome (mtDNA) of Arabidopsis specifies only 58 genes in 367 kb, whereas the 184 kb mtDNA in the liverwort Marchantia polymorpha codes for 66 genes, and the 58 kb genome in the green alga Prototheca wickerhamii encodes 63 genes. In Arabidopsis' mtDNA, genes for subunits of complex II, for several ribosomal proteins and for 16 tRNAs are missing, some of which have been transferred recently to the nuclear genome. Numerous integrated fragments originate from alien genomes, including 16 sequence stretches of plastid origin, 41 fragments of nuclear (retro)transposons and two fragments of fungal viruses. These immigrant sequences suggest that the large size of plant mitochondrial genomes is caused by secondary expansion as a result of integration and propagation, and is thus a derived trait established during the evolution of land plants.