Characterization of new chloroplast markers to determine biogeographical origin and crop type of Cannabis sativa

Characterization and evaluation of nine Cannabis sativa chloroplast SNP markers for crop type determination and biogeographical origin on European samples

Cannabis sativa can be classified in two main types, according to psychotropic cannabinoid ∆9-tetrahydrocannabinol (∆9-THC) content: the drug-type and the fiber-type. According to the European Monitoring Center for Drugs and Drug Addiction, most of the European Union countries consider the possession of cannabis, for personal use, a minor offense with possibility of incarceration. Despite of the model of legal supply (i.e., Spanish cannabis clubs, Netherlands coffee shops) or medical use (i.e., Italy), cannabis remains the most used and trafficked illicit plant in the European Union. Differentiating cannabis crops or tracing the biogeographical origin is crucial for law enforcement purposes. Chloroplast DNA (cpDNA) markers may assist to determine biogeographic origin and to differentiate hemp from marijuana. This research aims: to identify and to evaluate nine C. sativa cpDNA polymorphic SNP sites to differentiate crop type and to provide information about its biogeographical origin. Five SNaPshot™ assays for nine chloroplast markers were developed and conducted in marijuana samples seized in Chile, the USA-Mexico border and Spain, and hemp samples grown in Spain and in Italy. The SNapShot™ assays were tested on 122 cannabis samples, which included 16 blind samples, and were able to differentiate marijuana crop type from hemp crop type in all samples. Using phylogenetic analysis, genetic differences were observed between marijuana and hemp samples. Moreover, principal component analysis (PCA) supported the relationship among hemp samples, as well as for USA-Mexico border, Spanish, and Chilean marijuana samples. Genetic differences between groups based on the biogeographical origin and their crop type were observed. Increasing the number of genetic markers, including the most recently studied ones, and expanding the sample database will provide more accurate information about crop differentiation and biogeographical origin.

Evaluation of chloroplast DNA barcoding markers to individualize Papaver somniferum for forensic intelligence purposes

The opium poppy, Papaver somniferum L., is a forensically important plant due to the medicinal and illegal uses for the milky latex stored in the pods. This latex contains the alkaloids morphine, codeine, and thebaine that are used for their analgesic properties and/or for synthesizing other opioids. However, these compounds are highly addictive and have caused a national opioid epidemic. Two other Papaver species, P. setigerum DC. and P. bracteatum Lindl., are also of forensic interest because they pose both forensic and legal issues. They are largely uncontrolled under the Controlled Substances Act, making these species a common defense strategy. Current morphological and chemical identification methods have been moderately successful but have drawbacks. There is also a lack of sequencing data available. Therefore, exploiting the genome using chloroplast DNA barcoding markers could help to accurately identify these species of interest when plant material is taken. This study screened and assessed the genetic variation both between species and within populations of P. somniferum in nine cpDNA barcode regions (ndhF-rpl32, petA-psbJ, rpl32-trnL, rps16-trnQ, trnE-trnT, trnH-psbA, trnL-trnF, rpl16 intron, and psbE-petL). Published reference genomes from the NCBI GenBank database were aligned and compared for an initial in silico screening. Additionally, ten P. somniferum seed samples from various vendors were sequenced and compared across samples and to published reference data at the various barcode regions of interest. This study showed that the regions trnH-psbA and petA-psbJ have promise for utility in individualization for both inter- and intra-species individualization of P. somniferum.

Cannabis Seeds Authentication by Chloroplast and Nuclear DNA Analysis Coupled with High-Resolution Melting Method for Quality Control Purposes

Background: Cannabis plants and their seed have been used in many cultures as a source of medicine and feeding during history. Today, there is an increasing demand for cannabis seeds for medical use. Moreover, a seed sales market with no legal regulations has also grown. This may pose some issues if a quality control is not set in place. Identification of cannabis strains is important for quality control purposes in a nonregulated growing market and in cases of illegal traffic and medical use. Owing to the high price as a pharmacological drug, commercial products of cannabis plants and seeds for medical users are often subjected to adulterations, either when packing or distributing certified seeds in the market. Materials and Methods: Cannabis commercial seeds and cannabis seeds for medical use were analyzed with high-resolution melting (HRM) analysis using barcoding markers. Humulus lupulus L. plants from a local market were used as outgroup control. DNA barcoding uses specific regions of the genome to identify differences in the genetic sequence of conserved regions such as internal transcribed spacer (ITS) and rbcL. DNA barcoding data can be generated with real-time polymerase chain reaction combined with HRM analysis to distinguish specific conserved DNA regions of closely related species. HRM analysis is the method of choice for rapid analysis of sequence variation. Results: The melting temperature (Tm) of homogeneous packages was consistent with single genotypes. However, packages containing contaminating seeds showed Tm differences of 0.2°C on average. Conclusions: An effective, rapid, and low-cost method based on ITS nuclear DNA and on chloroplast rbcL regions for screening and detection of contamination in commercial cannabis seeds was developed and applied for the analysis of different samples. This approach can be used as a quality control tool for cannabis seeds or other plant material.

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.

Mitochondrial genomes do not appear to regulate flowering pattern/reproductive strategy in Cannabis sativa

Currently, the amount of genetic data for Cannabis is lacking due to the illegal nature of the plant. Our study used 73 Cannabis sativa whole-genome shotgun libraries to reveal eight different mtDNA haplotypes. The most common haplotype contained 60 of the 73 samples studied and was composed of only dioecious individuals. However, other haplotypes contained a mix of both mating strategies (i.e. monoecious and dioecious). From these haplotype groupings we further examined the fully annotated mitochondrial genomes of four hemp individuals with different mt haplotypes and recorded gene content, copy number variation and synteny. Our results revealed highly syntenic mitochondrial genomes that contained ~60 identifiable sequences for protein-coding genes, tRNAs and rRNAs and no obvious rearrangements or chimeric genes. We found no clear evidence that modern reproductive patterns are due to simple cytoplasmic male sterility mutations. It is likely the interaction between nuclear genetic components and the X/Y sex chromosomes that determines reproductive strategy. Additionally, we added 50 % more mitochondrial genomes to the publicly available repository.

Massively parallel sequencing of Cannabis sativa chloroplast hotspots for forensic typing

BACKGROUND

Marijuana (Cannabis sativa) is the most commonly used illicit drug in the USA, and the use of DNA barcodes could assist drug trafficking investigations by indicating the biogeographical origin and crop type of a sample and providing a means for linking cases. Additionally, the legality of marijuana in the USA remains complicated with some states fully legalizing marijuana for recreational use while federally marijuana remains completely illegal. Massively parallel sequencing (MPS) offers distinct advantages over capillary electrophoresis (CE), including more comprehensive coverage of target loci, analysis of hundreds of markers simultaneously, and high throughput capabilities.

METHODS

This study reports on the development of a MiSeq FGx® assay targeting seven "hotspot" regions in the Cannabis sativa chloroplast genome that are highly polymorphic and informative in attempts to determine biogeographical origin and distinguishing between marijuana and hemp. Sequencing results were compared to previous studies that used CE-based genotyping methods.

RESULTS

A total of 49 polymorphisms were observed, 16 of which have not been previously reported. Additionally, sequence data revealed isoalleles at one locus, which were able to differentiate two samples that had the same haplotype using CE-based methods. This study reports preliminary results from sequencing 14 hemp and marijuana samples from different countries using the developed MPS assay.

CONCLUSION

Future studies should genotype a more comprehensive sample set from around the world to build a haplotype database, which could be used to provide investigative leads for law enforcement agencies investigating marijuana trafficking.

Development and Validation of a Novel and Fast Detection Method for Cannabis sativa: A 19-Plex Short Tandem Repeat Typing System

In recent years, influenced by the legalization of Cannabis sativa in some countries and regions, the number of people who smoke or abuse C. sativa has continuously grown, cases of transnational C. sativa trafficking have also been increasing. Therefore, fast and accurate identification and source tracking of C. sativa have become urgent social needs. In this study, we developed a new 19-plex short tandem repeats (STRs) typing system for C. sativa, which includes 15 autosomal STRs (D02-CANN1, C11-CANN1, 4910, B01-CANN1, E07-CANN1, 9269, B05-CANN1, H06-CANN2, 5159, nH09, CS1, ANUCS 305, 3735, and ANUCS 302 and 9043), two X-chromosome STRs (ANUCS 501 and 1528), one sex-determining marker (DM016, on Y-chromosome), and a quality control marker (DM029, on autosome). The whole polymerase chain reaction (PCR) process could finish within 1 h, making the system suitable for fast detection. The PCR products were detected and separated with an Applied Biosystems 3500XL Genetic Analyser. Developmental validation studies indicated that the 19-plex typing system was accurate, reliable and sensitive, which could also deconvolute mixed C. sativa samples. Specifically, the sensitivity study showed that a full genotyping profile was obtainable with as low as 125 pg of C. sativa DNA. The species specificity study demonstrated that this multiplex has no cross-reactivity with common non C. sativa DNA. In the population study, a total of 162 alleles at 15 autosomal STRs and 14 alleles at two X-chromosome STRs were detected among 85 samples. The efficiency parameters, including the total discrimination power (TDP) and the combined power of exclusion (CPE) of the system, were calculated to exceed 0.999 999 999 999 988 and 0.998 455 889 684 078, respectively, further proving that the system could meet the needs of individual identification. To the extent of the known studies, this is the first study that included the C. sativa sex-determining marker. In conclusion, the developed new 19-plex STR typing system can successfully achieve the purposes of species identification, gender determination, and individual identification, which could be a powerful tool in tracing trade routes of particular drug syndicates or dealers or in linking certain C. sativa to a crime scene.

A Ge.F.I. - ISFG European collaborative study on DNA identification of Cannabis sativa samples using a 13-locus multiplex STR method

Cannabis sativa is the most used controlled substance in Europe. With the advent of new and less restrictive European laws on cannabis sale for recreational use (including in Italy), an increase in indoor cannabis crops were observed. This increase was possible due to the availability of cannabis seeds through the internet market. Genetic identification of cannabis can link seizures and if in possession then might aid in an investigation. A 13-locus multiplex STR method was previously developed and validated by Houston et al. A collaborative exercise was organized by the Italian Forensic Geneticists - International Society of Forensic Genetics (Ge.F.I. - ISFG) Working Group with the aim to test the reproducibility, reliability and robustness of this multiplex cannabis STR kit. Twenty-one laboratories from three European countries participated in the collaborative exercise and were asked to perform STR typing of two cannabis samples. Cannabis DNA samples and the multiplex STR kit were provided by the University of Barcelona and Sam Houston State University. Different platforms for PCR amplification, capillary electrophoresis (CE) and genotyping software were selected at the discretion of the participating laboratories. Although the participating laboratories used different PCR equipment, CE platforms and genotyping software, concordant results were obtained from the majority of the samples. The overall genotyping success ratio was 96%. Only minor artifacts were observed. The mean peak height ratio was estimated to be 76.3% and 78.1% for sample 1 and sample 2, respectively. The lowest amount of -1 / + 1 stutter percentage produced, when the height of the parent allele was higher than 8000 RFU, resulted to be less than 10% of the parent allele height. Few common issues were observed such as a minor peak imbalance in some heterozygous loci, some artifact peaks and few instances of allelic drop-out. The results of this collaborative exercise demonstrated the robustness and applicability of the 13-locus system for cannabis DNA profiling for forensic purposes.

Evaluation of the trnK-matK-trnK, ycf3, and accD-psal chloroplast regions to differentiate crop type and biogeographical origin of Cannabis sativa

Cannabis sativa (marijuana and hemp) is one of the most controversial crops worldwide. In the USA, the state-specific legalization of marijuana and recently legalized hemp pose a problem for law enforcement. This study seeks to utilize chloroplast hSTRs, INDEL, and SNPs markers to develop genotyping methods to aid in the differentiation of legal hemp from illicit marijuana and also for tracking the flow of trafficked marijuana. Three polymorphic regions: trnK-matK-trnK, ycf3, and accD-psal, of the C. sativa chloroplast genome were evaluated in order to distinguish crop type and biogeographic origin. A total of nine polymorphic sites were genotyped from five distinct populations (hemp from the USA and Canada, marijuana from Chile and USA-Mexico, and medical marijuana from Chile) with a custom fragment and SNaPshot^TM assay. The study also combined genotype results from the same sample set using 21 additional polymorphic markers from previous studies. The effectiveness of these multi-locus assays to distinguish sample groups was assessed using haplotype analysis, phylogenetic analysis, pairwise comparisons, and principal component analysis. Results indicated a clear separation of Canadian hemp using only the nine polymorphic sites developed in this study. The additional 21 markers were able to separate US hemp from both marijuana groups to a significant level (p < 0.05) when assessing average Fixation Indices (F_ST). This study demonstrated the applicability of these organelle markers for the determination of crop type and biogeographic origin of C. sativa. However, a more extensive database is needed to evaluate the true discriminatory power of these markers.

Investigation of chloroplast regions rps16 and clpP for determination of Cannabis sativa crop type and biogeographical origin

Cannabis sativa can be classified as either hemp (a legal crop containing less than 0.3% delta-9-tetrahydrocannabinol, THC) or marijuana (an illegal drug containing more than 0.3% THC). Despite its legalization in 33 states for medicinal or recreational use, marijuana remains the most commonly used illicit drug in the USA, and it is heavily trafficked into and within the country. Discriminating between marijuana and hemp is critical to the legal process. Genetic analysis provides a means of analyzing samples unsuitable for chemical analysis, and in addition to discriminating between crop types, DNA may be able to determine the biogeographical origin of samples. In addition, the sharing of rare haplotypes between different seizures may be useful for linking cases and providing investigative leads to law enforcement. This study evaluates the potential of two highly polymorphic regions of the chloroplast genome of C. sativa, rps16 and clpP, to be used for determination of crop type and biogeographical origin. Custom fragment analysis and SNaPshot™ assays were developed to genotype nine polymorphic loci in hemp samples from the USA and Canada, marijuana samples from USA-Mexico and Chile, and medical marijuana samples from Chile. Haplotype analysis revealed eight haplotypes. Only Canadian hemp could be completely differentiated from the other sample groups by haplotype. Phylogenetic analysis and principal component analysis suggested a closer relationship among USA-Mexico marijuana, Chilean marijuana and medical marijuana, and USA hemp. Genotyping additional polymorphisms in future studies is expected to reveal further differences between these sample groups.