- Department of Plant Physiology, Institute for Biological Research “Siniša Stanković” - National Institute of the Republic of Serbia, University of Belgrade, Belgrade, Serbia
Nepeta nuda L. shares a typical secondary chemistry with other Nepeta species (fam. Lamiaceae), characterized by the tendency to intensively produce monoterpenoid iridoids, whereas the phenylpropanoid chemistry is steered towards the production of a caffeic acid ester, rosmarinic acid. Combining complementary state-of-the-art analytical techniques, N. nuda metabolome was here comprehensively characterized in the quest for the organ-specific composition of phenolics and terpenoids that possess well-defined functions in plant-biotic interactions as well as therapeutic potential. N. nuda inflorescences showed generally higher constitutive levels of specialized metabolites, as compared to leaves, and the composition of major iridoids and phenolics in reproductive organs was found to be more conserved than in leaves across 13 populations from the Central Balkans. The results suggest that N. nuda plants most likely invest more in constitutive than inducible biosynthesis of functional metabolites in flowers, since they are of essential importance for both pollination and defense against herbivores and pathogens. Conversely, specialized metabolism of leaves is found to be more susceptible to reprograming in response to differential growth conditions. The defense strategy of leaves, primarily functioning in CO2 fixation during photosynthesis, more likely relies on the induction of metabolite levels following plant-environment interplay. Organ-specific biosynthesis of iridoids in N. nuda is found to be tightly regulated at the transcriptional level, and high constitutive levels of these compounds in inflorescences most likely result from the up-regulated expression of several key genes (NnG8H, NnNEPS1, NnNEPS2, and NnNEPS3) determining the metabolic flux through the pathway. The organ-specific content of rosmarinic acid and co-expression patterns of the corresponding biosynthetic genes were much less correlated, which suggests independent organ-specific transcriptional regulation of the iridoid and phenolic pathways. Knowledge gathered within the present study can assist growers to select productive genotypes and manipulate phenology of N. nuda towards maximizing yields and facilitating its integration into pest management systems and other applications related to human health.
1 Introduction
Plants produce a diverse array of specialized metabolites that have evolved specific physiological and ecological functions. These metabolites occur either as volatile or non-volatile compounds and play important roles in plant’s adaptation to different environmental conditions and defense against herbivores and microorganisms, as well as in processes such as pollination and seed dispersal (Wink, 2018). Phytochemical diversity varies across different levels of biological organization, i.e., among populations, among individuals within a population, and within an individual plant (among plant organs, during development, and across seasons) (Moore et al., 2014; Wetzel and Whitehead, 2020).
Variation in several categories of specialized metabolites is well studied at the interpopulation level of numerous species, including those belonging to the genus Nepeta. The interpopulation chemical polymorphism detected in Nepeta sessilifolia, N. heliotropifolia, and N. fissa was reported to be correlated with the elevation (Talebi et al., 2019). It has been suggested that the amounts of different types of monoterpenes decreased and the amounts of oxygenated compounds increased along the altitudinal gradient. Environmental conditions are highlighted as important factors influencing the yield and the chemical composition of essential oils of N. fissa (Talebi et al., 2017) and N. heliotropifolia (Yarmoohammadi et al., 2017). In another study, 29 accessions of N. kotschyi, N. menthoides, N. crassifolia, and N. cataria were cultivated under western Tehran (Iran) environmental conditions and evaluated for their phenolics composition (Hadi et al., 2017). Diverse composition of specialized metabolites in various plant tissues and organs is an important component of the overall plant phenotype, and changes in phytochemical profiles are one of the fastest phenotypic responses to the dynamic environmental conditions (Junker, 2016). These plant-part-specific differences may also be related to the ecological functions of the vegetative and reproductive tissues.
Nepeta nuda L. (subfam. Nepetoidae, fam. Lamiaceae) is an herbaceous perennial species widespread across Europe and Asia. Its individuals, characterized by numerous erect stems, oblong-lanceolate lower leaves, and ovate upper leaves can grow 50 to 100 cm tall (Aćimović et al., 2020). Flowers are organized in lax or dense spike-like verticillasters that bloom from June to August, depending on the altitude (Aćimović et al., 2020). Health benefits of N. nuda are usually associated with the accumulation of terpenoids and phenolics, two major groups of specialized metabolites of the genus Nepeta. Terpenoids, predominately represented by iridoid compounds, are produced in most plant tissues where they primarily act either as pollinator attractants (Schultz et al., 2004) or insect repellents (Bernier et al., 2005; Birkett et al., 2011; Sparks et al., 2018; Reichert et al., 2019; reviewed in Formisano et al., 2011 and Süntar et al., 2018), but can also function as insect pheromones (Glinwood et al., 1999; Birkett and Pickett, 2003) and are appealing to felines (Uenoyama et al., 2021 and 2022). Numerous studies on the chemical composition of N. nuda plant extracts and essential oils revealed a high level of chemodiversity, which may be ecologically driven, as it has been proven to be strongly influenced by various environmental factors, such as geographical origin, light conditions or altitude (Sharma et al., 2021; Petrović et al., 2024). Only few studies report on variations in iridoid and phenolic composition among populations of N. nuda within a specific geographical region (Narimani et al., 2017; Petrova et al., 2022; Petrović et al., 2024). In a recent study, we reported low differences in the iridoid and phenolic profiles in leaves among populations in the Central Balkans, since the majority of chemical and genetic variations were found within populations (Petrović et al., 2024). Moreover, other recent studies suggested organ-specific biosynthesis and accumulation of bioactive metabolites in N. nuda (Petrova et al., 2022; Zaharieva et al., 2023). Thus, the content of rosmarinic acid and 1,5,9-epi-deoxyloganic acid in N. nuda plants from Mt. Pirin (Bulgaria) was significantly higher in flowers than in leaves (Petrova et al., 2022).
The objective of this study was to comparatively investigate the patterns of metabolites’ variation in N. nuda vegetative (leaves) and reproductive (inflorescences) organs across 13 natural populations within the Central Balkans, adopting both targeted and untargeted metabolomics approaches. The determination of variation patterns in the plant metabolome in nature is crucial for predicting the persistence of populations and species in the future, understanding the species’ ecological interactions, and conserving their genetic diversity (Labarrere et al., 2019). Characterization of intraspecies chemodiversity in N. nuda within the Central Balkans can indirectly enable mapping of areas suitable for cultivation and intendedly assist farmers and other stakeholders to identify accessions with the highest crop potential and earmark them for future cultivation. N. nuda can be considered attractive medicinal crop for montane and subalpine regions (up to 2100 m a.s.l.), due to high adaptive success and production of essential oils with prominent bioactive properties (e.g., Gkinis et al., 2010; Formisano et al., 2011; Musso et al., 2017; reviewed in Sharma et al., 2021; Zaharieva et al., 2023). Considering differential but equally important ecophysiological roles of leaves and inflorescences, we intended to investigate which plant part is more susceptible to metabolic reprogramming in response to external influences and discuss the potential ecological implications arising from variations in the composition of terpenes and phenolics in these two organs of N. nuda. We further delved into the molecular background of the organ-specific differences in the biosynthesis and accumulation of iridoids and phenolics, thus providing evidence for their mutually independent transcriptional regulation of metabolic fluxes in different organs.
2 Materials and methods
2.1 Chemicals and reagents
Acetonitrile (Fisher Scientific UK, Leicestershire, UK) and formic acid (Merck, Darmstadt, Germany) were of MS grade. Ultra-pure deionized water was generated using a Water Purification System (New Human Power I Integrate, Human Corporation, Republic of Korea). Standards of 1,5,9-epi-deoxyloganic acid and 5,9-dehydronepetalactone were isolated from natural sources as previously described by Aničić et al. (2021). The standard of cis,trans-nepetalactone was a generous gift from Entomol Products LLC (San Francisco, CA, USA). Analytical standards of loganin, rosmarinic acid, caffeic acid, luteolin and apigenin were purchased from Sigma-Aldrich (Hamburg, Germany).
2.2 Collection of plant material
Aboveground parts of flowering N. nuda plants were collected in June-August 2022 from 13 populations across Serbia: 1. Mali Ljukten (43°32’35”N; 20°48’38”E); 2. Brodica (44°29’04”N; 21°50’27”E); 3. Debeli lug (44°21’45”N; 21°54’01”E); 4. Rtanj (43°44’13”N; 21°57’03”E); 5. Straža (43°50’46”N; 21°42’03”E); 6. Židilje (44°00’42”N; 21°38’45”E); 7. Vinatovača (44°04’14”N; 21°45’36”E); 8. Vlasina (42°41’39”N; 22°22’44”E); 9. Donji Krivodol (43°06’18”N; 22°55’41”E); 10. Topli Do (43°22’05”N; 22°37’58”E); 11. Janjska reka (43°25’33”N; 22°31’11”E); 12. Balta Berilovac (43°24’10”N; 22°30’43”E); 13. Gornje selo (42°11’21”N; 20°56’20”E). Plants were identified in the field by the authors, and representative specimens, assigned with the voucher numbers, have been deposited in the Herbarium of the Institute of Botany and Botanical Garden, University of Belgrade, BEOU (acronym follows Thiers, 2023). Plant material was dried in the shade at room temperature until constant mass and subsequently stored in paper bags. Samples were kept in the dark at room temperature until use. N. nuda aboveground parts were detached to leaves and inflorescences.
2.3 Preparation of N. nuda methanol extracts
Samples containing either leaves or inflorescences (50 mg) of N. nuda were ground in liquid nitrogen and extracted overnight at 4 °C with 1 mL of 96% methanol (w:v = 1:20). After vortexing for 1 min, samples were sonicated in an ultrasonic bath (RK100, Bandelin, Berlin, Germany) for 1 h at 4°C. Supernatants were separated after centrifugation at 10,000g for 10 min and subsequently filtered through 0.2 µm-pore size cellulose filters (Agilent Technologies, Santa Clara, CA, USA). Each population was represented by 3 biological replicates.
2.4 Characterization of metabolites in N. nuda leaves and inflorescences using UHPLC-ESI-QToF-MS analysis
Three randomly selected N. nuda samples (populations 3, 8, and 9) of both inflorescences and leaves were subjected to metabolic fingerprinting using an Agilent 1290 Infinity UHPLC system in combination with a quadrupole time-of-flight mass spectrometer (6530C QToF-MS, Agilent Technologies, Inc., Santa Clara, CA, USA). The mass detector, equipped with an electrospray ionization (ESI) source, was operated in both positive and negative ionization modes, in the mass range from 100 to 1000 m/z. Samples were chromatographically separated on a Zorbax C18 column (2.1 × 50 mm, particle size 1.8 µm; Agilent Technologies, Inc., Santa Clara, CA, USA) thermostated at 40°C. The gradient elution program, UHPLC and MS parameters as well as ion source settings were described by Petrović et al. (2024).
The Agilent MassHunter software was used for acquisition, control, and MS data collection from the 6530C Q-ToF-MS instrument. The R Studio platform (‘enviPick’ and ‘xcms’ R packages) was used for the evaluation and presentation of MS data (Zengin et al., 2021). The metabolites were identified based on their monoisotopic masses and MS2 fragmentation and confirmed using literature data (Alimpić Aradski et al., 2023; Aničić et al., 2021; Dienaitė et al., 2018; Goldansaz et al., 2019). Accurate masses of components and fragment ions were calculated using the ChemDraw software (version 12.0, CambridgeSoft, Cambridge, MA, USA). The CAS SciFinder-n database (https://scifinder-n.cas.org/) was used to search for chemical compounds by formulae and structures.
2.5 UHPLC/DAD/(±)HESI-MS2 targeted metabolic profiling of N. nuda leaves and inflorescences
A targeted UHPLC/DAD/(±)HESI/MS2 metabolomics approach was adopted to quantify individual phenolics (rosmarinic acid, caffeic acid, apigenin, and luteolin) and iridoids (cis,trans- and trans,trans-nepetalactone, 5,9-dehydronepetalactone, nepetanudoside A, and 1,5,9-epi-deoxyloganic acid) in N. nuda methanol extracts of leaves and inflorescences. Analyses were performed on a Dionex UltiMate 3000 UHPLC system coupled with a DAD detector and configured with a triple quadrupole mass spectrometer (TSQ Quantum Access Max, Thermo Fisher Scientific, Basel, Switzerland). Samples were chromatographically separated on a Hypersil gold C18 analytical column (50 × 2.1 mm, 1.9 µm particle size; Thermo Fisher Scientific, Waltham, MA, USA), using the elution gradient and the flow rate previously described by Aničić et al. (2021). The mobile phase consisted of (A) water + 0.1% formic acid and (B) acetonitrile + 0.1% formic acid. The injection volume was 10 µL. The selected reaction monitoring (SRM) mode of the instrument was used for the quantification of the targeted compounds by direct comparison with the commercial standards. Nepetanudoside A was quantified based on the calibration curve of loganin, while quantification of aucubin was performed using the calibration curve for nepetaside. Calibration curves of pure standards revealed good linearity, with r2 values exceeding 0.99 (peak areas vs. concentration). The total amount of each phenolic and iridoid compound was evaluated by the calculation of its peak area and is expressed as µg per 100 mg leaf/inflorescence DW.
2.6 GC/MS non-targeted metabolomics of methanol extracts of N. nuda leaves and inflorescences
Profiling of volatile compounds in methanol extracts of N. nuda leaves and inflorescences was performed using an Agilent 8890 gas chromatography (GC) system coupled with a Mass Selective Detector (5977B GC/MSD, Agilent Technologies, Santa Clara, CA, USA) and connected to an automated sample extraction and enrichment platform (Centri®, Markes International Ltd., Bridgend, UK). Chromatographic separation procedure is described in detail by Petrović et al. (2024). The constituents of the methanol extracts were identified by comparison of their mass spectra and retention times with those of the respective standards as well as by comparison with the NIST05 library.
2.7 Statistical analyses of metabolomics data
For hierarchical cluster analysis (HCA), the input variables were scaled to the [min, max] range. HCA was performed based on the Spearman’s method of cluster agglomeration, using the Morpheus software (https://software.broadinstitute.org/morpheus). Principal component analysis (PCA) was constructed using the Past 4 software (version 4.14; Hammer and Harper, 2001). For the comparison between populations, quantitative metabolomics data were subjected to the Tukey’s post hoc test (p < 0.05) of one-way ANOVA. To test the statistical difference between leaves and inflorescences, t test was adopted.
2.8 Comparative metabolic profiling and gene expression analysis in leaves and inflorescences of greenhouse-grown N. nuda
2.8.1 Plant material and metabolic profiling of iridoids and phenolics
Seeds of N. nuda, collected in June 2022 in the locality Debeli Lug (East Serbia) were germinated under in vitro conditions, as previously described (Aničić et al., 2024). Following clonal propagation of a single genotype under in vitro conditions, adopting previously optimized protocols for other Nepeta species (Aničić et al., 2024), plantlets were acclimatized in a greenhouse of the Institute for Biological Research “Siniša Stanković” - National Institute of the Republic of Serbia, University of Belgrade, Serbia. Samples of leaves and inflorescences were collected from fully-flowering plants in their third year of growth under greenhouse conditions.
Leaves and inflorescences were excised from plants and immediately frozen in liquid nitrogen (LN). Samples were mechanically grinded in LN until fine powder, weighted and aliquoted into batches of 100 mg, and subsequently stored at −80°C until use. The first group of batches containing three biological replicates, designated for metabolomics, was dried by lyophilization until constant mass, weighted, and extracted with 96% methanol (100:1 = w:v). Following vigorous vortexing for 1 min, the extraction was continued in an ultrasonic bath for 1 h, and supernatants were separated by centrifugation for 10 min at 10,000g. Supernatants were filtered through 0.2-mm cellulose filters (Agilent Technologies, Santa Clara, CA, USA) and stored at 4°C until use. Quantification of the targeted iridoids and phenolics in samples was performed using a UHPLC-DAD-(±)HESI/MS2 instrument as described in section 2.5.
2.8.2 RNA extraction and qPCR assessment of the expression of biosynthetic genes
The second batch of samples was used to extract RNA from N. nuda leaves and inflorescences, and all the extractions were performed in three biological replicates. The RNA extraction was performed applying a modified CTAB method (Gasic et al., 2004), which was followed by the quantification using a N60 Nano-Photometer® (Implen GmbH, Munich, Germany) and a Qubit 3.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA), while the RNA integrity was electrophoretically checked. The procedure of the DNase I (Termo Fisher Scientific, USA) treatment and the protocol for cDNA synthesis were described earlier by Aničić et al. (2024).
Gene expression analyses were performed by real-time PCR using QuantStudio™ 3 Real-Time PCR System (Thermo Fisher Scientific, Carlsbad, CA, USA). Thermocycler conditions were previously described by Aničić et al. (2018). The reactions were performed using Maxima SYBR Green/ROX Master Mix (2X) (Termo Fisher Scientific, USA), cDNA corresponding to 50 ng RNA, and 0.3 µM primers, according to the manufacturer’s recommendations. The expression levels of targeted candidates of iridoid- and phenolic-biosynthetic genes were calculated according to the 2–ΔΔCt method (Livak and Schmittgen, 2001) using GAPDH as a housekeeping gene. Primer pairs were designed based on candidate genes whose sequences and accession numbers are listed in Supplementary Table 1.
2.8.3 Statistical analysis of gene expression data
qPCR data were processed and visualized using the Past 4 software (version 4.16; Hammer and Harper, 2001). To test the statistical difference between leaves and inflorescences, t test was adopted. Correlation analysis, adopting the Ward’s algorithm, was performed to examine the co-expression patterns of targeted genes in leaves and inflorescences of N. nuda.
Linear regression analysis was performed to visualize relations between metabolomics and gene-expression data. Coefficients of determination (R2) were calculated using the Past 4 software (version 4.16; Hammer and Harper, 2001).
3 Results and discussion
In this study, we adopted two complementary analytical techniques to cover different parts of N. nuda metabolome. A UHPLC-ESI-QToF-MS instrument, operating in both positive and negative ionization modes, enabled the chromatographic separation and identification of polar molecules in methanol extracts, primarily phenolics. It enabled us to comparatively and simultaneously analyze both glycosides and aglycones of iridoids in leaves and inflorescences. The solvent of choice for the present study was methanol as being vastly used for the extraction of both phenolics (e.g., Mišić et al., 2015; Hadi et al., 2017; Dienaitė et al., 2018; Aničić et al., 2021) and terpenoids (Aničić et al., 2024), including iridoid glycosides and volatile iridoid aglycones found in Nepeta species (Takeda et al., 1995; Petrova et al., 2022; Zaharieva et al., 2023).
3.1 Comparison of qualitative non-targeted phenolic profiles in N. nuda leaves and inflorescences
Totally 38 phenolic compounds were identified, including derivatives of hydroxybenzoic and hydroxycinnamic acids as well as flavonoid glycosides and aglycones (Supplementary Table 2).
Nine derivatives of hydroxybenzoic acid were identified mainly as hexosides (Supplementary Table 2), giving a specific loss of 162 Da in MS2 spectra. All hydroxybenzoic acid derivatives were detected in leaves, while only one compound, dihydroxybenzoic acid (3), was not identified in the tested inflorescences (Supplementary Table 2). The majority of compounds belonging to the group of hydroxybenzoic acids were more abundant in inflorescences than in leaves of N. nuda except galloyl hexoside (1) and dihydroxybenzoic acid hexoside (4) (Figure 1A).
Figure 1. (A) Heatmap of the scaled QToF-MS data with the samples and compounds arranged according to the hierarchical cluster analysis (HCA), adopting the Spearman’s method of cluster agglomeration. (B) Chemical structures of the characteristic iridoid aglycones identified in Nepeta nuda methanol extracts of leaves and inflorescences. Iridoids predominating in inflorescences are represented in blue, while those prevalent in leaves are presented in green. (C) Representative TIC GC/MS chromatograms of methanol extracts of leaves and inflorescences. HCA plot of the scaled GC/MS data with populations and compounds arranged according to the Spearman’s method of cluster agglomeration. Plant photo denotes N. nuda individual captured on site (locality Straža, SE Serbia).
All nineteen compounds from the group of hydroxycinnamic acids were detected in inflorescences, whereas compounds 13, 24, 27, and 28 were not found in leaves (Supplementary Table 2). Interestingly, esters of caffeic acid (23), malic acid (25), tartaric acid (11 and 21), tartronic acid (12 and 16), and glycolic acid (14) were more abundant in leaves than in inflorescences. The majority of other identified hydroxycinnamic acids were more abundant in inflorescences (Figure 1A). Compound 24 was not previously detected in any Nepeta species, but its MS2 spectra showed fragments typical for rosmarinic acid, which is characteristic for some Nepeta taxa (Mišić et al., 2015). This compound showed a neutral loss of 162 Da leading to the formation of a fragment ion at 359 m/z that corresponds to the exact mass of deprotonated rosmarinic acid. By studying other MS2 fragments, this compound was marked as rosmarinic acid hexoside.
All identified flavonoid compounds belong to the subgroup of flavones, which is typical for the genus Nepeta (Jamzad et al., 2003). The identified flavonoid aglycones were generally more abundant in inflorescences than in leaves. Hexuronides of luteolin (4 compounds) and apigenin (1 compound) were also recorded in the analyzed samples of N. nuda (Supplementary Table 2). Inflorescences were especially rich in these compounds. Flavones (mainly apigenin and luteolin derivatives) glycosylated with hexuronic (glucuronic) acid, often with another acyl residue (acetyl, coumaroyl, caffeoyl or similar), have been previously reported as constituents of several Nepeta species (Rabee et al., 2020). Compound 59 was identified as luteolin 7-O-(2”-caffeoyl)-hexuronide. In its MS2 spectra, ions were detected at 443 and 337 m/z, formed by the neutral loss of caffeic acid and luteolin, respectively (Supplementary Figure 1). This compound showed an unusual fragmentation pattern as its MS2 base peak was not the mass of luteolin aglycone (285 m/z) but corresponded to the mass of caffeic acid-H-H2O and it was previously identified in Satureja biflora aerial parts (Moghadam et al., 2015).
In the present study, we adopted the combination of analytical techniques previously reported to be convenient for establishing chemical characterization of N. nuda populations from the Central Balkans (Petrović et al., 2024). The mentioned study was focused on the chemical profiling of leaves dried in silica gel immediately after harvesting, while the samples used in this study were air-dried at room temperature until they reached a constant mass. This enabled us to compare and discuss the phytochemical composition of methanol extracts as influenced by drying procedure. Thus, vanillic acid was recorded only in silica gel-dried samples of N. nuda leaves (Petrović et al., 2024). Within the group of hydroxycinnamic acids, caffeic acid hexoside (13), ferulic acid (27), and nepetoidin A or B (28) were not present in air-dried leaves but were detected in inflorescences. The qualitative composition of flavonoids was also slightly perturbed, as luteolin (61) was recorded only in air-dried leaves, while chrysoeriol and acacetin were identified in leaves dried on silica-gel (Petrović et al., 2024).
3.2 Comparison of qualitative non-targeted iridoid profiles of N. nuda leaves and inflorescences
Twenty-seven iridoid compounds were identified in N. nuda inflorescences and leaves, 11 out of them being glycosides and 16 aglycones (Supplementary Table 2). One nitrogen analogue of nepetalactone, nepetalactam, was also recorded in N. nuda inflorescences. The presence of iridoid compounds was confirmed by the comparison with the available literature on Nepeta species and through the evaluation of their MS data. Some of the major iridoids identified in leaves and inflorescences of N. nuda are presented within Figure 1B.
Iridoid aglycones were mostly visible in the positive ionization mode due to their polarity, except for 7-deoxyloganetin, nepetalic and nepetonic acids (compounds 47, 48, and 50, respectively), which gave more abundant peaks in the negative ionization mode. Among the identified iridoid aglycones, 7-deoxyloganetin (47) and dihydronepetalactone (51) were recorded in inflorescences only, while cis,trans-nepetalactone (52) and nepetalactol acetate (53) were characteristic for N. nuda leaves.
Nepetalactones and their derivatives were the most abundant iridoid aglycones in the analyzed N. nuda samples. This iridoid monoterpenoid can usually be found in the form of four diastereoisomers, three of them, trans,trans-, cis,trans-, and trans,cis-nepetalactone, being identified in the present study. Previous studies dealing with the phytochemical characterization of N. nuda reported all 4 diastereoisomers: trans,trans- (De Pooter et al., 1988; Handjieva et al., 1996; Petrović et al., 2024), cis,trans- (De Pooter et al., 1988; Kökdil et al., 1999; Kobaisy et al., 2005; Mancini et al., 2009; Narimani et al., 2017; Akbaba et al., 2021; Petrović et al., 2024), trans,cis- (Regnier et al., 1967; De Pooter et al., 1988; Mancini et al., 2009; Gkinis et al., 2010; Bozok et al., 2017; Narimani et al., 2017; Akbaba et al., 2021; Aćimović et al., 2022), and cis,cis-nepetalactone (De Pooter et al., 1988; Kökdil et al., 1999; Kobaisy et al., 2005; Mancini et al., 2009; Bozari et al., 2013; Gormez et al., 2013; Zaharieva et al., 2023). In our study, trans,trans-nepetalactone was more abundant in inflorescences than in leaves (Figure 1A), while cis,trans-nepetalactone was recorded in leaves only. trans,cis-Nepetalactone was recorded in trace amounts in both leaves and inflorescences (Supplementary Table 2).
5,9-Dehydronepetalactone (40) is a dehydrogenated nepetalactone, isolated for the first time from N. cataria (Sastry et al., 1972). In N. nuda, this compound is primarily accumulated in inflorescences, while leaves contain considerably lower amounts (Figure 1A). This compound was previously reported in N. nuda (De Pooter et al., 1988; Kökdil et al., 1999; Petrović et al., 2024) and it was quantified in significant amounts in methanol extracts of N. rtanjensis, N. argolica, and N. parnassica (Mišić et al., 2015; Aničić et al., 2018, 2020, 2021). The content of 5,9-dehydronepetalactone in N. teydea EO was higher when it was prepared from flowering plants than from plants in the vegetative growth stage (Velasco-Negueruela et al., 1989), which is consistent with the findings of the present study.
Dihydronepetalactone (51) was identified in N. nuda inflorescences only (Supplementary Table 2). It is formed by nepetalactone hydrogenation (Feaster et al., 2009) and is known to be a common constituent of N. cataria (Regnier et al., 1967; Sharma and Cannoo, 2013; Sengupta et al., 2018; Handjeva et al., 1996; Patel et al., 2022).
Nepetalic acid (syn. 3α-hydroxy-4aα,7α,7aα-dihydronepetalactone, 3R-hydroxy-4aR,7R,7aR-dihydronepetalactone) can be synthetically converted into nepetalactone (Chauhan and Zhang, 2008), which marks this compound and plants rich in it as potentially valuable industrial resources. N. nuda is a good candidate as it accumulates nepetalic acid (48) in high amounts, especially in inflorescences (Supplementary Table 2, Figure 1A).
A derivative of nepetalic acid, its methyl isomer, nepetonic acid (2‚-[(1-methyl-2-al)-ethyl-5R-methylcyclopenta‚-carboxylic acid) (50), was identified in N. nuda samples in the present study (Supplementary Table 2). Methyl acetal of nepetalic acid (54) can be synthesized by nepetalic acid acylation (Scialdone and Liauw, 2008) and it was identified exclusively in inflorescences of N. nuda, adopting the UHPLC/QToF-MS technique (Supplementary Table 2).
Nepetalactam (81), classified under the group of tetrahydropyridine metabolites (lactams or cyclic amides), is a nitrogen analogue of (4aS,7S,7aR; cis,trans-)-nepetalactone previously identified in the essential oil of N. cataria (Handjeva et al., 1996; Patel et al., 2022), but also synthetized from nepetalic acid and nepetalactone (Eisenbraun et al., 1988; Chauhan et al., 2014). Interestingly, this compound was accumulated only in N. nuda inflorescences (Supplementary Table 2), indicating its possible contribution in plant-biotic interactions and reproductive strategy of this plant species. This compound was not identified in our previous study adopting a similar analytical approach (Petrović et al., 2024), most certainly since inflorescences were not chemically characterized.
Derivatives of iridoid glycosides were identified in the negative ionization mode (Supplementary Table 2), mostly showing [M–H]– as a molecular ion, while adducts with formic acid ([M + HCOOH–H]–) were observed in three iridoid glycosides (31, 33, and 34). The majority of identified compounds from this group were present in both inflorescences and leaves. The exception was caffeoyl-1,5,9-epi-deoxyloganic acid (39), which was identified in N. nuda inflorescences only.
1,5,9-epi-Deoxyloganic acid was previously identified in many Nepeta species (Murai et al., 1984; Nagy et al., 1998; Aničić et al., 2021, including N. nuda (Kökdil et al., 1999; Petrova et al., 2022; Petrović et al., 2024). In the present study, 1,5,9- epi-deoxyloganic acid was detected in both leaves and inflorescences, being more abundant in the latter. Its three derivatives were also recorded. 1,5,9-epi-Deoxyloganic acid hexoside was very abundant in leaves and especially in inflorescences, while methyl-1,5,9-epi-deoxyloganic acid was found in leaves of N. nuda. To the best of our knowledge, caffeoyl-1,5,9-epi-deoxyloganic acid was not previously recorded in any Nepeta species, and in the present study it was identified in significant amounts only in inflorescences of N. nuda.
The most abundant compound from the subgroup of iridoid glycosides, compared by peak area, was nepetanudoside A (compound 34). This compound, typical for N. nuda (Takeda et al., 1995; Petrović et al., 2024), was detected in both inflorescences and leaves but was more abundant in the reproductive organs (Figure 1A). The same trend was recorded for nepetanudoside C (31), which was previously reported for N. nuda (Takeda et al., 1996; Petrović et al., 2024).
Aucubin (29), present as a minor component in methanol extracts, was more abundant in inflorescences than in leaves of N. nuda and was previously reported for other Nepeta species, including N. septemcrenata (El-Moaty, 2010), N. rtanjensis, and N. argolica (Aničić et al., 2021). Lamiol derivative, 5-deoxylamiol (32), is an iridoid glycoside typical for Lamium species (Alipieva et al., 2003 and 2007) and has only recently been identified in N. nuda (Petrović et al., 2024). On the other hand, 6-deoxylamioside (38), which is here identified in both inflorescences and leaves, was previously reported for L. amplexicaule (Guiso and Martino, 1983).
To visualize the relations in metabolic fingerprints of N. nuda leaves and inflorescences, HCA was performed based on the Spearman’s method of cluster agglomeration, utilizing the relative quantification data (peak areas) (Figure 1A). These two plant organs are clearly diversified based on the composition of phenolics and iridoids, as the majority of the identified compounds are more abundant in inflorescences than in leaves. Compounds, on the other hand, are divided into two major clusters, based on their predominance in leaves or inflorescences.
3.3 Comparison of qualitative non-targeted profiles of volatile organic compounds in N. nuda leaves and inflorescences
In total, 36 VOCs were identified in methanol extracts of N. nuda inflorescences and leaves (Table 1).
The representative GC/MS TIC chromatograms are shown in Figure 1C. These VOCs belong to the groups of monoterpenes (6 monoterpenes and 13 oxygenated monoterpenoids) and sesquiterpenes (7 sesquiterpenes and 5 oxygenated sesquiterpenoids). Nepetalactam (36G), recorded only in inflorescences, was not classified in any of the two groups. Monoterpenoids were more abundant in leaves, while inflorescences were richer in oxygenated monoterpenoids (Table 1). Among monoterpenoids, β-pinene (3G) and sabinene (2G) predominated in N. nuda leaves. 1,8-Cineole (7G) was more abundant in leaves, although a significant amount of this compound was recorded in inflorescences, as well. From the group of monoterpenoid iridoids, two nepetalactone stereoisomers were detected: trans,trans- (18G) and cis,trans-nepetalactone (19G), the former predominating in inflorescences and the latter in leaves. Nepetonic acid (35G) was also more abundant in inflorescences. Interestingly, trans,cis- isomer of nepetalactone and some of nepetalactone derivatives (5,9-dehydronepetalactone, dihydronepetalactone, nepetalic acid), previously identified by the UHPLC-ESI-QToF-MS technique, were not recorded in methanol extracts of N. nuda samples using the GC/MS analysis. This can be, at least partially, ascribed to differential chromatography conditions (UHPLC vs. GC) and huge differences in the sensitivities of mass spectrometers.
Sesquiterpenoids and oxygenated sesquiterpenoids were more abundant in leaves than in inflorescences of N. nuda (Table 1). Germacrene D (25G) and beta-caryophyllene (21G) were detected as the major sesquiterpenoids in leaves, followed by β-bourbonene (20G). In inflorescences, germacrene D (25G) was the most abundant compound as well but was followed by β-farnesene (23G) and beta-caryophyllene (21G). Oxygenated diterpene phytol (33G) and oxygenated triterpene squalene (34G) were by far more abundant in leaves of N. nuda than in inflorescences. In sum, inflorescences and leaves of N. nuda were clearly distinguished based on the content of teprenoids and related compounds and can therefore be considered as different chemotypes.
Only scarce literature data reported on the composition of N. nuda methanol extracts revealed by GC/MS characterization. Other studies primarily focused on analyzing essential oils composition and reported the following compounds as dominant components: caryophyllene oxide (Kökdil et al., 1998), camphor and 1,8-cineole (Kilic et al., 2011), caryophyllene (Alim et al., 2009), nepetalactone (Bozok et al., 2017; Handjieva et al., 1996; Kökdil et al., 1996), 1,8-cineole and 4aα,7β,7aα-nepetalactone (Mamadalieva et al., 2016), nepetalactone and germacrene (Gormez et al., 2013), 1,8-cineole (Chalchat et al., 1998), and 1,8-cineole and a mixture of nepetalactones and germacrene-D (De Pooter et al., 1988). It must be considered that these differences might arise from many factors, including plant origin, growth conditions, developmental stage, plant part analyzed, a subspecies or a cultivar in question, essential oils isolation procedure, and GC/MS conditions. A study by Zaharieva et al. (2023) revealed 4aα,7β,7aα-nepetalactone, 1,8-cineole/eucalyptol, and germacrene D as the most abundant volatile compounds in N. nuda plants grown under either in vitro or ex vitro conditions. Caryophyllene, β-ocimene, bicyclogermacrene, β-pinene, myrcene, and humulene were also found in quite fair quantities.
The relative GC/MS quantitative data (peak areas) of methanol extracts of N. nuda were subjected to HCA to visualize the phytochemical relations between inflorescences and leaves (Figure 1C). Leaves cluster separately from inflorescences on the HCA plot constructed based on the Spearman’s algorithm, confirming the clear diversification between these plant parts. Compounds are clustered based on their predominance in inflorescences or leaves.
3.4 Organ-specific variation in the quantitative content of major iridoids and phenolics among the Central Balkans’ populations of N. nuda
The UHPLC/DAD/(±)HESI−MS2 analysis was targeted towards totally 9 compounds belonging to the groups of iridoids (5 compounds) and phenolics (4 compounds), and the selection was based on previous studies reporting high abundance of these compounds in N. nuda (Mišić et al., 2015; Aras et al., 2016; Smiljković et al., 2018; Petrova et al., 2022; Petrović et al., 2024). From the group of iridoids, we selected three aglycones (cis,trans- and trans,cis-nepetalactone, and 5,9-dehydronepetalactone) and two glycosides (1,5,9-epi-deoxyloganic acid and nepetanudoside A). The most abundant phenolics in N. nuda were rosmarinic acid and caffeic acid (Mišić et al., 2015; Petrova et al., 2022), and these two phenolic acids were quantified along with two flavonoids from the group of flavones, luteolin and apigenin (Figure 2). The majority of the analyzed metabolites were more abundant in inflorescences than in leaves of N. nuda. Literature data support these findings, as inflorescences have generally higher constitutive levels of defensive metabolites than leaves (Brown et al., 2003; Damle et al., 2005; Smallegange et al., 2007), with some exceptions (Godschalx et al., 2016). On the other hand, the amounts of these compounds were variable among the 13 analyzed populations.
Figure 2. UHPLC/DAD/(±)HESI−MS2 quantification of targeted iridoids and phenolics in methanol extracts of Nepeta nuda leaves (green bars) and inflorescences (blue bars), in samples originating from 13 Central Balkan populations. SE of three biological replicates is presented. Values labeled with different letters, for leaves (green letters) and inflorescences (blue letters) independently, are significantly different (p < 0.05) according to the Tukey’s post hoc test of one way ANOVA. Orange asterisks above the bars denote significantly different values between leaves and inflorescences according to the t-test, p-values, *p < 0.05, **p < 0.01, ***p < 0.001. Abbreviations of N. nuda populations: 1-Mali Ljukten, 2-Brodica, 3-Debeli lug, 4-Rtanj, 5-Straža, 6-Židilje, 7-Vinatovača, 8-Vlasina, 9-Donji Krivodol, 10-Topli Do, 11-Janjska reka, 12-Balta Berilovac, 13-Gornje selo.
The results display obvious differentiation between leaves and inflorescences (Figure 2). The major iridoid glycoside, nepetanudoside A, was present in both inflorescences and leaves of N. nuda, and in populations 1, 4, 5, 7, 8, 10, and 13, the content of this compound was significantly higher in leaves than in inflorescences. Population 5 was especially rich in this compound. 1,5,9-epi-Deoxyloganic acid was more abundant in inflorescences than in leaves of N. nuda, and this implies to all analyzed populations (Figure 2). The highest amounts of this compound were recorded in samples belonging to the populations 5 and 13. Of the iridoid aglycones, trans,trans-nepetalactone and 5,9-dehydronepetalactone were more abundant in flowers than in leaves of N. nuda. The content of trans,trans-nepetalactone varied between populations, with the population 4 and 5 being especially rich in this compound. The highest amounts of 5,9-dehydronepetalactone were recorded in samples belonging to the populations 2, 3, and 4. Interestingly, cis,trans-nepetalactone was recorded in significant amounts in only 5 populations (2, 3, 6, 7, and 8), primarily in leaves (Figure 2).
To visualize variations in the quantitative composition of the major iridoids in inflorescences and leaves of 13 N. nuda populations, PCA was performed adopting the Ward’s method of data agglomeration (Figure 3A). The major contributors to the diversification of samples along the PC1 explaining 75.76% of the total variability are nepetanudoside A, cis,trans-nepetalactone, and 1,5,9-epi-deoxyloganic acid. As for the PC2, which contributes to the total variability with 17.43%, cis,trans-nepetalactone and 1,5,9-epi-deoxyloganic acid are the major contributors to the diversification along this component. Samples of N. nuda leaves are more scattered along PC1 than samples of inflorescences, implying higher uniformity in the accumulation of the targeted compounds by inflorescences.
Figure 3. Principal component analysis (PCA) biplots constructed based on the UHPLC/DAD/(±)HESI−MS2 quantitative data, separately for iridoids (A) and phenolics (B), and cumulatively for both groups of metabolites (C). Participation of the variables (metabolites) in the first two PCs is indicated by the corresponding loading plots. Inflorescences are labeled with blue symbols while leaves are colored green.
Among phenolics, rosmarinic acid was by far the most abundant compound in the analyzed samples. It was recorded in considerable amounts in both leaves and inflorescences of N. nuda. Its amounts in inflorescences reached ∼802 µg 100 mg-1 dry weight (DW) (population 6), while in leaves ∼678 µg 100 mg-1 DW was recorded in samples from the population 11 (Figure 2). In the majority of the analyzed populations, the amounts of this compound were slightly higher in inflorescences than in leaves. Caffeic acid, luteolin, and apigenin were also more abundant in inflorescences than in leaves in all analyzed populations. Population 9 was characterized by the highest amounts of luteolin (6.50 µg 100 mg-1 DW) and apigenin (7.19 µg 100 mg-1 DW), while the highest amounts of caffeic acid were recorded for the population 11 (Figure 2).
When PCA was performed with quantitative data for phenolics only, it was obvious that samples of leaves and inflorescences were clearly diversified along PC2 (Figure 3B). The major contributors to the diversification along PC1, which explains 99.98% of the total variability, and along PC2 (0.02%), are apigenin and luteolin. The third contributor is caffeic acid. Results imply that phenolics content is more conserved in leaves than in inflorescences.
PCA plot was further constructed based on the quantitative data for both iridoids and phenolics (Figure 3C). PC1 explains 70.52% of the total variability, while PC2 is responsible for 15.82% of the variability. The major contributors to the diversification along PC1 are nepetanudoside A, rosmarinic acid, and cis,trans-nepetalactone. The same compounds, but with perturbed order, are the major factors that diversify samples along PC2.
3.5 Molecular background of iridoid biosynthesis in N. nuda inflorescences and leaves
We further aimed to get deeper into the molecular background of the organ-specific iridoid diversity in N. nuda. To exclude the effect of genotype and environmental conditions on the chemical phenotype, we performed clonal propagation of a single genotype originating from the population 3 (abbreviation DL1). Plantlets having the uniform genetic background were acclimatized and grown under controlled greenhouse conditions. Samples of leaves and inflorescences were collected from fully-flowering plants, in three biological repetitions.
Targeted metabolic profiling of five iridoids (cis,trans-nepetalactol, trans,trans-nepetalactone, 5,9-dehydronepetalactone, 1,5,9-epi-deoxyloganic acid, and nepetanudoside A) in leaves and inflorescences of greenhouse-grown N. nuda in parallel with co-expression analysis of iridoid biosynthesis-related genes, revealed the molecular background of their organ-specific composition. As expected, all the quantified iridoids were more abundant in inflorescences than in the leaves of N. nuda plants (Figure 4A). The difference in the abundance of specific iridoids between the two plant organs was the most pronounced for 1,5,9-epi-deoxyloganic acid and cis,trans-nepetalactol, the latter being recorded only in inflorescences. The correlation analysis based on the quantitative data revealed significant positive correlations for the majority of comparisons (Figure 4B). These results are in accordance with the study of Gomes et al. (2024), which associated the highest concentrations of nepetalactones, nepetalic acid and nepetalactam with the floral-bud and partial-flowering stages of N. cataria.
Figure 4. (A) Box-plot diagram depicting the distribution of the content of three iridoid aglycones (c,t-NLL, t,t-NL, 5,9-DNL) and two iridoid glycosides (1,5,9-eDLA and NnuA) in leaves (green dots) and inflorescences (blue dots) of Nepeta nuda grown under greenhouse conditions. (B) Pearson’s correlations based on the content of five major iridoids in N. nuda leaves and inflorescences represented by the intensity of yellow (positive correlations) or blue (negative correlations), as indicated on the color scale. (C) A proposed iridoid metabolic pathway in N. nuda. Unknown enzymes of the biosynthetic pathway are marked with the question mark. (D) Box-plot diagram depicting the distribution of the transcript levels of nine iridoid biosynthesis-related genes in leaves (green dots) and inflorescences (blue dots) of N. nuda. (E) Pearson correlations based on the expression level of nine biosynthetic genes in N. nuda leaves and inflorescences. Non-significantly correlated values are crossed (p > 0.05). IPP, isopentenyl pyrophosphate; DMAPP, dimethylallyl pirophosphate; GPP, geranyl pyrophosphate; GE, geraniol; 8HG, 8-hydroxygeraniol; 8OG, 8-oxogeranial; 8CE, 8-oxocitronellyl enolate; NLL, nepetalactol; NL, nepetalactone; DNL, dehydronepetalactone; c,t-NLL, cis,trans-nepetalactol; t,t-NL, trans,trans-nepetalactone; 5,9-DNL, 5,9-dehydronepetalactone; 1,5,9-eDLA, 1,5,9-epideoxyloganic acid; NnuA, nepetanudoside A; NnGPPS, geranyl diphosphate synthase; NnGES, geraniol synthase; NnG8H, geraniol 8-hydroxylase; Nn8HGO, 8-hydroxygeraniol oxidoreductase; NnISY, iridoid synthase; NnNEPS, nepetalactol-related short-chain dehydrogenase; NnMLPL, major latex protein-like enzyme. Asterisks denote significantly different values according to the t test (*p < 0.05, **p < 0.01, ***p < 0.001).
The iridoid biosynthetic pathway in Nepeta species has been extensively researched in the past, and it is believed that nepetalactol is the common precursor of all other iridoids in these species. The biosynthesis of nepetalactol involves a series of enzymatic steps catalyzed by geranyl diphosphate synthase (GPPS), geraniol synthase (GES), geraniol 8-hydroxylase (G8H), 8-hydroxygeraniol oxidoreductase (8HGO), iridoid synthase (ISY), NAD-dependent nepetalactol-related short-chain-dehydrogenase/reductase (NEPS), and major latex protein-like (MLPL) enzymes (Figure 4C). In Nepeta species, ISYs are mainly responsible for the stereoselective 1,4-reduction of 8-oxogeranial (8OG) to uncyclized and reactive 8-oxocitronellyl enol (8CE) and for determining the stereochemistry of C7 (Sherden et al., 2018; Lichman et al., 2019b and 2020; Hernaíndez Lozada et al., 2022). Enolate intermediate 8-oxocitronellyl enol in Nepeta is cyclized by NEPSs (Lichman et al., 2019a, 2019), but can also undergo a spontaneous cyclization to produce predominately cis,trans- steroisomer of nepetalactol (Hernaíndez Lozada et al., 2022). Although the upstream parts of the pathway have been well described, the downstream parts that are leading to various derivatives of nepetalactone, as well as those of the iridoid glycosides branch, remain largely unclear.
By BLAST-searching for the available N. nuda transcriptomes, candidates for 9 iridoid biosynthesis-related genes (NnGPPS, NnGES, NnG8H, Nn8HGO, NnISY, NnNEPS1, NnNEPS2, NnNEPS3, and NnMLPL) were identified and their expression profiles were comparatively analyzed in N. nuda leaves and inflorensces to more deeply investigate the molecular background of their variation between the two organs and to point to the key determinants of the organ-specific iridoid metabolism in N. nuda. The expression levels of NnGPPS and Nn8HGO, two genes of the early biosynthetic pathway (EBGs), did not significantly differ between leaves and inflorescences (Figure 4D). However, NnGES, NnG8H, and NnISY displayed around 1.6-, 6-, and 2.2-fold higher transcript levels in inflorescences than in leaves, respectively. Interestingly, genes responsible for the steps directly preceding the formation of nepetalactol and nepetalactone (NnNEPS1, NnNEPS2, NnNEPS3, and NnMLPL) showed significantly higher expression level in inflorescences than in leaves (Figure 4D). In the case of NnNEPS1, which belongs to the subgroup of NEPS oxidases and is responsible for converting nepetalactol into nepetalactone, around 40-fold higher expression level was recorded in inflorescences than in leaves. The NEPS oxidase (NnNEPS1) identified in the transcriptomes of N. nuda leaves, that corresponded to the N. mussinii NEPS1, has been reported to have cis,trans-nepetalactol dehydrogenase activity (Lichman et al., 2019a). NnNEPS2 is a homologue of N. sibirica NEPS2 (NsNEPS2), which acts both as a cis-trans cyclase and an oxidase (Lichman et al., 2019a; Hernaíndez Lozada et al., 2022). NnNEPS3 is a homologue of N. mussinii NEPS3 displaying a specific 7S-cis,cis-nepetalactone cyclase activity (Lichman et al., 2019a). The composition, expression level, and ratio of NEPS enzymes are responsible for setting the stereochemistry of nepetalactones (Lichman et al., 2019a), which is most certainly the case in N. nuda as well. The results indicate that NnG8H and the 3 NEPSs might serve as the regulatory points, determining the flux through the iridoid pathway and thus the level of iridoid biosynthesis and accumulation in different organs. Hence, the high level of nepetalactone and other iridoids in inflorescences is, most likely, the result of upregulated expression of NnG8H, NnNEPS1, NnNEPS2, and NnNEPS3 in these organs. The majority of correlations observed between the analysed biosynthetic genes were positive and statistically significant, which suggests the regulation of iridoid biosynthesis in different organs at the transcriptional level (Figure 4E). The exception were NnGPPS and Nn8HGO, which displayed no significant correlations with other genes. Conversely, the expression of NnGES was significantly positively correlated with that of NnG8H, NnNEPS2, and NnNEPS3.
Some previous studies dealing with Nepeta species revealed organ-specific expression patterns of the iridoid biosynthesis-related genes (Aničić et al., 2018, and 2024, Lichman et al., 2020). Closed buds and open flowers of N. cataria and N. mussinii were more enriched with transcripts of NEPSs, G8H, and ISY than mature leaves (Lichman et al., 2020). Among NEPS enzymes, NEPS1 displayed the highest expression level in open flowers and closed buds of N. cataria, while NEPS5 transcripts were the most abundant in closed buds of N. mussinnii (Lichman et al., 2020). This is in accordance with the findings of the present study, as NnNEPS1 displayed a significantly higher expression level in inflorescences than in mature leaves of N. nuda and more transcript levels than the other two NEPSs (NnNEPS2 and NnNEPS3). Although we did not separately assessed closed buds from open flowers, the analyzed inflorescences contained flowers in all development stages.
The expression of iridoid biosynthesis-related genes is coordinately regulated and well correlated with metabolite pools in mature leaves and inflorescences (Table 2), supporting the hypothesis that the organ-specific biosynthesis of iridoids is controlled at the transcriptional level. In other words, gene transcripts of targeted genes were more abundant in inflorescences, where the majority of targeted iridoids accumulated. The exceptions are NnGPPS, NnGES, and Nn8HGO, as their expression levels were not significantly correlated with the content of targeted metabolites. The majority of significant correlations were observed between the three NEPSs and G8H on one hand and iridoid aglycones (5,9-dehydronepetalactone, cis,trans-nepetalactol, and trans,trans-nepetalactone) on the other. The most significant correlations were observed between the expression of NnNEPS1 and NnNEPS3 and the three iridoid aglycones. NnNEPS1 was also positively correlated with the content of quantified iridoid aglycones, nepetanudoside A and 1,5,9-epi-deoxyloganic acid (Table 2), indicating possible involvement of this NEPS oxidase in the formation of iridoid glycosides. These speculations should be further investigated.
Table 2. Correlations (R2) between quantitative metabolomics data (for 5 selected iridoid compounds) and qPCR expression data for 9 iridoid biosynthesis-related genes in the analyzed samples of Nepeta nuda (n= 324).
3.6 Possible ecological implications of iridoids’ composition in N. nuda inflorescences and leaves
Due to a variety of important roles that iridoids play in plant-biotic interactions, we here discuss possible ecological implications of their organ-specific composition in N. nuda. Data show that N. nuda leaves exhibit greater diversity in terms of the targeted phenolics’ and iridoids’ content, and that phytochemical composition of inflorescences is more conserved among the analyzed populations. Plants likely invest more resources in the constitutive production of floral metabolites functioning in plant-biotic interactions (i.e., plant-pollinator, plant-pathogen, and plant-herbivore interplay), rather than in inducible production. It was previously hypothesized that constitutive resistance predominates in flowers (McCall and Irwin, 2006; Zangerl, 2003). Floral chemistry influences the interactions of plants with pollinators and other species and maximizes the fitness of both plants and flower visitors (Palmer-Young et al., 2019; Kessler and Kalske, 2018). The defense mechanisms of reproductive tissues, such as flowers, are crucial for plant reproduction, and their defensive strategies may differ from those of leaves, which have the primary function of fixing CO2 during photosynthesis. According to Pichersky and Gershenzon (2002), floral volatiles serve as attractants for species-specific pollinators, whereas the volatiles emitted from vegetative parts, especially those released after herbivory, appear to protect plants by deterring herbivores and by attracting their enemies. Flowers are essentially heterotrophic as they heavily rely on organic molecules synthesized in leaves and roots. However, they biosynthesize and accumulate complex molecules that sustain plant interactions with animal pollinators and pathogens, play roles in plant fertilization, embryo development, and the initial stages of fruit and seed set.
The chemical composition of colorful long-lasting flower spikes of N. nuda is a key determinant of their color, odor, and overall function in plant-biotic interactions. Nepeta flowers with their visual, chemical, and nutritional traits are very attractive to honeybees and other pollinators. Honeybees are the main pollinators of N. nuda (Bozek, 2003). The most important visitors of flowers of N. cataria, a species sympatric to N. nuda, are honeybees, solitary bees, and bumblebees (Cruden and Miller-Ward, 1981; Sih and Baltus, 1987), but also flies and butterflies (Aizen and Harder, 2009), thus it may be presumed that N. nuda is visited by similar groups of pollinators, as these two species share similar phytochemical profiles.
The increased content of iridoids in the reproductive organs of N. nuda might indicate their functional involvement in complex plant-biotic interactions. It is well documented that iridoids function as general toxicants, growth and reproductive inhibitors, repellents or oviposition-deterrents (Bernier et al., 2005; Birkett et al., 2011; Sparks et al., 2018; Reichert et al., 2019; reviewed in Formisano et al., 2011; Süntar et al., 2018). Two groups of iridoids present in Nepeta plants, volatile aglycones and glycosylated iridoids, most likely provide a two-level chemical defense, the former acting as repellents and/or pollinator attractants, and the latter contributing as deterrents. Both subgroups of iridoids also provide antimicrobial protection.
Constituents of Nepeta species, such as nepetalactone, dihydronepetalactone, and nepetalactol, are proven repellents against insects, such as mosquitos (Aedes, Anopheles, Culex), cockroaches (Periplaneta, Blattella), ticks (Rhipicephalus, Ixodes), flies (Stomoxys, Musca), mites (Dermanyssus, Dermatophagoides), and termites (Reticulitermes) (Bernier et al., 2005; Birkett et al., 2011; reviewed in Formisano et al., 2011; Sparks et al., 2018; reviewed in Süntar et al., 2018; Reichert et al., 2019). Dihydronepetalactone is reported to be two times more active than commercial N,N-diethyl-3-methylbenzamide (DEET) (Feaster et al., 2009). Notably, nepetalactone and nepetalactol mediate plant-insect interactions, but are also produced by a number of insects, most notably aphids, acting as sex pheromones (Fernández-Grandon and Poppy, 2015). In addition, dihydronepetalactones are the components of defensive secretions of some ant species (Cavill and Clark, 1967). As reported by Harney et al. (1978), nepetalic acid, along with nepetalactone and catnip essential oil, had a depressant effect in mice when administered intraperitoneally.
On the other hand, compounds from the group of iridoid glycosides present in Nepeta species, less volatile and more stable than nepetalactones and their derivatives, have rarely been investigated for their interaction with insects and generally for their contribution to the overall defense strategies. 1,5,9-epi-Deoxyloganic acid has been scarcely investigated for its biological activities, thus, only studies dealing with its antimicrobial and anti-inflammatory potential (Aničić et al., 2021), as well as enzyme inhibition and DNA protection activities (Başar et al., 2023) are presented in the literature. To the best of our knowledge, no data about its role in plant-biotic interactions are available. Despite only scarce information about the bioactivities of the major iridoid glycosides of N. nuda, plenty is known about some iridoid glycosides present as minor constituents. Such an example is aucubin, which displays potent insecticidal activity against various insects (Vivekanandhan et al., 2022, 2023) and antifungal activity (Marak et al., 2002). Comprehensive literature search revealed that the roles of nepetenudoside A and C, as well as of nepetaside, 5-deoxylamiol, and 6-deoxylamioside in plant-biotic interactions have not been previously documented. Due to their physico-chemical properties, iridoid glycosides of Nepeta sp. might provide additional values in the development of novel biopesticides.
The composition of floral volatiles, a chemically diverse group of plant metabolites with multiple functions, is shaped by environmental, ecological and evolutionary factors (Dötterl and Gershenzon, 2023). Although iridoids are by far the most abundant group of terpenes in Nepeta species, their synergistic effects with other terpenes must not be ruled out. The unique metabolome of each individual plant constitutes the defense machinery against biotic environmental factors and determines its success in the changing environment. Several reports showed insecticidal and repellent activity of eucalyptol, which is abundant in both leaves and inflorescences of N. nuda. Eucalyptol induces intoxication and has a prominent effect on locomotor activity, knock-down, and repellence on nymphs of Triatoma infestans and Rhodnius prolixus (Moretti et al., 2015). It also displays repellence against a number of insects, including house fly (Sukontason et al., 2004), Periplaneta americana (Scriven and Meloan, 1984), and Aedes aegypti (Klocke et al., 1987). On the other hand, eucalyptol is attractive to honey bees and euglosine bees (Effmert et al., 2005).
It was previously reported that the content and ratio of sesquiterpenoids β-farnesene and germacrene D in flowers of Tanacetum vulgare modulated both defense and pollination (Li et al., 2019). Aphid alarm pheromone β-farnesene is the key volatile cue responsible for specific attraction of beetles and avoidance by aphids at the early flowering stages (Li et al., 2019). As these terpenoids are also found in N. nuda inflorescences and leaves, it could be hypothesized that they can play important roles in plant-biotic interactions.
3.7 Molecular background of the biosynthesis of phenolics in N. nuda inflorescences and leaves
As expected, rosmarinic acid was by far the most abundant phenolic compound in leaves and inflorescences of N. nuda grown under greenhouse conditions (Figures 5A, B). Other phenolics quantified in samples were present in much lower amounts. All analyzed compounds were found in significantly higher amounts in inflorescences than in leaves. The amount of rosmarinic acid was significantly negatively correlated with those of caffeic acid, luteolin, and apigenin (Figure 5C). The only non-significant correlation was observed between luteolin and apigenin.
Figure 5. (A) Box-plot diagram depicting the distribution of the content of one phenolic acid (CA) and two flavonoids (Ap and Lu) in leaves (green dots) and inflorescences (blue dots) of N. nuda grown under greenhouse conditions. (B) RA is presented by a separate box-plot diagram. (C) Pearson’s correlations based on the content of four major phenolics in N. nuda leaves and inflorescences, represented by the intensity of yellow (positive correlations) or blue (negative correlations), as indicated on the color scale. Non-significantly correlated values are crossed (p > 0.05). (D) A proposed biosynthetic pathway of rosmarinic acid in N. nuda. (E) Box-plot diagram depicting the distribution of the transcript levels of seven rosmarinic acid biosynthesis-related genes in N. nuda leaves and inflorescences, represented by the intensity of yellow (positive correlations) or blue (negative correlations), as indicated on the color scale. (F) Pearson’s correlations based on the expression level of seven genes of the rosmarinic acid biosynthetic pathway in N. nuda leaves (green dots) and inflorescences (blue dots). Non-significantly correlated values are crossed (p > 0.05). CA, caffeic acid; RA, rosmarinic acid; Ap, apigenin; Lu, luteolin; PAL, phenylalanine ammonia-lyase; C4H, cinnamic acid 4-hydroxylase; 4CL, 4-coumaric acid CoA-ligase; TAT, tyrosine aminotransferase; HPPR, hydroxyphenylpyruvate reductase; RAS, 4-coumaroyl-CoA:4′-hydroxyphenyllactic acid 4-coumaroyltransferase; CYP, cytochrome P450-dependent monooxygenase. Asterisks denote significantly different values according to the t test (*p < 0.05, ***p < 0.001).
The biosynthesis of rosmarinic acid (Figure 5D) involves phenylpropanoid- and tyrosine-derived branches that use the amino acids L-phenylalanine and L-tyrosine as precursors, respectively (Matsuno et al., 2002; Petersen et al., 2009; Wang et al., 2017; Trócsányi et al., 2020; Suthar et al., 2021) and it proceeds by activating the following enzymes: phenylalanine ammonia-lyase (PAL), cinnamic acid 4-hydroxylase (C4H), 4-coumaric acid CoA-ligase (4CL), tyrosine aminotransferase (TAT), hydroxyphenylpyruvate reductase (HPPR), “rosmarinic acid synthase” (RAS; 4-coumaroyl-CoA:4′-hydroxyphenyllactic acid 4-coumaroyltransferase), and cytochrome P450-dependent monooxygenase. We here analyzed the expression profiles of 7 genes related to the biosynthesis of rosmarinic acid (NnPAL, NnC4H, Nn4CL, NnTAT, NnHPPR, NnRAS, and NnCYP) in N. nuda leaves and inflorescences. Looking at their enzymatic activities involved in this biosynthetic pathway, only two of them might be specific for rosmarinic acid biosynthesis, NnRAS and NnCYP. These two genes showed significantly higher expression levels in inflorescences than in leaves (Figure 5E), which is in accordance with a higher content of rosmarinic acid in these organs. Enzymes belonging to the general phenylpropanoid pathway (NnPAL, NnC4H, and Nn4CL) provide precursors for the formation of other metabolites such as lignin and other groups of phenolic compounds (e.g., flavonoids, anthocyanins). Transcript levels of NnC4H and Nn4CL were higher in inflorescences than in leaves, while NnPAL did not significantly differ between the two organs (Figure 5E). TAT and HPPR are also primary enzymes, and their activity results in the formation of precursors in the biosynthesis of tocopherols and plastoquinones. These two genes did not significantly differentially express between leaves and inflorescences of N. nuda (Figure 5E). The only significant correlations were positive correlations between the expression patterns of the analyzed rosmarinic acid-related biosynthetic genes: Nn4CL and NnC4H, NnCYP and Nn4CL, and NnCYP and NnRAS (Figure 5F).
The results further suggest that the expression levels of NnPAL, NnC4H, and Nn4CL, which belong to the general phenylpropanoid pathway, as well as of NnRAS and NnCYP, are not significantly correlated with the rosmarinic acid content in the two plant organs (Table 3). This is not surprising, taking into account that the majority of analyzed genes are related to the general phenylpropanoid pathway and are thus involved in providing the precursors for the biosynthesis of various classes of phenolic compounds. Furthermore, the differences in the quantitative content of the major phenolic compounds between leaves and inflorescences are less prominent than those recorded for iridoids.
Table 3. Correlations (R2) between quantitative metabolomics data (for 4 selected phenolic compounds) and qPCR expression data for 7 rosmarinic acid biosynthesis-related genes in mature leaves and inflorescences of Nepeta nuda (n= 108).
3.8 Possible ecological implications of differential phenolics’ composition in N. nuda inflorescences and leaves
The involvement of phenolic compounds in plants ranges from growth and development (e.g., cell wall thickening, hormone production, and pigmentation), reproduction (e.g., pigmentation, fruit flavoring, and fruit protection), defense against stressors (e.g., osmoregulation, UV protection, anti-herbivory roles, and antimicrobial activity) (Dixon, 2001; Singh et al., 2021), and leaf and petal senescence (Tyagi et al., 2020). The content of phenolics in plants is determined by various endogenous and environmental factors. Thus, high constitutive and inducible leaf flavonoid content has been correlated with insect and pathogen resistance (Treutter, 2005). The function of phenolics in flowers is particularly versatile: from antioxidant activity, through the protection from UV radiation, to providing diverse colors. Phenolics accumulated in flowers primarily protect from nectar robbers and enhance memory and foraging efficiency of pollinators or, as in the case of flavonoids, exhibit high antioxidant activities supporting pollinator well-being (Slavković and Bendahmane, 2023). Moreover, flavonoids can regulate biotic interactions with mutualists and antagonists and thus reduce herbivore attacks and infection (Cushnie and Lamb, 2005; Palmer-Young et al., 2019).
Rosmarinic acid, which is by far the most abundant compound in N. nuda, is known for its antioxidant (Erkan et al., 2008; Muñoz-Muñoz et al., 2013), antimicrobial (Bais et al., 2002; Ivanov et al., 2022), insecticidal activities (Khan et al., 2019; Darrag et al., 2022), and various health promoting effects (Amoah et al., 2016; Alagawany et al., 2017; Nadeem et al., 2019). Rosmarinic acid displays insecticidal effects against the pea aphid, Acyrthosiphon pisum (Khan et al., 2019), and is a highly potent feeding deterrent for the tobacco hornworm (Manduca sexta) (Simmonds et al., 2019). Salvigenin, nepetoidin B, and rosmarinic acid display respectable insecticidal activities against the red palm weevil (Rhynchophorus ferrugineus) (Darrag et al., 2022). Conversely, this compound has been reported to reduce honey bee mortality caused by pesticides (Ulziibayar et al., 2022).
4 Conclusions
Inflorescences and leaves of Nepeta nuda showed a differential investment in the constitutive production of phenolics and terpenes, based on their chemical profiles, which emphasizes the importance of the chemical profiling of different plant organs and encompassing both groups of metabolites to properly characterize the species’ metabolome. The data strongly indicate that plants prioritize constitutive production of floral metabolites, whose qualitative and quantitative content is, to a large extent, conserved among different populations. NnG8H and NnNEPS1-3, are recognized as important regulatory points responsible for the increased content of iridoids in inflorescences.
The presented results lead to the unambiguous conclusion that specialized metabolism in leaves is more reprogrammable in response to differential growth conditions than that in inflorescences. Therefore, vegetative parts are highlighted as more convenient objects for studying the environmentally-driven chemical diversity among N. nuda populations. However, the specific mechanisms and the exact environmental cues associated with the organ-specific chemical composition in N. nuda remain largely unexplored. Due to its rich metabolite profile and a potential ornamental value, N. nuda can be a suitable alternative to N. cataria, which is the most intensively cultivated Nepeta species. The information about the intra-individual (organ-specific) chemical specificities of N. nuda provides strategic tools for growers to manipulate the phenology of this medicinal crop and achieve a higher productivity of metabolites of interest. Given the importance of iridoids as powerful insect repellents, future studies should be conducted towards optimization of agronomic practices to maximize yields in versatile geographic regions and to facilitate the integration of iridoid-containing biopesticides into pest management systems and other applications related to human health.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
LP: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – original draft. BF: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. MS: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Writing – original draft. BŠ: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – original draft. TB: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – review & editing. DMa: Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Writing – review & editing. JNŽ: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – review & editing. SD: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – review & editing. NA: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – review & editing. MM: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – review & editing. JB: Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing. UG: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Visualization, Writing – original draft. DMi: Conceptualization, Formal analysis, Funding acquisition, Investigation, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The research was financed by the Science Fund of the Republic of Serbia Grant No. 7749433, project acronym NEPETOME, and is also supported by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia [Contracts No. 451-03-47/2023-01/200007 and 451-03-66/2024-03/200007.
Acknowledgments
The authors would like to acknowledge David L. Hallahan, Ph.D. and Entomol Products LLC (GA, USA) for providing the standard of cis,trans-nepetalactone.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1452804/full#supplementary-material
References
Aćimović, M., Stanković-Jeremić, J., Cvetković, M. (2020). Phyto-pharmacological aspects of Nepeta nuda L.: A systematic review. Lek sirov. 40, 75–83. doi: 10.5937/leksir2040075A
Aćimović, M., Lončar, B., Pezo, M., Stanković Jeremić, J., Cvetković, M., Rat, M., et al. (2022). Volatile compounds of Nepeta nuda L. Rtanj Mountain (Serbia). Horticulturae. 8, 85. doi: 10.3390/horticulturae8020085
Aizen, M. A., Harder, L. D. (2009). The global stock of domesticated honey bees is growing slower than agricultural demand for pollination. Curr. Biol. 19, 915–918. doi: 10.1016/j.cub.2009.03.071
Akbaba, E., Bagci, E., Hritcu, L., Maniu, C. (2021). Improvement of memory deficits via acetylcholinesterase inhibitory activity of Nepeta nuda ssp. nuda essential oil in rats. Kuwait J. Sci. 24, 48(3). doi: 10.48129/kjs.v48i3.9217
Alagawany, M., Abd El-Hack, M. E., Farag, M. R., Gopi, M., Karthik, K., Malik, Y. S., et al. (2017). Rosmarinic acid: modes of action, medicinal values and health benefits. Anim. Health Res. Rev. 18, 167–176. doi: 10.1017/S1466252317000081
Alim, A., Goze, I., Cetin, A., Atas, A. D., Cetinus, S. A., Vural, N. (2009). Chemical composition and in vitro antimicrobial and antioxidant activities of the essential oil of Nepeta nuda L. subsp. Albiflora (Boiss.) gams. Afr J. Microbiol. Res. 3, 463–467. https://hdl.handle.net/20.500.12418/10083
Alimpić Aradski, A., Oalđe Pavlović, M., Janošević, D., Todorović, S., Gašić, U., Mišić, D., et al. (2023). Leaves micromorphology, chemical profile, and bioactivity of in vitro-propagated Nepeta cyrenaica (Lamiaceae). Phytochem. Anal. 34, 661–679. doi: 10.1002/pca.3257
Alipieva, K., Kokubun, T., Taskova, R., Evstatieva, L., Handjieva, N. (2007). LC–ESI-MS analysis of iridoid glucosides in Lamium species. Biochem. Syst. Ecol. 35, 17–22. doi: 10.1016/j.bse.2006.07.004
Alipieva, K. I., Taskova, R. M., Evstatieva, L. N., Handjieva, N. V., Popov, S. S. (2003). Benzoxazinoids and iridoid glucosides from four Lamium species. Phytochem. 64, 1413–1417. doi: 10.1016/j.phytochem.2003.08.001
Amoah, S., Sandjo, L., Kratz, J., Biavatti, M. (2016). Rosmarinic acid–pharmaceutical and clinical aspects. Planta medica. 82, 388–406. doi: 10.1055/s-0035-1568274
Aničić, N., Gašić, U., Lu, F., Ćirić, A., Ivanov, M., Jevtić, B., et al. (2021). Antimicrobial and immunomodulating activities of two endemic nepeta species and their major iridoids isolated from natural sources. Pharmaceuticals. 14, 414. doi: 10.3390/ph14050414
Aničić, N., Matekalo, D., Skorić, M., Gašić, U., Nestorović, Ž.J., Dmitrović, S., et al. (2024). Functional iridoid synthases from iridoid producing and non-producing Nepeta species (subfam. Nepetoidae, fam. Lamiaceae). Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1211453
Aničić, N., Matekalo, D., Skorić, M., Pećinar, I., Brkušanin, M., Živković, J. N., et al. (2018). Trichome-specific and developmentally regulated biosynthesis of nepetalactones in leaves of cultivated Nepeta rtanjensis plants. Ind. Crops Prod. 117, 347–358. doi: 10.1016/j.indcrop.2018.03.019
Aničić, N., Matekalo, D., Skorić, M., Živković, J. N., Petrović, L., Dragićević, M., et al. (2020). Alterations in nepetalactone metabolism during polyethylene glycol (PEG)-induced dehydration stress in two Nepeta species. Phytochem. 174, 112340. doi: 10.1016/j.phytochem.2020.112340
Aras, A., Bursal, E., Dogru, M. (2016). UHPLC-ESI-MS/MS analyses for quantification of phenolic compounds of Nepeta nuda subsp. Lydiae. J. App Pharm. Sci. 6, 009–013. doi: 10.7324/JAPS.2016.601102
Bais, H. P., Walker, T. S., Schweizer, H. P., Vivanco, J. M. (2002). Root specific elicitation and antimicrobial activity of rosmarinic acid in hairy root cultures of Ocimum basilicum. Plant Physiol. Biochem. 40, 983–995. doi: 10.1016/S0981-9428(02)01460-2
Başar, Y., Yenigün, S., İpek, Y., Behçet, L., Gül, F., Özen, T., et al. (2023). DNA protection, molecular docking, enzyme inhibition and enzyme kinetic studies of 1, 5, 9-epideoxyloganic acid isolated from Nepeta aristata with bio-guided fractionation. J. Biomol. Struct. Dyn. 42, 1–14. doi: 10.1080/07391102.2023.2250461
Bernier, U. R., Furman, K. D., Kline, D. L., Allan, S. A., Barnard, D. R. (2005). Comparison of contact and spatial repellency of catnip oil and N, N-diethyl-3-methylbenzamide (DEET) against mosquitoes. J. Med. Entomol. 42, 306–311. doi: 10.1603/0022-2585(2005)042[0306:cocasr]2.0.co;2
Birkett, M. A., Hassanali, A., Hoglund, S., Pettersson, J., Pickett, J. A. (2011). Repellent activity of catmint, Nepeta cataria, and iridoid nepetalactone isomers against Afro-tropical mosquitoes, ixodid ticks and red poultry mites. Phytochem. 72, 109–114. doi: 10.1016/j.phytochem.2010.09.016
Birkett, M. A., Pickett, J. A. (2003). Aphid sex pheromones: from discovery to commercial production. Phytochemistry 62, 651–656. doi: 10.1016/s0031-9422(02)00568-x
Bozari, S., Agar, G., Aksakal, O., Erturk, F. A., Yanmis, D. (2013). Determination of chemical composition and genotoxic effects of essential oil obtained from Nepeta nuda on Zea mays seedlings. Toxicol. Ind. Health 29, 339–348. doi: 10.1177/0748233711433939
Bozek, M. (2003). Pollen efficiency and foraging by insect pollinators in three catnip [Nepeta L.] species. J. Apic Sci. 47, 2. bwmeta1.element.agro-article-c15f86b6-2561-4147-b5ac-f924510b96c8.
Bozok, F., Cenet, M., Sezer, G., Ulukanli, Z. (2017). Essential oil and bioherbicidal potential of the aerial parts of Nepeta nuda subsp. albiflora (Lamiaceae). J. Essent. Oil-Bear. Plants. 20, 148–154. doi: 10.1080/0972060X.2016.1264279
Brown, P. D., Tokuhisa, J. G., Reichelt, M., Gershenzon, J. (2003). Variation of glucosinolate accumulation among different organs and developmental stages of Arabidopsis thaliana. Phytochem. 62, 471–481. doi: 10.1016/s0031-9422(02)00549-6
Cavill, G. W. K., Clark, D. V. (1967). Insect venoms, attractants, and repellents—VIII. Isodihydronepetalactone. J. Insect Physiol. 13, 131–135. doi: 10.1016/0022-1910(67)90009-1
Chalchat, J. C., Garry, R. P., Michet, A. (1998). Essential oils of myrtle (Myrtus communis L.) of the Mediterranean littoral. J. Essent. Oil Res. 10, 613–617. doi: 10.1080/10412905.1998.9700988
Chauhan, K. R., Khanna, H., Bathini, N. B., Le, T. C., Grieco, J. (2014). Biobased lactams as novel arthropod repellents. Nat. Prod. Commun. 9, 1934578X1400901201. doi: 10.1177/1934578X1400901201
Chauhan, K. R., Zhang, A. (2008). Methods of separating Z, E-nepetalactone and Z, E-nepetalactone from catnip oil.
Cruden, R. W., Miller-Ward, S. (1981). Pollen-ovule ratio, pollen size, and the ratio of stigmatic area to the pollen-bearing area of the pollinator: an hypothesis. Evolution. 35, 964–974. doi: 10.2307/2407867
Cushnie, T. T., Lamb, A. J. (2005). Antimicrobial activity of flavonoids. Int. J. Antimicrob. Agents. 26, 343–356. doi: 10.1016/j.ijantimicag.2005.09.002
Damle, M. S., Giri, A. P., Sainani, M. N., Gupta, V. S. (2005). Higher accumulation of proteinase inhibitors in flowers than leaves and fruits as a possible basis for differential feeding preference of Helicoverpa armigera on tomato (Lycopersicon esculentum Mill, Cv. Dhanashree). Phytochem. 66, 2659–2667. doi: 10.1016/j.phytochem.2005.09.006
Darrag, H. M., Almuhanna, H. T., Hakami, E. H. (2022). Secondary metabolites in basil, bio-insecticide, inhibition effect, and in silico molecular docking against proteolytic enzymes of the red palm weevil (Rhynchophorus ferrugineus). Plants. 11, 1087. doi: 10.3390/plants11081087
De Pooter, H. L., Nicolai, B., De Laet, J., De Buyck, L. F., Schamp, N. M., Goetghebeur, P. (1988). The essential oils of five Nepeta species. A preliminary evaluation of their use in chemotaxonomy by cluster analysis. Flavour Fragr J. 3, 155–159. doi: 10.1002/ffj.2730030404
Dienaitė, L., Pukalskienė, M., Matias, A. A., Pereira, C. V., Pukalskas, A., Venskutonis, P. R. (2018). Valorization of six Nepeta species by assessing the antioxidant potential, phytochemical composition and bioactivity of their extracts in cell cultures. J. Func. Foods. 45, 12–522. doi: 10.1016/j.jff.2018.04.004
Dixon, R. A. (2001). Natural products and plant disease resistance. Nature 411, 843–847. doi: 10.1038/35081178
Dötterl, S., Gershenzon, J. (2023). Chemistry, biosynthesis and biology of floral volatiles: roles in pollination and other functions. Nat. Prod. Rep. 40, 1901–1937. doi: 10.1039/D3NP00024A
Effmert, U., Große, J., Röse, U. S. R., Ehrig, F., Kägi, R., Piechulla, B. (2005). Volatile composition, emission pattern, and localization of floral scent emission in Mirabilis Jalapa (Nyctaginaceae). Am. J. Bot. 92, 2–12. doi: 10.3732/ajb.92.1.2
Eisenbraun, E. J., Sullins, D. W., Browne, C. E., Shoolery, J. N. (1988). (4aS, 7S, 7aR)-Nepetalactam and (4aS, 7S, 7aR)-2-[(3R, 4R, 4aR, 7S, 7aR)-octahydro-4, 7-dimethyl-1-oxocyclopenta [c] pyran-3-yl] nepetalactam. Nitrogen analogs of nepetalactone and nepetalic. psi.-anhydride. J. Org. Chem. 53, 3968–3972. doi: 10.1021/jo00252a016
El-Moaty, H. I. A. (2010). Essential oil and iridoid glycosides of Nepeta septemcrenata Erenb. J. Natural Products 3, 103–111. Available at: https://www.scinapse.io/papers/2188484169.
Erkan, N., Ayranci, G., Ayranci, E. (2008). Antioxidant activities of rosemary (Rosmarinus Officinalis L.) extract, blackseed (Nigella sativa L.) essential oil, carnosic acid, rosmarinic acid and sesamol. Food Chem. 110, 76–82. doi: 10.1016/j.foodchem.2008.01.058
Feaster, J. E., Scialdone, M. A., Todd, R. G., Gonzalez, Y. I., Foster, J. P., Hallahan, D. L. (2009). Dihydronepetalactones deter feeding activity by mosquitoes, stable flies, and deer ticks. J. Med. Entomol. 46, 832–840. doi: 10.1603/033.046.0413
Fernández-Grandon, G. M., Poppy, G. M. (2015). Response of Aphidius colemani to aphid sex pheromone varies depending on plant synergy and prior experience. Bull. Entom Res. 105, 507–514. doi: 10.1017/S0007485315000371
Formisano, C., Rigano, D., Senatore, F. (2011). Chemical constituents and biological activities of Nepeta species. Chem. Biodivers. 8, 1783–1818. doi: 10.1002/cbdv.201000191
Gasic, K., Hernandez, A., Korban, S. S. (2004). RNA extraction from different apple tissues rich in polyphenols and polysaccharides for cDNA library construction. Plant Mol. Biol. Rep. 22, 437–438. doi: 10.1007/BF02772687
Gkinis, G., Bozin, B., Mimica-Dukic, N., Tzakou, O. (2010). Antioxidant activity of Nepeta nuda L. ssp. nuda essential oil rich in nepetalactones from Greece. J. medicinal Food 13, 1176–1181. doi: 10.1089/jmf.2009.0218
Glinwood, R. T., Du, Y. J., Smiley, D. W. M., Powell, W. (1999). Comparative responses of parasitoids to synthetic and plant-extracted nepetalactone component of aphid sex pheromones. J. Chem. Ecol. 25, 1481–1488. doi: 10.1023/A:1020872412891
Godschalx, A. L., Stady, L., Watzig, B., Ballhorn, D. J. (2016). Is protection against florivory consistent with the optimal defense hypothesis? BMC Plant Biol. 16, 1–9. doi: 10.1186/s12870-016-0719-2
Goldansaz, S. M., Festa, C., Pagano, E., De Marino, S., Finamore, C., Parisi, O. A., et al. (2019). Phytochemical and biological studies of Nepeta asterotricha Rech. f. (Lamiaceae): Isolation of nepetamoside. Molecules 24, 1684. doi: 10.3390/molecules24091684
Gomes, E. N., Yuan, B., Patel, H. K., Lockhart, A., Wyenandt, C. A., Wu, Q., et al. (2024). Implications of the Propagation Method for the Phytochemistry of Nepeta cataria L. throughout a Growing Season. Molecules 29, 2001. doi: 10.3390/molecules29092001
Gormez, A., Bozari, S., Yanmis, D., Gulluce, M., Agar, G., Sahin, F. (2013). Antibacterial activity and chemical composition of essential oil obtained from Nepeta nuda against phytopathogenic bacteria. J. Essent. Oil Res. 25, 149–153. doi: 10.1080/10412905.2012.751060
Guiso, M., Martino, C. (1983). 6-Deoxylamioside, a new iridoid glucoside from Lamium amplexicaule. J. Nat. Prod. 46, 157–160. doi: 10.1021/np50026a002
Hadi, N., Sefidkon, F., Shojaeiyan, A., Šiler, B., Jafari, A. A., Aničić, N., et al. (2017). Phenolics’ composition in four endemic Nepeta species from Iran cultivated under experimental field conditions: The possibility of the exploitation of Nepeta germplasm. Ind. Crops Prod. 95, 475–484. doi: 10.1016/j.indcrop.2016.10.059
Hammer, Ø., Harper, D. A. (2001). Past: paleontological statistics software package for educaton and data anlysis. Palaeontologia electronica 4, 1. Available at: http://palaeo-electronica.org/2001_1/past/issue1_01.htm.
Handjieva, N. V., Popov, S. S., Evstatieva, L. N. (1996). Constituents of Essential Oils from Nepeta cataria L., N. grandiflora MB and N. nuda L. J. Essent. Oil Res. 8, 639–643. doi: 10.1080/10412905.1996.9701032
Harney, J. W., Barofsky, I. M., Leary, J. D. (1978). Behavioral and toxicological studies of cyclopentanoid monoterpenes from Nepeta cataria. Lloydia 41, 367–374.
Hernaíndez Lozada, N. J., Hong, B., Wood, J. C., Caputi, L., Basquin, J., Chuang, L., et al. (2022). Biocatalytic routes to stereo-divergent iridoids. Nat. Commun. 13, 1–13. doi: 10.1038/s41467-022-32414-w
Ivanov, M., Kostić, M., Stojković, D., Soković, M. (2022). Rosmarinic acid–modes of antimicrobial and antibiofilm activities of a common plant polyphenol. S AFR J. Bot. 146, 521–527. doi: 10.1016/j.sajb.2021.11.050
Jamzad, Z., Grayer, R. J., Kite, G. C., Simmonds, M. S., Ingrouille, M., Jalili, A. (2003). Leaf surface flavonoids in Iranian species of Nepeta (Lamiaceae) and some related genera. Biochem. Syst. Ecol. 31, 587–600. doi: 10.1016/S0305-1978(02)00221-1
Junker, R. R. (2016). “Multifunctional and diverse floral scents mediate biotic interactions embedded in communities,” in Deciphering Chemical Language of Plant Communication. Signaling and Communication in Plants. Eds. Blande, J., Glinwood, R. (Springer, Cham). doi: 10.1007/978-3-319-33498-1_11
Kessler, A., Kalske, A. (2018). Plant secondary metabolite diversity and species interactions. Annu. Rev. Ecol. Evol. Syst. 49, 115–138. doi: 10.1146/annurev-ecolsys-110617-062406
Khan, S., Taning, C. N. T., Bonneure, E., Mangelinckx, S., Smagghe, G., Ahmad, R., et al. (2019). Bioactivity-guided isolation of rosmarinic acid as the principle bioactive compound from the butanol extract of Isodon rugosus against the pea aphid, Acyrthosiphon pisum. PLoS One 14, e0215048. doi: 10.1371/journal.pone.0215048
Kilic, O., Hayta, S., Bagci, E. (2011). Chemical composition of essential oil of Nepeta nuda L. subsp. nuda (Lamiaceae) from Turkey. Asian J. Chem. 23, 2788.
Klocke, J. A., Darlington, M. V., Balandrin, M. F. (1987). 1, 8-Cineole (Eucalyptol), a mosquito feeding and ovipositional repellent from volatile oil of Hemizonia fitchii (Asteraceae). J. Chem. Ecol. 13, 2131–2141. doi: 10.1007/BF01012562
Kobaisy, M., Tellez, M. R., Dayan, F. E., Mamonov, L. K., Mukanova, G. S., Sitpaeva, G. T., et al. (2005). Composition and phytotoxic activity of Nepeta pannonica L. essential oil. J. Essent. Oil Res. 17, 704–707. doi: 10.1080/10412905.2005.9699037
Kökdil, G., Kurucu, S., Topçu, G. (1996). Composition of the essential oil of Nepeta nuda L. ssp. albiflora (Boiss.) Gams. Flavour Fragr J. 11, 167–169. doi: 10.1002/(SICI)1099-1026(199605)11:3%3C167::AID-FFJ567%3E3.0.CO,2-A
Kökdil, G., Kurucu, S., Yıldız, A. (1998). Essential oil composition of Nepeta nuda L. ssp. nuda. Flavour Fragr J. 13, 233–234. doi: 10.1002/(SICI)1099-1026(1998070)13:4%3C233::AID-FFJ730%3E3.0.CO,2-7
Kökdil, G., Yalcin, S. M., Topcu, G.Ü.L.A. Ç. T. I. (1999). Nepetalactones and other constituents of Nepeta nuda ssp. albiflora. Turk. J. Chem. 23, 99–104. https://journals.tubitak.gov.tr/chem/vol23/iss1/14
Labarrere, B., Prinzing, A., Dorey, T., Chesneau, E., Hennion, F. (2019). Variations of secondary metabolites among natural populations of sub-Antarctic ranunculus species suggest functional redundancy and versatility. Plants 8, 234. doi: 10.3390/plants8070234
Li, J., Hu, H., Mao, J., Yu, L., Stoopen, G., Wang, M., et al. (2019). Defense of pyrethrum flowers: repelling herbivores and recruiting carnivores by producing aphid alarm pheromone. New phytologist. 223, 1607–1620. doi: 10.1111/nph.15869
Lichman, B. R., Godden, G. T., Hamilton, J. P., Palmer, L., Kamileen, M. O., Zhao, D., et al. (2020). The evolutionary origins of the cat attractant nepetalactone in catnip. Sci. Adv. 6, eaba0721. doi: 10.1126/sciadv.aba0721
Lichman, B. R., Kamileen, M. O., Titchiner, G. R., Saalbach, G., Stevenson, C. E. M., Lawson, D. M., et al. (2019a). Uncoupled activation and cyclization in catmint reductive terpenoid biosynthesis. Nat. Chem. Biol. 15, 71–79. doi: 10.1038/s41589-018-0185-2
Lichman, B. R., O’Connor, S. E., Kries, H. (2019b). Biocatalytic strategies towards [4 + 2] cycloadditions. Chem. - A Eur. J. 25, 6864–6877. doi: 10.1002/chem.201805412
Livak, K. J., Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2– ΔΔCT method. methods 25, 402–408. doi: 10.1006/meth.2001.1262
Mamadalieva, N. Z., Sharopov, F. S., Satyal, P., Azimova, S. S., Wink, M. (2016). Chemical composition of the essential oils of some Central Asian Nepeta species (Lamiaceae) by GLC-MS. Nat. Prod. Commun. 11, 1934578X1601101229. doi: 10.1177/1934578X1601101229
Mancini, E., Apostolides Arnold, N., De Feo, V., Formisano, C., Rigano, D., Piozzi, F., et al. (2009). Phytotoxic effects of essential oils of Nepeta curviflora Boiss. and Nepeta nuda L. subsp. albiflora growing wild in Lebanon. J. Plant Interact. 4, 253–259. doi: 10.1080/17429140903225507
Marak, H. B., Biere, A., Van Damme, J. M. (2002). Two herbivore-deterrent iridoid glycosides reduce the in-vitro growth of a specialist but not of a generalist pathogenic fungus of Plantago lanceolata L. Chemoecology. 12, 185–192. Available at: https://www.jstor.org/stable/3448829.
Matsuno, M., Nagatsu, A., Ogihara, Y., Ellis, B. E., Mizukami, H. (2002). CYP98A6 from Lithospermum erythrorhizon encodes 4-coumaroyl-4′-hydroxyphenyllactic acid 3-hydroxylase involved in rosmarinic acid biosynthesis. FEBS Lett. 514, 219–224. doi: 10.1016/S0014-5793(02)02368-2
McCall, A. C., Irwin, R. E. (2006). Florivory: the intersection of pollination and herbivory. Ecol. 9, 1351–1365. doi: 10.1111/j.1461-0248.2006.00975.x
Mišić, D., Siler, B., Gašić, U., Avramov, S., Zivković, S., Nestorović Živković, J., et al. (2015). Simultaneous UHPLC/DAD/(+/-)HESI-MS/MS analysis of phenolic acids and nepetalactones in methanol extracts of Nepeta species: a possible application in chemotaxonomic studies. Phytochem. Anal. 26, 72–85. doi: 10.1002/pca.2538
Moghadam, S. E., Ebrahimi, S. N., Gafner, F., Ochola, J. B., Marubu, R. M., Lwande, W., et al. (2015). Metabolite profiling for caffeic acid oligomers in Satureja biflora. Ind. Crops Prod. 76, 892–899. doi: 10.1016/j.indcrop.2015.07.059
Moore, B. D., Andrew, R. L., Külheim, C., Foley, W. J. (2014). Explaining intraspecific diversity in plant secondary metabolites in an ecological context. New Phytol. 201, 733–750. doi: 10.1111/nph.12.526
Moretti, A. N., Zerba, E. N., Alzogaray, R. A. (2015). Lethal and sublethal effects of eucalyptol on T riatoma infestans and R hodnius prolixus, vectors of C hagas disease. Entomol. Exp. Appl. 154, 62–70. doi: 10.1111/eea.12256
Muñoz-Muñoz, J. L., Garcia-Molina, F., Ros, E., Tudela, J., García-Canovas, F., Rodriguez-Lopez, J. N. (2013). Prooxidant and antioxidant activities of rosmarinic acid. J. Food Biochem. 37, 396–408. doi: 10.1111/j.1745-4514.2011.00639.x
Murai, F., Tagawa, M., Damtoft, S., Jensen, S. R., Nielsen, B. J. (1984). (1R, 5R, 8S, 9S)-Deoxyloganic acid from Nepeta cataria. Chem. Pharm. Bull. 32, 2809–2814. doi: 10.1248/cpb.32.2809
Musso, L., Scaglia, B., Haj, G. A., Arnold, N. A., Adani, F., Scarì, G., et al. (2017). Chemical characterization and nematicidal activity of the essential oil of Nepeta nuda L. ssp. pubescens and Nepeta curviflora Boiss. from Lebanon. J. Essent. Oil-Bear. Plants. 20, 1424–1433. doi: 10.1080/0972060X.2017.1407678
Nadeem, M., Imran, M., Aslam Gondal, T., Imran, A., Shahbaz, M., Muhammad Amir, R., et al. (2019). Therapeutic potential of rosmarinic acid: A comprehensive review. Appl. Sci. 9, 3139. doi: 10.3390/app9153139
Nagy, T., Kocsis, Á., Morvai, M., Szabó, L., Podányi, B., Gergely, A., et al. (1998). 2′-, 4′-, and 6′-O-substituted 1, 5, 9-epideoxyloganic acids from Nepeta grandiflora. Phytochemistry. 47, 1067–1072. doi: 10.1016/S0031-9422(98)80074-5
Narimani, R., Moghaddam, M., Ghasemi Pirbalouti, A., Mojarab, S. (2017). Essential oil composition of seven populations belonging to two Nepeta species from Northwestern Iran. Int. J. Food Prop. 20, 2272–2279. doi: 10.1080/10942912.2017.1369104
Palmer-Young, E. C., Farrell, I. W., Adler, L. S., Milano, N. J., Egan, P. A., Junker, R. R., et al. (2019). Chemistry of Floral Rewards: Intra-and Interspecific Variability of Nectar and Pollen Secondary Metabolites across Taxa. Ecol. Monogr. 89, 1–23. Available at: https://www.jstor.org/stable/26641226.
Patel, H., Gomes, E. N., Yuan, B., Lyu, W., Wu, Q., Simon, J. E. (2022). Investigation of volatile iridoid terpenes in Nepeta cataria L. (Catnip) Genotypes. Molecules. 27, 7057. doi: 10.3390/molecules27207057
Petersen, M., Abdullah, Y., Benner, J., Eberle, D., Gehlen, K., Hücherig, S., et al. (2009). Evolution of rosmarinic acid biosynthesis. Phytochemistry. 70, 1663–1679. doi: 10.1016/j.phytochem.2009.05.010
Petrova, D., Gašić, U., Yocheva, L., Hinkov, A., Yordanova, Z., Chaneva, G., et al. (2022). Catmint (Nepeta nuda L.) phylogenetics and metabolic responses in variable growth conditions. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.866777
Petrović, L., Skorić, M., Šiler, B., Banjanac, T., Gašić, U., Matekalo, D., et al. (2024). Patterns of genetic variation of nepeta nuda L. from the central balkans: understanding drivers of chemical diversity. Plants. 13, 1483. doi: 10.3390/plants13111483
Pichersky, E., Gershenzon, J. (2002). The formation and function of plant volatiles: perfumes for pollinator attraction and defense. Curr. Opin. Plant Biol. 5, 237–243. doi: 10.1016/s1369-5266(02)00251-0
Rabee, M., Andersen, Ø.M., Fossen, T., Enerstvedt, K. H., Abu Ali, H., Rayyan, S. (2020). Acylated flavone O-glucuronides from the aerial parts of nepeta curviflora. Molecules. 25, 3782. doi: 10.3390/molecules25173782
Regnier, F. E., Waller, G. R., Eisenbraun, E. J. (1967). Studies on the composition of the essential oils of three Nepeta species. Phytochemistry. 6, 1281–1289. doi: 10.1016/S0031-9422(00)86089-6
Reichert, W., Ejercito, J., Guda, T., Dong, X., Wu, Q., Ray, A., et al. (2019). Repellency assessment of Nepeta cataria essential oils and isolated nepetalactones on Aedes aEgypti. Sci. Rep. 9, 1524. doi: 10.1038/s41598-018-36814-1
Sastry, S. D., Springstube, W. R., Waller, G. R. (1972). Identification of 5, 9-dehydronepetalactone, a new monoterpene from Nepeta cataria. Phytochem. 11, 453–455. doi: 10.1016/S0031-9422(00)90053-0
Schultz, G., Simbro, E., Belden, J., Zhu, J., Coats, J. (2004). Catnip, Nepeta cataria (Lamiales: Lamiaceae)—A closer look: Seasonal occurrence of nepetalactone isomers and comparative repellency of three terpenoids to insects. Environ. Entomol. 33, 1562–1569. doi: 10.1603/0046-225X-33.6.1562
Scriven, R., Meloan, C. E. (1984).Determining the active component in 1, 3, 3-trimethyl-2-oxabicyclo {2, 2, 2} octane (cineole) that repels the American cockroach, Periplaneta americana. Available online at: http://hdl.handle.net/1811/23001 (Accesed May 05, 2024).
Sengupta, S. K., Hutchenson, K. W., Hallahan, D. L., Gonzalez, Y. I., Manzer, L. E., Jackson, S. C., et al. (2018). Hydrogenation of naturally-derived nepetalactone as a topical insect repellent. ACS Sustain. Chem. Eng. 6, 9628–9639. doi: 10.1021/acssuschemeng.7b04521
Sharma, A., Cannoo, D. S. (2013). Phytochemical composition of essential oils isolated from different species of genus Nepeta of Labiatae family: a review. Pharmacophore 4 (6), 181–211. Available at: https://pharmacophorejournal.com/1sHUSCy.
Sharma, A., Cooper, R., Bhardwaj, G., Cannoo, D. S. (2021). The genus Nepeta: Traditional uses, phytochemicals and pharmacological properties. J. Ethnopharmacol. 268, 113679. doi: 10.1016/j.jep.2020.113679
Sherden, N. H., Lichman, B., Caputi, L., Zhao, D., Kamileen, M. O., Buell, C. R., et al. (2018). Identification of iridoid synthases from Nepeta species: Iridoid cyclization does not determine nepetalactone stereochemistry. Phytochem. 145, 48–56. doi: 10.1016/j.phytochem.2017.10.004
Sih, A., Baltus, M. S. (1987). Patch size, pollinator behavior, and pollinator limitation in catnip. Ecology. 68, 1679–1690. doi: 10.2307/1939860
Simmonds, M. S., Stevenson, P. C., Hanson, F. E. (2019). Rosmarinic acid in Canna generalis activates the medial deterrent chemosensory neurone and deters feeding in the tobacco hornworm Manduca sexta. Physiol. Entomol. 44, 140–147. doi: 10.1111/phen.12284
Singh, S., Kaur, I., Kariyat, R. (2021). The multifunctional roles of polyphenols in plant-herbivore interactions. Int. J. Mol. Sci. 22, 1442. doi: 10.3390/ijms22031442
Slavković, F., Bendahmane, A. (2023). Floral phytochemistry: impact of volatile organic compounds and nectar secondary metabolites on pollinator behavior and health. Chem. Biodivers, e202201139. doi: 10.1002/cbdv.202201139
Smallegange, R. C., van Loon, J. J. A., Blatt, S. E., Harvey, J. A., Agerbirk, N., Dicke, M. (2007). Flower vs. leaf feeding by Pieris brassicae: glucosinolate-rich flower tissues are preferred and sustain higher growth rate. J. Chem. Ecol. 33, 1831–1844. doi: 10.1007/s10886-007-9350-x
Smiljković, M., Dias, M. I., Stojkovic, D., Barros, L., Bukvički, D., Ferreira, I. C. F. R. (2018). Characterization of phenolic compounds in tincture of edible Nepeta nuda: development of antimicrobial mouthwash. Food Funct. 9, 5417–5425. doi: 10.1039/c8fo01466c
Sparks, J. T., Botsko, G., Swale, D. R., Boland, L. M., Patel, S. S., Dickens, J. C. (2018). Membrane proteins mediating reception and transduction in chemosensory neurons in mosquitoes. Front. Physiol. 9. doi: 10.3389/fphys.2018.01309
Sukontason, K. L., Boonchu, N., Sukontason, K., Choochote, W. (2004). Effects of eucalyptol on house fly (Diptera: Muscidae) and blow fly (Diptera: Calliphoridae). Rev. Inst Med. Trop. Sao Paulo. 46, 97–101. doi: 10.1590/S0036-46652004000200008
Süntar, I., Nabavi, S. M., Barreca, D., Fischer, N., Efferth, T. (2018). Pharmacological and chemical features of Nepeta L. genus: Its importance as a therapeutic agent. Phytother. Res. 32, 185–198. doi: 10.1002/ptr.5946
Suthar, M. K., Purohit, P. M., Saran, P. L. (2021). Molecular cloning, expression and in-silico characterization of rosmarinic acid synthase from Ocimum tenuiflorum L. J. Plant Biochem. Biotechnol. 30, 317–325. doi: 10.1007/s13562-020-00591-8
Takeda, Y., Morimoto, Y., Matsumoto, T., Honda, G., Tabata, M., Fujita, T., et al. (1995). Nepetanudoside, an iridoid glucoside with an unusual stereostructure from Nepeta nuda ssp. Albiflora. J. Nat. Prod. 58, 1217–1221. doi: 10.1021/np50122a009
Takeda, Y., Yagi, T., Matsumoto, T., Honda, G., Tabata, M., Fujita, T., et al. (1996). Nepetanudosides and iridoid glucosides having novel stereochemistry from Nepeta nuda ssp. albiflora. Phytochem. 42, 1085–1088. doi: 10.1016/0031-9422(96)00074-X
Talebi, S. M., Nohooji, M. G., Yarmohammadi, M. (2017). Infraspecific variations in essential oil compositions of Nepeta fissa from Iran. Nusantara. Bioscience. 9, 318–321. doi: 10.13057/nusbiosci/n090313
Talebi, S. M., Nohooji, M. G., Yarmohammadi, M., Khani, M., Matsyura, A. (2019). Effect of altitude on essential oil composition and on glandular trichome density in three Nepeta species (N. sessilifolia, N. heliotropifolia and N. fissa). Mediterr. Botany. 40, 81–93. doi: 10.5209/mbot.59730
Thiers, B. (2023). Index Herbariorum: A Global Directory of Public Herbaria and Associated Staff (New York Botanical Garden’s Virtual Herbarium).
Treutter, D. (2005). Significance of flavonoids in plant resistance and enhancement of their biosynthesis. Plant Biol. 7, 581–591. doi: 10.1055/s-2005-873009
Trócsányi, E., György, Z., Zámboriné-Németh, É. (2020). New insights into rosmarinic acid biosynthesis based on molecular studies. Curr. Plant Biol. 23, 100162. doi: 10.1016/j.cpb.2020.100162
Tyagi, K., Shukla, P., Rohela, G. K., Shabnam, A. A., Gautam, R. (2020). Plant phenolics: their biosynthesis, regulation, evolutionary significance, and role in senescence. Plant Phenolics Sustain. Agriculture. 1, 431–449. doi: 10.1007/978-981-19-6426-8_11
Uenoyama, R., Miyazaki, T., Adachi, M., Nishikawa, T., Hurst, J. L., Miyazaki, M. (2022). Domestic cats damage silver vine and catnip plants containing iridoids and enhance chemical pest defence. iScience. 25, 104455. doi: 10.1016/j.isci.2022.104455
Uenoyama, R., Miyazaki, T., Hurst, J. L., Beynon, R. J., Adachi, M., Murooka, T., et al. (2021). The characteristic response of domestic cats to plant iridoids allows them to gain chemical defense against mosquitoes. Sci. Adv. 7, eabd9135. doi: 10.1126/sciadv.abd9135
Ulziibayar, D., Begna, T., Bisrat, D., Jung, C. (2022). Longevity-enhancing Effects of Rosmarinic Acid Feeding on Honey bees (Apis mellifera L.) after Exposure to Some Pesticides Used in Strawberry Greenhouse. J. Apic. 2, 135–141. doi: 10.17519/apiculture.2022.06.37.2.135
Velasco-Negueruela, A., Perez-Alonso, M. J., Rodriquez, A. B. (1989). Essential oil analysis of Nepeta teydea Webb. & Berth. Flavour Fragr J. 4, 197–199. doi: 10.1002/ffj.2730040409
Vivekanandhan, P., Swathy, K., Bedini, S., Shivakumar, M. S. (2023). Bioprospecting of Metarhizium anisopliae derived crude extract: a ecofriendly insecticide against insect pest. Int. J. Trop. Insect Sci. 43, 1–12. doi: 10.1007/s42690-022-00935-y
Vivekanandhan, P., Swathy, K., Shivakumar, M. S. (2022). Identification of insecticidal molecule aucubin from Metarhizium anisopliae ethyl acetate crude extract against disease mosquito vector. Int. J. Trop. Insect Sci. 42, 3303–3318. doi: 10.1007/s42690-022-00828-0
Wang, G.-Q., Chen, J.-F., YI, B., Tan, H.-X., Zhang, L., Chen, W.-S. N. (2017). HPPR encodes the hydroxyphenylpyruvate reductase required for the biosynthesis of hydrophilic phenolic acids in Salvia miltiorrhiza. Chin. J. Nat. Med. 15, 917–927. doi: 10.1016/S1875-5364(18)30008-6
Wetzel, W. C., Whitehead, S. R. (2020). The many dimensions of phytochemical diversity: linking theory to practice. Ecol. Lett. 23, 16–32. doi: 10.1111/ele.13422
Wink, M. (2018). Plant secondary metabolites modulate insect behavior-steps toward addiction? Front. Physiol. 9. doi: 10.3389/fphys.2018.00364
Yarmoohammadi, M., Talebi, S. M., Nohooji, M. G. (2017). Infraspecific variations in essential oil and glandular trichomes in Nepeta heliotropifolia. Biodiversitas J. Biol. Diversity. 18, 964–970. doi: 10.13057/biodiv/d180314
Zaharieva, A., Rusanov, K., Rusanova, M., Paunov, M., Yordanova, Z., Mantovska, D., et al. (2023). Uncovering the interrelation between metabolite profiles and bioactivity of in vitro- and wild-grown catmint (Nepeta nuda L.). Metabolites. 13, 1099. doi: 10.3390/metabo13101099
Zangerl, A. R. (2003). Evolution of induced plant responses to herbivores. Basic Appl. Ecol. 4, 91–103. doi: 10.1078/1439-1791-00135
Keywords: metabolomics, UHPLC-ESI-QToF-MS, UHPLC/DAD/(±)HESI-MS2, GC/MS, phenolics, terpenes, iridoids, populations
Citation: Petrović L, Filipović B, Skorić M, Šiler B, Banjanac T, Matekalo D, Nestorović Živković J, Dmitrović S, Aničić N, Milutinović M, Božunović J, Gašić U and Mišić D (2024) Molecular background of the diverse metabolic profiles in leaves and inflorescences of naked catmint (Nepeta nuda L.). Front. Plant Sci. 15:1452804. doi: 10.3389/fpls.2024.1452804
Received: 21 June 2024; Accepted: 11 November 2024;
Published: 28 November 2024.
Edited by:
Moonhyuk Kwon, Gyeongsang National University, Republic of KoreaReviewed by:
Juan De Dios Franco-Navarro, Spanish National Research Council (CSIC), SpainJongmin Ahn, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Republic of Korea
Copyright © 2024 Petrović, Filipović, Skorić, Šiler, Banjanac, Matekalo, Nestorović Živković, Dmitrović, Aničić, Milutinović, Božunović, Gašić and Mišić. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Biljana Filipović, YmlsamFuYS5uaWtvbGljQGliaXNzLmJnLmFjLnJz; Marijana Skorić, bWRldmljQGliaXNzLmJnLmFjLnJz; Uroš Gašić, dXJvcy5nYXNpY0BpYmlzcy5iZy5hYy5ycw==; Danijela Mišić, ZG1pc2ljQGliaXNzLmJnLmFjLnJz