Identification of Chemotypic Markers in Three Chemotype Categories of Cannabis Using Secondary Metabolites Profiled in Inflorescences, Leaves, Stem Bark, and Roots.

Previous chemotaxonomic studies of cannabis only focused on tetrahydrocannabinol (THC) dominant strains while excluded the cannabidiol (CBD) dominant strains and intermediate strains (THC ≈ CBD). This study investigated the utility of the full spectrum of secondary metabolites in different plant parts in three cannabis chemotypes (THC dominant, intermediate, and CBD dominant) for chemotaxonomic discrimination. Hierarchical clustering, principal component analysis (PCA), and canonical correlation analysis assigned 21 cannabis varieties into three chemotypes using the content and ratio of cannabinoids, terpenoids, flavonoids, sterols, and triterpenoids across inflorescences, leaves, stem bark, and roots. The same clustering results were obtained using secondary metabolites, omitting THC and CBD. Significant chemical differences were identified in these three chemotypes. Cannabinoids, terpenoids, flavonoids had differentiation power while sterols and triterpenoids had none. CBD dominant strains had higher amounts of total CBD, cannabidivarin (CBDV), cannabichromene (CBC), α-pinene, ß-myrcene, (-)-guaiol, ß-eudesmol, α-eudesmol, α-bisabolol, orientin, vitexin, and isovitexin, while THC dominant strains had higher total THC, total tetrahydrocannabivarin (THCV), total cannabigerol (CBG), camphene, limonene, ocimene, sabinene hydrate, terpinolene, linalool, fenchol, α-terpineol, ß-caryophyllene, trans-ß-farnesene, α-humulene, trans-nerolidol, quercetin, and kaempferol. Compound levels in intermediate strains were generally equal to or in between those in CBD dominant and THC dominant strains. Overall, with higher amounts of ß-myrcene, (-)-guaiol, ß-eudesmol, α-eudesmol, and α-bisabolol, intermediate strains more resemble CBD dominant strains than THC dominant strains. The results of this study provide a comprehensive profile of bioactive compounds in three chemotypes for medical purposes. The simultaneous presence of a predominant number of identified chemotype markers (with or without THC and CBD) could be used as chemical fingerprints for quality standardization or strain identification for research, clinical studies, and cannabis product manufacturing.

Classification of cannabis strains in the Canadian market with discriminant analysis of principal components using genome-wide single nucleotide polymorphisms.

The cannabis community typically uses the terms "Sativa" and "Indica" to characterize drug strains with high tetrahydrocannabinol (THC) levels. Due to large scale, extensive, and unrecorded hybridization in the past 40 years, this vernacular naming convention has become unreliable and inadequate for identifying or selecting strains for clinical research and medicinal production. Additionally, cannabidiol (CBD) dominant strains and balanced strains (or intermediate strains, which have intermediate levels of THC and CBD), are not included in the current classification studies despite the increasing research interest in the therapeutic potential of CBD. This paper is the first in a series of studies proposing that a new classification system be established based on genome-wide variation and supplemented by data on secondary metabolites and morphological characteristics. This study performed a whole-genome sequencing of 23 cannabis strains marketed in Canada, aligned sequences to a reference genome, and, after filtering for minor allele frequency of 10%, identified 137,858 single nucleotide polymorphisms (SNPs). Discriminant analysis of principal components (DAPC) was applied to these SNPs and further identified 344 structural SNPs, which classified individual strains into five chemotype-aligned groups: one CBD dominant, one balanced, and three THC dominant clusters. These structural SNPs were all multiallelic and were predominantly tri-allelic (339/344). The largest portion of these SNPs (37%) occurred on the same chromosome containing genes for CBD acid synthases (CBDAS) and THC acid synthases (THCAS). The remainder (63%) were located on the other nine chromosomes. These results showed that the genetic differences between modern cannabis strains were at a whole-genome level and not limited to THC or CBD production. These SNPs contained enough genetic variation for classifying individual strains into corresponding chemotypes. In an effort to elucidate the confused genetic backgrounds of commercially available cannabis strains, this classification attempt investigated the utility of DAPC for classifying modern cannabis strains and for identifying structural SNPs.

A single nucleotide polymorphism assay sheds light on the extent and distribution of genetic diversity, population structure and functional basis of key traits in cultivated north American cannabis.

BACKGROUND: The taxonomic classification of Cannabis genus has been delineated through three main types: sativa (tall and less branched plant with long and narrow leaves), indica (short and highly branched plant with broader leaves) and ruderalis (heirloom type with short stature, less branching and small thick leaves). While still under discussion, particularly whether the genus is polytypic or monotypic, this broad classification reflects putative geographical origins of each group and putative chemotype and pharmacologic effect. METHODS: Here we describe a thorough investigation of cannabis accessions using a set of 23 highly informative and polymorphic SNP (Single Nucleotide Polymorphism) markers associated with important traits such as cannabinoid and terpenoid expression as well as fibre and resin production. The assay offers insight into cannabis population structure, phylogenetic relationship, population genetics and correlation to secondary metabolite concentrations. We demonstrate the utility of the assay for rapid, repeatable and cost-efficient genotyping of commercial and industrial cannabis accessions for use in product traceability, breeding programs, regulatory compliance and consumer education. RESULTS: We identified 5 clusters in the sample set, including industrial hemp (K5) and resin hemp, which likely underwent a bottleneck to stabilize cannabidiolic acid (CBDA) accumulation (K2, Type II & III). Tetrahydrocannabinolic acid (THCA) resin (Type I) makes up the other three clusters with terpinolene (K4 - colloquial "sativa" or "Narrow Leaflet Drug" (NLD), myrcene/pinene (K1) and myrcene/limonene/linalool (K3 - colloquial "indica", "Broad Leaflet Drug" (BLD), which also putatively harbour an active version of the cannabichrometic acid Synthase gene (CBCAS). CONCLUSION: The final chemical compositions of cannabis products have key traits related to their genetic identities. Our analyses in the context of the NCBI Cannabis sativa Annotation Release 100 allows for hypothesis testing with regards to secondary metabolite production. Genetic markers related to secondary metabolite production will be important in many sectors of the cannabis marketplace. For example, markers related to THC production will be important for adaptable and compliant large-scale seed production under the new US Domestic Hemp Production Program.

Novel genomic resources for a climate change sensitive mammal: characterization of the American pika transcriptome.

BACKGROUND: When faced with climate change, species must either shift their home range or adapt in situ in order to maintain optimal physiological balance with their environment. The American pika (Ochotona princeps) is a small alpine mammal with limited dispersal capacity and low tolerance for thermal stress. As a result, pikas have become an important system for examining biotic responses to changing climatic conditions. Previous research using amplified fragment length polymorphisms (AFLPs) has revealed evidence for environmental-mediated selection in O. princeps populations distributed along elevation gradients, yet the anonymity of AFLP loci and lack of available genomic resources precluded the identification of associated gene regions. Here, we harnessed next-generation sequencing technology in order to characterize the American pika transcriptome and identify a large suite of single nucleotide polymorphisms (SNPs), which can be used to elucidate elevation- and site-specific patterns of sequence variation. RESULTS: We constructed pooled cDNA libraries of O. princeps from high (1400 m) and low (300 m) elevation sites along a previously established transect in British Columbia. Transcriptome sequencing using the Roche 454 GS FLX titanium platform generated 780 million base pairs of data, which were assembled into 7,325 high coverage contigs. These contigs were used to identify 24,261 novel SNP loci. Using high resolution melt analysis, we developed 17 of these SNPs into genotyping assays, which were validated with independent DNA samples from British Columbia Canada and Oregon State USA. In addition, we detected haplotypes in the NADH dehydrogenase subunit 5 of the mitochondrial genome that were fixed and different among elevations, suggesting that this may be an informative target gene for studying the role of cellular respiration in local adaptation. We also identified contigs that were unique to each elevation, including a high elevation-specific contig that was a positive match with the hemoglobin alpha chain from the plateau pika, a species restricted to high elevation steppes in Asia. Elevation-specific contigs may represent candidate regions subject to differential levels of gene expression along this elevation gradient. CONCLUSIONS: To our knowledge, this is the first broad-scale, transcriptome-level study conducted within the Ochotonidae, providing novel genomic resources for studying pika ecology, behaviour and population history.

Genetic evidence for restricted dispersal along continuous altitudinal gradients in a climate change-sensitive mammal: the American Pika.

When faced with rapidly changing environments, wildlife species are left to adapt, disperse or disappear. Consequently, there is value in investigating the connectivity of populations of species inhabiting different environments in order to evaluate dispersal as a potential strategy for persistence in the face of climate change. Here, we begin to investigate the processes that shape genetic variation within American pika populations from the northern periphery of their range, the central Coast Mountains of British Columbia, Canada. At these latitudes, pikas inhabit sharp elevation gradients ranging from sea level to 1500 m, providing an excellent system for studying the effects of local environmental conditions on pika population genetic structure and gene flow. We found low levels of neutral genetic variation compared to previous studies from more southerly latitudes, consistent with the relatively recent post-glacial colonization of the study location. Moreover, significant levels of inbreeding and marked genetic structure were detected within and among sites. Although low levels of recent gene flow were revealed among elevations within a transect, potentially admixed individuals and first generation migrants were identified using discriminant analysis of principal components between populations separated by less than five kilometers at the same elevations. There was no evidence for historical population decline, yet there was signal for recent demographic contractions, possibly resulting from environmental stochasticity. Correlative analyses revealed an association between patterns of genetic variation and annual heat-to-moisture ratio, mean annual precipitation, precipitation as snow and mean maximum summer temperature. Changes in climatic regimes forecasted for the region may thus potentially increase the rate of population extirpation by further reducing dispersal between sites. Consequently, American pika may have to rely on local adaptations or phenotypic plasticity in order to survive predicted climate changes, although additional studies are required to investigate the evolutionary potential of this climate change sensitive species.

A noninvasive hair sampling technique to obtain high quality DNA from elusive small mammals.

Noninvasive genetic sampling approaches are becoming increasingly important to study wildlife populations. A number of studies have reported using noninvasive sampling techniques to investigate population genetics and demography of wild populations. This approach has proven to be especially useful when dealing with rare or elusive species. While a number of these methods have been developed to sample hair, feces and other biological material from carnivores and medium-sized mammals, they have largely remained untested in elusive small mammals. In this video, we present a novel, inexpensive and noninvasive hair snare targeted at an elusive small mammal, the American pika (Ochotona princeps). We describe the general set-up of the hair snare, which consists of strips of packing tape arranged in a web-like fashion and placed along travelling routes in the pikas' habitat. We illustrate the efficiency of the snare at collecting a large quantity of hair that can then be collected and brought back to the lab. We then demonstrate the use of the DNA IQ system (Promega) to isolate DNA and showcase the utility of this method to amplify commonly used molecular markers including nuclear microsatellites, amplified fragment length polymorphisms (AFLPs), mitochondrial sequences (800bp) as well as a molecular sexing marker. Overall, we demonstrate the utility of this novel noninvasive hair snare as a sampling technique for wildlife population biologists. We anticipate that this approach will be applicable to a variety of small mammals, opening up areas of investigation within natural populations, while minimizing impact to study organisms.