More Like this

1. A unifying framework for joint trait analysis under a non-infinitesimal model

   https://doi.org/10.1093/bioinformatics/bty254  

   Johnson, Ruth ; Shi, Huwenbo ; Pasaniuc, Bogdan ; Sankararaman, Sriram ( June 2018 , Bioinformatics)

   A large proportion of risk regions identified by genome-wide association studies (GWAS) are shared across multiple diseases and traits. Understanding whether this clustering is due to sharing of causal variants or chance colocalization can provide insights into shared etiology of complex traits and diseases.

   In this work, we propose a flexible, unifying framework to quantify the overlap between a pair of traits called UNITY (Unifying Non-Infinitesimal Trait analYsis). We formulate a Bayesian generative model that relates the overlap between pairs of traits to GWAS summary statistic data under a non-infinitesimal genetic architecture underlying each trait. We propose a Metropolis–Hastings sampler to compute the posterior density of the genetic overlap parameters in this model. We validate our method through comprehensive simulations and analyze summary statistics from height and body mass index GWAS to show that it produces estimates consistent with the known genetic makeup of both traits.

   The UNITY software is made freely available to the research community at: https://github.com/bogdanlab/UNITY.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
2. ACPred-FL: a sequence-based predictor using effective feature representation to improve the prediction of anti-cancer peptides

   https://doi.org/10.1093/bioinformatics/bty451  

   Wei, Leyi ; Zhou, Chen ; Chen, Huangrong ; Song, Jiangning ; Su, Ran ; Hancock, ed., John ( June 2018 , Bioinformatics)

   Anti-cancer peptides (ACPs) have recently emerged as promising therapeutic agents for cancer treatment. Due to the avalanche of protein sequence data in the post-genomic era, there is an urgent need to develop automated computational methods to enable fast and accurate identification of novel ACPs within the vast number of candidate proteins and peptides.

   To address this, we propose a novel predictor named Anti-Cancer peptide Predictor with Feature representation Learning (ACPred-FL) for accurate prediction of ACPs based on sequence information. More specifically, we develop an effective feature representation learning model, with which we can extract and learn a set of informative features from a pool of support vector machine-based models trained using sequence-based feature descriptors. By doing so, the class label information of data samples is fully utilized. To improve the feature representation, we further employ a two-step feature selection technique, resulting in a most informative five-dimensional feature vector for the final peptide representation. Experimental results show that such five features provide the most discriminative power for identifying ACPs than currently available feature descriptors, highlighting the effectiveness of the proposed feature representation learning approach. The developed ACPred-FL method significantly outperforms state-of-the-art methods.

   The web-server of ACPred-FL is available at http://server.malab.cn/ACPred-FL.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
3. QUBIC2: a novel and robust biclustering algorithm for analyses and interpretation of large-scale RNA-Seq data

   https://doi.org/10.1093/bioinformatics/btz692  

   Xie, Juan ; Ma, Anjun ; Zhang, Yu ; Liu, Bingqiang ; Cao, Sha ; Wang, Cankun ; Xu, Jennifer ; Zhang, Chi ; Ma, Qin ; Birol, ed., Inanc ( September 2019 , Bioinformatics)

   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.

   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.

   The source code of QUBIC2 is freely available at https://github.com/OSU-BMBL/QUBIC2.

   czhang87@iu.edu or qin.ma@osumc.edu

   Supplementary data are available at Bioinformatics online.

    

   more » « less
4. Vivarium: an interface and engine for integrative multiscale modeling in computational biology

   https://doi.org/10.1093/bioinformatics/btac049  

   Agmon, Eran ; Spangler, Ryan K. ; Skalnik, Christopher J. ; Poole, William ; Peirce, Shayn M. ; Morrison, Jerry H. ; Covert, Markus W. ; Valencia, ed., Alfonso ( February 2022 , Bioinformatics)

   This article introduces Vivarium—software born of the idea that it should be as easy as possible for computational biologists to define any imaginable mechanistic model, combine it with existing models and execute them together as an integrated multiscale model. Integrative multiscale modeling confronts the complexity of biology by combining heterogeneous datasets and diverse modeling strategies into unified representations. These integrated models are then run to simulate how the hypothesized mechanisms operate as a whole. But building such models has been a labor-intensive process that requires many contributors, and they are still primarily developed on a case-by-case basis with each project starting anew. New software tools that streamline the integrative modeling effort and facilitate collaboration are therefore essential for future computational biologists.

   Vivarium is a software tool for building integrative multiscale models. It provides an interface that makes individual models into modules that can be wired together in large composite models, parallelized across multiple CPUs and run with Vivarium’s discrete-event simulation engine. Vivarium’s utility is demonstrated by building composite models that combine several modeling frameworks: agent-based models, ordinary differential equations, stochastic reaction systems, constraint-based models, solid-body physics and spatial diffusion. This demonstrates just the beginning of what is possible—Vivarium will be able to support future efforts that integrate many more types of models and at many more biological scales.

   The specific models, simulation pipelines and notebooks developed for this article are all available at the vivarium-notebooks repository: https://github.com/vivarium-collective/vivarium-notebooks. Vivarium-core is available at https://github.com/vivarium-collective/vivarium-core, and has been released on Python Package Index. The Vivarium Collective (https://vivarium-collective.github.io) is a repository of freely available Vivarium processes and composites, including the processes used in Section 3. Supplementary Materials provide with an extensive methodology section, with several code listings that demonstrate the basic interfaces.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
5. ACTIVA : realistic single-cell RNA-seq generation with automatic cell-type identification using introspective variational autoencoders

   https://doi.org/10.1093/bioinformatics/btac095  

   Heydari, A. Ali ; Davalos, Oscar A. ; Zhao, Lihong ; Hoyer, Katrina K. ; Sindi, Suzanne S. ; Mathelier, ed., Anthony ( February 2022 , Bioinformatics)

   Single-cell RNA sequencing (scRNAseq) technologies allow for measurements of gene expression at a single-cell resolution. This provides researchers with a tremendous advantage for detecting heterogeneity, delineating cellular maps or identifying rare subpopulations. However, a critical complication remains: the low number of single-cell observations due to limitations by rarity of subpopulation, tissue degradation or cost. This absence of sufficient data may cause inaccuracy or irreproducibility of downstream analysis. In this work, we present Automated Cell-Type-informed Introspective Variational Autoencoder (ACTIVA): a novel framework for generating realistic synthetic data using a single-stream adversarial variational autoencoder conditioned with cell-type information. Within a single framework, ACTIVA can enlarge existing datasets and generate specific subpopulations on demand, as opposed to two separate models [such as single-cell GAN (scGAN) and conditional scGAN (cscGAN)]. Data generation and augmentation with ACTIVA can enhance scRNAseq pipelines and analysis, such as benchmarking new algorithms, studying the accuracy of classifiers and detecting marker genes. ACTIVA will facilitate analysis of smaller datasets, potentially reducing the number of patients and animals necessary in initial studies.

   We train and evaluate models on multiple public scRNAseq datasets. In comparison to GAN-based models (scGAN and cscGAN), we demonstrate that ACTIVA generates cells that are more realistic and harder for classifiers to identify as synthetic which also have better pair-wise correlation between genes. Data augmentation with ACTIVA significantly improves classification of rare subtypes (more than 45% improvement compared with not augmenting and 4% better than cscGAN) all while reducing run-time by an order of magnitude in comparison to both models.

   The codes and datasets are hosted on Zenodo (https://doi.org/10.5281/zenodo.5879639). Tutorials are available at https://github.com/SindiLab/ACTIVA.

   Supplementary data are available at Bioinformatics online.

    

   more » « less