Journal:Combined ambient ionization mass spectrometric and chemometric approach for the differentiation of hemp and marijuana varieties of Cannabis sativa

From LIMSWiki
Jump to navigationJump to search
Full article title Combined ambient ionization mass spectrometric and chemometric approach for the differentiation of hemp and marijuana varieties of Cannabis sativa
Journal Journal of Cannabis Research
Author(s) Chambers, Megan I.; Beyramysoltan, Samira; Garosi, Benedetta; Musah, Rabi A.
Author affiliation(s) State University of New York
Primary contact Email: rmusah at albany dot edu
Year published 2023
Volume and issue 5
Article # 5
DOI 10.1186/s42238-023-00173-0
ISSN 2522-5782
Distribution license Creative Commons Attribution 4.0 International
Website https://jcannabisresearch.biomedcentral.com/articles/10.1186/s42238-023-00173-0
Download https://jcannabisresearch.biomedcentral.com/counter/pdf/10.1186/s42238-023-00173-0.pdf (PDF)

Abstract

Background: Hemp and marijuana are the two major varieties of Cannabis sativa. While both contain Δ9-tetrahydrocannabinol (THC), the primary psychoactive component of C. sativa, they differ in the amount of THC that they contain. Presently, U.S. federal laws stipulate that C. sativa containing greater than 0.3% THC is classified as marijuana, while plant material that contains less than or equal to 0.3% THC is hemp. Current methods to determine THC content are chromatography-based, which requires extensive sample preparation to render the materials into extracts suitable for sample injection, for complete separation and differentiation of THC from all other analytes present. This can create problems for forensic laboratories due to the increased workload associated with the need to analyze and quantify THC in all C. sativa materials.

Method: The work presented herein combines direct analysis in real time high-resolution mass spectrometry (DART-HRMS) and advanced chemometrics to differentiate hemp and marijuana plant materials. Samples were obtained from several sources (e.g., commercial vendors, DEA-registered suppliers, and the recreational Cannabis market). DART-HRMS enabled the interrogation of plant materials with no sample pretreatment. Advanced multivariate data analysis approaches, including random forest and principal component analysis (PCA), were used to optimally differentiate these two varieties with a high level of accuracy.

Results: When PCA was applied to the hemp and marijuana data, distinct clustering that enabled their differentiation was observed. Furthermore, within the marijuana class, subclusters between recreational and DEA-supplied marijuana samples were observed. A separate investigation using the silhouette width index to determine the optimal number of clusters for the marijuana and hemp data revealed this number to be two. Internal validation of the model using random forest demonstrated an accuracy of 98%, while external validation samples were classified with 100% accuracy.

Discussion: The results show that the developed approach would significantly aid in the analysis and differentiation of C. sativa plant materials prior to launching painstaking confirmatory testing using chromatography. However, to maintain and/or enhance the accuracy of the prediction model and keep it from becoming outdated, it will be necessary to continue to expand it to include mass spectral data representative of emerging hemp and marijuana strains/cultivars.

Keywords: Cannabis sativa, ambient ionization mass spectrometry, direct analysis in real time—high-resolution mass spectrometry, multivariate data analysis, random forest, principal component analysis

Background

Among the greatest challenges to emerge for U.S. forensic laboratories in recent years are those attributed to the increased legalization and decriminalization of marijuana at the state level, in addition to the permitted production of hemp. The 2019 National Institute of Justice (NIJ) Report to Congress: Needs Assessment of Forensic Laboratories and Medical Examiner/Coroner Offices identified this area as requiring focused attention towards improving criminal justice practices in the USA.[1] The challenge that hemp and marijuana present is as follows: both are major varieties of the same species Cannabis sativa, often referred to as Cannabis. While they each contain Δ9-tetrahydrocannabinol (THC), which is the primary psychoactive component of C. sativa, marijuana and hemp differ in the amount of this molecule that is present. In 2018, the U.S. federal guidelines stipulated that C. sativa which contains greater than 0.3% THC is a scheduled controlled substance (i.e., marijuana), while plant material that contains less than or equal to 0.3% is a legal agricultural commodity (i.e., hemp).[2] This definition has imposed severe challenges on crime labs. Among them is the dramatic increase in workload that results from the need to analyze and quantify the THC content of all C. sativa samples so that seized material can be appropriately designated. This is a time-consuming and resource-intensive enterprise that to greater and greater extents is consuming even larger forensic lab resources. Furthermore, defining the error cutoff for the 0.3% designation presents a challenge for the analysis of samples whose THC level is at the threshold.

Traditionally, hemp and marijuana plant materials are differentiated by determining the THC content through chromatography-based approaches such as gas chromatography-flame ionization detection (GC-FID) and gas chromatography-mass spectrometry (GC–MS)[3], in addition to high-performance liquid chromatography (HPLC) coupled to ultraviolet (UV) detection.[4] However, to accurately determine the THC content with these approaches, THC must be separated from all other components in the material (i.e., cannabinoids, terpenes, etc.) prior to quantification. One way to achieve this is to extend run times to allow for baseline separation between cannabinoids and other analytes present. Another option is to introduce a chemical derivatization step into the sample preparation protocol (which can be time-consuming), to differentiate between cannabinoids and their corresponding cannabinoid acids (e.g., THC and tetrahydrocannabinolic acid [THCA]). Although many investigations have been successful at differentiating between hemp and marijuana varieties or strains[5][6][7][8], the methods are reliant upon chromatography and are therefore susceptible to the aforementioned delineated challenges that can arise using this technique (i.e., lengthy run times, column contamination, etc.). Research towards developing, optimizing, and validating methods suitable for field testing of Cannabis materials has also been investigated.

Colorimetric tests represent a large percentage of these methods, which yield a presumptive result (by producing a color change)[9] when Cannabis-related substances are present, without the need for additional instrumentation (i.e., it is visible to the naked eye). Some examples include the 4-aminophenol test[10][11], Fast Blue BB test[11][12], and Duquenois-Levine test.[13] Similar to chromatography-based methods, these tests all rely upon the detection of THC specifically, which can complicate analyses because both marijuana and hemp contain this compound. Thus, while the distinction between marijuana and hemp has been defined based on THC levels, this is accompanied by several analytical challenges (i.e., baseline separation of molecules by chromatography-based methods, lengthy sample preparation protocols, and presumptive tests that can yield false positives[14], etc.).

An alternative less arbitrary approach is to base the distinction between them on the genome-defined differences in their metabolome signatures (i.e., small-molecule profiles). Studies utilizing the genetic profiles of Cannabis, such as genotyping-by-sequencing (GBS) and single-nucleotide polymorphisms (SNPs), have shown that, although they represent the same species, hemp and marijuana differ at the genome-wide level.[15][16][17] However, in addition to the fact that many crime laboratories are not positioned to integrate these types of analyses into current workflows, one of the bottlenecks to the routine use of the genome-defined small-molecule profiles for species attribution is the challenge of accessing this information quickly and reliably. One way to rapidly reveal this information, and subsequently distinguish between hemp and marijuana, is to combine an ambient ionization mass spectrometric technique—e.g., direct analysis in real time high-resolution mass spectrometry (DART-HRMS)[18]—with advanced statistical analysis. Ambient ionization methods (e.g., DART-HRMS, desorption electrospray ionization [DESI-MS]) have proven successful at screening for cannabinoids in Cannabis plant materials[19][20][21] and Cannabis-derived products (e.g., edibles, personal-care products, vape products, concentrates).[19][21] The unique capabilities of DART-HRMS are well-suited for the analysis of complex plant materials; the results are characterized by having high chemical information content, and little to no sample preparation prior to interrogating the materials is required. When applied to DART-HRMS-derived spectra, statistical data processing has enabled the successful differentiation of psychoactive plant species[22] and their headspace chemical signatures.[23] A modified version of DART-MS analysis introduced thermal desorption (TD) into the methodology (TD-DART-MS). One study utilized TD-DART-MS data to differentiate four hemp cultivars using PCA and partial least squares discriminant analysis (PLS-DA).[24] Another found that the application of statistical analysis to DART-MS data derived from methanolic extracts of hemp and marijuana samples revealed the potential for utilizing this method for optimally differentiating hemp and marijuana varieties.[25]

The study presented here, which is summarized in the scheme presented in Fig. 1, utilized DART-HRMS, for the first time, to investigate the complex genome-defined chemical fingerprints of hemp and marijuana (with no sample pretreatment) for the purpose of distinguishing between these two C. sativa varieties using multivariate statistical approaches. Advanced chemometrics was applied to the DART-HRMS data derived from commercial hemp, recreational marijuana, and marijuana samples from Drug Enforcement Administration (DEA)-registered suppliers to develop a robust model by which they (i.e., hemp and marijuana) could be readily differentiated. The success rate of the developed model’s ability to predict external validation samples was 100%, indicating a high level of certainty. Importantly, the developed method circumvents the need to separate and differentiate cannabinoids by chromatography techniques (i.e., the traditional forensic approach for determining the THC concentration in a sample and which is used for differentiating between hemp and marijuana), in addition to bypassing all sample pretreatment steps.


Fig1 Chambers JofCannRes23 5.png

Fig. 1. Workflow for discrimination of hemp and marijuana samples.

Materials and methods

Cannabis sativa plant materials

Twenty-nine C. sativa flower samples of the hemp variety were purchased from three online vendors: (1) CBD Hemp Direct (Las Vegas, Nevada, USA), (2) Berkshire CBD (Brattleboro, Vermont, USA), and (3) Plain Jane (Berkeley, California, USA). These samples were used to build the model (i.e., training set). An additional 12 samples of hemp plant material were purchased from Plain Jane (Medford, Oregon, USA) at a later date to test the model (i.e., they were used for external validation). Additional information (e.g., cultivar/strain, vendor, batch number) for these hemp materials is provided (see Additional file 1).

C. sativa plant material of the marijuana variety was obtained from two DEA-registered sources. The National Institute on Drug Abuse (NIDA) (Research Triangle Park (RTP), North Carolina, USA) Drug Supply Program, which is part of the National Institutes of Health (NIH), provided the following four samples (i.e., cultivars) with varying levels of THC and cannabidiol (CBD) (the major non-psychoactive constituent in C. sativa): 1 g low THC cultivar (low THC/very high CBD), 1 g medium THC cultivar (medium THC/medium CBD), 1 g high THC cultivar (high THC/low CBD), and 1 g very high THC cultivar (very high THC/low CBD). The National Institute of Standards and Technology (NIST) (Gaithersburg, Maryland, USA) provided eight 0.5 g samples of marijuana that were confiscated by local law enforcement at different times over the past few years. Twenty-one strains of recreational marijuana were purchased from Garden Remedies Marijuana Dispensary (Melrose, Massachusetts, USA). Ten of the recreational samples were randomly selected for use in the development of the training model, while the remaining 11 samples were used to test the model (i.e., for external validation). Information for all marijuana samples (e.g., sample name, brand, supplier/vendor, batch number, etc.) is available (see Additional file 1).

Mass spectral acquisition and analysis of DART-HRMS-derived data

The collection of mass spectral data was achieved by employing DART-HRMS. Two DART-HRMS instruments were used: (1) mass spectral data collected for all hemp products and the marijuana samples from DEA-registered suppliers were analyzed using the DART-HRMS instrument at the University at Albany (UAlbany) (Albany, New York, USA) and were translated and calibrated prior to data processing; and (2) all recreational marijuana flower samples were analyzed at IonSense Inc. (Saugus, Massachusetts, USA), with the raw data files calibrated, processed, and evaluated at UAlbany. The DART SVP (simplified voltage and pressure) ion source at IonSense was coupled to a JEOL AccuTOF high-resolution time-of-flight (TOF) mass spectrometer (Peabody, Massachusetts, USA) with a resolving power of 6000 full width at half maximum (FWHM) and mass accuracy of 5 millimass units (mmu). Data were collected in positive-ion mode using a DART ion source grid voltage of 300 V with the following mass spectrometer settings: ring lens, 5 V; orifice 1, 20 V; orifice 2 voltage, 5 V; peak voltage, 600 V; and detector voltage, 2000 V. The DART SVP ion source at UAlbany was also coupled to a JEOL AccuTOF high-resolution TOF mass spectrometer. The only difference between the DART ion source settings used at the two facilities was that the grid voltage at UAlbany was 250 V instead of 300 V. All mass spectral data were collected at a DART gas temperature of 350 °C using ultra-high purity helium gas at a flow rate of 2 L/min. Mass spectra were collected at a rate of 1 spectrum per second over a mass range of m/z 60–1000. TSSPro 3.0 software from Shrader Software Solutions (Grosse Pointe, Michigan, USA) was used for the calibration, spectral averaging, background subtraction, and peak centroiding of mass spectral data. Polyethylene glycol (PEG 600) (Sigma Aldrich, St. Louis, Missouri, USA) was used as the mass calibrant for all samples. Processing of the mass spectra of hemp and marijuana samples was performed with the Mass Mountaineer software suite from RBC Software (Portsmouth, New Hampshire, USA).

Multivariate data analysis

The workflow which extended from DART-HRMS data collection to multivariate data analysis is displayed in Fig. 1. In Step 1, DART mass spectra of the C. sativa samples representing hemp and marijuana varieties were acquired. The spectra in the form of text files were imported into MATLAB 9.9.0, R2020b Software (The MathWorks, Inc., Natick, Massachusetts, USA) and R 3.5.1 (R Core Team 2018) for analysis. Each text file was comprised of a two-column matrix of m/z values and their corresponding abundances (i.e., ion counts). In Step 2, peaks were aligned along common m/z values by histogram estimation and nearest-neighbor correction methods using the “mspalign” function in MATLAB. The generated matrix contained the aligned spectra for the replicates of hemp and marijuana samples. The replicates for each sample were averaged, normalized, transformed (with log 10), and subjected to unsupervised (Step 3) and supervised analyses (Step 4). As shown in Step 3, PCA[26] and k-means[27][28] were used to recognize the similarity and dissimilarity patterns of the samples and to reveal possible clusters, respectively. Silhouette width indexes were calculated to indicate the optimal number of clusters characterized by k-means and to validate the goodness of the clustering results. The data matrix was analyzed using supervised random forest (RF)[29][30] (Step 4) to create a model for differentiating hemp and marijuana plant materials. RF is an ensemble of individual tree predictors, in which each tree in the forest is grown based on the independent replicas of training samples and variables. The samples not included in the replicates for a given tree (1/3 of the original dataset) are termed “out-of-bag” (OOB) for that tree. The overall accuracy and performance characteristics of the discrimination model were estimated based on the predictions of OOB observations and external validation samples.

Results

DART-HRMS analysis of Cannabis sativa plant material

Initial investigations of C. sativa plant material focused on obtaining the DART-HRMS chemical profiles for both hemp and marijuana flower samples. Detailed information about the samples, including variety, cultivar/strain, vendor, and the batch number (when available) is provided (see Additional file 1). All samples were analyzed by inserting the closed end of a glass melting point capillary tube into the material and presenting the coated surface into the DART gas stream for approximately five seconds. A total of 29 hemp strains (i.e., cultivars) were purchased from three vendors at the beginning of this study, which included 27 CBD flower products and two cannabigerol (CBG) flower products. CBD flower contains high levels of CBD and cannabidiolic acid (CBDA), while CBG flower contains high levels of CBG and cannabigerolic acid (CBGA). An additional 12 hemp samples were purchased at a later date to test the developed model. Utilizing DART-HRMS is optimal for analyzing hemp and marijuana samples in their native forms (i.e., with no sample pretreatment, such as a decarboxylation step) to rapidly obtain the small-molecule profiles (i.e., in under one minute). The DART-HR mass spectra of all hemp flower samples (training-set hemp and test-set hemp) collected in positive-ion mode under soft ionization conditions (20 V) are available (see Additional file 2).

Figure 2 shows representative DART-HR mass spectra acquired in positive-ion mode from analysis of C. sativa plant materials, including CBD (panel A) and CBG (panel D) hemp flower samples. The DART-HR mass spectra of all CBD hemp flower samples are very similar to one another; protonated masses consistent with CBD and CBDA were detected at m/z 315 and 359, respectively, in all samples. DART-HRMS analysis of the two CBG hemp flower samples also yielded these peaks, in addition to peaks at nominal m/z 317 and 361, which are consistent with the protonated masses of CBG and CBGA, respectively. The DART-HR mass spectra of the CBG hemp flower samples retained similarities with the CBD hemp flower profiles. However, indicative of the high CBG levels reported in the CBG flower samples, the relative intensities of the peaks attributed to CBG and CBGA were much higher in the DART-HR mass spectra of the CBG flower products.


Fig2 Chambers JofCannRes23 5.png

Fig. 2. Representative DART-HR mass spectra of commercial hemp flower (panels A and D), marijuana samples supplied by NIST (panel B) and NIDA (panel E), and recreational marijuana flower products (panels C and F). Peaks consistent with the protonated masses of THC/CBD, CBG, THCA/CBDA, and CBGA at nominal m/z 315, 317, 359, and 361, respectively, were detected in the various samples.

C. sativa plant material of the marijuana variety was acquired from two U.S. DEA-registered sources: (1) NIDA supplied four marijuana samples (approximately 1 g each) through the NIDA/NIH Drug Supply Program; and (2) NIST provided eight marijuana samples (0.5 g each). All 12 marijuana samples were received in powdered form and were analyzed by DART-HRMS in positive-ion mode using the capillary tube sampling technique. Figure 2 presents two spectra of representative NIST (panel B) and NIDA (panel E) marijuana materials. Commercially available recreational marijuana samples were also analyzed. The DART-HR mass spectra for all marijuana samples from these suppliers are available (see Additional file 2). In total, 21 recreational marijuana samples were purchased from the Garden Remedies Marijuana Dispensary Adult-Use Menu. These products spanned the various marijuana strain types available (i.e., indica-dominant, sativa-dominant, hybrid), which represent C. sativa subspecies. Figure 2 presents two representative DART-HR mass spectra for indica (panel C) and sativa (panel F) dominant flower samples. The mass spectral profiles of all recreational marijuana flower products are available (see Additional file 2). Ten of the samples were randomly selected for inclusion in the training model. The remaining 11 recreational flower samples were used to test the prediction ability of the model (i.e., for external validation).

Differentiation of hemp and marijuana varieties of C. sativa

The aim of this work was to accomplish the following: (1) develop a rapid, easy-to-use, and efficient means by which to differentiate hemp and marijuana varieties of C. sativa, and by extension, a method to identify C. sativa unknowns; and (2) circumvent some of the challenges typically encountered during the analysis of C. sativa materials when using chromatography-based methods. The approach is founded on the hypothesis that inherent in the small-molecule profiles of hemp and marijuana is the necessary information for the differentiation of these Cannabis varieties. Prior to the application of multivariate analysis methods to the features of the DART-HRMS-derived chemical profiles of hemp and marijuana, the spectra of all samples were binned to create a common m/z reference vector to ease their comparison. Accordingly, the “mspalign” function in MATLAB was performed with a hist resolution parameter of 0.01, while the peak relative abundance cutoff threshold was set to 0.1% of the maximum intensity to detect all potentially significant peaks. The marijuana samples provided by NIDA and NIST were packaged in plastic bags, the composition of which contributed to the DART-HRMS profiles of the samples. Thus, the m/z values derived from the packaging (e.g., nominal m/z 59, 75, 89, 107, 127) were removed from the data. Another m/z value that was removed was nominal m/z 371, which has been previously shown to be a plasticizer present on the capillary tubes used for sampling.[31] The resulting matrix had dimensions of 430 × 390 and contained the aligned spectra for the five replicates of each of the 41 hemp samples, the five replicates of each of the 21 recreational marijuana samples, and the 10 replicates of each of the 12 marijuana samples supplied by NIDA and NIST. The results of the preliminary PCA analysis were examined by Q residuals and Hotelling’s T2 statistic to detect any outliers, and this resulted in three spectra being removed from the data. Outlier spectra included those whose acquisition was accompanied by poor mass calibration or those that were not representative of a typical chemical profile. The averaging of sample replicates resulted in a matrix with dimensions of 74 × 390. Following logarithm transformation, the matrix was subjected to further analysis. Figure 3 panel A presents the PCA results as a 2-dimensional (2D) score plot, where the color-coded classes appear in the coordinate space represented by the first two principal components (PCs), which cover 41% of the data variance. While the recreational marijuana samples (cyan triangles) are located in close proximity to the NIDA-supplied marijuana sample that was reported to contain medium levels of both THC and CBD, they were distant from the other NIDA and NIST samples. These results support previous studies that indicated differences between marijuana sold at dispensaries, and that provided for research purposes by DEA-registered suppliers.[17][32] Clustering by k-means using one minus correlation metrics resulted in the categorization of the hemp samples into one cluster (magenta circles) and the marijuana samples into the other cluster (cyan circles).


Fig3 Chambers JofCannRes23 5.png

Fig. 3. 2D score plot resulting from PCA of hemp and marijuana sample spectra (panel A); 2D score plot of multidimensional scaling (MDS) analysis of the proximity matrix resulting from the application of supervised random forest (panel B). The magenta and cyan colors represent hemp and marijuana, respectively. The cyan triangles show the subset of recreational marijuana samples.

Even though the DART-HR mass spectra of hemp and marijuana plant materials are readily visually apparent, a more objective approach to the assessment of the identity of C. sativa material was devised, using the random forest algorithm. This was applied to the 74 × 390 matrix. A total of 33 flower samples (12 hemp and 11 marijuana) of the 74 total C. sativa samples were randomly selected for external validation to examine the ability of the model to accurately predict the class assignments for new sample unknowns. The number of variables (which were randomly sampled as candidates at each split), and the number of trees found to be optimal were 20 and 500, respectively. Figure 3, panel B displays the proximity matrix generated from using supervised RF with a multidimensional scaling (MDS) method to show the pairwise similarities in a 2D Cartesian space, with the magenta and cyan points corresponding to the hemp and marijuana samples, respectively. It demonstrates the number of times that observations ended up in the same leaf node. According to Figure 3, panel B, although the NIDA marijuana sample reported as low THC/very high CBD is located between the two groups, the samples belonging to each group are close together and separated from the samples of the other group.

The optimal number of clusters was estimated by computing the average silhouette (which measures the quality of the clustering) of observations for different numbers of clusters. Figure 4, panel A displays the average silhouette width over a range of the possible number of clusters. The optimal number of clusters is the one that maximizes the average silhouette width. Based on the information provided in Figure 4, panel A, the optimal number of clusters is two. The silhouette plot in Figure 4, panel B displays silhouette coefficients for each sample when the data are split into two clusters. The silhouette width of each sample is a measure of how similar each sample is to its respective cluster in comparison to the other cluster. As shown in Figure 4, the optimum number of clusters is two: cluster 1 (magenta) has 40 members with a mean width of 0.23, and cluster 2 (cyan) has 34 members with a mean width of 0.45. Cluster 1 and cluster 2 members correspond to the samples of hemp and marijuana, respectively. One hemp sample was falsely clustered with the marijuana samples. The average silhouette width for the cluster of marijuana samples is higher than the average silhouette width for the hemp samples. This demonstrates that the cluster of marijuana samples is denser and that the samples are more similar to one another.


Fig4 Chambers JofCannRes23 5.png

Fig. 4. The average silhouette width over a range of cluster numbers (2–6) reveals that the optimum number of clusters is 2 (panel A). A silhouette plot (i.e., the visualization of the silhouette width for each sample) reveals the results with two clusters (panel B). Cluster 1 contains 40 members and cluster 2 contains 34 members. Hemp samples are shown in magenta, while marijuana samples are shown in cyan.

To reveal the model’s ability to distinguish between hemp and marijuana samples, Table 1 presents the confusion matrix for the prediction of OOB samples, while Table 2 contains the performance characteristics of the model (accuracy, sensitivity, specificity, and precision) for predicting the OOB samples. According to this table, the model performed well and the accuracy for predicting OOB samples is 98%.

Table 1. Confusion matrix associated with the prediction of “out-of-bag” samples in the random forest model.
Confusion matrix Prediction
Hemp Marijuana
True Hemp (29) 1.00 0.00
Marijuana (22) 0.04 0.96
Table 2. Performance results of the random forest model for prediction of “out-of-bag” and external validation samples.
Out-of-bag samples
Accuracy: 0.98 (98%)
Sensitivity Specificity Precision
Hemp (29) 1.00 0.96 0.97
Marijuana (22) 0.96 1.00 1.00
External C. sativa plant materials
Accuracy: 1.00 (100%)
Sensitivity Specificity Precision
Hemp (12) 1.00 1.00 1.00
Marijuana (11) 1.00 1.00 1.00

Classification of external C. sativa plant materials

The remaining 11 recreational marijuana flower products that were not included in the training set, in addition to the 12 hemp products purchased after the model had been developed, were screened against the model to test its ability to classify samples that were unknown to the model. Table 3 shows the confusion matrix results for the prediction of the test samples (i.e., for external validation). In addition, Table 2 shows the performance characteristics of the model for predicting the external C. sativa samples, with all performance merits equal to 1 for both test sample sets (i.e., hemp and marijuana). The information presented in Tables 1, 2, and 3 reveal that the model is well-fitted for discriminating the two C. sativa varieties.

Table 3. Confusion matrix associated with the prediction of external validation samples using a random forest model.
Confusion matrix Prediction
Hemp Marijuana
True Hemp (12) 1.00 0.00
Marijuana (11) 0.00 1.00

Discussion

The most common methods for differentiating hemp and marijuana plant materials are chromatography-based approaches (e.g., GC-FID, GC–MS, HPLC–UV)[3][4], with the categorization based upon THC content. Several reports have emphasized the use of GC-FID[8][33][34][35][36][37] and GC–MS[33][36][38][39][40][41][42] methods for detection of natural cannabinoids (among other Cannabis-derived molecules) in various Cannabis plant materials. Modifications to standard GC-FID and GC–MS protocols include GC-vacuum UV (VUV) spectroscopy[43], two-dimensional GC-FID (GCxGC-FID)[44], and GCxGC-MS with multivariate curve resolution-alternating least squares (MCR-ALS).[45] However, these methods rely upon the quantification of THC, which can be plagued with a number of analytical challenges, such as baseline separation of peaks and lengthy sample preparation protocols.

In an effort to circumvent the need to extend run times or incorporate extra sample preparation steps, several studies have investigated alternative sample collection techniques coupled with chromatography-based methods to differentiate C. sativa varieties. One study demonstrated the use of capillary microextraction of volatiles (CMV) coupled with GC–MS to distinguish the headspace volatiles of marijuana and hemp products based on their apparently distinct volatile organic compound (VOC) profiles.[5] However, this report revealed that potential adulterants and inconsistent packaging of samples may have contributed to the observed distinctions.[5] Another study utilized GC–MS coupled with dispersive pipette extraction (DPX) to investigate forensic casework marijuana and donated hemp samples.[6] Although the approach was successful at differentiating the two varieties with greater than 98% accuracy, a significant reduction of THC stability after 48 hours indicated that the samples would need to be reanalyzed if there was a delay between sample preparation and instrumental analysis.[6] Another GC-based study sought to differentiate hemp and marijuana through their cannabinoid and terpene profiles using GC-FID and principal component analysis (PCA).[7] This study, which included two recreational cultivars and three pharmacy Cannabis samples, successfully distinguished between the two C. sativa varieties.[7] In this case, expanding the sample source diversity could strengthen the ability of the model to classify a wider range of Cannabis samples. Another study applied PCA algorithms to quantitative data acquired from high-performance liquid chromatography-mass spectrometry (HPLC–MS) analysis of Cannabis plant materials.[8] This study identified several cannabinoids essential for differentiating between Cannabis strain types[8] (i.e., strains within the marijuana variety) as opposed to specifically targeting the cannabinoids essential to differentiating C. sativa varieties (i.e., hemp and marijuana), which would be important for criminal justice purposes in the U.S. Although many of these investigations were successful at differentiating between hemp and marijuana varieties or strains, the methods are reliant upon chromatography and are therefore susceptible to the aforementioned delineated challenges that can arise using this technique (i.e., lengthy run times, column contamination, etc.).

Non-chromatographic approaches that circumvent the requirement to separate and/or differentiate between cannabinoids have also been investigated for distinguishing hemp and marijuana. A hand-held Raman spectrometer coupled with orthogonal partial least squares-discriminant analysis (OPLS-DA) tools proved successful in differentiating between the two C. sativa varieties.[46] However, “real” forensic casework samples are rarely received in pristine form, and as such, the Raman approach is susceptible to interferences from various components that may be associated with the complex matrix and interfere with the Raman signal. Another study utilized advanced statistical modeling of nuclear magnetic resonance (NMR) spectroscopy and mass spectral data of C. sativa extracts[47], which is unique in that it is typically difficult to utilize NMR for the analysis of complex matrices and mixtures. Although effective, this instrumentation is not commonly found in forensic or other Cannabis analysis laboratories due to expensive start-up and maintenance costs.

Colorimetric tests are also commonly used for differentiating between hemp and marijuana varieties of Cannabis, especially in forensic fieldwork, and these do not generally require instrumental analysis to arrive at a presumptive identification. A validated method utilizing the 4-aminophenol color test to differentiate hemp and marijuana revealed some degree of success.[10] However, this test can yield inconclusive results with samples that have THC and CBD levels that are within a factor of three of one another.[10] Another common color test for the identification of marijuana samples is the Fast Blue BB (FBBB) colorimetric test, which reacts with the cannabinoids present in Cannabis (primarily THC). A study utilizing this test found that hemp and marijuana plant materials could be classified correctly when linear discriminant analysis (LDA) was used to develop a model based on RGB (red, green, blue) numerical codes from both fluorescence and color images that resulted from the application of the FBBB color test.[12] Positive-ion mode electrospray ionization Fourier transform-ion cyclotron resonance mass spectrometry (ESI( +)FT-ICR MS, ESI( +)MS/MS, ultraviolet–visible (UV–Vis) spectroscopy, and thin-layer chromatography (TLC) techniques have been used to investigate the products (i.e., chromophores) resulting from the application of the FBBB test to marijuana samples.[48] In addition, direct analysis in real time-mass spectrometry (DART-MS) and 1H NMR techniques were coupled to identify the chromophores produced when various cannabinoids react with the FBBB reagent.[49] A third color test to identify marijuana through the presence of THC is the Duquenois-Levine test. Research has been conducted to characterize (by mass spectrometry) the chromophores formed when cannabinoids react with the Duquenois reagents.[13][50][51] Similar to the chromatography-based methods described, these tests all rely upon detection of THC specifically, which can complicate analyses because both marijuana and hemp contain this compound. Thus, while the distinction between marijuana and hemp has been defined based on THC levels, this is accompanied by the several aforementioned analytical challenges. By using the entire metabolomic profiles of hemp and marijuana acquired through ambient ionization mass spectrometry, the method presented here does not rely solely on the presence of any one molecule (or set of molecules), ratios of molecules to one another, or the ability to differentiate between cannabinoid isomers (i.e., THC and CBD).

The overall results of this study reveal that DART-HRMS yields consistent and unique chemical profiles for analyzed Cannabis materials that enable hemp and marijuana samples to be accurately differentiated, while circumventing challenges typically encountered with traditional chromatography methods (difficulties with cannabinoid separation and extensive sample preparation) and presumptive color tests (inconclusive or false positive results). Furthermore, this study utilized a sample set that demonstrates a balance between the total number of samples included, the number of replicates obtained, and a diversity in sources from which the C. sativa materials were acquired. This research provides a strong foundation upon which to develop a comprehensive mass spectral database for identifying unknown C. sativa variants through the acquisition of their DART-HR mass spectra. While the approach does not aim to replace confirmatory testing for THC concentrations, the model accomplishes the following: (1) bypasses the typical sample preparation steps required for analyzing materials by chromatography-based methods that seek to differentiate the samples through separation of their constituent cannabinoids; (2) reduces the chances for false positives that can result from presumptive color tests; and (3) serves as a supplementary tool for forensic investigators that enables more targeted confirmatory testing.

This is timely and highly relevant, given the introduction in the U.S. House of Representatives of the “H.R.6645 – Hemp Advancement Act of 2022” bill.[52] This act aims to amend the current federal ruling regarding hemp by: (1) changing the 0.3% [THC] designation to 1% and (2) replacing the word “delta-9” with the word “total” to include the various isomers of THC that have emerged in recent years.[52] The introduction of this bill underscores some of the disadvantages of utilizing THC cutoffs in particular as the sole means by which to identify hemp and marijuana. Among other issues, it upends well-established and long-standing practices in criminalistics in a fashion that is expensive to address, since it will require the development of an entirely new set of protocols and data processing steps. Furthermore, it may not stand the test of time, as the cutoff thresholds are subject to change in the future. A method such as the one presented here, and which does not solely rely upon a 0.3% THC cutoff, is not at risk of becoming outdated upon further advancements of this bill or others in the U.S. House and Senate.

Conclusions

A combined ambient ionization mass spectrometric (i.e., DART-HRMS) and chemometric approach was successfully used to create a prediction model that facilitated rapid high-accuracy differentiation of C. sativa hemp and marijuana plant materials obtained from multiple sources (i.e., commercial, DEA-registered, recreational). This method, which circumvents sample pretreatment steps (i.e., solvent extractions), addresses some of the difficulties encountered when analyzing samples using more conventional forensic analysis methodologies. A primary example of this is eliminating the need to separate and differentiate cannabinoids by chromatography techniques in order to determine the sample’s THC content, which is the primary basis for distinguishing between hemp and marijuana varieties of Cannabis for most methods. When new hemp and recreational marijuana flower products were screened against the model developed in this study, 100% accuracy in prediction was observed. The identities of m/z values that were determined to be important for the optimal differentiation of hemp and marijuana are the subject of continuing investigations. In addition, it is possible that C. sativa materials (of either the hemp or marijuana variety) with atypical levels of minor cannabinoids (such as CBN or isomers of THC) may respond differently in the DART gas stream and that this, in turn, may influence the results predicted by the model. Therefore, samples such as these will be investigated (as was done with the analysis of the two CBG hemp flower samples), along with new samples/strains from commercial and DEA-registered suppliers as they become available so that the model reflects ongoing changes in the chemical profiles of Cannabis products on the market.

Supplementary information

  • Additional file 1 (.docx): Supplementary Mass Spectral Data and Sample Information. (1) Information about C. sativa plant materials analyzed in this study.
  • Additional file 2 (.docx): Supplementary Mass Spectral Data for C. sativa Materials. (1) DART-HR mass spectra for hemp and marijuana materials.

Abbreviations, acronyms, and initialisms

  • 2D: two-dimensional
  • CBD: cannabidiol
  • CBDA: cannabidiolic acid
  • CBG: cannabigerol
  • CBGA: cannabigerolic acid
  • CMV: capillary microextraction of volatiles
  • DART: direct analysis in real-time
  • DEA: U.S. Drug Enforcement Administration
  • DESI: desorption electrospray ionization
  • DPX: dispersive pipette extraction
  • ESI: electrospray ionization
  • FBBB: Fast Blue BB
  • FID: flame ionization detection
  • FT: Fourier transform
  • FWHM: full width at half maximum
  • GBS: genotyping-by-sequencing
  • GC: gas chromatography
  • GCxGC: two-dimensional gas chromatography
  • HR: high-resolution
  • HRMS: high-resolution mass spectrometry
  • HPLC: high-performance liquid chromatography
  • ICR: ion cyclotron resonance
  • LDA: linear discriminant analysis
  • MCD-ALS: multivariate curve resolution-alternating least squares
  • MDS: multidimensional scaling
  • mmu: millimass unit
  • MS: mass spectrometry
  • NIDA: ational Institute on Drug Abuse
  • NIH: ational Institutes of Health
  • NIJ: ational Institute of Justice
  • NIST: ational Institute of Standards and Technology
  • NMR: uclear magnetic resonance
  • OOB: ut-of-bag
  • OPLS-DA: rthogonal partial least squares-discriminant analysis
  • PC: rincipal component
  • PCA: rincipal component analysis
  • PEG: olyethylene glycol
  • PLS-DA: artial least squares-discriminant analysis
  • RF: andom forest
  • RGB: ed, green, blue
  • RTP: Research Triangle Institute
  • SNP: single-nucleotide polymorphism
  • SUNY: State University of New York
  • SVP: simplified voltage and pressure
  • TD: thermal desorption
  • THCA: delta-9-tetrahydrocannabinolic acid
  • TLC: thin-layer chromatography
  • TOF: time-of-flight
  • UAlbany: The University at Albany
  • UV: ultraviolet
  • Vis: visible
  • VOC: volatile organic compounds
  • VUV: vacuum UV
  • 9-THC or THC: delta-9-tetrahydrocannabinol

Acknowledgements

Thanks are extended to the National Institute on Drug Abuse/National Institutes of Health (NIDA/NIH) and the National Institute of Standards and Technology (NIST) for supplying Cannabis sativa marijuana samples analyzed in this study. Thanks are extended to IonSense, Inc. for the analysis of recreational Cannabis flower products and to Dr. Brent Wilson (NIST) for helpful assistance.

Author contributions

RAM conceived of the project, data analysis, project design, and project management and drafted the manuscript; MIC contributed to the experimental work and data analysis and drafted the manuscript; SB contributed to the data processing and data analysis and drafted the manuscript; BG contributed to the experimental work and data analysis. The authors read and approved the final manuscript.

Funding

The financial support of the National Institute of Justice (NIJ), Office of Justice programs, U.S. Department of Justice (DOJ) under Grant Nos. 2015-DN-BX-K057, 2017-R2-CX-0020 and 2019-BU-DX-0026 to RAM; the U.S. National Science Foundation (NSF) under Grant No. 1429329 to RAM; the 2020 Northeastern Association of Forensic Scientists (NEAFS) Carol De Forest Research Grant to MIC; the Initiatives for Women Foundation (IFW) Karen R. Hitchcock New Frontiers award to MIC; and the Research Foundation of SUNY are gratefully acknowledged. The opinions, findings, and conclusions or recommendations expressed in this publication are those of the authors and do not necessarily reflect those of the DOJ and/or the NSF.

Availability of data and materials

The datasets analyzed in the current study are available upon request at the discretion of the corresponding author.

Competing interests

The authors declare that they have no competing interests.


References

  1. National Institute of Justice (April 2019). "Report to Congress: Needs Assessment of Forensic Laboratories and Medical Examiner/Coroner Offices" (PDF). U.S. Department of Justice. pp. 86–97. https://www.ojp.gov/pdffiles1/nij/253626.pdf. 
  2. Conaway, K.M. (20 December 2018). "H.R.2 - Agriculture Improvement Act of 2018". Congress.gov. Library of Congress. https://www.congress.gov/bill/115th-congress/house-bill/112. 
  3. 3.0 3.1 Pourseyed Lazarjani, Masoumeh; Torres, Stephanie; Hooker, Thom; Fowlie, Chris; Young, Owen; Seyfoddin, Ali (1 December 2020). "Methods for quantification of cannabinoids: a narrative review" (in en). Journal of Cannabis Research 2 (1): 35. doi:10.1186/s42238-020-00040-2. ISSN 2522-5782. PMC PMC7819317. PMID 33526084. https://jcannabisresearch.biomedcentral.com/articles/10.1186/s42238-020-00040-2. 
  4. 4.0 4.1 Laboratory and Scientific Section, United Nations Office on Drugs and Crime (2009). "Recommended methods for the identification and analysis of cannabis and cannabis products" (PDF). ISBN 978-92-1-148242-3. https://www.unodc.org/documents/scientific/ST-NAR-40-Ebook_1.pdf. 
  5. 5.0 5.1 5.2 Wiebelhaus, Nancy; Hamblin, D’Nisha; Kreitals, Natasha M.; Almirall, Jose R. (1 November 2016). "Differentiation of marijuana headspace volatiles from other plants and hemp products using capillary microextraction of volatiles (CMV) coupled to gas-chromatography–mass spectrometry (GC–MS)" (in en). Forensic Chemistry 2: 1–8. doi:10.1016/j.forc.2016.08.004. https://linkinghub.elsevier.com/retrieve/pii/S2468170916300285. 
  6. 6.0 6.1 6.2 Horne, Melissa; Mastrianni, Kaylee R.; Amick, Gray; Hardy, Rachel; Renneker, Elissa; Miller, Kevin W.P. (1 September 2020). "Fast Discrimination of Marijuana using Automated High‐throughput Cannabis Sample Preparation and Analysis by Gas Chromatography–Mass Spectrometry" (in en). Journal of Forensic Sciences 65 (5): 1709–1715. doi:10.1111/1556-4029.14525. ISSN 0022-1198. https://onlinelibrary.wiley.com/doi/10.1111/1556-4029.14525. 
  7. 7.0 7.1 7.2 Pacula, Rosalie Liccardo; Jacobson, Mireille; Maksabedian, Ervant J. (1 June 2016). "In the weeds: a baseline view of cannabis use among legalizing states and their neighbours: In the weeds: a baseline view of cannabis use among legalizing states and their neighbours" (in en). Addiction 111 (6): 973–980. doi:10.1111/add.13282. PMC PMC5216038. PMID 26687431. https://onlinelibrary.wiley.com/doi/10.1111/add.13282. 
  8. 8.0 8.1 8.2 8.3 Fischedick, Justin Thomas; Hazekamp, Arno; Erkelens, Tjalling; Choi, Young Hae; Verpoorte, Rob (1 December 2010). "Metabolic fingerprinting of Cannabis sativa L., cannabinoids and terpenoids for chemotaxonomic and drug standardization purposes" (in en). Phytochemistry 71 (17-18): 2058–2073. doi:10.1016/j.phytochem.2010.10.001. https://linkinghub.elsevier.com/retrieve/pii/S003194221000381X. 
  9. Philp, Morgan; Shimmon, Ronald; Tahtouh, Mark; Fu, Shanlin (5 February 2018). "Color Spot Test As a Presumptive Tool for the Rapid Detection of Synthetic Cathinones" (in en). Journal of Visualized Experiments (132): 57045. doi:10.3791/57045. ISSN 1940-087X. PMC PMC5912360. PMID 29443096. https://www.jove.com/t/57045/color-spot-test-as-a-presumptive-tool-for-the-rapid-detection-of-synthetic-cathinones. 
  10. 10.0 10.1 10.2 Lewis, Kenna; Wagner, Rebecca; Rodriguez‐Cruz, Sandra E.; Weaver, Michael J.; Dumke, Jonathan C. (1 January 2021). "Validation of the 4‐aminophenol color test for the differentiation of marijuana‐type and hemp‐type cannabis" (in en). Journal of Forensic Sciences 66 (1): 285–294. doi:10.1111/1556-4029.14562. ISSN 0022-1198. https://onlinelibrary.wiley.com/doi/10.1111/1556-4029.14562. 
  11. 11.0 11.1 Acosta, Alexander; Li, Li; Weaver, Mike; Capote, Ryan; Perr, Jeannette; Almirall, José (1 December 2022). "Validation of a combined Fast blue BB and 4-Aminophenol colorimetric test for indication of Hemp-type and Marijuana-type cannabis" (in en). Forensic Chemistry 31: 100448. doi:10.1016/j.forc.2022.100448. https://linkinghub.elsevier.com/retrieve/pii/S2468170922000510. 
  12. 12.0 12.1 Acosta, Alexander; Almirall, José (1 December 2021). "Differentiation between hemp-type and marijuana-type cannabis using the Fast Blue BB colorimetric test" (in en). Forensic Chemistry 26: 100376. doi:10.1016/j.forc.2021.100376. https://linkinghub.elsevier.com/retrieve/pii/S2468170921000722. 
  13. 13.0 13.1 Forrester, D.E. (15 April 1997). "The Duquenois Color Test for Marijuana: Spectroscopic and Chemical Studies". Georgetown University. https://www.proquest.com/openview/b2bc339ba2dc1b9f55ab83bb9a112ec9/1?pq-origsite=gscholar&cbl=18750&diss=y. 
  14. Gabrielson, R.; Sanders, T. (7 July 2016). "Busted: Tens of thousands of people every year are sent to jail based on the results of a $2 roadside drug test. Widespread evidence shows that these tests routinely produce false positives. Why are police departments and prosecutors still using them?". ProPublica. https://www.propublica.org/article/common-roadside-drug-test-routinely-produces-false-positives. 
  15. Sawler, Jason; Stout, Jake M.; Gardner, Kyle M.; Hudson, Darryl; Vidmar, John; Butler, Laura; Page, Jonathan E.; Myles, Sean (26 August 2015). Tinker, Nicholas A.. ed. "The Genetic Structure of Marijuana and Hemp" (in en). PLOS ONE 10 (8): e0133292. doi:10.1371/journal.pone.0133292. ISSN 1932-6203. PMC PMC4550350. PMID 26308334. https://dx.plos.org/10.1371/journal.pone.0133292. 
  16. Roman, Madeline G.; Houston, Rachel (1 November 2020). "Investigation of chloroplast regions rps16 and clpP for determination of Cannabis sativa crop type and biogeographical origin" (in en). Legal Medicine 47: 101759. doi:10.1016/j.legalmed.2020.101759. https://linkinghub.elsevier.com/retrieve/pii/S1344622320300936. 
  17. 17.0 17.1 Schwabe, Anna L.; Hansen, Connor J.; Hyslop, Richard M.; McGlaughlin, Mitchell E. (29 September 2021). "Comparative Genetic Structure of Cannabis sativa Including Federally Produced, Wild Collected, and Cultivated Samples". Frontiers in Plant Science 12: 675770. doi:10.3389/fpls.2021.675770. ISSN 1664-462X. PMC PMC8544287. PMID 34707624. https://www.frontiersin.org/articles/10.3389/fpls.2021.675770/full. 
  18. Cody, Robert B.; Laramée, James A.; Durst, H. Dupont (1 April 2005). "Versatile New Ion Source for the Analysis of Materials in Open Air under Ambient Conditions" (in en). Analytical Chemistry 77 (8): 2297–2302. doi:10.1021/ac050162j. ISSN 0003-2700. https://pubs.acs.org/doi/10.1021/ac050162j. 
  19. 19.0 19.1 Chambers, Megan I.; Musah, Rabi A. (1 March 2022). "DART-HRMS as a triage approach for the rapid analysis of cannabinoid-infused edible matrices, personal-care products and Cannabis sativa hemp plant material" (in en). Forensic Chemistry 27: 100382. doi:10.1016/j.forc.2021.100382. https://linkinghub.elsevier.com/retrieve/pii/S2468170921000783. 
  20. Rodriguez-Cruz, Sandra E. (15 January 2006). "Rapid analysis of controlled substances using desorption electrospray ionization mass spectrometry" (in en). Rapid Communications in Mass Spectrometry 20 (1): 53–60. doi:10.1002/rcm.2267. ISSN 0951-4198. https://onlinelibrary.wiley.com/doi/10.1002/rcm.2267. 
  21. 21.0 21.1 Chambers, Megan I.; Musah, Rabi A. (1 May 2023). "DART-HRMS triage approach part 2 – Application to the detection of cannabinoids and terpenes in recreational Cannabis products" (in en). Forensic Chemistry 33: 100469. doi:10.1016/j.forc.2023.100469. https://linkinghub.elsevier.com/retrieve/pii/S246817092300005X. 
  22. Beyramysoltan, Samira; Abdul-Rahman, Nana-Hawwa; Musah, Rabi A. (1 November 2019). "Call it a “nightshade”—A hierarchical classification approach to identification of hallucinogenic Solanaceae spp. using DART-HRMS-derived chemical signatures" (in en). Talanta 204: 739–746. doi:10.1016/j.talanta.2019.06.010. https://linkinghub.elsevier.com/retrieve/pii/S0039914019306290. 
  23. Appley, Meghan Grace; Beyramysoltan, Samira; Musah, Rabi Ann (24 September 2019). "Random Forest Processing of Direct Analysis in Real-Time Mass Spectrometric Data Enables Species Identification of Psychoactive Plants from Their Headspace Chemical Signatures" (in en). ACS Omega 4 (13): 15636–15644. doi:10.1021/acsomega.9b02145. ISSN 2470-1343. PMC PMC6761758. PMID 31572865. https://pubs.acs.org/doi/10.1021/acsomega.9b02145. 
  24. Dong, Wen; Liang, Jian; Barnett, Isabella; Kline, Paul C.; Altman, Elliot; Zhang, Mengliang (1 December 2019). "The classification of Cannabis hemp cultivars by thermal desorption direct analysis in real time mass spectrometry (TD-DART-MS) with chemometrics" (in en). Analytical and Bioanalytical Chemistry 411 (30): 8133–8142. doi:10.1007/s00216-019-02200-7. ISSN 1618-2642. http://link.springer.com/10.1007/s00216-019-02200-7. 
  25. Pieslak, J.R. (2021). "Analytical techniques for the differentiation of hemp and marijuana". OpenBU. Boston University Libraries. https://hdl.handle.net/2144/43518. 
  26. Jolliffe, Ian T.; Cadima, Jorge (13 April 2016). "Principal component analysis: a review and recent developments" (in en). Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 374 (2065): 20150202. doi:10.1098/rsta.2015.0202. ISSN 1364-503X. PMC PMC4792409. PMID 26953178. https://royalsocietypublishing.org/doi/10.1098/rsta.2015.0202. 
  27. Sammut, Claude; Webb, Geoffrey I., eds. (2010) (in en). Encyclopedia of Machine Learning. Boston, MA: Springer US. doi:10.1007/978-0-387-30164-8. ISBN 978-0-387-30768-8. http://link.springer.com/10.1007/978-0-387-30164-8. 
  28. Lloyd, S. (1 March 1982). "Least squares quantization in PCM" (in en). IEEE Transactions on Information Theory 28 (2): 129–137. doi:10.1109/TIT.1982.1056489. ISSN 0018-9448. http://ieeexplore.ieee.org/document/1056489/. 
  29. Liaw, A.; Wiener, M. (2002). "Classification and Regression by randomForest" (PDF). R News 2 (3): 18–22. ISSN 1609-3631. https://cogns.northwestern.edu/cbmg/LiawAndWiener2002.pdf. 
  30. Breiman, Leo (2001). "Random Forests". Machine Learning 45 (1): 5–32. doi:10.1023/A:1010933404324. http://link.springer.com/10.1023/A:1010933404324. 
  31. Beyramysoltan, Samira; Ventura, Mónica I.; Rosati, Jennifer Y.; Giffen-Lemieux, Justine E.; Musah, Rabi A. (7 April 2020). "Identification of the Species Constituents of Maggot Populations Feeding on Decomposing Remains—Facilitation of the Determination of Post Mortem Interval and Time Since Tissue Infestation through Application of Machine Learning and Direct Analysis in Real Time-Mass Spectrometry" (in en). Analytical Chemistry 92 (7): 5439–5446. doi:10.1021/acs.analchem.0c00199. ISSN 0003-2700. https://pubs.acs.org/doi/10.1021/acs.analchem.0c00199. 
  32. Vergara, Daniela; Bidwell, L. Cinnamon; Gaudino, Reggie; Torres, Anthony; Du, Gary; Ruthenburg, Travis C.; deCesare, Kymron; Land, Donald P. et al. (19 April 2017). "Compromised External Validity: Federally Produced Cannabis Does Not Reflect Legal Markets" (in en). Scientific Reports 7 (1): 46528. doi:10.1038/srep46528. ISSN 2045-2322. PMC PMC5395929. PMID 28422145. https://www.nature.com/articles/srep46528. 
  33. 33.0 33.1 Zekič, Jure; Križman, Mitja (11 December 2020). "Development of Gas-Chromatographic Method for Simultaneous Determination of Cannabinoids and Terpenes in Hemp" (in en). Molecules 25 (24): 5872. doi:10.3390/molecules25245872. ISSN 1420-3049. PMC PMC7763075. PMID 33322595. https://www.mdpi.com/1420-3049/25/24/5872. 
  34. Dussy, Franz E.; Hamberg, Cornelia; Luginbühl, Marco; Schwerzmann, Thomas; Briellmann, Thomas A. (1 April 2005). "Isolation of Δ9-THCA-A from hemp and analytical aspects concerning the determination of Δ9-THC in cannabis products" (in en). Forensic Science International 149 (1): 3–10. doi:10.1016/j.forsciint.2004.05.015. https://linkinghub.elsevier.com/retrieve/pii/S0379073804003408. 
  35. Fischedick, Justin; Van Der Kooy, Frank; Verpoorte, Robert (2010). "Cannabinoid Receptor 1 Binding Activity and Quantitative Analysis of Cannabis sativa L. Smoke and Vapor" (in en). Chemical and Pharmaceutical Bulletin 58 (2): 201–207. doi:10.1248/cpb.58.201. ISSN 0009-2363. http://www.jstage.jst.go.jp/article/cpb/58/2/58_2_201/_article. 
  36. 36.0 36.1 Hazekamp, Arno; Simons, Ruud; Peltenburg‐Looman, Anja; Sengers, Melvin; van Zweden, Rianne; Verpoorte, Robert (1 January 2004). "Preparative Isolation of Cannabinoids from Cannabis sativa by Centrifugal Partition Chromatography" (in en). Journal of Liquid Chromatography & Related Technologies 27 (15): 2421–2439. doi:10.1081/JLC-200028170. ISSN 1082-6076. https://www.tandfonline.com/doi/full/10.1081/JLC-200028170. 
  37. Hazekamp, A.; Fischedick, J. T. (1 July 2012). "Cannabis - from cultivar to chemovar: Towards a better definition of Cannabis potency" (in en). Drug Testing and Analysis 4 (7-8): 660–667. doi:10.1002/dta.407. https://onlinelibrary.wiley.com/doi/10.1002/dta.407. 
  38. Hazekamp, Arno; Peltenburg, Anja; Verpoorte, Rob; Giroud, Christian (1 September 2005). "Chromatographic and Spectroscopic Data of Cannabinoids from Cannabis sativa L." (in en). Journal of Liquid Chromatography & Related Technologies 28 (15): 2361–2382. doi:10.1080/10826070500187558. ISSN 1082-6076. https://www.tandfonline.com/doi/full/10.1080/10826070500187558. 
  39. Namdar, Dvory; Mazuz, Moran; Ion, Aurel; Koltai, Hinanit (1 March 2018). "Variation in the compositions of cannabinoid and terpenoids in Cannabis sativa derived from inflorescence position along the stem and extraction methods" (in en). Industrial Crops and Products 113: 376–382. doi:10.1016/j.indcrop.2018.01.060. https://linkinghub.elsevier.com/retrieve/pii/S092666901830061X. 
  40. Namdar, Dvory; Charuvi, Dana; Ajjampura, Vinayka; Mazuz, Moran; Ion, Aurel; Kamara, Itzhak; Koltai, Hinanit (1 June 2019). "LED lighting affects the composition and biological activity of Cannabis sativa secondary metabolites" (in en). Industrial Crops and Products 132: 177–185. doi:10.1016/j.indcrop.2019.02.016. https://linkinghub.elsevier.com/retrieve/pii/S0926669019301086. 
  41. Omar, Jone; Olivares, Maitane; Alzaga, Mikel; Etxebarria, Nestor (1 April 2013). "Optimisation and characterisation of marihuana extracts obtained by supercritical fluid extraction and focused ultrasound extraction and retention time locking GC-MS: Gas Chromatography" (in en). Journal of Separation Science 36 (8): 1397–1404. doi:10.1002/jssc.201201103. https://onlinelibrary.wiley.com/doi/10.1002/jssc.201201103. 
  42. Knight, Glenys; Hansen, Sean; Connor, Mark; Poulsen, Helen; McGovern, Catherine; Stacey, Janet (1 October 2010). "The results of an experimental indoor hydroponic Cannabis growing study, using the ‘Screen of Green’ (ScrOG) method—Yield, tetrahydrocannabinol (THC) and DNA analysis" (in en). Forensic Science International 202 (1-3): 36–44. doi:10.1016/j.forsciint.2010.04.022. https://linkinghub.elsevier.com/retrieve/pii/S0379073810001969. 
  43. Leghissa, Allegra; Smuts, Jonathan; Qiu, Changling; Hildenbrand, Zacariah L.; Schug, Kevin A. (1 January 2018). "Detection of cannabinoids and cannabinoid metabolites using gas chromatography with vacuum ultraviolet spectroscopy" (in en). Separation Science Plus 1 (1): 37–42. doi:10.1002/sscp.201700005. https://onlinelibrary.wiley.com/doi/10.1002/sscp.201700005. 
  44. Gröger, Th.; Schäffer, M.; Pütz, M.; Ahrens, B.; Drew, K.; Eschner, M.; Zimmermann, R. (1 July 2008). "Application of two-dimensional gas chromatography combined with pixel-based chemometric processing for the chemical profiling of illicit drug samples" (in en). Journal of Chromatography A 1200 (1): 8–16. doi:10.1016/j.chroma.2008.05.028. https://linkinghub.elsevier.com/retrieve/pii/S0021967308008297. 
  45. Omar, Jone; Olivares, Maitane; Amigo, José Manuel; Etxebarria, Nestor (1 April 2014). "Resolution of co-eluting compounds of Cannabis Sativa in comprehensive two-dimensional gas chromatography/mass spectrometry detection with Multivariate Curve Resolution-Alternating Least Squares" (in en). Talanta 121: 273–280. doi:10.1016/j.talanta.2013.12.044. https://linkinghub.elsevier.com/retrieve/pii/S0039914013010370. 
  46. Sanchez, Lee; Filter, Conor; Baltensperger, David; Kurouski, Dmitry (2020). "Confirmatory non-invasive and non-destructive differentiation between hemp and cannabis using a hand-held Raman spectrometer" (in en). RSC Advances 10 (6): 3212–3216. doi:10.1039/C9RA08225E. ISSN 2046-2069. PMC PMC9048763. PMID 35497720. http://xlink.rsc.org/?DOI=C9RA08225E. 
  47. Chen, Zewei; Harrington, Peter de Boves (19 November 2019). "Pipeline for High-Throughput Modeling of Marijuana and Hemp Extracts" (in en). Analytical Chemistry 91 (22): 14489–14497. doi:10.1021/acs.analchem.9b03290. ISSN 0003-2700. https://pubs.acs.org/doi/10.1021/acs.analchem.9b03290. 
  48. dos Santos, Nayara A.; Souza, Lindamara M.; Domingos, Eloilson; França, Hildegardo S.; Lacerda, Valdemar; Beatriz, Adilson; Vaz, Boniek G.; Rodrigues, Rayza R.T. et al. (1 August 2016). "Evaluating the selectivity of colorimetric test (Fast Blue BB salt) for the cannabinoids identification in marijuana street samples by UV–Vis, TLC, ESI(+)FT-ICR MS and ESI(+)MS/MS" (in en). Forensic Chemistry 1: 13–21. doi:10.1016/j.forc.2016.07.001. https://linkinghub.elsevier.com/retrieve/pii/S2468170916300297. 
  49. França, Hildegardo S.; Acosta, Alexander; Jamal, Adeel; Romao, Wanderson; Mulloor, Jerome; Almirall, Jose R. (1 March 2020). "Experimental and ab initio investigation of the products of reaction from Δ9-tetrahydrocannabinol (Δ9-THC) and the fast blue BB spot reagent in presumptive drug tests for cannabinoids" (in en). Forensic Chemistry 17: 100212. doi:10.1016/j.forc.2019.100212. https://linkinghub.elsevier.com/retrieve/pii/S2468170919301092. 
  50. Jacobs, Alexander D.; Steiner, Robert R. (1 June 2014). "Detection of the Duquenois–Levine chromophore in a marijuana sample" (in en). Forensic Science International 239: 1–5. doi:10.1016/j.forsciint.2014.02.031. https://linkinghub.elsevier.com/retrieve/pii/S0379073814000929. 
  51. Watanabe, Kazuhito; Honda, Go; Miyagi, Takeaki; Kanai, Masataka; Usami, Noriyuki; Yamaori, Satoshi; Iwamuro, Yoshiaki; Chinaka, Satoshi et al. (1 January 2017). "The Duquenois reaction revisited: mass spectrometric estimation of chromophore structures derived from major phytocannabinoids" (in en). Forensic Toxicology 35 (1): 185–189. doi:10.1007/s11419-016-0337-6. ISSN 1860-8965. http://link.springer.com/10.1007/s11419-016-0337-6. 
  52. 52.0 52.1 Pingree, C. (1 November 2022). "H.R.6645 - Hemp Advancement Act of 2022". Congress.gov. Library of Congress. https://www.congress.gov/bill/117th-congress/house-bill/6645. 

Notes

This presentation is faithful to the original, with only a few minor changes to presentation. Some grammar and punctuation was cleaned up to improve readability. In some cases important information was missing from the references, and that information was added. The original lists references in alphabetical order; they are listed by order of appearance for this version, by design.