Optimization of SPME–GC–MS and characterization of floral scents from Aquilegia japonica and A. amurensis flowers

Background The floral scents of plants play a key role in plant reproduction through the communication between plants and pollinators. Aquilegia as a model species for studying evolution, however, there have been few studies on the floral scents and relationships between floral scents and pollination for Aquilegia taxa. Methods In this study, three types of solid-phase micro-extraction (SPME) fiber coatings (DVB/PDMS, CAR/PDMS, DVB/CAR/PDMS) were evaluated for their performance in extracting volatile organic compounds (VOCs) from flowers of Aquilegia amurensis, which can contribute to the future studies of elucidating the role of floral scents in the pollination process. Results In total, 55 VOCs were identified, and among them, 50, 47 and 45 VOCs were extracted by the DVB/CAR/PDMS fiber, CAR/PDMS fiber and DVB/PDMS fibers, respectively. Only 30 VOCs were detected in A. japonica taxa. Furthermore, the relative contents of 8 VOCs were significant different (VIP > 1 and p < 0.05) between the A. amurensis and A. japonica. Conclusions The results can be applied in new studies of the relationships between the chemical composition of floral scents and the processes of attraction of pollinator. It may provide new ideas for rapid evolution and frequent interspecific hybridization of Aquilegia. Supplementary Information The online version contains supplementary material available at 10.1186/s13065-021-00754-1.


Background
Volatile organic compounds (VOCs), emitted by plant organs such as leaves, flowers and fruits, have served multiple biological functions, including defense against pathogens, parasites, herbivores and interactions with pollinators [1]. Among them, floral aromas are important in the reproductive processes of many plants by attracting pollinators. Traits with a large effect on pollinator preference could play an important role in the evolution of plant reproductive isolation and speciation [2][3][4]. In addition, it has been reported that diversification of the North American clade of Aquilegia (Columbines) was associated mainly with the difference in pollinators [5]. Researchers have studied the relationships between floral morphologies and pollinators. For example, the changes in nectar spur length and flower orientation are highly correlated with the shifts of pollinators from bee to hummingbird to hawkmoth [6]. Moreover, most attempts to classify interactions between insects and flowers have focused on floral odors [7]. For instance, Mimulus lewisii with three monoterpene volatiles can attract bumblebee pollinators, but due to the lack of the above three specific monoterpenes in its sister species M. cardinalis, the pollinator is not bumblebee, so the reproductive isolation between the two sister species can be maintained [4]. Similarly, a single volatile compound (indole) present in flowers of Ipomosis tenuituba but not its sister species I. aggregata, which can attract hawkmoths to flowers [8]. This information says little, however, about the relationships between floral scents and pollination, evolution, and phylogeny of Aquilegia taxa. Until now, approximately 1700 chemical compounds identified in floral scent have been isolated from more than 90 plant families [9]. Among these compounds, the monoterpenes limonene, (E)-β-ocimene, myrcene, linalool, α-and β-pinene, and the benzenoids benzaldehyde, methyl 2-hydroxybenzoate (methyl salicylate), benzyl alcohol, and 2-phenyl ethanol are most common [10].
In our study, using headspace solid-phase micro extraction coupled with gas chromatography-mass spectrometry (SPME-GC-MS), which is common method in the detection of VOCs, the floral scent characteristics of Aquilegia japonica and A. amurensis were evaluated. A. japonica populations are distributed in Northeast of China, North Korea, South Korea and Japan, while A. amurensis populations are restricted to the northern Greater Khingan Mountains of China, Siberia and Mongolia. A. japonica and A. amurensis are sister species, both of two species with different distribution areas are difficult to identify in nature because of their highly similar shape morphology traits. Therefore, Flora of China holds that both of two species are one species [11]. However, the analysis based on genome showed that the differentiation of the two species was obvious (unpublished). Thus, research focusing on the distribution and combination of floral scent compounds at species and subspecies levels may be of the utmost importance for understanding the molecules responsible for attracting pollinators and promoting adaptations and evolutionary processes in angiosperms.
In the analysis of the VOCs, the SPME technique is characterized by its simplicity, speed and sensitivity. It is a convenient sample preparation technique that can be followed by thermal desorption directly in an analytical instrument [12,13]. Recently, several types of SPME fiber coatings have become available for the extraction of analytes, such as nonpolar polydimethylsiloxane (PDMS) fibers, carboxen-polydimethylsiloxane (CAR-PDMS) fibers, polydimethylsiloxane/divinylbenzene (PDMS/ DVB) fibers and divinylbenzene/carboxen/polydimethyl siloxane (DVB/CAR/PDMS) fibers. Furthermore, due to the different compounds that make up the floral scents of different plant taxa, researchers use different types of fiber to study them, for example, Fan et al. [14] used PDMS/DVB fibers for Malus plants; Gao et al. [15] used CAR-PDMS fibers for Freesia × hybrid and Mohammed et al. [16] used DVB/CAR/PDMS fibers for Rose; Silva et al. [17] found that PDMS fiber in melon flowers has poor adsorption for polar compounds. In addition, previous studies have observed molecules with polarity such as protoanemonin, nonanal, dimethoxytoluene, 2-phenyl ethanol and phenyl acetaldehyde in Aquilegia's floral scents [18]. Therefore, in the present study, SPME fibers coated with PDMS/DVB (65 μm), CAR/PDMS (75 μm) and DVB/CAR/PDMS (50/30 μm) were used to identify fibers suitable for measuring the floral scents in Aquilegia. Consequently, our study has not only assessed the performance of different fibers in extracting the VOCs of Aquilegia flowers, but also evaluated the main differences in compounds among the two taxa and provided fundamental information for the scent traits of Aquilegia.

Plant material
The materials were cultivated in a garden from 2017 at Changchun, Jilin, China, including eight individuals of A. amurensis (52.308 N, 124.376 E) and five individuals of A. japonica (41.949 N, 127.925 E) collected from the wild. Fully expanded flowers of the same size, were collected at around between 9 and 10 a.m. In order to reduce the difference between individuals in intraspecies, the flowers collected from different individuals of A. amurensis should be mixed and then approximately 0.6 g was weighted and sealed into 20 mL solid-phase micro extraction (SPME) vials (Agilent Technologies, Germany) immediately for further analysis. A. japonica was also done the same treated. Nine samples of A. amurensis were selected to set three replicates for each fiber and explore the best coating of SPME fiber. Then, six samples of A. japonica were collected and extracted with the best coating of fiber to discriminate different scent intensities of Aquilegia taxa. In addition, an admixture of a certain number of accurately weighted n-alkanes (C7-C30) diluted with hexane (w = 5%) was used as a standard.

Gas chromatography-mass spectrometry experiments
To select an efficient type of fiber coating to extract volatile compounds from the flowers, SPME fibers with three different coatings were used: 65 μm DVB/PDMS (divinylbenzene/polydimethylsiloxane), 75 μm CAR/PDMS (carboxen/polydimethylsiloxane) and 50/30 μm DVB/ CAR/PDMS (divinylbenzene/carboxen/polydimethylsiloxane) (Supelco, Bellefonte, PA, USA). Prior to the analyses, fibers were conditioned for 30 min according to the temperature recommended by the manufacturer. After 10 min equilibration between the flower and the headspace, the SPME fiber was exposed to the headspace of the capped vial to absorb volatile compounds of each sample under heating at 60 °C for 30 min and for 10 min at room temperature. After extraction, the fiber was removed from the flask and immediately inserted into the gas chromatograph injector (GC-MS) for 3 min for thermal desorption at 240 °C.
The flower samples were analyzed and identified using a GC-MS Agilent 7890b gas chromatograph coupled with a 5977b mass spectrometer. Chromatographic separation (GC) was performed using a DB-5MS capillary column (30 m × 0.25 mm × 0.25 μm film thickness, Agilent Technologies, Wilmington, DE, USA). The analytical conditions used were as follows: spitless injection at 240 °C; helium as the carrier gas at a flow rate of 1.0 mL/ min; and GC column temperature program of GC was initially set at 40 °C for 2 min, then heated to 150 °C for 3 °C/min, maintained for 5 min, and finally increased to 250 °C at 20 °C/min and maintained for 8 min. For MS detection, an electron impact (EI) ionization system was used at 70 eV; the temperature of the transfer line and ionization source was 150 and 230 °C, respectively; and full-scan acquisition mode was performed with a mass range of 20-550 Da. Constituents were identified by comparing mass spectra with the National Institute of Standards and Technology (NIST) 14 library (similarity > 75%) and with published data (NIST, http:// webbo ok. nist. gov/ chemi stry/; PubChem, http:// pubch em. ncbi. nlm. nih. gov/). Moreover, the retention time of various compounds in the standard was measured according to the above experimental conditions. According to the retention time of compounds in the floral scents and n-alkanes in the standard, the retention index (RI) was calculated, and compared with the RI in the literature to further determine the components in the floral scents. In addition, relative amounts of compounds were calculated in relation to the total area of the chromatogram by normalizing the peak area (Chemstation B.07.05).

Comparison of compound extraction sensitivity
The extraction sensitivity of various compounds was evaluated by the cumulative area normalization value (CANV) [19]. The CANV is calculated in three steps as follows: In the equations: AVk is the average peak area of compound K measured by the three SPME fibers; Ak(X) is the absolute peak area of compound K extracted by the NAk(X) = Ak(X) AVk; NAn(X).
X SPME fiber, where X is any of the PDMS/DVB, CAR/ PDMS and DVB/CAR/PDMS SPME fibers; NAk(X) is the standardized value of peak area of compound K extracted by the X fiber; and CAk(X) is the cumulative area normalization value of one to more compounds extracted by the X fiber. When the CANV is larger, the sensitivity of the SPME fiber is considered to be higher.

Characterization of VOCs from A. japonica and A. amurensis flowers
One-way analysis of variance (ANOVA) using R software was performed to investigate the significant differences (p < 0.05) in the relative amounts of compounds between the two taxa. The GC-MS dataset was imported to SIMCA-P 14.1 software for statistical analysis. Principal component analysis (PCA) and partial least squares-discriminant analysis (PLS-DA) were used to differentiate the samples and identify marker metabolites. Afterwards, the variable influence on projection (VIP), which summarizes the importance of the X-variables in the PLS-DA model with many components, was used to illustrate the variables that contributed to the separation.

Fiber performance
Three kinds of SPME fibers were used for the SPME-GC/MS full scan analysis of A. amurensis samples. The total ion chromatogram was shown in Fig. 1 and clear ion spectrum was obtained. In our study, three types of fiber coatings (DVB/PDMS, CAR/PDMS, DVB/ CAR/PDMS) were evaluated for their performance in absorbing VOCs, which was determined based on the number of chromatographic peaks that they detected, from flowers of columbines. In total, 55 volatile compounds were identified through MS and RI analysis, belonging to the following different chemical classes: fatty acid derivatives (10), benzenoids (2), monoterpenoids (23) and sesquiterpenoids (20) (Additional file 2: Table S1). Among them, 50 volatile compounds were extracted by DVB/CAR/PDMS fiber, 47 volatile compounds were extracted by CAR/PDMS fiber and 45 volatile compounds were extracted by the DVB/PDMS fiber. The correlation between the three repetitions of each fiber in the detection of compounds was shown in Additional file 3: Table S2. The CANV of DVB/CAR/ PDMS, CAR/PDMS and DVB/PDMS fibers was 85.12, Fig. 1 Total ion chromatogram of VOCs collected from SPME fibers coating with DVB/PDMS (a), CAR/PDMS (b) and DVB/CAR/PDMS (c), respectively. The x axis represents retention time (min) and the y axis represents relative abundance 56.16 and 29.72, respectively. Therefore, the DVB/CAR/ PDMS fiber showed the best efficiency and was used to extract volatile compounds in A. japonica.

Discrimination of the different taxa
The identified compounds and their relative contents (%) in A. japonica flowers were analyzed using DVB/CAR/ PDMS-coated SPME fiber because this type of fiber was more efficient for the extraction of compounds. In order to ensure the accuracy, 6 repetitions were set. A total of 30 volatile compounds were putatively identified in this taxon through MS and RI analysis, including fatty acid derivatives (15), benzenoids (2) and monoterpenoids (13) ( Table 1). The correlation between the six replicates was shown in Additional file 4: Table S3. Furthermore, 12 analytes were not detected in A. amurensis taxa (15.75% of the total content in A. japonica), and 32 volatile compounds were not detected in A. japonica taxa (18.49% of the total content in A. amurensis).
In addition, the main floral scents in A. amurensis were d-limonene and 1R-α-pinene (46.51% and 10.91% of the total content, respectively) while the primary volatile components in A. japonica included d-limonene and 1-hexanol (constituting 55.72% and 9.19% of the total, respectively) ( Table 1). The relative contents of the different chemical classes (fatty acid derivatives, benzenoids, monoterpenoids and sesquiterpenoids) between the two taxa were calculated and compared (Additional file 1: Figure S1). The kinds of terpenes were more in A. amurensis than in A. japonica, and sesquiterpenoids were not detected in A. japonica. However, the kinds of fatty acid derivatives in A. japonica was more than that in another taxon (Additional file 1: Figure S1).

Fiber selection
The choice of the most appropriate fiber is made to cover as many metabolites as possible. To select the most efficient fiber coating for the extraction of VOCs in the Aquilegia taxa, three SPME fibers were used. In our study, the DVB/CAR/PDMS fiber exhibited better extraction efficiency than the DVB/PDMS and CAR/PDMS fibers, presenting the highest CANV (85.12) compared to the other fibers (56.16 and 29.72, respectively). The affinity of the fiber for an analyte depends on the principle of 'like dissolves like' . Previous studies have demonstrated that many polar molecules in the Aquilegia's flora scents [18] as well as the DVB/CAR/PDMS fiber had an intermediate polarity and some studies also confirmed its efficiency [20,21]. The high efficiency may be because the coating with three different components improves the ability to adsorb compounds [22]. The DVB/PDMS fiber is preferred for the extraction of analytes with higher molecular weights (MW 50-300), such as volatiles, amines, and nitroaromatic compounds. Specifically, 8 fatty acid derivatives, 1 benzenoids and 36 terpenes were identified using the DVB/PDMS fiber, fewer than those detected by the other two fiber types. However, the CAR/PDMS fiber is more efficient for the extraction of gasses and low molecular weight compounds (MW 30-225) [23]. Among the 5 compounds that just the CAR/PDMS fiber did not extracted, most have intermediate and higher molecular weights, which was consistent with the results of the Silva et al. [17].
Moreover, these adsorbent type coatings, carried out by sorption of analytes in internal pores, are formed by porous solids. Therefore, saturation of the surface available for adsorption occurs because of the limited thickness. Competition between compounds was more intense when used the CAR/PDMS fiber than used the DVB/CAR/PDMS. When considering repeatability, the CAR/PDMS fiber was better than the 50/30 m DVB/ CAR/PDMS fiber, and Kataoka et al. also reported this result [24]. However, when considering sensitivity, the DVB/CAR/PDMS fiber showed higher performance than that of CAR/PDMS and DVB/PDMS. Thus, the DVB/ CAR/PDMS fiber has been selected for use in the measurement of the floral scents of Aquilegia and six replicate A. japonica flowers were evaluated when using the DVB/ CAR/PDMS fiber.

Scent composition in relation to the pollinators of the two Aquilegia taxa
Speciation in radiating flowering plants is often accompanied by diversification of animal pollinators [24][25][26]. Perhaps the most well-known signal in Aquilegia is floral color, orientation and the structure of spurs [6,27]. Meanwhile, the roles of floral scents have been investigated in other systems [28][29][30], showing that the floral scents are important signals for communication between plants and pollinators, representing an important cue for pollinators [31,32]. Therefore, a prezygotic reproductive barrier is expected when the composition of the floral scent is different. For example, Huber et al. [33] proposed that flowers of two Gymnadenia species with different floral odors, as well as other floral traits such as color and spur length, attracted different pollinators, enhancing prezygotic isolation. The variability of floral scents among entomophilous plants has been reported to depend on the reliance on different pollinator groups with different olfactory preferences [34]. For example, the high relative content of the most volatile monoterpene alkenes (e.g. limonene) in the floral scent of Silene gallica and S. coeli-rosa pollinated by bees has suggested that these compounds are used as attractants of bees [35]. Jürgens and Dötterl investigated VOCs with significant differences are shown in bold RT retention time, ND not detected, RI retention index, VIP variable importance in projection *Represents significant differences between different taxa 0.01 < p < 0.05 **Represents significant differences between different taxa p < 0.01 Our study has identified that the floral scents of the two taxa are dominated by the same one compound (d-limonene), suggesting an adaptation to the same pollinator. Nevertheless, the low-abundance scent components may be effective specific attractants of potential pollinators and cannot be ignored [29]. For instance, the main floral scent compound of the floral four Aquilegia species that Jürgens and Dötterl studied was octanal Fig. 2 The PCA score plots (a) and PLS-DA score plots (b) for datasets of GC-MS from the two taxa (29.5-42%), but the pollinators for these species were varied. The visitation of A. chrysantha was visited by hawk moths may correlate with relatively high amount of 2-phenyl ethanol (13.5%) compared to that of the other three Aquilegia species [29]. Therefore, the fact that the two taxa share the same main floral scent components may be attributed to their closer phylogenetic relationship.
The notable differences between the taxa were the increase in the relative amounts of fatty acid derivatives and the decrease in the relative amounts of monoterpenoids in A. japonica and the detection of various sesquiterpenes only in A. amurensis. Among the fatty acid derivatives, the relative proportions of (Z)-3-hexen-1-ol acetate, (E)-3-hexen-1-ol and methyl decanoate (VIP > 1, p < 0.05) were significantly different between the two species, representing nearly 10% of the total floral scents of A. japonica but not detected in A. amurensis. However, (Z)-3-hexen-1-ol acetate is often released from vegetation rapidly after damage [36]. It can be hypothesized that this compound may have a defense function. The large number of low-abundance sesquiterpenoids in A. amurensis may represent biosynthetic byproducts, as the monoterpenes and sesquiterpenes are derived from the mevalonic acid pathway via farnesyl pyrophosphate [37]. Further experiments are necessary to draw conclusions regarding whether these sesquiterpenes are byproducts or serve critical functions in plant pollinator relationships, further experiments are necessary to draw conclusions.

Conclusions
In this study, by evaluating the properties of different coatings of SPME fibers, the method of extracting and identifying the VOCs of Aquilegia flowers can be optimized. The DVB/CAR/PDMS fiber had the good performance, including sensitivity and repeatability, which is suitable for the subsequent detection of Aquilegia floral scent compounds. In the flowers of two sister species of A. japonica and A. amurensis, except for the main component of VOCs was d-limonene, there were significant differences in the types and relative content of fatty acid derivatives and terpenoid. The types and relative content of fatty acid derivatives in A. japonica were higher than those of A. amurensis, while the types and relative content of monoterpenes were lower than those of A. amurensis, and no sesquiterpenes were detected in A. japonica; there were also significant differences in the contents of eight compounds, including 3-hexen-1-ol, (E)-,3-hexen-1-ol, acetate, (Z)-,methyl decanoate, 1R-α-pinene, (−)-β-pinene, 3-carene, γ-muurolene and α-muurolene. The result provides important information for the future studies involving the VOCs of Aquilegia flowers and can be applied to the new study of relationship between the chemical components of floral scents and the attraction process of pollinators.