Transcription Factor AsMYC2 Controls the Jasmonate-Responsive Expression of ASS1 Regulating Sesquiterpene Biosynthesis in Aquilaria sinensis (Lour.) Gilg

Identification of Transcription Factors of Santalene Synthase Gene Promoters and SaSSY Cis-Elements through Yeast One-Hybrid Screening in Santalum album L

The main components of sandalwood heartwood essential oil are terpenoids, approximately 80% of which are α-santalol and β-santalol. In the synthesis of the main secondary metabolites of sandalwood heartwood, the key gene, santalene synthase (SaSSY), can produce α-santalene and β-santalene by catalyzed (E, E)-FPP. Furthermore, santalene is catalyzed by the cytochrome monooxygenase SaCYP736A167 to form sandalwood essential oil, which then produces a fragrance. However, the upstream regulatory mechanism of the key gene santalene synthase remains unclear. In this study, SaSSY (Sal3G10690) promoter transcription factors and SaSSY cis-elements were screened. The results showed that the titer of the sandalwood cDNA library was 1.75 × 107 CFU/mL, 80% of the inserted fragments identified by PCR were over 750 bp in length, and the positivity rate of the library was greater than 90%. The promoter region of the SaSSY gene was shown to have the structural basis for potential regulatory factor binding. After sequencing and bioinformatics analysis, we successfully obtained 51 positive clones and identified four potential SaSSY transcriptional regulators. Sal6G03620 was annotated as the transcription factor MYB36-like, and Sal8G07920 was annotated as the small heat shock protein HSP20 in sandalwood. Sal1G00910 was annotated as a hypothetical protein of sandalwood. Sal4G10880 was annotated as a homeobox-leucine zipper protein (ATHB-15) in sandalwood. In this study, a cDNA library of sandalwood was successfully constructed using a yeast one-hybrid technique, and the transcription factors that might interact with SaSSY gene promoters were screened. This study provides a foundation for exploring the molecular regulatory mechanism involved in the formation of sandalwood heartwood.

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.

The dual role of GoPGF reveals that the pigment glands are synthetic sites of gossypol in aerial parts of cotton

Gossypol and the related terpenoids are stored in the pigment gland to protect cotton plants from biotic stresses, but little is known about the synthetic sites of these metabolites. Here, we showed that GoPGF, a key gene regulating gland formation, was expressed in gland cells and roots. The chromatin immunoprecipitation sequencing (ChIP-seq) analysis demonstrated that GoPGF targets GhJUB1 to regulate gland morphogenesis. RNA-sequencing (RNA-seq) showed high accumulation of gossypol biosynthetic genes in gland cells. Moreover, integrated analysis of the ChIP-seq and RNA-seq data revealed that GoPGF binds to the promoter of several gossypol biosynthetic genes. The cotton callus overexpressing GoPGF had dramatically increased the gossypol levels, indicating that GoPGF can directly activate the biosynthesis of gossypol. In addition, the gopgf mutant analysis revealed the existence of both GoPGF-dependent and -independent regulation of gossypol production in cotton roots. Our study revealed that the pigment glands are synthetic sites of gossypol in aerial parts of cotton and that GoPGF plays a dual role in regulating gland morphogenesis and gossypol biosynthesis. The study provides new insights for exploring the complex relationship between glands and the metabolites they store in cotton and other plant species.

Transcription factors TgbHLH95 and TgbZIP44 cotarget terpene biosynthesis gene TgGPPS in Torreya grandis nuts

Terpenes are volatile compounds responsible for aroma and the postharvest quality of commercially important xiangfei (Torreya grandis) nuts, and there is interest in understanding the regulation of their biosynthesis. Here, a transcriptomics analysis of xiangfei nuts after harvest identified 156 genes associated with the terpenoid metabolic pathway. A geranyl diphosphate (GPP) synthase (TgGPPS) involved in production of the monoterpene precursor GPP was targeted for functional characterization, and its transcript levels positively correlated with terpene levels. Furthermore, transient overexpression of TgGPPS in tobacco (Nicotiana tabacum) leaves or tomato (Solanum lycopersicum) fruit led to monoterpene accumulation. Analysis of differentially expressed transcription factors identified one basic helix-loop-helix protein (TgbHLH95) and one basic leucine zipper protein (TgbZIP44) as potential TgGPPS regulators. TgbHLH95 showed significant transactivation of the TgGPPS promoter, and its transient overexpression in tobacco leaves led to monoterpene accumulation, whereas TgbZIP44 directly bound to an ACGT-containing element in the TgGPPS promoter, as determined by yeast 1-hybrid test and electrophoretic mobility shift assay. Bimolecular fluorescence complementation, firefly luciferase complementation imaging, co-immunoprecipitation, and GST pull-down assays confirmed a direct protein-protein interaction between TgbHLH95 and TgbZIP44 in vivo and in vitro, and in combination these proteins induced the TgGPPS promoter up to 4.7-fold in transactivation assays. These results indicate that a TgbHLH95/TgbZIP44 complex activates the TgGPPS promoter and upregulates terpene biosynthesis in xiangfei nuts after harvest, thereby contributing to its aroma.

Unraveling the regulatory role of MYC2 on ASMT gene expression in wheat: Implications for melatonin biosynthesis and drought tolerance

Recognized for its multifaceted functions, melatonin is a hormone found in both animals and plants. In the plant kingdom, it plays diverse roles, regulating growth, development, and stress responses. Notably, melatonin demonstrates its significance by mitigating the effects of abiotic stresses like drought. However, understanding the precise regulatory mechanisms controlling melatonin biosynthesis genes, especially during monocots' response to stresses, requires further exploration. Seeking to understand the molecular basis of drought stress tolerance in wheat, we analyzed RNA-Seq libraries of wheat exposed to drought stress using bioinformatics methods. In light of our findings, we identified that the Myelocytomatosis oncogenes 2 (MYC2) transcription factor is a hub gene upstream of a main melatonin biosynthesis gene, N-acetylserotonin methyltransferase (ASMT), in the wheat drought response-gene network. Promoter analysis of the ASMT gene suggested that it might be a target gene of MYC2. We conducted a set of molecular and physiochemical assays along with robust machine learning approaches to elevate those findings further. MYC2 and ASMT were co-regulated under Jasmonate, drought, and a combination of them in the leaf tissues of wheat was detected. A meaningful correlation was observed among gene expression profiles, melatonin contents, photosynthetic activities, antioxidant enzyme activities, H2 O2 levels, and plasma membrane damage. The results indicated an evident relationship between jasmonic acid and the melatonin biosynthesis pathway. Moreover, it seems that the MYC2-ASMT module might contribute to wheat drought tolerance by regulating melatonin contents.

Wound-Induced Temporal Reprogramming of Gene Expression during Agarwood Formation in Aquilaria sinensis

Agarwood is a resinous heartwood of Aquilaria sinensis that is formed in response to mechanical wounding. However, the transcriptional response of A. sinensis to mechanical wounding during the agarwood formation process is still unclear. Here, three five-year-old A. sinensis trees were mechanically damaged by a chisel, and time-series transcriptomic analysis of xylem tissues in the treated area (TA) was performed at 15 (TA1), 70 (TA2) and 180 days after treatment (TA3). Samples from untreated areas at the corresponding time points (UA1, UA2, UA3, respectively) were collected as controls. A total of 1862 (TA1 vs. UA1), 961 (TA2 vs. UA2), 1370 (TA3 vs. UA3), 3305 (TA2 vs. TA1), 2625 (TA3 vs. TA1), 2899 (TA3 vs. TA2), 782 (UA2 vs. UA1), 4443 (UA3 vs. UA1) and 4031 (UA3 vs. UA2) genes were differentially expressed (DEGs). Functional enrichment analysis showed that DEGs were significantly enriched for secondary metabolic processes, signal transduction and transcriptional regulation processes. Most of the genes involved in lignin biosynthesis were more abundant in the TA groups, which included phenylalanine ammonia-lyase, 4-coumarate CoA ligase, cinnamate 4-hydroxylase, caffeoyl-CoA O-methyltransferase and cinnamoyl-CoA reductase. DEGs involved in sesquiterpene biosynthesis were also identified. Hydroxymethylglutaryl-CoA synthase, 3-hydroxy-3-methylglutaryl-coenzyme A reductase, phosphomevalonate kinase and terpene synthase genes were significantly increased in the TA groups, promoting sesquiterpene biosynthesis in the wounded xylem tissues. The TF-gene transcriptomic networks suggested that MYB DNA-binding, NAM, WRKY, HLH and AP2 TFs co-expressed with genes related to lignin and sesquiterpene synthesis, indicating their critical regulatory roles in the biosynthesis of these compounds. Overall, our study reveals a dynamic transcriptional response of A. sinensis to mechanical wounding, provides a resource for identifying candidate genes for molecular breeding of agarwood quality, and sheds light on the molecular mechanisms of agarwood formation in A. sinensis.

PnMYB4 negatively modulates saponin biosynthesis in Panax notoginseng through interplay with PnMYB1

Saponins are the main triterpenoid ingredients from Panax notoginseng, a well-known Chinese medicine, and are important sources for producing drugs to prevent and treat cerebrovascular and cardiovascular diseases. However, the transcriptional regulatory network of saponin biosynthesis in P. notoginseng is largely unknown. In the present study we demonstrated that one R2R3-MYB transcription factor, designated PnMYB4, acts as a repressor of saponin accumulation. Suppression of PnMYB4 in P. notoginseng calli significantly increased the saponin content and the expression level of saponin biosynthetic genes. PnMYB4 directly bound to the promoters of key saponin biosynthetic genes, including PnSS, PnSE, and PnDS, to repress saponin accumulation. PnMYB4 and the activator PnMYB1 could interacted with PnbHLH, which is a positive regulator of saponin biosynthesis, to modulate the biosynthesis of saponin. PnMYB4 competed with PnMYB1 for binding to PnbHLH, repressing activation of the promoters of saponin structural genes induced by the PnMYB1-PnbHLH complex. Our study reveals that a complex regulatory module of saponin biosynthesis is associated with positive and negative MYB transcriptional regulators and provides a theoretical basis for improving the content of saponins and efficacy of P. notoginseng.

TcbHLH14 a Jasmonate Associated MYC2-like Transcription Factor Positively Regulates Pyrethrin Biosynthesis in Tanacetum cinerariifolium

Natural pyrethrins have high application value, and are widely used as a green pesticide in crop pest prevention and control. Pyrethrins are mainly extracted from the flower heads of Tanacetum cinerariifolium; however, the natural content is low. Therefore, it is essential to understand the regulatory mechanisms underlying the synthesis of pyrethrins through identification of key transcription factors. We identified a gene encoding a MYC2-like transcription factor named TcbHLH14 from T. cinerariifolium transcriptome, which is induced by methyl jasmonate. In the present study, we evaluated the regulatory effects and mechanisms of TcbHLH14 using expression analysis, a yeast one-hybrid assay, electrophoretic mobility shift assay, and overexpression/virus-induced gene silencing experiments. We found that TcbHLH14 can directly bind to the cis-elements of the pyrethrins synthesis genes TcAOC and TcGLIP to activate their expression. The transient overexpression of TcbHLH14 enhanced expression of the TcAOC and TcGLIP genes. Conversely, transient silencing of TcbHLH14 downregulated the expression of TcAOC and TcGLIP and reduced the content of pyrethrins. In summary, these results indicate that the potential application of TcbHLH14 in improving the germplasm resources and provide a new insight into the regulatory network of pyrethrins biosynthesis of T. cinerariifolium to further inform the development of engineering strategies for increasing pyrethrins contents.

Biosynthesis and the Transcriptional Regulation of Terpenoids in Tea Plants (Camellia sinensis)

Terpenes, especially volatile terpenes, are important components of tea aroma due to their unique scents. They are also widely used in the cosmetic and medical industries. In addition, terpene emission can be induced by herbivory, wounding, light, low temperature, and other stress conditions, leading to plant defense responses and plant-plant interactions. The transcriptional levels of important core genes (including HMGR, DXS, and TPS) involved in terpenoid biosynthesis are up- or downregulated by the MYB, MYC, NAC, ERF, WRKY, and bHLH transcription factors. These regulators can bind to corresponding cis-elements in the promoter regions of the corresponding genes, and some of them interact with other transcription factors to form a complex. Recently, several key terpene synthesis genes and important transcription factors involved in terpene biosynthesis have been isolated and functionally identified from tea plants. In this work, we focus on the research progress on the transcriptional regulation of terpenes in tea plants (Camellia sinensis) and thoroughly detail the biosynthesis of terpene compounds, the terpene biosynthesis-related genes, the transcription factors involved in terpene biosynthesis, and their importance. Furthermore, we review the potential strategies used in studying the specific transcriptional regulation functions of candidate transcription factors that have been discriminated to date.

Exploring the metabolic basis of growth/defense trade-offs in complex environments with Nicotiana attenuata plants cosilenced in NaMYC2a/b expression

In response to challenges from herbivores and competitors, plants use fitness-limiting resources to produce (auto)toxic defenses. Jasmonate signaling, mediated by MYC2 transcription factors (TF), is thought to reconfigure metabolism to minimize these formal costs of defense and optimize fitness in complex environments. To study the context-dependence of this metabolic reconfiguration, we cosilenced NaMYC2a/b by RNAi in Nicotiana attenuata and phenotyped plants in the field and increasingly realistic glasshouse setups with competitors and mobile herbivores. NaMYC2a/b had normal phytohormonal responses, and higher growth and fitness in herbivore-reduced environments, but were devastated in high herbivore-load environments in the field due to diminished accumulations of specialized metabolites. In setups with competitors and mobile herbivores, irMYC2a/b plants had lower fitness than empty vector (EV) in single-genotype setups but increased fitness in mixed-genotype setups. Correlational analyses of metabolic, resistance, and growth traits revealed the expected defense/growth associations for most sectors of primary and specialized metabolism. Notable exceptions were some HGL-DTGs and phenolamides that differed between single-genotype and mixed-genotype setups, consistent with expectations of a blurred functional trichotomy of metabolites. MYC2 TFs mediate the reconfiguration of primary and specialized metabolic sectors to allow plants to optimize their fitness in complex environments.

RNA-Sequencing Reveals the Involvement of Sesquiterpene Biosynthesis Genes and Transcription Factors during an Early Response to Mechanical Wounding of Aquilaria sinensis

Plants respond to wounding by reprogramming the expression of genes involved in secondary metabolism. Aquilaria trees produce many bioactive secondary metabolites in response to wounding, but the regulatory mechanism of agarwood formation in the early response to mechanical wounding has remained unclear. To gain insights into the process of transcriptome changes and to determine the regulatory networks of Aquilaria sinensis to an early response (15 days) to mechanical wounding, we collected A. sinensis samples from the untreated (Asc1) and treated (Asf1) xylem tissues and performed RNA sequencing (RNA-seq). This generated 49,102,523 (Asc1) and 45,180,981 (Asf1) clean reads, which corresponded to 18,927 (Asc1) and 19,258 (Asf1) genes, respectively. A total of 1596 differentially expressed genes (DEGs) were detected in Asf1 vs. Asc1 (|log2 (fold change)| ≥ 1, Padj ≤ 0.05), of which 1088 were up-regulated and 508 genes were down-regulated. GO and KEGG enrichment analysis of DEGs showed that flavonoid biosynthesis, phenylpropanoid biosynthesis, and sesquiterpenoid and triterpenoid biosynthesis pathways might play important roles in wound-induced agarwood formation. Based on the transcription factor (TF)-gene regulatory network analysis, we inferred that the bHLH TF family could regulate all DEGs encoding for farnesyl diphosphate synthase, sesquiterpene synthase, and 1-deoxy-D-xylulose-5-phosphate synthase (DXS), which contribute to the biosynthesis and accumulation of agarwood sesquiterpenes. This study provides insight into the molecular mechanism regulating agarwood formation in A. sinensis, and will be helpful in selecting candidate genes for improving the yield and quality of agarwood.

Comparative transcriptome analysis reveals different defence responses during the early stage of wounding stress in Chi-Nan germplasm and ordinary Aquilaria sinensis

BACKGROUND

Agarwood is a valuable Chinese medicinal herb and spice that is produced from wounded Aquilaria spp., is widely used in Southeast Asia and is highly traded on the market. The lack of highly responsive Aquilaria lines has seriously restricted agarwood yield and the development of its industry. In this article, a comparative transcriptome analysis was carried out between ordinary A. sinensis and Chi-Nan germplasm, which is a kind of A. sinensis tree with high agarwood-producing capacity in response to wounding stress, to elucidate the molecular mechanism underlying wounding stress in different A. sinensis germplasm resources and to help identify and breed high agarwood-producing strains.

RESULTS

A total of 2427 and 1153 differentially expressed genes (DEGs) were detected in wounded ordinary A. sinensis and Chi-Nan germplasm compared with the control groups, respectively. KEGG enrichment analysis revealed that genes participating in starch metabolism, secondary metabolism and plant hormone signal transduction might play major roles in the early regulation of wound stress. 86 DEGs related to oxygen metabolism, JA pathway and sesquiterpene biosynthesis were identified. The majority of the expression of these genes was differentially induced between two germplasm resources under wounding stress. 13 candidate genes related to defence and sesquiterpene biosynthesis were obtained by WGCNA. Furthermore, the expression pattern of genes were verified by qRT-PCR. The candidate genes expression levels were higher in Chi-Nan germplasm than that in ordinary A. sinensis during early stage of wounding stress, which may play important roles in regulating high agarwood-producing capacity in Chi-Nan germplasm.

CONCLUSIONS

Compared with A. sinensis, Chi-Nan germplasm invoked different biological processes in response to wounding stress. The genes related to defence signals and sesquiterepene biosynthesis pathway were induced to expression differentially between two germplasm resources. A total of 13 candidate genes were identified, which may correlate with high agarwood-producting capacity in Chi-Nan germplasm during the early stage of wounding stress. These genes will contribute to the development of functional molecular markers and the rapid breeding highly of responsive Aquilaria lines.

Integrating Multiple Omics Identifies Phaeoacremonium rubrigenum Acting as Aquilaria sinensis Marker Fungus to Promote Agarwood Sesquiterpene Accumulation by Inducing Plant Host Phosphorylation

The present study aimed to explore the factors that promote persistent agarwood accumulation. To this end, we first investigated the morphological changes and volatile compound distribution in five layers of "Guan Xiang" agarwood. The agarwood-normal transition layer (TL), an essential layer of persistent agarwood accumulation, showed clear metabolic differences by microscopy and GC-MS analysis. Microbiome analysis revealed that Phaeocremonium rubrigenum was the predominant biomarker fungus in the TL of "Guan Xiang" agarwood samples. Among the seven isolated fungi, P. rubrigenum exhibited a significantly heightened ability to induce the production in Aquilaria sinensis seedlings, especially for sesquiterpene. Tracing the proteome profile changes in P. rubrigenum-induced A. sinensis calli for 18 ds showed that the fungus-induced sesquiterpene biosynthesis increased mainly through the mevalonate (MVA) pathway. Specifically, the phosphorylation modification level, instead of the protein abundance of transcription factors (TFs), showed corresponding changes during sesquiterpene biosynthesis, thus indicating that induced phosphorylation is the key reason for enhanced sesquiterpene production. IMPORTANCE Agarwood is an expensive resinous portion derived from Aquilaria plants and has been widely used as medicine, incense, and perfume. The factors involved in steady agarwood accumulation remain elusive. Our current study suggests that as a TL marker fungus, P. rubrigenum could persistently promote agarwood sesquiterpene accumulation by inducing phosphorylation of the TFs-MVA network in A. sinensis. Moreover, our work provides strategies to improve agarwood industry management and sheds light on the potential molecular mechanisms of plant adaptation to native microbial conditions.

GhMYC2 activates cytochrome P450 gene CYP71BE79 to regulate gossypol biosynthesis in cotton

GhMYC2 regulates the gossypol biosynthesis pathway in cotton through activation of the expression of gossypol synthesis gene CYP71BE79, CDNC, CYP706B1, DH1, and CYP82D113. Cotton is one of the main cash crops globally. Cottonseed contains fiber, fat, protein, and starch, and has important economic value. However, gossypol in cottonseed seriously affects the development and utilization of cottonseed. Nonetheless, gossypol has great application potential in agriculture, medicine, and industry. Therefore, it is very important to study gossypol biosynthesis and its upstream regulatory pathways. It has been reported that the content of gossypol in hairy roots of cotton is regulated through jasmonic acid signaling; however, the specific molecular mechanism has not been revealed yet. We found that the expression of basic helix-loop-helix family transcription factor GhMYC2 was significantly upregulated after exogenous administration of methyl jasmonate to cotton seedlings, and the content of gossypol changed significantly with the variation of GhMYC2 expression. Further studies revealed that GhMYC2 could specifically bind to the G-Box in the promoter region of CDNC, CYP706B1, DH1, CYP82D113, CYP71BE79 to activate its expression and regulate gossypol synthesis, and its activation of CYP71BE79 promoter was inhibited by GhJAZ2. Not only that GhMYC2 could also interact with GoPGF. In this work, the molecular mechanisms of gossypol biosynthesis regulated by GhMYC2 were analyzed. The results provide a theoretical basis for cultivating new varieties of low-gossypol or high-gossypol cotton and creating excellent germplasm resources.

Identification and Functional Analysis of SabHLHs in Santalum album L

Santalum album L., a semi-parasitic evergreen tree, contains economically important essential oil, rich in sesquiterpenoids, such as (Z) α- and (Z) β-santalol. However, their transcriptional regulations are not clear. Several studies of other plants have shown that basic-helix-loop-helix (bHLH) transcription factors (TFs) were involved in participating in the biosynthesis of sesquiterpene synthase genes. Herein, bHLH TF genes with similar expression patterns and high expression levels were screened by co-expression analysis, and their full-length ORFs were obtained. These bHLH TFs were named SaMYC1, SaMYC3, SaMYC4, SaMYC5, SabHLH1, SabHLH2, SabHLH3, and SabHLH4. All eight TFs had highly conserved bHLH domains and SaMYC1, SaMYC3, SaMYC4, and SaMYC5, also had highly conserved MYC domains. It was indicated that the eight genes belonged to six subfamilies of the bHLH TF family. Among them, SaMYC1 was found in both the nucleus and the cytoplasm, while SaMYC4 was only localized in the cytoplasm and the remaining six TFs were localized in nucleus. In a yeast one-hybrid experiment, we constructed decoy vectors pAbAi-SSy1G-box, pAbAi-CYP2G-box, pAbAi-CYP3G-box, and pAbAi-CYP4G-box, which had been transformed into yeast. We also constructed pGADT7-SaMYC1 and pGADT7-SabHLH1 capture vectors and transformed them into bait strains. Our results showed that SaMYC1 could bind to the G-box of SaSSy, and the SaCYP736A167 promoter, which SaSSy proved has acted as a key enzyme in the synthesis of santalol sesquiterpenes and SaCYP450 catalyzed the ligation of santalol sesquiterpenes into terpene. We have also constructed pGreenII 62-SK-SaMYC1, pGreenII 0800-LUC-SaSSy and pGreenII 0800-LUC-SaCYP736A167 via dual-luciferase fusion expression vectors and transformed them into Nicotiana benthamiana using an Agrobacterium-mediated method. The results showed that SaMYC1 was successfully combined with SaSSy or SaCYP736A167 promoter and the LUC/REN value was 1.85- or 1.55-fold higher, respectively, than that of the control group. Therefore, we inferred that SaMYC1 could activate both SaSSy and SaCYP736A167 promoters.

Identification of a diarylpentanoid-producing polyketide synthase revealing an unusual biosynthetic pathway of 2-(2-phenylethyl)chromones in agarwood

2-(2-Phenylethyl)chromones (PECs) are the principal constituents contributing to the distinctive fragrance of agarwood. How PECs are biosynthesized is currently unknown. In this work, we describe a diarylpentanoid-producing polyketide synthase (PECPS) identified from Aquilaria sinensis. Through biotransformation experiments using fluorine-labeled substrate, transient expression of PECPS in Nicotiana benthamiana, and knockdown of PECPS expression in A. sinensis calli, we demonstrate that the C6-C5-C6 scaffold of diarylpentanoid is the common precursor of PECs, and PECPS plays a crucial role in PECs biosynthesis. Crystal structure (1.98 Å) analyses and site-directed mutagenesis reveal that, due to its small active site cavity (247 Å3), PECPS employs a one-pot formation mechanism including a "diketide-CoA intermediate-released" step for the formation of the C6-C5-C6 scaffold. The identification of PECPS, the pivotal enzyme of PECs biosynthesis, provides insight into not only the feasibility of overproduction of pharmaceutically important PECs using metabolic engineering approaches, but also further exploration of how agarwood is formed.

TcMYC2 regulates Pyrethrin biosynthesis in Tanacetum cinerariifolium

Pyrethrins constitute a class of terpene derivatives with high insecticidal activity and are mainly synthesized in the capitula of the horticulturally important plant, Tanacetum cinerariifolium. Treatment of T. cinerariifolium with methyl jasmonate (MeJA) in the field induces pyrethrin biosynthesis, but the mechanism linking MeJA with pyrethrin biosynthesis remains unclear. In this study, we explored the transcription factors involved in regulating MeJA-induced pyrethrin biosynthesis. A single spray application of MeJA to T. cinerariifolium leaves rapidly upregulated the expression of most known pyrethrin biosynthesis genes and subsequently increased the total pyrethrin content in the leaf. A continuous 2-week MeJA treatment resulted in enhanced pyrethrin content and increased trichome density. TcMYC2, a key gene in jasmonate signaling, was screened at the transcriptome after MeJA treatment. TcMYC2 positively regulated expression of the pyrethrin biosynthesis genes TcCHS, TcAOC, and TcGLIP by directly binding to E-box/G-box motifs in the promoters. The stable overexpression of TcMYC2 in T. cinerariifolium hairy roots significantly increased the expression of TcAOC and TcGLIP. Further transient overexpression and viral-induced gene-silencing experiments demonstrated that TcMYC2 positively promoted pyrethrin biosynthesis. Collectively, the results reveal a novel molecular mechanism for MeJA-induced pyrethrin biosynthesis in T. cinerariifolium involving TcMYC2.

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.

Molecular Cloning and Functional Analysis of 1-Deoxy-D-Xylulose 5-Phosphate Reductoisomerase from Santalum album

Sandalwood (Santalum album L.) heartwood-derived essential oil contains a high content of sesquiterpenoids that are economically highly valued and widely used in the fragrance industry. Sesquiterpenoids are biosynthesized via the mevalonate acid and methylerythritol phosphate (MEP) pathways, which are also the sources of precursors for photosynthetic pigments. 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) is a secondary rate-limiting enzyme in the MEP pathway. In this paper, the 1416-bp open reading frame of SaDXR and its 897-bp promoter region, which contains putative conserved cis-elements involved in stress responsiveness (HSE and TC-rich repeats), hormone signaling (abscisic acid, gibberellin and salicylic acid) and light responsiveness, were cloned from 7-year-old S. album trees. A bioinformatics analysis suggested that SaDXR encodes a functional and conserved DXR protein. SaDXR was widely expressed in multiple tissues, including roots, twigs, stem sapwood, leaves, flowers, fruit and stem heartwood, displaying significantly higher levels in tissues with photosynthetic pigments, like twigs, leaves and flowers. SaDXR mRNA expression increased in etiolated seedlings exposed to light, and the content of chlorophylls and carotenoids was enhanced in all 35S::SaDXR transgenic Arabidopsis thaliana lines, consistent with the SaDXR expression level. SaDXR was also stimulated by MeJA and H₂O₂ in seedling roots. α-Santalol content decreased in response to fosmidomycin, a DXR inhibitor. These results suggest that SaDXR plays an important role in the biosynthesis of photosynthetic pigments, shifting the flux to sandalwood-specific sesquiterpenoids.

MYC2 Transcription Factors TwMYC2a and TwMYC2b Negatively Regulate Triptolide Biosynthesis in Tripterygium wilfordii Hairy Roots

Triptolide, an important bioactive diterpenoid extracted from the plant Tripterygium wilfordii, exhibits many pharmacological activities. MYC2 transcription factor (TF) plays an important role in the regulation of various secondary metabolites in plants. However, whether MYC2 TF could regulate the biosynthesis of triptolide in T. wilfordii is still unknown. In this study, two homologous MYC2 TF genes, TwMYC2a and TwMYC2b, were isolated from T. wilfordii hairy roots and functionally characterized. The analyses of the phylogenetic tree and subcellular localization showed that they were grouped into the IIIe clade of the bHLH superfamily with other functional MYC2 proteins and localized in the nucleus. Furthermore, yeast one-hybrid and GUS transactivation assays suggested that TwMYC2a and TwMYC2b inhibited the promoter activity of the miltiradiene synthase genes, TwTPS27a and TwTPS27b, by binding to the E-box (CACATG) and T/G-box (CACGTT) motifs in their promoters. Transgenic results revealed that RNA interference of TwMYC2a/b significantly enhanced the triptolide accumulation in hairy roots and liquid medium by upregulating the expression of several key biosynthetic genes, including TwMS (TwTPS27a/b), TwCPS (TwTPS7/9), TwDXR, and TwHMGR1. In summary, our findings show that TwMYC2a and TwMYC2b act as two negative regulators of triptolide biosynthesis in T. wilfordii hairy roots and also provide new insights on metabolic engineering of triptolide in the future.

Natural products in agarwood and Aquilaria plants: chemistry, biological activities and biosynthesis

Covering: Up to the end of 2019.Agarwood is a resinous portion of Aquilaria trees, which is formed in response to environmental stress factors such as physical injury or microbial attack. It is very sought-after among the natural incenses, as well as for its medicinal properties in traditional Chinese and Ayurvedic medicine. Interestingly, the chemical constituents of agarwood and healthy Aquilaria trees are quite different. Sesquiterpenes and 2-(2-phenethyl)chromones with diverse scaffolds commonly accumulate in agarwood. Similar structures have rarely been reported from the original trees that mainly contain flavonoids, benzophenones, xanthones, lignans, simple phenolic compounds, megastigmanes, diterpenoids, triterpenoids, steroids, alkaloids, etc. This review summarizes the chemical constituents and biological activities both in agarwood and Aquilaria trees, and their biosynthesis is discussed in order to give a comprehensive overview of the research progress on agarwood.

Chromosomal-level reference genome of the incense tree Aquilaria sinensis

Trees in the genus Aquilaria (Thymelaeaceae) are known as lign aloes, and are native to the forests of southeast Asia. Lign aloes produce agarwood as an antimicrobial defence. Agarwood has a long history of cultural and medicinal use, and is of considerable commercial value. However, due to habitat destruction and over collection, lign aloes are threatened in the wild. We present a chromosomal‐level assembly for Aquilaria sinensis, a lign aloe endemic to China known as the incense tree, based on Illumina short‐read, 10X Genomics linked‐read, and Hi‐C sequencing data. Our 783.8 Mbp A. sinensis genome assembly is of high physical contiguity, with a scaffold N50 of 87.6 Mbp, and high completeness, with a 95.8% BUSCO score for eudicotyledon genes. We include 17 transcriptomes from various plant tissues, providing a total of 35,965 gene models. We reveal the first complete set of genes involved in sesquiterpenoid production, plant defence, and agarwood production for the genus Aquilaria, including genes involved in the biosynthesis of sesquiterpenoids via the mevalonic acid (MVA), 1‐deoxy‐D‐xylulose‐5‐phosphate (DXP), and methylerythritol phosphate (MEP) pathways. We perform a detailed repeat content analysis, revealing that transposable elements account for ~61% of the genome, with major contributions from gypsy‐like and copia‐like LTR retroelements. We also provide a comparative analysis of repeat content across sequenced species in the order Malvales. Our study reveals the first chromosomal‐level genome assembly for a tree in the genus Aquilaria and provides an unprecedented opportunity to address a variety of applied, genomic and evolutionary questions in the Thymelaeaceae more widely.

Genome-wide survey of the bHLH super gene family in Brassica napus

BACKGROUND

The basic helix-loop-helix (bHLH) gene family is one of the largest transcription factor families in plants and is functionally characterized in diverse species. However, less is known about its functions in the economically important allopolyploid oil crop, Brassica napus.

RESULTS

We identified 602 potential bHLHs in the B. napus genome (BnabHLHs) and categorized them into 35 subfamilies, including seven newly separated subfamilies, based on phylogeny, protein structure, and exon-intron organization analysis. The intron insertion patterns of this gene family were analyzed and a total of eight types were identified in the bHLH regions of BnabHLHs. Chromosome distribution and synteny analyses revealed that hybridization between Brassica rapa and Brassica oleracea was the main expansion mechanism for BnabHLHs. Expression analyses showed that BnabHLHs were widely in different plant tissues and formed seven main patterns, suggesting they may participate in various aspects of B. napus development. Furthermore, when roots were treated with five different hormones (IAA, auxin; GA3, gibberellin; 6-BA, cytokinin; ABA, abscisic acid and ACC, ethylene), the expression profiles of BnabHLHs changed significantly, with many showing increased expression. The induction of five candidate BnabHLHs was confirmed following the five hormone treatments via qRT-PCR. Up to 246 BnabHLHs from nine subfamilies were predicted to have potential roles relating to root development through the joint analysis of their expression profiles and homolog function.

CONCLUSION

The 602 BnabHLHs identified from B. napus were classified into 35 subfamilies, and those members from the same subfamily generally had similar sequence motifs. Overall, we found that BnabHLHs may be widely involved in root development in B. napus. Moreover, this study provides important insights into the potential functions of the BnabHLHs super gene family and thus will be useful in future gene function research.

Genome-wide analysis of WRKY transcription factors in Aquilaria sinensis (Lour.) Gilg

The WRKY proteins are a superfamily of transcription factor that regulate diverse developmental and physiological processes in plants. Completion of the whole-genome sequencing of Aquilaria sinensis allowed us to perform a genome-wide investigation for WRKY proteins. Here, we predicted 70 WRKY genes from the A. sinensis genome and undertaken a comprehensive bioinformatic analysis. Due to their diverse structural features, the 70 AsWRKY genes are classified into three main groups (group I-III), with five subgroups (IIa-IIe) in group II, except two belong to none of them. Distinct expression profiles of AsWRKYs with RNA sequencing data revealed their diverse expression patterns among different tissues and in the process of whole-tree-inducing agarwood formation. Based on the expression characteristics, we predict some AsWRKYs are pseudogenes, and some may be involved in the biosynthesis of agarwood sesquiterpenes as activators or repressors. Among the tested genes treated with MeJA and H₂O₂, most of them are induced by H₂O₂, but downregulated by MeJA, implying the complexity of their involvement in signal transduction regulation. Our results not only provide a basic platform for functional identification of WRKYs in A. sinensis but important clues for further analysis their regulation role in agarwood formation.

Genome-wide identification and characterization of TALE superfamily genes in cotton reveals their functions in regulating secondary cell wall biosynthesis

BACKGROUND

Cotton fiber length and strength are both key traits of fiber quality, and fiber strength (FS) is tightly correlated with secondary cell wall (SCW) biosynthesis. The three-amino-acid-loop-extension (TALE) superclass homeoproteins are involved in regulating diverse biological processes in plants, and some TALE members has been identified to play a key role in regulating SCW formation. However, little is known about the functions of TALE members in cotton (Gossypium spp.).

RESULTS

In the present study, based on gene homology, 46, 47, 88 and 94 TALE superfamily genes were identified in G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. Phylogenetic and evolutionary analysis showed the evolutionary conservation of two cotton TALE families (including BEL1-like and KNOX families). Gene structure analysis also indicated the conservation of GhTALE members under selection. The analysis of promoter cis-elements and expression patterns suggested potential transcriptional regulation functions in fiber SCW biosynthesis and responses to some phytohormones for GhTALE proteins. Genome-wide analysis of colocalization of TALE transcription factors with SCW-related QTLs revealed that some BEL1-like genes and KNAT7 homologs may participate in the regulation of cotton fiber strength formation. Overexpression of GhKNAT7-A03 and GhBLH6-A13 significantly inhibited the synthesis of lignocellulose in interfascicular fibers of Arabidopsis. Yeast two-hybrid (Y2H) experiments showed extensive heteromeric interactions between GhKNAT7 homologs and some GhBEL1-like proteins. Yeast one-hybrid (Y1H) experiments identified the upstream GhMYB46 binding sites in the promoter region of GhTALE members and defined the downstream genes that can be directly bound and regulated by GhTALE heterodimers.

CONCLUSION

We comprehensively identified TALE superfamily genes in cotton. Some GhTALE members are predominantly expressed during the cotton fiber SCW thicking stage, and may genetically correlated with the formation of FS. Class II KNOX member GhKNAT7 can interact with some GhBEL1-like members to form the heterodimers to regulate the downstream targets, and this regulatory relationship is partially conserved with Arabidopsis. In summary, this study provides important clues for further elucidating the functions of TALE genes in regulating cotton growth and development, especially in the fiber SCW biosynthesis network, and it also contributes genetic resources to the improvement of cotton fiber quality.

Genome-wide analysis of basic helix-loop-helix superfamily members related to anthocyanin biosynthesis in eggplant (Solanum melongena L.)

The basic helix–loop–helix (bHLH) superfamily is considered the second largest transcription factor (TF) family. It plays regulatory roles in the developmental processes of plants and in their defense responses. In recent years, many bHLH superfamily genes have been identified and characterized in herbaceous and woody plants. However, the comprehensive genomic and functional analyses of these genes in eggplant (Solanum melongena L.) have not been reported. In this study, 121 bHLH TFs were identified in the recently released eggplant genome. The phylogeny, gene structure and conserved motifs of the SmbHLH gene were comprehensively studied. Subsequently, the phylogenetic relationship between the bHLH of eggplant and the bHLH of other species was analyzed, and the proteins were classified into 17 subfamilies. Among these protein sequences, 16 subgroups were clustered into the functional clades of Arabidopsis. Two candidate genes (SmbHLH1, SmbHLH117) that may be involved in anthocyanin biosynthesis were screened. The tissue specificity or differential expression of the bHLH genes in different tissues and under various light and temperature conditions suggested the differential regulation of tissue development and metabolism. This study not only provides a solid foundation for the functional dissection of the eggplant bHLH gene family but may also be useful for the future synthesis of anthocyanins in eggplant.

Characteristics and Expression Pattern of MYC Genes in Triticum aestivum, Oryza sativa, and Brachypodium distachyon

Myelocytomatosis oncogenes (MYC) transcription factors (TFs) belong to basic helix-loop-helix (bHLH) TF family and have a special bHLH_MYC_N domain in the N-terminal region. Presently, there is no detailed and systematic analysis of MYC TFs in wheat, rice, and Brachypodium distachyon. In this study, 26 TaMYC, 7 OsMYC, and 7 BdMYC TFs were identified and their features were characterized. Firstly, they contain a JAZ interaction domain (JID) and a putative transcriptional activation domain (TAD) in the bHLH_MYC_N region and a BhlH region in the C-terminal region. In some cases, the bHLH region is followed by a leucine zipper region; secondly, they display tissue-specific expression patterns: wheat MYC genes are mainly expressed in leaves, rice MYC genes are highly expressed in stems, and B. distachyon MYC genes are mainly expressed in inflorescences. In addition, three types of cis-elements, including plant development/growth-related, hormone-related, and abiotic stresses-related were identified in different MYC gene promoters. In combination with the previous studies, these results indicate that MYC TFs mainly function in growth and development, as well as in response to stresses. This study laid a foundation for the further functional elucidation of MYC genes.

The Scent of Stress: Evidence From the Unique Fragrance of Agarwood

Agarwood (Aquilaria spp.) fragrance and its origin in stress make it probably the most suitable model to study stress-induced aroma. Production being confined only to certain small pockets of South and Southeast Asia, agarwood is arguably the costliest wood in the world. Formation of fragrant agarwood resin is the outcome of complex biotic, abiotic, and physical stress on the Aquilaria trees. The intricate mechanism by which some 150 odd fragrant molecules that constitute agarwood aroma is formed is still not clearly understood. The present review therefore aims to bring to focus this less known but highly valuable stress-induced aroma from Asia. Discussions on agarwood species, occurrence, distribution, formation, and products have been included as foundation. Although global trade in agarwood and its products is estimated at US$6 billion to US$8 billion, no reliable data are readily available in literature. Therefore, an effort has been made to review the current status of agarwood trade. The element of stress and its correlation to agarwood aroma is discussed in the subsequent sections. Natural agarwood formation as well as technologies and interventions in agarwood induction are stress-based (natural and artificial injury, insect and fungal attack, chemical induction). The molecular triggers are gradually coming to light as new studies are implicating jasmonate, LOX signaling, and other stress reaction routes as the source of agarwood aroma. This review therefore has strived to compile the information that is scattered across scientific as well as other authentic literature and update the reader on the current status. More information about the specific roles of other vital stressors like insects, abiotic, and genetic factors is eagerly awaited from ongoing and future research to further understand the unique fragrance of agarwood.

Identification and functional characterization of three new terpene synthase genes involved in chemical defense and abiotic stresses in Santalum album

BACKGROUND

It is well known that aromatic essential oils extracted from the heartwood of Santalum album L. have wide economic value. However, little is known about the role of terpenoids in response to various adverse environmental stresses as other plants do in the form of signals during plant-environment interactions.

RESULTS

In this study, trace amounts of volatiles consisting of α-santalene, epi-β-santalene, β-santalene, α-santalol, β-santalol, (E)-α-bergamotene, (E)-β-farnesene and β-bisabolene were found in the leaves of mature S. album trees. We identified more than 40 candidate terpene synthase (TPS) unigenes by mining publicly-available RNA-seq data and characterized the enzymes encoded by three cDNAs: one mono-TPS catalyzes the formation of mostly α-terpineol, and two multifunctional sesqui-TPSs, one of which produces (E)-α-bergamotene and sesquisabinene as major products and another which catalyzes the formation of (E)-β-farnesene, (E)-nerolidol and (E,E)-farnesol as main products. Metabolite signatures and gene expression studies confirmed that santalol content is closely related with santalene synthase (SaSSY) transcripts in heartwood, which is key enzyme responsible for santalol biosynthesis. However, the expression of three new SaTPS genes differed significantly from SaSSY in the essential oil-producing heartwood. Increased activities of antioxidant enzymes, superoxide dismutase, catalase, peroxidase and ascorbate peroxidase, were detected in different tissues of S. album plants after applying 1 mM methyl jasmonate (MeJA) and 1 mM salicylic acid (SA), or exposure to 4°C, 38°C and high light intensity. MeJA and SA dramatically induced the expression of SaTPS1 and SaTPS2 in leaves. SaTPS1 to 3 transcripts were differentially activated among different tissues under adverse temperature and light stresses. In contrast, almost all SaSSY transcripts decreased in response to these environmental stresses, unlike SaTPS1 to 3.

CONCLUSIONS

Multifunctional enzymes were biochemically characterized, including one chloroplastic mono-TPS and two cytosolic sesqui-TPSs in sandalwood. Our results suggest the ecological importance of these three new SaTPS genes in defensive response to biotic attack and abiotic stresses in S. album.

Hydrogen peroxide burst triggers accumulation of jasmonates and salicylic acid inducing sesquiterpene biosynthesis in wounded Aquilaria sinesis

Agarwood, a non-timber fragrant wood, is produced in wounded Aquilaria trees and widely used in perfume, incense, and medicine. Sesquiterpene is one of its main active compounds. It has been demonstrated that hydrogen peroxide (H₂O₂) plays a role in promoting agarwood sesquiterpene biosynthesis, but little is known about its signaling pathway. In this study, the pruning of actively growing saplings of A. sinensis resulted in an H₂O₂ burst and the accumulation of jasmonic acid (JA), salicylic acid (SA), and ethylene (ET), which was followed by the up-regulation of sesquiterpene synthase and the production of sesquiterpene in the pruned stems. This process could be enhanced by absorbed H₂O₂ and inhibited by an H₂O₂ scavenger (ascorbate, AsA) in pruned stems, although the concentration of ET and transcription of ET-related synthase genes remained unaffected. These results confirmed that the H₂O₂ burst in wounded stems triggered JA and SA accumulation to promote agarwood sesquiterpene biosynthesis. ET was also activated by injury that was independent with H₂O₂. All results excavated a full-scale signaling transduction nets among multiple stress signals during wound-induced agarwood production in A. sinensis and provide a new insight into improving the artificial technology of agarwood production.

Agarwood Induction: Current Developments and Future Perspectives

Agarwood is a resinous part of the non-timber Aquilaria tree, which is a highly valuable product for medicine and fragrance purposes. To protect the endangered Aquilaria species, mass plantation of Aquilaria trees has become a sustainable way in Asian countries to obtain the highly valuable agarwood. As only physiologically triggered Aquilaria tree can produce agarwood, effective induction methods are long sought in the agarwood industry. In this paper, we attempt to provide an overview for the past efforts toward the understanding of agarwood formation, the evolvement of induction methods and their further development prospects by integrating it with high-throughput omics approaches.

Genome-wide identification of the GhARF gene family reveals that GhARF2 and GhARF18 are involved in cotton fibre cell initiation

Auxin signalling plays an essential role in regulating plant development. Auxin response factors (ARFs), which are critical components of auxin signalling, modulate the expression of early auxin-responsive genes by binding to auxin response factor elements (AuxREs). However, there has been no comprehensive characterization of this gene family in cotton. Here, we identified 56 GhARF genes in the assembled Gossypium hirsutum genome. This gene family was divided into 17 subfamilies, and 44 members of them were distributed across 21 chromosomes. GhARF6 and GhARF11 subfamily genes were predominantly expressed in vegetative tissues, whereas GhARF2 and GhARF18 subfamily genes were highly expressed during seed fibre cell initiation. GhARF2-1 and GhARF18-1 were exclusively expressed in trichomes, organs similar to cotton seed fibre cells, and overexpression of these genes in Arabidopsis enhances trichome initiation. Comparative transcriptome analysis combined with AuxRE prediction revealed 11 transcription factors as potential target genes of GhARF2 and GhARF18. Six of these genes were significantly expressed during seed fibre cell initiation and were bound by GhARF2-1 and GhARF18-1 in yeast one-hybrid assays. Our results suggest that GhARF2 and GhARF18 genes may be key regulators of cotton seed fibre initiation by regulating the expression of several transcription factor genes. This study deepens our understanding of auxin-mediated initiation of cotton seed fibre cells and helps us in breeding better cotton varieties in the future.

Genome-Wide Identification and Characterization of the Potato bHLH Transcription Factor Family

Plant basic/helix–loop–helix (bHLH) transcription factors participate in a number of biological processes, such as growth, development and abiotic stress responses. The bHLH family has been identified in many plants, and several bHLH transcription factors have been functionally characterized in Arabidopsis. However, no systematic identification of bHLH family members has been reported in potato (Solanum tuberosum). Here, 124 StbHLH genes were identified and named according to their chromosomal locations. The intron numbers varied from zero to seven. Most StbHLH proteins had the highly conserved intron phase 0, which accounted for 86.2% of the introns. According to the Neighbor-joining phylogenetic tree, 259 bHLH proteins acquired from Arabidopsis and potato were divided into 15 groups. All of the StbHLH genes were randomly distributed on 12 chromosomes, and 20 tandem duplicated genes and four pairs of duplicated gene segments were detected in the StbHLH family. The gene ontology (GO) analysis revealed that StbHLH mainly function in protein and DNA binding. Through the RNA-seq and quantitative real time PCR (qRT-PCR) analyses, StbHLH were found to be expressed in various tissues and to respond to abiotic stresses, including salt, drought and heat. StbHLH1, 41 and 60 were highly expressed in flower tissues, and were predicted to be involved in flower development by GO annotation. StbHLH45 was highly expressed in salt, drought and heat stress, which suggested its important role in abiotic stress response. The results provide comprehensive information for further analyses of the molecular functions of the StbHLH gene family.