CSE182-L16. LW statistics/assembly
|
|
- Alyson Hensley
- 6 years ago
- Views:
Transcription
1 CSE182-L16 LW statistics/assembly
2 Silly Quiz Who are these people, and what is the occasion?
3 Genome Sequencing and Assembly
4 Sequencing A break at T is shown here. Measuring the lengths using electrophoresis allows us to get the position of each T The same can be done with every nucleotide. Fluorescent labeling can help separate different nucleotides November 13 Bafna
5 Automated detectors read the terminating bases. The signal decays after 1000 bases. November 13 Bafna
6 Sequencing Genomes: Clone by Clone Clones are constructed to span the entire length of the genome. These clones are ordered and oriented correctly (Mapping) Each clone is sequenced individually November 13 Bafna
7 Shotgun Sequencing Shotgun sequencing of clones was considered viable However, researchers in 1999 proposed shotgunning the entire genome. November 13 Bafna
8 Library Create vectors of the sequence and introduce them into bacteria. As bacteria multiply you will have many copies of the same clone. November 13 Bafna
9 Whole Genome Shotgun Break up the entire genome into pieces Sequence ends, and assemble using a computer LW statistics & Repeats argue against the success of such an approach Alternative: build a roadmap of the genome, with physical clones mapped for each region. Sequence each of the clones, and put them together November 13 Bafna
10 Whole Genome Shotgun Break up the entire genome into pieces Sequence ends, and assemble using a computer LW statistics & Repeats argue against the success of such an approach Alternative: build a roadmap of the genome, with physical clones mapped for each region. Sequence each of the clones, and put them together
11 Shotgun Sequencing Shotgun sequencing of clones was considered viable for small genomes However, researchers in 1999 proposed shotgunning the entire genome.
12 Massively parallel sequencing Sanger sequencing allows us to sequence <=1000 bp in one lane, up to 96 lanes, in one run. Today, we can sequence many Mbp in a single run
13
14
15 Questions Algorithmic: How do you put the genome back together from the pieces? Statistical? How many pieces do you need to sequence, etc.? The answer to the statistical questions had already been given in the context of mapping, by Lander and Waterman.
16 Lander Waterman Statistics The fragments are falling randomly on the genome Overlapping fragments form islands of contiguous sequence. Ideally, we want one island for each chromosome. How many fragments should we sequence? L G
17 Lander Waterman Statistics G = Genome Length L = Fragment Length N = Number of Fragments T = Required Overlap c = Coverage = LN/G α = N/G θ = T/L σ = 1- θ L G
18 LW statistics: questions As the coverage c increases, more and more areas of the genome are likely to be covered. Ideally, you want to see 1 island of overlapping contigs. Q1: What is the expected number of islands? The number should increase at first, and gradually decrease.
19 Analysis: Expected Number Islands Computing Expected # islands. Let X i =1 if an island ends at position i, X i =0 otherwise. Number of islands = i X i Expected # islands = E( i X i ) = i E(X i )
20 Prob. of an island ending at i L T i E(X i ) = Prob (Island ends at pos. i) =Prob(clone began at position i-l+1 AND no clone began in the next L-T positions) E(X i ) = α( 1 α) L T = αe cσ Expected # islands = E(X i ) = Gαe cσ = Ne cσ i
21 Computing # islands As the coverage c increases, more and more areas of the genome are likely to be covered. Ideally, you want to see 1 island. Q1: What is the expected number of islands? Ans: N exp(-cσ) The number increases at first, and gradually decreases.
22 Expected # of clones in an island Expected # of clones in an island = e cσ Q: How? Why do we care? Often, at the beginning of a genome project, we do not know the length of the genome. This equation helps us determine the length.
23 Problem 1: size of contigs Islands might simply be too small in length σ = (1-T/L) = (1-50/500) = 0.9, c = 8. #Islands = N e -cσ = 36K Size of an island = 54K Not enough to make it an acceptable assembly! PLUS, there is the problem of Repeats, Chimerism etc.
24 Assembly Basics Three main components: Overlap Layout Consensus
25 Overlap Given a pair of fragments s 1 and s 2, do they belong together? Yes, if a prefix of s 2 matches a suffix of s 1 How would you compute such a match?
26 Overlap S[i,j] = optimum score of an alignment of s 1 [1..i] against a substring of that starts anywhere, but ends in j. s 2 [*..j] j i The best prefix-suffix alignment is given by: Max i {S[i,n] }
27 Overlap Detection Compute the best prefixsuffix alignments between each pair of fragments. Keep the high-scoring ones as evidence of true overlap. What is the problem?
28 Overlap detection problem Consider the number of fragments. The LW statistics say that we need good coverage (c=8, 10) to get most of the base-pairs. G = 3000Mb, L=500 Coverage LN/G = 10 N = 10*3*10 9 /500 = 6*10 7 Number of comparisons needed = 3.6 * Number of alignments per minute=6 Number of compute nodes = 100 Time needed (Number of years) = Not good! (Only a small fraction are true overlaps) ( ) =11M
29 k-mer based overlap (Piegeonhole principle again) Consider a 25bp sequence. Expected number of occurrences in the genome 3*10 9 *4-25 = 2*10-6 A 25-bp sequence appears is unique to the genome! Two overlapping sequences should share a 25-mer Two non-overlapping sequences should not! 25bp
30 Sorting k-mers Build a list of k-mers that appear in the sequences and their reverse complements Create a record with 4 entries: K-mer Sequence number Position in the sequence Reverse complementation flag Sort a vector of these according to k-mer How many records per k-mer are expected? If number of records exceeds threshold, discard (why?) K-mer S.id Pos.
31 Coalesce k-mer hits into longer, gap-free partial alignments. These extended k-mer hits are saved. For each pair of sequences, form a directed graph. For each maximal path in the graph, construct an alignment. Refine alignment via banded DP Alignment module
32 Problem2: Size Islands might simply be too small in length σ = (1-T/L) = (1-50/500) = 0.9, c = 8. #Islands = N e -cσ = 36K Size of an island = 54K Not enough to make it an acceptable assembly! PLUS, there is the problem of Repeats, Chimerism etc.
33 Solution 2: Clones can have matepairs Recall that we sequence about 1000bp of the end of a clone If we sequenced both ends, we get extra information, particularly if we know the length of the original clone.
34 Mate Pairs Mate-pairs allow you to merge islands (contigs) into super-contigs
35 Super-contigs are quite large Make clones of truly predictable length. EX: 3 sets can be used: 2Kb, 10Kb and 50Kb. The variance in these lengths should be small. Use the mate-pairs to order and orient the contigs, and make super-contigs.
36 Problem 3: Repeats
37 Repeats & Chimerisms 40-50% of the human genome is made up of repetitive elements. Repeats can cause great problems in the assembly! Chimerism causes a clone to be from two different parts of the genome. Can again give a completely wrong assembly
38 Repeat detection Lander Waterman strikes again! The expected number of clones in a Repeat containing island is MUCH larger than in a non-repeat containing island (contig). Thus, every contig can be marked as Unique, or non-unique. In the first step, throw away the non-unique islands. Repeat
39 Detecting Repeat Contigs 1: Read Density Compute the log-odds ratio of two hypotheses: H1: The contig is from a unique region of the genome. The contig is from a region that is repeated at least twice
40 Detecting Chimeric reads Chimeric reads: Reads that contain sequence from two genomic locations. Good overlaps: G(a,b) if a,b overlap with a high score Transitive overlap: T(a,c) if G(a,b), and G(b,c) Find a point x across which only transitive overlaps occur. X is a point of chimerism
41 Whole genome shotgun Input: Shotgun sequence fragments (reads) Mate pairs Output: A single sequence created by consensus of overlapping reads First generation of assemblers did not include mate-pairs (Phrap, CAP..) Second generation: CA, Arachne, Euler
42 Assembly Use k-mers to detect potential overlaps Use alignments to build contig graphs Decide the unique contigs based on LW statistics Discard repeat contigs Break chimeric contigs Use mate-pairs to build scaffolds Fill gaps
43 Assembly Use k-mers to detect potential overlaps Use alignments to build contig graphs Decide the unique contigs based on LW statistics Discard repeat contigs Break chimeric contigs Use mate-pairs to build scaffolds
44 Consensus Derivation Consensus sequence is created by converting pairwise read alignments into multiple-read alignments
45 Summary Whole genome shotgun is now routine: Human, Mouse, Rat, Dog, Chimpanzee.. Many Prokaryotes (One can be sequenced in a day) Plant genomes: Arabidopsis, Rice Model organisms: Worm, Fly, Yeast A lot is not known about genome structure, organization and function. Comparative genomics offers low hanging fruit
46 Final exam syllabus Take home The entire course, but emphasis will be given to post-midterm lectures HMMs, Gene-finding, mass spectrometry, Micro-array analysis, genome sequencing and assembly
47 What we did not cover
The first generation DNA Sequencing
The first generation DNA Sequencing Slides 3 17 are modified from faperta.ugm.ac.id/newbie/download/pak_tar/.../instrument20072.ppt slides 18 43 are from Chengxiang Zhai at UIUC. The strand direction http://en.wikipedia.org/wiki/dna
More informationSequence Assembly and Alignment. Jim Noonan Department of Genetics
Sequence Assembly and Alignment Jim Noonan Department of Genetics james.noonan@yale.edu www.yale.edu/noonanlab The assembly problem >>10 9 sequencing reads 36 bp - 1 kb 3 Gb Outline Basic concepts in genome
More informationBENG 183 Trey Ideker. Genome Assembly and Physical Mapping
BENG 183 Trey Ideker Genome Assembly and Physical Mapping Reasons for sequencing Complete genome sequencing!!! Resequencing (Confirmatory) E.g., short regions containing single nucleotide polymorphisms
More informationA shotgun introduction to sequence assembly (with Velvet) MCB Brem, Eisen and Pachter
A shotgun introduction to sequence assembly (with Velvet) MCB 247 - Brem, Eisen and Pachter Hot off the press January 27, 2009 06:00 AM Eastern Time llumina Launches Suite of Next-Generation Sequencing
More informationDe Novo Assembly of High-throughput Short Read Sequences
De Novo Assembly of High-throughput Short Read Sequences Chuming Chen Center for Bioinformatics and Computational Biology (CBCB) University of Delaware NECC Third Skate Genome Annotation Workshop May 23,
More informationMate-pair library data improves genome assembly
De Novo Sequencing on the Ion Torrent PGM APPLICATION NOTE Mate-pair library data improves genome assembly Highly accurate PGM data allows for de Novo Sequencing and Assembly For a draft assembly, generate
More informationDe novo genome assembly with next generation sequencing data!! "
De novo genome assembly with next generation sequencing data!! " Jianbin Wang" HMGP 7620 (CPBS 7620, and BMGN 7620)" Genomics lectures" 2/7/12" Outline" The need for de novo genome assembly! The nature
More informationSequence assembly. Jose Blanca COMAV institute bioinf.comav.upv.es
Sequence assembly Jose Blanca COMAV institute bioinf.comav.upv.es Sequencing project Unknown sequence { experimental evidence result read 1 read 4 read 2 read 5 read 3 read 6 read 7 Computational requirements
More informationEach cell of a living organism contains chromosomes
COVER FEATURE Genome Sequence Assembly: Algorithms and Issues Algorithms that can assemble millions of small DNA fragments into gene sequences underlie the current revolution in biotechnology, helping
More informationGenome Sequencing-- Strategies
Genome Sequencing-- Strategies Bio 4342 Spring 04 What is a genome? A genome can be defined as the entire DNA content of each nucleated cell in an organism Each organism has one or more chromosomes that
More informationMapping strategies for sequence reads
Mapping strategies for sequence reads Ernest Turro University of Cambridge 21 Oct 2013 Quantification A basic aim in genomics is working out the contents of a biological sample. 1. What distinct elements
More informationIntroduction to Bioinformatics. Genome sequencing & assembly
Introduction to Bioinformatics Genome sequencing & assembly Genome sequencing & assembly p DNA sequencing How do we obtain DNA sequence information from organisms? p Genome assembly What is needed to put
More information1. A brief overview of sequencing biochemistry
Supplementary reading materials on Genome sequencing (optional) The materials are from Mark Blaxter s lecture notes on Sequencing strategies and Primary Analysis 1. A brief overview of sequencing biochemistry
More informationMolecular Biology: DNA sequencing
Molecular Biology: DNA sequencing Author: Prof Marinda Oosthuizen Licensed under a Creative Commons Attribution license. SEQUENCING OF LARGE TEMPLATES As we have seen, we can obtain up to 800 nucleotides
More informationGenome Assembly Using de Bruijn Graphs. Biostatistics 666
Genome Assembly Using de Bruijn Graphs Biostatistics 666 Previously: Reference Based Analyses Individual short reads are aligned to reference Genotypes generated by examining reads overlapping each position
More informationHuman SNP haplotypes. Statistics 246, Spring 2002 Week 15, Lecture 1
Human SNP haplotypes Statistics 246, Spring 2002 Week 15, Lecture 1 Human single nucleotide polymorphisms The majority of human sequence variation is due to substitutions that have occurred once in the
More informationGenomics AGRY Michael Gribskov Hock 331
Genomics AGRY 60000 Michael Gribskov gribskov@purdue.edu Hock 331 Computing Essentials Resources In this course we will assemble and annotate both genomic and transcriptomic sequence assemblies We will
More informationReview of whole genome methods
Review of whole genome methods Suffix-tree based MUMmer, Mauve, multi-mauve Gene based Mercator, multiple orthology approaches Dot plot/clustering based MUMmer 2.0, Pipmaker, LASTZ 10/3/17 0 Rationale:
More informationFinishing Drosophila Ananassae Fosmid 2728G16
Finishing Drosophila Ananassae Fosmid 2728G16 Kyle Jung March 8, 2013 Bio434W Professor Elgin Page 1 Abstract For my finishing project, I chose to finish fosmid 2728G16. This fosmid carries a segment of
More informationChromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Supplementary Material
Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions Joshua N. Burton 1, Andrew Adey 1, Rupali P. Patwardhan 1, Ruolan Qiu 1, Jacob O. Kitzman 1, Jay Shendure 1 1 Department
More informationDNA sequencing. Course Info
DNA sequencing EECS 458 CWRU Fall 2004 Readings: Pevzner Ch1-4 Adams, Fields & Venter (ISBN:0127170103) Serafim Batzoglou s slides Course Info Instructor: Jing Li 509 Olin Bldg Phone: X0356 Email: jingli@eecs.cwru.edu
More informationOutline. DNA Sequencing. Whole Genome Shotgun Sequencing. Sequencing Coverage. Whole Genome Shotgun Sequencing 3/28/15
Outline Introduction Lectures 22, 23: Sequence Assembly Spring 2015 March 27, 30, 2015 Sequence Assembly Problem Different Solutions: Overlap-Layout-Consensus Assembly Algorithms De Bruijn Graph Based
More informationBioinformatics Support of Genome Sequencing Projects. Seminar in biology
Bioinformatics Support of Genome Sequencing Projects Seminar in biology Introduction The Big Picture Biology reminder Enzyme for DNA manipulation DNA cloning DNA mapping Sequencing genomes Alignment of
More informationGenome Sequence Assembly
Genome Sequence Assembly Learning Goals: Introduce the field of bioinformatics Familiarize the student with performing sequence alignments Understand the assembly process in genome sequencing Introduction:
More informationContact us for more information and a quotation
GenePool Information Sheet #1 Installed Sequencing Technologies in the GenePool The GenePool offers sequencing service on three platforms: Sanger (dideoxy) sequencing on ABI 3730 instruments Illumina SOLEXA
More informationAssembly of Ariolimax dolichophallus using SOAPdenovo2
Assembly of Ariolimax dolichophallus using SOAPdenovo2 Charles Markello, Thomas Matthew, and Nedda Saremi Image taken from Banana Slug Genome Project, S. Weber SOAPdenovo Assembly Tool Short Oligonucleotide
More informationRestriction Site Mapping:
Restriction Site Mapping: In making genomic library the DNA is cut with rare cutting enzymes and large fragments of the size of 100,000 to 1000, 000bp. They are ligated to vectors such as Pacmid or YAC
More informationDe novo whole genome assembly
De novo whole genome assembly Lecture 1 Qi Sun Bioinformatics Facility Cornell University Data generation Sequencing Platforms Short reads: Illumina Long reads: PacBio; Oxford Nanopore Contiging/Scaffolding
More informationWorkflow of de novo assembly
Workflow of de novo assembly Experimental Design Clean sequencing data (trim adapter and low quality sequences) Run assembly software for contiging and scaffolding Evaluation of assembly Several iterations:
More informationLectures 18, 19: Sequence Assembly. Spring 2017 April 13, 18, 2017
Lectures 18, 19: Sequence Assembly Spring 2017 April 13, 18, 2017 1 Outline Introduction Sequence Assembly Problem Different Solutions: Overlap-Layout-Consensus Assembly Algorithms De Bruijn Graph Based
More informationHigh throughput omics and BIOINFORMATICS
High throughput omics and BIOINFORMATICS Giuseppe D'Auria Seville, February 2009 Genomes from isolated bacteria $ $ $ $ $ $ $ $ $$ $ $ $ $ $ $ $ se q se uen q c se uen ing q c se uen ing qu c en ing c
More informationIntroduction to Molecular Biology
Introduction to Molecular Biology Bioinformatics: Issues and Algorithms CSE 308-408 Fall 2007 Lecture 2-1- Important points to remember We will study: Problems from bioinformatics. Algorithms used to solve
More informationA Brief Introduction to Bioinformatics
A Brief Introduction to Bioinformatics Dan Lopresti Associate Professor Office PL 404B dal9@lehigh.edu February 2007 Slide 1 Motivation Biology easily has 500 years of exciting problems to work on. Donald
More informationIllumina (Solexa) Throughput: 4 Tbp in one run (5 days) Cheapest sequencing technology. Mismatch errors dominate. Cost: ~$1000 per human genme
Illumina (Solexa) Current market leader Based on sequencing by synthesis Current read length 100-150bp Paired-end easy, longer matepairs harder Error ~0.1% Mismatch errors dominate Throughput: 4 Tbp in
More informationFinishing Fosmid DMAC-27a of the Drosophila mojavensis third chromosome
Finishing Fosmid DMAC-27a of the Drosophila mojavensis third chromosome Ruth Howe Bio 434W 27 February 2010 Abstract The fourth or dot chromosome of Drosophila species is composed primarily of highly condensed,
More informationAxiom mydesign Custom Array design guide for human genotyping applications
TECHNICAL NOTE Axiom mydesign Custom Genotyping Arrays Axiom mydesign Custom Array design guide for human genotyping applications Overview In the past, custom genotyping arrays were expensive, required
More informationMODULE 1: INTRODUCTION TO THE GENOME BROWSER: WHAT IS A GENE?
MODULE 1: INTRODUCTION TO THE GENOME BROWSER: WHAT IS A GENE? Lesson Plan: Title Introduction to the Genome Browser: what is a gene? JOYCE STAMM Objectives Demonstrate basic skills in using the UCSC Genome
More informationDe novo whole genome assembly
De novo whole genome assembly Lecture 1 Qi Sun Minghui Wang Bioinformatics Facility Cornell University DNA Sequencing Platforms Illumina sequencing (100 to 300 bp reads) Overlapping reads ~180bp fragment
More informationGenetics Lecture 21 Recombinant DNA
Genetics Lecture 21 Recombinant DNA Recombinant DNA In 1971, a paper published by Kathleen Danna and Daniel Nathans marked the beginning of the recombinant DNA era. The paper described the isolation of
More informationSequencing technologies. Jose Blanca COMAV institute bioinf.comav.upv.es
Sequencing technologies Jose Blanca COMAV institute bioinf.comav.upv.es Outline Sequencing technologies: Sanger 2nd generation sequencing: 3er generation sequencing: 454 Illumina SOLiD Ion Torrent PacBio
More informationFinished (Almost) Sequence of Drosophila littoralis Chromosome 4 Fosmid Clone XAAA73. Seth Bloom Biology 4342 March 7, 2004
Finished (Almost) Sequence of Drosophila littoralis Chromosome 4 Fosmid Clone XAAA73 Seth Bloom Biology 4342 March 7, 2004 Summary: I successfully sequenced Drosophila littoralis fosmid clone XAAA73. The
More informationSequencing technologies. Jose Blanca COMAV institute bioinf.comav.upv.es
Sequencing technologies Jose Blanca COMAV institute bioinf.comav.upv.es Outline Sequencing technologies: Sanger 2nd generation sequencing: 3er generation sequencing: 454 Illumina SOLiD Ion Torrent PacBio
More informationLawrence Berkeley National Laboratory Lawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Lawrence Berkeley National Laboratory Title Glomus intraradices: Initial Whole-Genome Shotgun Sequencing and Assembly Results Permalink https://escholarship.org/uc/item/4c13k1dh
More informationThis is a closed book, closed note exam. No calculators, phones or any electronic device are allowed.
MCB 104 MIDTERM #2 October 23, 2013 ***IMPORTANT REMINDERS*** Print your name and ID# on every page of the exam. You will lose 0.5 point/page if you forget to do this. Name KEY If you need more space than
More informationNext Gen Sequencing. Expansion of sequencing technology. Contents
Next Gen Sequencing Contents 1 Expansion of sequencing technology 2 The Next Generation of Sequencing: High-Throughput Technologies 3 High Throughput Sequencing Applied to Genome Sequencing (TEDed CC BY-NC-ND
More informationComputational Biology I LSM5191
Computational Biology I LSM5191 Lecture 5 Notes: Genetic manipulation & Molecular Biology techniques Broad Overview of: Enzymatic tools in Molecular Biology Gel electrophoresis Restriction mapping DNA
More informationCS 5984: Application of Basic Clustering Algorithms to Find Expression Modules in Cancer
CS 5984: Application of Basic Clustering Algorithms to Find Expression Modules in Cancer T. M. Murali January 31, 2006 Innovative Application of Hierarchical Clustering A module map showing conditional
More informationGenscan. The Genscan HMM model Training Genscan Validating Genscan. (c) Devika Subramanian,
Genscan The Genscan HMM model Training Genscan Validating Genscan (c) Devika Subramanian, 2009 96 Gene structure assumed by Genscan donor site acceptor site (c) Devika Subramanian, 2009 97 A simple model
More informationContigs Built with Fingerprints, Markers, and FPC V4.7
Methods Contigs Built with Fingerprints, Markers, and FPC V4.7 Carol Soderlund, 1,3 Sean Humphray, 2 Andrew Dunham, 2 and Lisa French 2 1 Clemson University Genomic Institute, Clemson, South Carolina 29634-5808,
More informationMultiple choice questions (numbers in brackets indicate the number of correct answers)
1 Multiple choice questions (numbers in brackets indicate the number of correct answers) February 1, 2013 1. Ribose is found in Nucleic acids Proteins Lipids RNA DNA (2) 2. Most RNA in cells is transfer
More informationExpressed Sequence Tags: Clustering and Applications
12 Expressed Sequence Tags: Clustering and Applications Anantharaman Kalyanaraman Iowa State University Srinivas Aluru Iowa State University 12.1 Introduction... 12-1 12.2 Sequencing ESTs... 12-2 12.3
More informationshort read genome assembly Sorin Istrail CSCI1820 Short-read genome assembly algorithms 3/6/2014
1 short read genome assembly Sorin Istrail CSCI1820 Short-read genome assembly algorithms 3/6/2014 2 Genomathica Assembler Mathematica notebook for genome assembly simulation Assembler can be found at:
More informationBiology 201 (Genetics) Exam #3 120 points 20 November Read the question carefully before answering. Think before you write.
Name KEY Section Biology 201 (Genetics) Exam #3 120 points 20 November 2006 Read the question carefully before answering. Think before you write. You will have up to 50 minutes to take this exam. After
More informationOutline. Gene Finding Questions. Recap: Prokaryotic gene finding Eukaryotic gene finding The human gene complement Regulation
Tues, Nov 29: Gene Finding 1 Online FCE s: Thru Dec 12 Thurs, Dec 1: Gene Finding 2 Tues, Dec 6: PS5 due Project presentations 1 (see course web site for schedule) Thurs, Dec 8 Final papers due Project
More informationWhole-Genome Sequencing and Assembly with High- Throughput, Short-Read Technologies
Whole-Genome Sequencing and Assembly with High- Throughput, Short-Read Technologies Andreas Sundquist 1 *, Mostafa Ronaghi 2, Haixu Tang 3, Pavel Pevzner 4, Serafim Batzoglou 1 1 Department of Computer
More informationGenetics and Genomics in Medicine Chapter 3. Questions & Answers
Genetics and Genomics in Medicine Chapter 3 Multiple Choice Questions Questions & Answers Question 3.1 Which of the following statements, if any, is false? a) Amplifying DNA means making many identical
More informationDe novo Genome Assembly
De novo Genome Assembly A/Prof Torsten Seemann Winter School in Mathematical & Computational Biology - Brisbane, AU - 3 July 2017 Introduction The human genome has 47 pieces MT (or XY) The shortest piece
More informationConnect-A-Contig Paper version
Teacher Guide Connect-A-Contig Paper version Abstract Students align pieces of paper DNA strips based on the distance between markers to generate a DNA consensus sequence. The activity helps students see
More informationCOPE: An accurate k-mer based pair-end reads connection tool to facilitate genome assembly
Bioinformatics Advance Access published October 8, 2012 COPE: An accurate k-mer based pair-end reads connection tool to facilitate genome assembly Binghang Liu 1,2,, Jianying Yuan 2,, Siu-Ming Yiu 1,3,
More informationMeasuring transcriptomes with RNA-Seq. BMI/CS 776 Spring 2016 Anthony Gitter
Measuring transcriptomes with RNA-Seq BMI/CS 776 www.biostat.wisc.edu/bmi776/ Spring 2016 Anthony Gitter gitter@biostat.wisc.edu Overview RNA-Seq technology The RNA-Seq quantification problem Generative
More informationCAP 5510: Introduction to Bioinformatics CGS 5166: Bioinformatics Tools. Giri Narasimhan
CAP 5510: Introduction to Bioinformatics CGS 5166: Bioinformatics Tools Giri Narasimhan ECS 254; Phone: x3748 giri@cis.fiu.edu www.cis.fiu.edu/~giri/teach/bioinfs15.html Gene Expression q Process of transcription
More informationHunting Down the Papaya Transgenes
Hunting Down the Papaya Transgenes Michael Schatz Center for Bioinformatics and Computational Biology University of Maryland January 16, 2008 PAG XVI Papaya Overview Carica papaya from the order Brassicales
More informationGene Identification in silico
Gene Identification in silico Nita Parekh, IIIT Hyderabad Presented at National Seminar on Bioinformatics and Functional Genomics, at Bioinformatics centre, Pondicherry University, Feb 15 17, 2006. Introduction
More informationNext-Generation Sequencing. Technologies
Next-Generation Next-Generation Sequencing Technologies Sequencing Technologies Nicholas E. Navin, Ph.D. MD Anderson Cancer Center Dept. Genetics Dept. Bioinformatics Introduction to Bioinformatics GS011062
More informationNCBI web resources I: databases and Entrez
NCBI web resources I: databases and Entrez Yanbin Yin Most materials are downloaded from ftp://ftp.ncbi.nih.gov/pub/education/ 1 Homework assignment 1 Two parts: Extract the gene IDs reported in table
More informationON USING DNA DISTANCES AND CONSENSUS IN REPEATS DETECTION
ON USING DNA DISTANCES AND CONSENSUS IN REPEATS DETECTION Petre G. POP Technical University of Cluj-Napoca, Romania petre.pop@com.utcluj.ro Abstract: Sequence repeats are the simplest form of regularity
More information3. human genomics clone genes associated with genetic disorders. 4. many projects generate ordered clones that cover genome
Lectures 30 and 31 Genome analysis I. Genome analysis A. two general areas 1. structural 2. functional B. genome projects a status report 1. 1 st sequenced: several viral genomes 2. mitochondria and chloroplasts
More informationLecture 8: Sequencing and SNP. Sept 15, 2006
Lecture 8: Sequencing and SNP Sept 15, 2006 Announcements Random questioning during literature discussion sessions starts next week for real! Schedule changes Moved QTL lecture up Removed landscape genetics
More informationGenome Assembly. J Fass UCD Genome Center Bioinformatics Core Friday September, 2015
Genome Assembly J Fass UCD Genome Center Bioinformatics Core Friday September, 2015 From reads to molecules What s the Problem? How to get the best assemblies for the smallest expense (sequencing) and
More informationGenomics and Gene Recognition Genes and Blue Genes
Genomics and Gene Recognition Genes and Blue Genes November 3, 2004 Eukaryotic Gene Structure eukaryotic genomes are considerably more complex than those of prokaryotes eukaryotic cells have organelles
More informationmeasuring gene expression December 5, 2017
measuring gene expression December 5, 2017 transcription a usually short-lived RNA copy of the DNA is created through transcription RNA is exported to the cytoplasm to encode proteins some types of RNA
More informationEcole de Bioinforma(que AVIESAN Roscoff 2014 GALAXY INITIATION. A. Lermine U900 Ins(tut Curie, INSERM, Mines ParisTech
GALAXY INITIATION A. Lermine U900 Ins(tut Curie, INSERM, Mines ParisTech How does Next- Gen sequencing work? DNA fragmentation Size selection and clonal amplification Massive parallel sequencing ACCGTTTGCCG
More informationMeasuring transcriptomes with RNA-Seq
Measuring transcriptomes with RNA-Seq BMI/CS 776 www.biostat.wisc.edu/bmi776/ Spring 2017 Anthony Gitter gitter@biostat.wisc.edu These slides, excluding third-party material, are licensed under CC BY-NC
More informationRecombinant DNA Technology. The Role of Recombinant DNA Technology in Biotechnology. yeast. Biotechnology. Recombinant DNA technology.
PowerPoint Lecture Presentations prepared by Mindy Miller-Kittrell, North Carolina State University C H A P T E R 8 Recombinant DNA Technology The Role of Recombinant DNA Technology in Biotechnology Biotechnology?
More informationComparative Bioinformatics. BSCI348S Fall 2003 Midterm 1
BSCI348S Fall 2003 Midterm 1 Multiple Choice: select the single best answer to the question or completion of the phrase. (5 points each) 1. The field of bioinformatics a. uses biomimetic algorithms to
More informationAnalysis of large deletions in human-chimp genomic alignments. Erika Kvikstad BioInformatics I December 14, 2004
Analysis of large deletions in human-chimp genomic alignments Erika Kvikstad BioInformatics I December 14, 2004 Outline Mutations, mutations, mutations Project overview Strategy: finding, classifying indels
More informationImproving the Accuracy of Base Calls and Error Predictions for GS 20 DNA Sequence Data
Improving the Accuracy of Base Calls and Error Predictions for GS 20 DNA Sequence Data Justin S. Hogg Department of Computational Biology University of Pittsburgh Pittsburgh, PA 15213 jsh32@pitt.edu Abstract
More informationHaploid Assembly of Diploid Genomes
Haploid Assembly of Diploid Genomes Challenges, Trials, Tribulations 13 October 2011 İnanç Birol Assembly By Short Sequencing IEEE InfoVis 2009 2 3 in Literature ~40 citations on tool comparisons ~20 citations
More informationDepartment of Computer Science and Engineering, University of
Optimizing the BAC-End Strategy for Sequencing the Human Genome Richard M. Karp Ron Shamir y April 4, 999 Abstract The rapid increase in human genome sequencing eort and the emergence of several alternative
More informationGenomics and Transcriptomics of Spirodela polyrhiza
Genomics and Transcriptomics of Spirodela polyrhiza Doug Bryant Bioinformatics Core Facility & Todd Mockler Group, Donald Danforth Plant Science Center Desired Outcomes High-quality genomic reference sequence
More informationMotif Discovery from Large Number of Sequences: a Case Study with Disease Resistance Genes in Arabidopsis thaliana
Motif Discovery from Large Number of Sequences: a Case Study with Disease Resistance Genes in Arabidopsis thaliana Irfan Gunduz, Sihui Zhao, Mehmet Dalkilic and Sun Kim Indiana University, School of Informatics
More informationLecture 11: Gene Prediction
Lecture 11: Gene Prediction Study Chapter 6.11-6.14 1 Gene: A sequence of nucleotides coding for protein Gene Prediction Problem: Determine the beginning and end positions of genes in a genome Where are
More informationSept 2. Structure and Organization of Genomes. Today: Genetic and Physical Mapping. Sept 9. Forward and Reverse Genetics. Genetic and Physical Mapping
Sept 2. Structure and Organization of Genomes Today: Genetic and Physical Mapping Assignments: Gibson & Muse, pp.4-10 Brown, pp. 126-160 Olson et al., Science 245: 1434 New homework:due, before class,
More informationBIOINFORMATICS 1 SEQUENCING TECHNOLOGY. DNA story. DNA story. Sequencing: infancy. Sequencing: beginnings 26/10/16. bioinformatic challenges
BIOINFORMATICS 1 or why biologists need computers SEQUENCING TECHNOLOGY bioinformatic challenges http://www.bioinformatics.uni-muenster.de/teaching/courses-2012/bioinf1/index.hbi Prof. Dr. Wojciech Makałowski"
More informationRestriction Enzymes (Site-Specific Endonuclease) Enzymes that recognize and cleave dsdna in a highly sequence specific manner.
Enzymes Restriction Enzymes (Site-Specific Endonuclease) Enzymes that recognize and cleave dsdna in a highly sequence specific manner. Generally recognize an inverted repeat sequence 4, 6, or 8 base pairs
More informationBasic Bioinformatics: Homology, Sequence Alignment,
Basic Bioinformatics: Homology, Sequence Alignment, and BLAST William S. Sanders Institute for Genomics, Biocomputing, and Biotechnology (IGBB) High Performance Computing Collaboratory (HPC 2 ) Mississippi
More informationBST227 Introduction to Statistical Genetics. Lecture 8: Variant calling from high-throughput sequencing data
BST227 Introduction to Statistical Genetics Lecture 8: Variant calling from high-throughput sequencing data 1 PC recap typical genome Differs from the reference genome at 4-5 million sites ~85% SNPs ~15%
More informationThe Genome Analysis Centre. Building Excellence in Genomics and Computational Bioscience
Building Excellence in Genomics and Computational Bioscience Wheat genome sequencing: an update from TGAC Sequencing Technology Development now Plant & Microbial Genomics Group Leader Matthew Clark matt.clark@tgac.ac.uk
More informationHigh Throughput Sequencing Technologies. UCD Genome Center Bioinformatics Core Monday 15 June 2015
High Throughput Sequencing Technologies UCD Genome Center Bioinformatics Core Monday 15 June 2015 Sequencing Explosion www.genome.gov/sequencingcosts http://t.co/ka5cvghdqo Sequencing Explosion 2011 PacBio
More informationGenome Reassembly From Fragments. 28 March 2013 OSU CSE 1
Genome Reassembly From Fragments 28 March 2013 OSU CSE 1 Genome A genome is the encoding of hereditary information for an organism in its DNA The mathematical model of a genome is a string of character,
More informationALGORITHMS IN BIO INFORMATICS. Chapman & Hall/CRC Mathematical and Computational Biology Series A PRACTICAL INTRODUCTION. CRC Press WING-KIN SUNG
Chapman & Hall/CRC Mathematical and Computational Biology Series ALGORITHMS IN BIO INFORMATICS A PRACTICAL INTRODUCTION WING-KIN SUNG CRC Press Taylor & Francis Group Boca Raton London New York CRC Press
More informationLecture Four. Molecular Approaches I: Nucleic Acids
Lecture Four. Molecular Approaches I: Nucleic Acids I. Recombinant DNA and Gene Cloning Recombinant DNA is DNA that has been created artificially. DNA from two or more sources is incorporated into a single
More informationAdvisors: Prof. Louis T. Oliphant Computer Science Department, Hiram College.
Author: Sulochana Bramhacharya Affiliation: Hiram College, Hiram OH. Address: P.O.B 1257 Hiram, OH 44234 Email: bramhacharyas1@my.hiram.edu ACM number: 8983027 Category: Undergraduate research Advisors:
More informationProblem Set 8. Answer Key
MCB 102 University of California, Berkeley August 11, 2009 Isabelle Philipp Online Document Problem Set 8 Answer Key 1. The Genetic Code (a) Are all amino acids encoded by the same number of codons? no
More informationFuzzy Genome Sequence Assembly for Single and Environmental Genomes
Fuzzy Genome Sequence Assembly for Single and Environmental Genomes Sara Nasser, Adrienne Breland, Frederick C. Harris Jr., Monica Nicolescu, and Gregory L. Vert Department of Computer Science & Engineering,
More informationFinishing of DELE Drosophila elegans has been sequenced using Roche 454 pyrosequencing and Illumina
Sarah Swiezy Dr. Elgin, Dr. Shaffer Bio 434W 27 February 2015 Finishing of DELE8596009 Abstract Drosophila elegans has been sequenced using Roche 454 pyrosequencing and Illumina technology. DELE8596009,
More informationGenome Assembly and Annotation of Isochrysis Galbana
Genome Assembly and Annotation of Isochrysis Galbana By: Yi Wang Institution: California State University San Marcos Date: May 14, 2014 Abstract Isochrysis Galbana is a species of cocoolithophores, which
More informationRNA-Sequencing analysis
RNA-Sequencing analysis Markus Kreuz 25. 04. 2012 Institut für Medizinische Informatik, Statistik und Epidemiologie Content: Biological background Overview transcriptomics RNA-Seq RNA-Seq technology Challenges
More informationMCB 102 University of California, Berkeley July 28-30, Problem Set 6
MCB 102 University of California, Berkeley July 28-30, 2009 Isabelle Philipp Handout Problem Set 6 The answer key will be posted by Tuesday July 28. Try to solve the problem sets always first without the
More information