scholarly journals QUBIC2: a novel and robust biclustering algorithm for analyses and interpretation of large-scale RNA-Seq data

2019 ◽  
Vol 36 (4) ◽  
pp. 1143-1149 ◽  
Author(s):  
Juan Xie ◽  
Anjun Ma ◽  
Yu Zhang ◽  
Bingqiang Liu ◽  
Sha Cao ◽  
...  
Keyword(s):  
Gene Expression ◽  
Gene Expression Data ◽  
Large Scale ◽  
Gaussian Model ◽  
Functional Gene ◽  
Superior Performance ◽  
Supplementary Information ◽  
Expression Data ◽  
Rna Seq ◽  
Gene Modules

Abstract Motivation The biclustering of large-scale gene expression data holds promising potential for detecting condition-specific functional gene modules (i.e. biclusters). However, existing methods do not adequately address a comprehensive detection of all significant bicluster structures and have limited power when applied to expression data generated by RNA-Sequencing (RNA-Seq), especially single-cell RNA-Seq (scRNA-Seq) data, where massive zero and low expression values are observed. Results We present a new biclustering algorithm, QUalitative BIClustering algorithm Version 2 (QUBIC2), which is empowered by: (i) a novel left-truncated mixture of Gaussian model for an accurate assessment of multimodality in zero-enriched expression data, (ii) a fast and efficient dropouts-saving expansion strategy for functional gene modules optimization using information divergency and (iii) a rigorous statistical test for the significance of all the identified biclusters in any organism, including those without substantial functional annotations. QUBIC2 demonstrated considerably improved performance in detecting biclusters compared to other five widely used algorithms on various benchmark datasets from E.coli, Human and simulated data. QUBIC2 also showcased robust and superior performance on gene expression data generated by microarray, bulk RNA-Seq and scRNA-Seq. Availability and implementation The source code of QUBIC2 is freely available at https://github.com/OSU-BMBL/QUBIC2. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals Inferring TF activities and activity regulators from gene expression data with constraints from TF perturbation data

2020 ◽  
Author(s):  
Cynthia Z Ma ◽  
Michael R Brent
Keyword(s):  
Gene Expression ◽  
Gene Expression Data ◽  
Large Scale ◽  
Supplementary Information ◽  
Activity Levels ◽  
Expression Data ◽  
Gene Expression Matrix ◽  
Perturbation Data ◽  
Carry Over ◽  
Necessary And Sufficient

Abstract Motivation The activity of a transcription factor (TF) in a sample of cells is the extent to which it is exerting its regulatory potential. Many methods of inferring TF activity from gene expression data have been described, but due to the lack of appropriate large-scale datasets, systematic and objective validation has not been possible until now. Results We systematically evaluate and optimize the approach to TF activity inference in which a gene expression matrix is factored into a condition-independent matrix of control strengths and a condition-dependent matrix of TF activity levels. We find that expression data in which the activities of individual TFs have been perturbed are both necessary and sufficient for obtaining good performance. To a considerable extent, control strengths inferred using expression data from one growth condition carry over to other conditions, so the control strength matrices derived here can be used by others. Finally, we apply these methods to gain insight into the upstream factors that regulate the activities of yeast TFs Gcr2, Gln3, Gcn4 and Msn2. Availability and implementation Evaluation code and data are available at https://doi.org/10.5281/zenodo.4050573. Contact [email protected] Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals LSTrAP-Kingdom: an automated pipeline to generate annotated gene expression atlases for kingdoms of life

2021 ◽  
Author(s):  
William Goh1 ◽  
Marek Mutwil1
Keyword(s):  
Gene Expression ◽  
Large Scale ◽  
Supplementary Information ◽  
Expression Data ◽  
Supplementary Data ◽  
Rna Seq ◽  
Analysis Pipeline ◽  
Study Gene Expression ◽  
Automated Pipeline ◽  
Bacteria And Fungi

Abstract Motivation There are now more than two million RNA sequencing experiments for plants, animals, bacteria and fungi publicly available, allowing us to study gene expression within and across species and kingdoms. However, the tools allowing the download, quality control and annotation of this data for more than one species at a time are currently missing. Results To remedy this, we present the Large-Scale Transcriptomic Analysis Pipeline in Kingdom of Life (LSTrAP-Kingdom) pipeline, which we used to process 134,521 RNA-seq samples, achieving ∼12,000 processed samples per day. Our pipeline generated quality-controlled, annotated gene expression matrices that rival the manually curated gene expression data in identifying functionally-related genes. Availability LSTrAP-Kingdom is available from: https://github.com/wirriamm/plants-pipeline and is fully implemented in Python and Bash. Supplementary information Supplementary data are available at Bioinformatics online.

scholarly journals runibic: a Bioconductor package for parallel row-based biclustering of gene expression data

2017 ◽  
Author(s):  
Patryk Orzechowski ◽  
Artur Pańszczyk ◽  
Xiuzhen Huang ◽  
Jason H. Moore
Keyword(s):  
Gene Expression ◽  
Gene Expression Data ◽  
Memory Management ◽  
Large Scale ◽  
R Package ◽  
Supplementary Information ◽  
Sequential Algorithm ◽  
Expression Data ◽  
Bioconductor Package ◽  
Parallel Version

AbstractMotivationBiclustering (called also co-clustering) is an unsupervised technique of simultaneous analysis of rows and columns of input matrix. From the first application to gene expression data, multiple algorithms have been proposed. Only a handful of them were able to provide accurate results and were fast enough to be suitable for large-scale genomic datasets.ResultsIn this paper we introduce a Bioconductor package with parallel version of UniBic biclustering algorithm: one of the most accurate biclustering methods that have been developed so far. For the convenience of usage, we have wrapped the algorithm in an R package called runibic. The package includes: (1) a couple of times faster parallel version of the original sequential algorithm,(2) muchmore efficient memory management, (3) modularity which allows to build new methods on top of the provided one, (4) integration with the modern Bioconductor packages such as SummarizedExperiment, ExpressionSetand biclust.AvailabilityThe package is implemented in R (3.4) and will be available in the new release of Bioconductor (3.6). Currently it could be downloaded from the following URL: http://github.com/athril/runibic/[email protected], [email protected] informationSupplementary informations are available in vignette of the package.

scholarly journals LSTrAP-Kingdom: an automated pipeline to generate annotated gene expression atlases for kingdoms of life

2021 ◽  
Author(s):  
William Goh ◽  
Marek Mutwil
Keyword(s):  
Gene Expression ◽  
Quality Control ◽  
Gene Expression Data ◽  
Large Scale ◽  
Expression Data ◽  
Rna Seq ◽  
Analysis Pipeline ◽  
Study Gene Expression ◽  
Automated Pipeline ◽  
Bacteria And Fungi

AbstractSummaryThere are now more than two million RNA sequencing experiments for plants, animals, bacteria and fungi publicly available, allowing us to study gene expression within and across species and kingdoms. However, the tools allowing the download, quality control and annotation of this data for more than one species at a time are currently missing. To remedy this, we present the Large-Scale Transcriptomic Analysis Pipeline in Kingdom of Life (LSTrAP-Kingdom) pipeline, which we used to process 134,521 RNA-seq samples, achieving ~12,000 processed samples per day. Our pipeline generated quality-controlled, annotated gene expression matrices that rival the manually curated gene expression data in identifying functionally-related genes.Availability and implementationLSTrAP-Kingdom is available from: https://github.com/wirriamm/plants-pipeline and is fully implemented in Python and Bash.

A gp78/AMFR protein-driven gene signature that predicts breast cancer outcome.

2021 ◽  
Vol 39 (15_suppl) ◽  
pp. 553-553
Author(s):  
Sandeep K. Singhal ◽  
Jung Byun ◽  
Samson Park ◽  
Tingfen Yan ◽  
Ryan Yancey ◽  
...  
Keyword(s):  
Breast Cancer ◽  
Gene Expression ◽  
Protein Expression ◽  
Gene Expression Data ◽  
Independent Predictor ◽  
Response To Therapy ◽  
Expression Data ◽  
Rna Seq ◽  
Gene Expression Signatures ◽  
Gene Modules

553 Background: gp78, also known as the autocrine motility factor receptor (AMFR) or RNF45, is a polytopic RING-type E3 ubiquitin ligase resident to the endoplasmic reticulum (ER) that plays major role in the cellular response to stress by regulating ER homeostasis and signaling through its participation in the unfolded protein response (UPR) and ER associated degradation. We used machine learning (ML) and statistical modeling (SM) to assess gp78 as a protein biomarker that is an independent predictor of breast cancer (bc) survival exclusively in women of self-reported African descent as opposed to European ancestry. Methods: We examined a cohort of racially diverse 555 BC bc patients who underwent surgery for their primary BC in Greenville, NC using ML and SM approach. We leveraged the availability of RNA-seq gene expression data on a portion of our bc cohort (N=136 of 555) to construct gene expression signatures. Results: Using antibodies developed in the Weissman lab and established methods for quantitative IHC, we have found that gp78 expression is significantly increased in the tumors of bc patients compared to normal breast epithelia. In addition, we found that gp78 is expressed at significantly higher levels in bc of non-Hispanic black women (NHB) compared to non-Hispanic white women (NHW) (p=0.0038), and that bc subtypes known to be more aggressive and associated with higher grades like, Basal (p=1.6e-12), Luminal B (p=2.3e-4) and HER2(8.3e-4), display significantly higher levels of gp78 compare to Luminal A. Moreover, Kaplan-Meier survival curve analyses show that gp78 protein expression is more significantly associated with poor survival in NHB women (HR:1.65, p=0.073) compared to NHW women (HR:2.01, p=0.004). Finally, multivariate analysis reveals that gp78 protein expression, based on quantitative IHC, is an independent predictor of poor bc survival exclusively in women of African (NHB) ancestry (HR:1.99, p=0.017). We leveraged the availability of RNA-seq gene expression data on a portion of our bc cohort to construct gene expression signatures or gene modules. An analysis of pooled publicly available data from 845 patients that underwent neoadjuvant chemotherapy for bc (primarily taxane and anthracycline based), reveals that gp78 gene modules are highly predictive of patient response to therapy. gp78-derived gene modules show both high fold difference and significance in predicting response to therapy (AUC:0.72) which is very similar to other multi-gene panels that are currently in clinical use including Prosigna, MammaPrint, and Oncotype Dx. Conclusions: Our results show that gp78/AMFR is an independent predictor of bc survival and response to therapy, based on race, thus implicating a role for this protein, and potentially the UPR, as underlying biological differences in tumor properties linked to genetic ancestry.

Gene Set Correlation Analysis and Visualization Using Gene Expression Data

2020 ◽  
Vol 15 ◽  
Author(s):  
Chen-An Tsai ◽  
James J. Chen
Keyword(s):  
Gene Expression ◽  
Correlation Analysis ◽  
Gene Expression Data ◽  
Differentially Expressed Gene ◽  
Differentially Expressed ◽  
Superior Performance ◽  
Expression Data ◽  
Gene Set ◽  
Gene Sets ◽  
Set Correlation

Background: Gene set enrichment analyses (GSEA) provide a useful and powerful approach to identify differentially expressed gene sets with prior biological knowledge. Several GSEA algorithms have been proposed to perform enrichment analyses on groups of genes. However, many of these algorithms have focused on identification of differentially expressed gene sets in a given phenotype. Objective: In this paper, we propose a gene set analytic framework, Gene Set Correlation Analysis (GSCoA), that simultaneously measures within and between gene sets variation to identify sets of genes enriched for differential expression and highly co-related pathways. Methods: We apply co-inertia analysis to the comparisons of cross-gene sets in gene expression data to measure the costructure of expression profiles in pairs of gene sets. Co-inertia analysis (CIA) is one multivariate method to identify trends or co-relationships in multiple datasets, which contain the same samples. The objective of CIA is to seek ordinations (dimension reduction diagrams) of two gene sets such that the square covariance between the projections of the gene sets on successive axes is maximized. Simulation studies illustrate that CIA offers superior performance in identifying corelationships between gene sets in all simulation settings when compared to correlation-based gene set methods. Result and Conclusion: We also combine between-gene set CIA and GSEA to discover the relationships between gene sets significantly associated with phenotypes. In addition, we provide a graphical technique for visualizing and simultaneously exploring the associations of between and within gene sets and their interaction and network. We then demonstrate integration of within and between gene sets variation using CIA and GSEA, applied to the p53 gene expression data using the c2 curated gene sets. Ultimately, the GSCoA approach provides an attractive tool for identification and visualization of novel associations between pairs of gene sets by integrating co-relationships between gene sets into gene set analysis.

scholarly journals Graph Convolutional Network for Drug Response Prediction Using Gene Expression Data

Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 772
Author(s):  
Seonghun Kim ◽  
Seockhun Bae ◽  
Yinhua Piao ◽  
Kyuri Jo
Keyword(s):  
Gene Expression ◽  
Gene Expression Data ◽  
Large Scale ◽  
Drug Response ◽  
Response Prediction ◽  
Biological Data ◽  
Expression Data ◽  
Convolutional Network ◽  
Essential Information ◽  
Protein Protein Interaction

Genomic profiles of cancer patients such as gene expression have become a major source to predict responses to drugs in the era of personalized medicine. As large-scale drug screening data with cancer cell lines are available, a number of computational methods have been developed for drug response prediction. However, few methods incorporate both gene expression data and the biological network, which can harbor essential information about the underlying process of the drug response. We proposed an analysis framework called DrugGCN for prediction of Drug response using a Graph Convolutional Network (GCN). DrugGCN first generates a gene graph by combining a Protein-Protein Interaction (PPI) network and gene expression data with feature selection of drug-related genes, and the GCN model detects the local features such as subnetworks of genes that contribute to the drug response by localized filtering. We demonstrated the effectiveness of DrugGCN using biological data showing its high prediction accuracy among the competing methods.

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.

Sign in / Sign up

Export Citation Format

Share Document