Lung cancer is the second most commonly occurring non-cutaneous cancer in the United States with the highest mortality rate among both men and women. In this study, we utilized three lung cancer microarray datasets generated by previous researchers to identify differentially expressed genes, altered signaling pathways, and assess the involvement of Hedgehog (Hh) pathway. The three datasets contain the expression levels of tens of thousands genes in normal lung tissues and squamous cell lung carcinoma. The datasets were combined and analyzed. The dysregulated genes and altered signaling pathways were identified using statistical methods. We then performed Fisher’s exact test on the significance of the association of Hh pathway downstream genes and squamous cell lung carcinoma.395 genes were found commonly differentially expressed in squamous cell lung carcinoma. The genes encoding fibrous structural protein keratins and cell cycle dependent genes encoding cyclin-dependent kinases were significantly up-regulated while the ones encoding LIM domains were down. Over 100 signaling pathways were implicated in squamous cell lung carcinoma, including cell cycle regulation pathway, p53 tumor-suppressor pathway, IL-8 signaling, Wnt-β-catenin pathway, mTOR signaling and EGF signaling. In addition, 37 out of 223 downstream molecules of Hh pathway were altered. The P-value from the Fisher’s exact test indicates that Hh signaling is implicated in squamous cell lung carcinoma.Numerous genes were altered and multiple pathways were dysfunctional in squamous cell lung carcinoma. Many of the altered genes have been implicated in different types of carcinoma while some are organ-specific. Hh signaling is implicated in squamous cell lung cancer, opening the door for exploring new cancer therapeutic treatment using GLI antagonist GANT 61.

The Cancer Genome Atlas (TCGA) provides a rich resource that can be used to understand how genes interact in cancer cells and has collected RNA-Seq gene expression data for many types of human cancer. However, mining the data to uncover the hidden gene-interaction patterns remains a challenge. Gaussian graphical model (GGM) is often used to learn genetic networks because it defines an undirected graphical structure, revealing the conditional dependences of genes. In this study, we focus on inferring gene interactions in 15 specific types of human cancer using RNA-Seq expression data and GGM with graphical lasso. We take advantage of the corresponding Kyoto Encyclopedia of Genes and Genomes pathway maps to define the subsets of related genes. RNA-Seq expression levels of the subsets of genes in solid cancerous tumor and normal tissues were extracted from TCGA. The gene expression data sets were cleaned and formatted, and the genetic network corresponding to each cancer type was then inferred using GGM with graphical lasso. The inferred networks reveal stable conditional dependences among the genes at the expression level and confirm the essential roles played by the genes that encode proteins involved in the two key signaling pathway phosphoinositide 3-kinase (PI3K)/AKT/mTOR and Ras/Raf/MEK/ERK in human carcinogenesis. These stable dependences elucidate the expression level interactions among the genes that are implicated in many different human cancers. The inferred genetic networks were examined to further identify and characterize a collection of gene interactions that are unique to cancer. The cross-cancer genetic interactions revealed from our study provide another set of knowledge for cancer biologists to propose strong hypotheses, so further biological investigations can be conducted effectively.

This article has been withdrawn by the authors. Evaluation by the journal with image analysis software determined that in Fig. 5A, lanes 1-3 of the RB1 immunoblot were duplicated in lanes 4 -6, lane 4 of the ERBB2 immunoblot was duplicated in lane 6, lane 5 of the STAT3 immunoblot was duplicated in lane 8, and lanes 1-3 of the actin immunoblot were flipped horizontally and reused in lanes 6 -8. The authors state that RB1, ERRB2, STAT3, and actin in Fig. 5A were created from phosphor-chemiluminescent digital imaging. The authors also state that they have replicate data supporting the conclusions of Fig. 5, A and B. In Fig. 6B, evaluation by the journal of the original data determined that single cell background fluorescence was duplicated. The authors maintain that the concern is about a single background cell not a positive cell to show transfection efficiency and is inconsequential to proving transfection.

No abstract

The PubMed database offers an extensive set of publication data that can be useful, yet inherently complex to use without automated computational techniques. Data repositories such as the Genomic Data Commons (GDC) and the Gene Expression Omnibus (GEO) offer experimental data storage and retrieval as well as curated gene expression profiles. Genetic interaction databases, including Reactome and Ingenuity Pathway Analysis, offer pathway and experiment data analysis using data curated from these publications and data repositories. We have created a method to generate and analyze consensus networks, inferring potential gene interactions, using large numbers of Bayesian networks generated by data mining publications in the PubMed database. Through the concept of network resolution, these consensus networks can be tailored to represent possible genetic interactions. We designed a set of experiments to confirm that our method is stable across variation in both sample and topological input sizes. Using gene product interactions from the KEGG pathway database and data mining PubMed publication abstracts, we verify that regardless of the network resolution or the inferred consensus network, our method is capable of inferring meaningful gene interactions through consensus Bayesian network generation with multiple, randomized topological orderings. Our method can not only confirm the existence of currently accepted interactions, but has the potential to hypothesize new ones as well. We show our method confirms the existence of known gene interactions such as JAK-STAT-PI3K-AKT-mTOR, infers novel gene interactions such as RAS- Bcl-2 and RAS-AKT, and found significant pathway-pathway interactions between the JAK-STAT signaling and Cardiac Muscle Contraction KEGG pathways.