Figures
Abstract
Pathogenesis related (PR) proteins are low molecular weight family of proteins induced in plants under various biotic and abiotic stresses. They play an important role in plant-defense mechanism. PRs have wide range of functions, acting as hydrolases, peroxidases, chitinases, anti-fungal, protease inhibitors etc. In the present study, an attempt has been made to analyze promoter regions of PR1, PR2, PR5, PR9, PR10 and PR12 of Arabidopsis thaliana and Oryza sativa. Analysis of cis-element distribution revealed the functional multiplicity of PRs and provides insight into the gene regulation. CpG islands are observed only in rice PRs, which indicates that monocot genome contains more GC rich motifs than dicots. Tandem repeats were also observed in 5’ UTR of PR genes. Thus, the present study provides an understanding of regulation of PR genes and their versatile roles in plants.
Citation: Kaur A, Pati PK, Pati AM, Nagpal AK (2017) In-silico analysis of cis-acting regulatory elements of pathogenesis-related proteins of Arabidopsis thaliana and Oryza sativa. PLoS ONE 12(9): e0184523. https://doi.org/10.1371/journal.pone.0184523
Editor: Rajesh Mehrotra, Birla Institute of Technology and Science, INDIA
Received: February 14, 2017; Accepted: August 27, 2017; Published: September 14, 2017
Copyright: © 2017 Kaur et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: Amritpreet Kaur received funding from university grant commission (UGC), New Delhi, India.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Plants are persistently under the threat of several pathogens like bacteria, viruses, fungi, nematodes and other threats. However plant pathogen interactions are extremely intricate and cause majority of plants to become impervious to the vast majority of pathogens. Further these interactions exhibit specific responses that permit only a few pathogens to colonize and spread disease [1–2]. Early recognition of a pathogen is an indispensable step for disease resistance in plants, which is followed by an activation of a series of defense responses during the interaction. During incompatible interactions, avirulence pathogen proteins (Avr) interact with host resistance (R) genes bringing about a series of defense responses such as: accumulation of Reactive Oxygen Species (ROS); enhancement in abscisic acid (ABA), salicylic acid (SA), jasmonic acid (JA), auxins and gibberellins; synthesis of pathogenesis related (PR) proteins; phytolexins accumulation; and hypersensitive response (HR) induction. Consequently, plants don’t develop disease symptoms and are safe. In case of compatible interactions, the virulent pathogen invades host machinery and results in systemic infection which prompts the development of disease symptoms [1–3].
PR proteins are a collective set of low molecular weight proteins which accumulate under various biotic and abiotic stresses and under specific physiological conditions like pollen development, leaf senescence, fruit development and ripening [4–7]. Such proteins have been considered to perform a number of functions, acting as transcription factors, protease inhibitors, enzymes involved in hydrolysis and many are associated with various metabolic pathways [4, 8]. PRs were isolated for the first time from tobacco leaves (Nicotiana tabacum) infected with tobacco mosaic virus [9] and subsequently reported from many other plant species including A. thaliana, alfalfa, barley, bean, carrot, chickpea, grape vine, maize, pepper, pearl millet, rice, rubber, soybean, sunflower, sorghum, tomato and wheat [10]. PR proteins have been characterized and classified into 17 families based on the sharing of amino acid sequences, serological relationships, and enzymatic or biological activity [11].
PR1 was the first PR protein to be discovered and has a molecular weight of 14 to 17 kD and acts as a molecular marker for systemic acquired resistance response. It has antifungal activity. PR2 proteins are β-1, 3-glucanases, and their molecular mass ranges from 33 to 44 kDa. They comprise large and highly complex gene families involved in pathogen defense as well as a wide range of normal developmental processes. They are induced in response to wounding or infection by viruses, bacteria and fungi. β-1,3-glucanases degrade pathogen’s cell walls by cleaving β-1,3-glucosidic bonds in β-1,3-glucan, a major component of fungal cell wall. PR3 proteins (chitinases) have molecular mass in the range of 15–43 kDa. They cleave the chitin polymers in fungal cell wall, resulting in a weakened cell wall and making fungal cells osmotically sensitive. PR4 proteins are chitin binding proteins having molecular mass between 9–30 kDa. These proteins bind to chitin, and play an important role in enhancing the chitinase activity [6, 12]. Thaumatin-like proteins (PR5) possess molecular mass between 18–25 kDa. These can act as antifungal; glucanase and xylanase inhibitors; and α-amylase and trypsin inhibitors. They are also known to be induced during wounding and by insect feeding; especially by phloem feeding insects [13]. Proteinase inhibitors (PR6) and endoproteinases (PR7) are highly stable defensive proteins of plant tissues that are both developmentally regulated and induced in response to insect and pathogen attacks. PR9 (peroxidase) catalyzes cross-linking of macromolecules in plant cell wall. It also produces free radical like H2O2 against a wide range of pathogens [6]. PR10 are ribosome inactivating proteins, known to inhibit translation in fungi. These proteins protect plant proteins and other cellular structures during dormancy, salinity or cold stress [14]. PR12 (defensin) are small cysteine rich peptides providing protection against a broad range of organisms. They are known to inhibit protein synthesis, enzyme activity and ion channel function [15]. PR15 and PR16 catalyze oxidation of oxalates by molecular oxygen, yield CO2 and H2O2. They have role in plant development, defense, signaling, differentiation and apoptosis [16].
For a number of PR proteins, activities are known or can be deduced. The majority of PRs (e.g. PR1, PR2, PR3, PR4, PR5, PR7, PR12, PR13 and PR14) possess antifungal activity, whereas, PR8 and PR11 are classified as endochitinases. PR15 and PR16 are oxalate oxidase and oxalate oxidase- like proteins, respectively [4, 16]. However, very little is known about molecular mechanism of gene expression of PR genes. In one study, Lodhi et al., 2008 [17] deduced a relationship among architecture of promoter sequence, positioning of nucleosome and expression of PR-1a in tobacco. Therefore, study of gene expression regulation of PR proteins is a crucial step in understanding the molecular mechanisms of plant defense response. Transcription regulation involves association between transcription factors and particular cis-acting regulatory elements (CAREs) of a specific gene involved in plant defense response [18]. CAREs are short regulatory motifs (5–20 bp) present in the promoter regions of target genes (typically, non-coding DNA). Promoters play an important role in controlling gene expression. Multiple CAREs such as TATA box, GC box, CAAT box contain coupling sites for transcription factors, enhancers and repressor elements required for proper spatiotemporal expression of genes [19]. Cis-acting regulatory elements are essential transcriptional gene regulatory units as they control various stress responses. Recent advancements in such experimental techniques as RNA interference, microarrays, RNAseq and others have allowed identification and investigation of promoter regions of target genes but these techniques are expensive and technically challenging. Therefore, computational methods are being used to search the promoter regions for different cis-elements responsible for the regulation of the genes [18]. Different computer programs can also be used to look for known cis-elements and to study their organization. Such web-based tools as PLACE [20], PlantCARE [21], AGRIS [22], TRANSFAC [23] and PlantPAN [24] have been developed for the analysis of cis regulatory elements in plant genes.
Examination of CAREs within the promoter sequences of PR genes as well as their combinatorial effects, will lead to better comprehension of regulation of their gene expression. Understanding of cis-elements can also allow us to effectively change the expression pattern of a gene in desired way, which further can provide new ways for the plant genetic engineering technology for protection of crops against biotic and abiotic stresses. To the best of our knowledge, no work has been reported on cis-elements of Arabidopsis thaliana and Oryza sativa PR. Therefore, the present study was planned, to characterize cis-acting regulatory elements (CAREs) of PR classes 1, 2, 5, 9, 10 and 12 with respect to their occurrence and putative role in model plants, Arabidopsis thaliana and Oryza sativa. We also tried to validate our in-silico work with wet lab studies wherever available.
Materials and methods
Search for PR genes of A. thaliana and O. sativa and their structure analysis
Gene sequences of PR-1, PR-2, PR-5, PR-9, PR-10 and PR-12 of A. thaliana and O. sativa were retrieved from The Arabidopsis Information Resource (TAIR) (https://www.arabidopsis.org/) and Rice Genome Annotation Project (RGAP) (http://rice.plantbiology.msu.edu/), respectively. All the gene sequences were verified against Phytozome v 11.0 (https://phytozome.jgi.doe.gov/pz/portal.html) and Plant Genome and System Biology (PGSB) (http://pgsb.helmholtz-muenchen.de/plant/plantsdb.jsp)databases. Comparative analysis of PR gene sequences was performed using MatGAT 2.02 (http://ww3.bergen.edu/faculty/jsmalley/matgat.html) [25] to find percentage similarity and NCBI-Conserved domain database (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) [26] to find conserved domains. Chromosome maps of A. thaliana and O. sativa PR genes were constructed by Chromosome Map Tools available at TAIR (https://www.arabidopsis.org/jsp/ChromosomeMap/tool.jsp) and Oryza base (http://viewer.shigen.info/oryzavw/maptool/MapTool.do), respectively. The intron/exon organization of splice variants of PR genes of both A. thaliana and O. sativa was retrieved from TAIR and Rice Genome Annotation Project respectively. Pearson correlation analysis was performed to find out the relationship between gene lengths of different PRs of A. thaliana and O. sativa. Gene Structure Display Server (http://gsds.cbi.pku.edu.cn/) [27] was used to visualize intron/exon organization of genes.
Retrieval of promoter regions and analysis of cis-regulatory elements
Promoter sequences (1.5 kbp upstream of translation start site) of each A. thaliana and O. sativa PR gene under consideration were retrieved from the respective databases i.e. TAIR and RGAP. The tools PlantCare (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [21] and AGRIS (http://arabidopsis.med.ohio-state.edu/AtcisDB/) [22] were used for scanning of cis-elements present in promoter regions of PR genes of A. thaliana whereas, PlantCare and PLACE (http://www.dna.affrc.go.jp/htdocs/PLACE/) [20] were used for identification of cis-elements in promoters of O. sativa PR sequences. The cis-elements so obtained were compared with each other and discussed in light of literature available. PlantPAN (http://plantpan2.itps.ncku.edu.tw/) [24] was used for the analysis of CpG/CpNpG islands and tandem repeats.
In silico analysis of PR genes expression
Genevestigator (https://genevestigator.com/gv/) [28] was used for performing gene expression analysis.
Results and discussion
Search for PR genes of A. thaliana and O. sativa and their structural analysis
Characteristics of 6 PR genes of A. thaliana and O. sativa retrieved from TAIR and RGAP databases are given in Tables 1 and 2. The gene size of 6 PR genes of A. thaliana varied from 535 bp in AtPR12 to around 1545 bp in AtPR9 and in O. sativa, gene size varied from 738 bp in OsPR12 to 2347 bp in OsPR2. Characteristics of each PR gene sequence are given in Tables 1 and 2 and include locus, chromosome number, strand, transcript count, transcript id, gene length, CDS length, number of exons, and protein length. MatGAT tool was used to compare PR gene sequences of A. thaliana and O. sativa, which revealed percentage similarity (evolutionary distance) among different PRs to range from 39.4% for PR10 to 67.3% for PR2 (S1 Table). PR10 proteins are coded by multigene families and shows higher inter specific variation. However, all of them possess conserved glycine rich loop in their sequence [29]. PR gene sequences from A. thaliana and O. sativa were subjected to comparative analysis for their domain architecture using CDD from NCBI. The analysis revealed the presence of SCP_PR1_like domain in PR1; glycol_hydro_1 domain in PR2; GH64-TLP-SF in PR5; plant_peroxidase_like domain in PR9; SRPBCC in PR10 and gamma-thionin in PR12 of A. thaliana and O. sativa, respectively. Outcome of this analysis is the presence of conserved domains in their sequences, which indicate that they are homologs. In-silico chromosome mapping of PR genes of A. thaliana and O. sativa is presented in Fig 1 (a and b respectively). PR genes of A. thaliana have been shown to be distributed on 4 out of 5 chromosomes and of O. sativa are distributed on 5 out of 12 chromosomes. In A. thaliana, chromosome 1 harbored 3 PR genes, whereas chromosome 5 showed no PR gene. In case of O. sativa, chromosome number 7 possessed 2 PR genes, while chromosome numbers 3, 5, 6, 8, 9, 10, 11 showed no PR genes.
The position of exons/introns and intronic phase distribution are important characteristics for gene structure analysis. Introns have been classified into three phases: phase 0 introns are present between two consecutive codons; phase 1 introns are present between the first and second nucleotide of a codon and phase2 introns are present between the second and third nucleotide of a codon [30]. Fig 2a shows the exon-intron and intron phase arrangement of PR genes of A. thaliana and O. sativa. PR genes of A. thaliana and O. sativa show almost same intronic phase distribution except PR9 gene. Exon count is same in the PR genes of A. thaliana and O. sativa. Table 3 shows the position and length of introns-exons organization in A. thaliana and O. sativa PR genes. OsPR2 contains intron of 998 nucleotides, whereas AtPR2 intron length is 94 nucleotides. AtPR1 and OsPR1 showed the absence of introns in their gene sequences. No variation has been observed in the exon lengths of A. thaliana and O. sativa PR12 gene i.e. in PR12 of both species exon 1 is 64 bp long, while exon 2 is 179 nucleotides long. The intron length however varies slightly like AtPR12 intron is of 107 bp while OsPR12 intron is 100 bp in length. The 5’ and 3’ UTR regions of AtPR12 and OsPR12 varied considerably which contributed to difference in overall length of PR12 gene in two species. Pearson correlation analysis for gene lengths of 5 PR genes (PR1, PR5, PR9, PR10 and PR12) of A. thaliana and O. sativa revealed a highly significant coefficient of correlation (r = 0.996 at p<0.001) (S1 Fig). Since the intron length of PR2 gene varied considerably among the two species studied (A. thaliana intron length 94 bp and O. sativa intron length 998 bp long), this gene was excluded from correlation analysis. The splice variants of PRs of A. thaliana and O. sativa were also analyzed as shown in Fig 2 (b and c); and except in AtPR10, no other splice variant has been observed in the above analysis.
(a) Blue boxes represent untranslated region. Exons are represented by green boxes and connecting lines represent introns. (b) exon-intron arrangement among splice forms in A. thaliana and O. sativa.
Retrieval of promoter regions and analysis of cis-regulatory elements
Promoter sequences up to1.5 kbp upstream from the translation start site of each PR gene of A. thaliana and O. sativa were scanned using PlantCare program for the identification of cis-acting regulatory elements (CAREs). The study revealed a total of 55 CAREs in A. thaliana whereas, 48 CAREs were identified in O. sativa. The length of cis-elements varied from 4–13 bp in A. thaliana and 4–10 bp in O. sativa (Fig 3).Majority of cis-elements possessed the length of 6 and 7 bp. Cis-elements are grouped into different functional categories as shown in Fig 4. Stress responsive cis-elements were found to be maximum, and out of those, light responsive elements constituted 34%. Fig 5 depicts the frequency of occurrence of different cis-elements at different positions in the 1.5 kbp of both reverse and forward strands of A. thaliana and O. sativa PR genes. In a case of O. sativa, most of the cis-elements lied between 51-100bp at the forward strand and 351–750 bp at reverse strand. In A. thaliana, maximum cis-elements lied within the range of 351–650 bp (forward strand) and 701–950 bp (reverse strand). The PR gene sequences of A. thaliana and O. sativa were also scanned for the presence of cis-elements by AGRIS and PLACE, respectively. Cis-elements so obtained were compared among themselves, and the results are presented in Table 4. The number of cis-elements including TATA and CAAT box was present in higher amount in A. thaliana PR genes as compared to that of O. sativa PR genes. Fig 6 shows the presence of different cis-elements along with their frequencies in A. thaliana only, O. sativa only and both A. thaliana and O. sativa PR genes. The role of different cis-elements in stress, hormonal regulation and cellular development are discussed in details below.
Stress responsive cis-elements.
Different elements associated with different stress responses such as oxidation, defense, light, cold, drought, dehydration, heat, low temperature, wound were observed in A. thaliana and O. sativa PR gene sequences. Oxidative stress is a component of many abiotic and biotic stress conditions such as drought, high temperature stress, UV-B radiation, salinity, metal toxicity, chilling and plant pathogen interactions. PR proteins are strongly induced in response to different types of stresses including wounding or infection by pathogens. Under normal conditions, production of Reactive oxygen species (ROS) in plants is very low and under various environmental stresses, ROS is drastically increased in plants disturbing the normal balance of superoxide anion (O2−), hydroxyl radical (·OH), singlet oxygen (1O2) and hydrogen peroxide (H2O2) in the intracellular environment [31, 32]. Microarray expression analysis of A. thaliana revealed a range of cis-elements responsive to different types of ROS. Such elements have been categorized into two different categories: common ROS-related e.g. TATCCAT/C-motif, GCN4_motif and G-box and ROS-specific element like W-box [33]. Some other elements like H-box, ethylene-responsive GCC elements, salicylic acid, ethylene, abscisic acid and calcium are known to contribute to the response to oxidative damage in A. thaliana [34].
TATCCAT/C-motif is an amylase element representing sugar repression responsiveness. It plays an important role in GA-regulated expression. AtPR9 (peroxidase) shows the presence of this motif. In A. thaliana, peroxidases play an important role in generating H2O2 during the defense response and also provide resistance against a wide range of pathogens. Peroxidases also play a vital role in leaf expansion [35].
GCN4_motif (TGTGTCA) is an essential cis-element required for an endosperm specific gene expression. AtPR12 shows the presence of GCN4_motif. AtPR12 has a role in protecting germinating seeds and developing seeds [36].
G-box (CACGTG) element is involved in response to light, abscisic acid, methyl-jasmonate and anaerobiosis and has a role in ethylene induction as well as in seed specific expression. It is also known as ABRE (ABA-responsive element) [37, 38]. It has been shown to be present in all A. thaliana PR genes except in AtPR2 and in all O. sativa PR genes (Fig 5a and 5b).
W box (TTGACC) is an elicitor responsive cis-element, present in AtPR 9 and AtPR 10. W-boxes are found to interact with transcription factors belonging to WRKY family. W box regulates the expression of defense-related (PR10) genes and has role in biotic and abiotic stresses; seed dormancy; senescence etc [39–40]. During stress response, AtPR10 is induced by ABA, ethylene, jasmonic acid and salicylic acid. This gene may be induced by ROS and may act as a protinase against cellulases and pectate lyases of the pathogen [29]. Increase in ROS level especially H2O2 has been shown to increase the PR10 in plants [41]. The presence of W box in AtPR9 indicates its role in senescence [40].
In addition to the above mentioned responsive elements, other oxidative stress responsive elements like AREs, ethylene-responsive GCC elements, ERE and H-box are also present in promoters of PR gene sequences (Fig 5a and 5b). AREs (Anaerobic responsive elements) are essential for anaerobic induction, present in AtPR5, AtPR9 and AtPR10. AREs are bipartite elements consisting of GC and GT motifs. GT motif resembles AtMYB2 transcription binding site, which is drought and low oxygen induced element [42].GCC-box, ethylene-responsive element is necessary for high-level jasmonate-mediated regulation of PR12 expression during plant defense response [43–44]. GCC element is also associated with the expression of many genes involved in different kinds of abiotic and biotic stresses. H-box is a root specific regulatory element present in AtPR9 and OsPR10. It regulates defense genes by elicitors and other stress stimuli [45]. The TGACG motif, also known as ‘as1 element’, is another well characterized cis-element present in plants. TGACG motif is methyl jasmonate responsive element present among A. thaliana and O. sativa PR gene sequences. The transcription of TGACG mediated PR sequences is regulated by binding of BZIP TGA factor to TGACG element [46]. TC- rich repeats are seen in AtPR1, AtPR2, OsPR2 and OsPR9 and have role in stress and defense responsiveness. Whereas, G-box [47] and TATCCAT/C-motif are also involved in regulating defense responses.
A number of cis-elements associated with light stress include ACE, AE-box, ATC-motif, ATCT motif, Box I, Box II, BoxW1, Box4, Box S, CATT motif, CG motif, Chs-CMA1a, Chs-unit 1 ml, G box, GA motif, GAG motif, GATA-motif, Gap box, GT1 motif, I-box, MNF1, Sp1, TCCC motif and TCT motif. Among these some elements such as ACE, AE-box, Box4, Box S, CATT motif, G box, I-box, Sp1, TCCC motif and TCT motif are present in both A. thaliana and O. sativa PRs, whereas Chs-unit 1 ml and Gap box (OsPR10); Box II (OsPR5); HSE (OsPR12); and C repeat/ DRE (OsPR1) are present only in O. sativa. GATA-motif and Chs-CMA1a motif (AtPR10); I-box and Sp1 (AtPR2); Box S (AtPR5); C repeat/ DRE (AtPR5); LTR (AtPR10); and MBS (AtPR1, AtPR2, AtPR5, AtPR9 and AtPR12) are present only in A. thaliana PRs.
Cis-elements in hormonal regulation.
Motifs involved in hormonal regulation were found to be second largest in number after stress responsive motifs present in PRs. Few motifs such as ABRE (abscisic acid), CGTCA and TGACG (methyl-jasmonate), GCC-box (ethylene), TCA element (salicylic acid) (Fig 5c) were present in both A. thaliana and O. sativa PRs. Some motifs were only limited to O. sativa PR genes, like TATC and GARE motifs are gibberellin-responsive elements present in OsPR2 and OsPR12, respectively. Abscisic acid responsive elements such as motif II b (OsPR2 and OsPR9), motif lib (OsPR5) and CE1 (OsPR12) were also observed only in O. sativa PRs. PLACE tool shows the presence of E-box and DPBF CORE DCDC3 motifs in OsPR5. Auxin responsive elements, AuxRR-core and TGA- element were found in AtPR9 and AtPR10, respectively. A MeJA responsive element, JERE motif was found in AtPR10. The AGRIS tool revealed the presence of CACATG motif, which is the binding site for MYC2 transcription element and has a role in jasmonic and abscisic acid signaling have been reported in all A. thaliana PRs [64].
Role of cis-elements in cellular development.
Motifs involved in cellular development are relatively less in number as compared to hormonal and stress responsive elements. Fifteen types of motifs related to cellular development were found. These cis-elements include AC-I, AC-II, 5’UTR Py-rich stretch, CAT box, F-box, GCN4_motif, Skn-1-motif, O2-site, circadian, CCGTCC box, motif I, Plant AP-2 Like, H-box, ARE, MBSI (Fig 6a and 6b). AC-I and AC-II are involved in xylem specific expression. AC-I is present in OsPR2, whereas, AC-II is found in AtPR1 and AtPR12. 5’ UTR Py-rich stretch (AtPR1, AtPR12 and OsPR2) and Plant AP-2 Like (OsPR2, OsPR5 and OsPR9) confer high transcription level. CAT box (AtPR2 and OsPR9) and CCGTCC box (OsPR2, OsPR9 and OsPR12) have role in meristem-specific activation, F-box (AtPR9) confers role in cell cycle regulation (Fig 5a and 5b). GCN4_motif (AtPR12) and Skn-1-motif (AtPR5, AtPR9, AtPR10, AtPR12 and in all O. sativa PRs except OsPR5) are involved in endosperm expression. O2-site (OsPR10, AtPR9 and AtPR10) involved in zein metabolism and circadian motif (present in all AtPRs except AtPR12, OsPR2 and OsPR10) has role in circadian control. Motif I (OsPR2) and H-box (OsPR10) are root specific regulatory motifs. ARE (AtPR5, AtPR9, AtPR10) is essential for anaerobic induction and MBSI (AtPR2) is involved in flavonoid biosynthesis gene regulation (Fig 5c).
Calcium responsive cis-elements.
Calcium (Ca2+) is an intracellular regulator, consequential for many plant biological functions. Ca2+signaling is a paramount mechanism evolved in plants to defend themselves against pathogens. It is required for inducing defense-related genes and hypersensitive cell death. Calmodulins (CaM) interact with specific TF (WRKY, MYB, and NAC) families, and regulate the expression of defense genes, but the direct role of CaM in regulating plant defense genes has not been studied so far [65,66]. There are few specific promoter motifs which are regulated by calcium. These include ABRE or ABRERATCAL (Abscisic Acid Responsive—Element), C-Repeat/ DRE (Drought-Responsive Element), Site II, CAM box, CRT and W-box [67, 68] ABRERATCAL (MACGYGB where M = C/A, Y = T/C, B = T/C/G) is the binding site harbored by ABA-induced gene promoter [69]. In the present study, presence of calcium responsive cis-elements was identified by using PlantPAN. We observed the presence of ABRE related elements in all the A. thaliana PR sequences (Fig 7), whereas, no ABRE like elements was observed in O. sativa PR sequences.
Analysis of conservation of cis-elements in promoter regions of PR genes.
Analysis of Phylogenetic conservation of sequences involves the identification of conserved motifs across the genes. The goal of this work was the identification of cis-regulatory sequences conserved in promoters of PR genes of A. thaliana and O. sativa. PlantCare data was analyzed to reveal the conserved sequence motifs in promoters of A. thaliana and O. sativa PR genes (S2 Table). CAAT and TATA box act as binding sites for transcription factors. CAAT box, important in core promoter activity is almost conserved in all PRs of A. thaliana and O. sativa, whereas, TATA box is conserved in most of the PRs except PR5 and PR9. Zuo and Li (2011) also showed the presence of TATA-less promoters in plant genome [70], which indicates that TATA box is not conserved in all the plant genomes. G-box is conserved in all PR promoter regions of A. thaliana and O. sativa except PR2. Ishige et al., [71] examined 11 different G-box tetramers in regulation of GUS gene expression and found each G-box sequence influenced gene expression in different ways. One of the G-box sequences, G-box 10 was shown to confer high level constitutive expression in roots, leaves and seeds. MBS cis-element is present in PR2, PR5, PR9 and PR12. MYB transcriptional factor requires MBS for the gene expression of drought inducible genes. CGTCA motif is also found to be conserved in PR5 and PR9 of A. thaliana and O. sativa, involved in methyl jasmonate (MeJA) responsiveness. It activates series of defense mechanisms in response to different abiotic stresses like drought, salinity and low temperature. MeJA motif in the 5’ UTR of PR genes infers a possible role in pathogen stress or wound responses. PR1 and PR5 show the presence of ABRE motif, a positive regulator of abscisic acid signaling under drought stress and high salt condition in the vegetative tissues of plants.
There are few CAREs which are unique to A.thaliana PR proteins (S2 Table). CAREs like A box and CCGTCC box are development related motifs involved in activation of meristem specific expression. GARE and TATC box are cis-acting regulatory elements involved in gibberellin responsiveness. CE1, Motif II b and motif lib are abscisic acid responsive elements unique to A. thaliana PRs. Cis-elements present only in PRs of O. sativa are mentioned in (S2 Table). F box is cis-element conserved in PR9; involved in regulating plant defense responses in response to biotic and abiotic stresses. GCN4 motif present in PR12 has role in endosperm specific gene expression.
Tandem repeats and CpG/CpNpG analysis by PlantPAN
The eukaryotic genome has a wide number of DNA repeats and these repeats have a role in genome evolution [72]. Repetitive DNA may be interspersed in a tandem configuration throughout the genome or may be restricted at some specific location. DNA tandem repeats according to their repeated unit length can be classified into three groups: (i) microsatellite—repeat unit less than 9 nucleotide in length (ii) minisatellite—with 6–100 bp (usually around 15 bp) repeats (iii) megasatellite—tandem repeats of longer units, with length more than 135 nucleotides [18, 73]. Microsatellites are codominant, abundant, multi-allelic, so can be used as molecular markers, in linkage mapping and gene tagging [74].The 1.5kbp upstream promoter region of PR genes revealed the presence of DNA tandem repeats. We found tandem repeat units in three AtPRs (AtPR1, AtPR9 and AtPR12) and one in OsPR (OsPR1) (Table 5). AtPR1 and AtPR9 contain mononucleotide repeats with a repeat size of 1 nucleotide. AtPR12 contain minisatellite. OsPR1 has tetranucleotide repeat with the repeat size of 4 nucleotides. Variation in length of tandem repeats in promoter region could be due to numerical changes like addition and deletion of transcription factor binding sites [75].
Epigenetic modifications like DNA methylation, chromatin remodeling and histone modification are heritable changes in gene expression which influence the phenotype [76]. Among these, DNA methylation is important and affects gene expression in plants and animals [77]. DNA methylation occurs at cytosine base, within CpG dinucleotide or may occur at CpNpG (N = A, C or T) sites [78]. CpG rich regions are named as CpG islands and to classify a genome region as CpG island, three conditions must be fulfilled (i) GC content should be above 50% (ii) length of CpG/CpNpG region should be greater than 200 bp (iii) ratio of observed-to-expected CpG dinucleotide number should be above 0.6 [18]. CpG islands are present at or near the gene’s transcription start site and they may regulate the tissue-specific gene expression [79]. DNA of plant species has been shown to contain more CpG dinucleotides than human DNA [80]. Methylation of cytosine at CpG islands has been shown to restrict the access of promoter region of genes to their transcription factors hence, preventing their expression [81]. Cytosine methylation patterns are not static; they change substantially with developmental state or with environmental conditions across the plant genome [82]. DNA methylation has also been shown to play an important role in plant embryogenesis, seed development, in regulating an immune response to infection by pathogens, environmental adaptation and stress resistance. Defects in methylation can cause defect in embryogenesis like abnormal cell division and seed viability reduction, developmental retardation, reduced plant size and partial sterility [77, 83–84]. CpG/CpNpG analysis revealed the occurrence of CpG/CpNpG islands in the second half of the promoter region (towards 3’ end) of all the OsPRs except in OsPR2 but none is identified in AtPRs (Table 6). The absence of CpG/CpNpG islands in the OsPR2 might be due to spontaneous deamination in the germline during evolution [85]. The study performed by Ferguson and Jiang [86] also showed that monocot genome contains more GC-rich content than dicots. The absence of CpG islands in A. thaliana as compared to O. sativa genome may be due to the difference in codon usage, GC-biased gene conversions and mutational biases prevalent in two species [87].
In silico analysis of PR genes expression
Meta-analysis of Genevestigator microarray datasets was performed on A. thaliana and O. sativa PR genes (Fig 8a and 8b). AtPR10 was a highly expressed gene in almost all stages of development except senescence, where, AtPR12 showed the highest expression during the senescence stage. AtPR1 and AtPR2 expression was maximal during developed rosette stage and minimal during the seedling stage. During the germinating seed stage, AtPR5, AtPR9, AtPR10 and AtPR12 showed same level of expression. Analysis of expression patterns of OsPRs revealed that OsPR10 showed the highest expression at dough stage whereas, OsPR12 appeared to be highly expressed at germination stage. OsPR9 and OsPR10 shared a similar kind of expression at many development stages except at heading, milk and dough stages. OsPR2 and OsPR12 also shared low expression during seedling, tillering and stem elongation stages.
Expression analysis of A. thaliana and O. sativa PR genes in response to hormonal treatments and abiotic stresses were also investigated by Genevestigator (Fig 9a and 9b). Based on microarray data available, it has been observed that in response to hormonal treatments, AtPR1, AtPR2 and AtPR5 were highly up-regulated by external application of salicylic acid, whereas; AtPR9 and AtPR10 were down-regulated. AtPR1 and AtPR2 were also up-regulated by IAA and ABA, respectively. In case of O. sativa PRs, OsPR9 and OsPR12 were minimally up-regulated by the application of jasmonic acid.
Furthermore, the expression profiles of AtPR genes in response to different abiotic stresses (cold, drought, heat, osmotic, salt and wounding) were also analyzed. Up-regulation was observed in AtPR1 and AtPR2 under drought; AtPR1, AtPR2, AtPR5 under osmotic stress and in wounding; and in AtPR2 under salt stress. Down-regulation was observed in AtPR1 and AtPR2 under heat stress; AtPR9 under salt stress and AtPR10 under cold stress. Expression profiles of OsPRs for cold, drought and salt were also retrieved. OsPR5 was highly and OsPR2 was minimally up-regulated, whereas OsPR1, OsPR9 and OsPR10 were down-regulated under drought stress. Under cold stress, OsPR10 was slightly up-regulated and OsPR1 was minimally down-regulated.
The present study, involved identification of different cis-elements associated with biotic and abiotic stresses present in 5’ UTR sequences of PR genes of A. thaliana and O. sativa. The effort has been made to validate the functions of PR’s cis-elements with the microarray data and with literature, wherever it is available. Microarray data indicated stage-specific expression of many AtPRs and OsPRs during different development stages. In case of O. sativa PRs, OsPR5 is induced under drought stress and in germinating stage; this may be due to the presence of AtMYC2 and RY-element cis motifs, respectively. Microarray data also revealed slight up-regulation of OsPR12 gene during the flowering stage, this can be linked up with the presence of A-box, cis-regulatory element involved in specific day-length to control flowering in O. sativa [88]. The oxidative stress responsive cis-element like ABRE has been shown to up-regulate OsPR2, OsPR5 and OsPR10 under stresses like drought and cold providing clue towards its diverse roles [89]. Various elicitor responsive elements like GCC box and W-box have been observed in AtPR5 and AtPR12. These were observed to be highly expressed during seed germination. AtPR1, AtPR2 and AtPR5 were up-regulated after wounding; this may due to the presence of wound responsive element (W-box) in their promoter regions. AtPR1 and AtPR2 are highly expressed under drought stress condition which maybe because of AtMYC2 element, which acts as a drought inducible element.
TGACG-motif present in AtPR5, AtPR9, AtPR10, OsPR2, OsPR5, OsPR9 and OsPR12 increases the production of secondary metabolites by arresting or delaying cell cycle in G1/S checkpoint. This motif tends to be present in the promoter region of Kip-related proteins (cell cycle inhibitor) and inhibits the active Cyclin-dependent kinase/Cyclin complex [90]. The combined data of cis-element and the meta-analysis of PRs provided insights into the role of many PRs, and hence this information can be employed in future studies of PRs from other plant species.
Conclusion
PR proteins play an important role in providing resistance to plants. Till now, no work has been reported on cis-elements present in A. thaliana and O. sativa PRs. Therefore, in the present work, an in-silico approach was followed to study the presence of cis-elements in PR genes. We also tried to validate our in-silico work with the wet lab studies wherever available. This work throws light on the promoter regions of PRs, which further can provide new ways for the plant genetic engineering technology for protection of crops against diseases.
Supporting information
S1 Fig. Pearson correlation analysis for gene lengths of AtPRs and OsPRs.
https://doi.org/10.1371/journal.pone.0184523.s001
(TIFF)
S1 Table. Percentage similarity between AtPRs and OsPRs.
https://doi.org/10.1371/journal.pone.0184523.s002
(DOCX)
S2 Table. CAREs along with their function conserved in both AtPRs and OsPRs; unique in AtPRs and unique in OsPRs.
https://doi.org/10.1371/journal.pone.0184523.s003
(DOCX)
Acknowledgments
The authors acknowledge Prof. A. K. Thukral of Guru Nanak Dev University, Amritsar (India) for helping in analyzing statistical data.
References
- 1. Che F, Entani T, Marumoto T, Taniguchi M, Takayama S, Isogai A. Identification of novel genes differentially expressed in compatible and incompatible interactions between rice and Pseudomonas avenae. Plant Sci. 2002; 162: 449–458.
- 2. Mazumder M, Das S, Saha U, Chatterjee M, Banerjee K, Basu D. Salicylic acid-mediated establishment of the compatibility between Alternaria brassicicola and Brassica juncea is mitigated by abscisic acid in Sinapis alba. Plant Physiol and Biochem. 2013; 70: 43–51.
- 3. Kundu S, Chakraborty D, Kundu A, Pal A. Proteomics approach combined with biochemical attributes to elucidate compatible and incompatible plant-virus interactions between Vigna mungo and Mungbean Yellow Mosaic India Virus. Proteome Sci. 2013:11–15.
- 4. Van loon LC, Van Strien EA. The families of pathogenesis-related proteins their activities and comparative analysis of PR-1 type proteins. Physiol and Mol Plant Pathology. 1999; 55: 85–97.
- 5. Graham MY, Weidner J, Wheeler K. Induced expression of pathogenesis-related protein genes in soybean by wounding and the Phytophthora sojae cell wall glucan elicitor. Physiol Mol Plant Pathol. 2003; 63: 141–149.
- 6. Hegde YR, Keshgond RS. Role of pathogenesis-related proteins in plant disease management—a review Agri Reviews. 2013; 34 (2): 145–151.
- 7. Jiang L, Wu J, Fan S, Li W, Dong L, Cheng Q, Xu P, Zhang S. Isolation and Characterization of a Novel Pathogenesis-Related Protein Gene (GmPRP) with Induced Expression in Soybean (Glycine max) during Infection with Phytophthora sojae. PLoS ONE. 2015; 10(6): e0129932 pmid:26114301
- 8. Scherer NM, Thompson CE, Freitas LB, Bonatto SL, Salzano FM. Patterns of molecular evolution in pathogenesis- related proteins. Genetics and Mol Bio. 2005; 4: 645–653.
- 9. Van Loon LC, Van Kammen A. Polyacrylamide disc electrophoresis of the soluble leaf proteins from Nicotiana tabaeum var. 'Samsun' and 'Samsun NN'. II. Changes in protein constitution after infection with tobacco mosaic virus. Virology. 1970; 40: 199–211.
- 10. Sudisha J, Sharathchandra RG, Amruthesh KN, Kumar A, Shetty HS. Plant Defence: Biological Control Progress in Biological Control. Mérillon JM. and Ramawat K.G. (eds.).2012; 12:378–403.
- 11. Christensen AB, Cho BH, Naesby M, Gregersen PL, Brandt T, Ordenna KM, et al. The molecular characterization of two barley proteins establishes the novel PR-17 family of pathogenesis related proteins. Mol Plant Pathol.2002; 3:135–144. pmid:20569319
- 12. Van Loon LC, Pierpoint WS, Boller Th, Conejero V. Recommendations for naming plant pathogenesis-related proteins. Plant Mol. Bio. Rep.1994; 12: 245–264.
- 13. Dafoe NJ, Gowen BE, Constabel P. Thaumatin-like proteins are differentially expressed and localized in phloem tissues of hybrid poplar. BMC Plant Bio. 2010; 10:191.
- 14. Liu J, Ekramoddoullah AKM. The family 10 of plant pathogenesis-related proteins: Their structure regulation and function in response to biotic and abiotic stresses. Physiol and Mol Plant Pathology 2006; 68: 3–13.
- 15. Vriens K, Cammue BPA, Thevissen K. Antifungal Plant Defensins: Mechanisms of Action and Production. Molecules. 2014; 19: 12280–12303. pmid:25153857
- 16. Caliskan M. Germin an Oxalate Oxidase has a function in many aspects of plant life. Turk J of Bio. 2000; 24: 717–724.
- 17. Lodhi N, Ranjan A, Singh M, Srivastava R, Singh SP, Chaturvedi CP, et al. Interactions between upstream and core promoter sequences determine gene expression and nucleosome positioning in tobacco PR-1a promoter. Biochimica et Biophysica Acta. 2008; 1779: 634–644. pmid:18723134
- 18. Kaur G,Pati PK.Analysis of cis-acting regulatory elements of Respiratory burst oxidase homolog (Rboh) gene families in Arabidopsis and rice provides clues for their diverse functions. Comput Bio and Chem. 2016; 62:104–18.
- 19. Wittkopp PJ, Kalay G. Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence. Nature Reviews Genetics. 2012; 13:59–69.
- 20. Higo K, Ugawa Y, Iwamoto M, Korenaga T. Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucl Acids Res. 1999; 27(1): 297–300. pmid:9847208
- 21. Lescot M, Dhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. PlantCARE a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002; 30 (1): 325–327. pmid:11752327
- 22. Davuluri RV, Sun H, Palaniswamy SK, Matthews N, Molina C, Kurtz M, Grotewold E. AGRIS: Arabidopsis Gene Regulatory Information Server an information resource of Arabidopsis cis-regulatory elements and transcription factors. BMC Bioinfo. 2003; 4:25.
- 23. Matys V, Fricke E, Geffers R, Gössling E, Haubrock M, Hehl R, et al.TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003; 31:374–8. pmid:12520026
- 24. Chow CN, Zheng HQ, Wu NY, Chien CH, Huang HD, Tzong-Yi Lee TY, et al. PlantPAN 2.0: an update of plant promoter analysis navigator for reconstructing transcriptional regulatory networks in plants. Nucl Acids Res. 2015: gkv1035v1-gkv1035.
- 25. Campanella JJ, Bitincka L, Smalley J. MatGAT: An application that generates similarity/identity matrices using protein or DNA sequences. BMC Bioinformatics.2003; 4:29–32. pmid:12854978
- 26. Marchler-Bauer A, Derbyshire DK, Gonzales NR, Lu S, Chitsaz F, Geer LY, Geer RC et al. CDD: NCBI’s conserved domain database. Nucleic Acids Research.2015; 43: 222–226.
- 27. Hu B, Jin J, Guo A, Zhang H, Luo J, Gao G. GSDS 2.0: an upgraded gene feature visualization server.Bioinfo.2015; 31(8):1296–1297.
- 28. Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W. Genevestigator. Arabidopsis Microarray Database and Analysis Toolbox. Plant Physiol. 2004; 136: 2621–2632. pmid:15375207
- 29. Jain S, Kumar A. The Pathogenesis Related Class 10 proteins in Plant Defense against Biotic and Abiotic Stresses. Adv Plants Agric Res.2015; 3(1): 00077.
- 30.
Patthy L. Introns: Phase Compatibility. eLS.2008.
- 31. Krishnamurthy A, Rathinasabapathi B. Oxidative stress tolerance in plants: Novel interplay between auxin and reactive oxygen species signaling. Plant Signaling & Behavior.2013; 8:e25761–5.
- 32. Sharma P, Jha AB, Dubey RS, Pessarakli M. Reactive oxygen species oxidative damage and antioxidative defense mechanism in plants under stressful conditions. J of Bot. 2012; 1: 27.
- 33. Petrov V, Vermeirssen V, De Clercq I, Breusegem FV, Minkov I, Vandepoele K, Gechev TS. Identification of cis-regulatory elements specific for different types of reactive oxygen species in Arabidopsis thaliana. Gene. 2012; 499: 52–60. pmid:22402413
- 34. Komarnytsky S, Borisjuk N. Functional analysis of promoters’ elements in plant. Genetic eng. 2003; 25:113–141.
- 35. Daudi A, Cheng Z, O'Brien JA, Mammarella N, Khan S, Ausubel FM, Bolwell GP. The apoplastic oxidative burst peroxidase in Arabidopsis is a major component of pattern-triggered immunity. Plant Cell. 2012; 24(1):275–87. pmid:22247251
- 36. Lacerda AF, Vasconcelos ÉAR, Pelegrini PB, Grossi de Sa MF. Antifungal defensins and their role in plant defense. Front. Microbiol. 2014; 5:116. pmid:24765086
- 37. Siberil Y, Doireau P, Gantet P. Plant bZIP G-box binding factors Modular structure and activation mechanisms. Eur J of Biochem. 2001; 268(22):5655–66.
- 38. Yoshihara T, Washida H, Takaiwa F.A 45-bp proximal region containing AACA and GCN4 motif is sufficient to confer endosperm-specific expression of the rice storage protein glutelin gene GluA-3.FEBS Letters. 1996: 383.
- 39. Choi C, Hwang SH, Fang IR, Kwon SI, Park SR, Ahn I, et al. Molecular characterization of Oryza sativa WRKY6 which binds to W-box-like element 1 of the Oryza sativa pathogenesis-related (PR) 10a promoter and confers reduced susceptibility to pathogens. New Phytol.2015; 208: 846–859. pmid:26083148
- 40. Rinerson CI, Scully ED, Palmer NA, Donze-Reiner T, Rabara RC, Tripathi P, et al. The WRKY transcription factor family and senescence in switchgrass. BMC Genomics.2015; 16: 912. pmid:26552372
- 41. Gadjev I, Vanderauwera S, Gechev TS, Laloi C, Minkov IN, Shulaev V, et al. Transcriptomic Footprints Disclose Specificity of Reactive Oxygen Species Signaling in Arabidopsis. Plant Physiol. 2006; 141(2): 436–445. pmid:16603662
- 42. Dolferus R, Klok EJ, Ismond K, Delessert C, Wilson S, Good A, et al. Molecular Basis of the Anaerobic Response in Plants. IUBMB Life.2000; 51(2):79–82.
- 43. Brown RL, Kazan K, McGrath KC, Maclean DJ, Manners JM. A Role for the GCC-Box in Jasmonate-Mediated Activation of the PDF1.2 Gene of Arabidopsis. Plant Physiol.2003; 132: 1020–1032. pmid:12805630
- 44. Grimmig B, Gonzalez-Perez MN, Leubner-Metzger G, Vogeli-Lange R, Meins F, Hain R, et al. Ozone-induced gene expression occurs via ethylene-dependentand -independent signalling. Plant Mol Biol. 2003; 51: 597–605.
- 45. Yu LM, Lamb CJ, Dixon RA. Purification and biochemical characterization of proteins which bind to the H-box cis-element implicated in transcriptional activation of plant defense genes. Plant J. 1993; 3(6):805–16. pmid:8401613
- 46. Wang P, Du Y, Zhao X, Miao Y, Song CP. The MPK6-ERF6-ROS-Responsive cis-Acting Element7/GCC Box Complex Modulates Oxidative Gene Transcription and the Oxidative Response in Arabidopsis. Plant Physiol.2013; 161(3):1392–408. pmid:23300166
- 47. Arias JA, Dixon RA, Lamb CJ. Dissection of the functional architecture of a plant defense gene promoter using a homologous in vitro transcription initiation system. Plant Cell.1993; 5: 485–496. pmid:8485404
- 48. Ezcurra I, Wycliffe P, Nehlin L, Ellerström M, Rask L. Transactivation of the Brassica napus napin promoter by ABI3 requires interaction of the conserved B2 and B3 domains of ABI3 with different cis-elements: B2 mediates activation through an ABRE whereas B3 interacts with an RY/G-box. The Plant J. 2000; 24: 57–66. pmid:11029704
- 49. Rushton PJ, Somssich IE, Ringler P, Shen QJ. WRKY transcription factors. Trends in Plant Science 2010; 15: 247–258. pmid:20304701
- 50. Le Gourrierec J, Li YF,Zhou DX. Transcriptional activation by Arabidopsis GT-1 may be through interaction with TFIIA—TBP—TATA complex. The Plant J. 1999; 18: 663–668. pmid:10417717
- 51. Chakravarthy S, Tuori RP, D'Ascenzo MD, Fobert PR, Despres C, Martin GB. The tomato transcription factor Pti4 regulates defense-related gene expression via GCC box and non-GCC box cis elements. Plant Cell. 2003; 15: 3033–3050. pmid:14630974
- 52. Giuliano G, Pichersky E, Malik VS, Timko MP, Scolnik PA, Cashmore AR. An evolutionarily conserved protein binding sequence upstream of a plant light-regulated gene. Proceedings of the National Academy of Sciences of the USA. 1988; 85(19): 7089–7093. pmid:2902624
- 53. Menkens AE, Schindler U, Cashmore AR. The G-box: a ubiquitous regulatory DNA element in plants bound by the GBF family of bZIP proteins. Trends in Biochem Sci.1995; 20: 506–510.
- 54. Liu J, Wang F, Yu G, Zhang X, Jia C, Qin J, Pan H. Functional Analysis of the Maize C-Repeat/DREMotif-Binding Transcription Factor CBF3Promoter in Response to Abiotic Stress. Int. J. Mol. Sci. 2015; 16:12131–12146. pmid:26030672
- 55. Nover L, Bharti K, Doring P, Mishra SK, Ganguli A, Scharf K. Arabidopsis and the heat stress transcription factor world: how many heat stress transcription factors do we need, Cell Stress & Chaperones.2001; 6 (3):177–189.
- 56. Hudson ME, Quail PH. Identification of Promoter Motifs Involved in the Network of Phytochrome A-Regulated Gene Expression by Combined Analysis of Genomic Sequence and Microarray Data. Plant Physiol. 2003; 133(4):1605–1616. pmid:14681527
- 57. Logemann E, Parniske M, Hahlbrock K. Modes of expression and common structural features of the complete phenylalanine ammonia-lyase gene family in parsley. Proc. Natl. Acad. Sci. U. S. A. 1995; 92: 5905–5909. pmid:7597051
- 58. Chan CS, Guo L, Shih MC. Promoter analysis of the nuclear gene encoding the chloroplast glyceraldehyde-3-phosphate dehydrogenase B subunit of Arabidopsis thaliana. Plant MolBiol.2001; 46: 131–141.
- 59. Chen W, Provart NJ, Glazebrook J, Katagiri F, Chang HS, et al. Expression profile matrix of Arabidopsis transcription factor genes suggests their putative functions in response to environmental stresses. Plant Cell. 2002; 14: 559–574. pmid:11910004
- 60. Abe H, Urao T, Ito T, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) Function as Transcriptional Activators in Abscisic Acid Signaling. The Plant Cell. 2003; 15(1): 63–78. pmid:12509522
- 61. Hattori T, Totsuka M, Hobo T, Kagaya Y, Yamamoto-Toyoda A. Experimentally determined sequence requirement of ACGT-containing abscisic acid response element. Plant Cell Physiol. 2002; 43:136–140. pmid:11828032
- 62. Kim SR, Kim Y, An G. Identification of methyl jasmonate and salicylic acid response elements from the nopaline synthase (nos) promoter. Plant Physiol.1993; 103: 97–103. pmid:8208860
- 63. Kim SY, Ma J, Perret P, Li Z, Thomas TL. Arabidopsis ABI5 Subfamily Members Have DistinctDNA-Binding and Transcriptional Activities. Plant Physiol.2002; 130: 688–697. pmid:12376636
- 64. Pozo MJ, Ent SVD, Van Loon LC, Pieterse CMJ. Transcription factor MYC2 involved in priming for enhanced defense during rhizobacteria-induced systemic resistance in Arabidopsis thaliana. New Phytologist. 2008; 180: 511–523. pmid:18657213
- 65. Ranty B, Aldon D, Galaud J. Plant Calmodulins and Calmodulin-Related Proteins. Plant Signaling & Behavior. 2006; 1: 96–104.
- 66. Reddy ASN, Ali GS, Celesnik H, Day IS. Coping with Stresses: Roles of Calcium- and Calcium/Calmodulin-Regulated Gene Expression. Plant Cell.2011; 23: 2010–2032. pmid:21642548
- 67. Whalley HJ, Sargeant AW, Steele JFC, Lacoere T, Lamb R, Saunders NJ, et al. Transcriptomic Analysis Reveals Calcium Regulation of Specific Promoter Motifs in Arabidopsis. The Plant Cell. 2011; 23: 4079–4095. pmid:22086087
- 68. Whalley HJ, Knight MR. Calcium signatures are decoded by plants to give specific gene responses. New Phytologist. 2013;197: 690–693. pmid:23190495
- 69. Huda KMK, Banu MSA, Pathi KM, Tuteja N. Reproductive Organ and Vascular Specific Promoter of the Rice Plasma Membrane Ca2+ATPase Mediates Environmental Stress Responses in Plants. PLoS ONE. 2013; 8(3): e57803. pmid:23469243
- 70. Zuo Y, Li Q. Identification of TATA and TATA-less promoters in plant genomes by integrating diversity measure, GC-Skew and DNA geometric flexibility. Genomics. 2011; 97: 112–120. pmid:21112384
- 71. Ishige F, Takaichi M, Foster R, Chua N, Oeda K. A G-box motif (GCCACGTGCC) tetramer confers high-level constitutive expression in dicot and monocot plants. The Plant Journal.1999; 18: 443–448.
- 72. Bell GI. Roles of repetitive sequences. Computers and Chemistry. 1992, 16: 135–143.
- 73. Mehrotra S, Goyal V. Repetitive sequences in plant nuclear DNA: Types Distribution Evolution and Function. Genomics Proteomics Bioinfo.2014; 12:164–171.
- 74. Liu F, Hu Z, Liu W, Li J, Wang W, Liang Z, et al. Distribution function and evolution characterization of microsatellite in Sargassum thunbergii (Fucales Phaeophyta) transcriptome and their application in marker development. Scientific reports. 2016; 6: 18947. pmid:26732855
- 75. Vinces DM, Legendre M, Caldara M, Hagihara M,Verstrepen KJ. Unstable Tandem Repeats in Promoters Confer Transcriptional Evolvability. Science. 2009; 324 (5931): 1213–1216. pmid:19478187
- 76. Barrera V, Peinado MA. Evaluation of single CpG sites are proxies of CpG island methylation states at the genome scale. Nucleic Acids Res. 2012; 40(22):11490–8. pmid:23066096
- 77. Wang P, Xia H, Zhang Y, Zhao S, Zhao C, et al. Genome-wide high-resolution mapping of DNA methylation identifies epigenetic variation across embryo and endosperm in Maize (Zea may). BMC Genomics.2015; 16:21. pmid:25612809
- 78. Rodrigues MGF, Soares MR, Gimenez DFJ, Fonseca LFS, et al.Detection of differentially methylated regions of irradiated fig tree selections. Sci. Agric. 2017; 74:285–293.
- 79. Vinson C, Chatterjee R. CG methylation. Epigenomics.2012; 6: 655–663.
- 80. Ashikawa I. Gene-associated CpG islands in plants as revealed by analyses of genomic sequences. The Plant J. 2001; 26(6): 617–625. pmid:11489175
- 81. Lim DHK, Maher ER. DNA methylation: a form of epigenetic control of gene expression. The Obstetrician & Gynaecologist. 2010; 12:37–42.
- 82. Grant-Downton RT, Dickinson HG. Epigenetics and its Implications for Plant Biology 2. The 'Epigenetic Epiphany': Epigenetics, Evolution and Beyond. Annals of Botany, 2006; 97: 11–27. pmid:16260442
- 83. Rambani A, Rice JH, Liu J, Lane T, et al. The Methylome of Soybean Roots during the Compatible Interaction with the Soybean Cyst Nematode. Plant Physiology.2015; 168:1364–1377. pmid:26099268
- 84. Zhang M, Kimatu JN, Xu K, Liu B. DNA cytosine methylation in plant development. J. Genet. Genomics. 2010; 37:1–12. pmid:20171573
- 85. Esteller M. CpG island hypermethylation and tumor suppressor genes: a booming present, a brighter future. Oncogene.2002; 21:5427–5440. pmid:12154405
- 86. Ferguson AA, Jiang N. Pack-MULEs-Recycling and reshaping genes through GC-biased acquisition. Mobile Genetic Elements. 2011; 1(2): 135–138. pmid:22016862
- 87. Singh R, Ming R, Yu Q. Nucleotide composition of the Nelumbo nucifera genome. Tropical Plant Bio. 2013; 6: 85–97.
- 88. Mongkolsiriwatana C, Pongtongkam P, Peyachoknagul S. In silico Promoter Analysis of Photoperiod-Responsive Genes Identified by DNA Microarray in Rice (Oryza sativa L.). Kasetsart J (Nat. Sci.). 2009; 43: 164–177.
- 89. Maruyama K, Todaka D, Mizoi J, Yoshida T, Kidokoro S, Matsukura S, et al. Identification of Cis-Acting Promoter Elements in Cold- and Dehydration- Induced Transcriptional Pathways in Arabidopsis Rice and Soybean. DNA Res. 2012; 19: 37–49. pmid:22184637
- 90. Nejad ES, Hossein AH, Soltani S. Regulatory TGACG-motif may elicit the secondary metabolite production through inhibition of active Cyclin-dependent kinase/Cyclin complex. Plant Omics J.2012; 5(6):553–558.