The complete chloroplast genome of the Syzygium buxifolium Hook. Et Arn.1833 and its phylogenetic analysis.

Syzygium buxifolium. Hook. Et Arn.1833 is a member of the Myrtaceae family. This species is used in traditional Chinese medicines. It possesses numerous synonyms, reflecting the ambiguity in its taxonomy. The chloroplast genome has been widely used for species identification and phylogenetic analysis. Regrettably, there is a lack of information regarding the chloroplast genome of S. buxifolium. Here, we intend to obtain the chloroplast genome of S. buxifolium to resolve its classification problems. In particular, we utilized Illumina sequencing technology to sequence, GetOrganelle to assemble, and CPGAVAS2 to characterize the chloroplast genome of S. buxifolium. The chloroplast genome of S. buxifolium had a length of 158,581 bp and consisted of 111 unique genes, comprising 78 protein-coding genes, 29 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. In addition, we identified 86 Simple Sequence Repeats, 345 tandem repetitive sequences, and 34 dispersed repetitive sequences using modules implemented in CPGAVAS2. Lastly, we carried out phylogenetic analysis using Phylosuite. The results indicated a close relationship between S. buxifolium and S. grijsii. This study offers novel genetic data for the molecular identification and subsequent phylogenetic analysis of the Syzygium genus.

Comparative chloroplast genomes study of five officinal Ardisia Species: Unraveling interspecific diversity and evolutionary insights in Ardisia.

Ardisia S.W. (Primulaceae), naturally distributed in tropical and subtropical regions, has edible and medicinal values and is prevalent in clinical and daily use in China. More genetic information for distinct species delineation is needed to support the development and utilization of the genus Ardisia. We sequenced, annotated, and compared the chloroplast genomes of five Ardisia species: A. brunnescens, A. pusilla, A. squamulosa, A. crenata, and A. brevicaulis in this study. We found a typical quadripartite structure in all five chloroplast genomes, with lengths ranging from 155,045 to 156,943 bp. Except for A. pusilla, which lacked the ycf15 gene, the other four Ardisia species contained 114 unique genes, including 79 protein-coding genes, 30 tRNAs, and four rRNAs. In addition, the rps19 pseudogene gene was present only in A. brunnescens. Five highly variable DNA barcodes were identified for five Ardisia species, including trnT-GGU-psbD, trnT-UGU-trnL-UAA, rps4-trnT-UGU, rpl32-trnL-UAG, and rpoB-trnC-GAA. The RNA editiing sites of protein-coding genes in the five Ardisia plastome were characterized and compared, and 274 (A. crenata)-288 (A. brevicaulis) were found. The results of the phylogenetic analysis were consistent with the morphological classification. Sequence alignment and phylogenetic analysis showed that ycf15 genes were highly divergent in Primulaceae. Reconstructions of ancestral character states indicated that leaf margin morphology is critical for classifying the genus Ardisia, with a rodent-like character being the most primitive. These results provide valuable information on the taxonomy and evolution of Ardisia plants.

Genome-wide Identification and Characterization of Natural Antisense Transcripts by Strand-specific RNA Sequencing in Ganoderma lucidum.

Ganoderma lucidum is a white-rot fungus best-known for its medicinal and ligninolytic activities. To discover the underlying genes responsible for these activities, we identified and characterized the natural antisense transcripts (NATs) using strand-specific (ss) RNA-seq data obtained from the mycelia, primordia and fruiting bodies. NATs were identified using a custom pipeline and then subjected to functional enrichment and differential expression analyses. A total of 1613 cis- and 244 trans- sense and antisense transcripts were identified. Mapping to GO terms and KEGG pathways revealed that NATs were frequently associated with genes of particular functional categories in particular stages. ssRT-qPCR experiments showed that the expression profiles of 30 of 50 (60%) transcripts were highly correlated with those of the RNA-seq results (r ≥ 0.9). Expression profiles of 22 of 25 (88%) pairs of NATs and STs were highly correlated (p ≤ 0.01), with 15 having r ≥ 0.8 and 4 having r ≤ -0.8. Six lignin-modifying genes and their NATs were analyzed in detail. Diverse patterns of differential expression among different stages and positive and negative correlations were observed. These results suggested that NATs were implicated in gene expression regulation in a function-group and developmental-stage specific manner through complex mechanisms.

Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa.

Laccases, as early as 1959, were proposed to catalyze the oxidative polymerization of monolignols. Genetic evidence in support of this hypothesis has been elusive due to functional redundancy of laccase genes. An Arabidopsis double mutant demonstrated the involvement of laccases in lignin biosynthesis. We previously identified a subset of laccase genes to be targets of a microRNA (miRNA) ptr-miR397a in Populus trichocarpa. To elucidate the roles of ptr-miR397a and its targets, we characterized the laccase gene family and identified 49 laccase gene models, of which 29 were predicted to be targets of ptr-miR397a. We overexpressed Ptr-MIR397a in transgenic P. trichocarpa. In each of all nine transgenic lines tested, 17 PtrLACs were down-regulated as analyzed by RNA-seq. Transgenic lines with severe reduction in the expression of these laccase genes resulted in an ∼40% decrease in the total laccase activity. Overexpression of Ptr-MIR397a in these transgenic lines also reduced lignin content, whereas levels of all monolignol biosynthetic gene transcripts remained unchanged. A hierarchical genetic regulatory network (GRN) built by a bottom-up graphic Gaussian model algorithm provides additional support for a role of ptr-miR397a as a negative regulator of laccases for lignin biosynthesis. Full transcriptome-based differential gene expression in the overexpressed transgenics and protein domain analyses implicate previously unidentified transcription factors and their targets in an extended hierarchical GRN including ptr-miR397a and laccases that coregulate lignin biosynthesis in wood formation. Ptr-miR397a, laccases, and other regulatory components of this network may provide additional strategies for genetic manipulation of lignin content.

High-throughput sequencing and characterization of the small RNA transcriptome reveal features of novel and conserved microRNAs in Panax ginseng.

microRNAs (miRNAs) play vital regulatory roles in many organisms through direct cleavage of transcripts, translational repression, or chromatin modification. Identification of miRNAs has been carried out in various plant species. However, no information is available for miRNAs from Panax ginseng, an economically significant medicinal plant species. Using the next generation high-throughput sequencing technology, we obtained 13,326,328 small RNA reads from the roots, stems, leaves and flowers of P. ginseng. Analysis of these small RNAs revealed the existence of a large, diverse and highly complicated small RNA population in P. ginseng. We identified 73 conserved miRNAs, which could be grouped into 33 families, and 28 non-conserved ones belonging to 9 families. Characterization of P. ginseng miRNA precursors revealed many features, such as production of two miRNAs from distinct regions of a precursor, clusters of two precursors in a transcript, and generation of miRNAs from both sense and antisense transcripts. It suggests the complexity of miRNA production in P. ginseng. Using a computational approach, we predicted for the conserved and non-conserved miRNA families 99 and 31 target genes, respectively, of which eight were experimentally validated. Among all predicted targets, only about 20% are conserved among various plant species, whereas the others appear to be non-conserved, indicating the diversity of miRNA functions. Consistently, many miRNAs exhibited tissue-specific expression patterns. Moreover, we identified five dehydration- and ten heat-responsive miRNAs and found the existence of a crosstalk among some of the stress-responsive miRNAs. Our results provide the first clue to the elucidation of miRNA functions in P. ginseng.

Genome-wide identification and characterization of novel genes involved in terpenoid biosynthesis in Salvia miltiorrhiza.

Terpenoids are the largest class of plant secondary metabolites and have attracted widespread interest. Salvia miltiorrhiza, belonging to the largest and most widely distributed genus in the mint family, is a model medicinal plant with great economic and medicinal value. Diterpenoid tanshinones are the major lipophilic bioactive components in S. miltiorrhiza. Systematic analysis of genes involved in terpenoid biosynthesis has not been reported to date. Searching the recently available working draft of the S. miltiorrhiza genome, 40 terpenoid biosynthesis-related genes were identified, of which 27 are novel. These genes are members of 19 families, which encode all of the enzymes involved in the biosynthesis of the universal isoprene precursor isopentenyl diphosphate and its isomer dimethylallyl diphosphate, and two enzymes associated with the biosynthesis of labdane-related diterpenoids. Through a systematic analysis, it was found that 20 of the 40 genes could be involved in tanshinone biosynthesis. Using a comprehensive approach, the intron/exon structures and expression patterns of all identified genes and their responses to methyl jasmonate treatment were analysed. The conserved domains and phylogenetic relationships among the deduced S. miltiorrhiza proteins and their homologues isolated from other plant species were revealed. It was discovered that some of the key enzymes, such as 1-deoxy-D-xylulose 5-phosphate synthase, 4-hydroxy-3-methylbut-2-enyl diphosphate reductase, hydroxymethylglutaryl-CoA reductase, and geranylgeranyl diphosphate synthase, are encoded by multiple gene members with different expression patterns and subcellular localizations, and both homomeric and heteromeric geranyl diphosphate synthases exist in S. miltiorrhiza. The results suggest the complexity of terpenoid biosynthesis and the existence of metabolic channels for diverse terpenoids in S. miltiorrhiza and provide useful information for improving tanshinone production through genetic engineering.

Comprehensive transcriptome analysis reveals novel genes involved in cardiac glycoside biosynthesis and mlncRNAs associated with secondary metabolism and stress response in Digitalis purpurea.

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.