Phiclust: a clusterability measure for single-cell transcriptomics reveals phenotypic subpopulations

The ability to discover new cell phenotypes by unsupervised clustering of single-cell transcriptomes has revolutionized biology. Currently, there is no principled way to decide whether a cluster of cells contains meaningful subpopulations that should be further resolved. Here, we present phiclust (ϕclust), a clusterability measure derived from random matrix theory that can be used to identify cell clusters with non-random substructure, testably leading to the discovery of previously overlooked phenotypes.

Genomic alterations and evolution of cell clusters in metastatic invasive micropapillary carcinoma of the breast

Invasive micropapillary carcinoma (IMPC) has very high rates of lymphovascular invasion and lymph node metastasis and has been reported in several organs. However, the genomic mechanisms underlying its metastasis are unclear. Here, we perform whole-genome sequencing of tumor cell clusters from primary IMPC and paired axillary lymph node metastases. Cell clusters in multiple lymph node foci arise from a single subclone of the primary tumor. We find evidence that the monoclonal metastatic ancestor in primary IMPC shares high frequency copy-number loss of PRDM16 and IGSF9 and the copy number gain of ALDH2. Immunohistochemistry analysis further shows that low expression of IGSF9 and PRDM16 and high expression of ALDH2 are associated with lymph node metastasis and poor survival of patients with IMPC. We expect these genomic and evolutionary profiles to contribute to the accurate diagnosis of IMPC.

A novel graph-based k-partitioning approach improves the detection of gene-gene correlations by single-cell RNA sequencing

Background Gene expression is regulated by transcription factors, cofactors, and epigenetic mechanisms. Coexpressed genes indicate similar functional categories and gene networks. Detecting gene-gene coexpression is important for understanding the underlying mechanisms of cellular function and human diseases. A common practice of identifying coexpressed genes is to test the correlation of expression in a set of genes. In single-cell RNA-seq data, an important challenge is the abundance of zero values, so-called “dropout”, which results in biased estimation of gene-gene correlations for downstream analyses. In recent years, efforts have been made to recover coexpressed genes in scRNA-seq data. Here, our goal is to detect coexpressed gene pairs to reduce the “dropout” effect in scRNA-seq data using a novel graph-based k-partitioning method by merging transcriptomically similar cells. Results We observed that the number of zero values was reduced among the merged transcriptomically similar cell clusters. Motivated by this observation, we leveraged a graph-based algorithm and develop an R package, scCorr, to recover the missing gene-gene correlation in scRNA-seq data that enables the reliable acquisition of cluster-based gene-gene correlations in three independent scRNA-seq datasets. The graphically partitioned cell clusters did not change the local cell community. For example, in scRNA-seq data from peripheral blood mononuclear cells (PBMCs), the gene-gene correlation estimated by scCorr outperformed the correlation estimated by the nonclustering method. Among 85 correlated gene pairs in a set of 100 clusters, scCorr detected 71 gene pairs, while the nonclustering method detected only 4 pairs of a dataset from PBMCs. The performance of scCorr was comparable to those of three previously published methods. As an example of downstream analysis using scCorr, we show that scCorr accurately identified a known cell type (i.e., CD4+ T cells) in PBMCs with a receiver operating characteristic area under the curve of 0.96. Conclusions Our results demonstrate that scCorr is a robust and reliable graph-based method for identifying correlated gene pairs, which is fundamental to network construction, gene-gene interaction, and cellular omic analyses. scCorr can be quickly and easily implemented to minimize zero values in scRNA-seq analysis and is freely available at https://github.com/CBIIT-CGBB/scCorr.

Social evolution of shared biofilm matrix components

Biofilm formation is an important and ubiquitous mode of growth among bacteria. Central to the evolutionary advantage of biofilm formation is cell-cell and cell-surface adhesion achieved by a variety of factors, some of which are diffusible compounds that may operate as classical public goods - factors that are costly to produce but may benefit other cells. An outstanding question is how diffusible matrix production, in general, can be stable over evolutionary timescales. In this work, using Vibrio cholerae as a model, we show that shared diffusible biofilm matrix proteins are indeed susceptible to cheater exploitation, and that the evolutionary stability of producing these matrix components fundamentally depends on biofilm spatial structure, intrinsic sharing mechanisms of these components, and flow conditions in the environment. We further show that exploitation of diffusible adhesion proteins is localized within a well-defined spatial range around cell clusters that produce them. Based on this exploitation range and the spatial distribution of cell clusters, we construct a model of costly diffusible matrix production and relate these length scales to the relatedness coefficient in social evolution theory. Our results show that production of diffusible biofilm matrix components is evolutionarily stable under conditions consistent with natural biofilm habitats and host environments. We expect the mechanisms revealed in this study to be relevant to other secreted factors that operate as cooperative public goods in bacterial communities, and the concept of exploitation range and the associated analysis tools to be generally applicable.