Simultaneous quantification of terpenes and cannabinoids by reversed-phase LC-APCI-MS/MS in Cannabis sativa L. samples combined with a subsequent chemometric analysis

Cannabis sativa L. has been the most discussed medicinal plant in recent years. In particular, the dynamic shift from a formerly illicit and tightly controlled substance to a plant recognized for both medicinal and recreational purposes has brought C. sativa into the global spotlight. Due to the ongoing international legalization processes, fast and convenient analytical methods for the quality control of C. sativa flowers for medicinal and recreational purposes are of tremendous interest. In this study, we report the development and validation of a liquid chromatography-tandem mass spectrometry (LC–MS/MS)-based method applying atmospheric pressure chemical ionization (APCI) to fully quantify 16 terpenes and 7 cannabinoids including their acidic forms by a single chromatographic method. The method presented here is unique and simple, as it eliminates the need for derivatization reactions and includes the unconventional analysis of volatile compounds by liquid chromatography. Samples were prepared by a simple and fast ethanolic extraction. Separation was accomplished within 25 min on a reversed-phase C18 column. Method validation was conducted according to international guidelines regarding selectivity, accuracy, precision, robustness, and linearity. Detection was done in multiple reaction monitoring, which allowed the simultaneous quantification of co-eluting analytes applying two selective mass transitions. In addition, due to reproducible in-source decarboxylation, the acidic forms of cannabinoids were reliably quantified using mass transitions of the neutral forms. The accuracy given as the bias was below 15% for all analytes. Matrix effects for cannabinoids were studied by spiking Humulus lupulus extracts with the analytes at varying concentrations. APCI did not show susceptibility toward ion suppression or enhancement. In addition, the recovery effect after spiking was between 80 and 120% for terpenes. Further, 55 authentic C. sativa extracts were fully quantified, and the obtained results for the terpene profiles were compared to state-of-the-art gas chromatography coupled to flame ionization detection. Comparable results were achieved, emphasizing the method’s applicability for cannabinoids and terpenes. Further, acquired metabolite patterns for C. sativa samples were studied, identifying a relationship between cannabinoid and terpene patterns, as well as the abundance of myrcene in CBD-dominant C. sativa strains. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1007/s00216-024-05349-y.


Introduction
In the past century, no other plant has experienced such a dramatic change in public perception as Cannabis sativa L. It has transitioned from being an illicit controlled substance to a promising drug candidate for a wide range of medical problems, such as spasticity in multiple sclerosis, neuropathic pain, anorexia, epilepsy, and chemotherapy-induced nausea [1][2][3][4].Although the recreational use of cannabis with a high THC content is still banned in many countries, the legal framework for its medical use is experiencing a dynamic transformation [4][5][6].Countries like Switzerland have introduced simplified prescription and dispensing regulations of medicinal cannabis by doctors and pharmacies.As a result, they no longer require an exemption permit from the Federal Office of Public Health (BAG) [7].As the demand for medicinal cannabis rises, so does the need for quality assurance standards.In response, the European Directorate for the Quality of Medicines and Health Care (EDQM) has initiated the addition of a new monograph for cannabis flowers in the European Pharmacopeia 11.5 (Ph.Eur.) [8].Although the exact nomenclature and classification of cannabis is still a matter of debate and distinctions between the species and subspecies Cannabis sativa L., ruderalis, and indica have been made and rejected, it is undisputed that cannabis is exceptionally chemical diverse, and classification should be based on its chemical composition [9][10][11].The most widely accepted classification method, which has also been incorporated into the cannabis flower monograph of the Ph.Eur.11.5, is based on subdividing cannabis into different chemotypes based on their tetrahydrocannabinol (THC) and cannabidiol (CBD) content [8].Although a distinction is made between THC dominant, THC and CBD balanced, and CBD dominant variants, other constituents such as cannabigerol (CBG) and minor cannabinoids are gaining interest [12][13][14].Both in vitro and in vivo studies have implied differences in the activity of isolated plant compounds and the whole extract.These observations have given rise to the "entourage effect," a synergistic interplay among phytochemicals, which manifests itself in enhancing and modulating effects through combinations of phytochemicals [13,[15][16][17].With the discovery of the entourage effect, the primary focus on cannabinoids has shifted towards other constituents, such as terpenes.To date, terpenes in cannabis have been mainly valued for their flavor and fragrance profile, but they have their own pharmacological effects.They can interact synergistically with each other and/or with cannabinoids [14,18].Methods to analyze the full composition of metabolites in cannabis have gained momentum.Metabolomics is being used as a tool to capture the full chemical diversity of cannabis and refine its classification into chemotypes [6,[19][20][21].As a result, it is also becoming increasingly evident that commercial labels on cannabis products, which make claims about a pharmacological effect and a strain, do not reliably match the chemical composition of the chemovar [20].The growing interest in the composition of cannabis flowers has also promoted the development of analytical methods.Gas chromatography (GC) hyphenated to flame ionization detection (FID) or mass spectrometry (MS) has been preferred for the analysis of terpenes in general.In contrast, liquid chromatography (LC) coupled with an ultraviolet (UV) detector or MS has been favored for cannabinoid analysis [16,[22][23][24].GC is well suited for analyzing terpenes, as the analytes are both thermostable and volatile.Moreover, GC provides high separation efficiency [9,14].Since terpenes consist of the same repeating isoprene building blocks, they differ only in the number of units and stereochemistry, making separation and mass spectra interpretation challenging [11,25,26].Although GC offers several advantages for terpene analysis, it encounters difficulties in analyzing cannabinoids.The acid forms of cannabinoids undergo partial or complete heat-induced decarboxylation into their neutral form during injection into the GC.To avoid this reaction, various derivatization strategies have been introduced to protect the carboxy groups and improve the volatility of the analytes.However, these derivatization reactions are time-consuming, and achieving complete derivatization is not always possible.For this reason, LC has become the preferred method for cannabinoid analysis [23,25].With the currently growing need for reliable characterized and quantified cannabis products, fast, convenient, and greener analytical methods are crucial.To date, no analytical method combines the analysis of cannabinoids and terpenes in a single step, removing the need for time-consuming derivatization reactions.Efforts have been made to analyze terpenes by LC-MS and atmospheric pressure chemical ionization (APCI).However, the study by Hyland et al. was restricted to seven terpenes and did not include cannabinoids [27].In this study, we developed and validated an LC-APCI-MS/MS method for simultaneously quantifying16 terpenes and 7 cannabinoids.The presented method does not require time and reactant-consuming derivatization and exhibits a shorter run time than most GC methods, leading to a fast and holistic analysis of the major components of cannabis flowers.The presented method takes the complex composition of cannabis flowers into account and allows a simplified chemovar determination.Well-characterized cannabis flowers allow for classification, which in return can be used for clinical research as a standardized and fully characterized phytopharmaceutical.Consequently, this would allow for an in depth understanding of deviating pharmacological properties between products.

LC-APCI-MS/MS method development
The LC-MS/MS system used in this study was a Sciex ExionLC HPLC system coupled to a TripleQuad 3500 MS equipped with an atmospheric pressure chemical ionization (APCI) source (AB Sciex, Redwood City, CA, USA).Chromatographic separation was achieved with a Symme-try® C18 column (4.6 × 100 mm, particle size 3.5 μm) from Waters (Milford, MA, USA) held at 45 °C.A pre-column was installed to prevent premature column degradation (UHPLC polar C18, 2.1 mm ID, SecurityGuard™, ULTRA Cartridges from Phenomenex).Mobile phases A and B consisted of 2 mM ammonium acetate with 0.1% (v/v) formic acid in water and 2 mM ammonium acetate with 0.1% (v/v) formic acid and 5% (v/v) water in methanol, respectively.The injection volume was set to 5 µL.The solvent gradient was as follows: starting conditions were set to 70% B from 0 to 1 min with a flow rate of 0.7 mL/min.The flow was then decreased to 0.6 mL/min, and the gradient was slowly increased to 98% B until minute 20.Conditions were held constant until minute 26.5.Thereafter, the gradient was rapidly decreased to 70% B and held until minute 28.The autosampler was cooled to 15 °C.MS data was acquired in positive mode with the following source conditions: curtain gas (CUR) 20.0 psi, collision gas (CAD) 9 psi, nebulizer current (NC) 5.0 μA, temperature (TEM) 500.0°C, and ion source gas 1 (GS1) 45.0 psi.Entrance potential (EP) was set to 10.0 V.The declustering potential (DP), collision energy (CE), and collision cell exit potential (CXP) were adjusted individually for each analyte, as specified in Table 1.The MRM was conducted as a scheduled experiment with a detection window of 200 s; the target cycle time across single MRM experiments was set to 1 s, which led to a total of 1680 cycles.Data acquisition was performed using Analyst® (AB Sciex, version 1.7.2.), and processing was carried out using Sciex OS (AB Sciex, version 2.0.0.45330).

Method validation
Method validation was conducted on the basis of the International Council for Harmonization of Technical Requirements for Registration of Pharmaceuticals for Human Use (ICH) Q2(R2) guidelines in regard to accuracy, precision, specificity, linearity, range, and robustness [28].Accuracy and precision were determined by preparing six calibration standards (Cals) and three quality control (QC) samples in the low, medium and high concentration range.The limit of quantification (LoQ) was placed at the lowest calibrator.Over the course of a week, the Cals and QCs were measured on three separate days.The accuracy is expressed as the bias.The concentrations of the duplicate QC samples were back-calculated using the calibration and the corresponding fitted regression function.The relative difference was calculated between the mean experimental and theoretical concentrations of the QC samples.The precision was investigated as the relative standard deviation and reported as precision for both intra-and inter-day performance in accordance to Peters et al. [29].The method's specificity was ensured by the presence of both quantifier and qualifier ions and by directly comparing the ion ratios of both MRM transitions for each analyte.In addition, the relative retention time of an analyte was not allowed to further deviate than + / − 0.02 from the internal standard.To rule out matrix effects, three Humulus lupulus samples were spiked at three different concentration levels along the high, medium, and low calibration range to provide an analyte-free matrix for the cannabinoids.Besides being free of cannabinoids, H. lupulus contains a terpene composition similar to that of C. sativa.These samples were analyzed in triplicate, and the recovery effect (RE) was determined according to Formula (1) [24].

Preparation of Cal and QC standards
C. sativa, depending on its chemotype, can be categorized into high THC, balanced THC:CBD, or low THC content varieties.While research and legislation have focused primarily on these two psychoactive compounds, minor compounds are gaining significance.To address this growing interest, two calibration and QC mixtures were prepared to cover a wide range of concentrations.To calibrate for minor compounds, two stock solutions with concentrations of 1000 μg/mL for terpenes and 40 μg/mL for cannabinoids were prepared, respectively.Direct dilutions were created by combining these stock solutions in varying proportions.
(1) On each validation day, 100 μL of the respective Cal or QC was combined with 20 μL of an IS stock solution, resulting in a final IS concentration of 20 μg/mL for propylparaben and 0.18 μg/mL for Δ 9 -THC-D3.Detailed concentrations for Cal and QC samples of minor components are presented in Table 2.For the detection of major compounds, individual stock solutions for Cals and QCs were prepared by direct dilution of each commercial 1000 μg/mL cannabinoid solution.On the corresponding validation day, each Cal and QC stock solution was further diluted 1:1 with an IS stock solution, resulting in a final IS concentration of 0.18 μg/mL for Δ 9 -THC-D3.The final QC and Cal concentrations can be viewed in Table 3.

Sample extraction
Extraction of authentic C. sativa flower samples and processed material was conducted by the Swiss Drug Testing Lab (Winterthur, Switzerland).For each sample, dried C. sativa flowers were ground and dissolved in EtOH to reach a final concentration of 20 mg/mL.The samples underwent a 5-min sonication, followed by centrifugation at 4000 RPM for an additional 15 min.In the final step, extracts were filtered using a 0.45-μm PTFE syringe filter.Filtered samples were stored in a freezer at − 18 °C until analysis.Plants were cultivated in Switzerland.Characteristics of samples can be found in the supplementary Table 1.

Analysis of authentic C. sativa samples
Authentic samples were analyzed at two different dilutions in order to cover both major and minor compounds.For minor compounds, the crude ethanolic extract was combined with an IS stock solution, resulting in a final concentration of 20 μg/mL for propylparaben and 0.18 μg/mL for Δ 9 -THC-D3.In the case of major compounds, the crude extract was initially diluted 50-fold with EtOH and then mixed in a 1:1 ratio with an IS stock solution to achieve a final concentration of 0.18 μg/mL Δ 9 -THC-D3.Between every five samples, a blank of EtOH, and after every 10 measurements, a system suitability test consisting of 200 μg/mL α-pinene, β-pinene, α-humulene, β-caryophyllene, and cis-and trans-citral was injected.Further, QC samples were analyzed after every 20 samples.

Comparison to state-of-the-art methods
In order to compare the efficiency of the LC-MS/MS quantification with that of conventional methods, 15 authentic cannabis samples were additionally analyzed and quantified by a fully validated GC-FID method published by our group previously.A GC Trace 1600 configured with an AI 3000 autosampler (Thermo Fisher Scientific, Waltham, MA, USA) was used for analysis and equipped with a BGB-wax column (30 m × 0.25 mm × 0.25 μm, BGB Analytik, Boeckten, Switzerland).Minor adjustments were made by setting the split ratio to 1:10 [24].As a carrier gas helium 6.0 (PanGas, Dagmersellen, Switzerland) was used.Data acquisition was performed using Chromeleon (Thermo Fisher Scientific, version 7.3.1.).The same minor calibration solutions were used in order to fully quantify authentic samples using GC-FID.

Data visualization and analysis
Chromatograms were generated using Python (version 3.11.6),which included the matplotlib and seaborn libraries.Pearson correlation and violin plots were created using GraphPad Prism (version: 10.1.2).Principal component analysis (PCA) was conducted using MATLAB (version: Release R2022b, version: 9.13.0;The MathWorks Inc.).
Prior to analysis, the data underwent pre-processing by excluding variables with zero entries and auto-scaling.Graphical abstract was created with biorender.com.

Method development
Complete chromatographic separation was achieved for almost all 23 analytes.An oven temperature of 45 °C allowed for optimal separation of these selected analytes (Supplementary Fig. 3).The target scan time was set to 1 s and provided a sufficient number of data points for quantification.Increasing the flow rate did not significantly alter the resolution between peaks (Supplementary Fig. 4).
Trans-citral (2) and carvacrol (3), as well as CBD (12) and CBG (11), overlapped; however, they were distinguishable by mass in MRM. Figure 1 provides a chromatographic overview of the total ion chromatogram (TIC).Varying intensities were observed for all terpenes and cannabinoids, although being present at the same concentration.Two diastereomers, namely menthone (4, 6) and citral (1, 3), were chromatographically separable.Diastereomer ratios can provide valuable information for authenticity and quality control of natural products due to stereoselective synthesis of terpenes in plants [24].Figure 2 provides an overview of the extracted ion chromatogram (XIC) of all analytes.The chromatographic analysis included the acid forms of THC and CBD, as their total content is relevant to legal limits, and the current chemotype classification is based on their ratio.Cannabinoid acids can be detected using APCI and the same mass transitions as their neutral forms due to reproducible decarboxylation in the ionization source.Decarboxylation of the cannabinoid acids in the APCI source was further observed during compound tuning and has been described in the literature [30,31].In positive mode APCI ionization, cannabinoid acids undergo a partial and consistent decarboxylation, as evidenced by the lower peak areas compared to their neutral counterparts.Although APCI in negative mode offers improved ionization efficiency by generating [M-H] + parent ions, the simultaneous analysis of both cannabinoids and terpenes requires a compromise [32].Therefore, the use of neutral transitions for the cannabinoid acids in positive mode APCI allows adequate ionization and quantification of all analytes.Since the acid forms are primarily found in the plant but not in processed cannabis products such as resins and oils, the focus of this study was also on CBDV, CBN, and CBG.Nonetheless, to further prove the broad applicability and flexibility of the method, a randomly selected authentic sample was screened for CBDVA, CBNA, and CBGA using transitions from the neutral forms.All three analytes were found to be present and were baseline-separated from their neutral form (Supplementary Fig. 5).From the calibration, we were able to deduce that the areas for THCA and CBDA at the same concentration were, on average, 11 and 10 times lower than their neutral counterparts, respectively.These consistent ratios between the areas/concentrations of the acid and neutral forms could be used as a correction factor for quantification using the neutral cannabinoids as a calibration standard.However, these correction factors are dependent on ionization source and MS settings and require individual assessment.This approach would require further investigation, and could, however, be interesting from an economic point of view, as the neutral forms are more cost-effective.

Method validation
In order to compensate for possible variations and measurement inaccuracies, we opted to select an internal standard (IS) as a correction factor.Cis-3-hexen-1-ol, propylparaben, and Δ 9 -THC-D3 were considered as potential IS candidates.In our in-house GC-FID method, cis-3-hexen-1-ol has successfully been used as an IS.However, during MRM tuning, it became apparent that very small and, therefore, unspecific fragments (< 50 Da) were selected for cis-3-hexen-1-ol.Thus, it was decided to continue with propylparaben and Δ 9 -THC-D3.Propylparaben proved suitable as an IS for cis-and trans-menthone, as well as citronellol.Δ 9 -THC-D3 served as an IS for all other analytes.In order to control for carry-over during method validation, a blank EtOH sample was injected after every five measurements.No carry-over was observed.Robustness was ensured by testing different oven temperatures and flow rates on critical pairs such as α-humulene and β-caryophyllene.Ultimately, two calibration ranges were chosen because the calibration models for the acids THCA and CBDA did not fully overlap in these ranges, and CBG showed saturation at higher concentrations.Additionally, a different weighting factor for CBD appeared better suited for both the high and low concentration ranges.Two dilution steps were necessary to cover both cannabinoids and terpenes fully.Major cannabinoids required a 100-fold dilution step as they oversaturated the detector when measured undiluted.In return, terpenes were only sufficiently detectable in undiluted samples.To compensate for this discrepancy, various injection volumes and dilutions of authentic samples were evaluated to ensure optimal coverage of all analytes.An injection volume of 5 μL provided symmetric peak shapes, and selecting two different sample preparations allowed an optimal coverage of all analytes.Unlike cannabinoids, terpenes did not ionize as efficiently, which increased their limit of detection.We further investigated whether it was possible to extend the minor calibration range to 60 μg/mL for CBD, CBDV, CBN, and Δ 9 -THC, thus eliminating the need for a second calibration.While feasible, this approach appeared impractical for our research question, as we aimed to retain the possibility of extending this method to resins or oils.In an industrial application, additional considerations would have to be made as to what type of cannabis products are to be analyzed and what the working calibration range is.Method validation took place on 3 separate days and resulted in a bias below 15% for all analytes in the minor and major calibration range.Inter-(RSD R ) and intraday (RSD T ) variability yielded similar results with, the exception of CBDA as a major component at the QC low level, which slightly exceeded 15%.If a value did not deviate by more than 20%, which was the case for all observed analytes, it was considered acceptable.The validation results are detailed in Table 4 for the minor components and Table 5 for the major components.Further, calibration models and details can be viewed in the supplementary materials (Tables 2-3 and Figs.1-2).Considering that terpenes are traditionally analyzed by GC, which has the advantage of high chromatographic resolution and separation performance, guaranteeing selectivity and specificity becomes even more important.Therefore, three criteria were selected to guarantee identification and quantification.To ensure analyte identity and, therefore, specificity, the relative retention time between the analyte and IS had to be within ± 0.02, both qualifier and quantifier ions had to be present, and the ion ratio had to closely match the expected ion ratio.Differences in the ion ratio < 20% were considered acceptable, < 40% marginal, and anything above 40% unacceptable.During validation all analytes met these criteria.Due to the minimal mass differences and structural similarities among terpenes, the resulting fragments alone were unreliable for identification.To exclude further interferences with other cannabinoids and their acids, a mixture consisting of 10 µg/mL THCV, THCVA, CBG, CBGA, CBDV, CBDVA, CBC, CBCA, CBNA, and Δ 8 -THC in acetonitrile was analyzed using the same method.No interferences were observed, and all cannabinoids, including their acids, were baseline separated (Supplementary Figs. 6 and 7).If required, the method presented here can be expanded to contain these 8 cannabinoids including their acids, without the need to acquire additional mass transitions.This further emphasizes the broad applicability of the method.
To assess potential matrix effects such as ion enhancement and suppression, ethanolic H. lupulus extracts were analyzed blank and spiked at three different concentrations in the high, medium, and low range of the minor calibration, and the RE was determined.As the H. lupulus extract is an analyte-free matrix for cannabinoids, it was possible to observe matrix effects in great detail.The results presented in Table 4 were determined for the minor components only, as the matrix is diluted 100-fold with EtOH for the major component determination.It was therefore assumed to have a minimal effect.The RE was in an acceptable range between 80 and 120% for all analytes.Minimal to no matrix effects were observed, which is consistent with the observation, that APCI is less affected by matrix effects compared to electrospray ionization (ESI) [33].

LC-APCI-MS/MS analysis compared to state-of-the art methods
LC is not the method of choice for the analysis of volatile compounds such as terpenes.This becomes even clearer when considering the structural similarities of terpenes and their widespread occurrence.Cannabis, in particular, is a multi-component mixture consisting of numerous cannabinoids, terpenes, and flavonoids.This special feature makes cannabis an interesting medicinal plant but also makes its analytical analysis a challenge.This multitude of chemical diverse classes often requires a combination of different analytical techniques to gain a complete understanding of the composition of cannabis, which can be both time-consuming and expensive [23,[34][35][36].A single method that is able to capture as much compositional variability with minimum sample preparation in flowers of C. sativa would be desirable.Our study showed that terpenes can be efficiently ionized by APCI and that volatility is not a deal-breaker for the analysis with LC.Furthermore, our chromatographic method was able to separate 16 terpenes and at least 15 cannabinoids within 25 min, which is significantly shorter than most GC methods.In addition, the derivatization reactions, a process that can have side reactions and be incomplete, can be omitted.To test the reliability of the method, we compared our LC-APCI-MS/MS method with a state-ofthe-art GC-FID method by randomly selecting and analyzing 15 authentic samples as a subset from the 55 authentic samples for their terpene composition.The same analytes were quantified in the samples by GC-FID, and the results of the two methods were compared using Bland-Altman plots.For reliable results, we only selected analytes that could be quantified with both methods above the LoQ in at least 10 of the samples.Samples in which the analytes were detected but were not within the calibration range were excluded from the method comparison.The calibration models, weights, and range for the analysis using GC-FID can be studied in the supplementary material Table 5.These criteria were only fulfilled by α-and β-pinene, myrcene, linalool, limonene, α-humulene, and β-caryophyllene, terpenes which are known to be major and ubiquitous components in cannabis.With the exception of limonene, all analytes were in good agreement between both methods.Figure 3 depicts the results from the method comparison using a Bland-Altman analysis.Results were considered acceptable if the percentage difference between the two methods was not greater or smaller than 20%.Limonene concentrations determined by LC-APCI-MS/ MS differed by more than double from those determined  by GC-FID.However, the limonene ratio between samples appeared consistent in both the GC-FID and LC-APCI-MS/MS method.As it was initially thought that another substance was co-eluting with limonene, samples were also investigated by GC-FID.Especially α-and γ-terpinene were suspected as co-eluting analytes since both analytes are isomeric forms of limonene, sharing identical masses and fragments.However, neither were detected in any of the samples, and co-elution of these two structural isomers of limonene could be excluded.In general, the concentrations determined by LC were slightly higher than those measured by GC, which may be related to the detection method and split-injection mode.Ultimately, the co-elution of isomeric analytes is far more likely using LC than GC due to lower resolution power.For the other terpenes studied, this seemed to not be a problem.In particular, coelutions of terpenes occurring at extremely low concentrations do not appear to have an influence on the quantification of major terpenes.This method allows faster and more reliable detection of the main constituents of C. sativa than traditional methods.It also has the advantage that the acid forms of the cannabinoids can be analyzed without derivatization.For the observation of extreme minor terpenes, it is still recommended to rely on the high separation power of GC.Since a minimum concentration of 0.05% v/w is required for terpenes to have a pharmacological effect, it is actually sufficient to consider and quantify only the major terpenes when cannabis is used for medicinal purposes [37].We propose that the analysis of minor terpenes can be of particular interest when, for example, the origin or environment of the plant is of interest.

Analysis of authentic C. sativa samples and chemometric analysis
A total of 55 authentic C. sativa flower samples were analyzed and quantified.All samples contained less than 1% total Δ 9 -THC and were therefore within the Swiss legal limit.The total CBD content varied in all samples from 7.5 to more than 30%.CBDV and CBN were detected at low levels in most samples.CBG content varied widely from trace levels to levels of around 1%.According to the Ph.Eur.11.5 chemotype classification, all samples were assigned to the CBD-dominant type [8].The terpene profile was dominated by myrcene, followed by β-caryophyllene, limonene, α-and β-pinene, and finally α-humulene.All other terpenes were either determined at low concentrations, below the LoQ or not detected.The complete profile of all analyzed compounds can be found in Table 5 in the Supplementary information.Figure 4A and B present the violin plots of the terpene profile and the main cannabinoid profile, respectively.Values below the LoQ were excluded from the analysis.As mentioned in the previous section, the limonene content was found to be higher than in the GC-FID measurements.As shown in Fig. 4A, limonene is one of the most important monoterpenes found in C. sativa flowers.This is in line with previous published reports [38].Also, the overall terpene concentrations determined are consistent Fig. 3 Bland-Altman plots in order to compare the concentration determination using LC-APCI-MS/MS with GC-FID of 15 authentic cannabis samples.Blue dotted lines indicate the cut-off for a 10% deviation; red dotted lines a cut-off for 20%.The average concentration [µg/ml] is depicted on the x-axis with the %-difference on the y-axis with the observation that the total terpene content is approximately 3-5% w/w of the flowers and that the CBD-dominant chemotype mostly expresses myrcene [39].The majority of CBD and Δ 9 -THC in the extracts are in their acid forms, CBDA and THCA, respectively (Fig. 4B).This is due to the plant's exclusive synthesis of cannabinoids in their acid forms, with decarboxylation occurring through non-enzymatic processes such as heating.Other processes, such as harvesting, storage, and processing, can also induce decarboxylation.Understanding the difference between acidic and non-acidic forms provides important information about storage and changes in cannabinoid profiles over time, which is particularly important for pharmaceutical use [10].To get a better understanding how cannabinoid levels correlate with terpene levels, a Pearson correlation was performed.The correlation matrix is shown in Fig. 5. Interestingly, the sesquiterpenes α-humulene and β-caryophyllene increased and decreased in concentration with cannabinoid levels.In contrast, monoterpene levels demonstrated an inverse relationship with cannabinoid concentrations, with the exception of CBDA and THCA.This observation could be explained by the fact that the cannabinoid acids are decarboxylated through processes such as heating and processing, resulting in acid-free forms.However, terpenes can also evaporate during such processes, with a more significant effect on monoterpenes than sesquiterpenes.Subsequently, the terpene pattern could also serve as an indicator for the aging and processing of cannabis flowers.In addition, the double bond isomers α-and β-pinene showed the same directional concentration change as the other monoterpenes, with the exception of linalool, which showed an inverse concentration change.Since cannabinoids and terpenes originate from related building blocks (geranyl diphosphate originating from the methylerythritol 4-phosphate pathway), there may be a shortage of precursors for further terpene diversification, leading to an increased cannabinoid synthesis and decreased and selective monoterpene synthesis [40].However, this is only speculation.A principal component analysis (PCA) was performed as an explorative analysis.The biplot, including the loading vectors, is depicted in Fig. 6. 50.4% of the total variability in the data could be explained by the first two principal components.Overall, the analysis reveals the presence of two, possibly three, distinct groups.In the two right quadrants, C. sativa samples appear to cluster into one group.Mostly, all samples of processed material (C.sativa 8, 16, 45, 47) can be found in this region.These samples are characterized by a CBD content of > 25%.It also seems to form mainly due to the contribution of minor cannabinoids such as CBG, CBDV, and CBN.Further, acidfree forms of Δ 9 -THC and CBD appear to correlate with the Notably, correlations between the terpene and cannabinoid patterns are observed.One limitation of this study was the exclusive analysis of CBD-dominant chemotypes; however even within this chemotype, distinct subclasses exist.Given that some of the terpenes are present in pharmacologically relevant concentrations, their analysis should not be overloo ked.

Conclusion
C. sativa is a chemical factory for a diverse array of phytochemicals, ranging from volatile terpenes to psychoactive cannabinoids.Given the current dynamically changing legal situation, the increasing interest in C. sativa as a medicinal plant, and its expanding usage, a thorough understanding of its chemical composition is vital.Moreover, as a potent medicinal plant, it is central investigating the effects of storage, aging, and processing.Only through this deeper understanding can we ensure the production of high-quality phytopharmaceuticals, which is amplified by its recent inclusion into the Ph.Eur.11.5.However, as for the pharmacopeias and legal reason, the main focus has been set on the presence of two cannabinoids.Unfortunately, this view point neglects C. sativa's unique property of being a multicomponent drug.The synthesis of secondary metabolites in C. sativa is not random-there appear to be correlations between terpenes and cannabinoids.A restructuring of the current chemotype classification has been discussed several times in the literature, and a large number of metabolic studies have suggested the inclusion of minor cannabinoids and terpenes in the classification.However, separate analyses of cannabinoids and terpenes are currently preferred.Due to the volatility and structural similarities of terpenes, successful analysis requires high-resolution separation methods.Additionally, cannabinoid acids such as CBDA and THCA are of interest, yet their analysis via GC requires prior derivatization, which can lead to side reactions and incompleteness.To our knowledge, we present the first LC-APCI-MS/ MS method combining the analysis of terpenes, cannabinoids, and cannabinoid acids without the need for prior derivatization reactions.Method validation exhibited reliable results with the bias being below 15% for all analytes.Due to its selected calibration range, the presented method is suitable for analyzing both C. sativa flowers and processed materials.No relevant matrix effects were observed, further supporting the usage of APCI over ESI due to significantly lower susceptibility to ion enhancement and suppression.Methanol, ethanol, and water were used as solvents, making the method greener and sustainable.A drawback, however, is the inability to reliably quantify limonene due to either suspected co-elution with other terpenes or increased sensitivity of the detector.This observation is currently under further investigation.Overall, however, common terpenes and trace terpenes in C. sativa can be separated and quantified with results comparable to those of GC-FID.The method is rapid, with a runtime of under 30 min, and encompasses the two major classes of phytochemicals in C. sativa, allowing chemotype classification beyond THC and CBD using LC.

Fig. 1 Fig. 2
Fig. 1 TIC of 23 analytes in a multi-component mixture for quantitative analysis using LC-APCI-MS/MS

Fig. 4 Fig. 5
Fig. 4 Violin plots showing the concentration of terpenes (A) and major cannabinoids (B) in % w/w of C. sativa flowers.The pharmacologically active threshold of 0.05% w/w for terpenes is indicated by the red dotted line[17]

Fig. 6
Fig. 6 Loading plot after PCA of 55 C. sativa samples.Vectors in blue indicate the directionality and influence of the variables, and red dots represent individual C. sativa samples

Table 1
Simultaneous quantification of terpenes and cannabinoids by reversed-phase LC-APCI-MS/MS…

Table 2
Concentration of Cals and QCs for the analysis of the minor calibration range of C. sativa Cal 1 Cal 2 Cal 3 Cal 4 Cal 5 Cal 6 QC High QC med QC low

Table 3
Concentration of Cals and QCs for the analysis of the major calibration range of C. sativa

Table 4
Results for method validation of terpenes and cannabinoids in the minor concentration range.Inter-(RSD Simultaneous quantification of terpenes and cannabinoids by reversed-phase LC-APCI-MS/MS…

Table 5
Results for method validation for cannabinoids in the major concentration range