- Research
- Open access
- Published:
Pungency related gene network in Allium sativum L., response to sulfur treatments
BMC Genomic Data volume 25, Article number: 35 (2024)
Abstract
Pungency of garlic (Allium sativum L.) is generated from breakdown of the alk(en)yl cysteine sulphoxide (CSO), alliin and its subsequent breakdown to allicin under the activity of alliinase (All). Based on recent evidence, two other important genes including Sulfite reductase (SiR) and Superoxide dismutase (SOD) are thought to be related to sulfur metabolism. These three gene functions are in sulfate assimilation pathway. However, whether it is involved in stress response in crops is largely unknown. In this research, the order and priority of simultaneous expression of three genes including All, SiR and SOD were measured on some garlic ecotypes of Iran, collected from Zanjan, Hamedan and Gilan, provinces under sulfur concentrations (0, 6, 12, 24 and 60 g/ per experimental unit: pot) using real-time quantitative PCR (RT-qPCR) analysis. For understanding the network interactions between studied genes and other related genes, in silico gene network analysis was constructed to investigate various mechanisms underlying stimulation of A. sativum L. to cope with imposed sulfur. Complicated network including TF-TF, miRNA-TF, and miRNA-TF-gene, was split into sub-networks to have a deeper insight. Analysis of q-RT-PCR data revealed the highest expression in All and SiR genes respectively. To distinguish and select significant pathways in sulfur metabolism, RESNET Plant database of Pathway Studio software v.10 (Elsevier), and other relative data such as chemical reactions, TFs, miRNAs, enzymes, and small molecules were extracted. Complex sub-network exhibited plenty of routes between stress response and sulfate assimilation pathway. Even though Alliinase did not display any connectivity with other stress response genes, it showed binding relation with lectin functional class, as a result of which connected to leucine zipper, exocellulase, peroxidase and ARF functional class indirectly. Integration network of these genes revealed their involvement in various biological processes such as, RNA splicing, stress response, gene silencing by miRNAs, and epigenetic. The findings of this research can be used to extend further research on the garlic metabolic engineering, garlic stress related genes, and also reducing or enhancing the activity of the responsible genes for garlic pungency for health benefits and industry demands.
Introduction
In biological ecosystems, along with numerous growth-limiting factors, abiotic and biotic stresses are important challenging physiological barriers that have major effects on plant growth and development. Faced with adverse environmental conditions, plants reprogram their cellular activities through several genes and minimize stress damage using regulatory mechanisms, including post-transcriptional regulation of gene expression [1]. Transcription Factors (TFs) and non-coding RNAs are the important regulatory elements in functional genomics [2]. All biological processes in the organisms are managed primarily by changes in the activity and expression of key genes. The ability of a cell to switch on and off a gene drives all biological function and activity. Because of key effect of genes on metabolism, researchers have focused much interest on gene expression profiling to identify those key genes and gene clusters whose expression changes and variation in different biological states [3,4,5,6]. Research studies have shown that changes in gene expression are the key basis of response to environmental factors, stresses, defense response against pathogens, regulation of metabolic pathways, regulation of photosynthesis or symbiosis, so plant by using this possibility and resourcefulness, regulates and manages its biological activities for survival and reproduction. At the level of transcription, specific transcription factors (TFs) bind DNA in order to activate or repress the expression of a gene. MiRNAs repress gene expression post-transcriptionally by interacting with complementary sequences located in the 3′UTR of their target mRNAs [7]. For example, co-expressed genes and genes encoding interacting proteins tend to be regulated by common TFs [8]. Although synthetic genetic interactions mostly occur between homologous genes, large gene families have been identified that complicate interactions between some important genes [9]. Genes encoding TFs that control miRNA expression have a higher chance to be post-transcriptionally repressed by the miRNA. Furthermore, genes co-regulated by miRNAs are less functionally linked than genes co-regulated by TFs. Therefore, different types of molecular interactions provide additional insights and new horizons of valuable information about gene regulation and cell function, which indicates the need and comprehensive scientific attention to integrating data and extracting existing connections and correlations [10, 11].
Gene network analysis (GNA) use computer and bioinformatic data tools to model and simulate genetic regulatory networks. The important aim of GNA is to help life science researchers build a model of a genetic regulatory network using knowledge and information related to regulatory interactions in combination with gene expression [6, 7, 9].
Garlic (Allium sativum L.), a diploid (2n = 2x = 16) plant species is an important economic and medicinal plant with more than 5,000 years of planting history on the planet. Pungency of garlic is related to organosulfur components including allicin. Allicin is considered responsible for most of the pharmacological activity of crushed raw garlic cloves [4]. When the tissues of any Allium species are disrupted, these amino acid derivatives are cleaved by the enzyme alliinase (EC 4.4.1.4) into their corresponding sulfenic acids, and volatile sulfur compounds are produced that give the characteristic flavor and bioactivity of the species [5, 6]. Despite its agronomic importance, garlic remains largely a wildl plant, hampering its commercial potential. Garlic is the first species within the Allium genus to be sequenced [11, 12]. According to these sequencing results, garlic experienced a recent occurrence of burst of transposable elements. Alliinase genes and content are thought to have rapidly expanded during the burst of transposable elements which helps explain the evolution of allicin biosynthesis-related genes [13]. Two other genes including sulfite reductase (SiR) and superoxide dismutase (SOD) are thought related to sulfate assimilation pathway. However, whether it is involved in stress response in crops is largely unknown. These genes are therefore potential candidates of the alliin biosynthesis pathway.
Pungent and weak odor in garlic has its own advantages depending on the purpose of garlic consumption. Certainly, food processing industries as well as spice uses require the maximum aroma and spiciness of garlic, if it is consumed fresh and raw, a little spiciness of garlic can be desirable. Genetic manipulation or environmental and agricultural interventions to achieve the various goals above will require detailed knowledge on the network of related genes involved in the biosynthesis of spiciness in garlic and onion.
Materials and methods
Plant materials
The academic permission for collect and research on medicinal plants was obtained from Head of Biotechnology, Department Research Institute of Modern Biological Techniques, University of Zanjan, Zanjan, Iran. The study complies with all relevant guidelines. Three different garlic clones (with high, medium and low odor) (12) that had been selected based on pungency spectrophotometric method [12] were cultured in pots (Fig. 1) containing different concentrations of sulfur (0, 6, 12, 18, 24, 60 g.). Thiobacillus bacteria were added to each pot containing 1 kg of soil. The pots were given a month to undergo for oxidation process of sulfur and the dissolved sulfur in the soil. After one month, 3 mentioned clones were planted in three replicates, in depths of 5 cm. After three months of cultivation and at the end of the growing season, the plants were harvested and the garlic cloves harvested were used for RNA extraction and gene expression studies.
Gene expression analysis
RNAs were extracted based on Ribo Spin Plant Kit, which prepared of GeneAll Oguem-dong, Songpa-gu, Seoul, Korea. The cDNAs were synthesized using Thermo kit (GeneAll Oguem-dong, Songpa-gu, Seoul, Korea) accordingly [14]. Expression of three pungency related genes including All, SiR and SOD (Table 1) were evaluated under sulfur concentrations (0, 6, 12, 18, 24 and 60 g) using q-RT-PCR. The characteristics of studied primers were shown in Table 1 [12].
Gene network analysis
To distinguish and select significant pathways in sulfur metabolism, RESNET Plant database of Pathway Studio software v.10 (Elsevier), and other relative data such as chemical reactions, TFs, miRNAs, enzymes, and small molecules were extracted. This database includes aliases for genes in the model plants and the other plants including tobacco, and tomato. To predict interaction between genes and sulfur molecule, various statistical tests such as, Fisher’s Exact Test were used [15, 16]. To make statistical network based on sulfur and candidate genes, union selected, physical and direct interaction algorithms were used.
Promoter analysis
Promoter analysis can provide valuable information about underlying regulatory mechanism and function of genes in response to lots of signals [17, 18]. The sequence of SIR (AT5G04590) of Arabidopsis (as a model plant) was downloaded from the NCBI database (www.ncbi.nlm.nih.gov). 1500 Kb upstream (from the start codon) was extracted as promoter sequences from Ensembl through the BioMart tool (http://plants.ensembl.org/biomart). To find the transcription factors and their binding sites across the promoter region Plantpan v.2 database was used (http://plantpan2.itps.ncku.edu.tw/promoter.php) [19].
Result and discussion
Some of used garlic clones for this research were present in Fig. 2. The sequence and priority of the simultaneous expression of the three candidate genes compared to each other is shown in Fig. 3. Gene expression analysis revealed the highest expression in All and SiR genes respectively.
Allinase gene is the first that was observed in this study in the first place in the amount and speed of expression. The next gene is the SiR and the third gene in terms of the speed of gene expression is SOD. This graph specifically shows the order and priority of expression of the three mentioned genes at the same time. Because the focus of the research is on the study of the in silico gene network involved in the absorption and transport of sulfur in the garlic plant, therefore, we have not provided the molecular data of gene expression changes here, and we only wanted to provide the necessary primary documentation give to the readers including experimental and field works (Fig. 1), providing the primers (Table 1) and the expression order of the three studied genes (Fig. 3). However, the changes in the expression of the three studied genes are shown in Fig. 4.
A gene network was constructed to investigate various mechanisms stimulation of A. sativum L. to cope with imposed stress (Fig. 5). To closer inspect, complicated network, including TF-TF, miRNA-TF, and miRNA-TF-gene, was split into sub-network. Figure 6, complicated sub- network, exhibited plenty of routes between stress response and sulfate assimilation pathway.
GO analysis is a strong approach in understanding the molecular mechanisms underpinning developmental and environmental processes and offer a reliable tool for GO gene selection [20]. Integration network revealed the genes are involved in various biological processes such as, RNA splicing, stress response, gene silencing by miRNAs, and epigenetic.
Ethylene-insensitive3-like 3, EIL3 (SLIM1) key regulator in sulfur assimilation and garlic pungency
EIL3 is considered as central hub in sulfur response and metabolism [21]. Network analysis showed EIL3 controls lots of genes involving not only in sulfur assimilation such as, miR395, SULTR1;2, and APS4, but also glucosinolate biosynthesis process such as MYB34 [22, 23]. Glucosinolates, secondary metabolites, has roles in plant defense and inducers of anticarcinogenic in human [24]. In addition, Allicin has a huge number of activities such as, immune response, inhibiting the proliferation of tumor cells and induced apoptosis in gastric epithelial, breast cancer cells (Fig. 5). Since abundance of Allinase was high under normal condition, and plant needs to assimilate sulfur, it could be concluded that the sulfur uptake and assimilation, the intensity of pungency and drug metabolic in garlic might be controlled by EIL3.
Alliinase interplay between oxidase and metabolic pathways
Even though Alliinase did not display any connectivity with other stress response genes, it showed binding relation with lectin functional class, as a result of which connected to leucine zipper, exocellulase, peroxidase and ARF functional class indirectly (Fig. 5). ARFs (Auxin Response Factors) families involved in auxin signaling are regulated by miR167 and miR160 under abiotic stress [15, 25] Maruyama-Nakashita [20] reported expression of SULTR1;1 was induced by ARFs and Aux/IAA proteins under sulfur deficiency.
Moreover, Because of having two domains, EGF_alliinase (PF04863) and Alliinase_C (PF04864), it is presumed that Alliinase might be a part of primitive plant defense response. Therefore, this gene is speculated to be a significant linker between secondary metabolic and abiotic stress tolerance pathways.
Nonexpresser of PR genes 1 (NPR1) prior regulator in response to sulfur stress
Network analysis of SIR showed interaction with histone functional class, and post-transcriptional gene silencing [26, 27], and its expression is controlled by NPR1, central key in response to salt and oxidative stress tolerance in Arabidopsis [28].
Sub-network showed the expression of NPR1 is affected by glutathione, redox signaling molecule in defense response [29, 30] and interacts with leucine zipper and bZIP transcription factor functional classes. Moreover, CSDs are regulated by miR398, DCL1, involved in various biological processes such as RNA interferences, gene silencing by miRNAs and production of miRNAs, and LSD, is regulated by NPR1; as a consequence, NPR1 regulate them indirectly (Fig. 6).
Since SIR induced earlier than CSDs during sulfur supply conditions, NPR1 might be considered as a pivotal factor in response to stress. In addition, promoter analysis of SIR revealed that most of transcription factors belong to bHLH and bZIP family such as REV, NF-YA10, NF-YB1, and NF-YA2, which produce the high number of alternative splicing variants in barely. Alternative splicing might act as regulatory link between miRNAs and stress response [31]. Therefore, SIR might be considered as a key mediator gene signaling pathway in response to stress in garlic.
Conclusion
For the modern biologist, there are numerous computational strategies that can be employed to assay gene expression. Many of these are based on utilizing collections of expressed sequence tags (ESTs), unique segments of cDNA with base sequences identical to at least part of the coding region of a gene [20, 32]. Gene expression network reconstruction and analysis are starting to be widely used to characterize and predict biosystem behavior, giving rise to a new branch of biological knowledge, ‘network genomics’ [33]. Until recently, such analyses have been limited to one level of manifestation of the genetic information, i.e. transcript networks [34] or metabolic networks (https://doi.org/10.1371/journal.pone.0058759). However, changes in transcript levels are transferred to changes in metabolite levels and thereby to physiological endpoints via adaptations of physiology and homeostasis [35]. Complicated sub- network, exhibited plenty of routes between stress response and sulfate assimilation pathway. Even though Alliinase (alliin lyase) did not display any connectivity with other stress response genes, it showed binding relation with lectin functional class, as a result of which connected to leucine zipper, exocellulase, peroxidase and ARF functional class indirectly. Integration network of these genes revealed they are involved in various biological processes such as, RNA splicing, stress response, and gene silencing by miRNAs. In this research, the order and priority of simultaneous expression of three important genes including All, SiR and SOD were measured on some garlic ecotypes of Iran. Also, based on total information resulted of simultaneous amplification curve of ALL, SiR and SOD genes, in silico gene network related to these three key genes were presented. The findings of this research can be used in further research on the garlic metabolic engineering, garlic stress related genes, and also reducing or enhancing the activity of the responsible genes for garlic pungency and health beneficial. This pungency related network analyzing that is first report on garlic pungency, let us total view for distribution and function of related genes.
Availability of data and materials
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.
References
Zhang Y, Xu J, Li R, Ge Y, Li Y, Li R. Plants’ Response to Abiotic Stress: Mechanisms and Strategies. Int J Mol Sci. 2023;24(13):10915. https://doi.org/10.3390/ijms241310915. PMID:37446089;PMCID:PMC10341657.
Ebrahimi Khaksefidi R, Mirlohi S, Khalaji F, Fakhari Z, Shiran B, Fallahi H, Rafiei F, Budak H, Ebrahimie E. Differential expression of seven conserved microRNAs in response to abiotic stress and their regulatory network in Helianthus annuus. Front Plant Sci. 2015;6:741. https://doi.org/10.3389/fpls.2015.00741.
Murray D, Doran P, MacMathuna P, et al. In silico gene expression analysis – an overview. Mol Cancer. 2007;6:50. https://doi.org/10.1186/1476-4598-6-50.
Lawson LD, Hunsaker SM. Allicin bioavailability and bioequivalence from garlic supplements and garlic foods. Nutrients. 2018;10(7):812. https://doi.org/10.3390/nu10070812. PMID: 29937536; PMCID: PMC6073756.
Eady CC, Kamoi T, Kato M, Porter NG, Davis S, Shaw M, Kamoi A, Amai S. Silencing onion lachrymatory factor synthase causes a significant change in the sulfur secondary metabolite profile. Am Soc Plant Biol. 2005;147:2096–106.
Yang X, Su Y, Wu J, et al. Parallel analysis of global garlic gene expression and alliin content following leaf wounding. BMC Plant Biol. 2021;21:174. https://doi.org/10.1186/s12870-021-02948-0.
Defoort J, Van de Peer Y, Vermeirssen V. Function, dynamics and evolution of network motif modules in integrated gene regulatory networks of worm and plant. Nucleic Acids Res. 2018;46(13):6480–503.
Manke T, Bringas R, Vingron M. Correlating protein-DNA and protein-protein interaction networks. J Mol Biol. 2003;333:75–85.
Costanzo M, Baryshnikova A, Bellay J, Kim Y, Spear ED, Sevier CS, Ding H, Koh JL, Toufighi K, Mostafavi S. The genetic landscape of a cell. Science. 2010;327:425–31.
Mitra K, Carvunis AR, Ramesh SK, Ideker T. Integrative approaches for finding modular structure in biological networks. Nat Rev Genet. 2013;14:719–32.
Yang X, Su Y, Wu J, et al. Parallel analysis of global garlic gene expression and alliin content following leaf wounding. BMC Plant Biol. 2021;21:174. https://doi.org/10.1186/s12870-021-02948-0.
Ammarellou A, Kazemitabar K, Najafi H. Study on morpho-physiological and molecular diversity and pungency related genes expression in garlic. Sari: Ph. D. Thesis, Sari Agricultural Sciences and Natural Resources University; 2015. p. 200.
Ketter CAT, Randle WM. Pungency assessment in onions. In: Karcher SJ, editor. Tested studies for laboratory teaching. New York: Chapter 11. Association for Biology Laboratory Education (ABLE); 1998. p. 177–96.
Wang G, Tian C, Wang Y, Wan F, Hu L, Xiong A, Tian J. Selection of reliable reference genes for quantitative RT-PCR in garlic under salt stress. PeerJ. 2019;7:e7319. https://doi.org/10.7717/peerj.7319.
Nikitin A, Egorov S, Daraselia N, Mazo I. Pathway studio—the analysis and navigation of molecular networks. Bioinformatics. 2003;19(16):2155–7.
Ebrahimie M, Esmaeili F, Cheraghi S, Houshmand F, Shabani L, Ebrahimie E. Efficient and simple production of insulin-producing cells from embryonal carcinoma stem cells using mouse neonate pancreas extract, as a natural inducer. PLoS One. 2014;9(3):e90885.
Deihimi T, Niazi A, Ebrahimi M, Kajbaf K, Fanaee S, Bakhtiarizadeh MR, Ebrahimie E. Finding the undiscovered roles of genes: an approach using mutual ranking of coexpressed genes and promoter architecture-case study: dual roles of thaumatin like proteins in biotic and abiotic stresses. Springerplus. 2012;1(1):30.
Babgohari MZ, Ebrahimie E, Niazi A. In silico analysis of high affinity potassium transporter (HKT) isoforms in different plants. Aquatic Biosystems. 2014;10:9–23.
Chow CN, Zheng HQ, Wu NY, Chien CH, Huang HD, Lee TY, . . . Chang WC. PlantPAN 2.0: an update of plant promoter analysis navigator for reconstructing transcriptional regulatory networks in plants. Nucleic Acids Res. 2016;44(D1):D1154-D1160.
Fruzangohar M, Ebrahimie E, Ogunniyi AD, Mahdi LK, Paton JC, et al. Comparative GO: A Web Application for Comparative Gene Ontology and Gene Ontology-Based Gene Selection in Bacteria. PLoS One. 2013;8(3):e58759.
Maruyama-Nakashita A, Nakamura Y, Tohge T, Saito K, Takahashi H. Arabidopsis SLIM1 is a central transcriptional regulator of plant sulfur response and metabolism. Plant Cell. 2006;18(11):3235–51.
Kawashima CG, Matthewman CA, Huang S, Lee BR, Yoshimoto N, Koprivova A, et al. Interplay of SLIM1 and miR395 in the regulation of sulfate assimilation in Arabidopsis. Plant J. 2011;66(5):863–76.
Wawrzyńska A, Sirko A. To control and to be controlled: understanding the Arabidopsis SLIM1 function in sulfur deficiency through comprehensive investigation of the EIL protein family. Front Plant Sci. 2014;5.
Gigolashvili T, Engqvist M, Yatusevich R, Müller C, Flügge UI. HAG2/MYB76 and HAG3/MYB29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. New Phytol. 2008;177(3):627–42.
Khraiwesh B, Zhu JK, Zhu J. Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants. Biochim Biophys Acta. 2012;1819(2):137–48.
Pipal A, Goralik-Schramel M, Lusser A, Lanzanova C, Sarg B, Loidl A, Lindner H, Rossi V, Loidl P. Regulation and processing of maize histone deacetylase Hda1 by limited proteolysis. Plant Cell. 2003;15:1904–17.
Takeda S, Tadele Z, Hofmann I, Probst AV, Angelis KJ, Kaya H, . . . Shibahara KI. BRU1, a novel link between responses to DNA damage and epigenetic gene silencing in Arabidopsis. Gen Dev 2004;18(7):782–793.
Jayakannan M, Bose J, Babourina O, Rengel Z, Shabala S. Salicylic acid in plant salinity stress signalling and tolerance. Plant Growth Regul. 2015. https://doi.org/10.1007/s10725-015-0028-z.
Nafisi M, Goregaoker S, Botanga CJ, Glawischnig E, Olsen CE, Halkier BA, Glazebrook J. Arabidopsis cytochrome P450 monooxygenase 71A13 catalyzes the conversion of indole-3-acetaldoxime in camalexin synthesis. Plant Cell. 2007;19:2039–52.
Kovács G, Kalmár T, Endreffy E, Ondrik Z, Iványi B, Rikker C, et al. Efficient Targeted Next Generation Sequencing-Based Workflow for Differential Diagnosis of Alport-Related Disorders. PLoS One. 2016;11(3):e0149241. https://doi.org/10.1371/journal.pone.0149241.
Panahi Y, Hosseini MS, Khalili N, Naimi E, Majeed M, Sahebkar A. Antioxidant and anti-inflammatory effects of curcuminoid-piperine combination in subjects with metabolic syndrome: A randomized controlled trial and an updated meta-analysis. Clin Nutr. 2015;34(6):1101–8. https://doi.org/10.1016/j.clnu.2014.12.019.
Forst CV. Network genomics: a novel approach for the analysis of biological systems in the post-genomic era. Molecular Biology Reporter. 2002;29:265–80.
Featherstone DE, Broadie K. Wrestling with pleiotropy: genomic and topological analysis of the yeast gene expression network. BioEssays. 2002;24:267–74.
Fiehn O, Weckwerth W. Deciphering metabolic networks. Eur J Biochem. 2003;270:579–88.
Nikiforova V. Gakie`reB, KempaS, AdamikM, WillmitzerL, HesseH, Hoefgen R Towards dissecting nutrient metabolism in plants:a systems biology case study on sulfur metabolism. J Exp Bot. 2004;55:1861–70.
Acknowledgements
The author is grateful for the spiritual and scientific support of Research Institute of Modern Biological Techniques, University of Zanjan. Zanjan, Iran.
Funding
Some experimental and study requirements such as collection site, analysis, and interpretation of data were supported by Research Institute of Modern Biological Techniques, University of Zanjan, Iran.
Author information
Authors and Affiliations
Contributions
A.A. carried out collection field and pot experiments, measured data, and performed data analysis. designed the experiment and then revised and edited the article. Author read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The use of plant parts in the present study complies with international, national and/or institutional guidelines. All methods were performed in accordance with the relevant guidelines/regulations/legislation.
Consent for publication
Not applicable.
Competing interests
Not applicable.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ammarellou, A. Pungency related gene network in Allium sativum L., response to sulfur treatments. BMC Genom Data 25, 35 (2024). https://doi.org/10.1186/s12863-024-01206-0
Published:
DOI: https://doi.org/10.1186/s12863-024-01206-0