Epigenomics and Single-Cell Sequencing Define a Developmental Hierarchy in Langerhans Cell Histiocytosis

Total Page:16

File Type:pdf, Size:1020Kb

Epigenomics and Single-Cell Sequencing Define a Developmental Hierarchy in Langerhans Cell Histiocytosis Published OnlineFirst July 25, 2019; DOI: 10.1158/2159-8290.CD-19-0138 RESEARCH ARTICLE Epigenomics and Single-Cell Sequencing Define a Developmental Hierarchy in Langerhans Cell Histiocytosis Florian Halbritter1,2, Matthias Farlik1,3, Raphaela Schwentner2, Gunhild Jug2, Nikolaus Fortelny1, Thomas Schnöller2, Hanja Pisa2, Linda C. Schuster1, Andrea Reinprecht4, Thomas Czech4, Johannes Gojo5, Wolfgang Holter2,5,6, Milen Minkov2,7, Wolfgang M. Bauer3, Ingrid Simonitsch-Klupp8, Christoph Bock1,9,10,11, and Caroline Hutter2,5,6 ABSTRACT Langerhans cell histiocytosis (LCH) is a rare neoplasm predominantly affecting children. It occupies a hybrid position between cancers and inflammatory dis- eases, which makes it an attractive model for studying cancer development. To explore the molecular mechanisms underlying the pathophysiology of LCH and its characteristic clinical heterogeneity, we investigated the transcriptomic and epigenomic diversity in primary LCH lesions. Using single-cell RNA sequencing, we identified multiple recurrent types of LCH cells within these biopsies, including puta- tive LCH progenitor cells and several subsets of differentiated LCH cells. We confirmed the presence of proliferative LCH cells in all analyzed biopsies using IHC, and we defined an epigenomic and gene- regulatory basis of the different LCH-cell subsets by chromatin-accessibility profiling. In summary, our single-cell analysis of LCH uncovered an unexpected degree of cellular, transcriptomic, and epigenomic heterogeneity among LCH cells, indicative of complex developmental hierarchies in LCH lesions. SIGNIFICANCE: This study sketches a molecular portrait of LCH lesions by combining single-cell tran- scriptomics with epigenome profiling. We uncovered extensive cellular heterogeneity, explained in part by an intrinsic developmental hierarchy of LCH cells. Our findings provide new insights and hypotheses for advancing LCH research and a starting point for personalizing therapy. See related commentary by Gruber et al., p. 1343. 1CeMM Research Center for Molecular Medicine of the Austrian Acad- Note: Supplementary data for this article are available at Cancer Discovery emy of Sciences, Vienna, Austria. 2St. Anna Children’s Cancer Research Online (http://cancerdiscovery.aacrjournals.org/). 3 Institute (CCRI), Vienna, Austria. Department of Dermatology, Medical F. Halbritter and M. Farlik contributed equally to this work. University of Vienna, Vienna, Austria. 4Department of Neurosurgery, Medi- cal University of Vienna, Vienna, Austria. 5Department of Pediatrics, Medi- Corresponding Authors: Caroline Hutter, St. Anna Children’s Cancer cal University of Vienna, Vienna, Austria. 6St. Anna Children’s Hospital, Research Institute (CCRI), St. Anna Kinderkrebsforschung, Zimmerman- St. Anna Kinderspital, Vienna, Austria. 7Department of Pediatrics, Adoles- nplatz 10, 1090 Vienna, Austria. Phone: 4314-0470-4043; E-mail: caro- cent Medicine and Neonatology, Rudolfstiftung Hospital, Vienna, Austria. [email protected]; and Christoph Bock, CeMM Research Center for 8Clinical Institute of Pathology, Medical University of Vienna, Vienna, Aus- Molecular Medicine of the Austrian Academy of Sciences, Lazarettgasse tria. 9Department of Laboratory Medicine, Medical University of Vienna, 14, 1090 Vienna, Austria. E-mail: [email protected] Vienna, Austria. 10Max Planck Institute for Informatics, Saarland Informat- Cancer Discov 2019;9:1–16 11 ics Campus, Saarbrücken, Germany. Ludwig Boltzmann Institute for Rare doi: 10.1158/2159-8290.CD-19-0138 and Undiagnosed Diseases, Vienna, Austria. ©2019 American Association for Cancer Research. OF1 | CANCER DISCOVERY OCTOBER 2019 www.aacrjournals.org Downloaded from cancerdiscovery.aacrjournals.org on September 25, 2021. © 2019 American Association for Cancer Research. Published OnlineFirst July 25, 2019; DOI: 10.1158/2159-8290.CD-19-0138 INTRODUCTION BRAFV600E mutation (6–8), but the pathophysiology of LCH is different from cancers with oncogenic BRAFV600E mutations, Langerhans cell histiocytosis (LCH) is a rare hematopoietic such as melanoma, non–small cell lung cancer, and colorec- neoplasm driven by gain-of-function mutations in the MAPK tal cancer. LCH cells carry few genetic alterations beyond pathway (1, 2). The name of the disease originates from the BRAFV600E or another oncogenic driver in the MAPK pathway characteristic expression of CD207, also known as langerin, (MAP2K1, ARAF), and they show no progressive accumula- which is conventionally associated with Langerhans cells—a tion of recurrent somatic mutations (9). Instead, high levels type of resident epidermal immune cell (3). However, the rela- of inflammatory infiltrate suggest that cell-extrinsic, immu- tionship between LCH cells and Langerhans cells remains con- nologic stimuli could be key contributors to the formation of troversial, and LCH has also been linked to dendritic cells and LCH lesions (10–12). LCH is lethal in up to 20% of patients other myeloid cell types (4–6). LCH is phenotypically character- with risk organ–positive, multisystem disease treated accord- ized by an accumulation of CD1A and CD207 double-positive ing to current standard of care. In contrast, the disease often cells in various tissues (Fig. 1A). These lesions can develop in resolves spontaneously following unknown mechanisms in almost any organ, but are most common in skin and bone (1). patients with single lesions (1, 12). LCH severity is clinically assessed by the extent of the dis- ease (involvement of one or more organs, i.e., single-system RESULTS vs. multisystem disease) and the involvement of known risk organs (liver, spleen, hematopoietic system). Treatment deci- Langerhans Cell Histiocytosis Lesions Are sions are currently based solely on clinical presentation, as Composed of Distinct Cell Populations no reliable molecular markers for risk stratification have LCH cells within the same lesion show wide variability in been identified yet. LCH is most commonly driven by the their CD1A and CD207 levels based on IHC (Fig. 1B) and flow OCTOber 2019 CANCER DISCOVERY | OF2 Downloaded from cancerdiscovery.aacrjournals.org on September 25, 2021. © 2019 American Association for Cancer Research. Published OnlineFirst July 25, 2019; DOI: 10.1158/2159-8290.CD-19-0138 RESEARCH ARTICLE Halbritter et al. A Langerhans cells histiocytosis (LCH) is a rare BELCH cells are Marker genes identify multiple cell types in single-cell transcriptome data of LCH lesions paediatric disease that affects various organs heterogeneous CD1A CD207 MAP2K1 CD33 BRAF Lesions formed by immune Bone lesion within biopsies Brain cells of unknown origin CD1A Cell surface markers: Lymph node CD1A, CD207 High Skin Genetic drivers: Lung Gain of MAPK signaling, CD3D IL32 FOXP3 CD19 CD79A V600E Liver 60% BRAF CD207 Spleen Incidence: Bone ∼4.6/million (age 0–15) Hematoxylin and eosin stain Skin lesion t-SNE dimension 2 IL3RACLEC4C TCF4 CD14 CD163 BRAFV600E Normalized UMI count Low t-SNE dimension 1 CDLCH biopsy samples analyzed using LCH biopsy samples analyzed using F Single-cell transcriptomics captures LCH cells single-cell transcriptome sequencing single-cell transcriptome sequencing and their microenvironment Identifier Biopsy Sex Age Extent Plasmacytoid CD1A+ CD207 + LCH_ELymph FemaleChild Multi- dendritic cells (LCH cells) node system LCH_H Skin Female Infant Multi- B cells system LCH_I Bone Male Child Single- system LCH_J Skin Female Aged Multi- system LCH_K Bone Male Child Single- system t-SNE dimension 2 LCH_N Skin Male Adult Multi- t-SNE dimension 2 system T cells LCH_O Bone Male Child Single- 19,044 single-cell Monocytes/ system transcriptomes macrophages t-SNE dimension 1 t-SNE dimension 1 Patient: EHIJKNO LCH cellImmune cell Figure 1. Single-cell transcriptome analysis captures cellular and molecular diversity of LCH lesions. A, Clinical presentation and incidence of LCH. Shown is an X-ray image of LCH bone lesions in the skull (top right), a photograph of extensive LCH skin lesions (bottom right), and a hematoxylin and eosin staining of an LCH biopsy (bottom left). B, Cellular heterogeneity in an LCH biopsy (lymph node), as revealed by immunostaining for CD1A, CD207, and an antibody that specifically detects mutated BRAFV600E protein. C, Overview of patient biopsy samples analyzed by single-cell RNA sequencing (RNA-seq) in this study. D, Low-dimensional projection (t-SNE plot) of the combined single-cell RNA-seq dataset across all analyzed biopsies, comprising a total of 19,044 single-cell transcriptome profiles. E, Molecular heterogeneity in LCH illustrated by low-dimensional projection (as in D) of all single-cell transcriptome profiles overlaid with the expression levels of selected marker genes (dark blue color indicates high expression levels). F, Cellular hetero- geneity in LCH illustrated by low-dimensional projection (as in D), annotated with cell types inferred from marker-gene expression. cytometry (Supplementary Fig. S1). To characterize the cel- Focusing initially on LCH-specific transcriptional profiles lular and molecular landscape of LCH lesions, we performed that were shared across cells and across patients, we com- droplet-based single-cell transcriptome sequencing of seven pared CD1A and CD207 marker–positive LCH cells with four LCH lesions. The biopsies were obtained from patients with immune-cell populations identified in all biopsies (Supple- multisystem disease (n = 4) and single-system disease (n = 3), mentary Fig. S2D). The LCH cells showed high expression of and they were collected
Recommended publications
  • Single Amino Acid Repeats in the Proteome World: Structural, Functional, and Evolutionary Insights
    RESEARCH ARTICLE Single Amino Acid Repeats in the Proteome World: Structural, Functional, and Evolutionary Insights Amitha Sampath Kumar, Divya Tej Sowpati, Rakesh K. Mishra* Centre for Cellular and Molecular Biology, Council of Scientific and Industrial Research, Uppal Road, Hyderabad, 500007, India * [email protected] a11111 Abstract Microsatellites or simple sequence repeats (SSR) are abundant, highly diverse stretches of short DNA repeats present in all genomes. Tandem mono/tri/hexanucleotide repeats in the coding regions contribute to single amino acids repeats (SAARs) in the proteome. While SSRs in the coding region always result in amino acid repeats, a majority of SAARs arise due to a combination of various codons representing the same amino acid and not as a con- OPEN ACCESS sequence of SSR events. Certain amino acids are abundant in repeat regions indicating a Citation: Kumar AS, Sowpati DT, Mishra RK (2016) Single Amino Acid Repeats in the Proteome positive selection pressure behind the accumulation of SAARs. By analysing 22 proteomes World: Structural, Functional, and Evolutionary including the human proteome, we explored the functional and structural relationship of Insights. PLoS ONE 11(11): e0166854. amino acid repeats in an evolutionary context. Only ~15% of repeats are present in any doi:10.1371/journal.pone.0166854 known functional domain, while ~74% of repeats are present in the disordered regions, sug- Editor: Eugene A. Permyakov, Russian Academy of gesting that SAARs add to the functionality of proteins by providing flexibility, stability and Medical Sciences, RUSSIAN FEDERATION act as linker elements between domains. Comparison of SAAR containing proteins across Received: August 28, 2016 species reveals that while shorter repeats are conserved among orthologs, proteins with Accepted: November 5, 2016 longer repeats, >15 amino acids, are unique to the respective organism.
    [Show full text]
  • Supplementary Table 1: Adhesion Genes Data Set
    Supplementary Table 1: Adhesion genes data set PROBE Entrez Gene ID Celera Gene ID Gene_Symbol Gene_Name 160832 1 hCG201364.3 A1BG alpha-1-B glycoprotein 223658 1 hCG201364.3 A1BG alpha-1-B glycoprotein 212988 102 hCG40040.3 ADAM10 ADAM metallopeptidase domain 10 133411 4185 hCG28232.2 ADAM11 ADAM metallopeptidase domain 11 110695 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 195222 8038 hCG40937.4 ADAM12 ADAM metallopeptidase domain 12 (meltrin alpha) 165344 8751 hCG20021.3 ADAM15 ADAM metallopeptidase domain 15 (metargidin) 189065 6868 null ADAM17 ADAM metallopeptidase domain 17 (tumor necrosis factor, alpha, converting enzyme) 108119 8728 hCG15398.4 ADAM19 ADAM metallopeptidase domain 19 (meltrin beta) 117763 8748 hCG20675.3 ADAM20 ADAM metallopeptidase domain 20 126448 8747 hCG1785634.2 ADAM21 ADAM metallopeptidase domain 21 208981 8747 hCG1785634.2|hCG2042897 ADAM21 ADAM metallopeptidase domain 21 180903 53616 hCG17212.4 ADAM22 ADAM metallopeptidase domain 22 177272 8745 hCG1811623.1 ADAM23 ADAM metallopeptidase domain 23 102384 10863 hCG1818505.1 ADAM28 ADAM metallopeptidase domain 28 119968 11086 hCG1786734.2 ADAM29 ADAM metallopeptidase domain 29 205542 11085 hCG1997196.1 ADAM30 ADAM metallopeptidase domain 30 148417 80332 hCG39255.4 ADAM33 ADAM metallopeptidase domain 33 140492 8756 hCG1789002.2 ADAM7 ADAM metallopeptidase domain 7 122603 101 hCG1816947.1 ADAM8 ADAM metallopeptidase domain 8 183965 8754 hCG1996391 ADAM9 ADAM metallopeptidase domain 9 (meltrin gamma) 129974 27299 hCG15447.3 ADAMDEC1 ADAM-like,
    [Show full text]
  • Regulation of Cancer Stemness in Breast Ductal Carcinoma in Situ by Vitamin D Compounds
    Author Manuscript Published OnlineFirst on May 28, 2020; DOI: 10.1158/1940-6207.CAPR-19-0566 Author manuscripts have been peer reviewed and accepted for publication but have not yet been edited. Analysis of the Transcriptome: Regulation of Cancer Stemness in Breast Ductal Carcinoma In Situ by Vitamin D Compounds Naing Lin Shan1, Audrey Minden1,5, Philip Furmanski1,5, Min Ji Bak1, Li Cai2,5, Roman Wernyj1, Davit Sargsyan3, David Cheng3, Renyi Wu3, Hsiao-Chen D. Kuo3, Shanyi N. Li3, Mingzhu Fang4, Hubert Maehr1, Ah-Ng Kong3,5, Nanjoo Suh1,5 1Department of Chemical Biology, Ernest Mario School of Pharmacy; 2Department of Biomedical Engineering, School of Engineering; 3Department of Pharmaceutics, Ernest Mario School of Pharmacy; 4Environmental and Occupational Health Sciences Institute and School of Public Health, 5Rutgers Cancer Institute of New Jersey, New Brunswick; Rutgers, The State University of New Jersey, NJ, USA Running title: Regulation of cancer stemness by vitamin D compounds Key words: Breast cancer, cancer stemness, gene expression, DCIS, vitamin D compounds Financial Support: This research was supported by the National Institutes of Health grant R01 AT007036, R01 AT009152, ES005022, Charles and Johanna Busch Memorial Fund at Rutgers University and the New Jersey Health Foundation. Corresponding author: Dr. Nanjoo Suh, Department of Chemical Biology, Ernest Mario School of Pharmacy, Rutgers, The State University of New Jersey, 164 Frelinghuysen Road, Piscataway, New Jersey 08854. Tel: 848-445-8030, Fax: 732-445-0687; e-mail: [email protected] Disclosure of Conflict of Interest: “The authors declare no potential conflicts of interest” 1 Downloaded from cancerpreventionresearch.aacrjournals.org on October 1, 2021.
    [Show full text]
  • Content Based Search in Gene Expression Databases and a Meta-Analysis of Host Responses to Infection
    Content Based Search in Gene Expression Databases and a Meta-analysis of Host Responses to Infection A Thesis Submitted to the Faculty of Drexel University by Francis X. Bell in partial fulfillment of the requirements for the degree of Doctor of Philosophy November 2015 c Copyright 2015 Francis X. Bell. All Rights Reserved. ii Acknowledgments I would like to acknowledge and thank my advisor, Dr. Ahmet Sacan. Without his advice, support, and patience I would not have been able to accomplish all that I have. I would also like to thank my committee members and the Biomed Faculty that have guided me. I would like to give a special thanks for the members of the bioinformatics lab, in particular the members of the Sacan lab: Rehman Qureshi, Daisy Heng Yang, April Chunyu Zhao, and Yiqian Zhou. Thank you for creating a pleasant and friendly environment in the lab. I give the members of my family my sincerest gratitude for all that they have done for me. I cannot begin to repay my parents for their sacrifices. I am eternally grateful for everything they have done. The support of my sisters and their encouragement gave me the strength to persevere to the end. iii Table of Contents LIST OF TABLES.......................................................................... vii LIST OF FIGURES ........................................................................ xiv ABSTRACT ................................................................................ xvii 1. A BRIEF INTRODUCTION TO GENE EXPRESSION............................. 1 1.1 Central Dogma of Molecular Biology........................................... 1 1.1.1 Basic Transfers .......................................................... 1 1.1.2 Uncommon Transfers ................................................... 3 1.2 Gene Expression ................................................................. 4 1.2.1 Estimating Gene Expression ............................................ 4 1.2.2 DNA Microarrays ......................................................
    [Show full text]
  • Short-Term RANKL Exposure Initiates a Neoplastic Transcriptional Program in the Basal Epithelium of the Murine Salivary Gland T
    Cytokine 123 (2019) 154745 Contents lists available at ScienceDirect Cytokine journal homepage: www.elsevier.com/locate/cytokine Short-term RANKL exposure initiates a neoplastic transcriptional program in the basal epithelium of the murine salivary gland T Lan Haia,b,1, Maria M. Szwarca,1, David M. Lonarda, Kimal Rajapakshea, Dimuthu Pereraa, ⁎ Cristian Coarfaa, Michael Ittmannc, Rodrigo Fernandez-Valdiviad, John P. Lydona, a Department of Molecular and Cellular Biology, Baylor College of Medicine, Houston, TX, USA b Reproductive Medicine Center of Henan Provincial People’s Hospital, Zhengzhou, Henan Province, PR China c Department of Pathology, Dan L. Duncan Comprehensive Cancer Center, Baylor College of Medicine, Houston, TX, USA d Department of Pathology, Wayne State University School of Medicine Detroit, MI, USA ARTICLE INFO ABSTRACT Keywords: Although salivary gland cancers comprise only ∼3–6% of head and neck cancers, treatment options for patients Mouse with advanced-stage disease are limited. Because of their rarity, salivary gland malignancies are understudied RANKL compared to other exocrine tissue cancers. The comparative lack of progress in this cancer field is particularly Salivary gland evident when it comes to our incomplete understanding of the key molecular signals that are causal for the Tumor development and/or progression of salivary gland cancers. Using a novel conditional transgenic mouse Transcriptome (K5:RANKL), we demonstrate that Receptor Activator of NFkB Ligand (RANKL) targeted to cytokeratin 5-posi- RNA-seq tive basal epithelial cells of the salivary gland causes aggressive tumorigenesis within a short period of RANKL exposure. Genome-wide transcriptomic analysis reveals that RANKL markedly increases the expression levels of numerous gene families involved in cellular proliferation, migration, and intra- and extra-tumoral commu- nication.
    [Show full text]
  • High Expression of the Sd Synthase B4GALNT2 Associates with Good
    cells Article High Expression of the Sda Synthase B4GALNT2 Associates with Good Prognosis and Attenuates Stemness in Colon Cancer Michela Pucci y, Inês Gomes Ferreira y, Martina Orlandani, Nadia Malagolini, Manuela Ferracin and Fabio Dall’Olio * Department of Experimental, Diagnostic and Specialty Medicine (DIMES), General Pathology Building, University of Bologna, Via San Giacomo 14, Via San Giacomo 14, 40126 Bologna, Italy; [email protected] (M.P.); [email protected] (I.G.F.); [email protected] (M.O.); [email protected] (N.M.); [email protected] (M.F.) * Correspondence: [email protected]; Tel.: +39-051-2094704 These authors have equal contribution. y Received: 25 February 2020; Accepted: 7 April 2020; Published: 11 April 2020 Abstract: Background: The carbohydrate antigen Sda and its biosynthetic enzyme B4GALNT2 are highly expressed in normal colonic mucosa but are down-regulated to a variable degree in colon cancer tissues. Here, we investigated the clinical and biological importance of B4GALNT2 in colon cancer. Methods: Correlations of B4GALNT2 mRNA with clinical data were obtained from The Cancer Genome Atlas (TCGA) database; the phenotypic and transcriptomic changes induced by B4GALNT2 were studied in LS174T cells transfected with B4GALNT2 cDNA. Results: TCGA data indicate that patients with high B4GALNT2 expression in cancer tissues display longer survival than non-expressers. In LS174T cells, expression of B4GALNT2 did not affect the ability to heal a scratch wound or to form colonies in standard growth conditions but markedly reduced the growth in soft agar, the tridimensional (3D) growth as spheroids, and the number of cancer stem cells, indicating a specific effect of B4GALNT2 on the growth in poor adherence and stemness.
    [Show full text]
  • Distinct Transcription Profiles of Primary and Secondary Glioblastoma Subgroups
    Research Article Distinct Transcription Profiles of Primary and Secondary Glioblastoma Subgroups Cho-Lea Tso,1,2,6 William A. Freije,1 Allen Day,1 Zugen Chen,1 Barry Merriman,1 Ally Perlina,1 Yohan Lee,1 Ederlyn Q. Dia,3 Koji Yoshimoto,3 Paul S. Mischel,3,6 Linda M. Liau,4,6 Timothy F. Cloughesy,5,6 and Stanley F. Nelson1,6 Departments of 1Human Genetics, 2Medicine/Hematology-Oncology, 3Pathology and Laboratory Medicine, 4Neurosurgery, and 5Neurology, David Geffen School of Medicine, and 6Jonsson Comprehensive Cancer Center, University of California at Los Angeles, Los Angeles, California Abstract is possible that the malignant end points that are ultimately reached will prove to be shared in common by many types of tumors (5). The Glioblastomas are invasive and aggressive tumors of the brain, generally considered to arise from glial cells. A subset of these identification and functional assessment of genes altered in the cancers develops from lower-grade gliomas and can thus be process of malignant transformation is essential for understating clinically classified as ‘‘secondary,’’ whereas some glioblasto- the mechanism of cancer development and should facilitate the mas occur with no prior evidence of a lower-grade tumor and development of more effective treatments. can be clinically classified as ‘‘primary.’’ Substantial genetic Infiltrative astrocytic neoplasms are the most common brain differences between these groups of glioblastomas have been tumors of central nervous system in adults. Glioblastoma multi- identified previously. We used large-scale expression analyses forme (WHO grade IV) remains a devastating disease, with a to identify glioblastoma-associated genes (GAG) that are median survival of <1 year after diagnosis (6).
    [Show full text]
  • Transcriptomic Analysis for Prognostic Value in Head and Neck Squamous Cell Carcinoma
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 7 July 2021 Article Transcriptomic Analysis for Prognostic Value in Head and Neck Squamous Cell Carcinoma Li-Hsing Chi 1,2,3 , Alexander TH Wu 1, Michael Hsiao 4* and Yu-Chuan (Jack) Li 1,5* 1 The Ph.D. Program for Translational Medicine, College of Medical Science and Technology, Taipei Medical University and Academia Sinica, Taipei, Taiwan 2 Division of Oral and Maxillofacial Surgery, Department of Dentistry, Wan Fang Hospital, Taipei Medical University 3 Division of Oral and Maxillofacial Surgery, Department of Dentistry, Taipei Medical University Hospital, Taipei Medical University 4 Genomics Research Center, Academia Sinica, Taipei, Taiwan 5 Graduate Institute of Biomedical Informatics, College of Medical Science and Technology, Taipei Medical University, No.172-1, Sec. 2, Keelung Rd., Taipei 106, Taiwan * Correspondence: Hsiao: [email protected]; Li: [email protected] 1 Abstract: The survival analysis of the Cancer Genome Atlas (TCGA) dataset is a well-known 2 method to discover the gene expression-based prognostic biomarkers of head and neck squamous 3 cell carcinoma (HNSCC). A cutoff point is usually used in survival analysis for the patients’ 4 dichotomization in the continuous gene expression. There is some optimization software for cutoff 5 determination. However, the software’s predetermined cutoffs are usually set at the median or 6 quantiles of gene expression value to perform the analyses. There are also few clinicopathological 7 features available on their pre-processed data sets. We applied an in-house workflow, including 8 data retrieving and pre-processing, feature selection, sliding-window cutoff selection, Kaplan- 9 Meier survival analysis, and Cox proportional hazard modeling for biomarker discovery.
    [Show full text]
  • Downloaded SRA Files from Which Are Associated with Pigmentation (Table 1)
    Samaniego Castruita et al. BMC Genomics (2020) 21:543 https://doi.org/10.1186/s12864-020-06940-0 RESEARCH ARTICLE Open Access Analyses of key genes involved in Arctic adaptation in polar bears suggest selection on both standing variation and de novo mutations played an important role Jose Alfredo Samaniego Castruita†, Michael V. Westbury*† and Eline D. Lorenzen* Abstract Background: Polar bears are uniquely adapted to an Arctic existence. Since their relatively recent divergence from their closest living relative, brown bears, less than 500,000 years ago, the species has evolved an array of novel traits suited to its Arctic lifestyle. Previous studies sought to uncover the genomic underpinnings of these unique characteristics, and disclosed the genes showing the strongest signal of positive selection in the polar bear lineage. Here, we survey a comprehensive dataset of 109 polar bear and 33 brown bear genomes to investigate the genomic variants within these top genes present in each species. Specifically, we investigate whether fixed homozygous variants in polar bears derived from selection on standing variation in the ancestral gene pool or on de novo mutation in the polar bear lineage. Results: We find that a large number of sites fixed in polar bears are biallelic in brown bears, suggesting selection on standing variation. Moreover, we uncover sites in which polar bears are fixed for a derived allele while brown bears are fixed for the ancestral allele, which we suggest may be a signal of de novo mutation in the polar bear lineage. Conclusions: Our findings suggest that, among other mechanisms, natural selection acting on changes in genes derived from a combination of variation already in the ancestral gene pool, and from de novo missense mutations in the polar bear lineage, may have enabled the rapid adaptation of polar bears to their new Arctic environment.
    [Show full text]
  • A Transcriptomic Analysis of Head and Neck Squamous Cell Carcinomas for Prognostic Indications
    Journal of Personalized Medicine Article A Transcriptomic Analysis of Head and Neck Squamous Cell Carcinomas for Prognostic Indications Li-Hsing Chi 1,2,3 , Alexander T. H. Wu 1 , Michael Hsiao 4,5,* and Yu-Chuan (Jack) Li 1,6,* 1 The Ph.D. Program for Translational Medicine, College of Medical Science and Technology, Taipei Medical University and Academia Sinica, Taipei 11031, Taiwan; [email protected] (L.-H.C.); [email protected] (A.T.H.W.) 2 Division of Oral and Maxillofacial Surgery, Department of Dentistry, Wan Fang Hospital, Taipei Medical University, Taipei 11600, Taiwan 3 Division of Oral and Maxillofacial Surgery, Department of Dentistry, Taipei Medical University Hospital, Taipei Medical University, Taipei 11031, Taiwan 4 Genomics Research Center, Academia Sinica, Taipei 115024, Taiwan 5 Department of Biochemistry, College of Medicine, Kaohsiung Medical University, Kaohsiung 807378, Taiwan 6 Graduate Institute of Biomedical Informatics, College of Medical Science and Technology, Taipei Medical University, No.172-1, Sec. 2, Keelung Rd., Taipei 106339, Taiwan * Correspondence: [email protected] (M.H.); [email protected] (Y.-C.L.) Abstract: Survival analysis of the Cancer Genome Atlas (TCGA) dataset is a well-known method for discovering gene expression-based prognostic biomarkers of head and neck squamous cell carcinoma (HNSCC). A cutoff point is usually used in survival analysis for patient dichotomization when using continuous gene expression values. There is some optimization software for cutoff determination. However, the software’s predetermined cutoffs are usually set at the medians or quantiles of gene expression values. There are also few clinicopathological features available in pre-processed datasets.
    [Show full text]
  • Genome-Wide DNA Methylation and Gene Expression Analyses In
    Translational Psychiatry www.nature.com/tp ARTICLE OPEN Genome-wide DNA methylation and gene expression analyses in monozygotic twins identify potential biomarkers of depression ✉ Weijing Wang1, Weilong Li2, Yili Wu1, Xiaocao Tian3, Haiping Duan3, Shuxia Li 4, Qihua Tan 4,5 and Dongfeng Zhang 1 © The Author(s) 2021 Depression is currently the leading cause of disability around the world. We conducted an epigenome-wide association study (EWAS) in a sample of 58 depression score-discordant monozygotic twin pairs, aiming to detect specific epigenetic variants potentially related to depression and further integrate with gene expression profile data. Association between the methylation level of each CpG site and depression score was tested by applying a linear mixed effect model. Weighted gene co-expression network analysis (WGCNA) was performed for gene expression data. The association of DNA methylation levels of 66 CpG sites with depression score reached the level of P <1×10−4. These top CpG sites were located at 34 genes, especially PTPRN2, HES5, GATA2, PRDM7, and KCNIP1. Many ontology enrichments were highlighted, including Notch signaling pathway, Huntington disease, p53 pathway by glucose deprivation, hedgehog signaling pathway, DNA binding, and nucleic acid metabolic process. We detected 19 differentially methylated regions (DMRs), some of which were located at GRIK2, DGKA, and NIPA2. While integrating with gene expression data, HELZ2, PTPRN2, GATA2, and ZNF624 were differentially expressed. In WGCNA, one specific module was positively correlated with depression score (r = 0.62, P = 0.002). Some common genes (including BMP2, PRDM7, KCNIP1, and GRIK2) and enrichment terms (including complement and coagulation cascades pathway, DNA binding, neuron fate specification, glial cell differentiation, and thyroid gland development) were both identified in methylation analysis and WGCNA.
    [Show full text]
  • Primepcr™Assay Validation Report
    PrimePCR™Assay Validation Report Gene Information Gene Name Fc fragment of IgG binding protein Gene Symbol FCGBP Organism Human Gene Summary Description Not Available Gene Aliases FC(GAMMA)BP RefSeq Accession No. NC_000019.9, NT_011109.16 UniGene ID Hs.111732 Ensembl Gene ID ENSG00000090920 Entrez Gene ID 8857 Assay Information Unique Assay ID qHsaCEP0041227 Assay Type Probe - Validation information is for the primer pair using SYBR® Green detection Detected Coding Transcript(s) ENST00000221347, ENST00000601713 Amplicon Context Sequence GACAGAAAGTTTCAGTCCTCCAGGGCTCCACGTGGCCTCCAGCCGCCTGGCAA GCACTCACGTAGGCATGGATGTTGCTGCAGAGAATGCTCAGGTTCCCACCACCC AGGCAGAGA Amplicon Length (bp) 86 Chromosome Location 19:40408254-40408369 Assay Design Exonic Purification Desalted Validation Results Efficiency (%) 96 R2 0.9993 cDNA Cq 24.94 cDNA Tm (Celsius) 86 gDNA Cq 23.06 Specificity (%) 100 Information to assist with data interpretation is provided at the end of this report. Page 1/4 PrimePCR™Assay Validation Report FCGBP, Human Amplification Plot Amplification of cDNA generated from 25 ng of universal reference RNA Melt Peak Melt curve analysis of above amplification Standard Curve Standard curve generated using 20 million copies of template diluted 10-fold to 20 copies Page 2/4 PrimePCR™Assay Validation Report Products used to generate validation data Real-Time PCR Instrument CFX384 Real-Time PCR Detection System Reverse Transcription Reagent iScript™ Advanced cDNA Synthesis Kit for RT-qPCR Real-Time PCR Supermix SsoAdvanced™ SYBR® Green Supermix Experimental Sample qPCR Human Reference Total RNA Data Interpretation Unique Assay ID This is a unique identifier that can be used to identify the assay in the literature and online. Detected Coding Transcript(s) This is a list of the Ensembl transcript ID(s) that this assay will detect.
    [Show full text]