Genome-wide identification, expression, and sequence analysis of CONSTANS-like gene family in cannabis reveals a potential role in plant flowering time regulation

Comprehensive analyses of the ARF gene family in cannabis reveals their potential roles in regulating cannabidiol biosynthesis and male flower development

Background

Cannabidiol (CBD), as an important therapeutic property of the cannabis plants, is mainly produced in the flower organs. Auxin response factors (ARFs) are play a crucial role in flower development and secondary metabolite production. However, the specific roles of ARF gene family in cannabis remain unknown.

Methods

In this study, various bioinformatics analysis of CsARF genes were conducted using online website and bioinformatics, quantitative real time PCR technology was used to investigate the expression patterns of the CsARF gene family in different tissues of different cannabis varieties, and subcellular localization analysis was performed in tobacco leaf.

Results

In this study, 22 CsARF genes were identified and found to be unevenly distributed across 9 chromosomes of the cannabis genome. Phylogenetic analysis revealed that the ARF proteins were divided into 4 subgroups. Duplication analysis identified one pair of segmental/whole-genome duplicated CsARF, and three pairs of tandemly duplicated CsARF. Collinearity analysis revealed that two CsARF genes, CsARF4 and CsARF19, were orthologous in both rice and soybean. Furthermore, subcellular localization analysis showed that CsARF2 was localized in the nucleus. Tissue-specific expression analysis revealed that six genes were highly expressed in cannabis male flowers, and among these genes, 3 genes were further found to be highly expressed at different developmental stages of male flowers. Meanwhile, correlation analysis between the expression level of CsARF genes and CBD content in two cultivars 'H8' and 'Y7' showed that the expression level of CsARF13 was negatively correlated with CBD content, while the expression levels of six genes were positively correlated with CBD content. In addition, most of CsARF genes were responsive to IAA treatment.

Conclusion

Our study laid a foundation for the further studies of CsARFs function in cannabis, and provides candidate genes for breeding varieties with high CBD yield in cannabis production.

CONSTANS, a HUB for all seasons: How photoperiod pervades plant physiology regulatory circuits

How does a plant detect the changing seasons and make important developmental decisions accordingly? How do they incorporate daylength information into their routine physiological processes? Photoperiodism, or the capacity to measure the daylength, is a crucial aspect of plant development that helps plants determine the best time of the year to make vital decisions, such as flowering. The protein CONSTANS (CO) constitutes the central regulator of this sensing mechanism, not only activating florigen production in the leaves but also participating in many physiological aspects in which seasonality is important. Recent discoveries place CO in the center of a gene network that can determine the length of the day and confer seasonal input to aspects of plant development and physiology as important as senescence, seed size, or circadian rhythms. In this review, we discuss the importance of CO protein structure, function, and evolutionary mechanisms that embryophytes have developed to incorporate annual information into their physiology.

Genome-wide identification of the Pinus tabuliformis CONSTANS-like gene family and their potential roles in reproductive cone development

The CONSTANS-like (COL) genes, as a core transcription factor in the photoperiod regulation pathway, play a key role in plant reproduction development. However, their molecular characterization has rarely been studied in Pinus tabuliformis. Here, 10 PtCOL genes were identified in the P. tabuliformis genome and multiple sequence alignments have indicated that the PtCOL proteins contained highly conserved B-BOX1 and CCT domains. Sequence similarity analysis showed that PtCOL1 and PtCOL3 had the higher similarity with Norway spruce COLs (PaCOL2 and PaCOL1) and Arabidopsis COLs (AtCOL3, 4 and 5), respectively. Phylogeny and gene structure analyses revealed that PtCOLs were divided into three subgroups, each with identical or similar distributions of exons, introns, and motifs. Moreover, 10 PtCOLs were distributed on 6 chromosomes and PtCOL9 has syntenic gene pairs in both Ginkgo biloba and Sequoiadendron giganteum. Interestingly, in transcriptome profiles, most PtCOLs exhibited a diurnal oscillation pattern under both long (LD) and short (SD) day conditions. Additionally, PtCOLs were highly expressed in needles and female cones, and showed different spatial expression patterns. Among the ten PtCOLs, PtCOL1/3 heterologous overexpression Arabidopsis displayed a delayed-flowering phenotype under SD, indicating that they are likely to play a crucial role in the reproductive development. Additionally, PtCOL1 and PtCOL3 were not only capable of interacting with each other, but they were each capable of interacting with themselves. Furthermore, PtCOL1 and PtCOL3 were also involved in the MADS-box protein-protein interaction (PPI) network in P. tabuliformis cone development. Direct interactions of PtDAL11 with PtCOL1/3 impeded PtCOL1/3 translocation into the nucleus. In summary, this study provided comprehensive understanding for the functions of the PtCOL gene family and revealed their biological roles in the photoperiod-dependent P. tabuliformis cone development.

Comparative genomics of flowering behavior in Cannabis sativa

Cannabis sativa L. is a phenotypically diverse and multi-use plant used in the production of fiber, seed, oils, and a class of specialized metabolites known as phytocannabinoids. The last decade has seen a rapid increase in the licit cultivation and processing of C. sativa for medical end-use. Medical morphotypes produce highly branched compact inflorescences which support a high density of glandular trichomes, specialized epidermal hair-like structures that are the site of phytocannabinoid biosynthesis and accumulation. While there is a focus on the regulation of phytocannabinoid pathways, the genetic determinants that govern flowering time and inflorescence structure in C. sativa are less well-defined but equally important. Understanding the molecular mechanisms that underly flowering behavior is key to maximizing phytocannabinoid production. The genetic basis of flowering regulation in C. sativa has been examined using genome-wide association studies, quantitative trait loci mapping and selection analysis, although the lack of a consistent reference genome has confounded attempts to directly compare candidate loci. Here we review the existing knowledge of flowering time control in C. sativa, and, using a common reference genome, we generate an integrated map. The co-location of known and putative flowering time loci within this resource will be essential to improve the understanding of C. sativa phenology.

Genome-Wide Identification, Characterization and Expression Profiling of the CONSTANS-like Genes in Potato (Solanum tuberosum L.)

CONSTANS-like (COL) genes play important regulatory roles in flowering, tuber formation and the development of the potato (Solanum tuberosum L.). However, the COL gene family in S. tuberosum has not been systematically identified, restricting our knowledge of the function of these genes in S. tuberosum. In our study, we identified 14 COL genes, which were unequally distributed among eight chromosomes. These genes were classified into three groups based on differences in gene structure characteristics. The COL proteins of S. tuberosum and Solanum lycopersicum were closely related and showed high levels of similarity in a phylogenetic tree. Gene and protein structure analysis revealed similarities in the exon-intron structure and length, as well as the motif structure of COL proteins in the same subgroup. We identified 17 orthologous COL gene pairs between S. tuberosum and S. lycopersicum. Selection pressure analysis showed that the evolution rate of COL homologs is controlled by purification selection in Arabidopsis, S. tuberosum and S. lycopersicum. StCOL genes showed different tissue-specific expression patterns. StCOL5 and StCOL8 were highly expressed specifically in the leaves of plantlets. StCOL6, StCOL10 and StCOL14 were highly expressed in flowers. Tissue-specific expression characteristics suggest a functional differentiation of StCOL genes during evolution. Cis-element analysis revealed that the StCOL promoters contain several regulatory elements for hormone, light and stress signals. Our results provide a theoretical basis for the understanding of the in-depth mechanism of COL genes in regulating the flowering time and tuber development in S. tuberosum.

Overexpression of two CONSTANS-like 2 (MiCOL2) genes from mango delays flowering and enhances tolerance to abiotic stress in transgenic Arabidopsis

The CO/COL gene family plays an important role in regulating photoperiod-dependent flowering time in plants. In this study, two COL2 gene homologs, MiCOL2A and MiCOL2B, were isolated from 'SiJiMi' mango, and their expression patterns and functions were characterized. The MiCOL2A and MiCOL2B genes both belonged to the group Ⅰ of CO/COL gene family. MiCOL2A and MiCOL2B exhibited distinct circadian rhythms and were highly expressed in leaves during the flowering induction period. Subcellular localization analysis revealed that MiCOL2A and MiCOL2B are localized in the nucleus. The overexpression of MiCOL2A and MiCOL2B significantly delayed flowering time in Arabidopsis under both long-day (LD) and short-day (SD) conditions. The MiCOL2A and MiCOL2B overexpression Arabidopsis plants exhibited more tolerance to slat and drought stress after abiotic stress treatments, with greater ROS scavenging capacity and protective enzyme activity, less cell damage and death and higher expression of stress response genes than wild type plants. Bimolecular fluorescence complementation (BiFC) analysis showed that MiCOL2A and MiCOL2B interacted with several stress-related proteins, including zinc finger protein 4 (MiZFP4), MYB30-INTERACTING E3 LIGASE 1 (MiMIEL1) and RING zinc finger protein 34 (MiRZFP34). The results indicate that MiCOL2A and MiCOL2B are not only involved in flowering time but also play a positive role in abiotic stress responses in plants.

Genome-Wide Identification and Expression Analysis of Wall-Associated Kinase (WAK) Gene Family in Cannabis sativa L

Wall-associated kinases (WAKs) are receptors that bind pectin or small pectic fragments in the cell wall and play roles in cell elongation and pathogen response. In the Cannabis sativa (Cs) genome, 53 CsWAK/CsWAKL (WAK-like) protein family members were identified and characterized; their amino acid lengths and molecular weights varied from 582 to 983, and from 65.6 to 108.8 kDa, respectively. They were classified into four main groups by a phylogenetic tree. Out of the 53 identified CsWAK/CsWAKL genes, 23 CsWAK/CsWAKL genes were unevenly distributed among six chromosomes. Two pairs of genes on chromosomes 4 and 7 have undergone duplication. The number of introns and exons among CsWAK/CsWAKL genes ranged from 1 to 6 and from 2 to 7, respectively. The promoter regions of 23 CsWAKs/CsWAKLs possessed diverse cis-regulatory elements that are involved in light, development, environmental stress, and hormone responsiveness. The expression profiles indicated that our candidate genes (CsWAK1, CsWAK4, CsWAK7, CsWAKL1, and CsWAKL7) are expressed in leaf tissue. These genes exhibit different expression patterns than their homologs in other plant species. These initial findings are useful resources for further research work on the potential roles of CsWAK/CsWAKL in cellular signalling during development, environmental stress conditions, and hormone treatments.

Genome-Wide Identification, Characterization, and Expression Profile Analysis of CONSTANS-like Genes in Woodland Strawberry (Fragaria vesca)

CONSTANS-like (CO-like) gene is one of the most important regulators in the flowering process of the plant, playing a core role in the photoperiodic flowering induction pathway. In this study, we identified 10 distinct CO-like genes (FveCOs) in woodland strawberry (Fragaria vesca). They were classified into three groups with specific gene structure characteristics or protein domains in each group. The effect of selection pressure on the FveCOs in the woodland strawberry was tested by Ka/Ks, and it was shown that the evolution rate of FveCOs was controlled by purification selection factors. Intraspecific synteny analysis of woodland strawberry FveCOs showed that at least one duplication event existed in the gene family members. Collinearity analysis of woodland strawberry genome with genomes of Arabidopsis, rice (Oryza sativa), and apple (Malus × domestica) showed that CO-like genes of F. vesca and Malus × domestica owned higher similarity for their similar genomes compared with those of other two species. The FveCOs showed different tissue-specific expression patterns. Moreover, real-time quantitative PCR results revealed that the expressions of the most FveCOs followed a 24-h rhythm oscillation under both long-day (LD) and short-day (SD) conditions. Further expression analysis showed that the individual expression changing profile of FveCO3 and FveCO5 was opposite to each other under both LD and SD conditions. Moreover, the expression of FveCO3 and FveCO5 was both negatively correlated with the flowering time variation of the woodland strawberry grown under LD and SD conditions, indicating their potential vital roles in the photoperiodic flowering regulation. Further protein interaction network analysis also showed that most of the candidate interaction proteins of FveCO3 and FveCO5 were predicted to be the flowering regulators. Finally, LUC assay indicated that both FveCO3 and FveCO5 could bind to the promoter of FveFT1, the key regulator of flowering regulation in the woodland strawberry, and thus activate its expression. Taken together, this study laid a foundation for understanding the exact roles of FveCOs in the reproductive development regulation of the woodland strawberry, especially in the photoperiodic flowering process.

Characterization of Phytohormones and Transcriptomic Profiling of the Female and Male Inflorescence Development in Manchurian Walnut (Juglans mandshurica Maxim.)

Flowers are imperative reproductive organs and play a key role in the propagation of offspring, along with the generation of several metabolic products in flowering plants. In Juglans mandshurica, the number and development of flowers directly affect the fruit yield and subsequently its commercial value. However, owing to the lack of genetic information, there are few studies on the reproductive biology of Juglans mandshurica, and the molecular regulatory mechanisms underlying the development of female and male inflorescence remain unclear. In this study, phytohormones and transcriptomic sequencing analyses at the three stages of female and male inflorescence growth were performed to understand the regulatory functions underlying flower development. Gibberellin is the most dominant phytohormone that regulates flower development. In total, 14,579 and 7188 differentially expressed genes were identified after analyzing the development of male and female flowers, respectively, wherein, 3241 were commonly expressed. Enrichment analysis for significantly enriched pathways suggested the roles of MAPK signaling, phytohormone signal transduction, and sugar metabolism. Genes involved in floral organ transition and flowering were obtained and analyzed; these mainly belonged to the M-type MADS-box gene family. Three flowering-related genes (SOC1/AGL20, ANT, and SVP) strongly interacted with transcription factors in the co-expression network. Two key CO genes (CO3 and CO1) were identified in the photoperiod pathway. We also identified two GA20xs genes, one SVP gene, and five AGL genes (AGL8, AGL9, AGL15, AGL19, and AGL42) that contributed to flower development. The findings are expected to provide a genetic basis for the studies on the regulatory networks and reproductive biology in inflorescence development for J. mandshurica.

Genome-Wide Identification and Genetic Variations of the Starch Synthase Gene Family in Rice

Starch is a major ingredient in rice, and the amylose content of starch significantly impacts rice quality. OsSS (starch synthase) is a gene family related to the synthesis of amylose and amylopectin, and 10 members have been reported. In the present study, a synteny analysis of a novel family member belonging to the OsSSIV subfamily that contained a starch synthase catalytic domain showed that three segmental duplications and multiple duplications were identified in rice and other species. Expression data showed that the OsSS gene family is involved in diverse expression patterns. The prediction of miRNA targets suggested that OsSS are possibly widely regulated by miRNA functions, with miR156s targeted to OsSSII-3, especially. Haplotype analysis exhibited the relationship between amylose content and diverse genotypes. These results give new insight and a theoretical basis for the improved amylose content and eating quality of rice.