Search

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analyzed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes had always higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25 and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.

Salvia miltiorrhiza is well-known for its clinical practice in treating heart and cardiovascular diseases. Its roots used for traditional Chinese medicine materials are usually brick-red due to accumulation of red pigments, such as tanshinone IIA and tanshinone I. Here we report a S. miltiorrhiza line (shh) with orange roots. Compared with the red roots of normal S. miltiorrhiza plants, the contents of tanshinones with a single bond at C-15,16 were increased, whereas those with a double bond at C-15,16 were significantly decreased in shh. We assembled a high-quality chromosome-level genome of shh. Phylogenomic analysis showed that the relationship between two S. miltiorrhiza lines with red roots were closer than their relationship with shh. It indicates that shh could be not the mutant of an extant S. miltiorrhiza line with red roots. Comparative genomic and transcriptomic analyses showed that a 1.0kb DNA fragment was deleted in shh Sm2OGD3m. Complementation assay showed that overexpression of intact Sm2OGD3 in shh hairy roots recovered furan D-ring tanshinone accumulation. Consistently, in vitro protein assay showed that Sm2OGD3 catalyzed the conversion of cyptotanshinone, 15,16-dihydrotanshinone I and 1,2,15,16-tetrahydrotanshinone I into tanshinone IIA, tanshinone I and 1,2-dihydrotanshinone I, respectively. Thus, Sm2OGD3 functions as tanshinone 15,16-dehydrogenase and is a key enzyme in tanshinone biosynthesis. The results provide novel insights into the metabolic network of medicinally important tanshinone compounds.

Salvia miltiorrhiza is one of the most well known species in the genus Salvia of the Lamiaceae with great economic, academic and medicinal value. It was recorded as a traditional Chinese medicine m...

Background: Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field. Objective: This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production. Results: So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production. Conclusion: The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Aristolochic acids (AAs) and their derivatives are present in multiple Aristolochiaceae species that have been or are being used as medicinal materials. During the past decades, AAs have received increasing attention because of their nephrotoxicity and carcinogenicity. Elimination of AAs from medicinal materials using biotechnological approaches is important for improving medication safety. However, it has not been achieved because of the limited information available on AA biosynthesis. Here, we report a high-quality, reference-grade genome assembly of the AA-containing vine Aristolochia contorta. The total size of the assembly is 209.27 Mb, and it is assembled into 7 pseudochromosomes. Synteny analysis, Ks distribution, and 4DTv suggest an absence of whole-genome duplication (WGD) events in Aristolochia contorta after the angiosperm-wide WGD. Based on genomic, transcriptomic, and metabolic data, pathways and candidate genes were proposed for benzylisoquinoline alkaloid (BIA) and AA biosynthesis in A. contorta. Five O-methyltransferase genes, including AcOMT1–3, AcOMT5, and AcOMT7, were cloned and functionally characterized. The results provide a high-quality reference genome for AA-containing species of Aristolochiaceae. They lay a solid foundation for further elucidation of AA biosynthesis and regulation and for the molecular breeding of Aristolochiaceae medicinal materials.

Tanshinones are important bioactive components in Salvia miltiorrhiza and mainly accumulate in the periderms of mature roots. Tanshinone biosynthesis is a complicated process, and little is known about the third stage of the pathway. To investigate potential genes that are responsible for tanshinone biosynthesis, we conducted transcriptome profiling analysis of two S. miltiorrhiza cultivars. Differential expression analysis provided 2,149 differentially expressed genes (DEGs) for further analysis. GO and KEGG analysis showed that the DEGs were mainly associated with the biosynthesis of secondary metabolites. Weighted gene coexpression network analysis (WGCNA) was further performed to identify a “cyan” module associated with tanshinone biosynthesis. In this module, 25 cytochromes P450 (CYPs), three 2-oxoglutarate-dependent dioxygenases (2OGDs), one short-chain alcohol dehydrogenases (SDRs) and eight transcription factors were found to be likely involved in tanshinone biosynthesis. Among these CYPs, 14 CYPs have been reported previously, and 11 CYPs were identified in this study. Expression analysis showed that four newly identified CYPs were upregulated upon application of MeJA, suggesting their possible roles in tanshinone biosynthesis. Overall, this study not only identified candidate genes involved in tanshinone biosynthesis but also provided a basis for characterization of genes involved in important active ingredients of other traditional Chinese medicinal plants.

In this chapter, the discovery of noncoding RNAs, including lincRNA, natural antisense noncoding RNAs, circular RNAs, and microRNAs in Ganoderma lucidum is described. In particular, we describe the methods for validation and discovery of these noncoding RNAs, their characteristics, and their potential functions.

Alkaline/neutral invertases (NINVs), which irreversibly catalyze the hydrolysis of sucrose into fructose and glucose, play crucial roles in carbohydrate metabolism and plant development. Comprehensive insights into NINV genes are lacking in Salvia miltiorrhiza, a well-known traditional Chinese medicinal (TCM) plant with significant medicinal and economic value. Through genome-wide prediction, nine putative SmNINV genes, termed SmNINV1-SmNINV9, were identified. Integrated analysis of gene structures, sequence features, conserved domains, conserved motifs and phylogenetic trees revealed the conservation and divergence of SmNINVs. The identified SmNINVs were differentially expressed in roots, stems, leaves, flowers, and different root tissues. They also responded to drought, salicylic acid, yeast extract, and methyl jasmonate treatments. More importantly, computational prediction and experimental validation showed that SmNINV3 and SmNINV4 were targets of Smi-miR399, a conserved miRNA previously shown to affect Pi uptake and translocation through the cleavage of PHOSPHATE2 (PHO2). Consistently, analysis of 43 NINV genes and 26 miR399 sequences from Arabidopsis thaliana, Populus trichocarpa, Manihot esculenta, and Solanum lycopersicum showed that various AtNINV, PtNINV, MeNINV, and SlNINV genes were regulated by miR399. It indicates that the miR399-NINV module exists widely in plants. Furthermore, Smi-miR399 also cleaved SmPHO2 transcripts in S. miltiorrhiza, suggesting the complexity of NINVs, PHO2, and miR399 networks.

SmPPT, which encodes 4-hydroxybenzoate polyprenyl diphosphate transferase involved in ubiquinone biosynthesis, confers salt tolerance to S. miltiorrhiza through enhancing the activities of POD and CAT to scavenge ROS. Ubiquinone (UQ), also known as coenzyme Q (CoQ), is a key electron transporter in the mitochondrial respiratory system. UQ is composed of a benzene quinone ring and a polyisoprenoid side chain. Attachment of polyisoprenoid side chain to the benzene quinone ring is a rate-limiting step catalyzed by 4-hydroxybenzoate polyprenyl diphosphate transferase (PPT). So far, only a few plant PPT-encoding genes have been functionally analyzed. Through genome-wide analysis and subsequent molecular cloning, a PPT-encoding gene, termed SmPPT, was identified from an economically and academically important medicinal model plant, Salvia miltiorrhiza. SmPPT contained many putative cis-elements associated with abiotic stresses in the promoter region and were responsive to PEG-6000 and methyl jasmonate treatments. The deduced SmPPT protein contains the PT_UbiA conserved domain of polyprenyl diphosphate transferase and an N-terminal mitochondria transit peptide. Transient expression assay of SmPPT-GFP fusion protein showed that SmPPT was mainly localized in the mitochondria. SmPPT could functionally complement coq2 mutation and catalyzed UQ6 production in yeast cells. Overexpression of SmPPT increased UQ production and enhanced salt tolerance in S. miltiorrhiza. Under salinity stress conditions, transgenic plants accumulated less H2O2 and malondialdehyde and exhibited higher peroxidase (POD) and catalase (CAT) activities compared with wild-type plants. It indicates that SmPPT confers salt tolerance to S. miltiorrhiza at least partially through enhancing the activities of POD and CAT to scavenge ROS.

MIRNA (MIR) gene origin and early evolutionary processes, such as hairpin precursor sequence origination, promoter activity acquirement and the sequence of these two processes, are fundamental and fascinating subjects. Three models, including inverted gene duplication, spontaneous evolution and transposon transposition, have been proposed for de novo origination of hairpin precursor sequence. However, these models still open to discussion. In addition, de novo origination of MIR gene promoters has not been well investigated. Here, I systematically investigated the origin of evolutionarily young polyphenol oxidase gene (PPO)-targeting MIRs, including MIR1444, MIR058 and MIR12112, and a genomic region termed AasPPO-as-hp, which contained a hairpin-forming sequence. I found that MIR058 precursors and the hairpin-forming sequence of AasPPO-as-hp originated in an ancient PPO gene through forming short inverted repeats. Palindromic-like sequences and imperfect inverted repeats in the ancient PPO gene contributed to initiate the generation of short inverted repeats probably by causing errors during DNA duplication. Analysis of MIR058 and AasPPO-as-hp promoters showed that they originated in the 3ʹ-flanking region of the ancient PPO gene. Promoter activities were gained by insertion of a CAAT-box and multiple-copper-response element (CuRE)-containing miniature inverted-repeat transposable element (MITE) in the upstream of AT-rich TATA-box-like sequence. Gain of promoter activities occurred before hairpin-forming sequence origination. Sequence comparison of MIR1444, MIR058 and MIR12112 promoters showed frequent birth and death of CuREs, indicating copper could be vital for the origination and evolution of PPO-targeting MIRs. Based on the evidence obtained, a novel model for plant MIR origination and evolution is proposed.

Salvia miltiorrhiza Bunge is a perennial plant in the genus Salvia of the family Lamiaceae. It has a long medicinal history with the first recorded in about the second century BC. S. miltiorrhiza has great economic and medicinal value. Over 10% of the traditional Chinese patent medicines and simple preparations contain S. miltiorrhiza. Annual consumption of S. miltiorrhiza has exceeded 16 million kilograms in China. So far, more than two hundred chemical compounds have been isolated from this plant. S. miltiorrhiza and its ingredients have been clinically used for managing vascular diseases and academically shown potential in treating various other diseases. Moreover, S. miltiorrhiza has great academic value and is emerging as a model system in medicinal plant biology. This plant has a relatively short life cycle and a relatively small genome size. It is easy to propagate and cultivate. There are whole genome and transcriptome sequences, in vitro tissue culture systems, hairy root and crown gall induction systems and genetic transformation systems available for the research community of S. miltiorrhiza. This chapter provides an overview of the taxonomy, distribution, morphology, propagation, growth requirements, chemical constituents, and medicinal and academic significance of S. miltiorrhiza.

MicroRNAs (miRNAs in animals and plants or milRNAs in fungi) are endogenous noncoding RNAs that can regulate gene expression. However, little information is known about milRNAs and their target genes in Ganoderma lucidum. Here, we systematically predicted and characterised the milRNAs and their target genes across the three developmental stages of G. lucidum. A total of 168 unique milRNAs were predicted using a small RNA sequencing method. For them, 1612 target sequences corresponding to 1311 unique genes were predicted by degradome sequencing. We selected 42 predicted milRNAs and performed RT-PCR amplification and Sanger sequencing of the products. Five products were found to have sequences similar to those predicted, confirming the presence of milRNAs in G. lucidum, and demonstrating the difficulty in their validation. Among the 168 milRNAs, 111 were found to be significantly differentially expressed across the three developmental stages (q ≤ 0.05). The expression levels of 12 milRNAs were measured by stem-loop quantitative real-time polymerase chain reaction. Eight of them were in line with the sequencing results (r ≥ 0.9, p ≤ 0.05). These 12 milRNAs and their target genes form 16 milRNA-target gene pairs. The expression profiles of 8 of these 16 miRNA-target pairs were negatively correlated, according to real-time quantitative analysis, whereas the other eight pairs were positively correlated. Furthermore, the results of functional enrichment analysis showed that the target genes of milRNAs mapped to the Gene Ontology terms ‘GTP binding’ and ‘FAD binding’ were enriched in specific developmental stages. These target genes were related to the biosynthesis of triterpenes and polysaccharides and lignin degradation pathway in G. lucidum. In summary, this study indicates that milRNAs may play crucial regulatory roles in various biological processes of G. lucidum and open up new avenues for research on milRNAs’ biosyntheses and functions in basidiomycetes.

Flavonoids are a class of important secondary metabolites with a broad spectrum of pharmacological functions. Salvia miltiorrhiza Bunge (Danshen) is a well-known traditional Chinese medicinal herb with a broad diversity of flavonoids. However, flavonoid biosynthetic enzyme genes have not been systematically and comprehensively analyzed in S. miltiorrhiza. Through genome-wide prediction and molecular cloning, twenty six flavonoid biosynthesis-related gene candidates were identified, of which twenty are novel. They belong to nine families potentially encoding chalcone synthase (CHS), chalcone isomerase (CHI), flavone synthase (FNS), flavanone 3-hydroxylase (F3H), flavonoid 3′-hydroxylase (F3′H), flavonoid 3′,5′-hydroxylase (F3′5′H), flavonol synthase (FLS), dihydroflavonol 4-reductase (DFR), and anthocyanidin synthase (ANS), respectively. Analysis of intron/exon structures, features of deduced proteins and phylogenetic relationships revealed the conservation and divergence of S. miltiorrhiza flavonoid biosynthesis-related proteins and their homologs from other plant species. These genes showed tissue-specific expression patterns and differentially responded to MeJA treatment. Through comprehensive and systematic analysis, fourteen genes most likely to encode flavonoid biosynthetic enzymes were identified. The results provide valuable information for understanding the biosynthetic pathway of flavonoids in medicinal plants.

Cytosine DNA methylation is highly conserved epigenetic modification involved in a wide range of biological processes in eukaryotes. It was established and maintained by cytosine-5 DNA methyltransferases (C5-MTases) in plants. Through genome-wide identification, eight putativeSmC5-MTasegenes were identified from the genome ofSalvia miltiorrhiza, a well-known traditional Chinese medicine material and an emerging model medicinal plant. Based on conserved domains and phylogenetic analysis, eightSmC5-MTasegenes were divided into four subfamilies, includingMET,CMT,DRMandDNMT2. Genome-wide comparative analysis of theC5-MTasegene family inS. miltiorrhizaandArabidopsis thaliana, including gene structure, sequence features, sequence alignment and conserved motifs, was carried out. The results showed conservation and divergence of the members of each subfamily in plants. The length ofSmC5-MTaseopen reading frames ranges widely from 1,152 (SmDNMT2) to 5,034 bp (SmMET1). The intron number ofSmC5-MTasesvaries between 7 (SmDRM1) and 20 (SmCMT1andSmCMT2b). These features were similar to their counterparts fromArabidopsis. Sequence alignment and conserved motif analysis showed the existence of highly conserved and subfamily-specific motifs in the C5-MTases analyzed. Differential transcript abundance was detected forSmC5-MTases, implying genome-wide variance of DNA methylation in different organs and tissues. Transcriptome-wide analysis showed that the transcript levels of allSmC5-MTasegenes was slightly changed under yeast extract and methyl jasmonate treatments. SixSmC5-MTases, includingSmMET1,SmCMT1,SmCMT2a,SmCMT2b,SmCMT3andSmDRM1, were salicylic acid-responsive, suggesting the involvement ofSmC5-MTases in salicylic acid-dependent immunity. These results provide useful information for demonstrating the role of DNA methylation in bioactive compound biosynthesis and Dao-di herb formation in medicinal plants.

Cytosine DNA methylation is highly conserved epigenetic modification involved in a wide range of biological processes in eukaryotes. It was established and maintained by cytosine-5 DNA methyltransferases (C5-MTases) in plants. Through genome-wide identification, eight putative SmC5-MTase genes were identified from the genome of Salvia miltiorrhiza, a well-known traditional Chinese medicine material and an emerging model medicinal plant. Based on conserved domains and phylogenetic analysis, eight SmC5-MTase genes were divided into four subfamilies, including MET, CMT, DRM and DNMT2. Genome-wide comparative analysis of the C5-MTase gene family in S. miltiorrhiza and Arabidopsis thaliana, including gene structure, sequence features, sequence alignment and conserved motifs, was carried out. The results showed conservation and divergence of the members of each subfamily in plants. The length of SmC5-MTase open reading frames ranges widely from 1,152 (SmDNMT2) to 5,034 bp (SmMET1). The intron number of SmC5-MTases varies between 7 (SmDRM1) and 20 (SmCMT1 and SmCMT2b). These features were similar to their counterparts from Arabidopsis. Sequence alignment and conserved motif analysis showed the existence of highly conserved and subfamily-specific motifs in the C5-MTases analyzed. Differential transcript abundance was detected for SmC5-MTases, implying genome-wide variance of DNA methylation in different organs and tissues. Transcriptome-wide analysis showed that the transcript levels of all SmC5-MTase genes was slightly changed under yeast extract and methyl jasmonate treatments. Six SmC5-MTases, including SmMET1, SmCMT1, SmCMT2a, SmCMT2b, SmCMT3 and SmDRM1, were salicylic acid-responsive, suggesting the involvement of SmC5-MTases in salicylic acid-dependent immunity. These results provide useful information for demonstrating the role of DNA methylation in bioactive compound biosynthesis and Dao-di herb formation in medicinal plants.

MicroRNAs (miRNAs) are small noncoding regulatory RNAs that play key roles in the process of plant development. To date, extensive studies of miRNAs have been performed in a few model plants, but few efforts have focused on small RNAs (sRNAs) in conifers because of the lack of reference sequences for their enormous genomes. In this study, Solexa sequencing of three sRNA libraries obtained from dormant, reactivating, and active vascular cambium in Chinese fir (Cunninghamia lanceolata) using tangential cryosectioning identified 20 known miRNA families and 18 novel potential miRNAs, of which nine novel miRNA precursors were validated by RT-PCR and sequencing. More than half of these novel miRNAs displayed stage-specific expression patterns in the vascular cambium. Furthermore, analysing the 103 miRNAs and their predicted targets indicated that about 70% appeared to negatively regulate their targets, of which two target genes involved in the regulation of cambial cell division were validated via RNA ligase-mediated rapid amplification of 5'-cDNA ends (RLM 5'-RACE) and transient co-expression in Nicotiana benthamiana leaves. Interestingly, miRNA156 and miRNA172 may regulate the phase transition in vascular cambium from dormancy to active growth. These results provide new insights into the important regulatory functions of miRNAs in vascular cambium development and wood formation in conifers.

DEMETER-like DNA glycosylases (DMLs) initiate the base excision repair-dependent DNA demethylation to regulate a wide range of biological processes in plants. Six putative SmDML genes, termed SmDML1–SmDML6, were identified from the genome of S. miltiorrhiza, an emerging model plant for Traditional Chinese Medicine (TCM) studies. Integrated analysis of gene structures, sequence features, conserved domains and motifs, phylogenetic analysis and differential expression showed the conservation and divergence of SmDMLs. SmDML1, SmDML2 and SmDML4 were significantly down-regulated by the treatment of 5Aza-dC, a general DNA methylation inhibitor, suggesting involvement of SmDMLs in genome DNA methylation change. SmDML1 was predicted and experimentally validated to be target of Smi-miR7972. Computational analysis of forty whole genome sequences and almost all of RNA-seq data from Lamiids revealed that MIR7972s were only distributed in some plants of the three orders, including Lamiales, Solanales and Boraginales, and the number of MIR7972 genes varied among species. It suggests that MIR7972 genes underwent expansion and loss during the evolution of some Lamiids species. Phylogenetic analysis of MIR7972s showed closer evolutionary relationships between MIR7972s in Boraginales and Solanales in comparison with Lamiales. These results provide a valuable resource for elucidating DNA demethylation mechanism in S. miltiorrhiza.

Background FL 32608 Terpenoids are hydrocarbon compounds derived from common fivecarbon isoprene (C5H8) building blocks. They are formed through the condensation and subsequent modification of isoprene units in various ways including - among others - cyclization and/or oxygenation. Their synthesis is localized either to the chloroplast and/or to the cytoplasm/peroxisome/ endoplasmic reticulum and mitochondrion. Terpenoids represent a very large and diverse class of metabolites and play important roles in plant growth and development. In addition, they have been intensively used in human health care, disease treatment and in dietary supplements. Approximately 60% of natural products known so far are terpenoids. Conclusion This review briefly summarizes the biosynthetic pathways of major plant terpenoids. Then, five well-known and medicinally important diterpenoids, including paclitaxel, tanshinone, ginkgolide, triptolide and oridonin are discussed in detail. Their structures, occurrence, extraction and identification methods, pharmacological properties and clinical uses are also reviewed. Finally, the prospects of using biotechnology to produce medicinally important terpenoids are also briefly discussed.

Small RNA-mediated gene silencing is a vital regulatory mechanism in eukaryotes that requires ARGONAUTE (AGO) proteins. Salvia miltiorrhiza is a well-known traditional Chinese medicinal plant. Therefore, it is important to characterize S. miltiorrhiza AGO family genes as they may be involved in multiple metabolic pathways. This chapter introduces the detailed protocol for SmAGO gene prediction and molecular cloning. In addition, an Agrobacterium-mediated genetic transformation method for S. miltiorrhiza is presented. These methodologies can be used to functionally study SmAGO genes as well as other genes of interest in S. miltiorrhiza.

miR1444s are functionally significant miRNAs targeting polyphenol oxidase (PPO) genes for cleavage. MIR1444 genes were reported only in Populus trichocarpa. Through the computational analysis of 215 RNA-seq data, four whole genome sequences of Salicaceae species and deep sequencing of six P. trichocarpa small RNA libraries, we investigated the origin and evolution history of MIR1444s. A total of 23 MIR1444s were identified. Populus and Idesia species contain two MIR1444 genes, while Salix includes only one. Populus and Idesia MIR1444b genes and Salix MIR1444s were phylogenetically separated from Populus and Idesia MIR1444a genes. Ptr-miR1444a and ptr-miR1444b showed sequence divergence. Compared with ptr-miR1444b, ptr-miR1444a started 2 nt upstream of precursor, resulting in differential regulation of PPO targets. Sequence alignments showed that MIR1444 genes exhibited extensive similarity to their PPO targets, the characteristics of MIRs originated from targets through an inverted gene duplication event. Genome sequence comparison showed that MIR1444 genes in Populus and Idesia were expanded through the Salicoid genome duplication event. A copy of MIR1444 gene was lost in Salix through DNA segment deletion during chromosome rearrangements. The results provide significant information for the origin of plant miRNAs and the mechanism of Salicaceae gene evolution and divergence.

Small RNA-mediated gene silencing is a vital regulatory mechanism in eukaryotes that requires ARGONAUTE (AGO) proteins. Salvia miltiorrhiza is a well-known traditional Chinese medicinal plant. Therefore, it is important to characterize S. miltiorrhiza AGO family genes as they may be involved in multiple metabolic pathways. This chapter introduces the detailed protocol for SmAGO gene prediction and molecular cloning. In addition, an Agrobacterium-mediated genetic transformation method for S. miltiorrhiza is presented. These methodologies can be used to functionally study SmAGO genes as well as other genes of interest in S. miltiorrhiza.

Increasing evidence suggests that long non-coding RNAs (lncRNAs) play significant roles in plants. However, little is known about lncRNAs in Panax ginseng C. A. Meyer, an economically significant medicinal plant species. A total of 3,688 mRNA-like non-coding RNAs (mlncRNAs), a class of lncRNAs, were identified in P. ginseng. Approximately 40% of the identified mlncRNAs were processed into small RNAs, implying their regulatory roles via small RNA-mediated mechanisms. Eleven miRNA-generating mlncRNAs also produced siRNAs, suggesting the coordinated production of miRNAs and siRNAs in P. ginseng. The mlncRNA-derived small RNAs might be 21-, 22-, or 24-nt phased and could be generated from both or only one strand of mlncRNAs, or from super long hairpin structures. A full-length mlncRNA, termed MAR (multiple-function-associated mlncRNA), was cloned. It generated the most abundant siRNAs. The MAR siRNAs were predominantly 24-nt and some of them were distributed in a phased pattern. A total of 228 targets were predicted for 71 MAR siRNAs. Degradome sequencing validated 68 predicted targets involved in diverse metabolic pathways, suggesting the significance of MAR in P. ginseng. Consistently, MAR was detected in all tissues analyzed and responded to methyl jasmonate (MeJA) treatment. It sheds light on the function of mlncRNAs in plants.

SQUAMOSA promoter binding protein-likes (SPLs) are plant-specific transcription factors playing vital regulatory roles in plant growth and development. There is no information about SPLs in Salvia miltiorrhiza (Danshen), a significant medicinal plant widely used in Traditional Chinese medicine (TCM) for >1,700 years and an emerging model plant for TCM studies. Through genome-wide identification and subsequent molecular cloning, we identified a total 15 SmSPLs with divergent sequence features, gene structures, and motifs. Comparative analysis showed sequence conservation between SmSPLs and their Arabidopsis counterparts. A phylogenetic tree clusters SmSPLs into six groups. Many of the motifs identified commonly exist in a group/subgroup, implying their functional redundancy. Eight SmSPLs were predicted and experimentally validated to be targets of miR156/157. SmSPLs were differentially expressed in various tissues of S. milltiorrhiza. The expression of miR156/157-targeted SmSPLs was increased with the maturation of S. miltiorrhiza, whereas the expression of miR156/157 was decreased, confirming the regulatory roles of miR156/157 in SmSPLs and suggesting the functions of SmSPLs in S. miltiorrhiza development. The expression of miR156/157 was negatively correlated with miR172 during the maturation of S. miltiorrhiza. The results indicate the significance and complexity of SmSPL-, miR156-, and miR172-mediated regulation of developmental timing in S. miltiorrhiza.

Salvia miltiorrhiza is a well-known material of traditional Chinese medicine. Understanding the regulatory mechanisms of phenolic acid biosynthesis and metabolism are important for S. miltiorrhiza quality improvement. We report here that S. miltiorrhiza contains 19 polyphenol oxidases (PPOs), forming the largest PPO gene family in plant species to our knowledge. Analysis of gene structures and sequence features revealed the conservation and divergence of SmPPOs. SmPPOs were differentially expressed in plant tissues and eight of them were predominantly expressed in phloem and xylem, indicating that some SmPPOs are functionally redundant, whereas the others are associated with different physiological processes. Expression patterns of eighteen SmPPOs were significantly altered under MeJA treatment, and twelve were yeast extract and Ag+-responsive, suggesting the majority of SmPPOs are stress-responsive. Analysis of high-throughput small RNA sequences and degradome data showed that miR1444-mediated regulation of PPOs existing in P. trichocarpa is absent from S. miltiorrhiza. Instead, a subset of SmPPOs was posttranscriptionally regulated by a novel miRNA, termed Smi-miR12112. It indicates the specificity and significance of miRNA-mediated regulation of PPOs. The results shed light on the regulation of SmPPO expression and suggest the complexity of SmPPO-associated phenolic acid biosynthesis and metabolism.

Reverse transcription-qPCR (RT-qPCR) has become a popular method for gene expression studies. Its results require data normalization by housekeeping genes. No single gene is proved to be stably expressed under all experimental conditions. Therefore, systematic evaluation of reference genes is necessary. With the aim to identify optimum reference genes for RT-qPCR analysis of gene expression in different tissues of Panax ginseng and the seedlings grown under heat stress, we investigated the expression stability of eight candidate reference genes, including elongation factor 1-beta (EF1-β), elongation factor 1-gamma (EF1-γ), eukaryotic translation initiation factor 3G1 (IF3G1), eukaryotic translation initiation factor 3B (IF3B), actin (ACT), actin11 (ACT11), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and cyclophilin ABH-like protein (CYC), using four widely used computational programs: geNorm, Normfinder, BestKeeper, and the comparative ΔCt method. The results were then integrated using the web-based tool RefFinder. As a result, EF1-γ, IF3G1, and EF1-β were the three most stable genes in different tissues of P. ginseng, while IF3G1, ACT11, and GAPDH were the top three-ranked genes in seedlings treated with heat. Using three better reference genes alone or in combination as internal control, we examined the expression profiles of MAR, a multiple function-associated mRNA-like non-coding RNA (mlncRNA) in P. ginseng. Taken together, we recommended EF1-γ/IF3G1 and IF3G1/ACT11 as the suitable pair of reference genes for RT-qPCR analysis of gene expression in different tissues of P. ginseng and the seedlings grown under heat stress, respectively. The results serve as a foundation for future studies on P. ginseng functional genomics.

Background Gibberellin (GA), a classical phytohormone, plays significant roles in plant growth and development. It shares the important intermediate diphosphate precursor, GGPP, with the main lipophilic bioactive components, diterpenoid tanshinones in Salvia miltiorrhiza Bunge, one of the most important Traditional Chinese Medicine materials and an emerging model medicinal plant. Analysis of GA metabolism and regulation may help to demonstrate the biological functions of GAs and the crosstalk between GA metabolism and tanshinone biosynthesis in S. miltiorrhiza. However, genes involved in the conversion of ent-kaurene to GAs have not been systematically studied. Results Through genome-wide prediction and molecular cloning, twenty two candidate gibberellin metabolism pathway genes were systematically identified for the first time. It includes a SmKO, two SmKAOs, six SmGA20oxs, two SmGA3oxs and eleven SmGA2oxs, of which twenty genes are novel. The deduced proteins showed sequence conservation and divergence. Gibberellin metabolism pathway genes exhibited tissue-specific expression patterns and responded differentially to exogenous GA3 treatment, indicating differential regulation of gibberellin metabolism in different tissue types in S. miltiorrhiza. SmKAO1, SmKAO2, SmGA2ox2, and SmGA2ox4–SmGA2ox7 were significantly up-regulated; SmGA20ox2, SmGA3ox1, SmGA2ox1, SmGA2ox8, SmGA2ox10 and SmGA2ox11 were significantly down-regulated; while the responses of many other genes varied among different tissue-types and time-points of GA3 treatment, suggesting the complexity of feedback regulation. Tanshinone biosynthesis-related genes, such as SmCPS1 and SmKSL1, were up-regulated in response to GA3 treatment. Among the 22 identified genes, nine responded to yeast extract and Ag+-treatment in S. miltiorrhiza hairy roots. Moreover, tissue-specifically expressed splice variants were identified for SmKO, SmGA20ox3, SmGA2ox3 and SmGA2ox11, of which SmKOv1, SmGA20ox3v and SmGA2ox11v1 were GA3-responsive, suggesting the importance of alternative splicing in regulating GA metabolism. Conclusions The results show tissue-specifically expressed, feedback-regulated, stress-responsive and alternatively spliced novel genes and reveal multiple layer regulation of GA metabolism and crosstalk between gibberellin metabolism and tanshinone biosynthesis in S. miltiorrhiza. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2315-5) contains supplementary material, which is available to authorized users.

Plants develop important regulatory networks to adapt to the frequently-changing availability of copper (Cu). However, little is known about miRNA-associated Cu-regulatory networks in plant species other than Arabidopsis. Here, we report that Cu-responsive miRNAs in Populus trichocarpa (Torr. & Gray) include not only conserved miR397, miR398 and miR408, but also Populus-specific miR1444, suggesting the conservation and diversity of Cu-responsive miRNAs in plants. Copper-associated suppression of mature miRNAs is in company with the up-regulation of their target genes encoding Cu-containing proteins in Populus. The targets include miR397-targeted PtLAC5, PtLAC6 and PtLAC110a, miR398-targeted PtCSD1, PtCSD2a and PtCSD2b, miR408-targeted PtPCL1, PtPCL2, PtPCL3 and PtLAC4, and miR1444-targeted PtPPO3 and PtPPO6. Consistently, P. trichocarpa miR408 promoter-directed GUS gene expression is down-regulated by Cu in transgenic tobacco plants. Cu-response elements (CuREs) are found in the promoters of Cu-responsive miRNA genes. We identified 34 SQUAMOSA-promoter binding protein-like (SPL) genes, of which 17 are full-length PtSPL proteins or partial sequences with at least 300 amino acids. Phylogenetic analysis indicates that PtSPL3 and PtSPL4 are CuRE-binding proteins controlling Cu-responsive gene expression. Cu appears to be not involved in the regulation of these transcription factors because neither PtSPL3 nor PtSPL4 is Cu-regulated and no CuRE exists in their promoters.

MicroRNAs (miRNAs), a group of small non-coding RNAs, have recently become the subject of intense study. They are a class of post-transcriptional negative regulators playing vital roles in plant development and growth. However, little is known about their regulatory roles in the responses of trees to the stressful environments incurred over their long-term growth. Here, we report the cloning of small RNAs from abiotic stressed tissues of Populus trichocarpa (Ptc) and the identification of 68 putative miRNA sequences that can be classified into 27 families based on sequence homology. Among them, nine families are novel, increasing the number of the known Ptc-miRNA families from 33 to 42. A total of 346 targets was predicted for the cloned Ptc-miRNAs using penalty scores of

Trees constitute the majority of lignocellulosic biomass existing on our planet. Trees also serve as important feedstock materials for various industrial products. However, little is known about the regulatory mechanisms of cellulose synthase (CesA) genes of trees. Here, the cloning and characterization of three CesA genes (EgraCesA1, EgraCesA2, and EgraCesA3) from an economically important tree species, Eucalyptus grandis, are reported. All three genes were specifically expressed in xylem cells of eucalyptus undergoing secondary cell wall biosynthesis. The GUS gene, expressed under the control of the EgraCesA2 or EgraCesA3 promoter, was also localized in the secondary xylem in transgenic tobacco stems. However, the EgraCesA1 promoter alone or along with its 5'-UTR introns was insufficient to direct appropriate GUS expression. EgraCesA2 and EgraCesA3 gene expression was up-regulated in tension-stressed eucalyptus xylem cells. Accordingly, GUS expression directed by the EgraCesA2 or EgraCesA3 promoter was also up-regulated. EgraCesA1 had no such response. Thus, it is most unlikely that EgraCesA1 is a subunit of the EgraCesA2-EgraCesA3 complex. The presence of at least two types of cellulose biosynthesis machinery in wood formation is an important clue in deciphering the underpinnings of the perennial growth of trees in various environmental conditions. By analysing GUS gene expression directed by the EgraCesA3 promoter or its deletions, several negative and positive regulatory regions controlling gene expression in xylem or phloem were identified. Also a region which is likely to contain mechanical stress-responsive elements was deduced. These results will guide further studies on identifying cis-regulatory elements directing CesA gene transcription and wood formation regulatory networks.

MicroRNAs (miRNAs) are endogenous small RNAs that can have large-scale regulatory effects on development and on stress responses in plants. The endemic rust fungus Cronartium quercuum f. sp. fusiforme causes fusiform rust disease in pines, resulting in the development of spindle-shaped galls (cankers) on branches or stems. This disease is the most destructive disease of pines in the southern USA. To test whether miRNAs play roles in fusiform rust gall development, we cloned and identified 26 miRNAs from stem xylem of loblolly pine (Pinus taeda), which belong to four conserved and seven loblolly pine-specific miRNA families. Forty-three targets for nine of these 11 families were experimentally validated in vivo. Sequence analysis suggested that the target cleavage site may be determined not only by the miRNA sequence but also by the target sequence. Members of three loblolly pine-specific miRNA families target a large number of non-protein coding transcripts, and one of these families could also initiate secondary phased production from its target of a putative trans-acting short interfering RNA (ta-siRNA). Expression of 10 of these 11 miRNA families was significantly repressed in the galled stem. PCR-based transcript quantification showed complex expression patterns of these miRNAs and their targets in the galled tissues and in tissues surrounding the gall. We further predict 82 plant disease-related transcripts that may also response to miRNA regulation in pine. These results reveal a new genetic basis for host-pathogen interactions in the development of fusiform rust gall.

Background MicroRNAs (miRNAs) are small, non-coding RNAs that are important regulators of gene expression, and play major roles in plant development and their response to the environment. Root extracts from Panax notoginseng contain triterpene saponins as their principal bioactive constituent, and demonstrate medicinal properties. To investigate the novel and conserved miRNAs in P. notoginseng, three small RNA libraries constructed from 1-, 2-, and 3-year-old roots in which root saponin levels vary underwent high-throughput sequencing. Methods P. notoginseng roots, purified from 1-, 2-, and 3-year-old roots, were extracted for RNA, respectively. Three small libraries were constructed and subjected to next generation sequencing. Results Sequencing of the three libraries generated 67,217,124 clean reads from P. notoginseng roots. A total of 316 conserved miRNAs (belonging to 67 miRNA families and one unclassified family) and 52 novel miRNAs were identified. MIR156 and MIR166 were the largest miRNA families, while miR156i and miR156g showed the highest abundance of miRNA species. Potential miRNA target genes were predicted and annotated using Cluster of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes. Comparing these miRNAs between root samples revealed 33 that were differentially expressed between 2- and 1-year-old roots (8 increased, 25 decreased), 27 differentially expressed between 3- and 1-year-old roots (7 increased, 20 decreased), and 29 differentially expressed between 3- and 2-year-old roots (8 increased, 21 decreased). Two significantly differentially expressed miRNAs and four miRNAs predicted to target genes involved in the terpenoid backbone biosynthesis pathway were selected and validated by quantitative reverse transcription PCR. Furthermore, the expression patterns of these six miRNAs were analyzed in P. notoginseng roots, stems, and leaves at different developmental stages. Conclusions This study identified a large number of P. notoginseng miRNAs and their target genes, functional annotations, and gene expression patterns. It provides the first known miRNA profiles of the P. notoginseng root development cycle. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2010-6) contains supplementary material, which is available to authorized users.

Fungi have evolved powerful genomic and chemical defense systems to protect themselves against genetic destabilization and other organisms. However, the precise molecular basis involved in fungal defense remain largely unknown in Basidiomycetes. Here the complete genome sequence, as well as DNA methylation patterns and small RNA transcriptomes, was analyzed to provide a holistic overview of secondary metabolism and defense processes in the model medicinal fungus, Ganoderma sinense. We reported the 48.96 Mb genome sequence of G. sinense, consisting of 12 chromosomes and encoding 15,688 genes. More than thirty gene clusters involved in the biosynthesis of secondary metabolites, as well as a large array of genes responsible for their transport and regulation were highlighted. In addition, components of genome defense mechanisms, namely repeat-induced point mutation (RIP), DNA methylation and small RNA-mediated gene silencing, were revealed in G. sinense. Systematic bioinformatic investigation of the genome and methylome suggested that RIP and DNA methylation combinatorially maintain G. sinense genome stability by inactivating invasive genetic material and transposable elements. The elucidation of the G. sinense genome and epigenome provides an unparalleled opportunity to advance our understanding of secondary metabolism and fungal defense mechanisms.

DCL1, the core component for miRNA biogenesis, is itself regulated by miR162 in Arabidopsis. MiRNA-mediated feedback regulation of AtDCL1 is important to maintain the proper level of DCL1 transcripts. However, it is unknown whether the miRNA-mediated regulation of DCL1 is conserved among plants. We analyzed the SmDCL gene family in Salvia miltiorrhiza, an emerging model plant for Traditional Chinese Medicine (TCM) studies, using a comprehensive approach integrating genome-wide prediction, molecular cloning, gene expression profiling and posttranscriptional regulation analysis. A total of five SmDCLs were identified. Comparative analysis of SmDCLs and AtDCLs showed an apparent enlargement of SmDCL introns in S. miltiorrhiza. The absence of miR162 in S. miltiorrhiza and the loss of miR162 target site in SmDCL1 were unexpectedly found. Further analysis showed that the miR162 target site was not present in DCL1 from ancient plants and was gained during plant evolution. The gained miR162 target site might be lost in a few modern plants through nucleotide mutations. Our results provide evidence for the gain and loss of miR162 and its target sites in Dicer-like genes during evolution. The data is useful for understanding the evolution of miRNA-mediated feedback regulation of DCLs in plants.

Background WRKY proteins comprise a large family of transcription factors and play important regulatory roles in plant development and defense response. The WRKY gene family in Salvia miltiorrhiza has not been characterized. Results A total of 61 SmWRKYs were cloned from S. miltiorrhiza. Multiple sequence alignment showed that SmWRKYs could be classified into 3 groups and 8 subgroups. Sequence features, the WRKY domain and other motifs of SmWRKYs are largely conserved with Arabidopsis AtWRKYs. Each group of WRKY domains contains characteristic conserved sequences, and group-specific motifs might attribute to functional divergence of WRKYs. A total of 17 pairs of orthologous SmWRKY and AtWRKY genes and 21 pairs of paralogous SmWRKY genes were identified. Maximum likelihood analysis showed that SmWRKYs had undergone strong selective pressure for adaptive evolution. Functional divergence analysis suggested that the SmWRKY subgroup genes and many paralogous SmWRKY gene pairs were divergent in functions. Various critical amino acids contributed to functional divergence among subgroups were detected. Of the 61 SmWRKYs, 22, 13, 4 and 1 were predominantly expressed in roots, stems, leaves, and flowers, respectively. The other 21 were mainly expressed in at least two tissues analyzed. In S. miltiorrhiza roots treated with MeJA, significant changes of gene expression were observed for 49 SmWRKYs, of which 26 were up-regulated, 18 were down-regulated, while the other 5 were either up-regulated or down-regulated at different time-points of treatment. Analysis of published RNA-seq data showed that 42 of the 61 identified SmWRKYs were yeast extract and Ag+-responsive. Through a systematic analysis, SmWRKYs potentially involved in tanshinone biosynthesis were predicted. Conclusion These results provide insights into functional conservation and diversification of SmWRKYs and are useful information for further elucidating SmWRKY functions. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1411-x) contains supplementary material, which is available to authorized users.