Comprehensive Analysis of Mrna Expression Profiles in Head and Neck Cancer by Using Robust Rank Aggregation and Weighted Gene Coexpression Network Analysis

Total Page:16

File Type:pdf, Size:1020Kb

Comprehensive Analysis of Mrna Expression Profiles in Head and Neck Cancer by Using Robust Rank Aggregation and Weighted Gene Coexpression Network Analysis Hindawi BioMed Research International Volume 2020, Article ID 4908427, 21 pages https://doi.org/10.1155/2020/4908427 Research Article Comprehensive Analysis of mRNA Expression Profiles in Head and Neck Cancer by Using Robust Rank Aggregation and Weighted Gene Coexpression Network Analysis Zaizai Cao , Yinjie Ao , Yu Guo , and Shuihong Zhou Department of Otolaryngology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Zhejiang Province 310003, China Correspondence should be addressed to Shuihong Zhou; [email protected] Received 3 June 2020; Revised 2 November 2020; Accepted 23 November 2020; Published 10 December 2020 Academic Editor: Hassan Dariushnejad Copyright © 2020 Zaizai Cao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Background. Head and neck squamous cell cancer (HNSCC) is the sixth most common cancer in the world; its pathogenic mechanism remains to be further clarified. Methods. Robust rank aggregation (RRA) analysis was utilized to identify the metasignature dysregulated genes, which were then used for potential drug prediction. Weighted gene coexpression network analysis (WGCNA) was performed on all metasignature genes to find hub genes. DNA methylation analysis, GSEA, functional annotation, and immunocyte infiltration analysis were then performed on hub genes to investigate their potential role in HNSCC. Result. A total of 862 metasignature genes were identified, and 6 potential drugs were selected based on these genes. Based on the result of WGCNA, six hub genes (ITM2A, GALNTL1, FAM107A, MFAP4, PGM5, and OGN) were selected (GS > 0:1, MM > 0:75,GSp value < 0.05, and MM p value < 0.05). All six genes were downregulated in tumor tissue (FDR < 0:01) and were related to the clinical stage and prognosis of HNSCC in different degrees. Methylation analysis showed that the dysregulation of ITM2A, GALNTL1, FAM107A, and MFAP4 may be caused by hypermethylation. Moreover, the expression level of all 6 hub genes was positively associated with immune cell infiltration, and the result of GSEA showed that all hub genes may be involved in the process of immunoregulation. Conclusion. All identified hub genes could be potential biomarkers for HNSCC and provide a new insight into the diagnosis and treatment of head and neck tumors. 1. Introduction different screening criteria have a great impact on the research results. To solve this problem and obtain stable bio- Head and neck squamous cell carcinoma (HNSC) is the sixth markers, researchers proposed a novel rank aggregation most common cancer in the world [1]. Worldwide, more method: robust rank aggregation (RRA) [5], which has been than 300000 patients die of HNSC every year [2]. Although implemented as an R package (RobustRankAggreg) [5], to many treatments for HNSC such as surgery, chemotherapy, identify the overlapping genes among ranked gene lists [6], and radiotherapy have obtained some success, the 5-year sur- thus making the result more reliable. vival rate is still only 40-50% [3]. The chances of survival for WGCNA [7] is an effective method to find the clusters of patients with HNSCC depend largely on the initial stage of highly correlated genes and identify the hub genes of each cancer. Therefore, early detection and accurate diagnosis cluster [7]. This method has been widely applied in various are crucial for patients with HNSCC to receive treatment. biological contexts. In our study, a total of 24 independent In the past twenty years, with the application of microar- gene datasets were included in RRA analysis to identify ray and next-generation sequencing technologies, a great robust DEGs. We used these DEGs to predict the potential number of novel diagnostic or therapeutic biomarkers have small molecular drugs. The coexpression network was then been identified in HNSCC [4]. However, small samples in established by WGCNA to identify hub genes in these robust independent research, different platform technologies, and DEGs. The role of all hub genes in HNSCC was then 2 BioMed Research International 24 HNSCC GEO datasets Differential analysis (p value < 0.05) Differentially expressed genes in each GEO datasets RRA (adj p value < 0.05) Screen related small GO and KEGG Metasignature DEGs molecule drugs (CMap) enrichment analysis WGCNA in GSE65858 Key module : blue GS and MM filtering Explore the correlation Function annotation by 6 hub genes between hub genes and GSEA immunocyte infiltration Explore the correlation (i) Validation of the diagnostic role of hub genes. between methylation and (TCGA and GSE30784) expression of hub genes (ii) Validation of the correlation between hub genes and clinical characteristic (TCGA) (iii) III Explore the expression of hub genes in dysplasia tissue (GSE30784) 4 Explore the prognostic role of hub genes in HNSCC (TCGA) Figure 1: Simple flow chart of the entire study. validated by other independent databases. Furthermore, we MAL (FDR = 8:50e − 35), CRISP3 (FDR = 1:60e − 34), also utilized multiple online tools such as DiseaseMeth [8], CRNN (FDR = 4:32e − 33), and KRT4 (FDR = 1:24e − 30) MEXPRESS [9], and MethSurv [10] to evaluate the methyla- were the most significant downregulated genes. The chromo- tion level of hub genes. TIMER was used to assess the associ- somal locations and expression patterns of the top 100 DEGs ation between immune infiltration and hub genes. GSEA [11] are visualized in Figure 2. Chromosome 1 contains most analysis was applied to explore the biological functions of metasignature DEGs while X and Y chromosome contains these hub genes. To the best of our knowledge, this is the first no DEGs. It is clear that almost all displayed genes have the time to utilize RRA and WGCNA simultaneously for screen- same expression pattern in most of the independent studies, ing biomarkers of HNSCC. which indicates the reliability of the RRA analysis result. 2. Result 2.2. Functional Enrichment Analysis. We select the top 300 DEGs to perform GO and KEGG enrichment analyses. 2.1. Metasignature DEGs Identified by RRA Analysis. The Among the KEGG pathway database, we can find that these workflow of our study is shown in Figure 1, 24 independent DEGs were enriched in multiple cancer-related pathways like studies were used in RRA analysis, and a total of 466 upreg- focal adhesion, PI3K-Akt signal pathway, pathway in cancer, ulated genes and 396 downregulated genes were identified. small-cell lung cancer, transcriptional misregulation in can- The top 5 upregulated genes in tumor tissue were MMP1 cer, and chemical carcinogenesis (Figure 3(a)). Furthermore, (FDR = 4:77e − 53), MMP10 (FDR = 8:25e − 40), PTHLH in all terms of KEGG and GO, we found that these meta- (FDR = 1:48e − 38), MMP3 (FDR = 5:38e − 38), and SPP1 signature genes mostly involved in pathways associated with (FDR = 2:79e − 37) while TMPRSS11B (FDR = 1:53e − 36), the construction of ECM such as ECM receptor interaction, BioMed Research International 3 ISG15 PDPN MFAP2 IFI6 CLCA4 TGFBR3 SLC16A1 CRNN S100A7A DPT FMO2 LAMC2 KIF14 IL24 Heatmap APOL1 LAMB3 MMP11 8.49 RSAD2 FAM3B MAL FAP HOXD10 COL5A2 TPX2 MYO1B TGM3 C2orf54 CEACAM1 CEACAM5 CEACAM7 BST2 FUT3 KAT2B GPD1L LAMA3 FAM4107A FAM3D MGLL ABCA8 RBP1 HLF COL1A1 KRT16 KRT13 MFAP4 TMPRSS11B CDH3 IL8 GDPD3 CXCL10 PPL CXCL11 TNFRSF12A SPP1 ADH7 TDO2 0 HPGD RHCG TRIP13 TGFBI SPINK5 PAX9 SPINK7 COL4A2 PTK7 COL4A1 SCEL PLA2G7 CRISP3 POSTIN SH3BGRL2 DCBLD1 WDR66 KRT78 KRT4 FSCN1 ENDOU DFNA5 PTHLH INHBA FOXM1 CYP3A5 SERPINE1 ATP6V0A4 MMP13 –8.49 MMP12 MMP3 MMP1 LOXL2 EPHX2 SNAI2 SLURP1 SERPINH1 MAMDC2 TNC IFT3 PLAU PPP1R3C Figure 2: Circular visualization of expression patterns and chromosomal positions of top 100 DEGs. Red indicates gene upregulation, blue indicates downregulation, and white indicates genes that do not exist in a given dataset. extracellular matrix organization, collagen catabolic process, expression data of metasignature DEGs. The correlation collagen binding, collagen trimer, extracellular region, and coefficients between each module and each clinical trait were extracellular exosome (Figure 3). calculated, and it is clear that only the blue module and gray module were significantly associated with T grade of HNSCC 2.3. Screen the Candidate Small Molecule Drugs for HNSCC. (Figure 4(e)). Because genes in the gray module are not sig- According to our screen criteria, 6 small molecule drugs (repa- nificantly coexpressed with each other, we only chose the glinide ES = −0:848, thiostrepton ES = −0:863,levamisoleES blue module as a key module. A total of 102 genes were = −0:75, cortisone ES = −0:866, zimeldine ES = −0:784,and included in blue modules, and the result of enrichment anal- cyproterone ES = −0:742) were identified (Table 1). Their 2D ysis for these genes showed that the most significant GO and structures are visualized in Supplementary Fig 1. These poten- KEGG terms were related to cell metabolism, chemokine tial drugs can to some extent reverse the robust dysregulated activity, and transmembrane transport (Supplementary Fig genes in HNSCC, thus providing suggestions for us to develop 2). According to the value of GS and MM (GS > 0:1, MM > targeted drugs. 0:75,GSp value < 0.05, and MM p value < 0.05), 6 genes (ITM2A, GALNTL1, OGN, FAM107A, MFAP4, and PGM5), 2.4. Identification of Hub Genes in HNSCC Patients. To iden- which were also significantly correlated with each other tify the hub genes, we performed WGCNA on the GSE65858, (Figure
Recommended publications
  • Screening and Identification of Key Biomarkers in Clear Cell Renal Cell Carcinoma Based on Bioinformatics Analysis
    bioRxiv preprint doi: https://doi.org/10.1101/2020.12.21.423889; this version posted December 23, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Screening and identification of key biomarkers in clear cell renal cell carcinoma based on bioinformatics analysis Basavaraj Vastrad1, Chanabasayya Vastrad*2 , Iranna Kotturshetti 1. Department of Biochemistry, Basaveshwar College of Pharmacy, Gadag, Karnataka 582103, India. 2. Biostatistics and Bioinformatics, Chanabasava Nilaya, Bharthinagar, Dharwad 580001, Karanataka, India. 3. Department of Ayurveda, Rajiv Gandhi Education Society`s Ayurvedic Medical College, Ron, Karnataka 562209, India. * Chanabasayya Vastrad [email protected] Ph: +919480073398 Chanabasava Nilaya, Bharthinagar, Dharwad 580001 , Karanataka, India bioRxiv preprint doi: https://doi.org/10.1101/2020.12.21.423889; this version posted December 23, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Abstract Clear cell renal cell carcinoma (ccRCC) is one of the most common types of malignancy of the urinary system. The pathogenesis and effective diagnosis of ccRCC have become popular topics for research in the previous decade. In the current study, an integrated bioinformatics analysis was performed to identify core genes associated in ccRCC. An expression dataset (GSE105261) was downloaded from the Gene Expression Omnibus database, and included 26 ccRCC and 9 normal kideny samples. Assessment of the microarray dataset led to the recognition of differentially expressed genes (DEGs), which was subsequently used for pathway and gene ontology (GO) enrichment analysis.
    [Show full text]
  • Anti-FAM107A Antibody
    D122337 Anti-FAM107A antibody Cat. No. D122337 Package 25 μl/100 μl/200 μl Storage -20°C, pH7.4 PBS, 0.05% NaN3, 40% Glycerol Product overview Description A n ti-FAM107A rabbit polyclonal antibody Applications ELISA, IHC Immunogen Fusion protein of human FAM107A Reactivity Human Content 0 .7 mg/ml Host species Rabbit Ig class Immunogen-specific rabbit IgG Purification Antigen affinity purification Target information Symbol FAM107A Full name family with sequence similarity 107, member A Synonyms DRR1; TU3A Swissprot O95990 Target Background FAM107B is a 131 amino acid protein that is encoded by a gene that maps to human chromosome 10, which contains over 800 genes and 135 million nucleotides, making up nearly 4.5% of the human genome. PTEN is an important tumor suppressor gene located on chromosome 10 and, when defective, causes a genetic predisposition to cancer development known as Cowden syndrome. The chromosome 10 encoded gene ERCC6 is important for DNA repair and is linked to Cockayne syndrome which is characterized by extreme photosensitivity and premature aging. Tetrahydrobiopterin deficiency and a number of syndromes involving defective skull and facial bone fusion are also linked to chromosoSme 10. angonAs with most trisomies, trisomy 10 is rarBie and is deloeteriouts. ech Applications Immunohistochemistry Predicted cell location: Nucleus and Cytoplasm Predicted cell location: Nucleus and Cytoplasm Positive control: Human liver cancer Positive control: Human colon cancer Recommended dilution: 50-200 Recommended dilution: 50-200 The image on the left is immunohistochemistry of paraffin-embedded The image on the left is immunohistochemistry of paraffin-embedded Human liver cancer tissue using D122337(FAM107A Antibody) at Human colon cancer tissue using D122337(FAM107A Antibody) at dilution 1/30, on the right is treated with fusion protein.
    [Show full text]
  • In Cancer Cells by Heat Shock Stimulation
    583-593.qxd 19/7/2010 11:13 Ì ™ÂÏ›‰·583 INTERNATIONAL JOURNAL OF ONCOLOGY 37: 583-593, 2010 583 Induction of HITS, a newly identified family with sequence similarity 107 protein (FAM107B), in cancer cells by heat shock stimulation HIDEO NAKAJIMA1, YASUHITO ISHIGAKI2, QI-SHENG XIA1, TAKAYUKI IKEDA3, YOSHINO YOSHITAKE3, HIDETO YONEKURA3, TAKAYUKI NOJIMA4, TAKUJI TANAKA4, HISANORI UMEHARA5, NAOHISA TOMOSUGI2, TAKANOBU TAKATA2, TAKEO SHIMASAKI1, NAOKI NAKAYA1, ITARU SATO1, KAZUYUKI KAWAKAMI6, KEITA KOIZUMI7, TOSHINARI MINAMOTO6 and YOSHIHARU MOTOO1 1Department of Medical Oncology, 2Medical Research Institute and Departments of 3Biochemistry, 4Pathology and 5Hematology and Immunology, Kanazawa Medical University, 1-1 Daigaku, Uchinada, Kahoku, Ishikawa 920-0293; 6Division of Translational and Clinical Oncology, Cancer Research Institute and 7Department of Biophysical Genetics, Graduate School of Medical Science, Kanazawa University, 13-1 Takara-machi, Kanazawa, Ishikawa 920-0934, Japan Received March 31, 2010; Accepted May 14, 2010 DOI: 10.3892/ijo_00000707 Abstract. The Family with sequence similarity 107 (FAM107) Introduction possesses an N-terminal domain of unknown function (DUF1151) that is highly conserved beyond species. In human, The Family with sequence similarity 107 (FAM107) possesses FAM107A termed TU3A/DRR1 has been reported as a can- an N-terminal domain of unknown function (DUF1151) that didate tumor suppressor gene which expression is down- is conserved beyond species including mammalian, Xenopus, regulated in several types of cancer, however no studies have fish and Drosophila, and shows no homology match to investigated the other family protein, FAM107B. In the other functional conserved domains (http://www.ncbi.nlm. present study, we designated FAM107B as heat shock- nih.gov/structure/cdd/cdd.shtml).
    [Show full text]
  • Supplementary Data
    Supplemental figures Supplemental figure 1: Tumor sample selection. A total of 98 thymic tumor specimens were stored in Memorial Sloan-Kettering Cancer Center tumor banks during the study period. 64 cases corresponded to previously untreated tumors, which were resected upfront after diagnosis. Adjuvant treatment was delivered in 7 patients (radiotherapy in 4 cases, cyclophosphamide- doxorubicin-vincristine (CAV) chemotherapy in 3 cases). 34 tumors were resected after induction treatment, consisting of chemotherapy in 16 patients (cyclophosphamide-doxorubicin- cisplatin (CAP) in 11 cases, cisplatin-etoposide (PE) in 3 cases, cisplatin-etoposide-ifosfamide (VIP) in 1 case, and cisplatin-docetaxel in 1 case), in radiotherapy (45 Gy) in 1 patient, and in sequential chemoradiation (CAP followed by a 45 Gy-radiotherapy) in 1 patient. Among these 34 patients, 6 received adjuvant radiotherapy. 1 Supplemental Figure 2: Amino acid alignments of KIT H697 in the human protein and related orthologs, using (A) the Homologene database (exons 14 and 15), and (B) the UCSC Genome Browser database (exon 14). Residue H697 is highlighted with red boxes. Both alignments indicate that residue H697 is highly conserved. 2 Supplemental Figure 3: Direct comparison of the genomic profiles of thymic squamous cell carcinomas (n=7) and lung primary squamous cell carcinomas (n=6). (A) Unsupervised clustering analysis. Gains are indicated in red, and losses in green, by genomic position along the 22 chromosomes. (B) Genomic profiles and recurrent copy number alterations in thymic carcinomas and lung squamous cell carcinomas. Gains are indicated in red, and losses in blue. 3 Supplemental Methods Mutational profiling The exonic regions of interest (NCBI Human Genome Build 36.1) were broken into amplicons of 500 bp or less, and specific primers were designed using Primer 3 (on the World Wide Web for general users and for biologist programmers (see Supplemental Table 2) [1].
    [Show full text]
  • Rank Aggregation Via the Cross Entropy Algorithm
    Finding cancer genes through meta-analysis of of the individual microarray studies were hybridized with the A®ymetrix GeneChip Human Genome HG-U133A and record expression levels for 22,283 Probe IDs. To make our ¯nal results meaningful and comprehensive at the same time, we decided to focus on microarray experiments that study di®erent types of cancer in humans. The goal of our meta-analysis is to identify genetic factors which are common across di®erent types and stages of cancer. For that purpose, we have selected 20 di®erent cancer-related microarray experiments which are included in the contest meta-dataset and have explicit cell type groupings necessary for detecting di®erentially expressed genes. Here, we list the selected experiment IDs along with the number of samples in each experiment in the parenthesis: E-MEXP-72 (20), E-MEXP-83 (22), E-MEXP-76 (17), E-MEXP-97 (24), E-MEXP-121 (30), E-MEXP-149 (20), E-MEXP-231 (58), E-MEXP-353 (96), E-TABM-26 (57), E-MEXP-669 (24), GSE4475 (221), GSE1456 (159), GSE5090 (17), GSE1420 (24), GSE1577 (29), GSE1729 (43), GSE2485 (18), GSE2603 (21), GSE3585 (12), GSE4127 (29). The total number of selected arrays is 941 (about 1/6 of the overall number of samples in the meta-dataset). One can refer to ArrayExpress database which provides public access to the microarray data from these experiments (http://www.ebi.ac.uk/arrayexpress/ ). 3 Methodology The proposed meta-analysis approach to microarray data is a two-step procedure: 1. Individual Analysis. By analyzing each microarray dataset individually, a set of \interesting" genes (top-50 Probe IDs) that exhibit the largest di®erences in terms of expression values between the groups is obtained for each dataset.
    [Show full text]
  • Downloaded Per Proteome Cohort Via the Web- Site Links of Table 1, Also Providing Information on the Deposited Spectral Datasets
    www.nature.com/scientificreports OPEN Assessment of a complete and classifed platelet proteome from genome‑wide transcripts of human platelets and megakaryocytes covering platelet functions Jingnan Huang1,2*, Frauke Swieringa1,2,9, Fiorella A. Solari2,9, Isabella Provenzale1, Luigi Grassi3, Ilaria De Simone1, Constance C. F. M. J. Baaten1,4, Rachel Cavill5, Albert Sickmann2,6,7,9, Mattia Frontini3,8,9 & Johan W. M. Heemskerk1,9* Novel platelet and megakaryocyte transcriptome analysis allows prediction of the full or theoretical proteome of a representative human platelet. Here, we integrated the established platelet proteomes from six cohorts of healthy subjects, encompassing 5.2 k proteins, with two novel genome‑wide transcriptomes (57.8 k mRNAs). For 14.8 k protein‑coding transcripts, we assigned the proteins to 21 UniProt‑based classes, based on their preferential intracellular localization and presumed function. This classifed transcriptome‑proteome profle of platelets revealed: (i) Absence of 37.2 k genome‑ wide transcripts. (ii) High quantitative similarity of platelet and megakaryocyte transcriptomes (R = 0.75) for 14.8 k protein‑coding genes, but not for 3.8 k RNA genes or 1.9 k pseudogenes (R = 0.43–0.54), suggesting redistribution of mRNAs upon platelet shedding from megakaryocytes. (iii) Copy numbers of 3.5 k proteins that were restricted in size by the corresponding transcript levels (iv) Near complete coverage of identifed proteins in the relevant transcriptome (log2fpkm > 0.20) except for plasma‑derived secretory proteins, pointing to adhesion and uptake of such proteins. (v) Underrepresentation in the identifed proteome of nuclear‑related, membrane and signaling proteins, as well proteins with low‑level transcripts.
    [Show full text]
  • Somamer Reagents Generated to Human Proteins Number Somamer Seqid Analyte Name Uniprot ID 1 5227-60
    SOMAmer Reagents Generated to Human Proteins The exact content of any pre-specified menu offered by SomaLogic may be altered on an ongoing basis, including the addition of SOMAmer reagents as they are created, and the removal of others if deemed necessary, as we continue to improve the performance of the SOMAscan assay. However, the client will know the exact content at the time of study contracting. SomaLogic reserves the right to alter the menu at any time in its sole discretion. Number SOMAmer SeqID Analyte Name UniProt ID 1 5227-60 [Pyruvate dehydrogenase (acetyl-transferring)] kinase isozyme 1, mitochondrial Q15118 2 14156-33 14-3-3 protein beta/alpha P31946 3 14157-21 14-3-3 protein epsilon P62258 P31946, P62258, P61981, Q04917, 4 4179-57 14-3-3 protein family P27348, P63104, P31947 5 4829-43 14-3-3 protein sigma P31947 6 7625-27 14-3-3 protein theta P27348 7 5858-6 14-3-3 protein zeta/delta P63104 8 4995-16 15-hydroxyprostaglandin dehydrogenase [NAD(+)] P15428 9 4563-61 1-phosphatidylinositol 4,5-bisphosphate phosphodiesterase gamma-1 P19174 10 10361-25 2'-5'-oligoadenylate synthase 1 P00973 11 3898-5 26S proteasome non-ATPase regulatory subunit 7 P51665 12 5230-99 3-hydroxy-3-methylglutaryl-coenzyme A reductase P04035 13 4217-49 3-hydroxyacyl-CoA dehydrogenase type-2 Q99714 14 5861-78 3-hydroxyanthranilate 3,4-dioxygenase P46952 15 4693-72 3-hydroxyisobutyrate dehydrogenase, mitochondrial P31937 16 4460-8 3-phosphoinositide-dependent protein kinase 1 O15530 17 5026-66 40S ribosomal protein S3 P23396 18 5484-63 40S ribosomal protein
    [Show full text]
  • Exposure of Mice to Secondhand Smoke Elicits Both Transient and Long-Lasting Transcriptional Changes in Cancer-Related Functional Networks
    IJC International Journal of Cancer Exposure of mice to secondhand smoke elicits both transient and long-lasting transcriptional changes in cancer-related functional networks Stella Tommasi, Albert Zheng and Ahmad Besaratinia Department of Preventive Medicine, USC Keck School of Medicine, University of Southern California, Los Angeles, CA Secondhand smoke (SHS) has long been linked to lung cancer and other diseases in nonsmokers. Yet, the underlying mecha- nisms of SHS carcinogenicity in nonsmokers remain to be elucidated. We investigated the immediate and long-lasting effects of SHS exposure on gene expression in mice in vivo. We exposed mice whole body to SHS for 5 h/day, 5 days/week for 4 months in exposure chambers of a microprocessor-controlled smoking machine. Subsequently, we performed microarray gene expression profiling, genome-wide, to construct the pulmonary transcriptome of SHS-exposed mice, immediately after discon- tinuation of exposure (T0) and following 1-month (T1) and 7-month (T2) recoveries in clean air. Sub-chronic exposure of mice to SHS elicited a robust transcriptomic response, including both reversible and irreversible changes in gene expression. There were 674 differentially expressed transcripts immediately after treatment (T0), of which the majority were involved in xenobi- otic metabolism, signaling, and innate immune response. Reduced, yet, substantial numbers of differentially expressed tran- scripts were detectable in mice after cessation of SHS-exposure (254 transcripts at T1 and 30 transcripts at T2). Top biofunctional networks disrupted in SHS-exposed mice, even after termination of exposure, were implicated in cancer, respira- tory disease, and inflammatory disease. Our data show that exposure of mice to SHS induces both transient and long-lasting changes in gene expression, which impact cancer-related functional networks.
    [Show full text]
  • Family with Sequence Similarity 107: a Family of Stress Responsive Small Proteins with Diverse Functions in Cancer and the Nervous System (Review)
    BIOMEDICAL REPORTS 2: 321-325, 2014 Family with sequence similarity 107: A family of stress responsive small proteins with diverse functions in cancer and the nervous system (Review) HIDEO NAKAJIMA1,2 and KEITA KOIZUMI3 1Department of Oncology, Ageo Central General Hospital, Ageo, Saitama 362-8588; 2Center for AIDS Research, Kumamoto University, Kumamoto 860-0811; 3Research Center for Child Mental Development, Kanazawa University, Kanazawa, Ishikawa 920-8640, Japan Received January 30, 2014; Accepted February 04, 2014 DOI: 10.3892/br.2014.243 Abstract. Under conditions of acute stress, rapid adaptation is Introduction crucial for maximizing biological survival. The responses to environmental stress are often complex, involving numerous Family with sequence similarity 107 (FAM107) members genes and integrating events at the cellular and organismal contain an N-terminal domain of unknown function levels. The heat shock proteins (HSPs) are a family of highly (DUF1151) that is conserved across species, including conserved proteins that play critical roles in maintaining cell mammalian, Xenopus, fish and Drosophila, with no homolo- homeostasis and protecting cells under chronic and acute gous matches to other functionally conserved domains (Fig. 1). stress conditions. The genes for these stress-responding This family includes several hypothetical eukaryotic proteins proteins are widely distributed in organisms, tissues and cells. with largely undetermined functions. Mammals have two HSPs participate in a variety of physiological processes and genes, FAM107A and FAM107B, which encode for proteins are associated with various types of disease. In this review, we of 144 and 131 amino acids (aa), respectively. The C-terminal focused on family with sequence similarity 107 (FAM107), a variable regions of FAM107 members have a coiled-coil novel unique protein family that exhibits functional similarity domain that has been identified in several nuclear proteins, with HSPs during the cellular stress response.
    [Show full text]
  • Comprehensive Analysis Reveals Novel Gene Signature in Head and Neck Squamous Cell Carcinoma: Predicting Is Associated with Poor Prognosis in Patients
    5892 Original Article Comprehensive analysis reveals novel gene signature in head and neck squamous cell carcinoma: predicting is associated with poor prognosis in patients Yixin Sun1,2#, Quan Zhang1,2#, Lanlin Yao2#, Shuai Wang3, Zhiming Zhang1,2 1Department of Breast Surgery, The First Affiliated Hospital of Xiamen University, School of Medicine, Xiamen University, Xiamen, China; 2School of Medicine, Xiamen University, Xiamen, China; 3State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Xiamen University, Xiamen, China Contributions: (I) Conception and design: Y Sun, Q Zhang; (II) Administrative support: Z Zhang; (III) Provision of study materials or patients: Y Sun, Q Zhang; (IV) Collection and assembly of data: Y Sun, L Yao; (V) Data analysis and interpretation: Y Sun, S Wang; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors. #These authors contributed equally to this work. Correspondence to: Zhiming Zhang. Department of Surgery, The First Affiliated Hospital of Xiamen University, Xiamen, China. Email: [email protected]. Background: Head and neck squamous cell carcinoma (HNSC) remains an important public health problem, with classic risk factors being smoking and excessive alcohol consumption and usually has a poor prognosis. Therefore, it is important to explore the underlying mechanisms of tumorigenesis and screen the genes and pathways identified from such studies and their role in pathogenesis. The purpose of this study was to identify genes or signal pathways associated with the development of HNSC. Methods: In this study, we downloaded gene expression profiles of GSE53819 from the Gene Expression Omnibus (GEO) database, including 18 HNSC tissues and 18 normal tissues.
    [Show full text]
  • Transcriptome Sequencing Uncovers Novel Long Non-Coding and Small Nucleolar Rnas Dysregulated in Head and Neck Squamous Cell Carcinoma
    Transcriptome sequencing uncovers novel long non-coding and small nucleolar RNAs dysregulated in head and neck squamous cell carcinoma Angela E. Zou1, Jonjei Ku1, Thomas K. Honda1, Vicky Yu1, Selena Z. Kuo1, Hao Zheng1, Yinan Xuan1, Andrew Hinton2, Kevin T. Brumund1, Jonathan H. Lin3, Jessica Wang-Rodriguez3, Weg M. Ongkeko1* 1Division of Otolaryngology-Head and Neck Surgery, Department of Surgery, University of California, San Diego, La Jolla, California, United States of America 2Department of Pediatrics, University of California, San Diego, La Jolla, California, United States of America 3Veterans Administration Medical Center and Department of Pathology, University of California San Diego, La Jolla, California, United States of America * Corresponding author Shortened Title: Novel non-coding RNAs in head and neck cancer Keywords: HNSCC; Long non-coding RNAs; RNA-sequencing; Small nucleolar RNAs Zou 2 Abstract Head and neck squamous cell carcinoma persists as one of the most common and deadly malignancies, with early detection and effective treatment still posing formidable challenges. To expand our currently sparse knowledge of the non-coding alterations involved in the disease and identify potential biomarkers and therapeutic targets, we globally profiled the dysregulation of small nucleolar and long non-coding RNAs in head and neck tumors. Using next-generation RNA-sequencing data from 40 pairs of tumor and matched normal tissues, we found 2,808 long non-coding RNA (lncRNA) transcripts significantly differentially expressed by a fold change magnitude ≥ 2. Meanwhile, RNA-sequencing analysis of 31 tumor-normal pairs yielded 33 significantly dysregulated small nucleolar RNAs (snoRNA). In particular, we identified 2 dramatically downregulated lncRNAs and 1 downregulated snoRNA whose expression levels correlated significantly with overall patient survival, suggesting their functional significance and clinical relevance in head and neck cancer pathogenesis.
    [Show full text]
  • WO 2014/028907 Al 20 February 2014 (20.02.2014) P O P C T
    (12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) (19) World Intellectual Property Organization International Bureau (10) International Publication Number (43) International Publication Date WO 2014/028907 Al 20 February 2014 (20.02.2014) P O P C T (51) International Patent Classification: (74) Agents: EVANS, Judith et al; P.O. Box 23 1, Manassas, C12Q 1/68 (2006.01) G01N 33/00 (2006.01) VA 20108 (US). (21) International Application Number: (81) Designated States (unless otherwise indicated, for every PCT/US20 13/055469 kind of national protection available): AE, AG, AL, AM, AO, AT, AU, AZ, BA, BB, BG, BH, BN, BR, BW, BY, (22) International Filing Date: BZ, CA, CH, CL, CN, CO, CR, CU, CZ, DE, DK, DM, 16 August 2013 (16.08.2013) DO, DZ, EC, EE, EG, ES, FI, GB, GD, GE, GH, GM, GT, (25) Filing Language: English HN, HR, HU, ID, IL, IN, IS, JP, KE, KG, KN, KP, KR, KZ, LA, LC, LK, LR, LS, LT, LU, LY, MA, MD, ME, (26) Publication Language: English MG, MK, MN, MW, MX, MY, MZ, NA, NG, NI, NO, NZ, (30) Priority Data: OM, PA, PE, PG, PH, PL, PT, QA, RO, RS, RU, RW, SA, 61/684,029 16 August 2012 (16.08.2012) US SC, SD, SE, SG, SK, SL, SM, ST, SV, SY, TH, TJ, TM, 61/718,468 25 October 2012 (25. 10.2012) US TN, TR, TT, TZ, UA, UG, US, UZ, VC, VN, ZA, ZM, 61/745,207 2 1 December 2012 (21. 12.2012) US ZW. (71) Applicant: THE TRUSTEES OF COLUMBIA UNI¬ (84) Designated States (unless otherwise indicated, for every VERSITY IN THE CITY OF NEW YORK [US/US]; kind of regional protection available): ARIPO (BW, GH, 412 Low Memorial Library, 535 West 116th Street, New GM, KE, LR, LS, MW, MZ, NA, RW, SD, SL, SZ, TZ, York, NY 10027 (US).
    [Show full text]