DOCSLIB.ORG

S41467-019-13840-9.Pdf

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
S41467-019-13840-9.Pdf ARTICLE https://doi.org/10.1038/s41467-019-13840-9 OPEN Accurate quantification of circular RNAs identifies extensive circular isoform switching events Jinyang Zhang 1,2, Shuai Chen1, Jingwen Yang1 & Fangqing Zhao1,2,3* Detection and quantification of circular RNAs (circRNAs) face several significant challenges, including high false discovery rate, uneven rRNA depletion and RNase R treatment efficiency, and underestimation of back-spliced junction reads. Here, we propose a novel algorithm, fi 1234567890():,; CIRIquant, for accurate circRNA quanti cation and differential expression analysis. By con- structing pseudo-circular reference for re-alignment of RNA-seq reads and employing sophisticated statistical models to correct RNase R treatment biases, CIRIquant can provide more accurate expression values for circRNAs with significantly reduced false discovery rate. We further develop a one-stop differential expression analysis pipeline implementing two independent measures, which helps unveil the regulation of competitive splicing between circRNAs and their linear counterparts. We apply CIRIquant to RNA-seq datasets of hepa- tocellular carcinoma, and characterize two important groups of linear-circular switching and circular transcript usage switching events, which demonstrate the promising ability to explore extensive transcriptomic changes in liver tumorigenesis. 1 Computational Genomics Lab, Beijing Institutes of Life Science, Chinese Academy of Sciences, 100101 Beijing, China. 2 University of Chinese Academy of Sciences, 100049 Beijing, China. 3 Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, 650223 Kunming, China. *email: [email protected] NATURE COMMUNICATIONS | (2020) 11:90 | https://doi.org/10.1038/s41467-019-13840-9 | www.nature.com/naturecommunications 1 ARTICLE NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-019-13840-9 ircular RNA (circRNA) is a large class of RNA molecules Additionally, we profile the switching events in circRNA junction that contain a covalent circular structure formed by non- ratio and circular transcript usage, and characterize these two C ′ ′ canonical 3 to 5 end-joining event called back-splicing. groups of circRNAs with potential biological functions. We Previous studies have shown that circRNAs are widely present believe that CIRIquant, which provides an accurate and efficient and some are conserved in eukaryotic organisms, and have a quantification approach to characterize circRNAs and perform relative low abundance compared with canonical linear mRNA differential expression analysis after correcting experimental and transcripts1. Although the exact function of most circRNAs is still computational biases, will greatly improve our understanding of ambiguous, studies have shown that circRNAs may function as circRNA diversities and functions. sponges to sequester miRNAs or RNA binding proteins (RBPs)2–4. More recently, several studies found that internal ribosome entry sites or N6-methyladenosine can promote the Results translation of circRNAs, which results in the biogenesis of Challenges in quantifying the expression of circRNAs. To rig- circRNA-derived proteins5–7. Thus, circRNAs have great poten- orously evaluate the challenges in current quantification of cir- tials to play important roles in cellular metabolic process, and the cRNA expression, we collected 63 transcriptomic samples from characterization and quantification of circRNAs from high- six previous studies16–22, including both RiboMinus and Ribo- throughput RNA-seq data has become an emerging problem in Minus/RNase R RNA-seq libraries of four species (human, circRNA studies. mouse, fly and roundworm, Supplementary Table 1). All these A number of computational methods have been developed for RNA-seq datasets were aligned to their reference genomes and characterizing circRNAs8. Most of these methods employ align- the ribosomal RNA (rRNA) sequences using HISAT223 to assess ment based strategies to recognize back-spliced junction (BSJ) the mapping rate and rRNA sequence fraction. As shown in reads from circRNAs, and have limited sensitivity and notable Fig. 1a, RNA-seq datasets from four species showed an extra- false-positive rate in circRNA identification9. Recently, a model- ordinarily high variance in the efficiency of rRNA sequence based strategy is employed by Sailfish-cir10, which uses a quasi- depletion, which is largely due to the limited specificity and mapping method to acquire direct estimation of circular tran- efficiency of current RiboMinus transcriptome isolation kit. For script expression. This statistical model depends on the unique each sample, the raw reads were further aligned to the reference sequence between circular and linear transcripts, which limits its genome using the BWA-MEM algorithm24, and then subjected to ability on quantifying exonic circRNAs. Moreover, the ratio of CIRI2 for circRNA detection and quantification. We chose CPM BSJ reads to canonical linear reads at the junction, which repre- (counts per million mapped reads) to represent circRNA sents the splicing preference in precursor mRNAs, is also an expression levels to remove the biases derived from different important factor in circRNA analysis. However, among all cur- library insert size and sequencing depth. For 34 pairs of Ribo- rently available circRNA detection tools, only CIRI211,12 can Minus and RiboMinus/RNase R samples, the overall number of output the junction ratio directly. DCC13 and Sailfish-cir estimate detected circRNAs ranged from several thousand to over 30 the expression values of both linear and circular transcripts, thousand, in which most of circRNAs detected in the RiboMinus which can be used to calculate circRNA’s junction ratio. Hence, a samples can be validated in the RiboMinus/RNase R samples reliable computational tool is urgently needed for accurate (Fig. 1b). Although the total number of identified circRNAs quantification of circRNAs and their parental linear transcripts. increased after RNase R enrichment, over 50–80% highly Differential expression analyses of circRNAs in different sam- expressed circRNAs could be detected in both RiboMinus and ples is a routine analysis in circRNA studies. Currently, simple RiboMinus/RNase R samples. Considering that RNase R treat- statistical tests (e.g., t-test) or differential expression analysis ment is widely used for circRNA enrichment, we compared pipelines designed for linear RNA transcripts (e.g., DESeq214) are expression levels of both gene and circRNA in RNase R-treated often used to evaluate the significance of differentially expressed and untreated samples to investigate whether RNase R treatment circRNAs. Since most of circRNAs are expressed at extremely low can introduce bias in transcript quantification. As shown in levels, RNase R treatment are usually performed to enrich cir- Fig. 1b (right), gene expression levels decreased over 2-fold, which cRNAs. For studies using RiboMinus/RNase R-treated RNA-seq is in concordance with the degradation effect of linear RNAs in libraries, the variation of enrichment coefficient in the RNase R RNase R treatment. However, the relative expression value of treatment step may result in biased estimation of circRNA circRNAs also showed a certain level of reduction. Considering expression levels. In addition, analysis of experimental replicates that the RNase R treatment is expected to increase the saturation has shown the poor reproducibility of circRNA identification15, level of circRNA detection, the more circRNAs are identified, the which further indicates that accurate characterization and quan- more circular reads and the lower the relative expression value for tification of circular transcripts is crucial in circRNA studies. each specific circRNA. To overcome these limitations, we propose CIRIquant for To further investigate the effect of RNase R treatment on accurate quantification of both circRNAs and their parental circRNA quantification, we divided the detected circRNAs into transcripts and filtration of false positives BSJ reads. CIRIquant two groups, highly expressed circRNAs and lowly expressed can employ currently widely used tools (e.g., CIRI211,12, CIR- circRNAs, according to the ranking of their expression values. Cexplorer216, find_circ2, etc.) for circRNA identification, then Then, we calculated the proportion of reads from the two groups generated pseudo-reference sequences for the identified circRNA of circRNAs after RNase R treatment. As shown in Fig. 1c and transcripts to re-align putative BSJ reads. Using both alignment Supplementary Fig. 1, we observed an increased proportion of BSJ results of reads against the reference genome and pseudo cir- reads for highly expressed circRNAs after RNase R treatment, cRNA transcripts, CIRIquant not only achieves more accurate indicating that highly expressed circRNAs tended to be enriched and sensitive identification of BSJ reads, but also enables reliable in RiboMinus/RNase R-treated data. In order to assess the quantification of junction ratio for circRNAs. Moreover, CIRI- efficiency of RNase R treatment, we calculated the enrichment quant provides a convenient function for one-stop circRNA coefficient for circRNAs in different RNA-seq datasets. Surpris- differential expression analysis. We apply CIRIquant to survey ingly, the enrichment coefficient exhibited distinct difference circRNA expression profiles between hepatocellular carcinoma among samples from different species or even within the same (HCC) tumor samples and their adjacent normal tissues, which
Load more

Recommended publications
  • A Flexible Microfluidic System for Single-Cell Transcriptome Profiling
    www.nature.com/scientificreports OPEN A fexible microfuidic system for single‑cell transcriptome profling elucidates phased transcriptional regulators of cell cycle Karen Davey1,7, Daniel Wong2,7, Filip Konopacki2, Eugene Kwa1, Tony Ly3, Heike Fiegler2 & Christopher R. Sibley 1,4,5,6* Single cell transcriptome profling has emerged as a breakthrough technology for the high‑resolution understanding of complex cellular systems. Here we report a fexible, cost‑efective and user‑ friendly droplet‑based microfuidics system, called the Nadia Instrument, that can allow 3′ mRNA capture of ~ 50,000 single cells or individual nuclei in a single run. The precise pressure‑based system demonstrates highly reproducible droplet size, low doublet rates and high mRNA capture efciencies that compare favorably in the feld. Moreover, when combined with the Nadia Innovate, the system can be transformed into an adaptable setup that enables use of diferent bufers and barcoded bead confgurations to facilitate diverse applications. Finally, by 3′ mRNA profling asynchronous human and mouse cells at diferent phases of the cell cycle, we demonstrate the system’s ability to readily distinguish distinct cell populations and infer underlying transcriptional regulatory networks. Notably this provided supportive evidence for multiple transcription factors that had little or no known link to the cell cycle (e.g. DRAP1, ZKSCAN1 and CEBPZ). In summary, the Nadia platform represents a promising and fexible technology for future transcriptomic studies, and other related applications, at cell resolution. Single cell transcriptome profling has recently emerged as a breakthrough technology for understanding how cellular heterogeneity contributes to complex biological systems. Indeed, cultured cells, microorganisms, biopsies, blood and other tissues can be rapidly profled for quantifcation of gene expression at cell resolution.
    [Show full text]
  • Transposon Activation Mutagenesis As a Screening Tool for Identifying
    Chen et al. BMC Cancer 2013, 13:93 http://www.biomedcentral.com/1471-2407/13/93 RESEARCH ARTICLE Open Access Transposon activation mutagenesis as a screening tool for identifying resistance to cancer therapeutics Li Chen1,2*, Lynda Stuart2, Toshiro K Ohsumi3, Shawn Burgess4, Gaurav K Varshney4, Anahita Dastur1, Mark Borowsky3, Cyril Benes1, Adam Lacy-Hulbert2 and Emmett V Schmidt2 Abstract Background: The development of resistance to chemotherapies represents a significant barrier to successful cancer treatment. Resistance mechanisms are complex, can involve diverse and often unexpected cellular processes, and can vary with both the underlying genetic lesion and the origin or type of tumor. For these reasons developing experimental strategies that could be used to understand, identify and predict mechanisms of resistance in different malignant cells would be a major advance. Methods: Here we describe a gain-of-function forward genetic approach for identifying mechanisms of resistance. This approach uses a modified piggyBac transposon to generate libraries of mutagenized cells, each containing transposon insertions that randomly activate nearby gene expression. Genes of interest are identified using next- gen high-throughput sequencing and barcode multiplexing is used to reduce experimental cost. Results: Using this approach we successfully identify genes involved in paclitaxel resistance in a variety of cancer cell lines, including the multidrug transporter ABCB1, a previously identified major paclitaxel resistance gene. Analysis of co-occurring transposons integration sites in single cell clone allows for the identification of genes that might act cooperatively to produce drug resistance a level of information not accessible using RNAi or ORF expression screening approaches. Conclusion: We have developed a powerful pipeline to systematically discover drug resistance in mammalian cells in vitro.
    [Show full text]
  • WO 2019/079361 Al 25 April 2019 (25.04.2019) W 1P O PCT
    (12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) (19) World Intellectual Property Organization I International Bureau (10) International Publication Number (43) International Publication Date WO 2019/079361 Al 25 April 2019 (25.04.2019) W 1P O PCT (51) International Patent Classification: CA, CH, CL, CN, CO, CR, CU, CZ, DE, DJ, DK, DM, DO, C12Q 1/68 (2018.01) A61P 31/18 (2006.01) DZ, EC, EE, EG, ES, FI, GB, GD, GE, GH, GM, GT, HN, C12Q 1/70 (2006.01) HR, HU, ID, IL, IN, IR, IS, JO, JP, KE, KG, KH, KN, KP, KR, KW, KZ, LA, LC, LK, LR, LS, LU, LY, MA, MD, ME, (21) International Application Number: MG, MK, MN, MW, MX, MY, MZ, NA, NG, NI, NO, NZ, PCT/US2018/056167 OM, PA, PE, PG, PH, PL, PT, QA, RO, RS, RU, RW, SA, (22) International Filing Date: SC, SD, SE, SG, SK, SL, SM, ST, SV, SY, TH, TJ, TM, TN, 16 October 2018 (16. 10.2018) TR, TT, TZ, UA, UG, US, UZ, VC, VN, ZA, ZM, ZW. (25) Filing Language: English (84) Designated States (unless otherwise indicated, for every kind of regional protection available): ARIPO (BW, GH, (26) Publication Language: English GM, KE, LR, LS, MW, MZ, NA, RW, SD, SL, ST, SZ, TZ, (30) Priority Data: UG, ZM, ZW), Eurasian (AM, AZ, BY, KG, KZ, RU, TJ, 62/573,025 16 October 2017 (16. 10.2017) US TM), European (AL, AT, BE, BG, CH, CY, CZ, DE, DK, EE, ES, FI, FR, GB, GR, HR, HU, ΓΕ , IS, IT, LT, LU, LV, (71) Applicant: MASSACHUSETTS INSTITUTE OF MC, MK, MT, NL, NO, PL, PT, RO, RS, SE, SI, SK, SM, TECHNOLOGY [US/US]; 77 Massachusetts Avenue, TR), OAPI (BF, BJ, CF, CG, CI, CM, GA, GN, GQ, GW, Cambridge, Massachusetts 02139 (US).
    [Show full text]
  • Algorithms for Transcriptome Quantification and Reconstruction from RNA-Seq Data
    Georgia State University ScholarWorks @ Georgia State University Computer Science Dissertations Department of Computer Science Fall 11-16-2012 Algorithms for Transcriptome Quantification and Reconstruction from RNA-Seq Data Serghei Mangul Follow this and additional works at: https://scholarworks.gsu.edu/cs_diss Recommended Citation Mangul, Serghei, "Algorithms for Transcriptome Quantification and Reconstruction from RNA-Seq Data." Dissertation, Georgia State University, 2012. https://scholarworks.gsu.edu/cs_diss/71 This Dissertation is brought to you for free and open access by the Department of Computer Science at ScholarWorks @ Georgia State University. It has been accepted for inclusion in Computer Science Dissertations by an authorized administrator of ScholarWorks @ Georgia State University. For more information, please contact [email protected]. ALGORITHMS FOR TRANSCRIPTOME QUANTIFICATION AND RECONSTRUCTION FROM RNA-SEQ DATA by SERGHEI MANGUL Under the Direction of Dr. Alexander Zelikovsky ABSTRACT Massively parallel whole transcriptome sequencing and its ability to generate full tran- scriptome data at the single transcript level provides a powerful tool with multiple inter- related applications, including transcriptome reconstruction, gene/isoform expression esti- mation, also known as transcriptome quantification. As a result, whole transcriptome se- quencing has become the technology of choice for performing transcriptome analysis, rapidly replacing array-based technologies. The most commonly used transcriptome sequencing pro- tocol, referred to as RNA-Seq, generates short (single or paired) sequencing tags from the ends of randomly generated cDNA fragments. RNA-Seq protocol reduces the sequencing cost and significantly increases data throughput, but is computationally challenging to re- construct full-length transcripts and accurately estimate their abundances across all cell types. We focus on two main problems in transcriptome data analysis, namely, transcriptome reconstruction and quantification.
    [Show full text]
  • The Function and Evolution of C2H2 Zinc Finger Proteins and Transposons
    The function and evolution of C2H2 zinc finger proteins and transposons by Laura Francesca Campitelli A thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Department of Molecular Genetics University of Toronto © Copyright by Laura Francesca Campitelli 2020 The function and evolution of C2H2 zinc finger proteins and transposons Laura Francesca Campitelli Doctor of Philosophy Department of Molecular Genetics University of Toronto 2020 Abstract Transcription factors (TFs) confer specificity to transcriptional regulation by binding specific DNA sequences and ultimately affecting the ability of RNA polymerase to transcribe a locus. The C2H2 zinc finger proteins (C2H2 ZFPs) are a TF class with the unique ability to diversify their DNA-binding specificities in a short evolutionary time. C2H2 ZFPs comprise the largest class of TFs in Mammalian genomes, including nearly half of all Human TFs (747/1,639). Positive selection on the DNA-binding specificities of C2H2 ZFPs is explained by an evolutionary arms race with endogenous retroelements (EREs; copy-and-paste transposable elements), where the C2H2 ZFPs containing a KRAB repressor domain (KZFPs; 344/747 Human C2H2 ZFPs) are thought to diversify to bind new EREs and repress deleterious transposition events. However, evidence of the gain and loss of KZFP binding sites on the ERE sequence is sparse due to poor resolution of ERE sequence evolution, despite the recent publication of binding preferences for 242/344 Human KZFPs. The goal of my doctoral work has been to characterize the Human C2H2 ZFPs, with specific interest in their evolutionary history, functional diversity, and coevolution with LINE EREs.
    [Show full text]
  • Muscle Enriched Lamin Interacting Protein (Mlip) Binds Chromatin and Is Required for Myoblast Differentiation
    cells Article Muscle Enriched Lamin Interacting Protein (Mlip) Binds Chromatin and Is Required for Myoblast Differentiation Elmira Ahmady 1,2, Alexandre Blais 1,3,4 and Patrick G. Burgon 5,* 1 Department of Biochemistry, Microbiology, and Immunology, University of Ottawa, Ottawa, ON K1H 8M5, Canada; [email protected] (E.A.); [email protected] (A.B.) 2 Molecular Signaling Laboratory, University of Ottawa Heart Institute, Ottawa, ON K1Y 4W7, Canada 3 Ottawa Institute of Systems Biology, University of Ottawa, Ottawa, ON K1H 8M5, Canada 4 University of Ottawa Centre for Infection, Immunity and Inflammation (CI3), Ottawa, ON K1H 8M5, Canada 5 Department of Chemistry and Earth Sciences, College of Arts and Sciences, Qatar University, Doha 999043, Qatar * Correspondence: [email protected] Abstract: Muscle-enriched A-type lamin-interacting protein (Mlip) is a recently discovered Amniota gene that encodes proteins of unknown biological function. Here we report Mlip’s direct interaction with chromatin, and it may function as a transcriptional co-factor. Chromatin immunoprecipitations with microarray analysis demonstrated a propensity for Mlip to associate with genomic regions in close proximity to genes that control tissue-specific differentiation. Gel mobility shift assays confirmed that Mlip protein complexes with genomic DNA. Blocking Mlip expression in C2C12 myoblasts down- regulates myogenic regulatory factors (MyoD and MyoG) and subsequently significantly inhibits myogenic differentiation and the formation of myotubes. Collectively our data demonstrate that Mlip is required for C2C12 myoblast differentiation into myotubes. Mlip may exert this role as a transcriptional regulator of a myogenic program that is unique to amniotes. Citation: Ahmady, E.; Blais, A.; Burgon, P.G.
    [Show full text]
  • Whole Genome Sequencing of Familial Non-Medullary Thyroid Cancer Identifies Germline Alterations in MAPK/ERK and PI3K/AKT Signaling Pathways
    biomolecules Article Whole Genome Sequencing of Familial Non-Medullary Thyroid Cancer Identifies Germline Alterations in MAPK/ERK and PI3K/AKT Signaling Pathways Aayushi Srivastava 1,2,3,4 , Abhishek Kumar 1,5,6 , Sara Giangiobbe 1, Elena Bonora 7, Kari Hemminki 1, Asta Försti 1,2,3 and Obul Reddy Bandapalli 1,2,3,* 1 Division of Molecular Genetic Epidemiology, German Cancer Research Center (DKFZ), D-69120 Heidelberg, Germany; [email protected] (A.S.); [email protected] (A.K.); [email protected] (S.G.); [email protected] (K.H.); [email protected] (A.F.) 2 Hopp Children’s Cancer Center (KiTZ), D-69120 Heidelberg, Germany 3 Division of Pediatric Neurooncology, German Cancer Research Center (DKFZ), German Cancer Consortium (DKTK), D-69120 Heidelberg, Germany 4 Medical Faculty, Heidelberg University, D-69120 Heidelberg, Germany 5 Institute of Bioinformatics, International Technology Park, Bangalore 560066, India 6 Manipal Academy of Higher Education (MAHE), Manipal, Karnataka 576104, India 7 S.Orsola-Malphigi Hospital, Unit of Medical Genetics, 40138 Bologna, Italy; [email protected] * Correspondence: [email protected]; Tel.: +49-6221-42-1709 Received: 29 August 2019; Accepted: 10 October 2019; Published: 13 October 2019 Abstract: Evidence of familial inheritance in non-medullary thyroid cancer (NMTC) has accumulated over the last few decades. However, known variants account for a very small percentage of the genetic burden. Here, we focused on the identification of common pathways and networks enriched in NMTC families to better understand its pathogenesis with the final aim of identifying one novel high/moderate-penetrance germline predisposition variant segregating with the disease in each studied family.
    [Show full text]
  • Supplementary Table 1 Genes Tested in Qrt-PCR in Nfpas
    Supplementary Table 1 Genes tested in qRT-PCR in NFPAs Gene Bank accession Gene Description number ABI assay ID a disintegrin-like and metalloprotease with thrombospondin type 1 motif 7 ADAMTS7 NM_014272.3 Hs00276223_m1 Rho guanine nucleotide exchange factor (GEF) 3 ARHGEF3 NM_019555.1 Hs00219609_m1 BCL2-associated X protein BAX NM_004324 House design Bcl-2 binding component 3 BBC3 NM_014417.2 Hs00248075_m1 B-cell CLL/lymphoma 2 BCL2 NM_000633 House design Bone morphogenetic protein 7 BMP7 NM_001719.1 Hs00233476_m1 CCAAT/enhancer binding protein (C/EBP), alpha CEBPA NM_004364.2 Hs00269972_s1 coxsackie virus and adenovirus receptor CXADR NM_001338.3 Hs00154661_m1 Homo sapiens Dicer1, Dcr-1 homolog (Drosophila) (DICER1) DICER1 NM_177438.1 Hs00229023_m1 Homo sapiens dystonin DST NM_015548.2 Hs00156137_m1 fms-related tyrosine kinase 3 FLT3 NM_004119.1 Hs00174690_m1 glutamate receptor, ionotropic, N-methyl D-aspartate 1 GRIN1 NM_000832.4 Hs00609557_m1 high-mobility group box 1 HMGB1 NM_002128.3 Hs01923466_g1 heterogeneous nuclear ribonucleoprotein U HNRPU NM_004501.3 Hs00244919_m1 insulin-like growth factor binding protein 5 IGFBP5 NM_000599.2 Hs00181213_m1 latent transforming growth factor beta binding protein 4 LTBP4 NM_001042544.1 Hs00186025_m1 microtubule-associated protein 1 light chain 3 beta MAP1LC3B NM_022818.3 Hs00797944_s1 matrix metallopeptidase 17 MMP17 NM_016155.4 Hs01108847_m1 myosin VA MYO5A NM_000259.1 Hs00165309_m1 Homo sapiens nuclear factor (erythroid-derived 2)-like 1 NFE2L1 NM_003204.1 Hs00231457_m1 oxoglutarate (alpha-ketoglutarate)
    [Show full text]
  • Population-Haplotype Models for Mapping and Tagging Structural Variation Using Whole Genome Sequencing
    Population-haplotype models for mapping and tagging structural variation using whole genome sequencing Eleni Loizidou Submitted in part fulfilment of the requirements for the degree of Doctor of Philosophy Section of Genomics of Common Disease Department of Medicine Imperial College London, 2018 1 Declaration of originality I hereby declare that the thesis submitted for a Doctor of Philosophy degree is based on my own work. Proper referencing is given to the organisations/cohorts I collaborated with during the project. 2 Copyright Declaration The copyright of this thesis rests with the author and is made available under a Creative Commons Attribution Non-Commercial No Derivatives licence. Researchers are free to copy, distribute or transmit the thesis on the condition that they attribute it, that they do not use it for commercial purposes and that they do not alter, transform or build upon it. For any reuse or redistribution, researchers must make clear to others the licence terms of this work 3 Abstract The scientific interest in copy number variation (CNV) is rapidly increasing, mainly due to the evidence of phenotypic effects and its contribution to disease susceptibility. Single nucleotide polymorphisms (SNPs) which are abundant in the human genome have been widely investigated in genome-wide association studies (GWAS). Despite the notable genomic effects both CNVs and SNPs have, the correlation between them has been relatively understudied. In the past decade, next generation sequencing (NGS) has been the leading high-throughput technology for investigating CNVs and offers mapping at a high-quality resolution. We created a map of NGS-defined CNVs tagged by SNPs using the 1000 Genomes Project phase 3 (1000G) sequencing data to examine patterns between the two types of variation in protein-coding genes.
    [Show full text]
  • Supplemental Information
    SUPPLEMENTARY APPENDIX Plerixafor and G-CSF combination mobilizes hematopoietic stem and progenitors cells with a distinct transcriptional profile and a reduced in vivo homing capacity compared to plerixafor alone Maria Rosa Lidonnici, 1,2 * Annamaria Aprile, 1,3* Marta Claudia Frittoli, 4 Giacomo Mandelli, 1 Ylenia Paleari, 1,2 Antonello Spinelli, 5 Bernhard Gentner, 4 Matilde Zambelli, 6 Cristina Parisi, 6 Laura Bellio, 6 Elena Cassinerio, 7 Laura Zanaboni, 7 Maria Domenica Cappellini, 7 Fabio Ciceri, 4 Sarah Marktel 4 and Giuliana Ferrari 1,2 *MRL and AA contributed equally to this work. 1San Raffaele-Telethon Institute for Gene Therapy (SR-TIGET), Division of Regenerative Medicine, Stem Cells and Gene Therapy, IRCCS San Raffaele Sci - entific Institute, Milan; 2Vita-Salute San Raffaele University, Milan; 3Università degli Studi di Roma “Tor Vergata”; 4Hematology and Bone Marrow Transplan - tation Unit, IRCCS San Raffaele Scientific Institute, Milan; 5Experimental Imaging Centre, San Raffaele Scientific Institute, Milan; 6Immunohematology and Transfusion Medicine Unit, IRCCS San Raffaele Scientific Institute, Milan and 7Università di Milano, Ca Granda Foundation IRCCS, Italy Correspondence: [email protected] doi:10.3324/haematol.2016.154740 SUPPLEMENTARY FIGURES Fig. S1. Gene Ontology analysis of differentially expressed genes in Plx PB. Functional annotation by Partek software showed that genes could be grouped into a limited number of biological categories. Most of probesets upregulated in Plerixafor mobilized cells were enriched in cellular organization and metabolic process (p value <0.05) (A-B), while most of the probesets down regulated in Plx PB are involved in cell cycle process (p value <0.001) (C-D). One-way ANOVA statistical analysis was performed.
    [Show full text]
  • Evolutionary Analysis of Viral Sequences in Eukaryotic Genomes
    Evolutionary analysis of viral sequences in eukaryotic genomes Sean Schneider A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy University of Washington 2014 Reading Committee: James H. Thomas, Chair Willie Swanson Phil Green Program Authorized to Offer Degree: Genome Sciences ©Copyright 2014 Sean Schneider University of Washington Abstract Evolutionary analysis of viral sequences in eukaryotic genomes Sean Schneider Chair of the supervisory committee: Professor James H. Thomas Genome Sciences The focus of this work is several evolutionary analyses of endogenous viral sequences in eukaryotic genomes. Endogenous viral sequences can provide key insights into the past forms and evolutionary history of viruses, as well as the responses of host organisms they infect. In this work I have examined viral sequences in a diverse assortment of eukaryotic hosts in order to study coevolution between hosts and the organisms that infect them. This research consisted of two major lines of investigation. In the first portion of this work, I outline the hypothesis that the C2H2 zinc finger gene family in vertebrates has evolved by birth-death evolution in response to sporadic retroviral infection. The hypothesis suggests an evolutionary model in which newly duplicated zinc finger genes are retained by selection in response to retroviral infection. This hypothesis is supported by a strong association (R2=0.67) between the number of endogenous retroviruses and the number of zinc fingers in diverse vertebrate genomes. Based on this and other evidence, the zinc finger gene family appears to act as a “genomic immune system” against retroviral infections. The other major line of investigation in this work examines endogenous virus sequences utilized by parasitic wasps to disable hosts that they infect.
    [Show full text]
  • Characterizing Genomic Duplication in Autism Spectrum Disorder by Edward James Higginbotham a Thesis Submitted in Conformity
    Characterizing Genomic Duplication in Autism Spectrum Disorder by Edward James Higginbotham A thesis submitted in conformity with the requirements for the degree of Master of Science Graduate Department of Molecular Genetics University of Toronto © Copyright by Edward James Higginbotham 2020 i Abstract Characterizing Genomic Duplication in Autism Spectrum Disorder Edward James Higginbotham Master of Science Graduate Department of Molecular Genetics University of Toronto 2020 Duplication, the gain of additional copies of genomic material relative to its ancestral diploid state is yet to achieve full appreciation for its role in human traits and disease. Challenges include accurately genotyping, annotating, and characterizing the properties of duplications, and resolving duplication mechanisms. Whole genome sequencing, in principle, should enable accurate detection of duplications in a single experiment. This thesis makes use of the technology to catalogue disease relevant duplications in the genomes of 2,739 individuals with Autism Spectrum Disorder (ASD) who enrolled in the Autism Speaks MSSNG Project. Fine-mapping the breakpoint junctions of 259 ASD-relevant duplications identified 34 (13.1%) variants with complex genomic structures as well as tandem (193/259, 74.5%) and NAHR- mediated (6/259, 2.3%) duplications. As whole genome sequencing-based studies expand in scale and reach, a continued focus on generating high-quality, standardized duplication data will be prerequisite to addressing their associated biological mechanisms. ii Acknowledgements I thank Dr. Stephen Scherer for his leadership par excellence, his generosity, and for giving me a chance. I am grateful for his investment and the opportunities afforded me, from which I have learned and benefited. I would next thank Drs.
    [Show full text]