KCNAB2 overexpression inhibits human non-small-cell lung cancer cell growth in vitro and in vivo.

Non-small-cell lung cancer (NSCLC) accounts for approximately 85% of all lung cancer cases. NSCLC patients often have poor prognosis demanding urgent identification of novel biomarkers and potential therapeutic targets. KCNAB2 (regulatory beta subunit2 of voltage-gated potassium channel), encoding aldosterone reductase, plays a pivotal role in regulating potassium channel activity. In this research, we tested the expression of KCNAB2 as well as its potential functions in human NSCLC. Bioinformatics analysis shows that expression of KCNAB2 mRNA is significantly downregulated in human NSCLC, correlating with poor overall survival. In addition, decreased KCNAB2 expression was detected in different NSCLC cell lines and local human NSCLC tissues. Exogenous overexpression of KCNAB2 potently suppressed growth, proliferation and motility of established human NSCLC cells and promoted NSCLC cells apoptosis. In contrast, CRISPR/Cas9-induced KCNAB2 knockout further promoted the malignant biological behaviors of NSCLC cells. Protein chip analysis in the KCNAB2-overexpressed NSCLC cells revealed that KCNAB2 plays a possible role in AKT-mTOR cascade activation. Indeed, AKT-mTOR signaling activation was potently inhibited following KCNAB2 overexpression in NSCLC cells. It was however augmented by KCNAB2 knockout. In vivo, the growth of subcutaneous KCNAB2-overexpressed A549 xenografts was significantly inhibited. Collectively, KCNAB2 could be a novel effective gene for prognosis prediction of NSCLC. Targeting KCNAB2 may lead to the development of advanced therapies.

Identification of matrix-remodeling associated 5 as a possible molecular oncotarget of pancreatic cancer.

Pancreatic cancer has an extremely poor prognosis. Here we examined expression, potential functions and underlying mechanisms of MXRA5 (matrix remodeling associated 5) in pancreatic cancer. Bioinformatics studies revealed that MXRA5 transcripts are significantly elevated in pancreatic cancer tissues, correlating with the poor overall survival, high T-stage, N1 and pathologic stage of the patients. MXRA5 mRNA and protein expression is significantly elevated in microarray pancreatic cancer tissues and different pancreatic cancer cells. In primary and immortalized (BxPC-3 and PANC-1 lines) pancreatic cancer cells, shRNA-induced MXRA5 silencing or CRISPR/Cas9-mediated MXRA5 knockout suppressed cell survival, proliferation, migration, invasion, and epithelial-mesenchymal transition (EMT), while provoking cell apoptosis. Conversely, forced overexpression of MXRA5 further promoted pancreatic cancer cell progression and EMT. Bioinformatics studies and the protein chip analyses revealed that differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) in MXRA5-overexpressed primary pancreatic cancer cells were enriched in the PI3K-Akt-mTOR cascade. Indeed, Akt-mTOR activation in primary human pancreatic cancer cells was inhibited by MXRA5 shRNA or knockout, but was augmented following MXRA5 overexpression. In vivo, the growth of MXRA5 KO PANC-1 xenografts was largely inhibited in nude mice. Moreover, intratumoral injection of adeno-associated virus-packed MXRA5 shRNA potently inhibited primary pancreatic cancer cell growth in nude mice. Akt-mTOR activation was also largely inhibited in the MXRA5-depleted pancreatic cancer xenografts. Contrarily MXRA5 overexpression promoted primary pancreatic cancer cell growth in nude mice. Together, overexpressed MXRA5 is important for pancreatic cancer cell growth possibly through promoting EMT and Akt-mTOR activation. MXRA5 could be a potential therapeutic oncotarget for pancreatic cancer.

bHLH010/089 Transcription Factors Control Pollen Wall Development via Specific Transcriptional and Metabolic Networks in Arabidopsis thaliana.

The pollen wall is a specialized extracellular cell wall that protects male gametophytes from various environmental stresses and facilitates pollination. Here, we reported that bHLH010 and bHLH089 together are required for the development of the pollen wall by regulating their specific downstream transcriptional and metabolic networks. Both the exine and intine structures of bhlh010 bhlh089 pollen grains were severely defective. Further untargeted metabolomic and transcriptomic analyses revealed that the accumulation of pollen wall morphogenesis-related metabolites, including polysaccharides, glyceryl derivatives, and flavonols, were significantly changed, and the expression of such metabolic enzyme-encoding genes and transporter-encoding genes related to pollen wall morphogenesis was downregulated in bhlh010 bhlh089 mutants. Among these downstream target genes, CSLB03 is a novel target with no biological function being reported yet. We found that bHLH010 interacted with the two E-box sequences at the promoter of CSLB03 and directly activated the expression of CSLB03. The cslb03 mutant alleles showed bhlh010 bhlh089-like pollen developmental defects, with most of the pollen grains exhibiting defective pollen wall structures.

Bioinformatics analysis and experimental verification of the prognostic and biological significance mediated by fatty acid metabolism related genes for hepatocellular carcinoma.

Background: Liver cancer is among the leading causes of death related to cancer around the world. The most frequent type of human liver cancer is hepatocellular carcinoma (HCC). Fatty acid (FA) metabolism is an emerging hallmark that plays a promoting role in numerous malignancies. This study aimed to discover a FA metabolism-related risk signature and formulate a better model for HCC patients' prognosis prediction. Methods: We collected mRNA expression data and clinical parameters of patients with HCC using the TCGA databases, and the differential FA metabolism-related genes were explored. To create a risk prognostic model, we carried out the consensus clustering as well as univariate and multivariate Cox regression analyses. 16 genes were used to establish a prognostic model, which was then validated in the ICGC dataset. The accuracy of the model was performed using receiver operating characteristic (ROC) analyses, decision curve analysis (DCA) and nomogram. The immune cell infiltration level of risk genes was evaluated with single-sample GSEA (ssGSEA) algorithm. To reflect the response to immunotherapy, immunophenoscore (IPS) was obtained from TCGA-LIHC. Then, the expression of the candidate risk genes (p < 0.05) was validated by qRT-PCR, Western blotting and single-cell transcriptomics. Cellular function assays were performed to revealed the biological function of HAVCR1. Results: According to the TCGA-LIHC cohort analysis, the majority of the FA metabolism-related genes were expressed differentially in the HCC and normal tissues. The prognosis of patients with high-risk scores was observed to be worse. Multivariate COX regression analysis confirmed that the model can be employed as an independent prognosis factor for HCC patients. Furthermore, ssGSEA analysis revealed a link between the model and the levels of immune cell infiltration. Our model scoring mechanism also provides a high predictive value in HCC patients receiving anti-PDL1 immunotherapy. One of the FA metabolism-related genes, HAVCR1, displays a significant differential expression between normal and HCC cell lines. Hepatocellular carcinoma cells (Huh7, and HepG2) proliferation, motility, and invasion were all remarkably inhibited by HAVCR1 siRNA. Conclusion: Our study identified a novel FA metabolism-related prognostic model, revealing a better potential treatment and prevention strategy for HCC.

HbTGA1, a TGA Transcription Factor From Hevea brasiliensis, Regulates the Expression of Multiple Natural Rubber Biosynthesis Genes.

The TGA transcription factors are known to modulate the biosynthesis of secondary metabolites in plants. However, their regulatory function in natural rubber (NR) biosynthesis was not revealed in the rubber tree (Hevea brasiliensis). Here, 14 genes encoding TGA transcription factors (name HbTGA1-HbTGA14) were identified in the rubber tree. HbTGAs were differentially expressed in different tissues. HbTGA1 was expressed at its highest level in latex. We found specific in vitro and in vivo binding of the HbTGA1 protein with promoters of multiple NR biosynthesis genes (HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2). The activation of the promoters of HbHMGS2 and HbCPT6 was significantly suppressed by HbTGA1, while the activities of promoters of HbHMGR2, HbCPT8, and HbSRPP2 were increased by HbTGA1. The promoter activities of HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2 were significantly increased by HbTGA1 under jasmonate stress, while the promoter activities of HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2 were also significantly increased by HbTGA1 under salicylic acid stress. The present study provides insights into the role of TGA transcription factors in regulating the expression of NR biosynthesis genes from H. brasiliensis.

Identification of phosphoenolpyruvate carboxykinase 1 as a potential therapeutic target for pancreatic cancer.

Pancreatic cancer is the third leading cause of cancer-related mortalities and is characterized by rapid disease progression. Identification of novel therapeutic targets for this devastating disease is important. Phosphoenolpyruvate carboxykinase 1 (PCK1) is the rate-limiting enzyme of gluconeogenesis. The current study tested the expression and potential functions of PCK1 in pancreatic cancer. We show that PCK1 mRNA and protein levels are significantly elevated in human pancreatic cancer tissues and cells. In established and primary pancreatic cancer cells, PCK1 silencing (by shRNA) or CRISPR/Cas9-induced PCK1 knockout potently inhibited cell growth, proliferation, migration and invasion, and induced robust apoptosis activation. Conversely, ectopic overexpression of PCK1 in pancreatic cancer cells accelerated cell proliferation and migration. RNA-seq analyzing of differentially expressed genes (DEGs) in PCK1-silenced pancreatic cancer cells implied that DEGs were enriched in the PI3K-Akt-mTOR cascade. In pancreatic cancer cells, Akt-mTOR activation was largely inhibited by PCK1 shRNA, but was augmented after ectopic PCK1 overexpression. In vivo, the growth of PCK1 shRNA-bearing PANC-1 xenografts was largely inhibited in nude mice. Akt-mTOR activation was suppressed in PCK1 shRNA-expressing PANC-1 xenograft tissues. Collectively, PCK1 is a potential therapeutic target for pancreatic cancer.

Genome-wide identification and expression analysis of terpene synthase gene family in Aquilaria sinensis.

Agarwood is the resinous portion of Aquilaria trees, and has been widely used as medicine and incense. Sesquiterpenes are the main chemical characteristic constituents of agarwood. Terpene synthase (TPS) is a critical enzyme responsible for biosynthesis of sesquiterpene compounds. However, limited information is available on genome-wide identification and characterization of the TPS family in Aquilaria trees. In this study, TPS gene family was identified and characterized in Aquilaria sinensis by bioinformatics methods. The expression of those genes was analyzed by RNA-seq and quantitative real-time PCR. Transcription factors regulating TPS gene expression were identified by yeast one-hybrid and dual-luciferase assay. In total, 26 AsTPS genes (AsTPS1-AsTPS26) were identified, which were classified into five subgroups. Many putative cis-elements putatively involved in stresses and phytohormones (especially jasmonic acid) were identified in the promoter regions of AsTPSs, suggesting that AsTPSs genes may be regulated by stresses and jasmonic acid. Expression analysis revealed seven TPS genes encoding sesquiterpene synthetases were induced by wounding and methyl jasmonic acid (MeJA), which may be related to sesquiterpene biosynthesis. By yeast one-hybrid screening, a ERF transcription factor AsERF1 was identified to interact with the AsTPS1 promoter. Subcellular localization analysis indicated AsERF1 was a nucleus-localized protein. Transient transfection of AsERF1 in leaves of Nicotiana benthamiana significantly enhanced the promoter activation of AsTPS1, suggesting AsERF1 may participate in sesquiterpene biosynthesis by regulating AsTPS1 expression. These data generated in this study provide a foundation for future studies on functional roles and regulation mechanisms of AsTPS in sesquiterpene biosynthesis and agarwood formation.

Transcriptomes analysis reveals novel insight into the molecular mechanisms of somatic embryogenesis in Hevea brasiliensis.

BACKGROUND: Somatic embryogenesis (SE) is a promising technology for plant vegetative propagation, which has an important role in tree breeding. Though rubber tree (Hevea brasiliensis Muell. Arg.) SE has been founded, few late SE-related genes have been identified and the molecular regulation mechanisms of late SE are still not well understood. RESULTS: In this study, the transcriptomes of embryogenic callus (EC), primary embryo (PE), cotyledonary embryo (CE), abnormal embryo (AE), mature cotyledonary embryo (MCE) and withered abnormal embryo (WAE) were analyzed. A total of 887,852,416 clean reads were generated, 85.92% of them were mapped to the rubber tree genome. The de novo assembly generated 36,937 unigenes. The differentially expressed genes (DEGs) were identified in the pairwise comparisons of CE vs. AE and MCE vs. WAE, respectively. The specific common DEGs were mainly involved in the phytohormones signaling pathway, biosynthesis of phenylpropanoid and starch and sucrose metabolism. Among them, hormone signal transduction related genes were significantly enriched, especially the auxin signaling factors (AUX-like1, GH3.1, SAUR32-like, IAA9-like, IAA14-like, IAA27-like, IAA28-like and ARF5-like). The transcription factors including WRKY40, WRKY70, MYBS3-like, MYB1R1-like, AIL6 and bHLH93-like were characterized as molecular markers for rubber tree late SE. CML13, CML36, CAM-7, SERK1 and LEAD-29-like were also related to rubber tree late SE. In addition, histone modification had crucial roles during rubber tree late SE. CONCLUSIONS: This study provides important information to elucidate the molecular regulation during rubber tree late SE.

Tobacco rattle virus-induced gene silencing in Hevea brasiliensis.

Virus-induced gene silencing (VIGS) is a powerful gene-silencing tool that has been intensively applied in plants. To data, the application of VIGS in rubber tree has not yet been reported. In this study, we described the efficient gene silencing in rubber tree by VIGS. The gene encoding Hevea brasiliensis phytoene desaturase (HbPDS) was identified in rubber tree genome. Small interfering RNAs from HbPDS and the silencing gene fragment were predicted and a length of 399 bp was selected to be tested. We showed that the tobacco rattle virus (TRV)-VIGS could induce effective HbPDS silencing in rubber tree. This study was the first to report VIGS in rubber tree. The present TRV-VIGS method could be used to perform reverse genetic approaches to identify unknown gene functions and might be further applied to produce gene silenced rubber tree plants, to advance functional gene of rubber tree.

Differential Expression of lncRNAs and miRNAs Between Self-Rooting Juvenile and Donor Clones Unveils Novel Insight Into the Molecular Regulation of Rubber Biosynthesis in Hevea brasiliensis.

The rubber tree (Hevea brasiliensis Muell. Arg.) is a tropical tree species that produce natural rubber. Self-rooted juvenile clones (SRJCs) are novel rubber tree planting materials developed through primary somatic embryogenesis. SRJCs have a higher rubber yield compared with donor clones (DCs). The molecular basis underlying increased rubber yield in SRJCs remains largely unknown. Here, the latex from SRJCs and DCs were collected for strand-specific and small RNA-seq methods. A total of 196 differentially expressed long noncoding RNAs (DELs), and 11 differentially expressed microRNAs were identified in latex between SRJCs and DCs. Targeted genes of DELs were markedly enriched for various biological pathways related to plant hormone signal transduction, photosynthesis, glutathione metabolism, and amino acids biosynthesis. DELs probably acted as cis-acting regulation was calculated, and these DELs relevant to potentially regulate rubber biosynthesis, reactive oxygen species metabolism, and epigenetic modification. Furthermore, the DELs acting as microRNA targets were studied. The interaction of microRNA and DELs might involve in the regulation of natural rubber biosynthesis.

HbWRKY27, a group IIe WRKY transcription factor, positively regulates HbFPS1 expression in Hevea brasiliensis.

Farnesyl pyrophosphate synthase (FPS) is a key enzyme that catalyzes the formation of farnesyl pyrophosphate, the main initiator for rubber chain initiation in Hevea brasiliensis Muell. Arg. The transcriptional regulatory mechanisms of the FPS gene still not well understood. Here, a WRKY transcription factor designated HbWRKY27 was obtained by screening the latex cDNA library applied the HbFPS1 promoter as bait. HbWRKY27 interacted with the HbFPS1 promoter was further identified by individual Y1H and EMSA assays. HbWRKY27 belongs to group IIe WRKY subfamily which contains a typical WRKY domain and C-X5-CX23-HXH motif. HbWRKY27 was localized to the nucleus. HbWRKY27 predominantly accumulated in latex. HbWRKY27 was up-regulated in latex by ethrel, salicylic acid, abscisic acid, and methyl jasmonate treatment. Transient expression of HbWRKY27 led to increasing the activity of the HbFPS1 promoter in tobacco plant, suggesting that HbWRKY27 positively regulates the HbFPS1 expression. Taken together, an upstream transcription factor of the key natural rubber biosynthesis gene HbFPS1 was identified and this study will provide novel transcriptional regulatory mechanisms of the FPS gene in Hevea brasiliensis.

Identification of histone methylation modifiers and their expression patterns during somatic embryogenesis in Hevea brasiliensis.

Histone methylation plays a crucial role in various biological processes, from heterochromatin formation to transcriptional regulation. Currently, no information is available regarding histone methylation modifiers in the important rubber-producing plant Hevea brasiliensis. Here, we identified 47 histone methyltransferase (HMT) genes and 25 histone demethylase (HDM) genes as possible members of the histone methylation modifiers in the rubber tree genome. According to the structural features of HMT and HDM, the HbHMTs were classified into two groups (HbPRMs and HbSDGs), the HbHDMs have two groups (HbLSDs and HbJMJs). Expression patterns were analyzed in five different tissues and at different phases of somatic embryogenesis. HbSDG10, 21, 25, 33, HbJMJ2, 18, 20 were with high expression at different phases of somatic embryogenesis. HbSDG10,14, 20, 21, 33 and HbPRMT4 were expressed highly in anther, HbSDG14, 20, 21, 22, 23, 33, 35 and HbPRMT1 HbJMJ7 and HbLSD1, 2, 3, 4 showed high expression levels in callus. HbSDG1, 7, 10, 13, 14, 18, 19, 21, 22, 23, 35, HbPRMT1, 8, HbJMJ5, 7, 11, 16, 20 and HbLSD2, 3, 4 were expressed highly in somatic embryo. HbSDG10, 21, 25, 33, HbLSD2, 3 were expressed highly in bud of regenerated plant. The analyses reveal that HbHMTs and HbHDMs exhibit different expression patterns at different phases during somatic embryogenesis, implying that some HbHMTs and HbHDMs play important roles during somatic embryogenesis. This study provide fundamental information for further studies on histone methylation in Hevea brasiliensis.

Identification and characterization of the MADS-box genes highly expressed in the laticifer cells of Hevea brasiliensis.

MADS-box transcription factors possess many functions in plant reproduction and development. However, few MADS-box genes related to secondary metabolites regulation have been identified. In Hevea brasiliensis, natural rubber is a representative cis-polyisoprenoids in secondary metabolism which occurs in the rubber laticifer cells, the molecular regulation basis of natural rubber biosynthesis is not clear. Here, a total of 24 MADS-box genes including 4 type I MADS-box genes and 20 type II MADS-box genes were identified in the transcriptome of rubber tree latex. The phylogenetic analysis was performed to clarify the evolutionary relationships of all the 24 rubber tree MADS-box proteins with MADS-box transcription factors from Arabidopsis thaliana and Oryza sativa. Four type I MADS-box genes were subdivided into Mα (3 genes) and Mß (1 gene). Twenty type II MADS-box genes were subclassified into MIKC* (8 genes) and MIKCc (12 genes). Eight MADS-box genes (HblMADS3, 5, 6, 7, 9, 13, 23, 24) were predominant expression in laticifers. ABA up-regulated the expression of HblMADS9, and the expression of HblMADS3, HblMADS5, HblMADS24 were up-regulated by MeJA. The function of HblMADS24 was elucidated. HblMADS24 bound HbFPS1 promoter in yeast and HblMADS24 activated HbFPS1 promoter in tobacco plants. Moreover, we proposed that HblMADS24 is a transcription activator of HbFPS1 which taking part in natural rubber biosynthesis.

The 14-3-3 protein HbGF14a interacts with a RING zinc finger protein to regulate expression of the rubber transferase gene in Hevea brasiliensis.

Hevea brasiliensis is a key commercial source of natural rubber (cis 1,4-polyisoprene). In H. brasiliensis, rubber transferase is responsible for cis-1,4-polymerization of isoprene units from isopentenyl diphosphate and thus affects the yield of rubber. Little is known about the regulatory mechanisms of the rubber transferase gene at a molecular level. In this study we show that the 5'UTR intron of the promoter of the rubber transferase gene (HRT2) suppresses the expression of HRT2. A H. brasiliensis RING zinc finger protein (designated as HbRZFP1) was able to interact specifically with the HRT2 promoter to down-regulate its transcription in vivo. A 14-3-3 protein (named as HbGF14a) was identified as interacting with HbRZFP1, both in yeast and in planta. Transient co-expression of HbGF14a and HbRZFP1-encoding cDNAs resulted in HbRZFP1-mediated HRT2 transcription inhibition being relieved. HbGF14a repressed the protein-DNA binding of HbRZFP1 with the HRT2 promoter in yeast. We propose a regulatory mechanism by which the binding of HbGF14a to HbRZFP1 interferes with the interaction of HbRZFP1 with the HRT2 promoter, thereby repressing the protein-DNA binding between them. This study provides new insights into the role of HbGF14a in mediating expression of the rubber transferase gene in Hevea brasiliensis.

Identification, characterization and expression analysis of genes involved in steroidal saponin biosynthesis in Dracaena cambodiana.

Dracaena cambodiana is a traditional medicinal plant used for producing dragon's blood. The plants and dragon's blood of D. cambodiana contain a rich variety of steroidal saponins. However, little is known about steroidal saponin biosynthesis and its regulation in D. cambodiana. Here, 122 genes encoding enzymes involved in steroidal saponin biosynthesis were identified based on transcriptome data, with 29 of them containing complete open reading frames (ORF). Transcript expression analysis revealed that several genes related to steroidal saponin biosynthesis showed distinct tissue-specific expression patterns; the expression levels of genes encoding the key enzymes involved in the biosynthesis and early modification of steroidal saponins were significantly down-regulated in the stems in response to the inducer of dragon's blood, exhibiting positive correlations with the content of steroidal saponins. These results provide insights on the steroidal saponins biosynthetic pathway and mechanisms underlying induced formation of dragon's blood in D. cambodiana.

Transcriptome-Wide Identification and Characterization of MYB Transcription Factor Genes in the Laticifer Cells of Hevea brasiliensis.

MYB transcription factors hold vital roles in the regulation of plant secondary metabolic pathways. Laticifers in rubber trees (Hevea brasiliensis) are of primary importance in natural rubber production because natural rubber is formed and stored within these structures. To understand the role of MYB transcription factors in the specialized cells, we identified 44 MYB genes (named HblMYB1 to HblMYB44) by using our previously obtained transcriptome database of rubber tree laticifer cells and the public rubber tree genome database. Expression profiles showed that five MYB genes were highly expressed in the laticifers. HblMYB19 and HblMYB44 were selected for further study. HblMYB19 and HblMYB44 bound the promoters of HbFDPS1, HbSRPP, and HRT1 in yeast. Furthermore, the transient overexpression of HblMYB19 and HblMYB44 in tobacco plants significantly increased the activity of the promoters of HbFDPS1, HbSRPP, and HRT1. Basing on this information, we proposed that HblMYB19 and HblMYB44 are the regulators of HbFDPS1, HbSRPP, and HRT1, which are involved in the biosynthesis pathway of natural rubber.

Identification and characterization of the abscisic acid (ABA) receptor gene family and its expression in response to hormones in the rubber tree.

Abscisic acid (ABA) is an essential phytohormone involved in diverse physiological processes. Although genome-wide analyses of the ABA receptor PYR/PYL/RCAR (PYL) protein/gene family have been performed in certain plant species, little is known about the ABA receptor protein/gene family in the rubber tree (Hevea brasiliensis). In this study, we identified 14 ABA receptor PYL proteins/genes (designated HbPYL1 through HbPYL14) in the most recent rubber tree genome. A phylogenetic tree was constructed, which demonstrated that HbPYLs can be divided into three subfamilies that correlate well with the corresponding Arabidopsis subfamilies. Eight HbPYLs are highly expressed in laticifers. Five of the eight genes are simultaneously regulated by ABA, jasmonic acid (JA) and ethylene (ET). The identification and characterization of HbPYLs should enable us to further understand the role of ABA signal in the rubber tree.

De Novo transcriptome characterization of Dracaena cambodiana and analysis of genes involved in flavonoid accumulation during formation of dragon's blood.

Dragon's blood is a red resin mainly extracted from Dracaena plants, and has been widely used as a traditional medicine in East and Southeast Asia. The major components of dragon's blood are flavonoids. Owing to a lack of Dracaena plants genomic information, the flavonoids biosynthesis and regulation in Dracaena plants remain unknown. In this study, three cDNA libraries were constructed from the stems of D. cambodiana after injecting the inducer. Approximately 266.57 million raw sequencing reads were de novo assembled into 198,204 unigenes, of which 34,873 unique sequences were annotated in public protein databases. Many candidate genes involved in flavonoid accumulation were identified. Differential expression analysis identified 20 genes involved in flavonoid biosynthesis, 27 unigenes involved in flavonoid modification and 68 genes involved in flavonoid transport that were up-regulated in the stems of D. cambodiana after injecting the inducer, consistent with the accumulation of flavonoids. Furthermore, we have revealed the differential expression of transcripts encoding for transcription factors (MYB, bHLH and WD40) involved in flavonoid metabolism. These de novo transcriptome data sets provide insights on pathways and molecular regulation of flavonoid biosynthesis and transport, and improve our understanding of molecular mechanisms of dragon's blood formation in D. cambodiana.

HbMADS4, a MADS-box Transcription Factor from Hevea brasiliensis, Negatively Regulates HbSRPP.

In plants MADS-box transcription factors (TFs) play important roles in growth and development. However, no plant MADS-box gene has been identified to have a function related to secondary metabolites regulation. Here, a MADS-box TF gene, designated as HbMADS4, was isolated from Hevea brasiliensis by the yeast one-hybrid experiment to screen the latex cDNA library using the promoter of the gene encoding H. brasiliensis small rubber particle protein (HbSRPP) as bait. HbMADS4 was 984-bp containing 633-bp open reading frame encoding a deduced protein of 230 amino acid residues with a typical conserved MADS-box motif at the N terminus. HbMADS4 was preferentially expressed in the latex, but little expression was detected in the leaves, flowers, and roots. Its expression was inducible by methyl jasmonate and ethylene. Furthermore, transient over-expression and over-expression of HbMADS4 in transgenic tobacco plants significantly suppressed the activity of the HbSRP promoter. Altogether, it is proposed that HbMADS4 is a negative regulator of HbSRPP which participates in the biosynthesis of natural rubber.

Comparative Transcriptome Analysis of Latex Reveals Molecular Mechanisms Underlying Increased Rubber Yield in Hevea brasiliensis Self-Rooting Juvenile Clones.

Rubber tree (Hevea brasiliensis) self-rooting juvenile clones (JCs) are promising planting materials for rubber production. In a comparative trial between self-rooting JCs and donor clones (DCs), self-rooting JCs exhibited better performance in rubber yield. To study the molecular mechanism associated with higher rubber yield in self-rooting JCs, we sequenced and comparatively analyzed the latex of rubber tree self-rooting JCs and DCs at the transcriptome level. Total raw reads of 34,632,012 and 35,913,020 bp were obtained from the library of self-rooting JCs and DCs, respectively, by using Illumina HiSeq 2000 sequencing technology. De novo assemblies yielded 54689 unigenes from the library of self-rooting JCs and DCs. Among 54689 genes, 1716 genes were identified as differentially expressed between self-rooting JCs and DCs via comparative transcript profiling. Functional analysis showed that the genes related to the mass of categories were differentially enriched between the two clones. Several genes involved in carbohydrate metabolism, hormone metabolism and reactive oxygen species scavenging were up-regulated in self-rooting JCs, suggesting that the self-rooting JCs provide sufficient molecular basis for the increased rubber yielding, especially in the aspects of improved latex metabolisms and latex flow. Some genes encoding epigenetic modification enzymes were also differentially expressed between self-rooting JCs and DCs. Epigenetic modifications may lead to gene differential expression between self-rooting JCs and DCs. These data will provide new cues to understand the molecular mechanism underlying the improved rubber yield of H. brasiliensis self-rooting clones.