Skip to main content

High-throughput interspecies profiling of acidic plant hormones using miniaturised sample processing



Acidic phytohormones are small molecules controlling many physiological functions in plants. A comprehensive picture of their profiles including the active forms, precursors and metabolites provides an important insight into ongoing physiological processes and is essential for many biological studies performed on plants.


A high-throughput sample preparation method for liquid chromatography–tandem mass spectrometry determination of 25 acidic phytohormones classed as auxins, jasmonates, abscisates and salicylic acid was optimised. The method uses a small amount of plant tissue (less than 10 mg fresh weight) and acidic extraction in 1 mol/L formic acid in 10% aqueous methanol followed by miniaturised purification on reverse phase sorbent accommodated in pipette tips organised in a 3D printed 96-place interface, capable of processing 192 samples in one run. The method was evaluated in terms of process efficiency, recovery and matrix effects as well as establishing validation parameters such as accuracy and precision. The applicability of the method in relation to the amounts of sample collected from distantly related plant species was evaluated and the results for phytohormone profiles are discussed in the context of literature reports.


The method developed enables high-throughput profiling of acidic phytohormones with minute amounts of plant material, and it is suitable for large scale interspecies studies.


The acidic phytohormones are structurally diverse signalling molecules regulating plant growth, development and defence. They comprise the phytohormonal classes of auxins (AUXs), jasmonates (JAs), abscisates (ABAs), salicylates (SAs) and gibberellins. AUXs are indolic compounds, precursors and metabolites of bioactive indole-3-acetic acid (IAA), which is crucial to the regulation of almost every aspect of growth and development throughout a plant life [1]. The catabolism of IAA can happen either reversibly through the formation of indole-3-acetyl-1-O-ß-D-glucose (IAA-glc), indole-3-acetyl aspartate (IAA-Asp) and indole-3-acetyl glutamate (IAA-Glu) [2] or irreversibly through the formation of 2-oxindole-3-acetic acid (oxIAA) [3] and, in turn, oxIAA can be further glycosylated to 2-oxindole-3-acetyl-1-O-ß-D-glucose (oxIAA-glc) [4,5,6]. The IAA conjugation and oxidation mechanisms appeared early during plant evolution and are evolutionarily conserved across the whole plant kingdom [7,8,9,10]. JAs are cyclopentenones or cyclopentanones induced by mechanical stresses or responding to wounding or necrotrophic pathogen and herbivore attacks and they contribute to the regulation of plant growth and developmental processes [11]. The pathway for jasmonate biosynthesis begins with the conversion of α-linolenic and roughanic acid to cis-( +)-12-oxo-phytodienoic acid (cis-OPDA) and dinor-oxo-phytodienoic acid (dn-OPDA) respectively. These two cyclopentenone-ring-containing molecules are subsequently subjected to reduction and β-oxidation reactions that lead to the formation of jasmonic acid (JA). JA can be conjugated to amino acids to produce the biologically active form ( +)-7-iso-jasmonyl-L-isoleucine (JA-Ile) [12] or hydroxylated to inactive 12-hydroxyjasmonic acid (12-OHJA) (reviewed in [13]). While jasmonates are widely distributed in seed plants, the presence of JA and its derivatives in lower plants and algae is less understood, as the evidence regarding function of JAs in extant plant species is contradictory (reviewed in [14]). ABAs are isoprenoids derived from abscisic acid (ABA), which is responsive to abiotic stresses such as drought, salinity, heat and cold [15]. There are three different ABA hydroxylation pathways which oxidise one of the methyl groups of the ring structure (C-7′, C-8′, and C-9′), and the hydroxylation of ABA triggers further inactivation steps (reviewed in [16]). Hydroxylation at the C-8′ position is commonly thought to be the predominant ABA catabolic pathway. Consistent with this, the hydroxylated compounds phaseic acid (PA) and dihydrophaseic acid (DPA) are the most widespread and abundant ABA catabolites [16]. The 7′-hydroxyABA (7´-OHABA) form is found in a variety of plant species as a minor catabolite and 9’-hydroxyABA and its isomer neophaseic acid (neoPA) have been identified in several angiosperms [16]. SAs are phenolic compounds that have emerged as key plant defence molecules with critical roles in different aspects of plant immunity (reviewed in [17]). Despite the importance of salicylic acid (SA) in plant defence, its metabolic and signalling pathways are not fully understood. This is mainly due to the fact that these phenolic compounds have been relegated for a long time to the category of ‘secondary metabolites’ as they were considered non-essential for critical processes [18, 19]. Gibberellins are diterpenoids involved in plant developmental processes such as organ growth, seed germination, maturation and flowering [20]. Due to sensitivity issues caused by extremely low endogenous concentrations in plants and lower ionization efficiency of their molecules using mass spectrometry (MS) detection [21,22,23], they were omitted from our method.

The endogenous levels of the signalling forms of hormones and hence their activities in plants are tightly modulated by a combination of biosynthesis, metabolism, transport or release from their conjugated forms [16, 24,25,26]. Data on their concentrations reflect processes occurring in plants and offer an important insight into plant physiology. Depending on the plant tissue or the plant species itself, phytohormones usually occur at very low concentrations. The analytical methods employed in the determination of plant hormone levels therefore need to be very sensitive and capable of covering wide ranges of concentrations in order to track both low basal levels and levels that are elevated in response to stimuli. Moreover, the plant matrix is quite challenging from the analytical point of view. Nowadays plant hormone quantification is mostly performed by means of liquid chromatography–tandem mass spectrometry (LC–MS/MS) systems with electrospray ionisation (ESI). The current trend is to analyse large numbers of phytohormones in decreasing amounts of plant material using simple, high-throughput sample preparation techniques to provide comprehensive information about phytohormonal levels with sufficient statistical significance.

Due to the different chemical natures of targeted analytes, sample preparation from plant material prior to multi-hormone LC–MS/MS analysis is based on one of two strategies: simple, more generic, or laborious, comprising several selective steps. The simplest way is to only employ an extraction step [27,28,29,30]. This approach usually requires larger amounts of plant material, hundreds of mg fresh weight (FW) or tens of mg dry weight (DW), reaching compromise between extraction efficiency and ion suppression caused by interfering matrix present in the extracted sample. Another approach is extraction combined with derivatisation, which selectively enhances ionisation of the targeted analytes and lowers the limits of MS detection in the presence of plant matrix. For this approach very small amounts (units of mg and lower) of plant tissue are sufficient [31]. The most common form of sample preparation is extraction followed by solid phase extraction (SPE) purification using ion exchange or reverse phase sorbents [32,33,34], which can be also carried out in a 96-well set up [35]. SPE purification reduces the amount of interfering matrix in the sample and enables pre-concentration of the analytes, and it usually requires tens of mg FW plant tissue.

Small-scale purification based on classical SPE takes advantage of sensitive MS/MS detection and the ability to pre-concentrate the extracted analytes, simultaneously reducing the consumption of sample, sorbents and solvents [36, 37]. Miniaturisation also facilitates automation, increases throughput in sample preparation and decreases the economic and environmental impact. The in-tip microSPE approach using small plates of Empore discs placed in ordinary pipette tips, originally developed for the purification of proteins [38], has been applied to the purification of phytohormone classes such as cytokinins [39] and AUXs [40] in units of mg FW Arabidopsis thaliana plant tissue.

The aim of this work was to develop and validate a simple, miniaturised, high-throughput method for acidic phytohormone profiling from small amounts of plant samples that would be widely applicable in large-scale interspecies studies. Here we use this approach to analyse different amounts of nine matrices from plant species ranging from algae to land plants. The selected plant species are models often used in plant research. Moreover, they were chosen with respect to the representation of various groups within the phylogenetic tree. Determination of endogenous phytohormone profiles from evolutionarily related landmark species provides further substantial tiles to identify patterns that may have contributed to the diversification of extant plant lineages.

Results and discussion

LC–MS/MS optimisation

LC–MS analysis of plant hormones predominantly employs reverse stationary phases, particularly columns with C18 chemistry [28, 29, 31, 32, 34, 35]. To retain acidic phytohormones, taking into account the requirements for MS detection and the type of stationary phase, the mobile phases typically consist of volatile additives (formic or acetic acid) in water and acidified or pure methanol or acetonitrile. In this study, to optimise separation conditions two reverse phase C18 columns were tested: UPLC CSH C18 (2.1 ×100 mm, 1.7 μm, Waters) and Kinetex Evo C18 (2.1 ×150 mm, 2.6 μm, Phenomenex). The analysis performed on the CSH column (1.7 μm particle size) under conditions using 10 mmol/L formic acid and acetonitrile as mobile phases [33] generated backpressures close to the limits of the 1290 Infinity LC system (Agilent), and the elution of polar analytes was concentrated at the very beginning of the linear gradient. Using a Kinetex Evo column with a larger particle size (2.6 μm) lowered the system backpressure, and methanol instead of acetonitrile increased retention of the most polar compounds (e.g. oxIAA-glc), improved the overall chromatographic separation of the analytes and shortened the analytical run to 20 min (Table 1, Fig. 1).

Table 1 LC–MS/MS method parameters
Fig. 1
figure 1

LC–MS/MS analysis of 25 acidic phytohormones, their precursors and metabolites

Using the optimised chromatographic method was not possible to distinguish between 11-OHJA and 12-OHJA, because both compounds have the same retention times and fragmentation pattern. As in [33], the sum of unresolved 11-OHJA and 12-OHJA has been used in this study. The 11-OHJA and 12-OHJA can be separated by a gas chromatography MS method [41] or by using LC–MS with isocratic elution (Acquity UPLC BEH C18 column, 92% of 0.1% formic acid in water and 8% of 0.1% formic acid acetonitrile) [42].

High-throughput protocol optimisation

An in-tip microSPE stagetip assembled from SDB-XC (poly(styrenedivinylbenzene)) and C18 sorbents (Fig. 2) was used to optimise sample preparation conditions. Type and amount of the sorbent were set as previously published for AUXs [40]. Both microSPEs made up of SDB-XC/C18 and SPE Oasis® HLB (poly(divinylbenzene-co–N-vinylpyrrolidone) copolymer) sorbents provide reverse phase retention of analytes. The predicted pKa values for acidic phytohormones range from 2.8 to 5.1 (Table 1). In order to improve their retention on reverse phase sorbents, analytes need to be in an uncharged form at a pH below their pKa, which in this case requires acidic conditions. This corresponds to the following SPE purification methods for individual acidic phytohormone classes. For auxin analysis, plant extracts in sodium phosphate buffer (50 mmol/L, pH 7) were adjusted to pH 2.7 with hydrochloric acid before loading on an Oasis® HLB SPE column [43]. A similar approach was chosen for auxin metabolite profiling using a microSPE made from SDB-XC and C18 sorbents [40]. Purification of stress hormones (JAs, ABA, SA and IAA) on an Oasis® HLB SPE column included a preconditioning step with 0.1% formic acid in water before loading plant extracts in 10% aqueous methanol [33]. ABA and its metabolites were extracted in 10% aqueous methanol containing 1% acetic acid before being purified on an Oasis® HLB SPE column [44].

Fig. 2
figure 2

High-throughput sample preparation workflow

The impact of sample acidification on the retention of standards on microSPE sorbent (SDB-XC/C18) was tested. Authentic standards (1 pmol) in 10% aqueous methanol and in 1 mol/L formic acid in 10% aqueous methanol were loaded on microSPE tips and processed as described in Fig. 2. The mean purification recovery of all phytohormones was 50% for standards loaded in 10% aqueous methanol and 64% for standards in 1 mol/L formic acid in 10% aqueous methanol. Acidification was beneficial for the retention and subsequent elution of all phytohormones except for cis-OPDA and OPC-8 (Additional file 1: Fig. S1). Next, the concentration of acid in the extraction solvent and the method of acidification were investigated. Process efficiency (PE) was assessed using 2 mg FW of Arabidopsis thaliana 10-day-old seedling matrix spiked with internal standard (IS; 1 pmol of all AUXs and JAs detected in ESI + ; 2 pmol for JAs detected in ESI-, ABAs and SA) extracted in 10% aqueous methanol containing different concentrations of formic acid (0.1–0.25–0.5–1 mol/L) and 10% aqueous methanol acidified to pH 2.7 (HCl) before purification. The mean PEs for extracts in 0.1–0.25–0.5–1 mol/L formic acid in 10% aqueous methanol and 10% aqueous methanol acidified to pH 2.7 (HCl) before purification were 55%, 59%, 58%, 57% and 60%, respectively. The results thus show no significant trend regardless of formic acid concentration or whether the sample was extracted in acidified solution or acidified before purification (Fig. 3).

Fig. 3
figure 3

The impact of acidification on process efficiency (PE). PE (mean area of IS spike before extraction and purification / mean area of IS spike without extraction and purification injected in the same amount × 100) was estimated using 2 mg FW of 10-day-old Arabidopsis thaliana seedlings. Error bars represent ± SD, n = 4

A time-dependent increase in estimated cis-OPDA levels was previously reported in ground Arabidopsis thaliana plant material kept on ice [27]. This phenomenon was explained by possible residual activity of lipases releasing cis-OPDA from galactolipids—Arabidopsides [45]. Likewise to cis-OPDA, dn-OPDA is a constituent of the Arabidopsides [46]. When the influence of extraction approach on measured endogenous levels of phytohormones in Arabidopsis thaliana 10-day-old seedlings was studied (acidic extraction in 1 mol/L formic acid in 10% aqueous methanol and extraction in 10% aqueous methanol with pH adjusted before purification), endogenous cis-OPDA and dn-OPDA levels were found to be hugely dependent on the type of extraction (Additional file 1: Fig. S2). From these results we concluded that extraction in 1 mol/L formic acid can prevent potential enzyme activity and provide stable extraction conditions. On the other hand, the acidic conditions may decrease the stability of some analytes, especially JA or IAA conjugates with amino acids and glucose [34]. Regardless of stability, acidic conditions are crucial for the retention of the majority of analytes on reverse phase sorbent. The impact of analyte stability on precise quantification is compensated by appropriate internal standards, which are added to the plant material before extraction and get exposed to acidic conditions in the same way as the endogenous analytes.

To elute acidic plant hormones retained on the SPE reverse phase sorbent, a high percentage of organic solvent (e.g. 80% aqueous methanol) is utilised [33, 40, 43, 44]. In reverse phase chromatography of acidic phytohormones the gradients usually start at a low organic content. In order to achieve high quality of chromatographic peak shapes while injecting larger sample volumes (e.g. 10 μL), it is beneficial to include an evaporation step, which removes the excess organic solvent, and reconstitute the samples in a solvent with a low organic content similar to the composition of mobile phases at the beginning of the gradient [47]. This step also provides an opportunity to pre-concentrate the analytes. The optimum means of sample evaporation after elution (50 μL of 80% aqueous methanol) was investigated. Direct lyophilisation of 96-well plates was the most high-throughput approach available. The best recoveries were obtained using evaporation in vacuo (94%), followed by lyophilisation (93%); evaporation under a stream of nitrogen (89%) gave the lowest recovery (Additional file 1: Fig. S3). The most problematic compound was MeJA, in all probability due to its volatility [48].

Method validation

The optimised high-throughput method was then validated. With the exception of OPC-8, the accuracies fell in the range 85–115% (Table 2). Despite the very similar structure and physico-chemical properties of OPC-8 and cis-OPDA, which differ only in the presence of a double bound on the cyclopentane ring (cis-OPDA), the internal standard, [2H5]-OPDA, was not able to fully compensate for losses during sample preparation, as can be seen from the different recoveries of these compounds (OPC-8 23% compared to cis-OPDA 43%) (Table 2). Precisions were lower than 15% for all analytes (Table 2).

Table 2 Validation parameters

Matrix effects (ME) express the signal suppression (ME below 100%) or enhancement (ME above 100%) caused by constituents of the sample matrix co-eluting with analytes; this occurs frequently when using MS/MS detection [49]. The ME for 2 mg FW of 10-day-old Arabidopsis thaliana seedlings ranged from 40 to 114%. The most pronounced ME was for MeJA (40%). Rather than ion suppression, the result is a consequence of the known volatility of MeJA [48] manifested during sample evaporation/lyophilisation, as the calculation of ME and also PE is related to neat standard solution, which was not subjected to sample preparation including evaporation (lyophilisation).

The recovery (RE) represents the losses of analyte in the course of sample preparation processes [49]. In our experiment the recoveries ranged from 22 to 105%. The process efficiency (PE) reflects both the influence of sample loss during sample preparation (RE) and ME [49]. The PE values were 27–107%. The lowest RE and PE results were obtained for oxIAA-glc (22 and 28% respectively). The low stability of oxIAA-glc contributed to these results, since spiking before extraction/purification resulted in longer exposure to acidic conditions compared to spiking after extraction/purification or neat standard solution.

Phytohormone profiles of plant species representatives

The optimised high-throughput sample preparation protocol was applied to profile phytohormones in nine matrices of different plant species across the Plantae (Fig. 4). To evaluate the influence of sample amount on the ability of the method to quantify the acidic phytohormones, different amounts of these matrices (1–2–4–6–8 mg of FW and 0.1–0.2–0.4–0.6–0.8 mg of DW, Stichococcus bacillaris) were processed. The endogenous levels and the minimal sample weight necessary for detection and quantification of acidic phytohormones are listed in Table 3. Importantly, the sample weight and the levels of phytohormones found in that sample showed a linear relationship expressed as R2 ranging from 0.9820 to 0.9999 (Additional file 1: Table S1). Finally, we investigated how AUXs, JAs and ABAs contributed to the total pool in several distantly related species (Fig. 4).

Fig. 4
figure 4

Schematic illustration of endogenous contents of AUXs, JAs and ABAs in selected species of land plants and algae assigned to groups. Left panel—Proportion of IAA, amide-linked IAA conjugates (IAA-Asp and IAA-Glu), IAA-glc, oxIAA and oxIAA-glc expressed as percentage of the total auxin pool. Middle panel–Proportion of JA, sum of 11-OHJA and 12-OHJA, 9,10-dhJA, JA-Ile, cis-OPDA and dn-OPDA expressed as percentage of the total jasmonate pool. Right panel–Proportion of ABA, PA, DPA, neoPA and 7´-OHABA expressed as percentage of the total abscisate pool. Results reported in Table 3 were used to draw this chart, when more species per group, the proportion was constructed based on the mean of results

Table 3 Phytohormone levels in different plant species

Free IAA was quantified in all the plant species analysed (Table 3). While free IAA was the only auxin detected in the chlorophyte S. bacillaris, in the monocot T. aestivum and in the two asterids, N. tabacum and S. lycopersicum, oxIAA followed by the glucosylesters of IAA and oxIAA were the most abundant auxin forms in the moss Physcomitrium patens, confirming previous findings reported for microalgae and bryophytes [50,51,52]. Among the seed plants, the proportion of ester conjugates in the total auxin pool was high in the conifer P. abies and in the two Brassicaceae species, A. thaliana and B. rapa. IAA-glc was the major glucosylester form in gymnosperms, while oxIAA-glc represented the majority of the total auxin pool in Brassicaceae, in accordance with published data [3, 10, 53]. The two amide conjugates of IAA, IAA-Asp and IAA-Glu, made up the majority of the total auxin pool in the woody eudicot hybrid poplar, confirming a previous finding reported in poplar where the IAA conjugate concentration was found to be higher than that of free IAA [54]. Interestingly, oxIAA was not detected in hybrid poplar leaves, while oxIAA was previously found to contribute to the asymmetric IAA distribution in bent roots of black poplar [55]. These differences in the auxin pool in poplar species might reflect the different tissues analyzed here (leaves, this study; roots, [55]) and/or species-specific metabolic pathways (hybrid poplar, this study; black poplar, [55]). Our results seem to corroborate the finding that universal elements exist among land plants that appear to form a metabolome in order to efficiently modulate auxin levels via conjugation and/or oxidation.

In jasmonate profiles, JA was detected in all the species analysed, with the exception of the moss P. patens and the chlorophyte S. bacillaris. Interestingly, the dihydro-oxylipin 9,10-dhJA and JA-Ile were the only JAs detected in the chlorophyte S. bacillaris, and to our knowledge this is the first report of a study where jasmonates were analyzed in this species. 9,10-dhJA was also found in fungi and plants, but descriptions of the presence of this dihydro-pentanone in plants are still limited. However, it has been reported that cis-OPDA was detected in the charophyte Klebsormidium flaccidum, but not in the charophyte Chara braunii (reviewed in [14]). cis-OPDA was the only jasmonate detected in the moss P. patens and this finding is in agreement with previous studies in the two bryophytes P. patens and Marchantia polymorpha, where JA and JA-Ile were not detected [56,57,58], whereas these metabolites could be found in a wide range of other bryophytes [51]. Literature on the distribution of JA and its derivatives in the lower plants, bryophytes and algae, appears to be contradictory and the presence or absence of jasmonates in bryophytes and algae requires further studies for confirmation [14]. Nonetheless, currently available evidence suggests that cis-OPDA is not only the precursor of JA biosynthesis but has distinct physiological roles from those of JA-Ile, acting through a JA-independent signalling pathway, and that its biosynthesis originated in the algal lineage, before the emergence of land plants [14, 25, 59, 60]. For gymnosperms, cis-OPDA was found to be the most abundant jasmonate in the conifer P. abies, followed by JA and JA-Ile (Fig. 4) and this is in agreement with findings described in [61]. In the grass T. aestivum, cis-OPDA and, to a lesser extent, 12-OHJA, were the two most abundant jasmonates detected (Fig. 4), while JA, JA-Ile and 9,10-dhJA were found only at very low concentrations (Table 3). With respect to the dicotyledonous species, JA and JA-Ile similarly contributed to the total jasmonate pool, whereas apparent dissimilarities were found for the other jasmonates profiled (Fig. 4). For instance, dn-OPDA was detected only in Arabidopsis samples and, interestingly, it was at a similar concentration to that of cis-OPDA (Table 3). The dihydro-pentanone 9,10-dhJA was found only in the two asterids and B. rapa, not in A. thaliana or hybrid poplar (Table 3).

ABA was detected in all the species analyzed. ABA, followed by neoPA, contributed the most to the total abscisate pool in the chlorophyte S. bacillaris, confirming the presence of ABA in this species (reviewed in [62]). In the bryophytes, ABA was the only abscisate found in the moss P. patens. However, various ABA metabolites were detected in several other bryophytes, as reported in [51]. The conifer P. abies seems to prefer to catabolise ABA through the formation of neoPA, as this metabolite was found to be more abundant than DPA and 7´-OHABA. This is in contrast with findings reported in [63] where the main hydroxylated ABA catabolite was PA, followed by 9′-hydroxyABA and DPA, whereas neoPA and 7′-OHABA were below detection limits. These differences might be due to the age of the plants analyzed (2 week-old, this study; 6 week-old, [63]) and/or on the growth conditions used (vermiculite-grown seedlings, this study; hydroponic-grown seedlings, [63]). Among the angiosperms, the grass and the asterids appear to preferentially catabolise ABA through the formation of PA, while hybrid poplar catabolises ABA through DPA formation (Fig. 4). ABA was the only abscisate detected in A. thaliana, whereas DPA was found to be the most abundant ABA catabolite, followed by PA, 7´-OHABA and neoPA in B. rapa (Table 3). It should be noted that the catabolic dynamics that regulate the ABA content in plants can be subject to significant changes upon stresses that trigger the ABA response.


The acidic plant hormones are important regulators of plant processes, and comprehensive data on their profiles is often fundamental for plant biology research. A simple method for acidic plant hormone profiling comprising acidic extraction, high-throughput reverse phase in-tip microSPE sample purification using a 3D printed device and sensitive LC–MS/MS analysis was optimised and validated. The method allows simultaneous screening of 25 phytohormones and their precursors and metabolites and makes it possible to track their presence and the dynamics of their turnover in various species, requiring only units of mg of plant samples. Our protocol currently represents a fast and simple approach to phytohormone profiling that facilitates the processing of large sets of diverse plant samples.



Authentic and stable isotope labelled standards (Table 1) were provided from in-house standard library of JAs: cis-OPDA, OPC-8, OPC-6, OPC-4, dn-OPDA, JA, MeJA, 9,10-dhJA, 12-OHJA, JA-Val, JA-Ile, JA-Phe and [2H5]-OPDA, [2H6]-JA, [2H2]-( − )-JA-Ile were purchased from OlChemIm Ltd. (Olomouc, Czech Republic), JA-Trp was synthetised as described in [64], 11-OHJA was kindly provided by Dr. Otto Miersch, [2H6]-( ±)-MeJA was synthesized by acetyl chloride–methanol-d4 esterification; AUXs: IAA was purchased from Sigma Aldrich (St. Louis, MO, USA), oxIAA, IAA-Asp, IAA-Glu, IAA-glc, oxIAA-glc and [13C6]-IAA, [13C6]-IAA-Asp, [13C6]-IAA-Glu were purchased from OlChemIm Ltd. (Olomouc, Czech Republic), [13C6]-oxIAA was synthesized as described in [65], [13C6]-IAA-glc, [13C6]-oxIAA-glc were synthesized as described in [66, 67]; ABAs: ABA, PA, DPA, neoPA, 7′-OHABA and [2H6]-ABA, [2H3]-PA, [2H3]-DPA, [2H3]-neoPA, [2H4]-7′-OHABA were purchased from National Research Council Canada (Saskatoon, Canada); SAs: SA and [2H4]-SA were purchased from Sigma Aldrich (St. Louis, MO, USA). Methanol, acetonitrile, formic acid, all LC–MS grade, acetic acid of gradient grade and hydrochloric acid were purchased from Sigma-Aldrich (St. Luis, MO, USA). Purified Milli-Q water from a Simplicity 185 System was used for preparation of aqueous solutions.

Plant material and growth conditions

Arabidopsis thaliana — seedlings of 10-day-old Arabidopsis thaliana ecotype ‘Col-0’ grown in vitro in Murashige and Skoog medium Petri dishes in a growth chamber at 20 °C, under long-day photoperiod (16 h light, 8 h dark) with a light intensity of 90 µmol photons m−2 s−1, were used for the optimisation of the extraction and purification protocol, method validation and determination of phytohormones in a series of matrices from different specimens.

Populus tremula x alba – the hybrid poplar seedlings were propagated vegetatively and grown in vitro in jars containing 20 mL of media (2.183 g/L Murashige and Skoog 1A media, 1 g/L glucose, and 8 g/L plant agar enriched with 112 mg/L Gamborg’s B5 vitamin mixture, L-Glutamine 20 mg/L, pantothenate 0.1 mg/L, L-cysteine chlorhydrate 0.1 mg/L and biotin 0.001 mg/L) under long-day conditions (16 h light at 22 °C, 8 h dark at 18 °C). The leaves from 90-day-old seedlings were cut at petiole and used in this study for phytohormone quantification.

Brassica rapa — leaves from 4-week-old Brassica rapa (cultivar ‘Pekinensis’) seedlings, germinated as described in [68] and cultivated in soil with vermiculite (2:1, v/v, respectively) in a growth chamber under long-day conditions at 21 °C, were used for phytohormone quantification.

Solanum lycopersicum — leaves from 8-week-old tomato seedlings (cultivar ‘Tornado’) cultivated in soil under environmentally controlled greenhouse conditions (natural photoperiod) during spring 2020 (from the middle of March to the middle of May) (Olomouc, Czech Republic), were used for phytohormone quantification.

Nicotiana tabacum — leaves from 8-week-old tobacco seedlings, cultivated in soil in a growth chamber under long-day conditions (16 h light at 20 °C, 8 h dark at 18 °C) were used for phytohormone quantification.

Triticum aestivum — leaves from 7-day-old seedlings of winter-type wheat (cultivar ‘Turandot’) were used in this study for phytohormone quantification. Seedlings were obtained from dormant seeds that were soaked for 5 h in tap water to germinate, sown into wet vermiculite, cultivated under long-day conditions at 21 °C, and irrigated every second day with 0.5 L of Hoagland’s nutrient solution [69].

Picea abies — 14-day-old spruce whole seedlings, grown in a growth chamber under long-day conditions (16 h light at 22 °C, 8 h dark at 18 °C) as described in [10], were used in this study for phytohormone quantification.

Physcomitrium patens — 3-week-old whole gametophores from Physcomitrium patens (Hedw.) Bruch & Schimp (ecotype ‘Gransden 2004′) wild type were used in this study for phytohormone quantification. The plants were grown axenically in 9-cm Petri dishes on Knop medium (100 mg/L Ca(NO3)2·4H2O, 25 mg/L KCl, 25 mg/L KH2PO4, 25 mg/L MgSO4·7H2O and 1.25 mg/L FeSO4·7H2O, pH 5.8). Knop medium was supplemented with 92.05 mg/L ammonium tartrate, 0.5 mg/L nicotinic acid, 0.125 mg/L p-amino benzoic acid, 2.5 mg/L thiamine HCl, trace-element solution (0.614 mg/L H3BO3, 0.389 mg/L MnCl2·4H2O, 0.059 mg/L NiCl2·6H2O, 0.055 mg/L CoCl2·6H2O, 0.055 mg/L CuSO4·5H2O, 0.055 mg/L ZnSO4·7H2O, 0.0386 mg/L Al2(SO4)3·18H2O, 0.028 mg/L KBr, 0.028 mg/L KI, 0.028 mg/L LiCl, 0.028 mg/L SnCl2·2H2O) and 200 mg/L glucose. Medium was solidified with 1.5% (w/v) plant agar. Plants were cultured under standard conditions in a growth chamber at 20 ± 1 °C, under a 16/8 h light/dark photoperiod with a light intensity of 90 µmol photons m−2 s−1. Plants were subcultured onto fresh medium every 3 weeks.

Stichococcus bacillaris — green algae were sampled in Nové Hrady (Czech Republic) in summer of 2013. Briefly, green algae were first pre-cultivated in batch cultures in shaking 100 mL Erlenmeyer flasks that were irradiated by 40 µmol photon m−2 s−1 from fluorescence tubes (Osram T8 Biolux). Late exponential growth was achieved within 5 days following the inoculation. Then, the cultures were transferred to glass tubes (V = 300 mL) in which the cultures were aerated and cultivated in a stronger light (120 µmol photon m−2 s−1) and in Bold´s basal medium [70]. Sample material was lyophilized, homogenized under liquid nitrogen and then used in this study for phytohormone quantification.

When harvested, all plant tissues were immediately frozen in liquid nitrogen and stored at − 80 °C.

Sample extraction

Plant tissue was homogenised and powdered with a mortar and pestle under liquid nitrogen and reweighed into samples of approximately 10 mg FW or 1 mg of DW (Stichococcus bacillaris). The samples were extracted in cold 1 mol/L formic acid in 10% aqueous methanol so that 1 mL of the extraction solvent contained 10 mg FW or 1 mg DW. Internal standards were added to the samples to give 5 pmols of [13C6]-IAA, [13C6]-oxIAA, [13C6]-IAA-Asp, [13C6]-IAA-Glu, [13C6]-IAA-glc, [13C6]-oxIAA-glc, [2H5]-OPDA, [2H2]-( − )-JA-Ile, [2H6]-( ±)-MeJA and 10 pmols of [2H6]-JA, [2H6]-ABA, [2H3]-PA, [2H3]-DPA, [2H3]-neoPA, [2H4]-7′-OHABA, [2H4]-SA per 1 mL. To each of the samples, 4 ceria-stabilised zirconium oxide 2 mm beads (Retsch GmbH, Haan, Germany) were added, the samples were placed on an MM 400 mixer mill (Retsch GmbH, Haan, Germany) (29 Hz, 10 min, precooled holders), and centrifuged (25 200 g, 20 min, 8 °C), and 100 μL (1 mg of FW, 0.1 mg DW) or 200 μL (2 mg of FW, 0.2 mg of DW) of supernatant was loaded on a stagetip (Fig. 2A). The samples were handled at 8 °C in a CoolRack® in a CoolBox (Biocision, Larkspur, CA, USA) in the course of the sample preparation steps.

In-tip microSPE high-throughput purification

For in-tip microSPE columns, the stagetips [38] were self-assembled using ordinary pipette tips (2–200 μL, Eppendorf) and 3 layers of SDB-XC and C18 sorbent (3 M Empore, St. Paul, MN, USA) similarly to the method published for auxin metabolite profiling [40] (Fig. 2B).

A 96-place tip-holder was projected in 123D Design (AutoDesk Inc., San Rafael, CA, USA) and printed by DeltiX (Trilab, Czech Republic) using polylactic acid filament. This holder was designed to be compatible with regular 96-well plates and a centrifuge rotor for 96-well plates (M-20, 75003624, Heraeus Megafuge 16 K centrifuge, Thermo Fisher Scientific, Waltham, MA, USA) (Fig. 2B). The design and dimensions of the device are provided in Additional file 1: Fig. S4. Depending on number of samples, the stagetips were accommodated in 3D-printed 96-place tip-holders adaptable to process different number of large sets of samples, up to 192 samples in one run.

The purification protocol was adopted from [40] with modifications (Fig. 2A). The stagetips, located in a centrifuge, were conditioned with 50 μL of acetone (320 g, 10 min, 8 °C), 50 μL of methanol (320 g, 10 min, 8 °C), and 50 μL of water (627 g, 15 min, 8 °C). 100 or 200 μL of sample extract (1742 g, 20 min; 8 °C) was applied, washed with 50 μL of 0.1% aqueous acetic acid in water (819 g, 20 min, 8 °C) and eluted with 80% aqueous methanol (1280 g, 20 min, 8 °C) into a 96-well plate suitable for direct injection onto an LC–MS/MS system. Apart from the sample loads, all liquids were applied with an eight-channel pipette (Fig. 2B). Finally, the eluents in the 96-well plate were lyophilized (Labconco, Kansas city, MO, USA) and the plate stored at − 20 °C until required for analysis. When evaporation techniques tested TurboVap LV system (Caliper Life Sciences, Hopkinton, MA, USA) for evaporation under gentle stream of nitrogen or Acid-Resistant CentriVaps benchtop concentrator (Labconco, Kansas city, MO, USA) for evaporation in vacuo were used.

LC–MS/MS analysis

All analyses were carried out on an Agilent 6490 Triple Quadrupole LC/MS system coupled to a 1290 Infinity LC system (Agilent Technologies, Santa Clara, CA, USA). The data were processed in the MassHunter Quantitative software package version B.09.00 (Agilent Technologies, Santa Clara, CA, USA).

Purified and freeze-dried extracts were each dissolved in 40 μL of 20% aqueous acetonitrile and injected (10 μL) onto a Kinetex Evo C18 reverse phase column 2.1 ×150 mm, particle size 2.6 μm (Phenomenex, Torrance, CA, USA) protected by an inlet filter, using 10 mmol/L formic acid in water and methanol as mobile phases at a flow rate of 0.3 mL/min and a temperature of 50 °C. Gradient elution started at 2 min, going from 20 to 90% of methanol at 13.5 min, continuing within 0.5 min to 100% methanol. After 1 min of 100% methanol, the methanol content was decreased to 20% over the next 0.5 min and the column was equilibrated with 20% methanol for 3.5 min before the next injection.

The MS system was operated in dynamic multiple reaction monitoring mode in ESI positive and negative ionisation mode simultaneously (Table 1). The MS settings, multiple reaction monitoring transitions, and collision energies were optimised to previous standards. The nozzle voltage was set to 0 V, the capillary voltage to 2800/3000 V positive/negative mode, and the drying gas was at 130 °C with a flow rate of 14 L/min. The sheath gas was heated to 400 °C and its flow rate was set to 12 L/min.

Method validation

The method was validated in terms of accuracy and precision at five concentration levels using 2 mg FW of 10-day-old Arabidopsis thaliana seedlings spiked with authentic (0.1, 0.5, 1.0, 5.0 and 10.0 pmol) and internal standards (1 pmol for all AUXs, JAs detected in ESI + , 2 pmol for JAs, ABAs, SA detected in ESI-) (Table 1) with four replicates. The accuracy was calculated using the determined levels of the analytes (pmols) with the endogenous levels subtracted divided by the nominal level of the spike (pmol) and expressed as a percentage of the nominal level. Precision was calculated as the relative standard deviation of determined levels (%). The PE, RE and ME were assessed as in [49] using 2 mg FW of plant matrix spiked before or after extraction and the purification procedure at one concentration level (4 pmol) with four replicates. Briefly, to evaluate PE, the mean peak area of the IS spiked before sample preparation was divided by the mean peak area of the neat solution of IS without any extraction and purification, expressed as a percentage. RE was calculated as the percentage of the mean peak rate of authentic standards spiked before and after the extraction and purification process. ME was calculated by dividing the mean peak area of authentic standards spiked after sample preparation by the peak area of the neat solution with the same amount of analyte. All validation samples were processed by the optimised extraction and purification protocol (Fig. 2).

To evaluate the influence of the method of drying the samples after purification, 50 μL aliquots of solutions of all authentic standards (at 0.5 pmol) in 80% aqueous methanol were processed using evaporation under a stream of nitrogen, in vacuo drying and lyophilisation, with four replicates of each. The dried samples were redissolved in 50 μL 20% aqueous acetonitrile mimicking the conditions immediately before LC–MS/MS analysis. The recovery was calculated by dividing the peak area of a dried standard by that of the original authentic standard solution, expressed as a percentage.

Calibration and quantification

The contents of all samples were quantified using logarithmically transformed calibration curves constructed by plotting the responses of calibration standards (analyte peak area divided by IS area multiplied by IS concentration) against their known concentrations. The calibrations spanned from six to three orders of magnitude (Table 1).

Amount of plant matrix required for phytohormone profiling

The weight of sample necessary for quantification of the analytes was evaluated using 1–2–4–6–8 mg FW or 0.1–0.2–0.4–0.6–0.8 mg DW of plant matrix. Either 1 or 2 mg FW (0.1 or 0.2 mg of DW) plant matrix were loaded on a stagetip, and eluates of 2 mg (0.2 mg of DW) were combined to obtain samples of 4 to 8 mg FW (0.4–0.8 mg DW). Each sample amount was analysed in either triplicate, or, for 1 mg FW (0.1 mg DW), in duplicate. The mean concentration of each analyte was calculated from all samples in which the analyte was quantified, e.g. an analyte quantified at 1 mg (n = 2) was also quantified at 2 mg (n = 3), 4 mg (n = 3), 6 mg (n = 3) and 8 mg (n = 3), giving a total of 14 samples (n = 14).

Availability of data and materials

The data used for this manuscript are available from the corresponding author on request.



7′-Hydroxyabscisic acid


9,10-Dihydrojasmonic acid


11-Hydroxyjasmonic acid


12-Hydroxyjasmonic acid


Abscisic acid






cis-( +)-12-Oxo-phytodienoic acid


Dinor-oxo-phytodienoic acid


Dihydrophaseic acid


Dry weight


Electrospray ionisation


Fresh weight


Indole-3-acetic acid


Indole-3-acetyl aspartate


Indole-3-acetyl glutamate




Internal standard


Jasmonic acid






Liquid chromatography


Matrix effects


Methyl jasmonate


Mass spectrometry


Tandem mass spectrometry


Neophaseic acid


2-Oxindole-3-acetic acid




(Z)-4-[3-oxo-2-(pent-2-enyl)cyclopentyl]butanoic acid


(Z)-6-[3-Oxo-2-(pent-2-enyl)cyclopentyl]hexanoic acid


(Z)-8-[3-Oxo-2-(pent-2-enyl)cyclopentyl]octanoic acid


Phaseic acid


Process efficiency




Salicylic acid




Standard deviation


Solid phase extraction


  1. Casanova-Saez R, Mateo-Bonmati E, Ljung K. Auxin metabolism in plants. Cold Spring Harb Perspect Biol. 2021;13: a039867.

    Article  CAS  PubMed  Google Scholar 

  2. Hayashi K, Arai K, Aoi Y, Tanaka Y, Hira H, Guo RP, et al. The main oxidative inactivation pathway of the plant hormone auxin. Nat Commun. 2021;12:6752.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Porco S, Pencik A, Rashed A, Voss U, Casanova-Saez R, Bishopp A, et al. Dioxygenase-encoding AtDAO1 gene controls IAA oxidation and homeostasis in Arabidopsis. Proc Natl Acad Sci. 2016;113:11016–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Tanaka K, Hayashi K, Natsume M, Kamiya Y, Sakakibara H, Kawaide H, et al. UGT74D1 catalyzes the glucosylation of 2-Oxindole-3-acetic acid in the auxin metabolic pathway in arabidopsis. Plant Cell Physiol. 2014;55:218–28.

    Article  CAS  PubMed  Google Scholar 

  5. Brunoni F, Collani S, Simura J, Schmid M, Bellini C, Ljung K. A bacterial assay for rapid screening of IAA catabolic enzymes. Plant Methods. 2019;15(1):1–10.

    Article  CAS  Google Scholar 

  6. Muller K, Dobrev PI, Pencik A, Hosek P, Vondrakova Z, Filepova R, et al. dioxygenase for auxin oxidation 1 catalyzes the oxidation of IAA amino acid conjugates. Plant Physiol. 2021;187:103–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Ludwig-Muller J. Auxin conjugates: their role for plant development and in the evolution of land plants. J Exp Bot. 2011;62:1757–73.

    Article  PubMed  Google Scholar 

  8. Cooke TJ, Poli D, Sztein AE, Cohen JD. Evolutionary patterns in auxin action. Plant Mol Biol. 2002;49:319–38.

    Article  CAS  PubMed  Google Scholar 

  9. Zhang J, Peer WA. Auxin homeostasis: the DAO of catabolism. J Exp Bot. 2017;68:3145–54.

    Article  CAS  PubMed  Google Scholar 

  10. Brunoni F, Collani S, Casanova-Saez R, Simura J, Karady M, Schmid M, et al. Conifers exhibit a characteristic inactivation of auxin to maintain tissue homeostasis. New Phytol. 2020;226:1753–65.

    Article  CAS  PubMed  Google Scholar 

  11. Wasternack C, Hause B. Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. an update to the review in annals of botany. Ann Bot. 2007;2013(111):1021–58.

    Google Scholar 

  12. Fonseca S, Chini A, Hamberg M, Adie B, Porzel A, Kramell R, et al. (+)-7-iso-Jasmonoyl-L-isoleucine is the endogenous bioactive jasmonate. Nat Chem Biol. 2009;5:344–50.

    Article  CAS  PubMed  Google Scholar 

  13. Wasternack C. Jasmonates: an update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann Bot. 2007;100:681–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Han GZ. Evolution of jasmonate biosynthesis and signaling mechanisms. J Exp Bot. 2017;68:1323–31.

    CAS  PubMed  Google Scholar 

  15. Chen K, Li GJ, Bressan RA, Song CP, Zhu JK, Zhao Y. Abscisic acid dynamics, signaling, and functions in plants. J Integr Plant Biol. 2020;62:25–54.

    Article  CAS  PubMed  Google Scholar 

  16. Nambara E, Marion-Poll A. Abscisic acid biosynthesis and catabolism. Ann Rev Plant Biol. 2005;56:165–85.

    Article  CAS  Google Scholar 

  17. Zhang YL, Li X. Salicylic acid: biosynthesis, perception, and contributions to plant immunity. Curr Opin Plant Biol. 2019;50:29–36.

    Article  CAS  PubMed  Google Scholar 

  18. Lefevere H, Bauters L, Gheysen G. Salicylic acid biosynthesis in plants. Front Plant Sci. 2020;11:338.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Dempsey MA, Vlot AC, Wildermuth MC, Klessig DF. Salicylic acid biosynthesis and metabolism. Arabidopsis Book. 2011;9: e0156.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Hedden P. The current status of research on gibberellin biosynthesis. Plant Cell Physiol. 2020;61:1832–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Deng T, Wu DP, Duan CF, Yan XH, Du Y, Zou J, et al. Spatial profiling of gibberellins in a single leaf based on microscale matrix solid-phase dispersion and precolumn derivatization coupled with ultraperformance liquid chromatography-tandem mass spectrometry. Anal Chem. 2017;89:9537–43.

    Article  CAS  PubMed  Google Scholar 

  22. Hao YH, Zhang Z, Wang L, Liu C, Lei AW, Yuan BF, et al. Stable isotope labeling assisted liquid chromatography-electrospray tandem mass spectrometry for quantitative analysis of endogenous gibberellins. Talanta. 2015;144:341–8.

    Article  CAS  PubMed  Google Scholar 

  23. Urbanova T, Tarkowska D, Novak O, Hedden P, Strnad M. Analysis of gibberellins as free acids by ultra performance liquid chromatography-tandem mass spectrometry. Talanta. 2013;112:85–94.

    Article  CAS  PubMed  Google Scholar 

  24. Ljung K. Auxin metabolism and homeostasis during plant development. Development. 2013;140:943–50.

    Article  CAS  PubMed  Google Scholar 

  25. Wasternack C, Strnad M. Jasmonate signaling in plant stress responses and development—active and inactive compounds. New Biotechnol. 2016;33:604–13.

    Article  CAS  Google Scholar 

  26. Maruri-Lopez I, Yaniri Aviles-Baltazar N, Buchala A, Serrano M. Intra and extracellular journey of the phytohormone salicylic acid. Front Plant Sci. 2019;10:423.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Trapp MA, De Souza GD, Rodrigues E, Boland W, Mithofer A. Validated method for phytohormone quantification in plants. Front Plant Sci. 2014;5:417.

    Google Scholar 

  28. Erland LAE, Shukla MR, Glover WB, Saxena PK. A simple and efficient method for analysis of plant growth regulators: a new tool in the chest to combat recalcitrance in plant tissue culture. Plant Cell Tissue Organ Cult. 2017;131:459–70.

    Article  CAS  Google Scholar 

  29. Sheflin AM, Kirkwood JS, Wolfe LM, Jahn CE, Broeckling CD, Schachtman DP, et al. High-throughput quantitative analysis of phytohormones in sorghum leaf and root tissue by ultra-performance liquid chromatography-mass spectrometry. Anal Bioanal Chem. 2019;411:4839–48.

    Article  CAS  PubMed  Google Scholar 

  30. Yonny ME, Ballesteros-Gomez A, Adamo MLT, Torresi AR, Nazareno MA, Rubio S. Supramolecular solvent-based high-throughput sample treatment for monitoring phytohormones in plant tissues. Talanta. 2020;219: 121249.

    Article  CAS  PubMed  Google Scholar 

  31. Cai WJ, Yu L, Wang W, Sun MX, Feng YQ. Simultaneous determination of multiclass phytohormones in submilligram plant samples by one-pot multifunctional derivatization-assisted liquid chromatography-tandem mass spectrometry. Anal Chem. 2019;91:3492–9.

    Article  CAS  PubMed  Google Scholar 

  32. Cao ZY, Sun LH, Mou RX, Zhang LP, Lin XY, Zhu ZW, et al. Profiling of phytohormones and their major metabolites in rice using binary solid-phase extraction and liquid chromatography-triple quadrupole mass spectrometry. J Chromatogr A. 2016;1451:67–74.

    Article  CAS  PubMed  Google Scholar 

  33. Flokova K, Tarkowska D, Miersch O, Strnad M, Wasternack C, Novak O. UHPLC-MS/MS based target profiling of stress-induced phytohormones. Phytochemistry. 2014;105:147–57.

    Article  CAS  PubMed  Google Scholar 

  34. Simura J, Antoniadi I, Siroka J, Tarkowska D, Strnad M, Ljung K, et al. Plant hormonomics: multiple phytohormone profiling by targeted metabolomics. Plant Physiol. 2018;177:476–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Balcke GU, Handrick V, Bergau N, Fichtner M, Henning A, Stellmach H, et al. An UPLC-MS/MS method for highly sensitive high-throughput analysis of phytohormones in plant tissues. Plant Methods. 2012;8:47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Miggiels P, Wouters B, van Westen GJP, Dubbelman AC, Hankemeier T. Novel technologies for metabolomics: more for less. Trends Anal Chem. 2019;120: 115323.

    Article  Google Scholar 

  37. Burato JSD, Medina DAV, de Toffoli AL, Maciel EVS, Lancas FM. Recent advances and trends in miniaturized sample preparation techniques. J Sep Sci. 2020;43:202–25.

    Article  Google Scholar 

  38. Rappsilber J, Mann M, Ishihama Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using stagetips. Nat Protoc. 2007;2:1896–906.

    Article  CAS  PubMed  Google Scholar 

  39. Svacinova J, Novak O, Plackova L, Lenobel R, Holik J, Strnad M, et al. A new approach for cytokinin isolation from Arabidopsis tissues using miniaturized purification: pipette tip solid-phase extraction. Plant Methods. 2012;8:17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Pencik A, Casanova-Saez R, Pilarova V, Zukauskaite A, Pinto R, Micol JL, et al. Ultra-rapid auxin metabolite profiling for high-throughput mutant screening in Arabidopsis. J Exp Bot. 2018;69:2569–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Miersch O, Neumerkel J, Dippe M, Stenzel I, Wasternack C. Hydroxylated jasmonates are commonly occurring metabolites of jasmonic acid and contribute to a partial switch-off in jasmonate signaling. New Phytol. 2008;177:114–27.

    Article  CAS  PubMed  Google Scholar 

  42. Glauser G, Grata E, Dubugnon L, Rudaz S, Farmer EE, Wolfender JL. Spatial and temporal dynamics of jasmonate synthesis and accumulation in Arabidopsis in response to wounding. J Biol Chem. 2008;283:16400–7.

    Article  CAS  PubMed  Google Scholar 

  43. Novak O, Henykova E, Sairanen I, Kowalczyk M, Pospisil T, Ljung K. Tissue-specific profiling of the Arabidopsis thaliana auxin metabolome. Plant J. 2012;72:523–36.

    Article  CAS  PubMed  Google Scholar 

  44. Tureckova V, Novak O, Strnad M. Profiling ABA metabolites in Nicotiana tabacum L. leaves by ultra-performance liquid chromatography-electrospray tandem mass spectrometry. Talanta. 2009;80:390–9.

    Article  CAS  PubMed  Google Scholar 

  45. Stelmach BA, Muller A, Hennig P, Gebhardt S, Schubert-Zsilavecz M, Weiler EW. A novel class of oxylipins, sn1-O-(12-oxophytodienoyl)-sn2-O-(hexadecatrienoyl)-monogalactosyl diglyceride, from Arabidopsis thaliana. J Biol Chem. 2001;276:12832–8.

    Article  CAS  PubMed  Google Scholar 

  46. Genva M, Akong FO, Andersson MX, Deleu M, Lins L, Fauconnier ML. New insights into the biosynthesis of esterified oxylipins and their involvement in plant defense and developmental mechanisms. Phytochem Rev. 2019;18:343–58.

    Article  CAS  Google Scholar 

  47. Keunchkarian S, Reta M, Romero L, Castells C. Effect of sample solvent on the chromatographic peak shape of analytes eluted under reversed-phase liquid chromatogaphic conditions. J Chromatogr A. 2006;1119:20–8.

    Article  CAS  PubMed  Google Scholar 

  48. Jang G, Shim JS, Jung C, Song JT, Lee HY, Chung PJ, et al. Volatile methyl jasmonate is a transmissible form of jasmonate and its biosynthesis is involved in systemic jasmonate response in wounding. Plant Biotechnol Rep. 2014;8:409–19.

    Article  Google Scholar 

  49. Matuszewski BK, Constanzer ML, Chavez-Eng CM. Strategies for the assessment of matrix effect in quantitative bioanalytical methods based on HPLC-MS/MS. Anal Chem. 2003;75:3019–30.

    Article  CAS  PubMed  Google Scholar 

  50. Stirk WA, Ordog V, Novak O, Rolcik J, Strnad M, Balint P, et al. Auxin and cytokinin relationships in 24 microalgal strains. J Phycol. 2013;49:459–67.

    Article  CAS  PubMed  Google Scholar 

  51. Zaveska Drabkova L, Dobrev PI, Motyka V. Phytohormone profiling across the bryophytes. PLoS ONE. 2015;10: e0125411.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Kosakivska IV, Voytenko LV, Likhnyovskiy RV, Ustinova AY. Effect of temperature on accumulation of abscisic acid and indole-3-acetic acid in Triticum aestivum L seedlings. Genet Plant Physiol. 2014;4(3–4):201–8.

    Google Scholar 

  53. Pavlovic I, Petrik I, Tarkowska D, Lepedus H, Vujcic Bok V, Radic Brkanac S, et al. Correlations between phytohormones and drought tolerance in selected Brassica crops: Chinese cabbage, white cabbage and kale. Int J Mol Sci. 2018;19:2866.

    Article  PubMed Central  Google Scholar 

  54. Junghans U, Polle A, Duchting P, Weiler E, Kuhlman B, Gruber F, et al. Adaptation to high salinity in poplar involves changes in xylem anatomy and auxin physiology. Plant Cell Environ. 2006;29:1519–31.

    Article  CAS  PubMed  Google Scholar 

  55. De Zio E, Trupiano D, Karady M, Antoniadi I, Montagnoli A, Terzaghi M, et al. Tissue-specific hormone profiles from woody poplar roots under bending stress. Physiol Plant. 2019;165:101–13.

    Article  PubMed  Google Scholar 

  56. Stumpe M, Gobel C, Faltin B, Beike AK, Hause B, Himmelsbach K, et al. The moss Physcomitrella patens contains cyclopentenones but no jasmonates: mutations in allene oxide cyclase lead to reduced fertility and altered sporophyte morphology. New Phytol. 2010;188:740–9.

    Article  CAS  PubMed  Google Scholar 

  57. De Leon IP, Schmelz EA, Gaggero C, Castro A, Alvarez A, Montesano M. Physcomitrella patens activates reinforcement of the cell wall, programmed cell death and accumulation of evolutionary conserved defence signals, such as salicylic acid and 12-oxo-phytodienoic acid, but not jasmonic acid, upon Botrytis cinerea infection. Mol Plant Pathol. 2012;13:960–74.

    Article  Google Scholar 

  58. Yamamoto Y, Ohshika J, Takahashi T, Ishizaki K, Kohchi T, Matusuura H, et al. Functional analysis of allene oxide cyclase, MpAOC, in the liverwort Marchantia polymorpha. Phytochemistry. 2015;116:48–56.

    Article  CAS  PubMed  Google Scholar 

  59. Stintzi A, Weber H, Reymond P, Browse J, Farmer EE. Plant defense in the absence of jasmonic acid: The role of cyclopentenones. Proc Natl Acad Sci. 2001;98:12837–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Monte I, Kneeshaw S, Franco-Zorrilla JM, Chini A, Zamarreno AM, Garcia-Mina JM, et al. An ancient COI1-independent function for reactive electrophilic oxylipins in thermotolerance. Curr Biol. 2020;30:962–71.

    Article  CAS  PubMed  Google Scholar 

  61. Alallaq S, Ranjan A, Brunoni F, Novak O, Lakehal A, Bellini C. Red light controls adventitious root regeneration by modulating hormone homeostasis in Picea abies seedlings. Front Plant Sci. 2020;11: 586140.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Hartung W. The evolution of abscisic acid (ABA) and ABA function in lower plants, fungi and lichen. Funct Plant Biol. 2010;37:806–12.

    Article  CAS  Google Scholar 

  63. Pashkovskiy PP, Vankova R, Zlobin IE, Dobrev P, Ivanov YV, Kartashov AV, et al. Comparative analysis of abscisic acid levels and expression of abscisic acid-related genes in scots pine and Norway spruce seedlings under water deficit. Plant Physiol Biochem. 2019;140:105–12.

    Article  CAS  PubMed  Google Scholar 

  64. Kramell R, Schmidt J, Schneider G, Sembdner G, Schreiber K. Synthesis Of N-(Jasmonoyl)amino acid conjugates. Tetrahedron. 1988;44:5791–807.

    Article  CAS  Google Scholar 

  65. Kowalczyk M, Sandberg G. Quantitative analysis of indole-3-acetic acid metabolites in Arabidopsis. Plant Physiol. 2001;127:1845–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Kai K, Horita J, Wakasa K, Miyagawa H. Three oxidative metabolites of indole-3-acetic acid from Arabidopsis thaliana. Phytochemistry. 2007;68:1651–63.

    Article  CAS  PubMed  Google Scholar 

  67. Kai K, Nakamura S, Wakasa K, Miyagawa H. Facile preparation of deuterium-labeled standards of indole-3-acetic acid (IAA) and its metabolites to quantitatively analyze the disposition of exogenous IAA in Arabidopsis thaliana. Biosci Biotechnol Biochem. 2007;71:1946–54.

    Article  CAS  PubMed  Google Scholar 

  68. Pavlovic I, Pencik A, Novak O, Vujcic V, Brkanac SR, Lepedus H, et al. Short-term salt stress in Brassica rapa seedlings causes alterations in auxin metabolism. Plant Physiol Biochem. 2018;125:74–84.

    Article  CAS  PubMed  Google Scholar 

  69. Hoagland DR, Arnon DI. The water-culture method for growing plants without soil. 2nd ed. California: California agricultural experiment station; 1950.

    Google Scholar 

  70. Bischoff HW, Bold HC. Phycological studies IV some soil algae from enchanted rock and related algal species. Texas: Univ Texas Publ; 1963.

    Google Scholar 

Download references


The authors thank Hana Svobodová and Miroslava Špičáková for technical support, and David Friedecký and Lukáš Kouřil for 3D printing.


This work was financially supported by Czech Science Foundation project No. 19-10464Y and by the Palacký University Olomouc Young Researcher grant no. JG_2020_002.

Author information

Authors and Affiliations



JŠ designed and performed the analytical experiments, wrote and organised the manuscript. FB wrote the biological part of the background and the results and discussion sections of the manuscript. ON contributed to manuscript writing. AP supervised development of the sample preparation protocol. VM and AŽ synthesised analytical standards. KF, MS and ON supervised the research, discussed results and reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jitka Široká.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Fig. S1. Purification recoveries of standards in acidified and non-acidified solution. Calculated as mean area for standard solution (1 pmol) processed by microSPE divided by the area for the corresponding amount of standard without microSPE × 100. Error bars represent ± SD, n = 3. Fig. S2. Endogenous levels of dn-OPDA and cis-OPDA found in 10-day-old Arabidopsis thaliana seedlings using different extraction conditions over time. Error bars represent ± SD, n = 4. Fig. S3. Influence of the method of drying samples after purification using evaporation under a stream of nitrogen, in vacuo drying and lyophilisation, expressed as recovery after evaporation. Calculated as mean peak area for dried and re-dissolved standard divided by area for the standard solution at the same concentration without evaporation, expressed as a percentage. Error bars represent ± SD, n = 4. Fig. S4. Design and dimensions of 3D printed 96-place microSPE holder. Table S1. Range of sample weights (1–2–4–6–8 mg of FW or 0.1–0.2–0.4–0.6–0.8 mg of DW) in which the analyte was quantified, R2 of regression of sample weight and analyte level found in sample.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Široká, J., Brunoni, F., Pěnčík, A. et al. High-throughput interspecies profiling of acidic plant hormones using miniaturised sample processing. Plant Methods 18, 122 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Plant hormones
  • Miniaturisation
  • In-tip microSPE
  • High-throughput
  • 3D printing
  • Liquid chromatography
  • Mass spectrometry
  • Evolutionarily distant plant species