Next Article in Journal
Genetic Diversity and Population Structure of Common Bean (Phaseolus vulgaris L.) Landraces in the Lazio Region of Italy
Next Article in Special Issue
Willow (Salix babylonica) Extracts Can Act as Biostimulants for Enhancing Salinity Tolerance of Maize Grown in Soilless Culture
Previous Article in Journal
Domestication Potential of Garcinia kola Heckel (Clusiaceae): Searching for Diversity in South Cameroon
Previous Article in Special Issue
Strawberry Biostimulation: From Mechanisms of Action to Plant Growth and Fruit Quality
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

‘Garlic-lipo’4Plants: Liposome-Encapsulated Garlic Extract Stimulates ABA Pathway and PR Genes in Wheat (Triticum aestivum)

1
Department of Plant Physiology and Plant Ecology, Georgikon Campus, Institute of Agronomy, Hungarian University of Agriculture and Life Sciences, Festetics Str. 7, 8360 Keszthely, Hungary
2
Research Institute for Medicinal Plants and Herbs Ltd., Lupaszigeti Str. 4, 2011 Budakalász, Hungary
3
Department of Biochemistry and Medical Chemistry, University of Pécs Medical School, Szigeti Str. 12, 7633 Pécs, Hungary
4
Department of Information Technology and Its Applications, Faculty of Information Technology, University of Pannonia, Gasparich Str. 18, 8900 Zalaegerszeg, Hungary
5
EduCoMat Ltd., Iskola Str. 12/A, 8360 Keszthely, Hungary
6
Institute of Metagenomics, University of Debrecen, Egyetem Square 1, 4032 Debrecen, Hungary
7
Department of Environmental Biology, Georgikon Campus, Hungarian University of Agriculture and Life Sciences, Festetics Str. 7, 8360 Keszthely, Hungary
8
Department of Molecular Biotechnology and Microbiology, Institute of Biotechnology, Faculty of Science and Technology, University of Debrecen, Egyetem Square 1, 4132 Debrecen, Hungary
*
Author to whom correspondence should be addressed.
Plants 2023, 12(4), 743; https://doi.org/10.3390/plants12040743
Submission received: 8 January 2023 / Revised: 28 January 2023 / Accepted: 31 January 2023 / Published: 7 February 2023
(This article belongs to the Special Issue Plant Biostimulation)

Abstract

:
Recently, environmentally friendly crop improvements using next-generation plant biostimulants (PBs) come to the forefront in agriculture, regardless of whether they are used by scientists, farmers, or industries. Various organic and inorganic solutions have been investigated by researchers and producers, focusing on tolerance to abiotic and biotic stresses, crop quality, or nutritional deficiency. Garlic has been considered a universal remedy ever since antiquity. A supercritical carbon dioxide garlic extract encapsulated in nanoscale liposomes composed of plant-derived lipids was examined as a possible PB agent. The present study focused on the characterization of the genes associated with the pathways involved in defense response triggered by the liposome nanoparticles that were loaded with supercritical garlic extracts. This material was applied to Triticum aestivum in greenhouse experiments using foliar spraying. The effects were examined in a large-scale genome-wide transcriptional profiling experiment by collecting the samples four times (0 min, used as a control, and 15 min, 24 h, and 48 h after spraying). Based on a time-course expression analysis, the dynamics of the cellular response were determined by examining differentially expressed genes and applying a cluster analysis. The results suggested an enhanced expression of abscisic acid (ABA) pathway and pathogenesis-related (PR) genes, of which positive regulation was found for the AP2-, C2H2-, HD-ZIP-, and MYB-related transcription factor families.

Graphical Abstract

1. Introduction

Garlic (Allium sativum L., Amaryllidaceae) is a well-known medicinal herb with a remarkable therapeutic repute of providing many beneficial bioactive compounds [1,2,3]. Among the phyto-components of garlic, there is a high level of fructans [4], flavonoids [5], phenols, and polyphenol content [6,7,8]. The antimicrobial effect of garlic was reported in in vitro experiments that showed an effective response against many plant-pathogenic fungi and bacteria [9,10,11,12,13]. Plant-based extracts have several advantages for crops, such as increased resistance to abiotic and biotic stress and the possibility of reduced application of pesticides and mineral fertilizers [14,15]. Plant biostimulants (PBs) could be effective for achieving these aims by using diverse natural substances that generate physiological and molecular processes in plants to enhance nutritional efficiency, abiotic stress tolerance, and crop quality attributes [16,17]. There are several studies using biostimulants with different chemical and biological ingredients, but there is llimited literature about the use of extracts from garlic as PBs. Hayat et al. used aqueous garlic extracts (AGEs) as a biostimulant for foliar spray and fertilizationto improve crop quality and to promote the physiological potential of cucumber, tomato, eggplant, and pepper seedlings. Moreover, the AGEs exerted priming effects and defense responses against pathogenic fungal infections. This defense mechanism probably activated reactive oxygen species and antioxidant enzymes, such as superoxide dismutase (SOD), and peroxidase (POD) [18,19,20,21]. The application of AGEs has improved the growth and physiological efficiency of fava bean [22] and snap bean plants, and the level of endogenous phytohormones, particularly auxins, gibberellins, and salicylates [23].
Hormonal balance is another factor that has large effects on dormancy release, organ formation, and development. The role of hormones in the dormancy release and sprouting of garlic cloves has already been demonstrated [24,25]. These studies suggested that abscisic acid (ABA) played a major role in the whole differentiation/growth process of garlic plants. Garlic contains large amounts of the compounds described above; therefore, we hypothesize that garlic may have on significant effect when applied to wheat plants and stimulate the physiological processes associated with these compounds. This assumption was examined in this study, which aimed to analyze these relevant processes.
Extraction techniques using supercritical fluids may create greener processes and end products than conventional alternatives [26], leading to efficient solvents with better transport properties (diffusivity, mass transfer coefficient, and penetration ability) than many liquid organic solvents [27]. Supercritical extracts of algae as biostimulants were comprehensively investigated in winter wheat, rapeseed, and mustard [28,29]. Hincapié et al. (2008) [30] reported that using supercritical carbon dioxide (SC-CO2) as a solvent produced garlic extracts, which were tested on Tetranychus urticae. The SC-CO2 garlic extracts do not contain any toxic extraction solvents and preserve oxidation-prone substances because they are not exposed to oxygen or high temperatures during the extraction process [31]. To date, the use of SC-CO2 garlic extracts as PBs in wheat has not been reported.
The nanoscale drug delivery system is a technology used extensively in medicine to increase drug efficiency and to reach target tissues. These modern nanotechnological solutions are considered in agronomy as well [32,33,34]. Karny et al. (2018) [35] reported the use of liposomes composed of plant-derived lipids that penetrated the leaves and translocated to other leaves and to the roots of tomato plants, where active compounds were released. The efficiency of penetration of the applied nanoparticles was up to 33 percent, compared to 0.1 percent of crop-protection agents in a similar treatment. Moreover, airborne nanoparticles disintegrated into safe molecular building blocks (phospholipids) over longer airborne distances to reduce the environmental toll. An example is the liposome-encapsulated plant conditioner product, Elice16indures® of RIMPH Ltd., Budakalász, Hungary (containing 11 herbal CO2 extracts involving garlic), and in both phytotron and earlier field experiments, it was found that this agent triggered the plant defense system through hormonal pathways in Hordeum vulgare [36] and Glycine max [37,38].
In this study, the effect of SC-CO2 garlic extract encapsulated in nanoscale liposomes (preliminarily named as ‘Garlic-lipo’) as a PB was examined. The mechanism of action was investigated by using genome-wide transcriptional profiling of genes associated with the defense response triggered by the ‘Garlic-lipo’. The results presented here are based on the QuantSeq 3′ mRNA sequencing experiments, which data have been previously reported by Kutasy et al. (2022) [39].

2. Materials and Methods

2.1. Garlic-Lipo Stability and Particle Size Measurements of Liposomes

We tested the ‘Garlic lipo’ agent as an active component of ELICE Vakcina®, a member of the Elice16 product line developed in the Research Institute for Medicinal Plants and Herbs (RIMPH) Ltd., Hungary (https://gynki.hu/en/rimph-botanicals/products/ accessed on 1 January 2023), and which formulation and biological effect have already been characterized [36,38]. According to the formulation technology of ELICE Vakcina, the supercritical fluid extraction of fine-grated fresh garlic with natural carbon dioxide (FLAVEX Naturextrakte GmbH, Rehlingen-Siersburg, Germany) was encapsulated in sunflower lecithin-based liposomes.
The electrostatic stabilization, dispersion properties, and nanoparticle size distribution were measured by dynamic light scattering using the Zetasizer Ultra instrument at room temperature (Malvern Panalytica, Malvern, UK). The polydispersity index (PI) and Zeta potential were measured in four technical repeats. The dilution of the samples was 10−1 each time using Milli-Q water. The diluted suspension was placed in a disposable, low-volume cuvette with a path length of 10 mm (Malvern Instruments, Malvern, UK). The size distribution of each liposome sample and the particle size mode were measured using detectors at 3 different angles to account for the front (13°), side (90°), and back (173°) scatter of light at room temperature. Each sample was measured at each angle in triplicate, with adaptive correlation applied to each, to improve overall data quality. The data were processed using the ZS Xplorer Software V.2.3.1 The particle size of the liposomes is given from the cumulative analysis of four measurements of the average Z-diameter and PI.

2.2. Plant Materials and Treatment with ‘Garlic-Lipo’

The seeds of T. aestivum cultivar ‘Cellule’ (France, 2012) were germinated in Petri dishes at 30 °C/20 °C with a 16 h photoperiod, and they were planted in pots and grown under controlled greenhouse conditions. The plants were sprayed with ‘Garlic-lipo’ in a dosage of 240 g·ha−1 using a Euro Pulvé plot sprayer at the BBCH12 stage [40]. Fresh leaves were collected before the treatment as control (0) and at 15 min, 24 h, and 48 h after treatment (Figure 1). The samples were stored in a RNA Shield (Zymo Research, Irvine, CA, USA) at −25 °C until RNA analysis. Three biological repeats were used for further analyses.

2.3. Sequencing, Data Processing, Gene-Level Quantification, and DEG Determination

The Illumina Gene Expression Profiling (GEx) library preparation, RNA sequencing, sequence read pre-processing, gene-level quantification, and differentially expressed gene (DEG) analysis were performed as described in our previous papers [39,41]. The RNA libraries were sequenced with a final output single-end, with 14–26 M x 85 bp long reads. The raw reads were deposited in the National Center for Biotechnology Information (NCBI) database (as the BioProject PRJNA808851) and the Sequence Read Archive (SRA) database (SRR18107544, SRR18107543, SRR18107542, and SRR18107541). A de novo reference transcript was created using these SRA datasets to reconstruct transcripts that represent the experimental samples. This transcript dataset was deposited in the Transcriptome Shotgun Assembly (TSA) database at DBJ/EMBL/GenBank under the accession GJUY00000000.1. To estimate transcript abundances, CountTable was created [39], which is accessible in Mendeley Data (https://data.mendeley.com/datasets/p66v4yxbtp/1, accessed on 17 March 2022). DEGs were calculated based on the CountTable by applying a pairwise differential expression analysis as described by Hegedűs et al. (2022) [41].

2.4. Time-Course Expression Analysis

A time-course expression analysis of the count data arising from the RNA-seq was performed to detect genomic features with significant temporal expression changes and with significant differences between the experimental groups [39]. The genes with significantly different expression levels were divided into nine clusters. The TimeCourse Table is accessible in Mendeley Data (https://data.mendeley.com/datasets/xvvscxpz6w/1, accessed on 18 March 2022). According to the expression dynamics of the changes of the clustered genes, these nine clusters were classified into four groups: (1) expression level gradually increased; (2) expression level gradually decreased; (3) expression level initially increased and then reduced; and (4) expression level initially decreased and then returned to its original level. The clusters were annotated with Gene Ontology (GO) terms and visualized [39]. Group 1 genes with gradually increased expression levels were further examined and discussed in this study.

2.5. Functional Annotation

The methods of functional annotation and GO analysis were described in a previous publication [37]. The annotation of the entire transcriptome is presented in the AnnotationTable that is accessible in Mendeley Data (https://data.mendeley.com/datasets/p66v4yxbtp/1, accessed on 17 March 2022) and described in Kutasy et al. (2022) [39].

2.6. Real-Time PCR Analysis

Total RNA was quantified using a Nanodrop 2000c spectrophotometer and a Qubit 2.0 fluorometer (Thermo Fischer Scientific, Waltham, MA, USA). Five hundred nanograms (ng) of total RNA was reverse transcribed using the M-MuLV RT (Maxima First Strand cDNA Synthesis Kit, Thermo Fischer Scientific, Waltham, MA, USA). A total of 50 ng of cDNA was used in 10 μL reactions for real-time PCR using the Xceed qPCR SG 2× Mix (Institute of Applied Biotechnologies, Praha-Strašnice, Czech Republic) and a CFX384 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). The data (Tables S3 and S4) were analyzed using the ΔΔCt method [42]. Tubulin and ubiquitin expressions were used as a reference for gene expression.

3. Results

To date, the effect of garlic extract was investigated in vegetable plants in studies that focused on yield quality and physiological function, as well as stress responses. However, it has not yet been examined in monocotyledonous plants. Moreover, the effect of garlic extract on changes in gene expression levels has not been studied.
High levels of flavonoids, phenolics, fructans, and abscisic acid hormone (ABA) were found in garlic. Therefore, in this study, the processes related to these compounds were examined, which might be stimulated by garlic extract.
This study presents the effects of garlic SC-CO2 extract on wheat plants by investigating ‘Garlic-lipo’. The formulation of this agent was carried out with a focus on better uptake by plants after treatment, with a strict analytical control measuring the size of the nanoparticle and the Zeta potential. The mechanism of this effect was investigated using transcriptome analysis, such as mRNA-sequencing, determination of DEGs, and time-course expression analysis, showing high expression changes in response to the treatment. The in silico analysis of mRNA sequencing and gene expression was validated using RT-qPCR methods.

3.1. Dynamic Light Scattering

The size determination parameters are summarized in Table 1, and representative figures of the measurements are shown in Figure 2. The Zeta potential analysis is summarized in Table 2 and Figure 3. The liposome properties of the investigated samples indicated that either the particle size or the polydispersity index (p > 0.05, Pearson) was relatively constant (149.2 nm in size, polydispersity index = 0.13035, and n = 4) in the four measured samples.

3.2. Pairwise Expression Analysis and Determination of Top 50 DEGs

The expression levels of DEGs measured at three sampling times (15 min/24 h/48 h) after the ‘Garlic-lipo’ treatment were compared to the control (0 min) for the largest difference in the top 50 genes presented in the heatmaps (Figure 4, Figures S1 and S2). Many differentially expressed (DE) features were detected in the 15 min sample [(probability > 0.9): 1712; up-regulated (M > 0): 919; down-regulated (M < 0): 793]; in the 24 h sample [(probability > 0.9): 2485; up-regulated (M > 0): 1279; down-regulated (M < 0): 1206]; and in the 48 h sample [(probability > 0.9): 2578; up-regulated (M > 0): 865; down-regulated (M < 0): 1713].
Proteins involved in stress response and the ABA pathway showed up-regulated levels in the 15 min sample (Figure S1), such as ABA-responsive protein HVA22, late embryogenesis abundant protein (LEA), 9-cis-epoxycarotenoid dioxygenase (NCED3), and U-box domain-containing protein 19 (PUB19). The down-regulated genes were found to be associated with photosynthesis, as well as ribosomal RNA genes.
The top 50 DEGs were mostly down-regulated in the 24 h sample (Figure S2), except for pathogenesis-related protein 1 (PR1) and aldehyde dehydrogenase family 7 member A1 (ALDH7). Photosynthesis-related and ribosome-related genes were shown to be down-regulated. The largest difference was found in the 48 h sample, which is presented in Figure 4.
While there were several ABA-related genes determined in the 15 min sample and the 24 h sample, only one was found in the 48 h sample. As for up-regulated genes, the expression level of flavanone 3-hydroxylase (F3H) gene associated with the ABA pathway was increased. Moreover, PR1-5 genes and other stress-induced genes were proven to be up-regulated. The up-regulated PR genes, such as PR1, endo-beta-1,3-glucanase (PR2), chitinas (PR3), wheatwin-2 protein (PR4), and thaumatin protein (PR5), were significantly different in their expressions. One of the PRs, the Bowman–Birk-type trypsin inhibitor (PR6), showed a down-regulated expression level.
Biotic and abiotic stress response genes, such as alpha-amylase/trypsin inhibitor (ATI), caseinolytic protease B1 (ClpB1), wheat-induced resistance 1B (WIR1B), and heavy metal-associated isoprenylated plant protein 27 (HIPP27), were up-regulated in the sample of 48 h. Different transferases, such as UDP-glucose flavonoid 3-O-glucosyltransferase 7 (GT7), tricetin 3′,4′,5′-O-trimethyltransferase (OMT2), and acetylserotonin O-methyltransferase 1 (ASMT), were shown to be up-regulated level in the same sample.
Moreover, biotic and abiotic stress-responding transcription factors (TFs) that consist of HD-ZIP, C2H2, and AP2 families, such as zinc finger protein 12 (ZAT12—C2H2 family), homeobox-leucine zipper transcription factors 22 and 24 (HOX22 and HOX24—HD-ZIP family), dehydration-responsive element-binding protein 1G (DREB1G—AP2 family), and MYBS1—MYB-related family, were observed in each sample evaluated at all of the three timepoints of this study.
The determination of the top 50 DEGs demonstrated the changes in the expression levels of ABA pathway-related genes, PRs, and TF genes (Figure 5). To obtain more information about the physiological processes in wheat after the biostimulation treatment applied here, the samples were examined using a time-course expression analysis.

3.3. Time-Course Expression Analysis

A time-course expression analysis was an appropriate tool for evaluating the data arising from the RNA-seq technology to represent genes with significant changes in their temporal expressions and with significant differences between the experimental groups. Using TSA and CountTable data, 5287 contigs were tested and filtered, which resulted in nine clusters. The clusters were classified into four groups according to their dynamics of changes: those of which the expression level (i) gradually increased, (ii) gradually decreased, (iii) initially increased and then reduced, or (iv) initially decreased and then returned to its original level. The annotated genes were examined by using the GO terms and biological processes that were already presented in a previous paper (Kutasy 2022). The group of genes with expression levels that gradually increased ((i) cluster 1) was shown mainly to be associated with defense response processes, while the groups of genes with expression levels that gradually decreased ((ii) clusters 3, 4, 6, and 7) were shown to have down-regulated cellular biosynthetic and metabolic processes. Moreover, the groups of genes with expression levels that initially increased and then reduced ((iii) clusters 2, 8, and 9) were mostly shown to be involved in different transport and translation processes. The group of genes with expression levels that initially decreased and then returned to their original levels ((iv) cluster5) was associated with RNA-related (transcription) and metabolic processes (Kutasy 2022). The genes of the nine clusters were also classified using the Kyoto Encyclopedia for Genes and Genomes (KEGG) analysis (Table S1). In this study, the genes with expression levels that gradually increased showed defense response associations (Figure 6).
The time-course analysis identified 24 pathways, with a total of 201 sequences (Tables S1 and S2) in (i) cluster 1, and enzymes that may be involved in primary defense reactions. Regarding defense response, plant–pathogen interaction; mitogen-activated protein kinase (MAPK) signaling pathway; plant hormone signal transduction; flavonoid, terpenoid, and phenylpropanoid biosynthesis; glutathione metabolism; and metabolism by cytochrome P450 were found to be the most important pathways (Figure 6). The disease-resistance protein PR1 was involved in the pathways of plant–pathogen interaction (map04626), plant hormone signal transduction (map04075), and MAPK signaling (map04016), while the RPM1 protein was found to be associated with the plant–pathogen interaction (map04626) pathway. The mono- and diterpenoid biosynthesis (map00902, 00904) pathways were identified in all of the samples, involving the proteins of 7-deoxyloganetin glucosyltransferase (UGT85A23-24), ent-kaurene oxidase (CYP701A), and ent-sandaracopimaradiene 3-hydroxylase (CYP701A8). In the phenylpropanoid biosynthesis (map00940) pathway, the proteins of shikimate O-hydroxycinnamoyl-transferase (HST), caffeic acid 3-O-methyltransferase/acetylserotonin O-methyltransferase (COMT), and peroxidase were found. The cytochrome P450 and glutathione metabolic (map00980, 00480) pathways were also identified, involving the proteins of glutathione S-transferases (GST23 and GSTU6) and L-ascorbate peroxidase (APX).
Based on the time-course expression analyses of (ii) clusters 3, 4, 6, and 7, the down-regulated cellular biosynthetic and metabolic processes may be involved in photosynthesis-associated processes; the processes of peroxisome, glyoxylate, and dicarboxylate metabolisms; linoleic acid metabolism; oxidative phosphorylation; and ribosomal processes. Regarding the plant–pathogen interaction pathway, heat shock proteins were observed.
The time-course analysis of (iii) clusters 2, 8, and 9 showed enzymes that may be involved in process of endocytosis; MAPK signaling pathway; and cysteine and methionine metabolisms. These clusters are also associated with ribosomal processes in the cytosol.
Moreover, the time-course analysis of (iv) cluster 5 was associated with mRNA processing, RNA phosphodiester-bound hydrolysis, and nucleic acid metabolic processes. The KEGG analysis of cluster 5 gave no results.

3.4. Real-Time PCR Analysis

The most important defense response-relevant genes of the top 50 DEGs were examined using real-time PCR. PR1-5 showed an increased expression level during the whole experimental period as it was the highest in the 48 h sample—as it was observed in the transcriptome analysis. The genes associated with the ABA pathway (NCED3, HVA22, PUB19, and LEA31) and TFs (HOX24, ZAT12, and MYBS1) also showed similar changes in their expression levels as revealed in the transcriptome analysis. Typically, these genes showed maximal expression levels in the 15 min sample and the 24 h sample. Moreover, a TF (DREB1G) was prognosticated to have the highest expression level in the 48 h sample in the DEG analysis, and this was confirmed by the real-time PCR (Figure 7 and Tables S3 and S4).

4. Discussion

Worldwide, the scientific community and commercial enterprises have focused on the development of environmentally safe methods that facilitate plant growth and crop protection. PBs could be effective chievinging these aims as they are diverse natural substances that generate physiological and molecular processes in plants to enhance nutritional efficiency, abiotic and biotic stress tolerances, and crop quality attributes [16,17,43]. Using plant-based extracts has several advantages, such as enhanced yields and quality of crops, low toxicity to humans and the environment, increased resistance of crops to abiotic and biotic stresses, and reduction in the application of pesticides and mineral fertilizers [14,15]. Aqueous garlic extract as a plant biostimulant was recently reported as a promising agent to trigger plant immunity [18]. Focusing on low-input sustainable agriculture, we aimed to reinforce the efficacy of garlic extract using SC-CO2 extraction and liposome nanotechnology for better uptake in and delivery to plants. A novel ‘Garlic-lipo’ production and its treatment on wheat plants were designed. The time-course effect of ‘Garlic-lipo’ was tested by using whole-genome transcriptional profiling and enzyme activity measurements. The up-regulated genes were validated by using RT-qPCR experiments.
The liposome size determination and the monitoring of colloid characteristics during the ‘Garlic-lipo’ production were controlled using dynamic light scattering. The liposome sizes were in the range between 134.2 and 164.2 nm, and the colloid stability proved to be high, showing a Zeta potential between −40.66 and −40.53 mV. These physicochemical properties of the nanoliposomes provide the product a high stability and great bioactivity [44]. However, there are only a few reports on the application of such a product in agricultural practice.
This product was tested using time-course experiments. The QuantSeq 3′ mRNA sequencing for RNA quantification in the samples had been performed earlier (data reported in Kutasy 2022 [39]). The determination of DEGs showed an increase in the expression levels of genes associated with the ABA pathway, which was observed especially 15 min after treatment. A parallel can be drawn with the remarkable amounts of endogenous phytohormones described in garlic clove extracts, among which a high level of ABA was found by Arguello (1991) [24]. The cultivar, photoperiod, and temperature affected the hormone levels when ABA showed up to 10–25-fold higher levels than cytokinins, auxins, and jasmonic acid [45]. A treatment with a high ABA-content garlic extract, therefore, may stimulate plant responses to biotic or abiotic stresses and increase yield levels, as was shown by previous exogenous ABA treatments [46,47,48]. Based on this information, it was hypothesized that a garlic extract with a high ABA content might stimulate ABA pathway-related genes. This theory was proven here in these experiments. In these experiments, we observed a rapid up-regulation of these genes at 15 min after the ‘Garlic-lipo’ treatment, but their expressions gradually decreased when the samples were taken 24 h after the treatment. Among the top 50 DEGs (15 vs. 0 min), the homologous protein of stress-and ABA-induced HVA22, LEA, and NCED3 is the key rate-limiting enzyme in the ABA biosynthetic pathway regulating multi-abiotic stress tolerance, as was found in [49,50,51,52,53]. The regulators of salt-tolerant [54] and negative regulators of drought stress-response protein [55] and PUB19 [56] were also up-regulated, which were also known to be sensitive to ABA in Arabidopsis thaliana [57]. The down-regulated genes were associated with photosynthesis, as well as ribosomal RNA genes. A down-regulation of photosynthetic genes by ABA signaling was observed by Staneloni et al. (2008) [58] and Gao et al. (2016) [59], and this phenomenon might be associated with ABA-mediated chlorophyll degradation and leaf senescence.
The regulation of ABA-signaling genes was associated with the overexpression of TFs, with their highest values observed at 15 min. TFs are members of the HD-ZIP, bZIP, C2H2, and AP2 families, such as ZAT12 (C2H2 family), that participate in the induction and repression of cold-responsive and oxidative stress genes [60]; HOX22 (HD-ZIP family) is involved in the defense mechanism of pathogen-related stress response of fungi [61], and drought- [62] and heat-stress-tolerance-induced HOX24 [63] and MYBS1 [64] are the essential components of the sugar signaling pathway regulating the α-amylase gene promoter [65]. A high expression level was found for DREB1G (AP2 family), which is a typical C-repeat-binding factor (CBF/DREB1) whose specific functions was found to be related to cold stress response [66,67]. Interestingly, this TF was the only one present when the sample was taken 48 h after the treatment. A microarray analysis using Arabidopsis cultured cells revealed that several ABA-inducible genes and the TF ZAT12 are induced by oxidative stresses, which might be a regulator in the reactive oxygen species (ROS) scavenging mechanism that is involved in multiple stresses, such as wounding, pathogen infection, and abiotic stresses [68].
The CBF cold-responsive pathway includes the three key TFs of CBF1, CBF2, and CBF3, and this pathway also plays a significant role in the cold acclimation of Arabidopsis [69]. The CBF proteins as transcriptional activators recognize both the C-repeat (CRT) and the DREB. It was shown that the overexpression of CBF2 and CBF3 can induce several cold-responsive (COR) genes [70,71].
Vogel et al. (2005) [60] reported that ZAT12 was also induced with CBFs, which regulated the expression of 24 cold-responsive genes, indicating that ZAT12 might be involved in cold acclimation via a novel cold-response pathway [72,73].
The up-regulated PR genes were found both in the 24 h sample and in the 48 h sample after the ‘Garlic-lipo’ treatment, of which they showed a higher expression level in the 48 h sample. The top 50 DEGs of 48 vs. 0 h indicated significantly higher expressions of PR genes, such as PR1, which are known as the hallmarks of defense pathways in wheat [74,75]. PR2 is involved in the resistance response of wheat to stripe rust (Puccinia striiformis f. sp. tritici) [76,77]; PR3 constitutes the second largest group of antifungal proteins [78]; PR4 has antifungal activity against several pathogenic fungi (Fusarium sp., Botrytis sp.) [79,80]; and PR5 is mostly associated with responses to biotic stresses, especially drought and osmotic stresses [81]. Additionally, F3H, which is an ABA-induced abiotic stress gene that catalyzes the third step of the central flavonoid biosynthetic pathway [82,83], was found to be up-regulated in this sample. One of the PRs, the PR6, showed a down-regulated expression level [84,85].
To summarize, the majority of the top 50 DEGs might be grouped into (i) ABA-pathway and (ii) PR genes; however, their induction was strongly associated with time. A rapid response to the drug was observed in ABA signaling, but PR induction was detected after 24–48 h of the treatment. It was presumed that the regulation of these genes was under the control of (iii) TFs, which were not under strict time control, such as ABA and PR. The observed time-shift effect of PR genes is characteristic of priming inducers that might be triggered by fungus or chemical treatments [85,86]. The expression dynamics of PR1, PR3, PR4, PR5, PR6, NCED3, HVA22, PUB19, LEA 31, HOX24, DREB1G, MYBS1, and ZAT12 were validated by the RT-qPCR as well.
Additionally, biotic and abiotic stress-response genes were up-regulated in the 48 h sample, such as ATI, which has a defense mechanism against insect attack [87]; ClpB1, which is a high-molecular-weight chaperon protein that is part of the heat shock protein 100 (HSP 100) family [88]; WIR1B, which is a small glycine- and proline-rich protein having a key role in defense reactions against fusarium head blight (FHB) infection caused by the fungus Fusarium graminearum [89] and in the plant–pathogen interactions of Blumeria graminis [90]; and HIPP27, which is a metallochaperone that has been reported to play an important role in cadmium detoxification [91] and in response to nematode infection [92]. Different transferases showed up-regulated levels in the 48 h sample, such as GT7, which primarily transports sugar molecules to aglycones to enhance solubility and stability during flavonol biosynthesis [93]; OMT, which catalyzes three sequential methylations of flavone tricetin; and ASMT1, which is responsible for the syntheses of N-acetylserotonin and melatonin. Melatonin is synthesized from tryptophan and plays an important role in enhancing the resistance of plants against biotic and abiotic stresses [94]. Melatonin treatment significantly enhanced the drought tolerance of wheat seedlings by decreasing membrane damage [95]. The down-regulated genes were genes associated with photosynthesis and ribosomal RNA genes.
Since the effects of ‘Garlic-lipo’ showed a strong time factor relation, a time-course gene expression analysis was performed to determine other gene clusters induced by the treatment effects. This investigation resulted in nine clusters, in which the genes were annotated by determining their GO terms and using a KEGG analysis (Table S1). Both the GO and KEGG analyses showed an increase in the expression levels of genes associated with defense response processes. During the experimental period, the genes involved in the biosynthesis of PR1-5 genes, terpenoids, and phenylpropanoids operated at a consistently high level. PRs have been found to be activated in response to different biotic and abiotic stresses, so PR1-5 proteins have a wide range of functions, such as protease inhibitor, glucanase, chitinase, and thaumatin protein [96]. Furthermore, the CYP450 and glutathione transferase genes involved in stress response showed a high level of expression. The GST genes are highly inducible by a wide range of stress conditions, including biotic stress [97].
A down-regulation of genes involved in photosynthesis was observed based on a heat map and KEGG analysis. For the entire period, the down-expression of the processes involved in photosynthesis was detected using time-course GO analysis, but it was not one of the most prominent processes (Figure 2 in Kutasy et al. (2022) [39]). Moreover, the KEGG analysis showed that the process of restoring photosynthesis had already begun. In the porphyrin and chlorophyll metabolism pathway (map00860) and the carbon fixation pathway in photosynthetic organisms (map00710), the expression of chlorophyllase (CLH) and malate dehydrogenase (MDH) were found to be up-regulated.
To summarize, the main part of the top 50 up-regulated DEGs affecting plant development and pathogen defense mechanism was found after the ‘Garlic-lipo’ treatment. These genes suggest a strong regulation of the ABA pathway and resistance mechanism. This phenomenon indicates physiological relations, with the synthesis of the plant hormone ABA having a key role in plant development, control of organ size, stomal closure, and signaling under different plant stresses. The ‘Garlic-lipo’ may rapidly affect the ABA signaling pathway in plants in response to nanoparticle exposure, leading to enhanced ABA-related gene activity and defense system. The assumed phytopathogen activity of this agent was also demonstrated to trigger PR genes involving systemic acquired resistance (PR1, PR2, and PR5) and local acquired resistance (PR3, PR4). Additionally, the common upregulated DEGs of the investigated samples found in the KEGG pathway analysis were involved in monoterpenoid, phenylpropanoid, diterpenoid, and stilbenoid biosynthesis, and glutathione and phenylalanine metabolism. These pathways have a key role in plant development, cell cycle, stem growth regulation, induced germination, defense mechanisms, and disease resistance (Figure 8). Based on all of these results, it can be concluded that the transcriptomic changes, which were examined using several approaches, unanimously indicate that the use of ‘Garlic-lipo’ may improve the health and development of plants by strengthening their response to biotic and abiotic stresses involving an emerging pathogen attack.

5. Conclusions

These results highlight the significant effects of the SC-CO2 garlic extract as the up-regulations of several defense-response genes were demonstrated. The increase in the expressions of the ABA pathway and PR genes was shown to be associated mainly with defense mechanisms. In conclusion, the plant’s innate immunity and defense response mechanisms might be triggered systemically by this agent, suggesting a potential priming active substance. Due to its extraction process and nanotechnological formulation, ‘Garlic-lipo’ may be a promising PB that can support sustainable agriculture. To substantiate the effects described here, we are planning to conduct further infection tests and field experiments.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants12040743/s1, Figure S1: Heat map of the top 50 DEGs based on the comparison of the 15 min vs. control samples; Figure S2: Heat map of the top 50 DEGs based on the comparison of the 24 h vs. control samples; Table S1: The significantly expressed genes based on time-course analysis; Table S2: The result of time-course expression and KEGG analysis in nine clusters; Table S3: Primers of real-time PCR; Table S4: Real-time PCR datasets.

Author Contributions

Conceptualization, E.V.; methodology, B.K., K.D., G.H. and E.V; preparation of the investigated material ‘Garlic-lipo’, M.K.; software, G.H. and E.V.; validation, B.K. and K.D.; formal analysis, B.K., G.H. and E.V.; project administration, E.V.; RT-qPCR investigation, N.K.; data curation, B.K.; writing—original draft preparation, B.K.; writing—review and editing, Z.P.A. and E.V.; visualization, B.K.; supervision, E.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work was founded by the KFI_16-1-2017-0457—Development and production of a plant-based pesticide-plant conditioner for use in organic farming—a project of the Hungarian Government. This work was also supported by the GINOP-2.3.4-15-2016-00002.

Data Availability Statement

Plant sample from SC-CO2 garlic extract-treated Triticum aestivum. Published: Mendeley data, 17 March 2022|Version 1| DOI: 10.17632/p66v4yxbtp.1 Contributors: Barbara Kutasy, Kincső Decsi and Eszter Virág https://data.mendeley.com/datasets/p66v4yxbtp/1 (accessed on 17 March 2022). TimeCourse Table of plant sample from SC-CO2 garlic extract-treated Triticum aestivum. Published: Mendeley data, 18 March 2022|Version 1| DOI: 10.17632/xvvscxpz6w.1 Contributors: Barbara Kutasy, Kincső Decsi and Eszter Virág https://data.mendeley.com/datasets/xvvscxpz6w/1 (accessed on 18 March 2022). Triticum aestivum cultivar ‘Cellule’ (bread wheat). Published: National Center for Biotechnology Information (NCBI) database as BioProject PRJNA808851 (Sequence Read Archive (SRA) database: SRR18107544, SRR18107543, SRR18107542, and SRR18107541; Transcriptome Shotgun Assembly (TSA) database: GJUY00000000.1.); Registration date: 21 February 2022 Contributors: Barbara Kutasy, Kincső Decsi and Eszter Virág. https://www.ncbi.nlm.nih.gov/bioproject?LinkName=biosample_bioproject&from_uid=26113471 (accessed on 21 February 2022).

Acknowledgments

We express our thanks to József Péter Pallos, RIMPH Ltd. Budakalász, Hungary, for project financing and to Zsófia Thomas-Nyári for administration. We express our thanks to Plant-Art Research Ltd. Tata, Hungary for providing the plants grown in this study and the treatments and to Xenovea Ltd., Szeged, Hungary, for performing the NGS sequencing. We are grateful to the editor and to the reviewers for their valuable comments that have helped improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lawson, L.D. Garlic: A review of its medicinal effects and indicated active compounds. Blood 1998, 179, 62. [Google Scholar] [CrossRef]
  2. Diretto, G.; Rubio-Moraga, A.; Argandoña, J.; Castillo, P.; Gómez-Gómez, L.; Ahrazem, O. Tissue-specific accumulation of sulfur compounds and saponins in different parts of garlic cloves from purple and white ecotypes. Molecules 2017, 22, 1359. [Google Scholar] [CrossRef]
  3. Divya, B.; Suman, B.; Venkataswamy, M.; Thyagaraju, K. A study on phytochemicals, functional groups and mineral composition of Allium sativum (garlic) cloves. Int. J. Curr. Pharm. Res. 2017, 9, 42–45. [Google Scholar] [CrossRef]
  4. Muir, J.G.; Shepherd, S.J.; Rosella, O.; Rose, R.; Barrett, J.S.; Gibson, P.R. Fructan and free fructose content of common Australian vegetables and fruit. J. Agric. Food Chem. 2007, 55, 6619–6627. [Google Scholar] [CrossRef] [PubMed]
  5. Miean, K.H.; Mohamed, S. Flavonoid (myricetin, quercetin, kaempferol, luteolin, and apigenin) content of edible tropical plants. J. Agric. Food Chem. 2001, 49, 3106–3112. [Google Scholar] [CrossRef]
  6. Lanzotti, V. The analysis of onion and garlic. J. Chromatogr. A 2006, 1112, 3–22. [Google Scholar] [CrossRef]
  7. Gorinstein, S.; Drzewiecki, J.; Leontowicz, H.; Leontowicz, M.; Najman, K.; Jastrzebski, Z.; Zachwieja, Z.; Barton, H.; Shtabsky, B.; Katrich, E. Comparison of the bioactive compounds and antioxidant potentials of fresh and cooked Polish, Ukrainian, and Israeli garlic. J. Agric. Food Chem. 2005, 53, 2726–2732. [Google Scholar] [CrossRef]
  8. Nencini, C.; Menchiari, A.; Franchi, G.G.; Micheli, L. In vitro antioxidant activity of aged extracts of some Italian Allium species. Plant Foods Hum. Nutr. 2011, 66, 11–16. [Google Scholar] [CrossRef]
  9. Perello, A.E.; Noll, U.; Slusarenko, A.J. In vitro efficacy of garlic extract to control fungal pathogens of wheat. J. Med. Plants Res. 2013, 7, 1809–1817. [Google Scholar] [CrossRef]
  10. Mougou, I.; Boughalleb-M’hamdi, N. Biocontrol of Pseudomonas syringae pv. syringae affecting citrus orchards in Tunisia by using indigenous Bacillus spp. and garlic extract. Egypt. J. Biol. Pest Control 2018, 28, 60. [Google Scholar] [CrossRef] [Green Version]
  11. Sarfraz, M.; Nasim, M.J.; Jacob, C.; Gruhlke, M.C. Efficacy of allicin against plant pathogenic fungi and unveiling the underlying mode of action employing yeast based chemogenetic profiling approach. Appl. Sci. 2020, 10, 2563. [Google Scholar] [CrossRef]
  12. Slusarenko, A.J.; Patel, A.; Portz, D. Control of plant diseases by natural products: Allicin from garlic as a case study. Sustain. Dis. Manag. A Eur. Context 2008, 121, 313–322. [Google Scholar] [CrossRef]
  13. Curtis, H.; Noll, U.; Störmann, J.; Slusarenko, A.J. Broad-spectrum activity of the volatile phytoanticipin allicin in extracts of garlic (Allium sativum L.) against plant pathogenic bacteria, fungi and Oomycetes. Physiol. Mol. Plant Pathol. 2004, 65, 79–89. [Google Scholar] [CrossRef]
  14. Godlewska, K.; Ronga, D.; Michalak, I. Plant extracts-importance in sustainable agriculture. Ital. J. Agron. 2021, 16. [Google Scholar] [CrossRef]
  15. Rouphael, Y.; Colla, G. Biostimulants in agriculture. Front. Plant Sci. 2020, 11, 40. [Google Scholar] [CrossRef]
  16. Du Jardin, P. Plant biostimulants: Definition, concept, main categories and regulation. Sci. Hortic. 2015, 196, 3–14. [Google Scholar] [CrossRef]
  17. Calvo, P.; Nelson, L.; Kloepper, J.W. Agricultural uses of plant biostimulants. Plant Soil 2014, 383, 3–41. [Google Scholar] [CrossRef]
  18. Hayat, S.; Ahmad, H.; Ali, M.; Hayat, K.; Khan, M.A.; Cheng, Z. Aqueous garlic extract as a plant biostimulant enhances physiology, improves crop quality and metabolite abundance, and primes the defense responses of receiver plants. Appl. Sci. 2018, 8, 1505. [Google Scholar] [CrossRef]
  19. Hayat, S.; Ahmad, H.; Ali, M.; Ren, K.; Cheng, Z. Aqueous garlic extract stimulates growth and antioxidant enzymes activity of tomato (Solanum lycopersicum). Sci. Hortic. 2018, 240, 139–146. [Google Scholar] [CrossRef]
  20. Hayat, S.; Cheng, Z.; Ahmad, H.; Ali, M.; Chen, X.; Wang, M. Garlic, from remedy to stimulant: Evaluation of antifungal potential reveals diversity in phytoalexin allicin content among garlic cultivars; allicin containing aqueous garlic extracts trigger antioxidants in cucumber. Front. Plant Sci. 2016, 7, 1235. [Google Scholar] [CrossRef] [Green Version]
  21. Olusanmi, M.; Amadi, J. Studies on the antimicrobial properties and phytochemical screening of garlic (Allium sativum) extracts. Ethnobot. Leafl. 2010, 2009, 10. [Google Scholar]
  22. Mohamed, M.H.; Badr, E.A.; Sadak, M.S.; Khedr, H.H. Effect of garlic extract, ascorbic acid and nicotinamide on growth, some biochemical aspects, yield and its components of three faba bean (Vicia faba L.) cultivars under sandy soil conditions. Bull. Natl. Res. Cent. 2020, 44, 100. [Google Scholar] [CrossRef]
  23. Elzaawely, A.A.; Ahmed, M.E.; Maswada, H.F.; Al-Araby, A.A.; Xuan, T.D. Growth traits, physiological parameters and hormonal status of snap bean (Phaseolus vulgaris L.) sprayed with garlic cloves extract. Arch. Agron. Soil Sci. 2018, 64, 1068–1082. [Google Scholar] [CrossRef]
  24. Arguello, J.; Ledesma, A.; Bottini, R. Hormonal regulation of dormancy in garlic (Allium sativum L.) cv Rosado Paraguayo. Agriscientia 1991, 8, 9–14. [Google Scholar] [CrossRef]
  25. Rohkin Shalom, S.; Gillett, D.; Zemach, H.; Kimhi, S.; Forer, I.; Zutahy, Y.; Tam, Y.; Teper-Bamnolker, P.; Kamenetsky, R.; Eshel, D. Storage temperature controls the timing of garlic bulb formation via shoot apical meristem termination. Planta 2015, 242, 951–962. [Google Scholar] [CrossRef] [PubMed]
  26. Beckman, E.J. Supercritical and near-critical CO2 in green chemical synthesis and processing. J. Supercrit. Fluids 2004, 28, 121–191. [Google Scholar] [CrossRef]
  27. Brunner, G. Supercritical fluids: Technology and application to food processing. J. Food Eng. 2005, 67, 21–33. [Google Scholar] [CrossRef]
  28. Michalak, I.; Chojnacka, K.; Dmytryk, A.; Wilk, R.; Gramza, M.; Rój, E. Evaluation of supercritical extracts of algae as biostimulants of plant growth in field trials. Front. Plant Sci. 2016, 7, 1591. [Google Scholar] [CrossRef]
  29. Michalak, I.; Chojnacka, K.; Saeid, A. Plant growth biostimulants, dietary feed supplements and cosmetics formulated with supercritical CO2 algal extracts. Molecules 2017, 22, 66. [Google Scholar] [CrossRef]
  30. Hincapié, C.A.L.; López, G.E.P.; Torres, R.C. Comparison and characterization of garlic (Allium sativum L.) bulbs extracts and their effect on mortality and repellency of Tetranychus urticae Koch (Acari: Tetranychidae). Chil. J. Agric. Res. 2008, 68, 317–327. [Google Scholar] [CrossRef]
  31. Bravi, E.; Perretti, G.; Falconi, C.; Marconi, O.; Fantozzi, P. Antioxidant effects of supercritical fluid garlic extracts in sunflower oil. J. Sci. Food Agric. 2017, 97, 102–107. [Google Scholar] [CrossRef] [PubMed]
  32. Ye, Y.; Cota-Ruiz, K.; Hernández-Viezcas, J.A.; Valdés, C.; Medina-Velo, I.A.; Turley, R.S.; Peralta-Videa, J.R.; Gardea-Torresdey, J.L. Manganese nanoparticles control salinity-modulated molecular responses in Capsicum annuum L. through priming: A sustainable approach for agriculture. ACS Sustain. Chem. Eng. 2020, 8, 1427–1436. [Google Scholar] [CrossRef]
  33. Latef, A.A.H.A.; Alhmad, M.F.A.; Abdelfattah, K.E. The possible roles of priming with ZnO nanoparticles in mitigation of salinity stress in lupine (Lupinus termis) plants. J. Plant Growth Regul. 2017, 36, 60–70. [Google Scholar] [CrossRef]
  34. Dutta, P. Seed priming: New vistas and contemporary perspectives. In Advances in Seed Priming; Springer: Singapore, 2018; pp. 3–22. [Google Scholar]
  35. Karny, A.; Zinger, A.; Kajal, A.; Shainsky-Roitman, J.; Schroeder, A. Therapeutic nanoparticles penetrate leaves and deliver nutrients to agricultural crops. Sci. Rep. 2018, 8, 7589. [Google Scholar] [CrossRef] [PubMed]
  36. Hegedűs, G.; Kutasy, B.; Kiniczky, M.; Decsi, K.; Juhász, Á.; Nagy, Á.; Pallos, J.P.; Virág, E. Liposomal Formulation of Botanical Extracts may Enhance Yield Triggering PR Genes and Phenylpropanoid Pathway in Barley (Hordeum vulgare). Plants 2022, 11, 2969. [Google Scholar] [CrossRef]
  37. Decsi, K.; Kutasy, B.; Kiniczky, M.; Hegedűs, G.; Virág, E. RNA-seq datasets of field soybean cultures conditioned by Elice16Indures® biostimulator. Data Brief 2022, 42, 108182. [Google Scholar] [CrossRef]
  38. Decsi, K.; Kutasy, B.; Hegedűs, G.; Alföldi, Z.P.; Kálmán, N.; Nagy, Á.; Virág, E. Natural immunity stimulation using ELICE16INDURES® plant conditioner in field culture of soybean. Heliyon 2023, 9, e12907. [Google Scholar] [CrossRef] [PubMed]
  39. Kutasy, B.; Decsi, K.; Kiniczky, M.; Hegedűs, G.; Virág, E. Time-course gene expression profiling data of Triticum aestivum treated by supercritical CO2 garlic extract encapsulated in nanoscale liposomes. Data Brief 2022, 42, 108287. [Google Scholar] [CrossRef]
  40. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  41. Hegedűs, G.; Nagy, Á.; Decsi, K.; Kutasy, B.; Virág, E. Transcriptome datasets of β-Aminobutyric acid (BABA)-primed mono-and dicotyledonous plants, Hordeum vulgare and Arabidopsis thaliana. Data Brief 2022, 41, 107983. [Google Scholar] [CrossRef] [PubMed]
  42. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  43. Beckers, G.J.; Conrath, U. Priming for stress resistance: From the lab to the field. Curr. Opin. Plant Biol. 2007, 10, 425–431. [Google Scholar] [CrossRef] [PubMed]
  44. Hanachi, A.; Bianchi, A.; Kahn, C.J.; Velot, E.; Arab-Tehrany, E.; Cakir-Kiefer, C.; Linder, M. Encapsulation of Salmon Peptides in Marine Liposomes: Physico-Chemical Properties, Antiradical Activities and Biocompatibility Assays. Mar. Drugs 2022, 20, 249. [Google Scholar] [CrossRef]
  45. Atif, M.J.; Amin, B.; Ghani, M.I.; Ali, M.; Liu, X.; Zhang, Y.; Cheng, Z. Allium sativum L.(Garlic) bulb enlargement as influenced by differential combinations of photoperiod and temperature. Food Chem. 2021, 338, 127991. [Google Scholar] [CrossRef] [PubMed]
  46. Travaglia, C.; Reinoso, H.; Cohen, A.; Luna, C.; Tommasino, E.; Castillo, C.; Bottini, R. Exogenous ABA increases yield in field-grown wheat with moderate water restriction. J. Plant Growth Regul. 2010, 29, 366–374. [Google Scholar] [CrossRef]
  47. Aroca, R.; del Mar Alguacil, M.; Vernieri, P.; Ruiz-Lozano, J.M. Plant responses to drought stress and exogenous ABA application are modulated differently by mycorrhization in tomato and an ABA-deficient mutant (sitiens). Microb. Ecol. 2008, 56, 704–719. [Google Scholar] [CrossRef] [PubMed]
  48. Chan, Z. Expression profiling of ABA pathway transcripts indicates crosstalk between abiotic and biotic stress responses in Arabidopsis. Genomics 2012, 100, 110–115. [Google Scholar] [CrossRef] [PubMed]
  49. Huang, Y.; Guo, Y.; Liu, Y.; Zhang, F.; Wang, Z.; Wang, H.; Wang, F.; Li, D.; Mao, D.; Luan, S. 9-cis-epoxycarotenoid dioxygenase 3 regulates plant growth and enhances multi-abiotic stress tolerance in rice. Front. Plant Sci. 2018, 9, 162. [Google Scholar] [CrossRef]
  50. Hao, G.P.; Zhang, X.H.; Wang, Y.Q.; Wu, Z.Y.; Huang, C.L. Nucleotide variation in the NCED3 region of Arabidopsis thaliana and its association study with abscisic acid content under drought stress. J. Integr. Plant Biol. 2009, 51, 175–183. [Google Scholar] [CrossRef]
  51. Zhang, W.; Yang, H.; You, S.; Xu, Y.; Ran, K.; Fan, S. Cloning, characterization and functional analysis of the role MhNCED3, a gene encoding 9-cis-epoxycarotenoid dioxygenase in Malus hupehensis Rehd., plays in plant tolerance to osmotic and Cd2+ stresses. Plant Soil 2014, 381, 143–160. [Google Scholar] [CrossRef]
  52. Iuchi, S.; Kobayashi, M.; Taji, T.; Naramoto, M.; Seki, M.; Kato, T.; Tabata, S.; Kakubari, Y.; Yamaguchi-Shinozaki, K.; Shinozaki, K. Regulation of drought tolerance by gene manipulation of 9-cis-epoxycarotenoid dioxygenase, a key enzyme in abscisic acid biosynthesis in Arabidopsis. Plant J. 2001, 27, 325–333. [Google Scholar] [CrossRef] [Green Version]
  53. Shen, Q.; Chen, C.-N.; Brands, A.; Pan, S.-M.; Tuan-Hua, D.H. The stress-and abscisic acid-induced barley gene HVA22: Developmental regulation and homologues in diverse organisms. Plant Mol. Biol. 2001, 45, 327–340. [Google Scholar] [CrossRef] [PubMed]
  54. Bergler, J.; Hoth, S. Plant U-box armadillo repeat proteins AtPUB18 and AtPUB19 are involved in salt inhibition of germination in Arabidopsis. Plant Biol. 2011, 13, 725–730. [Google Scholar] [CrossRef]
  55. Seo, D.H.; Ahn, M.Y.; Park, K.Y.; Kim, E.Y.; Kim, W.T. The N-terminal UND motif of the Arabidopsis U-box E3 ligase PUB18 is critical for the negative regulation of ABA-mediated stomatal movement and determines its ubiquitination specificity for exocyst subunit Exo70B1. Plant Cell 2016, 28, 2952–2973. [Google Scholar] [CrossRef]
  56. Azevedo, C.; Santos-Rosa, M.J.; Shirasu, K. The U-box protein family in plants. Trends Plant Sci. 2001, 6, 354–358. [Google Scholar] [CrossRef]
  57. Liu, Y.-C.; Wu, Y.-R.; Huang, X.-H.; Sun, J.; Xie, Q. AtPUB19, a U-box E3 ubiquitin ligase, negatively regulates abscisic acid and drought responses in Arabidopsis thaliana. Mol. Plant 2011, 4, 938–946. [Google Scholar] [CrossRef] [PubMed]
  58. Staneloni, R.J.; Rodriguez-Batiller, M.J.; Casal, J.J. Abscisic acid, high-light, and oxidative stress down-regulate a photosynthetic gene via a promoter motif not involved in phytochrome-mediated transcriptional regulation. Mol. Plant 2008, 1, 75–83. [Google Scholar] [CrossRef]
  59. Gao, S.; Gao, J.; Zhu, X.; Song, Y.; Li, Z.; Ren, G.; Zhou, X.; Kuai, B. ABF2, ABF3, and ABF4 promote ABA-mediated chlorophyll degradation and leaf senescence by transcriptional activation of chlorophyll catabolic genes and senescence-associated genes in Arabidopsis. Mol. Plant 2016, 9, 1272–1285. [Google Scholar] [CrossRef] [PubMed]
  60. Vogel, J.T.; Zarka, D.G.; Van Buskirk, H.A.; Fowler, S.G.; Thomashow, M.F. Roles of the CBF2 and ZAT12 transcription factors in configuring the low temperature transcriptome of Arabidopsis. Plant J. 2005, 41, 195–211. [Google Scholar] [CrossRef]
  61. González-Diéguez, D.O.; Orozco, C. Identificación de marcadores moleculares (SNP) y genes candidatos asociados a la tolerancia genética al complejo mancha de asfalto en variedades nativas de maíz (Zea mays L.). Cienc. Tecnol. Salud 2020, 7, 218–235. [Google Scholar] [CrossRef]
  62. Pasini, L.; Bergonti, M.; Fracasso, A.; Marocco, A.; Amaducci, S. Microarray analysis of differentially expressed mRNAs and miRNAs in young leaves of sorghum under dry-down conditions. J. Plant Physiol. 2014, 171, 537–548. [Google Scholar] [CrossRef]
  63. Chaudhary, C.; Sharma, N.; Khurana, P. Decoding the wheat awn transcriptome and overexpressing Ta Rca1β in rice for heat stress tolerance. Plant Mol. Biol. 2021, 105, 133–146. [Google Scholar] [CrossRef]
  64. Mohanty, B. Genomic architecture of promoters and transcriptional regulation of candidate genes in rice involved in tolerance to anaerobic germination. Curr. Plant Biol. 2022, 29, 100236. [Google Scholar] [CrossRef]
  65. Lu, C.-A.; Lin, C.-C.; Lee, K.-W.; Chen, J.-L.; Huang, L.-F.; Ho, S.-L.; Liu, H.-J.; Hsing, Y.-I.; Yu, S.-M. The SnRK1A protein kinase plays a key role in sugar signaling during germination and seedling growth of rice. Plant Cell 2007, 19, 2484–2499. [Google Scholar] [CrossRef]
  66. Sakuma, Y.; Liu, Q.; Dubouzet, J.G.; Abe, H.; Shinozaki, K.; Yamaguchi-Shinozaki, K. DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration-and cold-inducible gene expression. Biochem. Biophys. Res. Commun. 2002, 290, 998–1009. [Google Scholar] [CrossRef]
  67. Moon, S.-J.; Min, M.K.; Kim, J.; Kim, D.Y.; Yoon, I.S.; Kwon, T.R.; Byun, M.O.; Kim, B.-G. Ectopic expression of OsDREB1G, a member of the OsDREB1 subfamily, confers cold stress tolerance in rice. Front. Plant Sci. 2019, 10, 297. [Google Scholar] [CrossRef] [PubMed]
  68. Fujita, M.; Fujita, Y.; Noutoshi, Y.; Takahashi, F.; Narusaka, Y.; Yamaguchi-Shinozaki, K.; Shinozaki, K. Crosstalk between abiotic and biotic stress responses: A current view from the points of convergence in the stress signaling networks. Curr. Opin. Plant Biol. 2006, 9, 436–442. [Google Scholar] [CrossRef] [PubMed]
  69. Medina, J.; Catalá, R.; Salinas, J. The CBFs: Three Arabidopsis transcription factors to cold acclimate. Plant Sci. 2011, 180, 3–11. [Google Scholar] [CrossRef] [PubMed]
  70. Guo, X.; Liu, D.; Chong, K. Cold signaling in plants: Insights into mechanisms and regulation. J. Integr. Plant Biol. 2018, 60, 745–756. [Google Scholar] [CrossRef] [PubMed]
  71. Gilmour, S.J.; Fowler, S.G.; Thomashow, M.F. Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol. Biol. 2004, 54, 767–781. [Google Scholar] [CrossRef] [PubMed]
  72. Xie, M.; Sun, J.; Gong, D.; Kong, Y. The roles of Arabidopsis C1-2i subclass of C2H2-type zinc-finger transcription factors. Genes 2019, 10, 653. [Google Scholar] [CrossRef] [Green Version]
  73. Maruyama, K.; Sakuma, Y.; Kasuga, M.; Ito, Y.; Seki, M.; Goda, H.; Shimada, Y.; Yoshida, S.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Identification of cold-inducible downstream genes of the Arabidopsis DREB1A/CBF3 transcriptional factor using two microarray systems. Plant J. 2004, 38, 982–993. [Google Scholar] [CrossRef] [PubMed]
  74. Lu, S.; Friesen, T.L.; Faris, J.D. Molecular characterization and genomic mapping of the pathogenesis-related protein 1 (PR-1) gene family in hexaploid wheat (Triticum aestivum L.). Mol. Genet. Genom. 2011, 285, 485–503. [Google Scholar] [CrossRef]
  75. Naz, R.; Bano, A.; Wilson, N.L.; Guest, D.; Roberts, T.H. Pathogenesis-related protein expression in the apoplast of wheat leaves protected against leaf rust following application of plant extracts. Phytopathology 2014, 104, 933–944. [Google Scholar] [CrossRef]
  76. Kemp, G.; Botha, A.-M.; Kloppers, F.; Pretorius, Z. Disease development and β-1, 3-glucanase expression following leaf rust infection in resistant and susceptible near-isogenic wheat seedlings. Physiol. Mol. Plant Pathol. 1999, 55, 45–52. [Google Scholar] [CrossRef]
  77. Liu, B.; Xue, X.; Cui, S.; Zhang, X.; Han, Q.; Zhu, L.; Liang, X.; Wang, X.; Huang, L.; Chen, X. Cloning and characterization of a wheat β-1, 3-glucanase gene induced by the stripe rust pathogen Puccinia striiformis f. sp. tritici. Mol. Biol. Rep. 2010, 37, 1045. [Google Scholar] [CrossRef]
  78. Kasprzewska, A. Plant chitinases-regulation and function. Cell. Mol. Biol. Lett. 2003, 8, 809–824. [Google Scholar] [PubMed]
  79. Bertini, L.; Caporale, C.; Testa, M.; Proietti, S.; Caruso, C. Structural basis of the antifungal activity of wheat PR4 proteins. FEBS Lett. 2009, 583, 2865–2871. [Google Scholar] [CrossRef] [PubMed]
  80. Caporale, C.; Di Berardino, I.; Leonardi, L.; Bertini, L.; Cascone, A.; Buonocore, V.; Caruso, C. Wheat pathogenesis-related proteins of class 4 have ribonuclease activity. FEBS Lett. 2004, 575, 71–76. [Google Scholar] [CrossRef]
  81. De Jesús-Pires, C.; Ferreira-Neto, J.R.; Pacifico Bezerra-Neto, J.; Kido, E.A.; de Oliveira Silva, R.L.; Pandolfi, V.; Wanderley-Nogueira, A.C.; Binneck, E.; da Costa, A.F.; Pio-Ribeiro, G. Plant thaumatin-like proteins: Function, evolution and biotechnological applications. Curr. Protein Pept. Sci. 2020, 21, 36–51. [Google Scholar] [CrossRef]
  82. Owens, D.K.; Crosby, K.C.; Runac, J.; Howard, B.A.; Winkel, B.S. Biochemical and genetic characterization of Arabidopsis flavanone 3β-hydroxylase. Plant Physiol. Biochem. 2008, 46, 833–843. [Google Scholar] [CrossRef]
  83. Han, Y.; Huang, K.; Liu, Y.; Jiao, T.; Ma, G.; Qian, Y.; Wang, P.; Dai, X.; Gao, L.; Xia, T. Functional analysis of two flavanone-3-hydroxylase genes from Camellia sinensis: A critical role in flavonoid accumulation. Genes 2017, 8, 300. [Google Scholar] [CrossRef]
  84. Rakwal, R.; Agrawal, G.K.; Jwa, N.-S. Characterization of a rice (Oryza sativa L.) Bowman–Birk proteinase inhibitor: Tightly light regulated induction in response to cut, jasmonic acid, ethylene and protein phosphatase 2A inhibitors. Gene 2001, 263, 189–198. [Google Scholar] [CrossRef]
  85. Wu, S.-W.; Wang, H.-W.; Yang, Z.-D.; Kong, L.-R. Expression comparisons of pathogenesis-related (PR) genes in wheat in response to infection/infestation by Fusarium, Yellow dwarf virus (YDV) aphid-transmitted and hessian fly. J. Integr. Agric. 2014, 13, 926–936. [Google Scholar] [CrossRef]
  86. Pasquer, F. Effects of Plant Protection Compounds on Wheat Gene Expression. Ph.D. Thesis, University of Zurich, Zürich, Switzerland, 2005. [Google Scholar]
  87. Carbonero, P.; Díaz, I.; Vicente-Carbajosa, J.; Alfonso-Rubi, J.; Gaddour, K.; Lara, P. Cereal α-amylase/trypsin inhibitors and transgenic insect resistance. In Genetics and Breeding for Crop Quality and Resistance; Springer: Ordrecht, The Netherlands, 1999; pp. 147–158. [Google Scholar] [CrossRef]
  88. Erdayani, E.; Nagarajan, R.; Grant, N.P.; Gill, K.S. Genome-wide analysis of the HSP101/CLPB gene family for heat tolerance in hexaploid wheat. Sci. Rep. 2020, 10, 3948. [Google Scholar] [CrossRef] [PubMed]
  89. Biselli, C.; Bagnaresi, P.; Faccioli, P.; Hu, X.; Balcerzak, M.; Mattera, M.G.; Yan, Z.; Ouellet, T.; Cattivelli, L.; Valè, G. Comparative transcriptome profiles of near-isogenic hexaploid wheat lines differing for effective alleles at the 2DL FHB resistance QTL. Front. Plant Sci. 2018, 9, 37. [Google Scholar] [CrossRef] [PubMed]
  90. Douchkov, D.; Johrde, A.; Nowara, D.; Himmelbach, A.; Lueck, S.; Niks, R.; Schweizer, P. Convergent evidence for a role of WIR1 proteins during the interaction of barley with the powdery mildew fungus Blumeria graminis. J. Plant Physiol. 2011, 168, 20–29. [Google Scholar] [CrossRef]
  91. Zhao, J.; Zhou, H.; Li, X. Ubiquitin-specific protease16 interacts with a heavy metal associated isoprenylated plant protein27 and modulates cadmium tolerance. Plant Signal. Behav. 2013, 8, e25680. [Google Scholar] [CrossRef]
  92. De Abreu-Neto, J.B.; Turchetto-Zolet, A.C.; de Oliveira, L.F.V.; Bodanese Zanettini, M.H.; Margis-Pinheiro, M. Heavy metal-associated isoprenylated plant protein (HIPP): Characterization of a family of proteins exclusive to plants. FEBS J. 2013, 280, 1604–1616. [Google Scholar] [CrossRef]
  93. Dooner, H.K.; Robbins, T.P.; Jorgensen, R.A. Genetic and developmental control of anthocyanin biosynthesis. Annu. Rev. Genet. 1991, 25, 173–199. [Google Scholar] [CrossRef]
  94. Tan, D.-X.; Hardeland, R.; Manchester, L.C.; Korkmaz, A.; Ma, S.; Rosales-Corral, S.; Reiter, R.J. Functional roles of melatonin in plants, and perspectives in nutritional and agricultural science. J. Exp. Bot. 2012, 63, 577–597. [Google Scholar] [CrossRef]
  95. Cui, G.; Zhao, X.; Liu, S.; Sun, F.; Zhang, C.; Xi, Y. Beneficial effects of melatonin in overcoming drought stress in wheat seedlings. Plant Physiol. Biochem. 2017, 118, 138–149. [Google Scholar] [CrossRef]
  96. Jain, D.; Khurana, J.P. Role of pathogenesis-related (PR) proteins in plant defense mechanism. In Molecular Aspects of Plant-Pathogen Interaction; Springer: Singapore, 2018; pp. 265–281. [Google Scholar] [CrossRef]
  97. Gullner, G.; Komives, T.; Király, L.; Schröder, P. Glutathione S-transferase enzymes in plant-pathogen interactions. Front. Plant Sci. 2018, 9, 1836. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Sample collection at three timepoints of the nanoscale liposome-encapsulated SC-CO2 garlic extract (‘Garlic-lipo’) on Triticum aestivum in greenhouse experiments. The effects were examined in a large-scale genome-wide transcriptional profiling experiment collecting samples at various timepoints (0 min is the control, and 15 min, 24 h, and 48 h after spaying) and were evaluated by using DEG and time-course expression analysis.
Figure 1. Sample collection at three timepoints of the nanoscale liposome-encapsulated SC-CO2 garlic extract (‘Garlic-lipo’) on Triticum aestivum in greenhouse experiments. The effects were examined in a large-scale genome-wide transcriptional profiling experiment collecting samples at various timepoints (0 min is the control, and 15 min, 24 h, and 48 h after spaying) and were evaluated by using DEG and time-course expression analysis.
Plants 12 00743 g001
Figure 2. Size determination diagram of liposomes from the ’Garlic-lipo’ colloid. Liposome sizes are in the range between 134.2 and 164.2 nm.
Figure 2. Size determination diagram of liposomes from the ’Garlic-lipo’ colloid. Liposome sizes are in the range between 134.2 and 164.2 nm.
Plants 12 00743 g002
Figure 3. Determination of the Zeta potential based on the electrophoretic mobility of liposomes. (a) Zeta potential distribution: the colloid stability is proven to be high, showing the Zeta potential between −40.66 and −40.53 mV. (b) Phase plot to estimate the electrophoretic mobility of particles: Phase difference is proportional to the electrophoretic mobility of the nanomaterial and, thus to the Zeta potential of the aqueous medium.
Figure 3. Determination of the Zeta potential based on the electrophoretic mobility of liposomes. (a) Zeta potential distribution: the colloid stability is proven to be high, showing the Zeta potential between −40.66 and −40.53 mV. (b) Phase plot to estimate the electrophoretic mobility of particles: Phase difference is proportional to the electrophoretic mobility of the nanomaterial and, thus to the Zeta potential of the aqueous medium.
Plants 12 00743 g003
Figure 4. Heat map of the top 50 differentially expressed genes (DEGs) based on the comparison of the 48 h vs. control samples. The heatmap color range is from red for positive Z-score values to green for Z-score negative values. The number of differentially expressed (DE) features (Probability > 0.9): 2578. Up-regulated (M > 0): 865. Down-regulated (M < 0): 1713. For the annotation of the transcript IDs, see the Annotation Table (DOI: 10.17632/p66v4yxbtp.1).
Figure 4. Heat map of the top 50 differentially expressed genes (DEGs) based on the comparison of the 48 h vs. control samples. The heatmap color range is from red for positive Z-score values to green for Z-score negative values. The number of differentially expressed (DE) features (Probability > 0.9): 2578. Up-regulated (M > 0): 865. Down-regulated (M < 0): 1713. For the annotation of the transcript IDs, see the Annotation Table (DOI: 10.17632/p66v4yxbtp.1).
Plants 12 00743 g004
Figure 5. The determination of the top 50 differentially expressed genes (DEGs) in silico is represented in three sampling times. The numerical analysis of DEGs in a pairwise comparison of two times was examined using OmixBox.BioBam based on the RSEM and edgeR programs. These quantitative statistical methods to evaluate the significance of individual genes were implemented to examine the changes in the expression levels of (a) PR genes (control vs. 48 h); (b) ABA pathway-related genes (control vs. 15 min); and (c) TF genes (control vs. 15 min).
Figure 5. The determination of the top 50 differentially expressed genes (DEGs) in silico is represented in three sampling times. The numerical analysis of DEGs in a pairwise comparison of two times was examined using OmixBox.BioBam based on the RSEM and edgeR programs. These quantitative statistical methods to evaluate the significance of individual genes were implemented to examine the changes in the expression levels of (a) PR genes (control vs. 48 h); (b) ABA pathway-related genes (control vs. 15 min); and (c) TF genes (control vs. 15 min).
Plants 12 00743 g005
Figure 6. KEGG analysis of cluster 1 at the three investigated times, relative to the control (0 min). The expression levels of these genes gradually increase. The graph shows the numbers of transcripts that are annotated (DOI:10.17632/p66v4yxbtp.1) and can be placed into the KEGG pathway (Tables S1 and S2) using the TimeCourseTable (DOI: 10.17632/xvvscxpz6w.1).
Figure 6. KEGG analysis of cluster 1 at the three investigated times, relative to the control (0 min). The expression levels of these genes gradually increase. The graph shows the numbers of transcripts that are annotated (DOI:10.17632/p66v4yxbtp.1) and can be placed into the KEGG pathway (Tables S1 and S2) using the TimeCourseTable (DOI: 10.17632/xvvscxpz6w.1).
Plants 12 00743 g006
Figure 7. Results of the real-time PCR reactions performed at three sampling timepoints compared to control. The changes in the expression levels of (a) PR genes; (b) ABA pathway-related genes; and (c) TF genes. The values are represented as mean ± SEM (n = 4). p < 0.05 vs. Control, p < 0.05 vs. treated group. One asterisk indicates statistically significant difference between the means of a treated sample set compared to the mean of the control sample set to 5%; two asterisks indicate statistically significant difference to 1%; three asterisks indicate statistically significant difference to 0.1%.
Figure 7. Results of the real-time PCR reactions performed at three sampling timepoints compared to control. The changes in the expression levels of (a) PR genes; (b) ABA pathway-related genes; and (c) TF genes. The values are represented as mean ± SEM (n = 4). p < 0.05 vs. Control, p < 0.05 vs. treated group. One asterisk indicates statistically significant difference between the means of a treated sample set compared to the mean of the control sample set to 5%; two asterisks indicate statistically significant difference to 1%; three asterisks indicate statistically significant difference to 0.1%.
Plants 12 00743 g007
Figure 8. Summary of up-regulated DEGs (yellow) found in the KEGG pathways (not colored), with the patway-related secondary metabolites (blue), and their physiological roles (red) in defense mechanism and plant development as a response to ‘Garlic-lipo’ treatment.
Figure 8. Summary of up-regulated DEGs (yellow) found in the KEGG pathways (not colored), with the patway-related secondary metabolites (blue), and their physiological roles (red) in defense mechanism and plant development as a response to ‘Garlic-lipo’ treatment.
Plants 12 00743 g008
Table 1. Particle size distribution (PSD) by intensity statistics.
Table 1. Particle size distribution (PSD) by intensity statistics.
NameMeanMinimumMaximum
Z-Average (nm)149.2134.2164.2
Polydispersity Index (PI)0.130350.13030.1304
Intercept0.9510.9490.953
Peak One Mean by Intensity (nm)160.6129.6191.6
Peak One Area by Intensity (%)100100100
Derived Mean Count Rate (kcps)171,3001.71 × 1051.71 × 105
Table 2. Measurement of Zeta potential and quality of individual liposomes of the ’Garlic-lipo’ samples.
Table 2. Measurement of Zeta potential and quality of individual liposomes of the ’Garlic-lipo’ samples.
NameMeanMinimumMaximum
Derived Mean Count Rate (kcps)171,3001.71 × 1051.71 × 105
Zeta Potential (mV)−40.595−40.66−40.53
Zeta Peak One Mean−41.18−41.43−40.93
Zeta Peak Two Mean−9.2−9.25−9.15
Conductivity (mS/cm)0.0074390.0073990.007479
Wall Zeta Potential (mV)−66.28−66.33−66.23
Zeta Deviation (mV)8.6328.6078.657
Derived Mean Count Rate (kcps)85,1508.52 × 1048.52 × 104
Reference Beam Count Rate (kcps)3671.536453698
Quality Factor4.5094.5084.51
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kutasy, B.; Kiniczky, M.; Decsi, K.; Kálmán, N.; Hegedűs, G.; Alföldi, Z.P.; Virág, E. ‘Garlic-lipo’4Plants: Liposome-Encapsulated Garlic Extract Stimulates ABA Pathway and PR Genes in Wheat (Triticum aestivum). Plants 2023, 12, 743. https://doi.org/10.3390/plants12040743

AMA Style

Kutasy B, Kiniczky M, Decsi K, Kálmán N, Hegedűs G, Alföldi ZP, Virág E. ‘Garlic-lipo’4Plants: Liposome-Encapsulated Garlic Extract Stimulates ABA Pathway and PR Genes in Wheat (Triticum aestivum). Plants. 2023; 12(4):743. https://doi.org/10.3390/plants12040743

Chicago/Turabian Style

Kutasy, Barbara, Márta Kiniczky, Kincső Decsi, Nikoletta Kálmán, Géza Hegedűs, Zoltán Péter Alföldi, and Eszter Virág. 2023. "‘Garlic-lipo’4Plants: Liposome-Encapsulated Garlic Extract Stimulates ABA Pathway and PR Genes in Wheat (Triticum aestivum)" Plants 12, no. 4: 743. https://doi.org/10.3390/plants12040743

APA Style

Kutasy, B., Kiniczky, M., Decsi, K., Kálmán, N., Hegedűs, G., Alföldi, Z. P., & Virág, E. (2023). ‘Garlic-lipo’4Plants: Liposome-Encapsulated Garlic Extract Stimulates ABA Pathway and PR Genes in Wheat (Triticum aestivum). Plants, 12(4), 743. https://doi.org/10.3390/plants12040743

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop