The application of RNA-seq to the comprehensive analysis of plant mitochondrial transcriptomes

Systematic identification of reference genes for qRT-PCR of Ardisia kteniophylla A. DC under different experimental conditions and for anthocyanin-related genes studies

Ardisia kteniophylla A. DC, widely known as folk medicinal herb and ornamental plant, has been extensively investigated due to its unique leaf color, anti-cancer and other pharmacological activities. The quantitative real-time PCR (qRT-PCR) was an excellent tool for the analysis of gene expression with its high sensitivity and quantitative properties. Normalizing gene expression with stable reference genes was essential for qRT-PCR accuracy. In addition, no studies have yet been performed on the selection, verification and stability of internal reference genes suitable for A. kteniophylla, which has greatly hindered the biomolecular researches of this species. In this study, 29 candidate genes were successfully screened according to stable expression patterns of large-scale RNA seq data that from a variety of tissues and the roots of different growth stages in A. kteniophylla. The candidates were then further determined via qRT-PCR in various experimental samples, including MeJA, ABA, SA, NaCl, CuSO4, AgNO3, MnSO4, CoCl2, drought, low temperature, heat, waterlogging, wounding and oxidative stress. To assess the stability of the candidates, five commonly used strategies were employed: delta-CT, geNorm, BestKeeper, NormFinder, and the comprehensive tool RefFinder. In summary, UBC2 and UBA1 were found to be effective in accurately normalizing target gene expression in A. kteniophella regardless of experimental conditions, while PP2A-2 had the lowest stability. Additionally, to verify the reliability of the recommended reference genes under different colored leaf samples, we examined the expression patterns of six genes associated with anthocyanin synthesis and regulation. Our findings suggested that PAP1 and ANS3 may be involved in leaf color change in A. kteniphella. This study successfully identified the ideal reference gene for qRT-PCR analysis in A. kteniphella, providing a foundation for future research on gene function, particularly in the biosynthesis of anthocyanins.

Evolution of cox2 introns in angiosperm mitochondria and efficient splicing of an elongated cox2i691 intron

In angiosperms, the mitochondrial cox2 gene harbors up to two introns, commonly referred to as cox2i373 and cox2i691. We studied the cox2 from 222 fully-sequenced mitogenomes from 30 angiosperm orders and analyzed the evolution of their introns. Unlike cox2i373, cox2i691 shows a distribution among plants that is shaped by frequent intron loss events driven by localized retroprocessing. In addition, cox2i691 exhibits sporadic elongations, frequently in domain IV of introns. Such elongations are poorly related to repeat content and two of them showed the presence of LINE transposons, suggesting that increasing intron size is very likely due to nuclear intracelular DNA transfer followed by incorporation into the mitochondrial DNA. Surprisingly, we found that cox2i691 is erroneously annotated as absent in 30 mitogenomes deposited in public databases. Although each of the cox2 introns is ∼1.5 kb in length, a cox2i691 of 4.2 kb has been reported in Acacia ligulata (Fabaceae). It is still unclear whether its unusual length is due to a trans-splicing arrangement or the loss of functionality of the interrupted cox2. Through analyzing short-read RNA sequencing of Acacia with a multi-step computational strategy, we found that the Acacia cox2 is functional and its long intron is spliced in cis in a very efficient manner despite its length.

Evolution of mitochondrial RNA editing in extant gymnosperms

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported.

Total and Mitochondrial Transcriptomic and Proteomic Insights into Regulation of Bioenergetic Processes for Shoot Fast-Growth Initiation in Moso Bamboo

As a fast-growing, woody grass plant, Moso bamboo (Phyllostachys edulis) can supply edible shoots, building materials, fibrous raw material, raw materials for crafts and furniture and so on within a relatively short time. Rapid growth of Moso bamboo occurs after the young bamboo shoots are covered with a shell and emerge from the ground. However, the molecular reactions of bioenergetic processes essential for fast growth remain undefined. Herein, total and mitochondrial transcriptomes and proteomes were compared between spring and winter shoots. Numerous key genes and proteins responsible for energy metabolism were significantly upregulated in spring shoots, including those involved in starch and sucrose catabolism, glycolysis, the pentose phosphate pathway, the tricarboxylic acid cycle and oxidative phosphorylation. Accordingly, significant decreases in starch and soluble sugar, higher ATP content and higher rates of respiration and glycolysis were identified in spring shoots. Further, the upregulated genes and proteins related to mitochondrial fission significantly increased the number of mitochondria, indirectly promoting intracellular energy metabolism. Moreover, enhanced alternate-oxidase and uncoupled-protein pathways in winter shoots showed that an efficient energy-dissipating system was important for winter shoots to adapt to the low-temperature environment. Heterologous expression of PeAOX1b in Arabidopsis significantly affected seedling growth and enhanced cold-stress tolerance. Overall, this study highlights the power of comparing total and mitochondrial omics and integrating physiochemical data to understand how bamboo initiates fast growth through modulating bioenergetic processes.

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

In land plant mitochondria, C-to-U RNA editing converts cytidines into uridines at highly specific RNA positions called editing sites. This editing step is essential for the correct functioning of mitochondrial proteins. When using sequence homology information, edited positions can be computationally predicted with high precision. However, predictions based on the sequence contexts of such edited positions often result in lower precision, which is limiting further advances on novel genetic engineering techniques for RNA regulation. Here, a deep convolutional neural network called Deepred-Mt is proposed. It predicts C-to-U editing events based on the 40 nucleotides flanking a given cytidine. Unlike existing methods, Deepred-Mt was optimized by using editing extent information, novel strategies of data augmentation, and a large-scale training dataset, constructed with deep RNA sequencing data of 21 plant mitochondrial genomes. In comparison to predictive methods based on sequence homology, Deepred-Mt attains significantly better predictive performance, in terms of average precision as well as F1 score. In addition, our approach is able to recognize well-known sequence motifs linked to RNA editing, and shows that the local RNA structure surrounding editing sites may be a relevant factor regulating their editing. These results demonstrate that Deepred-Mt is an effective tool for predicting C-to-U RNA editing in plant mitochondria. Source code, datasets, and detailed use cases are freely available at https://github.com/aedera/deepredmt.

Horizontal gene transfers dominate the functional mitochondrial gene space of a holoparasitic plant

Although horizontal gene transfer (HGT) is common in angiosperm mitochondrial DNAs (mtDNAs), few cases of functional foreign genes have been identified. The one outstanding candidate for large-scale functional HGT is the holoparasite Lophophytum mirabile, whose mtDNA has lost most native genes but contains intact foreign homologs acquired from legume host plants. To investigate the extent to which this situation results from functional replacement of native by foreign genes, functional mitochondrial gene transfer to the nucleus, and/or loss of mitochondrial biochemical function in the context of extreme parasitism, we examined the Lophophytum mitochondrial and nuclear transcriptomes by deep paired-end RNA sequencing. Most foreign mitochondrial genes in Lophophytum are highly transcribed, accurately spliced, and efficiently RNA edited. By contrast, we found no evidence for functional gene transfer to the nucleus or loss of mitochondrial functions in Lophophytum. Many functional replacements occurred via the physical replacement of native genes by foreign genes. Some of these events probably occurred as the final act of HGT itself. Lophophytum mtDNA has experienced an unprecedented level of functional replacement of native genes by foreign copies. This raises important questions concerning population-genetic and molecular regimes that underlie such a high level of foreign gene takeover.

Differentially Expressed Genes Shared by Two Distinct Cytoplasmic Male Sterility (CMS) Types of Silene vulgaris Suggest the Importance of Oxidative Stress in Pollen Abortion

Cytoplasmic male sterility (CMS), encoded by the interacting mitochondrial and nuclear genes, causes pollen abortion or non-viability. CMS is widely used in agriculture and extensively studied in crops. Much less is known about CMS in wild species. We performed a comparative transcriptomic analysis of male sterile and fertile individuals of Silene vulgaris, a model plant for the study of gynodioecy, to reveal the genes responsible for pollen abortion in this species. We used RNA-seq datasets previously employed for the analysis of mitochondrial and plastid transcriptomes of female and hermaphrodite flower buds, making it possible to compare the transcriptomes derived from three genomes in the same RNA specimen. We assembled de novo transcriptomes for two haplotypes of S. vulgaris and identified differentially expressed genes between the females and hermaphrodites, associated with stress response or pollen development. The gene for alternative oxidase was downregulated in females. The genetic pathways controlling CMS in S. vulgaris are similar to those in crops. The high number of the differentially expressed nuclear genes contrasts with the uniformity of organellar transcriptomes across genders, which suggests these pathways are evolutionarily conserved and that selective mechanisms may shield organellar transcription against changes in the cytoplasmic transcriptome.

Shared up-regulation and contrasting down-regulation of gene expression distinguish desiccation-tolerant from intolerant green algae

Among green plants, desiccation tolerance is common in seeds and spores but rare in leaves and other vegetative green tissues. Over the last two decades, genes have been identified whose expression is induced by desiccation in diverse, desiccation-tolerant (DT) taxa, including, e.g., late embryogenesis abundant proteins (LEA) and reactive oxygen species scavengers. This up-regulation is observed in DT resurrection plants, mosses, and green algae most closely related to these Embryophytes. Here we test whether this same suite of protective genes is up-regulated during desiccation in even more distantly related DT green algae, and, importantly, whether that up-regulation is unique to DT algae or also occurs in a desiccation-intolerant relative. We used three closely related aquatic and desert-derived green microalgae in the family Scenedesmaceae and capitalized on extraordinary desiccation tolerance in two of the species, contrasting with desiccation intolerance in the third. We found that during desiccation, all three species increased expression of common protective genes. The feature distinguishing gene expression in DT algae, however, was extensive down-regulation of gene expression associated with diverse metabolic processes during the desiccation time course, suggesting a switch from active growth to energy-saving metabolism. This widespread downshift did not occur in the desiccation-intolerant taxon. These results show that desiccation-induced up-regulation of expression of protective genes may be necessary but is not sufficient to confer desiccation tolerance. The data also suggest that desiccation tolerance may require induced protective mechanisms operating in concert with massive down-regulation of gene expression controlling numerous other aspects of metabolism.

Identification and Selection of Reference Genes for Quantitative Transcript Analysis in Corydalis yanhusuo

Corydalis yanhusuo is a medicinal plant frequently used in traditional Chinese medicine, which has effective medical effects in many aspects. Real-time polymerase chain reaction (RT-PCR) has been one of the most widely used methods in biosynthesis research due to its high sensitivity and quantitative properties in gene expression analysis. To obtain accurate normalization, reference genes are often selected in advance; however, no reference genes are available in C. yanhusuo. Herein, 12 reference gene candidates, named cyclophilin 2 (CYP2), elongation factor 1-α (EF1-α), protein phosphatase 2 (PP2A), SAND protein family (SAND), polypyrimidine tract-binding protein (PTBP), TIP41-like protein (TIP41), lyceraldehyde-3-phosphate hydrogenase (GAPDH), ubiquitin-conjugating enzyme 9 (UBC9), cyclophilin 1 (CYP1), tubulin beta (TUBA), thioredoxin (YLS8), and polyubiquitin 10 (UBQ10), were selected for stability analysis. After being treated with hormone, UV, salt, metal, oxidative, drought, cold (4 °C), and hot stresses (40 °C), the qRT-PCR data of the selected genes was analyzed with NormFinder, geNorm, and BestKeeper. The result indicated that GAPDH, SNAD, and PP2A were the top three most stable reference genes under most treatments. This study selected and validated reliable reference genes in C. yanhusuo under various environmental conditions, which can provide great help for future research on gene expression normalization in C. yanhusuo.

Selection and validation of appropriate reference genes for real-time quantitative PCR analysis in Momordica charantia

Real time quantitative reverse transcription PCR (RT-qPCR) has been attracting more attention for its high sensitivity in gene expression analysis. Given the widely use of RT-qPCR in normalization, it is playing a pivotal role for seeking suitable reference genes in different species. In current work, 12 candidate reference genes including Actin 2 (ACT2), Cyclophilin 2 (CYP2), Glyceraldehyde-3-phosphate dehydrogenase C2 (GAPC2), Elongation factor 1-α (EF1-α), Nuclear cap binding protein 20 (NCBP20), Serine/threonine-protein phosphatase PP2A (PP2A), Polypyrimidine tract-binding protein 1 (PTBP1), SAND family protein (SNAD), TIP41-like protein (TIP41), Tubulin beta-6 (TUB6), Ubiquitin-conjugating enzyme 9 (UBC9) and Glyceraldehyde-3-phosphatedehydrogenase (GAPDH) were screened from the transcriptome datasets of M. charantia. Afterwards, GeNorm, NormFinder and BestKeeper algorithms were applied to assess the expression stability of these 12 genes under different abiotic stresses including drought, cold, high-salt, hormone, UV, oxidative and metal stress. The results indicated that 12 selected genes exhibited various stability across the samples under different external stress conditions, but TIP41, PTBP1 and PP2A presented high stability among all the reference genes. To validate the suitability of the identified reference genes, the results of hormone subset were compared with RNA sequencing (RNA-seq) data, and the relative abundance of Ascorbate peroxidase 1（APX1）was used to confirm the reliability of the results. This work assesses the stability of reference genes in M. charantia under different abiotic stress conditions, which will be beneficent for accurate normalization of target genes in M. charantia.

Arms races with mitochondrial genome soft sweeps in a gynodioecious plant, Plantago lanceolata

Biological situations involving conflict can create arms race situations with repeated fixations of different functional variants, producing selective sweeps and lowering neutral diversity in genome regions linked to the functional locus. However, they can sometimes lead to balancing selection, potentially creating long coalescent times for sites with functionally different variants, and, if recombination occurs rarely, for extended haplotypes carrying such variants. We tested between these possibilities in a gynodioecious plant, Plantago lanceolata, in which cytoplasmic male-sterility factors conflict with nuclear restorers of male fertility. We find low mitochondrial diversity, which does not support very long-term coexistence of highly diverged mitochondrial haplotypes. Interestingly, however, we found a derived haplotype that is associated with male fertility in a restricted geographic region, and that has fixed differences from the ancestral sequence in several genes, suggesting that it did not arise very recently. Taken together, the results suggest arms race events that involved "soft" selective sweeps involving a moderately old-established haplotype, consistent with the frequency fluctuations predicted by theoretical models of gynodioecy.

Mitochondrial genome and transcriptome analysis of five alloplasmic male-sterile lines in Brassica juncea

Background

Alloplasmic lines, in which the nuclear genome is combined with wild cytoplasm, are often characterized by cytoplasmic male sterility (CMS), regardless of whether it was derived from sexual or somatic hybridization with wild relatives. In this study, we sequenced and analyzed the mitochondrial genomes of five such alloplasmic lines in Brassica juncea.

Results

The assembled and annotated mitochondrial genomes of the five alloplasmic lines were found to have virtually identical gene contents. They preserved most of the ancestral mitochondrial segments, and the same candidate male sterility gene (orf108) was found harbored in mitotype-specific sequences. We also detected promiscuous sequences of chloroplast origin that were conserved among plants of the Brassicaceae, and found the RNA editing profiles to vary across the five mitochondrial genomes.

Conclusions

On the basis of our characterization of the genetic nature of five alloplasmic mitochondrial genomes, we speculated that the putative candidate male sterility gene orf108 may not be responsible for the CMS observed in Brassica oxyrrhina and Diplotaxis catholica. Furthermore, we propose the potential coincidence of CMS in alloplasmic lines. Our findings lay the foundation for further elucidation of male sterility gene.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5721-2) contains supplementary material, which is available to authorized users.

High Level of Conservation of Mitochondrial RNA Editing Sites Among Four Populus Species

RNA editing occurs in the endosymbiont organelles of higher plants as C-to-U conversions of defined nucleotides. The availability of large quantities of RNA sequencing data makes it possible to identify RNA editing sites and to quantify their editing extent. We have investigated RNA editing in 34 protein-coding mitochondrial transcripts of four Populus species, a genus noteworthy for its remarkably small number of RNA editing sites compared to other angiosperms. 27 of these transcripts were subject to RNA editing in at least one species. In total, 355 RNA editing sites were identified with high confidence, their editing extents ranging from 10 to 100%. The most heavily edited transcripts were ccmB with the highest density of RNA editing sites (53.7 sites / kb) and ccmFn with the highest number of sites (39 sites). Most of the editing events are at position 1 or 2 of the codons, usually altering the encoded amino acid, and are highly conserved among the species, also with regard to their editing extent. However, one SNP was found in the newly sequenced and annotated mitochondrial genome of P. alba resulting in the loss of an RNA editing site compared to P. tremula and P. davidiana This SNP causes a C-to-T transition and an amino acid exchange from Ser to Phe, highlighting the widely discussed role of RNA editing in compensating mutations.

Toxicity of Melaleuca alternifolia essential oil to the mitochondrion and NAD+/NADH dehydrogenase in Tribolium confusum

BACKGROUND

In our previous study, Melaleuca alternifolia essential oil (EO) was considered to have an insecticidal effect by acting on the mitochondrial respiratory chain in insects. However, the mode of action is not fully understood.

METHODS

In this study, we investigated the insecticidal efficacy of the M. alternifolia EO against another major stored-product pest, Tribolium confusum Jacquelin du Val. Rarefaction and vacuolization of the mitochondrial matrix were evident in oil-fumigated T. confusum adults.

RESULTS

Alterations to the mitochondria confirmed the insecticidal effect of the M. alternifolia EO. Furthermore, comparative transcriptome analysis of T. confusum using RNA-seq indicated that most of the differentially expressed genes were involved in insecticide detoxification and mitochondrial function. The biochemical analysis showed that the intracellular NAD+/NADH ratio is involved in the differential effect of the M. alternifolia EO.

DISCUSSION

These results led us to conclude that NAD+/NADH dehydrogenase may be the prime target site for the M. alternifolia EO in insects, leading to blocking of the mitochondrial respiratory chain.

The S. pombe mitochondrial transcriptome

Mitochondrial gene expression is largely controlled through post-transcriptional processes including mitochondrial RNA (mt-RNA) processing, modification, decay, and quality control. Defective mitochondrial gene expression results in mitochondrial oxidative phosphorylation (OXPHOS) deficiency and has been implicated in human disease. To fully understand mitochondrial transcription and RNA processing, we performed RNA-seq analyses of mt-RNAs from the fission yeast Schizosaccharomyces pombe RNA-seq analyses show that the abundance of mt-RNAs vary greatly. Analysis of data also reveals mt-RNA processing sites including an unusual RNA cleavage event by mitochondrial tRNA (mt-tRNA) 5'-end processing enzyme RNase P. Additionally, this analysis reveals previously unknown mitochondrial transcripts including the rnpB-derived fragment, mitochondrial small RNAs (mitosRNAs) such as mt-tRNA-derived fragments (mt-tRFs) and mt-tRNA halves, and mt-tRNAs marked with 3'-CCACCA/CCACC in S. pombe Finally, RNA-seq reveals that inactivation of trz2 encoding S. pombe mitochondrial tRNA 3'-end processing enzyme globally impairs mt-tRNA 3'-end processing, inhibits mt-mRNA 5'-end processing, and causes accumulation of unprocessed transcripts, demonstrating the feasibility of using RNA-seq to examine the protein known or predicted to be involved in mt-RNA processing in S. pombe Our work uncovers the complexity of a fungal mitochondrial transcriptome and provides a framework for future studies of mitochondrial gene expression using S. pombe as a model system.

Plant organelle RNA editing and its specificity factors: enhancements of analyses and new database features in PREPACT 3.0

Background

Gene expression in plant chloroplasts and mitochondria is affected by RNA editing. Numerous C-to-U conversions, accompanied by reverse U-to-C exchanges in some plant clades, alter the genetic information encoded in the organelle genomes. Predicting and analyzing RNA editing, which ranges from only few sites in some species to thousands in other taxa, is bioinformatically demanding.

Results

Here, we present major enhancements and extensions of PREPACT, a WWW-based service for analysing, predicting and cataloguing plant-type RNA editing. New features in PREPACT’s core include direct GenBank accession query input and options to restrict searches to candidate U-to-C editing or to sites where editing has been documented previously in the references. The reference database has been extended by 20 new organelle editomes. PREPACT 3.0 features new modules “EdiFacts” and “TargetScan”. EdiFacts integrates information on pentatricopeptide repeat (PPR) proteins characterized as site-specific RNA editing factors. PREPACT’s editome references connect into EdiFacts, linking editing events to specific co-factors where known. TargetScan allows position-weighted querying for sequence motifs in the organelle references, optionally restricted to coding regions or sequences around editing sites, or in queries uploaded by the user. TargetScan is mainly intended to evaluate and further refine the proposed PPR-RNA recognition code but may be handy for other tasks as well. We present an analysis for the immediate sequence environment of more than 15,000 documented editing sites finding strong and different bias in the editome data sets.

Conclusions

We exemplarily present the novel features of PREPACT 3.0 aimed to enhance the analyses of plant-type RNA editing, including its new modules EdiFacts integrating information on characterized editing factors and TargetScan aimed to analyse RNA editing site recognition specificities.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2244-9) contains supplementary material, which is available to authorized users.

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

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

The Role of Non-Coding RNAs in Cytoplasmic Male Sterility in Flowering Plants

The interactions between mitochondria and nucleus substantially influence plant development, stress response and morphological features. The prominent example of a mitochondrial-nuclear interaction is cytoplasmic male sterility (CMS), when plants produce aborted anthers or inviable pollen. The genes responsible for CMS are located in mitochondrial genome, but their expression is controlled by nuclear genes, called fertility restorers. Recent explosion of high-throughput sequencing methods enabled to study transcriptomic alterations in the level of non-coding RNAs under CMS biogenesis. We summarize current knowledge of the role of nucleus encoded regulatory non-coding RNAs (long non-coding RNA, microRNA as well as small interfering RNA) in CMS. We also focus on the emerging data of non-coding RNAs encoded by mitochondrial genome and their possible involvement in mitochondrial-nuclear interactions and CMS development.

De novo Transcriptome Sequencing of MeJA-Induced Taraxacum koksaghyz Rodin to Identify Genes Related to Rubber Formation

Increase in the consumption of natural rubber (NR) has necessitated the identification of alternative sources of NR. The quality of NR produced by Taraxacum koksaghyz Rodin (TKS) is comparable to that from Hevea brasiliensis (H.brasiliensis), and therefore, TKS is being considered as an alternative source of NR. Here, we sequenced the TKS root transcriptome after wild TKS seedlings were treated with methyl jasmonate (MeJA) for 0, 6, and 24 h. The clean reads generated for each experimental line were assembled into 127,833 unigenes. The Kyoto encyclopedia of genes and genomes pathway prediction suggested that methyl jasmonate regulated secondary metabolism in TKS. Differential expression analysis showed that the expression of HMGCR, FPPS, IDI, GGPPS, and REF/SRPP increased with methyl jasmonate treatment. Interestingly, differential expression analysis of the jasmonate (JA)-related transcription factors (TFs), indicated that certain genes encoding these transcription factors (namely, bHLH, MYB, AP2/EREBP, and WRKY) showed the same expression pattern in the lines treated for 6 h and 24 h. Moreover, HMGCR was up-regulated in the transgenic seedlings overexpressing DREB. We predicted that methyl jasmonate regulated secondary metabolism and affected rubber biosynthesis via the interaction of the JA-related TFs with genes associated with rubber biosynthesis in TKS.

The mitochondrial genome of the terrestrial carnivorous plant Utricularia reniformis (Lentibulariaceae): Structure, comparative analysis and evolutionary landmarks

The carnivorous plants of the family Lentibulariaceae have attained recent attention not only because of their interesting lifestyle, but also because of their dynamic nuclear genome size. Lentibulariaceae genomes span an order of magnitude and include species with the smallest genomes in angiosperms, making them a powerful system to study the mechanisms of genome expansion and contraction. However, little is known about mitochondrial DNA (mtDNA) sequences of this family, and the evolutionary forces that shape this organellar genome. Here we report the sequencing and assembly of the complete mtDNA from the endemic terrestrial Brazilian species Utricularia reniformis. The 857,234bp master circle mitochondrial genome encodes 70 transcriptionaly active genes (42 protein-coding, 25 tRNAs and 3 rRNAs), covering up to 7% of the mtDNA. A ltrA-like protein related to splicing and mobility and a LAGLIDADG homing endonuclease have been identified in intronic regions, suggesting particular mechanisms of genome maintenance. RNA-seq analysis identified properties with putative diverse and important roles in genome regulation and evolution: 1) 672kbp (78%) of the mtDNA is covered by full-length reads; 2) most of the 243kbp intergenic regions exhibit transcripts; and 3) at least 69 novel RNA editing sites in the protein-coding genes. Additional genomic features are hypothetical ORFs (48%), chloroplast insertions, including truncated plastid genes that have been lost from the chloroplast DNA (5%), repeats (5%), relics of transposable elements mostly related to LTR retrotransposons (5%), and truncated mitovirus sequences (0.4%). Phylogenetic analysis based on 32 different Lamiales mitochondrial genomes corroborate that Lentibulariaceae is a monophyletic group. In summary, the U. reniformis mtDNA represents the eighth largest plant mtDNA described to date, shedding light on the genomic trends and evolutionary characteristics and phylogenetic history of the family Lentibulariaceae.

Non-coding RNA may be associated with cytoplasmic male sterility in Silene vulgaris

Cytoplasmic male sterility (CMS) is a widespread phenomenon in flowering plants caused by mitochondrial (mt) genes. CMS genes typically encode novel proteins that interfere with mt functions and can be silenced by nuclear fertility-restorer genes. Although the molecular basis of CMS is well established in a number of crop systems, our understanding of it in natural populations is far more limited. To identify CMS genes in a gynodioecious plant, Silene vulgaris, we constructed mt transcriptomes and compared transcript levels and RNA editing patterns in floral bud tissue from female and hermaphrodite full siblings. The transcriptomes from female and hermaphrodite individuals were very similar overall with respect to variation in levels of transcript abundance across the genome, the extent of RNA editing, and the order in which RNA editing and intron splicing events occurred. We found only a single genomic region that was highly overexpressed and differentially edited in females relative to hermaphrodites. This region is not located near any other transcribed elements and lacks an open-reading frame (ORF) of even moderate size. To our knowledge, this transcript would represent the first non-coding mt RNA associated with CMS in plants and is, therefore, an important target for future functional validation studies.

Stable Internal Reference Genes for Normalizing Real-Time Quantitative PCR in Baphicacanthus cusia under Hormonal Stimuli and UV Irradiation, and in Different Plant Organs

Baphicacanthus cusia (Nees) Bremek, the plant source for many kinds of drugs in traditional Chinese medicine, is widely distributed in South China, especially in Fujian. Recent studies about B. cusia mainly focus on its chemical composition and pharmacological effects, but further analysis of the plant's gene functions and expression is required to better understand the synthesis of its effective compounds. Real-time quantitative polymerase chain reaction (RT-qPCR) is a powerful method for gene expression analysis. It is necessary to select a suitable reference gene for expression normalization to ensure the accuracy of RT-qPCR results. Ten candidate reference genes were selected from the transcriptome datasets of B. cusia in this study, and the expression stability was assessed across 60 samples representing different tissues and organs under various conditions, including ultraviolet (UV) irradiation, hormonal stimuli (jasmonic acid methyl ester and abscisic acid), and in different plant organs. By employing different algorithms, such as geNorm, NormFinder, and BestKeeper, which are complementary approaches based on different statistical procedures, 18S rRNA was found to be the most stable gene under UV irradiation and hormonal stimuli, whereas ubiquitin-conjugating enzyme E2 was the best suitable gene for different plant organs. This novel study aimed to screen for suitable reference genes and corresponding primer pairs specifically designed for gene expression studies in B. cusia, in particular for RT-qPCR analyses.

Identification by the DArTseq method of the genetic origin of the Coffea canephora cultivated in Vietnam and Mexico

BACKGROUND

The coffee species Coffea canephora is commercially identified as "Conilon" when produced in Brazil, or "Robusta" when produced elsewhere in the world. It represents approximately 40 % of coffee production worldwide. While the genetic diversity of wild C. canephora has been well studied in the past, only few studies have addressed the genetic diversity of currently cultivated varieties around the globe. Vietnam is the largest Robusta producer in the world, while Mexico is the only Latin American country, besides Brazil, that has a significant Robusta production. Knowledge of the genetic origin of Robusta cultivated varieties in countries as important as Vietnam and Mexico is therefore of high interest.

RESULTS

Through the use of Sequencing-based diversity array technology-DArTseq method-on a collection of C. canephora composed of known accessions and accessions cultivated in Vietnam and Mexico, 4,021 polymorphic SNPs were identified. We used a multivariate analysis using SNP data from reference accessions in order to confirm and further fine-tune the genetic diversity of C. canephora. Also, by interpolating the data obtained for the varieties from Vietnam and Mexico, we determined that they are closely related to each other, and identified that their genetic origin is the Robusta Congo - Uganda group.

CONCLUSIONS

The genetic characterization based on SNP markers of the varieties grown throughout the world, increased our knowledge on the genetic diversity of C. canephora, and contributed to the understanding of the genetic background of varieties from very important coffee producers. Given the common genetic origin of the Robusta varieties cultivated in Vietnam, Mexico and Uganda, and the similar characteristics of climatic areas and relatively high altitude where they are grown, we can state that the Vietnamese and the Mexican Robusta have the same genetic potential to produce good cup quality.

ChloroSeq, an Optimized Chloroplast RNA-Seq Bioinformatic Pipeline, Reveals Remodeling of the Organellar Transcriptome Under Heat Stress

Although RNA-Seq has revolutionized transcript analysis, organellar transcriptomes are rarely assessed even when present in published datasets. Here, we describe the development and application of a rapid and convenient method, ChloroSeq, to delineate qualitative and quantitative features of chloroplast RNA metabolism from strand-specific RNA-Seq datasets, including processing, editing, splicing, and relative transcript abundance. The use of a single experiment to analyze systematically chloroplast transcript maturation and abundance is of particular interest due to frequent pleiotropic effects observed in mutants that affect chloroplast gene expression and/or photosynthesis. To illustrate its utility, ChloroSeq was applied to published RNA-Seq datasets derived from Arabidopsis thaliana grown under control and abiotic stress conditions, where the organellar transcriptome had not been examined. The most appreciable effects were found for heat stress, which induces a global reduction in splicing and editing efficiency, and leads to increased abundance of chloroplast transcripts, including genic, intergenic, and antisense transcripts. Moreover, by concomitantly analyzing nuclear transcripts that encode chloroplast gene expression regulators from the same libraries, we demonstrate the possibility of achieving a holistic understanding of the nucleus-organelle system. ChloroSeq thus represents a unique method for streamlining RNA-Seq data interpretation of the chloroplast transcriptome and its regulators.

Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae

MYB transcription factors play an important role in regulating key plant developmental processes involving defense, cell shape, pigmentation, and root formation. Within this gene family, sequences containing an R2R3 MYB domain are the most abundant type and exhibit a wide diversity of functions. In this study, we identify 559 R2R3 MYB genes using whole genome data from four species of Solanaceae and reconstruct their evolutionary relationships. We compare the Solanaceae R2R3 MYBs to the well-characterized Arabidopsis thaliana sequences to estimate functional diversity and to identify gains and losses of MYB clades in the Solanaceae. We identify numerous R2R3 MYBs that do not appear closely related to Arabidopsis MYBs, and thus may represent clades of genes that have been lost along the Arabidopsis lineage or gained after the divergence of Rosid and Asterid lineages. Despite differences in the distribution of R2R3 MYBs across functional subgroups and species, the overall size of the R2R3 subfamily has changed relatively little over the roughly 50 million-year history of Solanaceae. We added our information regarding R2R3 MYBs in Solanaceae to other data and performed a meta-analysis to trace the evolution of subfamily size across land plants. The results reveal many shifts in the number of R2R3 genes, including a 54 % increase along the angiosperm stem lineage. The variation in R2R3 subfamily size across land plants is weakly positively correlated with genome size and strongly positively correlated with total number of genes. The retention of such a large number of R2R3 copies over long evolutionary time periods suggests that they have acquired new functions and been maintained by selection. Discovering the nature of this functional diversity will require integrating forward and reverse genetic approaches on an -omics scale.

Selection and Validation of Appropriate Reference Genes for Quantitative Real-Time PCR Analysis of Gene Expression in Lycoris aurea

Lycoris aurea (L' Hér.) Herb, a perennial grass species, produces a unique variety of pharmacologically active Amaryllidaceae alkaloids. However, the key enzymes and their expression pattern involved in the biosynthesis of Amaryllidaceae alkaloids (especially for galanthamine) are far from being fully understood. Quantitative real-time polymerase chain reaction (qRT-PCR), a commonly used method for quantifying gene expression, requires stable reference genes to normalize its data. In this study, to choose the appropriate reference genes under different experimental conditions, 14 genes including YLS8 (mitosis protein YLS8), CYP2 (Cyclophilin 2), CYP 1 (Cyclophilin 1), TIP41 (TIP41-like protein), EXP2 (Expressed protein 2), PTBP1 (Polypyrimidine tract-binding protein 1), EXP1 (Expressed protein 1), PP2A (Serine/threonine-protein phosphatase 2A), β-TUB (β-tubulin), α-TUB (α-tubulin), EF1-α (Elongation factor 1-α), UBC (Ubiquitin-conjugating enzyme), ACT (Actin) and GAPDH (Glyceraldehyde 3-phosphate dehydrogenase) were selected from the transcriptome datasets of L. aurea. And then, expressions of these genes were assessed by qRT-PCR in various tissues and the roots under different treatments. The expression stability of the 14 candidates was analyzed by three commonly used software programs (geNorm, NormFinder, and BestKeeper), and their results were further integrated into a comprehensive ranking based on the geometric mean. The results show the relatively stable genes for each subset as follows: (1) EXP1 and TIP41 for all samples; (2) UBC and EXP1 for NaCl stress; (3) PTBP1 and EXP1 for heat stress, polyethylene glycol (PEG) stress and ABA treatment; (4) UBC and CYP2 for cold stress; (5) PTBP1 and PP2A for sodium nitroprusside (SNP) treatment; (6) CYP1 and TIP41 for methyl jasmonate (MeJA) treatment; and (7) EXP1 and TIP41 for various tissues. The reliability of these results was further enhanced through comparison between part qRT-PCR result and RNA sequencing (RNA-seq) data. In summary, our results identified appropriate reference genes for qRT-PCR in L. aurea, and will facilitate gene expression studies under these conditions.

High transcript abundance, RNA editing, and small RNAs in intergenic regions within the massive mitochondrial genome of the angiosperm Silene noctiflora

BACKGROUND

Species within the angiosperm genus Silene contain the largest mitochondrial genomes ever identified. The enormity of these genomes (up to 11 Mb in size) appears to be the result of increased non-coding DNA, which represents >99 % of the genome content. These genomes are also fragmented into dozens of circular-mapping chromosomes, some of which contain no identifiable genes, raising questions about if and how these 'empty' chromosomes are maintained by selection. To assess the possibility that they contain novel and unannotated functional elements, we have performed RNA-seq to analyze the mitochondrial transcriptome of Silene noctiflora.

RESULTS

We identified regions of high transcript abundance in almost every chromosome in the mitochondrial genome including those that lack any annotated genes. In some cases, these transcribed regions exhibited higher expression levels than some core mitochondrial protein-coding genes. We also identified RNA editing sites throughout the genome, including 97 sites that were outside of protein-coding gene sequences and found in pseudogenes, introns, UTRs, and transcribed intergenic regions. Unlike in protein-coding sequences, however, most of these RNA editing sites were only edited at intermediate frequencies. Finally, analysis of mitochondrial small RNAs indicated that most were likely degradation products from longer transcripts, but we did identify candidates for functional small RNAs that mapped to intergenic regions and were not associated with longer RNA transcripts.

CONCLUSIONS

Our findings demonstrate transcriptional activity in many localized regions within the extensive intergenic sequence content in the S. noctiflora mitochondrial genome, supporting the possibility that the genome contains previously unidentified functional elements. However, transcription by itself is not proof of functional importance, and we discuss evidence that some of the observed transcription and post-transcriptional modifications are non-adaptive. Therefore, further investigations are required to determine whether any of the identified transcribed regions have played a functional role in the proliferation and maintenance of the enormous non-coding regions in Silene mitochondrial genomes.

Variable frequency of plastid RNA editing among ferns and repeated loss of uridine-to-cytidine editing from vascular plants

The distinct distribution and abundance of C-to-U and U-to-C RNA editing among land plants suggest that these two processes originated and evolve independently, but the paucity of information from several key lineages limits our understanding of their evolution. To examine the evolutionary diversity of RNA editing among ferns, we sequenced the plastid transcriptomes from two early diverging species, Ophioglossum californicum and Psilotum nudum. Using a relaxed automated approach to minimize false negatives combined with manual inspection to eliminate false positives, we identified 297 C-to-U and three U-to-C edit sites in the O. californicum plastid transcriptome but only 27 C-to-U and no U-to-C edit sites in the P. nudum plastid transcriptome. A broader comparison of editing content with the leptosporangiate fern Adiantum capillus-veneris and the hornwort Anthoceros formosae uncovered large variance in the abundance of plastid editing, indicating that the frequency and type of RNA editing is highly labile in ferns. Edit sites that increase protein conservation among species are more abundant and more efficiently edited than silent and non-conservative sites, suggesting that selection maintains functionally important editing. The absence of U-to-C editing from P. nudum plastid transcripts and other vascular plants demonstrates that U-to-C editing loss is a recurrent phenomenon in vascular plant evolution.

Selection of appropriate reference genes for quantitative real-time PCR in Oxytropis ochrocephala Bunge using transcriptome datasets under abiotic stress treatments

BACKGROUND

Oxytropis ochrocephala Bunge, an indigenous locoweed species in China, poses great threats to livestock on grasslands. There is a need for further genetic study in the plants per se, for understanding the basis of its acclimation mechanism in various unfavorable environmental conditions and to implement effective control measures. Quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) is the most commonly used method for gene expression analysis. To facilitate gene expression studies and obtain more accurate qRT-PCR data, normalization relative to stable reference genes is required. The aim of this study was to select the most stable reference genes for transcriptional analysis in O. ochrocephala.

RESULTS

We selected 12 candidate reference genes, 18S ribosomal RNA (18S RNA), actin2/7 (ACT7), β-actin (ACTB), actin101 (ACT101), actin11 (ACT11), β-tubulin (TUB), α-tubulin (TUA), glyceraldehyde-3-phosphate dehydrogenase-1 (GAPDH1), GAPDH2, metallothionein-like protein (MET), fructose-bisphosphate aldolase (FBA) and histone H3 (HIS), from the transcriptome datasets of O. ochrocephala and determined the suitability by analyzing their expression levels when exposed to a range of abiotic stress conditions. By employing software packages including geNorm, NormFinder and BestKeeper, HIS, ACT7, and ACT101 were assessed as the most suitable set for normalization in all samples. When normalized with the most stable reference genes, the expression patterns of the three target genes were in accordance with those in the transcriptome data, indicating that the reference genes selected in this study are suitable.

CONCLUSIONS

The study provided appropriate reference genes for accurate normalization in qRT-PCR analysis in O. ochrocephala and emphasized the importance of validating reference genes for gene expression analysis under specific experimental condition. The usage of inappropriate reference gene would cause misinterpretation.