Microarray meta-analysis to explore abiotic stress-specific gene expression patterns in Arabidopsis

BackgroundAbiotic stresses are the major limiting factors that affect plant growth, development, yield and final quality. Deciphering the underlying mechanisms of plants’ adaptations to stresses using few datasets might overlook the different aspects of stress tolerance in plants, which might be simultaneously and consequently operated in the system. Fortunately, the accumulated microarray expression data offer an opportunity to infer abiotic stress-specific gene expression patterns through meta-analysis. In this study, we propose to combine microarray gene expression data under control, cold, drought, heat, and salt conditions and determined modules (gene sets) of genes highly associated with each other according to the observed expression data.ResultsBy analyzing the expression variations of the Eigen genes from different conditions, we had identified two, three, and five gene modules as cold-, heat-, and salt-specific modules, respectively. Most of the cold- or heat-specific modules were differentially expressed to a particular degree in shoot samples, while most of the salt-specific modules were differentially expressed to a particular degree in root samples. A gene ontology (GO) analysis on the stress-specific modules suggested that the gene modules exclusively enriched stress-related GO terms and that different genes under the same GO terms may be alternatively disturbed in different conditions. The gene regulatory events for two genes, DREB1A and DEAR1, in the cold-specific gene module had also been validated, as evidenced through the literature search.ConclusionsOur protocols study the specificity of the gene modules that were specifically activated under a particular type of abiotic stress. The biplot can also assist to visualize the stress-specific gene modules. In conclusion, our approach has the potential to further elucidate mechanisms in plants and beneficial for future experiments design under different abiotic stresses.

[1]  Kazuo Shinozaki,et al.  Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. , 2006, Annual review of plant biology.

[2]  Anthony Trewavas,et al.  Aspects of plant intelligence. , 2003, Annals of botany.

[3]  Todd P. Michael,et al.  Analysis of Global Gene Expression in Brachypodium distachyon Reveals Extensive Network Plasticity in Response to Abiotic Stress , 2014, PloS one.

[4]  Li-Yu Daisy Liu.  Coefficient of intrinsic dependence: a new measure of association , 2005 .

[5]  Terence P. Speed,et al.  A comparison of normalization methods for high density oligonucleotide array data based on variance and bias , 2003, Bioinform..

[6]  K. Gabriel,et al.  The biplot graphic display of matrices with application to principal component analysis , 1971 .

[7]  R. Kream,et al.  Comparing Bioinformatic Gene Expression Profiling Methods: Microarray and RNA-Seq , 2014, Medical science monitor basic research.

[8]  J. Young,et al.  Plantæ , 1871, Transactions of the Glasgow Geological Society.

[9]  K. Shinozaki,et al.  Toward understanding transcriptional regulatory networks in abiotic stress responses and tolerance in rice , 2012, Rice.

[10]  P. Urwin,et al.  The interaction of plant biotic and abiotic stresses: from genes to the field. , 2012, Journal of experimental botany.

[11]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[12]  Reinhard Schneider,et al.  Using graph theory to analyze biological networks , 2011, BioData Mining.

[13]  Chien-Yu Chen,et al.  Statistical identification of gene association by CID in application of constructing ER regulatory network , 2009, BMC Bioinformatics.

[14]  Wen Huang,et al.  The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant , 2001, Nucleic Acids Res..

[15]  K. Shinozaki,et al.  The transcriptional regulatory network in the drought response and its crosstalk in abiotic stress responses including drought, cold, and heat , 2014, Front. Plant Sci..

[16]  Kazuo Shinozaki,et al.  DEAR1, a transcriptional repressor of DREB protein that mediates plant defense and freezing stress responses in Arabidopsis , 2009, Journal of Plant Research.

[17]  P. Ronald.  Plant Genetics, Sustainable Agriculture and Global Food Security , 2011, Genetics.

[18]  Douglas G Altman,et al.  Key Issues in Conducting a Meta-Analysis of Gene Expression Microarray Datasets , 2008, PLoS medicine.

[19]  Burkhard Morgenstern,et al.  Meta-Analysis of Pathway Enrichment: Combining Independent and Dependent Omics Data Sets , 2014, PloS one.

[20]  Manjit S. Kang,et al.  GGE Biplot Analysis: A Graphical Tool for Breeders, Geneticists, and Agronomists , 2002 .

[21]  S. Horvath,et al.  A General Framework for Weighted Gene Co-Expression Network Analysis , 2005, Statistical applications in genetics and molecular biology.

[22]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[23]  T. Speed,et al.  Summaries of Affymetrix GeneChip probe level data. , 2003, Nucleic acids research.

[24]  M. Sachs,et al.  Alteration of Gene Expression During Environmental Stress in Plants , 1986 .

[25]  Karl Pearson F.R.S..  LIII. On lines and planes of closest fit to systems of points in space , 1901 .

[26]  S. Rasmussen,et al.  Transcriptome Responses to Combinations of Stresses in Arabidopsis1[W][OA] , 2013, Plant Physiology.

[27]  R. Albert.  Scale-free networks in cell biology , 2005, Journal of Cell Science.

[28]  Ya-Chun Hsiao,et al.  A Stepwise Approach of Finding Dependent Variables via Coefficient of Intrinsic Dependence , 2016, J. Comput. Biol..

[29]  L. Liu,et al.  Identifying Gene Set Association Enrichment Using the Coefficient of Intrinsic Dependence , 2013, PloS one.

[30]  Steve Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[31]  Jessica Gurevitch,et al.  Meta‐analysis and meta‐regression of transcriptomic responses to water stress in Arabidopsis , 2016, The Plant journal : for cell and molecular biology.

[32]  Burkhard Morgenstern,et al.  MarVis-Pathway: integrative and exploratory pathway analysis of non-targeted metabolomics data , 2014, Metabolomics.

[33]  Piero Fariselli,et al.  Blurring contact maps of thousands of proteins: what we can learn by reconstructing 3D structure , 2011, BioData Mining.

[34]  Zhou Du,et al.  agriGO: a GO analysis toolkit for the agricultural community , 2010, Nucleic Acids Res..

[35]  King-Jen Chang,et al.  Major Functional Transcriptome of an Inferred Center Regulator of an ER(−) Breast Cancer Model System , 2012, Cancer informatics.

[36]  H. Hotelling.  Analysis of a complex of statistical variables into principal components. , 1933 .

[37]  Arnold M Saxton,et al.  Comparison of threshold selection methods for microarray gene co-expression matrices , 2009, BMC Research Notes.

[38]  Joachim Kilian,et al.  Plant Core Environmental Stress Response Genes Are Systemically Coordinated during Abiotic Stresses , 2013, International journal of molecular sciences.

[39]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[40]  Jean Yee Hwa Yang,et al.  Comparison study of microarray meta-analysis methods , 2010, BMC Bioinformatics.

[41]  Robert H. Halstead,et al.  Matrix Computations , 2011, Encyclopedia of Parallel Computing.

[42]  Lin Song,et al.  Comparison of co-expression measures: mutual information, correlation, and model based indices , 2012, BMC Bioinformatics.