More Like this

1. Estimating error models for whole genome sequencing using mixtures of Dirichlet-multinomial distributions

   https://doi.org/10.1093/bioinformatics/btx133  

   Wu, Steven H. ; Schwartz, Rachel S. ; Winter, David J. ; Conrad, Donald F. ; Cartwright, Reed A. ; Stegle, ed., Oliver ( March 2017 , Bioinformatics)

   Accurate identification of genotypes is an essential part of the analysis of genomic data, including in identification of sequence polymorphisms, linking mutations with disease and determining mutation rates. Biological and technical processes that adversely affect genotyping include copy-number-variation, paralogous sequences, library preparation, sequencing error and reference-mapping biases, among others.

   We modeled the read depth for all data as a mixture of Dirichlet-multinomial distributions, resulting in significant improvements over previously used models. In most cases the best model was comprised of two distributions. The major-component distribution is similar to a binomial distribution with low error and low reference bias. The minor-component distribution is overdispersed with higher error and reference bias. We also found that sites fitting the minor component are enriched for copy number variants and low complexity regions, which can produce erroneous genotype calls. By removing sites that do not fit the major component, we can improve the accuracy of genotype calls.

   Methods and data files are available at https://github.com/CartwrightLab/WuEtAl2017/ (doi:10.5281/zenodo.256858).

   Supplementary data is available at Bioinformatics online.

    

   more » « less
2. lordFAST: sensitive and Fast Alignment Search Tool for LOng noisy Read sequencing Data

   https://doi.org/10.1093/bioinformatics/bty544  

   Haghshenas, Ehsan ; Sahinalp, S. Cenk ; Hach, Faraz ; Berger, ed., Bonnie ( July 2018 , Bioinformatics)

   Recent advances in genomics and precision medicine have been made possible through the application of high throughput sequencing (HTS) to large collections of human genomes. Although HTS technologies have proven their use in cataloging human genome variation, computational analysis of the data they generate is still far from being perfect. The main limitation of Illumina and other popular sequencing technologies is their short read length relative to the lengths of (common) genomic repeats. Newer (single molecule sequencing – SMS) technologies such as Pacific Biosciences and Oxford Nanopore are producing longer reads, making it theoretically possible to overcome the difficulties imposed by repeat regions. Unfortunately, because of their high sequencing error rate, reads generated by these technologies are very difficult to work with and cannot be used in many of the standard downstream analysis pipelines. Note that it is not only difficult to find the correct mapping locations of such reads in a reference genome, but also to establish their correct alignment so as to differentiate sequencing errors from real genomic variants. Furthermore, especially since newer SMS instruments provide higher throughput, mapping and alignment need to be performed much faster than before, maintaining high sensitivity.

   We introduce lordFAST, a novel long-read mapper that is specifically designed to align reads generated by PacBio and potentially other SMS technologies to a reference. lordFAST not only has higher sensitivity than the available alternatives, it is also among the fastest and has a very low memory footprint.

   lordFAST is implemented in C++ and supports multi-threading. The source code of lordFAST is available at https://github.com/vpc-ccg/lordfast.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
3. MAGUS+eHMMs: improved multiple sequence alignment accuracy for fragmentary sequences

   https://doi.org/10.1093/bioinformatics/btab788  

   Shen, Chengze ; Zaharias, Paul ; Warnow, Tandy ; Boeva, ed., Valentina ( November 2021 , Bioinformatics)

   Multiple sequence alignment is an initial step in many bioinformatics pipelines, including phylogeny estimation, protein structure prediction and taxonomic identification of reads produced in amplicon or metagenomic datasets, etc. Yet, alignment estimation is challenging on datasets that exhibit substantial sequence length heterogeneity, and especially when the datasets have fragmentary sequences as a result of including reads or contigs generated by next-generation sequencing technologies. Here, we examine techniques that have been developed to improve alignment estimation when datasets contain substantial numbers of fragmentary sequences. We find that MAGUS, a recently developed MSA method, is fairly robust to fragmentary sequences under many conditions, and that using a two-stage approach where MAGUS is used to align selected ‘backbone sequences’ and the remaining sequences are added into the alignment using ensembles of Hidden Markov Models further improves alignment accuracy. The combination of MAGUS with the ensemble of eHMMs (i.e. MAGUS+eHMMs) clearly improves on UPP, the previous leading method for aligning datasets with high levels of fragmentation.

   UPP is available on https://github.com/smirarab/sepp, and MAGUS is available on https://github.com/vlasmirnov/MAGUS. MAGUS+eHMMs can be performed by running MAGUS to obtain the backbone alignment, and then using the backbone alignment as an input to UPP.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
4. TopHap: rapid inference of key phylogenetic structures from common haplotypes in large genome collections with limited diversity

   https://doi.org/10.1093/bioinformatics/btac186  

   Caraballo-Ortiz, Marcos A. ; Miura, Sayaka ; Sanderford, Maxwell ; Dolker, Tenzin ; Tao, Qiqing ; Weaver, Steven ; Pond, Sergei L. K. ; Kumar, Sudhir ; Schwartz, ed., Russell ( March 2022 , Bioinformatics)

   Building reliable phylogenies from very large collections of sequences with a limited number of phylogenetically informative sites is challenging because sequencing errors and recurrent/backward mutations interfere with the phylogenetic signal, confounding true evolutionary relationships. Massive global efforts of sequencing genomes and reconstructing the phylogeny of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains exemplify these difficulties since there are only hundreds of phylogenetically informative sites but millions of genomes. For such datasets, we set out to develop a method for building the phylogenetic tree of genomic haplotypes consisting of positions harboring common variants to improve the signal-to-noise ratio for more accurate and fast phylogenetic inference of resolvable phylogenetic features.

   We present the TopHap approach that determines spatiotemporally common haplotypes of common variants and builds their phylogeny at a fraction of the computational time of traditional methods. We develop a bootstrap strategy that resamples genomes spatiotemporally to assess topological robustness. The application of TopHap to build a phylogeny of 68 057 SARS-CoV-2 genomes (68KG) from the first year of the pandemic produced an evolutionary tree of major SARS-CoV-2 haplotypes. This phylogeny is concordant with the mutation tree inferred using the co-occurrence pattern of mutations and recovers key phylogenetic relationships from more traditional analyses. We also evaluated alternative roots of the SARS-CoV-2 phylogeny and found that the earliest sampled genomes in 2019 likely evolved by four mutations of the most recent common ancestor of all SARS-CoV-2 genomes. An application of TopHap to more than 1 million SARS-CoV-2 genomes reconstructed the most comprehensive evolutionary relationships of major variants, which confirmed the 68KG phylogeny and provided evolutionary origins of major and recent variants of concern.

   TopHap is available at https://github.com/SayakaMiura/TopHap.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
5. ISSRseq: An extensible method for reduced representation sequencing

   https://doi.org/10.1111/2041-210X.13784  

   Sinn, Brandon T. ; Simon, Sandra J. ; Santee, Mathilda V. ; DiFazio, Stephen P. ; Fama, Nicole M. ; Barrett, Craig F. ( December 2021 , Methods in Ecology and Evolution)

   The capability to generate densely sampled single nucleotide polymorphism (SNP) data is essential in diverse subdisciplines of biology, including crop breeding, pathology, forensics, forestry, ecology, evolution and conservation. However, the wet‐laboratory expertise and bioinformatics training required to conduct genome‐scale variant discovery remain limiting factors for investigators with limited resources.

   Here we present ISSRseq, a PCR‐based method for reduced representation of genomic variation using simple sequence repeats as priming sites to sequence inter simple sequence repeat (ISSR) regions. Briefly, ISSR regions are amplified with single primers, pooled, used to construct sequencing libraries with a commercially available kit, and sequenced on the Illumina platform. We also present a flexible bioinformatic pipeline that assembles ISSR loci, calls and hard filters variants, outputs data matrices in common formats, and conducts population analyses using R.

   Using three angiosperm species as case studies, we demonstrate that ISSRseq is highly repeatable, necessitates only simple wet‐laboratory skills and commonplace instrumentation, is flexible in terms of the number of single primers used, and can generate genomic‐scale variant discovery on par with existing RRS methods which require more complex wet‐laboratory procedures.

   ISSRseq represents a straightforward approach to SNP genotyping in any organism, and we predict that this method will be particularly useful for those studying population genomics and phylogeography of non‐model organisms. Furthermore, the ease of ISSRseq relative to other RRS methods should prove useful to those lacking advanced expertise in wet‐laboratory methods or bioinformatics.

    

   more » « less