A Database of Human Immune Receptor Alleles Recovered from Population Sequencing Data

Total Page:16

File Type:pdf, Size:1020Kb

A Database of Human Immune Receptor Alleles Recovered from Population Sequencing Data A Database of Human Immune Receptor Alleles Recovered from Population Sequencing Data This information is current as Yaxuan Yu, Rhodri Ceredig and Cathal Seoighe of September 26, 2021. J Immunol published online 23 January 2017 http://www.jimmunol.org/content/early/2017/01/20/jimmun ol.1601710 Downloaded from Supplementary http://www.jimmunol.org/content/suppl/2017/01/20/jimmunol.160171 Material 0.DCSupplemental Why The JI? Submit online. http://www.jimmunol.org/ • Rapid Reviews! 30 days* from submission to initial decision • No Triage! Every submission reviewed by practicing scientists • Fast Publication! 4 weeks from acceptance to publication *average by guest on September 26, 2021 Subscription Information about subscribing to The Journal of Immunology is online at: http://jimmunol.org/subscription Permissions Submit copyright permission requests at: http://www.aai.org/About/Publications/JI/copyright.html Email Alerts Receive free email-alerts when new articles cite this article. Sign up at: http://jimmunol.org/alerts Errata An erratum has been published regarding this article. Please see next page or: /content/198/9/3758.full.pdf The Journal of Immunology is published twice each month by The American Association of Immunologists, Inc., 1451 Rockville Pike, Suite 650, Rockville, MD 20852 Copyright © 2017 by The American Association of Immunologists, Inc. All rights reserved. Print ISSN: 0022-1767 Online ISSN: 1550-6606. Published January 23, 2017, doi:10.4049/jimmunol.1601710 The Journal of Immunology A Database of Human Immune Receptor Alleles Recovered from Population Sequencing Data Yaxuan Yu,* Rhodri Ceredig,† and Cathal Seoighe* High-throughput sequencing data from TCRs and Igs can provide valuable insights into the adaptive immune response, but bioinformatics pipelines for analysis of these data are constrained by the availability of accurate and comprehensive repositories of TCR and Ig alleles. We have created an analytical pipeline to recover immune receptor alleles from genome sequencing data. Applying this pipeline to data from the 1000 Genomes Project we have created Lym1K, a collection of immune receptor alleles that combines known, well-supported alleles with novel alleles found in the 1000 Genomes Project data. We show that Lym1K leads to a significant improvement in the alignment of short read sequences from immune receptors and that the addition of novel alleles discovered from genome sequence data are likely to be particularly significant for comprehensive analysis of populations that are not currently well represented in existing repositories of immune alleles. The Journal of Immunology, 2017, 198: 000–000. Downloaded from he adaptive immune system provides protection from input sequences with the reference V, D (for Ig H chain and TCR disease by initiating specific immune responses to specific b-chain), and J genes, determine the CDR3 region, and estimate T foreign Ags following their recognition by the clonally the repertoire diversity. Two main factors that directly influence distributed surface-bound receptors of T and B lymphocytes, namely the quality of these analyses are the alignment algorithm and the TCRs for T cells and Igs for B cells. In its secreted form, the B cell’s completeness and accuracy of the reference database. In the past 3 http://www.jimmunol.org/ Ig molecule is called an Ab. The variability in molecular shapes of y, several software packages have been developed with different foreign Ags represents a significant challenge for the immune alignment algorithms for analyzing next-generation sequencing system, requiring development of an extensive receptor repertoire sequence data of TCR and/or Ig (11–19). These studies all focused for Ag recognition. In response to this challenge, TCRs and Igs on improving the accuracy, effectiveness, and completeness of achieve vast diversity by a unique mechanism called VDJ recom- their alignment algorithms. Most of the software packages use bination, which is a largely stochastic process of gene rearrange- the reference V(D)J sequences provided by the international ment encoding the variable parts of TCRs and Igs, respectively. ImMunoGeneTics information system (IMGT) (20), which is cur- Additionally, genes encoding Igs undergo somatic hypermutation rently the most comprehensive collection of TCR and Ig germline during T cell–dependent B cell responses, resulting in further ex- gene sequences available. by guest on September 26, 2021 pansion in Ab diversity and subsequent selection of B cell clones However, previous studies have reported numerous novel V generating Igs with higher affinity for Ag. Quantitative analysis of alleles that are not found in the IMGT database (21–25), suggesting the diversity of TCR and Ig repertoires can shed light on receptor– that this database is incomplete. Additionally, a study of 226 Ag interactions and reveal insights into the state of the immune IGHV alleles in IMGT reported a high rate of false positives (of system under diverse conditions, such as aging (1–3), chronic viral 226 alleles, 104 were classified as very likely to contain se- infections (4–6), and autoimmune diseases (7–10). quencing errors) (26). This is because most of the IMGT TRBV Next-generation sequencing technology revolutionized the and IGHV alleles were included from early studies, conducted analysis of TCR/Ig diversity by providing vast amounts of data at between 1984 and 1996 (27, 28). After that, there were only two much higher resolution, leading also to a need for effective and major updates on IGHV alleles in 2002 and 2009 introducing in robust bioinformatics pipelines. Generally, these pipelines align the total 36 alleles, and no updates on TRBV alleles have been sub- mitted in the past 20 y. Furthermore, many alleles were reported in only one study, and multiple allelic variants of the same gene were *School of Mathematics, Statistics and Applied Mathematics, National University of † derived from the same individual. For instance, an individual can Ireland Galway, Galway H91 HX31, Ireland; and Discipline of Psychology, National University of Ireland Galway, Galway H91 HX31, Ireland have at most two alleles of any gene; however, seven IGHV3-15 ORCIDs: 0000-0002-4393-727X (R.C.); 0000-0003-0496-4029 (C.S.). alleles were derived from one individual. This was pointed out by Received for publication October 5, 2016. Accepted for publication January 3, 2017. Wang et al. (26), but the issue remains unresolved in the current version of IMGT. Additionally, many TRBV and TRAV alleles in This work was supported by the European Union Seventh Framework Programme (FP7/2007–2013) under Grant 317040 (QuanTI). Funding for open access was by the IMGT are annotated as having come from rearranged sequences, European Union Seventh Framework Programme. which may contain somatic mutations and partially truncated 39 Address correspondence and reprint requests to Prof. Cathal Seoighe, National ends resulting from VDJ recombination. University of Ireland Galway, School of Mathematics, Statistics and Applied Since the development of IMGT, a large and ever-increasing Mathematics, University Road, Galway H91 HX31, Ireland. E-mail address: [email protected] number of full human genomes have been sequenced. The chro- The online version of this article contains supplemental material. mosomal regions containing the immune receptor alleles have both Abbreviations used in this article: G1K, 1000 Genomes Project; IMGT, international a high level of intrachromosomal duplication and a high level of ImMunoGeneTics information system; SNP, single nucleotide polymorphism; SRA, diversity and thus present challenges for genome assembly. National Center for Biotechnology Information Sequence Read Archive; VCF, var- However, improvement of the reference genome over time has the iant cell format. potential to enable recovery of a comprehensive set of immune gene Copyright Ó 2017 by The American Association of Immunologists, Inc. 0022-1767/17/$30.00 alleles found in human populations from population resequencing www.jimmunol.org/cgi/doi/10.4049/jimmunol.1601710 2 TCR/Ig ALLELES RECOVERED FROM POPULATION SEQUENCING DATA studies. We have developed a bioinformatics pipeline that can be we randomly chose a V, D (only for the H chain of Igs and the b-chain of used to generate a TCR/Ig reference database using data from TCRs), and J gene assuming uniform gene usage. We then decided the population resequencing studies. We applied this pipeline to data particular allele of the chosen V gene based on the frequency of its oc- currence in the chosen population. Second, we “mutated” the chosen V and from the 1000 Genomes Project (G1K) (29), which contains the J segments by introducing 0–10 mismatches on the V region and 0–5 most comprehensive collection of global human variation currently mismatches on the J region. To simulate junction diversity, we inserted available in the public domain, to create the Lym1K database. We zero to six randomly generated nucleotides into the V-D and D-J junction a aligned short read sequence data from human immune receptors to (V-J junction only for the Ig L chain and TCR -chain). We applied our simulation pipeline for TCRB, TCRA, IGH, and IGL Lym1K and obtained significantly improved alignments compared genes separately. In all, we produced 27 datasets, with separate datasets with what could be achieved using the IMGT database as reference. corresponding to each of the 26 populations included in G1K and 1 dataset We also report evidence that the diversity of immune receptor al- derived from the pool of all populations. Each dataset consisted of 20 leles observed in non-African populations is particularly under- simulated samples, and each sample contained 200,000 simulated represented in the existing database. sequences. G1K data Materials and Methods Immune receptor sequencing dataset We retrieved variant cell format (VCF) files corresponding to the TCRB, TCRA, IGH, and IGL regions from phase 3 of G1K. We used the R We obtained two datasets from the National Center for Biotechnology In- package SNPRelate (30) to perform principal component analysis sepa- formation Sequence Read Archive (SRA).
Recommended publications
  • Whole Genome Sequencing Data of Multiple Individuals of Pakistani
    www.nature.com/scientificdata oPeN Whole genome sequencing data DAtA DeScriptor of multiple individuals of Pakistani descent Shahid Y. Khan1, Muhammad Ali1, Mei-Chong W. Lee2, Zhiwei Ma3, Pooja Biswas4, Asma A. Khan5, Muhammad Asif Naeem5, Saima Riazuddin 6, Sheikh Riazuddin5,7,8, Radha Ayyagari4, J. Fielding Hejtmancik 3 & S. Amer Riazuddin1 ✉ Here we report whole genome sequencing of four individuals (H3, H4, H5, and H6) from a family of Pakistani descent. Whole genome sequencing yielded 1084.92, 894.73, 1068.62, and 1005.77 million mapped reads corresponding to 162.73, 134.21, 160.29, and 150.86 Gb sequence data and 52.49x, 43.29x, 51.70x, and 48.66x average coverage for H3, H4, H5, and H6, respectively. We identifed 3,529,659, 3,478,495, 3,407,895, and 3,426,862 variants in the genomes of H3, H4, H5, and H6, respectively, including 1,668,024 variants common in the four genomes. Further, we identifed 42,422, 39,824, 28,599, and 35,206 novel variants in the genomes of H3, H4, H5, and H6, respectively. A major fraction of the variants identifed in the four genomes reside within the intergenic regions of the genome. Single nucleotide polymorphism (SNP) genotype based comparative analysis with ethnic populations of 1000 Genomes database linked the ancestry of all four genomes with the South Asian populations, which was further supported by mitochondria based haplogroup analysis. In conclusion, we report whole genome sequencing of four individuals of Pakistani descent. Background & Summary The completion of Human Genome Project ignited several large scale efforts to characterize variations in the human genome, which led to a comprehensive catalog of the common variants including single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels), across the entire human genome1,2.
    [Show full text]
  • Genomic Data Standards Resources and Initiatives Cited in the Supplemental Information to the Genomic Data Sharing Policy
    Genomic Data Standards Resources and Initiatives Cited in the Supplemental Information to the Genomic Data Sharing Policy IMPORTANT NOTE: The National Institutes of Health makes no endorsement of the non-NIH-funded genomic data standards resources and/or initiatives included in this document. Categories Resources and Initiatives Common Data Element Resource Portal: http://www.nlm.nih.gov/cde/ NIH encourages the use of common data elements (CDEs) in clinical research, patient registries, and other human subject research in order to improve data quality and opportunities for comparison and combination of data from multiple studies and with electronic health records. This portal provides access to NIH-supported CDE initiatives and other tools and resources that can assist investigators developing protocols for data collection. Clinical Genome Resource (ClinGen): http://www.nih.gov/news/health/sep2013/nhgri-25.htm In 2013, the NIH National Human Genome Research Institute (NHGRI) and the Eunice Kennedy Shriver National Institute NIH-Funded of Child Health and Human Development (NICHD) awarded three grants totalling over $25 million to support a consortium of research groups to design and implement a framework for evaluating variants that are relevant to patient care (e.g., play a role in disease). International Collaboration for Clinical Genomics (ICCG): http://www.iccg.org/about-the-iccg/clingen/ ICCG was awarded as a grant under ClinGen. ICCG was charged with developing standard formats for gathering and depositing data in ClinVar (http://www.ncbi.nlm.nih.gov/clinvar/). It will work with a variety of different stakeholder groups, including clinical laboratories and existing locus-specific databases, to obtain robust data sets on genomic variants and disease associations.
    [Show full text]
  • Comparative Analysis of Pacbio and Oxford Nanopore Sequencing Technologies for Transcriptomic Landscape Identification of Penaeus Monodon
    life Brief Report Comparative Analysis of PacBio and Oxford Nanopore Sequencing Technologies for Transcriptomic Landscape Identification of Penaeus monodon Zulema Udaondo 1 , Kanchana Sittikankaew 2, Tanaporn Uengwetwanit 2 , Thidathip Wongsurawat 1,3 , Chutima Sonthirod 4, Piroon Jenjaroenpun 1,3 , Wirulda Pootakham 4, Nitsara Karoonuthaisiri 2 and Intawat Nookaew 1,* 1 Department of Biomedical Informatics, University of Arkansas for Medical Sciences, Little Rock, AR 72205, USA; [email protected] (Z.U.); [email protected] (T.W.); [email protected] (P.J.) 2 National Center for Genetic Engineering and Biotechnology (BIOTEC), National Science and Technology Development Agency, Pathum Thani 12120, Thailand; [email protected] (K.S.); [email protected] (T.U.); [email protected] (N.K.) 3 Division of Bioinformatics and Data Management for Research, Department of Research and Development, Faculty of Medicine, Siriraj Hospital, Mahidol University, Bangkok 10700, Thailand 4 National Omics Center (NOC), National Science and Technology Development Agency (NSTDA), Pathum Thani 12120, Thailand; [email protected] (C.S.); [email protected] (W.P.) * Correspondence: [email protected]; Tel.: +1-501-686-6025; Fax: +1-501-603-1766 Citation: Udaondo, Z.; Sittikankaew, K.; Uengwetwanit, T.; Wongsurawat, Abstract: With the advantages that long-read sequencing platforms such as Pacific Biosciences T.; Sonthirod, C.; Jenjaroenpun, P.; (Menlo Park, CA, USA) (PacBio) and Oxford Nanopore Technologies (Oxford, UK) (ONT) can offer, Pootakham, W.; Karoonuthaisiri, N.; various research fields such as genomics and transcriptomics can exploit their benefits. Selecting an Nookaew, I. Comparative Analysis of appropriate sequencing platform is undoubtedly crucial for the success of the research outcome, PacBio and Oxford Nanopore thus there is a need to compare these long-read sequencing platforms and evaluate them for specific Sequencing Technologies for research questions.
    [Show full text]
  • Globalfungi, a Global Database of Fungal Occurrences from High
    www.nature.com/scientificdata OPEN GlobalFungi, a global database DATA DEscrIPTor of fungal occurrences from high-throughput-sequencing metabarcoding studies Tomáš Větrovský1,6, Daniel Morais1,6, Petr Kohout1,6, Clémentine Lepinay1,6, Camelia Algora1, Sandra Awokunle Hollá1, Barbara Doreen Bahnmann1, Květa Bílohnědá1, Vendula Brabcová1, Federica D’Alò2, Zander Rainier Human1, Mayuko Jomura 3, Miroslav Kolařík1, Jana Kvasničková1, Salvador Lladó1, Rubén López-Mondéjar1, Tijana Martinović1, Tereza Mašínová1, Lenka Meszárošová1, Lenka Michalčíková1, Tereza Michalová1, Sunil Mundra4,5, Diana Navrátilová1, Iñaki Odriozola 1, Sarah Piché-Choquette 1, Martina Štursová1, Karel Švec1, Vojtěch Tláskal 1, Michaela Urbanová1, Lukáš Vlk1, Jana Voříšková1, Lucia Žifčáková1 & Petr Baldrian 1 ✉ Fungi are key players in vital ecosystem services, spanning carbon cycling, decomposition, symbiotic associations with cultivated and wild plants and pathogenicity. The high importance of fungi in ecosystem processes contrasts with the incompleteness of our understanding of the patterns of fungal biogeography and the environmental factors that drive those patterns. To reduce this gap of knowledge, we collected and validated data published on the composition of soil fungal communities in terrestrial environments including soil and plant-associated habitats and made them publicly accessible through a user interface at https://globalfungi.com. The GlobalFungi database contains over 600 million observations of fungal sequences across > 17 000 samples with geographical locations and additional metadata contained in 178 original studies with millions of unique nucleotide sequences (sequence variants) of the fungal internal transcribed spacers (ITS) 1 and 2 representing fungal species and genera. The study represents the most comprehensive atlas of global fungal distribution, and it is framed in such a way that third-party data addition is possible.
    [Show full text]
  • Closure of the NCBI SRA and Implications for the Long-Term Future of Genomics Data Storage GB Editorial Team*
    GB Editorial Team Genome Biology 2011, 12:402 http://genomebiology.com/2011/12/3/402 EDITORIAL Closure of the NCBI SRA and implications for the long-term future of genomics data storage GB Editorial Team* e National Center for Biotechnology Information accept submissions doesn’t meant that the SRA is closing, (NCBI) in the US recently announced that, as a result of merely changing and the European Nucleotide Archive budgetary constraints, it would no longer be accepting (ENA) at EMBL-EBI will remain. e NCBI’s decision submissions to its Sequence Read Archive (SRA) and that was based on budgetary constraints. It should be noted over the course of the next year or so it would slowly that most people don’t realize that storage space is only a phase out support for this database (http://www.ncbi. minor fraction of the budget of the database; the bulk of nlm.nih.gov/sra). ere seems to be a certain amount of the cost is associated with the staff who maintain the confusion in the community about what effect this database, process the submissions, develop the software decision will have. At Genome Biology we feel that the and so on. free availability of data is an important concept for science, SS: From the outside, it appears that the SRA is closing so we asked the views of various interested people on what because of NIH budgetary considerations. One problem the short-term implications of this announcement will is that the amount of sequence being generated is be, and also how they envisaged the future of data storage growing at an extraordinary rate, probably faster than in the long term.
    [Show full text]
  • L3: Short Read Alignment to a Reference Genome
    L3: Short Read Alignment to a Reference Genome Shamith Samarajiwa CRUK Autumn School in Bioinformatics Cambridge, September 2017 Where to get help! http://seqanswers.com http://www.biostars.org http://www.bioconductor.org/help/mailing-list Read the posting guide before sending email! Overview ● Understand the difference between reference genome builds ● Introduction to Illumina sequencing ● Short read aligners ○ BWA ○ Bowtie ○ STAR ○ Other aligners ● Coverage and Depth ● Mappability ● Use of decoy and sponge databases ● Alignment Quality, SAMStat, Qualimap ● Samtools and Picard tools, ● Visualization of alignment data ● A very brief look at long reads, graph genome aligners and de novo genome assembly Reference Genomes ● A haploid representation of a species genome. ● The human genome is a haploid mosaic derived from 13 volunteer donors from Buffalo, NY. USA. ● In regions where there is known large scale population variation, sets of alternate loci (178 in GRCh38) are assembled alongside the reference locus. ● The current build has around 500 gaps, whereas the first version had ~150,000 gaps. Genome Reference Consortium: https://www.ncbi.nlm.nih.gov/grc £600 30X Sequencers http://omicsmap.com Illumina sequencers Illumina Genome Analyzer Illumina sequencing technology ● Illumina sequencing is based on the Solexa technology developed by Shankar Balasubramanian and David Klenerman (1998) at the University of Cambridge. ● Multiple steps in “Sequencing by synthesis” (explained in next slide) ○ Library Preparation ○ Bridge amplification and Cluster generation ○ Sequencing using reversible terminators ○ Image acquisition and Fastq generation ○ Alignment and data analysis Illumina Flowcell Sequencing By Synthesis technology Illumina Sequencing Bridge Amplification Sequencing Cluster growth Incorporation of fluorescence, reversibly terminated tagged nt Multiplexing • Multiplexing gives the ability to sequence multiple samples at the same time.
    [Show full text]
  • An Introduction to Next-Generation Sequencing Technology
    An introduction to Next-Generation Sequencing Technology www.illumina.com/technology/next-generation-sequencing.html For Research Use Only. Not for use in diagnostic procedures. Table of Contents Table of Contents 2 I. Welcome to Next-Generation Sequencing 3 a. The Evolution of Genomic Science 3 b. The Basics of NGS Chemistry 4 c. Advances in Sequencing Technology 5 Paired-End Sequencing 5 Tunable Coverage and Unlimited Dynamic Range 6 Advances in Library Preparation 6 Multiplexing 7 Flexible, Scalable Instrumentation 7 II. NGS Methods 8 a. Genomics 8 Whole-Genome Sequencing 8 Exome Sequencing 8 De novo Sequencing 9 Targeted Sequencing 9 b. Transcriptomics 11 Total RNA and mRNA Sequencing 11 Targeted RNA Sequencing 11 Small RNA and Noncoding RNASequencing 11 c. Epigenomics 12 Methylation Sequencing 12 ChIP Sequencing 12 Ribosome Profiling 12 III. Illumina DNA-to-Data NGS Solutions 13 a. The Illumina NGS Workflow 13 b. Integrated Data Analysis 13 IV. Glossary 14 V. References 15 For Research Use Only. Not for use in diagnostic procedures. I. Welcome to Next-Generation Sequencing a. The Evolution of Genomic Science DNA sequencing has come a long way since the days of two-dimensional chromatography in the 1970s. With the advent of the Sanger chain termination method1 in 1977, scientists gained the ability to sequence DNA in a reliable, reproducible manner. A decade later, Applied Biosystems introduced the first automated, capillary electrophoresis (CE)-based sequencing instruments,the AB370 in 1987 and the AB3730xl in 1998,instruments that became the primary workhorses for the NIH-led and Celera-led Human Genome Projects.2 While these “first-generation” instruments were considered high throughput for their time, the Genome Analyzer emerged in 2005 and took sequencing runs from 84 kilobase (kb) per run to 1 gigabase (Gb) per run.3 The short read, massively parallel sequencing technique was a fundamentally different approach that revolutionized sequencing capabilities and launched the “next generation” in genomic science.
    [Show full text]
  • Sampling Complementary RNA Reads from the Sequence Read Archive Mario Stanke1,2* , Willy Bruhn1, Felix Becker1,2 and Katharina J
    Stanke et al. BMC Bioinformatics (2019) 20:558 https://doi.org/10.1186/s12859-019-3182-x SOFTWARE Open Access VARUS: sampling complementary RNA reads from the sequence read archive Mario Stanke1,2* , Willy Bruhn1, Felix Becker1,2 and Katharina J. Hoff1,2 Abstract Background: Vast amounts of next generation sequencing RNA data has been deposited in archives, accompanying very diverse original studies. The data is readily available also for other purposes such as genome annotation or transcriptome assembly. However, selecting a subset of available experiments, sequencing runs and reads for this purpose is a nontrivial task and complicated by the inhomogeneity of the data. Results: This article presents the software VARUS that selects, downloads and aligns reads from NCBI’s Sequence Read Archive, given only the species’ binomial name and genome. VARUS automatically chooses runs from among all archived runs to randomly select subsets of reads. The objective of its online algorithm is to cover a large number of transcripts adequately when network bandwidth and computing resources are limited. For most tested species VARUS achieved both a higher sensitivity and specificity with a lower number of downloaded reads than when runs were manually selected. At the example of twelve eukaryotic genomes, we show that RNA-Seq that was sampled with VARUS is well-suited for fully-automatic genome annotation with BRAKER. Conclusions: With VARUS, genome annotation can be automatized to the extent that not even the selection and quality control of RNA-Seq has to be done manually. This introduces the possibility to have fully automatized genome annotation loops over potentially many species without incurring a loss of accuracy over a manually supervised annotation process.
    [Show full text]
  • SRA Handbook
    SRA Handbook Last Updated: 2016 Jan 14 National Center for Biotechnology Information (US) Bethesda (MD) National Center for Biotechnology Information (US), Bethesda (MD) NLM Citation: SRA Handbook [Internet]. Bethesda (MD): National Center for Biotechnology Information (US); 2010-. iii This documentation provides an overview and help manual for the Sequence Read Archive (SRA) at the National Center for Biotechnology Information. iv SRA Handbook Table of Contents Introduction to SRA................................................................................................................................. 1 Concepts ............................................................................................................................................................. 3 Overview......................................................................................................................................................... 3 Concepts.......................................................................................................................................................... 3 Overview........................................................................................................................................................... 9 An Introduction to the Sequence Read Archive...................................................................................... 9 Terms of Usage .............................................................................................................................................
    [Show full text]
  • Training Modules Documentation Release 1
    training_modules Documentation Release 1 Marc B Rossello Sep 24, 2021 ENA Data Submission 1 General Guide on ENA Data Submission3 2 How to Register a Study 29 3 How to Register Samples 35 4 Preparing Files for Submission 47 5 How to Submit Raw Reads 83 6 How to Submit Assemblies 103 7 How to Submit Targeted Sequences 145 8 How to Submit Other Analyses 161 9 General Guide on ENA Data Retrieval 183 10 Updating Metadata Objects 223 11 Updating Assemblies 231 12 Updating Annotated Sequences 233 13 Data Release Policies 235 14 Common Run Submission Errors 237 15 Common Sample Submission Errors 243 16 Tips for Sample Taxonomy 245 17 Requesting New Taxon IDs 249 18 Metagenome Submission Queries 255 19 Locus Tag Prefixes 261 20 Archive Generated FASTQ Files 263 i 21 Third Party Tools 269 22 Introductory Webinar 271 ii training_modules Documentation, Release 1 Welcome to the guidelines for submission and retrieval for the European Nucleotide Archive. Please use the links to find instructions specific to your needs. If you’re completely new to ENA, you can see an introductory webinar at the bottom of the page. ENA Data Submission 1 training_modules Documentation, Release 1 2 ENA Data Submission CHAPTER 1 General Guide on ENA Data Submission Welcome to the general guide for the European Nucleotide Archive submission. Please take a moment to view this introduction and consider the options available to you before you begin your submission. ENA allows submissions via three routes, each of which is appropriate for a different set of submission types.
    [Show full text]
  • Assessing Metagenomic Signals Recovered from Lyuba, a 42,000-Year-Old Permafrost-Preserved Woolly Mammoth Calf
    G C A T T A C G G C A T genes Article Assessing Metagenomic Signals Recovered from Lyuba, a 42,000-Year-Old Permafrost-Preserved Woolly Mammoth Calf Giada Ferrari 1,2,* ID , Heidi E. L. Lischer 3,4, Judith Neukamm 1,5 ID , Enrique Rayo 1, Nicole Borel 6, Andreas Pospischil 6, Frank Rühli 1, Abigail S. Bouwman 1 and Michael G. Campana 1,7,* ID 1 Institute of Evolutionary Medicine, University of Zurich, 8057 Zurich, Switzerland; [email protected] (J.N.); [email protected] (E.R.); [email protected] (F.R.); [email protected] (A.S.B.) 2 Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, 0316 Oslo, Norway 3 Institute of Evolutionary Biology and Environmental Studies, University of Zurich, 8057 Zurich, Switzerland; [email protected] 4 Swiss Institute of Bioinformatics (SIB), 1015 Lausanne, Switzerland 5 Institute for Archaeological Sciences, University of Tübingen, 72070 Tübingen, Germany 6 Institute of Veterinary Pathology, University of Zurich, 8057 Zurich, Switzerland; [email protected] (N.B.); [email protected] (A.P.) 7 Center for Conservation Genomics, Smithsonian Conservation Biology Institute, Washington, DC 20008, USA * Correspondence: [email protected] (G.F.); [email protected] (M.G.C.); Tel.: +47-22-85-50-65 (G.F.); +1-202-633-4183 (M.G.C.) Received: 2 August 2018; Accepted: 30 August 2018; Published: 31 August 2018 Abstract: The reconstruction of ancient metagenomes from archaeological material, and their implication in human health and evolution, is one of the most recent advances in paleomicrobiological studies.
    [Show full text]
  • A Thesis Entitled Human Genome and Transcriptome Analysis with Next
    A Thesis entitled Human Genome and Transcriptome Analysis with Next-Generation Sequencing by Basil Khuder Submitted to the Graduate Faculty as partial fulfillment of the requirements for the Master of Science in Biomedical Sciences Degree in Bioinformatics, Proteomics, and Genomics ___________________________________________ Dr. Alexi Fedorov, Committee Chair ___________________________________________ Dr. David Kennedy, Committee Member ___________________________________________ Dr. Robert Blumenthal, Committee Member ___________________________________________ Dr. Amanda Bryant-Friedrich, Dean College of Graduate Studies The University of Toledo August 2017 Copyright 2017, Basil Khuder This document is copyrighted material. Under copyright law, no parts of this document may be reproduced without the expressed permission of the author. An Abstract of Human Genome and Transcriptome Analysis with Next-Generation Sequencing by Basil Khuder Submitted to the Graduate Faculty as partial fulfillment of the requirements for the Masters of Science in Biomedical Sciences Degree in Bioinformatics, Proteomics, Genomics The University of Toledo August 2017 Advancements in Next-Generation Sequencing technologies and steep declines in costs have enabled sequencing to occur at astronomical rates. With this technology, researchers have made great strides in progressing our understanding of the human genome. Additionally, this surplus of data has also opened the doors to many popular databases, such as such as NCBI’s Sequencing-Read Archive, to help host the enormity of files. The onslaught of data has also, however, created predicaments for scientists, as researchers are still trying to find the most optimal methods for using and processing these data to answer some of their most challenging questions. Here we present two approaches to normalize and analyze high-throughput data, that can help respond to questions about the human genome and transcriptome and demonstrate how constructing sequence analysis methods can produce substantial biological implications.
    [Show full text]