More Like this

1. MIDAS2: Metagenomic Intra-species Diversity Analysis System

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

   Zhao, Chunyu ; Dimitrov, Boris ; Goldman, Miriam ; Nayfach, Stephen ; Pollard, Katherine S. ; Schwartz, ed., Russell ( November 2022 , Bioinformatics)

   The Metagenomic Intra-Species Diversity Analysis System (MIDAS) is a scalable metagenomic pipeline that identifies single nucleotide variants (SNVs) and gene copy number variants in microbial populations. Here, we present MIDAS2, which addresses the computational challenges presented by increasingly large reference genome databases, while adding functionality for building custom databases and leveraging paired-end reads to improve SNV accuracy. This fast and scalable reengineering of the MIDAS pipeline enables thousands of metagenomic samples to be efficiently genotyped.

   The source code is available at https://github.com/czbiohub/MIDAS2.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
2. High-performance single-cell gene regulatory network inference at scale: the Inferelator 3.0

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

   Skok Gibbs, Claudia ; Jackson, Christopher A. ; Saldi, Giuseppe-Antonio ; Tjärnberg, Andreas ; Shah, Aashna ; Watters, Aaron ; De Veaux, Nicholas ; Tchourine, Konstantine ; Yi, Ren ; Hamamsy, Tymor ; et al ( February 2022 , Bioinformatics)

   Gene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above.

   In this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type-specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data.

   The inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/).

   Supplementary data are available at Bioinformatics online.

    

   more » « less
3. CellMeSH: probabilistic cell-type identification using indexed literature

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

   Mao, Shunfu ; Zhang, Yue ; Seelig, Georg ; Kannan, Sreeram ( February 2022 , Bioinformatics)

   Birol, Inanc (Ed.)

   Single-cell RNA sequencing (scRNA-seq) is widely used for analyzing gene expression in multi-cellular systems and provides unprecedented access to cellular heterogeneity. scRNA-seq experiments aim to identify and quantify all cell types present in a sample. Measured single-cell transcriptomes are grouped by similarity and the resulting clusters are mapped to cell types based on cluster-specific gene expression patterns. While the process of generating clusters has become largely automated, annotation remains a laborious ad hoc effort that requires expert biological knowledge.

   Here, we introduce CellMeSH—a new automated approach to identifying cell types for clusters based on prior literature. CellMeSH combines a database of gene–cell-type associations with a probabilistic method for database querying. The database is constructed by automatically linking gene and cell-type information from millions of publications using existing indexed literature resources. Compared to manually constructed databases, CellMeSH is more comprehensive and is easily updated with new data. The probabilistic query method enables reliable information retrieval even though the gene–cell-type associations extracted from the literature are noisy. CellMeSH is also able to optionally utilize prior knowledge about tissues or cells for further annotation improvement. CellMeSH achieves top-one and top-three accuracies on a number of mouse and human datasets that are consistently better than existing approaches.

   Web server at https://uncurl.cs.washington.edu/db_query and API at https://github.com/shunfumao/cellmesh.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
4. Scalable preprocessing for sparse scRNA-seq data exploiting prior knowledge

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

   Mukherjee, Sumit ; Zhang, Yue ; Fan, Joshua ; Seelig, Georg ; Kannan, Sreeram ( June 2018 , Bioinformatics)

   Single cell RNA-seq (scRNA-seq) data contains a wealth of information which has to be inferred computationally from the observed sequencing reads. As the ability to sequence more cells improves rapidly, existing computational tools suffer from three problems. (i) The decreased reads-per-cell implies a highly sparse sample of the true cellular transcriptome. (ii) Many tools simply cannot handle the size of the resulting datasets. (iii) Prior biological knowledge such as bulk RNA-seq information of certain cell types or qualitative marker information is not taken into account. Here we present UNCURL, a preprocessing framework based on non-negative matrix factorization for scRNA-seq data, that is able to handle varying sampling distributions, scales to very large cell numbers and can incorporate prior knowledge.

   We find that preprocessing using UNCURL consistently improves performance of commonly used scRNA-seq tools for clustering, visualization and lineage estimation, both in the absence and presence of prior knowledge. Finally we demonstrate that UNCURL is extremely scalable and parallelizable, and runs faster than other methods on a scRNA-seq dataset containing 1.3 million cells.

   Source code is available at https://github.com/yjzhang/uncurl_python.

   Supplementary data are available at Bioinformatics online.

    

   more » « less
5. The Tomato Expression Atlas

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

   Fernandez-Pozo, Noe ; Zheng, Yi ; Snyder, Stephen I. ; Nicolas, Philippe ; Shinozaki, Yoshihito ; Fei, Zhangjun ; Catala, Carmen ; Giovannoni, James J. ; Rose, Jocelyn K. C. ; Mueller, Lukas A. ; et al ( March 2017 , Bioinformatics)

   With the development of new high-throughput DNA sequencing technologies and decreasing costs, large gene expression datasets are being generated at an accelerating rate, but can be complex to visualize. New, more interactive and intuitive tools are needed to visualize the spatiotemporal context of expression data and help elucidate gene function. Using tomato fruit as a model, we have developed the Tomato Expression Atlas to facilitate effective data analysis, allowing the simultaneous visualization of groups of genes at a cell/tissue level of resolution within an organ, enhancing hypothesis development and testing in addition to candidate gene identification. This atlas can be adapted to different types of expression data from diverse multicellular species.

   The Tomato Expression Atlas is available at http://tea.solgenomics.net/. Source code is available at https://github.com/solgenomics/Tea.

   Supplementary data are available at Bioinformatics online.

    

   more » « less