Figure S1 the Pipeline of Calling and Filtering Rare Variants from RNA-Seq Data in This Study

Total Page:16

File Type:pdf, Size:1020Kb

Figure S1 the Pipeline of Calling and Filtering Rare Variants from RNA-Seq Data in This Study BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S1 The pipeline of calling and filtering rare variants from RNA-seq data in this study. 1 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S2 Validation of variants called from RNA-seq by Sanger sequencing in an independent cohort of three samples. Sanger sequencing chromatograms of the validated variants among 10 variants randomly selected from each samples. Sample1 Sample2 2 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Sample3 3 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S3 The sample quality control filtering based on the number of rare variants carried by each sample. The number of rare variants called from RNA-seq of each sample was plotted and those outside of the linear phase of the distribution were filtered out. The number of variants in the linear phase is shown in orange and the number of variants which is either too low or too high in the exponential phase is shown in grey, indicating the excluded samples. 250 200 150 100 variant number variant 50 0 0 100 200 300 400 500 600 sample ID 4 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S4 The distribution of mutation types among the 63 recurrent rare variants found in our study. SNVs 2% 5% missense_NA 30% missense_D missense_T 63% stoploss/stopgain_NA missense_NA: rare missense SNVs with SIFT functional prediction not available; missense_D: rare missense SNVs with SIFT functional prediction to be deleterious; missense_T: rare missense SNVs with SIFT functional prediction to be tolerant; stoploss/stopgain_NA: SNVs with a gene truncating effect (stop gain or stop loss) with SIFT functional prediction not available. 5 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S5. The expression level of genes with recurrent deleterious rare variants based on data in DICE. (https://dice-database.org/) A. polyarticular JIA 6 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis B. oligoarticular JIA C. systemic JIA DICE:Database of Immune Cell Expression; TPM:Transcripts per million. 7 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S6 The expression level of genes with recurrent deleterious rare variants based on in HPA (https://www.proteinatlas.org/) A. polyarticular JIA 8 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis B. oligoarticular JIA C. systemic JIA HPA: Human protein Atlas; pTPM: Protein- transcripts per million. 9 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S7 The differentially expressed genes in the meta-analysis of polyarticular and oligoarticular subtypes of three RNA-seq datasets (|FC|>2, adjusted P-value <0.05). Volcano plot showing the gene differential expression analysis of three RNA-seq datasets. The vertical lines correspond to 2.0-fold up and down, respectively, and the horizontal line represents adjust-p value of 0.05. The red points in the plot represent the 43 significantly up-regulated genes; and the blue points represent the 4 significantly down-regulated genes. 10 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Figure S8 Volcano plots showing the results of pathway over-representation analysis. It shows the log of the FDR versus the enrichment ratio for all the pathways in the KEGG database, highlighting the degree by which the significant pathways stand out from the background. The size and color of the dot is proportional to the number of input genes falling into each pathway. (A) Pathway analysis of 63 genes with recurrent rare variants; (B) Pathway analysis of 138 differentially expressed genes of polyarticular and oligoarticular subtypes of JIA. A B 11 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis Table S1 The fraction of rare variants called from RNA-sequencing data were detected in whole exome sequencing (WES) data of each of the 100 randomly selected GTEx samples. RNA-seq WES Fraction RNA-seq WES Fraction SRR1069024 SRR2167400 72.70% SRR1361368 SRR2167124 60.80% SRR1071931 SRR2165980 66.20% SRR1363611 SRR2165743 66.70% SRR1072345 SRR2168332 77.40% SRR1366543 SRR2165999 80.40% SRR1072529 SRR2166291 88.00% SRR1368529 SRR2167121 75.80% SRR1072871 SRR2167072 83.10% SRR1370742 SRR2165890 83.30% SRR1074646 SRR2167359 62.70% SRR1378026 SRR2167957 82.10% SRR1076511 SRR2165120 69.20% SRR1381435 SRR2168744 82.20% SRR1077503 SRR2166156 80.90% SRR1381481 SRR2166915 87.80% SRR1079259 SRR2165448 60.40% SRR1382615 SRR2168325 72.20% SRR1079660 SRR2167414 79.70% SRR1383855 SRR2156637 78.30% SRR1080998 SRR2167470 80.40% SRR1386967 SRR2166500 76.50% SRR1083432 SRR1343762 76.70% SRR1388081 SRR2170921 64.90% SRR1083657 SRR2167152 62.20% SRR1389133 SRR1310297 74.00% SRR1084577 SRR2165457 90.20% SRR1392547 SRR1348058 76.60% SRR1084720 SRR2164914 77.10% SRR1396191 SRR2167760 85.50% SRR1087654 SRR2167068 68.80% SRR1397673 SRR2165117 85.70% SRR1090508 SRR2166372 91.30% SRR1400032 SRR2163570 82.20% SRR1091452 SRR2164811 62.70% SRR1400054 SRR2163573 82.50% SRR1092422 SRR2166301 100% SRR1407758 SRR2157513 78.80% SRR1096571 SRR1305587 75.80% SRR1409510 SRR2167950 84.20% SRR1099260 SRR2167942 90.70% SRR1413892 SRR2166259 86.50% SRR1099648 SRR2167756 91.40% SRR1416560 SRR2167440 85.50% SRR1099744 SRR2165963 72.70% SRR1420669 SRR2170534 69.70% SRR1101324 SRR2166998 77.90% SRR1433044 SRR2167335 78.70% SRR1101787 SRR2167293 70.50% SRR1433860 SRR2166939 86.80% SRR1102031 SRR2166948 70.20% SRR1435833 SRR2164488 68.90% SRR1310412 SRR2165954 87.10% SRR1440224 SRR2165753 77.80% SRR1311419 SRR2166173 88.50% SRR1440603 SRR2167421 91.70% SRR1311644 SRR2167963 85.30% SRR1444852 SRR2165960 78.60% SRR1312598 SRR2170783 60.70% SRR1446244 SRR2170487 65.90% SRR1314916 SRR2167408 66.70% SRR1446850 SRR2167131 75.00% SRR1320137 SRR2166827 71.20% SRR1447523 SRR2167496 90.40% SRR1322619 SRR2165570 64.10% SRR1457574 SRR2167277 84.60% SRR1326536 SRR2167467 85.50% SRR1459739 SRR2166269 70.30% SRR1328407 SRR2164694 94.40% SRR1464618 SRR2165024 69.10% SRR1329047 SRR2167309 89.60% SRR1468850 SRR2163078 76.90% SRR1331311 SRR2165452 67.90% SRR1474839 SRR2165523 84.60% SRR1331684 SRR2170126 86.00% SRR1475930 SRR2166052 72.60% SRR1332007 SRR2164492 84.80% SRR1476207 SRR1306141 66.00% SRR1333570 SRR2165222 93.50% SRR1476300 SRR2166515 81.50% SRR1337540 SRR2166369 78.30% SRR1477991 SRR2167387 69.00% 12 Meng X, et al. Ann Rheum Dis 2021; 80:626–631. doi: 10.1136/annrheumdis-2020-218359 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) Ann Rheum Dis SRR1338880 SRR2167707 62.50% SRR1480789 SRR2167332 90.50% SRR1339619 SRR1310396 60.00% SRR1481060 SRR2164704 88.00% SRR1345266 SRR2156686 66.00% SRR1481351 SRR1308426 78.80% SRR1345541 SRR2155763 81.20% SRR1486021 SRR2164933 62.70% SRR1345847 SRR2170211 63.20% SRR1488261 SRR1307843 84.10% SRR1349032 SRR2167437 62.70% SRR1489116 SRR2157311 78.90% SRR1349276 SRR2166865 60.90% SRR1489926 SRR2164748 83.70% SRR1352530 SRR1307722 79.40% SRR1491332 SRR2165363 88.60% SRR1355102 SRR2156394 73.00% SRR1498313 SRR2170789 88.90% RNA-seq= the SRR ID of the RNA-seq file; WES= the SRR ID of the WES file for the same sample as the RNA-seq file; fraction= the fraction of rare variants called from RNA-seq data was detected by the WES data 13 Meng X, et al.
Recommended publications
  • Investigation of the Fatty Acid Transporter-Encoding Genes
    www.nature.com/scientificreports OPEN Investigation of the fatty acid transporter-encoding genes SLC27A3 and SLC27A4 in autism Received: 13 February 2015 1 1 1 1 Accepted: 12 October 2015 Motoko Maekawa , Yoshimi Iwayama , Tetsuo Ohnishi , Manabu Toyoshima , 1 1 1 1 2,3 Published: 09 November 2015 Chie Shimamoto , Yasuko Hisano , Tomoko Toyota , Shabeesh Balan , Hideo Matsuzaki , Yasuhide Iwata3, Shu Takagai3, Kohei Yamada3, Motonori Ota4, Satoshi Fukuchi5, Yohei Okada6,7, Wado Akamatsu6,8, Masatsugu Tsujii3,9, Nobuhiko Kojima10, Yuji Owada11, Hideyuki Okano6, Norio Mori3 & Takeo Yoshikawa1 The solute carrier 27A (SLC27A) gene family encodes fatty acid transport proteins (FATPs) and includes 6 members. During fetal and postnatal periods of development, the growing brain requires a reliable supply of fatty acids. Because autism spectrum disorders (ASD) are now recognized as disorders caused by impaired early brain development, it is possible that functional abnormalities of SLC27A genes may contribute to the pathogenesis of ASD. Here, we confirmed the expression of SLC27A3 and SLC27A4 in human neural stem cells derived from human induced pluripotent stem cells, which suggested their involvement in the developmental stage of the central nervous system. Additionally, we resequenced the SLC27A3 and SLC27A4 genes using 267 ASD patient and 1140 control samples and detected 47 (44 novel and 29 nonsynonymous) and 30 (17 novel and 14 nonsynonymous) variants for the SLC27A3 and SLC27A4, respectively, revealing that they are highly polymorphic with multiple rare variants. The SLC27A4 Ser209 allele was more frequently represented in ASD samples. Furthermore, we showed that a SLC27A4 Ser209 mutant resulted in significantly higher fluorescently-labeled fatty acid uptake into bEnd3 cells, a mouse brain capillary-derived endothelial cell line, compared with SLC27A4 Gly209, suggesting that the functional change may contribute to ASD pathophysiology.
    [Show full text]
  • Identification of the Binding Partners for Hspb2 and Cryab Reveals
    Brigham Young University BYU ScholarsArchive Theses and Dissertations 2013-12-12 Identification of the Binding arP tners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non- Redundant Roles for Small Heat Shock Proteins Kelsey Murphey Langston Brigham Young University - Provo Follow this and additional works at: https://scholarsarchive.byu.edu/etd Part of the Microbiology Commons BYU ScholarsArchive Citation Langston, Kelsey Murphey, "Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non-Redundant Roles for Small Heat Shock Proteins" (2013). Theses and Dissertations. 3822. https://scholarsarchive.byu.edu/etd/3822 This Thesis is brought to you for free and open access by BYU ScholarsArchive. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of BYU ScholarsArchive. For more information, please contact [email protected], [email protected]. Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactions and Non-Redundant Roles for Small Heat Shock Proteins Kelsey Langston A thesis submitted to the faculty of Brigham Young University in partial fulfillment of the requirements for the degree of Master of Science Julianne H. Grose, Chair William R. McCleary Brian Poole Department of Microbiology and Molecular Biology Brigham Young University December 2013 Copyright © 2013 Kelsey Langston All Rights Reserved ABSTRACT Identification of the Binding Partners for HspB2 and CryAB Reveals Myofibril and Mitochondrial Protein Interactors and Non-Redundant Roles for Small Heat Shock Proteins Kelsey Langston Department of Microbiology and Molecular Biology, BYU Master of Science Small Heat Shock Proteins (sHSP) are molecular chaperones that play protective roles in cell survival and have been shown to possess chaperone activity.
    [Show full text]
  • UNIVERSITY of CALIFORNIA Los Angeles
    UNIVERSITY OF CALIFORNIA Los Angeles Dissecting transcriptional control by Klf4 in somatic cell reprogramming A dissertation submitted in partial satisfaction of the requirements for the degree Doctor of Philosophy in Biological Chemistry by Huajun Zhou 2017 ABSTRACT OF THE DISSERTATION Dissecting transcriptional control by Klf4 in somatic cell reprogramming by Huajun Zhou Doctor of Philosophy in Biological Chemistry University of California, Los Angeles, 2017 Professor Gregory S. Payne, Chair Ectopic expression of four transcription factors, Oct4, Sox2, Klf4, and c-Myc, coverts somatic cells directly into induced pluripotent stem cells (iPSCs), which are functionally equivalent to embryonic stem cells (ESCs). The discovery of iPSC has been reshaping the methodology of disease modeling and drug screening in the past decade, and provides tremendous promise for regenerative medicine. However, the mechanism underlying this conversion process, reprogramming, is not yet fully understood. I aimed to dissect the reprogramming process, by characterizing the functional domains of one reprogramming factor Klf4. The transcriptional activation domain (TAD) of Klf4 was revealed to be critical for reprogramming. To search for the factors that mediates the functionality of Klf4 TAD, I identified transcriptional coactivators CBP/p300 and Mediator complex as the physical interaction partners of Klf4 TAD, and further showed that this interaction is functionally required ii for Klf4 mediated transcriptional activation in reprograming. Clathrin heavy chain, initially identified as a physical interaction partner of Klf4 TAD, was shown to be not required for Klf4 transcriptional activation. Clathrin heavy chain was furthered characterized for its potential transcriptional activation activity in CHC-TFE3, a chromosomal fusion discovered in renal cell carcinoma.
    [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]
  • Views of the NIH
    CLINICAL EPIDEMIOLOGY www.jasn.org Genetic Variants Associated with Circulating Fibroblast Growth Factor 23 Cassianne Robinson-Cohen ,1 Traci M. Bartz,2 Dongbing Lai,3 T. Alp Ikizler,1 Munro Peacock,4 Erik A. Imel,4 Erin D. Michos,5 Tatiana M. Foroud,3 Kristina Akesson,6,7 Kent D. Taylor,8 Linnea Malmgren,6,7 Kunihiro Matsushita,5,9,10 Maria Nethander,11 Joel Eriksson,12 Claes Ohlsson,12 Daniel Mellström,12 Myles Wolf,13 Osten Ljunggren,14 Fiona McGuigan,6,7 Jerome I. Rotter,8 Magnus Karlsson,6,7 Michael J. Econs,3,4 Joachim H. Ix,15,16 Pamela L. Lutsey,17 Bruce M. Psaty,18,19 Ian H. de Boer ,20 and Bryan R. Kestenbaum 20 Due to the number of contributing authors, the affiliations are listed at the end of this article. ABSTRACT Background Fibroblast growth factor 23 (FGF23), a bone-derived hormone that regulates phosphorus and vitamin D metabolism, contributes to the pathogenesis of mineral and bone disorders in CKD and is an emerging cardiovascular risk factor. Central elements of FGF23 regulation remain incompletely under- stood; genetic variation may help explain interindividual differences. Methods We performed a meta-analysis of genome-wide association studies of circulating FGF23 con- centrations among 16,624 participants of European ancestry from seven cohort studies, excluding par- ticipants with eGFR,30 ml/min per 1.73 m2 to focus on FGF23 under normal conditions. We evaluated the association of single-nucleotide polymorphisms (SNPs) with natural log–transformed FGF23 concentra- tion, adjusted for age, sex, study site, and principal components of ancestry.
    [Show full text]
  • Supplemental Table 1. Complete Gene Lists and GO Terms from Figure 3C
    Supplemental Table 1. Complete gene lists and GO terms from Figure 3C. Path 1 Genes: RP11-34P13.15, RP4-758J18.10, VWA1, CHD5, AZIN2, FOXO6, RP11-403I13.8, ARHGAP30, RGS4, LRRN2, RASSF5, SERTAD4, GJC2, RHOU, REEP1, FOXI3, SH3RF3, COL4A4, ZDHHC23, FGFR3, PPP2R2C, CTD-2031P19.4, RNF182, GRM4, PRR15, DGKI, CHMP4C, CALB1, SPAG1, KLF4, ENG, RET, GDF10, ADAMTS14, SPOCK2, MBL1P, ADAM8, LRP4-AS1, CARNS1, DGAT2, CRYAB, AP000783.1, OPCML, PLEKHG6, GDF3, EMP1, RASSF9, FAM101A, STON2, GREM1, ACTC1, CORO2B, FURIN, WFIKKN1, BAIAP3, TMC5, HS3ST4, ZFHX3, NLRP1, RASD1, CACNG4, EMILIN2, L3MBTL4, KLHL14, HMSD, RP11-849I19.1, SALL3, GADD45B, KANK3, CTC- 526N19.1, ZNF888, MMP9, BMP7, PIK3IP1, MCHR1, SYTL5, CAMK2N1, PINK1, ID3, PTPRU, MANEAL, MCOLN3, LRRC8C, NTNG1, KCNC4, RP11, 430C7.5, C1orf95, ID2-AS1, ID2, GDF7, KCNG3, RGPD8, PSD4, CCDC74B, BMPR2, KAT2B, LINC00693, ZNF654, FILIP1L, SH3TC1, CPEB2, NPFFR2, TRPC3, RP11-752L20.3, FAM198B, TLL1, CDH9, PDZD2, CHSY3, GALNT10, FOXQ1, ATXN1, ID4, COL11A2, CNR1, GTF2IP4, FZD1, PAX5, RP11-35N6.1, UNC5B, NKX1-2, FAM196A, EBF3, PRRG4, LRP4, SYT7, PLBD1, GRASP, ALX1, HIP1R, LPAR6, SLITRK6, C16orf89, RP11-491F9.1, MMP2, B3GNT9, NXPH3, TNRC6C-AS1, LDLRAD4, NOL4, SMAD7, HCN2, PDE4A, KANK2, SAMD1, EXOC3L2, IL11, EMILIN3, KCNB1, DOK5, EEF1A2, A4GALT, ADGRG2, ELF4, ABCD1 Term Count % PValue Genes regulation of pathway-restricted GDF3, SMAD7, GDF7, BMPR2, GDF10, GREM1, BMP7, LDLRAD4, SMAD protein phosphorylation 9 6.34 1.31E-08 ENG pathway-restricted SMAD protein GDF3, SMAD7, GDF7, BMPR2, GDF10, GREM1, BMP7, LDLRAD4, phosphorylation
    [Show full text]
  • Supplemental Information
    Supplemental information Dissection of the genomic structure of the miR-183/96/182 gene. Previously, we showed that the miR-183/96/182 cluster is an intergenic miRNA cluster, located in a ~60-kb interval between the genes encoding nuclear respiratory factor-1 (Nrf1) and ubiquitin-conjugating enzyme E2H (Ube2h) on mouse chr6qA3.3 (1). To start to uncover the genomic structure of the miR- 183/96/182 gene, we first studied genomic features around miR-183/96/182 in the UCSC genome browser (http://genome.UCSC.edu/), and identified two CpG islands 3.4-6.5 kb 5’ of pre-miR-183, the most 5’ miRNA of the cluster (Fig. 1A; Fig. S1 and Seq. S1). A cDNA clone, AK044220, located at 3.2-4.6 kb 5’ to pre-miR-183, encompasses the second CpG island (Fig. 1A; Fig. S1). We hypothesized that this cDNA clone was derived from 5’ exon(s) of the primary transcript of the miR-183/96/182 gene, as CpG islands are often associated with promoters (2). Supporting this hypothesis, multiple expressed sequences detected by gene-trap clones, including clone D016D06 (3, 4), were co-localized with the cDNA clone AK044220 (Fig. 1A; Fig. S1). Clone D016D06, deposited by the German GeneTrap Consortium (GGTC) (http://tikus.gsf.de) (3, 4), was derived from insertion of a retroviral construct, rFlpROSAβgeo in 129S2 ES cells (Fig. 1A and C). The rFlpROSAβgeo construct carries a promoterless reporter gene, the β−geo cassette - an in-frame fusion of the β-galactosidase and neomycin resistance (Neor) gene (5), with a splicing acceptor (SA) immediately upstream, and a polyA signal downstream of the β−geo cassette (Fig.
    [Show full text]
  • Comprehensive Array CGH of Normal Karyotype Myelodysplastic
    Leukemia (2011) 25, 387–399 & 2011 Macmillan Publishers Limited All rights reserved 0887-6924/11 www.nature.com/leu LEADING ARTICLE Comprehensive array CGH of normal karyotype myelodysplastic syndromes reveals hidden recurrent and individual genomic copy number alterations with prognostic relevance A Thiel1, M Beier1, D Ingenhag1, K Servan1, M Hein1, V Moeller1, B Betz1, B Hildebrandt1, C Evers1,3, U Germing2 and B Royer-Pokora1 1Institute of Human Genetics and Anthropology, Medical Faculty, Heinrich Heine University, Duesseldorf, Germany and 2Department of Hematology, Oncology and Clinical Immunology, Heinrich Heine University, Duesseldorf, Germany About 40% of patients with myelodysplastic syndromes (MDSs) 40–50% of MDS cases have a normal karyotype. MDS patients present with a normal karyotype, and they are facing different with a normal karyotype and low-risk clinical parameters are courses of disease. To advance the biological understanding often assigned into the IPSS low and intermediate-1 risk groups. and to find molecular prognostic markers, we performed a high- resolution oligonucleotide array study of 107 MDS patients In the absence of genetic or biological markers, prognostic (French American British) with a normal karyotype and clinical stratification of these patients is difficult. To better prognosticate follow-up through the Duesseldorf MDS registry. Recurrent these patients, new parameters to identify patients at higher risk hidden deletions overlapping with known cytogenetic aberra- are urgently needed. With the more recently introduced modern tions or sites of known tumor-associated genes were identi- technologies of whole-genome-wide surveys of genetic aberra- fied in 4q24 (TET2, 2x), 5q31.2 (2x), 7q22.1 (3x) and 21q22.12 tions, it is hoped that more insights into the biology of disease (RUNX1, 2x).
    [Show full text]
  • Supplementary Table S4. FGA Co-Expressed Gene List in LUAD
    Supplementary Table S4. FGA co-expressed gene list in LUAD tumors Symbol R Locus Description FGG 0.919 4q28 fibrinogen gamma chain FGL1 0.635 8p22 fibrinogen-like 1 SLC7A2 0.536 8p22 solute carrier family 7 (cationic amino acid transporter, y+ system), member 2 DUSP4 0.521 8p12-p11 dual specificity phosphatase 4 HAL 0.51 12q22-q24.1histidine ammonia-lyase PDE4D 0.499 5q12 phosphodiesterase 4D, cAMP-specific FURIN 0.497 15q26.1 furin (paired basic amino acid cleaving enzyme) CPS1 0.49 2q35 carbamoyl-phosphate synthase 1, mitochondrial TESC 0.478 12q24.22 tescalcin INHA 0.465 2q35 inhibin, alpha S100P 0.461 4p16 S100 calcium binding protein P VPS37A 0.447 8p22 vacuolar protein sorting 37 homolog A (S. cerevisiae) SLC16A14 0.447 2q36.3 solute carrier family 16, member 14 PPARGC1A 0.443 4p15.1 peroxisome proliferator-activated receptor gamma, coactivator 1 alpha SIK1 0.435 21q22.3 salt-inducible kinase 1 IRS2 0.434 13q34 insulin receptor substrate 2 RND1 0.433 12q12 Rho family GTPase 1 HGD 0.433 3q13.33 homogentisate 1,2-dioxygenase PTP4A1 0.432 6q12 protein tyrosine phosphatase type IVA, member 1 C8orf4 0.428 8p11.2 chromosome 8 open reading frame 4 DDC 0.427 7p12.2 dopa decarboxylase (aromatic L-amino acid decarboxylase) TACC2 0.427 10q26 transforming, acidic coiled-coil containing protein 2 MUC13 0.422 3q21.2 mucin 13, cell surface associated C5 0.412 9q33-q34 complement component 5 NR4A2 0.412 2q22-q23 nuclear receptor subfamily 4, group A, member 2 EYS 0.411 6q12 eyes shut homolog (Drosophila) GPX2 0.406 14q24.1 glutathione peroxidase
    [Show full text]
  • Mouse Epb42 Knockout Project (CRISPR/Cas9)
    https://www.alphaknockout.com Mouse Epb42 Knockout Project (CRISPR/Cas9) Objective: To create a Epb42 knockout Mouse model (C57BL/6J) by CRISPR/Cas-mediated genome engineering. Strategy summary: The Epb42 gene (NCBI Reference Sequence: NM_013513 ; Ensembl: ENSMUSG00000023216 ) is located on Mouse chromosome 2. 13 exons are identified, with the ATG start codon in exon 1 and the TAA stop codon in exon 13 (Transcript: ENSMUST00000102490). Exon 2~5 will be selected as target site. Cas9 and gRNA will be co-injected into fertilized eggs for KO Mouse production. The pups will be genotyped by PCR followed by sequencing analysis. Note: Homozygotes for a targeted null mutation exhibit erythrocytic abnormalities including mild spherocytosis, altered ion transport, and dehydration. Exon 2 starts from about 0.53% of the coding region. Exon 2~5 covers 31.07% of the coding region. The size of effective KO region: ~5547 bp. The KO region does not have any other known gene. Page 1 of 9 https://www.alphaknockout.com Overview of the Targeting Strategy Wildtype allele 5' gRNA region gRNA region 3' 1 2 3 4 5 13 Legends Exon of mouse Epb42 Knockout region Page 2 of 9 https://www.alphaknockout.com Overview of the Dot Plot (up) Window size: 15 bp Forward Reverse Complement Sequence 12 Note: The 1985 bp section upstream of Exon 2 is aligned with itself to determine if there are tandem repeats. No significant tandem repeat is found in the dot plot matrix. So this region is suitable for PCR screening or sequencing analysis. Overview of the Dot Plot (down) Window size: 15 bp Forward Reverse Complement Sequence 12 Note: The 1230 bp section downstream of Exon 5 is aligned with itself to determine if there are tandem repeats.
    [Show full text]
  • Supplementary Material
    BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) J Neurol Neurosurg Psychiatry Page 1 / 45 SUPPLEMENTARY MATERIAL Appendix A1: Neuropsychological protocol. Appendix A2: Description of the four cases at the transitional stage. Table A1: Clinical status and center proportion in each batch. Table A2: Complete output from EdgeR. Table A3: List of the putative target genes. Table A4: Complete output from DIANA-miRPath v.3. Table A5: Comparison of studies investigating miRNAs from brain samples. Figure A1: Stratified nested cross-validation. Figure A2: Expression heatmap of miRNA signature. Figure A3: Bootstrapped ROC AUC scores. Figure A4: ROC AUC scores with 100 different fold splits. Figure A5: Presymptomatic subjects probability scores. Figure A6: Heatmap of the level of enrichment in KEGG pathways. Kmetzsch V, et al. J Neurol Neurosurg Psychiatry 2021; 92:485–493. doi: 10.1136/jnnp-2020-324647 BMJ Publishing Group Limited (BMJ) disclaims all liability and responsibility arising from any reliance Supplemental material placed on this supplemental material which has been supplied by the author(s) J Neurol Neurosurg Psychiatry Appendix A1. Neuropsychological protocol The PREV-DEMALS cognitive evaluation included standardized neuropsychological tests to investigate all cognitive domains, and in particular frontal lobe functions. The scores were provided previously (Bertrand et al., 2018). Briefly, global cognitive efficiency was evaluated by means of Mini-Mental State Examination (MMSE) and Mattis Dementia Rating Scale (MDRS). Frontal executive functions were assessed with Frontal Assessment Battery (FAB), forward and backward digit spans, Trail Making Test part A and B (TMT-A and TMT-B), Wisconsin Card Sorting Test (WCST), and Symbol-Digit Modalities test.
    [Show full text]
  • Nuclear PTEN Safeguards Pre-Mrna Splicing to Link Golgi Apparatus for Its Tumor Suppressive Role
    ARTICLE DOI: 10.1038/s41467-018-04760-1 OPEN Nuclear PTEN safeguards pre-mRNA splicing to link Golgi apparatus for its tumor suppressive role Shao-Ming Shen1, Yan Ji2, Cheng Zhang1, Shuang-Shu Dong2, Shuo Yang1, Zhong Xiong1, Meng-Kai Ge1, Yun Yu1, Li Xia1, Meng Guo1, Jin-Ke Cheng3, Jun-Ling Liu1,3, Jian-Xiu Yu1,3 & Guo-Qiang Chen1 Dysregulation of pre-mRNA alternative splicing (AS) is closely associated with cancers. However, the relationships between the AS and classic oncogenes/tumor suppressors are 1234567890():,; largely unknown. Here we show that the deletion of tumor suppressor PTEN alters pre-mRNA splicing in a phosphatase-independent manner, and identify 262 PTEN-regulated AS events in 293T cells by RNA sequencing, which are associated with significant worse outcome of cancer patients. Based on these findings, we report that nuclear PTEN interacts with the splicing machinery, spliceosome, to regulate its assembly and pre-mRNA splicing. We also identify a new exon 2b in GOLGA2 transcript and the exon exclusion contributes to PTEN knockdown-induced tumorigenesis by promoting dramatic Golgi extension and secretion, and PTEN depletion significantly sensitizes cancer cells to secretion inhibitors brefeldin A and golgicide A. Our results suggest that Golgi secretion inhibitors alone or in combination with PI3K/Akt kinase inhibitors may be therapeutically useful for PTEN-deficient cancers. 1 Department of Pathophysiology, Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Shanghai Jiao Tong University School of Medicine (SJTU-SM), Shanghai 200025, China. 2 Institute of Health Sciences, Shanghai Institutes for Biological Sciences of Chinese Academy of Sciences and SJTU-SM, Shanghai 200025, China.
    [Show full text]