DOCSLIB.ORG

Whole Genome Sequencing, De Novo Assembly and Phenotypic Profiling

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Whole Genome Sequencing, De Novo Assembly and Phenotypic Profiling INVESTIGATION Whole Genome Sequencing, de Novo Assembly and Phenotypic Profiling for the New Budding Yeast Species Saccharomyces jurei Samina Naseeb,1,2 Haya Alsammar,1 Tim Burgis, Ian Donaldson, Norman Knyazev, Christopher Knight, and Daniela Delneri2 Manchester Institute of Biotechnology, The University of Manchester, UK, M1 7DN ORCID IDs: 0000-0003-3599-5813 (S.N.); 0000-0001-9815-4267 (C.K.); 0000-0001-8070-411X (D.D.) ABSTRACT Saccharomyces sensu stricto complex consist of yeast species, which are not only important in KEYWORDS the fermentation industry but are also model systems for genomic and ecological analysis. Here, we present evolution the complete genome assemblies of Saccharomyces jurei, a newly discovered Saccharomyces sensu stricto fitness species from high altitude oaks. Phylogenetic and phenotypic analysis revealed that S. jurei is more closely PacBio related to S. mikatae, than S. cerevisiae, and S. paradoxus. The karyotype of S. jurei presents two reciprocal translocation chromosomal translocations between chromosome VI/VII and I/XIII when compared to the S. cerevisiae Saccharomyces genome. Interestingly, while the rearrangement I/XIII is unique to S. jurei, the other is in common with S. mikatae strain IFO1815, suggesting shared evolutionary history of this species after the split between S. cerevisiae and S. mikatae. The number of Ty elements differed in the new species, with a higher number of Ty elements present in S. jurei than in S. cerevisiae. Phenotypically, the S. jurei strain NCYC 3962 has relatively higher fitness than the other strain NCYC 3947T under most of the environmental stress conditions tested and showed remarkably increased fitness in higher concentration of acetic acid compared to the other sensu stricto species. Both strains were found to be better adapted to lower temperatures compared to S. cerevisiae. Saccharomyces sensu stricto yeasts, currently comprise eight species: surrounding soil at an altitude of 1000 m above sea level in Saint Auban, S. cerevisiae, S. paradoxus, S. uvarum, S. mikatae, S. kudriavzevii, France (Naseeb et al. 2017b). It is known that species within the sensu S. arboricola, S. eubayanus, S. jurei (Martini and Martini 1987; Wang stricto group are reproductively isolated and possess post- zygotic bar- and Bai 2008; Naumov et al. 2000; Naumov et al. 1995a; Naumov et al. riers (Naumov 1987). Moreover, yeasts within this group exhibit almost 1995b; Libkind et al. 2011; Naseeb et al. 2017b) and two natural hy- identical karyotypes with 16 chromosomes (Cardinali and Martini brids: S. pastorianus (Masneuf et al. 1998; Querol and Bond 2009) and 1994; Carle and Olson 1985; Naumov et al. 1996). S. bayanus (Nguyen et al. 2011). Saccharomyces jurei is the latest ad- In the modern era of yeast genetics, the advances in sequencing dition to the sensu stricto clade and was isolated from oak tree bark and technology have lead to the whole genome sequencing of many Saccharomyces sensu stricto species (S. cerevisiae, S. bayanus var. Copyright © 2018 Naseeb et al. uvarum, S. kudriavzevii, S. mikatae, S. paradoxus, S. eubayanus doi: https://doi.org/10.1534/g3.118.200476 and S. arboricola)(Libkindet al. 2011; Liti et al. 2013; Cliften Manuscript received June 4, 2018; accepted for publication July 11, 2018; et al. 2003; Kellis et al. 2003; Scannell et al. 2011; Casaregola et al. published Early Online August 14, 2018. 2000). To date, more than 1000 S. cerevisiae strains belonging to This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/ different geographical and environmental origins have been se- licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction quenced and assembled (Engel and Cherry 2013; Peter et al. in any medium, provided the original work is properly cited. 2018). The availability of sequencing data from multiple strains of Supplemental material available at Figshare: https://figshare.com/s/ Saccharomycotina yeast species has enhanced our understanding of 60bbbc1e98886077182a. 1These authors contributed equally in this work biological mechanisms and comparative genomics. Researchers are 2Corresponding authors: Samina Naseeb: [email protected]; now combining comparative genomics with population ecology to Daniela Delneri: [email protected] better understand the genetic variations, taxonomy, evolution and Volume 8 | September 2018 | 2967 n Table 1 Strains used in this study Species Strain number References S. jurei NCYC 3947T (Naseeb et al. 2017b) NCYC 3962 S. cerevisiae NCYC 505T (Vaughan Martini and Kurtzman 1985) S. paradoxus CBS 432T (Naumov 1987) S. mikatae NCYC 2888T (IFO 1815T) (Yamada et al. 1993) S. kudriavzevii NCYC 2889T (IFO 1802T) (Yamada et al. 1993) S. arboricola CBS 10644T (Wang and Bai 2008) S. eubayanus PYCC 6148T (CBS 12357T) (Libkind et al. 2011) S. uvarum NCYC 2669 (CBS 7001) (Pulvirenti et al. 2000) S. pastorianus NCYC 329T (CBS 1538T) (Martini and Martini 1987) speciation of yeast strains in nature. Genome variation provides the EDTA and1M DTT) containing 1 mg/ml lyticase. The cells were in- raw material for evolution, and may arise by various mechanisms cubated at 37° for 1 hr. Following incubation, the cells were mixed with including gene duplication, horizontal gene transfer, hybridization stop solution (3M NaCl, 100mM Tris pH 7.5 and 20mM EDTA) and and micro and macro rearrangements (Fischer et al. 2001; Seoighe 60 ml of 10% SDS. The cell suspension was vortexed and mixed with et al. 2000; Lynch 2002; Hall et al. 2005; Naseeb et al. 2017a; Naseeb 500 ml phenol-chloroform. The samples were centrifuged at 13000 rpm et al. 2016; Naseeb and Delneri 2012). Synteny conservation studies for 2 min to separate the aqueous phase from the organic phase. The have shown highly variable rates of genetic rearrangements between upper aqueous phase was transferred to a clean 1.5 ml tube and phenol- individual lineages both in vertebrates and in yeasts (Bourque et al. chloroform step was repeated twice until a white interface was no 2005; Fischer et al. 2006; Vakirlis et al. 2016). This genome variation longer present. The aqueous phase was washed with 1 ml absolute is a means of evolutionary adaptation to environmental changes. An ethanol by centrifugation at 13000 rpm for 10 min. The pellet was air understanding of the genetic machinery linked to phenotypic var- dried and resuspended in 30 ml of sterile milliQ water. iation provides knowledge of the distribution of Saccharomyces spe- Genomic DNA for PacBio sequencing was extracted using Qiagen cies in different environments, and their ability to withstand specific Genomic-tip 20/G kit (cat. No. 10223) following manufacturer’srec- conditions (Goddard and Greig 2015; Jouhten et al. 2016; Brice et al. ommended instructions. The yield of all DNA samples was assessed by 2018; Peter et al. 2018). the nanodrop spectrophotometer (ND-1000) and by Qubit 2.0 fluo- Recently, we isolated two strains (NCYC 3947T and NCYC 3962) of rometer (catalog no. Q32866). Purity and integrity of DNA was Saccharomyces jurei from Quercus robur bark and surrounding soil checked by electrophoresis on 0.8% (w/v) agarose gel and by calculating (Naseeb et al. 2017b). The initial sequencing of ITS1, D1/D2 and seven the A260/A280 ratios. other nuclear genes showed that both strains of S. jurei were closely related to S. mikatae and S. paradoxus and grouped in Saccharomyces Library preparation for Illumina and PacBio sequencing sensu stricto complex. We also showed that S. jurei can readily hybridize Paired end whole-genomesequencingwasperformedusingthe Illumina with other sensu stricto species but the resulting hybrids were sterile HiSeq platform. FastQC (Babraham Bioinformatics) was used to apply (Naseeb et al. 2017b). Here, we represent high quality de novo sequence quality control to sequence reads, alignment of the reads was performed and assembly of both strains (NCYC 3947T and NCYC 3962) of S. jurei. using BOWTIE2 (Langmead and Salzberg 2012) and post-processed The phylogenetic analysis placed S. jurei in the sensu stricto clade, in a using SAMTOOLS (Li et al. 2009). small monophyletic group with S. mikatae. By combining Illumina For Pacbio sequencing, genomic DNA (10 mg) of NCYC 3947T and HiSeq and PacBio data, we were able to assemble full chromosomes NCYC 3962 strains was first DNA damage repaired, sheared with and carry out synteny analysis. Moreover, we show that S. jurei NCYC Covaris G-tube, end repaired and exonuclease treated. SMRTbell li- 3962 had higher fitness compared to NCYC 3947T under different brary (10-20kb size) was prepared by ligation of hairpin adaptors at environmental conditions. Fitness of S. jurei strains at different tem- both ends according to PacBio recommended procedure (Pacific Bio- peratures showed that it was able to grow at wider range of tempera- science, No: 100-259-100). The resulting library was then size selected tures (12°-37°). using Blue Pippin with 7-10kb cut-off. Sequencing run was performed on PacBio RS II using P6/C4 chemistry for 4 hr. The genome was MATERIAL AND METHODS assembled using SMRT analysis and HGAP3 pipeline was made using default settings. Yeast strains Strains used in this study are presented in Table 1. All strains were Genome assembly, annotation, orthology and grown and maintained on YPDA (1% w/v yeast extract, 2% w/v Bacto- chromosomal structural plots peptone, 2% v/v glucose and 2% w/v agar). Species names and strains The PacBio sequences were assembled using hierarchical genome- number are stated in Table 1. assembly process (HGAP) (Chin et al. 2013). Protein coding gene models were predicted using Augustus (Stanke and Morgenstern DNA Extraction 2005) and the Yeast Genome Annotation Pipeline (Byrne and Wolfe For Illumina Hiseq, the total DNA was extracted from an overnight 2005). In addition, protein sequences from other Saccharomyces species grown culture of yeast strains by using the standard phenol/chloroform were aligned to the genome assembly using tblastn (Gertz et al.
Load more

Recommended publications
  • De Novo Genomic Analyses for Non-Model Organisms: an Evaluation of Methods Across a Multi-Species Data Set
    De novo genomic analyses for non-model organisms: an evaluation of methods across a multi-species data set Sonal Singhal [email protected] Museum of Vertebrate Zoology University of California, Berkeley 3101 Valley Life Sciences Building Berkeley, California 94720-3160 Department of Integrative Biology University of California, Berkeley 1005 Valley Life Sciences Building Berkeley, California 94720-3140 June 5, 2018 Abstract High-throughput sequencing (HTS) is revolutionizing biological research by enabling sci- entists to quickly and cheaply query variation at a genomic scale. Despite the increasing ease of obtaining such data, using these data effectively still poses notable challenges, especially for those working with organisms without a high-quality reference genome. For every stage of analysis – from assembly to annotation to variant discovery – researchers have to distinguish technical artifacts from the biological realities of their data before they can make inference. In this work, I explore these challenges by generating a large de novo comparative transcriptomic dataset data for a clade of lizards and constructing a pipeline to analyze these data. Then, using a combination of novel metrics and an externally validated variant data set, I test the efficacy of my approach, identify areas of improvement, and propose ways to minimize these errors. I arXiv:1211.1737v1 [q-bio.GN] 8 Nov 2012 find that with careful data curation, HTS can be a powerful tool for generating genomic data for non-model organisms. Keywords: de novo assembly, transcriptomes, suture zones, variant discovery, annotation Running Title: De novo genomic analyses for non-model organisms 1 Introduction High-throughput sequencing (HTS) is poised to revolutionize the field of evolutionary genetics by enabling researchers to assay thousands of loci for organisms across the tree of life.
    [Show full text]
  • MATCH-G Program
    MATCH-G Program The MATCH-G toolset MATCH-G (Mutational Analysis Toolset Comparing wHole Genomes) Download the Toolset The toolset is prepared for use with pombe genomes and a sample genome from Sanger is included. However, it is written in such a way as to allow it to utilize any set or number of chromosomes. Simply follow the naming convention "chromosomeX.contig.embl" where X=1,2,3 etc. For use in terminal windows in Mac OSX or Unix-like environments where Perl is present by default. Toolset Description Version History 2.0 – Bug fixes related to scaling to other genomes. First version of GUI interface. 1.0 – Initial Release. Includes build, alignment, snp, copy, and gap resolution routines in original form. Please note this project in under development and will undergo significant changes. Project Description The MATCH-G toolset has been developed to facilitate the evaluation of whole genome sequencing data from yeasts. Currently, the toolset is written for pombe strains, however the toolset is easily scalable to other yeasts or other sequenced genomes. The included tools assist in the identification of SNP mutations, short (and long) insertion/deletions, and changes in gene/region copy number between a parent strain and mutant or revertant strain. Currently, the toolset is run from the command line in a unix or similar terminal and requires a few additional programs to run noted below. Installation It is suggested that a separate folder be generated for the toolset and generated files. Free disk space proportional to the original read datasets is recommended. The toolset utilizes and requires a few additional free programs: Bowtie rapid alignment software: http://bowtie-bio.sourceforge.net/index.shtml (Langmead B, Trapnell C, Pop M, Salzberg SL.
    [Show full text]
  • Large Scale Genomic Rearrangements in Selected Arabidopsis Thaliana T
    bioRxiv preprint doi: https://doi.org/10.1101/2021.03.03.433755; this version posted March 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. 1 Large scale genomic rearrangements in selected 2 Arabidopsis thaliana T-DNA lines are caused by T-DNA 3 insertion mutagenesis 4 5 Boas Pucker1,2+, Nils Kleinbölting3+, and Bernd Weisshaar1* 6 1 Genetics and Genomics of Plants, Center for Biotechnology (CeBiTec), Bielefeld University, 7 Sequenz 1, 33615 Bielefeld, Germany 8 2 Evolution and Diversity, Department of Plant Sciences, University of Cambridge, Cambridge, 9 United Kingdom 10 3 Bioinformatics Resource Facility, Center for Biotechnology (CeBiTec, Bielefeld University, 11 Sequenz 1, 33615 Bielefeld, Germany 12 + authors contributed equally 13 * corresponding author: Bernd Weisshaar 14 15 BP: [email protected], ORCID: 0000-0002-3321-7471 16 NK: [email protected], ORCID: 0000-0001-9124-5203 17 BW: [email protected], ORCID: 0000-0002-7635-3473 18 19 20 21 22 23 24 25 26 27 28 29 30 page 1 bioRxiv preprint doi: https://doi.org/10.1101/2021.03.03.433755; this version posted March 7, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
    [Show full text]
  • An Open-Sourced Bioinformatic Pipeline for the Processing of Next-Generation Sequencing Derived Nucleotide Reads
    bioRxiv preprint doi: https://doi.org/10.1101/2020.04.20.050369; this version posted May 28, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. An open-sourced bioinformatic pipeline for the processing of Next-Generation Sequencing derived nucleotide reads: Identification and authentication of ancient metagenomic DNA Thomas C. Collin1, *, Konstantina Drosou2, 3, Jeremiah Daniel O’Riordan4, Tengiz Meshveliani5, Ron Pinhasi6, and Robin N. M. Feeney1 1School of Medicine, University College Dublin, Ireland 2Division of Cell Matrix Biology Regenerative Medicine, University of Manchester, United Kingdom 3Manchester Institute of Biotechnology, School of Earth and Environmental Sciences, University of Manchester, United Kingdom [email protected] 5Institute of Paleobiology and Paleoanthropology, National Museum of Georgia, Tbilisi, Georgia 6Department of Evolutionary Anthropology, University of Vienna, Austria *Corresponding Author Abstract The emerging field of ancient metagenomics adds to these Bioinformatic pipelines optimised for the processing and as- processing complexities with the need for additional steps sessment of metagenomic ancient DNA (aDNA) are needed in the separation and authentication of ancient sequences from modern sequences. Currently, there are few pipelines for studies that do not make use of high yielding DNA cap- available for the analysis of ancient metagenomic DNA ture techniques. These bioinformatic pipelines are tradition- 1 4 ally optimised for broad aDNA purposes, are contingent on (aDNA) ≠ The limited number of bioinformatic pipelines selection biases and are associated with high costs.
    [Show full text]
  • Sequence Alignment/Map Format Specification
    Sequence Alignment/Map Format Specification The SAM/BAM Format Specification Working Group 3 Jun 2021 The master version of this document can be found at https://github.com/samtools/hts-specs. This printing is version 53752fa from that repository, last modified on the date shown above. 1 The SAM Format Specification SAM stands for Sequence Alignment/Map format. It is a TAB-delimited text format consisting of a header section, which is optional, and an alignment section. If present, the header must be prior to the alignments. Header lines start with `@', while alignment lines do not. Each alignment line has 11 mandatory fields for essential alignment information such as mapping position, and variable number of optional fields for flexible or aligner specific information. This specification is for version 1.6 of the SAM and BAM formats. Each SAM and BAMfilemay optionally specify the version being used via the @HD VN tag. For full version history see Appendix B. Unless explicitly specified elsewhere, all fields are encoded using 7-bit US-ASCII 1 in using the POSIX / C locale. Regular expressions listed use the POSIX / IEEE Std 1003.1 extended syntax. 1.1 An example Suppose we have the following alignment with bases in lowercase clipped from the alignment. Read r001/1 and r001/2 constitute a read pair; r003 is a chimeric read; r004 represents a split alignment. Coor 12345678901234 5678901234567890123456789012345 ref AGCATGTTAGATAA**GATAGCTGTGCTAGTAGGCAGTCAGCGCCAT +r001/1 TTAGATAAAGGATA*CTG +r002 aaaAGATAA*GGATA +r003 gcctaAGCTAA +r004 ATAGCT..............TCAGC -r003 ttagctTAGGC -r001/2 CAGCGGCAT The corresponding SAM format is:2 1Charset ANSI X3.4-1968 as defined in RFC1345.
    [Show full text]
  • Assembly Exercise
    Assembly Exercise Turning reads into genomes Where we are • 13:30-14:00 – Primer Design to Amplify Microbial Genomes for Sequencing • 14:00-14:15 – Primer Design Exercise • 14:15-14:45 – Molecular Barcoding to Allow Multiplexed NGS • 14:45-15:15 – Processing NGS Data – de novo and mapping assembly • 15:15-15:30 – Break • 15:30-15:45 – Assembly Exercise • 15:45-16:15 – Annotation • 16:15-16:30 – Annotation Exercise • 16:30-17:00 – Submitting Data to GenBank Log onto ILRI cluster • Log in to HPC using ILRI instructions • NOTE: All the commands here are also in the file - assembly_hands_on_steps.txt • If you are like me, it may be easier to cut and paste Linux commands from this file instead of typing them in from the slides Start an interactive session on larger servers • The interactive command will start a session on a server better equipped to do genome assembly $ interactive • Switch to csh (I use some csh features) $ csh • Set up Newbler software that will be used $ module load 454 A norovirus sample sequenced on both 454 and Illumina • The vendors use different file formats unknown_norovirus_454.GACT.sff unknown_norovirus_illumina.fastq • I have converted these files to additional formats for use with the assembly tools unknown_norovirus_454_convert.fasta unknown_norovirus_454_convert.fastq unknown_norovirus_illumina_convert.fasta Set up and run the Newbler de novo assembler • Create a new de novo assembly project $ newAssembly de_novo_assembly • Add read data to the project $ addRun de_novo_assembly unknown_norovirus_454.GACT.sff
    [Show full text]
  • Supplemental Material Nanopore Sequencing of Complex Genomic
    Supplemental Material Nanopore sequencing of complex genomic rearrangements in yeast reveals mechanisms of repeat- mediated double-strand break repair McGinty, et al. Contents: Supplemental Methods Supplemental Figures S1-S6 Supplemental References Supplemental Methods DNA extraction: The DNA extraction protocol used is a slight modification of a classic ethanol precipitation. First, yeast cultures were grown overnight in 2 ml of complete media (YPD) and refreshed for 4 hours in an additional 8 ml YPD to achieve logarithmic growth. The cultures were spun down, and the pellets were resuspended in 290 µl of solution containing 0.9M Sorbitol and 0.1M EDTA at pH 7.5. 10 µl lyticase enzyme was added, and the mixture was incubated for 30 minutes at 37oC in order to break down the yeast cell wall. The mixture was centrifuged at 8,000 rpm for two minutes, and the pellet was resuspended in 270 µl of solution containing 50mM Tris 20mM EDTA at pH 7.5, and 30 µl of 10% SDS. Following five minutes incubation at room temperature, 150 µl of chilled 5M potassium acetate solution was added. This mixture was incubated for 10 min at 4oC and centrifuged for 10 min at 13,000 rpm. The supernatant was then combined with 900 µl of pure ethanol, which had been chilled on ice. This solution was stored overnight at -20oC, and then centrifuged for 20 minutes at 4,000 rpm in a refrigerated centrifuge. The pellet was then washed twice with 70% ethanol, allowed to dry completely, and then resuspended in 50 µl TE (10mM Tris, 1mM EDTA, pH 8).
    [Show full text]
  • Providing Bioinformatics Services on Cloud
    Providing Bioinformatics Services on Cloud Christophe Blanchet, Clément Gauthey InfrastructureC. Blanchet Distributed and C. Gauthey for Biology IDB-IBCPEGI CF13, CNRS Manchester, FR3302 - LYON 9 April - FRANCE2013 http://idee-b.ibcp.fr Infrastructure Distributed for Biology - IDB CNRS-IBCP FR3302, Lyon, FRANCE IDB acknowledges co-funding by the European Community's Seventh Framework Programme (INFSO-RI-261552) and the French National Research Agency's Arpege Programme (ANR-10-SEGI-001) Bioinformatics Today • Biological data are big data • 1512 online databases (NAR Database Issue 2013) • Institut Sanger, UK, 5 PB • Beijing Genome Institute, China, 4 sites, 10 PB ➡ Big data in lot of places • Analysing such data became difficult • Scale-up of the analyses : gene/protein to complete genome/ proteome, ... • Lot of different daily-used tools • That need to be combined in workflows • Usual interfaces: portals, Web services, federation,... ➡ Datacenters with ease of access/use • Distributed resources ADN Experimental platforms: NGS, imaging, ... • ADN • Bioinformatics platforms BI Federation of datacenters M ➡ BI CC ADN ADN BI ADN EGI CF13, Manchester, 9 April 2013 ADN ADN BI ADN ADN M M Sequencing Genomes Complete genome sequencing become a lab commodity with NGS (cheap and efficient) source: www.genomesonline.org source: www.politigenomics.com/next-generation-sequencing-informatics EGI CF13, Manchester, 9 April 2013 Infrastructures in Biology Lot of tools and web services to treat and vizualize lot of data EGI CF13, Manchester, 9 April 2013
    [Show full text]
  • Syntax Highlighting for Computational Biology Artem Babaian1†, Anicet Ebou2, Alyssa Fegen3, Ho Yin (Jeffrey) Kam4, German E
    bioRxiv preprint doi: https://doi.org/10.1101/235820; this version posted December 20, 2017. The copyright holder has placed this preprint (which was not certified by peer review) in the Public Domain. It is no longer restricted by copyright. Anyone can legally share, reuse, remix, or adapt this material for any purpose without crediting the original authors. bioSyntax: Syntax Highlighting For Computational Biology Artem Babaian1†, Anicet Ebou2, Alyssa Fegen3, Ho Yin (Jeffrey) Kam4, German E. Novakovsky5, and Jasper Wong6. 5 10 15 Affiliations: 1. Terry Fox Laboratory, BC Cancer, Vancouver, BC, Canada. [[email protected]] 2. Departement de Formation et de Recherches Agriculture et Ressources Animales, Institut National Polytechnique Felix Houphouet-Boigny, Yamoussoukro, Côte d’Ivoire. [[email protected]] 20 3. Faculty of Science, University of British Columbia, Vancouver, BC, Canada [[email protected]] 4. Faculty of Mathematics, University of Waterloo, Waterloo, ON, Canada. [[email protected]] 5. Department of Medical Genetics, University of British Columbia, Vancouver, BC, 25 Canada. [[email protected]] 6. Genome Science and Technology, University of British Columbia, Vancouver, BC, Canada. [[email protected]] Correspondence†: 30 Artem Babaian Terry Fox Laboratory BC Cancer Research Centre 675 West 10th Avenue Vancouver, BC, Canada. V5Z 1L3. 35 Email: [[email protected]] bioRxiv preprint doi: https://doi.org/10.1101/235820; this version posted December 20, 2017. The copyright holder has placed this preprint (which was not certified by peer review) in the Public Domain. It is no longer restricted by copyright. Anyone can legally share, reuse, remix, or adapt this material for any purpose without crediting the original authors.
    [Show full text]
  • Software List for Biology, Bioinformatics and Biostatistics CCT
    Software List for biology, bioinformatics and biostatistics v CCT - Delta Software Version Application short read assembler and it works on both small and large (mammalian size) ALLPATHS-LG 52488 genomes provides a fast, flexible C++ API & toolkit for reading, writing, and manipulating BAMtools 2.4.0 BAM files a high level of alignment fidelity and is comparable to other mainstream Barracuda 0.7.107b alignment programs allows one to intersect, merge, count, complement, and shuffle genomic bedtools 2.25.0 intervals from multiple files Bfast 0.7.0a universal DNA sequence aligner tool analysis and comprehension of high-throughput genomic data using the R Bioconductor 3.2 statistical programming BioPython 1.66 tools for biological computation written in Python a fast approach to detecting gene-gene interactions in genome-wide case- Boost 1.54.0 control studies short read aligner geared toward quickly aligning large sets of short DNA Bowtie 1.1.2 sequences to large genomes Bowtie2 2.2.6 Bowtie + fully supports gapped alignment with affine gap penalties BWA 0.7.12 mapping low-divergent sequences against a large reference genome ClustalW 2.1 multiple sequence alignment program to align DNA and protein sequences assembles transcripts, estimates their abundances for differential expression Cufflinks 2.2.1 and regulation in RNA-Seq samples EBSEQ (R) 1.10.0 identifying genes and isoforms differentially expressed EMBOSS 6.5.7 a comprehensive set of sequence analysis programs FASTA 36.3.8b a DNA and protein sequence alignment software package FastQC
    [Show full text]
  • Computational Protocol for Assembly and Analysis of SARS-Ncov-2 Genomes
    PROTOCOL Computational Protocol for Assembly and Analysis of SARS-nCoV-2 Genomes Mukta Poojary1,2, Anantharaman Shantaraman1, Bani Jolly1,2 and Vinod Scaria1,2 1CSIR Institute of Genomics and Integrative Biology, Mathura Road, Delhi 2Academy of Scientific and Innovative Research (AcSIR) *Corresponding email: MP [email protected]; AS [email protected]; BJ [email protected]; VS [email protected] ABSTRACT SARS-CoV-2, the pathogen responsible for the ongoing Coronavirus Disease 2019 pandemic is a novel human-infecting strain of Betacoronavirus. The outbreak that initially emerged in Wuhan, China, rapidly spread to several countries at an alarming rate leading to severe global socio-economic disruption and thus overloading the healthcare systems. Owing to the high rate of infection of the virus, as well as the absence of vaccines or antivirals, there is a lack of robust mechanisms to control the outbreak and contain its transmission. Rapid advancement and plummeting costs of high throughput sequencing technologies has enabled sequencing of the virus in several affected individuals globally. Deciphering the viral genome has the potential to help understand the epidemiology of the disease as well as aid in the development of robust diagnostics, novel treatments and prevention strategies. Towards this effort, we have compiled a comprehensive protocol for analysis and interpretation of the sequencing data of SARS-CoV-2 using easy-to-use open source utilities. In this protocol, we have incorporated strategies to assemble the genome of SARS-CoV-2 using two approaches: reference-guided and de novo. Strategies to understand the diversity of the local strain as compared to other global strains have also been described in this protocol.
    [Show full text]
  • Alignment-Free Sequence Comparison: Benefits, Applications, and Tools Andrzej Zielezinski1, Susana Vinga2, Jonas Almeida3 and Wojciech M
    Zielezinski et al. Genome Biology (2017) 18:186 DOI 10.1186/s13059-017-1319-7 REVIEW Open Access Alignment-free sequence comparison: benefits, applications, and tools Andrzej Zielezinski1, Susana Vinga2, Jonas Almeida3 and Wojciech M. Karlowski1* Abstract Alignment-free sequence analyses have been applied to problems ranging from whole-genome phylogeny to the classification of protein families, identification of horizontally transferred genes, and detection of recombined sequences. The strength of these methods makes them particularly useful for next-generation sequencing data processing and analysis. However, many researchers are unclear about how these methods work, how they compare to alignment-based methods, and what their potential is for use for their research. We address these questions and provide a guide to the currently available alignment-free sequence analysis tools. Introduction most alignment programs also model inserted/deleted The 1980s and 1990s were a flourishing time not only states (gaps). However, as our understanding of complex for pop music but also for bioinformatics, where the evolutionary scenarios and our knowledge about the pat- emergence of sequence comparison algorithms revolu- terns and properties of biological sequences advanced, tionized the computational and molecular biology fields. we gradually uncovered some downsides of sequence At that time, many computational biologists quickly be- comparisons based solely on alignments. came stars in the field by developing programs for se- quence alignment, which is a method that positions the Five cases where alignment-based sequence ana- biological sequences’ building blocks to identify regions lysis might be troublesome of similarity that may have consequences for functional, First, alignment-producing programs assume that hom- structural, or evolutionary relationships.
    [Show full text]