In silico identification of myb and bhlh families reveals candidate transcription factors for secondary metabolic pathways in cannabis sativa L

Plant secondary metabolic pathways are finely regulated by the activity of transcription factors, among which members of the bHLH and MYB subfamilies play a main role. Cannabis sativa L. is a unique officinal plant species with over 600 synthesized phytochemicals having diverse scale-up industrial and pharmaceutical usage. Despite comprehensive knowledge of cannabinoids’ metabolic pathways, very little is known about their regulation, while the literature on flavonoids’ metabolic pathways is still scarce. In this study, we provide the first genome-wide analysis of bHLH and MYB families in C. sativa reference cultivar CBDRx and identification of candidate coding sequences for these transcription factors. Cannabis sativa bHLHs and MYBs were then classified into functional subfamilies through comparative phylogenetic analysis with A. thaliana transcription factors. Analyses of gene structure and motif distribution confirmed that CsbHLHs and CsMYBs belonging to the same evolutionary clade share common features at both gene and amino acidic level. Candidate regulatory genes for key metabolic pathways leading to flavonoid and cannabinoid synthesis in Cannabis were also retrieved. Furthermore, a candidate gene approach was used to identify structural enzyme-coding genes for flavonoid and cannabinoid synthesis. Taken as a whole, this work represents a valuable resource of candidate genes for further investigation of the C. sativa cannabinoid and flavonoid metabolic pathways for genomic studies and breeding programs.