Reprogramming of DNA methylation and changes of gene expression in grafted Hevea brasiliensis.

Rubber tree (Hevea brasiliensis) is reproduced by bud grafting for commercial planting, but significant intraclonal variations exist in bud-grafted clones. DNA methylation changes related to grafting may be partly responsible for intraclonal variations. In the current study, whole-genome DNA methylation profiles of grafted rubber tree plants (GPs) and their donor plants (DPs) were evaluated by whole-genome bisulfite sequencing. Data showed that DNA methylation was downregulated and DNA methylations in CG, CHG, and CHH sequences were reprogrammed in GPs, suggesting that grafting induced the reprogramming of DNA methylation. A total of 5,939 differentially methylated genes (DMGs) were identified by comparing fractional methylation levels between GPs and DPs. Transcriptional analysis revealed that there were 9,798 differentially expressed genes (DEGs) in the DP and GP comparison. A total of 1,698 overlapping genes between DEGs and DMGs were identified. These overlapping genes were markedly enriched in the metabolic pathway and biosynthesis of secondary metabolites by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Global DNA methylation and transcriptional analyses revealed that reprogramming of DNA methylation is correlated with gene expression in grafted rubber trees. The study provides a whole-genome methylome of rubber trees and an insight into the molecular mechanisms underlying the intraclonal variations existing in the commercial planting of grafted rubber trees.

Molecular evolution and characterization of type III polyketide synthase gene family in Aquilaria sinensis.

2-(2-Phenylethyl) chromone (PEC) and its derivatives are markers of agarwood formation and are also related to agarwood quality. However, the biosynthetic and regulatory mechanisms of PECs still remain mysterious. Several studies suggested that type III polyketide synthases (PKSs) contribute to PEC biosynthesis in Aquilaria sinensis. Furthermore, systematic studies on the evolution of PKSs in A. sinensis have rarely been reported. Herein, we comprehensively analyzed PKS genes from 12 plant genomes and characterized the AsPKSs in detail. A unique branch contained only AsPKS members was identified through evolutionary analysis, including AsPKS01 that was previously indicated to participate in PEC biosynthesis. AsPKS07 and AsPKS08, two tandem-duplicated genes of AsPKS01 and lacking orthologous genes in evolutionary models, were selected for their transient expression in the leaves of Nicotiana benthamiana. Subsequently, PECs were detected in the extracts of N. benthamiana leaves, suggesting that AsPKS07 and AsPKS08 promote PEC biosynthesis. The interaction between the promoters of AsPKS07, AsPKS08 and five basic leucine zippers (bZIPs) from the S subfamily indicated that their transcripts could be regulated by these transcription factors (TFs) and might further contribute to PECs biosynthesis in A. sinensis. Our findings provide valuable insights into the molecular evolution of the PKS gene family in A. sinensis and serve as a foundation for advancing PEC production through the bioengineering of gene clusters. Ultimately, this contribution is expected to shed light on the mechanism underlying agarwood formation.

Revealing the role of CCoAOMT1: fine-tuning bHLH transcription factors for optimal anther development.

The tapetum, a crucial innermost layer encompassing male reproductive cells within the anther wall, plays a pivotal role in normal pollen development. The transcription factors (TFs) bHLH010/089/091 redundantly facilitate the rapid nuclear accumulation of DYSFUNCTIONAL TAPETUM 1, a gatekeeper TF in the tapetum. Nevertheless, the regulatory mechanisms governing the activity of bHLH010/089/091 remain unknown. In this study, we reveal that caffeoyl coenzyme A O-methyltransferase 1 (CCoAOMT1) is a negative regulator affecting the nuclear localization and function of bHLH010 and bHLH089, probably through their K259 site. Our findings underscore that CCoAOMT1 promotes the nuclear export and degradation of bHLH010 and bHLH089. Intriguingly, elevated CCoAOMT1 expression resulted in defective pollen development, mirroring the phenotype observed in bhlh010 bhlh089 mutants. Moreover, our investigation revealed that the K259A mutation in the bHLH089 protein disrupted its translocation from the nucleus to the cytosol and impeded its degradation induced by CCoAOMT1. Importantly, transgenic plants with the probHLH089::bHLH089K259A construct failed to rescue proper pollen development or gene expression in bhlh010 bhlh089 mutants. Collectively, these findings emphasize the need to maintain balanced TF homeostasis for male fertility. They firmly establish CCoAOMT1 as a pivotal regulator that is instrumental in achieving equilibrium between the induction of the tapetum transcriptional network and ensuring appropriate anther development.

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.

Systematic evolution of bZIP transcription factors in Malvales and functional exploration of AsbZIP14 and AsbZIP41 in Aquilaria sinensis.

Introduction: Agarwood, the dark-brown resin produced by Aquilaria trees, has been widely used as incense, spice, perfume or traditional medicine and 2-(2-phenethyl) chromones (PECs) are the key markers responsible for agarwood formation. But the biosynthesis and regulatory mechanism of PECs were still not illuminated. The transcription factor of basic leucine zipper (bZIP) presented the pivotal regulatory roles in various secondary metabolites biosynthesis in plants, which might also contribute to regulate PECs biosynthesis. However, molecular evolution and function of bZIP are rarely reported in Malvales plants, especially in Aquilaria trees. Methods and results: Here, 1,150 bZIPs were comprehensively identified from twelve Malvales and model species genomes and the evolutionary process were subsequently analyzed. Duplication types and collinearity indicated that bZIP is an ancient or conserved TF family and recent whole genome duplication drove its evolution. Interesting is that fewer bZIPs in A. sinensis than that species also experienced two genome duplication events in Malvales. 62 AsbZIPs were divided into 13 subfamilies and gene structures, conservative domains, motifs, cis-elements, and nearby genes of AsbZIPs were further characterized. Seven AsbZIPs in subfamily D were significantly regulated by ethylene and agarwood inducer. As the typical representation of subfamily D, AsbZIP14 and AsbZIP41 were localized in nuclear and potentially regulated PECs biosynthesis by activating or suppressing type III polyketide synthases (PKSs) genes expression via interaction with the AsPKS promoters. Discussion: Our results provide a basis for molecular evolution of bZIP gene family in Malvales and facilitate the understanding the potential functions of AsbZIP in regulating 2-(2-phenethyl) chromone biosynthesis and agarwood formation.

Genomic insight into domestication of rubber tree.

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

Systematic investigation of the R2R3-MYB gene family in Aquilaria sinensis reveals a transcriptional repressor AsMYB054 involved in 2-(2-phenylethyl)chromone biosynthesis.

Trees in the genus Aquilaria produce agarwood, a valuable resin used in medicine, perfumes, and incense. 2-(2-Phenethyl)chromones (PECs) are characteristic components of agarwood; however, molecular mechanisms underlying PEC biosynthesis and regulation remain largely unknown. The R2R3-MYB transcription factors play important regulatory roles in the biosynthesis of various secondary metabolites. In this study, 101 R2R3-MYB genes in Aquilaria sinensis were systematically identified and analyzed at the genome-wide level. Transcriptomic analysis revealed that 19 R2R3-MYB genes were significantly regulated by an agarwood inducer, and showed significant correlations with PEC accumulation. Expression and evolutionary analyses revealed that AsMYB054, a subgroup 4 R2R3-MYB, was negatively correlated with PEC accumulation. AsMYB054 was located in the nucleus and functioned as a transcriptional repressor. Moreover, AsMYB054 could bind to the promoters of the PEC biosynthesis related genes AsPKS02 and AsPKS09, and inhibit their transcriptional activity. These findings suggested that AsMYB054 functions as a negative regulator of PEC biosynthesis via the inhibition of AsPKS02 and AsPKS09 in A. sinensis. Our results provide a comprehensive understanding of the R2R3-MYB subfamily in A. sinensis and lay a foundation for further functional analyses of R2R3-MYB genes in PEC biosynthesis.

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.

Endothelial cell-specific molecule 1 drives cervical cancer progression.

The expression, biological functions and underlying molecular mechanisms of endothelial cell-specific molecule 1 (ESM1) in human cervical cancer remain unclear. Bioinformatics analysis revealed that ESM1 expression was significantly elevated in human cervical cancer tissues, correlating with patients' poor prognosis. Moreover, ESM1 mRNA and protein upregulation was detected in local cervical cancer tissues and various cervical cancer cells. In established and primary cervical cancer cells, ESM1 shRNA or CRISPR/Cas9-induced ESM1 KO hindered cell proliferation, cell cycle progression, in vitro cell migration and invasion, and induced significant apoptosis. Whereas ESM1 overexpression by a lentiviral construct accelerated proliferation and migration of cervical cancer cells. Further bioinformatics studies and RNA sequencing data discovered that ESM1-assocaited differentially expressed genes (DEGs) were enriched in PI3K-Akt and epithelial-mesenchymal transition (EMT) cascades. Indeed, PI3K-Akt cascade and expression of EMT-promoting proteins were decreased after ESM1 silencing in cervical cancer cells, but increased following ESM1 overexpression. Further studies demonstrated that SYT13 (synaptotagmin 13) could be a primary target gene of ESM1. SYT13 silencing potently inhibited ESM1-overexpression-induced PI3K-Akt cascade activation and cervical cancer cell migration/invasion. In vivo, ESM1 knockout hindered SiHa cervical cancer xenograft growth in mice. In ESM1-knockout xenografts tissues, PI3K-Akt inhibition, EMT-promoting proteins downregulation and apoptosis activation were detected. In conclusion, overexpressed ESM1 is important for cervical cancer growth in vitro and in vivo, possibly by promoting PI3K-Akt activation and EMT progression. ESM1 represents as a promising diagnostic marker and potential therapeutic target of cervical 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.

Identification and Characterization of Chalcone Isomerase Genes Involved in Flavonoid Production in Dracaena cambodiana.

Dragon's blood is a traditional medicine in which flavonoids are the main bioactive compounds; however, the underlying formation mechanism of dragon's blood remains largely poorly understood. Chalcone isomerase (CHI) is the key enzyme in the flavonoid biosynthesis pathway. However, CHI family genes are not well understood in Dracaena cambodiana Pierre ex Gagnep, an important source plant of dragon's blood. In this study, 11 CHI family genes were identified from D. cambodiana, and they were classified into three types. Evolutionary and transcriptional profiling analysis revealed that DcCHI1 and DcCHI4 might be involved in flavonoid production. Both DcCHI1 and DcCHI4 displayed low expression levels in stem under normal growth conditions and were induced by methyl jasmonate (MeJA), 6-benzyl aminopurine (6-BA, synthetic cytokinin), ultraviolet-B (UV-B), and wounding. The recombinant proteins DcCHI1 and DcCHI4 were expressed in Escherichia coli and purified by His-Bind resin chromatography. Enzyme activity assay indicated that DcCHI1 catalyzed the formation of naringenin from naringenin chalcone, while DcCHI4 lacked this catalytic activity. Overexpression of DcCHI1 or DcCHI4 enhanced the flavonoid production in D. cambodiana and tobacco. These findings implied that DcCHI1 and DcCHI4 play important roles in flavonoid production. Thus, our study will not only contribute to better understand the function and expression regulation of CHI family genes involved in flavonoid production in D. cambodiana but also lay the foundation for developing the effective inducer of dragon's blood.

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.

Genome sequence of the agarwood tree Aquilaria sinensis (Lour.) Spreng: the first chromosome-level draft genome in the Thymelaeceae family.

BACKGROUD: Aquilaria sinensis (Lour.) Spreng is one of the important plant resources involved in the production of agarwood in China. The agarwood resin collected from wounded Aquilaria trees has been used in Asia for aromatic or medicinal purposes from ancient times, although the mechanism underlying the formation of agarwood still remains poorly understood owing to a lack of accurate and high-quality genetic information. FINDINGS: We report the genomic architecture of A. sinensis by using an integrated strategy combining Nanopore, Illumina, and Hi-C sequencing. The final genome was ∼726.5 Mb in size, which reached a high level of continuity and a contig N50 of 1.1 Mb. We combined Hi-C data with the genome assembly to generate chromosome-level scaffolds. Eight super-scaffolds corresponding to the 8 chromosomes were assembled to a final size of 716.6 Mb, with a scaffold N50 of 88.78 Mb using 1,862 contigs. BUSCO evaluation reveals that the genome completeness reached 95.27%. The repeat sequences accounted for 59.13%, and 29,203 protein-coding genes were annotated in the genome. According to phylogenetic analysis using single-copy orthologous genes, we found that A. sinensis is closely related to Gossypium hirsutum and Theobroma cacao from the Malvales order, and A. sinensis diverged from their common ancestor ∼53.18-84.37 million years ago. CONCLUSIONS: Here, we present the first chromosome-level genome assembly and gene annotation of A. sinensis. This study should contribute to valuable genetic resources for further research on the agarwood formation mechanism, genome-assisted improvement, and conservation biology of Aquilaria species.