Comparative genomics analysis of bHLH genes in cucurbits identifies a novel gene regulating cucurbitacin biosynthesis

Abstract

The basic helix–loop–helix (bHLH) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, and have undergone significant expansion during land plant evolution by gene duplications. In cucurbit crops, several bHLH genes have been found to be responsible for agronomic traits such as bitterness. However, the characterization of bHLH genes across the genomes of cucurbit species has not been reported, and how they have evolved and diverged remains largely unanswered. Here we identified 1,160 bHLH genes in seven cucurbit crops and performed a comprehensive comparative genomics analysis. We determined orthologous and paralogous bHLH genes across cucurbit crops by syntenic analysis between or within species. Orthology and phylogenetic analysis of the tandem-duplicated bHLH genes in the Bt cluster, which regulate the biosynthesis of cucurbitacins, suggest that this cluster is derived from three ancestral genes after the cucurbit-common tetraploidization event. Interestingly, we identified a new conserved cluster paralogous to the Bt cluster that includes two tandem bHLH genes, and the evolutionary history and expression profiles of these two genes in the new cluster suggest the involvement of one gene (Brp) in the regulation of cucurbitacin biosynthesis in roots. Further biochemical and transgenic assays in melon hairy roots supported the function of Brp. This study provides useful information for further investigating the functions of bHLH TFs and novel insights into the regulation of cucurbitacin biosynthesis in cucurbit crops and other plants.

References

[1] Riechmann JL, Heard J, Martin G et al. Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000; 290: 2105.

[2] Feller A, Machemer K, Braun EL et al. Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant J. 2011; 66: 94-116.

[3] Carretero-Paulet L, Galstyan A, Roig-Villanova I et al. Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in Arabidopsis, poplar, rice, moss, and algae. Plant Physiol. 2010; 153: 1398-412.

[4] Pires N, Dolan L . Origin and diversification of basic-helix-loop-helix proteins in plants. Mol Biol Evol. 2010; 27: 862-74.

[5] Toledo-Ortiz G, Huq E, Quail PH . The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003; 15: 1749-70.

[6] Li X, Duan X, Jiang H et al. Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis. Plant Physiol. 2006; 141: 1167-84.

[7] Wang J, Sun P, Li Y et al. An overlooked Paleotetraploidization in Cucurbitaceae. Mol Biol Evol. 2018; 35: 16-26.

[8] Sun H, Wu S, Zhang G et al. Karyotype stability and unbiased fractionation in the paleo-allotetraploid Cucurbita genomes. Mol Plant. 2017; 10: 1293-306.

[9] Dong H, Zhao H, Li S et al. Genome-wide association studies reveal that members of bHLH subfamily 16 share a conserved function in regulating flag leaf angle in rice (Oryza sativa). PLoS Genet. 2018; 14: e1007323.

[10] Yang F, Wang Q, Schmitz G et al. The bHLH protein ROX acts in concert with RAX1 and LAS to modulate axillary meristem formation in Arabidopsis. Plant J. 2012; 71: 61-70.

[11] Paik I, Kathare PK, Kim J-I et al. Expanding roles of PIFs in signal integration from multiple processes. Mol Plant. 2017; 10: 1035-46.

[12] Radoeva T, Lokerse AS, Llavata-Peris CI et al. A robust auxin response network controls embryo and suspensor development through a basic helix loop helix transcriptional module. Plant Cell. 2019; 31: 52-67.

[13] Ren Y, Zhao Q, Zhao X et al. Expression analysis of the MdCIbHLH1 gene in apple flower buds and seeds in the process of dormancy. Hortic Plant J. 2016; 2: 61-6.

[14] Zhang HB, Bokowiec MT, Rushton PJ et al. Tobacco transcription factors NtMYC2a and NtMYC2b form nuclear complexes with the NtJAZ1 repressor and regulate multiple jasmonate-inducible steps in nicotine biosynthesis. Mol Plant. 2012; 5: 73-84.

[15] Shen Q, Lu X, Yan T et al. The jasmonate-responsive AaMYC2 transcription factor positively regulates artemisinin biosynthesis in Artemisia annua. New Phytol. 2016; 210: 1269-81.

[16] Ribeiro B, Lacchini E, Bicalho KU et al. A seed-specific regulator of triterpene saponin biosynthesis in Medicago truncatula. Plant Cell. 2020; 32: 2020-42.

[17] Shang Y, Ma Y, Zhou Y et al. Biosynthesis, regulation, and domestication of bitterness in cucumber. Science. 2014; 346: 1084-8.

[18] Zhou Y, Ma Y, Zeng J et al. Convergence and divergence of bitterness biosynthesis and regulation in Cucurbitaceae. Nat Plants. 2016; 2: 16183.

[19] Xi H, He Y, Chen H . Functional characterization of SmbHLH13 in anthocyanin biosynthesis and flowering in eggplant. Hortic Plant J. 2021; 7: 73-80.

[20] Massari ME, Murre C . Helix-loop-helix proteins: regulators of transcription in eukaryotic organisms. Mol Cell Biol. 2000; 20: 429-40.

[21] Sun H, Fan HJ, Ling HQ . Genome-wide identification and characterization of the bHLH gene family in tomato. BMC Genomics. 2015; 16: 9.

[22] Heim MA, Jakoby M, Werber M et al. 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: 735-47.

[23] Kane J, Freeling M, Lyons E . The evolution of a high copy gene array in Arabidopsis. J Mol Evol. 2010; 70: 531-44.

[24] Fang L, Cheng F, Wu J et al. The impact of genome triplication on tandem gene evolution in Brassica rapa. Front Plant Sci. 2012; 3: 261.

[25] Wang N, Cui Y, Liu Y et al. Requirement and functional redundancy of Ib subgroup bHLH proteins for iron deficiency responses and uptake in Arabidopsis thaliana. Mol Plant. 2013; 6: 503-13.

[26] Cui Y, Chen C-L, Cui M et al. Four IVa bHLH transcription factors are novel interactors of FIT and mediate JA inhibition of iron uptake in Arabidopsis. Mol Plant. 2018; 11: 1166-83.

[27] Yu J, Ke T, Tehrim S et al. PTGBase: an integrated database to study tandem duplicated genes in plants. Database (Oxford). 2015; 2015: bav017.

[28] Qiao X, Li Q, Yin H et al. Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants. Genome Biol. 2019; 20: 38.

[29] Frey M, Chomet P, Glawischnig E et al. Analysis of a chemical plant defense mechanism in grasses. Science. 1997; 277: 696-9.

[30] Qi X, Bakht S, Leggett M et al. A gene cluster for secondary metabolism in oat: implications for the evolution of metabolic diversity in plants. Proc Natl Acad Sci USA. 2004; 101: 8233-8.

[31] Hen-Avivi S, Savin O, Racovita RC et al. A metabolic gene cluster in the wheat W1 and the barley Cer-cqu loci determines β-diketone biosynthesis and glaucousness. Plant Cell. 2016; 28: 1440-60.

[32] Winzer T, Gazda V, He Z et al. A Papaver somniferum 10-gene cluster for synthesis of the anticancer alkaloid noscapine. Science. 2012; 336: 1704-8.

[33] Shoji T, Yuan L . ERF gene clusters: working together to regulate metabolism. Trends Plant Sci. 2021; 26: 23-32.

[34] Nutzmann HW, Huang A, Osbourn A . Plant metabolic clusters - from genetics to genomics. New Phytol. 2016; 211: 771-89.

[35] Nutzmann HW, Doerr D, Ramírez-Colmenero A et al. Active and repressed biosynthetic gene clusters have spatially distinct chromosome states. Proc Natl Acad Sci USA. 2020; 117: 13800-9.

[36] Van Moerkercke A, Steensma P, Schweizer F et al. The bHLH transcription factor BIS1 controls the iridoid branch of the monoterpenoid indole alkaloid pathway in Catharanthus roseus. Proc Natl Acad Sci USA 2015; 112: 8130-5.

[37] Van Moerkercke A, Steensma P, Gariboldi I et al. The basic helix-loop-helix transcription factor BIS2 is essential for monoterpenoid indole alkaloid production in the medicinal plant Catharanthus roseus. Plant J. 2016; 88: 3-12.

[38] Singh SK, Patra B, Paul P et al. BHLH IRIDOID SYNTHESIS 3 is a member of a bHLH gene cluster regulating terpenoid indole alkaloid biosynthesis in Catharanthus roseus. Plant Direct. 2021; 5: e00305.

[39] Kumar S, Stecher G, Tamura K . MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016; 33: 1870-4.

[40] Kent WJ . BLAT-the BLAST-like alignment tool. Genome Res. 2002; 12: 656-64.

[41] Hu B, Jin J, Guo AY et al. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2015; 31: 1296-7.

[42] Bailey TL, Williams N, Misleh C et al. MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006; 34: W369-73.

[43] Wang Y, Tang H, DeBarrry J et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49- e49.

[44] Chen C, Chen H, Zhang Y et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020; 13: 1194-202.

[45] Cheng F, Wu J, Fang L et al. Syntenic gene analysis between Brassica rapa and other Brassicaceae species. Front Plant Sci. 2012; 3: 198.

[46] Wang D, Zhang Y, Zhang Z et al. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics. 2010; 8: 77-80.

[47] Li Z, Zhang Z, Yan P et al. RNA-Seq improves annotation of protein-coding genes in the cucumber genome. BMC Genomics. 2011; 12: 540.

[48] Yano R, Nonaka S, Ezura H . Melonet-DB, a grand RNA-Seq gene expression atlas in melon (Cucumis melo L.). Plant Cell Physiol. 2018; 59: e4.

[49] Guo S, Zhao S, Sun H et al. Resequencing of 414 cultivated and wild watermelon accessions identifies selection for fruit quality traits. Nat Genet. 2019; 51: 1616-23.

[50] Wu S, Shamimuzzaman M, Sun H et al. The bottle gourd genome provides insights into Cucurbitaceae evolution and facilitates mapping of a Papaya ring-spot virus resistance locus. Plant J. 2017; 92: 963-75.

[51] Jiang B, Xie D, Liu W et al. De novo assembly and characterization of the transcriptome, and development of SSR markers in wax gourd (Benicasa hispida). PLoS One. 2013; 8: e71054.

[52] Cui J, Yang Y, Luo S et al. Whole-genome sequencing provides insights into the genetic diversity and domestication of bitter gourd (Momordica spp.). Hortic Res. 2020; 7: 85.

[53] Pertea M, Kim D, Pertea GM et al. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. 2016; 11: 1650-67.

[54] Wang S, Yang X, Xu M et al. A rare SNP identified a TCP transcription factor essential for tendril development in cucumber. Mol Plant. 2015; 8: 1795-808.