Research & Reviews: Journal of Botanical Sciences
MenuISSN: 2320-0189
ISSN: 2320-0189
+1-845-458-6882Department of Biochemistry and Molecular Biology, Shahjalal University of Science and Technology, Sylhet-3114, Bangladesh
Received date:29/05/2016;Accepted date:02/07/2016;Published date:04/07/2016
Visit for more related articles atraybet36
Detoxificationof a highly cytotoxic metabolite methylglyoxal (MG) is carried out through the sequential action of two glyoxalase enzymes, glyoxalase I (GLYI) and glyoxalase II (GLYII). Physiological and metabolic role of glyoxalase proteins have been comprehensively investigated in various plant and animal species. Previously genome-wide analysis of soybean indicates the presence of fortyone GLYI and twenty-three GLYII proteins encoded by twenty-four and twelve genes, respectively. Expression of GmGLYI and GmGLYII genes was found to be modulated by different developmental stimuli and abiotic stresses. In the present study, expression of soybean glyoxalase genes has been analyzed at five distinct developmental stages and eight tissues as well as in response to various biotic stresses using publicly available microarray data. GmGLYI-7, GmGLYI-10/21, GmGLYII-1/8, and GmGLYII-6 showed high level of expression at all the developmental stages and tissues with few exceptions; while GmGLYI-6/9, GmGLYI-20, GmGLYII-6 and GmGLYII-10 showed most biotic stress specific modulation. However, the observed alteration is highly complex that depends on the duration and type of infected organisms.
Glyoxalase, Glycine max, Biotic stress, Microarray.
MG: Methylglyoxal; GLYI: Glyoxalase I; GLYII: Glyoxalase II; GSH: Reduced Glutathione; DPI: Day post infection; HPI: Hour Post Infection
乙二醛酶系统是一种进化的强烈反对ved pathway, consists of two thiol-dependent enzymes; glyoxalase І (GLYI) and glyoxalase II (GLYII). Glyoxalase pathway detoxifies a cytotoxic metabolite, methylglyoxal (MG) to its non-toxic form, D-lactate with the help of reduced glutathione (GSH). Free MG and GSH forms non-enzymatic hemithioacetal adduct, which isomerized into S-D-lactoyl-glutathione (SLG) by the action of GLYI. GLYII further hydrolyzes SLG into D-lactate and regenerates GSH [1]. Glyoxalase pathway was first discovered more than hundred years back (1913) in two different organisms, rabbit and dog [2]. Later, glyoxalase genes/proteins have been reported from a wide range of species includingEscherichia coli,yeast, fungus, plants, mammals,Homo sapiensand found to be highly conserved [2,3].
Chemically methylglyoxal (MG) is a reactive α,β-dicarbonyl aldehyde. It can disrupts cellular functions by forming adduct with cellular macromolecules such as DNA, RNA and proteins; and forms advance glycation end products (AGEs) [2]. Intracellular MG level was found to be increased under various adverse conditions in all living organisms and highly correlated with stress induced pathogenesis. Increased level of MG found to be connected with various diseases; diabetes, cardiovascular diseases, cancer, nephropathy, retinopathy, neuropathy [4]. Due to reduced GLYI activity, MG gets accumulated in diabetes patients, that promotes endothelial cells inflammation, dysfunction and vascular damage [5]. Enhanced accumulation of MG has been reported in various plant species under stresses that lead to stress induced senescence [6,7] . Accumulation of SLG is also toxic to living system as it could inhibit DNA synthesis [8]. Thus glyoxalase system has played a vital role to combat MG accumulation under adverse conditions. There are numerous studies in plants that have shown the ability of glyoxalase enzymes in conferring tolerance against multiple stresses by maintaining intracellular MG level [6,9,10] . Moreover, glyoxalase enzymes have been reported to be involved in different other important cellular functions, such as cell division and proliferation, microtubule assembly and protection against oxoaldehydes toxicity [11]. Thus, glyoxalase pathway has been considered as “marker for cell growth and division”. Over-expression ofGLYIand/orGLYIIgenes in transgenic plants enhanced tolerance against multiple abiotic stresses [6,9,12,13] . For this, both accumulation of MG and glyoxalase gene expression/enzyme activity are considered as potential plant stress biomarker [14].
Glyoxalase genes/proteins have been extensively studied in various micro-organisms and mammals, but limited information available about plant glyoxalases. The first genome wide analysis of glyoxalase genes have been reported from two model plants-ArabidopsisthalianaandOryza sativa[10]. Authors predicted the presence of 11 potentialGLYIand 3GLYIIgenes in rice; and 11 putativeGLYIand 5GLYIIgenes in Arabidopsis [10]. Expression analysis ofAtGLYI/II andOsGLYI/II genes showed their developmental and tissue specific regulation. Further, qRT-PCR analysis performed in two contrasting rice genotypes, i.e., IR64 and Pokkali identified two highly stress inducible members in rice-OsGLYI-6andOsGLYI-11. In soybean, the genome wide analysis of glyoxalase genes has been reported recently [15]. Comprehensive soybean genome analysis identified 41 putative GLYI and 23 putative GLYII proteins encoded by 24 and 12 genes, respectively. Expression profiling of these genes was carried out at different tissues and developmental stages as well as under two major abiotic stresses- salinity and drought using publicly available RNAseq data [15]. Present study revealed the biotic stress inducible expression ofGmGLYIandGmGLYIIgenes.
Identification of Corresponding ProbesetforGmglyiandGmglyiiGenes
To retrieve the expression data of soybean glyoxalase genes in various developmental stages, tissues and biotic stresses from Genevestigator, the corresponding probesets for all GmGLYI and GmGLYII genes are required.
Probesets为被确定using NetAffx Analysis Center (http://www.affymetrix.com/analysis/index.affx?navMode=cat530006&aId=netaffxNav) online Probe Match tool using respective CDS sequences as query. Probe with the highest value was considered in case of genes with more than one probeset; while more than one gene showing same probeset was considered as same transcriptional profile.
Microarray-Based Expression Analysis of Soybean Glyoxalase Genes
Normalized and curated expression data at various developmental stages, tissues and biotic stresses were retrieved from publicly available Affymetrix soybean genome array using Genevestigator (https://www.genevestigator.com/gv/plant.jsp) [16,17] . Genevestigator processed a large set of manually curated and quality-controlled microarrays and RNAseq data from a large variety of samples, tissues and conditions. Various developmental stages of soybean such as germination (61 samples), seed growth (504 samples), flowering (3 samples), fruit formation (45 samples), and bean development (89 samples) are analyzed. Similarly, expression was checked at eight tissues such as root (12 samples), seedlings (2 samples), hypocotyl (14 samples), leaf (3 samples), inflorescence (2 samples), seed (39 samples), embryo (8 samples), and endosperm (8 samples). Detailed information about the corresponding experiments has been provided asAdditional File-1.To analyze the expression in response to different microorganisms viz.Bradyrhizobiumjaponicum,Phakopsorapachyrhizi,Heteroderaglycines,Spodopteralitura, andRhizophagusirregularis,相对信号比的值被检索using Genevestigator with default parameters (using experiments with id: GM-00028, GM-00033, GM-00036, GM-00046, and GM-00048) and the log2 transformed fold change values were used to generate heap map. Detailed information of all these experiments has been included asAdditional File-2. All these experiments were replicated atleast thrice and fold change values were calculated by comparing with the respective control with same number of replicates. Heat maps with hierarchical clustering were generated using Multi Experiment Viewer software (http://www.tm4.org/mev/)with Manhattan distance metric method [18].
Expression analysis ofGmGLYIandGmGLYIIgenes at various developmental stages and tissues
植物是一个高度复杂的过程发展that needs integrated action of several pathways and genes [19]. Different stages of development represent different sub-sets of gene population. Soybean has five distinct developmental stages-germination, shoot growth, flowering, fruit formation and bean development. Expression of allGmGLYIandGmGLYIIgenes were analyzed using genevestigator (https://genevestigator.com/gv/doc/intro_plant.jsp) based micro-array data and plotted in a scatter diagram with standard deviation. Some genes share same probe set which were considered as same transcript profile; while some genes did not have a respective probe, and thus no gene expression was available. Among the total of 24GmGLYIand 12GmGLYIIgenes; expression data was found for 19GmGLYIand 8GmGLYIIgenes. This analysis revealed thatGmGLYI-7 andGmGLYI-10/21 maintained high transcript abundance at all the developmental stages except flowering, whereGmGLYI-20 showed maximum expression(Figure 1A). Similarly,GmGLYII-1/8 andGmGLYII-6 showed higher level of expression at all the developmental stages except flowering, whereGmGLYII-4/5 showed highest abundance(Figure 1B). Further, expression of allGmGLYIandGmGLYIIgenes were analyzed at eight distinct soybean tissues, such as- root, seedling, hypocotyl, leaf, inflorescence, seed, embryo, and endosperm.GmGLYI-7 andGmGLYI-10/21 maintained high transcript abundance at all the analyzed tissues except leaf, whereGmGLYI-23 showed the highest peak(Figure 1C). In case ofGmGLYIIfamily members,GmGLYII-1/8 showed similar pattern with maximum level of expression at all the analyzed tissues except root, whereGmGLYII-6 showed the highest expression(Figure 1D). This indicates the development and tissue specific alteration ofGmGLYIandGmGLYIItranscripts. Genes with very high level of expression at all these stages and tissues indicates the constitutive expression due to their imperative role in regular cellular/metabolic processes.
Figure 1:Expression Profiling Of Gmglyi andGmglyii Genes Under Different Developmental Stages and Tissues Of Soybean.
Expression analysis ofGmGLYIandGmGLYIIgenes in response to various biotic stresses
Abiotic stress specific response ofGmGLYIandGmGLYIIgenes was analyzed previously [15]. Apart from plant development, tissue specificity and abiotic stress response, expression of glyoxalase genes have been reported to be modulated by biotic stress too [14]. Expression of riceGLYIgene was found to be down-regulated in response toXanthomonasoryzaepv. oryzae orPyriculariagriseainfection [20]. Proteomic comparison ofAspergillusflavus-resistant and -susceptible maize kernel embryo showed drastic up-regulation of GLYI protein in the resistant embryos [21]. Fungal infection significantly enhances intercellular methylglyoxal (MG) level in susceptible genotypes and thus high GLYI activity is crucial. Similarly, induction of GLYI activity was observed in rice andBrassicain response to brown planthopper (Nilaparvatalugens) andSclerotiniasclerotioruminfection, respectively [21,22] . All these studies convincingly indicate the role/involvement of glyoxalase genes/proteins in biotic stress modulation pathways.
To gain deep insights into the regulation of glyoxalase genes under biotic stress, the expression profiles of allGmGLYIandGmGLYIIgenes were analyzed in response to various types of microorganisms using publicly available microarray data. Expression ofGmGLYIandGmGLYIIgenes were checked withBradyrhizobiumjaponicum,Phakopsorapachyrhizi,Heteroderaglycines,Spodopteralitura, andRhizophagusirregularisinfection at different time interval.B. japonicumis a legume-root nodulating nitrogen-fixing species which improves crop yields.P. pachyrhizicauses asian soybean rust disease that creates prominent dark to reddish-brown lesions.H. glycinesare devastating soybean pest that infects mainly in the roots. It causes suppression of growth, root and stem necrosis and leaf chlorosis, and ultimately loss in the total yield.S. liturais a noctuid moth which is considered a major pest of many crops including soybean.R. irregularisis an arbuscularmycorrhizal fungus that is commonly used in scientific studies to check the effects of arbuscularmycorrhizal fungi on plant growth and development. In response toB. japonicuminfection (6 to 48 hpi)GmGLYIgenes showed scattered pattern up/down-regulation of expression; while allGmGLYIIgenes showed down-regulation till 18 hpi and up-regulation on/after 24 hpi(Figure 2A and 2B). In response toP. pachyrhiziinfection for 12 to 248 hpi, the level of modulation was the highest(Figure 2A and 2B).GmGLYI-6/9, GmGLYI-19/24andGmGLYI-20showed high up-regulation at the early stage (12 hpi) as well as late stage (216 and 248 hpi) of infection; whileGmGLYI-10/21, GmGLYI-18andGmGLYI-23showed strong down-regulation(Figure 2A). In case ofGmGLYII,GmGLYII-7/9, GmGLYII-6andGmGLYII-10high level of down-regulation at the early stage (12 hpi) as well as late stage (216 and 248 hpi) of infection, and the rest members showed medium to low level of alterations(Figure 2B). Similarly, in response toH. glycines(2-10 dpi), most of theGmGLYItranscripts level were found to be down-regulated, exceptGmGLYI-20that showed strong induction at later stages (5 and 10 dpi)(Figure 2A).然而,大多数的GmGLYIItranscripts level were found to be up-regulated, whileGmGLYII-6andGmGLYII-10showed strong reduction at early stage (2 dpi)(Figure 2B). In response toS. lituraandR. irregularis,GmGLYI-17/22andGmGLYI-2/13showed medium to high level of up-regulation, while rest of theGmGLYIandGmGLYIIgenes showed minimum alterations(Figure 2A and 2B). This indicates the diversified effect of different pathogenic and symbiotic micro-organisms on the expression ofGmGLYIandGmGLYIItranscripts.
Figure 2:Expression Analyses of Soybean Glyoxalase Genes in Response to Different Pathogenic and Symbiotic Microorganism Infection.
Taken together, this study reveals the strong correlation between biotic stress and alteration of glyoxalase transcripts. As generation of MG is a common phenomenon of stressed plant, transcript of glyoxalase genes (MG metabolizing pathway) is highly modulated by all stresses. Accumulation of MG has been reported previously in maize susceptible plants in response toAspergillusflavusinfection [21]. Thus, the observed information will encourage researcher to carry out further functional characterization of glyoxalase genes from various plant species in response to biotic stress.
Author acknowledges Shahjalal University of Science and Technology, Sylhet, Bangladesh for providing the logistic support and Department of Biochemistry and Molecular Biology of the same University for providing the laboratory space. Author acknowledges the support of Genevestigator.
