Integration of Transcriptomic Data Identifies Global and Cell-Specific Asthma-Related Gene Expression Signatures

Over 140,000 transcriptomic studies performed in healthy and diseased cell and tissue types, at baseline and after exposure to various agents, are available in public repositories. Integrating results of transcriptomic datasets has been an attractive approach to identify gene expression signatures that are more robust than those obtained for individual datasets, especially datasets with small sample size. We developed Reproducible Analysis and Validation of Expression Data (RAVED), a pipeline that facilitates the creation of R Markdown reports detailing reproducible analysis of publicly available transcriptomic data, and used it to analyze asthma and glucocorticoid response microarray and RNA-Seq datasets. Subsequently, we used three approaches to integrate summary statistics of these studies and identify cell/tissue-specific and global asthma and glucocorticoid-induced gene expression changes. Transcriptomic integration methods were incorporated into an online app called REALGAR, where end-users can specify datasets to integrate and quickly obtain results that may facilitate design of experimental studies.

[1]  Audrey Kauffmann,et al.  Bioinformatics Applications Note Arrayqualitymetrics—a Bioconductor Package for Quality Assessment of Microarray Data , 2022 .

[2]  Bart De Moor,et al.  BioMart and Bioconductor: a powerful link between biological databases and microarray data analysis , 2005, Bioinform..

[3]  Francesco Vallania,et al.  Methods to increase reproducibility in differential gene expression via meta-analysis , 2016, Nucleic acids research.

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

[5]  Alexander A. Morgan,et al.  Expression-based genome-wide association study links the receptor CD44 in adipose tissue with type 2 diabetes , 2012, Proceedings of the National Academy of Sciences.

[6]  Manuel A. R. Ferreira,et al.  Multiancestry association study identifies new asthma risk loci that colocalize with immune cell enhancer marks , 2017, Nature Genetics.

[7]  P. O'Byrne,et al.  Corticosteroid‐induced gene expression in allergen‐challenged asthmatic subjects taking inhaled budesonide , 2012, British journal of pharmacology.

[8]  B. Zuraw,et al.  The anti-inflammatory effect of glucocorticoids is mediated by glucocorticoid-induced leucine zipper in epithelial cells. , 2007, The Journal of allergy and clinical immunology.

[9]  Wolfgang Huber,et al.  RNA-Seq workflow: gene-level exploratory analysis and differential expression , 2015, F1000Research.

[10]  Rui Hong,et al.  Disease-Specific Integration of Omics Data to Guide Functional Validation of Genetic Associations , 2017, AMIA.

[11]  Hideaki Sugawara,et al.  The Sequence Read Archive , 2010, Nucleic Acids Res..

[12]  Sangsoo Kim,et al.  Combining multiple microarray studies and modeling interstudy variation , 2003, ISMB.

[13]  M. Vijayan,et al.  The Transcriptomics of Glucocorticoid Receptor Signaling in Developing Zebrafish , 2013, PloS one.

[14]  A. Litonjua,et al.  Vitamin D Modulates Expression of the Airway Smooth Muscle Transcriptome in Fatal Asthma , 2015, PloS one.

[15]  Manuel A. R. Ferreira,et al.  Shared genetic origin of asthma, hay fever and eczema elucidates allergic disease biology , 2017, Nature Genetics.

[16]  Andrew E. Jaffe,et al.  Bioinformatics Applications Note Gene Expression the Sva Package for Removing Batch Effects and Other Unwanted Variation in High-throughput Experiments , 2022 .

[17]  M. Loza,et al.  Regulation of T cells in airway disease by beta-agonist. , 2010, Frontiers in bioscience.

[18]  Y. Bossé,et al.  Future clinical implications emerging from recent genome-wide expression studies in asthma , 2014, Expert review of clinical immunology.

[19]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[20]  Andrea Rau,et al.  Differential meta-analysis of RNA-seq data from multiple studies , 2013, BMC Bioinformatics.

[21]  J. Cidlowski,et al.  Antiinflammatory action of glucocorticoids--new mechanisms for old drugs. , 2005, The New England journal of medicine.

[22]  Scott T. Weiss,et al.  RNA-Seq Transcriptome Profiling Identifies CRISPLD2 as a Glucocorticoid Responsive Gene that Modulates Cytokine Function in Airway Smooth Muscle Cells , 2014, PloS one.

[23]  P. Khatri,et al.  Robust classification of bacterial and viral infections via integrated host gene expression diagnostics , 2016, Science Translational Medicine.

[24]  Y. Prakash Airway smooth muscle in airway reactivity and remodeling: what have we learned? , 2013, American journal of physiology. Lung cellular and molecular physiology.

[25]  B. Himes,et al.  Using omics approaches to understand pulmonary diseases , 2017, Respiratory Research.

[26]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[27]  Rainer Breitling,et al.  A comparison of meta-analysis methods for detecting differentially expressed genes in microarray experiments , 2008, Bioinform..

[28]  Lior Pachter,et al.  Near-optimal probabilistic RNA-seq quantification , 2016, Nature Biotechnology.

[29]  Paul Theodor Pyl,et al.  HTSeq – A Python framework to work with high-throughput sequencing data , 2014 .

[30]  J. Thèze,et al.  Interleukin-7 Compartmentalizes Its Receptor Signaling Complex to Initiate CD4 T Lymphocyte Response* , 2010, The Journal of Biological Chemistry.

[31]  Lior Pachter,et al.  Differential analysis of RNA-seq incorporating quantification uncertainty , 2016, Nature Methods.

[32]  Adnan Custovic,et al.  Asthma endotypes: a new approach to classification of disease entities within the asthma syndrome. , 2011, The Journal of allergy and clinical immunology.

[33]  Rafael A. Irizarry,et al.  A framework for oligonucleotide microarray preprocessing , 2010, Bioinform..

[34]  S. Wenzel Asthma phenotypes: the evolution from clinical to molecular approaches , 2012, Nature Medicine.

[35]  Purvesh Khatri,et al.  A comprehensive time-course–based multicohort analysis of sepsis and sterile inflammation reveals a robust diagnostic gene set , 2015, Science Translational Medicine.

[36]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[37]  N. Laird,et al.  Meta-analysis in clinical trials. , 1986, Controlled clinical trials.

[38]  John P A Ioannidis,et al.  Improving Validation Practices in “Omics” Research , 2011, Science.

[39]  Dennis B. Troup,et al.  NCBI GEO: mining millions of expression profiles—database and tools , 2004, Nucleic Acids Res..

[40]  H. Hammad,et al.  The airway epithelium in asthma , 2012, Nature Medicine.

[41]  M. Pufall,et al.  Transcriptional Regulation of Human Dual Specificity Protein Phosphatase 1 (DUSP1) Gene by Glucocorticoids , 2010, PloS one.

[42]  Rainer Breitling,et al.  Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments , 2004, FEBS letters.

[43]  Richard M Myers,et al.  Genomic determination of the glucocorticoid response reveals unexpected mechanisms of gene regulation. , 2009, Genome research.