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.

A putative E3 ubiquitin ligase substrate receptor degrades transcription factor SmNAC to enhance bacterial wilt resistance in eggplant.

Bacterial wilt caused by Ralstonia solanacearum is a severe soil-borne disease globally, limiting the production in Solanaceae plants. SmNAC negatively regulated eggplant resistance to Bacterial wilt (BW) though restraining salicylic acid (SA) biosynthesis. However, other mechanisms through which SmNAC regulates BW resistance remain unknown. Here, we identified an interaction factor, SmDDA1b, encoding a substrate receptor for E3 ubiquitin ligase, from the eggplant cDNA library using SmNAC as bait. SmDDA1b expression was promoted by R. solanacearum inoculation and exogenous SA treatment. The virus-induced gene silencing of the SmDDA1b suppressed the BW resistance of eggplants; SmDDA1b overexpression enhanced the BW resistance of tomato plants. SmDDA1b positively regulates BW resistance by inhibiting the spread of R. solanacearum within plants. The SA content and the SA biosynthesis gene ICS1 and signaling pathway genes decreased in the SmDDA1b-silenced plants but increased in SmDDA1b-overexpression plants. Moreover, SmDDB1 protein showed interaction with SmCUL4 and SmDDA1b and protein degradation experiments indicated that SmDDA1b reduced SmNAC protein levels through proteasome degradation. Furthermore, SmNAC could directly bind the SmDDA1b promoter and repress its transcription. Thus, SmDDA1b is a novel regulator functioning in BW resistance of solanaceous crops via the SmNAC-mediated SA pathway. Those results also revealed a negative feedback loop between SmDDA1b and SmNAC controlling BW resistance.

Identification and characterization analysis of candidate genes controlling mushroom leaf development in Chinese kale by BSA-seq.

Mushroom leaves (MLs) are malformed leaves that develop from the leaf veins in some of Chinese kale genotypes. To study the genetic model and molecular mechanism of ML development in Chinese kale, the F2 segregation population was constructed by two inbred lines, genotype Boc52 with ML and genotype Boc55 with normal leaves (NL). In the present study, we have identified for the first time that the development of mushroom leaves may be affected by the change of adaxial-abaxial polarity of leaves. Examination of the phenotypes of F1 and F2 segregation populations suggested that ML development is controlled by two dominant major genes inherited independently. BSA-seq analysis showed that a major quantitative trait locus (QTL) qML4.1 that controls ML development is located within 7.4 Mb on chromosome kC4. The candidate region was further narrowed to 255 kb by linkage analysis combined with insertion/deletion (InDel) markers, and 37 genes were predicted in this region. According to the expression and annotation analysis, a B3 domain-containing transcription factor NGA1-like gene, BocNGA1, was identified as a key candidate gene for controlling ML development in Chinese kale. Fifteen single nucleotide polymorphisms (SNPs) were found in coding sequences and 21 SNPs and 3 InDels found in the promoter sequences of BocNGA1 from the genotype Boc52 with ML. The expression levels of BocNGA1 in ML genotypes are significantly lower than in the NL genotypes, which suggests that BocNGA1 may act as a negative regulator for ML genesis in Chinese kale. This study provides a new foundation for Chinese kale breeding and for the study of the molecular mechanism of plant leaf differentiation. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01364-6.

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.

BocODD1 and BocODD2 Regulate the Biosynthesis of Progoitrin Glucosinolate in Chinese Kale.

Progoitrin (2-hydroxy-3-butenyl glucosinolate, PRO) is the main source of bitterness of Brassica plants. Research on the biosynthesis of PRO glucosinolate can aid the understanding of the nutritional value in Brassica plants. In this study, four ODD genes likely involved in PRO biosynthesis were cloned from Chinese kale. These four genes, designated as BocODD1-4, shared 75-82% similarities with the ODD sequence of Arabidopsis. The sequences of these four BocODDs were analyzed, and BocODD1 and BocODD2 were chosen for further study. The gene BocODD1,2 showed the highest expression levels in the roots, followed by the leaves, flowers, and stems, which is in accordance with the trend of the PRO content in the same tissues. Both the expression levels of BocODD1,2 and the content of PRO were significantly induced by high- and low-temperature treatments. The function of BocODDs involved in PRO biosynthesis was identified. Compared with the wild type, the content of PRO was increased twofold in the over-expressing BocODD1 or BocODD2 plants. Meanwhile, the content of PRO was decreased in the BocODD1 or BocODD2 RNAi lines more than twofold compared to the wildtype plants. These results suggested that BocODD1 and BocODD2 may play important roles in the biosynthesis of PRO glucosinolate in Chinese kale.

Genome-wide analysis of histone acetyltransferase and histone deacetylase families and their expression in fruit development and ripening stage of pepper (Capsicum annuum).

The fruit development and ripening process involve a series of changes regulated by fine-tune gene expression at the transcriptional level. Acetylation levels of histones on lysine residues are dynamically regulated by histone acetyltransferases (HATs) and histone deacetylases (HDACs), which play an essential role in the control of gene expression. However, their role in regulating fruit development and ripening process, especially in pepper (Capsicum annuum), a typical non-climacteric fruit, remains to understand. Herein, we performed genome-wide analyses of the HDAC and HAT family in the pepper, including phylogenetic analysis, gene structure, encoding protein conserved domain, and expression assays. A total of 30 HAT and 15 HDAC were identified from the pepper genome and the number of gene differentiation among species. The sequence and phylogenetic analysis of CaHDACs and CaHATs compared with other plant HDAC and HAT proteins revealed gene conserved and potential genus-specialized genes. Furthermore, fruit developmental trajectory expression profiles showed that CaHDAC and CaHAT genes were differentially expressed, suggesting that some are functionally divergent. The integrative analysis allowed us to propose CaHDAC and CaHAT candidates to be regulating fruit development and ripening-related phytohormone metabolism and signaling, which also accompanied capsaicinoid and carotenoid biosynthesis. This study provides new insights into the role of histone modification mediate development and ripening in non-climacteric fruits.

The 3D architecture of the pepper genome and its relationship to function and evolution.

The organization of chromatin into self-interacting domains is universal among eukaryotic genomes, though how and why they form varies considerably. Here we report a chromosome-scale reference genome assembly of pepper (Capsicum annuum) and explore its 3D organization through integrating high-resolution Hi-C maps with epigenomic, transcriptomic, and genetic variation data. Chromatin folding domains in pepper are as prominent as TADs in mammals but exhibit unique characteristics. They tend to coincide with heterochromatic regions enriched with retrotransposons and are frequently embedded in loops, which may correlate with transcription factories. Their boundaries are hotspots for chromosome rearrangements but are otherwise depleted for genetic variation. While chromatin conformation broadly affects transcription variance, it does not predict differential gene expression between tissues. Our results suggest that pepper genome organization is explained by a model of heterochromatin-driven folding promoted by transcription factories and that such spatial architecture is under structural and functional constraints.

Overexpression of EgrIAA20 from Eucalyptus grandis, a Non-Canonical Aux/IAA Gene, Specifically Decouples Lignification of the Different Cell-Types in Arabidopsis Secondary Xylem.

Wood (secondary xylem) formation is regulated by auxin, which plays a pivotal role as an integrator of developmental and environmental cues. However, our current knowledge of auxin-signaling during wood formation is incomplete. Our previous genome-wide analysis of Aux/IAAs in Eucalyptus grandis showed the presence of the non-canonical paralog member EgrIAA20 that is preferentially expressed in cambium. We analyzed its cellular localization using a GFP fusion protein and its transcriptional activity using transactivation assays, and demonstrated its nuclear localization and strong auxin response repressor activity. In addition, we functionally tested the role of EgrIAA20 by constitutive overexpression in Arabidopsis to investigate for phenotypic changes in secondary xylem formation. Transgenic Arabidopsis plants overexpressing EgrIAA20 were smaller and displayed impaired development of secondary fibers, but not of other wood cell types. The inhibition in fiber development specifically affected their cell wall lignification. We performed yeast-two-hybrid assays to identify EgrIAA20 protein partners during wood formation in Eucalyptus, and identified EgrIAA9A, whose ortholog PtoIAA9 in poplar is also known to be involved in wood formation. Altogether, we showed that EgrIAA20 is an important auxin signaling component specifically involved in controlling the lignification of wood fibers.

The Capsicum MYB31 regulates capsaicinoid biosynthesis in the pepper pericarp.

Pepper (Capsicum) are consumed worldwide as vegetables and food additives due to their pungent taste. Capsaicinoids are the bioactive compounds that confer the desired pungency to pepper fruits. Capsaicinoid biosynthesis was thought to occur exclusively in fruit placenta. Recently, biosynthesis in the pericarp of extremely pungent varieties was discovered, however, the mechanism of capsaicinoid biosynthesis regulation in the pericarp remains largely unknown. Here, the capsaicinoid contents of placenta and pericarp were analyzed. The results indicated that the Capsicum chinense pericarp accumulated a vast amount of capsaicinoids. Expression of the master regulator MYB31 and capsaicinoid biosynthesis genes (CBGs) were significantly upregulated in the pericarp in C. chinense accessions compared to accessions in other tested species. Moreover, in fruit of extremely-pungent 'Trinidad Moruga Scorpion' (C. chinense) and low-pungent '59' inbred line (C. annuum), the capsaicinoid accumulation patterns in the pericarp were consistent with expression levels of CBGs and MYB31. Silencing MYB31 in 'Trinidad Moruga Scorpion' pericarp leads to a significantly decreased CBGs transcription level and capsaicinoids content. Taken together, our results provide insights into the molecular mechanism arising from the expression of MYB31 in the pericarp that results in exceedingly hot peppers.

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.

Specific-Locus Amplified Fragment Sequencing (SLAF-Seq) as High-Throughput SNP Genotyping Methods.

Most plant agronomic traits are quantitatively inherited. Identification of quantitative trait loci (QTL) is a challenging target for most scientists and crop breeders as large-scale genotyping is difficult. Molecular marker technology has continuously evolved from hybridization-based technology to PCR-based technology, and finally, sequencing-based high-throughput single-nucleotide polymorphisms (SNPs). High-throughput sequencing technologies can provide strategies for sequence-based SNP genotyping. Here we describe the SLAF-seq that can be applied as the SNP genotyping approach. The high-throughput SNP genotyping methods will prove useful for the construction of high-density genetic maps and identification of QTLs for their deployment in plant breeding and facilitate genome-wide selection (GWS) and genome-wide association studies (GWAS).

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

Plant biosynthesis involves numerous specialized metabolites with diverse chemical natures and biological activities. The biosynthesis of metabolites often exclusively occurs in response to tissue-specific combinatorial developmental cues that are controlled at the transcriptional level. Capsaicinoids are a group of specialized metabolites that confer a pungent flavor to pepper fruits. Capsaicinoid biosynthesis occurs in the fruit placenta and combines its developmental cues. Although the capsaicinoid biosynthetic pathway has been largely characterized, the regulatory mechanisms that control capsaicinoid metabolism have not been fully elucidated. In this study, we combined fruit placenta transcriptome data with weighted gene coexpression network analysis (WGCNA) to generate coexpression networks. A capsaicinoid-related gene module was identified in which the MYB transcription factor CaMYB48 plays a critical role in regulating capsaicinoid in pepper. Capsaicinoid biosynthetic gene (CBG) and CaMYB48 expression primarily occurs in the placenta and is consistent with capsaicinoid biosynthesis. CaMYB48 encodes a nucleus-localized protein that primarily functions as a transcriptional activator through its C-terminal activation motif. CaMYB48 regulates capsaicinoid biosynthesis by directly regulating the expression of CBGs, including AT3a and KasIa. Taken together, the results of this study indicate ways to generate robust networks optimized for the mining of CBG-related regulators, establishing a foundation for future research elucidating capsaicinoid regulation.

Systematic analysis of the Capsicum ERF transcription factor family: identification of regulatory factors involved in the regulation of species-specific metabolites.

BACKGROUND: ERF transcription factors (TFs) belong to the Apetala2/Ethylene responsive Factor (AP2/ERF) TF family and play a vital role in plant growth and development processes. Capsorubin and capsaicinoids have relatively high economic and nutritional value, and they are specifically found in Capsicum. However, there is little understanding of how ERFs participate in the regulatory networks of capsorubin and capsaicinoids biosynthesis. RESULTS: In this study, a total of 142 ERFs were identified in the Capsicum annuum genome. Subsequent phylogenetic analysis allowed us to divide ERFs into DREB (dehydration responsive element binding proteins) and ERF subfamilies, and further classify them into 11 groups with several subgroups. Expression analysis of biosynthetic pathway genes and CaERFs facilitated the identification of candidate genes related to the regulation of capsorubin and capsaicinoids biosynthesis; the candidates were focused in cluster C9 and cluster C10, as well as cluster L3 and cluster L4, respectively. The expression patterns of CaERF82, CaERF97, CaERF66, CaERF107 and CaERF101, which were found in cluster C9 and cluster C10, were consistent with those of accumulating of carotenoids (ß-carotene, zeaxanthin and capsorubin) in the pericarp. In cluster L3 and cluster L4, the expression patterns of CaERF102, CaERF53, CaERF111 and CaERF92 were similar to those of the accumulating capsaicinoids. Furthermore, CaERF92, CaERF102 and CaERF111 were found to be potentially involved in temperature-mediated capsaicinoids biosynthesis. CONCLUSION: This study will provide an extremely useful foundation for the study of candidate ERFs in the regulation of carotenoids and capsaicinoids biosynthesis in peppers.

Jasmonate-Inducible R2R3-MYB Transcription Factor Regulates Capsaicinoid Biosynthesis and Stamen Development in Capsicum.

Jasmonates (JAs) play an important role in plant developmental processes and regulate the biosynthesis of various specialized metabolites, and transcription factors are crucial in mediating JA signaling to regulate these processes. Capsaicinoids (Caps) are intriguing specialized metabolites produced uniquely by Capsicum species that give their fruits a pungent flavor to defend against herbivory and pathogens. In this study, we identify a R2R3-MYB transcription factor CaMYB108 and demonstrate its roles in regulating the biosynthesis of Caps and stamen development. Transcriptional analysis indicated that CaMYB108 was preferentially expressed in the flower and fruit, while the subcellular localization of CaMYB108 was shown to be the nucleus. Virus-induced gene silencing of CaMYB108 led to the expression of capsaicinoid biosynthetic genes (CBGs), and the contents of Caps dramatically reduce. Moreover, the CaMYB108-silenced plants showed delayed anther dehiscence and reduced pollen viability. Transient overexpression of CaMYB108 caused the expression of CBGs to be upregulated, and the Caps content significantly increased. The results of dual-luciferase reporter assays showed that CaMYB108 targeted CBG promoters. In addition, the expression of CaMYB108 and CBGs was inducible by methyl jasmonate and was consistent with the increased content of Caps. Overall, our results indicate that CaMYB108 is involved in the regulation of Caps biosynthesis and stamen development.

Natural variations in the MYB transcription factor MYB31 determine the evolution of extremely pungent peppers.

Plants produce countless specialized metabolites crucial for their development and fitness, and many are useful bioactive compounds. Capsaicinoids are intriguing genus-specialized metabolites that confer a pungent flavor to Capsicum fruits, and they are widely applied in different areas. Among the five domesticated Capsicum species, Capsicum chinense has a high content of capsaicinoids, which results in an extremely hot flavor. However, the species-specific upregulation of capsaicinoid-biosynthetic genes (CBGs) and the evolution of extremely pungent peppers are not well understood. We conducted genetic and functional analyses demonstrating that the quantitative trait locus Capsaicinoid1 (Cap1), which is identical to Pun3 contributes to the level of pungency. The Cap1/Pun3 locus encodes the Solanaceae-specific MYB transcription factor MYB31. Capsicum species have evolved placenta-specific expression of MYB31, which directly activates expression of CBGs and results in genus-specialized metabolite production. The capsaicinoid content depends on MYB31 expression. Natural variations in the MYB31 promoter increase MYB31 expression in C. chinense via the binding of the placenta-specific expression of transcriptional activator WRKY9 and augmentation of CBG expression, which promotes capsaicinoid biosynthesis. Our findings provide insights into the evolution of extremely pungent C. chinense, which is due to natural variations in the master regulator, and offers targets for engineering or selecting flavor in Capsicum.

Construction of a high density genetic map of an interspecific cross of Capsicum chinense and Capsicum annuum and QTL analysis of floral traits.

The yield of pepper plants (Capsicum spp.) is their most important trait and is affected by the flower number and flowering time. Capsicum annuum produces a single flower per node and has an early flowering habit. By contrast, Capsicum chinense yields multiple flowers per node and has a late flowering character. However, the genetic mechanism underlying the control of these floral traits remains largely unknown. In this study, 150 F2 populations from an interspecific cross between the inbred lines 740 (C. chinense) and CA1 (C. annuum) and their parents were used to construct a molecular genetic linkage map using the specific length amplified fragment sequencing (SLAF-seq) technique. This linkage map, spanning 1,586.78 cM in length, contained 9,038 markers on 12 chromosomes, with a mean marker distance of 0.18 cM. Phenotypic data on the flowering time and flower number per node were collected over multiple years, and QTL analysis identified 6 QTLs for the flowering time and flower number per node by composite interval mapping (CIM) and genome-wide composite interval mapping (GCIM) methods at least in two environments. The candidate genes within the major QTL were predicted. In the major flowering time QTL, the candidate gene Capana02g000700, which encodes the homeotic protein APETALA2, was identified. Quantitative reverse-transcription PCR (qRT-PCR) analysis indicated that its expression level in 740 was higher than that in CA1. Gene expression analysis indicated that the expression of Capana02g000700 was significantly upregulated in flowers, and many floral development-related genes were found to be coexpressed with Capana02g000700, supporting the function of this gene in association with flowering time in C. chinense and C. annuum species.

Overexpression of AtEDT1/HDG11 in Chinese Kale (Brassica oleracea var. alboglabra) Enhances Drought and Osmotic Stress Tolerance.

Plants are constantly challenged by environmental stresses, including drought and high salinity. Improvement of drought and osmotic stress tolerance without yield decrease has been a great challenge in crop improvement. The Arabidopsis ENHANCED DROUGHT TOLERANCE1/HOMEODOMAIN GLABROUS11 (AtEDT1/HDG11), a protein of the class IV HD-Zip family, has been demonstrated to significantly improve drought tolerance in Arabidopsis, rice, and pepper. Here, we report that AtEDT1/HDG11 confers drought and osmotic stress tolerance in the Chinese kale. AtEDT1/HDG11-overexpression lines exhibit auxin-overproduction phenotypes, such as long hypocotyls, tall stems, more root hairs, and a larger root system architecture. Compared with the untransformed control, transgenic lines have significantly reduced stomatal density. In the leaves of transgenic Chinese kale plants, proline (Pro) content and reactive oxygen species-scavenging enzyme activity was significantly increased after drought and osmotic stress, particularly compared to wild kale. More importantly, AtEDT1/HDG11-overexpression leads to abscisic acid (ABA) hypersensitivity, resulting in ABA inhibitor germination and induced stomatal closure. Consistent with observed phenotypes, the expression levels of auxin, ABA, and stress-related genes were also altered under both normal and/or stress conditions. Further analysis showed that AtEDT1/HDG11, as a transcription factor, can target the auxin biosynthesis gene YUCC6 and ABA response genes ABI3 and ABI5. Collectively, our results provide a new insight into the role of AtEDT1/HDG11 in enhancing abiotic stress resistance through auxin- and ABA-mediated signaling response in Chinese kale.

Purple foliage coloration in tea (Camellia sinensis L.) arises from activation of the R2R3-MYB transcription factor CsAN1.

Purple foliage always appears in Camellia sinensis families; however, the transcriptional regulation of anthocyanin biosynthesis is unknown. The tea bud sport cultivar 'Zijuan' confers an abnormal pattern of anthocyanin accumulation, resulting in a mutant phenotype that has a striking purple color in young foliage and in the stem. In this study, we aimed to unravel the underlying molecular mechanism of anthocyanin biosynthetic regulation in C. sinensis. Our results revealed that activation of the R2R3-MYB transcription factor (TF) anthocyanin1 (CsAN1) specifically upregulated the bHLH TF CsGL3 and anthocyanin late biosynthetic genes (LBGs) to confer ectopic accumulation of pigment in purple tea. We found CsAN1 interacts with bHLH TFs (CsGL3 and CsEGL3) and recruits a WD-repeat protein CsTTG1 to form the MYB-bHLH-WDR (MBW) complex that regulates anthocyanin accumulation. We determined that the hypomethylation of a CpG island in the CsAN1 promoter is associated with the purple phenotype. Furthermore, we demonstrated that low temperature and long illumination induced CsAN1 promoter demethylation, resulting in upregulated expression to promote anthocyanin accumulation in the foliage. The successful isolation of CsAN1 provides important information on the regulatory control of anthocyanin biosynthesis in C. sinensis and offers a genetic resource for the development of new varieties with enhanced anthocyanin content.