A versatile, rapid Agrobacterium-mediated transient expression system for functional genomics studies in cannabis seedling.

MAIN CONCLUSION: We have developed and optimized a rapid, versatile Agrobacterium-mediated transient expression system for cannabis seedlings that can be used in functional genomics studies of both hemp-type and drug-type cannabis. Cannabis (Cannabis sativa L.) holds great promise in the medical and food industries due to its diverse chemical composition, including specialized cannabinoids. However, the study of key genes involved in various biological processes, including secondary metabolite biosynthesis, has been hampered by the lack of efficient in vivo functional analysis methods. Here, we present a novel, short-cycle, high-efficiency transformation method for cannabis seedlings using Agrobacterium tumefaciens. We used the RUBY reporter system to monitor transformation results without the need for chemical treatments or specialized equipment. Four strains of A. tumefaciens (GV3101, EHA105, LBA4404, and AGL1) were evaluated for transformation efficiency, with LBA4404 and AGL1 showing superior performance. The versatility of the system was further demonstrated by successful transformation with GFP and GUS reporter genes. In addition, syringe infiltration was explored as an alternative to vacuum infiltration, offering simplicity and efficiency for high-throughput applications. Our method allows rapid and efficient in vivo transformation of cannabis seedlings, facilitating large-scale protein expression and high-throughput characterization studies.

Intra-leaf modeling of Cannabis leaflet shape produces leaf models that predict genetic and developmental identities.

The iconic, palmately compound leaves of Cannabis have attracted significant attention in the past. However, investigations into the genetic basis of leaf shape or its connections to phytochemical composition have yielded inconclusive results. This is partly due to prominent changes in leaflet number within a single plant during development, which has so far prevented the proper use of common morphometric techniques. Here, we present a new method that overcomes the challenge of nonhomologous landmarks in palmate, pinnate, and lobed leaves, using Cannabis as an example. We model corresponding pseudo-landmarks for each leaflet as angle-radius coordinates and model them as a function of leaflet to create continuous polynomial models, bypassing the problems associated with variable number of leaflets between leaves. We analyze 341 leaves from 24 individuals from nine Cannabis accessions. Using 3591 pseudo-landmarks in modeled leaves, we accurately predict accession identity, leaflet number, and relative node number. Intra-leaf modeling offers a rapid, cost-effective means of identifying Cannabis accessions, making it a valuable tool for future taxonomic studies, cultivar recognition, and possibly chemical content analysis and sex identification, in addition to permitting the morphometric analysis of leaves in any species with variable numbers of leaflets or lobes.

Comparative proteomic and metabolomic analyses reveal stress responses of hemp to salinity.

KEY MESSAGE: Integrated omics analyses outline the cellular and metabolic events of hemp plants in response to salt stress and highlight several photosynthesis and energy metabolism related pathways as key regulatory points. Soil salinity affects many physiological processes of plants and leads to crop yield losses worldwide. For hemp, a crop that is valued for multiple aspects, such as its medical compounds, fibre, and seed, a comprehensive understanding of its salt stress responses is a prerequisite for resistance breeding and tailoring its agronomic performance to suit certain industrial applications. Here, we first observed the phenotype of salt-stressed hemp plants and found that under NaCl treatment, hemp plants displayed pronounced growth defects, as indicated by the significantly reduced average height, number of leaves, and chlorophyll content. Next, we conducted comparative proteomics and metabolomics to dissect the complex salt-stress response mechanisms. A total of 314 proteins and 649 metabolites were identified to be differentially behaving upon NaCl treatment. Functional classification and enrichment analysis unravelled that many differential proteins were proteases associated with photosynthesis. Through metabolic pathway enrichment, several energy-related pathways were found to be altered, such as the biosynthesis and degradation of branched-chain amino acids, and our network analysis showed that many ribosomal proteins were involved in these metabolic adaptations. Taken together, for hemp plants, influences on chloroplast function probably represent a major toxic effect of salinity, and modulating several energy-producing pathways possibly through translational regulation is presumably a key protective mechanism against the negative impacts. Our data and analyses provide insights into our understanding of hemp's stress biology and may lay a foundation for future functional genomics studies.

Integrated Transcriptomic and Metabolomic Analysis Reveal the Dynamic Process of Bama Hemp Seed Development and the Accumulation Mechanism of α-Linolenic Acid and Linoleic Acid.

Bama County is a world-famous longevity county in the Guangxi Province, China. Bama hemp is a traditional seed used in hemp cultivation in the Bama County. The seeds contain abundant unsaturated fatty acids, particularly linoleic acid (LA) and linolenic acid in the golden ratio. These two substances have been proven to be related to human health and the prevention of various diseases. However, the seed development and seed oil accumulation mechanisms remain unclear. This study employed a combined analysis of physiological, transcriptomic, and metabolomic parameters to elucidate the fatty acid formation patterns in Bama hemp seeds throughout development. We found that seed oil accumulated at a late stage in embryo development, with seed oil accumulation following an "S″-shaped growth curve, and positively correlated with seed size, sugar content, protein content, and starch content. Transcriptome analysis identified genes related to the metabolism of LA, α-linolenic acid (ALA), and jasmonic acid (JA). We found that the FAD2 gene was upregulated 165.26 folds and the FAD3 gene was downregulated 6.15 folds at day 21. Metabolomic changes in LA, ALA, and JA compounds suggested a competitive relationship among these substances. Our findings indicate that the peak period of substance accumulation and nutrient accumulation in Bama hemp seeds occurs during the midstage of seed development (day 21) rather than in the late stage (day 40). The results of this research will provide a theoretical basis for local cultivation and deep processing of Bama hemp.

Quantitative Estimation of Promoter Activity in Cannabis sativa Using Agroinfiltration-Based Transient Gene Expression.

To properly assess promoter activity, which is critical for understanding biosynthetic pathways in different plant species, we use agroinfiltration-based transient gene expression assay. We compare the activity of several known promoters in Nicotiana benthamiana with their activity in Cannabis sativa (both hemp and medicinal cannabis), which has attracted much attention in recent years for its industrial, medicinal, and recreational properties. Here we describe an optimized protocol for transient expression in Cannabis combined with a ratiometric GUS reporter system that allows more accurate evaluation of promoter activity and reduces the effects of variable infiltration efficiency.

Bacterial endophytic community composition varies by hemp cultivar in commercially sourced seed.

The seed-endophytic bacterial community is a potentially beneficial and heritable fraction of the plant microbiome. Its utilization as a sustainable crop improvement strategy could be especially valuable for species such as hemp, where production is being scaled up and new challenges will be faced in managing crop productivity and health. However, little is known about the makeup and variation of the hemp seed microbiome. This study profiled the endophytic bacterial communities harboured by 16 hemp cultivars sourced from commercial suppliers in Europe. A 16S rDNA amplicon sequencing approach identified 917 amplicon sequence variants across samples. Taxonomic classification of sequences revealed 4 phyla and 87 genera to be represented in the dataset. Several genera were widespread while some were specific to one or a few cultivars. Flavobacterium, Pseudomonas, and Pantoea were notable in their high overall abundance and prevalence, but community composition was variable and no one taxon was universally abundant, suggesting a high degree of flexibility in community assembly. Taxonomic composition and alpha diversity differed among cultivars, though further work is required to understand the relative influence of hemp genetic factors on community structure. The taxonomic profiles presented here can be used to inform further work investigating the functional characteristics and potential plant-growth-promoting traits of seed-borne bacteria in hemp.

Genetic insights into agronomic and morphological traits of drug-type cannabis revealed by genome-wide association studies.

Cannabis sativa L., previously concealed by prohibition, is now a versatile and promising plant, thanks to recent legalization, opening doors for medical research and industry growth. However, years of prohibition have left the Cannabis research community lagging behind in understanding Cannabis genetics and trait inheritance compared to other major crops. To address this gap, we conducted a comprehensive genome-wide association study (GWAS) of nine key agronomic and morphological traits, using a panel of 176 drug-type Cannabis accessions from the Canadian legal market. Utilizing high-density genotyping-by-sequencing (HD-GBS), we successfully generated dense genotyping data in Cannabis, resulting in a catalog of 800 K genetic variants, of which 282 K common variants were retained for GWAS analysis. Through GWAS analysis, we identified 18 markers significantly associated with agronomic and morphological traits. Several identified markers exert a substantial phenotypic impact, guided us to putative candidate genes that reside in high linkage-disequilibrium (LD) with the markers. These findings lay a solid foundation for an innovative cannabis research, leveraging genetic markers to inform breeding programs aimed at meeting diverse needs in the industry.

Low UV radiation influenced DNA methylation, gene regulation, cell proliferation, viability, and biochemical differentiation in the cell suspension cultures of Cannabis indica.

The effect of low artificial Ultraviolet (UV) on the DNA methylation remains controversial. This study addresses how differential photoperiods of UV radiation affect the biochemical and molecular behaviors of Cannabis indica cell suspension cultures. The cell suspensions were illuminated with the compact fluorescent lamps (CFL), emitting a combination of 10% UVB, 30% UVA, and the rest visible wavelengths for 0, 4, 8, and 16 h. The applied photoperiods influenced cell morphological characteristics. The 4 h photoperiod was the most effective treatment for improving biomass, growth index and cell viability percentage while these indices remained non-significant in the 16 h treatment. The methylation-sensitive amplified polymorphism (MASP) assay revealed that the UV radiation was epigenetically accompanied by DNA hypermethylation. The light-treated cells significantly displayed higher relative expression of the cannabidiolic| acid synthase (CBDAS) and delta9-tetrahydrocannabinolic acid synthase (THCAS) genes about 4-fold. The expression of the olivetolic acid cyclase (OAC) and olivetol synthase (OLS) genes exhibited an upward trend in response to the UV radiation. The light treatments also enhanced the proline content and protein concentration. The 4 h illumination was significantly capable of improving the cannabidiol (CBD) and delta-9-tetrahydrocannabinol (THC) concentrations, in contrast with 16 h. By increasing the illumination exposure time, the activity of the phenylalanine ammonia-lyase (PAL) enzyme linearly upregulated. The highest amounts of the phenylpropanoid derivatives were observed in the cells cultured under the radiation for 4 h. Taken collective, artificial UV radiation can induce DNA methylation modifications and impact biochemical and molecular differentiation in the cell suspensions in a photoperiod-dependent manner.

Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L.

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

Genome-wide identification and analysis of ascorbate peroxidase (APX) gene family in hemp (Cannabis sativa L.) under various abiotic stresses.

Ascorbate peroxidase (APX) plays a critical role in molecular mechanisms such as plant development and defense against abiotic stresses. As an important economic crop, hemp (Cannabis sativa L.) is vulnerable to adverse environmental conditions, such as drought, cold, salt, and oxidative stress, which lead to a decline in yield and quality. Although APX genes have been characterized in a variety of plants, members of the APX gene family in hemp have not been completely identified. In this study, we (1) identified eight members of the CsAPX gene family in hemp and mapped their locations on the chromosomes using bioinformatics analysis; (2) examined the physicochemical characteristics of the proteins encoded by these CsAPX gene family members; (3) investigated their intraspecific collinearity, gene structure, conserved domains, conserved motifs, and cis-acting elements; (4) constructed a phylogenetic tree and analyzed interspecific collinearity; and (5) ascertained expression differences in leaf tissue subjected to cold, drought, salt, and oxidative stresses using quantitative real-time-PCR (qRT-PCR). Under all four stresses, CsAPX6, CsAPX7, and CsAPX8 consistently exhibited significant upregulation, whereas CsAPX2 displayed notably higher expression levels under drought stress than under the other stresses. Taken together, the results of this study provide basic genomic information on the expression of the APX gene family and pave the way for studying the role of APX genes in abiotic stress.

Loss of daylength sensitivity by splice site mutation in Cannabis pseudo-response regulator.

Photoperiod insensitivity (auto-flowering) in drug-type Cannabis sativa circumvents the need for short day (SD) flowering requirements making outdoor cultivation in high latitudes possible. However, the benefits of photoperiod insensitivity are counterbalanced by low cannabinoid content and poor flower quality in auto-flowering genotypes. Despite recent studies in cannabis flowering, a mechanistic understanding of photoperiod insensitivity is still lacking. We used a combination of genome-wide association study and genetic fine-mapping to identify the genetic cause of auto-flowering in cannabis. We then used gene expression analyses and transient transformation assays to characterize flowering time control. Herein, we identify a splice site mutation within circadian clock gene PSEUDO-RESPONSE REGULATOR 37 (CsPRR37) in auto-flowering cannabis. We show that CsPRR37 represses FT expression and its circadian oscillations transition to a less repressive state during SD as compared to long days (LD). We identify several key circadian clock genes whose expression is altered in auto-flowering cannabis, particularly under non-inductive LD. Research into the pervasiveness of this mutation and others affecting flowering time will help elucidate cannabis domestication history and advance cannabis breeding toward a more sustainable outdoor cultivation system.

Deep learning-based quantification and transcriptomic profiling reveal a methyl jasmonate-mediated glandular trichome formation pathway in Cannabis sativa.

Cannabis glandular trichomes (GTs) are economically and biotechnologically important structures that have a remarkable morphology and capacity to produce, store, and secrete diverse classes of secondary metabolites. However, our understanding of the developmental changes and the underlying molecular processes involved in cannabis GT development is limited. In this study, we developed Cannabis Glandular Trichome Detection Model (CGTDM), a deep learning-based model capable of differentiating and quantifying three types of cannabis GTs with a high degree of efficiency and accuracy. By profiling at eight different time points, we captured dynamic changes in gene expression, phenotypes, and metabolic processes associated with GT development. By integrating weighted gene co-expression network analysis with CGTDM measurements, we established correlations between phenotypic variations in GT traits and the global transcriptome profiles across the developmental gradient. Notably, we identified a module containing methyl jasmonate (MeJA)-responsive genes that significantly correlated with stalked GT density and cannabinoid content during development, suggesting the existence of a MeJA-mediated GT formation pathway. Our findings were further supported by the successful promotion of GT development in cannabis through exogenous MeJA treatment. Importantly, we have identified CsMYC4 as a key transcription factor that positively regulates GT formation via MeJA signaling in cannabis. These findings provide novel tools for GT detection and counting, as well as valuable information for understanding the molecular regulatory mechanism of GT formation, which has the potential to facilitate the molecular breeding, targeted engineering, informed harvest timing, and manipulation of cannabinoid production.

The development of a next-generation sequencing panel targeting cannabinoid synthase genes to distinguish between marijuana and hemp.

Hemp and marijuana, both derived from Cannabis sativa L. (C. sativa), are subject to divergent legal regulations due to their different Δ9-tetrahydrocannabinol (Δ9-THC) contents. Cannabinoid synthase genes are considered the key enzymes that determine the chemical composition or chemotype of a particular cultivar. However, existing methods for crop type differentiation based on previous synthase gene theories have limitations in terms of precision and specificity, and a wider range of cannabis varieties must be considered when examining cannabis-based genetic markers. A custom next-generation sequencing (NGS) panel was developed targeting all synthase genes, including Δ9-THC acid synthase, cannabidiolic acid synthase, and cannabichromenic acid synthase, as well as the pseudogenes across diverse C. sativa samples, spanning reference hemp and marijuana, commercial hemp derivatives, and seized marijuana extracts. Interpretation of NGS data revealed a relationship between genotypes and underlying chemotypes, with the principal component analysis indicating a clear distinction between hemp and marijuana clusters. This differentiation was attributed to variations in both synthase genes and pseudogene variants. Finally, this study proposes a genetic cannabis classification method using a differentiation flow chart with novel synthase markers. The flow chart successfully differentiated hemp from marijuana with a 1.3% error rate (n = 147).

Characterization of indigenous populations of cannabis in Iran: a morphological and phenological study.

BACKGROUND: Cannabis is a historically, culturally, and economically significant crop in human societies, owing to its versatile applications in both industry and medicine. Over many years, native cannabis populations have acclimated to the various environments found throughout Iran, resulting in rich genetic and phenotypic diversity. Examining phenotypic diversity within and between indigenous populations is crucial for effective plant breeding programs. This study aimed to classify indigenous cannabis populations in Iran to meet the needs of breeders and breeding programs in developing new cultivars. RESULTS: Here, we assessed phenotypic diversity in 25 indigenous populations based on 12 phenological and 14 morphological traits in male and female plants. The extent of heritability for each parameter was estimated in both genders, and relationships between quantitative and time-based traits were explored. Principal component analysis (PCA) identified traits influencing population distinctions. Overall, populations were broadly classified into early, medium, and late flowering groups. The highest extent of heritability of phenological traits was found in Start Flower Formation Time in Individuals (SFFI) for females (0.91) Flowering Time 50% in Individuals (50% of bracts formed) (FT50I) for males (0.98). Populations IR7385 and IR2845 exhibited the highest commercial index (60%). Among male plants, the highest extent of Relative Growth Rate (RGR) was observed in the IR2845 population (0.122 g.g- 1.day- 1). Finally, populations were clustered into seven groups according to the morphological traits in female and male plants. CONCLUSIONS: Overall, significant phenotypic diversity was observed among indigenous populations, emphasizing the potential for various applications. Early-flowering populations, with their high RGR and Harvest Index (HI), were found as promising options for inclusion in breeding programs. The findings provide valuable insights into harnessing the genetic diversity of indigenous cannabis for diverse purposes.

Leaf Spot Disease Development and Its Effect on Yield of Essential Oil-Producing Hemp Cultivars in Kentucky.

Leaf spot diseases are common on field hemp, causing defoliation, and blighting. The most prevalent leaf spot diseases of hemp in Kentucky are Bipolaris leaf spot (Bipolaris gigantea), Cercospora leaf spot (Cercospora flagellaris), and Septoria leaf spot (Septoria cannabis). In this study, disease progression, cultivar susceptibility, and yield loss were examined using cultivars from four relatedness groups at two locations and in two growing seasons. Septoria leaf spot was the first leaf spot disease to be observed in the field, followed by Bipolaris leaf spot. Both diseases reduced canopy density. Cercospora leaf spot developed in the late reproductive stages as harvest approached. A wide range of susceptibility was documented, suggesting genetic variability across cultivars. Trump group cultivars were the most susceptible, while Otto II group cultivars were the least susceptible. Most importantly, leaf spot diseases had minimal influence on floral biomass and no effect on CBD yield, suggesting that, regardless of disease severity, leaf spot diseases may seldom warrant management. While the importance of foliar disease and corresponding yield loss can shift over time, variation in disease progress among leaf spot diseases and susceptibility of hemp cultivars documented in this study suggest potential disease management through cultural practices such as cultivar and planting date selections.

Development of a real-time PCR assay for detection of hemp russet mite (Aculops cannabicola).

Of the many arthropod species affecting hemp (Cannabis sativa L.) cultivation in the United States, one species of particular importance is the hemp russet mite (Aculops cannabicola, HRM). Hemp russet mite is a microscopic arthropod which feeds on all parts of hemp plants. Due to its minute size, HRM can proliferate undetected for a long time, complicating management efforts and causing serious economic losses. DNA sequencing and PCR assays can facilitate accurate identification and early detection of HRM in infested-plants. Therefore, a real-time SYBR Green based species-specific PCR assay (quantitative PCR, qPCR) was developed for the identification of HRM DNA by amplification of a 104 bp Internal Transcribed Spacer 1 (ITS1) sequence. The detection limit was estimated to be approximately 48 copies of the HRM marker gene sequence. The real-time-PCR assay is rapid, detects all life stages of mite under 2 hours. A 10-fold serial dilution of the plasmid DNA containing the ITS1 insert were used as standards in the real-time PCR assay. The quantification cycle (Cq) value of the assay showed a strong linear relationship with HRM DNA with R2 of 0.96. The assay was tested against several commonly found hemp pests including two-spotted spider mite and western flower thrips to determine specificity of the assay and to show that no non-target species DNA was amplified. The outcomes of this research will have important applications for agricultural biosecurity through accurate identification of HRM, early detection and timely deployment of management tactics to manage and prevent pest outbreaks.

Problematic Cannabis Use Is Associated with Reduced Rectal Microbial Species Richness and Diversity Among a Pilot Sample of Young Sexual and Gender Minorities.

Compared to young heterosexual men, young sexual and gender minorities (YSGM) have elevated systemic inflammation and unique intestinal microbial profiles, influenced by HIV infection and substance use. However, links between cannabis use and microbial dysbiosis in this population have not been well described. In this pilot study, we aimed to characterize the complex interrelationships between cannabis use and microbial community structure in YSGM in relationship to HIV status. Cannabis use was assessed by self-administered Cannabis Use Disorder Identification Test (CUDIT) questionnaires and rectal microbial community alpha-diversity metrics were assessed via 16S ribosomal ribonucleic acid (rRNA) sequencing in a subset of YSGM (n = 42) in the RADAR cohort (aged 16-29) in Chicago. Multivariable regression models were used to assess the relationship between cannabis use and microbiome alpha-diversity metrics, adjusting for HIV status and other risk characteristics, including inflammation, which was evaluated by plasma levels of C-reactive protein (CRP). Problematic cannabis use, but not general use, was significantly inversely associated with microbial community richness (Adj. Beta = -8.13; 95% confidence interval [CI]: -15.68 to -0.59) and Shannon diversity (Adj. Beta = -0.04; 95% CI: -0.07 to 0.009). No significant association was observed between CUDIT score and community evenness, nor was any significant moderation observed by HIV status. We observed that problematic cannabis use was associated with reduced microbial community richness and Shannon diversity, adjusting for within population differences in inflammation and HIV status. Future research should aim to assess how cannabis use contributes to microbiome-related health factors among YSGM and if decreasing cannabis use can restore gut microbial community structure.

Genetic Mapping, Identification, and Characterization of a Candidate Susceptibility Gene for Powdery Mildew in Cannabis sativa.

Powdery mildew (PM) in Cannabis sativa is most frequently caused by the biotrophic fungus Golovinomyces ambrosiae. Based on previously characterized variation in susceptibility to PM, biparental populations were developed by crossing the most resistant cultivar evaluated, 'FL 58', with a susceptible cultivar, 'TJ's CBD'. F1 progeny were evaluated and displayed a range of susceptibility, and two were self-pollinated to generate two F2 populations. In 2021, the F2 populations (n = 706) were inoculated with PM and surveyed for disease severity. In both F2 populations, 25% of the progeny were resistant, while the remaining 75% showed a range of susceptibility. The F2 populations, as well as selected F1 progeny and the parents, were genotyped with a single-nucleotide polymorphism array, and a consensus genetic map was produced. A major effect quantitative trait locus on C. sativa chromosome 1 (Chr01) and other smaller-effect quantitative trait loci (QTL) on four other chromosomes were identified. The most associated marker on Chr01 was located near CsMLO1, a candidate susceptibility gene. Genomic DNA and cDNA sequencing of CsMLO1 revealed a 6.8-kb insertion in FL 58, relative to TJ's CBD, of which 846 bp are typically spliced into the mRNA transcript encoding a premature stop codon. Molecular marker assays were developed using CsMLO1 sequences to distinguish PM-resistant and PM-susceptible genotypes. These data support the hypothesis that a mutated MLO susceptibility gene confers resistance to PM in C. sativa and provides new genetic resources to develop resistant cultivars. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.