Coexpression network analysis reveals an MYB transcriptional activator involved in capsaicinoid biosynthesis in hot peppers

Two telomere-to-telomere gapless genomes reveal insights into Capsicum evolution and capsaicinoid biosynthesis

Chili pepper (Capsicum) is known for its unique fruit pungency due to the presence of capsaicinoids. The evolutionary history of capsaicinoid biosynthesis and the mechanism of their tissue specificity remain obscure due to the lack of high-quality Capsicum genomes. Here, we report two telomere-to-telomere (T2T) gap-free genomes of C. annuum and its wild nonpungent relative C. rhomboideum to investigate the evolution of fruit pungency in chili peppers. We precisely delineate Capsicum centromeres, which lack high-copy tandem repeats but are extensively invaded by CRM retrotransposons. Through phylogenomic analyses, we estimate the evolutionary timing of capsaicinoid biosynthesis. We reveal disrupted coding and regulatory regions of key biosynthesis genes in nonpungent species. We also find conserved placenta-specific accessible chromatin regions, which likely allow for tissue-specific biosynthetic gene coregulation and capsaicinoid accumulation. These T2T genomic resources will accelerate chili pepper genetic improvement and help to understand Capsicum genome evolution.

Haplotype analysis, regulatory elements and docking simulation of structural models of different AT3 copies in the genus Capsicum

Capsaicinoids are responsible for the pungency in Capsicum species. These are synthesized by the Capsaicin synthase (CS) encoded by the AT3 gene, which catalyzes the transference of an acyl moiety from a branched-chain fatty acid-CoA ester to the vanillylamine to produce capsaicinoids. Some AT3 gene copies have been identified on the Capsicum genome. The absence of capsaicinoid in some nonpungent accessions is related to mutant AT3 alleles. The differences between CS protein copies can affect the tridimensional structure of the protein and the affinity for its substrates, and this could affect fruit pungency. This study characterized 32 AT3 sequences covering Capsicum pungent and non-pungent accessions. These were clustered in AT3-D1 and AT3-D2 groups and representative sequences were analyzed. Genomic upstream analysis shows different regulatory elements, mainly responsive to light and abiotic stress. AT3-D1 and AT3-D2 gene expression was confirmed in fruit tissues of C. annuum. Amino acid substitutions close to the predictable HXXXD and DFGWG motifs were also identified. AT3 sequences were modeled showing a BAHD acyltransferase structure with two connected domains. A pocket with different shape, size and composition between AT3 models was found inside the protein, with the conserved motif HXXXD exposed to it, and a channel for their accessibility. CS substrates exhibit high interaction energies with the His and Asp conserved residues. AT3 models have different interaction affinities with the (E)-8-methylnon-6-enoyl-CoA, 8-methylnonanoyl-CoA and vanillylamine substrates. These results suggested that AT3-D1 and AT3-D2 sequences encode CS enzymes with different regulatory factors and substratum affinities.Communicated by Ramaswamy H. Sarma.

Biochemical Aspects of Putative Aminotransferase Responsible for Converting Vanillin to Vanillylamine in the Capsaicinoid Biosynthesis Pathway in Capsicum Plants

The biosynthesis pathway of capsaicinoids includes the conversion of vanillin to vanillylamine, where putative aminotransferase (pAMT) is thought to be the enzyme responsible in Capsicum plants. The objectives of this study were to prove that pAMT is the enzyme responsible for this conversion in plants and to clarify its catalytic properties using biochemical methods. Both an extract of habanero placenta and recombinant pAMT (rpAMT) constructed by using an Escherichia coli expression system were able to convert vanillin to vanillylamine in the presence of γ-aminobutyric acid as an amino donor and pyridoxal phosphate as a coenzyme. Conversion from vanillin to vanillylamine by the placenta extract was significantly attenuated by adding an anti-pAMT antibody to the reaction system. The amino donor specificity and affinity for vanillin of rpAMT were similar to those of the placenta extract. We thus confirmed that pAMT is the enzyme responsible for the conversion of vanillin to vanillylamine in capsaicinoid synthesis in Capsicum fruits. Therefore, we propose that pAMT should henceforth be named vanillin aminotransferase (VAMT).

Systematic analysis of the UDP-glucosyltransferase family: discovery of a member involved in rutin biosynthesis in Solanum melongena

Eggplant (Solanum melongena) is an economically important crop and rich in various nutrients, among which rutin that has positive effects on human health is found in eggplant. Glycosylation mediated by UDP-glycosyltransferases (UGTs) is a key step in rutin biosynthesis. However, the UGT gene has not been reported in eggplant to date. Herein, 195 putative UGT genes were identified in eggplant by genome-wide analysis, and they were divided into 17 subgroups (Group A-P and Group R) according to the phylogenetic evolutionary tree. The members of Groups A, B, D, E and L were related to flavonol biosynthesis, and rutin was the typical flavonol. The expression profile showed that the transcriptional levels of SmUGT genes in Clusters 7-10 were closely related to those of rutin biosynthetic pathway genes. Notably, SmUGT89B2 was classified into Cluster 7 and Group B; its expression was consistent with rutin accumulation in different tissues and different leaf stages of eggplant. SmUGT89B2 was located in the nucleus and cell membrane. Virus-induced gene silencing (VIGS) and transient overexpression assays showed that SmUGT89B2 can promote rutin accumulation in eggplant. These findings provide new insights into the UGT genes in eggplant, indicating that SmUGT89B2 is likely to encode the final enzyme in rutin biosynthesis.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

Stressing the importance of plant specialized metabolites: omics-based approaches for discovering specialized metabolism in plant stress responses

Plants produce a diverse range of specialized metabolites that play pivotal roles in mediating environmental interactions and stress adaptation. These unique chemical compounds also hold significant agricultural, medicinal, and industrial values. Despite the expanding knowledge of their functions in plant stress interactions, understanding the intricate biosynthetic pathways of these natural products remains challenging due to gene and pathway redundancy, multifunctionality of proteins, and the activity of enzymes with broad substrate specificity. In the past decade, substantial progress in genomics, transcriptomics, metabolomics, and proteomics has made the exploration of plant specialized metabolism more feasible than ever before. Notably, recent advances in integrative multi-omics and computational approaches, along with other technologies, are accelerating the discovery of plant specialized metabolism. In this review, we present a summary of the recent progress in the discovery of plant stress-related specialized metabolites. Emphasis is placed on the application of advanced omics-based approaches and other techniques in studying plant stress-related specialized metabolism. Additionally, we discuss the high-throughput methods for gene functional characterization. These advances hold great promise for harnessing the potential of specialized metabolites to enhance plant stress resilience in the future.

CaHMA1 promotes Cd accumulation in pepper fruit

The main planting areas for pepper (Capsicum sp.) are high in cadmium (Cd), which is the most prevalent heavy metal pollutant worldwide. Breeding pepper cultivars with low Cd levels can promote sustainable agricultural production and ensure the safety of pepper products. To identify breeding targets for reducing Cd accumulation in pepper fruits, we performed a genome-wide association study on 186 accessions. Polymorphisms were associated with fruit Cd content in a genomic region containing a homolog of Arabidopsis (Arabidopsis thaliana) Heavy metal-transporting ATPase 1 (HMA1) encoding a P-type ATPase. In two cultivars with contrasting Cd accumulation, transcriptome analysis revealed differentially expressed genes enriched for carbohydrate metabolism and photosynthesis in fruits with high Cd accumulation, and a Cd2+/Zn2+-exporting ATPase gene (HMA). Heterologous expression of CaHMA1 in yeast increases Cd sensitivity. Overexpression of CaHMA1 conferred a severe increase in Cd content in Arabidopsis plants, whereas reduced CaHMA1 expression in pepper fruits decreased Cd content. We propose that CaHMA1 expression may be an important component of the high Cd accumulation in pepper plants.

Virus-Induced Gene Silencing in the Tea Plant (Camellia sinensis)

The recent availability of a number of tea plant genomes has sparked substantial interest in using reverse genetics to explore gene function in tea (Camellia sinensis). However, a hurdle to this is the absence of an efficient transformation system, and virus-induced gene silencing (VIGS), a transient transformation system, could be an optimal choice for validating gene function in the tea plant. In this study, phytoene desaturase (PDS), a carotenoid biosynthesis gene, was used as a reporter to evaluate the VIGS system. The injection sites of the leaves (leaf back, petiole, and stem) for infiltration were tested, and the results showed that petiole injection had the most effective injection, without leading to necrotic lesions that cause the leaves to drop. Tea leaves were inoculated with Agrobacterium harboring a tobacco rattle virus plasmid (pTRV2) containing a CsPDS silencing fragment. The tea leaves exhibited chlorosis symptoms 7-14 days after inoculation, depending on the cultivar. In the chlorosis plants, the coat protein (CP) of tobacco rattle virus (TRV) was detected and coincided with the lower transcription of CsPDS and reduced chlorophyll content compared with the empty vector control, with 81.82% and 54.55% silencing efficiency of 'LTDC' and 'YSX', respectively. These results indicate that the VIGS system with petiole injection could quickly and effectively silence a gene in tea plants.

ERF5 enhances protocorm-like body regeneration via enhancement of STM expression in Dendrobium orchid

Orchid plants develop protocorms upon germination and produce protocorm-like structures called protocorm-like bodies (PLBs) from protocorms and somatic cells via tissue culture. Protocorm-like bodies have broad technical application potential in the orchid industry and their regeneration is a distinct developmental process in the plant kingdom. However, little is known about this unparalleled developmental program. In this study, we identified a PLB-abundant gene, ethylene response factor (ERF), and a transcription factor named DoERF5, and determined its important role in PLB regeneration in Dendrobium orchid. Overexpression of DoERF5 in Dendrobium greatly enhanced the PLB regeneration from PLB and stem explants, and upregulated the expression of WOUND-INDUCED DEDIFFERENTIATION (DoWIND) homologs and SHOOT MERISTEMLESS (DoSTM), as well as the genes involved in cytokinin biosynthesis (DoIPT) and the cytokinin response factors (DoARRs). However, silencing DoERF5 reduced the regeneration rate of PLBs, and downregulated the expression of DoWIND homologs, DoSTM and DoARRs. We demonstrated that DoERF5 is directly bound to the DoSTM promoter and regulates its expression. In addition, overexpression of DoSTM in Dendrobium orchid resulted in favorable regeneration of PLBs. Our results clarify that DoERF5 regulates the regeneration of PLB by enhancing DoSTM expression. Our findings provide new insights into how DoERF5 mediates PLB regeneration and offers technical potential in improving clonal propagation, preservation, and the bioengineering of orchids.

An R-R-type MYB transcription factor promotes non-climacteric pepper fruit carotenoid pigment biosynthesis

Carotenoids are major accessory pigments in the chloroplast, and they also act as phytohormones and volatile compound precursors to influence plant development and confer characteristic colours, affecting both the aesthetic and nutritional value of fruits. Carotenoid pigmentation in ripening fruits is highly dependent on developmental trajectories. Transcription factors incorporate developmental and phytohormone signalling to regulate the biosynthesis process. By contrast to the well-established pathways regulating ripening-related carotenoid biosynthesis in climacteric fruit, carotenoid regulation in non-climacteric fruit is poorly understood. Capsanthin is the primary carotenoid of non-climacteric pepper (Capsicum) fruit; its biosynthesis is tightly associated with fruit ripening, and it confers red pigmentation to the ripening fruit. In the present study, using a coexpression analysis, we identified an R-R-type MYB transcription factor, DIVARICATA1, and demonstrated its role in capsanthin biosynthesis. DIVARICATA1 encodes a nucleus-localised protein that functions primarily as a transcriptional activator. Functional analyses showed that DIVARICATA1 positively regulates carotenoid biosynthetic gene (CBG) transcript levels and capsanthin levels by directly binding to and activating CBG promoter transcription. Furthermore, an association analysis revealed a significant positive association between DIVARICATA1 transcription level and capsanthin content. ABA promotes capsanthin biosynthesis in a DIVARICATA1-dependent manner. Comparative transcriptomic analysis of DIVARICATA1 in Solanaceae plants showed that its function likely differs among species. Moreover, the pepper DIVARICATA1 gene could be regulated by the ripening regulator MADS-RIN. The present study illustrates the transcriptional regulation of capsanthin biosynthesis and offers a target for breeding peppers with high red colour intensity.

Single-cell transcriptome atlas reveals spatiotemporal developmental trajectories in the basal roots of moso bamboo (Phyllostachys edulis)

Roots are essential for plant growth and development. Bamboo is a large Poaceae perennial with 1642 species worldwide. However, little is known about the transcriptional atlas that underpins root cell-type differentiation. Here, we set up a modified protocol for protoplast preparation and report single-cell transcriptomes of 14 279 filtered single cells derived from the basal root tips of moso bamboo. We identified four cell types and defined new cell-type-specific marker genes for the basal root. We reconstructed the developmental trajectories of the root cap, epidermis, and ground tissues and elucidated critical factors regulating cell fate determination. According to in situ hybridization and pseudotime trajectory analysis, the root cap and epidermis originated from a common initial cell lineage, revealing the particularity of bamboo basal root development. We further identified key regulatory factors for the differentiation of these cells and indicated divergent root developmental pathways between moso bamboo and rice. Additionally, PheWOX13a and PheWOX13b ectopically expressed in Arabidopsis inhibited primary root and lateral root growth and regulated the growth and development of the root cap, which was different from WOX13 orthologs in Arabidopsis. Taken together, our results offer an important resource for investigating the mechanism of root cell differentiation and root system architecture in perennial woody species of Bambusoideae.

Genome-Wide Identification of WRKY Gene Family and Functional Characterization of CcWRKY25 in Capsicum chinense

Pepper is renowned worldwide for its distinctive spicy flavor. While the gene expression characteristics of the capsaicinoid biosynthesis pathway have been extensively studied, there are already a few reports regarding transcriptional regulation in capsaicin biosynthesis. In this study, 73 WRKYs were identified in the genome of Capsicum chinense, and their physicochemical traits, DNA, and protein sequence characteristics were found to be complex. Combining RNA-seq and qRT-PCR data, the WRKY transcription factor CA06g13580, which was associated with the accumulation tendency of capsaicinoid, was screened and named CcWRKY25. CcWRKY25 was highly expressed in the placenta of spicy peppers. The heterologous expression of CcWRKY25 in Arabidopsis promoted the expression of genes PAL, 4CL1, 4CL2, 4CL3, CCR, and CCoAOMT and led to the accumulation of lignin and flavonoids. Furthermore, the expression of the capsaicinoid biosynthesis pathway genes (CBGs) pAMT, AT3, and KAS was significantly reduced in CcWRKY25-silenced pepper plants, resulting in a decrease in the amount of capsaicin. However, there was no noticeable difference in lignin accumulation. The findings suggested that CcWRKY25 could be involved in regulating capsaicinoid synthesis by promoting the expression of genes upstream of the phenylpropanoid pathway and inhibiting CBGs' expression. Moreover, the results highlighted the role of CcWRKY25 in controlling the pungency of pepper and suggested that the competitive relationship between lignin and capsaicin could also regulate the spiciness of the pepper.

Comparative transcriptome analysis reveals the involvement of an MYB transcriptional activator, SmMYB108, in anther dehiscence in eggplant

Male sterility is a highly attractive agronomic trait as it effectively prevents self-fertilization and facilitates the production of high-quality hybrid seeds in plants. Timely release of mature pollen following anther dehiscence is essential for stamen development in flowering plants. Although several theories have been proposed regarding this, the specific mechanism of anther development in eggplant remains elusive. In this study, we selected an R2R3-MYB transcription factor gene, SmMYB108, that encodes a protein localized primarily in the nucleus by comparing the transcriptomics of different floral bud developmental stages of the eggplant fertile line, F142. Quantitative reverse transcription polymerase chain reaction revealed that SmMYB108 was preferentially expressed in flowers, and its expression increased significantly on the day of flowering. Overexpression of SmMYB108 in tobacco caused anther dehiscence. In addition, we found that SmMYB108 primarily functions as a transcriptional activator via C-terminal activation (amino acid 262-317). Yeast one-hybrid and dual-luciferase reporter assays revealed that genes (SmMYB21, SmARF6, and SmARF8) related to anther development targeted the SmMYB108 promoter. Overall, our results provide insights into the molecular mechanisms involved in the regulation of anther development by SmMYB108.

Genetic Regulation, Environmental Cues, and Extraction Methods for Higher Yield of Secondary Metabolites in Capsicum

Capsicum (chili pepper) is a widely popular and highly consumed fruit crop with beneficial secondary metabolites such as capsaicinoids, carotenoids, flavonoids, and polyphenols, among others. Interestingly, the secondary metabolite profile is a dynamic function of biosynthetic enzymes, regulatory transcription factors, developmental stage, abiotic and biotic environment, and extraction methods. We propose active manipulable genetic, environmental, and extraction controls for the modulation of quality and quantity of desired secondary metabolites in Capsicum species. Specific biosynthetic genes such as Pun (AT3) and AMT in the capsaicinoids pathway and PSY, LCY, and CCS in the carotenoid pathway can be genetically engineered for enhanced production of capsaicinoids and carotenoids, respectively. Generally, secondary metabolites increase with the ripening of the fruit; however, transcriptional regulators such as MYB, bHLH, and ERF control the extent of accumulation in specific tissues. The precise tuning of biotic and abiotic factors such as light, temperature, and chemical elicitors can maximize the accumulation and retention of secondary metabolites in pre- and postharvest settings. Finally, optimized extraction methods such as ultrasonication and supercritical fluid method can lead to a higher yield of secondary metabolites. Together, the integrated understanding of the genetic regulation of biosynthesis, elicitation treatments, and optimization of extraction methods can maximize the industrial production of secondary metabolites in Capsicum.

Weighted gene co-expression network analysis of nitrogen (N)-responsive genes and the putative role of G-quadruplexes in N use efficiency (NUE) in rice

Rice is an important target to improve crop nitrogen (N) use efficiency (NUE), and the identification and shortlisting of the candidate genes are still in progress. We analyzed data from 16 published N-responsive transcriptomes/microarrays to identify, eight datasets that contained the maximum number of 3020 common genes, referred to as N-responsive genes. These include different classes of transcription factors, transporters, miRNA targets, kinases and events of post-translational modifications. A Weighted gene co-expression network analysis (WGCNA) with all the 3020 N-responsive genes revealed 15 co-expression modules and their annotated biological roles. Protein-protein interaction network analysis of the main module revealed the hub genes and their functional annotation revealed their involvement in the ubiquitin process. Further, the occurrences of G-quadruplex sequences were examined, which are known to play important roles in epigenetic regulation but are hitherto unknown in N-response/NUE. Out of the 3020 N-responsive genes studied, 2298 contained G-quadruplex sequences. We compared these N-responsive genes containing G-quadruplex sequences with the 3601 genes we previously identified as NUE-related (for being both N-responsive and yield-associated). This analysis revealed 389 (17%) NUE-related genes containing G-quadruplex sequences. These genes may be involved in the epigenetic regulation of NUE, while the rest of the 83% (1811) genes may regulate NUE through genetic mechanisms and/or other epigenetic means besides G-quadruplexes. A few potentially important genes/processes identified as associated with NUE were experimentally validated in a pair of rice genotypes contrasting for NUE. The results from the WGCNA and G4 sequence analysis of N-responsive genes helped identify and shortlist six genes as candidates to improve NUE. Further, the hitherto unavailable segregation of genetic and epigenetic gene targets could aid in informed interventions through genetic and epigenetic means of crop improvement.

MYB24 Negatively Regulates the Biosynthesis of Lignin and Capsaicin by Affecting the Expression of Key Genes in the Phenylpropanoid Metabolism Pathway in Capsicum chinense

The wide application of pepper is mostly related to the content of capsaicin, and phenylpropanoid metabolism and its branch pathways may play an important role in the biosynthesis of capsaicin. The expression level of MYB24, a transcription factor screened from the transcriptome data of the pepper fruit development stage, was closely related to the spicy taste. In this experiment, CcMYB24 was cloned from Hainan Huangdenglong pepper, a hot aromatic pepper variety popular in the world for processing, and its function was confirmed by tissue expression characteristics, heterologous transformation in Arabidopsis thaliana, and VIGS technology. The results showed that the relative expression level of CcMYB24 was stable in the early stage of pepper fruit development, and increased significantly from 30 to 50 days after flowering. Heterologous expression led to a significant increase in the expression of CcMYB24 and decrease in lignin content in transgenic Arabidopsis thaliana plants. CcMYB24 silencing led to a significant increase in the expression of phenylpropanoid metabolism pathway genes PAL, 4CL, and pAMT; lignin branch CCR1 and CAD; and capsaicin pathway CS, AT3, and COMT genes in the placenta of pepper, with capsaicin content increased by more than 31.72% and lignin content increased by 20.78%. However, the expression of PAL, pAMT, AT3, COMT, etc., in the corresponding pericarps did not change significantly. Although CS, CCR1, and CAD increased significantly, the relative expression amount was smaller than that in placental tissue, and the lignin content did not change significantly. As indicated above, CcMYB24 may negatively regulate the formation of capsaicin and lignin by regulating the expression of genes from phenylpropanoid metabolism and its branch pathways.

MYB30 and MYB14 form a repressor-activator module with WRKY8 that controls stilbene biosynthesis in grapevine

When exposed to pathogen infection or ultraviolet (UV) radiation, grapevine (Vitis vinifera) plants rapidly accumulate the stilbenoid resveratrol (Res) with concomitant increase of stilbene synthase (STS), the key enzyme in stilbene biosynthesis. Although a few transcription factors have been shown to regulate STSs, the molecular mechanism governing the regulation of STSs is not well elucidated. Our previous work showed that a VvMYB14-VvWRKY8 regulatory loop fine-tunes stilbene biosynthesis in grapevine through protein-protein interaction; overexpression of VvWRKY8 down-regulates VvMYB14 and VvSTS15/21; and application of exogenous Res up-regulates WRKY8 expression. Here, we identified an R2R3-MYB repressor, VvMYB30, which competes with the activator VvMYB14 for binding to the common binding sites in the VvSTS15/21 promoter. Similar to VvMYB14, VvMYB30 physically interacts with VvWRKY8 through their N-termini, forming a complex that does not bind DNA. Exposure to UV-B/C stress induces VvMYB14, VvWRKY8, and VvSTS15/21, but represses VvMYB30 in grapevine leaves. In addition, MYB30 expression is up-regulated by VvWRKY8-overexpression or exogenous Res. These findings suggest that the VvMYB14-VvWRKY8-VvMYB30 regulatory circuit allows grapevine to respond to UV stress by producing Res and prevents over-accumulation of Res to balance metabolic costs. Our work highlights the stress-mediated induction and feedback inhibition of stilbene biosynthesis through a complex regulatory network involving multiple positive and negative transcriptional regulators.

Metabolomic and transcriptomic analyses of the flavonoid biosynthetic pathway in blueberry (Vaccinium spp.)

As a highly economic small fruit crop, blueberry is enjoyed by most people in terms of color, taste, and rich nutrition. To better understand its coloring mechanism on the process of ripening, an integrative analysis of the metabolome and transcriptome profiles was performed in three blueberry varieties at three developmental stages. In this study, 41 flavonoid metabolites closely related to the coloring in blueberry samples were analyzed. It turned out that the most differential metabolites in the ripening processes were delphinidin-3-O-arabinoside (dpara), peonidin-3-O-glucoside (pnglu), and delphinidin-3-O-galactoside (dpgal), while the most differential metabolites among different varieties were flavonols. Furthermore, to obtain more accurate and comprehensive transcripts of blueberry during the developmental stages, PacBio and Illumina sequencing technology were combined to obtain the transcriptome of the blueberry variety Misty, for the very first time. Finally, by applying the gene coexpression network analysis, the darkviolet and bisque4 modules related to flavonoid synthesis were determined, and the key genes related to two flavonoid 3', 5'-hydroxylase (F3'5'H) genes in the darkviolet module and one bHLH transcription factor in the bisque4 module were predicted. It is believed that our findings could provide valuable information for the future study on the molecular mechanism of flavonoid metabolites and flavonoid synthesis pathways in blueberries.

G2-LIKE CAROTENOID REGULATOR (SlGCR) is a positive regulator of lutein biosynthesis in tomato

Lutein is an oxygen-containing carotenoid synthesized in plant chloroplasts and chromoplasts. It plays an indispensable role in promoting plant growth and maintaining eye health in humans. The rate-limiting step of lutein biosynthesis is catalyzed by the lycopene ε-cyclase enzyme (LCYE). Although great progress has been made in the identification of transcription factors involved in the lutein biosynthetic pathway, many systematic molecular mechanisms remain to be elucidated. Here, using co-expression analysis, we identified a gene, G2-LIKE CAROTENOID REGULATOR (SlGCR), encoding a GARP G2-like transcription factor, as the potential regulator of SlLCYE in tomato. Silencing of SlGCR reduced the expression of carotenoid biosynthetic genes and the accumulation of carotenoids in tomato leaves. By contrast, overexpression of SlGCR in tomato fruit significantly increased the expression of relevant genes and enhanced the accumulation of carotenoids. SlGCR can directly bind to the SlLCYE promoter and activate its expression. In addition, we also discovered that expression of SlGCR was negatively regulated by the master regulator SlRIN, thereby inhibiting lutein synthesis during tomato fruit ripening. Taken together, we identified SlGCR as a novel regulator involved in tomato lutein biosynthesis, elucidated the regulatory mechanism, and provided a potential tool for tomato lutein metabolic engineering.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-022-00088-z.

Gas Chromatography-Mass Spectrometry Metabolite Analysis Combined with Transcriptomic and Proteomic Provide New Insights into Revealing Cuticle Formation during Pepper Development

The cuticle plays an important role for the quality of pepper fruit. However, the molecular mechanism of cuticle formation in pepper fruit remains unclear. Our results showed that the wax was continuously accumulated during pepper development, while the cutin monomer first increased and then decreased. Hexadecanoic acid and 10,16-hydroxyhexadecanoic acid were the main components of wax and cutin, respectively. Combined with transcriptome and proteome, the formation patterns of wax and cutin polyester network for pepper cuticle was proposed. The 18 pairs of consistent expression genes and proteins involved in cuticle formation were revealed. Meanwhile, 12 key genes were screened from fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, fatty acid elongation, cutin, suberine, and wax biosynthesis, glycerolipid metabolism, and transport pathway. This study would provide important candidate genes and theoretical basis for the molecular mechanism of cuticle formation, which is essential for the breeding of peppers.

Genome-Wide Analysis of the MYB-Related Transcription Factor Family in Pepper and Functional Studies of CaMYB37 Involvement in Capsaicin Biosynthesis

Chili pepper is an important economic vegetable worldwide. MYB family gene members play an important role in the metabolic processes in plant growth and development. In this study, 103 pepper MYB-related members were identified and grouped into nine subfamilies according to phylogenetic relationships. Additionally, a total of 80, 20, and 37 collinear gene pairs were identified between pepper and tomato, pepper and Arabidopsis, and tomato and Arabidopsis, respectively. We performed promoter cis-element analysis and showed that CaMYB-related members may be involved in multiple biological processes such as growth and development, secondary metabolism, and circadian rhythm regulation. Expression pattern analysis indicated that CaMYB37 is significantly more enriched in fruit placenta, suggesting that this gene may be involved in capsaicin biosynthesis. Through VIGS, we confirmed that CaMYB37 is critical for the biosynthesis of capsaicin in placenta. Our subcellular localization studies revealed that CaMYB37 localized in the nucleus. On the basis of yeast one-hybrid and dual-luciferase reporter assays, we found that CaMYB37 directly binds to the promoter of capsaicin biosynthesis gene AT3 and activates its transcription, thereby regulating capsaicin biosynthesis. In summary, we systematically identified members of the CaMYB-related family, predicted their possible biological functions, and revealed that CaMYB37 is critical for the transcriptional regulation of capsaicin biosynthesis. This work provides a foundation for further studies of the CaMYB-related family in pepper growth and development.

Integrated omics analysis identified genes and their splice variants involved in fruit development and metabolites production in Capsicum species

To date, several transcriptomic studies during fruit development have been reported; however, no comprehensive integrated study on expression diversity, alternative splicing, and metabolomic profiling was reported in Capsicum. This study analyzed RNA-seq data and untargeted metabolomic profiling from early green (EG), mature green (MG), and breaker (Br) fruit stages from two Capsicum species, i.e., C. annuum (Cann) and C. frutescens (Cfrut) from Northeast India. A total of 117,416 and 96,802 alternatively spliced events (AltSpli-events) were identified from Cann and Cfrut, respectively. Among AltSpli-events, intron retention (IR; 32.2% Cann and 25.75% Cfrut) followed by alternative acceptor (AA; 15.4% Cann and 18.9% Cfrut) were the most abundant in Capsicum. Around 7600 genes expressed in at least one fruit stage of Cann and Cfrut were AltSpli. The study identified spliced variants of genes including transcription factors (TFs) potentially involved in fruit development/ripening (Aux/IAA 16-like, ETR, SGR1, ARF, CaGLK2, ETR, CaAGL1, MADS-RIN, FUL1, SEPALLATA1), carotenoid (PDS, CA1, CCD4, NCED3, xanthoxin dehydrogenase, CaERF82, CabHLH100, CaMYB3R-1, SGR1, CaWRKY28, CaWRKY48, CaWRKY54), and capsaicinoids or flavonoid biosynthesis (CaMYB48, CaWRKY51), which were significantly differentially spliced (DS) between consecutive Capsicum fruit stages. Also, this study observed that differentially expressed isoforms (DEiso) from 38 genes with differentially spliced events (DSE) were significantly enriched in various metabolic pathways such as starch and sucrose metabolism, amino acid metabolism, cysteine cutin suberin and wax biosynthesis, and carotenoid biosynthesis. Furthermore, the metabolomic profiling revealed that metabolites from aforementioned pathways such as carbohydrates (mainly sugars such as D-fructose, D-galactose, maltose, and sucrose), organic acids (carboxylic acids), and peptide groups significantly altered during fruit development. Taken together, our findings could help in alternative splicing-based targeted studies of candidate genes involved in fruit development and ripening in Capsicum crop.

Genomic selection with fixed-effect markers improves the prediction accuracy for Capsaicinoid contents in Capsicum annuum

Capsaicinoids provide chili peppers (Capsicum spp.) with their characteristic pungency. Several structural and transcription factor genes are known to control capsaicinoid contents in pepper. However, many other genes also regulating capsaicinoid contents remain unknown, making it difficult to develop pepper cultivars with different levels of capsaicinoids. Genomic selection (GS) uses genome-wide random markers (including many in undiscovered genes) for a trait to improve selection efficiency. In this study, we predicted the capsaicinoid contents of pepper breeding lines using several GS models trained with genotypic and phenotypic data from a training population. We used a core collection of 351 Capsicum accessions and 96 breeding lines as training and testing populations, respectively. To obtain the optimal number of single nucleotide polymorphism (SNP) markers for GS, we tested various numbers of genome-wide SNP markers based on linkage disequilibrium. We obtained the highest mean prediction accuracy (0.550) for different models using 3294 SNP markers. Using this marker set, we conducted GWAS and selected 25 markers that were associated with capsaicinoid biosynthesis genes and quantitative trait loci for capsaicinoid contents. Finally, to develop more accurate prediction models, we obtained SNP markers from GWAS as fixed-effect markers for GS, where 3294 genome-wide SNPs were employed. When four to five fixed-effect markers from GWAS were used as fixed effects, the RKHS and RR-BLUP models showed accuracies of 0.696 and 0.689, respectively. Our results lay the foundation for developing pepper cultivars with various capsaicinoid levels using GS for capsaicinoid contents.

Whole-Genome Identification and Analysis of Multiple Gene Families Reveal Candidate Genes for Theasaponin Biosynthesis in Camellia oleifera

Camellia oleifera is an economically important oilseed tree. Seed meals of C. oleifera have a long history of use as biocontrol agents in shrimp farming and as cleaning agents in peoples’ daily lives due to the presence of theasaponins, the triterpene saponins from the genus Camellia. To characterize the biosynthetic pathway of theasaponins in C. oleifera, members of gene families involved in triterpenoid biosynthetic pathways were identified and subjected to phylogenetic analysis with corresponding members in Arabidopsis thaliana, Camellia sinensis, Actinidia chinensis, Panax ginseng, and Medicago truncatula. In total, 143 triterpenoid backbone biosynthetic genes, 1169 CYP450s, and 1019 UGTs were identified in C. oleifera. The expression profiles of triterpenoid backbone biosynthetic genes were analyzed in different tissue and seed developmental stages of C. oleifera. The results suggested that MVA is the main pathway for triterpenoid backbone biosynthesis. Moreover, the candidate genes for theasaponin biosynthesis were identified by WGCNA and qRT-PCR analysis; these included 11 CYP450s, 14 UGTs, and eight transcription factors. Our results provide valuable information for further research investigating the biosynthetic and regulatory network of theasaponins.

Universal gene co-expression network reveals receptor-like protein genes involved in broad-spectrum resistance in pepper (Capsicum annuum L.)

Receptor-like proteins (RLPs) on plant cells have been implicated in immune responses and developmental processes. Although hundreds of RLP genes have been identified in plants, only a few RLPs have been functionally characterized in a limited number of plant species. Here, we identified RLPs in the pepper (Capsicum annuum) genome and performed comparative transcriptomics coupled with the analysis of conserved gene co-expression networks (GCNs) to reveal the role of core RLP regulators in pepper-pathogen interactions. A total of 102 RNA-seq datasets of pepper plants infected with four pathogens were used to construct CaRLP-targeted GCNs (CaRLP-GCNs). Resistance-responsive CaRLP-GCNs were merged to construct a universal GCN. Fourteen hub CaRLPs, tightly connected with defense-related gene clusters, were identified in eight modules. Based on the CaRLP-GCNs, we evaluated whether hub CaRLPs in the universal GCN are involved in the biotic stress response. Of the nine hub CaRLPs tested by virus-induced gene silencing, three genes (CaRLP264, CaRLP277, and CaRLP351) showed defense suppression with less hypersensitive response-like cell death in race-specific and non-host resistance response to viruses and bacteria, respectively, and consistently enhanced susceptibility to Ralstonia solanacearum and/or Phytophthora capsici. These data suggest that key CaRLPs are involved in the defense response to multiple biotic stresses and can be used to engineer a plant with broad-spectrum resistance. Together, our data show that generating a universal GCN using comprehensive transcriptome datasets can provide important clues to uncover genes involved in various biological processes.

Ethylene-Inducible AP2/ERF Transcription Factor Involved in the Capsaicinoid Biosynthesis in Capsicum

Ethylene is very important in the process of plant development and regulates the biosynthesis of many secondary metabolites. In these regulatory mechanisms, transcription factors (TFs) that mediate ethylene signals play a very important role. Capsaicinoids (CAPs) are only synthesized and accumulated in Capsicum species, causing their fruit to have a special pungent taste, which can protect against attack from herbivores and pathogens. In this study, we identified the TF CcERF2, which is induced by ethylene, and demonstrated its regulatory effect on CAPs biosynthesis. Transcriptome sequencing analysis revealed that the expression patterns of CcERF2 and multiple genes associated with CAPs biosynthesis were basically the same. The spatiotemporal expression results showed CcERF2 was preferentially expressed in the placenta of the spicy fruit. Ethylene can induce the expression of CcERF2 and CAPs biosynthesis genes (CBGs). CcERF2 gene silencing and 1-methylcyclopropene (1-MCP) and pyrazinamide (PZA) treatments caused a decrease in expression of CBGs and a sharp decrease in content of CAPs. The results indicated that CcERF2 was indeed involved in the regulation of structural genes of the CAPs biosynthetic pathway.

Transcriptome Profiles of Leaves and Roots of Goldenrain Tree (Koelreuteria paniculata Laxm.) in Response to Cadmium Stress

Cadmium (Cd) pollution is a widespread environmental problem. In this study, we explored the transcriptome and biochemical responses of goldenrain tree (Koelreuteria paniculata Laxm.) leaves and roots to Cd stress. Leaf and root growth decreased substantially under Cd stress (50 mg/L CdCl₂), but leaf and root antioxidant mechanisms were significantly activated. In RNA-seq analysis, roots treated with 25 mg/L CdCl₂ featured enriched GO terms in cellular components related to intracellular ribonucleoprotein complex, ribonucleoprotein complex, and macromolecular complex. In leaves under Cd stress, most differentially expressed genes were enriched in the cellular component terms intrinsic component of membrane and membrane part. Weighted gene co-expression network analysis and analysis of module–trait relations revealed candidate genes associated with superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities and malondialdehyde (MDA). Ten transcription factors responded to Cd stress expression, including those in C2H2, MYB, WRKY, and bZIP families. Transcriptomic analysis of goldenrain tree revealed that Cd stress rapidly induced the intracellular ribonucleoprotein complex in the roots and the intrinsic component of membrane in the leaves. The results also indicate directions for further analyses of molecular mechanisms of Cd tolerance and accumulation in goldenrain tree.

Transcriptional Regulation of Ripening in Chili Pepper Fruits (Capsicum spp.)

Chili peppers represent a very important horticultural crop that is cultivated and commercialized worldwide. The ripening process makes the fruit palatable, desirable, and attractive, thus increasing its quality and nutritional value. This process includes visual changes, such as fruit coloration, flavor, aroma, and texture. Fruit ripening involves a sequence of physiological, biochemical, and molecular changes that must be finely regulated at the transcriptional level. In this review, we integrate current knowledge about the transcription factors involved in the regulation of different stages of the chili pepper ripening process.

Capsicum chinense MYB Transcription Factor Genes: Identification, Expression Analysis, and Their Conservation and Diversification With Other Solanaceae Genomes

Myeloblastosis (MYB) genes are important transcriptional regulators of plant growth, development, and secondary metabolic biosynthesis pathways, such as capsaicinoid biosynthesis in Capsicum. Although MYB genes have been identified in Capsicum annuum, no comprehensive study has been conducted on other Capsicum species. We identified a total of 251 and 240 MYB encoding genes in Capsicum chinense MYBs (CcMYBs) and Capsicum baccatum MYBs (CbMYBs). The observation of twenty tandem and 41 segmental duplication events indicated expansion of the MYB gene family in the C. chinense genome. Five CcMYB genes, i.e., CcMYB101, CcMYB46, CcMYB6, CcPHR8, and CcRVE5, and two CaMYBs, i.e., CaMYB3 and CaHHO1, were found within the previously reported capsaicinoid biosynthesis quantitative trait loci. Based on phylogenetic analysis with tomato MYB proteins, the Capsicum MYBs were classified into 24 subgroups supported by conserved amino acid motifs and gene structures. Also, a total of 241 CcMYBs were homologous with 225 C. annuum, 213 C. baccatum, 125 potato, 79 tomato, and 23 Arabidopsis MYBs. Synteny analysis showed that all 251 CcMYBs were collinear with C. annuum, C. baccatum, tomato, potato, and Arabidopsis MYBs spanning over 717 conserved syntenic segments. Using transcriptome data from three fruit developmental stages, a total of 54 CcMYBs and 81 CaMYBs showed significant differential expression patterns. Furthermore, the expression of 24 CcMYBs from the transcriptome data was validated by quantitative real-time (qRT) PCR analysis. Eight out of the 24 CcMYBs validated by the qRT-PCR were highly expressed in fiery hot C. chinense than in the lowly pungent C. annuum. Furthermore, the co-expression analysis revealed several MYB genes clustered with genes from the capsaicinoid, anthocyanin, phenylpropanoid, carotenoid, and flavonoids biosynthesis pathways, and related to determining fruit shape and size. The homology modeling of 126 R2R3 CcMYBs showed high similarity with that of the Arabidopsis R2R3 MYB domain template, suggesting their potential functional similarity at the proteome level. Furthermore, we have identified simple sequence repeat (SSR) motifs in the CcMYB genes, which could be used in Capsicum breeding programs. The functional roles of the identified CcMYBs could be studied further so that they can be manipulated for Capsicum trait improvement.

Genome-wide identification of the Capsicum bHLH transcription factor family: discovery of a candidate regulator involved in the regulation of species-specific bioactive metabolites

BACKGROUND

The basic helix-loop-helix (bHLH) transcription factors (TFs) serve crucial roles in regulating plant growth and development and typically participate in biological processes by interacting with other TFs. Capsorubin and capsaicinoids are found only in Capsicum, which has high nutritional and economic value. However, whether bHLH family genes regulate capsorubin and capsaicinoid biosynthesis and participate in these processes by interacting with other TFs remains unknown.

RESULTS

In this study, a total of 107 CabHLHs were identified from the Capsicum annuum genome. Phylogenetic tree analysis revealed that these CabHLH proteins were classified into 15 groups by comparing the CabHLH proteins with Arabidopsis thaliana bHLH proteins. The analysis showed that the expression profiles of CabHLH009, CabHLH032, CabHLH048, CabHLH095 and CabHLH100 found in clusters C1, C2, and C3 were similar to the profile of carotenoid biosynthesis in pericarp, including zeaxanthin, lutein and capsorubin, whereas the expression profiles of CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 found in clusters L5, L6 and L9 were consistent with the profile of capsaicinoid accumulation in the placenta. Moreover, CabHLH007, CabHLH009, CabHLH026 and CabHLH086 also might be involved in temperature-mediated capsaicinoid biosynthesis. Yeast two-hybrid (Y2H) assays demonstrated that CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 could interact with MYB31, a master regulator of capsaicinoid biosynthesis.

CONCLUSIONS

The comprehensive and systematic analysis of CabHLH TFs provides useful information that contributes to further investigation of CabHLHs in carotenoid and capsaicinoid biosynthesis.

Identification of Chlorophyll Metabolism- and Photosynthesis-Related Genes Regulating Green Flower Color in Chrysanthemum by Integrative Transcriptome and Weighted Correlation Network Analyses

Green chrysanthemums are difficult to breed but have high commercial value. The molecular basis for the green petal color in chrysanthemum is not fully understood. This was investigated in the present study by RNA sequencing analysis of white and green ray florets collected at three stages of flower development from the F₁ progeny of the cross between Chrysanthemum × morifolium “Lüdingdang” with green-petaled flowers and Chrysanthemum vistitum with white-petaled flowers. The chlorophyll content was higher and chloroplast degradation was slower in green pools than in white pools at each developmental stage. Transcriptome analysis revealed that genes that were differentially expressed between the two pools were enriched in pathways related to chlorophyll metabolism and photosynthesis. We identified the transcription factor genes CmCOLa, CmCOLb, CmERF, and CmbHLH as regulators of the green flower color in chrysanthemum by differential expression analysis and weighted gene co-expression network analysis. These findings can guide future efforts to improve the color palette of chrysanthemum flowers through genetic engineering.

Genome-Wide Identification and Analysis of the MYB Transcription Factor Gene Family in Chili Pepper (Capsicum spp.)

The MYB transcription factor family is very large and functionally diverse in plants, however, only a few members of this family have been reported and characterized in chili pepper (Capsicum spp.). In the present study, we performed genome-wide analyses of the MYB family in Capsicum annuum, including phylogenetic relationships, conserved domain, gene structure organization, motif protein arrangement, chromosome distribution, chemical properties predictions, RNA-seq expression, and RT-qPCR expression assays. A total of 235 non-redundant MYB proteins were identified from C. annuum, including R2R3-MYB, 3R-MYB, atypical MYB, and MYB-related subclasses. The sequence analysis of CaMYBs compared with other plant MYB proteins revealed gene conservation, but also potential specialized genes. Tissue-specific expression profiles showed that CaMYB genes were differentially expressed, suggesting that they are functionally divergent. Furthermore, the integration of our data allowed us to propose strong CaMYBs candidates to be regulating phenylpropanoid, lignin, capsaicinoid, carotenoid, and vitamin C biosynthesis, providing new insights into the role of MYB transcription factors in secondary metabolism. This study adds valuable knowledge about the functions of CaMYB genes in various processes in the Capsicum genus.