"The UCSC Genome Browser"
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
ENCODE Genome-Wide Data on the UCSC Genome Browser Melissa Cline ENCODE Data Coordination Center (DCC) UC Santa Cruz
ENCODE Genome-Wide Data on the UCSC Genome Browser Melissa Cline ENCODE Data Coordination Center (DCC) UC Santa Cruz http://encodeproject.org Slides at http://genome-preview.ucsc.edu/ What is ENCODE? • International consortium project with the goal of cataloguing the functional regions of the human genome GTTTGCCATCTTTTG! CTGCTCTAGGGAATC" CAGCAGCTGTCACCA" TGTAAACAAGCCCAG" GCTAGACCAGTTACC" CTCATCATCTTAGCT" GATAGCCAGCCAGCC" ACCACAGGCATGAGT" • A gold mine of experimental data for independent researchers with available disk space ENCODE covers diverse regulatory processes ENCODE experiments are planned for integrative analysis Example of ENCODE data Genes Dnase HS Chromatin Marks Chromatin State Transcription Factor Binding Transcription RNA Binding Translation ENCODE tracks on the UCSC Genome Browser ENCODE tracks marked with the NHGRI helix There are currently 2061 ENCODE experiments at the ENCODE DCC How to find the data you want Finding ENCODE tracks the hard way A better way to find ENCODE tracks Finding ENCODE metadata descriptions Visualizing: Genome Browser tricks that every ENCODE user should know Turning ENCODE subtracks and views on and off Peaks Signal View on/off Subtrack on/off Right-click to the subtrack display menu Subtrack Drag and Drop Sessions: the easy way to save and share your work Downloading data with less pain 1. Via the Downloads button on the track details page 2. Via the File Selection tool Publishing: the ENCODE data release policy Every ENCODE subtrack has a “Restricted Until” date Key points of the ENCODE data release policy • Anyone is free to download and analyze data. • One cannot submit publications involving ENCODE data unless – the data has been at the ENCODE DCC for at least nine months, or – the data producers have published on the data, or – the data producers have granted permission to publish. -
Sequence Alignment/Map) Is a Text Format for Storing Sequence Alignment Data in a Series of Tab Delimited ASCII Columns
NGS FILE FORMATS SEQUENCE FILE FORMATS FASTA FORMAT FASTA Single sequence example: >HWI-ST398_0092:1:1:5372:2486#0/1 TTTTTCGTTCTTTTCATGTACCGCTTTTTGTTCGGTTAGATCGGAAGAGCGGTTCAGCAGGAATGCCGAGACCGAT ACGTAGCAGCAGCATCAGTACGACTACGACGACTAGCACATGCGACGATCGATGCTAGCTGACTATCGATG Multiple sequence example: >Sequence Name 1 TTTTTCGTTCTTTTCATGTACCGCTTTTTGTTCGGTTAGATCGGAAGAGCGGTTCAGCAGGAATGCCGAGACCGAT ACGTAGCAGCAGCATCAGTACGACTACGACGACTAGCACATGCGACGATCGATGCTAGCTGACTATCGATG >Sequence Name 2 ACGTAGACACGACTAGCATCAGCTACGCATCGATCAGCATCGACTAGCATCACACATCGATCAGCATCACGACTAGCAT AGCATCGACTACACTACGACTACGATCCACGTACGACTAGCATGCTAGCGCTAGCTAGCTAGCTAGTCGATCGATGAGT AGCTAGCTAGCTAGC >Sequence Name 3 ACTCAGCATGCATCAGCATCGACTACGACTACGACATCGACTAGCATCAGCAT SEQUENCE FILE FORMATS FASTQ FORMAT FASTQ Text based format for storing sequence data and corresponding quality scores for each base. To enable a one-one correspondence between the base sequence and the quality score the score is stored as a single one letter/number code using an offset of the standard ASCII code. Quality scores range from 0 to 40 and represent a log10 score for the probability of being wrong. E.g. score of 30 => 1:1000 chance of error SEQUENCE FILE FORMATS FASTQ FORMAT FASTQ Each fastq file contain multiple entries and each entry consists of 4 lines: 1. header line beginning with “@“ and sequence name 2. sequence line 3. header line beginning with “+” which can have the name but rarely does 4. quality score line SEQUENCE FILE FORMATS FASTQ FORMAT FASTQ @HWI-ST398_0092:6:73:5372:2486#0/1 TTTTTCGTTCTTTTCATGTACCGCTTTTTGTTCGGTTAGATCGGAAGAGCGGTTCAGCAGGAATGCCGAGACCGAT -
No Evidence for Recent Selection at FOXP2 Among Diverse Human Populations
Article No Evidence for Recent Selection at FOXP2 among Diverse Human Populations Graphical Abstract Authors Elizabeth Grace Atkinson, Amanda Jane Audesse, Julia Adela Palacios, Dean Michael Bobo, Ashley Elizabeth Webb, Sohini Ramachandran, Brenna Mariah Henn Correspondence [email protected] (E.G.A.), [email protected] (B.M.H.) In Brief An in-depth examination of diverse sets of human genomes argues against a recent selective evolutionary sweep of FOXP2, a gene that was believed to be critical for speech evolution in early hominins. Highlights d No support for positive selection at FOXP2 in large genomic datasets d Sample composition and genomic scale significantly affect selection scans d An intronic ROI within FOXP2 is expressed in human brain cells and cortical tissue d This ROI contains a large amount of constrained, human- specific polymorphisms Atkinson et al., 2018, Cell 174, 1424–1435 September 6, 2018 ª 2018 Elsevier Inc. https://doi.org/10.1016/j.cell.2018.06.048 Article No Evidence for Recent Selection at FOXP2 among Diverse Human Populations Elizabeth Grace Atkinson,1,8,9,10,* Amanda Jane Audesse,2,3 Julia Adela Palacios,4,5 Dean Michael Bobo,1 Ashley Elizabeth Webb,2,6 Sohini Ramachandran,4 and Brenna Mariah Henn1,7,* 1Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY, USA 2Department of Molecular Biology, Cell Biology and Biochemistry, Brown University, Providence, RI 02912, USA 3Neuroscience Graduate Program, Brown University, Providence, RI 02912, USA 4Department of Ecology and Evolutionary -
Databases/Resources on the Web
Jon K. Lærdahl, Structural Bioinforma�cs Databases/Resources on the web Jon K. Lærdahl [email protected] Jon K. Lærdahl, A lot of biological databases Structural Bioinforma�cs available on the web... MetaBase, the database of biological bioinforma�cs.ca – links directory databases (1801 entries) (620 databases) -‐ h�p://metadatabase.org -‐ h�p://bioinforma�cs.ca/links_directory Jon K. Lærdahl, Structural Bioinforma�cs btw, the bioinforma�cs.ca links directory is an excellent resource bioinforma�cs.ca – links directory h�p://bioinforma�cs.ca/links_directory Currently 1459 tools 620 databases 164 “resources” The problem is not to find a tool or database, but to know what is “gold” and what is “junk” Jon K. Lærdahl, Some important centres for Structural Bioinforma�cs bioinforma�cs Na�onal Center for Biotechnology Informa�on (NCBI) – part of the US Na�onal Library of Medicine (NLM), a branch of the Na�onal Ins�tutes of Health – located in Bethesda, Maryland European Bioinforma�cs Ins�tute (EMBL-‐EBI) – part of part of European Molecular Biology Laboratory (EMBL) – located in Hinxton, Cambridgeshire, UK Jon K. Lærdahl, NCBI databases Structural Bioinforma�cs Provided the GenBank DNA sequence database since 1992 Online Mendelian Inheritance in Man (OMIM) -‐ known diseases with a gene�c component and links to genes – started early 1960s as a book – online version, OMIM, since 1987 – on the WWW by NCBI in 1995 – currently >22,000 entries (14,400 genes) EST -‐ nucleo�de database subset that contains only Expressed Sequence Tag -
The UCSC Genome Browser Database: 2021 Update Jairo Navarro Gonzalez 1,*, Ann S
D1046–D1057 Nucleic Acids Research, 2021, Vol. 49, Database issue Published online 22 November 2020 doi: 10.1093/nar/gkaa1070 The UCSC Genome Browser database: 2021 update Jairo Navarro Gonzalez 1,*, Ann S. Zweig1, Matthew L. Speir1, Daniel Schmelter1, Kate R. Rosenbloom1, Brian J. Raney1, Conner C. Powell1, Luis R. Nassar1, Nathan D. Maulding1, Christopher M. Lee 1, Brian T. Lee 1,AngieS.Hinrichs1, Alastair C. Fyfe1, Jason D. Fernandes1, Mark Diekhans 1, Hiram Clawson1, Jonathan Casper1, Anna Benet-Pages` 1,2, Galt P. Barber1, David Haussler1,3, Robert M. Kuhn1, Maximilian Haeussler1 and W. James Kent1 Downloaded from https://academic.oup.com/nar/article/49/D1/D1046/5998393 by guest on 27 September 2021 1Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA, 2Medical Genetics Center (MGZ), Munich, Germany and 3Howard Hughes Medical Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA Received September 21, 2020; Revised October 19, 2020; Editorial Decision October 20, 2020; Accepted November 18, 2020 ABSTRACT pires to quickly incorporate and contextualize vast amounts of genomic information. For more than two decades, the UCSC Genome Apart from incorporating data from researchers and con- Browser database (https://genome.ucsc.edu)has sortia, the Browser also provides tools available for users to provided high-quality genomics data visualization view and compare their own data with ease. Custom tracks and genome annotations to the research community. allow users to quickly view a dataset, and track hubs allow As the field of genomics grows and more data be- users to extensively organize their data and share it privately come available, new modes of display are required using a URL. -
Alternate-Locus Aware Variant Calling in Whole Genome Sequencing Marten Jäger1,2, Max Schubach1, Tomasz Zemojtel1,Knutreinert3, Deanna M
Jäger et al. Genome Medicine (2016) 8:130 DOI 10.1186/s13073-016-0383-z RESEARCH Open Access Alternate-locus aware variant calling in whole genome sequencing Marten Jäger1,2, Max Schubach1, Tomasz Zemojtel1,KnutReinert3, Deanna M. Church4 and Peter N. Robinson1,2,3,5,6* Abstract Background: The last two human genome assemblies have extended the previous linear golden-path paradigm of the human genome to a graph-like model to better represent regions with a high degree of structural variability. The new model offers opportunities to improve the technical validity of variant calling in whole-genome sequencing (WGS). Methods: We developed an algorithm that analyzes the patterns of variant calls in the 178 structurally variable regions of the GRCh38 genome assembly, and infers whether a given sample is most likely to contain sequences from the primary assembly, an alternate locus, or their heterozygous combination at each of these 178 regions. We investigate 121 in-house WGS datasets that have been aligned to the GRCh37 and GRCh38 assemblies. Results: We show that stretches of sequences that are largely but not entirely identical between the primary assembly and an alternate locus can result in multiple variant calls against regions of the primary assembly. In WGS analysis, this results in characteristic and recognizable patterns of variant calls at positions that we term alignable scaffold-discrepant positions (ASDPs). In 121 in-house genomes, on average 51.8 ± 3.8 of the 178 regions were found to correspond best to an alternate locus rather than the primary assembly sequence, and filtering these genomes with our algorithm led to the identification of 7863 variant calls per genome that colocalized with ASDPs. -
The UCSC Genome Browser Database: Update 2011 Pauline A
D876–D882 Nucleic Acids Research, 2011, Vol. 39, Database issue Published online 18 October 2010 doi:10.1093/nar/gkq963 The UCSC Genome Browser database: update 2011 Pauline A. Fujita1,*, Brooke Rhead1, Ann S. Zweig1, Angie S. Hinrichs1, Donna Karolchik1, Melissa S. Cline1, Mary Goldman1, Galt P. Barber1, Hiram Clawson1, Antonio Coelho1, Mark Diekhans1, Timothy R. Dreszer1, Belinda M. Giardine2, Rachel A. Harte1, Jennifer Hillman-Jackson1, Fan Hsu1, Vanessa Kirkup1, Robert M. Kuhn1, Katrina Learned1, Chin H. Li1, Laurence R. Meyer1, Andy Pohl1,3, Brian J. Raney1, Kate R. Rosenbloom1, Kayla E. Smith1, David Haussler1,4 and W. James Kent1 1Center for Biomolecular Science and Engineering, School of Engineering, University of California Santa Cruz Downloaded from (UCSC), Santa Cruz, CA 95064, 2Center for Comparative Genomics and Bioinformatics, Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA 16802, USA, 3Centre for Genomic Regulation (CRG), Barcelona, Spain and 4Howard Hughes Medical Institute, UCSC, Santa Cruz, CA 95064, USA Received September 15, 2010; Accepted September 30, 2010 nar.oxfordjournals.org ABSTRACT differs among species, with recent assemblies of the human The University of California, Santa Cruz Genome genome being the most richly annotated. The Genome Browser (http://genome.ucsc.edu) offers online Browser contains mapping and sequencing annotation tracks describing assembly, gap and GC percent details access to a database of genomic sequence and at Health Sciences Center Library on February 4, 2011 annotation data for a wide variety of organisms. for all assemblies. Most organisms also have tracks con- taining alignments of RefSeq genes (3,4), mRNAs and The Browser also has many tools for visualizing, ESTs from GenBank (5) as well as gene and gene predic- comparing and analyzing both publicly available tion tracks such as Ensembl Genes (6). -
Galaxy Platform for NGS Data Analyses
Galaxy Platform For NGS Data Analyses Weihong Yan [email protected] Collaboratory Web Site http://qcb.ucla.edu/collaboratory Collaboratory Workshops Workshop Outline ü Day 1 § UCLA galaxy and user account § Galaxy web interface and management § Tools for NGS analyses and their application § Data formats § Build/share workflow and history § Q and A ü Day 2 § Galaxy Tools for RNA-seq analysis § Galaxy Tools for ChIP-seq analysis § Galaxy Tools for annotation. § Q and A *** Published datasets/results will be used in the tutorial UCLA Galaxy http://galaxy.hoffman2.idre.ucla.edu ü Hardware – Headnode (1) 96Gb memory, 12 core – Computing nodes (8) 48Gb memory, 12 core – Storage 100 Tb disk space ü Galaxy Resource Management - Hoffman2 grid engine Default: 1 core/job bowtie, bwa, tophat, cuffdiff, cufflinks, gatk programs: 4 core/job UCLA Galaxy http://galaxy.hoffman2.idre.ucla.edu ü galaxy login account: login: your email associated with ucla ü Disk quota: 1 Tb/user Galaxy Account Management Installed tools Launch analysis and view result History of execu7on and results Raw Reads *_qseq.txt, *.fastq Upload to Galaxy File transfer protocol (ftp) deMultiplex Barcode splitter, deMultiplex workflow fastqc, compute quality statistics, Quality Assessment draw quality score boxplot, draw nuclotides distribution Process Reads Trim sequences, sickle, scythe Alignment to bwa, bowtie, bowtie2, tophat Reference Format Conversion Text manipulation toolkit, BEDTools, SAM Results (sam/bam) Tools, java genomics toolkit, picard toolkit Downstream Analyses BS-Seeker2, cufflinks, cuffdiff, macs, macs2, GATK, CEAS Visualization Genome browser, IGV Repositories of Galaxy Tools https://toolshed.g2.bx.psu.edu ü History panel contains all datasets that are uploaded and results derived from certain analyses ü A history can be organized, annotated, and managed as a project ü History is sharable. -
Genetic Analyses in a Bonobo (Pan Paniscus) with Arrhythmogenic Right Ventricular Cardiomyopathy Received: 19 September 2017 Patrícia B
www.nature.com/scientificreports OPEN Genetic analyses in a bonobo (Pan paniscus) with arrhythmogenic right ventricular cardiomyopathy Received: 19 September 2017 Patrícia B. S. Celestino-Soper1, Ty C. Lynnes1, Lili Zhang1, Karen Ouyang1, Samuel Wann2, Accepted: 21 February 2018 Victoria L. Clyde2 & Matteo Vatta1,3 Published: xx xx xxxx Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a disorder that may lead to sudden death and can afect humans and other primates. In 2012, the alpha male bonobo of the Milwaukee County Zoo died suddenly and histologic evaluation found features of ARVC. This study sought to discover a possible genetic cause for ARVC in this individual. We sequenced our subject’s DNA to search for deleterious variants in genes involved in cardiovascular disorders. Variants found were annotated according to the human genome, following currently available classifcation used for human diseases. Sequencing from the DNA of an unrelated unafected bonobo was also used for prediction of pathogenicity. Twenty-four variants of uncertain clinical signifcance (VUSs) but no pathogenic variants were found in the proband studied. Further familial, functional, and bonobo population studies are needed to determine if any of the VUSs or a combination of the VUSs found may be associated with the clinical fndings. Future genotype-phenotype establishment will be benefcial for the appropriate care of the captive zoo bonobo population world-wide as well as conservation of the bobono species in its native habitat. Cardiovascular disease is a leading cause of death for both human and non-human primates, who share similar genomes and many environmental and lifestyle characteristics1. However, while humans most ofen die due to coronary artery disease (CAD), non-human primates living in captivity are more ofen afected by hypertensive cardiomyopathy1. -
Cgb − a Unix Shell Program to Create Custom Instances of the Ucsc Genome Browser
CGB − A UNIX SHELL PROGRAM TO CREATE CUSTOM INSTANCES OF THE UCSC GENOME BROWSER Vincenzo Forgetta1* & Ken Dewar1 1Department of Human Genetics, McGill University, Montreal, Quebec, Canada. *Corresponding author: Lady Davis Institute for Medical Research, Jewish General Hospital, 3755 Côte Ste-Catherine Road, H-465 (Pavilion H), Montreal, Québec, H3T 1E2, Tel: +1 (514) 340-8222 ext.3982, Fax: +1 (514) 340-7564, Email: [email protected] Keywords: UCSC Genome Browser, genomics, custom mirror procedure. Abstract The UCSC Genome Browser is a popular tool for the exploration and analysis of reference genomes. Mirrors of the UCSC Genome Browser and its contents exist at multiple geographic locations, and this mirror procedure has been modified to support genome sequences not maintained by UCSC and generated by individual researchers. While straightforward, this procedure is lengthy and tedious and would benefit from automation, especially when processing many genome sequences. We present a Unix shell program that facilitates the creation of custom instances of the UCSC Genome Browser for genome sequences not being maintained by UCSC. It automates many steps of the browser creation process, provides password protection for each browser instance, and automates the creation of basic annotation tracks. As an example we generate a custom UCSC Genome Browser for a bacterial genome obtained from a massively parallel sequencing platform. Introduction In the past, large institutions sequenced de novo the genome of organisms such as human (Lander et al., 2001), mouse (Waterston et al., 2002) and fly (Adams et al., 2000), and bioinformatics tools were created to provide the scientific research community with access to analyzing these reference genome resources. -
UCSC Genome Browser Outline
UCSC Genome Browser Introduction to ab initio and evidence-based gene finding Wilson Leung 06/2006 Outline Introduction to annotation ab initio gene finding Basics of the UCSC Browser Evidence-based gene finding using the UCSC Browser 1 AAACAACAATCATAAATAGAGGAAGTTTTCGGAATATACGATAAGTGAAAT ATCGTTCTTAAAAAAGAGCAAGAACAGTTTAACCATTGAAAACAAGATTATT CCAATAGCCGTAAGAGTTCATTTAATGACAATGACGATGGCGGCAAAGTCGA TGAAGGACTAGTCGGAACTGGAAATAGGAATGCGCCAAAAGCTAGTGCAGC TAAACATCAATTGAAACAAGTTTGTACATCGATGCGCGGAGGCGCTTTTCTC TCAGGATGGCTGGGGATGCCAGCACGTTAATCAGGATACCAATTGAGGAGGT GCCCRAGCTCACCEAGAGCCGGCCAATAAGGACCCATCGGGGGGGCCGCTTA TGTGGAAGCCAAACATTAAACCATAGGCAACCGATTTGTGGGAATCGAATTT AAGAAACGGCGGTCAGCCACCCGCTCAACAAGTGCCAAAGCCATCTTGGGG GCATACGCCTTCATCAAATTTGGGCGGAACTTGGGGCGAGGACGATGATGGC GCCGATAGCACCAGCGTTTGGACGGGTCAGTCATTCCAHATAEGCACAACGT CTGGTGTTGCAGTCGGTGCCATAGCGCCTGGCEGTTGGCGCCGCTGCTGGTCC CCAATGGGGACAGGCTGTTGCTGTTGGTGTTGGAGTCGGAGTTGCCTTAAAC TCGACTGGAAATAACAATGCGCCGGCAACAGGAGCCCTGCCTGCCGTGGCTC CCGCGGCGCCGCGGAGAATAGTNGTNATAANAGNAGAAOATOATAATACTACAAAATATTTTGTGGCGIGCAIGCAACOAAOAATAGTNGGNCGACCAGCT GTCCGAAATGTGGGGACATCATCCTCAGATTGCTCACAATCAGATTCAGAAT ATTAACAAAATGGTCGGCCCCGTTGTTAGTGCAACAGGGTCAAATATCGCAA GCTCAAATATTGGCCCAAGCGGTGTTGGTTCCGTATCCGGTAATGTCGGGGC ACAATGGGGAGCCACACAGGCCGCGTTGGGGCCCCAAGGTATTTCCAAGCA AATCACTGGATGGGAGGAACCACAATCAGATTCAGAATATTAACAAAATGG TCGGCCCCGTTGTTATGGATAAAAAATTTGTGTCTTCGTACGGAGATTATGTT GTTAATCAATTTTATTAAGATATTTAAATAAATATGTGTACCTTTCACGAGAA ATTTGCTTACCTTTTCGACACACACACTTATACAGACAGGTAATAATTACCTT TTGAGCAATTCGATTTTCATAAAATATACCTAAATCGCATCGTCTATGAATCT TTGTAATACTTTCGAATTTAATTATTAGTCTACATTAATATTGATACCGTTTTT -
BMC Bioinformatics Biomed Central
BMC Bioinformatics BioMed Central Database Open Access Atlas – a data warehouse for integrative bioinformatics Sohrab P Shah, Yong Huang, Tao Xu, Macaire MS Yuen, John Ling and BF Francis Ouellette* Address: UBC Bioinformatics Centre, University of British Columbia, Vancouver, BC, Canada Email: Sohrab P Shah - [email protected]; Yong Huang - [email protected]; Tao Xu - [email protected]; Macaire MS Yuen - [email protected]; John Ling - [email protected]; BF Francis Ouellette* - [email protected] * Corresponding author Published: 21 February 2005 Received: 04 September 2004 Accepted: 21 February 2005 BMC Bioinformatics 2005, 6:34 doi:10.1186/1471-2105-6-34 This article is available from: http://www.biomedcentral.com/1471-2105/6/34 © 2005 Shah et al; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract Background: We present a biological data warehouse called Atlas that locally stores and integrates biological sequences, molecular interactions, homology information, functional annotations of genes, and biological ontologies. The goal of the system is to provide data, as well as a software infrastructure for bioinformatics research and development. Description: The Atlas system is based on relational data models that we developed for each of the source data types. Data stored within these relational models are managed through Structured Query Language (SQL) calls that are implemented in a set of Application Programming Interfaces (APIs).