scholarly journals ReQTL – an allele-level measure of variation-expression genomic relationships

2018 ◽  
Author(s):  
Liam Spurr ◽  
Nawaf Alomran ◽  
Piotr Słowiński ◽  
Muzi Li ◽  
Pavlos Bousounis ◽  
...  
Keyword(s):  
Gene Expression ◽  
Genetic Variation ◽  
Rna Sequencing ◽  
High Performance ◽  
Supplementary Information ◽  
Phenotypic Trait ◽  
Sequencing Data ◽  
Single Nucleotide Variants ◽  
Genomic Relationships ◽  
Supplementary Material

MotivationBy testing for association of DNA genotypes with 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 can be measured at expressed genome regions, and differs from the DNA genotype in sites subjected to regulatory forces. Therefore, assessment of correlation between RNA genetic variation and gene expression can reveal regulatory genomic relationships in addition to eQTLs.ResultsWe introduce ReQTL, an eQTL modification which substitutes the DNA allele count for the variant allele frequency (VAF) at expressed SNV loci in the transcriptome. We exemplify the method on sets of RNA-sequencing data from human tissues obtained though the Genotype-Tissue Expression Project (GTEx) and demonstrate that ReQTL analyses show consistently high performance and sufficient power to identify both previously known and novel molecular associations. The majority of the SNVs implicated in significant cis-ReQTLs identified by our analysis were previously reported as significant cis-eQTL loci. Notably, trans ReQTL loci in our data were substantially enriched in RNA-editing sites. In summary, ReQTL analyses are computationally feasible and do not require matched DNA data, hence they have a high potential to facilitate the discovery of novel molecular interactions through exploration of the increasingly accessible RNA-sequencing datasets.Availability and implementationSample scripts used in our ReQTL analyses are available with the Supplementary Material (ReQTL_sample_code)[email protected] or [email protected] InformationRe_QTL_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 GeneQC: A quality control tool for gene expression estimation based on RNA-sequencing reads mapping

2018 ◽  
Author(s):  
Adam McDermaid ◽  
Xin Chen ◽  
Yiran Zhang ◽  
Juan Xie ◽  
Cankun Wang ◽  
...  
Keyword(s):  
Gene Expression ◽  
Rna Sequencing ◽  
Model Fitting ◽  
Functional Enrichment ◽  
Supplementary Information ◽  
Rna Seq ◽  
Quality Control Tool ◽  
Differential Gene ◽  
Supplementary Material ◽  
The Impact

AbstractMotivationOne of the main benefits of using modern RNA-sequencing (RNA-Seq) technology is the more accurate gene expression estimations compared with previous generations of expression data, such as the microarray. However, numerous issues can result in the possibility that an RNA-Seq read can be mapped to multiple locations on the reference genome with the same alignment scores, which occurs in plant, animal, and metagenome samples. Such a read is so-called a multiple-mapping read (MMR). The impact of these MMRs is reflected in gene expression estimation and all downstream analyses, including differential gene expression, functional enrichment, etc. Current analysis pipelines lack the tools to effectively test the reliability of gene expression estimations, thus are incapable of ensuring the validity of all downstream analyses.ResultsOur investigation into 95 RNA-Seq datasets from seven species (totaling 1,951GB) indicates an average of roughly 22% of all reads are MMRs for plant and animal species. Here we present a tool called GeneQC (Gene expression Quality Control), which can accurately estimate the reliability of each gene’s expression level. The underlying algorithm is designed based on extracted genomic and transcriptomic features, which are then combined using elastic-net regularization and mixture model fitting to provide a clearer picture of mapping uncertainty for each gene. GeneQC allows researchers to determine reliable expression estimations and conduct further analysis on the gene expression that is of sufficient quality. This tool also enables researchers to investigate continued re-alignment methods to determine more accurate gene expression estimates for those with low reliability.AvailabilityGeneQC is freely available at http://bmbl.sdstate.edu/GeneQC/[email protected] informationSupplementary data are available at Bioinformatics online.

scholarly journals scGEApp: a Matlab app for feature selection on single-cell RNA sequencing data

2019 ◽  
Author(s):  
James J. Cai
Keyword(s):  
Gene Expression ◽  
Feature Selection ◽  
Single Cell ◽  
Rna Sequencing ◽  
User Interfaces ◽  
Dropout Rate ◽  
Supplementary Information ◽  
Sequencing Data ◽  
Full Spectrum ◽  
Single Cell Rna Sequencing

AbstractMotivationThe recent development of single-cell technologies, especially single-cell RNA sequencing (scRNA-seq), provides an unprecedented level of resolution to the cell type heterogeneity. It also enables the study of gene expression variability across individual cells within a homogenous cell population. Feature selection algorithms have been used to select biologically meaningful genes while controlling for sampling noise. An easy-to-use application for feature selection on scRNA-seq data requires integration of functions for data filtering, normalization, visualization, and enrichment analyses. Graphic user interfaces (GUIs) are desired for such an application.ResultsWe used native Matlab and App Designer to develop scGEApp for feature selection on singlecell gene expression data. We specifically designed a new feature selection algorithm based on the 3D spline fitting of expression mean (μ), coefficient of variance (CV), and dropout rate (rdrop), making scGEApp a unique tool for feature selection on scRNA-seq data. Our method can be applied to single-sample or two-sample scRNA-seq data, identify feature genes, e.g., those with unexpectedly high CV for given μ and rdrop of those genes, or genes with the most feature changes. Users can operate scGEApp through GUIs to use the full spectrum of functions including normalization, batch effect correction, imputation, visualization, feature selection, and downstream analyses with GSEA and GOrilla.Availabilityhttps://github.com/jamesjcai/scGEAppContact:[email protected] informationSupplementary data are available at Bioinformatics online.

scholarly journals Receptor-driven invasion profiles in diffuse intrinsic pontine glioma (DIPG)

2021 ◽  
Author(s):  
Anju Karki ◽  
Noah E Berlow ◽  
Jin-Ah Kim ◽  
Esther Hulleman ◽  
Qianqian Liu ◽  
...  
Keyword(s):  
Gene Expression ◽  
Growth Factor ◽  
Protein Expression ◽  
Cell Viability ◽  
Rna Sequencing ◽  
Cell Invasion ◽  
Growth Factor Receptor ◽  
Diffuse Intrinsic Pontine Glioma ◽  
Sequencing Data ◽  
Pontine Glioma

Abstract Background Diffuse intrinsic pontine glioma (DIPG) is a devastating pediatric cancer with unmet clinical need. DIPG is invasive in nature, where tumor cells interweave into the fiber nerve tracts of the pons making the tumor unresectable. Accordingly, novel approaches in combating the disease is of utmost importance and receptor-driven cell invasion in the context of DIPG is under-researched area. Here we investigated the impact on cell invasion mediated by PLEXINB1, PLEXINB2, platelet growth factor receptor (PDGFR)α, PDGFRβ, epithelial growth factor receptor (EGFR), activin receptor 1 (ACVR1), chemokine receptor 4 (CXCR4) and NOTCH1. Methods We used previously published RNA-sequencing data to measure gene expression of selected receptors in DIPG tumor tissue versus matched normal tissue controls (n=18). We assessed protein expression of the corresponding genes using DIPG cell culture models. Then, we performed cell viability and cell invasion assays of DIPG cells stimulated with chemoattractants/ligands. Results RNA-sequencing data showed increased gene expression of receptor genes such as PLEXINB2, PDGFRα, EGFR, ACVR1, CXCR4 and NOTCH1 in DIPG tumors compared to the control tissues. Representative DIPG cell lines demonstrated correspondingly increased protein expression levels of these genes. Cell viability assays showed minimal effects of growth factors/chemokines on tumor cell growth in most instances. Recombinant SEMA4C, SEM4D, PDGF-AA, PDGF-BB, ACVA, CXCL12 and DLL4 ligand stimulation altered invasion in DIPG cells. Conclusions We show that no single growth factor-ligand pair universally induces DIPG cell invasion. However, our results reveal a potential to create a composite of cytokines or anti-cytokines to modulate DIPG cell invasion.

FC 011KIDNEYNETWORK: USING KIDNEY DERIVED GENE EXPRESSION DATA TO PREDICT AND PRIORITIZE NOVEL GENES INVOLVED IN KIDNEY DISEASE

2021 ◽  
Vol 36 (Supplement_1) ◽  
Author(s):  
Floranne Boulogne ◽  
Laura Claus ◽  
Henry Wiersma ◽  
Roy Oelen ◽  
Floor Schukking ◽  
...  
Keyword(s):  
Gene Expression ◽  
Kidney Disease ◽  
Candidate Gene ◽  
Exome Sequencing ◽  
Rna Sequencing ◽  
Expression Patterns ◽  
Genetic Diagnosis ◽  
Specific Gene ◽  
Sequencing Data ◽  
Exome Sequencing Data

Abstract Background and Aims Genetic testing in patients with suspected hereditary kidney disease does not always reveal the genetic cause for the patient's disorder. Potentially pathogenic variants can reside in genes that are not known to be involved in kidney disease, which makes it difficult to prioritize and interpret the relevance of these variants. As such, there is a clear need for methods that predict the phenotypic consequences of gene expression in a way that is as unbiased as possible. To help identify candidate genes we have developed KidneyNetwork, in which tissue-specific expression is utilized to predict kidney-specific gene functions. Method We combined gene co-expression in 878 publicly available kidney RNA-sequencing samples with the co-expression of a multi-tissue RNA-sequencing dataset of 31,499 samples to build KidneyNetwork. The expression patterns were used to predict which genes have a kidney-related function, and which (disease) phenotypes might be caused when these genes are mutated. By integrating the information from the HPO database, in which known phenotypic consequences of disease genes are annotated, with the gene co-expression network we obtained prediction scores for each gene per HPO term. As proof of principle, we applied KidneyNetwork to prioritize variants in exome-sequencing data from 13 kidney disease patients without a genetic diagnosis. Results We assessed the prediction performance of KidneyNetwork by comparing it to GeneNetwork, a multi-tissue co-expression network we previously developed. In KidneyNetwork, we observe a significantly improved prediction accuracy of kidney-related HPO-terms, as well as an increase in the total number of significantly predicted kidney-related HPO-terms (figure 1). To examine its clinical utility, we applied KidneyNetwork to 13 patients with a suspected hereditary kidney disease without a genetic diagnosis. Based on the HPO terms “Renal cyst” and “Hepatic cysts”, combined with a list of potentially damaging variants in one of the undiagnosed patients with mild ADPKD/PCLD, we identified ALG6 as a new candidate gene. ALG6 bears a high resemblance to other genes implicated in this phenotype in recent years. Through the 100,000 Genomes Project and collaborators we identified three additional patients with kidney and/or liver cysts carrying a suspected deleterious variant in ALG6. Conclusion We present KidneyNetwork, a kidney specific co-expression network that accurately predicts what genes have kidney-specific functions and may result in kidney disease. Gene-phenotype associations of genes unknown for kidney-related phenotypes can be predicted by KidneyNetwork. We show the added value of KidneyNetwork by applying it to exome sequencing data of kidney disease patients without a molecular diagnosis and consequently we propose ALG6 as a promising candidate gene. KidneyNetwork can be applied to clinically unsolved kidney disease cases, but it can also be used by researchers to gain insight into individual genes to better understand kidney physiology and pathophysiology. Acknowledgments This research was made possible through access to the data and findings generated by the 100,000 Genomes Project; http://www.genomicsengland.co.uk.

scholarly journals LSTrAP-Crowd: Prediction of novel components of bacterial ribosomes with crowd-sourced analysis of RNA sequencing data

2020 ◽  
Author(s):  
Benedict Hew ◽  
Qiao Wen Tan ◽  
William Goh ◽  
Jonathan Wei Xiong Ng ◽  
Kenny Koh ◽  
...  
Keyword(s):  
Gene Expression ◽  
Protein Synthesis ◽  
Rna Sequencing ◽  
Gene Expression Data ◽  
Large Scale ◽  
Bacterial Resistance ◽  
Expression Data ◽  
Sequencing Data ◽  
Novel Proteins ◽  
Novel Antibiotics

AbstractBacterial resistance to antibiotics is a growing problem that is projected to cause more deaths than cancer in 2050. Consequently, novel antibiotics are urgently needed. Since more than half of the available antibiotics target the bacterial ribosomes, proteins that are involved in protein synthesis are thus prime targets for the development of novel antibiotics. However, experimental identification of these potential antibiotic target proteins can be labor-intensive and challenging, as these proteins are likely to be poorly characterized and specific to few bacteria. In order to identify these novel proteins, we established a Large-Scale Transcriptomic Analysis Pipeline in Crowd (LSTrAP-Crowd), where 285 individuals processed 26 terabytes of RNA-sequencing data of the 17 most notorious bacterial pathogens. In total, the crowd processed 26,269 RNA-seq experiments and used the data to construct gene co-expression networks, which were used to identify more than a hundred uncharacterized genes that were transcriptionally associated with protein synthesis. We provide the identity of these genes together with the processed gene expression data. The data can be used to identify other vulnerabilities or bacteria, while our approach demonstrates how the processing of gene expression data can be easily crowdsourced.

scholarly journals Single-Cell Transcriptome Analysis Reveals Dynamic Cell Populations and Differential Gene Expression Patterns in Control and Aneurysmal Human Aortic Tissue

Circulation ◽  
2020 ◽  
Vol 142 (14) ◽  
pp. 1374-1388
Author(s):  
Yanming Li ◽  
Pingping Ren ◽  
Ashley Dawson ◽  
Hernan G. Vasquez ◽  
Waleed Ageedi ◽  
...  
Keyword(s):  
Gene Expression ◽  
Single Cell ◽  
Rna Sequencing ◽  
Aortic Wall ◽  
Genome Wide Association ◽  
Aortic Tissue ◽  
Sequencing Data ◽  
Genome Wide ◽  
Single Cell Rna Sequencing ◽  
Differential Gene

Background: Ascending thoracic aortic aneurysm (ATAA) is caused by the progressive weakening and dilatation of the aortic wall and can lead to aortic dissection, rupture, and other life-threatening complications. To improve our understanding of ATAA pathogenesis, we aimed to comprehensively characterize the cellular composition of the ascending aortic wall and to identify molecular alterations in each cell population of human ATAA tissues. Methods: We performed single-cell RNA sequencing analysis of ascending aortic tissues from 11 study participants, including 8 patients with ATAA (4 women and 4 men) and 3 control subjects (2 women and 1 man). Cells extracted from aortic tissue were analyzed and categorized with single-cell RNA sequencing data to perform cluster identification. ATAA-related changes were then examined by comparing the proportions of each cell type and the gene expression profiles between ATAA and control tissues. We also examined which genes may be critical for ATAA by performing the integrative analysis of our single-cell RNA sequencing data with publicly available data from genome-wide association studies. Results: We identified 11 major cell types in human ascending aortic tissue; the high-resolution reclustering of these cells further divided them into 40 subtypes. Multiple subtypes were observed for smooth muscle cells, macrophages, and T lymphocytes, suggesting that these cells have multiple functional populations in the aortic wall. In general, ATAA tissues had fewer nonimmune cells and more immune cells, especially T lymphocytes, than control tissues did. Differential gene expression data suggested the presence of extensive mitochondrial dysfunction in ATAA tissues. In addition, integrative analysis of our single-cell RNA sequencing data with public genome-wide association study data and promoter capture Hi-C data suggested that the erythroblast transformation-specific related gene( ERG ) exerts an important role in maintaining normal aortic wall function. Conclusions: Our study provides a comprehensive evaluation of the cellular composition of the ascending aortic wall and reveals how the gene expression landscape is altered in human ATAA tissue. The information from this study makes important contributions to our understanding of ATAA formation and progression.

scholarly journals neoepiscope improves neoepitope prediction with multivariant phasing

2019 ◽  
Vol 36 (3) ◽  
pp. 713-720 ◽  
Author(s):  
Mary A Wood ◽  
Austin Nguyen ◽  
Adam J Struck ◽  
Kyle Ellrott ◽  
Abhinav Nellore ◽  
...  
Keyword(s):  
False Negative ◽  
Supplementary Information ◽  
Supplementary File ◽  
Sequencing Data ◽  
Single Nucleotide Variants ◽  
Single Nucleotide ◽  
Somatic Variant ◽  
Negative Results ◽  
Multiple Datasets ◽  
False Negative Results

Abstract Motivation The vast majority of tools for neoepitope prediction from DNA sequencing of complementary tumor and normal patient samples do not consider germline context or the potential for the co-occurrence of two or more somatic variants on the same mRNA transcript. Without consideration of these phenomena, existing approaches are likely to produce both false-positive and false-negative results, resulting in an inaccurate and incomplete picture of the cancer neoepitope landscape. We developed neoepiscope chiefly to address this issue for single nucleotide variants (SNVs) and insertions/deletions (indels). Results Herein, we illustrate how germline and somatic variant phasing affects neoepitope prediction across multiple datasets. We estimate that up to ∼5% of neoepitopes arising from SNVs and indels may require variant phasing for their accurate assessment. neoepiscope is performant, flexible and supports several major histocompatibility complex binding affinity prediction tools. Availability and implementation neoepiscope is available on GitHub at https://github.com/pdxgx/neoepiscope under the MIT license. Scripts for reproducing results described in the text are available at https://github.com/pdxgx/neoepiscope-paper under the MIT license. Additional data from this study, including summaries of variant phasing incidence and benchmarking wallclock times, are available in Supplementary Files 1, 2 and 3. Supplementary File 1 contains Supplementary Table 1, Supplementary Figures 1 and 2, and descriptions of Supplementary Tables 2–8. Supplementary File 2 contains Supplementary Tables 2–6 and 8. Supplementary File 3 contains Supplementary Table 7. Raw sequencing data used for the analyses in this manuscript are available from the Sequence Read Archive under accessions PRJNA278450, PRJNA312948, PRJNA307199, PRJNA343789, PRJNA357321, PRJNA293912, PRJNA369259, PRJNA305077, PRJNA306070, PRJNA82745 and PRJNA324705; from the European Genome-phenome Archive under accessions EGAD00001004352 and EGAD00001002731; and by direct request to the authors. Supplementary information Supplementary data are available at Bioinformatics online.

Sign in / Sign up

Export Citation Format

Share Document