High-Throughput Chemotyping of Cannabis and Hemp Extracts Using an Ultraviolet Microplate Reader and Multivariate Classifiers

Profiling alkaloids in Aconitum pendulum N. Busch collected from different elevations of Qinghai province using widely targeted metabolomics

Aconitum pendulum N. Busch (Ranunculaceae) is rich in alkaloids with anti-inflammatory and analgesic activities. Many studies have focused on the identification or quantification of alkaloid components using low-throughput tests. However, the metabolic differences of plants from environmentally distinct regions remain unclear. The present study profiled alkaloid chemical compounds in the rhizomes of A. pendulum from different regions. A total of 80 chemical compounds were identified using a widely targeted ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) approach. Principal component, hierarchical clustering, and orthogonal partial least squares-discriminant analyses of the chemical compounds indicated that the plants from 6 regions clearly separated into distinct groups. A total of 19 compounds contributed the most to the metabolite differences between collection areas and were identified as potential metabolic markers. The anti-inflammatory activities of the A. pendulum extracts were also evaluated and the potential environmental effects on the regulation of metabolite composition and bioactivity were explored. These results improve our understanding of the variation in chemical composition of plants from different regions and will serve as a reference for evaluating the medicinal value of A. pendulum in different environments.

In Situ Determination of Cannabidiol in Hemp Oil by Near-Infrared Spectroscopy

Cannabidiol (CBD, 1) is an active component of hemp oil and many other products that offers diverse health benefits. Near-infrared spectroscopy (NIRS) coupled with chemometrics was utilized to quantify the CBD (1) concentration in the hemp oil through the containing glass vial. NIRS provided a fast and cost-effective tool to measure chemical profiles for the hemp oil samples with various concentrations of CBD (1) and its acid precursor, i.e., cannabidiolic acid (CBDA, 2). The measured NIR spectra were transformed by using a Savitzky-Golay first-derivative filter to remove baseline drift. Two self-optimizing chemometric methods, super partial least-squares regression (sPLSR) and self-optimizing support vector elastic net (SOSVEN), were applied to construct automatically multivariate models that predict the concentrations of CBD (1) and total CBD (sum of 1 and 2 concentrations) of the hemp oil samples. The SOSVEN had validation errors of 6.4 mg/mL for the prediction of CBD (1) concentration and 6.6 mg/mL for the prediction of total CBD concentration, which are significantly lower than the errors given by sPLSR. Other than the lower validation errors, SOSVEN has another advantage over sPLSR in that it builds a multivariate model while selecting spectral features at the same time. These results demonstrated that NIR spectroscopy combined with chemometrics can be used as a rapid and cost-effective approach to determine the CBD (1) and total CBD concentrations in hemp oil. Manufacturers would benefit from the fast and reliable approach in quality assurance.

Self-Optimizing Support Vector Elastic Net

Chemometrics is widely used to solve various quantitative and qualitative problems in analytical chemistry. A self-optimizing chemometrics method facilitates scientists to exploit the advantages of chemometrics. In this report, a parameter-free support vector elastic net that self-optimizes two key regularization constants, i.e., λ for L2 regularization and t for L1 regularization, is developed and referred to as self-optimizing support vector elastic net (SOSVEN). Response surface modeling (RSM) and bootstrapped Latin partitions (BLPs) are incorporated for the optimization. Responses at a set of design points over the ranges of the two factors are evaluated with an internal BLP validation using a calibration set. A 2-dimensional interpolation with a cubic spline fits a response surface to determine the best condition that gives the best-estimated response. The SOSVEN with RSM had comparable performances with the one tuned by grid search, while the RSM is more efficient. The developed SOSVEN was compared with two parameter-free chemometrics methods, super partial least-squares regression (sPLSR) and super support vector regression (sSVR) for calibration, and sPLS-discriminant analysis (sPLS-DA) and support vector classification (SVC) for classification. For calibration, the SOSVEN with RSM worked equivalently well or better than the other two self-optimizing methods for the evaluations using meat and hemp oil data sets. For classification, a reference wine data set and mass spectra of different marijuana extracts were used. The three classifiers had similar performances to identify the cultivars of wines with nearly 98% of accuracy. The SOSVEN significantly outperformed sPLS-DA and SVC to classify the mass spectra of marijuana extracts with an overall accuracy of 97%. These results demonstrated excellent abilities of SOSVEN for classification and calibration.

Automatic soft independent modeling for class analogies

Soft independent modeling of class analogy (SIMCA) is an important method for authentication. The key parameters for SIMCA, the number of principal components and the decision threshold, determine the model's performance. In this report, a self-optimizing SIMCA that automatically determines these two parameters is devised and referred to as automatic SIMCA (aSIMCA). An efficient optimization is obtained by incorporating response surface modeling (RSM) and bootstrapped Latin partitions with the model-building dataset. A set of design points over the ranges of the two parameters are evaluated with respect to sensitivity and specificity by using the model-building data from target and non-target classes. Averages of the sensitivity and specificity are used as responses for the design points. A 2-dimensional interpolation and a bivariate cubic polynomial were used to model the response surface. As a control method, a grid search that evaluates all combinations of the two parameters over the same ranges was performed in parallel to determine the best conditions for SIMCA and the modeling performance was compared to aSIMCA with RSM. The developed aSIMCA methods were evaluated by authenticating two botanical extracts sets, i.e., marijuana and hemp, with spectral datasets collected from various spectroscopic techniques, including nuclear magnetic resonance, high-resolution mass, and ultraviolet spectrometry. Results of a paired t-test indicated that the aSIMCA with the RSM had similar performance with the one optimized by the grid search for modeling marijuana and hemp, while the RSM was more computationally efficient. The 2-dimensional interpolation is preferred because the better efficiency and the fit to the response surface is more precise.

Pipeline for High-Throughput Modeling of Marijuana and Hemp Extracts

Authentication of Cannabis products is important for assuring the quality of manufacturing, with the increasing consumption and regulation. In this report, a two-stage pipeline was developed for high-throughput screening and chemotyping the spectra from two sets of botanical extracts from the Cannabis genus. The first set contains different marijuana samples with higher concentrations of tetrahydrocannabinol (THC). The other set includes samples from hemp, a variety of Cannabis sativa with the THC concentration below 0.3%. The first stage applies the technique of class modeling to determine whether spectra belong to marijuana or hemp and reject novel spectra that may be neither marijuana nor hemp. An automatic soft independent modeling of class analogy (aSIMCA) that self-optimizes the number of principal components and the decision threshold is utilized in the first pipeline process to achieve excellent efficiency and efficacy. Once these spectra are recognized by aSIMCA as marijuana or hemp, they are then routed to the appropriate classifiers in the second stage for chemotyping the spectra, i.e., identifying these spectra into different chemotypes so that the pharmacological properties and cultivars of the spectra can be recognized. Three multivariate classifiers, a fuzzy rule building expert system (FuRES), super partial least-squares-discriminant analysis (sPLS-DA), and support vector machine tree type entropy (SVMtreeH), are employed for chemotyping. The discriminant ability of the pipeline was evaluated with different spectral data sets of these two groups of botanical samples, including proton nuclear magnetic resonance, mass, and ultraviolet spectra. All evaluations gave good results with accuracies greater than 95%, which demonstrated promising application of the pipeline for automated high-throughput screening and chemotyping marijuana and hemp, as well as other botanical products.