Skip to main content
Download PDF

Four distinct types of dehydration stress memory genes in Arabidopsis thaliana

BMC Plant Biology volume 13, Article number: 229 (2013) Cite this article

Abstract

Background

How plants respond to dehydration stress has been extensively researched. However, how plants respond to multiple consecutive stresses is virtually unknown. Pre-exposure to various abiotic stresses (including dehydration) may alter plants’ subsequent responses by improving resistance to future exposures. These observations have led to the concept of ‘stress memory’ implying that during subsequent exposures plants provide responses that are different from those during their first encounter with the stress. Genes that provide altered responses in a subsequent stress define the ‘memory genes’ category; genes responding similarly to each stress form the ‘non-memory’ category.

Results

Using a genome-wide RNA-Seq approach we determine the transcriptional responses of Arabidopsis plants that have experienced multiple exposures to dehydration stress and compare them with the transcriptional behavior of plants encountering the stress for the first time. The major contribution of this study is the revealed existence of four distinct, previously unknown, transcription memory response patterns of dehydration stress genes in A.thaliana. The biological relevance for each of the four memory types is considered in the context of four overlapping strategies employed by a plant to improve its stress tolerance and/or survival: 1) increased synthesis of protective, damage-repairing, and detoxifying functions; 2) coordinating photosynthesis and growth under repetitive stress; 3) re-adjusting osmotic and ionic equilibrium to maintain homeostasis; and 4) re-adjusting interactions between dehydration and other stress/hormone regulated pathways.

Conclusions

The results reveal the unknown, hitherto, existence of four distinct transcription memory response types in a plant and provide genome-wide characterization of memory and non-memory dehydration stress response genes in A.thaliana. The transcriptional responses during repeated exposures to stress are different from known responses occurring during a single exposure. GO analyses of encoded proteins suggested implications for the cellular/organismal protective, adaptive, and survival functions encoded by the memory genes. The results add a new dimension to our understanding of plants’ responses to dehydration stress and to current models for interactions between different signaling systems when adjusting to repeated spells of water deficits.

Background

Environmental stresses exert evolutionary pressure on organisms, which, in turn, have developed sophisticated responses to cope and to survive. Plants react to unfavorable conditions by dynamically changing physiological behavior and expression levels of implicated genes. Pre-exposure to various abiotic stresses (salinity, cold, high temperature) may alter plants’ subsequent responses by improving resistance to future exposures [1, 2]. Pre-treatment (priming) with hormones (jasmonic acid, JA, abscisic acid, ABA, salicylic acid, SA) increased the systemic immunity and induced stronger responses from responding genes upon subsequent treatments, relative to non-primed plants [3–7]. These observations have led to the concept of ‘stress memory’, implying that during subsequent exposures plants provide responses that are different from their responses during the first encounter with the stress [8].

Drought-triggered dehydration stress is one of the most common environmental stresses endured by plants and their responses to dehydration stress are extensively researched at organismal, cellular, and genome levels. Recently, we demonstrated that Arabidopsis plants subjected to several cycles of dehydration/water recovery treatments maintained higher relative water content than plants experiencing dehydration stress for the first time [9]. Dehydration stress ‘memory’ also affected gene expression from response genes. Analysis of 14 dehydration response genes after multiple stress treatments revealed the existence of two different types of transcriptional responses: genes that produced transcripts at similar levels during each stress and genes that in a subsequent response significantly increased their transcript levels and transcription rates relative to the first response [9]. Accordingly, genes that provide altered responses in a subsequent stress were referred to as ‘memory genes’ to distinguish them from ‘non-memory’ genes responding similarly to each stress.

Here, we explore how broadly memory genes are distributed within the dehydration stress responsive fraction of the Arabidopsis thaliana genome. Comprehensive RNA-Seq analyses of the transcriptomes of plants prior to stress, during a first stress, and during a third dehydration stress exposure revealed an unexpected diversity of memory-type responses at the transcriptional level. Analyses of the memory genes by GO functional categories indicated a biased functional distribution among the memory response patterns. Results are interpreted in the context of the possible biological relevance of these memory responses.

Methods

Plant growth and treatments

Arabidopsis thaliana (Col-0) plants were grown in potting soil in growth rooms at 22°C with a 12-h light photoperiod and light intensity of 180 μmol m-2 s-1. Repeated dehydration stresses were performed by air-drying for 2 h followed by a 22 h period of full re-hydration recovery as described [9]. RNA-Seq analyses were performed on rosette leaves from pre-stressed (W) plants, from plants exposed to the first dehydration treatment (S1), and during a third stress (S3) following two stress/recovery cycles. Plants from two independent biological samples were used. After removal from soil and before initiating the stress/recovery cycle, experimental plants were conditioned overnight in humid chambers and used to establish the basal transcript levels (W) for the transcriptome analyses. The transcriptional behavior of the marker memory genes RD29B and RAB18 were monitored in leaf samples taken from fresh plants in soil, from overnight conditioned plants after their removal from soil, and from plants kept under watered conditions throughout the stress/recovery as indicators of the normal dehydration stress behavior of plants in S1 and S2 and internal controls for the whole genome transcriptome analyses.

RNA extraction and RNA-Seq library construction

Leaf tissues were collected and immediately frozen in liquid nitrogen. Total RNA was extracted with Trizol (Invitrogen Inc. Carlsbad, CA, USA), treated with DNase I (Qiagen, Valencia, CA), and purified using Qiagen RNeasy Mini Kit. RNA integrity was confirmed on a Bioanalyzer 2100 using Nano 6000 LabChip (Agilent Technologies, Santa Clara, CA). Complementary DNA sequencing library was prepared from the total RNA using the mRNA-Seq Sample Preparation Kit (Illumina, San Diego, CA). Briefly, poly-adenylated RNA was isolated from 10 μg total RNA by Sera-Mag Magnetic Oligo-dT beads (ThermoFisher Scientific, Waltham, MA). RNA Purified mRNA was fragmented, annealed to high concentrations of random hexamers, and reverse transcribed. Following second strand cDNA synthesis, end repair, and A-tailing, Oligo adapters complementary to sequencing primers were ligated to cDNA fragment ends. Resultant cDNA libraries were size fractionated on an agarose gel, 200 bp fragments excised, and amplified by 15 cycles of polymerase chain reaction. Clusters were generated from the cDNA sequencing library on the surface of a flowcell in the Cluster Station (Illumina) by so-called bridge amplification. Replicates for the watered, S1 and S3 sample libraries were each run on a single lane in a flowcell on an Illumina GAIIx at the Genomics Core Facility at the University of Nebraska-Lincoln.

Reverse transcription and real-time PCR

Total RNA isolation and reverse transcription with oligo(dT) (18418–012, Invitrogen) were performed as described previously [10]. The amounts of individual genes were measured with gene-specific primers by real-time PCR analysis with a CyclerIQ real-time PCR Instrument (Bio-Rad) and SYBR Green mixture (Bio-Rad). The relative expression or amount of specific genes was quantitated with the 2-ΔΔCt calculation [11], according to the manufacturer’s software (Bio-Rad), where the reference gene was ubiquitin. Primers used in real-time RT-PCR are in Additional file 1: Table S1.

Bioinformatics analysis

Transcriptome sequencing of the watered, S1, and S3 samples yielded a total of 53.4, 76.8, and 78.4 million reads, respectively, summed over the two biological replicates per sample (see Additional file 2: Table S2). The read length for S1 and S3 is 75 bases, while for watered it is 101 bases. To determine the quality of the replicates we performed a least-square simple linear regression for each of the three samples. We calculated the R2 statistic (0.96 ≤ R2 ≤ 0.99) and slope (1.07 ≤ b ≤ 1.09), which provide measures of goodness-of-fit and correlation, respectively, using the regress function in MATLAB® (version 8.0.0.783 [R2012b]; The MathWorks™) (Additional file 3: Table S3). For use in all further analyses, the Arabidopsis thaliana genome and gene models were downloaded from Illumina’s iGenomes (Ensembl, TAIR10), and gene functional descriptions and Gene Ontology assignments from the Arabidopsis Information Resource (release TAIR10, [12]).

The bowtie (version 2.1.0; [13]) and tophat (version 2.0.8; [14]) packages were used with default parameters to map the RNA sequence reads from watered, S1, and S3 to the genome (see Additional file 3: Table S3), and to determine the expression quantity of known transcripts in each sample. The cuffdiff tool from the cufflinks package (version 2.0.2; [15]) was run with default parameters to calculate expression changes and associated q-values (False Discovery Rate adjusted p-values) for each gene, between the samples S1 and water, and S3 and S1. We further classify genes as being significantly differentially expressed when all three of the following conditions are met: q ≤ 0.05; | log2(fold change) | ≥ 1; and the FPKM-normalized expression value of at least one sample out of the two needs to be larger than the 25th percentile. The output files of cuffdiff are further annotated (in-house Perl script) by adding gene functional descriptions and GO classifications, and merged into a master file containing all data for S1 versus water and S3 versus S1 (Additional file 3: Table S3, Additional file 4: Table S4, and Additional file 5: Table S5).

From that master file we determined the 6579 significant drought-responsive genes (S1 versus water), and using that initial set we then looked at the significant responses in S3 versus S1 (1963 genes) and the remainder that did not respond (4616 genes). We assigned simple classifications to the types of response during the first stress (+ or -) and the second stress (+, -, or =), combining them into six classes: [++], [--], [+-], [-+], [+=], and [-=] (Table 1; Table 2; Additional file 4: Table S4, Additional file 5: Table S5). Two additional classes [=/+] and [=/-] contain genes that in S1 did not change significantly expression (according to the three criteria for significance, above) compared to pre-stressed levels in W but significantly changed transcription in S2. Formally, these genes do not belong to the S1 dehydration-stress responding fraction.

Table 1 Dehydrations response and memory genes in Arabidopsis thaliana
Full size table
Table 2 Distribution of dehydration stress responding genes according to GO function
Full size table

The raw transcriptome sequence files for watered, S1, and S3 have been uploaded, together with gene expression result files, to NCBI’s Gene Expression Omnibus under sequence number GSE48235.

Results

Dehydration stress memory response genes of Arabidopsis thaliana

Watered non-stressed plants and plants subjected to one or three dehydration stresses by exposure to dry air for 2 hours followed by 22 h periods of watered recovery intervals [9] were analyzed. The levels of mRNAs from the leaves of these plants were determined by RNA-Seq. Transcripts from 33,555 A. thaliana genes were identified (Table 1; Additional file 3: Table S3). Of these, 6579 genes (~20% of the genes in the genome) significantly increased or decreased their transcript levels during the first dehydration stress (S1), compared to their transcript levels in watered plants (W), representing the dehydration stress-responsive genes of A. thaliana. Of the genes responding in S1, 3396 genes were up-regulated and 3183 were down-regulated (Table 1; Additional file 4: Table S4). A comparison of the transcript levels during the first (S1) and third dehydration stress (S3) identified 1,963 dehydration responsive genes displaying significantly different amounts of transcripts in S3 versus S1. Our operational criterion for transcriptional memory is that the transcriptional responses to similar stress conditions must be different. Accordingly, these 1,963 genes display dehydration stress transcriptional memory (Table 1; Additional file 4: Table S4).

Four distinct transcriptional memory response-patterns were recognized within the S1-responding fraction. We designate these transcription patterns as [+/+], [-/-], [+/-] and [-/+] memory responses, where the first sign indicates higher (+) or lower (-) transcript levels in S1 relative to the levels in pre-stressed watered plants (W). The second sign indicates higher (+) or lower (-) transcripts in S3 relative to the levels in S1. Non-memory genes are denoted as [+/=] or [-/=] indicating that up-regulated or down-regulated transcript levels occurring in S1, respectively, have similar levels in S3. In addition, 798 genes were up-regulated in S2 ([=/+] W = S1 < S3) and 573 were down-regulated ([=/-] W = S1 > S3). These delayed response patterns are different from the memory responses represented by the four S1-responding memory types. Here, we focus mainly on the transcriptional patterns of the S1-responding gene fraction of Arabidopsis.

Of the four different memory classes, 362 genes denoted as [+/+] memory genes were induced in S1 and induced to higher levels in S3 (Table 1; Figure 1A). All 7 memory genes described in our earlier study [9] were present in the [+/+] memory category in the genome-wide RNA-Seq data (Additional file 4: Table S4). The [-/-] class contained 310 genes with decreased transcript levels in S1 and further reduced levels in S3 (Table 1; Figure 1B). The [-/+] class contained 434 memory genes that were down-regulated in S1, but produced significantly higher transcript levels in S3 (Figure 1C); conversely, 857 [+/-] memory genes, were induced in S1 but had lower levels in S3 (Figure 1D). The latter two categories ‘revise’ their transcriptional behavior in a subsequent stress: after robustly responding in S1, these genes show weaker/no responses in S3, producing transcripts at levels close to their initial watered (pre-stressed) levels. These genes are referred to as ‘revised response’ memory genes. The distribution of transcript levels for the dehydration stress-responsive memory genes in S1 and S3 illustrating the four memory categories is shown (Figure 2A, B). The two non-memory gene classes ([+/=] and [-/=]) consistently provide transcriptional responses in S3 similar to those in S1 (Figure 1E, F).

Figure 1
figure 1

Transcription memory response patterns of dehydration stress responding genes. Three genes from each memory category were randomly chosen from the RNA-Seq datasets to illustrate the diverse transcriptional response patterns. A) [+/+] memory genes; B) [-/-] memory genes; C) [-/+] memory genes; D) [+/-] memory genes; E) [+/=] non-memory genes; F) [-/=] non-memory genes. Transcript levels measured by real-time qRT-PCR are under initial pre-stressed watered (W) conditions, after the first (S1), and the third (S3) exposures to dehydration stress. Measurements were performed on three biological samples, each with three technical replicates. Bars are the mean ± s.e.m., n = 3 replicates. For interpretation of the response signs see text.

Full size image
Figure 2
figure 2

Distribution of dehydration stress responding genes in A. thaliana in S1 and S3. Left panel). Transcript levels from dehydration stress responding genes that are up-regulated or down-regulated during S1 (color key at the bottom) plotted by the log2 of their S1 levels along the x-axis, and by log2 of their S1/watered ratio along the y-axis. Right panel) Transcript levels of dehydration stress memory genes: as in (A) except the y-axis is the log2 of the S3/watered ratio. The clustering of the four colors illustrates the distribution of the four distinct memory response types: revised response [+/-] and [-/+] memory genes clustering closer to their pre-stressed (W) levels, while the [+/+] and [-/-] increasing separation from these levels.

Full size image

Functional distribution of Arabidopsis memory genes

The large number of the dehydration stress-memory responsive genes and the diversity of their transcriptional responses raised questions of their functions and of whether there was a preferential association of particular cellular functions with any of the four memory types. To address these questions, genes displaying transcriptional memory were analyzed for biological function according to their Gene Ontology (GO) classification (Table 2; Additional file 5: Table S5).

Genes involved in ABA/abiotic stresses responses are signature for the [+/+] memory genes

Genes implicated in responses to salt, salinity, cold/heat acclimation, and abscisic acid (ABA) constitute about a quarter of the [+/+] memory genes. LEA genes represent 3% among them (Table 2; Additional file 5: Table S5; Figure 3).

Figure 3
figure 3

Distribution of dehydration stress memory genes in A. thaliana according to GO functions. The functional GO categories represented by the more abundant (5% or more of the genes within a memory response type) GO groups of responding genes. The four memory and two non-memory types are shown in a Venn diagram format: A) the two memory and one non-memory types up-regulated in S1; B) the two memory and one non-memory types down-regulated in S1.

Full size image

Ribosome/protein synthesis, chloroplast, and thylakoid membrane associated functions are signature for [-/-] memory genes

Memory genes encoding ribosomal, as well as chloroplast and photosynthesis proteins are the major constituents of the [-/-] memory subgroup (Additional file 4: Table S4; Additional file 5: Table S5; Figure 3). Encoded proteins are involved in ribosome structure and amino acid biosynthesis, as well as in photosynthesis, in the light harvesting complex (photosystem II), and in responses to light (red/far red light, in particular). The largest number of thylakoid membrane-associated proteins is encoded by [-/-] memory genes.

Functions encoded by [-/+] revised response memory genes

There are no particular types of cellular functions specifically enriched within the [-/+] group of memory genes, although the largest percentages of genes encode chloroplast and membrane (plasma, organellar, and thylakoid membrane) related functions. Many of these are for electron transport, photosystem II assembly, chloroplast organization and re-location, functions similar to those encoded by [-/-] memory genes (Table 2; Additional file 5: Table S5; Figure 3).

Dehydration stress [+/-] memory genes are shared with multiple signaling pathways

The signature [+/-] dehydration memory genes are regulated by multiple signaling pathways, including the ABA, ethylene, auxin (IAA), gibberrellic acid (GA), jasmonic acid (JA), and salicylic acid (SA) pathways (Table 2; Additional file 5: Table S5; Figure 3). A prominent functional group represented by [+/-] memory genes encodes proteins associated with membranes. Although membrane (plasma and organellar)-related genes are highly represented within the entire drought-responding fraction (including the non-memory and all four memory types), genes for tonoplast intrinsic proteins regulating water transport (GAMMA-TIP, PIP2B, TMP-C, PIP2A, TIP2, PIP1A, RD28) and the inward K+ channel proteins (KAT1 and KAT2) display, exclusively, [+/-] memory responses (Table 2; Additional file 5: Table S5).

Dehydration stress memory of transcription factor genes

Among the ~1500 Arabidopsis genes encoding transcription factors (TFs) [16], members of about two dozen families have been implicated in responses to drought [17–19]. Among these families, transcriptional memory behavior was found for members of ten of the families (Additional file 6: Table S6). Five families (AP2/ERF, bHLH, Homeodomain-Zip, MYB/Myb-like, and ZF) are represented in all memory categories; GRAS and HSF members display only [+/+] and [+/-] memory, and the CAAT family is represented only by [+/+] memory genes. TF genes are least represented in the [-/-] memory subgroup (Table 2; Additional file 6: Table S6).

[+/+] memory genes for MYB2, MYB112, MYB13, MYB47, ATHB-7, ATHB-14, and ABF2 specifically cluster with ABA/abiotic stress responding genes but are rarely found among other hormonally responsive pathways (Figure 4A); six putative TF genes display [-/-] memory type responses and only three are shared with other response pathways (Figure 4B). TFs with [+/-] or [-/+] revised response memory behavior are broadly shared with multiple abiotic and hormone response networks, outlining distinct clusters (Figure 4C, D).

Figure 4
figure 4

Dehydration stress memory genes encoding TFs. Hierarchical clustering of dehydration stress memory genes encoding transcription factors implicated in the crosstalk between various abiotic stress and hormone regulated pathways. Presence (red, value 1) or non-presence (black, value 0) of annotation in stress and hormone regulated pathways is determined for each transcription factor, and then clustered along both rows (genes) and columns (annotation) via unweighted average Euclidean distance. A) TF genes displaying [+/+] memory; B) [-/-]; C [+/-]; D [-/+].

Full size image

Discussion

Biological relevance of the transcription memory behavior

Plants experience more dehydration stress during the day when transpiration exceeds the rate of water uptake and fully, or partially, recover during the night when the rate of water uptake exceeds the greatly diminished loss of water from leaves [20]. This creates a diurnal cycle of oscillating leaf water potential that is probably most pronounced during periods of moderate drought conditions [21, 22]. The repetitive dehydration stress system employed in this study has similarities to this natural diurnal stress. We hypothesize that dehydration stress memory helps plants to prepare for the next day’s stress if they were stressed the day before, despite alleviated stress signals during the night [9].

Genes displaying transcriptional memory are among the genes responding to the first dehydration stress. Altering their expression levels in subsequent stresses, presumably, allows the plant to finely tune its responses to ongoing/recurring dehydration stress. The possible biological relevance of the genes displaying transcriptional memory is considered in the context of four overlapping strategies generally employed by a plant to improve its stress tolerance and/or survival: 1) increased synthesis of membrane protecting, damage-repairing, and detoxifying functions; 2) coordinating photosynthesis and growth under repetitive stress; 3) re-adjusting osmotic and ionic equilibrium to maintain homeostasis [23]; and 4) re-adjusting interactions between dehydration and other stress/hormone regulated pathways [24, 25].

Increased synthesis of protective functions

In general, the proteins encoded by [+/+] memory genes are involved in cell protective roles. These include heat shock proteins and chaperons, proteins involved in the repair of stress-induced damages to membranes, in modifying membrane lipid composition, in regulating membrane fluidity and permeability to toxic ions, dehydrins (RD29B, RAB18, LTi30), and lipid transfer proteins (LTP2, LTP3, LTP4) [26, 27]. LEA (late embryogenesis abundant) proteins functioning as molecular chaperones to maintain membrane structures, ion balance and homeostasis [28, 29] and enzymes for the synthesis of isoleucine (toxin degradation), serine (redox responses), and proline (an osmolyte) are encoded by [+/+] memory genes. Thereby, [+/+] memory genes ensure elevated synthesis of factors critical for cell survival under multiple dehydration stresses (Table 2; Figure 3).

Abscisic acid (ABA) is the key mediator of dehydration signaling and is involved in responses to other abiotic stress response pathways [30, 31]. ABA responsive genes that respond to dehydration are often points of overlap in various abiotic stress pathways [32, 33] (Table 2). The [+/+] memory and [+/=] non-memory genes, generally, encode proteins with similar protective functions. However, the memory component, presumably, allows the plant to optimize its stress tolerance within the context of protection. Only a few [-/-] memory genes and a smaller number of down-regulated [-/=] genes are shared between dehydration/ABA and other abiotic stress responding systems (Table 2; Additional file 5: Table S5; Figure 3).

Photosynthesis and growth during repeated stresses

More than 40% of the [-/-] memory genes encode ribosome/protein synthesis, chloroplast, and thylakoid membrane related proteins indicating repetitive stress results in a stronger attenuation of the pathways for photosystem II electron transport, small RuBisCo subunit and sugar biosynthesis, as well as for ribosome assembly, protein translation, DNA replication and histones (Additional file 4: Table S4; Additional file 5: Table S5; Figure 3). None of the [+/+] memory genes encodes a thylakoid membrane function. Of note, 12 [+/+] memory genes associated with chloroplasts encode functions that are different, even opposite to, the functions encoded by the [-/-] memory genes. Chlorophyll and toxin catabolic enzymes, enzymes for the synthesis of isoleucine (toxin degradation), serine (redox responses), and proline (an osmolyte) are encoded by [+/+] memory genes. Photosynthetic and starch synthesizing activities, in addition to proteins involved in chloroplast organization and relocation, are encoded by [-/-] memory genes. Down-regulated [-/=] non-memory genes encode ribosome- and chloroplast-associated proteins as well (Table 2; Additional file 5: Table S3; Figure 3). The functional implications of the [-/-] transcriptional response patterns (decreased protein synthesis, cell growth and photosynthesis) and the protective functions by the [+/+] memory genes are in agreement with known responses of plants enduring dehydration stress [34–37]. Of note, on the background of a rather even distribution of GO functions encoded by the late responding [=/+] and [=/-] genes, the fraction of genes encoding chloroplast and thylakoid membrane proteins constitute ~ 20% among the [=/-] genes (Table 2; Additional file 5: Table S5).

Re-adjusting cellular homeostasis

The [-/+] and [+/-] revised memory genes respond robustly to the initial stress, but in a subsequent exposure provide responses at levels closer to their non-stressed levels. Presumably, these genes help in restoring homeostasis as the plant adjusts to dehydration stress. Thus, [-/+] memory genes for chloroplast organization, relocation, photosynthetic, light harvesting, and metabolic functions are similar to the functions encoded by the [-/-] memory genes (Additional file 5: Table S5; Figure 3). However, by producing significantly more transcripts in S3 than in S1([-/+]), or returning to pre-stressed (W) levels in S3 ([+/-]), these genes support the ability of plants to restore photosynthesis after temporary water withdrawal [38]. Presumably, the [-/+] genes contribute to an initial decrease in chloroplast/photosynthetic activity in S1, but by reverting to pre-stressed transcription levels in S3, contribute to a subsequent re-adjustment. The [-/+] memory genes GLK1 and GLK2, implicated in regulating the photosynthetic apparatus and chloroplast development in a cell-autonomous manner [39, 40] are candidate genes for a role in such an adjustment.

Within the same GO category, [-/+] and [+/-] memory genes may encode opposing biochemical activities. For example, a number of the [-/+] memory genes encode chlorophyll and carbohydrate (starch) synthesizing enzymes, while chlorophyll and starch degrading activities are encoded by [+/-] memory genes (Additional file 5: Table S5). Juxtaposing these functions would suggest complementary roles during both the initial and repeated stresses. Acting together, these genes are likely to decrease overall chlorophyll and starch biosynthesis during the first stress but to partially restore it during subsequent stresses. These complementary transcriptional patterns and functional roles support the fine coordination occurring between metabolic and energy adjustments during adaptation to drought [41].

Membrane-associated memory genes present another paradigm of biological relevance for the dehydration memory genes. As dehydration may damage membrane integrity, it is not surprising that a large number of genes (more than 750 of memory and non-memory) encode functions related to membranes; more than 260 of the late responding genes encode membrane-associated proteins (Table 2; Additional file 5: Table S5; Figure 3). Proteins regulating osmotic pressure, water balance, and wall modifications have been implicated in plants’ stress responses and environmental adaptation [42]. These functions are encoded by genes of all four-memory response types, suggesting that after providing an up/down transcriptional response in S1, a large number of Arabidopsis genes alter their transcript outputs (further up- or down-regulating or reversing towards pre-stressed levels) in a subsequent exposure to more finely modify cell membrane structure, ion balance and homeostasis.

Dehydration stress memory genes in crosstalk with other response pathways

Individual genes that are regulated by multiple different stimuli represent overlapping points of these plant’s signaling networks [43–48]. The subset of the dehydration stress memory genes that are at these convergent points revealed new aspects of plants’ stress responses. For example the dehydration/JA pathways share the largest number of both memory and non-memory genes co-regulated by these pathways, consistent with reported crosstalk between jasmonic acid (JA) and ABA regulated signaling networks during dehydration [45, 47, 49], salt [50], and cold stresses [51]. More than 200 JA-responsive genes are up-regulated in S1 (including memory and non-memory genes) and about 80 are down-regulated (Table 2; Additional file 5: Table S5). In S3 stress, however, 89 [+/=] and 46 [-/=] genes continued to provide similar transcript levels, seven produced higher ([+/+]), and 16 lower ([-/-]), transcript levels, but 140 ([+/-]) genes did not respond in the subsequent stress. Thereby, in addition to a finer tuning, the altered transcript levels from the memory genes suggest an altered crosstalk between these pathways in S3 compared to S1.

Non-memory ([+/=] and [-/=]) and unidirectional memory ([+/+] and [-/-]) dehydration stress responding genes respond every time the plant experiences dehydration stress. Concurrent with amplified responses from the [+/+] and [-/-] memory genes, however, engagement of the [+/-] and [-/+] genes lessens in S3, suggesting the nature and the dynamics of the interactions between dehydration and other signaling pathways are different during repeated exposures than the interactions occurring in a single exposure.

Lastly, we note that TF encoded by memory genes can be critical for the expression of dependent genes in S1, but that they do not necessarily determine the memory behavior of regulated genes in S3. For example, the [+/-] memory of the MYC2 gene, identified as a master regulator of the crosstalk between the ABA, SA, GA, JA, and auxin signaling pathways [52–54], correlates with the [+/-] memory of a large number of MYC2-dependent genes common for these signaling pathways (Figure 4C; Additional file 6: Table S6). However, the signature gene, RD22, directly regulated by MYC2 [55, 56] is a non-memory gene that is up-regulated in both S1 and in S3 (Additional file 4: Table S4; Additional file 6: Table S6) despite the fairly low levels of MYC2 transcripts in S3. Evidently, another TF activates RD22 in S3. The results suggest that diverse gene-specific mechanisms are involved in regulating the behavior of dehydration stress responding genes and that the memory behavior of individual TFs alone does not necessarily determine or predict the memory or non-memory behavior of their targets.

Conclusion

The genome-wide response of Arabidopsis genes to dehydration stress revealed the existence of four distinct transcriptional memory response patterns. By altering transcript levels, and presumably the levels of encoded proteins, memory genes are likely to alter the cellular responses and the crosstalk between overlapping pathways. Adjustments of expression of memory genes, together with the consistent responses from the non-memory genes, allow the plant to optimize its responses and the interactions between various signaling systems. Transcriptional memory, like defense gene priming, can provide the benefits of a more robust or modified stress response while reducing the costs of the state of preparedness [25]. The behavior of transcriptional memory genes adds a new dimension to our understanding of plants’ responses to dehydration stress and to current models for interactions between different signaling systems. Revealing the molecular mechanisms of transcriptional memory responses may be critical for understanding how plants’ adapt to changing environments and is emerging as a new area in plant abiotic and biotic stress response research.

References

  1. Goh CH, Nam HG, Park YS: Stress memory in plants: a negative regulation of stomatal response and transient induction of rd22 gene to light in abscisic acid-entrained Arabidopsis plants. Plant J. 2003, 36: 240-255. 10.1046/j.1365-313X.2003.01872.x.

    Article  PubMed  CAS  Google Scholar 

  2. Sung DY, Kaplan F, Lee KJ, Guy CL: Acquired tolerance to temperature extremes. Trends Plant Sci. 2003, 8: 179-187. 10.1016/S1360-1385(03)00047-5.

    Article  PubMed  CAS  Google Scholar 

  3. Baldwin IT, Schmelz E: Immunological “memory” in the induced accumulation of nicotine in wild tobacco. Ecology. 1996, 77: 236-246. 10.2307/2265673.

    Article  Google Scholar 

  4. Conrath U, Beckers GJ, Flors V, Garcia-Agustin P, Jakab G, Mauch F, Newman MA, Pieterse CM, Poinssot B, Pozo MJ, et al: Priming: getting ready for battle. Mol Plant Microbe Interact. 2006, 19: 1062-1071. 10.1094/MPMI-19-1062.

    Article  PubMed  CAS  Google Scholar 

  5. Ton J, D’Alessandro M, Jourdie V, Jakab G, Karlen D, Held M, Mauch-Mani B, Turlings TC: Priming by airborne signals boosts direct and indirect resistance in maize. Plant J. 2007, 49: 16-26.

    Article  PubMed  CAS  Google Scholar 

  6. Ton J, Jakab G, Toquin V, Flors V, Iavicoli A, Maeder MN, Metraux JP, Mauch-Mani B: Dissecting the beta-aminobutyric acid-induced priming phenomenon in Arabidopsis. Plant Cell. 2005, 17: 987-999. 10.1105/tpc.104.029728.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Jaskiewicz M, Conrath U, Peterhansel C: Chromatin modification acts as a memory for systemic acquired resistance in the plant stress response. EMBO Rep. 2011, 12: 50-55. 10.1038/embor.2010.186.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  8. Bruce TJ, Matthes MC, Napier JA, Pickett JA: Stressful “memories” of plants: evidence for possible mechanisms. Plant Sci. 2007, 173: 603-608. 10.1016/j.plantsci.2007.09.002.

    Article  CAS  Google Scholar 

  9. Ding Y, Fromm M, Avramova Z: Multiple exposures to drought ‘train’ transcriptional responses in Arabidopsis. Nat Commun. 2012, 3: 740.

    Article  PubMed  Google Scholar 

  10. Ding Y, Avramova Z, Fromm M: Two distinct roles of ARABIDOPSIS HOMOLOG OF TRITHORAX1 (ATX1) at promoters and within transcribed regions of ATX1-regulated genes. Plant Cell. 2011, 23: 350-363. 10.1105/tpc.110.080150.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  PubMed  CAS  Google Scholar 

  12. Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M, et al: The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012, 40: D1202-D1210. 10.1093/nar/gkr1090.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9: 357-359. 10.1038/nmeth.1923.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14: R36. 10.1186/gb-2013-14-4-r36.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L: Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013, 31: 46-53.

    Article  PubMed  CAS  Google Scholar 

  16. Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, et al: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290: 2105-2110. 10.1126/science.290.5499.2105.

    Article  PubMed  CAS  Google Scholar 

  17. Shinozaki K, Yamaguchi-Shinozaki K, Seki M: Regulatory network of gene expression in the drought and cold stress responses. Curr Opin Plant Biol. 2003, 6: 410-417. 10.1016/S1369-5266(03)00092-X.

    Article  PubMed  CAS  Google Scholar 

  18. Zhang JZ, Creelman RA, Zhu JK: From laboratory to field. Using information from Arabidopsis to engineer salt, cold, and drought tolerance in crops. Plant Physiol. 2004, 135: 615-621. 10.1104/pp.104.040295.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Fujita Y, Fujita M, Shinozaki K, Yamaguchi-Shinozaki K: ABA-mediated transcriptional regulation in response to osmotic stress in plants. J Plant Res. 2011, 124: 509-525. 10.1007/s10265-011-0412-3.

    Article  PubMed  CAS  Google Scholar 

  20. Ruggiero C, Angelino G, Maggio A: Developmental regulation of water uptake in wheat. J of plant Physiol. 2007, 164: 1170-1178. 10.1016/j.jplph.2006.06.017.

    Article  CAS  Google Scholar 

  21. Tardieu F, Parent B, Simonneau T: Control of leaf growth by abscisic acid: hydraulic or non-hydraulic processes?. Plant Cell Environ. 2010, 33: 636-647. 10.1111/j.1365-3040.2009.02091.x.

    Article  PubMed  Google Scholar 

  22. Mane SP, Vasquez-Robinet C, Sioson AA, Heath LS, Grene R: Early PLDalpha-mediated events in response to progressive drought stress in Arabidopsis: a transcriptome analysis. J Exp Bot. 2007, 58: 241-252.

    Article  PubMed  CAS  Google Scholar 

  23. Mahajan S, Tuteja N: Cold, salinity and drought stresses: an overview. Arch Biochem Biophys. 2005, 444: 139-158. 10.1016/j.abb.2005.10.018.

    Article  PubMed  CAS  Google Scholar 

  24. De Vos M, Van Oosten VR, Van Poecke RM, Van Pelt JA, Pozo MJ, Mueller MJ, Buchala AJ, Metraux JP, Van Loon LC, Dicke M, et al: Signal signature and transcriptome changes of Arabidopsis during pathogen and insect attack. Mol Plant Microbe Interact. 2005, 18: 923-937. 10.1094/MPMI-18-0923.

    Article  PubMed  CAS  Google Scholar 

  25. van Hulten M, Pelser M, van Loon LC, Pieterse CM, Ton J: Costs and benefits of priming for defense in Arabidopsis. Proc Natl Acad Sci USA. 2006, 103: 5602-5607. 10.1073/pnas.0510213103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Torres-Schumann S, Godoy JA, Pintor-Toro JA: A probable lipid transfer protein gene is induced by NaCl in stems of tomato plants. Plant Mol Biol. 1992, 18: 749-757. 10.1007/BF00020016.

    Article  PubMed  CAS  Google Scholar 

  27. Kader JC: Lipid-transfer proteins in plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 627-654. 10.1146/annurev.arplant.47.1.627.

    Article  PubMed  CAS  Google Scholar 

  28. Battaglia M, Olvera-Carrillo Y, Garciarrubio A, Campos F, Covarrubias AA: The enigmatic LEA proteins and other hydrophilins. Plant Physiol. 2008, 148: 6-24. 10.1104/pp.108.120725.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Huang D, Wu W, Abrams SR, Cutler AJ: The relationship of drought-related gene expression in Arabidopsis thaliana to hormonal and environmental factors. J Exp Bot. 2008, 59: 2991-3007. 10.1093/jxb/ern155.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Thomashow MF: PLANT COLD ACCLIMATION: Freezing tolerance genes and regulatory mechanisms. Annu Rev Plant Physiol Plant Mol Biol. 1999, 50: 571-599. 10.1146/annurev.arplant.50.1.571.

    Article  PubMed  CAS  Google Scholar 

  31. Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D’Angelo C, Bornberg-Bauer E, Kudla J, Harter K: The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J. 2007, 50: 347-363. 10.1111/j.1365-313X.2007.03052.x.

    Article  PubMed  CAS  Google Scholar 

  32. Kreps JA, Wu Y, Chang HS, Zhu T, Wang X, Harper JF: Transcriptome changes for Arabidopsis in response to salt, osmotic, and cold stress. Plant Physiol. 2002, 130: 2129-2141. 10.1104/pp.008532.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. Seki M, Ishida J, Narusaka M, Fujita M, Nanjo T, Umezawa T, Kamiya A, Nakajima M, Enju A, Sakurai T, et al: Monitoring the expression pattern of around 7,000 Arabidopsis genes under ABA treatments using a full-length cDNA microarray. Funct Integr Genomics. 2002, 2: 282-291. 10.1007/s10142-002-0070-6.

    Article  PubMed  CAS  Google Scholar 

  34. Muller B, Pantin F, Genard M, Turc O, Freixes S, Piques M, Gibon Y: Water deficits uncouple growth from photosynthesis, increase C content, and modify the relationships between C and growth in sink organs. J Exp Bot. 2011, 62: 1715-1729. 10.1093/jxb/erq438.

    Article  PubMed  CAS  Google Scholar 

  35. Pinheiro C, Chaves MM: Photosynthesis and drought: can we make metabolic connections from available data?. J Exp Bot. 2011, 62: 869-882. 10.1093/jxb/erq340.

    Article  PubMed  CAS  Google Scholar 

  36. Chaves MM, Flexas J, Pinheiro C: Photosynthesis under drought and salt stress: regulation mechanisms from whole plant to cell. Ann Bot. 2009, 103: 551-560.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Harb A, Krishnan A, Ambavaram MM, Pereira A: Molecular and physiological analysis of drought stress in Arabidopsis reveals early responses leading to acclimation in plant growth. Plant Physiol. 2010, 154: 1254-1271. 10.1104/pp.110.161752.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  38. Hayano-Kanashiro C, Calderon-Vazquez C, Ibarra-Laclette E, Herrera-Estrella L, Simpson J: Analysis of gene expression and physiological responses in three Mexican maize landraces under drought stress and recovery irrigation. PLoS One. 2009, 4: e7531. 10.1371/journal.pone.0007531.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Waters MT, Moylan EC, Langdale JA: GLK transcription factors regulate chloroplast development in a cell-autonomous manner. Plant J. 2008, 56: 432-444. 10.1111/j.1365-313X.2008.03616.x.

    Article  PubMed  CAS  Google Scholar 

  40. Waters MT, Wang P, Korkaric M, Capper RG, Saunders NJ, Langdale JA: GLK transcription factors coordinate expression of the photosynthetic apparatus in Arabidopsis. Plant Cell. 2009, 21: 1109-1128. 10.1105/tpc.108.065250.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  41. Sperdouli I, Moustakas M: Interaction of proline, sugars, and anthocyanins during photosynthetic acclimation of Arabidopsis thaliana to drought stress. J Plant Physiol. 2012, 169: 577-585. 10.1016/j.jplph.2011.12.015.

    Article  PubMed  CAS  Google Scholar 

  42. Moore JP, Vicre-Gibouin M, Farrant JM, Driouich A: Adaptations of higher plant cell walls to water loss: drought vs desiccation. Physiol Plant. 2008, 134: 237-245. 10.1111/j.1399-3054.2008.01134.x.

    Article  PubMed  CAS  Google Scholar 

  43. Ma S, Bohnert HJ: Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell-specific expression. Genome Biol. 2007, 8: R49. 10.1186/gb-2007-8-4-r49.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Swindell WR: The association among gene expression responses to nine abiotic stress treatments in Arabidopsis thaliana. Genetics. 2006, 174: 1811-1824. 10.1534/genetics.106.061374.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  45. Adie BA, Perez-Perez J, Perez-Perez MM, Godoy M, Sanchez-Serrano JJ, Schmelz EA, Solano R: ABA is an essential signal for plant resistance to pathogens affecting JA biosynthesis and the activation of defenses in Arabidopsis. Plant Cell. 2007, 19: 1665-1681. 10.1105/tpc.106.048041.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  46. Suhita D, Raghavendra AS, Kwak JM, Vavasseur A: Cytoplasmic alkalization precedes reactive oxygen species production during methyl jasmonate- and abscisic acid-induced stomatal closure. Plant Physiol. 2004, 134: 1536-1545. 10.1104/pp.103.032250.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  47. Jung C, Lyou SH, Yeu S, Kim MA, Rhee S, Kim M, Lee JS, Choi YD, Cheong JJ: Microarray-based screening of jasmonate-responsive genes in Arabidopsis thaliana. Plant Cell Rep. 2007, 26: 1053-1063. 10.1007/s00299-007-0311-1.

    Article  PubMed  CAS  Google Scholar 

  48. Anderson JP, Badruzsaufari E, Schenk PM, Manners JM, Desmond OJ, Ehlert C, Maclean DJ, Ebert PR, Kazan K: Antagonistic interaction between abscisic acid and jasmonate-ethylene signaling pathways modulates defense gene expression and disease resistance in Arabidopsis. Plant Cell. 2004, 16: 3460-3479. 10.1105/tpc.104.025833.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Munemasa S, Oda K, Watanabe-Sugimoto M, Nakamura Y, Shimoishi Y, Murata Y: The coronatine-insensitive 1 mutation reveals the hormonal signaling interaction between abscisic acid and methyl jasmonate in Arabidopsis guard cells. Specific impairment of ion channel activation and second messenger production. Plant Physiol. 2007, 143: 1398-1407. 10.1104/pp.106.091298.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  50. Walia H, Wilson C, Condamine P, Liu X, Close TJ: Role of jasmonic acid in improving salinity tolerance of barley revealed by physiological and large scale gene expression studies. Plant, Cell and Environ. 2007, 30: 410-421. 10.1111/j.1365-3040.2006.01628.x.

    Article  CAS  Google Scholar 

  51. Wong CE, Li Y, Labbe A, Guevara D, Nuin P, Whitty B, Diaz C, Golding GB, Gray GR, Weretilnyk EA, et al: Transcriptional profiling implicates novel interactions between abiotic stress and hormonal responses in Thellungiella, a close relative of Arabidopsis. Plant Physiol. 2006, 140: 1437-1450. 10.1104/pp.105.070508.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  52. Abe H, Yamaguchi-Shinozaki K, Urao T, Iwasaki T, Hosokawa D, Shinozaki K: Role of arabidopsis MYC and MYB homologs in drought- and abscisic acid-regulated gene expression. Plant Cell. 1997, 9: 1859-1868.

    PubMed  CAS  PubMed Central  Google Scholar 

  53. Dombrecht B, Xue GP, Sprague SJ, Kirkegaard JA, Ross JJ, Reid JB, Fitt GP, Sewelam N, Schenk PM, Manners JM, et al: MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell. 2007, 19: 2225-2245. 10.1105/tpc.106.048017.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  54. Kazan K, Manners JM: MYC2: the master in action. Mol Plant. 2013, 6: 686-703. 10.1093/mp/sss128.

    Article  PubMed  CAS  Google Scholar 

  55. Yamaguchi-Shinozaki K, Shinozaki K: The plant hormone abscisic acid mediates the drought-induced expression but not the seed-specific expression of rd22, a gene responsive to dehydration stress in Arabidopsis thaliana. Mol Gen Genet. 1993, 238: 17-25.

    PubMed  CAS  Google Scholar 

  56. Boter M, Ruiz-Rivero O, Abdeen A, Prat S: Conserved MYC transcription factors play a key role in jasmonate signaling both in tomato and Arabidopsis. Genes Dev. 2004, 18: 1577-1591. 10.1101/gad.297704.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors are grateful to Stephanie Canny for producing Figure 2. Partial support was provided by NSF award MCB-1121898 to Z.A., M.F., and J.-J.R., and by the Nebraska Research Initiative to J.-J.R.

Author information

Author notes

    Authors and Affiliations

    1. University of Science & Technology of China, 443 Huangshang Road, Hefei, Anhui, 230027, China

      Yong Ding

    2. University of Nebraska School of Biological Sciences, 1901 Vine Street, Lincoln, 68588, USA

      Yong Ding, Ning Liu & Zoya Avramova

    3. University of Nebraska Center for Biotechnology and Center for Plant Science Innovation, 1901 Vine Street, Lincoln, 68588, USA

      Laetitia Virlouvet, Jean-Jack Riethoven & Michael Fromm

    Authors
    1. Yong Ding
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Ning Liu
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Laetitia Virlouvet
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Jean-Jack Riethoven
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Michael Fromm
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Zoya Avramova
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding author

    Correspondence to Zoya Avramova.

    Additional information

    Authors’ contributions

    YD, NL and LV performed experiments, J-JR performed bioinformatics analyses. MF and ZA conceived the study and interpreted results. ZA wrote the paper. All authors read and approved the final manuscript.

    Yong Ding, Ning Liu contributed equally to this work.

    Electronic supplementary material

    Authors’ original submitted files for images

    Rights and permissions

    Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

    Reprints and Permissions

    About this article

    Cite this article

    Ding, Y., Liu, N., Virlouvet, L. et al. Four distinct types of dehydration stress memory genes in Arabidopsis thaliana. BMC Plant Biol 13, 229 (2013). https://doi.org/10.1186/1471-2229-13-229

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/1471-2229-13-229

    Keywords

    • Gene Ontology
    • Thylakoid Membrane
    • Transcriptional Response
    • Dehydration Stress
    • Memory Response

    BMC Plant Biology

    ISSN: 1471-2229

    Contact us