agronomy Pilot Cultivation of the Local Endemic Cretan Marjoram Origanum microphyllum (Benth.) Vogel (Lamiaceae): Effect of Fertilizers on Growth and Herbal Quality Features

: Wild phytogenetic resources are threatened by overexploitation. This pressure on species and natural ecosystems can be alleviated in part by recruiting and domesticating wild-growing species under pilot cultivation with tailor-made fertilization schemes. This study focused on the pilot cultivation of Origanum microphyllum —a critically endangered local endemic plant of Crete, Greece—investigating the effect of conventional and integrated nutrient management (INM) fertilizers by foliar or root application and biostimulant. Above-ground biomass together with leaf chlorophyll ﬂuorescence and color (SPAD meter, DA meter, Chroma Meter) were determined. Leaf chlorophyll, antioxidant compound (carotenoids, phenols, ﬂavonoids), and nutrient contents were also assessed. The results showed that fertilization did not signiﬁcantly affect plant growth and leaf nutrient content. Root fertilization was associated with greener leaves compared to foliar. The same trend was generally evident for antioxidant compound content. The small size of leaves may have impeded the efﬁciency of the foliar application. In conclusion, root application of conventional or INM fertilizers seems more suitable to promote visual quality and herbal antioxidant proﬁle of O. microphyllum , than the foliar one.


Introduction
Conventional and complementary herbal medicines are traded globally both in developing and developed countries [1,2]. In the pharmaceutical industry, a considerable portion of target drugs still relies directly or indirectly on herbal material [2]. Herbal raw materials are also essential in other sectors (e.g., perfumery, cosmetics, and food industry) with a growing respective demand [3][4][5]. Although it is largely unnoticed by the public, a sizeable merit of the traded herbal material is yet derived from wild-growing plant populations of aromatic-medicinal plants (MAPs), threatening them at local scales and often across their ranges [6]. Many MAPs are nowadays threatened with extinction due to unsustainable practices, over-exploitation, and limited or compromised awareness [7]. In this perspective, considerable ecological risks arise not only for the range-restricted species, but also for the whole habitat in which these MAPs thrive in the wild [6]. Given the growing demand for herbal material in various economic sectors and the accelerating depletion of wild phytogenetic resources due to over-collection for direct use, the domestication of targeted species of interest and their pilot cultivation including the development of fertilization protocols appears to be a sustainable alternative.
Apart from developing new varieties of existing crops with superior characteristics in the face of climate change impacts, there is also a modern tendency in developing new crops to support livelihoods with a clear focus on the sustainable exploitation of neglected and underutilized species [8][9][10][11][12]. In agricultural crops, the nutrient requirements are conventionally satisfied by applying inorganic fertilizers. While stimulating plant growth and productivity when applied, their excessive use in several cases comes at the cost of a sizeable environmental impact [13]. For instance, the N and P surplus not recovered by the soil-plant system may drive eutrophication events in both freshwater systems and coastal zones [13,14]. In addition, the production of conventional fertilizers contributes to greenhouse gas emissions [15]. By contrast, organic fertilizers have been associated with low environmental footprint and decreased risks due to the slow release of plant nutrients to soil [15]. However, because of the latter the beneficial effect of organic fertilizers on crop productivity is sometimes deemed ambiguous. Therefore, replacing conventional with organic fertilizers is aimed to create a more eco-friendly and sustainable farming system. These environmental benefits can be further enhanced by the use of biostimulants [16]. Although biostimulants do not conventionally include nutrients and therefore cannot substitute fertilizers, they can decrease their utilization by stimulating increased root nutrient absorption [16]. Towards that direction, the combined use of conventional, organic fertilizers and biostimulants (the so-called integrated nutrient management; INM) represents a strategy that could sustain high yields and limit undesirable impact [17,18].
It is known that diverse fertilization strategies may differentially affect not only plant growth and yield, but also product quality [15]. Herbal material quality encompasses several features, the relative importance of which differs based on the intended use [4,5]. A critical indicator of the herbal material quality is the content of antioxidant compounds determining high-added value [14,19]. These secondary metabolites act protectively towards the influence of several external stressors and limit many diseases in humans [20,21]. Among antioxidants, carotenoids, flavonoids, and phenolic compounds are of particular interest due to their strong potential to reduce stress events both at cellular and organismal levels [14,19]. There is already a lot of evidence that the herbal material produced with the use of organic and INM fertilizers has elevated the content of antioxidants as compared to conventional fertilizers [20,[22][23][24]. In this respect, organic and INM fertilizers are an effective means of upgrading herbal material quality as related to health promoting compound content. There is also evidence that schemes with organic and INM fertilizers may stimulate leaf chlorophyll content, thus resulting to dark green color in leaf appearance which is often associated with high ornamental value [25]. The level and uniformity of leaf greenness is an additional quality indicator, which is conventionally employed from harvest to sale, involving growers, processors, distributors, and consumers [4].
This study deals with the development of a pilot and tailor-made cultivation and fertilization protocol for the Cretan marjoram Origanum microphyllum (Benth.) Vogel (Lamiaceae), aiming to optimize both yield and key herbal quality traits (color, antioxidant compound content). O. microphyllum is a local endemic plant of western and eastern Crete (Greece) [26]; this Cretan marjoram is well appreciated by local people, especially in the western part of the island where it is often directly harvested from the wild and traded at local markets in dried form particularly for its calming and anti-spasmodic properties and pleasant spicy scent [12]. Our study herein represents the first attempt to cultivate this Critically Endangered Cretan endemic plant [27] in its original native range (Crete), serving both conservation purposes and sustainable exploitation of this neglected and underutilized species. The latter is envisaged due to the multiple interest of the Cretan marjoram in three economic sectors, i.e., the ornamental [10], the agro-alimentary [11], and the medicinal-cosmetic sector [12].

Origin of Plant Material
In 2019, botanical expeditions were organized to explore different areas for wildgrowing populations of O. microphyllum with natural vigorous growth in rocky wild habitats of western and eastern Crete ( The herbarium samples of O. microphyllum were taxonomically identified, and a unique IPEN (International Plant Exchange Network) accession number was assigned to them by the Institute of Plant Breeding and Genetic Resources (IPBGR) of the Hellenic Agricultural Organization Demeter (ELGO Demeter). New plants were initially raised ex-situ through seed germination trials (pretreatment by seed immersing in 50 ppm GA3 solution for 24 h) and pilot rooting of cuttings of the IPEN accession number GR-1-BBGK-19,31 assigned to the original plant collection from Mt. Dikti, eastern Crete at altitude of 1000 m. These pilot propagation trials resulted in adequate number of ex situ raised plants for further field experimentation which were transplanted in 2 L plastic pots from the company of AFI GLAVAKI KE SIA OE, Tree and Plant Nurseries, Aridea PELLAS, GR-58400, Greece.

Establishment of Field Experiment
The field experiment was located in a research field of 20 × 25 m (35 • 19 N, 25 • 6 E, altitude 60 m above sea level) of the Hellenic Mediterranean University, Heraklion, Crete (Greece). On 21 January 2021, potted plants of O. microphyllum were established in field at 80 cm distance between rows and at 40 cm between plant individuals. The length of each row was 20 m with east-west direction. The experiment also included as "guard plants" other local Cretan endemics, i.e., Carlina diae (Rech. f.) Meusel and Kästner, Sideritis syriaca L. subsp. syriaca, Origanum dictamnus L., and Verbascum arcturus L.) established in equal numbers as the species under study, thus reducing heterogeneity in air velocity between inner and outer rows.
The experimental design involved completely randomized blocks (three blocks per treatment; see below) and 10 O. microphyllum plants per block, which were randomly planted in different rows. Six fertilization treatments were introduced starting from 11 March 2021 and they were weekly performed till the final harvest. Irrigation was supplied through an automatic system using 2 L h −1 adjustable drippers, supplying the plants with water three times per week. No visible pests and diseases were necessary to control, but the elimination of weeds was frequently required and it was manually achieved. The final harvest of O. microphyllum was completed at the end of June 2021. situ through seed germination trials (pretreatment by seed immersing in 50 ppm GA3 solution for 24 h) and pilot rooting of cuttings of the IPEN accession number GR-1-BBGK-19,31 assigned to the original plant collection from Mt. Dikti, eastern Crete at altitude of 1000 m. These pilot propagation trials resulted in adequate number of ex situ raised plants for further field experimentation which were transplanted in 2 L plastic pots from the company of AFI GLAVAKI KE SIA OE, Tree and Plant Nurseries, Aridea PELLAS, GR-58400, Greece.

Establishment of Field Experiment
The field experiment was located in a research field of 20 × 25 m (35°19′ N, 25°6′ E, altitude 60 m above sea level) of the Hellenic Mediterranean University, Heraklion, Crete (Greece). On 21 January 2021, potted plants of O. microphyllum were established in field at 80 cm distance between rows and at 40 cm between plant individuals. The length of each row was 20 m with east-west direction. The experiment also included as "guard plants"

Soil Analysis
A composite surface soil sample, constituted of 10 sub-samples, was obtained from the research field; it was fined using a 2 mm sieve and was analyzed in triplicate for the properties described hereafter. Particle size distribution was determined by the hydrometer method [28]. Organic C was determined by the wet oxidation method [29]. The Kjeldahl method was used to estimate the total N content [30] and CaCO 3 was measured with a calcimeter. The pH was determined in a 1:2 (w/v) water suspension, the electrical conductivity was estimated in the saturation extract (EC se ) and the sodium absorption ratio (SAR) was evaluated by the concentrations of water-soluble Na, Ca, and Mg [31], and the cation exchange capacity (CEC) was calculated by the [Co(NH 3 ) 6 ]Cl 3 method (ISO 23470).
Phosphorus was extracted using 0.5 M NaHCO 3 , pH 8.5 and was measured by the molybdenum blue-ascorbic acid method [32]. Both NO 3 and NH 4 ions available in soil were extracted with 1 M KCl and they were measured using UV-Vis spectrometry and the sodium salicylate-sodium nitroprusside method, respectively [33]. Exchangeable K, Ca, and Mg were extracted with 1 M CH 3 COONH 4 , pH 7 [34]; K was estimated with flame photometry, while Ca and Mg using atomic absorption spectrometry. The micronutrients Cu, Zn, Fe, and Mn available in soil were extracted with DTPA [35] and were measured using atomic absorption spectrometry as well. Soil available B was extracted with hot water and was estimated with the azomethine-H method using UV-Vis spectrometry [36].

Fertilization Treatments
During the experimentation with O. microphyllum fertilization, treatments were monitored weekly using INM and conventional fertilizers in liquid or soluble granule form or water (control) administered by soil and foliar applications. The foliar treatments were performed using a 5 L plastic handheld sprayer (working in low pressure) until apparent wetness, and the soil applications were manually performed (100 mL nutrient solution per plant). The INM involved four commercial semi-organic fertilizers from THEOFRASTOS company (Industrial area of Korinthos, GR-20100 Korinthos, Greece), which were made from high quality edible raw materials and limited amounts of chemical solvents. The conventional fertilizers were supplied in soluble powder or granule form, except the liquid fertilizer for micronutrients (Plex Mix). The composition and commercial details of fertilizers are demonstrated in Table S1

Plant Measurements
Periodical plant and leaf level measurements were conducted until the end of flowering of O. microphyllum at the experimental setting (25 May 2021). The selected leaf samples for measurements were fully expanded with no apparent symptoms of either insect damage or pathogen infection. In all cases, 15 min was the maximum time between sampling and the initiation of the evaluation; otherwise, leaf samples were placed in vials, flash-frozen in liquid nitrogen, and were transferred for storage in a freezer (−80 • C). Replicate leaves were sampled from separate plant individuals.

Non-Invasive Evaluation of Leaf Coloration at Three Growth Stages
Since leaf coloration is known to be affected by fertilization with implications for both visually perceived quality and photosynthetic capacity [37], three different methods were employed in the current study to evaluate leaf coloration: (i) leaf SPAD value estimating chlorophyll content with the use of a SPAD-502 meter (Konica Minolta Corp., Solna, Sweden); (ii) the index of absorbance difference (I AD ) was estimated in leaves by using the DA meter (tr DA Meter, T.R. Turoni, Italy); (iii) leaf color was calculated by using a Chroma Meter (Model CR-400, Minolta Corp., Japan). CIE L*a*b* parameters were recorded using D65 illuminants and a 10 • Standard Observer as a reference system. L* was used to specify lightness ranging from 0 (black) to 100 (white); a* was used to define the intensity in a range of green to red (negative values indicate green and positive values red); b* was used to define the intensity in a range from blue to yellow (negative values specify blue and positive values yellow [5]). All these measurements were performed in situ during three different growth stages of O. microphyllum (vegetative, early flowering, and full flowering). Three points were recorded per replicate leaf, which were further averaged, and three replicate leaves were assessed per treatment.

Non-Invasive Evaluation of Photosynthetic Performance at Three Growth Stages
The leaf photosynthetic performance was estimated by calculating the ratio of variable to maximum chlorophyll fluorescence (F v /F m ) [38,39], using a chlorophyll fluorometer (OS-30P, Op-tiSciences, Hudson, NH, USA). Leaves were first dark-adapted (≥20 min) by employing leaf clips and then the F v /F m was estimated by applying a saturated photosynthetic photon flux density of 3000 µmol m −2 s −1 . These measurements were performed in situ at three different developmental stages of O. microphyllum (vegetative, early flowering, and full flowering). Similarly, three points were recorded per replicate leaf and further averaged, and three replicate leaves were assessed per treatment.

Plant Growth
The growth of O. microphyllum plant individuals was evaluated using the aboveground part of the plants (fresh and dry weight) (± 0.01 g; MXX-412; Denver Instruments, Bohemia, NY, USA). The samples were placed in a forced-air drying oven for 72 h at 80 • C following the measurements of the fresh weight of the plants, with six replicate plants per treatment.

Leaf Photosynthetic Pigment Content
A total of 0.5 g for each sample was fine chopped and homogenized with the addition of 10 mL of 80% acetone and was then filtered. The filtered extract was diluted by adding 2 mL of 80% acetone per mL of extract. The whole process of extraction took place in a dark room since chlorophyll is light sensitive [40,41]. The final extract was subjected to reading on a spectrophotometer (Mapada UV-1800; Shanghai. Mapada Instruments Co., Ltd., Shanghai, China) and the total chlorophyll and carotenoid contents were estimated [42]. Three leaves were assessed per treatment, four samples were pooled for each replicate (collected from different plant individuals), and the assay was performed twice.

Leaf Total Phenolic and Total Flavonoid Contents
Leaf samples at 0.1 g were extracted with 1 mL of 80% aqueous methanol in an ultrasonic bath for 10 min and then were centrifuged at 15,000× g for 10 min. Following the protocol of Chen et al. (2021) [19], the content of total phenols and total flavonoids Agronomy 2022, 12, 94 7 of 17 was calculated using the Folin-Ciocalteu assay and aluminum chloride colorimetric assay, respectively. The absorbance against prepared reagent blank was referenced using a microplate reader (Infinite 200 PRO, TECAN, Männedorf, Switzerland). The total phenolic content was measured using gallic acid as the standard reference and the gallic acid equivalent (GAE) was expressed as mg per g fresh mass. The total flavonoid content was measured using rutin as the standard reference and the rutin equivalent (RUE) was expressed as mg per g fresh mass. Three leaves were assessed per treatment, four samples were pooled for each replicate (collected from different plant individuals) and the assay was performed twice.

Leaf Soluble Sugar Content
Leaf samples at 0.1 g were fine chopped and homogenized in 1 mL deionized water and incubated in a water bath at 100 • C for 30 min, following centrifugation at 15,000× g for 15 min at room temperature (25 • C). Then, 0.1 mL of the previous solution was mixed with anthrone ethyl acetate and H 2 SO 4 and the soluble sugar content was assayed in the supernatant by measuring the absorbance at 630 nm using a spectrophotometer (Mapada UV-1800; Shanghai. Mapada Instruments Co., Ltd., Shanghai, China), according to Dubois et al. (1956) [43]. Soluble sugar content was expressed per fresh weight basis (mg g −1 ). Three leaves were assessed per treatment, four samples were pooled for each replicate (collected from different plant individuals), and the assay was performed twice.

Leaf Nutrient Analyses
Leaf nutrient analysis was conducted in order to assess the role of fertilization on plant nutrition. The samples of O. microphyllum were washed with distilled water and then dried at 70 • C, weighed, ground, and analyzed for total N by the Kjeldahl method [30]. In addition, sub-samples were ashed at 500 • C for 4 h minimum [44] and then the ash was dissolved in 2 M HCl, following filtration. The filtrate was used to determine the P, K, Ca, Mg, Cu, Zn, Fe, Mn, and B content in the samples by employing the analytical methods described previously for soil analysis and the nutrient content was expressed per dry weight basis. Three replicates were assessed per treatment, four samples were pooled for each replicate (collected from different plant individuals), and the assay was performed twice.

Statistical Analysis
The data of each property were tested for homogeneity of variances (Duncan's test) and analysis of variance (ANOVA) was conducted. Then, the protected least significant differences test (LSD-test) was used for means comparison, at p ≤ 0.05 For all statistical analyses, the SAS software was used (SAS Institute, Cary, NC, USA).

Soil Properties of the Pilot Field
The soil of the pilot field had sandy-loam texture and alkaline reaction. In addition, it was calcareous with high content of CaCO 3 and low organic C content, EC se , and SAR. Moreover, the soil had high content of available NO 3 -N [45] and sufficient concentrations of available P and K [46,47]. As for the soil available B and metallic micronutrients, they ranged at levels higher or similar to the sufficiency levels reported [48] (Supplementary  Table S2).

Leaf Color
The effect of fertilization on leaf color (a critical herb quality trait) was evaluated at three growth stages (vegetative, early flowering, full flowering) by employing three techniques (SPAD meter, DA meter, Chroma Meter). At vegetative and early flowering stages, leaf SPAD value differences among treatments were not significant (Figure 2A,B). At full flowering stage, plants treated with INM-sa had higher leaf SPAD value as compared to the ones receiving foliar fertilization (i.e., INM-fa and ChF-fa; Figure 2C). However, despite the differences observed among certain fertilization treatments, in all cases the SPAD values were similar to that of control. particular growth stage, in all cases of fertilization treatments the leaf IAD was similar to that of control (Supplementary Figure S1C).
At all three growth stages, the leaf L* value (expressing lightness) of plants receiving ChF-sa or MPE-sa was lower in comparison to the rest of the treatments (Supplementary Figure S2).
At all three growth stages, the leaf a* value (expressing intensity in the green to red range) was higher in plants receiving foliar fertilization (i.e., INM-fa and ChF-fa) compared to control plants or those receiving ChF-sa and MPE-sa (Supplementary Figure S3).
At all three growth stages, the leaf b* value (expressing intensity in the blue to yellow range) was lower in plants receiving foliar fertilization (i.e., INM-fa and ChF-fa) and ΙΝΜsa compared to the rest of the treatments (Supplementary Figure S4).

Leaf Photosynthetic Performance
At three growth stages, the effect of fertilization on the overall photosynthetic efficiency was evaluated by determining Fv/Fm values (Figure 3). At vegetative stage, all fertilization treatments had Fv/Fm values similar to control ( Figure 3A). At early flowering stage, the leaves of plants receiving MPE-sa had the highest Fv/Fm value compared to the  Figure S1A). At early flowering stage, differences among all treatments were not significant (Supplementary Figure S1B), whereas at full flowering stage, plants receiving foliar fertilization (i.e., INM-fa and ChF-fa) had higher leaf I AD compared to plants receiving fertilization through the root (i.e., INM-sa, ChF-sa, and MPE-sa). Once more, at the particular growth stage, in all cases of fertilization treatments the leaf I AD was similar to that of control (Supplementary Figure S1C).

At vegetative stage, leaf I AD (index of absorbance difference) of plants receiving ChF-sa was higher compared to control plants, and those receiving INM-sa and MPE-sa (Supplementary
At all three growth stages, the leaf L* value (expressing lightness) of plants receiving ChF-sa or MPE-sa was lower in comparison to the rest of the treatments (Supplementary Figure S2).
At all three growth stages, the leaf a* value (expressing intensity in the green to red range) was higher in plants receiving foliar fertilization (i.e., INM-fa and ChF-fa) compared to control plants or those receiving ChF-sa and MPE-sa (Supplementary Figure S3).
At all three growth stages, the leaf b* value (expressing intensity in the blue to yellow range) was lower in plants receiving foliar fertilization (i.e., INM-fa and ChF-fa) and INM-sa compared to the rest of the treatments (Supplementary Figure S4).

Leaf Photosynthetic Performance
At three growth stages, the effect of fertilization on the overall photosynthetic efficiency was evaluated by determining

Leaf Photosynthetic Performance
At three growth stages, the effect of fertilization on the overall photosynthetic ciency was evaluated by determining Fv/Fm values (Figure 3). At vegetative stage, al tilization treatments had Fv/Fm values similar to control ( Figure 3A). At early flowe stage, the leaves of plants receiving MPE-sa had the highest Fv/Fm value compared t rest of the treatments, including the control ( Figure 3B). At full flowering stage, p receiving ChF-sa and MPE-sa had the highest Fv/Fm values, while the ones receiving I fa had the lowest ( Figure 3C).

Leaf Chlorophyll Content and Antioxidant Compound Content
Plants receiving ChF-sa showed the highest leaf chlorophyll content, followed b ones receiving ΙΜΝ-sa ( Figure 4A). The lowest leaf chlorophyll content was not plants receiving INM-fa ( Figure 4A).
The fertilization affected the content of leaves in important antioxidants (caroten phenols, and flavonoids). Plants receiving ChF-sa and INM-sa had the highest leaf ca noid content, followed by the ones receiving MPE-sa ( Figure 4B). Plants receiving IN had the lowest leaf carotenoid content ( Figure 4B).
Independently of the applied scheme, plants receiving fertilization showed h leaf phenolic content compared to control plants ( Figure 4C), while those receiving sa also presented higher leaf phenolic content in comparison to ChF-fa ( Figure 4C).

Leaf Chlorophyll Content and Antioxidant Compound Content
Plants receiving ChF-sa showed the highest leaf chlorophyll content, followed by the ones receiving IMN-sa ( Figure 4A). The lowest leaf chlorophyll content was noted in plants receiving INM-fa ( Figure 4A).

Leaf Soluble Sugar Content
Independently of the applied scheme, plants receiving fertilization presented higher leaf soluble sugar content compared to control and those receiving ChF-sa had lower leaf soluble sugar content compared to the ones receiving ChF-fa, INM-fa, and INM-sa ( Figure  5).  The fertilization affected the content of leaves in important antioxidants (carotenoids, phenols, and flavonoids). Plants receiving ChF-sa and INM-sa had the highest leaf carotenoid content, followed by the ones receiving MPE-sa ( Figure 4B). Plants receiving INM-fa had the lowest leaf carotenoid content ( Figure 4B).
Independently of the applied scheme, plants receiving fertilization showed higher leaf phenolic content compared to control plants ( Figure 4C), while those receiving ChF-sa also presented higher leaf phenolic content in comparison to ChF-fa ( Figure 4C).
Regardless of the fertilization scheme, plants receiving fertilization had higher leaf flavonoid content compared to control ( Figure 4D), and those receiving ChF-fa showed lower leaf flavonoid content compared to the other fertilization schemes ( Figure 4D).

Leaf Soluble Sugar Content
Independently of the applied scheme, plants receiving fertilization presented higher leaf soluble sugar content compared to control and those receiving ChF-sa had lower leaf soluble sugar content compared to the ones receiving ChF-fa, INM-fa, and INM-sa ( Figure 5).

Leaf Soluble Sugar Content
Independently of the applied scheme, plants receiving fertilization presented higher leaf soluble sugar content compared to control and those receiving ChF-sa had lower leaf soluble sugar content compared to the ones receiving ChF-fa, INM-fa, and INM-sa ( Figure  5).

Plant Growth and Leaf Nutrients Content
Fertilization was not found to affect either the above-ground fresh and dry plant mass or the water content (data not shown).
The effect of fertilization scheme on nutrient absorption by plants was performed by analyzing leaves for plant essential macronutrients and micronutrients (Tables 1 and 2, respectively). As far as the macronutrients are concerned, Mg concentration was not affected by fertilization. For the rest macronutrients, significant differences among all treatments were observed; however, except for P, the N, K, and Ca concentrations in leaves of plants receiving fertilization of any kind were similar or lower than those of control plants. In the case of P, enhanced leaf content was apparent in plants receiving ChF-fa compared to the remaining treatments, including control (Table 1).    As far as the micronutrients are concerned, significant differences among all treatments were obtained for all elements. However, Cu, Zn, and Fe concentrations in leaves of plants receiving fertilization of any kind were similar or lower than those of control plants.
Regarding the remaining elements, plants receiving foliar fertilization (i.e., INM-fa or ChF-fa) presented significantly higher leaf Mn content, while MPE-sa led to significantly higher leaf B content.

Discussion
This study represents the first step to bring in cultivation the Cretan marjoram O. microphyllum (Critically Endangered Cretan local endemic), evaluating the fertilization which is optimal for both plant growth and herbal quality. In this direction, the possibility of using organic fertilizer in place of conventional ones was pursued owing to both the lower environmental impact, and the possibility of future certification as 'organic' or 'environmentally friendly' product, either of which is valuable and well-appreciated in the value chains of herbal material markets [49][50][51]. Cretan marjoram is currently valued locally by people in western Crete for its strong sweet scent, calming and anti-spasmodic properties and therefore it is occasionally traded at local markets in dried form derived from wild-harvested mountainous populations [12]. O. microphyllum has a strong potential in the medicinal-cosmetic sector due to the fact that relevant phytochemical studies showed that its essential oil contains high levels of bioactive molecules such as terpinen-4-ol, sabinene, and linalool, exhibiting antioxidant and antiproliferative activities against human cancer cell lines [52]. Furthermore, extracts from O. microphyllum aerial parts are effective against visceral leishmaniasis, a common disease in the world today [53]. Previous studies suggest the domestication of this neglected and underutilized plant as a promising case with agro-alimentary potential due to its domestic use in culinary preparations and for tea making [11]. In fact, these uses of O. microphyllum are in the same line with those of the globally well-appreciated sweet marjoram (Origanum majorana L.); the latter is considered nowadays as a traditional herbal medicine with approved indications by the European Medicines Agency [54]. This actually reflects health claims for O. majorana and safety of use based on traditional ethnobotanical knowledge for decades or centuries, and is associated with extant marketing and distribution channels in ornamental plants which can be profited for the trade of O. microphyllum as well at local and global scales [10].
Our study offered insight regarding the fertilization of O. microphyllum under pilot cultivation. In general, fertilization of O. microphyllum exerted limited effects on plant growth, and nutrient content, because soil was adequately rich in nutrients. However, it is worth noting that foliar application of both kinds of fertilizers enhanced Mn in leaves and biostimulant application to soil increased B. In general, the nutrients' concentrations in leaves of O. microphyllum ranged at levels almost similar to values reported in literature for sweet marjoram (O. majorana) of northern Greece [55] and Turkey [56], except for Mg, Zn, and B, with the latter being higher herein compared to both studies. As far as herbal material quality is concerned, it was indeed affected by fertilization. The quality criteria of the traded herbal material largely depend on the intended use, and they usually incorporate visual quality aspects [4,5]. For example, the intensity and uniformity of greenness of plant material is typically employed as a quality rating trait throughout the production and distribution chain [4]. At harvested full flowering stage, the leaf SPAD value of plants receiving foliar fertilization (i.e., INM-fa and ChF-fa) tended to be lower than the ones receiving root fertilization (i.e., INM-sa, ChF-sa, and MPE-sa); however, the SPAD values of all fertilization treatments were similar to control ( Figure 2C). In line with these results, plants receiving both kinds of fertilization by foliar application were less green (i.e., higher I AD value) compared to plants receiving root fertilization (Supplementary Figure S1C); however, again the I AD values of all fertilization treatments were similar to control. Accordingly, leaves of plants receiving ChF-sa and MPE-sa were darker (i.e., lower L* value) than the ones of plants receiving foliar fertilization and those of control Figure S2C). In agreement with this trend, leaves of plants receiving foliar fertilization, but also INM-sa were yellower (i.e., lower b* value) compared to the rest of the treatments (Supplementary Figure S4C). When taken together, these results indicate that both kinds of fertilization applied by soil were associated with greener leaves in contrast to the foliar one.
Other studies dealing with different species report that leaves were greener when plants received organic compared to conventional fertilizers [25].
An important quality aspect is the level of antioxidant compounds, including carotenoids, phenols, and flavonols [14]. In this line, previous studies suggest that organic fertilizers often promote higher content of antioxidant compounds [20,[22][23][24]. In our study, O. microphyllum plants receiving root fertilization showed higher carotenoid content than the ones receiving foliar fertilization with INM-fa being less effective ( Figure 4B). In all cases, fertilization stimulated total phenolic and flavonoid contents, with ChF-fa being less effective in the latter case ( Figure 4C,D). Considering all three antioxidants' contents of O. microphyllum together, ChF-sa and INM-sa stand out as the most stimulatory fertilization schemes applied preferably through the root.

Conclusions
In this pilot field cultivation of Origanum microphyllum, the effect of fertilization on plant growth and herbal quality was evaluated for the first time in an attempt to develop knowledge on cultivation techniques regarding this promising new MAP crop.
In general, the fertilization schemes examined have exerted limited effects on O. microphyllum plant growth and nutrient content; the fertilization when applied through the root system was associated with enhanced visually perceived quality traits (increased green color intensity) and higher medicinal value (increased antioxidant compounds' content). The small size of leaves may have impeded the efficiency of the foliar application. Consequently, the results of this study suggest that the use of soil application upgraded visual quality and herbal antioxidant profile in O. microphyllum. Although the present field study is limited to a single cultivation cycle performed in a pilot way, it however presents first-time baseline values for a range of herbal quality features related to O. microphyllum, thus allowing for the consideration of an enhanced value chain for this neglected and underutilized species.  Table S1: Product details of purchased fertilizers used in the pilot cultivation of Origanum microphyllum in Heraklion, Crete, Greece, Table S2: Experimental field soil properties of the research garden of the Hellenic Mediterranean University used for the fertilization trials on Origanum microphyllum.