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

{{DISPLAYTITLE:Journal:Combined ambient ionization mass spectrometric and chemometric approach for the differentiation of hemp and marijuana varieties of ''Cannabis sativa''}}
{{Infobox journal article
|name         = 
|image        = 
|alt          = <!-- Alternative text for images -->
|caption      = 
|title_full   = Combined ambient ionization mass spectrometric and chemometric approach for the differentiation of hemp and marijuana varieties of ''Cannabis sativa''
|journal      = ''Journal of Cannabis Research''
|authors      = Chambers, Megan I.; Beyramysoltan, Samira; Garosi, Benedetta; Musah, Rabi A.
|affiliations = State University of New York
|contact      = Email: rmusah at albany dot edu
|editors      = 
|pub_year     = 2023
|vol_iss      = '''5'''
|at           = 5
|doi          = [https://doi.org/10.1186/s42238-023-00173-0 10.1186/s42238-023-00173-0]
|issn         = 2522-5782
|license      = [https://creativecommons.org/licenses/by/4.0/ Creative Commons Attribution 4.0 International]
|website      = [https://jcannabisresearch.biomedcentral.com/articles/10.1186/s42238-023-00173-0 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 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 [[Tetrahydrocannabinol|Δ<sup>9</sup>-tetrahydrocannabinol]] (THC), the primary [[Psychoactive drug|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 (material)|sample]] preparation to render the materials into [[Cannabis concentrate|extracts]] suitable for sample injection, for complete separation and differentiation of THC from all other analytes present. This can create problems for [[Forensic science|forensic]] [[Laboratory|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]] [[Resolution (mass spectrometry)|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 [[Cannabis strains|strains]]/[[cultivar]]s.

'''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 science|forensic]] [[Laboratory|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.<ref>{{Cite web |last=National Institute of Justice |date=April 2019 |title=Report to Congress: Needs Assessment of Forensic Laboratories and Medical Examiner/Coroner Offices |url=https://www.ojp.gov/pdffiles1/nij/253626.pdf |format=PDF |publisher=U.S. Department of Justice |pages=86–97}}</ref> 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 [[Tetrahydrocannabinol|Δ<sup>9</sup>-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).<ref>{{Cite web |last=Conaway, K.M. |date=20 December 2018 |title=H.R.2 - Agriculture Improvement Act of 2018 |work=Congress.gov |url=https://www.congress.gov/bill/115th-congress/house-bill/112 |publisher=Library of Congress}}</ref> 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'' [[Sample (material)|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]]-[[Chromatography detector|flame ionization detection]] (GC-FID) and [[gas chromatography-mass spectrometry]] (GC–MS)<ref name=":0">{{Cite journal |last=Pourseyed Lazarjani |first=Masoumeh |last2=Torres |first2=Stephanie |last3=Hooker |first3=Thom |last4=Fowlie |first4=Chris |last5=Young |first5=Owen |last6=Seyfoddin |first6=Ali |date=2020-12 |title=Methods for quantification of cannabinoids: a narrative review |url=https://jcannabisresearch.biomedcentral.com/articles/10.1186/s42238-020-00040-2 |journal=Journal of Cannabis Research |language=en |volume=2 |issue=1 |pages=35 |doi=10.1186/s42238-020-00040-2 |issn=2522-5782 |pmc=PMC7819317 |pmid=33526084}}</ref>, in addition to [[high-performance liquid chromatography]] (HPLC) coupled to [[Ultraviolet–visible spectroscopy|ultraviolet]] (UV) detection.<ref name=":1">{{Cite web |last=Laboratory and Scientific Section, United Nations Office on Drugs and Crime |date=2009 |title=Recommended methods for the identification and analysis of cannabis and cannabis products |url=https://www.unodc.org/documents/scientific/ST-NAR-40-Ebook_1.pdf |format=PDF |isbn=978-92-1-148242-3}}</ref> However, to accurately determine the THC content with these approaches, THC must be separated from all other components in the material (i.e., [[cannabinoid]]s, [[terpene]]s, 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 [[Cannabis strains|strains]]<ref name=":2">{{Cite journal |last=Wiebelhaus |first=Nancy |last2=Hamblin |first2=D’Nisha |last3=Kreitals |first3=Natasha M. |last4=Almirall |first4=Jose R. |date=2016-11 |title=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) |url=https://linkinghub.elsevier.com/retrieve/pii/S2468170916300285 |journal=Forensic Chemistry |language=en |volume=2 |pages=1–8 |doi=10.1016/j.forc.2016.08.004}}</ref><ref name=":3">{{Cite journal |last=Horne |first=Melissa |last2=Mastrianni |first2=Kaylee R. |last3=Amick |first3=Gray |last4=Hardy |first4=Rachel |last5=Renneker |first5=Elissa |last6=Miller |first6=Kevin W.P. |date=2020-09 |title=Fast Discrimination of Marijuana using Automated High‐throughput Cannabis Sample Preparation and Analysis by Gas Chromatography–Mass Spectrometry |url=https://onlinelibrary.wiley.com/doi/10.1111/1556-4029.14525 |journal=Journal of Forensic Sciences |language=en |volume=65 |issue=5 |pages=1709–1715 |doi=10.1111/1556-4029.14525 |issn=0022-1198}}</ref><ref name=":4">{{Cite journal |last=Pacula |first=Rosalie Liccardo |last2=Jacobson |first2=Mireille |last3=Maksabedian |first3=Ervant J. |date=2016-06 |title=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 |url=https://onlinelibrary.wiley.com/doi/10.1111/add.13282 |journal=Addiction |language=en |volume=111 |issue=6 |pages=973–980 |doi=10.1111/add.13282 |pmc=PMC5216038 |pmid=26687431}}</ref><ref name=":5">{{Cite journal |last=Fischedick |first=Justin Thomas |last2=Hazekamp |first2=Arno |last3=Erkelens |first3=Tjalling |last4=Choi |first4=Young Hae |last5=Verpoorte |first5=Rob |date=2010-12 |title=Metabolic fingerprinting of Cannabis sativa L., cannabinoids and terpenoids for chemotaxonomic and drug standardization purposes |url=https://linkinghub.elsevier.com/retrieve/pii/S003194221000381X |journal=Phytochemistry |language=en |volume=71 |issue=17-18 |pages=2058–2073 |doi=10.1016/j.phytochem.2010.10.001}}</ref>, 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)<ref>{{Cite journal |last=Philp |first=Morgan |last2=Shimmon |first2=Ronald |last3=Tahtouh |first3=Mark |last4=Fu |first4=Shanlin |date=2018-02-05 |title=Color Spot Test As a Presumptive Tool for the Rapid Detection of Synthetic Cathinones |url=https://www.jove.com/t/57045/color-spot-test-as-a-presumptive-tool-for-the-rapid-detection-of-synthetic-cathinones |journal=Journal of Visualized Experiments |language=en |issue=132 |pages=57045 |doi=10.3791/57045 |issn=1940-087X |pmc=PMC5912360 |pmid=29443096}}</ref> 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<ref name=":6">{{Cite journal |last=Lewis |first=Kenna |last2=Wagner |first2=Rebecca |last3=Rodriguez‐Cruz |first3=Sandra E. |last4=Weaver |first4=Michael J. |last5=Dumke |first5=Jonathan C. |date=2021-01 |title=Validation of the 4‐aminophenol color test for the differentiation of marijuana‐type and hemp‐type cannabis |url=https://onlinelibrary.wiley.com/doi/10.1111/1556-4029.14562 |journal=Journal of Forensic Sciences |language=en |volume=66 |issue=1 |pages=285–294 |doi=10.1111/1556-4029.14562 |issn=0022-1198}}</ref><ref name=":7">{{Cite journal |last=Acosta |first=Alexander |last2=Li |first2=Li |last3=Weaver |first3=Mike |last4=Capote |first4=Ryan |last5=Perr |first5=Jeannette |last6=Almirall |first6=José |date=2022-12 |title=Validation of a combined Fast blue BB and 4-Aminophenol colorimetric test for indication of Hemp-type and Marijuana-type cannabis |url=https://linkinghub.elsevier.com/retrieve/pii/S2468170922000510 |journal=Forensic Chemistry |language=en |volume=31 |pages=100448 |doi=10.1016/j.forc.2022.100448}}</ref>, Fast Blue BB test<ref name=":7" /><ref name=":8">{{Cite journal |last=Acosta |first=Alexander |last2=Almirall |first2=José |date=2021-12 |title=Differentiation between hemp-type and marijuana-type cannabis using the Fast Blue BB colorimetric test |url=https://linkinghub.elsevier.com/retrieve/pii/S2468170921000722 |journal=Forensic Chemistry |language=en |volume=26 |pages=100376 |doi=10.1016/j.forc.2021.100376}}</ref>, and Duquenois-Levine test.<ref name=":9">{{Cite web |last=Forrester, D.E. |date=15 April 1997 |title=The Duquenois Color Test for Marijuana: Spectroscopic and Chemical Studies |url=https://www.proquest.com/openview/b2bc339ba2dc1b9f55ab83bb9a112ec9/1?pq-origsite=gscholar&cbl=18750&diss=y |publisher=Georgetown University}}</ref> 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<ref>{{Cite web |last=Gabrielson, R.; Sanders, T. |date=07 July 2016 |title=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? |work=ProPublica |url=https://www.propublica.org/article/common-roadside-drug-test-routinely-produces-false-positives}}</ref>, 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.<ref>{{Cite journal |last=Sawler |first=Jason |last2=Stout |first2=Jake M. |last3=Gardner |first3=Kyle M. |last4=Hudson |first4=Darryl |last5=Vidmar |first5=John |last6=Butler |first6=Laura |last7=Page |first7=Jonathan E. |last8=Myles |first8=Sean |date=2015-08-26 |editor-last=Tinker |editor-first=Nicholas A. |title=The Genetic Structure of Marijuana and Hemp |url=https://dx.plos.org/10.1371/journal.pone.0133292 |journal=PLOS ONE |language=en |volume=10 |issue=8 |pages=e0133292 |doi=10.1371/journal.pone.0133292 |issn=1932-6203 |pmc=PMC4550350 |pmid=26308334}}</ref><ref>{{Cite journal |last=Roman |first=Madeline G. |last2=Houston |first2=Rachel |date=2020-11 |title=Investigation of chloroplast regions rps16 and clpP for determination of Cannabis sativa crop type and biogeographical origin |url=https://linkinghub.elsevier.com/retrieve/pii/S1344622320300936 |journal=Legal Medicine |language=en |volume=47 |pages=101759 |doi=10.1016/j.legalmed.2020.101759}}</ref><ref name=":10">{{Cite journal |last=Schwabe |first=Anna L. |last2=Hansen |first2=Connor J. |last3=Hyslop |first3=Richard M. |last4=McGlaughlin |first4=Mitchell E. |date=2021-09-29 |title=Comparative Genetic Structure of Cannabis sativa Including Federally Produced, Wild Collected, and Cultivated Samples |url=https://www.frontiersin.org/articles/10.3389/fpls.2021.675770/full |journal=Frontiers in Plant Science |volume=12 |pages=675770 |doi=10.3389/fpls.2021.675770 |issn=1664-462X |pmc=PMC8544287 |pmid=34707624}}</ref> However, in addition to the fact that many crime laboratories are not positioned to integrate these types of analyses into current [[workflow]]s, 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]] [[Resolution (mass spectrometry)|high-resolution]] [[mass spectrometry]] (DART-HRMS)<ref>{{Cite journal |last=Cody |first=Robert B. |last2=Laramée |first2=James A. |last3=Durst |first3=H. Dupont |date=2005-04-01 |title=Versatile New Ion Source for the Analysis of Materials in Open Air under Ambient Conditions |url=https://pubs.acs.org/doi/10.1021/ac050162j |journal=Analytical Chemistry |language=en |volume=77 |issue=8 |pages=2297–2302 |doi=10.1021/ac050162j |issn=0003-2700}}</ref>—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<ref name=":11">{{Cite journal |last=Chambers |first=Megan I. |last2=Musah |first2=Rabi A. |date=2022-03 |title=DART-HRMS as a triage approach for the rapid analysis of cannabinoid-infused edible matrices, personal-care products and Cannabis sativa hemp plant material |url=https://linkinghub.elsevier.com/retrieve/pii/S2468170921000783 |journal=Forensic Chemistry |language=en |volume=27 |pages=100382 |doi=10.1016/j.forc.2021.100382}}</ref><ref>{{Cite journal |last=Rodriguez-Cruz |first=Sandra E. |date=2006-01-15 |title=Rapid analysis of controlled substances using desorption electrospray ionization mass spectrometry |url=https://onlinelibrary.wiley.com/doi/10.1002/rcm.2267 |journal=Rapid Communications in Mass Spectrometry |language=en |volume=20 |issue=1 |pages=53–60 |doi=10.1002/rcm.2267 |issn=0951-4198}}</ref><ref name=":12">{{Cite journal |last=Chambers |first=Megan I. |last2=Musah |first2=Rabi A. |date=2023-05 |title=DART-HRMS triage approach part 2 – Application to the detection of cannabinoids and terpenes in recreational Cannabis products |url=https://linkinghub.elsevier.com/retrieve/pii/S246817092300005X |journal=Forensic Chemistry |language=en |volume=33 |pages=100469 |doi=10.1016/j.forc.2023.100469}}</ref> and ''Cannabis''-derived products (e.g., edibles, personal-care products, vape products, [[Cannabis concentrate|concentrates]]).<ref name=":11" /><ref name=":12" /> 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<ref>{{Cite journal |last=Beyramysoltan |first=Samira |last2=Abdul-Rahman |first2=Nana-Hawwa |last3=Musah |first3=Rabi A. |date=2019-11 |title=Call it a “nightshade”—A hierarchical classification approach to identification of hallucinogenic Solanaceae spp. using DART-HRMS-derived chemical signatures |url=https://linkinghub.elsevier.com/retrieve/pii/S0039914019306290 |journal=Talanta |language=en |volume=204 |pages=739–746 |doi=10.1016/j.talanta.2019.06.010}}</ref> and their headspace chemical signatures.<ref>{{Cite journal |last=Appley |first=Meghan Grace |last2=Beyramysoltan |first2=Samira |last3=Musah |first3=Rabi Ann |date=2019-09-24 |title=Random Forest Processing of Direct Analysis in Real-Time Mass Spectrometric Data Enables Species Identification of Psychoactive Plants from Their Headspace Chemical Signatures |url=https://pubs.acs.org/doi/10.1021/acsomega.9b02145 |journal=ACS Omega |language=en |volume=4 |issue=13 |pages=15636–15644 |doi=10.1021/acsomega.9b02145 |issn=2470-1343 |pmc=PMC6761758 |pmid=31572865}}</ref> A modified version of DART-MS analysis introduced [[Thermal desorption spectroscopy|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).<ref>{{Cite journal |last=Dong |first=Wen |last2=Liang |first2=Jian |last3=Barnett |first3=Isabella |last4=Kline |first4=Paul C. |last5=Altman |first5=Elliot |last6=Zhang |first6=Mengliang |date=2019-12 |title=The classification of Cannabis hemp cultivars by thermal desorption direct analysis in real time mass spectrometry (TD-DART-MS) with chemometrics |url=http://link.springer.com/10.1007/s00216-019-02200-7 |journal=Analytical and Bioanalytical Chemistry |language=en |volume=411 |issue=30 |pages=8133–8142 |doi=10.1007/s00216-019-02200-7 |issn=1618-2642}}</ref> 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.<ref>{{Cite web |last=Pieslak, J.R. |date=2021 |title=Analytical techniques for the differentiation of hemp and marijuana |work=OpenBU |url=https://hdl.handle.net/2144/43518 |publisher=Boston University Libraries}}</ref>

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.


[[File:Fig1 Chambers JofCannRes23 5.png|800px]]
{{clear}}
{| 
 | style="vertical-align:top;" |
{| border="0" cellpadding="5" cellspacing="0" width="800px"
 |-
  | style="background-color:white; padding-left:10px; padding-right:10px;" |<blockquote>'''Fig. 1.''' Workflow for discrimination of hemp and marijuana samples.</blockquote>
 |- 
|}
|}

==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 mass spectrometry|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<ref>{{Cite journal |last=Jolliffe |first=Ian T. |last2=Cadima |first2=Jorge |date=2016-04-13 |title=Principal component analysis: a review and recent developments |url=https://royalsocietypublishing.org/doi/10.1098/rsta.2015.0202 |journal=Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences |language=en |volume=374 |issue=2065 |pages=20150202 |doi=10.1098/rsta.2015.0202 |issn=1364-503X |pmc=PMC4792409 |pmid=26953178}}</ref> and k-means<ref>{{Cite book |date=2010 |editor-last=Sammut |editor-first=Claude |editor2-last=Webb |editor2-first=Geoffrey I. |title=Encyclopedia of Machine Learning |url=http://link.springer.com/10.1007/978-0-387-30164-8 |language=en |publisher=Springer US |place=Boston, MA |doi=10.1007/978-0-387-30164-8 |isbn=978-0-387-30768-8}}</ref><ref>{{Cite journal |last=Lloyd |first=S. |date=1982-03 |title=Least squares quantization in PCM |url=http://ieeexplore.ieee.org/document/1056489/ |journal=IEEE Transactions on Information Theory |language=en |volume=28 |issue=2 |pages=129–137 |doi=10.1109/TIT.1982.1056489 |issn=0018-9448}}</ref> 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)<ref>{{Cite journal |last=Liaw, A.; Wiener, M. |year=2002 |title=Classification and Regression by randomForest |url=https://cogns.northwestern.edu/cbmg/LiawAndWiener2002.pdf |format=PDF |journal=R News |volume=2 |issue=3 |pages=18–22 |issn=1609-3631}}</ref><ref>{{Cite journal |last=Breiman |first=Leo |date=2001 |title=Random Forests |url=http://link.springer.com/10.1023/A:1010933404324 |journal=Machine Learning |volume=45 |issue=1 |pages=5–32 |doi=10.1023/A:1010933404324}}</ref> (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.


[[File:Fig2 Chambers JofCannRes23 5.png|900px]]
{{clear}}
{| 
 | style="vertical-align:top;" |
{| border="0" cellpadding="5" cellspacing="0" width="900px"
 |-
  | style="background-color:white; padding-left:10px; padding-right:10px;" |<blockquote>'''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.</blockquote>
 |- 
|}
|}

''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., ''[[Cannabis indica|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.<ref>{{Cite journal |last=Beyramysoltan |first=Samira |last2=Ventura |first2=Mónica I. |last3=Rosati |first3=Jennifer Y. |last4=Giffen-Lemieux |first4=Justine E. |last5=Musah |first5=Rabi A. |date=2020-04-07 |title=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 |url=https://pubs.acs.org/doi/10.1021/acs.analchem.0c00199 |journal=Analytical Chemistry |language=en |volume=92 |issue=7 |pages=5439–5446 |doi=10.1021/acs.analchem.0c00199 |issn=0003-2700}}</ref> 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 T<sup>2</sup> 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.<ref name=":10" /><ref>{{Cite journal |last=Vergara |first=Daniela |last2=Bidwell |first2=L. Cinnamon |last3=Gaudino |first3=Reggie |last4=Torres |first4=Anthony |last5=Du |first5=Gary |last6=Ruthenburg |first6=Travis C. |last7=deCesare |first7=Kymron |last8=Land |first8=Donald P. |last9=Hutchison |first9=Kent E. |last10=Kane |first10=Nolan C. |date=2017-04-19 |title=Compromised External Validity: Federally Produced Cannabis Does Not Reflect Legal Markets |url=https://www.nature.com/articles/srep46528 |journal=Scientific Reports |language=en |volume=7 |issue=1 |pages=46528 |doi=10.1038/srep46528 |issn=2045-2322 |pmc=PMC5395929 |pmid=28422145}}</ref> 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).


[[File:Fig3 Chambers JofCannRes23 5.png|900px]]
{{clear}}
{| 
 | style="vertical-align:top;" |
{| border="0" cellpadding="5" cellspacing="0" width="900px"
 |-
  | style="background-color:white; padding-left:10px; padding-right:10px;" |<blockquote>'''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.</blockquote>
 |- 
|}
|}

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.


[[File:Fig4 Chambers JofCannRes23 5.png|700px]]
{{clear}}
{| 
 | style="vertical-align:top;" |
{| border="0" cellpadding="5" cellspacing="0" width="700px"
 |-
  | style="background-color:white; padding-left:10px; padding-right:10px;" |<blockquote>'''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.</blockquote>
 |- 
|}
|}

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

{| 
 | style="vertical-align:top;" |
{| class="wikitable" border="1" cellpadding="5" cellspacing="0" width="80%"
 |-
  | colspan="4" style="background-color:white; padding-left:10px; padding-right:10px;" |'''Table 1.''' Confusion matrix associated with the prediction of “out-of-bag” samples in the random forest model.
 |- 
  ! colspan="2" rowspan="2" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Confusion matrix
  ! colspan="2" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Prediction
 |- 
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Hemp
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Marijuana
 |-
  | rowspan="2" style="background-color:white; padding-left:10px; padding-right:10px;" |True
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Hemp (29)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.00
 |-  
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Marijuana (22)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.04
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.96
 |-
|}
|}

{| 
 | style="vertical-align:top;" |
{| class="wikitable" border="1" cellpadding="5" cellspacing="0" width="80%"
 |-
  | colspan="4" style="background-color:white; padding-left:10px; padding-right:10px;" |'''Table 2.''' Performance results of the random forest model for prediction of “out-of-bag” and external validation samples.
 |- 
  ! rowspan="3" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |
  ! colspan="3" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Out-of-bag samples
 |- 
  ! colspan="3" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Accuracy: 0.98 (98%)
 |- 
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Sensitivity
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Specificity
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Precision
 |-
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Hemp (29)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.96
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.97
 |-  
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Marijuana (22)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.96
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
 |-  
  ! rowspan="3" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |
  ! colspan="3" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |External ''C. sativa'' plant materials
 |- 
  ! colspan="3" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Accuracy: 1.00 (100%)
 |- 
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Sensitivity
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Specificity
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Precision
 |- 
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Hemp (12)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
 |-
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Marijuana (11)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |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.

{| 
 | style="vertical-align:top;" |
{| class="wikitable" border="1" cellpadding="5" cellspacing="0" width="80%"
 |-
  | colspan="4" style="background-color:white; padding-left:10px; padding-right:10px;" |'''Table 3.''' Confusion matrix associated with the prediction of external validation samples using a random forest model.
 |- 
  ! colspan="2" rowspan="2" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Confusion matrix
  ! colspan="2" style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Prediction
 |- 
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Hemp
  ! style="background-color:#e2e2e2; padding-left:10px; padding-right:10px;" |Marijuana
 |-
  | rowspan="2" style="background-color:white; padding-left:10px; padding-right:10px;" |True
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Hemp (12)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |1.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.00
 |-  
  | style="background-color:white; padding-left:10px; padding-right:10px;" |Marijuana (11)
  | style="background-color:white; padding-left:10px; padding-right:10px;" |0.00
  | style="background-color:white; padding-left:10px; padding-right:10px;" |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)<ref name=":0" /><ref name=":1" />, with the categorization based upon THC content. Several reports have emphasized the use of GC-FID<ref name=":5" /><ref name=":13">{{Cite journal |last=Zekič |first=Jure |last2=Križman |first2=Mitja |date=2020-12-11 |title=Development of Gas-Chromatographic Method for Simultaneous Determination of Cannabinoids and Terpenes in Hemp |url=https://www.mdpi.com/1420-3049/25/24/5872 |journal=Molecules |language=en |volume=25 |issue=24 |pages=5872 |doi=10.3390/molecules25245872 |issn=1420-3049 |pmc=PMC7763075 |pmid=33322595}}</ref><ref>{{Cite journal |last=Dussy |first=Franz E. |last2=Hamberg |first2=Cornelia |last3=Luginbühl |first3=Marco |last4=Schwerzmann |first4=Thomas |last5=Briellmann |first5=Thomas A. |date=2005-04 |title=Isolation of Δ9-THCA-A from hemp and analytical aspects concerning the determination of Δ9-THC in cannabis products |url=https://linkinghub.elsevier.com/retrieve/pii/S0379073804003408 |journal=Forensic Science International |language=en |volume=149 |issue=1 |pages=3–10 |doi=10.1016/j.forsciint.2004.05.015}}</ref><ref>{{Cite journal |last=Fischedick |first=Justin |last2=Van Der Kooy |first2=Frank |last3=Verpoorte |first3=Robert |date=2010 |title=Cannabinoid Receptor 1 Binding Activity and Quantitative Analysis of Cannabis sativa L. Smoke and Vapor |url=http://www.jstage.jst.go.jp/article/cpb/58/2/58_2_201/_article |journal=Chemical and Pharmaceutical Bulletin |language=en |volume=58 |issue=2 |pages=201–207 |doi=10.1248/cpb.58.201 |issn=0009-2363}}</ref><ref name=":14">{{Cite journal |last=Hazekamp |first=Arno |last2=Simons |first2=Ruud |last3=Peltenburg‐Looman |first3=Anja |last4=Sengers |first4=Melvin |last5=van Zweden |first5=Rianne |last6=Verpoorte |first6=Robert |date=2004-01 |title=Preparative Isolation of Cannabinoids from Cannabis sativa by Centrifugal Partition Chromatography |url=https://www.tandfonline.com/doi/full/10.1081/JLC-200028170 |journal=Journal of Liquid Chromatography & Related Technologies |language=en |volume=27 |issue=15 |pages=2421–2439 |doi=10.1081/JLC-200028170 |issn=1082-6076}}</ref><ref>{{Cite journal |last=Hazekamp |first=A. |last2=Fischedick |first2=J. T. |date=2012-07 |title=Cannabis - from cultivar to chemovar: Towards a better definition of Cannabis potency |url=https://onlinelibrary.wiley.com/doi/10.1002/dta.407 |journal=Drug Testing and Analysis |language=en |volume=4 |issue=7-8 |pages=660–667 |doi=10.1002/dta.407}}</ref> and GC–MS<ref name=":13" /><ref name=":14" /><ref>{{Cite journal |last=Hazekamp |first=Arno |last2=Peltenburg |first2=Anja |last3=Verpoorte |first3=Rob |last4=Giroud |first4=Christian |date=2005-09 |title=Chromatographic and Spectroscopic Data of Cannabinoids from Cannabis sativa L. |url=https://www.tandfonline.com/doi/full/10.1080/10826070500187558 |journal=Journal of Liquid Chromatography & Related Technologies |language=en |volume=28 |issue=15 |pages=2361–2382 |doi=10.1080/10826070500187558 |issn=1082-6076}}</ref><ref>{{Cite journal |last=Namdar |first=Dvory |last2=Mazuz |first2=Moran |last3=Ion |first3=Aurel |last4=Koltai |first4=Hinanit |date=2018-03 |title=Variation in the compositions of cannabinoid and terpenoids in Cannabis sativa derived from inflorescence position along the stem and extraction methods |url=https://linkinghub.elsevier.com/retrieve/pii/S092666901830061X |journal=Industrial Crops and Products |language=en |volume=113 |pages=376–382 |doi=10.1016/j.indcrop.2018.01.060}}</ref><ref>{{Cite journal |last=Namdar |first=Dvory |last2=Charuvi |first2=Dana |last3=Ajjampura |first3=Vinayka |last4=Mazuz |first4=Moran |last5=Ion |first5=Aurel |last6=Kamara |first6=Itzhak |last7=Koltai |first7=Hinanit |date=2019-06 |title=LED lighting affects the composition and biological activity of Cannabis sativa secondary metabolites |url=https://linkinghub.elsevier.com/retrieve/pii/S0926669019301086 |journal=Industrial Crops and Products |language=en |volume=132 |pages=177–185 |doi=10.1016/j.indcrop.2019.02.016}}</ref><ref>{{Cite journal |last=Omar |first=Jone |last2=Olivares |first2=Maitane |last3=Alzaga |first3=Mikel |last4=Etxebarria |first4=Nestor |date=2013-04 |title=Optimisation and characterisation of marihuana extracts obtained by supercritical fluid extraction and focused ultrasound extraction and retention time locking GC-MS: Gas Chromatography |url=https://onlinelibrary.wiley.com/doi/10.1002/jssc.201201103 |journal=Journal of Separation Science |language=en |volume=36 |issue=8 |pages=1397–1404 |doi=10.1002/jssc.201201103}}</ref><ref>{{Cite journal |last=Knight |first=Glenys |last2=Hansen |first2=Sean |last3=Connor |first3=Mark |last4=Poulsen |first4=Helen |last5=McGovern |first5=Catherine |last6=Stacey |first6=Janet |date=2010-10 |title=The results of an experimental indoor hydroponic Cannabis growing study, using the ‘Screen of Green’ (ScrOG) method—Yield, tetrahydrocannabinol (THC) and DNA analysis |url=https://linkinghub.elsevier.com/retrieve/pii/S0379073810001969 |journal=Forensic Science International |language=en |volume=202 |issue=1-3 |pages=36–44 |doi=10.1016/j.forsciint.2010.04.022}}</ref> 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<ref>{{Cite journal |last=Leghissa |first=Allegra |last2=Smuts |first2=Jonathan |last3=Qiu |first3=Changling |last4=Hildenbrand |first4=Zacariah L. |last5=Schug |first5=Kevin A. |date=2018-01 |title=Detection of cannabinoids and cannabinoid metabolites using gas chromatography with vacuum ultraviolet spectroscopy |url=https://onlinelibrary.wiley.com/doi/10.1002/sscp.201700005 |journal=Separation Science Plus |language=en |volume=1 |issue=1 |pages=37–42 |doi=10.1002/sscp.201700005}}</ref>, two-dimensional GC-FID (GCxGC-FID)<ref>{{Cite journal |last=Gröger |first=Th. |last2=Schäffer |first2=M. |last3=Pütz |first3=M. |last4=Ahrens |first4=B. |last5=Drew |first5=K. |last6=Eschner |first6=M. |last7=Zimmermann |first7=R. |date=2008-07 |title=Application of two-dimensional gas chromatography combined with pixel-based chemometric processing for the chemical profiling of illicit drug samples |url=https://linkinghub.elsevier.com/retrieve/pii/S0021967308008297 |journal=Journal of Chromatography A |language=en |volume=1200 |issue=1 |pages=8–16 |doi=10.1016/j.chroma.2008.05.028}}</ref>, and GCxGC-MS with multivariate curve resolution-alternating least squares (MCR-ALS).<ref>{{Cite journal |last=Omar |first=Jone |last2=Olivares |first2=Maitane |last3=Amigo |first3=José Manuel |last4=Etxebarria |first4=Nestor |date=2014-04 |title=Resolution of co-eluting compounds of Cannabis Sativa in comprehensive two-dimensional gas chromatography/mass spectrometry detection with Multivariate Curve Resolution-Alternating Least Squares |url=https://linkinghub.elsevier.com/retrieve/pii/S0039914013010370 |journal=Talanta |language=en |volume=121 |pages=273–280 |doi=10.1016/j.talanta.2013.12.044}}</ref> 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.<ref name=":2" /> However, this report revealed that potential adulterants and inconsistent packaging of samples may have contributed to the observed distinctions.<ref name=":2" /> Another study utilized GC–MS coupled with dispersive pipette extraction (DPX) to investigate forensic casework marijuana and donated hemp samples.<ref name=":3" /> 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.<ref name=":3" /> Another GC-based study sought to differentiate hemp and marijuana through their cannabinoid and terpene profiles using GC-FID and [[principal component analysis]] (PCA).<ref name=":4" /> This study, which included two recreational cultivars and three pharmacy ''Cannabis'' samples, successfully distinguished between the two ''C. sativa'' varieties.<ref name=":4" /> 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.<ref name=":5" /> This study identified several cannabinoids essential for differentiating between ''Cannabis'' strain types<ref name=":5" /> (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.<ref>{{Cite journal |last=Sanchez |first=Lee |last2=Filter |first2=Conor |last3=Baltensperger |first3=David |last4=Kurouski |first4=Dmitry |date=2020 |title=Confirmatory non-invasive and non-destructive differentiation between hemp and cannabis using a hand-held Raman spectrometer |url=http://xlink.rsc.org/?DOI=C9RA08225E |journal=RSC Advances |language=en |volume=10 |issue=6 |pages=3212–3216 |doi=10.1039/C9RA08225E |issn=2046-2069 |pmc=PMC9048763 |pmid=35497720}}</ref> 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<ref>{{Cite journal |last=Chen |first=Zewei |last2=Harrington |first2=Peter de Boves |date=2019-11-19 |title=Pipeline for High-Throughput Modeling of Marijuana and Hemp Extracts |url=https://pubs.acs.org/doi/10.1021/acs.analchem.9b03290 |journal=Analytical Chemistry |language=en |volume=91 |issue=22 |pages=14489–14497 |doi=10.1021/acs.analchem.9b03290 |issn=0003-2700}}</ref>, 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.<ref name=":6" /> 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.<ref name=":6" /> 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.<ref name=":8" /> 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.<ref>{{Cite journal |last=dos Santos |first=Nayara A. |last2=Souza |first2=Lindamara M. |last3=Domingos |first3=Eloilson |last4=França |first4=Hildegardo S. |last5=Lacerda |first5=Valdemar |last6=Beatriz |first6=Adilson |last7=Vaz |first7=Boniek G. |last8=Rodrigues |first8=Rayza R.T. |last9=Carvalho |first9=Verônica V. |last10=Merlo |first10=Bianca B. |last11=Kuster |first11=Ricardo M. |date=2016-08 |title=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 |url=https://linkinghub.elsevier.com/retrieve/pii/S2468170916300297 |journal=Forensic Chemistry |language=en |volume=1 |pages=13–21 |doi=10.1016/j.forc.2016.07.001}}</ref> 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.<ref>{{Cite journal |last=França |first=Hildegardo S. |last2=Acosta |first2=Alexander |last3=Jamal |first3=Adeel |last4=Romao |first4=Wanderson |last5=Mulloor |first5=Jerome |last6=Almirall |first6=Jose R. |date=2020-03 |title=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 |url=https://linkinghub.elsevier.com/retrieve/pii/S2468170919301092 |journal=Forensic Chemistry |language=en |volume=17 |pages=100212 |doi=10.1016/j.forc.2019.100212}}</ref> 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.<ref name=":9" /><ref>{{Cite journal |last=Jacobs |first=Alexander D. |last2=Steiner |first2=Robert R. |date=2014-06 |title=Detection of the Duquenois–Levine chromophore in a marijuana sample |url=https://linkinghub.elsevier.com/retrieve/pii/S0379073814000929 |journal=Forensic Science International |language=en |volume=239 |pages=1–5 |doi=10.1016/j.forsciint.2014.02.031}}</ref><ref>{{Cite journal |last=Watanabe |first=Kazuhito |last2=Honda |first2=Go |last3=Miyagi |first3=Takeaki |last4=Kanai |first4=Masataka |last5=Usami |first5=Noriyuki |last6=Yamaori |first6=Satoshi |last7=Iwamuro |first7=Yoshiaki |last8=Chinaka |first8=Satoshi |last9=Aramaki |first9=Hironori |last10=Yamamoto |first10=Ikuo |date=2017-01 |title=The Duquenois reaction revisited: mass spectrometric estimation of chromophore structures derived from major phytocannabinoids |url=http://link.springer.com/10.1007/s11419-016-0337-6 |journal=Forensic Toxicology |language=en |volume=35 |issue=1 |pages=185–189 |doi=10.1007/s11419-016-0337-6 |issn=1860-8965}}</ref> 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.<ref name=":15">{{Cite web |last=Pingree, C. |date=01 November 2022 |title=H.R.6645 - Hemp Advancement Act of 2022 |work=Congress.gov |url=https://www.congress.gov/bill/117th-congress/house-bill/6645 |publisher=Library of Congress}}</ref> 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.<ref name=":15" /> 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==

*[https://static-content.springer.com/esm/art%3A10.1186%2Fs42238-023-00173-0/MediaObjects/42238_2023_173_MOESM1_ESM.docx Additional file 1] (.docx): Supplementary Mass Spectral Data and Sample Information. (1) Information about ''C. sativa'' plant materials analyzed in this study.

*[https://static-content.springer.com/esm/art%3A10.1186%2Fs42238-023-00173-0/MediaObjects/42238_2023_173_MOESM2_ESM.docx 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
*'''∆<sup>9</sup>-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==
{{Reflist|colwidth=30em}}

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

<!--Place all category tags here-->
[[Category:LIMSwiki journal articles (added in 2023)]]
[[Category:LIMSwiki journal articles (all)]]
[[Category:LIMSwiki journal articles on cannabis research]]
[[Category:LIMSwiki journal articles on cannabis testing]]