Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes

The spatio-temporal biosynthesis of floral flavonols is controlled by differential phylogenetic MYB regulators in Freesia hybrida

Floral flavonols play specific pivotal roles in pollinator attraction, pollen germination and fertility, in addition to other functions in vegetative organs. For many plants, the process of flavonol biosynthesis in late flower development stages and in mature flower tissues is poorly understood, in contrast to early flower development stages. It is thought that this process may be regulated independently of subgroup 7 R2R3 MYB (SG7 MYB) transcription factors. In this study, two FLS genes were shown to be expressed synchronously with the flower development-specific and tissue-specific biosynthesis of flavonols in Freesia hybrida. FhFLS1 contributed to flavonol biosynthesis in early flower buds, toruses and calyxes, and was regulated by four well-known SG7 MYB proteins, designated as FhMYBFs, with at least partial regulatory redundancy. FhFLS2 accounted for flavonols in late developed flowers and in the petals, stamens and pistils, and was targeted directly by non SG7 MYB protein FhMYB21L2. In parallel, AtMYB21 and AtMYB24 also activated AtFLS1, a gene highly expressed in Arabidopsis anthers and pollen, indicating the conserved regulatory roles of MYB21 against FLS genes in these two evolutionarily divergent angiosperm plants. Our results reveal a novel regulatory and synthetic mechanism underlying flavonol biosynthesis in floral organs and tissues which may be exploited to investigate supplementary roles of flavonols in flowers.

Tomato transcriptional repressor MYB70 directly regulates ethylene-dependent fruit ripening

Ethylene is a key plant hormone controlling the ripening of climacteric fruits, and several transcription factors acting as important regulators of fruit ripening have been identified in tomato (Solanum lycopersicum), a model for climacteric fruits. The vast majority of these transcription factors are transcriptional activators, however, and the associated transcriptional regulatory mechanisms of most regulators are unclear. Here, we report on a tomato transcriptional repressor (termed SlMYB70) that negatively regulates fruit ripening by directly modulating ethylene biosynthesis. As an EAR motif-containing MYB transcription factor-encoding gene, SlMYB70 displayed a ripening-associated expression pattern and was responsive to ethylene. RNA interference (RNAi)-mediated repression of SlMYB70 accelerated fruit ripening, but overexpression of SlMYB70 delayed fruit ripening. Ethylene production was noticeably increased and decreased in SlMYB70-RNAi and SlMYB70-overexpressing lines, respectively, compared with wild-type tomatoes. SlMYB70 was proven to be a transcriptional repressor, dependent on the EAR repression motif, and to repress the transcription of two ethylene biosynthesis genes in fruit ripening, namely SlACS2 and SlACO3. The promoters of SlACS2 and SlACO3 are directly bound by SlMYB70, which was verified using a combination of yeast one-hybrid chromatin immunoprecipitation quantitative polymerase chain reaction and electrophoretic mobility shift assays. These results suggest that SlMYB70 negatively regulates fruit ripening via the direct transcriptional repression of ethylene biosynthesis genes, which provides insights into the ethylene-mediated key regulatory hierarchy in climacteric fruit ripening, and also highlights different types of transcriptional regulation of fruit ripening.

Advances in Biosynthesis and Biological Functions of Proanthocyanidins in Horticultural Plants

Proanthocyanidins are colorless flavonoid polymers condensed from flavan-3-ol units. They are essential secondary plant metabolites that contribute to the nutritional value and sensory quality of many fruits and the related processed products. Mounting evidence has shown that the accumulation of proanthocyanidins is associated with the resistance of plants against a broad spectrum of abiotic and biotic stress conditions. The biosynthesis of proanthocyanidins has been examined extensively, allowing for identifying and characterizing the key regulators controlling the biosynthetic pathway in many plants. New findings revealed that these specific regulators were involved in the proanthocyanidins biosynthetic network in response to various environmental conditions. This paper reviews the current knowledge regarding the control of key regulators in the underlying proanthocyanidins biosynthetic and molecular mechanisms in response to environmental stress. Furthermore, it discusses the directions for future research on the metabolic engineering of proanthocyanidins production to improve food and fruit crop quality.

Transcriptome Dynamics during Black and White Sesame (Sesamum indicum L.) Seed Development and Identification of Candidate Genes Associated with Black Pigmentation

Seed coat color is a crucial agronomic trait in sesame (Sesamum indicum L.) since it is strongly linked to seed oil, proteins, and lignans contents, and also influences consumer preferences. In East Asia, black sesame seed is used in the treatment and the prevention of various diseases. However, in sesame, little is known about the establishment of the seed coat color, and only one gene has been reported to control black pigmentation. This study provides an overview of developing seeds transcriptome of two varieties of sesame "Zhongfengzhi No.1" (white seed) and "Zhongzhi No.33" (black seed) and shed light on genes involving in black seed formation. Until eight days post-anthesis (DPA), both the seeds of the two varieties were white. The black sesame seed turned to yellow between 9 and 11 DPA and then black between 12 and 14 DPA. The black and white sesame showed similar trend-expressed genes with the numbers increased at the early stages of seed development. The differentially expressed genes (DEGs) number increased with seed development in the two sesame varieties. We examined the DEGs and uncovered that more were up-regulated at the early stages. The DEGs between the black and white sesame were mainly enriched in 37 metabolic pathways, among which the flavonoid biosynthesis and biosynthesis of secondary metabolites were dominants. Furthermore, we identified 20 candidate genes associated with pigment biosynthesis in black sesame seed, among which 10 were flavonoid biosynthesis and regulatory genes. These genes also include isochorismate and polyphenol oxidase genes. By comparing the phenotypes and genes expressions of the black and white sesame seed at different development stages, this work revealed the important role of 8-14 DPA in black pigment biosynthesis and accumulation. Moreover, it unfolded candidate genes associated with black pigmentation in sesame. These findings provide a vast transcriptome dataset and list of genes that will be targeted for functional studies related to the molecular mechanism involved in biosynthesis and regulation of seed coat color in sesame.

Screening of key genes responsible for Pennisetum setaceum 'Rubrum' leaf color using transcriptome sequencing

Pennisetum setaceum 'Rubrum' is an ornamental grass plant that produces purple leaves in high-light environments and light purple or green leaves in low-light environments, the latter of which greatly reduces its aesthetic appeal. Therefore, we aimed to identify the key genes associated with leaf coloration and elucidate the molecular mechanisms involved in the color changes in P. setaceum 'Rubrum' leaves. We performed transcriptome sequencing of P. setaceum 'Rubrum' leaves before and after shading. A total of 19,043 differentially expressed genes were identified, and the numbers of upregulated and downregulated genes at T1 stage, when compared with their expression at the T0 stage, were 10,761 and 8,642, respectively. The possible pathways that determine P. setaceum 'Rubrum' leaf color included flavonoid biosynthesis, flavone and flavonol biosynthesis, and carotenoid biosynthesis. There were 31 differentially expressed genes related to chlorophyll metabolism, of which 21 were related to chlorophyll biosynthesis and 10 to chlorophyll degradation, as well as three transcription factors that may be involved in the regulation of chlorophyll degradation. There were 31 key enzyme genes involved in anthocyanin synthesis and accumulation in P. setaceum 'Rubrum' leaves, with four transcription factors that may be involved in the regulation of anthocyanin metabolism. The transcriptome data were verified and confirmed reliable by real-time fluorescence quantitative PCR analysis. These findings provide a genetic basis for improving leaf color in P. setaceum 'Rubrum.'

A full-length transcriptome and gene expression analysis reveal genes and molecular elements expressed during seed development in Gnetum luofuense

BACKGROUND

Gnetum is an economically important tropical and subtropical gymnosperm genus with various dietary, industrial and medicinal uses. Many carbohydrates, proteins and fibers accumulate during the ripening of Gnetum seeds. However, the molecular mechanisms related to this process remain unknown.

RESULTS

We therefore assembled a full-length transcriptome from immature and mature G. luofuense seeds using PacBio sequencing reads. We identified a total of 5726 novel genes, 9061 alternative splicing events, 3551 lncRNAs, 2160 transcription factors, and we found that 8512 genes possessed at least one poly(A) site. In addition, gene expression comparisons of six transcriptomes generated by Illumina sequencing showed that 14,323 genes were differentially expressed from an immature stage to a mature stage with 7891 genes upregulated and 6432 genes downregulated. The expression of 14 differentially expressed transcription factors from the MADS-box, Aux/IAA and bHLH families was validated by qRT-PCR, suggesting that they may have important roles in seed ripening of G. luofuense.

CONCLUSIONS

These findings provide a valuable molecular resource for understanding seed development of gymnosperms.

Characterization of Metabolites and Transcripts Involved in Flower Pigmentation in Primula vulgaris

Primula vulgaris exhibits a wide range of flower colors and is a valuable ornamental plant. The combination of flavonols/anthocyanins and carotenoids provides various colorations ranging from yellow to violet-blue. However, the complex metabolic networks and molecular mechanisms underlying the different flower colors of P. vulgaris remain unclear. Based on comprehensive analysis of morphological anatomy, metabolites, and gene expression in different-colored flowers of P. vulgaris, the mechanisms relating color-determining compounds to gene expression profiles were revealed. In the case of P. vulgaris flower color, hirsutin, rosinin, petunidin-, and cyanidin-type anthocyanins and the copigment herbacetin contributed to the blue coloration, whereas peonidin-, cyandin-, and delphinidin-type anthocyanins showed high accumulation levels in pink flowers. The color formation of blue and pink were mainly via the regulation of F3'5'H (c53168), AOMT (c47583, c44905), and 3GT (c50034). Yellow coloration was mainly due to gossypetin and carotenoid, which were regulated by F3H (c43100), F3 1 (c53714), 3GT (c53907) as well as many carotenoid biosynthetic pathway-related genes. Co-expression network and transient expression analysis suggested a potential direct link between flavonoid and carotenoid biosynthetic pathways through MYB transcription factor regulation. This work reveals that transcription changes influence physiological characteristics, and biochemistry characteristics, and subsequently results in flower coloration in P. vulgaris.

Identification of Anthocyanins-Related Glutathione S-Transferase (GST) Genes in the Genome of Cultivated Strawberry (Fragaria × ananassa)

Anthocyanins are responsible for the red color of strawberry, they are a subclass of flavonoids synthesized in cytosol and transferred to vacuole to form the visible color. Previous studies in model and ornamental plants indicated members of the glutathione S-transferase (GST) gene family were involved in vacuolar accumulation of anthocyanins. In the present study, a total of 130 FaGST genes were identified in the genome of cultivated strawberry (Fragaria × ananassa), which were unevenly distributed across the 28 chromosomes from the four subgenomes. Evolutionary analysis revealed the expansion of FaGST family was under stable selection and mainly drove by WGD/segmental duplication event. Classification and phylogenetic analysis indicated that all the FaGST genes were clarified into seven subclasses, among which FaGST1, FaGST37, and FaGST97 belonging to Phi class were closely related to FvRAP, an anthocyanin-related GST of wildwood strawberry, and this clade was clustered with other known anthocyanin-related GSTs. RNAseq-based expression analysis at different developmental stages of strawberry revealed that the expression of FaGST1, FaGST37, FaGST39, FaGST73, and FaGST97 was gradually increased during the fruit ripening, consistent with the anthocyanins accumulation. These expression patterns of those five FaGST genes were also significantly correlated with those of other anthocyanin biosynthetic genes such as FaCHI, FaCHS, and FaANS, as well as anthocyanin regulatory gene FaMYB10. These results indicated FaGST1, FaGST37, FaGST39, FaGST73, and FaGST97 may function in vacuolar anthocyanin accumulation in cultivated strawberry.

The strawberry transcription factor FaRAV1 positively regulates anthocyanin accumulation by activation of FaMYB10 and anthocyanin pathway genes

The RAV (related to ABI3/viviparous 1) group of transcription factors (TFs) play multifaceted roles in plant development and stress responses. Here, we show that strawberry (Fragaria × ananassa) FaRAV1 positively regulates anthocyanin accumulation during fruit ripening via a hierarchy of activation processes. Dual-luciferase assay screening of all fruit-expressed AP2/ERFs showed FaRAV1 had the highest transcriptional activation of the promoter of FaMYB10, a key activator of anthocyanin biosynthesis. Yeast one-hybrid and electrophoretic mobility shift assays indicated that FaRAV1 could directly bind to the promoter of FaMYB10. Transient overexpression of FaRAV1 in strawberry fruit increased FaMYB10 expression and anthocyanin production significantly. Correspondingly, transient RNA interference-induced silencing of FaRAV1 led to decreases in FaMYB10 expression and anthocyanin content. Transcriptome analysis of FaRAV1-overexpressing strawberry fruit revealed that transcripts of phenylpropanoid and flavonoid biosynthesis pathway genes were up-regulated. Luciferase assays showed that FaRAV1 could also activate the promoters of strawberry anthocyanin biosynthetic genes directly, revealing a second level of FaRAV1 action in promoting anthocyanin accumulation. These results show that FaRAV1 stimulates anthocyanin accumulation in strawberry both by direct activation of anthocyanin pathway gene promoters and by up-regulation of FaMYB10, which also positively regulates these genes.

Postharvest seed coat darkening in pinto bean (Phaseolus vulgaris) is regulated by Psd , an allele of the basic helix-loop-helix transcription factor P

Pinto bean (Phaseolus vulgaris) is one of the leading market classes of dry beans that is most affected by postharvest seed coat darkening. The process of seed darkening poses a challenge for bean producers and vendors as they encounter significant losses in crop value due to decreased consumer preference for darker beans. Here, we identified a novel allele of the P gene, Psd , responsible for the slow darkening seed coat in pintos, and identified trait-specific sequence polymorphisms which are utilized for the development of new gene-specific molecular markers for breeding. These tools can be deployed to help tackle this economically important issue for bean producers.

SUMMARY

Postharvest seed coat darkening in pinto bean is an undesirable trait that reduces the market value of the stored crop. Regular darkening (RD) pintos darken faster after harvest and accumulate higher level of proanthocyanidins (PAs) compared to slow darkening (SD) cultivars. Although the markers cosegregating with the SD trait have been known for some time, the SLOW DARKENING (Sd) gene identity had not been proven.Here, we identified Psd as a candidate for controlling the trait. Genetic complementation, transcript abundance, metabolite analysis, and inheritance study confirmed that Psd is the Sd gene. Psd is another allele of the P (Pigment) gene, whose loss-of-function alleles result in a white seed coat. Psd encodes a bHLH transcription factor with two transcript variants but only one is involved in PA biosynthesis. An additional glutamate residue in the activation domain, and/or an arginine to histidine substitution in the bHLH domain of the Psd-1 transcript in the SD cultivar is likely responsible for the reduced activity of this allele compared to the allele in a RD cultivar, leading to reduced PA accumulation.Overall, we demonstrate that a novel allele of P, Psd , is responsible for the SD phenotype, and describe the development of new, gene-specific, markers that could be utilized in breeding to resolve an economically important issue for bean producers.

Metabolomics analysis reveals Embden Meyerhof Parnas pathway activation and flavonoids accumulation during dormancy transition in tree peony

BACKGROUND

Bud dormancy is a sophisticated strategy which plants evolve to survive in tough environments. Endodormancy is a key obstacle for anti-season culture of tree peony, and sufficient chilling exposure is an effective method to promote dormancy release in perennial plants including tree peony. However, the mechanism of dormancy release is still poorly understood, and there are few systematic studies on the metabolomics during chilling induced dormancy transition.

RESULTS

The tree peony buds were treated with artificial chilling, and the metabolmics analysis was employed at five time points after 0-4 °C treatment for 0, 7, 14, 21 and 28 d, respectively. A total of 535 metabolites were obtained and devided into 11 groups including flavonoids, amino acid and its derivatives, lipids, organic acids and its derivates, nucleotide and its derivates, alkaloids, hydroxycinnamoyl derivatives, carbohydrates and alcohols, phytohormones, coumarins and vitamins. Totally, 118 differential metabolites (VIP ≥ 1, P < 0.05) during chilling treatment process were detected, and their KEGG pathways involved in several metabolic pathways related to dormancy. Sucrose was the most abundant carbohydrate in peony bud. Starch was degradation and Embden Meyerhof Parnas (EMP) activity were increased during the dormancy release process, according to the variations of sugar contents, related enzyme activities and key genes expression. Flavonoids synthesis and accumulation were also promoted by prolonged chilling. Moreover, the variations of phytohormones (salicylic acid, jasmonic acid, abscisic acid, and indole-3-acetic acid) indicated they played different roles in dormancy transition.

CONCLUSION

Our study suggested that starch degradation, EMP activation, and flavonoids accumulation were crucial and associated with bud dormancy transition in tree peony.

Transcriptomic analysis of tea plant (Camellia sinensis) revealed the co-expression network of 4111 paralogous genes and biosynthesis of quality-related key metabolites under multiple stresses

The tea plant is an essential economic plant in many countries. However, its growing season renders them vulnerable to stresses. To understand the transcriptomic influences of these stresses on tea plants, we sequenced and analyzed the transcriptomes under drought, high-temperature, and pest. Paralogs were identified by comparing 14 evolutionarily close genomes. The differentially expressed paralog (DEPs) genes were analyzed regarding single or multiple stresses, and 1075 of the 4111 DEPs were commonly found in all the stresses. The co-expression network of the DEPs and TFs indicated that genes of catechin biosynthesis were associated with most transcription factors specific to each stress. The genes playing a significant role in the late response to drought and pest stress mainly functioned in the early response to high-temperature. This study revealed the relationship between stress and regulation of QRM synthesis and the role of QRMs in response to these (a)biotic stresses.

Transcriptomic and metabolomic profiling reveals the effect of LED light quality on morphological traits, and phenylpropanoid-derived compounds accumulation in Sarcandra glabra seedlings

BACKGROUND

Sarcandra glabra is an evergreen and traditional Chinese herb with anti-oxidant, anti-bacterial, anti-inflammatory, and anti-tumor effects. Light is one of the most influential factor affecting the growth and quality of herbs. In recent times, the introduction of Light Emission Diode (LED) technology has been widely used for plants in greenhouse. However, the impact of such lights on plant growth and the regulatory mechanism of phenylpropanoid-derived compounds in S. glabra remain unclear.

RESULTS

The red LED light (RL) substantially increased the plant height and decreased the stem diameter and leaf area relative to the white LED light (WL), while the blue LED light (BL) significantly reduced the height and leaf area of S. glabra. According to transcriptomic profiling, 861, 378, 47, 10,033, 7917, and 6379 differentially expressed genes (DEGs) were identified among the groups of leaf tissue under BL (BY) vs. leaf tissue under RL (RY), BY vs. leaf tissue under WL (WY), RY vs. WY, root tissue under WL (WG) vs. WY, stem tissue under WL (WJ) vs. WG, and WJ vs. WY, respectively. We identified 46 genes encoding for almost all known enzymes involved in phenylpropanoid biosynthesis, e.g., phenylalanine ammonia lyase (PAL), chalcone synthase (CHS), and flavonol synthase (FLS). We found 53 genes encoding R2R3-MYB proteins and bHLH proteins, respectively, where several were related to flavonoids biosynthesis. A total of 454 metabolites were identified based on metabolomic profiling, of which 44, 87, and 296 compounds were differentially produced in WY vs. RY, WY vs. BY, and WY vs. WG. In BY there was a substantial reduction in the production of esculetin, caffeic acid, isofraxidin, and fraxidin, while the yields of quercitrin and kaempferol were significantly up-regulated. In RY, the contents of cryptochlorogenic acid, cinnamic acid, and kaempferol decreased significantly. Besides, in WG, the production of metabolites (e.g. chlorogenic acid, cryptochlorogenic acid, and scopolin) declined, while their yields increased significantly (e.g. esculetin, fraxetin, isofraxidin, and fraxidin).

CONCLUSION

These results provide further insight into the regulatory mechanism of accumulation patterns of phenylpropanoid-derived compounds in S. glabra under various light conditions, allowing optimum breeding conditions to be developed for this plant.

Dynamic Changes of Phenolic Compounds and Their Associated Gene Expression Profiles Occurring during Fruit Development and Ripening of the Donghong Kiwifruit

The newly released Donghong kiwifruit is a promising commercial cultivar. The dynamic changes of major phenolic compounds (flavonols, flavanols, phenolic acids, and anthocyanins) during the representative stages of fruit development and ripening of the Donghong kiwifruit were determined by high-performance liquid chromatography. The corresponding time-course transcriptional changes were evaluated using the combined analysis of RNA-Seq and qRT-PCR. The most predominant phenolic compound in the Donghong kiwifruit was epicatechin. Cyanidin 3-O-[2-O-(β-xylosyl)-β-galactoside] and cyanidin 3-O-β-galactoside were two essential anthocyanins detected. Candidate genes and pathways involved in phenolic compounds biosynthesis were highlighted. The structural genes (AcLDOX2, Ac5GGT1, and Ac5AT2) and the transcription factor (bHLH74-2) were strongly associated with anthocyanin biosynthesis. AcMYB4-1 may be a novel transcription factor that reduces anthocyanin accumulation. Results from the study may be a very useful supplement to current knowledge of molecular mechanisms to elucidate coloration in the red-fleshed kiwifruit and could help breeders modify the kiwifruit germplasm.

The Regulation of Floral Colour Change in Pleroma raddianum (DC.) Gardner

Floral colour change is a widespread phenomenon in angiosperms, but poorly understood from the genetic and chemical point of view. This article investigates this phenomenon in Pleroma raddianum, a Brazilian endemic species whose flowers change from white to purple. To this end, flavonoid compounds and their biosynthetic gene expression were profiled. By using accurate techniques (Ultra Performance Liquid Chromatography-High-Resolution Mass Spectrometry (UPLC-HRMS)), thirty phenolic compounds were quantified. Five key genes of the flavonoid biosynthetic pathway were partially cloned, sequenced, and the mRNA levels were analysed (RT-qPCR) during flower development. Primary metabolism was also investigated by gas chromatography coupled to mass spectrometry (GC-EIMS), where carbohydrates and organic acids were identified. Collectively, the obtained results suggest that the flower colour change in P. raddianum is determined by petunidin and malvidin whose accumulation coincides with the transcriptional upregulation of early and late biosynthetic genes of the flavonoid pathway, mainly CHS and ANS, respectively. An alteration in sugars, organic acids and phenolic co-pigments is observed together with the colour change. Additionally, an increment in the content of Fe³⁺ ions in the petals, from the pink to purple stage, seemed to influence the saturation of the colour.

Integrated transcriptomic and metabolomic data reveal the flavonoid biosynthesis metabolic pathway in Perilla frutescens (L.) leaves

Perilla frutescens (L.) is an important medicinal and edible plant in China with nutritional and medical uses. The extract from leaves of Perilla frutescens contains flavonoids and volatile oils, which are mainly used in traditional Chinese medicine. In this study, we analyzed the transcriptomic and metabolomic data of the leaves of two Perilla frutescens varieties: JIZI 1 and JIZI 2. A total of 9277 differentially expressed genes and 223 flavonoid metabolites were identified in these varieties. Chrysoeriol, apigenin, malvidin, cyanidin, kaempferol, and their derivatives were abundant in the leaves of Perilla frutescens, which were more than 70% of total flavonoid contents. A total of 77 unigenes encoding 15 enzymes were identified as candidate genes involved in flavonoid biosynthesis in the leaves of Perilla frutescens. High expression of the CHS gene enhances the accumulation of flavonoids in the leaves of Perilla frutescens. Our results provide valuable information on the flavonoid metabolites and candidate genes involved in the flavonoid biosynthesis pathways in the leaves of Perilla frutescens.

FtMYB18 acts as a negative regulator of anthocyanin/proanthocyanidin biosynthesis in Tartary buckwheat

KEY MESSAGE

FtMYB18 plays a role in the repression of anthocyanins and proanthocyanidins accumulation by strongly down-regulating the CHS and DFR genes in Tartary buckwheat, and the C5 motif plays an important role in this process. Anthocyanins and proanthocyanidins (PAs) are important flavonoids in Tartary buckwheat (Fagopyrum tataricum Gaertn.), which provides various vibrant color and stronge abiotic stress resistance. Their synthesis is generally regulated by MYB transcription factors at transcription level. However, the negative regulations of MYB and their effects on flavonol metabolism are poorly understood. A SG4-like MYB subfamily TF, FtMYB18, containing C5 motif was identified from Tartary buckwheat. The expression of FtMYB18 was not only showed a negative correlation with anthocyanins and PAs content but also strongly respond to MeJA and ABA. As far as the transgenic lines with FtMYB18 overexpression, anthocyanins and PAs accumulations were decreased through down-regulating expression levels of NtCHS and NtDFR in tobacco, AtDFR and AtTT12 in Arabidopsis, FtCHS, FtDFR and FtANS in Tartary buckwheat hairy roots, respectively. However, FtMYB18 showed no effect on the FLS gene expression and the metabolites content in flavonol synthesis branch. The further molecular interaction analysis indicated FtMYB18 could mediate the inhibition of anthocyanins and PAs synthesis by forming MBW transcriptional complex with FtTT8 and FtTTG1, or MYB-JAZ complex with FtJAZ1/-3/-4/-7. Importantly, in FtMYB18 mutant lines with C5 motif deletion (FtMYB18-C), both of anthocyanins and PAs accumulations had recovered to the similar level as that in wild type, which was attributed to the weakened MBW complex activity or the deficient molecular interaction between FtMYB18ΔC5 with FtJAZ3/-4. The results showed that FtMYB18 could suppress anthocyanins and PAs synthesis at transcription level through the specific interaction of C5 motif with other proteins in Tartary buckwheat.

Utilization of microRNAs and their regulatory functions for improving biotic stress tolerance in tea plant [Camellia sinensis (L.) O. Kuntze]

MicroRNAs play a central role in responses to biotic stressors through their interactions with their target mRNAs. Tea plant (Camellia sinensis L.), an important beverage crop, is vulnerable to tea geometrid and anthracnose disease that causes considerable crop loss and tea production worldwide. Sustainable production of tea in the current scenario to biotic factors is major challenges. To overcome the problem of biotic stresses, high-throughput sequencing (HTS) with bioinformatics analyses has been used as an effective approach for the identification of stress-responsive miRNAs and their regulatory functions in tea plant. These stress-responsive miRNAs can be utilized for miRNA-mediated gene silencing to enhance stress tolerance in tea plant. Therefore, this review summarizes the current understanding of miRNAs regulatory functions in tea plant responding to Ectropis oblique and Colletotrichum gloeosporioides attacks for future miRNA research. Also, it highlights the utilization of miRNA-mediated gene silencing strategies for developing biotic stress-tolerant tea plant.

Genome-wide interologous interactome map (TeaGPIN) of Camellia sinensis

Tea, prepared from the young leaves of Camellia sinensis, is a non-alcoholic beverage globally consumed due to its antioxidant properties, strong taste and aroma. Although, the genomic data of this medicinally and commercially important plant is available, studies related to its sub-cellular interactomic maps are less explored. In this work, we propose a genome-wide interologous protein-protein interaction (PPI) network of tea, termed as TeaGPIN, consisting of 12,033 nodes and 216,107 interactions, developed using draft genome of tea and known PPIs exhaustively collected from 49 template plants. TeaGPIN interactions are prioritized using domain-domain interactions along with the interolog information. A high-confidence TeaGPIN consisting of 5983 nodes and 58,867 edges is reported and its interactions are further evaluated using protein co-localization similarities. Based on three network centralities (degree, betweenness and eigenvector), 1302 key proteins are reported in tea to have p-value <0.01 by comparing the TeaGPIN with 10,000 realizations of Erdős-Rényi and Barabási-Albert based corresponding random network models. Functional content of TeaGPIN is assessed using KEGG and GO annotations and its modular architecture is explored. Network based characterization is carried-out on the transcription factors, and proteins involved flavonoid biosynthesis and photosynthesis pathways to find novel candidates involved in various regulatory processes. We believe the proposed TeaGPIN will impart useful insights in understanding various mechanisms related to growth and development as well as defence against biotic and abiotic perturbations.

Phenotypic, Phytochemical, and Transcriptomic Analysis of Black Sorghum (Sorghum bicolor L.) Pericarp in Response to Light Quality

Black sorghum [Sorghum bicolor (L.) Moench] is characterized by the black appearance of the pericarp and production of 3-deoxyanthocyanidins (3-DOA), which are valued for their cytotoxicity to cancer cells and as natural food colorants and antioxidant additives. The black pericarp phenotype is not fully penetrant in all environments, which implicates the light spectrum and/or photoperiod as the critical factor for trait expression. In this study, black- or red-pericarp genotypes were grown under regimes of visible light, visible light supplemented with UVA or supplemented with UVA plus UVB (or dark control). Pericarp 3-DOAs and pericarp pigmentation were maximized in the black genotype exposed to a light regime supplemented with UVB. Changes in gene expression during black pericarp development revealed that ultraviolet light activates genes related to plant defense, reactive oxygen species, and secondary metabolism, suggesting that 3-DOA accumulation is associated with activation of flavonoid biosynthesis and several overlapping defense and stress signaling pathways.

Transcriptomic Analysis Reveals Important Roles of Lignin and Flavonoid Biosynthetic Pathways in Rice Thermotolerance During Reproductive Stage

Rice is one of the major staple cereals in the world, but heat stress is increasingly threatening its yield. Analyzing the thermotolerance mechanism from new thermotolerant germplasms is very important for rice improvement. Here, physiological and transcriptome analyses were used to characterize the difference between two germplasms, heat-sensitive MH101 and heat-tolerant SDWG005. Two genotypes exhibited diverse heat responses in pollen viability, pollination characteristics, and antioxidant enzymatic activity in leaves and spikelets. Through cluster analysis, the global transcriptomic changes indicated that the ability of SDWG005 to maintain a steady-state balance of metabolic processes played an important role in thermotolerance. After analyses of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, we found that the thermotolerance mechanism in SDWG00 was associated with reprogramming the cellular activities, such as response to abiotic stress and metabolic reorganization. In contrast, the down-regulated genes in MH101 that appeared to be involved in DNA replication and DNA repair proofreading, could cause serious injury to reproductive development when exposed to high temperature during meiosis. Furthermore, we identified 77 and 11 differentially expressed genes (DEGs) involved in lignin and flavonoids biosynthetic pathways, respectively. Moreover, we found that more lignin deposition and flavonoids accumulation happened in SDWG005 than in MH101 under heat stress. The results indicated that lignin and flavonoid biosynthetic pathways might play important roles in rice heat resistance during meiosis.

De novo transcriptome sequencing of Rhododendron molle and identification of genes involved in the biosynthesis of secondary metabolites

BACKGROUND

Rhododendron molle (Ericaceae) is a traditional Chinese medicinal plant, its flower and root have been widely used to treat rheumatism and relieve pain for thousands of years in China. Chemical studies have revealed that R. molle contains abundant secondary metabolites such as terpenoinds, flavonoids and lignans, some of which have exhibited various bioactivities including antioxidant, hypotension and analgesic activity. In spite of immense pharmaceutical importance, the mechanism underlying the biosynthesis of secondary metabolites remains unknown and the genomic information is unavailable.

RESULTS

To gain molecular insight into this plant, especially on the information of pharmaceutically important secondary metabolites including grayanane diterpenoids, we conducted deep transcriptome sequencing for R. molle flower and root using the Illumina Hiseq platform. In total, 100,603 unigenes were generated through de novo assembly with mean length of 778 bp, 57.1% of these unigenes were annotated in public databases and 17,906 of those unigenes showed significant match in the KEGG database. Unigenes involved in the biosynthesis of secondary metabolites were annotated, including the TPSs and CYPs that were potentially responsible for the biosynthesis of grayanoids. Moreover, 3376 transcription factors and 10,828 simple sequence repeats (SSRs) were also identified. Additionally, we further performed differential gene expression (DEG) analysis of the flower and root transcriptome libraries and identified numerous genes that were specifically expressed or up-regulated in flower.

CONCLUSIONS

To the best of our knowledge, this is the first time to generate and thoroughly analyze the transcriptome data of both R. molle flower and root. This study provided an important genetic resource which will shed light on elucidating various secondary metabolite biosynthetic pathways in R. molle, especially for those with medicinal value and allow for drug development in this plant.

Comprehensive Transcriptomic Analysis for Developing Seeds of a Synthetic Brassica Hexaploid

Polyploidization is a universal phenomenon in plants and plays a crucial role in evolution. In this study, the transcriptomes of developing seeds of a synthetic Brassica hexaploid and its parents (B. rapa and B. carinata) were analyzed to find the gene expression changes in hexaploid seeds. There were 3166 and 3893 DEGs between the Brassica hexaploid and its parents at the full-size stage and mature stage, respectively, most of which were upregulated in hexaploid seeds compared to its parents. At the mature stage, the hexaploid seeds showed a greater difference from its parents. These DEGs had a wide range of functions, which may account for the physiological and morphological differences between the Brassica hexaploid and its parents. The KEGG pathway analysis revealed that hexaploid seeds had higher levels of expression of genes involved in metabolic pathways, RNA transport and biosynthesis of secondary metabolites, and the expression levels in the photosynthesis-related pathways were significantly higher than those in B. rapa. Transgressive expression was the main non-additive expression pattern of the Brassica hexaploid. The gene expression difference between the Brassica hexaploid and its paternal parent was more significant than that with its maternal parent, which may be due in part to the cytoplasmic and maternal effects. Moreover, transcription factor genes, such as G2-like, MYB and mTERF, were highly expressed in hexaploid seeds, possibly promoting their resistance to stress. Our results may provide valuable insights into the adaptation mechanisms of polyploid plants.

PsbHLH1, a novel transcription factor involved in regulating anthocyanin biosynthesis in tree peony (Paeonia suffruticosa)

Flower color is one of the most important features of ornamental plants. Anthocyanin composition and concentration are usually closely related to flower color formation. The biosynthesis of anthocyanin is regulated by a series of structural genes and regulatory genes. The basic helix-loop-helix proteins (bHLHs) are considered as one of the key transcription factors known as the regulators of anthocyanin biosynthesis. However, the bHLH transcription factor family of tree peony (Paeonia suffruticosa) has not been systematically studied in previous studies, especially for the regulation of petal pigmentation. The aim of this study was to identify bHLH genes and unravel their underlying molecular mechanism involved in the regulation of anthocyanin biosynthesis in tree peony. Based on transcriptome profiling analysis, we identified three bHLHs candidate anthocyanin regulators, PsbHLH1, PsbHLH2, and PsbHLH3. PsbHLH1-3 were phylogenetically clustered in the IIIf bHLH subgroup, which is involved in anthocyanin biosynthesis in other plant species. In addition, three bHLH proteins were localized in the nucleus and displayed transcriptional activation activity in a yeast hybrid system. Through a series of functional experiments, we further demonstrated that PsbHLH1 could transcriptionally activate the expression of PsDFR and PsANS via directly binding to their promoters. These results laid a solid foundation to better understand the regulatory mechanisms of anthocyanin biosynthesis in P. suffruticosa and to benefit molecular breeding of tree peony cultivars with novel color.

MYB5-like and bHLH influence flavonoid composition in pomegranate

The fruit of the pomegranate (Punica granatum L.) is an important nutraceutical food rich in polyphenolic compounds, including hydrolysable tannins, anthocyanins and flavonols. Their composition varies according to cultivar, tissue and fruit development stage and is probably regulated by a combination of MYB and bHLH type transcription factors (TFs). In this study, metabolomics analysis during fruit developmental stages in the main pomegranate cultivars, Wonderful and Valenciana with contrasting colour of their ripe fruits, showed that flavonols were mostly present in flowers while catechins were highest in unripe fruits and anthocyanins in late fruit maturation stages. A novel MYB TF, PgMYB5-like, was identified, which differs from previously isolated pomegranate TFs by unique C-terminal protein motifs and lack of the amino-acid residues conserved among anthocyanins promoting MYBs. In both pomegranate cultivars the expression of PgMYB5-like was high at flowering stage, while it decreased during fruit ripening. A previously identified bHLH-type TF, PgbHLH, also showed high transcript levels at flowering stage in both cultivars, while it showed a decrease in expression during fruit ripening in cv. Valenciana, but not in cv. Wonderful. Functional analysis of both TFs was performed by agro-infiltration into Nicotiana benthamiana leaves. Plants infiltrated with the PgMYB5-like+PgbHLH combined construct showed a specific and significant accumulation of intermediates of the flavonoid pathway, especially dihydroflavonols, while anthocyanins were not produced. Thus, we propose a role for PgMYB5-like and PgbHLH in the first steps of flavonoid production in flowers and in unripe fruits. The expression patterns of these two TFs may be key in determining the differential flavonoid composition in both flowers and fruits of the pomegranate varieties Wonderful and Valenciana.

Gibberellins modulate local auxin biosynthesis and polar auxin transport by negatively affecting flavonoid biosynthesis in the root tips of rice

As critical signalling molecules, both gibberellin (GA) and auxin play essential roles in regulating root elongation, and many studies have been shown that auxin influences GA biosynthesis and signalling. However, the mechanism by which GA affects auxin in root elongation is still unknown. In this study, root elongation and DR5-GUS activity were analyzed in rice seedlings. Paclobutrazol-induced short root phenotypes could be partially reversed by co-treatment with IAA, and the inhibition of root elongation caused by naphthylphthalamic acid could be partially reversed when plants were co-treated with GA. DR5-GUS activity was increased in the presence of GA and was reduced at the root tip of paclobutrazol-treated seedlings, indicating that GA could regulate local auxin biosynthesis and polar auxin transport (PAT) in rice root tips. Our RNA-seq analysis showed that GA was involved in the regulation of flavonoid biosynthesis. Flavonoid accumulation level in ks1 root tips was significantly increased and negatively correlated with GA content in GA- and PAC-treated seedlings. GA also rescued the decreased DR5-GUS activity induced by quercetin in rice root tips, confirming that flavonoids act as an intermediary in GA-mediated auxin biosynthesis and PAT. Based on RNA-seq and qPCR analyses, we determined that GA regulates local auxin biosynthesis and polar auxin transport by modulating the expression of OsYUCCA6 and PIN. Our findings provide valuable new insights into the interactions between GA and auxin in the root tips of rice.

Integrated Metabolomics and Transcriptomics Analyses Reveal the Molecular Mechanisms Underlying the Accumulation of Anthocyanins and Other Flavonoids in Cowpea Pod (Vigna unguiculata L.)

As an important vegetable crop of the legume family, cowpea (Vigna unguiculata L.) is grown widely for its tender pod with good taste and nutrition. The purple cowpea pods attract more attention mainly for the eye-catching color and health-promoting ingredients. Initially, large quantities of two major anthocyanins (delphinidin 3-O-glucoside and cyanidin 3-O-glucoside) and nine kinds of flavonoids (most are quercetin-based flavonol glycosides) were separated and identified from purple cowpea pod by ultra-high performance liquid chromatography coupled with quadrupole Orbitrap high-resolution mass spectrometry. To study them systematically, two representative cowpea cultivars with a drastic difference in anthocyanin accumulation were further analyzed by the integration of metabolomics and transcriptomics. A total of 56 differentially accumulated metabolites and 4142 differentially expressed genes were identified, respectively. On the basis of the comprehensive analysis of multiomic data, it was shown that VuMYB90-1, VuMYB90-2, VuMYB90-3, VuCPC, VuMYB4, and endogenous bHLH and WD40 proteins coordinately control anthocyanin and flavonoid accumulation via transcriptional regulation of structural genes in purple cowpea pod.

Metabolome and transcriptome analyses reveal chlorophyll and anthocyanin metabolism pathway associated with cucumber fruit skin color

BACKGROUND

Fruit skin color play important role in commercial value of cucumber, which is mainly determined by the content and composition of chlorophyll and anthocyanins. Therefore, understanding the related genes and metabolomics involved in composition of fruit skin color is essential for cucumber quality and commodity value.

RESULTS

The results showed that chlorophyll a, chlorophyll b and carotenoid content in fruit skin were higher in Lv (dark green skin) than Bai (light green skin) on fruit skin. Cytological observation showed more chloroplast existed in fruit skin cells of Lv. A total of 162 significantly different metabolites were found between the fruit skin of the two genotypes by metabolome analysis, including 40 flavones, 9 flavanones, 8 flavonols, 6 anthocyanins, and other compounds. Crucial anthocyanins and flavonols for fruit skin color, were detected significantly decreased in fruit skin of Bai compared with Lv. By RNA-seq assay, 4516 differentially expressed genes (DEGs) were identified between two cultivars. Further analyses suggested that low expression level of chlorophyll biosynthetic genes, such as chlM, por and NOL caused less chlorophylls or chloroplast in fruit skin of Bai. Meanwhile, a predicted regulatory network of anthocyanin biosynthesis was established to illustrate involving many DEGs, especially 4CL, CHS and UFGT.

CONCLUSIONS

This study uncovered significant differences between two cucumber genotypes with different fruit color using metabolome and RNA-seq analysis. We lay a foundation to understand molecular regulation mechanism on formation of cucumber skin color, by exploring valuable genes, which is helpful for cucumber breeding and improvement on fruit skin color.

Isolation and Analysis of Anthocyanin Pathway Genes from Ribes Genus Reveals MYB Gene with Potent Anthocyanin-Inducing Capabilities

Horticultural crops of the Ribes genus are valued for their anthocyanin-rich fruits, but until now, there were no data about the genes and regulation of their flavonoid pathway. In this study, the coding sequences of flavonoid pathway enzymes and their putative regulators MYB10, bHLH3 and WD40 were isolated, and their expression analyzed in fruits with varying anthocyanin levels from different cultivars of four species belonging to the Ribes genus. Transcription levels of anthocyanin synthesis enzymes and the regulatory gene RrMYB10 correlated with fruit coloration and anthocyanin quantities of different Ribes cultivars. Regulatory genes were tested for the ability to modulate anthocyanin biosynthesis during transient expression in the leaves of two Nicotiana species and to activate Prunus avium promoters of late anthocyanin biosynthesis genes in N. tabacum. Functional tests showed a strong capability of RrMyb10 to induce anthocyanin synthesis in a heterologous system, even without the concurrent expression of any heterologous bHLH, whereas RrbHLH3 enhanced MYB-induced anthocyanin synthesis. Data obtained in this work facilitate further analysis of the anthocyanin synthesis pathway in key Ribes species, and potent anthocyanin inducer RrMyb10 can be used to manipulate anthocyanin expression in heterologous systems.

The identification of an R2R3-MYB transcription factor involved in regulating anthocyanin biosynthesis in Primulina swinglei flowers

Primulina genus is an ideal wild ornamental flower and emerging model for studying biosynthesis, diversity, and evolution of flower pigment. However, the molecular mechanism underlying anthocyanin biosynthesis and regulation in Primulina remains unknown. Here, changes in anthocyanin content and the expression profiles of anthocyanin biosynthetic structural genes were examined in developing Primulina swinglei flowers and three other organs. Seventy-three R2R3-MYB transcription factor genes were identified from transcriptome of P. swinglei flowers, two of which, PsMYB1 and PsMYB2, are candidate regulators of anthocyanin biosynthesis according to clustering analysis. Furthermore, transient over-expression studies using tobacco leaves showed distinct pigment accumulation following co-infection with PsMYB1 and MrbHLH1 (a previously confirmed anthocyanin regulator from Morella rubra). Additionally, dual luciferase assays showed that PsMYB1 trans-activated the PsANS promoter, with the addition of MrbHLH1 resulting in a 5-fold increase in the intensity of this interaction. PsMYB1 did not, however, have any effect on the PsF3H promoter. The expression profile and dual luciferase assays showed that PsMYB2 plays no roles in anthocyanin regulation. Therefore, PsMYB1 is proposed to be the transcription factor gene regulating anthocyanin biosynthesis in P. swinglei.

Alternatively Spliced BnaPAP2.A7 Isoforms Play Opposing Roles in Anthocyanin Biosynthesis of Brassica napus L

Brassica napus L. (rapeseed, oilseed rape, and canola) and varieties of its two diploid parents, B. oleracea and B. rapa, display a large amount of variation in anthocyanin pigmentation of the leaf, stem, and fruit. Here, we demonstrate that BnaPAP2.A7, an ortholog of the B. oleracea anthocyanin activator BoMYB2 that confers purple traits, positively regulates anthocyanin biosynthesis in leaves of B. napus. Sequencing of BnaPAP2.A7 and transgenic analysis suggests that activation of this gene in purple rapeseed may result from a single nucleotide and/or 2bp insertion in its promoter region. BnaPAP2.A7 gives rise to three splice variants, designated BnaPAP2.A7-744, BnaPAP2.A7-910, and BnaPAP2.A7-395 according to the length of the transcripts. While BnaPAP2.A7-744 encodes a full-length R2R3-MYB, both BnaPAP2.A7-910 and BnaPAP2.A7-395 encode truncated proteins that lack both a partial R3 repeat and the complete C terminal domain, and so in vitro are unable to interact with the Arabidopsis bHLH protein AtTT8. Although expression of either BnaPAP2.A7-910 or BnaPAP2.A7-395 in green rapeseed does not result in purple leaves, both genes do modify genome-wide gene expression, with a strong repression of anthocyanin-related genes. We have demonstrated that BnaPAP.A7 regulates anthocyanin accumulation in leaves of B. napus and propose a potential mechanism for modulation of anthocyanin biosynthesis by alternative splicing.

Genomic Characterization and Expression Analysis of Basic Helix-Loop-Helix (bHLH) Family Genes in Traditional Chinese Herb Dendrobium officinale

Dendrobium officinale Kimura et Migo is of great importance as a traditional Chinese herb due to its abundant metabolites. The family of basic helix-loop-helix (bHLH) transcription factors widely exists in plants and plays an essential role in plant growth and development, secondary metabolism as well as responses to environmental changes. However, there is limited information on bHLH genes in D. officinale. In the present study, a total of 98 putative DobHLH genes were identified at the genomic level, which could be classified into 18 clades. Gene structures and conserved motifs in DobHLH genes showed high conservation during their evolution. The conserved amino acids and DNA bindings of DobHLH proteins were predicted, both of which are pivotal for their function. Furthermore, gene expression from eight tissues showed that some DobHLH genes were ubiquitously expressed while other DobHLH genes were expressed in the specific tissues. Expressional changes of DobHLH genes under MeJA and ABA treatments were detected by qRT-PCR. The protein-protein interactions between DobHLHs were predicted and several interactions were confirmed by yeast two hybrid. Therefore, our results here contribute to the understanding of bHLH genes in D. officinale and lay a foundation for the further functional study of its biological processes.

Genome-wide identification and characterization of bHLH family genes from Ginkgo biloba

Basic helix–loop–helix (bHLH) proteins, one of the most important and largest transcription factor family in plants, play important roles in regulating growth and development, stress response. In recent years, many bHLH family genes have been identified and characterized in woody plants. However, a systematic analysis of the bHLH gene family has not been reported in Ginkgo biloba, the oldest relic plant species. In this study, we identifed a total of 85 GbbHLH genes from the genomic and transcriptomic databases of G. biloba, which were classified into 17 subfamilies based on the phylogenetic analysis. Gene structures analysis indicated that the number of exon–intron range in GbbHLHs from 0 to 12. The MEME analysis showed that two conserved motifs, motif 1 and motif 2, distributed in most GbbHLH protein. Subcellular localization analysis exhibited that most GbbHLHs located in nucleus and a few GbbHLHs were distributed in chloroplast, plasma membrane and peroxisome. Promoter cis-element analysis revealed that most of the GbbHLH genes contained abundant cis-elements that involved in plant growth and development, secondary metabolism biosynthesis, various abiotic stresses response. In addition, correlation analysis between gene expression and flavonoid content screened seven candidate GbbHLH genes involved in flavonoid biosynthesis, providing the targeted gene encoding transcript factor for increase the flavonoid production through genetic engineering in G. biloba.

LbNR-Derived Nitric Oxide Delays Lycium Fruit Coloration by Transcriptionally Modifying Flavonoid Biosynthetic Pathway

Anthocyanin-derived fleshy fruit pigmentation has become an excellent system for studying the regulatory network underlying fruit ripening and quality. The transcriptional control of anthocyanin biosynthesis by MYB-bHLH-WDR complexes has been well established, but the intermediate signals through which the environmental or developmental cues regulate these transcription factors remain poorly understood. Here we found that nitric oxide (NO) production during Lycium fruit ripening decreased progressively presenting a negative relationship with anthocyanins. After cloning of the nitric reductase (NR) gene from Lycium barbarum (LbNR) plants, we demonstrated that LbNR-derived NO partially inhibited anthocyanin biosynthesis but enhanced proanthocyanidin (PA) accumulation, and delayed fruit coloration. Application of the NO donor, sodium nitroprusside (SNP), produced a similar effect. The endogenous or exogenous NO downregulated the transcripts both of the regulatory genes and the structural genes that related to anthocyanin biosynthesis, while upregulated both of those genes that related to PA biosynthesis. Given there is a significant negative relationship between the levels of anthocyanins and PAs during Lycium fruit ripening, NO not only inhibited anthocyanin de novo biosynthesis but redirected the flavonoid biosynthetic pathway from anthocyanins to PA production. Two types of LrMYB transcription factors of opposite nature, namely anthocyanin-specific and PA-specific, which belong to the R2R3-MYB subfamily and 1R-MYB subfamily, respectively, were identified from L. ruthenicum fruits. It was further found that NO acts by antagonizing the ABA signaling, a phytohormone we have previously shown playing a positive role in Lycium fruit coloration. Our results provided particularly novel information about NO-ABA-anthocyanin interplay during Lycium fruit development and ripening, which may fill a gap between the developmental cues and the transcriptional regulation of anthocyanin biosynthesis.

Integrative analyses of transcriptomics and metabolomics upon seed germination of foxtail millet in response to salinity

Salinity stress has become an expanding threat to food security worldwide. Revealing the mechanisms of salinity tolerance in plants has immense significance. Foxtail millet (Setaria italica L.) has been regarded as a model crop for exploring mechanisms under stress, considering its extreme adaptation abilities to adverse ecologies. In present study, two foxtail millet cultivars of Yugu2 and An04 with contrasting salt tolerance properties were investigated through integrative analyses of transcriptomics and metabolomics. In the transcriptomics results, 8887 and 12,249 DEGs were identified in Yugu2 and An04 in response to salinity, respectively, and 3149 of which were overlapped between two varieties. These salinity-responsive genes indicated that ion transport, redox homeostasis, phytohormone metabolism, signaling and secondary metabolism were enriched in Yugu2 by GO and KEGG analyses. The integrative omics analysis implied that phenylpropanoid, flavonoid and lignin biosynthesis pathways, and lysophospholipids were vital in determining the foxtail millet salinity tolerance. Importantly, the tolerance of Yugu2 attributed to higher efficiencies of ion channel and antioxidant system. All these provide a comprehensive regulatory network of foxtail millet to cope with salinity, and shed some lights on salt tolerance which is relevant for other cereal crops.

Molecular characterization of an anthocyanin-related glutathione S-transferase gene in Japanese gentian with the CRISPR/Cas9 system

BACKGROUND

The blue pigmentation of Japanese gentian flowers is due to a polyacylated anthocyanin, gentiodelphin, and all associated biosynthesis genes and several regulatory genes have been cloned and characterized. However, the final step involving the accumulation of anthocyanins in petal vacuoles remains unclear. We cloned and analyzed the glutathione S-transferases (GSTs) in Japanese gentian that are known to be involved in anthocyanin transport in other plant species.

RESULTS

We cloned GST1, which is expressed in gentian flower petals. Additionally, this gene belongs to the Phi-type GST clade related to anthocyanin biosynthesis. We used the CRISPR/Cas9-mediated genome editing system to generate loss-of-function GST1 alleles. The edited alleles were confirmed by Sanger and next-generation sequencing analyses. The GST1 genome-edited lines exhibited two types of mutant flower phenotypes, severe (almost white) and mild (pale blue). The phenotypes were associated with decreased anthocyanin accumulation in flower petals. In the GST1 genome-edited lines, sugar-induced stress conditions inhibited the accumulation of anthocyanins in stems and leaves, suggestvhing that GST1 is necessary for stress-related anthocyanin accumulation in organs other than flowers. These observations clearly demonstrate that GST1 is the gene responsible for anthocyanin transport in Japanese gentian, and is necessary for the accumulation of gentiodelphin in flowers.

CONCLUSIONS

In this study, an anthocyanin-related GST gene in Japanese gentian was functionally characterized. Unlike other biosynthesis genes, the functions of GST genes are difficult to examine in in vitro studies. Thus, the genome-editing strategy described herein may be useful for in vivo investigations of the roles of transport-related genes in gentian plants.

Methylation of MdMYB1 locus mediated by RdDM pathway regulates anthocyanin biosynthesis in apple

Methylation at the MdMYB1 promoter in apple sports has been reported as a regulator of the anthocyanin pathway, but little is known about how the locus is recognized by the methylation machinery to regulate anthocyanin accumulation. In this study, we analysed three differently coloured 'Fuji' apples and found that differences in the transcript levels of MdMYB1, which encodes a key regulator of anthocyanin biosynthesis, control the anthocyanin content (and therefore colour) in fruit skin. The CHH methylation levels in the MR3 region (-1246 to -780) of the MdMYB1 promoter were found to be negatively correlated with MdMYB1 expression. Thus, they were ideal materials to study DNA methylation in apple sports. The protein of RNA-directed DNA methylation (RdDM) pathway responsible for CHH methylation, MdAGO4, was found to interact with the MdMYB1 promoter. MdAGO4s can interact with MdRDM1 and MdDRM2s to form an effector complex, fulfilling CHH methylation. When MdAGO4s and MdDRM2s were overexpressed in apple calli and Arabidopsis mutants, those proteins increase the CHH methylation of AGO4-binding sites. In electrophoretic mobility shift assays, MdAGO4s were found to specifically bind to sequence containing ATATCAGA. Knockdown of MdNRPE1 did not affect the binding of MdAGO4s to the c3 region of the MdMYB1 promoter in 35S::AGO4 calli. Taken together, our data show that the MdMYB1 locus is methylated through binding of MdAGO4s to the MdMYB1 promoter to regulate anthocyanin biosynthesis by the RdDM pathway.

A minor QTL, SG3, encoding an R2R3-MYB protein, negatively controls grain length in rice

SG3, an R2R3 MYB protein coding gene that tightly linked to a major QTLGS3, negatively regulates grain length while dependent on the status ofGS3in rice. It is often very difficult to isolate a minor QTL that is closely linked to a major QTL in rice. In this study, we focused on the isolation of a minor grain length QTL, small grain 3 (SG3), which is closely linked to the major QTL grain size 3 (GS3). The genetic effect of SG3 on grain length was dependent on GS3 status. Its genetic effect was larger in the presence of nonfunctional sg3 than functional SG3. A large number of near-isogenic F2 plants in which GS3 was fixed with nonfunctional alleles were developed to clone SG3. A gene encoding an R2R3 MYB domain transcriptional regulator was identified as the candidate gene for SG3. SG3 overexpression and knockdown plants showed shortened and elongated grains, respectively, which demonstrated that SG3 acts as a negative regulator of grain length. SG3 was preferentially expressed in panicles after flowering, and SG3 acted as a transcription activator. Comparative sequencing analysis identified a 12-bp insertion in the third exon of NYZ that led to a frameshift and resulted in a premature stop codon. The insertion/deletion was associated with grain length in the presence of functional GS3 in the indica subspecies. SG3 and GS3 were frequently in coupling phase in indica rice, making them good targets for the breeding of cultivars with short or long grains. The isolation of the SG3 MYB gene provides new gene resource and contributes to the regulatory network of grain length in rice.

Outstanding questions in flower metabolism

The great diversity of flowers, their color, odor, taste, and shape, is mostly a result of the metabolic processes that occur in this reproductive organ when the flower and its tissues develop, grow, and finally die. Some of these metabolites serve to advertise flowers to animal pollinators, other confer protection towards abiotic stresses, and a large proportion of the molecules of the central metabolic pathways have bioenergetic and signaling functions that support growth and the transition to fruits and seeds. Although recent studies have advanced our general understanding of flower metabolism, several questions still await an answer. Here, we have compiled a list of open questions on flower metabolism encompassing molecular aspects, as well as topics of relevance for agriculture and the ecosystem. These questions include the study of flower metabolism through development, the biochemistry of nectar and its relevance to promoting plant-pollinator interaction, recycling of metabolic resources after flowers whiter and die, as well as the manipulation of flower metabolism by pathogens. We hope with this review to stimulate discussion on the topic of flower metabolism and set a reference point to return to in the future when assessing progress in the field.

MiR156 regulates anthocyanin biosynthesis through SPL targets and other microRNAs in poplar

Anthocyanins biosynthesized from the flavonoid pathway are types of pigments that are involved in the protection of poplar from biotic and abiotic stresses. Previous researchers studying anthocyanin-related transcription factors and structural genes in poplar have made significant discoveries. However, little is known about the regulatory role of microRNAs in anthocyanin biosynthesis in poplar. Here, we overexpressed miR156 in poplar to study the comprehensive effects of the miR156-SPL module on the biosynthesis of anthocyanins. Small RNA sequencing analysis revealed 228 microRNAs differentially expressed in transgenic poplar plants with dramatically increased miR156 levels. Furthermore, integrated microRNAomic and transcriptomic analysis suggested that two microRNAs, miR160h, and miR858, have the potential to affect anthocyanin accumulation in poplar by regulating auxin response factors and MYB transcription factors, respectively. Additionally, the accumulation of miR160h and miR858 displayed a positive correlation with miR156 levels, suggesting a possible interaction between the miR156-SPL module and these microRNAs in poplar. Last, metabolomics analysis revealed that the levels of anthocyanins, flavones, and flavonols were substantially elevated in transgenic poplar plants overexpressing miR156 compared with the wild type, whereas the total lignin content was reduced in the transgenic plants. Taken together, our results indicate that miR156 can fine tune the anthocyanin biosynthetic pathway via multiple factors, including microRNAs, transcription factors, and the levels of structural genes, in poplar. This provides additional clues for understanding the complex regulatory network of anthocyanin biosynthesis in woody plants.

MYB repressors and MBW activation complex collaborate to fine-tune flower coloration in Freesia hybrida

Floral anthocyanin has multiple ecological and economic values, its biosynthesis largely depends on the conserved MYB-bHLH-WD40 (MBW) activation complex and MYB repressors hierarchically with the MBW complex. In contrast to eudicots, the MBW regulatory network model has not been addressed in monocots because of the lack of a suitable system, as grass plants exhibit monotonous floral pigmentation patterns. Presently, the MBW regulatory network was investigated in a non-grass monocot plant, Freesia hybrida. FhMYB27 and FhMYBx with different functional manners were confirmed to be anthocyanin related R2R3 and R3 MYB repressors, respectively. Particularly, FhMYBx could obstruct the formation of positive MBW complex by titrating bHLH proteins, whereas FhMYB27 mainly defected the activator complex into suppressor via its repression domains in C-terminus. Furthermore, the hierarchical and feedback regulatory loop was verified, indicating the synergistic and sophisticated regulatory network underlying Freesia anthocyanin biosynthesis was quite similar to that reported in eudicot plants.

Yueqing Li, Xiaotong Shan, et al. study the MYB-bHLH-WD40 (MBW) regulatory network in a non-grass monocot plant, Freesia hybrida. They report two anthocyanin related MYB repressors FhMYB27 and FhMYBx and verified their involvement in a functional feedback loop with MBW to regulate anthocyanin biosynthesis.

Genome assembly of wild tea tree DASZ reveals pedigree and selection history of tea varieties

Wild teas are valuable genetic resources for studying domestication and breeding. Here we report the assembly of a high-quality chromosome-scale reference genome for an ancient tea tree. The further RNA sequencing of 217 diverse tea accessions clarifies the pedigree of tea cultivars and reveals key contributors in the breeding of Chinese tea. Candidate genes associated with flavonoid biosynthesis are identified by genome-wide association study. Specifically, diverse allelic function of CsANR, CsF3'5'H and CsMYB5 is verified by transient overexpression and enzymatic assays, providing comprehensive insights into the biosynthesis of catechins, the most important bioactive compounds in tea plants. The inconspicuous differentiation between ancient trees and cultivars at both genetic and metabolic levels implies that tea may not have undergone long-term artificial directional selection in terms of flavor-related metabolites. These genomic resources provide evolutionary insight into tea plants and lay the foundation for better understanding the biosynthesis of beneficial natural compounds.

Full-length sequencing of ginkgo transcriptomes for an in-depth understanding of flavonoid and terpenoid trilactone biosynthesis

Ginkgo biloba L. is regarded as the most ancient living tree, and its kernel has been used as a traditional Chinese medicine for more than 2,000 years. The leaf extracts of this tree have been among the bestselling herbal remedies in Western countries since the last century. To understand the biosynthesis of the pharmacologically active ingredients in G. biloba, flavonoids and terpenoid trilactones (TTLs), we sequenced the transcriptomes of G. biloba leaves, kernels and testae with Iso-Seq and RNA-Seq technologies and obtained 152,524 clean consensus reads. When these reads were used to improve the annotation of the G. biloba genome, 4,856 novel genes, 25,583 new isoforms of previously annotated genes and 4,363 lncRNAs were discovered. Gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that genes involved in growth, regulation and response to stress were more likely to be regulated by alternative splicing (AS) or alternative polyadenylation (APA), which represent the two most important posttranscriptional regulation mechanisms. It was found that some of the characterized genes involved in the biosynthesis of flavonoids and TTLs were also possibly regulated by AS and APA. Using phylogenetic and gene expression pattern analyses, some candidate genes for the biosynthesis of flavonoids and TTLs were screened. After qRT-PCR validation, the final candidate genes for flavonoid biosynthesis included three UDP-glycosyltransferases and one MYB transcription factor, while the candidate genes for TTL biosynthesis included two cytochrome P450 and one WRKY transcription factor. Our study suggested that Iso-Seq may play an important role in improving genome annotation, elucidating AS and APA mechanisms and discovering candidate genes involved in the biosynthesis of some secondary metabolites.

TRANSPARENT TESTA GLABRA1, a Key Regulator in Plants with Multiple Roles and Multiple Function Mechanisms

TRANSPARENT TESTA GLABRA1 (TTG1) is a WD40 repeat protein. The phenotypes caused by loss-of-function of TTG1 were observed about half a century ago, but the TTG1 gene was identified only about twenty years ago. Since then, TTG1 has been found to be a plant-specific regulator with multiple roles and multiple functional mechanisms. TTG1 is involved in the regulation of cell fate determination, secondary metabolisms, accumulation of seed storage reserves, plant responses to biotic and abiotic stresses, and flowering time in plants. In some processes, TTG1 may directly or indirectly regulate the expression of downstream target genes via forming transcription activator complexes with R2R3 MYB and bHLH transcription factors. Whereas in other processes, TTG1 may function alone or interact with other proteins to regulate downstream target genes. On the other hand, the studies on the regulation of TTG1 are very limited. So far, only the B3-domain family transcription factor FUSCA3 (FUS3) has been found to regulate the expression of TTG1, phosphorylation of TTG1 affects its interaction with bHLH transcription factor TT2, and TTG1 proteins can be targeted for degradation by the 26S proteasome. Here, we provide an overview of TTG1, including the identification of TTG1, the functions of TTG1, the possible function mechanisms of TTG1, and the regulation of TTG1. We also proposed potential research directions that may shed new light on the regulation and functional mechanisms of TTG1 in plants.

MYB21 interacts with MYC2 to control the expression of terpene synthase genes in flowers of Freesia hybrida and Arabidopsis thaliana

Previously, linalool was found to be the most abundant component among the cocktail of volatiles released from flowers of Freesia hybrida. Linalool formation is catalysed by monoterpene synthase TPS1. However, the regulatory network developmentally modulating the expression of the TPS1 gene in Freesia hybrida remains unexplored. In this study, three regulatory genes, FhMYB21L1, FhMYB21L2, and FhMYC2, were screened from 52 candidates. Two MYB transcription factor genes were synchronously expressed with FhTPS1 and could activate its expression significantly when overexpressed, and the binding of FhMYB21L2 to the MYBCORE sites in the FhTPS1 promoter was further confirmed, indicating a direct role in activation. FhMYC2 showed an inverse expression pattern compared with FhTPS1; its expression led to a decreased binding of FhMYB21 to the FhTPS1 promoter to reduce its activation capacity when co-expressed, suggesting a role for an MYB-bHLH complex in the regulation of the FhTPS1 gene. In Arabidopsis, both MYB21 and MYC2 regulators were shown to activate the expression of sesquiterpene synthase genes, and the regulatory roles of AtMYB21 and AtMYC2 in the expression of the linalool synthase gene were also confirmed, implying conserved functions of the MYB-bHLH complex in these two evolutionarily divergent plants. Moreover, the expression ratio between MYB21 and MYC2 orthologues might be a determinant factor in floral linalool emission.

Transcriptome Sequencing and Chemical Analysis Reveal the Formation Mechanism of White Florets in Carthamus tinctorius L

Carthamus tinctorius L. (safflower), an economic crop and herb, has been extensively studied for its diverse chemical constituents and pharmacological effects, but the mechanism of safflower pigments (SP) leading to different colors of florets has not been clarified. In the present study, we compared the contents of SP in two varieties of safflower with white and red florets, named Xinhonghua No. 7 (WXHH) and Yunhong No. 2 (RYH). The results showed the contents of SP in RYH were higher than WXHH. To investigate genes related to SP, we obtained six cDNA libraries of florets from the two varieties by transcriptome sequencing. A total of 225,008 unigenes were assembled and 40 unigenes related to safflower pigment biosynthesis were annotated, including 7 unigenes of phenylalanine ammonia-lyase (PAL), 20 unigenes of 4-coumarate-CoA ligase (4CL), 1 unigene of trans-cinnamate 4-monooxygenase (C4H), 7 unigenes of chalcone synthase (CHS), 4 unigenes of chalcone isomerase (CHI), and 1 unigene of flavanone 3-hydroxylase (F3H). Based on expression levels we selected 16 differentially expressed unigenes (DEGs) and tested them using reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR), which was consistent with the sequencing results. Consequently, we speculated that in WXHH, 3 PALs, 3 4CLs, 1 C4H, 1 CHS, and 1 CHI, which were down-regulated, and 1 F3H, which was up-regulated, may play a key role in the formation of white florets.

Specialized phenolic compounds in seeds: structures, functions, and regulations

Plants produce a huge diversity of specialized metabolites (SM) throughout their life cycle that play important physiological and ecological functions. SM can protect plants and seeds against diseases, predators, and abiotic stresses, or support their interactions with beneficial or symbiotic organisms. They also have strong impacts on human nutrition and health. Despite this importance, the biosynthesis and biological functions of most of the SM remain elusive and their diversity and/or quantity have been reduced in most crops during domestication. Seeds present a large number of SM that are important for their physiological, agronomic, nutritional or industrial qualities and hence, provide interesting models for both studying biosynthesis and producing large amounts of specialized metabolites. For instance, phenolics are abundant and widely distributed in seeds. More specifically, flavonoid pathway has been instrumental for understanding environmental or developmental regulations of specialized metabolic pathways, at the molecular and cellular levels. Here, we summarize current knowledge on seed phenolics as model, and discuss how recent progresses in omics approaches could help to further characterize their diversity, regulations, and the underlying molecular mechanisms involved.

Isolation and identification of MYB transcription factors (MYB19Long and MYB19Short) involved in raised spot anthocyanin pigmentation in lilies (Lilium spp.)

Although anthocyanin color patterns on flowers are among the most attractive characteristics, the genetic mechanisms through which color patterns are developed are not well understood, especially for color patterns associated with altered petal structure. Lilium species and cultivars often develop raised spots, where the interior surfaces of tepals increase to develop bumps with accompanying anthocyanin accumulation. The aim of this study was to identify transcription factors regulating pigmentation of the bumps. We identified two R2R3-MYB genes, MYB19Long and MYB19Short, in Lilium leichtlinii, L. lancifolium, and Asiatic hybrid lily cultivars. Their amino acid sequences were similar; however, part of the C-terminal region was triplicated in MYB19Long. Spatial and temporal expression profiles in lilies were strongly associated with anthocyanin biosynthesis gene expression in the bumps, and some defects were found in these genes in L. lancifolium 'Pure Gold' that developed colorless bumps. Thus, both MYB19Long and MYB19Short were likely to be involved in the bump pigmentation. MYB19Long had a stronger ability to stimulate target gene expression than MYB19Short, and expression levels of MYB19Long were greater than those of MYB19Short in lily tepals; thus, the ability to biosynthesize anthocyanin pigments was greater for MYB19Long than for MYB19Short. Among the F₁ population, MYB19Short expression was found only in the tepals of F₁ plants that developed bumps, although all of the F₁ plants possessed the MYB19Short gene, indicating that MYB19 expression followed bump development. These findings helped to elucidate the genetic mechanisms underlying raised spot development.