Non‑Small‑Cell Lung Cancer Pathological Subtype‑Related Gene

Total Page:16

File Type:pdf, Size:1020Kb

Non‑Small‑Cell Lung Cancer Pathological Subtype‑Related Gene 356 MOLECULAR AND CLINICAL ONCOLOGY 8: 356-361, 2018 Non‑small‑cell lung cancer pathological subtype‑related gene selection and bioinformatics analysis based on gene expression profiles JIANGPENG CHEN1, XIAOQI DONG2, XUN LEI1, YINYIN XIA1, QING ZENG1, PING QUE1, XIAOYAN WEN1, SHAN HU1 and BIN PENG1 1School of Public Health and Management, Chongqing Medical University, Chongqing 400016; 2Department of Respiratory Diseases, The First Affiliated Hospital of Zhejiang University, Hangzhou, Zhejiang 310003, P.R. China Received March 21, 2016; Accepted November 21, 2017 DOI: 10.3892/mco.2017.1516 Abstract. Lung cancer is one of the most common malignant Introduction diseases and a major threat to public health on a global scale. Non-small-cell lung cancer (NSCLC) has a higher degree of Lung cancer is one of the most common malignant diseases malignancy and a lower 5-year survival rate compared with that and a major threat to public health on a global scale. The main of small-cell lung cancer. NSCLC may be mainly divided into types of lung cancer are small-cell lung cancer (SCLC) and non- two pathological subtypes, adenocarcinoma and squamous cell small-cell lung cancer (NSCLC). NSCLC has a higher degree carcinoma. The aim of the present study was to identify disease of malignancy and a lower 5-year survival rate compared with genes based on the gene expression profile and the shortest path SCLC, and may be divided into two major histopathological analysis of weighted functional protein association networks subtypes, namely adenocarcinoma (ADC) and squamous cell with the existing protein-protein interaction data from the Search carcinoma (SCC). As NSCLC is a complex disease, the iden- Tool for the Retrieval of Interacting Genes. The gene expression tification of disease genes has been among the main goals of profile (GSE10245) was downloaded from the National Center biologists and clinicians. Although our understanding of lung for Biotechnology Information Gene Expression Omnibus data- cancer, particularly NSCLC, has improved significantly in base, including 40 lung adenocarcinoma and 18 lung squamous recent years, several questions require further elucidation. cell carcinoma tissues. A total of 8 disease genes were identi- With the completion of the Human Genome Project, scien- fied using Naïve Bayesian Classifier based on the Maximum tific research has entered the post-genome era. Furthermore, Relevance Minimum Redundancy feature selection method with the advances in biological research and the development following preprocessing. An additional 21 candidate genes of high throughput biotechnologies, the quantities of candidate were selected using the shortest path analysis with Dijkstra's genes and numerous genomic loci have been reported to play a algorithm. The AURKA and SLC7A2 genes were selected three vital role in the pathogenesis of NSCLC. However, the mecha- and two times in the shortest path analysis, respectively. All nism underlying this disease has not yet been fully elucidated. those genes participate in a number of important pathways, As is well-known, different analytical methods may result in such as oocyte meiosis, cell cycle and cancer pathways with different consequences due to the ̔small N large P̓ problem in Kyoto Encyclopedia of Genes and Genomes pathway enrich- microarray data, i.e., the gene expression data usually come with ment analysis. The present findings may provide novel insights only dozens of tissue samples, but with up to tens of thousands into the pathogenesis of NSCLC and enable the development of gene features. This extreme sparseness is considered to of novel therapeutic strategies. However, further investigation is significantly compromise the performance of a classifier (1). As required to confirm these findings. a result, the ability to extract a subset of informative genes while removing irrelevant or redundant genes is crucial for accurate classification (2). Hence, the methodology research is of great importance and the reuse of microarray data is feasible and extremely valuable. It is meaningful to identify disease-specific genes to help provide a novel therapeutic target for tumor treat- ment and elucidate the mechanism of the disease. Correspondence to: Dr Bin Peng, School of Public Health and Management, Chongqing Medical University, 1 Yixueyuan Road, The aim of the present study was to identify the NSCLC Chongqing 400016, P.R. China subtype-related genes based on the classification method E-mail: [email protected] and the shortest path analysis based on the protein-protein interaction (PPI) network, and to investigate their underlying Key words: non-small-cell lung cancer, microarray, bioinformatics functions by Gene Ontology (GO) and Kyoto Encyclopedia of analysis, feature selection, pathway Genes and Genomes (KEGG) pathway enrichment analysis with the Database for Annotation, Visualization and Integrated CHEN et al: NON-SMALL CELL LUNG CANCER RELATED GENE SELECTION AND BIOINFORMATICS ANALYSIS 357 Discovery (DAVID). This workflow may be of value for Let S denote the subset of features we are seeking. The identifying disease genes and may provide insight into the maximum relevance condition is pathogenesis of malignant tumors. To the best of our knowl- edge, this investigation is the first of its type in lung cancer (2) studies, and it may be useful for further lung cancer candidate gene verification, thus enabling a better understanding of the molecular mechanisms of NSCLC and the development more The minimum redundancy condition is effective diagnostic and intervention methods. (3) Materials and methods Data source. The gene expression dataset (GDS3627, acces- sion no. GSE10245) was downloaded from the National Center where S is used to represent the number of features in S, for Biotechnology Information Gene Expression Omnibus and c represents the targeted classes. (GEO) from a research conducted by Kuner et al (3) on lung The mutual information quotient criterion may be cancer. The data were processed with the Affymatrix GPL570 described as follows. platform and generated from 58 tissue samples (40 ADC and 18 SCC). All the tissue samples underwent careful analysis of maxФ (D,R), Ф = D/R (4) the histopathology, estimation of tumor and stromal content and subsequent macro-dissection of the tumor regions prior to This optimization may be computed efficiently by heuristic enrollment in the microarray experiment. No ethics committee algorithm approval is required to obtain these data, since they are publicly available. In addition, only expression data are presented, and (5) no personal patient information is revealed. Preprocessing of microarray data. Data preprocessing, including where xj ∈ XF -Sm-1, XF represents the original feature set. quantile normalization and log2 transformation was completed by the GEOquery package (4). In order to minimize false posi- Naïve Bayes Classifier.The Naïve Bayes Classifier is a simple tives, the GeneSelector package in R programme was used to probabilistic classifier based on applying Bayes' theorem with rank gene lists. Briefly, selection of differential expressed genes strong independence assumptions between features. For a may be highly affected by the statistical methods used, and the sample s with m gene expression levels {g1, g2,....., gm} for the overlap of genes selected using the same significance value may m features, the posterior probability that s belongs to class ck is be quite low. GeneSelector ranks genes on the basis of how well they perform in a selected number of statistical tests and mini- mize the number of false positives (5). Ranked lists based on (6) ordinary Student's t-test, significance analysis of microarrays, and fold-change were first generated for each of the criteria. The Markov chain model within the GeneSelector package was used where p(gt|ck) are conditional tables estimated from training to aggregate the ranked lists of all criteria into a list of genes that examples. are significant, independently of the method used. To reduce The Naïve Bayes Classifier was used after all the features noise and due to the fact that prediction methods such as Naïve ranked by the mRMR. A series of classifiers were then Bayes Classifier prefer categorical data, the gene expression data constructed when the features were added into the classifier were discretized using the respective µ (mean) and σ (standard one by one. For example, the first classifier was constructed by deviation): Any data larger than µ + σ/2 were transformed to the first feature, the second classifier by the first two features, state 1; any data between µ - σ/2 and µ + σ/2 were transformed and the rest in the same manner. to state 0; and any data smaller than µ ‑ σ/2 were transformed Using incremental feature selection (IFS), the number N to state -1. These three states correspond to the overexpression, may be determined. Its idea is to compare prediction accuracy baseline, and underexpression of genes, respectively. According defined in the following selection among different Ns, and to the previous studies, discretization of gene expression data select the one with the highest accuracy. generally leads to better prediction accuracy. The prediction accuracy was formulated by Maximum relevance minimum redundancy (mRMR) feature (7) selection method. The mRMR feature selection method was proposed by Ding et al (6). For discrete variables, the mutual information I of two variables x and y is defined based on where TP represents the number of true positives, TN
Recommended publications
  • Genome-Wide Analysis of 5-Hmc in the Peripheral Blood of Systemic Lupus Erythematosus Patients Using an Hmedip-Chip
    INTERNATIONAL JOURNAL OF MOLECULAR MEDICINE 35: 1467-1479, 2015 Genome-wide analysis of 5-hmC in the peripheral blood of systemic lupus erythematosus patients using an hMeDIP-chip WEIGUO SUI1*, QIUPEI TAN1*, MING YANG1, QIANG YAN1, HUA LIN1, MINGLIN OU1, WEN XUE1, JIEJING CHEN1, TONGXIANG ZOU1, HUANYUN JING1, LI GUO1, CUIHUI CAO1, YUFENG SUN1, ZHENZHEN CUI1 and YONG DAI2 1Guangxi Key Laboratory of Metabolic Diseases Research, Central Laboratory of Guilin 181st Hospital, Guilin, Guangxi 541002; 2Clinical Medical Research Center, the Second Clinical Medical College of Jinan University (Shenzhen People's Hospital), Shenzhen, Guangdong 518020, P.R. China Received July 9, 2014; Accepted February 27, 2015 DOI: 10.3892/ijmm.2015.2149 Abstract. Systemic lupus erythematosus (SLE) is a chronic, Introduction potentially fatal systemic autoimmune disease characterized by the production of autoantibodies against a wide range Systemic lupus erythematosus (SLE) is a typical systemic auto- of self-antigens. To investigate the role of the 5-hmC DNA immune disease, involving diffuse connective tissues (1) and modification with regard to the onset of SLE, we compared is characterized by immune inflammation. SLE has a complex the levels 5-hmC between SLE patients and normal controls. pathogenesis (2), involving genetic, immunologic and envi- Whole blood was obtained from patients, and genomic DNA ronmental factors. Thus, it may result in damage to multiple was extracted. Using the hMeDIP-chip analysis and valida- tissues and organs, especially the kidneys (3). SLE arises from tion by quantitative RT-PCR (RT-qPCR), we identified the a combination of heritable and environmental influences. differentially hydroxymethylated regions that are associated Epigenetics, the study of changes in gene expression with SLE.
    [Show full text]
  • Bioinformatic Analysis of Structure and Function of LIM Domains of Human Zyxin Family Proteins
    International Journal of Molecular Sciences Article Bioinformatic Analysis of Structure and Function of LIM Domains of Human Zyxin Family Proteins M. Quadir Siddiqui 1,† , Maulik D. Badmalia 1,† and Trushar R. Patel 1,2,3,* 1 Alberta RNA Research and Training Institute, Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive, Lethbridge, AB T1K 3M4, Canada; [email protected] (M.Q.S.); [email protected] (M.D.B.) 2 Department of Microbiology, Immunology and Infectious Disease, Cumming School of Medicine, University of Calgary, 3330 Hospital Drive, Calgary, AB T2N 4N1, Canada 3 Li Ka Shing Institute of Virology, University of Alberta, Edmonton, AB T6G 2E1, Canada * Correspondence: [email protected] † These authors contributed equally to the work. Abstract: Members of the human Zyxin family are LIM domain-containing proteins that perform critical cellular functions and are indispensable for cellular integrity. Despite their importance, not much is known about their structure, functions, interactions and dynamics. To provide insights into these, we used a set of in-silico tools and databases and analyzed their amino acid sequence, phylogeny, post-translational modifications, structure-dynamics, molecular interactions, and func- tions. Our analysis revealed that zyxin members are ohnologs. Presence of a conserved nuclear export signal composed of LxxLxL/LxxxLxL consensus sequence, as well as a possible nuclear localization signal, suggesting that Zyxin family members may have nuclear and cytoplasmic roles. The molecular modeling and structural analysis indicated that Zyxin family LIM domains share Citation: Siddiqui, M.Q.; Badmalia, similarities with transcriptional regulators and have positively charged electrostatic patches, which M.D.; Patel, T.R.
    [Show full text]
  • Association of Gene Ontology Categories with Decay Rate for Hepg2 Experiments These Tables Show Details for All Gene Ontology Categories
    Supplementary Table 1: Association of Gene Ontology Categories with Decay Rate for HepG2 Experiments These tables show details for all Gene Ontology categories. Inferences for manual classification scheme shown at the bottom. Those categories used in Figure 1A are highlighted in bold. Standard Deviations are shown in parentheses. P-values less than 1E-20 are indicated with a "0". Rate r (hour^-1) Half-life < 2hr. Decay % GO Number Category Name Probe Sets Group Non-Group Distribution p-value In-Group Non-Group Representation p-value GO:0006350 transcription 1523 0.221 (0.009) 0.127 (0.002) FASTER 0 13.1 (0.4) 4.5 (0.1) OVER 0 GO:0006351 transcription, DNA-dependent 1498 0.220 (0.009) 0.127 (0.002) FASTER 0 13.0 (0.4) 4.5 (0.1) OVER 0 GO:0006355 regulation of transcription, DNA-dependent 1163 0.230 (0.011) 0.128 (0.002) FASTER 5.00E-21 14.2 (0.5) 4.6 (0.1) OVER 0 GO:0006366 transcription from Pol II promoter 845 0.225 (0.012) 0.130 (0.002) FASTER 1.88E-14 13.0 (0.5) 4.8 (0.1) OVER 0 GO:0006139 nucleobase, nucleoside, nucleotide and nucleic acid metabolism3004 0.173 (0.006) 0.127 (0.002) FASTER 1.28E-12 8.4 (0.2) 4.5 (0.1) OVER 0 GO:0006357 regulation of transcription from Pol II promoter 487 0.231 (0.016) 0.132 (0.002) FASTER 6.05E-10 13.5 (0.6) 4.9 (0.1) OVER 0 GO:0008283 cell proliferation 625 0.189 (0.014) 0.132 (0.002) FASTER 1.95E-05 10.1 (0.6) 5.0 (0.1) OVER 1.50E-20 GO:0006513 monoubiquitination 36 0.305 (0.049) 0.134 (0.002) FASTER 2.69E-04 25.4 (4.4) 5.1 (0.1) OVER 2.04E-06 GO:0007050 cell cycle arrest 57 0.311 (0.054) 0.133 (0.002)
    [Show full text]
  • Supplementary Materials
    Supplementary materials Supplementary Table S1: MGNC compound library Ingredien Molecule Caco- Mol ID MW AlogP OB (%) BBB DL FASA- HL t Name Name 2 shengdi MOL012254 campesterol 400.8 7.63 37.58 1.34 0.98 0.7 0.21 20.2 shengdi MOL000519 coniferin 314.4 3.16 31.11 0.42 -0.2 0.3 0.27 74.6 beta- shengdi MOL000359 414.8 8.08 36.91 1.32 0.99 0.8 0.23 20.2 sitosterol pachymic shengdi MOL000289 528.9 6.54 33.63 0.1 -0.6 0.8 0 9.27 acid Poricoic acid shengdi MOL000291 484.7 5.64 30.52 -0.08 -0.9 0.8 0 8.67 B Chrysanthem shengdi MOL004492 585 8.24 38.72 0.51 -1 0.6 0.3 17.5 axanthin 20- shengdi MOL011455 Hexadecano 418.6 1.91 32.7 -0.24 -0.4 0.7 0.29 104 ylingenol huanglian MOL001454 berberine 336.4 3.45 36.86 1.24 0.57 0.8 0.19 6.57 huanglian MOL013352 Obacunone 454.6 2.68 43.29 0.01 -0.4 0.8 0.31 -13 huanglian MOL002894 berberrubine 322.4 3.2 35.74 1.07 0.17 0.7 0.24 6.46 huanglian MOL002897 epiberberine 336.4 3.45 43.09 1.17 0.4 0.8 0.19 6.1 huanglian MOL002903 (R)-Canadine 339.4 3.4 55.37 1.04 0.57 0.8 0.2 6.41 huanglian MOL002904 Berlambine 351.4 2.49 36.68 0.97 0.17 0.8 0.28 7.33 Corchorosid huanglian MOL002907 404.6 1.34 105 -0.91 -1.3 0.8 0.29 6.68 e A_qt Magnogrand huanglian MOL000622 266.4 1.18 63.71 0.02 -0.2 0.2 0.3 3.17 iolide huanglian MOL000762 Palmidin A 510.5 4.52 35.36 -0.38 -1.5 0.7 0.39 33.2 huanglian MOL000785 palmatine 352.4 3.65 64.6 1.33 0.37 0.7 0.13 2.25 huanglian MOL000098 quercetin 302.3 1.5 46.43 0.05 -0.8 0.3 0.38 14.4 huanglian MOL001458 coptisine 320.3 3.25 30.67 1.21 0.32 0.9 0.26 9.33 huanglian MOL002668 Worenine
    [Show full text]
  • Follow-Up of Loci from the International Genomics of Alzheimer’S Disease Project Identifies TRIP4 As a Novel Susceptibility Gene
    OPEN Citation: Transl Psychiatry (2014) 4, e358; doi:10.1038/tp.2014.2 © 2014 Macmillan Publishers Limited All rights reserved 2158-3188/14 www.nature.com/tp ORIGINAL ARTICLE Follow-up of loci from the International Genomics of Alzheimer’s Disease Project identifies TRIP4 as a novel susceptibility gene A Ruiz1,32, S Heilmann2,3,32, T Becker4,5,32, I Hernández1, H Wagner6, M Thelen6, A Mauleón1, M Rosende-Roca1, C Bellenguez7,8,9, JC Bis10, D Harold11, A Gerrish11, R Sims11, O Sotolongo-Grau1, A Espinosa1, M Alegret1, JL Arrieta12, A Lacour4, M Leber4, J Becker6, A Lafuente1, S Ruiz1, L Vargas1, O Rodríguez1, G Ortega1, M-A Dominguez1, IGAP33, R Mayeux13,14, JL Haines15,16, MA Pericak-Vance17,18, LA Farrer19,20,21,22,23, GD Schellenberg24, V Chouraki23, LJ Launer25, C van Duijn26,27,28, S Seshadri23, C Antúnez29, MM Breteler4, M Serrano-Ríos30, F Jessen4,6, L Tárraga1, MM Nöthen2,3, W Maier4,6, M Boada1,31 and A Ramírez2,6 To follow-up loci discovered by the International Genomics of Alzheimer’s Disease Project, we attempted independent replication of 19 single nucleotide polymorphisms (SNPs) in a large Spanish sample (Fundació ACE data set; 1808 patients and 2564 controls). Our results corroborate association with four SNPs located in the genes INPP5D, MEF2C, ZCWPW1 and FERMT2, respectively. Of these, ZCWPW1 was the only SNP to withstand correction for multiple testing (P = 0.000655). Furthermore, we identify TRIP4 (rs74615166) as a novel genome-wide significant locus for Alzheimer’s disease risk (odds ratio = 1.31; confidence interval 95% (1.19–1.44); P = 9.74 × 10−9).
    [Show full text]
  • Genomic and Transcriptome Analysis Revealing an Oncogenic Functional Module in Meningiomas
    Neurosurg Focus 35 (6):E3, 2013 ©AANS, 2013 Genomic and transcriptome analysis revealing an oncogenic functional module in meningiomas XIAO CHANG, PH.D.,1 LINGLING SHI, PH.D.,2 FAN GAO, PH.D.,1 JONATHAN RUssIN, M.D.,3 LIYUN ZENG, PH.D.,1 SHUHAN HE, B.S.,3 THOMAS C. CHEN, M.D.,3 STEVEN L. GIANNOTTA, M.D.,3 DANIEL J. WEISENBERGER, PH.D.,4 GAbrIEL ZADA, M.D.,3 KAI WANG, PH.D.,1,5,6 AND WIllIAM J. MAck, M.D.1,3 1Zilkha Neurogenetic Institute, Keck School of Medicine, University of Southern California, Los Angeles, California; 2GHM Institute of CNS Regeneration, Jinan University, Guangzhou, China; 3Department of Neurosurgery, Keck School of Medicine, University of Southern California, Los Angeles, California; 4USC Epigenome Center, Keck School of Medicine, University of Southern California, Los Angeles, California; 5Department of Psychiatry, Keck School of Medicine, University of Southern California, Los Angeles, California; and 6Division of Bioinformatics, Department of Preventive Medicine, Keck School of Medicine, University of Southern California, Los Angeles, California Object. Meningiomas are among the most common primary adult brain tumors. Although typically benign, roughly 2%–5% display malignant pathological features. The key molecular pathways involved in malignant trans- formation remain to be determined. Methods. Illumina expression microarrays were used to assess gene expression levels, and Illumina single- nucleotide polymorphism arrays were used to identify copy number variants in benign, atypical, and malignant me- ningiomas (19 tumors, including 4 malignant ones). The authors also reanalyzed 2 expression data sets generated on Affymetrix microarrays (n = 68, including 6 malignant ones; n = 56, including 3 malignant ones).
    [Show full text]
  • Mouse Trip4 Conditional Knockout Project (CRISPR/Cas9)
    https://www.alphaknockout.com Mouse Trip4 Conditional Knockout Project (CRISPR/Cas9) Objective: To create a Trip4 conditional knockout Mouse model (C57BL/6J) by CRISPR/Cas-mediated genome engineering. Strategy summary: The Trip4 gene (NCBI Reference Sequence: NM_019797 ; Ensembl: ENSMUSG00000032386 ) is located on Mouse chromosome 9. 16 exons are identified, with the ATG start codon in exon 4 and the TGA stop codon in exon 16 (Transcript: ENSMUST00000119245). Exon 5 will be selected as conditional knockout region (cKO region). Deletion of this region should result in the loss of function of the Mouse Trip4 gene. To engineer the targeting vector, homologous arms and cKO region will be generated by PCR using BAC clone RP23-363H13 as template. Cas9, gRNA and targeting vector will be co-injected into fertilized eggs for cKO Mouse production. The pups will be genotyped by PCR followed by sequencing analysis. Note: Exon 5 starts from about 5.85% of the coding region. The knockout of Exon 5 will result in frameshift of the gene. The size of intron 4 for 5'-loxP site insertion: 3867 bp, and the size of intron 5 for 3'-loxP site insertion: 1682 bp. The size of effective cKO region: ~670 bp. The cKO region does not have any other known gene. Page 1 of 8 https://www.alphaknockout.com Overview of the Targeting Strategy Wildtype allele gRNA region 5' gRNA region 3' 1 5 6 16 Targeting vector Targeted allele Constitutive KO allele (After Cre recombination) Legends Exon of mouse Trip4 Homology arm cKO region loxP site Page 2 of 8 https://www.alphaknockout.com Overview of the Dot Plot Window size: 10 bp Forward Reverse Complement Sequence 12 Note: The sequence of homologous arms and cKO region is aligned with itself to determine if there are tandem repeats.
    [Show full text]
  • Differential Expression Profile Prioritization of Positional Candidate Glaucoma Genes the GLC1C Locus
    LABORATORY SCIENCES Differential Expression Profile Prioritization of Positional Candidate Glaucoma Genes The GLC1C Locus Frank W. Rozsa, PhD; Kathleen M. Scott, BS; Hemant Pawar, PhD; John R. Samples, MD; Mary K. Wirtz, PhD; Julia E. Richards, PhD Objectives: To develop and apply a model for priori- est because of moderate expression and changes in tization of candidate glaucoma genes. expression. Transcription factor ZBTB38 emerges as an interesting candidate gene because of the overall expres- Methods: This Affymetrix GeneChip (Affymetrix, Santa sion level, differential expression, and function. Clara, Calif) study of gene expression in primary cul- ture human trabecular meshwork cells uses a positional Conclusions: Only1geneintheGLC1C interval fits our differential expression profile model for prioritization of model for differential expression under multiple glau- candidate genes within the GLC1C genetic inclusion in- coma risk conditions. The use of multiple prioritization terval. models resulted in filtering 7 candidate genes of higher interest out of the 41 known genes in the region. Results: Sixteen genes were expressed under all condi- tions within the GLC1C interval. TMEM22 was the only Clinical Relevance: This study identified a small sub- gene within the interval with differential expression in set of genes that are most likely to harbor mutations that the same direction under both conditions tested. Two cause glaucoma linked to GLC1C. genes, ATP1B3 and COPB2, are of interest in the con- text of a protein-misfolding model for candidate selec- tion. SLC25A36, PCCB, and FNDC6 are of lesser inter- Arch Ophthalmol. 2007;125:117-127 IGH PREVALENCE AND PO- identification of additional GLC1C fami- tential for severe out- lies7,18-20 who provide optimal samples for come combine to make screening candidate genes for muta- adult-onset primary tions.7,18,20 The existence of 2 distinct open-angle glaucoma GLC1C haplotypes suggests that muta- (POAG) a significant public health prob- tions will not be limited to rare descen- H1 lem.
    [Show full text]
  • Uncovering Cancer Gene Regulation by Accurate Regulatory Network Inference from Uninformative Data
    www.nature.com/npjsba ARTICLE OPEN Uncovering cancer gene regulation by accurate regulatory network inference from uninformative data Deniz Seçilmiş 1, Thomas Hillerton1, Daniel Morgan 1, Andreas Tjärnberg 2, Sven Nelander3, Torbjörn E. M. Nordling 4 and ✉ Erik L. L. Sonnhammer 1 The interactions among the components of a living cell that constitute the gene regulatory network (GRN) can be inferred from perturbation-based gene expression data. Such networks are useful for providing mechanistic insights of a biological system. In order to explore the feasibility and quality of GRN inference at a large scale, we used the L1000 data where ~1000 genes have been perturbed and their expression levels have been quantified in 9 cancer cell lines. We found that these datasets have a very low signal-to-noise ratio (SNR) level causing them to be too uninformative to infer accurate GRNs. We developed a gene reduction pipeline in which we eliminate uninformative genes from the system using a selection criterion based on SNR, until reaching an informative subset. The results show that our pipeline can identify an informative subset in an overall uninformative dataset, allowing inference of accurate subset GRNs. The accurate GRNs were functionally characterized and potential novel cancer-related regulatory interactions were identified. npj Systems Biology and Applications (2020) 6:37 ; https://doi.org/10.1038/s41540-020-00154-6 1234567890():,; INTRODUCTION where the main aim is to improve the SNR of the dataset by Living organisms are orchestrated by the biochemical reactions permanently removing the least informative genes and their that occur as a result of the interactions between biomolecules.
    [Show full text]
  • The Human Gene Connectome As a Map of Short Cuts for Morbid Allele Discovery
    The human gene connectome as a map of short cuts for morbid allele discovery Yuval Itana,1, Shen-Ying Zhanga,b, Guillaume Vogta,b, Avinash Abhyankara, Melina Hermana, Patrick Nitschkec, Dror Friedd, Lluis Quintana-Murcie, Laurent Abela,b, and Jean-Laurent Casanovaa,b,f aSt. Giles Laboratory of Human Genetics of Infectious Diseases, Rockefeller Branch, The Rockefeller University, New York, NY 10065; bLaboratory of Human Genetics of Infectious Diseases, Necker Branch, Paris Descartes University, Institut National de la Santé et de la Recherche Médicale U980, Necker Medical School, 75015 Paris, France; cPlateforme Bioinformatique, Université Paris Descartes, 75116 Paris, France; dDepartment of Computer Science, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel; eUnit of Human Evolutionary Genetics, Centre National de la Recherche Scientifique, Unité de Recherche Associée 3012, Institut Pasteur, F-75015 Paris, France; and fPediatric Immunology-Hematology Unit, Necker Hospital for Sick Children, 75015 Paris, France Edited* by Bruce Beutler, University of Texas Southwestern Medical Center, Dallas, TX, and approved February 15, 2013 (received for review October 19, 2012) High-throughput genomic data reveal thousands of gene variants to detect a single mutated gene, with the other polymorphic genes per patient, and it is often difficult to determine which of these being of less interest. This goes some way to explaining why, variants underlies disease in a given individual. However, at the despite the abundance of NGS data, the discovery of disease- population level, there may be some degree of phenotypic homo- causing alleles from such data remains somewhat limited. geneity, with alterations of specific physiological pathways under- We developed the human gene connectome (HGC) to over- come this problem.
    [Show full text]
  • A Network View of Microrna and Gene Interactions in Different Pathological Stages of Colon Cancer Jia Wen, Benika Hall and Xinghua Shi*
    Wen et al. BMC Medical Genomics 2019, 12(Suppl 7):158 https://doi.org/10.1186/s12920-019-0597-1 Research Open Access A network view of microRNA and gene interactions in different pathological stages of colon cancer Jia Wen, Benika Hall and Xinghua Shi* From 14th International Symposium on Bioinformatics Research and Applications (ISBRA’18) Beijing, China. 8–11 June 2018 Abstract Background: Colon cancer is one of the common cancers in human. Although the number of annual cases has decreased drastically, prognostic screening and translational methods can be improved. Hence, it is critical to understand the molecular mechanisms of disease progression and prognosis. Results: In this study, we develop a new strategy for integrating microRNA and gene expression profiles together with clinical information toward understanding the regulation of colon cancer. Particularly, we use this approach to identify microRNA and gene expression networks that are specific to certain pathological stages. To demonstrate the application of our method, we apply this approach to identify microRNA and gene interactions that are specific to pathological stages of colon cancer in The Cancer Genome Atlas (TCGA) datasets. Conclusions: Our results show that there are significant differences in network connections between miRNAs and genes in different pathological stages of colon cancer. These findings point to a hypothesis that these networks signify different roles of microRNA and gene regulation in the pathogenesis and tumorigenesis of colon cancer. Background survival rate is significantly decreased to approximately Colon cancer, which is reported to be one of the few cur- 12%. The drastic decrease in survival rate in colon can- able cancers, is one of the most common cancers around cer speaks to the need for better early diagnostic and the world.
    [Show full text]
  • Profiling and Validation of the Circular RNA Repertoire in Adult Murine Hearts
    Accepted Manuscript Profiling and Validation of the Circular RNA Repertoire in Adult Murine Hearts Tobias Jakobi, Lisa F. Czaja-Hasse, Richard Reinhardt, Christoph Dieterich PII: S1672-0229(16)30033-X DOI: http://dx.doi.org/10.1016/j.gpb.2016.02.003 Reference: GPB 200 To appear in: Genomics, Proteomics & Bioinformatics Received Date: 20 November 2015 Revised Date: 19 January 2016 Accepted Date: 15 February 2016 Please cite this article as: T. Jakobi, L.F. Czaja-Hasse, R. Reinhardt, C. Dieterich, Profiling and Validation of the Circular RNA Repertoire in Adult Murine Hearts, Genomics, Proteomics & Bioinformatics (2016), doi: http:// dx.doi.org/10.1016/j.gpb.2016.02.003 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. Profiling and Validation of the Circular RNA Repertoire in Adult Murine Hearts Tobias Jakobi1,2,*,a, Lisa F. Czaja-Hasse3,b , Richard Reinhardt3,c, Christoph Dieterich1,2,d 1Section of Bioinformatics and Systems Cardiology, Department of Internal Medicine III, University Hospital Heidelberg, Heidelberg D-69120, Germany 2 German Centre for Cardiovascular Research (DZHK), partner site Heidelberg / Mannheim, Heidelberg D-69120, Germany 3Max Planck-Genome-Centre Cologne, Cologne D-50829, Germany * Corresponding author. E-mail: [email protected] (Jakobi T).
    [Show full text]