Characterization of the polyphenol oxidase gene family reveals a novel microRNA involved in posttranscriptional regulation of PPOs in Salvia miltiorrhiza

Genome-Wide Identification and Functional Analysis of Salvia miltiorrhiza MicroRNAs Reveal the Negative Regulatory Role of Smi-miR159a in Phenolic Acid Biosynthesis

MicroRNAs (miRNAs) are a group of endogenous small non-coding RNAs in plants. They play critical functions in various biological processes during plant growth and development. Salvia miltiorrhiza is a well-known traditional Chinese medicinal plant with significant medicinal, economic, and academic values. In order to elucidate the role of miRNAs in S. miltiorrhiza, six small RNA libraries from mature roots, young roots, stems, mature leaves, young leaves and flowers of S. miltiorrhiza and one degradome library from mixed tissues were constructed. A total of 184 miRNA precursors, generating 137 known and 49 novel miRNAs, were genome-widely identified. The identified miRNAs were predicted to play diversified regulatory roles in plants through regulating 891 genes. qRT-PCR and 5' RLM-RACE assays validated the negative regulatory role of smi-miR159a in SmMYB62, SmMYB78, and SmMYB80. To elucidate the function of smi-miR159a in bioactive compound biosynthesis, smi-miR159a transgenic hairy roots were generated and analyzed. The results showed that overexpression of smi-miR159a caused a significant decrease in rosmarinic acid and salvianolic acid B contents. qRT-PCR analysis showed that the targets of smi-miR159a, including SmMYB62, SmMYB78, and SmMYB80, were significantly down-regulated, accompanied by the down-regulation of SmPAL1, SmC4H1, Sm4CL1, SmTAT1, SmTAT3, SmHPPR1, SmRAS, and SmCYP98A14 genes involved in phenolic acid biosynthesis. It suggests that smi-miR159a is a significant negative regulator of phenolic acid biosynthesis in S. miltiorrhiza.

The Smi-miR858a-SmMYB module regulates tanshinone and phenolic acid biosynthesis in Salvia miltiorrhiza

Tanshinones and phenolic acids are two major classes of bioactive compounds in Salvia miltiorrhiza. Revealing the regulatory mechanism of their biosynthesis is crucial for quality improvement of S. miltiorrhiza medicinal materials. Here we demonstrated that Smi-miR858a-Smi-miR858c, a miRNA family previously known to regulate flavonoid biosynthesis, also played critical regulatory roles in tanshinone and phenolic acid biosynthesis in S. miltiorrhiza. Overexpression of Smi-miR858a in S. miltiorrhiza plants caused significant growth retardation and tanshinone and phenolic acid reduction. Computational prediction and degradome and RNA-seq analyses revealed that Smi-miR858a could directly cleave the transcripts of SmMYB6, SmMYB97, SmMYB111, and SmMYB112. Yeast one-hybrid and transient transcriptional activity assays showed that Smi-miR858a-regulated SmMYBs, such as SmMYB6 and SmMYB112, could activate the expression of SmPAL1 and SmTAT1 involved in phenolic acid biosynthesis and SmCPS1 and SmKSL1 associated with tanshinone biosynthesis. In addition to directly activating the genes involved in bioactive compound biosynthesis pathways, SmMYB6, SmMYB97, and SmMYB112 could also activate SmAOC2, SmAOS4, and SmJMT2 involved in the biosynthesis of methyl jasmonate, a significant elicitor of plant secondary metabolism. The results suggest the existence of dual signaling pathways for the regulation of Smi-miR858a in bioactive compound biosynthesis in S. miltiorrhiza.

Structural and expression analysis of polyphenol oxidases potentially involved in globe artichoke (C. cardunculus var. scolymus L.) tissue browning

Globe artichoke capitula are susceptible to browning due to oxidation of phenols caused by the activity of polyphenol oxidases (PPOs), this reduces their suitability for fresh or processed uses. A genome-wide analysis of the globe artichoke PPO gene family was performed. Bioinformatics analyses identified eleven PPOs and their genomic and amino acidic features were annotated. Cis-acting element analysis identified a gene regulatory and functional profile associated to plant growth and development as well as stress response. For some PPOs, phylogenetic analyses revealed a structural and functional conservation with different Asteraceae PPOs, while the allelic variants of the eleven PPOs investigated across four globe artichoke varietal types identified several SNP/Indel variants, some of which having impact on gene translation. By RTqPCR were assessed the expression patterns of PPOs in plant tissues and in vitro calli characterized by different morphologies. Heterogeneous PPO expression profiles were observed and three of them (PPO6, 7 and 11) showed a significant increase of transcripts in capitula tissues after cutting. Analogously, the same three PPOs were significantly up-regulated in calli showing a brown phenotype due to oxidation of phenols. Our results lay the foundations for a future application of gene editing aimed at disabling the three PPOs putatively involved in capitula browning.

Recent Advances of Polyphenol Oxidases in Plants

Polyphenol oxidase (PPO) is present in most higher plants, but also in animals and fungi. PPO in plants had been summarized several years ago. However, recent advances in studies of PPO in plants are lacking. This review concludes new researches on PPO distribution, structure, molecular weights, optimal temperature, pH, and substrates. And, the transformation of PPO from latent to active state was also discussed. This state shift is a vital reason for elevating PPO activity, but the activation mechanism in plants has not been elucidated. PPO has an important role in plant stress resistance and physiological metabolism. However, the enzymatic browning reaction induced by PPO is a major problem in the production, processing, and storage of fruits and vegetables. Meanwhile, we summarized various new methods that had been invented to decrease enzymatic browning by inhibiting PPO activity. In addition, our manuscript included information on several important biological functions and the transcriptional regulation of PPO in plants. Furthermore, we also prospect some future research areas of PPO and hope they will be useful for future research in plants.

Genome-Wide Identification and Characterization of the PPO Gene Family in Cotton (Gossypium) and Their Expression Variations Responding to Verticillium Wilt Infection

Polyphenol oxidases (PPOs) are copper-binding metalloproteinases encoded by nuclear genes, ubiquitously existing in the plastids of microorganisms, plants, and animals. As one of the important defense enzymes, PPOs have been reported to participate in the resistant processes that respond to diseases and insect pests in multiple plant species. However, PPO gene identification and characterization in cotton and their expression patterns under Verticillium wilt (VW) treatment have not been clearly studied. In this study, 7, 8, 14, and 16 PPO genes were separately identified from Gossypium arboreum, G. raimondii, G. hirsutum, and G. barbadense, respectively, which were distributed within 23 chromosomes, though mainly gathered in chromosome 6. The phylogenetic tree manifested that all the PPOs from four cotton species and 14 other plants were divided into seven groups, and the analyses of the conserved motifs and nucleotide sequences showed highly similar characteristics of the gene structure and domains in the cotton PPO genes. The dramatically expressed differences were observed among the different organs at various stages of growth and development or under the diverse stresses referred to in the published RNA-seq data. Quantitative real-time PCR (qRT-PCR) experiments were also performed on the GhPPO genes in the roots, stems, and leaves of VW-resistant MBI8255 and VW-susceptible CCRI36 infected with Verticillium dahliae V991, proving the strong correlation between PPO activity and VW resistance. A comprehensive analysis conducted on cotton PPO genes contributes to the screening of the candidate genes for subsequent biological function studies, which is also of great significance for the in-depth understanding of the molecular genetic basis of cotton resistance to VW.

MicroRNAs in Medicinal Plants

Medicinal plant microRNAs (miRNAs) are an endogenous class of small RNA central to the posttranscriptional regulation of gene expression. Biosynthetic research has shown that the mature miRNAs in medicinal plants can be produced from either the standard messenger RNA splicing mechanism or the pre-ribosomal RNA splicing process. The medicinal plant miRNA function is separated into two levels: (1) the cross-kingdom level, which is the regulation of disease-related genes in animal cells by oral intake, and (2) the intra-kingdom level, which is the participation of metabolism, development, and stress adaptation in homologous or heterologous plants. Increasing research continues to enrich the biosynthesis and function of medicinal plant miRNAs. In this review, peer-reviewed papers on medicinal plant miRNAs published on the Web of Science were discussed, covering a total of 78 species. The feasibility of the emerging role of medicinal plant miRNAs in regulating animal gene function was critically evaluated. Staged progress in intra-kingdom miRNA research has only been found in a few medicinal plants, which may be mainly inhibited by their long growth cycle, high demand for growth environment, immature genetic transformation, and difficult RNA extraction. The present review clarifies the research significance, opportunities, and challenges of medicinal plant miRNAs in drug development and agricultural production. The discussion of the latest results furthers the understanding of medicinal plant miRNAs and helps the rational design of the corresponding miRNA/target genes functional modules.

Genome-wide investigation and expression profiling of polyphenol oxidase (PPO) family genes uncover likely functions in organ development and stress responses in Populus trichocarpa

Background

Trees such as Populus are planted extensively for reforestation and afforestation. However, their successful establishment greatly depends upon ambient environmental conditions and their relative resistance to abiotic and biotic stresses. Polyphenol oxidase (PPO) is a ubiquitous metalloproteinase in plants, which plays crucial roles in mediating plant resistance against biotic and abiotic stresses. Although the whole genome sequence of Populus trichocarpa has long been published, little is known about the PPO genes in Populus, especially those related to drought stress, mechanical damage, and insect feeding. Additionally, there is a paucity of information regarding hormonal responses at the whole genome level.

Results

A genome-wide analysis of the poplar PPO family was performed in the present study, and 18 PtrPPO genes were identified. Bioinformatics and qRT-PCR were then used to analyze the gene structure, phylogeny, chromosomal localization, gene replication, cis-elements, and expression patterns of PtrPPOs. Sequence analysis revealed that two-thirds of the PtrPPO genes lacked intronic sequences. Phylogenetic analysis showed that all PPO genes were categorized into 11 groups, and woody plants harbored many PPO genes. Eighteen PtrPPO genes were disproportionally localized on 19 chromosomes, and 3 pairs of segmented replication genes and 4 tandem repeat genomes were detected in poplars. Cis-acting element analysis identified numerous growth and developmental elements, secondary metabolism processes, and stress-related elements in the promoters of different PPO members. Furthermore, PtrPPO genes were expressed preferentially in the tissues and fruits of young plants. In addition, the expression of some PtrPPOs could be significantly induced by polyethylene glycol, abscisic acid, and methyl jasmonate, thereby revealing their potential role in regulating the stress response. Currently, we identified potential upstream TFs of PtrPPOs using bioinformatics.

Conclusions

Comprehensive analysis is helpful for selecting candidate PPO genes for follow-up studies on biological function, and progress in understanding the molecular genetic basis of stress resistance in forest trees might lead to the development of genetic resources.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08028-9.

Identification and Characterization of Verticillium nonalfalfae-Responsive MicroRNAs in the Roots of Resistant and Susceptible Hop Cultivars

MicroRNAs are 21- to 24-nucleotide-long, non-coding RNA molecules that regulate gene expression at the post-transcriptional level. They can modulate various biological processes, including plant response and resistance to fungal pathogens. Hops are grown for use in the brewing industry and, recently, also for the pharmaceutical industry. Severe Verticillium wilt caused by the phytopathogenic fungus Verticillium nonalfalfae, is the main factor in yield loss in many crops, including hops (Humulus lupulus L.). In our study, we identified 56 known and 43 novel miRNAs and their expression patterns in the roots of susceptible and resistant hop cultivars after inoculation with V. nonalfalfae. In response to inoculation with V. nonalfalfae, we found five known and two novel miRNAs that are differentially expressed in the susceptible cultivar and six known miRNAs in the resistant cultivar. Differentially expressed miRNAs target 49 transcripts involved in protein localization and pigment synthesis in the susceptible cultivar, whereas they are involved in transcription factor regulation and hormone signalling in the resistant cultivar. The results of our study suggest that the susceptible and resistant hop cultivars respond differently to V. nonalfalfae inoculation at the miRNA level and that miRNAs may contribute to the successful defence of the resistant cultivar.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

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.

Genome-Wide Investigation of the NF-X1 Gene Family in Populus trichocarpa Expression Profiles during Development and Stress

Poplar are planted extensively in reforestation and afforestation. However, their successful establishment largely depends on the environmental conditions of the newly established plantation and their resistance to abiotic as well as biotic stresses. NF-X1, a widespread transcription factor in plants, plays an irreplaceable role in plant growth, development, and stress tolerance. Although the whole genome sequence of Populus trichocarpa has been published for a long time, little is known about the NF-X1 genes in poplar, especially those related to drought stress, mechanical damage, insect feeding, and hormone response at the whole genome level. In this study, whole genome analysis of the poplar NF-X1 family was performed, and 4 PtrNF-X1 genes were identified. Then, bioinformatics analysis and qRT-PCR were applied to analyze the gene structure, phylogeny, chromosomal localization, gene replication, Cis-elements, and expression patterns of PtrNF-X1genes. Sequence analysis revealed that one-quarter of the PtrNF-X1 genes did not contain introns. Phylogenetic analysis revealed that all NF-X1 genes were split into three subfamilies. The number of two pairs of segmented replication genes were detected in poplars. Cis-acting element analysis identified a large number of elements of growth and development and stress-related elements on the promoters of different NF-X1 members. In addition, some PtrNF-X1 could be significantly induced by polyethylene glycol (PEG) and abscisic acid (ABA), thus revealing their potential role in regulating stress response. Comprehensive analysis is helpful in selecting candidate NF-X1 genes for the follow-up study of the biological function, and molecular genetic progress of stress resistance in forest trees provides genetic resources.

Genome-wide identification and functional characterization of natural antisense transcripts in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most widely used traditional medicines. Natural antisense transcripts (NATs) are a class of long noncoding RNAs that can regulate gene expression. Here, we identified 812 NATs, including 168 cis-NATs and 644 trans-NATs from twelve root, flower, and leaf samples of S. miltiorrhiza using RNA-seq. The expression profiles for 41 of 50 NATs and their sense transcripts (STs) obtained from RNA-Seq were validated using qRT-PCR. The expression profiles of 17 NATs positively correlated with their STs. GO and KEGG pathway analyses mapped the STs for cis-NATs to pathways for biosynthesis of secondary metabolites. We characterized four NATs in detail, including NAT0001, NAT0002, NAT0004, and NAT00023. Their STs are kaurene synthase-like 1 and the homologs of UDP-glucose flavonoid 3-O-glucosyltransferase 6, UDP-glycosyltransferase 90A1, and beta-glucosidase 40, respectively. The first gene is involved in the biosynthesis of bioactive tanshinones, the next two are involved in anthocyanin biosynthesis, whereas the last is involved in phenylpropanoid biosynthesis. Besides, we found seven STs that are potential targets of miRNAs. And we found two miRNAs including miR156a and miR7208, might originate from NATs, NAT0112 and NAT0086. The results suggest that S. miltiorrhiza NATs might interact with STs, produce miRNAs, and be regulated by miRNAs. They potentially play significant regulatory roles in the biosynthesis of bioactive compounds.

Recent advances in polyphenol oxidase-mediated plant stress responses

Plant polyphenol oxidases (PPOs) are ubiquitous copper metalloenzymes with a biochemistry that has been known for more than a century. By the 1990s, biologists began to recognize the importance of PPOs in plant response to the infestation of herbivores and pathogens; ideas concerning a defensive role for PPOs arose to address observed evidence, and several testable hypotheses were suggested. Two pivotal discoveries in tomato (Lycopersicon esculentum Miller) plants, an inverse correlation between PPO levels and insect growth and PPO induction by defence signals, have driven many studies of PPO defence functions in the context of abiotic and biotic stresses. During the past three decades, extensive molecular research in transgenic and non-transgenic systems has partly revealed the sophisticated mechanisms underlying PPO defence against herbivores and pathogens. These understandings, rather than theoretical predictions, have driven the development of new hypotheses and advanced PPO-related studies. Here, we review progress in PPO family features, expression regulation and the defensive role of PPOs in plants. We propose assumptions of an extended range of co- and post-transcriptional processes to the regulation of unexplored PPO expression. In addition, the identification of endogenous PPO substrates and downstream targets of PPO action will be useful for elucidating PPO defensive roles. The potential effects of PPO-mediated oxidative defences on herbivore performance ultimately needs to be further investigated. Therefore, expanding multidisciplinary approaches to unexplored dimensions of PPO defence function should be a future priority.

Canavanine Increases the Content of Phenolic Compounds in Tomato (Solanum lycopersicum L.) Roots

Canavanine (CAN) is a nonproteinogenic amino acid, and its toxicity comes from its utilization instead of arginine in many cellular processes. As presented in previous experiments, supplementation of tomato (Solanum lycopersicum L.) with CAN led to decreased nitric oxide (NO) level and induced secondary oxidative stress. CAN improved total antioxidant capacity in roots, with parallel inhibition of enzymatic antioxidants. The aim of this work was to determine how CAN-dependent limitation of NO emission and reactive oxygen species overproduction impact content, localization, and metabolism of phenolic compounds (PCs) in tomato roots. Tomato seedlings were fed with CAN (10 and 50 µM) for 24 or 72 h. Inhibition of root growth due to CAN supplementation correlated with increased concentration of total PCs; CAN (50 µM) led to the homogeneous accumulation of PCs all over the roots. CAN increased also flavonoids content in root tips. The activity of polyphenol oxidases and phenylalanine ammonia-lyase increased only after prolonged treatment with 50 µM CAN, while expressions of genes encoding these enzymes were modified variously, irrespectively of CAN dosage and duration of the culture. PCs act as the important elements of the cellular antioxidant system under oxidative stress induced by CAN.

smi-miR396b targeted SmGRFs, SmHDT1, and SmMYB37/4 synergistically regulates cell growth and active ingredient accumulation in Salvia miltiorrhiza hairy roots

MIR396b had been cloned and overexpressed in Salvia miltiorrhiza hairy roots. MiR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to regulate cell growth and secondary metabolism in S. miltiorrhiza hairy roots. Danshen (Salvia miltiorrhiza Bunge) is a valuable medicinal herb with two kinds of clinically used natural products, salvianolic acids and tanshinones. miR396 is a conserved microRNA and plays extensive roles in plants. However, it is still unclear how miR396 works in S. miltiorrhiza. In this study, an smi-MIR396b has been cloned from S. miltiorrhiza. Overexpression of miR396b in danshen hairy roots inhibited hairy root growth, reduced salvianolic acid concentration, but enhanced tanshinone accumulation, resulting in the biomass and total salvianolic acids respectively reduced to 55.5 and 72.1% of the control and total tanshinones increased up to 1.91-fold of the control. Applied degradome sequencing, 5'RLM-RACE, and qRT-PCR, 13 targets for miR396b were identified including seven conserved SmGRF1-7 and six novel ones. Comparative transcriptomics and microRNomics analysis together with qRT-PCR results confirmed that miR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to mediate the phytohormone, especially gibberellin signaling pathways and consequentially resulted in the phenotype variation of miR396b-OE hairy roots. Furthermore, miR396b could be activated by methyl jasmonate, abscisic acid, gibberellin, salt, and drought stresses. The findings in this study indicated that smi-miR396b acts as an upstream and central regulator in cell growth and the biosynthesis of tanshinones and salvianolic acids, shedding light on the coordinated regulation of plant growth and biosynthesis of active ingredients in S. miltiorrhiza.

L-DOPA synthesis in Mucuna pruriens (L.) DC. is regulated by polyphenol oxidase and not CYP 450/tyrosine hydroxylase: An analysis of metabolic pathway using biochemical and molecular markers

Mucuna pruriens L., commonly known as velvetbean or cow-itch, is a self-pollinated tropical legume of the family Fabaceae, known for its medicinal properties. The active principle L-DOPA extracted from the plant is a potent drug used in the treatment of Parkinson's disease. Although, it is hypothesized that a single step reaction can produce L-DOPA, the presence of optional routes makes the pathway more intricate. For instance, the catecholamine biosynthetic pathway, which leads to L-DOPA production, could occur by hydroxylation of tyrosine to L-DOPA either by polyphenol oxidase (PPO) or tyrosine hydroxylase (TH). Furthermore, Cytochrome P450 (CYP) enzymes can also cause hydroxylation of tyrosine, resulting in L-DOPA synthesis. Therefore, the present investigation was focused on validating the step, which catalyzes the synthesis of L-DOPA, at the biochemical and molecular levels. Enzyme inhibitor studies showed significant inhibition of PPO enzyme with corresponding decrease in L-DOPA synthesis while TH and CYP inhibition had no effect on L-DOPA synthesis. Activity staining of non-denaturing PAGE gel for PPO and TH showed activity only to PPO enzyme. Following in-gel assay and tryptic digestion of the excised stained gel portion, peptide recovery and LC-MS/MS analysis were performed. Degenerate primers based on peptide sequence resulted in an 800bp amplicon. The subsequent sub-cloning, RACE analysis and BLAST search resulted in the isolation of full-length PPO coding sequence of 1800 bp. Structure prediction and phylogenetic analysis of the obtained sequence revealed strong similarity to other plant PPO's like Glycine max, Vigna radiata and Vicia faba of the same family.

Analysis of the laccase gene family and miR397-/miR408-mediated posttranscriptional regulation in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most commonly used traditional Chinese medicine materials. It contains important bioactive phenolic compounds, such as salvianolic acids, flavonoids and anthocyanins. Elucidation of phenolic compound biosynthesis and its regulatory mechanism is of great significance for S. miltiorrhiza quality improvement. Laccases (LACs) are multicopper-containing enzymes potentially involved in the polymerization of phenolic compounds. So far, little has been known about LAC genes in S. miltiorrhiza. Through systematic investigation of the whole genome sequence and transcriptomes of S. miltiorrhiza, we identified 65 full-length SmLAC genes (SmLAC1–SmLAC65). Phylogenetic analysis showed that 62 of the identified SmLACs clustered with LACs from Arabidopsis and Populus trichocarpa in seven clades (C1–C7), whereas the other three fell into one S. miltiorrhiza-specific clade (C8). All of the deduced SmLAC proteins contain four conserved signature sequences and three typical Cu-oxidase domains, and gene structures of most LACs from S. miltiorrhiza, Arabidopsis and P. trichocarpa were highly conserved, however SmLACs encoding C8 proteins showed distinct intron-exon structures. It suggests the conservation and diversity of plant LACs in gene structures. The majority of SmLACs exhibited tissue-specific expression patterns, indicates manifold functions of SmLACs played in S. miltiorrhiza. Analysis of high-throughput small RNA sequences and degradome data and experimental validation using the 5′ RACE method showed that 23 SmLACs were targets of Smi-miR397. Among them, three were also targeted by Smi-miR408. It suggests the significance of miR397 and miR408 in posttranscriptional regulation of SmLAC genes. Our results provide a foundation for further demonstrating the functions of SmLACs in the production of bioactive phenolic compounds in S. miltiorrhiza.

De novo origination of MIRNAs through generation of short inverted repeats in target genes

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.

Two New Polyphenol Oxidase Genes of Tea Plant (Camellia sinensis) Respond Differentially to the Regurgitant of Tea Geometrid, Ectropis obliqua

Polyphenol oxidases (PPOs) have been reported to play an important role in protecting plants from attacks by herbivores. Though PPO genes in other plants have been extensively studied, research on PPO genes in the tea plant (Camellia sinensis) is lacking. In particular, which members of the PPO gene family elicit the defense response of the tea plant are as yet unknown. Here, two new PPO genes, CsPPO1 and CsPPO2, both of which had high identity with PPOs from other plants, were obtained from tea leaves. The full length of CsPPO1 contained an open reading frame (ORF) of 1740 bp that encoded a protein of 579 amino acids, while CsPPO2 contained an ORF of 1788 bp that encoded a protein of 595 amino acids. The deduced CsPPO1 and CsPPO2 proteins had calculated molecular masses of 64.6 and 65.9 kDa; the isoelectric points were 6.94 and 6.48, respectively. The expression products of recombinant CsPPO1 and CsPPO2 in Escherichia coli were about 91 and 92 kDa, respectively, but the recombinant proteins existed in the form of an inclusion body. Whereas CsPPO1 is highly expressed in stems, CsPPO2 is highly expressed in roots. Further results showed that the expression of CsPPO1 and CsPPO2 was wound- and Ectropis obliqua-induced, and that regurgitant, unlike treatment with wounding plus deionized water, significantly upregulated the transcriptional expression of CsPPO2 but not of CsPPO1. The difference between regurgitant and wounding indicates that CsPPO2 may play a more meaningful defensive role against E. obliqua than CsPPO1. Meanwhile, we found the active component(s) of the regurgitant elicited by the expression of CsPPO may contain small molecules (under 3-kDa molecular weight). These conclusions advance the understanding of the biological function of two new PPO genes and show that one of these, CsPPO2, may be a promising gene for engineering tea plants that are resistant to E. obliqua.

Genomic Variation Among and Within Six Juglans Species

Genomic analysis in Juglans (walnuts) is expected to transform the breeding and agricultural production of both nuts and lumber. To that end, we report here the determination of reference sequences for six additional relatives of Juglans regia: Juglans sigillata (also from section Dioscaryon), Juglans nigra, Juglans microcarpa, Juglans hindsii (from section Rhysocaryon), Juglans cathayensis (from section Cardiocaryon), and the closely related Pterocarya stenoptera. While these are ‘draft’ genomes, ranging in size between 640Mbp and 990Mbp, their contiguities and accuracies can support powerful annotations of genomic variation that are often the foundation of new avenues of research and breeding. We annotated nucleotide divergence and synteny by creating complete pairwise alignments of each reference genome to the remaining six. In addition, we have re-sequenced a sample of accessions from four Juglans species (including regia). The variation discovered in these surveys comprises a critical resource for experimentation and breeding, as well as a solid complementary annotation. To demonstrate the potential of these resources the structural and sequence variation in and around the polyphenol oxidase loci, PPO1 and PPO2 were investigated. As reported for other seed crops variation in this gene is implicated in the domestication of walnuts. The apparently Juglandaceae specific PPO1 duplicate shows accelerated divergence and an excess of amino acid replacement on the lineage leading to accessions of the domesticated nut crop species, Juglans regia and sigillata.

The Pentatricopeptide Repeat Gene Family in Salvia miltiorrhiza: Genome-Wide Characterization and Expression Analysis

The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in plants and plays important roles in posttranscriptional regulation. In this study, we combined whole genome sequencing and transcriptomes to systematically investigate PPRs in Salvia miltiorrhiza, which is a well-known material of traditional Chinese medicine and an emerging model system for medicinal plant studies. Among 562 identified SmPPRs, 299 belong to the P subfamily while the others belong to the PLS subfamily. The majority of SmPPRs have only one exon and are localized in the mitochondrion or chloroplast. As many as 546 SmPPRs were expressed in at least one tissue and exhibited differential expression patterns, which indicates they likely play a variety of functions in S. miltiorrhiza. Up to 349 SmPPRs were salicylic acid-responsive and 183 SmPPRs were yeast extract and Ag⁺-responsive, which indicates these genes might be involved in S. miltiorrhiza defense stresses and secondary metabolism. Furthermore, 23 salicylic acid-responsive SmPPRs were co-expressed with phenolic acid biosynthetic enzyme genes only while 16 yeast extract and Ag⁺-responsive SmPPRs were co-expressed with tanshinone biosynthetic enzyme genes only. Two SmPPRs were co-expressed with both phenolic acid and tanshinone biosynthetic enzyme genes. The results provide a useful platform for further investigating the roles of PPRs in S. miltiorrhiza.

Characterization of Conserved and Novel microRNAs in Lilium lancifolium Thunb. by High-Throughput Sequencing

MicroRNAs (miRNAs) are among the class of noncoding small RNA molecules and play a crucial role in post-transcriptional regulation in plants. Although Lilium is one of the most popular ornamental flowers worldwide, however, there is no report on miRNAs identification. In the present study, therefore, miRNAs and their targets were identified from flower, leaf, bulblet and bulb of Lilium lancifolium Thunb. by high-throughput sequencing and bioinformatics analysis. In this study, a total of 38 conserved miRNAs belonging to 17 miRNA families and 44 novel miRNAs were identified. In total, 366 target genes for conserved miRNAs and 415 target genes for novel miRNAs were predicted. The majority of the target genes for conserved miRNAs were transcriptional factors and novel miRNAs targeted mainly protein coding genes. A total of 53 cleavage sites belonging to 6 conserved miRNAs families and 14 novel miRNAs were identified using degradome sequencing. Twenty-three miRNAs were randomly selected, then, their credibility was confirmed using northern blot or stem-loop qRT-PCR. The results from qRT-PCR analysis showed the expression pattern of 4 LL-miRNAs was opposite to their targets. Therefore, our finding provides an important basis to understand the biological functions of miRNAs in Lilium.