scholarly journals Construction of whole genomes from scaffolds using single cell strand-seq data

2018 ◽  
Author(s):  
Mark Hills ◽  
Ester Falconer ◽  
Kieran O’Neil ◽  
Ashley D. Sanders ◽  
Kerstin Howe ◽  
...  
Keyword(s):  
Single Cell ◽  
Sequence Data ◽  
Model Organisms ◽  
Tasmanian Devil ◽  
Template Strand ◽  
Modern Molecular Biology ◽  
Whole Genomes ◽  
Dna Strand ◽  
Genome Assemblies

Accurate reference genome sequences provide the foundation for modern molecular biology and genomics as the interpretation of sequence data to study evolution, gene expression and epigenetics depends heavily on the quality of the genome assembly used for its alignment. Correctly organising sequenced fragments such as contigs and scaffolds in relation to each other is a critical and often challenging step in the construction of robust genome references. We previously identified misoriented regions in the mouse and human reference assemblies using Strand-seq, a single cell sequencing technique that preserves DNA directionality1, 2. Here we demonstrate the ability of Strand-seq to build and correct full-length chromosomes, by identifying which scaffolds belong to the same chromosome and determining their correct order and orientation, without the need for overlapping sequences. We demonstrate that Strand-seq exquisitely maps assembly fragments into large related groups and chromosome-sized clusters without using new assembly data. Using template strand inheritance as a bi-allelic marker, we employ genetic mapping principles to cluster scaffolds that are derived from the same chromosome and order them within the chromosome based solely on directionality of DNA strand inheritance. We prove the utility of our approach by generating improved genome assemblies for several model organisms including the ferret, pig, Xenopus, zebrafish, Tasmanian devil and the Guinea pig.

scholarly journals Construction of Whole Genomes from Scaffolds Using Single Cell Strand-Seq Data

2021 ◽  
Vol 22 (7) ◽  
pp. 3617
Author(s):  
Mark Hills ◽  
Ester Falconer ◽  
Kieran O’Neill ◽  
Ashley D. Sanders ◽  
Kerstin Howe ◽  
...  
Keyword(s):  
Single Cell ◽  
Sequence Data ◽  
Model Organisms ◽  
Tasmanian Devil ◽  
Template Strand ◽  
Modern Molecular Biology ◽  
Whole Genomes ◽  
Dna Strand ◽  
Genome Assemblies

Accurate reference genome sequences provide the foundation for modern molecular biology and genomics as the interpretation of sequence data to study evolution, gene expression, and epigenetics depends heavily on the quality of the genome assembly used for its alignment. Correctly organising sequenced fragments such as contigs and scaffolds in relation to each other is a critical and often challenging step in the construction of robust genome references. We previously identified misoriented regions in the mouse and human reference assemblies using Strand-seq, a single cell sequencing technique that preserves DNA directionality Here we demonstrate the ability of Strand-seq to build and correct full-length chromosomes by identifying which scaffolds belong to the same chromosome and determining their correct order and orientation, without the need for overlapping sequences. We demonstrate that Strand-seq exquisitely maps assembly fragments into large related groups and chromosome-sized clusters without using new assembly data. Using template strand inheritance as a bi-allelic marker, we employ genetic mapping principles to cluster scaffolds that are derived from the same chromosome and order them within the chromosome based solely on directionality of DNA strand inheritance. We prove the utility of our approach by generating improved genome assemblies for several model organisms including the ferret, pig, Xenopus, zebrafish, Tasmanian devil and the Guinea pig.

scholarly journals VGEA: an RNA viral assembly toolkit

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12129
Author(s):  
Paul E. Oluniyi ◽  
Fehintola Ajogbasile ◽  
Judith Oguzie ◽  
Jessica Uwanibe ◽  
Adeyemi Kayode ◽  
...  
Keyword(s):  
De Novo ◽  
Sequence Data ◽  
Workflow Management ◽  
Viral Population ◽  
Lassa Virus ◽  
Viral Genomes ◽  
Bioinformatics Tools ◽  
Reference Sequences ◽  
Genome Assemblies

Next generation sequencing (NGS)-based studies have vastly increased our understanding of viral diversity. Viral sequence data obtained from NGS experiments are a rich source of information, these data can be used to study their epidemiology, evolution, transmission patterns, and can also inform drug and vaccine design. Viral genomes, however, represent a great challenge to bioinformatics due to their high mutation rate and forming quasispecies in the same infected host, bringing about the need to implement advanced bioinformatics tools to assemble consensus genomes well-representative of the viral population circulating in individual patients. Many tools have been developed to preprocess sequencing reads, carry-out de novo or reference-assisted assembly of viral genomes and assess the quality of the genomes obtained. Most of these tools however exist as standalone workflows and usually require huge computational resources. Here we present (Viral Genomes Easily Analyzed), a Snakemake workflow for analyzing RNA viral genomes. VGEA enables users to map sequencing reads to the human genome to remove human contaminants, split bam files into forward and reverse reads, carry out de novo assembly of forward and reverse reads to generate contigs, pre-process reads for quality and contamination, map reads to a reference tailored to the sample using corrected contigs supplemented by the user’s choice of reference sequences and evaluate/compare genome assemblies. We designed a project with the aim of creating a flexible, easy-to-use and all-in-one pipeline from existing/stand-alone bioinformatics tools for viral genome analysis that can be deployed on a personal computer. VGEA was built on the Snakemake workflow management system and utilizes existing tools for each step: fastp (Chen et al., 2018) for read trimming and read-level quality control, BWA (Li & Durbin, 2009) for mapping sequencing reads to the human reference genome, SAMtools (Li et al., 2009) for extracting unmapped reads and also for splitting bam files into fastq files, IVA (Hunt et al., 2015) for de novo assembly to generate contigs, shiver (Wymant et al., 2018) to pre-process reads for quality and contamination, then map to a reference tailored to the sample using corrected contigs supplemented with the user’s choice of existing reference sequences, SeqKit (Shen et al., 2016) for cleaning shiver assembly for QUAST, QUAST (Gurevich et al., 2013) to evaluate/assess the quality of genome assemblies and MultiQC (Ewels et al., 2016) for aggregation of the results from fastp, BWA and QUAST. Our pipeline was successfully tested and validated with SARS-CoV-2 (n = 20), HIV-1 (n = 20) and Lassa Virus (n = 20) datasets all of which have been made publicly available. VGEA is freely available on GitHub at: https://github.com/pauloluniyi/VGEA under the GNU General Public License.

scholarly journals From single nuclei to whole genome assemblies

2019 ◽  
Author(s):  
Merce Montoliu-Nerin ◽  
Marisol Sánchez-García ◽  
Claudia Bergin ◽  
Manfred Grabherr ◽  
Barbara Ellis ◽  
...  
Keyword(s):  
Single Cell ◽  
Large Scale ◽  
Genomic Data ◽  
Life Cycles ◽  
Genomic Research ◽  
Metagenomic Data ◽  
Model Organisms ◽  
Genomic Study ◽  
And Function ◽  
Genome Assemblies

SummaryA large proportion of Earth's biodiversity constitutes organisms that cannot be cultured, have cryptic life-cycles and/or live submerged within their substrates1–4. Genomic data are key to unravel both their identity and function5. The development of metagenomic methods6,7 and the advent of single cell sequencing8–10 have revolutionized the study of life and function of cryptic organisms by upending the need for large and pure biological material, and allowing generation of genomic data from complex or limited environmental samples. Genome assemblies from metagenomic data have so far been restricted to organisms with small genomes, such as bacteria11, archaea12 and certain eukaryotes13. On the other hand, single cell technologies have allowed the targeting of unicellular organisms, attaining a better resolution than metagenomics8,9,14–16, moreover, it has allowed the genomic study of cells from complex organisms one cell at a time17,18. However, single cell genomics are not easily applied to multicellular organisms formed by consortia of diverse taxa, and the generation of specific workflows for sequencing and data analysis is needed to expand genomic research to the entire tree of life, including sponges19, lichens3,20, intracellular parasites21,22, and plant endophytes23,24. Among the most important plant endophytes are the obligate mutualistic symbionts, arbuscular mycorrhizal (AM) fungi, that pose an additional challenge with their multinucleate coenocytic mycelia25. Here, the development of a novel single nuclei sequencing and assembly workflow is reported. This workflow allows, for the first time, the generation of reference genome assemblies from large scale, unbiased sorted, and sequenced AM fungal nuclei circumventing tedious, and often impossible, culturing efforts. This method opens infinite possibilities for studies of evolution and adaptation in these important plant symbionts and demonstrates that reference genomes can be generated from complex non-model organisms by isolating only a handful of their nuclei.

scholarly journals qc3C: reference-free quality control for Hi-C sequencing data

2021 ◽  
Author(s):  
Matthew Z. DeMaere ◽  
Aaron E. Darling
Keyword(s):  
High Throughput Sequencing ◽  
Software Tool ◽  
Model Organisms ◽  
Sequencing Data ◽  
Proximity Ligation ◽  
Genome Wide ◽  
Costly Failure ◽  
Wet Lab ◽  
Genome Assemblies

AbstractHi-C is a sample preparation method that enables high-throughput sequencing to capture genome-wide spatial interactions between DNA molecules. The technique has been successfully applied to solve challenging problems such as 3D structural analysis of chromatin, scaffolding of large genome assemblies and more recently the accurate resolution of metagenome-assembled genomes (MAGs). Despite continued refinements, however, Hi-C library preparation remains a complex laboratory protocol and diligent quality management is recommended to avoid costly failure. Current wet-lab protocols for Hi-C library QC provide only a crude assay, while commonly used sequence-based QC methods demand a reference genome; the quality of which can skew results. We propose a new, reference-free approach for Hi-C library quality assessment that requires only a modest amount of sequencing data. The algorithm builds upon the observation that proximity ligation events are likely to create k -mers that would not naturally occur in the sample. Our software tool (qc3C) is to our knowledge the first to implement a reference-free Hi-C QC tool, and also provides reference-based QC, enabling Hi-C to be more easily applied to non-model organisms and environmental samples. We characterise the accuracy of the new algorithm on simulated and real datasets and compare it to reference-based methods.

scholarly journals mtDNAcombine: tools to combine sequences from multiple studies

2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Eleanor F. Miller ◽  
Andrea Manica
Keyword(s):  
Sequence Data ◽  
Data Extraction ◽  
Bayesian Skyline Plot ◽  
Model Organisms ◽  
Data Sets ◽  
Data Handling ◽  
Online Database ◽  
Genetic Studies ◽  
Wide Range ◽  
Existing Data

Abstract Background Today an unprecedented amount of genetic sequence data is stored in publicly available repositories. For decades now, mitochondrial DNA (mtDNA) has been the workhorse of genetic studies, and as a result, there is a large volume of mtDNA data available in these repositories for a wide range of species. Indeed, whilst whole genome sequencing is an exciting prospect for the future, for most non-model organisms’ classical markers such as mtDNA remain widely used. By compiling existing data from multiple original studies, it is possible to build powerful new datasets capable of exploring many questions in ecology, evolution and conservation biology. One key question that these data can help inform is what happened in a species’ demographic past. However, compiling data in this manner is not trivial, there are many complexities associated with data extraction, data quality and data handling. Results Here we present the mtDNAcombine package, a collection of tools developed to manage some of the major decisions associated with handling multi-study sequence data with a particular focus on preparing sequence data for Bayesian skyline plot demographic reconstructions. Conclusions There is now more genetic information available than ever before and large meta-data sets offer great opportunities to explore new and exciting avenues of research. However, compiling multi-study datasets still remains a technically challenging prospect. The mtDNAcombine package provides a pipeline to streamline the process of downloading, curating, and analysing sequence data, guiding the process of compiling data sets from the online database GenBank.

scholarly journals Hawkeye and AMOS: visualizing and assessing the quality of genome assemblies

2011 ◽  
Vol 14 (2) ◽  
pp. 213-224 ◽  
Author(s):  
M. C. Schatz ◽  
A. M. Phillippy ◽  
D. D. Sommer ◽  
A. L. Delcher ◽  
D. Puiu ◽  
...  
Keyword(s):  
Genome Assemblies

scholarly journals A new Plasmodium vivax reference sequence with improved assembly of the subtelomeres reveals an abundance of pir genes

2016 ◽  
Vol 1 ◽  
pp. 4 ◽  
Author(s):  
Sarah Auburn ◽  
Ulrike Böhme ◽  
Sascha Steinbiss ◽  
Hidayat Trimarsanto ◽  
Jessica Hostetler ◽  
...  
Keyword(s):  
South America ◽  
Plasmodium Vivax ◽  
Sequence Data ◽  
Ex Vivo ◽  
Vital Role ◽  
Asia Pacific ◽  
Reference Sequence ◽  
Manual Curation ◽  
Illumina Sequence

Plasmodium vivax is now the predominant cause of malaria in the Asia-Pacific, South America and Horn of Africa. Laboratory studies of this species are constrained by the inability to maintain the parasite in continuous ex vivo culture, but genomic approaches provide an alternative and complementary avenue to investigate the parasite’s biology and epidemiology. To date, molecular studies of P. vivax have relied on the Salvador-I reference genome sequence, derived from a monkey-adapted strain from South America. However, the Salvador-I reference remains highly fragmented with over 2500 unassembled scaffolds.  Using high-depth Illumina sequence data, we assembled and annotated a new reference sequence, PvP01, sourced directly from a patient from Papua Indonesia. Draft assemblies of isolates from China (PvC01) and Thailand (PvT01) were also prepared for comparative purposes. The quality of the PvP01 assembly is improved greatly over Salvador-I, with fragmentation reduced to 226 scaffolds. Detailed manual curation has ensured highly comprehensive annotation, with functions attributed to 58% core genes in PvP01 versus 38% in Salvador-I. The assemblies of PvP01, PvC01 and PvT01 are larger than that of Salvador-I (28-30 versus 27 Mb), owing to improved assembly of the subtelomeres.  An extensive repertoire of over 1200 Plasmodium interspersed repeat (pir) genes were identified in PvP01 compared to 346 in Salvador-I, suggesting a vital role in parasite survival or development. The manually curated PvP01 reference and PvC01 and PvT01 draft assemblies are important new resources to study vivax malaria. PvP01 is maintained at GeneDB and ongoing curation will ensure continual improvements in assembly and annotation quality.

scholarly journals The universally-conserved transcription factor RfaH is recruited to a hairpin structure of the non-template DNA strand

eLife ◽  
2018 ◽  
Vol 7 ◽  
Author(s):  
Philipp K Zuber ◽  
Irina Artsimovitch ◽  
Monali NandyMazumdar ◽  
Zhaokun Liu ◽  
Yuri Nedialkov ◽  
...  
Keyword(s):  
Virulence Genes ◽  
Complex Structure ◽  
Transcription Regulator ◽  
Transcription Bubble ◽  
Template Strand ◽  
Domains Of Life ◽  
Dna Strand ◽  
Functional Analyses ◽  
Template Dna

RfaH, a transcription regulator of the universally conserved NusG/Spt5 family, utilizes a unique mode of recruitment to elongating RNA polymerase to activate virulence genes. RfaH function depends critically on an ops sequence, an exemplar of a consensus pause, in the non-template DNA strand of the transcription bubble. We used structural and functional analyses to elucidate the role of ops in RfaH recruitment. Our results demonstrate that ops induces pausing to facilitate RfaH binding and establishes direct contacts with RfaH. Strikingly, the non-template DNA forms a hairpin in the RfaH:ops complex structure, flipping out a conserved T residue that is specifically recognized by RfaH. Molecular modeling and genetic evidence support the notion that ops hairpin is required for RfaH recruitment. We argue that both the sequence and the structure of the non-template strand are read out by transcription factors, expanding the repertoire of transcriptional regulators in all domains of life.

scholarly journals Apoptotic endonuclease EndoG regulates alternative splicing of human telomerase catalytic subunit hTERT

2016 ◽  
Vol 62 (5) ◽  
pp. 544-554 ◽  
Author(s):  
D.D. Zhdanov ◽  
D.A. Vasina ◽  
E.V. Orlova ◽  
V.S. Orlova ◽  
M.V. Pokrovskaya ◽  
...  
Keyword(s):  
Alternative Splicing ◽  
Telomerase Activity ◽  
Catalytic Subunit ◽  
Over Expression ◽  
Htert Gene ◽  
Template Strand ◽  
Non Coding Rna ◽  
Dna Strand ◽  
Human Telomerase ◽  
Long Non Coding Rna

Human telomerase catalytic subunit hTERT is subjected to alternative splicing results in loss of its function and leads to decrease of telomerase activity. However, very little is known about the mechanism of hTERT pre-mRNA alternative splicing. Apoptotic endonuclease EndoG is known to participate this process. The aim of this study was to determine the role of EndoG in regulation of hTERT alternative splicing. Increased expression of b-deletion splice variant was determined during EndoG over-expression in CaCo-2 cell line, after EndoG treatment of cell cytoplasm and nuclei and after nuclei incubation with EndoG digested cell RNA. hTERT alternative splicing was induced by 47-mer RNA oligonucleotide in naked nuclei and in cells after transfection. Identified long non-coding RNA, that is the precursor of 47-mer RNA oligonucleotide. Its size is 1754 nucleotides. Based on the results the following mechanism was proposed. hTERT pre-mRNA is transcribed from coding DNA strand while long non-coding RNA is transcribed from template strand of hTERT gene. EndoG digests long non-coding RNA and produces 47-mer RNA oligonucleotide complementary to hTERT pre-mRNA exon 8 and intron 8 junction place. Interaction of 47-mer RNA oligonucleotide and hTERT pre-mRNA causes alternative splicing.

Sign in / Sign up

Export Citation Format

Share Document