PRISM: Recovering cell type specific expression profiles from individual composite RNA-seq samples.
Ontology highlight
ABSTRACT: Motivation
A major challenge in analyzing cancer patient transcriptomes is that the tumors are inherently heterogeneous and evolving. We analyzed 214 bulk RNA samples of a longitudinal, prospective ovarian cancer cohort and found that the sample composition changes systematically due to chemotherapy and between the anatomical sites, preventing direct comparison of treatment-naive and treated samples.Results
To overcome this, we developed PRISM, a latent statistical framework to simultaneously extract the sample composition and cell type specific whole-transcriptome profiles adapted to each individual sample. Our results indicate that the PRISM-derived composition-free transcriptomic profiles and signatures derived from them predict the patient response better than the composite raw bulk data. We validated our findings in independent ovarian cancer and melanoma cohorts, and verified that PRISM accurately estimates the composition and cell type specific expression through whole-genome sequencing and RNA in situ hybridization experiments.Availability
https://bitbucket.org/anthakki/prism.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Hakkinen A
PROVIDER: S-EPMC8479664 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Conventionally, overall gene expressions from microarrays are used to infer gene networks, but it is challenging to account splicing isoforms. High-throughput RNA Sequencing has made splice variant profiling practical. However, its true merit in quantifying splicing isoforms and isoform-specific exon expressions is not well explored in inferring gene networks. This study demonstrates SpliceNet, a method to infer isoform-specific co-expression networks from exon-level RNA-Seq data, using large dimensional trace. It goes beyond differentially expressed genes and infers splicing isoform network changes between normal and diseased samples. It eases the sample size bottleneck; evaluations on simulated data and lung cancer-specific ERBB2 and MAPK signaling pathways, with varying number of samples, evince the merit in handling high exon to sample size ratio datasets. Inferred network rewiring of well established Bcl-x and EGFR centered networks from lung adenocarcinoma expression data is in good agreement with literature. Gene level evaluations demonstrate a substantial performance of SpliceNet over canonical correlation analysis, a method that is currently applied to exon level RNA-Seq data. SpliceNet can also be applied to exon array data. SpliceNet is distributed as an R package available at http://www.jjwanglab.org/SpliceNet.
| S-EPMC4150760 | biostudies-literature
Project description:Cell-type specific gene expression profiles are needed for many computational methods operating on bulk RNA-Seq samples, such as deconvolution of cell-type fractions and digital cytometry. However, the gene expression profile of a cell type can vary substantially due to both technical factors and biological differences in cell state and surroundings, reducing the efficacy of such methods. Here, we investigated which factors contribute most to this variation. We evaluated different normalization methods, quantified the variance explained by different factors, evaluated the effect on deconvolution of cell type fractions, and examined the differences between UMI-based single-cell RNA-Seq and bulk RNA-Seq. We investigated a collection of publicly available bulk and single-cell RNA-Seq datasets containing B and T cells, and found that the technical variation across laboratories is substantial, even for genes specifically selected for deconvolution, and this variation has a confounding effect on deconvolution. Tissue of origin is also a substantial factor, highlighting the challenge of using cell type profiles derived from blood with mixtures from other tissues. We also show that much of the differences between UMI-based single-cell and bulk RNA-Seq methods can be explained by the number of read duplicates per mRNA molecule in the single-cell sample. Our work shows the importance of either matching or correcting for technical factors when creating cell-type specific gene expression profiles that are to be used together with bulk samples.
| S-EPMC7505444 | biostudies-literature
Project description:Small RNAs (sRNAs) are bioactive molecules that can be detected in biofluids, reflecting physiological and pathological states. In plasma, sRNAs are found within extracellular vesicles (EVs) and in extravesicular compartments, offering potential sources of highly sensitive biomarkers. Deep sequencing strategies to profile sRNAs favor the detection of microRNAs (miRNAs), the best-known class of sRNAs. Phospho-RNA-seq, through the enzymatic treatment of sRNAs with T4 polynucleotide kinase (T4-PNK), has been recently developed to increase the detection of thousands of previously inaccessible RNAs. In this study, we investigated the value of phospho-RNA-seq on both the EVs and extravesicular plasma subfractions. Phospho-RNA-seq increased the proportion of sRNAs used for alignment and highlighted the diversity of the sRNA transcriptome. Unsupervised clustering analysis using sRNA counts matrices correctly classified the EVs and extravesicular samples only in the T4-PNK treated samples, indicating that phospho-RNA-seq stresses the features of sRNAs in each plasma subfraction. Furthermore, T4-PNK treatment emphasized specific miRNA variants differing in the 5'-end (5'-isomiRs) and certain types of tRNA fragments in each plasma fraction. Phospho-RNA-seq increased the number of tissue-specific messenger RNA (mRNA) fragments in the EVs compared with the extravesicular fraction, suggesting that phospho-RNA-seq favors the discovery of tissue-specific sRNAs in EVs. Overall, the present data emphasizes the value of phospho-RNA-seq in uncovering RNA-based biomarkers in EVs.
| S-EPMC10380198 | biostudies-literature
Project description:In this study, mutations present in a series of human melanomas (stage IV disease) will be determined, using autologous blood cells to obtain a reference genome. From each of the samples that are analyzed, tumour-infiltrating T lymphocytes have also been isolated. This offers a unique opportunity to determine which (fraction of) mutations in human cancer leads to epitopes that are recognized by T cells. The resulting information is likely to be of value to understand how T cell activating drugs exert their action.
| S-ECPF-ERAD-61 | biostudies-other
Project description:MotivationHigh-throughput RNA sequencing has revolutionized the scope and depth of transcriptome analysis. Accurate reconstruction of a phenotype-specific transcriptome is challenging due to the noise and variability of RNA-seq data. This requires computational identification of transcripts from multiple samples of the same phenotype, given the underlying consensus transcript structure.ResultsWe present a Bayesian method, integrated assembly of phenotype-specific transcripts (IntAPT), that identifies phenotype-specific isoforms from multiple RNA-seq profiles. IntAPT features a novel two-layer Bayesian model to capture the presence of isoforms at the group layer and to quantify the abundance of isoforms at the sample layer. A spike-and-slab prior is used to model the isoform expression and to enforce the sparsity of expressed isoforms. Dependencies between the existence of isoforms and their expression are modeled explicitly to facilitate parameter estimation. Model parameters are estimated iteratively using Gibbs sampling to infer the joint posterior distribution, from which the presence and abundance of isoforms can reliably be determined. Studies using both simulations and real datasets show that IntAPT consistently outperforms existing methods for the IntAPT. Experimental results demonstrate that, despite sequencing errors, IntAPT exhibits a robust performance among multiple samples, resulting in notably improved identification of expressed isoforms of low abundance.Availability and implementationThe IntAPT package is available at http://github.com/henryxushi/IntAPT.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC8097681 | biostudies-literature
Project description:Avian pathogenic Escherichia coli (APEC) causes great economic loss to the poultry industry worldwide. Chicken type II pneumocytes (CP II cells) secrete surfactants and modulate lung immunity to decrease the infection of the invading pathogen. Nevertheless, the pathogenesis of CP II cells to APEC infection remains poorly understood. Therefore, we conducted global gene expression profiling of CP II cells after APEC-O78 infection to explore the host-pathogen interaction. The differentially expressed genes of CP II cells to APEC infection were characterized by RNA-seq with EB-seq algorithm. In consequence, the mRNA of 18996 genes was identified, and CP II cells responded to APEC infection with marked changes in the expression of 1390 genes. Among them, there are 803 down-regulated mRNAs and 587 up-regulated mRNAs. The KEGG prediction and Gene Ontology terms analysis revealed that the major enriched pathways were related to NF-κB signaling pathway, apoptosis pathway, tight junction, and cytokine-cytokine receptor interaction and other pathways. We adopted qRT-PCR to verify the validity of the selected gene expression. The fold induction of qPCR was similar to the RNA-seq results. These results provide a better understanding of the pathogenesis of APEC, especially apoptosis pathway involved in APEC infection.
| S-EPMC6550405 | biostudies-literature
Project description:In order to study molecular changes in the stroma from tissue samples it is recommended to separate tumor tissue from stromal tissue. This is particularly relevant to mouse tumor xenograft models where tumor, particularly metastatic tumors, can be small and difficult to separate from the host tissue. In our research we compared qualitatively the ability of high-throughput mRNA sequencing, RNA-Seq, and microarrays to detect tumor (human) and stromal (mouse) expression from mixed tumor-stromal samples in terms of the genes and pathways that are involved in cross-alignment (RNA-Seq) and cross-hybridization (microarrays). Human samples consisted of total RNA obtained from MDA-MB-231 human breast carcinoma cell line and isolated from three independent cultures of sub-confluent MDA-MB-231 cell lines in exponential phase of growth. Mouse samples were obtained from NOD scid gamma mice, and normal lung tissue was harvested from three independent age-matched mice.
| S-ECPF-GEOD-40890 | biostudies-other
Project description:The L1000 technology, a cost-effective high-throughput transcriptomics technology, has been applied to profile a collection of human cell lines for their gene expression response to > 30,000 chemical and genetic perturbations. In total, there are currently over 3 million available L1000 profiles. Such a dataset is invaluable for the discovery of drug and target candidates and for inferring mechanisms of action for small molecules. The L1000 assay only measures the mRNA expression of 978 landmark genes while 11,350 additional genes are computationally reliably inferred. The lack of full genome coverage limits knowledge discovery for half of the human protein coding genes, and the potential for integration with other transcriptomics profiling data. Here we present a Deep Learning two-step model that transforms L1000 profiles to RNA-seq-like profiles. The input to the model are the measured 978 landmark genes while the output is a vector of 23,614 RNA-seq-like gene expression profiles. The model first transforms the landmark genes into RNA-seq-like 978 gene profiles using a modified CycleGAN model applied to unpaired data. The transformed 978 RNA-seq-like landmark genes are then extrapolated into the full genome space with a fully connected neural network model. The two-step model achieves 0.914 Pearson's correlation coefficients and 1.167 root mean square errors when tested on a published paired L1000/RNA-seq dataset produced by the LINCS and GTEx programs. The processed RNA-seq-like profiles are made available for download, signature search, and gene centric reverse search with unique case studies.
| S-EPMC9472394 | biostudies-literature
Project description:The importance of cell type-specific gene expression in disease-relevant tissues is increasingly recognized in genetic studies of complex diseases. However, most gene expression studies are conducted on bulk tissues, without examining cell type-specific expression profiles. Several computational methods are available for cell type deconvolution (i.e. inference of cellular composition) from bulk RNA-Seq data, but few of them impute cell type-specific expression profiles. We hypothesize that with external prior information such as single cell RNA-seq and population-wide expression profiles, it can be computationally tractable to estimate both cellular composition and cell type-specific expression from bulk RNA-Seq data. Here we introduce CellR, which addresses cross-individual gene expression variations to adjust the weights of cell-specific gene markers. It then transforms the deconvolution problem into a linear programming model while taking into account inter/intra cellular correlations and uses a multi-variate stochastic search algorithm to estimate the cell type-specific expression profiles. Analyses on several complex diseases such as schizophrenia, Alzheimer's disease, Huntington's disease and type 2 diabetes validated the efficiency of CellR, while revealing how specific cell types contribute to different diseases. In summary, CellR compares favorably against competing approaches, enabling cell type-specific re-analysis of gene expression data on bulk tissues in complex diseases.
| S-EPMC8219045 | biostudies-literature
Project description:Liquid biopsy is expected to be a promising cancer screening method because of its low invasiveness and the possibility of detecting multiple types in a single test. In the last decade, many studies on cancer detection using small RNAs in blood have been reported. To put small RNA tests into practical use as a multiple cancer type screening test, it is necessary to develop a method that can be applied to multiple facilities. We collected samples of eight cancer types and healthy controls from 20 facilities to evaluate the performance of cancer type classification. A total of 2,475 cancer samples and 496 healthy control samples were collected using a standardized protocol. After obtaining a small RNA expression profile, we constructed a classification model and evaluated its performance. First, we investigated the classification performance using samples from five single facilities. Each model showed areas under the receiver curve (AUC) ranging from 0.67 to 0.89. Second, we performed principal component analysis (PCA) to examine the characteristics of the facilities. The degree of hemolysis and the data acquisition period affected the expression profiles. Finally, we constructed the classification model by reducing the influence of these factors, and its performance had an AUC of 0.76. The results reveal that small RNA can be used for the classification of cancer types in samples from a single facility. However, interfacility biases will affect the classification of samples from multiple facilities. These findings will provide important insights to improve the performance of multiple cancer type classifications using small RNA expression profiles acquired from multiple facilities.
| S-EPMC9207371 | biostudies-literature