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.

Immunogenicity and safety of heterologous immunisation with Ad5-nCOV in healthy adults aged 60 years and older primed with an inactivated SARS-CoV-2 vaccine (CoronaVac): a phase 4, randomised, observer-blind, non-inferiority trial.

Background: People over 60 have been found to develop less protection after two doses of inactivated COVID-19 vaccines than younger people. Heterologous immunisation could potentially induce more robust immune responses compared to homologous immunisation. We aimed to assess the immunogenicity and safety of a heterologous immunisation with an adenovirus type 5-vectored vaccine (Ad5-nCOV, Convidecia) among elderly who were primed with an inactivated vaccine (CoronaVac) previously. Methods: We did a randomised, observer-blinded, non-inferiority trial in healthy adults aged 60 years and older in Lianshui County (Jiangsu, China) between August 26, 2021 and May 15, 2022. 199 eligible participants who had received two doses of CoronaVac in the past 3-6 months were randomised (1:1) to receive a third dose of Convidecia (group A, n = 99) or CoronaVac (group B, n = 100), while 100 participants primed with one dose of CoronaVac in the past 1-2 months were randomised equally to receive a second dose of Convidecia (group C, n = 50) or CoronaVac (group D, n = 50). Participants and investigators were masked to the vaccine received. Primary outcomes were the geometric mean titers (GMTs) of neutralising antibodies against live SARS-CoV-2 virus 14 days after boosting and 28-day adverse reactions. This study was registered with ClinicalTrials.govNCT04952727. Findings: A heterologous third dose of Convidecia resulted in a 6.2-fold (GMTs: 286.4 vs 48.2), 6.3-fold (45.9 vs 7.3) and 7.5-fold (32.9 vs 4.4) increase in neutralising antibodies against SARS-CoV-2 wild-type, delta (B.1.617.2) and omicron (BA.1.1) 14 days post boosting, respectively, compared with the homologous boost. The heterologous booster with Convidecia induced significantly higher neutralsing activities, with up to 91% inhibition in binding of Spike to ACE2 for BA.4 and BA.5 variants, compared with 35% inhibition induced by three doses of CoronaVac. For participants primed with one dose of CoronaVac, a heterologous dose of Convidecia induced higher neutralising antibodies against wild-type than two doses of CoronaVac (GMTs: 70.9 vs 9.3, p < 0.0001), but not for that against variants of concern (GMTs against delta: 5.0 vs 4.0, p = 0.4876; GMTs against omicron: 4.8 vs 3.7, p = 0.4707). Adverse reactions were reported by 8 (8.1%) participants in group A and 4 (4.0%) in group B (p ＞ 0.05), and 8 (16.0%) in group C and 1 (2.0%) in group D (p = 0.031). Interpretation: In elderly individuals primed with two doses of CoronaVac, the heterologous immunisation with Convidecia induced strong antibodies against SARS-CoV-2 wildtype and variants of concern, which could be an alternative regimen for enhancing protection in this vulnerable population. Funding: National Natural Science Foundation of China, Jiangsu Provincial Key Research and Development Program, and Jiangsu Science Fund for Distinguished Young Scholars Program.

Safety and immunogenicity of aerosolised Ad5-nCoV, intramuscular Ad5-nCoV, or inactivated COVID-19 vaccine CoronaVac given as the second booster following three doses of CoronaVac: a multicentre, open-label, phase 4, randomised trial.

BACKGROUND: Aerosolised Ad5-nCoV is the first approved mucosal respiratory COVID-19 vaccine to be used as a booster after the primary immunisation with COVID-19 vaccines. This study aimed to evaluate the safety and immunogenicity of aerosolised Ad5-nCoV, intramuscular Ad5-nCoV, or inactivated COVID-19 vaccine CoronaVac given as the second booster. METHODS: This is an open-label, parallel-controlled, phase 4 randomised trial enrolling healthy adult participants (≥18 years) who had completed a two-dose primary immunisation and a booster immunisation with inactivated COVID-19 vaccines (CoronaVac only) at least 6 months before, in Lianshui and Donghai counties, Jiangsu Province, China. We recruited eligible participants from previous trials in China (NCT04892459, NCT04952727, and NCT05043259) as cohort 1 (with the serum before and after the first booster dose available), and from eligible volunteers in Lianshui and Donghai counties, Jiangsu Province, as cohort 2. Participants were randomly assigned at a ratio of 1:1:1, using a web-based interactive response randomisation system, to receive the fourth dose (second booster) of aerosolised Ad5-nCoV (0·1 mL of 1·0 × 1011 viral particles per mL), intramuscular Ad5-nCoV (0·5 mL of 1·0 × 1011 viral particles per mL), or inactivated COVID-19 vaccine CoronaVac (0·5 mL), respectively. The co-primary outcomes were safety and immunogenicity of geometric mean titres (GMTs) of serum neutralising antibodies against prototype live SARS-CoV-2 virus 28 days after the vaccination, assessed on a per-protocol basis. Non-inferiority or superiority was achieved when the lower limit of the 95% CI of the GMT ratio (heterologous group vs homologous group) exceeded 0·67 or 1·0, respectively. This study was registered with ClinicalTrials.gov, NCT05303584 and is ongoing. FINDINGS: Between April 23 and May 23, 2022, from 367 volunteers screened for eligibility, 356 participants met eligibility criteria and received a dose of aerosolised Ad5-nCoV (n=117), intramuscular Ad5-nCoV (n=120), or CoronaVac (n=119). Within 28 days of booster vaccination, participants in the intramuscular Ad5-nCoV group reported a significantly higher frequency of adverse reactions than those in the aerosolised Ad5-nCoV and intramuscular CoronaVac groups (30% vs 9% and 14%, respectively; p<0·0001). No serious adverse events related to the vaccination were reported. The heterologous boosting with aerosolised Ad5-nCoV triggered a GMT of 672·4 (95% CI 539·7-837·7) and intramuscular Ad5-nCoV triggered a serum neutralising antibody GMT of 582·6 (505·0-672·2) 28 days after the booster dose, both of which were significantly higher than the GMT in the CoronaVac group (58·5 [48·0-71·4]; p<0·0001). INTERPRETATION: A heterologous fourth dose (second booster) with either aerosolised Ad5-nCoV or intramuscular Ad5-nCoV was safe and highly immunogenic in healthy adults who had been immunised with three doses of CoronaVac. FUNDING: National Natural Science Foundation of China, Jiangsu Provincial Science Fund for Distinguished Young Scholars, and Jiangsu Provincial Key Project of Science and Technology Plan.

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.

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.

Prolonged interval time between blastocyst biopsy and vitrification compromised the outcomes in preimplantation genetic testing.

This study aimed to evaluate to what extent the different interval times between trophectoderm (TE) biopsy and vitrification influence the clinical outcomes in preimplantation genetic testing (PGT) cycles. Patients who underwent frozen embryo transfer (FET) after PGT between 2015 and 2019 were recruited. In total, 297 cycles with single day 5 euploid blastocyst transfer were included. These cycles were divided into three groups according to the interval times: <1 h group, 1-2 h group, and ≥2 h group. Blastocyst survival, clinical pregnancy, miscarriage, and ongoing pregnancy rates were compared. The results showed that, in PGT-SR cycles, survival rate in the ≥2 h group (96.72%) was significantly lower than in the <1 h group (100%, P = 0.047). The clinical pregnancy rate in the ≥2 h group was 55.93%, significantly lower than in the <1 h group (74.26%, P = 0.017). The ongoing pregnancy rates in the 1-2 h group and the ≥2 h group were 48.28% and 47.46%, respectively, significantly lower than that in the <1 h group (67.33%, P < 0.05). The miscarriage rate in the 1-2 h group was 18.42%, significantly higher than that in the <1 h group (5.33%, P = 0.027). In PGT-A cycles, the clinical pregnancy and ongoing pregnancy rates in the <1 h group were 67.44% and 53.49%, respectively, higher than that in the 1-2 h group (52.94%, 47.06%, P > 0.05) and the ≥2 h group (52.63%, 36.84%, P > 0.05). In conclusion, vitrification of blastocysts beyond 1 h after biopsy significantly influences embryo survival and clinical outcomes and is therefore not recommended.

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.

Tritium diffusion in a Li2TiO3 crystal terminated with the (001) surface from first-principles calculations.

The tritium release behavior of the Li2TiO3 crystal has become an important index to evaluate its comprehensive performance as a solid breeder material in nuclear fusion reactors. The tritium diffusion on the surface (surface diffusion) and diffusion from the inside to the surface (hopping diffusion) in Li2TiO3 crystals with a 1/3-Li(001) surface are systematically investigated by the first-principles method. Possible adsorption sites, diffusion pathways and energy barriers of surface diffusion and hopping diffusion have been calculated and analyzed, respectively. Tritium atoms are found to diffuse preferentially along the [100] direction on the surface and two equivalent pathways across the surface were identified. The obtained activation energies are about 0.50 eV for surface diffusion and 1.56 eV for hopping diffusion. The local density of states and Bader charge for typical surface diffusion and hopping diffusion pathways are calculated and analyzed. The results reveal that the tritium (T) atom bonds with neighboring oxygen (O) atoms during the surface diffusion, while the T-O interaction is significantly weakened in the hopping diffusion which results in the higher activation energy than that of surface diffusion. In combination with our previous work, a complete tritium diffusion model for the Li2TiO3 crystal is proposed and the corresponding tritium diffusion coefficients are obtained. Our obtained activation energies are in the same range as previous experimental data and could provide theoretical support for the future related experiments.

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.