De Novo Transcriptome Assembly of Cedar (Cedrela odorata L.) and Differential Gene Expression Involved in Herbivore Resistance.

Timber trees are targets of herbivorous attacks. The identification of genes associated with pest resistance can be accomplished through differential expression analysis using transcriptomes. We reported the de novo assembly of cedar (Cedrela odorata L.) transcriptome and the differential expression of genes involved in herbivore resistance. The assembly and annotation of the transcriptome were obtained using RNAseq from healthy cedar plants and those infested with Chrysobothris yucatanensis. A total of 325.6 million reads were obtained, and 127,031 (97.47%) sequences were successfully assembled. A total of 220 herbivory-related genes were detected, of which 170 genes were annotated using GO terms, and 161 genes with 245 functions were identified-165, 75, and 5 were molecular functions, biological processes, and cellular components, respectively. To protect against herbivorous infestation, trees produce toxins and volatile compounds which are modulated by signaling pathways and gene expression related to molecular functions and biological processes. The limited number of genes identified as cellular components suggests that there are minimal alterations in cellular structure in response to borer attack. The chitin recognition protein, jasmonate ZIM-domain (JAZ) motifs, and response regulator receiver domain were found to be overexpressed, whereas the terpene synthase, cytochrome P450, and protein kinase domain gene families were underexpressed. This is the first report of a cedar transcriptome focusing on genes that are overexpressed in healthy plants and underexpressed in infested plants. This method may be a viable option for identifying genes associated with herbivore resistance.

Functional characterization of a bark-specific monoterpene synthase potentially involved in wounding- and methyl jasmonate-induced linalool emission in rubber (Hevea brasiliensis).

Rubber (Hevea brasiliensis) is a latex-producing plant that often encounters mechanical wounding, as well as pathogen and pest attacks through wound sites during and after tapping. Terpenoids play an important role in the ecological interactions of many plant species, and their diversity is mainly generated by enzymes known as terpene synthases (TPS). In this study, one cDNA sequence encoding a putative terpene synthase, HbTPS20, was obtained from the bark tissues of H. brasiliensis. The encoded protein contains 610 amino acids with a putative N-terminal plastid transit peptide of approximately 70 residues. It belongs to the TPS-b subfamily. Further phylogenetic analysis showed that HbTPS20 formed a separate branch that diverged from the progenitor of all other potentially functional terpene synthases of the rubber TPS-b subfamily. The truncated HbTPS20 without the signal peptide coding sequence was successfully expressed in E. coli and in vitro enzymatic assays with geranyl diphosphate (GPP) or neryl diphosphate (NPP) as a substrate defined HbTPS20 as an active linalool synthase (HbLIS) with the ability to produce linalool as the principal product. RT-qPCR analysis showed that the highest transcript levels of HbTPS20 were found in barks, and this gene was expressed at 2.26- and 250-fold greater levels in the bark tissues of wounded and MeJA-treated plants, respectively, than in those of the control plants. This indicates that this gene may be involved in the induced stress responses of rubber.

In Silico Genome-Wide Mining and Analysis of Terpene Synthase Gene Family in Hevea Brasiliensis.

Terpene synthases (TPSs) catalyze terpenoid synthesis and affect the intracellular isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) concentration. In this study, we mined the in silico genome-wide TPS genes of Hevea brasiliensis and identified 47 full-length TPS genes. They had DDXXD, DXDD, NSE/DTE, RR(X)8 W, EA(X)W, and other conserved motifs. The phylogenetic tree analysis revealed that the TPSs of H.brasiliensis (HbTPSs) were divided into five subfamilies, TPS-a, TPS-b, TPS-c, TPS-e/f, and TPS-g. HbTPSs were predicted to have functions in the cellular components, molecular functions, and biological processes. HbTPSs were involved in seven pathways, which were K14173, K14175, K15803, K04120, K04121, K17982, and K12742 in the secondary metabolite pathway prediction. Three-dimensional structures of HbTPSs of 7 pathways were predicted, and DDXXD, NSE/DTE, and EA(X)W conserved motifs near the binding sites were found. Cis-acting elements analysis showed that they had more cis-acting elements related to phytohormone responsiveness, which indicated that terpenoid biosynthesis might be related to phytohormone regulation. RNA-Seq analysis showed that different HbTPSs were expressed differentially in different tissues. This study's results help reveal the role of HbTPSs and their molecular mechanism and help resolve the regulatory mechanism of terpenoid biosynthesis in H.brasiliensis.

A Bioinformatics Tool for Efficient Retrieval of High-Confidence Terpene Synthases (TPS) and Application to the Identification of TPS in Coffea and Quillaja.

Terpenoids are a class of compounds that are found in all living organisms. In plants, some terpenoids are part of primary metabolism, but most terpenes found in plants are classified as specialized metabolites, encoded by terpene synthases (TPS). It is not obvious how to assign the putative product of a given TPS using bioinformatics tools. Phylogenetic analyses easily assign TPS into families; however members of the same TPS family can synthetize more than one terpenoid-and, in many biotechnological applications, researchers are more interested in the product of a given TPS rather than its phylogenetic profile. Automated protein annotation can be used to classify TPS based on their products, despite the family they belong to. Here, we implement an automated bioinformatics method, search_TPS, to identify TPS proteins that synthesize mono, sesqui and diterpenes in Angiosperms. We verified the applicability of the method by classifying wet lab validated TPS and applying it to find TPS proteins in Coffea arabica, C. canephora, C. eugenioides, and Quillaja saponaria. Search_TPS is a computational tool based on PERL scripts that carries out a series of HMMER searches against a curated database of TPS profile hidden Markov models. The tool is freely available at https://github.com/liliane-sntn/TPS .

Tropane alkaloids and terpenes synthase genes of Datura stramonium (Solanaceae).

BACKGROUND: Plants have evolved physical-chemical defense to prevent/diminish damage by their enemies. Chemical defense involves the synthesis' pathways of specialized toxic, repellent, or anti-nutritive metabolites to herbivores. Molecular evolutionary studies have revealed the origin of new genes, acquisition and functional diversification along time in different plant lineages. METHODS: Using bioinformatic tools we analyze gene divergence of tropane alkaloids (TAs) and terpene synthases (TPSs) in Datura stramonium and other species of Solanaceae; compared gene and amino acids sequence of TAs and TPSs on genomes, cDNA and proteins sequences of Viridiplantae. We analyzed two recently assembled genomes of D. stramonium (Ticumán and Teotihuacán), transcriptomes of Datura metel and genomes of other Solanaceae. Hence, we analyzed variation of TAs and TPSs to infer genes involved in plant defense and plant responses before stress. We analyzed protein modeling and molecular docking to predict interactions between H6H and ligand; we translated the sequences (Teo19488, Tic8550 and Tic8549) obtained from the two genomes of D. stramonium by using Swiss-Model and Ramachandran plot and MolProbity structure validation of Teo19488 protein model. RESULTS: For TAs, we detected an expansion event in the tropinone reductase II (TRII) and the ratio synonymous/non-synonymous substitutions indicate positive selection. In contrast, a contraction event and negative selection was detected in tropinone reductase I (TRI). In Hy-oscyamine 6 b-hydroxylase (H6H), enzyme involved in the production of tropane alkaloids atropine and scopolamine, the synonymous/non-synonymous substitution ratio in its dominion indicates positive selection. For terpenes (TPS), we found 18 DsTPS in D. stramomiun and seven in D. metel; evolutionary analyses detected positive selection in TPS10.1 and TPS10.2 of D. stramonium and D. metel. Comparison of copies of TPSs in D. stramonium detected variation among them in the binding site. Duplication events and differentiation of TAs and TPSs of D. stramonium, as compared to other Solanaceae, suggest their possible involvement on adaptive evolution of defense to herbivores. Protein modeling and docking show that the three protein structures obtained of DsH6H from Teo19488, Tic-8550 and Tic8549 maintain the same interactions and the union site of 2OG-FeII_Oxy with the Hy-o ligand as in 6TTM of D. metel. CONCLUSION: Our results indicate differences in the number of gene copies involved in the synthesis of tropane alkaloids, between the genomes of D. stramonium from two Mexican populations. More copies of genes related to the synthesis of tropane alkaloids (TRI, TRII, H6H, PMT) are found in D. stramonium as compared to Viridiplantae. Likewise, for terpene synthases (TPS), TPS-10 is duplicated in D. stramonium and D. metel. Further studies should be directed to experimentally assess gain (overexpression) or loss (silencing) of function of duplicated genes.

Corrigendum: Reference Genes for qPCR Analysis in Resin-Tapped Adult Slash Pine As a Tool to Address the Molecular Basis of Commercial Resinosis.

[This corrects the article on p. 849 in vol. 7, PMID: 27379135.].

Reference Genes for qPCR Analysis in Resin-Tapped Adult Slash Pine As a Tool to Address the Molecular Basis of Commercial Resinosis.

Pine oleoresin is a major source of terpenes, consisting of turpentine (mono- and sesquiterpenes) and rosin (diterpenes) fractions. Higher oleoresin yields are of economic interest, since oleoresin derivatives make up a valuable source of materials for chemical industries. Oleoresin can be extracted from living trees, often by the bark streak method, in which bark removal is done periodically, followed by application of stimulant paste containing sulfuric acid and other chemicals on the freshly wounded exposed surface. To better understand the molecular basis of chemically-stimulated and wound induced oleoresin production, we evaluated the stability of 11 putative reference genes for the purpose of normalization in studying Pinus elliottii gene expression during oleoresinosis. Samples for RNA extraction were collected from field-grown adult trees under tapping operations using stimulant pastes with different compositions and at various time points after paste application. Statistical methods established by geNorm, NormFinder, and BestKeeper softwares were consistent in pointing as adequate reference genes HISTO3 and UBI. To confirm expression stability of the candidate reference genes, expression profiles of putative P. elliottii orthologs of resin biosynthesis-related genes encoding Pinus contorta ß-pinene synthase [PcTPS-(-)ß-pin1], P. contorta levopimaradiene/abietadiene synthase (PcLAS1), Pinus taeda α-pinene synthase [PtTPS-(+)αpin], and P. taeda α-farnesene synthase (PtαFS) were examined following stimulant paste application. Increased oleoresin yields observed in stimulated treatments using phytohormone-based pastes were consistent with higher expression of pinene synthases. Overall, the expression of all genes examined matched the expected profiles of oleoresin-related transcript changes reported for previously examined conifers.

De novo assembly of Eugenia uniflora L. transcriptome and identification of genes from the terpenoid biosynthesis pathway.

Pitanga (Eugenia uniflora L.) is a member of the Myrtaceae family and is of particular interest due to its medicinal properties that are attributed to specialized metabolites with known biological activities. Among these molecules, terpenoids are the most abundant in essential oils that are found in the leaves and represent compounds with potential pharmacological benefits. The terpene diversity observed in Myrtaceae is determined by the activity of different members of the terpene synthase and oxidosqualene cyclase families. Therefore, the aim of this study was to perform a de novo assembly of transcripts from E. uniflora leaves and to annotation to identify the genes potentially involved in the terpenoid biosynthesis pathway and terpene diversity. In total, 72,742 unigenes with a mean length of 1048bp were identified. Of these, 43,631 and 36,289 were annotated with the NCBI non-redundant protein and Swiss-Prot databases, respectively. The gene ontology categorized the sequences into 53 functional groups. A metabolic pathway analysis with KEGG revealed 8,625 unigenes assigned to 141 metabolic pathways and 40 unigenes predicted to be associated with the biosynthesis of terpenoids. Furthermore, we identified four putative full-length terpene synthase genes involved in sesquiterpenes and monoterpenes biosynthesis, and three putative full-length oxidosqualene cyclase genes involved in the triterpenes biosynthesis. The expression of these genes was validated in different E. uniflora tissues.