Performance of a Genomic Sequencing Classifier for the Preoperative Diagnosis of Cytologically Indeterminate Thyroid Nodules

Total Page:16

File Type:pdf, Size:1020Kb

Performance of a Genomic Sequencing Classifier for the Preoperative Diagnosis of Cytologically Indeterminate Thyroid Nodules Supplementary Online Content Patel KN, Angell TE, Babiarz J, et al. Performance of a genomic sequencing classifier for the preoperative diagnosis of cytologically indeterminate thyroid nodules. JAMA Surg. Published online May 23, 2018. doi:10.1001/jamasurg.2018.1153 eMethods 1. Training cohort eMethods 2. Reference methods eMethods 3. RNA purification eMethods 4. Library preparation eMethods 5. Next-generation sequencing eMethods 6. RNA sequencing pipeline, feature extraction, and quality control eAppendix. Genomic sequence classifier to gene expression classifier comparison on a per- samples basis eTable 1. Blinding of the independent test set eTable 2. Composition of the core ensemble model training set eTable 3. Feature sets used in each classifier within the final ensemble model eTable 4. List of 1115 core genes deriving the ensemble model prediction eTable 5. Histology subtype comparison between validation cohorts eTable 6. Prevalence of malignancy between validation cohorts eTable 7. Performance comparison between the genomic sequence classifier and gene expression classifier eFigure 1. Afirma gene expression classifier system eFigure 2. Standards for Reporting of Diagnostic Accuracy Studies diagram of sample flow through the study eReferences. This supplementary material has been provided by the authors to give readers additional information about their work. © 2018 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 eMethods 1. Training cohort The following thyroid nodule FNA samples were included in the training set (eTable 2): ENHANCE Arm 1 A dedicated molecular sample was obtained when the cytology specimen was collected from a nodule ≥ 1 cm during clinical care, and subsequently underwent surgery. Two central histopathologists reviewed the histopathology slides while blinded to all molecular information and blinded to the local histopathology diagnosis. A similarly blinded third central histopathologist reviewed the slides if agreement was not achieved between the first two blinded histopathologists. If two out of three agreement was not achieved while blinded to each other’s diagnosis, then all 3 histopathologists re-viewed digitized slide images in open discussion to achieve a diagnosis while blinded to all molecular information. ENHANCE Arm 2 A dedicated molecular sample was obtained when the cytology specimen was collected from a nodule ≥ 1 cm during clinical care. Arm 2 samples were all unoperated, Bethesda II, or Bethesda III/IV and GEC benign, and lacked 2015 American Thyroid Association high suspicion sonographic pattern findings. Additionally, they had clinical follow-up (mean 23 months, range 17-32) and either a repeat FNA that was cytology benign, or had no growth (< 50% increase in volume or < 20% increase in 2 or more dimensions) or development of high suspicion ultrasound findings after the initial FNA. Nodules were excluded from Arm 2 if repeat FNA was Bethesda V or VI, GEC suspicious, or they underwent surgery. Arm 2 nodules served as truly benign samples, recognizing that GEC benign samples were underrepresented among operated Arm 1 samples. VERA-CVP (non Cyto-I) samples Samples described in the clinical validation of the Afirma GEC1 with sufficient materials remaining. Only Bethesda II, V, and VI samples with histopathology labels defined by an expert panel of pathologists were allowed in the training set. 60% of these samples were randomly chosen into the training set. VERA-Train Samples used in the training set of the Afirma GEC.1 VERA-Extra Collected and associated with histopathology labels identically to VERA-CVP, but these samples were not used in the training or validation of the Afirma GEC. CLIA-GEC B Samples from the CLIA stream that are GEC Benign. These samples do not have long term follow-up or a histopathology label. Their benign GEC prediction is used as a surrogate label in algorithm training. © 2018 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 eMethods 2. Reference methods BRAF V600E status BRAF V600E status was determined from genomic DNA using Competitive Allele Specific Taqman PCR (castPCR™, Thermo Fisher, Waltham, MA) for BRAF 1799T>A mutation, as previously described.2,3 Briefly, genomic DNA was purified with the AllPrep Micro Kit (Qiagen, Hilden, Germany) and quantified with Quanti-iT PicoGreen dsDNA Assay Kit (Thermo Fisher, Waltham, MA). Five ng of DNA was tested with wild-type and mutant assays on an ABI7900HT. Samples were labelled BRAF V600E positive if the variant allele frequency was ≥5% and wild type if the allele frequency was <5%. Medullary Thyroid Cancer Histopathology diagnoses, including medullary thyroid cancer, were previously established by an expert panel of thyroid histopathologists while blinded to all clinical and molecular data.1 © 2018 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 eMethods 3. RNA purification RNA was purified with the AllPrep Micro kit (Qiagen, Hilden, Germany) as previously described.1 RNA was quantified using the QuantiFluor RNA System (Promega, Madison, WI). Fluorescence was read with a Tecan Infinite 200 Pro plate reader (Tecan, Männedorf, Switzerland). RNA Integrity Number was determined with the Bioanalyzer 2100 (Agilent, Santa Clara, CA). © 2018 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 eMethods 4. Library preparation Samples were randomized and plated into 96 well plates according to their random order. Each plate contained Universal Human Reference RNA (Agilent, Santa Clara, CA), a benign thyroid tissue control sample, a malignant thyroid tissue control sample, a medullary thyroid carcinoma tissue control sample and 6 FNAs that were run on every plate in the study. Additionally, 3 samples from each plate were randomly selected to be included as technical replicates. 15 ng of total RNA was transferred to a 96 well plate. The TruSeq RNA Access Library Preparation Kit (Illumina, San Diego, CA) was adapted for use on the Microlab STAR robotics platform (Hamilton, Reno, NV). During library preparation, total RNA is fragmented, reverse transcribed, end-repaired, A- tailed, and Illumina adapters with individual indexes are ligated. Following PCR and AMpure XP (Beckman Coulter, Indianapolis, IN) cleanup, library size and quantity was determined with the Fragment Analyzer (Advanced Analytical, Ankeny, IA). 250 ng of 4 libraries were combined and sequentially captured with the human exome to remove ribosomal RNA, intronic, and intergenic sequences. Following PCR and AMpure XP (Beckman Coulter, Indianapolis, IN) cleanup, library size and quantity were determined with the Bioanalyzer 2100 (Agilent, Santa Clara, CA). © 2018 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 eMethods 5. Next-generation sequencing Libraries were normalized to 2 nM, pooled to 16 samples per sequencing run, and denatured according to the manufacturer’s instructions. 1% phiX library (Illumina, San Diego, CA) was spiked into each sequencing run. Denatured and diluted libraries were loaded onto NextSeq 500 machines (Illumina, San Diego, CA) and sequenced with a NextSeq v2 High Output 150 cycle kit (Illumina, San Diego, CA) for paired end 2x76 cycle sequencing. Sequencing runs were required to have >75% of bases ≥Q30 and <1% phiX error rate. © 2018 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 eMethods 6. RNA sequencing pipeline, feature extraction, and quality control RNA-seq data was used to generate gene expression counts, identify variants, detect fusion-pairs, and calculate loss of heterozygosity (LOH) statistics. Raw sequencing data (FASTQ file) was aligned to human reference genome assembly 37 (Genome Reference Consortium) using STAR RNA-seq aligner.4 Expression counts were obtained by HTSeq5 and normalized using DESeq26 accounting for sequencing depth and gene-wise variability. Variants were identified using GATK variant calling pipeline,7,8 and fusion-pairs detected using STAR-Fusion.9 We developed a loss of heterozygosity (LOH) statistic at chromosome and genome level using variants identified genome-wide. The statistic quantifies the magnitude of LOH by calculating the proportion of variants that have a variant allele frequency (VAF; fraction of reads carrying the alternative allele) away from 0.5 (<0.2 or >0.8) after pre-filtering of variants that has a VAF exactly at zero or one, or is located in cytoband regions exhibiting abnormal excess of LOH signatures across all training samples. To exclude low quality samples from downstream analysis, quality metrics were evaluated against pre- specified acceptance metrics for total numbers of sequenced and uniquely mapped reads, the overall proportion of exonic reads among mapped, the mean per-base coverage, the uniformity of base coverage, and base duplication and mismatch rates. All these QC metrics were generated using RNA-SeQC.10 Any sample that failed a QC metric was reprocessed from total RNA through library preparation and sequencing if sufficient RNA was available. Only samples passing the quality criteria were used for downstream analysis. © 2018 American Medical Association. All rights reserved. Downloaded From: https://jamanetwork.com/ on 09/28/2021 eAppendix. Genomic sequence classifier to gene expression classifier comparison on
Recommended publications
  • Combined Analysis of Chip-Seq and Gene Microarray Datasets Identify the E2-Mediated Genes in Erα-Dependent Manner in Osteosarcoma
    ONCOLOGY REPORTS 38: 2335-2342, 2017 Combined analysis of ChIP-seq and gene microarray datasets identify the E2-mediated genes in ERα-dependent manner in osteosarcoma KANGSONG TIAN1, WEI QI1, QIAN YAN1, FEnG ZHanG1, DELEI SONG1, HaIyanG ZHanG2 and MING LV1 1Trauma Department of Orthopedics, 2Microscopic Department of Orthopedics, Zibo Central Hospital, Zibo, Shandong 255036, P.R. China Received March 14, 2017; Accepted August 11, 2017 DOI: 10.3892/or.2017.5914 Abstract. Osteosarcoma is a common bone tumor which is Introduction affected by E2, the most representative estrogen. Gene regula- tion function of E2 is highly dependent on estrogen receptor. ESR1, also known as ER-alpha or ERα, is an important The purpose of this study was to explore the gene regulation estrogen receptor (ER) and involves in the gene regula- patterns of E2 through estrogen receptor α (ESR1) in osteo- tion in various diseases and biological processes, including sarcoma based on the combined analysis of ChIP-seq and breast cancer (1,2), osteosarcoma (3) and cell growth (4). gene microarray. All of the datasets were downloaded from 17β-estradiol (E2) is one of the most representative estrogens the Gene Expression Omnibus (GEO). Differential expression responsible for the development and reproductive capability genes (DEGs) in E2 treated U2OS cells expressing ESR1 of female characteristics (5). Besides, it participates in the (U2OS-ERα) compared with those treated with vehicle were progression of many diseases, for example, E2 could regulate obtained based on R programming software. ESR1-specific the proliferation of breast cancer cells through focal adhesion binding sites (peaks) in E2 treated U2OS cells were identified and chemokine signaling pathways (6); through autophagy through MACS.
    [Show full text]
  • Prospective Isolation of NKX2-1–Expressing Human Lung Progenitors Derived from Pluripotent Stem Cells
    The Journal of Clinical Investigation RESEARCH ARTICLE Prospective isolation of NKX2-1–expressing human lung progenitors derived from pluripotent stem cells Finn Hawkins,1,2 Philipp Kramer,3 Anjali Jacob,1,2 Ian Driver,4 Dylan C. Thomas,1 Katherine B. McCauley,1,2 Nicholas Skvir,1 Ana M. Crane,3 Anita A. Kurmann,1,5 Anthony N. Hollenberg,5 Sinead Nguyen,1 Brandon G. Wong,6 Ahmad S. Khalil,6,7 Sarah X.L. Huang,3,8 Susan Guttentag,9 Jason R. Rock,4 John M. Shannon,10 Brian R. Davis,3 and Darrell N. Kotton1,2 2 1Center for Regenerative Medicine, and The Pulmonary Center and Department of Medicine, Boston University School of Medicine, Boston, Massachusetts, USA. 3Center for Stem Cell and Regenerative Medicine, Brown Foundation Institute of Molecular Medicine, University of Texas Health Science Center, Houston, Texas, USA. 4Department of Anatomy, UCSF, San Francisco, California, USA. 5Division of Endocrinology, Diabetes and Metabolism, Beth Israel Deaconess Medical Center and Harvard Medical School, Boston, Massachusetts, USA. 6Department of Biomedical Engineering and Biological Design Center, Boston University, Boston, Massachusetts, USA. 7Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts, USA. 8Columbia Center for Translational Immunology & Columbia Center for Human Development, Columbia University Medical Center, New York, New York, USA. 9Department of Pediatrics, Monroe Carell Jr. Children’s Hospital, Vanderbilt University, Nashville, Tennessee, USA. 10Division of Pulmonary Biology, Cincinnati Children’s Hospital, Cincinnati, Ohio, USA. It has been postulated that during human fetal development, all cells of the lung epithelium derive from embryonic, endodermal, NK2 homeobox 1–expressing (NKX2-1+) precursor cells.
    [Show full text]
  • GJA4/Connexin 37 Mutations Correlate with Secondary Lymphedema Following Surgery in Breast Cancer Patients
    biomedicines Article GJA4/Connexin 37 Mutations Correlate with Secondary Lymphedema Following Surgery in Breast Cancer Patients Mahrooyeh Hadizadeh 1,2, Seiied Mojtaba Mohaddes Ardebili 1, Mansoor Salehi 2, Chris Young 3, Fariborz Mokarian 4, James McClellan 5, Qin Xu 6, Mohammad Kazemi 2, Elham Moazam 4, Behzad Mahaki 7 ID and Maziar Ashrafian Bonab 8,* 1 Department of Medical Genetics, Faculty of Medicine, Tabriz University of Medical Sciences, Tabriz 5166614766, Iran; [email protected] (M.H.); [email protected] (S.M.M.A.) 2 Department of Genetics and Molecular Biology, Isfahan University of Medical Sciences, Isfahan 81746753461, Iran; [email protected] (M.S.); [email protected] (M.K.) 3 School of Allied Health Sciences, Faculty of Health and Life Sciences, De Montfort University, Leicester LE1 9BH, UK; [email protected] 4 Cancer Prevention Research Centre, Isfahan University of Medical Sciences, Isfahan 8184917911, Iran; [email protected] (F.M.); [email protected] (E.M.) 5 School of Biological Sciences, University of Portsmouth, Portsmouth PO1 2DY, UK; [email protected] 6 School of Pharmacy, Faculty of Health and Life Sciences, De Montfort University, Leicester LE1 9BH, UK; [email protected] 7 Department of Occupational Health Engineering, School of Health, Isfahan University of Medical Sciences, Isfahan 8174673461, Iran; [email protected] 8 Department of Biological Sciences, University of Chester, Chester CH1 4BJ, UK * Correspondence: [email protected]; Tel.: +44-(0)1244-513-056 Received: 31 December 2017; Accepted: 13 February 2018; Published: 22 February 2018 Abstract: Lymphedema is a condition resulting from mutations in various genes essential for lymphatic development and function, which leads to obstruction of the lymphatic system.
    [Show full text]
  • Supporting Information
    Supporting Information Mulders et al. 10.1073/pnas.0905780106 SI Materials and Methods in Opti-MEM (Invitrogen) to myoblasts, in both cases at a final Animals. Hemizygous DM500 mice, derived from the DM300– oligo concentration of 200 nM. Fresh medium was supplemented 328 line (1), express a transgenic human DM1 locus, which bears after 4 hours. After 24 h, medium was changed. RNA was a repeat segment that has expanded to Ϸ500 CTG triplets, isolated 48 h after transfection. because of intergenerational triplet-repeat instability. For the isolation of immortal DM500 myoblasts, DM500 mice were RNA Isolation. Typically, RNA from cultured cells was isolated crossed with H-2Kb-tsA58 transgenic mice (2). Homozygous using the Aurum Total RNA mini kit (guanidine-HCl/ HSALR20b mice express a (CUG)250 segment in the context of mercaptoethanol-based lysis, silica membrane binding; Bio- a human skeletal actin transcript (3). All animal experiments Rad), according to the manufacturer’s protocol. RNA from were approved by the Institutional Animal Care and Use Com- muscle tissue was isolated using TRIzol reagent (Invitrogen). mittees of the Radboud University Nijmegen and the University Alternative methods to isolate RNA from cultured cells in- of Rochester Medical Center. volved: (i) use of the TRIzol reagent, according to the manu- facturer’s protocol; (ii) an oligo(dT) affinity column (Nucleo- Cell Culture. An immortal mouse myoblast cell culture expressing Trap mRNA mini kit; Macherey-Nagel) for the isolation of hDMPK (CUG)500 transcripts was derived from GPS tissue poly(A) RNA; or (iii) a TRIzol procedure preceded by a isolated from DM500 mice additionally expressing 1 copy of the proteinase K digestion of the whole cell lysate (7): in short, cells H-2Kb-tsA58 allele (4).
    [Show full text]
  • MOZ and MORF, Two Large Mystic Hats in Normal and Cancer Stem Cells
    Oncogene (2007) 26, 5408–5419 & 2007 Nature Publishing Group All rights reserved 0950-9232/07 $30.00 www.nature.com/onc REVIEW MOZ and MORF, two large MYSTic HATs in normal and cancer stem cells X-J Yang and M Ullah Molecular Oncology Group, Department of Medicine, McGill University Health Center, Montre´al, Que´bec, Canada Genes of the human monocytic leukemia zinc-finger protein pattern. For cancer biology, it is thus important to MOZ (HUGO symbol, MYST3) and its paralog MORF understand the fundamental mechanisms whereby chro- (MYST4) are rearranged in chromosome translocations matin structure and function are regulated. In the associated withacute myeloid leukemia and/or benign past two decades, our knowledge about regulation in uterine leiomyomata. Both proteins have intrinsic histone this field has exploded. Known regulatory mechanisms acetyltransferase activity and are components of quartet include chromatin assembly, ATP-dependent remodeling, complexes withnoncatalytic subunits containing thebromo- covalent modification, condensin-mediated condensation, domain, plant homeodomain-linked (PHD) finger and replacement with histone variants, and association of proline-tryptophan-tryptophan-proline (PWWP)-containing noncoding RNA (reviewed by Horn and Peterson, 2002; domain, three types of structural modules characteristic of Khorasanizadeh, 2004; Li et al., 2007). Covalent chromatin regulators. Although leukemia-derived fusion pro- modification can occur at both the DNA and histone teins suchas MOZ-TIF2 promote self-renewal of leukemic levels. With histones, modifications include acetylation, stem cells, recent studies indicate that murine MOZ and phosphorylation, methylation, ubiquitination, and many MORF are important for proper development of hema- others (for reviews, see Spencer and Davie, 1999; Strahl topoietic and neurogenic progenitors, respectively, thereby and Allis, 2000; Berger, 2002; Jason et al., 2002; highlighting the importance of epigenetic integrity in Kouzarides, 2007).
    [Show full text]
  • Seq2pathway Vignette
    seq2pathway Vignette Bin Wang, Xinan Holly Yang, Arjun Kinstlick May 19, 2021 Contents 1 Abstract 1 2 Package Installation 2 3 runseq2pathway 2 4 Two main functions 3 4.1 seq2gene . .3 4.1.1 seq2gene flowchart . .3 4.1.2 runseq2gene inputs/parameters . .5 4.1.3 runseq2gene outputs . .8 4.2 gene2pathway . 10 4.2.1 gene2pathway flowchart . 11 4.2.2 gene2pathway test inputs/parameters . 11 4.2.3 gene2pathway test outputs . 12 5 Examples 13 5.1 ChIP-seq data analysis . 13 5.1.1 Map ChIP-seq enriched peaks to genes using runseq2gene .................... 13 5.1.2 Discover enriched GO terms using gene2pathway_test with gene scores . 15 5.1.3 Discover enriched GO terms using Fisher's Exact test without gene scores . 17 5.1.4 Add description for genes . 20 5.2 RNA-seq data analysis . 20 6 R environment session 23 1 Abstract Seq2pathway is a novel computational tool to analyze functional gene-sets (including signaling pathways) using variable next-generation sequencing data[1]. Integral to this tool are the \seq2gene" and \gene2pathway" components in series that infer a quantitative pathway-level profile for each sample. The seq2gene function assigns phenotype-associated significance of genomic regions to gene-level scores, where the significance could be p-values of SNPs or point mutations, protein-binding affinity, or transcriptional expression level. The seq2gene function has the feasibility to assign non-exon regions to a range of neighboring genes besides the nearest one, thus facilitating the study of functional non-coding elements[2]. Then the gene2pathway summarizes gene-level measurements to pathway-level scores, comparing the quantity of significance for gene members within a pathway with those outside a pathway.
    [Show full text]
  • Defining Functional Interactions During Biogenesis of Epithelial Junctions
    ARTICLE Received 11 Dec 2015 | Accepted 13 Oct 2016 | Published 6 Dec 2016 | Updated 5 Jan 2017 DOI: 10.1038/ncomms13542 OPEN Defining functional interactions during biogenesis of epithelial junctions J.C. Erasmus1,*, S. Bruche1,*,w, L. Pizarro1,2,*, N. Maimari1,3,*, T. Poggioli1,w, C. Tomlinson4,J.Lees5, I. Zalivina1,w, A. Wheeler1,w, A. Alberts6, A. Russo2 & V.M.M. Braga1 In spite of extensive recent progress, a comprehensive understanding of how actin cytoskeleton remodelling supports stable junctions remains to be established. Here we design a platform that integrates actin functions with optimized phenotypic clustering and identify new cytoskeletal proteins, their functional hierarchy and pathways that modulate E-cadherin adhesion. Depletion of EEF1A, an actin bundling protein, increases E-cadherin levels at junctions without a corresponding reinforcement of cell–cell contacts. This unexpected result reflects a more dynamic and mobile junctional actin in EEF1A-depleted cells. A partner for EEF1A in cadherin contact maintenance is the formin DIAPH2, which interacts with EEF1A. In contrast, depletion of either the endocytic regulator TRIP10 or the Rho GTPase activator VAV2 reduces E-cadherin levels at junctions. TRIP10 binds to and requires VAV2 function for its junctional localization. Overall, we present new conceptual insights on junction stabilization, which integrate known and novel pathways with impact for epithelial morphogenesis, homeostasis and diseases. 1 National Heart and Lung Institute, Faculty of Medicine, Imperial College London, London SW7 2AZ, UK. 2 Computing Department, Imperial College London, London SW7 2AZ, UK. 3 Bioengineering Department, Faculty of Engineering, Imperial College London, London SW7 2AZ, UK. 4 Department of Surgery & Cancer, Faculty of Medicine, Imperial College London, London SW7 2AZ, UK.
    [Show full text]
  • Anti-GJA4 / Connexin 37 Antibody (ARG58815)
    Product datasheet [email protected] ARG58815 Package: 50 μg anti-GJA4 / Connexin 37 antibody Store at: -20°C Summary Product Description Rabbit Polyclonal antibody recognizes GJA4 / Connexin 37 Tested Reactivity Hu, Ms, Rat Predict Reactivity Hm Tested Application ICC, IHC-Fr, WB Host Rabbit Clonality Polyclonal Isotype IgG Target Name GJA4 / Connexin 37 Species Human Immunogen Synthetic peptide corresponding to aa. 3-17 of Human Connexin 37 (DWGFLEKLLDQVQEH). Conjugation Un-conjugated Alternate Names Connexin-37; Gap junction alpha-4 protein; CX37; Cx37 Application Instructions Application table Application Dilution ICC 0.5 - 1 µg/ml IHC-Fr 1:200 - 1:1000 WB 0.1 - 0.5 µg/ml Application Note * The dilutions indicate recommended starting dilutions and the optimal dilutions or concentrations should be determined by the scientist. Properties Form Liquid Purification Affinity purification with immunogen. Buffer 0.9% NaCl, 0.2% Na2HPO4, 0.05% Thimerosal, 0.05% Sodium azide and 5% BSA. Preservative 0.05% Thimerosal and 0.05% Sodium azide Stabilizer 5% BSA Concentration 0.5 mg/ml Storage instruction For continuous use, store undiluted antibody at 2-8°C for up to a week. For long-term storage, aliquot and store at -20°C or below. Storage in frost free freezers is not recommended. Avoid repeated freeze/thaw cycles. Suggest spin the vial prior to opening. The antibody solution should be gently mixed before use. www.arigobio.com 1/3 Note For laboratory research only, not for drug, diagnostic or other use. Bioinformation Gene Symbol GJA4 Gene Full Name gap junction protein, alpha 4, 37kDa Background This gene encodes a member of the connexin gene family.
    [Show full text]
  • A Computational Approach for Defining a Signature of Β-Cell Golgi Stress in Diabetes Mellitus
    Page 1 of 781 Diabetes A Computational Approach for Defining a Signature of β-Cell Golgi Stress in Diabetes Mellitus Robert N. Bone1,6,7, Olufunmilola Oyebamiji2, Sayali Talware2, Sharmila Selvaraj2, Preethi Krishnan3,6, Farooq Syed1,6,7, Huanmei Wu2, Carmella Evans-Molina 1,3,4,5,6,7,8* Departments of 1Pediatrics, 3Medicine, 4Anatomy, Cell Biology & Physiology, 5Biochemistry & Molecular Biology, the 6Center for Diabetes & Metabolic Diseases, and the 7Herman B. Wells Center for Pediatric Research, Indiana University School of Medicine, Indianapolis, IN 46202; 2Department of BioHealth Informatics, Indiana University-Purdue University Indianapolis, Indianapolis, IN, 46202; 8Roudebush VA Medical Center, Indianapolis, IN 46202. *Corresponding Author(s): Carmella Evans-Molina, MD, PhD ([email protected]) Indiana University School of Medicine, 635 Barnhill Drive, MS 2031A, Indianapolis, IN 46202, Telephone: (317) 274-4145, Fax (317) 274-4107 Running Title: Golgi Stress Response in Diabetes Word Count: 4358 Number of Figures: 6 Keywords: Golgi apparatus stress, Islets, β cell, Type 1 diabetes, Type 2 diabetes 1 Diabetes Publish Ahead of Print, published online August 20, 2020 Diabetes Page 2 of 781 ABSTRACT The Golgi apparatus (GA) is an important site of insulin processing and granule maturation, but whether GA organelle dysfunction and GA stress are present in the diabetic β-cell has not been tested. We utilized an informatics-based approach to develop a transcriptional signature of β-cell GA stress using existing RNA sequencing and microarray datasets generated using human islets from donors with diabetes and islets where type 1(T1D) and type 2 diabetes (T2D) had been modeled ex vivo. To narrow our results to GA-specific genes, we applied a filter set of 1,030 genes accepted as GA associated.
    [Show full text]
  • Genome-Wide Profiling of Active Enhancers in Colorectal Cancer
    Genome-wide proling of active enhancers in colorectal cancer Min Wu ( [email protected] ) Wuhan University https://orcid.org/0000-0003-1372-4764 Qinglan Li Wuhan University Xiang Lin Wuhan University Ya-Li Yu Zhongnan Hospital, Wuhan University Lin Chen Wuhan University Qi-Xin Hu Wuhan University Meng Chen Zhongnan Hospital, Wuhan University Nan Cao Zhongnan Hospital, Wuhan University Chen Zhao Wuhan University Chen-Yu Wang Wuhan University Cheng-Wei Huang Wuhan University Lian-Yun Li Wuhan University Mei Ye Zhongnan Hospital, Wuhan University https://orcid.org/0000-0002-9393-3680 Article Keywords: Colorectal cancer, H3K27ac, Epigenetics, Enhancer, Transcription factors Posted Date: December 10th, 2020 DOI: https://doi.org/10.21203/rs.3.rs-119156/v1 License: This work is licensed under a Creative Commons Attribution 4.0 International License. Read Full License Genome-wide profiling of active enhancers in colorectal cancer Qing-Lan Li1, #, Xiang Lin1, #, Ya-Li Yu2, #, Lin Chen1, #, Qi-Xin Hu1, Meng Chen2, Nan Cao2, Chen Zhao1, Chen-Yu Wang1, Cheng-Wei Huang1, Lian-Yun Li1, Mei Ye2,*, Min Wu1,* 1 Frontier Science Center for Immunology and Metabolism, Hubei Key Laboratory of Cell Homeostasis, Hubei Key Laboratory of Developmentally Originated Disease, Hubei Key Laboratory of Intestinal and Colorectal Diseases, College of Life Sciences, Wuhan University, Wuhan, Hubei 430072, China 2Division of Gastroenterology, Department of Geriatrics, Hubei Clinical Centre & Key Laboratory of Intestinal and Colorectal Diseases, Zhongnan Hospital, Wuhan University, Wuhan, Hubei 430072, China #Equal contribution to the study. Contact information *Correspondence should be addressed to Dr. Min Wu, Email: [email protected], Tel: 86-27-68756620, or Dr.
    [Show full text]
  • The Biomarkers of Key Mirnas and Target Genes Associated with Acute Myocardial Infarction
    The biomarkers of key miRNAs and target genes associated with acute myocardial infarction Qi Wang1, Bingyan Liu2,3, Yuanyong Wang4, Baochen Bai1, Tao Yu3 and Xian–ming Chu1,5 1 Department of Cardiology, The Affiliated hospital of Qingdao University, Qingdao, China 2 School of Basic Medicine, Qingdao University, Qingdao, China 3 Institute for Translational Medicine, Qingdao University, Qingdao, China 4 Department of Thoracic Surgery, Affiliated Hospital of Qingdao University, Qingdao, China 5 Department of Cardiology, The Affiliated Cardiovascular Hospital of Qingdao University, Qingdao, China ABSTRACT Background. Acute myocardial infarction (AMI) is considered one of the most prominent causes of death from cardiovascular disease worldwide. Knowledge of the molecular mechanisms underlying AMI remains limited. Accurate biomarkers are needed to predict the risk of AMI and would be beneficial for managing the incidence rate. The gold standard for the diagnosis of AMI, the cardiac troponin T (cTnT) assay, requires serial testing, and the timing of measurement with respect to symptoms affects the results. As attractive candidate diagnostic biomarkers in AMI, circulating microRNAs (miRNAs) are easily detectable, generally stable and tissue specific. Methods. The Gene Expression Omnibus (GEO) database was used to compare miRNA expression between AMI and control samples, and the interactions between miRNAs and mRNAs were analysed for expression and function. Furthermore, a protein-protein interaction (PPI) network was constructed. The miRNAs identified in the bioinformatic analysis were verified by RT-qPCR in an H9C2 cell line. The miRNAs in plasma samples from patients with AMI (n D 11) and healthy controls (n D 11) were used to construct Submitted 23 December 2019 receiver operating characteristic (ROC) curves to evaluate the clinical prognostic value Accepted 14 April 2020 of the identified miRNAs.
    [Show full text]
  • Downloaded from [266]
    Patterns of DNA methylation on the human X chromosome and use in analyzing X-chromosome inactivation by Allison Marie Cotton B.Sc., The University of Guelph, 2005 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in The Faculty of Graduate Studies (Medical Genetics) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) January 2012 © Allison Marie Cotton, 2012 Abstract The process of X-chromosome inactivation achieves dosage compensation between mammalian males and females. In females one X chromosome is transcriptionally silenced through a variety of epigenetic modifications including DNA methylation. Most X-linked genes are subject to X-chromosome inactivation and only expressed from the active X chromosome. On the inactive X chromosome, the CpG island promoters of genes subject to X-chromosome inactivation are methylated in their promoter regions, while genes which escape from X- chromosome inactivation have unmethylated CpG island promoters on both the active and inactive X chromosomes. The first objective of this thesis was to determine if the DNA methylation of CpG island promoters could be used to accurately predict X chromosome inactivation status. The second objective was to use DNA methylation to predict X-chromosome inactivation status in a variety of tissues. A comparison of blood, muscle, kidney and neural tissues revealed tissue-specific X-chromosome inactivation, in which 12% of genes escaped from X-chromosome inactivation in some, but not all, tissues. X-linked DNA methylation analysis of placental tissues predicted four times higher escape from X-chromosome inactivation than in any other tissue. Despite the hypomethylation of repetitive elements on both the X chromosome and the autosomes, no changes were detected in the frequency or intensity of placental Cot-1 holes.
    [Show full text]