Application of Weighted Gene Co-expression Network Analysis for Data from Paired Design

Investigating how genes jointly affect complex human diseases is important, yet challenging. The network approach (e.g., weighted gene co-expression network analysis (WGCNA)) is a powerful tool. However, genomic data usually contain substantial batch effects, which could mask true genomic signals. Paired design is a powerful tool that can reduce batch effects. However, it is currently unclear how to appropriately apply WGCNA to genomic data from paired design. In this paper, we modified the current WGCNA pipeline to analyse high-throughput genomic data from paired design. We illustrated the modified WGCNA pipeline by analysing the miRNA dataset provided by Shiah et al. (2014), which contains forty oral squamous cell carcinoma (OSCC) specimens and their matched non-tumourous epithelial counterparts. OSCC is the sixth most common cancer worldwide. The modified WGCNA pipeline identified two sets of novel miRNAs associated with OSCC, in addition to the existing miRNAs reported by Shiah et al. (2014). Thus, this work will be of great interest to readers of various scientific disciplines, in particular, genetic and genomic scientists as well as medical scientists working on cancer.

[1]  Edwin Wang,et al.  Understanding genomic alterations in cancer genomes using an integrative network approach. , 2013, Cancer letters.

[2]  J. Castle,et al.  An integrative genomics approach to infer causal associations between gene expression and disease , 2005, Nature Genetics.

[3]  E. Davidson,et al.  The hardwiring of development: organization and function of genomic regulatory systems. , 1997, Development.

[4]  Edwin Wang,et al.  Signaling network analysis of ubiquitin-mediated proteins suggests correlations between the 26S proteasome and tumor progression. , 2009, Molecular bioSystems.

[5]  Wei Shi,et al.  Comprehensive MicroRNA Profiling for Head and Neck Squamous Cell Carcinomas , 2010, Clinical Cancer Research.

[6]  Hongwei Li,et al.  Differential mRNA expression profiling of oral squamous cell carcinoma by high-throughput RNA sequencing , 2015, Journal of biomedical research.

[7]  Y. Leea,et al.  Analysis of oncogenic signaling networks in glioblastoma identifies ASPM as a molecular target , 2006 .

[8]  S. Horvath,et al.  Weighted gene coexpression network analysis strategies applied to mouse weight , 2007, Mammalian Genome.

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

[10]  Rajiv Sarin,et al.  Mutational landscape of gingivo-buccal oral squamous cell carcinoma reveals new recurrently-mutated genes and molecular subgroups , 2013, Nature Communications.

[11]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[12]  Alex E. Lash,et al.  Gene Expression Omnibus: NCBI gene expression and hybridization array data repository , 2002, Nucleic Acids Res..

[13]  Edwin Wang,et al.  Network Analysis Reveals A Signaling Regulatory Loop in the PIK3CA-mutated Breast Cancer Predicting Survival Outcome , 2017, Genom. Proteom. Bioinform..

[14]  S. Jiang,et al.  Downregulated miR329 and miR410 promote the proliferation and invasion of oral squamous cell carcinoma by targeting Wnt-7b. , 2014, Cancer research.

[15]  M. Manikandan,et al.  Oral squamous cell carcinoma: microRNA expression profiling and integrative analyses for elucidation of tumourigenesis mechanism , 2016, Molecular Cancer.

[16]  C. K. Hsiao,et al.  miRSystem: An Integrated System for Characterizing Enriched Functions and Pathways of MicroRNA Targets , 2012, PloS one.

[17]  Stefan Joos,et al.  Activation of MAP kinase signaling through ERK5 but not ERK1 expression is associated with lymph node metastases in oral squamous cell carcinoma (OSCC). , 2008, Neoplasia.

[18]  G. Rubin,et al.  The Role of the Genome Project in Determining Gene Function: Insights from Model Organisms , 1996, Cell.

[19]  Edwin Wang,et al.  Signaling network assessment of mutations and copy number variations predict breast cancer subtype-specific drug targets. , 2013, Cell reports.

[20]  S. Horvath,et al.  Variations in DNA elucidate molecular networks that cause disease , 2008, Nature.

[21]  Jinfeng Zou,et al.  Identification and Construction of Combinatory Cancer Hallmark-Based Gene Signature Sets to Predict Recurrence and Chemotherapy Benefit in Stage II Colorectal Cancer. , 2016, JAMA oncology.

[22]  S. Shenouda,et al.  MicroRNA function in cancer: oncogene or a tumor suppressor? , 2009, Cancer and Metastasis Reviews.

[23]  Sean R. Davis,et al.  NCBI GEO: archive for functional genomics data sets—update , 2012, Nucleic Acids Res..

[24]  A. Kolokythas,et al.  microRNA from brush biopsy to characterize oral squamous cell carcinoma epithelium , 2016, Cancer medicine.

[25]  E. Wang,et al.  MicroRNA Regulatory Patterns on the Human Metabolic Network~!2008-10-09~!2008-10-22~!2008-11-28~! , 2008 .

[26]  J. Ferlay,et al.  Estimates of the worldwide mortality from eighteen major cancers in 1985. Implications for prevention and projections of future burden , 1993, International journal of cancer.

[27]  J. Picquenot,et al.  Expression of p53 family members and CD44 in oral squamous cell carcinoma (OSCC) in relation to tumorigenesis. , 2010, Histology and histopathology.

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

[29]  Joshua M. Stuart,et al.  A Gene-Coexpression Network for Global Discovery of Conserved Genetic Modules , 2003, Science.