scholarly journals RsQTL: correlation of expressed SNVs with splicing using RNA-sequencing data

2019 ◽  
Author(s):  
Justin Sein ◽  
Liam F. Spurr ◽  
Pavlos Bousounis ◽  
N M Prashant ◽  
Hongyu Liu ◽  
...  
Keyword(s):  
Rna Sequencing ◽  
Tissue Expression ◽  
Supplementary Information ◽  
Rna Seq ◽  
Sequencing Data ◽  
Dynamic Nature ◽  
Dynamic Variation ◽  
Exon Junctions ◽  
Variant Allele Fraction ◽  
Allele Fraction

SummaryRsQTL is a tool for identification of splicing quantitative trait loci (sQTLs) from RNA-sequencing (RNA-seq) data by correlating the variant allele fraction at expressed SNV loci in the transcriptome (VAFRNA) with the proportion of molecules spanning local exon-exon junctions at loci with differential intron excision (percent spliced in, PSI). We exemplify the method on sets of RNA-seq data from human tissues obtained though the Genotype-Tissue Expression Project (GTEx). RsQTL does not require matched DNA and can identify a subset of expressed sQTL loci. Due to the dynamic nature of VAFRNA, RsQTL is applicable for the assessment of conditional and dynamic variation-splicing relationships.Availability and implementationhttps://github.com/HorvathLab/[email protected] or [email protected] InformationRsQTL_Supplementary_Data.zip

scholarly journals ReQTL: identifying correlations between expressed SNVs and gene expression using RNA-sequencing data

2019 ◽  
Vol 36 (5) ◽  
pp. 1351-1359 ◽  
Author(s):  
Liam F Spurr ◽  
Nawaf Alomran ◽  
Pavlos Bousounis ◽  
Dacian Reece-Stremtan ◽  
N M Prashant ◽  
...  
Keyword(s):  
Gene Expression ◽  
Genetic Variation ◽  
Rna Sequencing ◽  
Variant Allele ◽  
Supplementary Information ◽  
Phenotypic Trait ◽  
Continuous Measure ◽  
Rna Seq ◽  
Variant Allele Fraction ◽  
Allele Fraction

Abstract Motivation By testing for associations between DNA genotypes and gene expression levels, expression quantitative trait locus (eQTL) analyses have been instrumental in understanding how thousands of single nucleotide variants (SNVs) may affect gene expression. As compared to DNA genotypes, RNA genetic variation represents a phenotypic trait that reflects the actual allele content of the studied system. RNA genetic variation at expressed SNV loci can be estimated using the proportion of alleles bearing the variant nucleotide (variant allele fraction, VAFRNA). VAFRNA is a continuous measure which allows for precise allele quantitation in loci where the RNA alleles do not scale with the genotype count. We describe a method to correlate VAFRNA with gene expression and assess its ability to identify genetically regulated expression solely from RNA-sequencing (RNA-seq) datasets. Results We introduce ReQTL, an eQTL modification which substitutes the DNA allele count for the variant allele fraction at expressed SNV loci in the transcriptome (VAFRNA). We exemplify the method on sets of RNA-seq data from human tissues obtained though the Genotype-Tissue Expression (GTEx) project and demonstrate that ReQTL analyses are computationally feasible and can identify a subset of expressed eQTL loci. Availability and implementation A toolkit to perform ReQTL analyses is available at https://github.com/HorvathLab/ReQTL. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals SPsimSeq: semi-parametric simulation of bulk and single cell RNA sequencing data

2019 ◽  
Author(s):  
Alemu Takele Assefa ◽  
Jo Vandesompele ◽  
Olivier Thas
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Rna Sequencing ◽  
Empirical Distribution ◽  
Supplementary Information ◽  
Rna Seq ◽  
Sequencing Data ◽  
Actual Distribution ◽  
Wide Range ◽  
Single Cell Rna Sequencing

SummarySPsimSeq is a semi-parametric simulation method for bulk and single cell RNA sequencing data. It simulates data from a good estimate of the actual distribution of a given real RNA-seq dataset. In contrast to existing approaches that assume a particular data distribution, our method constructs an empirical distribution of gene expression data from a given source RNA-seq experiment to faithfully capture the data characteristics of real data. Importantly, our method can be used to simulate a wide range of scenarios, such as single or multiple biological groups, systematic variations (e.g. confounding batch effects), and different sample sizes. It can also be used to simulate different gene expression units resulting from different library preparation protocols, such as read counts or UMI counts.Availability and implementationThe R package and associated documentation is available from https://github.com/CenterForStatistics-UGent/SPsimSeq.Supplementary informationSupplementary data are available at bioRχiv online.

scholarly journals Arnav: Site specific error models to identify variants in RNA

2018 ◽  
Author(s):  
Avinash Ramu ◽  
Donald F. Conrad
Keyword(s):  
Statistical Method ◽  
Rna Sequencing ◽  
High Specificity ◽  
Sequencing Data ◽  
Site Specific ◽  
Error Models ◽  
True Variant ◽  
Specific Error ◽  
Variant Allele Fraction ◽  
Allele Fraction

AbstractSummaryWe present Arnav (Analysis of RNA variants) a lightweight and easy-to-use statistical method for detecting mutations from RNA sequencing data. Site-specific error models allow Arnav to call variants with high specificity when the true variant allele fraction is unknown. We show the utility of Arnav by identifying variants using RNA sequencing data from the GTEx project.Availability and ImplementationArnav is implemented in C++ and is distributed under the GPL license at https://github.com/gatoravi/[email protected]

scholarly journals Rigorous Benchmarking of HLA Callers for RNA Sequencing Data

2021 ◽  
Author(s):  
Ram Ayyala ◽  
Junghyun Jung ◽  
Sergey Knyazev ◽  
SERGHEI MANGUL
Keyword(s):  
Rna Sequencing ◽  
Gold Standard ◽  
Large Scale ◽  
Human Leukocyte ◽  
Tissue Expression ◽  
Hla Typing ◽  
Read Length ◽  
Evaluation Metrics ◽  
Rna Seq ◽  
Sequencing Data

Although precise identification of the human leukocyte antigen (HLA) allele is crucial for various clinical and research applications, HLA typing remains challenging due to high polymorphism of the HLA loci. However, with Next-Generation Sequencing (NGS) data becoming widely accessible, many computational tools have been developed to predict HLA types from RNA sequencing (RNA-seq) data. However, there is a lack of comprehensive and systematic benchmarking of RNA-seq HLA callers using large-scale and realist gold standards. In order to address this limitation, we rigorously compared the performance of 12 HLA callers over 50,000 HLA tasks including searching 30 pairwise combinations of HLA callers and reference in over 1,500 samples. In each case, we produced evaluation metrics of accuracy that is the percentage of correctly predicted alleles (two and four-digit resolution) based on six gold standard datasets spanning 650 RNA-seq samples. To determine the influence of the relationship of the read length over the HLA region on prediction quality using each tool, we explored the read length effect by considering read length in the range 37-126 bp, which was available in our gold standard datasets. Moreover, using the Genotype-Tissue Expression (GTEx) v8 data, we carried out evaluation metrics by calculating the concordance of the same HLA type across different tissues from the same individual to evaluate how well the HLA callers can maintain consistent results across various tissues of the same individual. This study offers crucial information for researchers regarding appropriate choices of methods for an HLA analysis.

Curse: building expression atlases and co-expression networks from public RNA-Seq data

2018 ◽  
Vol 35 (16) ◽  
pp. 2880-2881 ◽  
Author(s):  
Dries Vaneechoutte ◽  
Klaas Vandepoele
Keyword(s):  
Rna Sequencing ◽  
State Of The Art ◽  
The State ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Rna Seq ◽  
Valuable Resource ◽  
Sequencing Data ◽  
Imperfect Quality ◽  
Sequence Read Archive

Abstract Summary Public RNA-Sequencing (RNA-Seq) datasets are a valuable resource for transcriptome analyses, but their accessibility is hindered by the imperfect quality and presentation of their metadata and by the complexity of processing raw sequencing data. The Curse suite was created to alleviate these problems. It consists of an online curation tool named Curse to efficiently build compendia of experiments hosted on the Sequence Read Archive, and a lightweight pipeline named Prose to download and process the RNA-Seq data into expression atlases and co-expression networks. Curse networks showed improved linking of functionally related genes compared to the state-of-the-art. Availability and implementation Curse, Prose and their manuals are available at http://bioinformatics.psb.ugent.be/webtools/Curse/. Prose was implemented in Java. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals NeoFuse: predicting fusion neoantigens from RNA sequencing data

2019 ◽  
Vol 36 (7) ◽  
pp. 2260-2261 ◽  
Author(s):  
Georgios Fotakis ◽  
Dietmar Rieder ◽  
Marlene Haider ◽  
Zlatko Trajanoski ◽  
Francesca Finotello
Keyword(s):  
Rna Sequencing ◽  
Human Leukocyte ◽  
Computational Prediction ◽  
Supplementary Information ◽  
Rna Seq ◽  
Sequencing Data ◽  
Leukocyte Antigen ◽  
Computational Workflow ◽  
Hla Binding ◽  
Proteins And Peptides

Abstract Summary Gene fusions can generate immunogenic neoantigens that mediate anticancer immune responses. However, their computational prediction from RNA sequencing (RNA-seq) data requires deep bioinformatics expertise to assembly a computational workflow covering the prediction of: fusion transcripts, their translated proteins and peptides, Human Leukocyte Antigen (HLA) types, and peptide-HLA binding affinity. Here, we present NeoFuse, a computational pipeline for the prediction of fusion neoantigens from tumor RNA-seq data. NeoFuse can be applied to cancer patients’ RNA-seq data to identify fusion neoantigens that might expand the repertoire of suitable targets for immunotherapy. Availability and implementation NeoFuse source code and documentation are available under GPLv3 license at https://icbi.i-med.ac.at/NeoFuse/. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals Putatively cancer-specific exon–exon junctions are shared across patients and present in developmental and other non-cancer cells

NAR Cancer ◽  
2020 ◽  
Vol 2 (1) ◽  
Author(s):  
Julianne K David ◽  
Sean K Maden ◽  
Benjamin R Weeder ◽  
Reid F Thompson ◽  
Abhinav Nellore
Keyword(s):  
Rna Splicing ◽  
Tissue Expression ◽  
The Body ◽  
The Cancer Genome Atlas ◽  
Cancer Type ◽  
Rna Seq ◽  
Normal Tissues ◽  
Cancer Types ◽  
Exon Junctions ◽  
Cancer Tissues

Abstract This study probes the distribution of putatively cancer-specific junctions across a broad set of publicly available non-cancer human RNA sequencing (RNA-seq) datasets. We compared cancer and non-cancer RNA-seq data from The Cancer Genome Atlas (TCGA), the Genotype-Tissue Expression (GTEx) Project and the Sequence Read Archive. We found that (i) averaging across cancer types, 80.6% of exon–exon junctions thought to be cancer-specific based on comparison with tissue-matched samples (σ = 13.0%) are in fact present in other adult non-cancer tissues throughout the body; (ii) 30.8% of junctions not present in any GTEx or TCGA normal tissues are shared by multiple samples within at least one cancer type cohort, and 87.4% of these distinguish between different cancer types; and (iii) many of these junctions not found in GTEx or TCGA normal tissues (15.4% on average, σ = 2.4%) are also found in embryological and other developmentally associated cells. These findings refine the meaning of RNA splicing event novelty, particularly with respect to the human neoepitope repertoire. Ultimately, cancer-specific exon–exon junctions may have a substantial causal relationship with the biology of disease.

schex avoids overplotting for large single-cell RNA-sequencing datasets

2019 ◽  
Vol 36 (7) ◽  
pp. 2291-2292 ◽  
Author(s):  
Saskia Freytag ◽  
Ryan Lister
Keyword(s):  
Single Cell ◽  
Rna Sequencing ◽  
R Package ◽  
Supplementary Information ◽  
Supplementary Data ◽  
Sequencing Data ◽  
Single Cell Rna Sequencing

Abstract Summary Due to the scale and sparsity of single-cell RNA-sequencing data, traditional plots can obscure vital information. Our R package schex overcomes this by implementing hexagonal binning, which has the additional advantages of improving speed and reducing storage for resulting plots. Availability and implementation schex is freely available from Bioconductor via http://bioconductor.org/packages/release/bioc/html/schex.html and its development version can be accessed on GitHub via https://github.com/SaskiaFreytag/schex. Supplementary information Supplementary data are available at Bioinformatics online.

A Two-Stage Poisson Model for Testing RNA-Seq Data

2011 ◽  
Vol 10 (1) ◽  
Author(s):  
Paul L. Auer ◽  
Rebecca W Doerge
Keyword(s):  
Rna Sequencing ◽  
Statistical Approach ◽  
Poisson Model ◽  
Real Data ◽  
Rna Seq ◽  
Sequencing Data ◽  
Sequencing Technology ◽  
Two Stage ◽  
Individual Gene ◽  
Unique Nature

RNA sequencing technology is providing data of unprecedented throughput, resolution, and accuracy. Although there are many different computational tools for processing these data, there are a limited number of statistical methods for analyzing them, and even fewer that acknowledge the unique nature of individual gene transcription. We introduce a simple and powerful statistical approach, based on a two-stage Poisson model, for modeling RNA sequencing data and testing for biologically important changes in gene expression. The advantages of this approach are demonstrated through simulations and real data applications.

RNA sequencing analysis for profiling activation of cancer-associated molecular pathways.

2019 ◽  
Vol 37 (15_suppl) ◽  
pp. e13032-e13032 ◽  
Author(s):  
Anton Buzdin ◽  
Andrew Garazha ◽  
Maxim Sorokin ◽  
Alex Glusker ◽  
Alexey Aleshin ◽  
...  
Keyword(s):  
Gene Expression ◽  
Original Data ◽  
Tissue Expression ◽  
Molecular Pathways ◽  
Sequencing Analysis ◽  
Rna Seq ◽  
Sequencing Data ◽  
Healthy Human ◽  
Tissue Samples ◽  
Normal Tissues

e13032 Background: Intracellular molecular pathways (IMPs) control all major events in the living cell. They are considered hotspots in contemporary oncology because knowledge of IMPs activation is essential for understanding mechanisms of molecular pathogenesis in oncology. Profiling IMPs requires RNA-seq data for tumors and for a collection of reference normal tissues. However, there is a shortage now in such profiles for normal tissues from healthy human donors, uniformly profiled in a single series of experiments. Access to the largest dataset of normal profiles GTEx is only partly available through the dbGaP. In TCGA database, norms are adjacent to surgically removed tumors and may be affected by tumor-linked growth factors, inflammation and altered vascularization. ENCODE datasets were for the autopsies of normal tissues, but they can’t form statistically significant reference groups. Methods: Tissue samples representing 20 organs were taken from post-mortal human healthy donors killed in road accidents no later than 36 hours after death, blood samples were taken from healthy volunteers. Gene expression was profiled in RNA-seq experiments using the same reagents, equipment and protocols. Bioinformatic algorithms for IMP analysis were developed and validated using experimental and public gene expression datasets. Results: From original sequencing data we constructed the biggest fully open reference expression database of normal human tissues including 465 profiles termed Oncobox Atlas of Normal Tissue Expression (ANTE, original data: GSE120795). We next developed a method termed Oncobox for interrogating activation of IMPs in human cancers. It includes modules of expression data harmonization and comparison and an algorithm for automatic annotation of molecular pathways. The Oncobox system enables accurate scoring of thousands molecular pathways using RNA-seq data. Oncobox pathway analysis is also applicable for quantitative proteomics and microRNA data in oncology. Conclusions: The Oncobox system can be used for a plethora of applications in cancer research including finding differentially regulated genes and IMPs, and for discovery of new pathway-related diagnostic and prognostic biomarkers.

Sign in / Sign up

Export Citation Format

Share Document