scholarly journals AGOUTI: improving genome assembly and annotation using transcriptome data

2015 ◽  
Author(s):  
Simo V. Zhang ◽  
Luting Zhuo ◽  
Matthew W. Hahn
Keyword(s):  
Genome Assembly ◽  
Gene Annotation ◽  
Supplementary Information ◽  
Gene Identification ◽  
Supplementary Data ◽  
Transcriptome Data ◽  
Rna Seq ◽  
Separate Gene ◽  
Gene Models ◽  
Genome Assemblies

AbstractSummaryCurrent genome assemblies consist of thousands of contigs. These incomplete and fragmented assemblies lead to errors in gene identification, such that single genes spread across multiple contigs are annotated as separate gene models. We present AGOUTI (Annotated Genome Optimization Using Transcriptome Information), a tool that uses RNA-seq data to simultaneously combine contigs into scaffolds and fragmented gene models into single models. We show that AGOUTI improves both the contiguity of genome assemblies and the accuracy of gene annotation, providing updated versions of each as output.AvailabilityThe software is implemented in python and is available from github.com/svm-zhang/[email protected] informationSupplementary data are available at Bioinformatics online.

scholarly journals ARBitR: an overlap-aware genome assembly scaffolder for linked reads

2020 ◽  
Author(s):  
Markus Hiltunen ◽  
Martin Ryberg ◽  
Hanna Johannesson
Keyword(s):  
Genome Assembly ◽  
General Public ◽  
Source Code ◽  
Draft Genome ◽  
Supplementary Information ◽  
Genomic Sequencing ◽  
Supplementary Data ◽  
Genome Assemblies ◽  
General Public License

Abstract Summary Linked genomic sequencing reads contain information that can be used to join sequences together into scaffolds in draft genome assemblies. Existing software for this purpose performs the scaffolding by joining sequences with a gap between them, not considering potential overlaps of contigs. We developed ARBitR to create scaffolds where overlaps are taken into account and show that it can accurately recreate regions where draft assemblies are broken. Availability and implementation ARBitR is written and implemented in Python3 for Unix-based operative systems. All source code is available at https://github.com/markhilt/ARBitR under the GNU General Public License v3. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals DEsingle for detecting three types of differential expression in single-cell RNA-seq data

2017 ◽  
Author(s):  
Zhun Miao ◽  
Ke Deng ◽  
Xiaowo Wang ◽  
Xuegong Zhang
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Negative Binomial ◽  
Single Cells ◽  
R Package ◽  
Supplementary Information ◽  
Binomial Model ◽  
Supplementary Data ◽  
Rna Seq ◽  
Real Zeros

AbstractSummaryThe excessive amount of zeros in single-cell RNA-seq data include “real” zeros due to the on-off nature of gene transcription in single cells and “dropout” zeros due to technical reasons. Existing differential expression (DE) analysis methods cannot distinguish these two types of zeros. We developed an R package DEsingle which employed Zero-Inflated Negative Binomial model to estimate the proportion of real and dropout zeros and to define and detect 3 types of DE genes in single-cell RNA-seq data with higher accuracy.Availability and ImplementationThe R package DEsingle is freely available at https://github.com/miaozhun/DEsingle and is under Bioconductor’s consideration [email protected] informationSupplementary data are available at bioRxiv online.

scholarly journals SVIM-asm: Structural variant detection from haploid and diploid genome assemblies

2020 ◽  
Author(s):  
David Heller ◽  
Martin Vingron
Keyword(s):  
Genetic Information ◽  
Source Code ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Diploid Genome ◽  
Insertions And Deletions ◽  
Structural Variant ◽  
Sequencing Technologies ◽  
Variant Detection ◽  
Genome Assemblies

AbstractMotivationWith the availability of new sequencing technologies, the generation of haplotype-resolved genome assemblies up to chromosome scale has become feasible. These assemblies capture the complete genetic information of both parental haplotypes, increase structural variant (SV) calling sensitivity and enable direct genotyping and phasing of SVs. Yet, existing SV callers are designed for haploid genome assemblies only, do not support genotyping or detect only a limited set of SV classes.ResultsWe introduce our method SVIM-asm for the detection and genotyping of six common classes of SVs from haploid and diploid genome assemblies. Compared against the only other existing SV caller for diploid assemblies, DipCall, SVIM-asm detects more SV classes and reached higher F1 scores for the detection of insertions and deletions on two recently published assemblies of the HG002 individual.Availability and ImplementationSVIM-asm has been implemented in Python and can be easily installed via bioconda. Its source code is available at github.com/eldariont/[email protected] informationSupplementary data are available online.

scholarly journals LIONS: analysis suite for detecting and quantifying transposable element initiated transcription from RNA-seq

2019 ◽  
Vol 35 (19) ◽  
pp. 3839-3841 ◽  
Author(s):  
Artem Babaian ◽  
I Richard Thompson ◽  
Jake Lever ◽  
Liane Gagnier ◽  
Mohammad M Karimi ◽  
...  
Keyword(s):  
Transposable Elements ◽  
Transposable Element ◽  
Test Data ◽  
Source Code ◽  
Supplementary Information ◽  
Transcriptional Networks ◽  
Supplementary Data ◽  
Rna Seq ◽  
Transcriptional Initiation ◽  
Instruction Manual

Abstract Summary Transposable elements (TEs) influence the evolution of novel transcriptional networks yet the specific and meaningful interpretation of how TE-derived transcriptional initiation contributes to the transcriptome has been marred by computational and methodological deficiencies. We developed LIONS for the analysis of RNA-seq data to specifically detect and quantify TE-initiated transcripts. Availability and implementation Source code, container, test data and instruction manual are freely available at www.github.com/ababaian/LIONS. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals MAJIQ-SPEL: web-tool to interrogate classical and complex splicing variations from RNA-Seq data

2017 ◽  
Vol 34 (2) ◽  
pp. 300-302 ◽  
Author(s):  
Christopher J Green ◽  
Matthew R Gazzara ◽  
Yoseph Barash
Keyword(s):  
Experimental Validation ◽  
Ucsc Genome Browser ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Rna Seq ◽  
Web Tool ◽  
Rt Pcr ◽  
Design Algorithm ◽  
Gene Isoforms ◽  
Downstream Analysis

Abstract Summary Analysis of RNA sequencing (RNA-Seq) data have highlighted the fact that most genes undergo alternative splicing (AS) and that these patterns are tightly regulated. Many of these events are complex, resulting in numerous possible isoforms that quickly become difficult to visualize, interpret and experimentally validate. To address these challenges we developed MAJIQ-SPEL, a web-tool that takes as input local splicing variations (LSVs) quantified from RNA-Seq data and provides users with visualization and quantification of gene isoforms associated with those. Importantly, MAJIQ-SPEL is able to handle both classical (binary) and complex, non-binary, splicing variations. Using a matching primer design algorithm it also suggests to users possible primers for experimental validation by RT-PCR and displays those, along with the matching protein domains affected by the LSV, on UCSC Genome Browser for further downstream analysis. Availability and implementation Program and code will be available athttp://majiq.biociphers.org/majiq-spel. Supplementary information Supplementary data are available atBioinformatics online.

scholarly journals LRez: C ++ API and toolkit for analyzing and managing Linked-Reads data

2021 ◽  
Author(s):  
Pierre Morisse ◽  
Claire Lemaitre ◽  
Fabrice Legeai
Keyword(s):  
Genome Assembly ◽  
Low Cost ◽  
Variant Calling ◽  
Supplementary Information ◽  
Supplementary Data ◽  
High Quality ◽  
Dna Molecule ◽  
Sequencing Technologies ◽  
Wide Range ◽  
Genomic Regions

Abstract Motivation Linked-Reads technologies combine both the high-quality and low cost of short-reads sequencing and long-range information, through the use of barcodes tagging reads which originate from a common long DNA molecule. This technology has been employed in a broad range of applications including genome assembly, phasing and scaffolding, as well as structural variant calling. However, to date, no tool or API dedicated to the manipulation of Linked-Reads data exist. Results We introduce LRez, a C ++ API and toolkit which allows easy management of Linked-Reads data. LRez includes various functionalities, for computing numbers of common barcodes between genomic regions, extracting barcodes from BAM files, as well as indexing and querying BAM, FASTQ and gzipped FASTQ files to quickly fetch all reads or alignments containing a given barcode. LRez is compatible with a wide range of Linked-Reads sequencing technologies, and can thus be used in any tool or pipeline requiring barcode processing or indexing, in order to improve their performances. Availability and implementation LRez is implemented in C ++, supported on Unix-based platforms, and available under AGPL-3.0 License at https://github.com/morispi/LRez, and as a bioconda module. Supplementary information Supplementary data are available at Bioinformatics Advances

scholarly journals Per-sample standardization and asymmetric winsorization lead to accurate clustering of RNA-seq expression profiles

2021 ◽  
Author(s):  
Davide Risso ◽  
Stefano Maria Pagnotta
Keyword(s):  
Single Cell ◽  
Expression Profiles ◽  
Unsupervised Clustering ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Rna Seq ◽  
Data Transformations ◽  
The Impact

Abstract Motivation Data transformations are an important step in the analysis of RNA-seq data. Nonetheless, the impact of transformation on the outcome of unsupervised clustering procedures is still unclear. Results Here, we present an Asymmetric Winsorization per Sample Transformation (AWST), which is robust to data perturbations and removes the need for selecting the most informative genes prior to sample clustering. Our procedure leads to robust and biologically meaningful clusters both in bulk and in single-cell applications. Availability The AWST method is available at https://github.com/drisso/awst. The code to reproduce the analyses is available at https://github.com/drisso/awst\_analysis. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals GPress: a framework for querying general feature format (GFF) files and expression files in a compressed form

2020 ◽  
Vol 36 (18) ◽  
pp. 4810-4812
Author(s):  
Qingxi Meng ◽  
Idoia Ochoa ◽  
Mikel Hernaez
Keyword(s):  
Single Cell ◽  
Data Streams ◽  
General Feature ◽  
Supplementary Information ◽  
Storage Space ◽  
Supplementary Data ◽  
Rna Seq ◽  
Sequencing Data ◽  
General Feature Format ◽  
Original File

Abstract Motivation Sequencing data are often summarized at different annotation levels for further analysis, generally using the general feature format (GFF) or its descendants, gene transfer format (GTF) and GFF3. Existing utilities for accessing these files, like gffutils and gffread, do not focus on reducing the storage space, significantly increasing it in some cases. We propose GPress, a framework for querying GFF files in a compressed form. GPress can also incorporate and compress expression files from both bulk and single-cell RNA-Seq experiments, supporting simultaneous queries on both the GFF and expression files. In brief, GPress applies transformations to the data which are then compressed with the general lossless compressor BSC. To support queries, GPress compresses the data in blocks and creates several index tables for fast retrieval. Results We tested GPress on several GFF files of different organisms, and showed that it achieves on average a 61% reduction in size with respect to gzip (the current de facto compressor for GFF files) while being able to retrieve all annotations for a given identifier or a range of coordinates in a few seconds (when run in a common laptop). In contrast, gffutils provides faster retrieval but doubles the size of the GFF files. When additionally linking an expression file, we show that GPress can reduce its size by more than 68% when compared to gzip (for both bulk and single-cell RNA-Seq experiments), while still retrieving the information within seconds. Finally, applying BSC to the data streams generated by GPress instead of to the original file shows a size reduction of more than 44% on average. Availability and implementation GPress is freely available at https://github.com/qm2/gpress. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals Lost in Translation: The Pitfalls of Ensembl Gene Annotations Between Human Genome Assemblies and Their Impact on Diagnostics

2020 ◽  
Author(s):  
Mohammed O.E Abdallah ◽  
Mahmoud Koko ◽  
Raj Ramesar
Keyword(s):  
Human Genome ◽  
Genome Assembly ◽  
Evolutionary Constraint ◽  
Clinical Genetics ◽  
Ensembl Gene ◽  
Protein Coding ◽  
Gene Annotations ◽  
Human Genome Assembly ◽  
Gene Models ◽  
Genome Assemblies

Abstract Background:The GRCh37 human genome assembly is still widely used in genomics despite the fact an updated human genome assembly (GRCh38) has been available for many years. A particular issue with relevant ramifications for clinical genetics currently is the case of the GRCh37 Ensembl gene annotations which has been archived, and thus not updated, since 2013. These Ensembl GRCh37 gene annotations are just as ubiquitous as the former assembly and are the default gene models used and preferred by the majority of genomic projects internationally. In this study, we highlight the issue of genes with discrepant annotations, that have been recognized as protein coding in the new but not the old assembly. These genes are ignored by all genomic resources that still rely on the archived and outdated gene annotations. Moreover, the majority if not all of these discrepant genes (DGs) are automatically discarded and ignored by all variant prioritization tools that rely on the GRCh37 Ensembl gene annotations.Methods:We performed bioinformatics analysis identifying Ensembl genes with discrepant annotations between the two most recent human genome assemblies, hg37, hg38, respectively. Clinical and phenotype gene curations have been obtained and compared for this gene set. Furthermore, matching RefSeq transcripts have also been collated and analyzed. ٌResults:We found hundreds of genes (N=267) that were reclassified as “protein-coding” in the new hg38 assembly. Notably, 169 of these genes also had a discrepant HGNC gene symbol between the two assemblies.Most genes had RefSeq matches (N=199/267) including all the genes with defined phenotypes in Ensembl genes GRCh38 assembly (N=10). However, many protein-coding genes remain missing from the current known RefSeq gene models (N=68)Conclusion: We found many clinically relevant genes in this group of neglected genes and we anticipate that many more will be found relevant in the future. For these genes, the inaccurate label of “non-protein-coding” hinders the possibility of identifying any causal sequence variants that overlap them. In addition, Important additional annotations such as evolutionary constraint metrics are also not calculated for these genes for the same reason, further relegating them into oblivion.

scholarly journals APAlyzer: a bioinformatics package for analysis of alternative polyadenylation isoforms

2020 ◽  
Vol 36 (12) ◽  
pp. 3907-3909 ◽  
Author(s):  
Ruijia Wang ◽  
Bin Tian
Keyword(s):  
Gene Expression ◽  
Alternative Polyadenylation ◽  
Supplementary Information ◽  
Human Tissues ◽  
Bioconductor Package ◽  
Supplementary Data ◽  
Rna Seq ◽  
Eukaryotic Genes ◽  
Polyadenylation Sites

Abstract Summary Most eukaryotic genes produce alternative polyadenylation (APA) isoforms. APA is dynamically regulated under different growth and differentiation conditions. Here, we present a bioinformatics package, named APAlyzer, for examining 3′UTR APA, intronic APA and gene expression changes using RNA-seq data and annotated polyadenylation sites in the PolyA_DB database. Using APAlyzer and data from the GTEx database, we present APA profiles across human tissues. Availability and implementation APAlyzer is freely available at https://bioconductor.org/packages/release/bioc/html/APAlyzer.html as an R/Bioconductor package. Supplementary information Supplementary data are available at Bioinformatics online.

Sign in / Sign up

Export Citation Format

Share Document