Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2017 Aug 18:6:e29655.
doi: 10.7554/eLife.29655.

Temporal network analysis identifies early physiological and transcriptomic indicators of mild drought in Brassica rapa

Affiliations

Temporal network analysis identifies early physiological and transcriptomic indicators of mild drought in Brassica rapa

Kathleen Greenham et al. Elife. .

Abstract

The dynamics of local climates make development of agricultural strategies challenging. Yield improvement has progressed slowly, especially in drought-prone regions where annual crop production suffers from episodic aridity. Underlying drought responses are circadian and diel control of gene expression that regulate daily variations in metabolic and physiological pathways. To identify transcriptomic changes that occur in the crop Brassica rapa during initial perception of drought, we applied a co-expression network approach to associate rhythmic gene expression changes with physiological responses. Coupled analysis of transcriptome and physiological parameters over a two-day time course in control and drought-stressed plants provided temporal resolution necessary for correlation of network modules with dynamic changes in stomatal conductance, photosynthetic rate, and photosystem II efficiency. This approach enabled the identification of drought-responsive genes based on their differential rhythmic expression profiles in well-watered versus droughted networks and provided new insights into the dynamic physiological changes that occur during drought.

Keywords: Brassica rapa; abiotic stress; daily rhythms; drought; photosynthesis; plant biology; transcriptomic network analysis.

PubMed Disclaimer

Conflict of interest statement

No competing interests declared.

Figures

Figure 1.
Figure 1.. Experimental design for evaluation of early effects of mild drought treatment on Brassica rapa R500 plants and soil water potential during the experiment.
(A) Schematic illustration of the experimental design with days of the experiment and collection/measurement times. Top: The time of the experiment in days after sowing (DAS) is plotted with the orange shade representing the days of mild drought treatment starting at the end of 15 DAS. Center: Purple shade (between 18 and 20 DAS) indicates the 2 days of continuous collection/measurements. The x-axis represents the 48 hr of the experiment in Zeitgeber Time (ZT0 = lights on in the growth chambers). The dotted yellow lines represent the ZT times for collection and measurements of both well-watered and droughted plants. Bottom: White and black bars represent the photoperiod in the growth chambers (14 hr/10 hr; Day/Night) for Days 3 and 4. (B) Soil water potential (Ψs) progression for well-watered (solid symbols) and mild droughted (open symbols) B. rapa R500 plants (n = at least 6) after the start of the drought treatment during the physiology (triangles) and RNA-seq (circles) time-course experiments.
Figure 2.
Figure 2.. Despite a significant (p<0.001) decrease in dry biomass accumulation, Brassica rapa R500 plants did not reach the wilting point in mild drought.
(A) Grams of dry above-ground biomass in well-watered (dark gray) and mild droughted (light gray) B. rapa R500 plants (n = 12) collected on Day 3 and Day 4. On both days, above-ground biomass was collected at ZT9. (B) Top and lateral view of well-watered and droughted B. rapa R500 plants on Day 3 and Day 4 of the experiment.
Figure 3.
Figure 3.. Specific time-of-day differences in physiology and in dynamics of Non-Structural Carbohydrates (NSC) were observed in B. rapa R500 plants subjected to mild drought relative to control plants.
(A) Photosynthetic rate, A (μmol CO2 m−2s−1) of well-watered controls (green) and mild droughted (orange) plants. (B) Stomatal conductance, gs (mmol H2O m−2s−1) of well-watered controls (green) and mild droughted (orange) plants. (C) Maximum efficiency of PSII in light adapted conditions, Fv’/Fm’ of well-watered controls (green) and mild droughted (orange) plants. (D) NSC accumulation in well-watered (green) and mild droughted (orange) plants. NSC content is expressed as dry weight and refers to the sum of leaves and cotyledon biomass. Tissue collection was carried out every 4 hr during the 2 day time course except for ZT9 when no sugar extraction was performed. White and black bars and gray shading represent the dark period in the growth chambers (14 hr/10 hr; Day/Night). All represented values are averages of at least eight replicates ± SE, asterisks represent a significant difference (p<0.01) between treatments. See Figure 3—figure supplement 1.
Figure 3—figure supplement 1.
Figure 3—figure supplement 1.. Water use efficiency (WUE) is not an accurate measure of plant response to drought.
Variation of intrinsic WUE (as A/gs) for well-watered controls (green) and mildly droughted (orange) plants. Gray shading indicates dark period.
Figure 4.
Figure 4.. Co-expression network analysis identifies modules of temporally regulated genes with similar phase patterns.
(A and B) The well-watered (A) and droughted (B) networks are shown with gene nodes colored by module membership. Network visualization was done in Cytoscape using a ForceDirected layout with an edge threshold cutoff of 0.1. Modules with significant overlap of gene membership between the well-watered and drought networks are similarly colored. (C and D) The module eigengene (ME) plots for the well-watered (C) and droughted (D) modules reveal temporal regulation of the top eight modules in each network containing ~80–85% of expressed genes. The ME plots are colored based on the modules in the network visualizations. White and black bars along the x-axis represent the day and night time points respectively.
Figure 5.
Figure 5.. Variation in gene expression over a day is greater than the effects of mild drought.
(A) Hierarchical clustering was performed on the log2 FPKM values for the Day 4 time points for well-watered and droughted plants. (B) Consensus network matrix for the well-watered (top) and droughted (bottom) networks. Consensus modules are listed along the bottom of both matrices with the well-watered and droughted modules along the side. Numbers beside the module number are the total gene counts for that module. Numbers inside the matrix are the number of genes in common between the two modules. Red numbers have significant overlap in gene count based on Fisher’s exact test with the –log(p) of the p-value encoding the coloring (0 = -log(1) and 50 = -log(1E-50)).
Figure 6.
Figure 6.. Expansion of gene expression patterns uncovered in the droughted network.
(A) The wM5 module from the well-watered network is shown with the mean log2 FPKM expression profiles for the genes in the wM5 network that overlap with dM1, dM5, and dM10 modules. (B) The nodes in the wM5 module are highlighted based on their membership in the dM1, dM5, and dM10 modules from the droughted network. The mean log2 FPKM expression values for the genes in common between the wM5 module and each droughted module is shown. Grey shading in the expression profile plots represents the standard error of the mean log2 FPKM expression levels. White and black bars along the x-axis represent the day and night time points, respectively.
Figure 7.
Figure 7.. Dynamic changes in physiology are correlated with co-expression networks.
(A and B) Correlations between module eigengenes and physiology measurements for well-watered (A) and droughted (B) plants. The numbers within the heat map represent correlations and P values (in parentheses; red, positively correlated, blue, negatively correlated) for the module – trait associations. (C) Scatterplots of gene significance versus module membership for the photosynthetic rate with dM1 (red), PSII efficiency with dM6 (blue), and stomatal conductance with dM13 (black). Significant correlations imply that hub genes within the module are also highly correlated with the physiology measure.
Figure 8.
Figure 8.. Modules positively correlated with physiology are associated with metabolism and light harvesting processes.
(A) Mean log2 FPKM expression profiles of genes in the dM1, dM5, dM10, and dM19 modules that are positively correlated with stomatal conductance and identified by JTK-CYCLE as having an amplitude change between well-watered and droughted samples. Grey lines are always well-watered. White and black bars along the x-axis represent the day and night time points, respectively. (B) Top 10 associated GO terms for dM1 (top), dM5 (middle), and dM10 (bottom). (C) Genes identified in the dM1, dM5, dM10, and dM19 modules with known Arabidopsis orthologs are shown in the network view. Nodes are colored based on module and size based on module membership. The larger nodes are highly connected within a module and have greater ‘hubness’. The nodes circled in purple were validated by Nanostring.
Figure 9.
Figure 9.. Modules negatively correlated with physiology are associated with abiotic stress and light response.
(A) Mean log2 FPKM expression profiles and of genes in the dM2, dM6, dM7, dM8, dM11, and dM20 modules that are negatively correlated with stomatal conductance and identified by JTK-CYCLE as having an amplitude change between well-watered and droughted samples. Grey lines are always well-watered. White and black bars along the x-axis represent the day and night time points, respectively. (B) Top 10 associated GO terms for dM2 (top), dM6 (middle), and dM11 (bottom). (C) Genes identified in the dM2, dM6, dM7, dM8, dM11, and dM20 modules with known Arabidopsis orthologs are shown in the network view. Nodes are colored based on module and size based on module membership. The larger nodes are highly connected within a module and have greater ‘hubness’. The nodes circled in magenta were validated by Nanostring.
Figure 10.
Figure 10.. Validation of drought-responsive genes supports alterations of light harvesting processes and metabolism.
Gene expression for a set of drought-responsive genes identified from the JTK-CYCLE filter. NanoString was performed on RNA isolated from leaf tissue harvested during the RNA-seq time course experiment. Each data point represents the mean of 5 individual plants. Gene counts were normalized to Bra021441, Bra014841, and Bra020305. Error bars represent SE. Grey shading indicates the dark period. Asterisks indicate significant difference (p<0.05). See Figure 10—figure supplement 1.
Figure 10—figure supplement 1.
Figure 10—figure supplement 1.. Validation of drought responsive genes.
Gene expression for a set of drought-responsive genes identified from the JTK-CYCLE filter. NanoString was performed on RNA isolated from leaf tissue harvested during the RNA-seq time course experiment. Each data point represents the mean of 5 individual plants. Gene counts were normalized to Bra021441, Bra014841, and Bra020305. Error bars represent SE. Black bars and grey shading indicates the dark period. Asterisks indicate significant difference (p<0.05).
Figure 10—figure supplement 2.
Figure 10—figure supplement 2.. RNA-seq expression of drought responsive genes assessed by NanoString.
FPKM expression values for the drought responsive genes analyzed by Nanostring and shown in Figure 10. Each data point represents the mean of two biological replicates. Black bars and grey shading indicates the dark period.

Comment in

  • Into the fourth dimension.
    Des Marais DL. Des Marais DL. Elife. 2017 Oct 4;6:e31328. doi: 10.7554/eLife.31328. Elife. 2017. PMID: 28976303 Free PMC article.

References

    1. Anderegg WR, Berry JA, Field CB. Linking definitions, mechanisms, and modeling of drought-induced tree death. Trends in Plant Science. 2012;17:693–700. doi: 10.1016/j.tplants.2012.09.006. - DOI - PubMed
    1. Ashraf M, Mehmood S. Response of four Brassica species to drought stress. Environmental and Experimental Botany. 1990;30:93–100. doi: 10.1016/0098-8472(90)90013-T. - DOI
    1. Attia Z, Domec JC, Oren R, Way DA, Moshelion M. Growth and physiological responses of isohydric and anisohydric poplars to drought. Journal of Experimental Botany. 2015;66:4373–4381. doi: 10.1093/jxb/erv195. - DOI - PMC - PubMed
    1. Blum A. Drought resistance, water-use efficiency, and yield potential—are they compatible, dissonant, or mutually exclusive? Australian Journal of Agricultural Research. 2005;56:1159–1168. doi: 10.1071/AR05069. - DOI
    1. Blum A. Effective use of water (EUW) and not water-use efficiency (WUE) is the target of crop yield improvement under drought stress. Field Crops Research. 2009;112:119–123. doi: 10.1016/j.fcr.2009.03.009. - DOI

Publication types

LinkOut - more resources