Search

Abstract Background Small non-coding RNAs (sRNAs) play key roles in plant development, growth and responses to biotic and abiotic stresses. At least four classes of sRNAs have been well characterized in plants, including repeat-associated siRNAs (rasiRNAs), microRNAs (miRNAs), trans-acting siRNAs (tasiRNAs) and natural antisense transcript-derived siRNAs. Chinese fir (Cunninghamia lanceolata) is one of the most important coniferous evergreen tree species in China. No sRNA from Chinese fir has been described to date. Results To obtain sRNAs in Chinese fir, we sequenced a sRNA library generated from seeds, seedlings, leaves, stems and calli, using Illumina high throughput sequencing technology. A comprehensive set of sRNAs were acquired, including conserved and novel miRNAs, rasiRNAs and tasiRNAs. With BLASTN and MIREAP we identified a total of 115 conserved miRNAs comprising 40 miRNA families and one novel miRNA with precursor sequence. The expressions of 16 conserved and one novel miRNAs and one tasiRNA were detected by RT-PCR. Utilizing real time RT-PCR, we revealed that four conserved and one novel miRNAs displayed developmental stage-specific expression patterns in Chinese fir. In addition, 209 unigenes were predicted to be targets of 30 Chinese fir miRNA families, of which five target genes were experimentally verified by 5' RACE, including a squamosa promoter-binding protein gene, a pentatricopeptide (PPR) repeat-containing protein gene, a BolA-like family protein gene, AGO1 and a gene of unknown function. We also demonstrated that the DCL3-dependent rasiRNA biogenesis pathway, which had been considered absent in conifers, existed in Chinese fir. Furthermore, the miR390-TAS3-ARF regulatory pathway was elucidated. Conclusions We unveiled a complex population of sRNAs in Chinese fir through high throughput sequencing. This provides an insight into the composition and function of sRNAs in Chinese fir and sheds new light on land plant sRNA evolution.

Abstract Background MicroRNAs (miRNAs) play key roles in diverse developmental processes, nutrient homeostasis and responses to biotic and abiotic stresses. The biogenesis and regulatory functions of miRNAs have been intensively studied in model angiosperms, such as Arabidopsis thaliana, Oryza sativa and Populus trichocarpa. However, global identification of Pinus densata miRNAs has not been reported in previous research. Results Here, we report the identification of 34 conserved miRNAs belonging to 25 miRNA families from a P. densata mRNA transcriptome database using local BLAST and MIREAP programs. The primary and/or precursor sequences of 29 miRNAs were further confirmed by RT-PCR amplification and subsequent sequencing. The average value of the minimal folding free energy indexes of the 34 miRNA precursors was 0.92. Nineteen (58%) mature miRNAs began with a 5' terminal uridine residue. Analysis of miRNA precursors showed that 19 mature miRNAs were novel members of 14 conserved miRNA families, of which 17 miRNAs were further validated by subcloning and sequencing. Using real-time quantitative RT-PCR, we found that the expression levels of 7 miRNAs were more than 2-fold higher in needles than in stems. In addition, 72 P. densata mRNAs were predicted to be targets of 25 miRNA families. Four target genes, including a nodal modulator 1-like protein gene, two GRAS family transcription factor protein genes and one histone deacetylase gene, were experimentally verified to be the targets of 3 P. densata miRNAs, pde-miR162a, pde-miR171a and pde-miR482a, respectively. Conclusions This study led to the discovery of 34 conserved miRNAs comprising 25 miRNA families from Pinus densata. These results lay a solid foundation for further studying the regulative roles of miRNAs in the development, growth and responses to environmental stresses in P. densata.

Abstract Background Digitalis purpurea is an important ornamental and medicinal plant. There is considerable interest in exploring its transcriptome. Results Through high-throughput 454 sequencing and subsequent assembly, we obtained 23532 genes, of which 15626 encode conserved proteins. We determined 140 unigenes to be candidates involved in cardiac glycoside biosynthesis. It could be grouped into 30 families, of which 29 were identified for the first time in D. purpurea. We identified 2660 mRNA-like npcRNA (mlncRNA) candidates, an emerging class of regulators, using a computational mlncRNA identification pipeline and 13 microRNA-producing unigenes based on sequence conservation and hairpin structure-forming capability. Twenty five protein-coding unigenes were predicted to be targets of these microRNAs. Among the mlncRNA candidates, only 320 could be grouped into 140 families with at least two members in a family. The majority of D. purpurea mlncRNAs were species-specific and many of them showed tissue-specific expression and responded to cold and dehydration stresses. We identified 417 protein-coding genes with regions significantly homologous or complementary to 375 mlncRNAs. It includes five genes involved in secondary metabolism. A positive correlation was found in gene expression between protein-coding genes and the homologous mlncRNAs in response to cold and dehydration stresses, while the correlation was negative when protein-coding genes and mlncRNAs were complementary to each other. Conclusions Through comprehensive transcriptome analysis, we not only identified 29 novel gene families potentially involved in the biosynthesis of cardiac glycosides but also characterized a large number of mlncRNAs. Our results suggest the importance of mlncRNAs in secondary metabolism and stress response in D. purpurea.

Aristolochia contorta Bunge is an academically and medicinally important plant species. It belongs to the magnoliids, with an uncertain phylogenetic position, and is one of the few plant species lacking a whole-genome duplication (WGD) event after the angiosperm-wide WGD. A. contorta has been an important traditional Chinese medicine material. Since it contains aristolochic acids (AAs), chemical compounds with nephrotoxity and carcinogenicity, the utilization of this plant has attracted widespread attention. Great efforts are being made to increase its bioactive compounds and reduce or completely remove toxic compounds. MicroRNAs (miRNAs) and natural antisense transcripts (NATs) are two classes of regulators potentially involved in metabolism regulation. Here, we report the identification and characterization of 223 miRNAs and 363 miRNA targets. The identified miRNAs include 51 known miRNAs belonging to 20 families and 172 novel miRNAs belonging to 107 families. A negative correlation between the expression of miRNAs and their targets was observed. In addition, we identified 441 A. contorta NATs and 560 NAT-sense transcript (ST) pairs, of which 12 NATs were targets of 13 miRNAs, forming 18 miRNA-NAT-ST modules. Various miRNAs and NATs potentially regulated secondary metabolism through the modes of miRNA-target gene–enzyme genes, NAT-STs, and NAT-miRNA-target gene–enzyme genes, suggesting the complexity of gene regulatory networks in A. contorta. The results lay a solid foundation for further manipulating the production of its bioactive and toxic compounds. [ABSTRACT FROM AUTHOR]

SUMMARY: Prunella vulgaris is one of the bestselling and widely used medicinal herbs. It is recorded as an ace medicine for cleansing and protecting the liver in Chinese Pharmacopoeia and has been used as the main constitutions of many herbal tea formulas in China for centuries. It is also a traditional folk medicine in Europe and other countries of Asia. Pentacyclic triterpenoids are a major class of bioactive compounds produced in P. vulgaris. However, their biosynthetic mechanism remains to be elucidated. Here, we report a chromosome‐level reference genome of P. vulgaris using an approach combining Illumina, ONT, and Hi‐C technologies. It is 671.95 Mb in size with a scaffold N50 of 49.10 Mb and a complete BUSCO of 98.45%. About 98.31% of the sequence was anchored into 14 pseudochromosomes. Comparative genome analysis revealed a recent WGD in P. vulgaris. Genome‐wide analysis identified 35 932 protein‐coding genes (PCGs), of which 59 encode enzymes involved in 2,3‐oxidosqualene biosynthesis. In addition, 10 PvOSC, 358 PvCYP, and 177 PvUGT genes were identified, of which five PvOSCs, 25 PvCYPs, and 9 PvUGTs were predicted to be involved in the biosynthesis of pentacyclic triterpenoids. Biochemical activity assay of PvOSC2, PvOSC4, and PvOSC6 recombinant proteins showed that they were mixed amyrin synthase (MAS), lupeol synthase (LUS), and β‐amyrin synthase (BAS), respectively. The results provide a solid foundation for further elucidating the biosynthetic mechanism of pentacyclic triterpenoids in P. vulgaris. Significance Statement: The first chromosome‐level reference genome of Prunella vulgaris and the Prunella genus was reported. Fifty‐nine 2,3‐oxidosqualene biosynthesis‐related genes and 39 pentacyclic triterpenoid biosynthesis‐related PvOSCs, PvCYPs, and PvUGTs genes were systematically analyzed. PvOSC2, PvOSC4, and PvOSC6 were experimentally verified to be mixed amyrin synthase, lupeol synthase, and β‐amyrin synthase, respectively. The results provide a solid foundation for further elucidating the biosynthetic mechanism of bioactive compounds in the widely used medicinal plant, P. vulgaris. [ABSTRACT FROM AUTHOR]