- Research article
- Open Access
Genome-wide identification of the Capsicum bHLH transcription factor family: discovery of a candidate regulator involved in the regulation of species-specific bioactive metabolites
BMC Plant Biology volume 21, Article number: 262 (2021)
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.
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.
The comprehensive and systematic analysis of CabHLH TFs provides useful information that contributes to further investigation of CabHLHs in carotenoid and capsaicinoid biosynthesis.
Peppers (Capsicum spp.), which include sweet and hot chili varieties, are important economic crops in the world . Five domesticated Capsicum species, namely, Capsicum annuum, Capsicum baccatum, Capsicum chinense, Capsicum frutescens and Capsicum pubescens, corresponding to a total of 38.1 million tons of hot peppers, are consumed by three-quarters of the world’s population (FAO, see www.fao.org) . Capsorubin and capsanthin can be described as carotenoids and are exclusively biosynthesized in peppers. Peppers are beneficial for preventing various diseases, such as eye diseases, cardiovascular diseases and certain cancers [2,3,4], and provide excellent natural colourants . The accumulation of carotenoid pigments results in colour variations in the ripe pepper fruit. The carotenoids biosynthetic pathway is currently known (Figure S1A) [6, 7]. Capsorubin or capsanthin is transformed by geranylgeranyl pyrophosphate (GGPP) in a series of enzymatic reactions with phytoene synthase (PSY), phytoene desaturase (PDS), lycopene β-cyclase (LCYB), etc. [6, 7]. Moreover, capsaicinoids, the alkaloids specific to Capsicum, are responsible for pungency. Capsaicinoids are produced in the fruit placenta of hot peppers . Capsaicinoids mainly contain capsaicin, dihydrocapsaicin and several analogues [9, 10]. Capsaicin and dihydrocapsaicin comprise approximately 91% of the total capsaicinoid content in Capsicum species . Numerous studies have demonstrated that capsaicinoids are produced after the condensation of phenylalanine and in fatty acid chain biosynthetic pathways (Figure S1B) [1, 12, 13].
The bHLH family is one of the largest transcription factors (TFs) in plants and is important for plant growth and development [14, 15]. bHLH proteins were originally found in eukaryotic species  and contain a basic and a helix–loop–helix (HLH) region. The two regions possess DNA-binding and protein-interacting abilities . The basic region contains 13–17 basic amino acids in the N-terminal domain and provides a DNA-binding region to bind to consensus hexanucleotide E-box (CANNTG) [17, 18], whereas the HLH region includes approximately 40 amino acids in the C-terminal domain with two alpha helixes separated by a loop of variable length [16, 19]. Furthermore, the HLH region promotes the interaction with other bHLH proteins and the formation of homodimers and heterodimeric complexes . Except for the two conserved regions, the bHLH protein sequences are dissimilar, revealing protein evolution .
The bHLH TFs have been identified in the majority of plants, including Arabidopsis thaliana, rice, and potato [15, 21,22,23,24,25,26,27,28,29,30]. WRKY, NAC and APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) TF families have been analysed systematically in peppers [31,32,33,34]. The bHLH TFs are closely related to the primary and specialized metabolites of the plant. PHYTOCHROME INTERACTING FACTOR 3 (PIF3), a bHLH TF, positively regulates anthocyanin biosynthesis by activating the transcription of anthocyanin biosynthetic genes in Arabidopsis thaliana . SlAN1 and SlPIF1a regulate anthocyanin and carotenoid biosynthesis in tomato, whereas SlPRE2 negatively regulates pigment accumulation in fruit [36,37,38]. The bHLH proteins typically regulate secondary metabolite biosynthesis by interacting with TFs such as MYB, WD40-repeat (WDR), and ETHYLENE RESPONSIVE FACTOR (ERF) [39,40,41,42,43]. bHLH TFs physically interact with MYB proteins and facilitate their regulation of anthocyanin biosynthesis in plants. PacMYBA physically interacts with several anthocyanin-related bHLH TFs controlling anthocyanin biosynthesis in sweet cherry by activating the promoters of PacDFR, PacANS and PacUFGT . AtTT8 interacts with MYB proteins PRODUCTION OF ANTHOCYANIN PIGMENT1/PRODUCTION OF ANTHOCYANIN PIGMENT (PAP1/PAP) and TRANSPARENT TESTA2 (TT2) to regulate anthocyanin and proanthocyanidin biosynthesis in Arabidopsis thaliana [44, 45]. Rosea1 (ROS1, an MYB type) and Delila (DEL, a bHLH type) collectively regulate anthocyanin biosynthesis in tomato . In peppers, CaMYC combined with CaMYB and CaWD40 regulates anthocyanin biosynthesis through the regulation of the transcription of a synthetic gene . The physical interaction and regulatory synergy between MYB and bHLH TFs are mediated by the R3 domain in MYB proteins and the N-terminal region in bHLH proteins [45, 47,48,49]. Additionally, bHLH-mediated ERF TFs regulate synthesis of alkaloids, such as terpenoid indoles alkaloids , glycoside alkaloids , and nicotine . Capsorubin and capsaicinoids are strictly biosynthesized in Capsicum, which possesses high economic and nutritional value. Our previous studies confirmed that MYB31, MYB108 and MYB48 were involved in capsaicinoid biosynthesis [1, 12, 50, 51]. Specifically, natural variations of the master regulator MYB31 promoter can affect capsaicinoid contents among different Capsicum species [1, 50]. However, whether bHLH TFs regulate capsorubin and capsaicinoid biosynthesis in peppers and orchestrate the regulation of capsaicinoid biosynthesis by interacting with MYB or other TFs needs to be addressed.
In this study, bHLH family genes were identified in the Capsicum annuum genome. A reference genome of C. annuum was sequenced (2n = 2x = 24), and its genome size was estimated to be 3.48 Gb by 19-mer analysis . Characteristic analysis of bHLH family members was systematically performed for better understanding the potential biological function of candidate bHLH TFs in regulating species-specific metabolite biosynthesis in peppers. Capsaicinoids biosynthesis is affected by environmental factors, such as water, temperature and light . Higher temperatures are beneficial to the accumulation of total capsaicinoids . bHLH TFs respond to temperature in Camellia sinensis  and Arabidopsis thaliana . Thus, the candidate bHLH TFs were also studied in response to different temperatures. Moreover, the interaction between bHLH TFs and MYB31 was analysed by Y2H. This study provides candidate bHLH TFs in carotenoid and capsaicinoid biosynthetic pathways in peppers.
Identification and chromosomal localizations of bHLHs in pepper
After excluding redundant sequences, a total of 107 bHLH proteins were identified from the C. annuum genome using the Hidden Markov Model (HMM) profile of the HLH domain (PF00010). As shown in Table S1, all identified CabHLH proteins encoded 190–940 amino residues. The molecular weight (Mw) of these proteins ranged from 21.61 kDa to 106.70 kDa, and the theoretical pI ranged from 4.60 to 9.91.
The 107 CabHLHs were renamed as CabHLH001 to CabHLH107 and mapped to the pepper chromosomes based on the chromosomal positions (Figure S2). They were distributed on 12 chromosomes. Chromosome 01 contained the largest number of CabHLH members (18). Chromosomes 05 and 09 included three CabHLH members, separately. However, the genes ranging from CabHLH092 to CabHLH107 were not located on any chromosomes. Genes located in some scaffolds were not assembled to chromosomes due to sequencing and assembly technology limitations.
Conserved amino acid residues in the bHLH domain
The amino acid sequences of the bHLH domain were used to perform multiple alignment analyses (Figure S3). The results indicated that bHLH family proteins possessed a conserved bHLH domain, which contained the basic, first helix, loop and second helix regions. As shown in Fig. 1, 20 amino acid residues were conserved with a consensus ratio greater than 50%, and six amino acid residues were conserved with a consensus ratio greater than 75% among the conserved bHLH domains. Five residues (His-5, Glu-9, Arg-10, Arg-12 and Arg-13) were conserved in the basic region. Four residues (Leu-23, Leu-26, Val-27 and Pro-28) were conserved in the first helix region. Lys-31 and Asp-34 were conserved in the loop region, and nine residues (Lys-35, Ala-36, Ser-37, Leu-39, Ala-42, Ile-43, Tyr-45, leu-49 and leu-56) were conserved in the second helix region. The residues Leu-23 and Leu-49 were extremely conserved among the 107 bHLH proteins in pepper.
Phylogenetic analysis of the bHLH family proteins
To classify the CabHLH proteins, a phylogenetic tree that contained all of the identified bHLH protein sequences in peppers and those in Arabidopsis thaliana was constructed with the neighbour-joining method (Fig. 2). According to the classification of AtbHLHs in a previous study , the CabHLH proteins were divided into 15 subfamilies and named groups I to XII. Group II contained the largest numbers of CabHLHs (25) and AtbHLHs (13), whereas group VII contained only one CabHLH and three AtbHLH proteins. The different number of CabHLHs and AtbHLHs proteins in the same group might arise from unequal duplication of the bHLH family during the plant's evolutionary process.
Members of the same group might possess similar biological functions. To preliminarily speculate on the biological functions of CabHLHs, another neighbour-joining phylogenetic tree was constructed using all bHLH proteins in Arabidopsis thaliana, tomato, rice and pepper (Figure S4). The functional characteristics of AtbHLHs and SlbHLHs have been reported and summarized in the literature, and these functional characteristics were used to evaluate the potential function of CabHLH in the same group (Table S2). SlPIF1a included in group VI could regulate carotenoid biosynthesis by a light-dependent mechanism in tomato . SlPIF1a mapped to CabHLH051 in peppers and was also classified into group VI. These results indicated that the members of group VI might be involved in carotenoid biosynthesis.
Analysis of conserved motifs in CabHLH proteins
To investigate the structural features of CabHLH proteins, the conserved amino acid motifs were analysed and identified using the Multiple EM for Motif Elicitation (MEME) Suite. A total of 15 conserved motifs containing 21–100 residues were found in motifs 1 to 15. The motifs information is provided in Table S3. Motifs 1 and 2 were located in the bHLH domain region, which appeared in all proteins. Motifs 3 to 15 were distributed outside the bHLH domain region. Motifs 6, 7, 8 and 11 were primarily restricted to group II. Motifs 9, 12 and 15 were specifically identified within group IV. Motifs 10 and 13 were found in group Id, and motif 14 was observed in group V. Generally, most of the same group’s proteins possessed common motifs in terms of alignment and position (Fig. 3).
Expression profiles of CabHLHs at different developmental stages in pericarp and placenta
Capsorubin and capsaicinoids are found in pepper fruit and closely associated with the biosynthetic genes transcriptional level during developmental stages. To obtain further insight into the potential functions of CabHLHs during capsorubin and capsaicinoid biosynthesis, the expression profiles of CabHLHs in different developmental pericarp and placenta were investigated. RNA-Seq raw data were obtained from Kim et al.  and included 6, 16, 25, 36, 38, 43 and 48 days post-anthesis (DPA) stages (Figs. 4 and 5). All the raw reads were spliced and remapped to version 2.0 of the C. annuum genome.
As shown in Fig. 4b, the expression levels of capsorubin biosynthetic genes gradually increased at 36 DPA, which was consistent with the accumulation profile of capsorubin in pericarp tissue. A total of 20 expressed CabHLH genes were not detectable. These genes were probably transcribed at a low level in different developmental pericarp stage. According to the similarity of the expression profiles, all CabHLH expression profiles in different developmental pericarp stages were hierarchically clustered and classified into 10 clusters (Fig. 4a). The expression profiles of CabHLHs in clusters C1 to C4 maintained good agreement with the expression profiles of capsorubin biosynthetic genes. The members of clusters C1 to C4 might be associated with capsorubin biosynthesis.
Capsaicinoids were abundantly produced from 13 to 25 DPA in the placenta's developmental stages, and the expression levels of capsaicinoid biosynthetic genes were high during this stage (Fig. 5b). Based on the similarity of the expression profiles, the expression of all the CabHLHs in different developmental placenta stages were hierarchically clustered into 10 clusters (Fig. 5a). A total of 16 genes were expressed at a low level that could not be detected. The expression profiles of CabHLHs in clusters L5, L6, L8 and L9 were similar to the expression profiles of capsaicinoid biosynthetic genes. Therefore, these results indicated that the members of clusters L5, L6, L8 and L9 might be associated with capsaicinoids biosynthesis.
Additionally, capsorubin and capsaicinoids are produced mainly in the pericarp and placental tissues of peppers. To confirm whether CabHLHs were specifically expressed in pericarp and placental tissues, the expression profiles of all identified CabHLHs in different tissues, including the leaf, root, stem, pericarp and placenta, were investigated. However, the RNA-Seq raw data uploaded by Kim et al. did not contain leaf, root or stem tissues from peppers . The RPKM values of these tissues, which were mapped to version 1.5 of the C. annuum genome, were directly published online. The heat map indicated that CabHLHs were not specifically expressed in certain tissues (Figure S5). Presumably, these TFs orchestrated functions in addition to regulating capsorubin and capsaicinoid biosynthesis.
Validation of candidate bHLH genes involved in capsorubin and capsaicinoids biosynthesis
To further verify the expression profiles of CabHLHs in the pericarp and placenta, ten CabHLHs from candidate clusters that might be related to capsorubin and capsaicinoid biosynthesis were selected for qRT-PCR analysis. The selected genes were relatively highly expressed in different developmental pericarp or placental tissues. As shown in Fig. 6a, the contents of zeaxanthin and capsorubin increased from the MG stage in pericarp tissue, whereas the lutein content from the branch of the non-synthetic capsorubin decreased. No significant difference in β-carotene content was noted in different developmental pericarp tissues. CabHLH032, CabHLH048, CabHLH095 and CabHLH100 expression was consistent with the accumulation profile of carotenoids (zeaxanthin and capsorubin) in pericarp, whereas CabHLH009 expression was similar to the accumulation profile of lutein in pericarp. However, these genes were also highly expressed in other tissues (roots, flowers, stems, placentas, leaves and seeds) (Fig. 6b). Thus, it was likely that the members of clusters C1, C2 and C3 were associated with capsorubin biosynthesis but also possessed other specific functions in other tissues.
The capsaicin and dihydrocapsaicin were initially produced at 10 DPA, peaked at 25 DPA in placental tissue, and then gradually decreased (Fig. 7a). CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 expression resembled the accumulation profile of capsaicin in different developmental placenta stages. A high expression level of CabHLH026 was observed in pericarp, seed and placental tissues, and high expression of CabHLH063 was evident in stem and leaf. CabHLH007, CabHLH009 and CabHLH086 were highly expressed in certain tissues (Fig. 7b). Therefore, CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 in clusters L5, L6 and L9 might be associated with capsaicinoid biosynthesis.
The expression of candidate CabHLHs associated with capsaicinoid biosynthesis in response to different temperatures
To obtain a preliminary understanding of whether capsaicinoid biosynthesis is regulated by CabHLH genes in response to different temperatures, the expression of five candidate CabHLHs and five important capsaicinoid biosynthetic genes at different temperatures was measured. As shown in Fig. 8a, the capsaicin and dihydrocapsaicin content tended to increase with increasing temperature from T15 to T25. The placenta capsaicin content with T25 treatment was significantly increased compared with that in peppers with T33 treatment. CabHLH007, CabHLH009 and CabHLH086 expression increased with increasing temperatures, which was similar to the accumulation of dihydrocapsaicin (Fig. 8b). CabHLH026 was highly expressed in T25, which was consistent with the accumulation of capsaicin and the expression of capsaicinoid biosynthetic genes (AT3, AMT, BCKDH and KasIa) (Fig. 8b; c). In contrast, CabHLH063 expression decreased with increasing temperature, which maintained consistency with the expression profile of the capsaicinoid biosynthetic gene Acl (Fig. 8b; c). Thus, CabHLH007, CabHLH009 and CabHLH086 expression was positively associated with an increase in dihydrocapsaicin content and temperature. CabHLH063 expression was negatively related to the increase in the capsaicin content and temperature, whereas CabHLH026 expression was positively related to these factors. These candidate genes might be related to capsaicinoid biosynthesis in response to different temperatures by regulating the transcription of capsaicinoid biosynthetic genes.
The interaction of candidate CabHLHs and MYB31 in yeast and identified bHLH binding sites
bHLH proteins typically bind to E-box binding sites of gene promoter and regulate the transcription. To characterize the potential association of bHLH TFs and the pathway biosynthetic genes, one thousand five hundred base pair nucleotide sequences upstream of the start codon (ATG) from capsorubin and capsaicin biosynthetic genes, including CCS, PSY, β-CH, β-LCY, Acl, AMT, AT3, BACT, BCKDH, COMTa, FatA and KasIa, were analysed using the PlantCARE database . As shown in Table S7, multiple bHLH DNA binding sites were detected in the promoters of capsorubin and capsaicin biosynthetic genes, except for CCS, β-CH and KasIa.
The bHLH protein executes function always by interacting with other TFs, such as MYB. We performed Y2H assay to verify the interaction between candidate CabHLHs and CaMYB31. The results indicated that these bHLHs interacted with MYB31 in a gene-dependent manner. CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 could interact with MYB31 in a heterologous system. CabHLH026 strongly interacts with CaMYB31 in yeast, whereas only weak interactions were observed in the group of CabHLH063-CaMYB31 and CabHLH086-CaMYB31 (Fig. 9). Therefore, it was likely that CabHLH could regulate capsaicinoid biosynthesis by interacting with CaMYB31, a master regulator of capsaicinoid biosynthesis.
The bHLH family has emerged as the second-largest TF family in plants . This family has been successfully identified and investigated in many plants, including Solanum tuberosum L. (124) , Camellia sinensis (L.) O. Ktze. (120) , Solanum lycopersicum (159) , Zea mays L. (208) , Brassica rapa pekinensis (230) , Arabidopsis thaliana (147) , Glycine max (L) Merr. (319) , Malus pumila (188)  and Oryza sativa L. (167) . In this study, 107 CabHLH genes were identified in the pepper genome. The conserved bHLH domain consisted of the basic, first helix, loop and second helix regions, which contained approximately 60 amino acids ; this structure was also observed in peppers (Figure S3). Twenty amino acid residues were conserved with a consensus ratio greater than 50% (Fig. 1), consistent with a previous study . Glu-9, Arg-12 and Arg-13 were conserved with a consensus ratio greater than 80%, whereas Leu-23 and Leu-49 were extremely conserved among 107 CabHLH proteins (Fig. 1). A previous study had confirmed that Glu-13 and Arg-16/Arg-17 (corresponding to Glu-9 and Arg-12/Arg-13 in this study; Fig. 1) played an important role in DNA binding . Leu-27 and Leu-61 (corresponding to Leu-23 and Leu-49 in this study; Fig. 1) are necessary for dimerization . Glu-13 and Arg-16 (corresponding to Glu-9 and Arg-12 in this study; Fig. 1) recognize the E-box (5’-CANNTG-3’), whereas the His/Lys-9, Glu-13 and Arg-17 (corresponding to His-5, Glu-9 and Arg-13 in this study; Fig. 1) were required for binding G-Box (5’-CACGTG-3’) [19, 59].
According to a previous study , the phylogenetic tree was divided into 15 subgroups (Fig. 2), which was similar to the classification of potato . Motifs 1 and 2 located in the bHLH domain were included in all CabHLH proteins (Fig. 3), which was consistent with the observation of MdbHLH proteins . Motifs 1 and 2 containing some conserved amino residues played important roles in DNA binding and protein dimerization [23, 58]. The motifs from 3 to 15 were distributed outside the bHLH domain and randomly arranged in CabHLH proteins (Fig. 3). The results indicated that the members of the same group had the same motif arrangements in most case.
Homologous genes generally possessed a similar function in plant growth and development processes . Some SlbHLHs and AtbHLHs included in groups II and VI were functionally identified (Table S2). MYC2 interacted with MYB TFs to regulate glucosinolate biosynthesis  and respond to abscisic acid , jasmonic acid , and light signalling  in Arabidopsis thaliana. TRANSPARENT TESTA8 (TT8), GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) combined with MYB and WDR TFs could regulate anthocyanin biosynthesis and trichome development in Arabidopsis thaliana [65,66,67]. These proteins were in group II (Figure S4). Some SlbHLHs in group II were also involved in anthocyanin biosynthesis  and trichome development  and responded to cold, osmotic and salt stress . Furthermore, MYC2 interacted with ERF TFs to regulate the biosynthesis of alkaloids such as steroidal glycoalkaloids , terpenoid indole alkaloids , and nicotine biosynthesis . The homologous genes of MYC2 in pepper were mapped to CabHLH006 and CabHLH018, and all of the proteins were classified into group II (Figure S4). These results indicated that group II might be related to secondary metabolites, development and signal response processes in plants. In addition, PIF3 combined with bZIP TFs could regulate anthocyanin biosynthesis . PIFs served as cellular signalling hubs that function by integrating multiple signals to regulate the transcriptional network during plant growth and development . SlPIF1a was included in group VI and regulated carotenoid biosynthesis in tomato . Homologous genes of these PIFs were mapped to CabHLH063 and CabHLH051 in pepper, and the proteins were located in group VI. The CabHLH protein members in groups II and VI probably regulated secondary metabolite biosynthesis by interacting with other TFs.
The expression profiles of CabHLHs in clusters C1, C2, C3 and C4 maintained good agreement with the expression of capsorubin biosynthetic genes (Fig. 4). The expression profiles of CabHLH009, CabHLH032, CabHLH048, CabHLH095 and CabHLH100 from these clusters aligned with the accumulation profile of carotenoids, including zeaxanthin, lutein and capsorubin, in pericarp (Fig. 6). Therefore, it was likely that the members of these clusters were related to capsorubin biosynthesis. The expression profiles of CabHLHs in clusters L5, L6, L8 and L9 were similar to the expression profiles of capsaicinoid biosynthetic genes (Fig. 5). CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 expression profiles were consistent with the accumulation profile of capsaicinoids in placenta (Fig. 7). These results indicated that the members of clusters L5, L6 and L9 might regulate capsaicinoid biosynthesis. Moreover, the similar transcript levels of the genes indicated that they might perform similar functions in some cases. CaMYC combined with CaMYB and CaWD40 to regulate anthocyanin biosynthesis  and was mapped to CabHLH0016 in cluster C1 (Fig. 4). Numerous studies have reported that MYC2 typically regulated alkaloid biosynthesis by interacting with ERF TFs [41,42,43], and ERF TFs were associated with capsaicinoid biosynthesis in Capsicum [34, 71]. MYC2 was mapped to two homologous genes CabHLH006 and CabHLH018 that were classified into clusters L8 and L7, respectively. The expression profile of cluster L8 maintained good agreement with the profile of capsaicinoid biosynthesis. SlPIF1a was related to carotenoid biosynthesis in tomato  and mapped to CabHLH051 in the pepper, which was classified into cluster C4. The members of cluster C4 were candidate regulators for carotenoid biosynthesis in this study. The members of clusters C1, C2, C3, C4, L5, L6, L8 and L9 were probably related to carotenoid and capsaicinoid biosynthesis.
Environmental factors are essential for secondary metabolite biosynthesis. Previous studies have demonstrated that bHLH TFs were responsive to a heat response in Camellia sinensis and Arabidopsis thaliana [26, 72]. PIFs acted as a signalling centre that integrated multiple signals to regulate plant growth and development, including high temperature and light . The expression of CabHLH063, a homologous gene of PIF, and Acl decreased with increasing temperature. The expression of CabHLH026 and four capsaicinoid biosynthetic genes (i.e., AT3, AMT, BCKDH and KasIa) peaked in T25, which was similar to capsaicin accumulation (Fig. 8c). These results indicated that the members of clusters L5, L6 and L9 might possess different functions in regulating capsaicinoid biosynthesis in response to different temperatures. In addition, bHLH proteins were typically involved in plant metabolic pathways by interacting with other TFs [39, 41,42,43,44,45,46]. MYB TFs played a key role in capsaicinoid biosynthesis as evidenced in previous studies [1, 12, 50, 51]. Y2H assays revealed that CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 could interact with MYB31 in a heterologous system. These results implied that CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 might be involved in capsaicinoid biosynthesis by interacting with MYB31.
A total of 107 CabHLH proteins were identified in peppers. These proteins were divided into 15 groups according to the classification of Arabidopsis thaliana. The bHLH conserved domain containing motifs 1 and 2 appeared in all CabHLHs. The expression profiles showed that clusters C1, C2, C3, C4, L5, L6, L8 and L9 were candidates for the regulation of carotenoid and capsaicinoid biosynthesis. CabHLH009, CabHLH032, CabHLH048, CabHLH095, CabHLH100, CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 were selected from the candidate clusters because they might contribute to carotenoid and capsaicinoid biosynthesis. CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 might be responsive to different temperatures to mediate capsaicinoid biosynthesis. CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 likely control capsaicinoid biosynthesis by interacting with MYB31. However, further studies showing how these candidate bHLHs regulate carotenoid and capsaicinoid biosynthesis are required.
Identification of bHLHs in peppers and their chromosomal locations
The bHLH protein sequences were retrieved from C. annuum ‘CM334’ genome  using the HMM profile of the HLH domain (PF00010) obtained from the PFAM database . Redundant sequences were filtered with HMMER software , SMART database  and NCBI Conserved Domain Search Service (CD Search) . The ExPASy server  was used to predict the full length of amino acid sequences, MW, PI and instability index of the proteins. The corresponding information is provided in Table S1. The chromosomal locations of gene loci were acquired from version 2.0 of the pepper genome. According to the chromosomal positions, all the identified CabHLHs were renamed consecutively.
Multiple alignments and phylogenetic analysis
Multiple sequence alignments were performed using conserved bHLH domain sequences of four plant species, including peppers, rice, tomato and Arabidopsis thaliana. The alignments were performed with ClustalX 2.1 with the default parameters. The bHLH protein sequences in Arabidopsis thaliana (131) were obtained from The Arabidopsis Information Resource (TAIR) database , whereas the bHLH protein sequences in rice (144) and tomato (145) were acquired from Plant TF Database version 4.0 . All of these bHLH protein sequences were renamed and listed in Table S4. The neighbour-joining phylogenetic trees were generated using MEGA X with 1000 bootstrap replications  and visualized using EvolView v3 .
Protein motif analysis
RNA-Seq data and expression profiles of CabHLHs
FPKM values of candidate CabHLHs were obtained from the pepper RNA-Seq raw data (GenBank: AYRZ00000000) [84, 85]. Heat maps showing the gene expression profiles in seven different developmental placenta and pericarp tissues (6 DPA to 48 DPA) were drawn by R language.
Plant materials and temperature treatments
The 59 inbred line (Capsicum annuum) authenticated by Professor Jianjun Lei (College of Horticulture, South China Agricultural University) was used for the experiments . The variety rights of this breed belong to our lab. If somebody needs this inbred line, they can acquire from the corresponding authors. All plant materials used in this study were provided by South China Agricultural University. According to institutional, national and international guidelines, these samples do not require specific permissions for research purposes. The seeds were sown using mixture substrate (peat, perlite and coir pith) in the plastic greenhouse under the following environment conditions: 35℃/23℃(day/night); 12 h/12 h (light/dark cycle), and 300–1000 μmol·m−2·s−1 of photosynthetic photon flux density (PPFD) at 12:00 AM. The sample used for RNA extraction was collected according to the methods of our previous study .
Peppers used for temperature treatment experiment were also cultivated using mixture substrate as described above. The temperature treatment and sample collection were performed as previously described . All the samples were performed with three biological replicates and one technical replicate for each biological replicate.
Quantitative Real-time PCR (qRT-PCR) analysis
The isolation and extraction of total RNA and RNA reverse transcription were performed according to Song et al. . The primers for qPCR were designed using Primer 5.0 (Table S5). Ploy (A)-binding protein (CA00g52140) and ubiquitin extension protein (CA12g20490) were used as housekeeping genes, given that transcriptome studies revealed that their transcription levels were fairly stable among different tissues and plant different development stages in Capsicum . qRT-PCR was performed according to Wang et al. . The relative expression of each CabHLH was calculated using the 2−ΔΔCt method . Three technical replicates of each biological replicate and three biological replicates were performed for all samples. The Dunnett’s t-test was used to determine the results of significant differences by using SPSS 22.
Yeast two-hybrid assay
The full-length cDNAs of CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 were cloned into a pMD19-T vector (6013, TaKaRa, China). The amplified full-length fragments of CabHLH007, CabHLH009, CabHLH026, CabHLH063 and CabHLH086 were ligated into the pGADT7 and pGBKT7 vectors using one-step cloning (C112, Vazyme, China), separately. Negative and positive plasmids and plasmids containing different CabHLH were transformed into the AH109 yeast strain (YC1010, Weidi Biotechnology, China). The Y2H assay was performed according to the manufacturer’s instructions (Clontech). Image processing was performed using Adobe Illustrator CS2020. All the primers were listed in Table S6.
Promoter cis-element discovery
The basis region of the bHLH protein could bind to E-box (CANNTG). One thousand five hundred base pair nucleotide sequences upstream of the start codon in the biosynthetic genes of capsorubin (CCS, PSY, β-CH, β-LCY) and capsaicin (Acl, AMT, AT3, BACT, BCKDH, CoMTa, FatA and KasIa) were retrieved from C. annuum genome and applied to E-box discovery. The cis-elements were identified by using the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).
Availability of data and materials
Most data generated or analysed during this study are included in this article and its supplemental files. The sequencing data (GenBank: AYRZ00000000) used and analysed during this study are available in the NCBI database (doi:10.1038/ng.2877).
Fragments per kilobase of transcript per million mapped reads
Multiple EM for Motif Elicitation
Hidden Markov model
Caffeic acid O-methyltransferase
A putative acyltransferase
p-Coumaroyl shikimate/quinate 3-hydroxylase
Hydroxycinnamoyl-CoA hydratase lyase
Branched-chain amino acid aminotransferase
Acyl carrier protein
Mature green stage
- B + 7:
Breaker plus 7 days
- B + 14:
Breaker plus 14 days
Ethylene responsive factor
Production of anthocyanin pigment1/production of anthocyanin pigment
Apetala2/ethylene response factor
Phytochrome interacting factor 3
Enhancer of glabra3
Yeast two-hybrid assays
photosynthetic photon flux density
Zhu Z, Sun B, Cai W, Zhou X, Mao Y, Chen C, Wei J, Cao B, Chen C, Chen G, Lei J. Natural variations in the MYB transcription factor MYB31 determine the evolution of extremely pungent peppers. New Phytol. 2019;223(2):922–38.
Sharma SK, Vij AS, Sharma M. Mechanisms and clinical uses of capsaicin. Eur J Pharmacol. 2013;720(1–3):55–62.
Spiller F, Alves MK, Vieira SM, Carvalho TA, Leite CE, Lunardelli A, Poloni JA, Cunha FQ, de Oliveira JR. Anti-inflammatory effects of red pepper (Capsicum baccatum) on carrageenan- and antigen-induced inflammation. J Pharm Pharmacol. 2008;60(4):473–8.
Melendez-Martinez AJ. An overview of carotenoids, apocarotenoids, and vitamin A in agro-food, nutrition, health, and disease. Mol Nutr Food Res. 2019;63(15):e1801045.
Gomes LM, Petito N, Costa VG, Falcao DQ, de Lima AK. Inclusion complexes of red bell pepper pigments with beta-cyclodextrin: preparation, characterisation and application as natural colorant in yogurt. Food Chem. 2014;148:428–36.
Berry HM, Rickett DV, Baxter CJ, Enfissi E, Fraser PD. Carotenoid biosynthesis and sequestration in red chilli pepper fruit and its impact on colour intensity traits. J Eep Bot. 2019;70(10):2637–50.
Guzman I, Hamby S, Romero J, Bosland PW, O’Connell MA. Variability of carotenoid biosynthesis in orange colored Capsicum spp. Plant Sci. 2010;179(1–2):49–59.
Fujiwake H, Suzuki T, Iwai K. Intracellular localization of capsaicin and its analogues in Capsicum fruit II. The vacuole as the intracellular accumulation site of capsaicinoid in the protoplast of Capsicum fruit. Plant Cell Physiol. 2006;21(6):1023–30.
Reeves G, Bosland PW, Coon D. ‘Trinidad Moruga Scorpion’ pepper is the world’s hottest measured chile pepper at more than two million scoville heat units. Horttechnology. 2012;22(4):534–8.
Carle R, Schieber A, Schweiggert U. Characterization of major and minor capsaicinoids and related compounds in chili pods (Capsicum frutescens L.) by high-performance liquid chromatography/atmospheric pressure chemical ionization mass spectrometry. Anal Chim Acta. 2006;557(1/2):236–44.
CJ Stewart BC Kang K Liu M Mazourek SL Moore EY Yoo BD Kim I Paran MM Jahn 2005 The Pun1 gene for pungency in pepper encodes a putative acyltransferase Plant J 42 5 675 688
Sun B, Zhu Z, Chen C, Chen G, Cao B, Chen C, Lei J. Jasmonate-inducible R2R3-MYB transcription factor regulates capsaicinoid biosynthesis and stamen development in Capsicum. J Agric Food Chem. 2019;67(39):10891–903.
Choi H, Yu Y, Jo YD, Ahn HK, Park B, Yeam I, Kim K, Park J, Choi YD, Park M. Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat Genet. 2014;46(3):270–8.
Ledent V, Vervoort M. The basic helix-loop-helix protein family: comparative genomics and phylogenetic analysis. Genome Res. 2001;11(5):754–70.
Toledo-Ortiz G, Huq E, Quail PH. The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003;15(8):1749–70.
Murre C, McCaw PS, Baltimore D. A new DNA binding and dimerization motif in immunoglobulin enhancer binding, daughterless, MyoD, and myc proteins. Cell. 1989;56(5):777–83.
Atchley WR, Terhalle W, Dress A. Positional dependence, cliques, and predictive motifs in the bHLH protein domain. J Mol Evol. 1999;48(5):501–16.
Massari ME, Murre C. Helix-loop-helix proteins: regulators of transcription in eucaryotic organisms. Mol Cell Biol. 2000;20(2):429–40.
Ferre-D’Amare AR, Pognonec P, Roeder RG, Burley SK. Structure and function of the b/HLH/Z domain of USF. Embo J. 1994;13(1):180–9.
Morgenstern B, Atchley WR. Evolution of bHLH transcription factors: modular evolution by domain shuffling? Mol Biol Evol. 1999;16(12):1654–63.
Heim MA, Jakoby M, Werber M, Martin C, Weisshaar B, Bailey PC. The basic helix-loop-helix transcription factor family in plants: a genome-wide study of protein structure and functional diversity. Mol Biol Evol. 2003;20(5):735–47.
Li X, Duan X, Jiang H, Sun Y, Tang Y, Yuan Z, Guo J, Liang W, Chen L, Yin J. Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis. Plant Physiol. 2006;141(4):1167–84.
Carretero-Paulet L, Galstyan A, Roig-Villanova I, Martinez-Garcia JF, Bilbao-Castro JR, Robertson DL. Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in Arabidopsis, poplar, rice, moss, and algae. Plant Physiol. 2010;153(3):1398–412.
Wang R, Zhao P, Kong N, Lu R, Pei Y, Huang C, Ma H, Chen Q. Genome-wide identification and characterization of the potato bHLH transcription factor family. Genes. 2018;9(1):54.
Wang P, Su L, Gao H, Jiang X, Wu X, Li Y, Zhang Q, Wang Y, Ren F. Genome-wide characterization of bHLH genes in grape and analysis of their potential relevance to abiotic stress tolerance and secondary metabolite biosynthesis. Front Plant Sci. 2018;9:64.
Cui X, Wang YX, Liu ZW, Wang WL, Li H, Zhuang J. Transcriptome-wide identification and expression profile analysis of the bHLH family genes in Camellia sinensis. Funct Integr Genomics. 2018;18(5):489–503.
Song XM, Huang ZN, Duan WK, Ren J, Liu TK, Li Y, Hou XL. Genome-wide analysis of the bHLH transcription factor family in Chinese cabbage (Brassica rapa ssp. pekinensis). Mol Genet Genomics. 2014;289(1):77–91.
Zhang T, Lv W, Zhang H, Ma L, Li P, Ge L, Li G. Genome-wide analysis of the basic Helix-Loop-Helix (bHLH) transcription factor family in maize. BMC Plant Biol. 2018;18(1):235.
Wang J, Hu Z, Zhao T, Yang Y, Chen T, Yang M, Yu W, Zhang B. Genome-wide analysis of bHLH transcription factor and involvement in the infection by yellow leaf curl virus in tomato (Solanum lycopersicum). BMC Genomics. 2015;16(1):39.
Sun H, Fan HJ, Ling HQ. Genome-wide identification and characterization of the bHLH gene family in tomato. BMC Genomics. 2015;16(1):9.
Diao W, Snyder JC, Wang S, Liu J, Pan B, Guo G, Ge W, Dawood M. Genome-wide analyses of the NAC transcription factor gene family in pepper (Capsicum annuum L.): chromosome location, phylogeny, structure, expression profiles, cis-elements in the promoter, and interaction network. Int J Mol Sci. 2018;19(4):1048–61.
Jin JH, Wang M, Zhang HX, Khan A, Wei AM, Luo DX, Gong ZH. Genome-wide identification of the AP2/ERF transcription factor family in pepper (Capsicum annuum L.). Genome. 2018;61(9):663–74.
Cheng Y, Yao ZP, Ruan MY, Ye QJ, Wang RQ, Zhou GZ, Luo J. In silico identification and characterization of the WRKY gene superfamily in pepper (Capsicum annuum L.). Genet Mol Res. 2016;15(3):1–12.
Song J, Chen C, Zhang S, Wang J, Huang Z, Chen M, Cao B, Zhu Z, Lei J. Systematic analysis of the Capsicum ERF transcription factor family: identification of regulatory factors involved in the regulation of species-specific metabolites. BMC Genomics. 2020;21(1):573.
Shin J, Park E, Choi G. PIF3 regulates anthocyanin biosynthesis in an HY5-dependent manner with both factors directly binding anthocyanin biosynthetic gene promoters in Arabidopsis. Plant J. 2007;49(6):981–94.
Park J, Nou I, Kim H, Afrin KS, Kang S, Rahim MA, Jung H. Expression of anthocyanin biosynthesis-related genes reflects the peel color in purple tomato. Hortic Environ Biote. 2018;59(3):435–45.
Zhu Z, Chen G, Guo X, Yin W, Yu X, Hu J, Hu Z. Overexpression of SlPRE2, an atypical bHLH transcription factor, affects plant morphology and fruit pigment accumulation in tomato. Sci Rep. 2017;7(1):5786.
Llorente B, D’Andrea L, Ruiz-Sola MA, Botterweg E, Pulido P, Andilla J, Loza-Alvarez P, Rodriguez-Concepcion M. Tomato fruit carotenoid biosynthesis is adjusted to actual ripening progression by a light-dependent mechanism. Plant J. 2016;85(1):107–19.
Shen X, Zhao K, Liu L, Zhang K, Yuan H, Liao X, Wang Q, Guo X, Li F, Li T. A role for PacMYBA in ABA-regulated anthocyanin biosynthesis in red-colored sweet cherry cv. Hong Deng (Prunus avium L.). Plant Cell Physiol. 2014;55(5):862–80.
Lu B, Cheng G, Zhang Z, Sun J, Ali M, Jia Q, Luo D, Gong Z, Li D. CaMYC, a novel transcription factor, regulates anthocyanin biosynthesis in color-leaved pepper (Capsicum annuum L.). J Plant Growth Regul. 2019;38(2):574–85.
Cardenas PD, Sonawane PD, Pollier J, Vanden BR, Dewangan V, Weithorn E, Tal L, Meir S, Rogachev I, Malitsky S. GAME9 regulates the biosynthesis of steroidal alkaloids and upstream isoprenoids in the plant mevalonate pathway. Nat Commun. 2016;7:10654.
Zhang H, Hedhili S, Montiel G, Zhang Y, Chatel G, Pre M, Gantet P, Memelink J. The basic helix-loop-helix transcription factor CrMYC2 controls the jasmonate-responsive expression of the ORCA genes that regulate alkaloid biosynthesis in Catharanthus roseus. Plant J. 2011;67(1):61–71.
De Boer K, Tilleman S, Pauwels L, Vanden BR, De Sutter V, Vanderhaeghen R, Hilson P, Hamill JD, Goossens A. APETALA2/ETHYLENE RESPONSE FACTOR and basic helix-loop-helix tobacco transcription factors cooperatively mediate jasmonate-elicited nicotine biosynthesis. Plant J. 2011;66(6):1053–65.
Baudry A, Heim MA, Dubreucq B, Caboche M, Weisshaar B, Lepiniec L. TT2, TT8, and TTG1 synergistically specify the expression of BANYULS and proanthocyanidin biosynthesis in Arabidopsis thaliana. Plant J. 2004;39(3):366–80.
Zimmermann IM, Heim MA, Weisshaar B, Uhrig JF. Comprehensive identification of Arabidopsis thaliana MYB transcription factors interacting with R/B-like BHLH proteins. Plant J. 2004;40(1):22–34.
Butelli E, Titta L, Giorgio M, Mock HP, Matros A, Peterek S, Schijlen EG, Hall RD, Bovy AG, Luo J. Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors. Nat Biotechnol. 2008;26(11):1301–8.
Grotewold E, Sainz MB, Tagliani L, Hernandez JM, Bowen B, Chandler VL. Identification of the residues in the Myb domain of maize C1 that specify the interaction with the bHLH cofactor R. Proc Natl Acad Sci U S A. 2000;97(25):13579–84.
Feller A, Machemer K, Braun EL, Grotewold E. Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant J. 2011;66(1):94–116.
Goff SA, Cone KC, Chandler VL. Functional analysis of the transcriptional activator encoded by the maize B gene: evidence for a direct functional interaction between two classes of regulatory proteins. Genes Dev. 1992;6(5):864–75.
Arce-Rodriguez ML, Ochoa-Alejo N. An R2R3-MYB transcription factor regulates capsaicinoid biosynthesis. Plant Physiol. 2017;174(3):1359–70.
Sun B, Zhou X, Chen C, Chen C, Chen K, Chen M, Liu S, Chen G, Cao B, Cao F, Lei J. Coexpression network analysis reveals an MYB transcriptional activator involved in capsaicinoid biosynthesis in hot peppers. Hortic Res-England. 2020;7(162):1–14.
Kim S, Park M, Yeom SI, Kim YM, Lee JM, Lee HA, Seo E, Choi J, Cheong K, Kim KT, Jung K, Lee GW, Oh SK, Bae C, Kim SB, Lee HY, Kim SY, Kim MS, Kang BC, Jo YD, Yang HB, Jeong HJ, Kang WH, Kwon JK, Shin C, Lim JY, Park JH, Huh JH, Kim JS, Kim BD, Cohen O, Paran I, Suh MC, Lee SB, Kim YK, Shin Y, Noh SJ, Park J, Seo YS, Kwon SY, Kim HA, Park JM, Kim HJ, Choi SB, Bosland PW, Reeves G, Jo SH, Lee BW, Cho HT, Choi HS, Lee MS, Yu Y, Do Choi Y, Park BS, van Deynze A, Ashrafi H, Hill T, Kim WT, Pai HS, Ahn HK, Yeam I, Giovannoni JJ, Rose JK, Sørensen I, Lee SJ, Kim RW, Choi IY, Choi BS, Lim JS, Lee YH, Choi D. Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat Genet. 2014;46(3):270–8.
Naves ER, de Avila SL, Sulpice R, Araujo WL, Nunes-Nesi A, Peres L, Zsogon A. Capsaicinoids: pungency beyond Capsicum. Trends Plant Sci. 2019;24(2):109–20.
Dong X, Yan Y, Jiang B, Shi Y, Jia Y, Cheng J, Shi Y, Kang J, Li H, Zhang D, Qi L, Han R, Zhang S, Zhou Y, Wang X, Terzaghi W, Gu H, Kang D, Yang S, Li J. The cold response regulator CBF1 promotes Arabidopsis hypocotyl growth at ambient temperatures. EMBO J. 2020;39(13):e103630.
Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouze P, Rombauts S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7.
Hudson KA, Hudson ME. A classification of basic helix-loop-helix transcription factors of soybean. Int J Genomics. 2015;2015(3):603182.
Mao K, Dong Q, Li C, Liu C, Ma F. Genome wide identification and characterization of apple bHLH transcription factors and expression analysis in response to drought and salt stress. Front Plant Sci. 2017;8:480.
Atchley WR, Fitch WM. A natural classification of the basic helix-loop-helix class of transcription factors. Proc Natl Acad Sci U S A. 1997;94(10):5172–6.
Shimizu T, Toumoto A, Ihara K, Shimizu M, Kyogoku Y, Ogawa N, Oshima Y, Hakoshima T. Crystal structure of PHO4 bHLH domain-DNA complex: flanking base recognition. EMBO J. 1997;16(15):4689–97.
Yue H, Wang M, Liu S, Du X, Song W, Nie X. Transcriptome-wide identification and expression profiles of the WRKY transcription factor family in Broomcorn millet (Panicum miliaceum L.). BMC Genomics. 2016;17:343.
Frerigmann H, Berger B, Gigolashvili T. bHLH05 is an interaction partner of MYB51 and a novel regulator of glucosinolate biosynthesis in Arabidopsis. Plant Physiol. 2014;166(1):349–69.
Abe H, Urao T, Ito T, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell. 2003;15(1):63–78.
Lorenzo O, Chico JM, Sanchez-Serrano JJ, Solano R. JASMONATE-INSENSITIVE1 encodes a MYC transcription factor essential to discriminate between different jasmonate-regulated defense responses in Arabidopsis. Plant Cell. 2004;16(7):1938–50.
Yadav V, Mallappa C, Gangappa SN, Bhatia S, Chattopadhyay S. A basic helix-loop-helix transcription factor in Arabidopsis, MYC2, acts as a repressor of blue light-mediated photomorphogenic growth. Plant Cell. 2005;17(7):1953–66.
Nesi N, Debeaujon I, Jond C, Pelletier G, Caboche M, Lepiniec L. The TT8 gene encodes a basic helix-loop-helix domain protein required for expression of DFR and BAN genes in Arabidopsis siliques. Plant Cell. 2000;12(10):1863–78.
Payne CT, Zhang F, Lloyd AM. GL3 encodes a bHLH protein that regulates trichome development in Arabidopsis through interaction with GL1 and TTG1. Genetics. 2000;156(3):1349–62.
Bernhardt C, Lee MM, Gonzalez A, Zhang F, Lloyd A, Schiefelbein J. The bHLH genes GLABRA3 (GL3) and ENHANCER OF GLABRA3 (EGL3) specify epidermal cell fate in the Arabidopsis root. Development. 2003;130(26):6431–9.
Xu J, van Herwijnen ZO, Drager DB, Sui C, Haring MA, Schuurink RC. SlMYC1 regulates type VI glandular trichome formation and terpene biosynthesis in tomato glandular cells. Plant Cell. 2018;30(12):2988–3005.
Feng HL, Ma NN, Meng X, Zhang S, Wang JR, Chai S, Meng QW. A novel tomato MYC-type ICE1-like transcription factor, SlICE1a, confers cold, osmotic and salt tolerance in transgenic tobacco. Plant Physiol Biochem. 2013;73:309–20.
Leivar P, Quail PH. PIFs: pivotal components in a cellular signaling hub. Trends Plant Sci. 2011;16(1):19–28.
Keyhaninejad N, Curry J, Romero J, O’Connell MA. Fruit specific variability in capsaicinoid accumulation and transcription of structural and regulatory genes in Capsicum fruit. Plant Sci. 2014;215–216:59–68.
Koini MA, Alvey L, Allen T, Tilley CA, Harberd NP, Whitelam GC, Franklin KA. High temperature-mediated adaptations in plant architecture require the bHLH transcription factor PIF4. Curr Biol. 2009;19(5):408–13.
Qin C, Yu C, Shen Y, Fang X, Chen L, Min J, Cheng J, Zhao S, Xu M, Luo Y, Yang Y, Wu Z, Mao L, Wu H, Ling-Hu C, Zhou H, Lin H, González-Morales S, Trejo-Saavedra DL, Tian H, Tang X, Zhao M, Huang Z, Zhou A, Yao X, Cui J, Li W, Chen Z, Feng Y, Niu Y, Bi S, Yang X, Li W, Cai H, Luo X, Montes-Hernández S, Leyva-González MA, Xiong Z, He X, Bai L, Tan S, Tang X, Liu D, Liu J, Zhang S, Chen M, Zhang L, Zhang L, Zhang Y, Liao W, Zhang Y, Wang M, Lv X, Wen B, Liu H, Luan H, Zhang Y, Yang S, Wang X, Xu J, Li X, Li S, Wang J, Palloix A, Bosland PW, Li Y, Krogh A, Rivera-Bustamante RF, Herrera-Estrella L, Yin Y, Yu J, Hu K, Zhang Z. Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization. Proc Natl Acad Sci U S A. 2014;111(14):5135–40.
El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, Qureshi M, Richardson LJ, Salazar GA, Smart A, Sonnhammer E, Hirsh L, Paladin L, Piovesan D, Tosatto S, Finn R. The Pfam protein families database in 2019. Nucleic Acids Res. 2019;47(D1):D427–32.
Potter SC, Luciani A, Eddy SR, Park Y, Lopez R, Finn RD. HMMER web server: 2018 update. Nucleic Acids Res. 2018;46(W1):W200–4.
Letunic I, Bork P. 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2018;46(D1):D493–6.
Lu S, Wang J, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Marchler GH, Song JS, Thanki N, Yamashita RA, Yang MZ, Zhang DC, Zheng CJ, Lanczycki CJ, Marchler-Bauer A. CDD/SPARCLE: the conserved domain database in 2020. Nucleic Acids Res. 2020;48(D1):D265–8.
Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003;31(13):3784–8.
Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M, Karthikeyan AS, Lee CH, Nelson WD, Ploetz L, Singh S, Wensel A, Hula E. The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012;40(D1):D1202–10.
Tian F, Yang DC, Meng YQ, Jin J, Gao G. PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res. 2020;48(D1):D1104–13.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.
Subramanian B, Gao S, Lercher MJ, Hu S, Chen WH. Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees. Nucleic Acids Res. 2019;47(W1):W270–5.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server issue):W202–8.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.
Wang Y, Chen G, Lei J, Cao B, Chen C. Identification and characterization of a LEA-like gene, CaMF5, specifically expressed in the anthers of male-fertile Capsicum annuum. Horticultural Plant J. 2020;6(1):39–48.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
The authors thank to lab members for assistance.
The work was supported by grants from the National Natural Science Foundation of China (32070331, 3207180173 and 32072580), China Postdoctoral Science Foundation (2020M682733) and National Key Research and Development Program (2018YFD1000800). The funding organizations played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Carotenoid (A) and capsaicinoid (B) biosynthetic pathways. Figure S2. Chromosomal localizations of CabHLHs. Figure S3. Multiple alignment analysis of the bHLH domain in pepper bHLH proteins. Figure S4. Phylogenetic tree of the bHLH family in Capsicum annuum, Arabidopsis thaliana, tomato and rice. Figure S5. The expression profiles of CabHLH genes in different tissues. Figure S6. Diagram of yeast two-hybrid assays experiments
The corresponding information of 107 CabHLHs. Table S2. Functionally characterized partial bHLH proteins from tomato and Arabidopsis thaliana. Table S3. Features of the CabHLH proteins motifs. Table S4. List of the bHLH genes in rice, tomato and Arabidopsis thaliana. Table S5. List of the primer for real-time quantitative PCR. Table S6. List of the primer for yeast two-hybrid assays. Table S7. Putative cis-elements of capsorubin and capsaicin biosynthetic genes promoters.
About this article
Cite this article
Liu, R., Song, J., Liu, S. et al. Genome-wide identification of the Capsicum bHLH transcription factor family: discovery of a candidate regulator involved in the regulation of species-specific bioactive metabolites. BMC Plant Biol 21, 262 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-021-03004-7
- Yeast two-hybrid assays