A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

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.

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.

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).

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.

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.

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.

NLR- and mlo-Based Resistance Mechanisms against Powdery Mildew in Cannabis sativa

Powdery mildew (PM) is one of the most common Cannabis sativa diseases. In spite of this, very few documented studies have characterized the resistance genes involved in PM defense mechanisms, or sources of natural genetic resistance in cannabis. The focus of the present work is on the two primary mechanisms for qualitative resistance against PM. The first is based on resistance (R) genes characterized by conserved nucleotide-binding site and/or leucine-rich repeat domains (NLRs). The second one involves susceptibility (S) genes, and particularly mildew resistance locus o (MLO) genes, whose loss-of-function mutations seem to be a reliable way to protect plants from PM infection. Cannabis defenses against PM are thus discussed, mainly detailing the strategies based on these two mechanisms. Emerging studies about this research topic are also reported and, based on the most significant results, a potential PM resistance model in cannabis plant-pathogen interactions is proposed. Finally, innovative approaches, based on the pyramiding of multiple R genes, as well as on genetic engineering and genome editing methods knocking out S genes, are discussed, to obtain durable PM-resistant cannabis cultivars with a broad-spectrum resistance range.

Transcriptomic and metabolomic analyses reveal the differential accumulation of phenylpropanoids and terpenoids in hemp autotetraploid and its diploid progenitor

BACKGROUND

Cannabis sativa, a dioecious plant that has been cultivated worldwide for thousands of years, is known for its secondary metabolites, especially cannabinoids, which possess several medicinal effects. In this study, we investigated the autopolyploidization effects on the biosynthesis and accumulation of these metabolites, transcriptomic and metabolomic analyses were performed to explore the gene expression and metabolic variations in industrial hemp autotetraploids and their diploid progenitors.

RESULTS

Through these analyses, we obtained 1,663 differentially expressed metabolites and 1,103 differentially expressed genes. Integrative analysis revealed that phenylpropanoid and terpenoid biosynthesis were regulated by polyploidization. No substantial differences were found in the cannabidiol or tetrahydrocannabinol content between tetraploids and diploids. Following polyploidization, some transcription factors, including nine bHLH and eight MYB transcription factors, affected the metabolic biosynthesis as regulators. Additionally, several pivotal catalytic genes, such as flavonol synthase/flavanone 3-hydroxylase, related to the phenylpropanoid metabolic pathway, were identified as being modulated by polyploidization.

CONCLUSIONS

This study enhances the overall understanding of the impact of autopolyploidization in C. sativa and the findings may encourage the application of polyploid breeding for increasing the content of important secondary metabolites in industrial hemp.

A forensic application of genetic markers for distinction between drug-type and fiber-type Cannabis sativa L

Genetic markers can represent a valuable tool for forensic purposes in discriminating between fiber-type and drug-type cannabis. The aim of this research was to evaluate developed genetic markers for tetrahydrocannabinolic acid synthase (THCAS) when applied on certified hemp (14 varieties) and forensic casework samples of four chemotypes (40 seizures). Chemotype-associated PCR-based markers did not enable reliable selective amplification despite the difference in cannabinoid composition. In order to characterize forensic samples of unknown origin, THCAS sequencing was performed. The comparison of THCAS sequences, including additional accessions, indicated high genetic similarity of hemp varieties. Confiscated samples of intermediate, THC, CBD and CBG type were clearly separated from fiber-type accessions and assigned to drug-type cluster. Despite the unknown origin, their position on the tree support the notion that they are more related to drug-type accessions than to the fiber-type. However, no clear distinction between chemotypes was found. Furthermore, 26 amino acid substitutions were revealed in THCAS that clearly separate hemp varieties and neither of them cluster with any other tested sample.

Comprehensive Phenotypic Characterization of Diverse Drug-Type Cannabis Varieties from the Canadian Legal Market

Cannabis (Cannabis sativa L.) stands as a historically significant and culturally important plant, embodying economic, social, and medicinal relevance for human societies. However, years of prohibition and stigmatization have hindered the cannabis research community, which is hugely undersized and suffers from a scarcity of understanding of cannabis genetics and how key traits are expressed or inherited. In this study, we conducted a comprehensive phenotypic characterization of 176 drug-type cannabis accessions, representative of Canada's legal market. We assessed germination methods, evaluated various traits including agronomic, morphological, and cannabinoid profiles, and uncovered significant variation within this population. Notably, the yield displayed a negative correlation with maturity-related traits but a positive correlation with the fresh biomass. Additionally, the potential THC content showed a positive correlation with maturity-related traits but a negative correlation with the yield. Significant differences were observed between the plants derived from regular female seeds and feminized seeds, as well as between the plants derived from cuttings and seeds for different traits. This study advances our understanding of cannabis cultivation, offering insights into germination practices, agronomic traits, morphological characteristics, and biochemical diversity. These findings establish a foundation for precise breeding and cultivar development, enhancing cannabis's potential in the legal market.

Cannabinoids function in defense against chewing herbivores in Cannabis sativa L

In the decades since the first cannabinoids were identified by scientists, research has focused almost exclusively on the function and capacity of cannabinoids as medicines and intoxicants for humans and other vertebrates. Very little is known about the adaptive value of cannabinoid production, though several hypotheses have been proposed including protection from ultraviolet radiation, pathogens, and herbivores. To test the prediction that genotypes with greater concentrations of cannabinoids will have reduced herbivory, a segregating F2 population of Cannabis sativa was leveraged to conduct lab- and field-based bioassays investigating the function of cannabinoids in mediating interactions with chewing herbivores. In the field, foliar cannabinoid concentration was inversely correlated with chewing herbivore damage. On detached leaves, Trichoplusia ni larvae consumed less leaf area and grew less when feeding on leaves with greater concentrations of cannabinoids. Scanning electron and light microscopy were used to characterize variation in glandular trichome morphology. Cannabinoid-free genotypes had trichomes that appeared collapsed. To isolate cannabinoids from confounding factors, artificial insect diet was amended with cannabinoids in a range of physiologically relevant concentrations. Larvae grew less and had lower rates of survival as cannabinoid concentration increased. These results support the hypothesis that cannabinoids function in defense against chewing herbivores.

Cannabis sativa: origin and history, glandular trichome development, and cannabinoid biosynthesis

Is Cannabis a boon or bane? Cannabis sativa has long been a versatile crop for fiber extraction (industrial hemp), traditional Chinese medicine (hemp seeds), and recreational drugs (marijuana). Cannabis faced global prohibition in the twentieth century because of the psychoactive properties of ∆9-tetrahydrocannabinol; however, recently, the perspective has changed with the recognition of additional therapeutic values, particularly the pharmacological potential of cannabidiol. A comprehensive understanding of the underlying mechanism of cannabinoid biosynthesis is necessary to cultivate and promote globally the medicinal application of Cannabis resources. Here, we comprehensively review the historical usage of Cannabis, biosynthesis of trichome-specific cannabinoids, regulatory network of trichome development, and synthetic biology of cannabinoids. This review provides valuable insights into the efficient biosynthesis and green production of cannabinoids, and the development and utilization of novel Cannabis varieties.

Cultivar-dependent phenotypic and chemotypic responses of drug-type Cannabis sativa L. to polyploidization

Cannabis sativa L. is a plant with a wide range of potential medicinal applications. In recent years, polyploidy has gained attention as a potential strategy for rapidly improving C. sativa, which, unlike other modern crops, has not yet benefitted from this established biotechnological application. Currently, no reports on high THCA and CBDA drug-type polyploid cultivars have been published. Moreover, it still needs to be clarified if different cultivars react similarly to polyploidization. For these reasons, we set out to evaluate and compare the phenotype and chemotype of three high Δ9-tetrahydrocannabinolic acid (THCA) and one high cannabidiolic acid (CBDA) drug-type cultivars in their diploid, triploid and tetraploid state through agronomic and metabolomic approaches. Our observations on plant morphology revealed a significant increase in plant height and leaf size with increasing ploidy levels in a cultivar-dependent manner. In contrast, cannabinoids were negatively affected by polyploidization, with the concentration of total cannabinoids, THCA, CBDA and cannabigerolic acid (CBGA) decreasing significantly in higher ploidy levels across all four cultivars. Headspace analysis of volatiles revealed that total volatile content decreased in triploids. On the other hand, tetraploids reacted differently depending on the cultivars. Two THCA dominant cultivars showed an increase in concentrations, while in the other two cultivars, concentrations decreased. Additionally, several rare compounds not present in diploids appeared in higher ploidy levels. Moreover, in one high THCA cultivar, a couple of elite tetraploid genotypes for cannabinoid and volatile production were identified, highlighting the role of cultivar and genotypic variability as an important factor in Cannabis sativa L. polyploids. Overall, our observations on plant morphology align with the giga phenotype observed in polyploids of other plant species. The decrease in cannabinoids and volatiles production in triploids have relevant implications regarding their commercial use. On the other hand, this study found that tetraploidization is a suitable approach to improve Cannabis sativa L. medicinal potential, although the response is cultivar and genotype-dependent. This work lays the ground for further improving, evaluating and harnessing Cannabis sativa L. chemical diversity by the breeding, biotechnological and pharmaceutical sectors.

Characterization of Thermo-Mechanical and Chemical Properties of Polypropylene/Hemp Fiber Biocomposites: Impact of Maleic Anhydride Compatibilizer and Fiber Content

This article presents a comprehensive study on the physical, mechanical, thermal, and chemical properties of polypropylene (PP) composites reinforced with hemp fibers (HF) and compatibilized with maleic anhydride (MAPP). The composites were processed using a twin-screw extruder, followed by hot compression at 190 °C. Subsequently, the composites were analyzed using Izod impact and Shore D hardness tests to evaluate their mechanical properties. Thermal properties were investigated through differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA), while X-ray diffraction (XRD) and Fourier-transform infrared spectroscopy (FTIR) were employed to study their chemical properties. Additionally, a statistical analysis was conducted to compare the average results of the impact and hardness tests. XRD analysis revealed that the addition of HF and MAPP led to the disappearance of peaks corresponding to the beta phase in pure PP. Hemp fibers exhibited an impressive crystallinity of 82.10%, surpassing other natural fibers, and had a significant molecular orientation angle (MFA) of 6.06°, making them highly desirable for engineering applications. The crystallite size was observed to be relatively large, at 32.49 nm. FTIR analysis demonstrated strong interactions between the fiber, compatibilizing agent, and polymer matrix. TGA tests showed that the addition of 5 and 10 wt.% MAPP resulted in complete degradation of the composites, similar to pure PP. DSC analyses indicated a reduction in crystallinity (Xc) due to the incorporation of HF and MAPP. Shore D hardness tests revealed an increase in hardness with the addition of 5 wt.% MAPP, while a steep decline in this property was observed with 10 wt.% MAPP. In terms of impact resistance, fractions of 3 and 5 wt.% MAPP in the composites exhibited improved performance compared to the pure polymer. Analysis of variance (ANOVA) was employed to ensure the statistical reliability of the mechanical test results. This comprehensive study sheds light on the diverse properties of PP composites reinforced with hemp fibers and compatibilized with MAPP, emphasizing their potential as sustainable materials for engineering applications. The results contribute to the understanding of the structural and functional aspects of these composites, guiding future research and developments in the field.

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.

A Guide to Cannabis Virology: From the Virome Investigation to the Development of Viral Biotechnological Tools

Cannabis sativa cultivation is experiencing a period of renewed interest due to the new opportunities for its use in different sectors including food, techno-industrial, construction, pharmaceutical and medical, cosmetics, and textiles. Moreover, its properties as a carbon sequestrator and soil improver make it suitable for sustainable agriculture and climate change mitigation strategies. The increase in cannabis cultivation is generating conditions for the spread of new pathogens. While cannabis fungal and bacterial diseases are better known and characterized, viral infections have historically been less investigated. Many viral infection reports on cannabis have recently been released, highlighting the increasing threat and spread of known and unknown viruses. However, the available information on these pathogens is still incomplete and fragmentary, and it is therefore useful to organize it into a single structured document to provide guidance to growers, breeders, and academic researchers. This review aims to present the historical excursus of cannabis virology, from the pioneering descriptions of virus-like symptoms in the 1940s/50s to the most recent high-throughput sequencing reports. Each of these viruses detected in cannabis will be categorized with an increasing degree of threat according to its potential risk to the crop. Lastly, the development of viral vectors for functional genetics studies will be described, revealing how cannabis virology is evolving not only for the characterization of its virome but also for the development of biotechnological tools for the genetic improvement of this crop.

Challenges and potentials of new breeding techniques in Cannabis sativa

Cannabis sativa L. is an ancient crop used for fiber and seed production and not least for its content of cannabinoids used for medicine and as an intoxicant drug. Due to the psychedelic effect of one of the compounds, tetrahydrocannabinol (THC), many countries had regulations or bands on Cannabis growing, also as fiber or seed crop. Recently, as many of these regulations are getting less tight, the interest for the many uses of this crop is increasing. Cannabis is dioecious and highly heterogenic, making traditional breeding costly and time consuming. Further, it might be difficult to introduce new traits without changing the cannabinoid profile. Genome editing using new breeding techniques might solve these problems. The successful use of genome editing requires sequence information on suitable target genes, a genome editing tool to be introduced into plant tissue and the ability to regenerate plants from transformed cells. This review summarizes the current status of Cannabis breeding, uncovers potentials and challenges of Cannabis in an era of new breeding techniques and finally suggests future focus areas that may help to improve our overall understanding of Cannabis and realize the potentials of the plant.

Improving CBCA synthase activity through rational protein design

Global interest for the minor cannabinoid cannabichromene (CBC) is growing steadily, as potential pharmaceutical applications continue to emerge. Due to low-yielding and unspecific extraction processes from its plant host Cannabis sativa, a biotechnological production is desirable. The complete heterologous biosynthesis of several other cannabinoids has recently been demonstrated as an accessible platform. However, the enzyme involved in the biosynthesis of CBC precursor cannabichromenic acid (CBCA) suffers from comparatively low catalytic efficiency, has not been crystallized, and remains poorly characterized. This study contributes to overcoming these challenges in three unique aspects. A deep‑learning‑assisted prediction of the CBCA synthase crystal structure using DeepMinds AlphaFold is performed and evaluated. The predicted CBCA synthase structure scored considerably higher in various quality assessments than the alternative template‑based homology modeling approach. A robust and practical understanding of crucial structure-function relationships for CBCA synthase is provided and a new binding mode for the substrate uncovered. Rational design approaches and computational analyses to suggest CBCAS variants with facilitated activity are applied. Through subsequent screening the substrate conversion of those variants is compared to the native enzyme. The best variant presented in this study increases CBCA production from crude lysate 22-fold and is one of five positions where substitutions had a significantly favorable impact on product formation.

Emerging mechanistic insights into the regulation of specialized metabolism in plants

Plants biosynthesize a broad range of natural products through specialized and species-specific metabolic pathways that are fuelled by core metabolism, together forming a metabolic network. Specialized metabolites have important roles in development and adaptation to external cues, and they also have invaluable pharmacological properties. A growing body of evidence has highlighted the impact of translational, transcriptional, epigenetic and chromatin-based regulation and evolution of specialized metabolism genes and metabolic networks. Here we review the forefront of this research field and extrapolate to medicinal plants that synthetize rare molecules. We also discuss how this new knowledge could help in improving strategies to produce useful plant-derived pharmaceuticals.

Current status and future prospects in cannabinoid production through in vitro culture and synthetic biology

For centuries, cannabis has been a rich source of fibrous, pharmaceutical, and recreational ingredients. Phytocannabinoids are the most important and well-known class of cannabis-derived secondary metabolites and display a broad range of health-promoting and psychoactive effects. The unique characteristics of phytocannabinoids (e.g., metabolite likeness, multi-target spectrum, and safety profile) have resulted in the development and approval of several cannabis-derived drugs. While most work has focused on the two main cannabinoids produced in the plant, over 150 unique cannabinoids have been identified. To meet the rapidly growing phytocannabinoid demand, particularly many of the minor cannabinoids found in low amounts in planta, biotechnology offers promising alternatives for biosynthesis through in vitro culture and heterologous systems. In recent years, the engineered production of phytocannabinoids has been obtained through synthetic biology both in vitro (cell suspension culture and hairy root culture) and heterologous systems. However, there are still several bottlenecks (e.g., the complexity of the cannabinoid biosynthetic pathway and optimizing the bioprocess), hampering biosynthesis and scaling up the biotechnological process. The current study reviews recent advances related to in vitro culture-mediated cannabinoid production. Additionally, an integrated overview of promising conventional approaches to cannabinoid production is presented. Progress toward cannabinoid production in heterologous systems and possible avenues for avoiding autotoxicity are also reviewed and highlighted. Machine learning is then introduced as a powerful tool to model, and optimize bioprocesses related to cannabinoid production. Finally, regulation and manipulation of the cannabinoid biosynthetic pathway using CRISPR- mediated metabolic engineering is discussed.

Genetic Variants Associated with Long-Terminal Repeats Can Diagnostically Classify Cannabis Varieties

Cannabis sativa (Cannabis) has recently been legalized in multiple countries globally for either its recreational or medicinal use. This, in turn, has led to a marked increase in the number of Cannabis varieties available for use in either market. However, little information currently exists on the genetic distinction between adopted varieties. Such fundamental knowledge is of considerable value and underpins the accelerated development of both a nascent pharmaceutical industry and the commercial recreational market. Therefore, in this study, we sought to assess genetic diversity across 10 Cannabis varieties by undertaking a reduced representation shotgun sequencing approach on 83 individual plants to identify variations which could be used to resolve the genetic structure of the assessed population. Such an approach also allowed for the identification of the genetic features putatively associated with the production of secondary metabolites in Cannabis. Initial analysis identified 3608 variants across the assessed population with phylogenetic analysis of this data subsequently enabling the confident grouping of each variety into distinct subpopulations. Within our dataset, the most diagnostically informative single nucleotide polymorphisms (SNPs) were determined to be associated with the long-terminal repeat (LTRs) class of retroelements, with 172 such SNPs used to fully resolve the genetic structure of the assessed population. These 172 SNPs could be used to design a targeted resequencing panel, which we propose could be used to rapidly screen different Cannabis plants to determine genetic relationships, as well as to provide a more robust, scientific classification of Cannabis varieties as the field moves into the pharmaceutical sphere.

Profiling Cannabinoid Contents and Expression Levels of Corresponding Biosynthetic Genes in Commercial Cannabis (Cannabis sativa L.) Cultivars

Cannabis (Cannabis sativa L.) is widely cultivated and studied for its psychoactive and medicinal properties. As the major cannabinoids are present in acidic forms in Cannabis plants, non-enzymatic processes, such as decarboxylation, are crucial for their conversion to neutral active cannabinoid forms. Herein, we detected the levels of cannabidivarin (CBDV), cannabidiol (CBD), cannabichromene (CBC), and Δ9-tetrahydrocannabinol (Δ9-THC) in the leaves and vegetative shoots of five commercial Cannabis cultivars using a combination of relatively simple extraction, decarboxylation, and high-performance liquid chromatography analyses. The CBDV, CBC, and Δ9-THC levels were 6.3-114.9, 34.4-187.2, and 57.6-407.4 μg/g, respectively, and the CBD levels were the highest, ranging between 1.2-8.9 μg/g in leaf and vegetative shoot tissues of Cannabis cultivars. Additionally, correlations were observed between cannabinoid accumulation and transcription levels of genes encoding key enzymes for cannabinoid biosynthesis, including CsCBGAS, CsCBDAS, CsCBCAS, and CsTHCAS. These data suggest that the high accumulation of cannabinoids, such as CBC, Δ9-THC, and CBD, might be derived from the transcriptional regulation of CsCBGAS and CsCBDAS in Cannabis plants.

Emerging trends in genomic and epigenomic regulation of plant specialised metabolism

Regulation of specialised metabolism genes is multilayered and complex, influenced by an array of genomic, epigenetic and epigenomic mechanisms. Here, we review the most recent knowledge in this field, drawing from discoveries in several plant species. Our aim is to improve understanding of how plant genome structure and function influence specialised metabolism. We also highlight key areas for future exploration. Gene regulatory mechanisms influencing specialised metabolism include gene duplication and neo-functionalization, conservation of operon-like clusters of specialised metabolism genes, local chromatin modifications, and the organisation of higher order chromatin structures within the nucleus. Genomic and epigenomic research to-date in the discipline have focused on a relatively small number of plant species, primarily at whole organ or tissue level. This is largely due to the technical demands of the experimental methods needed. However, a high degree of cell-type specificity of function exists in specialised metabolism, driven by similarly specific gene regulation. In this review we focus on the genomic characteristics of genes that are found in different types of clusters within the genome. We propose that acquisition of cell-resolution epigenomic datasets in emerging models, such as the glandular trichomes of Cannabis sativa, will yield important advances. Data such as chromatin accessibility and histone modification profiles can pinpoint which regulatory sequences are active in individual cell types and at specific times in development. These could provide fundamental biological insight as well as novel targets for genetic engineering and crop improvement.

Origin, Early History, Cultivation, and Characteristics of the Traditional Varieties of Moroccan Cannabis sativa L

Background: Cannabis has been cultivated and used for centuries in the north Moroccan Rif (local name is kif). However, its history is poorly known and the date of its first introduction and dispersal in Morocco is still difficult to be precise. Aim: The purpose of the present work is to review the literature on the origin, history, and cultivation of Cannabis in Morocco, as well as data on the morphological, genetic, and phytochemical characteristics of local cultivated varieties. Discussion: Considering the importance of preserving the fragile environment of the Rif and the future development of the Moroccan medical Cannabis market, which will require authentication of the raw material, the use of local strains which are well adapted to the particular environment of the Rif is highly recommended. However, there is no document that summarizes and clarifies the nomenclature and the characteristics of local Moroccan Cannabis. In addition, the recent adoption by Rif growers of improved hybrid cultivars is obliterating the traits and peculiarities of Moroccan Cannabis through genetic introgression. Conclusion: Summarizing and discussing the data from the literature on the characteristics of local Moroccan Cannabis varieties may be useful for their identification and the localization of the areas of the Rif region where their cultivation is still practiced.

Biosynthetic origins of unusual cannabimimetic phytocannabinoids in Cannabis sativa L: A review

Plants of Cannabis sativa L. (Cannabaceae) produce an array of more than 160 isoprenylated resorcinyl polyketides, commonly referred to as phytocannabinoids. These compounds represent molecules of therapeutic importance due to their modulation of the human endocannabinoid system (ECS). While understanding of the biosynthesis of the major phytocannabinoids Δ⁹-tetrahydrocannabinol (Δ⁹-THC) and cannabidiol (CBD) has grown rapidly in recent years, the biosynthetic origin and genetic regulation of many potentially therapeutically relevant minor phytocannabinoids remain unknown, which limits the development of chemotypically elite varieties of C. sativa. This review provides an up-to-date inventory of unusual phytocannabinoids which exhibit cannabimimetic-like activities and proposes putative metabolic origins. Metabolic branch points exploitable for combinatorial biosynthesis and engineering of phytocannabinoids with augmented therapeutic activities are also described, as is the role of phytocannabinoid remodelling to accelerate the therapeutic portfolio expansion in C. sativa.

Multi-Omics Approaches to Study Molecular Mechanisms in Cannabis sativa

Cannabis (Cannabis sativa L.), also known as hemp, is one of the oldest cultivated crops, grown for both its use in textile and cordage production, and its unique chemical properties. However, due to the legislation regulating cannabis cultivation, it is not a well characterized crop, especially regarding molecular and genetic pathways. Only recently have regulations begun to ease enough to allow more widespread cannabis research, which, coupled with the availability of cannabis genome sequences, is fuelling the interest of the scientific community. In this review, we provide a summary of cannabis molecular resources focusing on the most recent and relevant genomics, transcriptomics and metabolomics approaches and investigations. Multi-omics methods are discussed, with this combined approach being a powerful tool to identify correlations between biological processes and metabolic pathways across diverse omics layers, and to better elucidate the relationships between cannabis sub-species. The correlations between genotypes and phenotypes, as well as novel metabolites with therapeutic potential are also explored in the context of cannabis breeding programs. However, further studies are needed to fully elucidate the complex metabolomic matrix of this crop. For this reason, some key points for future research activities are discussed, relying on multi-omics approaches.

Chemical and genetic variation in feral Cannabis sativa populations across the Nebraska climate gradient

Cannabis sativa is a versatile crop that can be cultivated for fiber, seed, or phytochemicals. To take advantage of this versatility and the potential of Cannabis as a feedstock for the bioeconomy, genomics-enabled breeding programs must be strengthened and expanded. This work contributes to the foundation for such by investigating the phytochemistry and genomics of feral Cannabis populations collected from seventeen counties across the climate gradient of Nebraska. Flower tissue from male and female plants (28 total) was studied using (i) gas chromatography-mass spectrometry to assess cannabinoid profiles and (ii) RNA sequencing to determine transcript abundances. Both male and female flower tissues produced cannabinoids, and, though the compounds were more abundant in female flower tissue, the primary cannabinoid in both was usually cannabidiol. The expression of genes that mediate early steps on the cannabinoid biosynthetic pathway were upregulated in female relative to male flowers, suggesting that female versus male flower tissue cannabinoid abundance may be controlled at least in part at the transcriptional level. DNA sequencing was used to place feral Cannabis plants from Nebraska into a previously described genomic context, revealing that all the plants studied here are much more similar to previously characterized hemp-type Cannabis plants than to drug-type Cannabis plants, at least at the genetic level. This work provides foundational phytochemical knowledge and a large set of high-quality single nucleotide polymorphism markers for future studies of feral Nebraska Cannabis.

In planta Female Flower Agroinfiltration Alters the Cannabinoid Composition in Industrial Hemp (Cannabis sativa L.)

Industrial hemp is a diploid (2n = 20), dioecious plant, and an essential source of various phytochemical productions. More than 540 phytochemicals have been described, some of which proved helpful in the remedial treatment of human diseases. Therefore, further study of hemp phytochemicals in medicine is highly anticipated. Previously, we developed the vacuum agroinfiltration method, which allows the transient gene expression in hemp tissues including female flowers, where cannabinoids are produced and accumulated. In this study, we attempted to alter the composition of total CBD and THC. The RT-PCR and sanger sequence identified eleven copies of the CBDAS gene, two copies of the THCAS gene, and one CBCAS gene. Binary vectors were constructed to overexpress the CBDAS gene and silence the THCAS gene via RNA interference. The Transcript level of the CBDAS gene was increased by more than 10 times than the plants used as a control, which led to a 54% higher total CBD content. The silencing of the THCAS gene led to downregulation of the THCAS gene, with an 80% reduction in transcript levels, and total THC content was reduced to 43% compared with mock plant. These results suggest that hemp vacuum infiltration is highly effective for metabolic engineering of cannabinoids in hemp.

Evaluation of tetrahydrocannabinolic acid (THCA) synthase polymorphisms for distinguishing between marijuana and hemp

The Controlled Substances Act (CSA) classifies marijuana (Cannabis sativa) as a Schedule I illicit drug. However, the recent Agriculture Improvement Act of 2018 (U.S. Farm Bill) removed hemp from the definition of marijuana in the CSA, making it a legal crop. As a result, many hemp products are now available, including strains of hemp buds high in other cannabinoids such as cannabidiol (CBD) or cannabigerol (CBG). The genetic inheritance of chemical phenotype (chemotype) has been widely studied, with the tetrahydrocannabinolic acid (THCA) synthase gene at the forefront. Previous studies have speculated that there are two forms of the THCA gene, one that produces an active enzyme (present in marijuana) and one that cannot produce a functional enzyme (present in hemp). A DNA analysis method is desirable for determining crop type in sample types inconducive to chemical analysis, such as immature crops, trace residues, small leaf fragments, seeds, and root material. This study optimized and evaluated a previously reported single nucleotide polymorphism (SNP) assay for determining C. sativa crop type. Furthermore, the presence or absence of 15 cannabinoids, including THC and THCA, was reported in cannabis reference materials and 15 legal hemp flower samples. The SNP assay correctly identified crop type in most samples. However, several marijuana samples were classified as hemp, and several hemp seeds were classified as marijuana. Two strains of legal CBG hemp flowers were also classified as marijuana, indicating that factors other than the genetic variation of the THCA synthase gene should be considered when determining crop type.

Combined use of specific length amplified fragment sequencing (SLAF-seq) and bulked segregant analysis (BSA) for rapid identification of genes influencing fiber content of hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.), an ancient crop, is a significant source of high-quality fiber that primarily caters to the textile industry worldwide. Fiber content is a crucial quantitative trait for evaluating fiber yield in hemp. Understanding the genetic mechanisms involved in hemp breeding is essential for improving yield. In this study, we developed 660 F1 plants from a cross between Jindao-15 (high fiber content fiber-use variety) and Fire No.1 (low fiber content fiber-use variety), and thirty plants each with high and low fiber content were selected from 305 monoecious plants of this population according to 5%-10% of population size for quantitative traits. The DNA from these plants was extracted to establish two bulk DNA pools and then subjected to the restriction digestion by the enzymes RsaI and HaeIII to obtain 314-364 bp digestion fragments and subjected to sequencing using specific length amplified fragment sequencing (SLAF-seq). Finally, we successfully developed 368,404 SLAF tags, which led to the detection of 25,133 high-quality SNPs. Combing with the resequencing results of parents, the SNPs of mixed pools were then subjected to the SNP-Index correlation algorithm, which revealed four candidate regions related to fiber content traits on Chromosome 1, with a length of 8.68 Mb and containing 389 annotated genes. The annotation information and the comparison results identified 15 genes that were highly likely to modulate the fiber content of hemp. Further, qPCR validation identified six genes (LOC115705530, LOC115705875, LOC115704794, LOC115705371, LOC115705688 and LOC115707511) that were highly positively correlated with influencing the hemp fiber content. These genes were involved in the transcription regulation, auxin and water transportion, one carbon and sugar metabolism. And non-synnoumous mutation SNPs which may play vital role in influencing the fiber content were detected in LOC115705875, LOC115704794, LOC115705688 and LOC115707511. Thus, our study highlights the importance of the combined use of SLAF-Seq and Bulked Segregant analysis (BSA) to locate genes related to hemp fiber content rapidly. Hence, our study provides novel mechanistic inputs for the fast identification of genes related to important agronomic traits of hemp and other crops catering to the textile industry.

Genome-Wide Identification and Expression Analysis of Homeodomain Leucine Zipper Subfamily IV (HD-ZIP IV) Gene Family in Cannabis sativa L

The plant-specific homeodomain zipper family (HD-ZIP) of transcription factors plays central roles in regulating plant development and environmental resistance. HD-ZIP transcription factors IV (HDZ IV) have been involved primarily in the regulation of epidermal structure development, such as stomata and trichomes. In our study, we identified nine HDZ IV-encoding genes in Cannabis sativa L. by conducting a computational analysis of cannabis genome resources. Our analysis suggests that these genes putatively encode proteins that have all the conserved domains of HDZ IV transcription factors. The phylogenetic analysis of HDZ IV gene family members of cannabis, rice (Oryza sativa), and Arabidopsis further implies that they might have followed distinct evolutionary paths after divergence from a common ancestor. All the identified cannabis HDZ IV gene promoter sequences have multiple regulation motifs, such as light- and hormone-responsive elements. Furthermore, experimental evidence shows that different HDZ IV genes have different expression patterns in root, stem, leaf, and flower tissues. Four genes were primarily expressed in flowers, and the expression of CsHDG5 (XP_030501222.1) was also correlated with flower maturity. Fifty-nine genes were predicted as targets of HDZ IV transcription factors. Some of these genes play central roles in pathogen response, flower development, and brassinosteroid signaling. A subcellular localization assay indicated that one gene of this family is localized in the Arabidopsis protoplast nucleus. Taken together, our work lays fundamental groundwork to illuminate the function of cannabis HDZ IV genes and their possible future uses in increasing cannabis trichome morphogenesis and secondary metabolite production.

Standardized genome-wide function prediction enables comparative functional genomics: a new application area for Gene Ontologies in plants

Background

Genome-wide gene function annotations are useful for hypothesis generation and for prioritizing candidate genes potentially responsible for phenotypes of interest. We functionally annotated the genes of 18 crop plant genomes across 14 species using the GOMAP pipeline.

Results

By comparison to existing GO annotation datasets, GOMAP-generated datasets cover more genes, contain more GO terms, and are similar in quality (based on precision and recall metrics using existing gold standards as the basis for comparison). From there, we sought to determine whether the datasets across multiple species could be used together to carry out comparative functional genomics analyses in plants. To test the idea and as a proof of concept, we created dendrograms of functional relatedness based on terms assigned for all 18 genomes. These dendrograms were compared to well-established species-level evolutionary phylogenies to determine whether trees derived were in agreement with known evolutionary relationships, which they largely are. Where discrepancies were observed, we determined branch support based on jackknifing then removed individual annotation sets by genome to identify the annotation sets causing unexpected relationships.

Conclusions

GOMAP-derived functional annotations used together across multiple species generally retain sufficient biological signal to recover known phylogenetic relationships based on genome-wide functional similarities, indicating that comparative functional genomics across species based on GO data holds promise for generating novel hypotheses about comparative gene function and traits.

Tobacco Rattle Virus as a Tool for Rapid Reverse-Genetics Screens and Analysis of Gene Function in Cannabis sativa L

Medical cannabis (Cannabis sativa L.) is quickly becoming a central agricultural crop as its production has continued to increase globally. The recent release of the cannabis reference genomes provides key genetic information for the functional analysis of cannabis genes. Currently, however, the established tools for in vivo gene functional analysis in cannabis are very limited. In this study, we investigated the use of the tobacco rattle virus (TRV) as a possible tool for virus-induced gene silencing (VIGS) and virus-aided gene expression (VAGE). Using leaf photobleaching as a visual marker of PHYTOENE DESATURASE (PDS) silencing, we found that VIGS was largely restricted to the agro-infiltrated leaves. However, when agro-infiltration was performed under vacuum, VIGS increased dramatically, which resulted in intense PDS silencing and an increased photobleaching phenotype. The suitability of TRV as a vector for virus-aided gene expression (VAGE) was demonstrated by an analysis of DsRed fluorescence protein. Interestingly, a DsRed signal was also observed in glandular trichomes in TRV₂-DsRed-infected plants, which suggests the possibility of trichome-related gene function analysis. These results indicate that TRV, despite its limited spread, is an attractive vector for rapid reverse-genetics screens and for the analysis of gene function in cannabis.

Industry-Based Misconceptions Regarding Cross-Pollination of Cannabis spp

Cross-pollination of commercial crops has been an ongoing issue in many species. Cannabis spp. encompasses the classifications of marijuana [high in Δ⁹-tetrahydrocannabinol (THC)] and hemp (below 0.3% THC). As such, cannabis is the most recent crop facing the dilemma of cross-pollination and is leading to litigation. These litigations are driven by the large misunderstanding of the impacts of cross-pollination within the cannabis industry. The misconception is that if hemp is cross-pollinated by high THC cannabis, the hemp will become "hot" (high in THC) thereby rendering the crop illegal under the 2018 Farm Bill. However, there are many factors that contribute to the amount of THC a plant may produce. This article examines and refutes the misconception of cross-pollination increasing THC levels by highlighting several methods of how THC may become high in a given hemp crop.

Cannabis sativa research trends, challenges, and new-age perspectives

Cannabis sativa L. has been one of the oldest medicinal plants cultivated for 10,000 years for several agricultural and industrial applications. However, the plant became controversial owing to some psychoactive components that have adverse effects on human health. In this review, we analyzed the trends in cannabis research for the past two centuries. We discussed the historical transitions of cannabis from the category of herbal medicine to an illicit drug and back to a medicinal product post-legalization. In addition, we address the new-age application of immuno-suppressive and anti-inflammatory extracts for the treatment of COVID-19 inflammation. We further address the influence of the legal aspects of cannabis cultivation for medicinal, pharmaceutical, and biotechnological research. We reviewed the up-to-date cannabis genomic resources and advanced technologies for their potential application in genomic-based cannabis improvement. Overall, this review discusses the diverse aspects of cannabis research developments ranging from traditional use as herbal medicine to the latest potential in COVID-19, legal practices with updated patent status, and current state of art genetic and genomic tools reshaping cannabis biotechnology in modern age agriculture and pharmaceutical industry.

Fertilization Following Pollination Predominantly Decreases Phytocannabinoids Accumulation and Alters the Accumulation of Terpenoids in Cannabis Inflorescences

In the last decades, growing evidence showed the therapeutic capabilities of Cannabis plants. These capabilities were attributed to the specialized secondary metabolites stored in the glandular trichomes of female inflorescences, mainly phytocannabinoids and terpenoids. The accumulation of the metabolites in the flower is versatile and influenced by a largely unknown regulation system, attributed to genetic, developmental and environmental factors. As Cannabis is a dioecious plant, one main factor is fertilization after successful pollination. Fertilized flowers are considerably less potent, likely due to changes in the contents of phytocannabinoids and terpenoids; therefore, this study examined the effect of fertilization on metabolite composition by crossbreeding (-)-Δ⁹-trans-tetrahydrocannabinol (THC)- or cannabidiol (CBD)-rich female plants with different male plants: THC-rich, CBD-rich, or the original female plant induced to develop male pollen sacs. We used advanced analytical methods to assess the phytocannabinoids and terpenoids content, including a newly developed semi-quantitative analysis for terpenoids without analytical standards. We found that fertilization significantly decreased phytocannabinoids content. For terpenoids, the subgroup of monoterpenoids had similar trends to the phytocannabinoids, proposing both are commonly regulated in the plant. The sesquiterpenoids remained unchanged in the THC-rich female and had a trend of decrease in the CBD-rich female. Additionally, specific phytocannabinoids and terpenoids showed an uncommon increase in concentration followed by fertilization with particular male plants. Our results demonstrate that although the profile of phytocannabinoids and their relative ratios were kept, fertilization substantially decreased the concentration of nearly all phytocannabinoids in the plant regardless of the type of fertilizing male. Our findings may point to the functional roles of secondary metabolites in Cannabis.

Genome-Wide Analysis of Alternative Splicing and Non-Coding RNAs Reveal Complicated Transcriptional Regulation in Cannabis sativa L

It is of significance to mine the structural genes related to the biosynthetic pathway of fatty acid (FA) and cellulose as well as explore the regulatory mechanism of alternative splicing (AS), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) in the biosynthesis of cannabinoids, FA and cellulose, which would enhance the knowledge of gene expression and regulation at post-transcriptional level in Cannabis sativa L. In this study, transcriptome, small RNA and degradome libraries of hemp 'Yunma No.1' were established, and comprehensive analysis was performed. As a result, a total of 154, 32 and 331 transcripts encoding key enzymes involved in the biosynthesis of cannabinoids, FA and cellulose were predicted, respectively, among which AS occurred in 368 transcripts. Moreover, 183 conserved miRNAs, 380 C. sativa-specific miRNAs and 7783 lncRNAs were predicted. Among them, 70 miRNAs and 17 lncRNAs potentially targeted 13 and 17 transcripts, respectively, encoding key enzymes or transporters involved in the biosynthesis of cannabinoids, cellulose or FA. Finally, the crosstalk between AS and miRNAs or lncRNAs involved in cannabinoids and cellulose was also predicted. In summary, all these results provided insights into the complicated network of gene expression and regulation in C. sativa.

Cannabis labelling is associated with genetic variation in terpene synthase genes

Analysis of over 100 Cannabis samples quantified for terpene and cannabinoid content and genotyped for over 100,000 single nucleotide polymorphisms indicated that Sativa- and Indica-labelled samples were genetically indistinct on a genome-wide scale. Instead, we found that Cannabis labelling was associated with variation in a small number of terpenes whose concentrations are controlled by genetic variation at tandem arrays of terpene synthase genes.

In silico analysis enabling informed design for genome editing in medicinal cannabis; gene families and variant characterisation

BACKGROUND

Cannabis has been used worldwide for centuries for industrial, recreational and medicinal use, however, to date no successful attempts at editing genes involved in cannabinoid biosynthesis have been reported. This study proposes and develops an in silico best practices approach for the design and implementation of genome editing technologies in cannabis to target all genes involved in cannabinoid biosynthesis.

RESULTS

A large dataset of reference genomes was accessed and mined to determine copy number variation and associated SNP variants for optimum target edit sites for genotype independent editing. Copy number variance and highly polymorphic gene sequences exist in the genome making genome editing using CRISPR, Zinc Fingers and TALENs technically difficult. Evaluation of allele or additional gene copies was determined through nucleotide and amino acid alignments with comparative sequence analysis performed. From determined gene copy number and presence of SNPs, multiple online CRISPR design tools were used to design sgRNA targeting every gene, accompanying allele and homologs throughout all involved pathways to create knockouts for further investigation. Universal sgRNA were designed for highly homologous sequences using MultiTargeter and visualised using Sequencher, creating unique sgRNA avoiding SNP and shared nucleotide locations targeting optimal edit sites.

CONCLUSIONS

Using this framework, the approach has wider applications to all plant species regardless of ploidy number or highly homologous gene sequences.

SIGNIFICANCE STATEMENT

Using this framework, a best-practice approach to genome editing is possible in all plant species, including cannabis, delivering a comprehensive in silico evaluation of the cannabinoid pathway diversity from a large set of whole genome sequences. Identification of SNP variants across all genes could improve genome editing potentially leading to novel applications across multiple disciplines, including agriculture and medicine.

Genome-Wide Characterization of the MLO Gene Family in Cannabis sativa Reveals Two Genes as Strong Candidates for Powdery Mildew Susceptibility

Cannabis sativa is increasingly being grown around the world for medicinal, industrial, and recreational purposes. As in all cultivated plants, cannabis is exposed to a wide range of pathogens, including powdery mildew (PM). This fungal disease stresses cannabis plants and reduces flower bud quality, resulting in significant economic losses for licensed producers. The Mildew Locus O (MLO) gene family encodes plant-specific proteins distributed among conserved clades, of which clades IV and V are known to be involved in susceptibility to PM in monocots and dicots, respectively. In several studies, the inactivation of those genes resulted in durable resistance to the disease. In this study, we identified and characterized the MLO gene family members in five different cannabis genomes. Fifteen Cannabis sativa MLO (CsMLO) genes were manually curated in cannabis, with numbers varying between 14, 17, 19, 18, and 18 for CBDRx, Jamaican Lion female, Jamaican Lion male, Purple Kush, and Finola, respectively (when considering paralogs and incomplete genes). Further analysis of the CsMLO genes and their deduced protein sequences revealed that many characteristics of the gene family, such as the presence of seven transmembrane domains, the MLO functional domain, and particular amino acid positions, were present and well conserved. Phylogenetic analysis of the MLO protein sequences from all five cannabis genomes and other plant species indicated seven distinct clades (I through VII), as reported in other crops. Expression analysis revealed that the CsMLOs from clade V, CsMLO1 and CsMLO4, were significantly upregulated following Golovinomyces ambrosiae infection, providing preliminary evidence that they could be involved in PM susceptibility. Finally, the examination of variation within CsMLO1 and CsMLO4 in 32 cannabis cultivars revealed several amino acid changes, which could affect their function. Altogether, cannabis MLO genes were identified and characterized, among which candidates potentially involved in PM susceptibility were noted. The results of this study will lay the foundation for further investigations, such as the functional characterization of clade V MLOs as well as the potential impact of the amino acid changes reported. Those will be useful for breeding purposes in order to develop resistant cultivars.

Analysis of Sequence Variability and Transcriptional Profile of Cannabinoid synthase Genes in Cannabis sativa L. Chemotypes with a Focus on Cannabichromenic acid synthase

Cannabis sativa L. has been long cultivated for its narcotic potential due to the accumulation of tetrahydrocannabinolic acid (THCA) in female inflorescences, but nowadays its production for fiber, seeds, edible oil and bioactive compounds has spread throughout the world. However, some hemp varieties still accumulate traces of residual THCA close to the 0.20% limit set by European Union, despite the functional gene encoding for THCA synthase (THCAS) is lacking. Even if some hypotheses have been produced, studies are often in disagreement especially on the role of the cannabichromenic acid synthase (CBCAS). In this work a set of European Cannabis genotypes, representative of all chemotypes, were investigated from a chemical and molecular point of view. Highly specific primer pairs were developed to allow an accurate distinction of different cannabinoid synthases genes. In addition to their use as markers to detect the presence of CBCAS at genomic level, they allowed the analysis of transcriptional profiles in hemp or marijuana plants. While the high level of transcription of THCAS and cannabidiolic acid synthase (CBDAS) clearly reflects the chemical phenotype of the plants, the low but stable transcriptional level of CBCAS in all genotypes suggests that these genes are active and might contribute to the final amount of cannabinoids.

Emerging diseases of Cannabis sativa and sustainable management

Cultivation of cannabis plants (Cannabis sativa L., marijuana) has taken place worldwide for centuries. In Canada, legalization of cannabis in October 2018 for the medicinal and recreational markets has spurned interest in large-scale growing. This increased production has seen a rise in the incidence and severity of plant pathogens, causing a range of previously unreported diseases. The objective of this review is to highlight the important diseases currently affecting the cannabis and hemp industries in North America and to discuss various mitigation strategies. Progress in molecular diagnostics for pathogen identification and determining inoculum sources and methods of pathogen spread have provided useful insights. Sustainable disease management approaches include establishing clean planting stock, modifying environmental conditions to reduce pathogen development, implementing sanitation measures, and applying fungal and bacterial biological control agents. Fungicides are not currently registered for use and hence there are no published data on their efficacy. The greatest challenge remains in reducing microbial loads (colony-forming units) on harvested inflorescences (buds). Contaminating microbes may be introduced during the cultivation and postharvest phases, or constitute resident endophytes. Failure to achieve a minimum threshold of microbes deemed to be safe for utilization of cannabis products can arise from conventional and organic cultivation methods, or following applications of beneficial biocontrol agents. The current regulatory process for approval of cannabis products presents a challenge to producers utilizing biological control agents for disease management. © 2021 The Author. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Genome-wide development of insertion-deletion (InDel) markers for Cannabis and its uses in genetic structure analysis of Chinese germplasm and sex-linked marker identification

Background

Cannabis sativa L., a dioecious plant derived from China, demonstrates important medicinal properties and economic value worldwide. Cannabis properties have been usually harnessed depending on the sex of the plant. To analyse the genetic structure of Chinese Cannabis and identify sex-linked makers, genome-wide insertion-deletion (InDel) markers were designed and used.

Results

In this study, a genome-wide analysis of insertion-deletion (InDel) polymorphisms was performed based on the recent genome sequences. In total, 47,558 InDels were detected between the two varieties, and the length of InDels ranged from 4 bp to 87 bp. The most common InDels were tetranucleotides, followed by pentanucleotides. Chromosome 5 exhibited the highest number of InDels among the Cannabis chromosomes, while chromosome 10 exhibited the lowest number. Additionally, 31,802 non-redundant InDel markers were designed, and 84 primers evenly distributed in the Cannabis genome were chosen for polymorphism analysis. A total of 38 primers exhibited polymorphisms among three accessions, and of the polymorphism primers, 14 biallelic primers were further used to analyse the genetic structure. A total of 39 fragments were detected, and the PIC value ranged from 0.1209 to 0.6351. According to the InDel markers and the flowering time, the 115 Chinese germplasms were divided into two subgroups, mainly composed of cultivars obtained from the northernmost and southernmost regions, respectively. Additional two markers, “Cs-I1–10” and “Cs-I1–15”, were found to amplify two bands (398 bp and 251 bp; 293 bp and 141 bp) in the male plants, while 389-bp or 293-bp bands were amplified in female plants. Using the two markers, the feminized and dioecious varieties could also be distinguished.

Conclusion

Based on the findings obtained herein, we believe that this study will facilitate the genetic improvement and germplasm conservation of Cannabis in China, and the sex-linked InDel markers will provide accurate sex identification strategies for Cannabis breeding and production.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07883-w.

Origin and Evolution of the Cannabinoid Oxidocyclase Gene Family

Cannabis is an ancient crop representing a rapidly increasing legal market, especially for medicinal purposes. Medicinal and psychoactive effects of Cannabis rely on specific terpenophenolic ligands named cannabinoids. Recent whole-genome sequencing efforts have uncovered variation in multiple genes encoding the final steps in cannabinoid biosynthesis. However, the origin, evolution, and phylogenetic relationships of these cannabinoid oxidocyclase genes remain unclear. To elucidate these aspects, we performed comparative genomic analyses of Cannabis, related genera within the Cannabaceae family, and selected outgroup species. Results show that cannabinoid oxidocyclase genes originated in the Cannabis lineage from within a larger gene expansion in the Cannabaceae family. Localization and divergence of oxidocyclase genes in the Cannabis genome revealed two main syntenic blocks, each comprising tandemly repeated cannabinoid oxidocyclase genes. By comparing these blocks with those in genomes from closely related species, we propose an evolutionary model for the origin, neofunctionalization, duplication, and diversification of cannabinoid oxidocycloase genes. Based on phylogenetic analyses, we propose a comprehensive classification of three main clades and seven subclades that are intended to aid unequivocal referencing and identification of cannabinoid oxidocyclase genes. Our data suggest that cannabinoid phenotype is primarily determined by the presence/absence of single-copy genes. Although wild populations of Cannabis are still unknown, increased sampling of landraces and wild/feral populations across its native geographic range is likely to uncover additional cannabinoid oxidocyclase sequence variants.

Genomic analyses provide comprehensive insights into the domestication of bast fiber crop ramie (Boehmeria nivea)

Ramie (Boehmeria nivea) is an economically important natural fiber-producing crop that has been cultivated for thousands of years in China; however, the evolution of this crop remains largely unknown. Here, we report a ramie domestication analysis based on genome assembly and resequencing of cultivated and wild accessions. Two chromosome-level genomes representing wild and cultivated ramie were assembled de novo. Numerous structural variations between two assemblies, together with the genetic variations from population resequencing, constituted a comprehensive genomic variation map for ramie. Domestication analysis identified 71 high-confidence selective sweeps comprising 320 predicted genes, and 29 genes from sweeps were associated with fiber growth in the expression. In addition, we identified seven genetic loci associated with the fiber yield trait in the segregated population derived from the crossing of two assembled accessions, and two of which showed an overlap with the selective sweeps. These findings indicated that bast fiber traits were focused on during the domestication history of ramie. This study sheds light on the domestication of ramie and provides a valuable resource for biological and breeding studies of this important crop.

Response of essential oil hemp (Cannabis sativa L.) growth, biomass, and cannabinoid profiles to varying fertigation rates

Five essential oil hemp (Cannabis sativa L.) cultivars (Cherry Blossom, Cherry Blossom (Tuan), Berry Blossom, Cherry Wine, and Cherry Blossom × Trump) were treated with six fertigation treatments to quantify the effects of synthetic fertilizer rates and irrigation electrical conductivity on plant growth, biomass accumulation, and cannabinoid profiles. Irrigation water was injected with a commercial 20-20-20 fertilizer at rates of 0, 50, 150, 300, 450, and 600 ppm nitrogen equating to 0.33 (control), 0.54, 0.96, 1.59, 2.22, and 2.85 dS m-1, respectively. Plants were grown under artificial lighting (18 hr) to maintain vegetative growth for eight weeks, followed by an eight-week flowering period. High linear relationship between chlorophyll concentrations and SPAD-502 measurements validated the utilization of SPAD meters to rapidly identify nutrient deficiency in essential oil hemp. Cultivars expressed significant variation in plant height and cannabinoid profiles (% dry mass), in concurrence with limited biomass and cannabinoid (g per plant) yield variation. Cherry Blossom was the best performing cultivar and Cherry Wine was the least productive. Variation in plant growth, biomass, and cannabinoid concentrations were affected to a greater extent by fertilizer rates. Optimal fertilizer rates were observed at 50 ppm N, while increased fertilizer rates significantly reduced plant growth, biomass accumulation, and cannabinoid concentrations. Increased fertilizer rates (> 300 ppm N) resulted in compliant THC levels (< 0.3%), although when coupled with biomass reductions resulted in minimal cannabinoid yields. Additionally, CBD concentration demonstrated higher sensitivity to increased fertilizer rates (> 300 ppm N) compared to THC and CBG (> 450 ppm N). The results of this study can serve as a guide when using fertigation methods on essential oil hemp cultivars; although results may differ with cultivar selection, environmental conditions, and management practices.