Combinatorial prediction of marker panels from single-cell transcriptomic data

Programs, origins and immunomodulatory functions of myeloid cells in glioma

Gliomas are incurable malignancies notable for having an immunosuppressive microenvironment with abundant myeloid cells, the immunomodulatory phenotypes of which remain poorly defined1. Here we systematically investigate these phenotypes by integrating single-cell RNA sequencing, chromatin accessibility, spatial transcriptomics and glioma organoid explant systems. We discovered four immunomodulatory expression programs: microglial inflammatory and scavenger immunosuppressive programs, which are both unique to primary brain tumours, and systemic inflammatory and complement immunosuppressive programs, which are also expressed by non-brain tumours. The programs are not contingent on myeloid cell type, developmental origin or tumour mutational state, but instead are driven by microenvironmental cues, including tumour hypoxia, interleukin-1β, TGFβ and standard-of-care dexamethasone treatment. Their relative expression can predict immunotherapy response and overall survival. By associating the respective programs with mediating genomic elements, transcription factors and signalling pathways, we uncover strategies for manipulating myeloid-cell phenotypes. Our study provides a framework to understand immunomodulation by myeloid cells in glioma and a foundation for the development of more-effective immunotherapies.

snRNA-seq of long-preserved FFPE samples from colorectal liver metastasis lesions with diverse prognoses

Differences in prognostic outcomes are prevalent in patients with colorectal cancer liver metastases. Comparative analysis of tissue samples, particularly applying single-cell transcriptome sequencing technology, can provide a deeper understanding of potential impacting factors. However, long-term monitoring for prognosis determination necessitates extended preservation of tissue samples using formalin-fixed and paraffin-embedded (FFPE) treatments, which can cause substantial RNA degradation, presenting challenges to single-cell or single-nucleus sequencing. In this study, employing snRandom-seq, a single-nucleus RNA sequencing (snRNA-seq) technology specifically for FFPE samples, we tested multiple lesion samples from 18 distinctive colorectal cancer liver metastasis cases with diverse prognostic outcomes that have been preserved for at least three years (mostly over five years). The process yielded expression data from 82,285 cells. The high-quality snRNA-seq data demonstrate the feasibility of single-nucleus sequencing in long-term preserved FFPE samples, offering potential insights into the heterogeneity between different prognoses of colorectal cancer liver metastases, and the relationship between the heterogeneity within different lesions of the same patient and prognosis.

Probe set selection for targeted spatial transcriptomics

Targeted spatial transcriptomic methods capture the topology of cell types and states in tissues at single-cell and subcellular resolution by measuring the expression of a predefined set of genes. The selection of an optimal set of probed genes is crucial for capturing the spatial signals present in a tissue. This requires selecting the most informative, yet minimal, set of genes to profile (gene set selection) for which it is possible to build probes (probe design). However, current selections often rely on marker genes, precluding them from detecting continuous spatial signals or new states. We present Spapros, an end-to-end probe set selection pipeline that optimizes both gene set specificity for cell type identification and within-cell type expression variation to resolve spatially distinct populations while considering prior knowledge as well as probe design and expression constraints. We evaluated Spapros and show that it outperforms other selection approaches in both cell type recovery and recovering expression variation beyond cell types. Furthermore, we used Spapros to design a single-cell resolution in situ hybridization on tissues (SCRINSHOT) experiment of adult lung tissue to demonstrate how probes selected with Spapros identify cell types of interest and detect spatial variation even within cell types.

Hierarchical marker genes selection in scRNA-seq analysis

When analyzing scRNA-seq data containing heterogeneous cell populations, an important task is to select informative marker genes to distinguish various cell clusters and annotate the clusters with biologically meaningful cell types. In existing analysis methods and pipelines, marker genes are typically identified using a one-vs-all strategy, examining differential expression between one cell cluster versus the combination of all other cell clusters. However, this strategy applied to cell clusters belonging to closely related cell types often generates overlapping marker genes, which capture the common signature of closely related cell clusters but provide limited information for distinguishing them. To address the limitations of the one-vs-all strategy, we propose a hierarchical marker gene selection strategy that groups similar cell clusters and selects marker genes in a hierarchical manner. This strategy is able to improve the accuracy and interpretability of cell type identification in single-cell RNA-seq data.

Computational identification of surface markers for isolating distinct subpopulations from heterogeneous cancer cell populations

Intratumor heterogeneity reduces treatment efficacy and complicates our understanding of tumor progression and there is a pressing need to understand the functions of heterogeneous tumor cell subpopulations within a tumor, yet systems to study these processes in vitro are limited. Single-cell RNA sequencing (scRNA-seq) has revealed that some cancer cell lines include distinct subpopulations. Here, we present clusterCleaver, a computational package that uses metrics of statistical distance to identify candidate surface markers maximally unique to transcriptomic subpopulations in scRNA-seq which may be used for FACS isolation. With clusterCleaver, ESAM and BST2/tetherin were experimentally validated as surface markers which identify and separate major transcriptomic subpopulations within MDA-MB-231 and MDA-MB-436 cells, respectively. clusterCleaver is a computationally efficient and experimentally validated workflow for identification of surface markers for tracking and isolating transcriptomically distinct subpopulations within cell lines. This tool paves the way for studies on coexisting cancer cell subpopulations in well-defined in vitro systems.

scPanel: a tool for automatic identification of sparse gene panels for generalizable patient classification using scRNA-seq datasets

Single-cell RNA sequencing (scRNA-seq) technologies can generate transcriptomic profiles at a single-cell resolution in large patient cohorts, facilitating discovery of gene and cellular biomarkers for disease. Yet, when the number of biomarker genes is large, the translation to clinical applications is challenging due to prohibitive sequencing costs. Here, we introduce scPanel, a computational framework designed to bridge the gap between biomarker discovery and clinical application by identifying a sparse gene panel for patient classification from the cell population(s) most responsive to perturbations (e.g. diseases/drugs). scPanel incorporates a data-driven way to automatically determine a minimal number of informative biomarker genes. Patient-level classification is achieved by aggregating the prediction probabilities of cells associated with a patient using the area under the curve score. Application of scPanel to scleroderma, colorectal cancer, and COVID-19 datasets resulted in high patient classification accuracy using only a small number of genes (<20), automatically selected from the entire transcriptome. In the COVID-19 case study, we demonstrated cross-dataset generalizability in predicting disease state in an external patient cohort. scPanel outperforms other state-of-the-art gene selection methods for patient classification and can be used to identify parsimonious sets of reliable biomarker candidates for clinical translation.

An intestinal TH17 cell-derived subset can initiate cancer

Approximately 25% of cancers are preceded by chronic inflammation that occurs at the site of tumor development. However, whether this multifactorial oncogenic process, which commonly occurs in the intestines, can be initiated by a specific immune cell population is unclear. Here, we show that an intestinal T cell subset, derived from interleukin-17 (IL-17)-producing helper T (TH17) cells, induces the spontaneous transformation of the intestinal epithelium. This subset produces inflammatory cytokines, and its tumorigenic potential is not dependent on IL-17 production but on the transcription factors KLF6 and T-BET and interferon-γ. The development of this cell type is inhibited by transforming growth factor-β1 (TGFβ1) produced by intestinal epithelial cells. TGFβ signaling acts on the pretumorigenic TH17 cell subset, preventing its progression to the tumorigenic stage by inhibiting KLF6-dependent T-BET expression. This study therefore identifies an intestinal T cell subset initiating cancer.

CellCover Captures Neural Stem Cell Progression in Mammalian Neocortical Development

Definition of cell classes across the tissues of living organisms is central in the analysis of growing atlases of single-cell RNA sequencing (scRNA-seq) data across biomedicine. Marker genes for cell classes are most often defined by differential expression (DE) methods that serially assess individual genes across landscapes of diverse cells. This serial approach has been extremely useful, but is limited because it ignores possible redundancy or complementarity across genes that can only be captured by analyzing multiple genes simultaneously. We aim to identify discriminating panels of genes. To efficiently explore the vast space of possible marker panels, leverage the large number of cells often sequenced, and overcome zero-inflation in scRNA-seq data, we propose viewing gene panel selection as a variation of the "minimal set-covering problem" in combinatorial optimization. We show that this new method, CellCover, captures cell-class-specific signals in the developing mouse neocortex that are distinct from those defined by DE methods. Transfer learning experiments across mouse, primate, and human data demonstrate that CellCover identifies markers of conserved cell classes in neurogenesis, as well as temporal progression in both progenitors and neurons. Exploring markers of human outer radial glia (oRG, or basal RG) across mammals, we show that transcriptomic elements of this key cell type in the expansion of the human cortex appeared in gliogenic precursors of the rodent before the full program emerged in the primate lineage. We have assembled the public datasets we use in this report at NeMO analytics where the expression of individual genes {NeMO Individual Genes} and marker gene panels can be freely explored {NeMO: Telley 3 Sets Covering Panels}, {NeMO: Telley 12 Sets Covering Panels}, and {NeMO: Sorted Brain Cell Covering Panels}. CellCover is available in {CellCover R} and {CellCover Python}.

CellMarkerPipe: cell marker identification and evaluation pipeline in single cell transcriptomes

Assessing marker genes from all cell clusters can be time-consuming and lack systematic strategy. Streamlining this process through a unified computational platform that automates identification and benchmarking will greatly enhance efficiency and ensure a fair evaluation. We therefore developed a novel computational platform, cellMarkerPipe ( https://github.com/yao-laboratory/cellMarkerPipe ), for automated cell-type specific marker gene identification from scRNA-seq data, coupled with comprehensive evaluation schema. CellMarkerPipe adaptively wraps around a collection of commonly used and state-of-the-art tools, including Seurat, COSG, SC3, SCMarker, COMET, and scGeneFit. From rigorously testing across diverse samples, we ascertain SCMarker's overall reliable performance in single marker gene selection, with COSG showing commendable speed and comparable efficacy. Furthermore, we demonstrate the pivotal role of our approach in real-world medical datasets. This general and opensource pipeline stands as a significant advancement in streamlining cell marker gene identification and evaluation, fitting broad applications in the field of cellular biology and medical research.

Computational identification of surface markers for isolating distinct subpopulations from heterogeneous cancer cell populations

Intratumor heterogeneity reduces treatment efficacy and complicates our understanding of tumor progression. There is a pressing need to understand the functions of heterogeneous tumor cell subpopulations within a tumor, yet biological systems to study these processes in vitro are limited. With the advent of single-cell RNA sequencing (scRNA-seq), it has become clear that some cancer cell line models include distinct subpopulations. Heterogeneous cell lines offer a unique opportunity to study the dynamics and evolution of genetically similar cancer cell subpopulations in controlled experimental settings. Here, we present clusterCleaver, a computational package that uses metrics of statistical distance to identify candidate surface markers maximally unique to transcriptomic subpopulations in scRNA-seq which may be used for FACS isolation. clusterCleaver was experimentally validated using the MDA-MB-231 and MDA-MB-436 breast cancer cell lines. ESAM and BST2/tetherin were experimentally confirmed as surface markers which identify and separate major transcriptomic subpopulations within MDA-MB-231 and MDA-MB-436 cells, respectively. clusterCleaver is a computationally efficient and experimentally validated workflow for identification and enrichment of distinct subpopulations within cell lines which paves the way for studies on the coexistence of cancer cell subpopulations in well-defined in vitro systems.

Direct comparison of mass cytometry and single-cell RNA sequencing of human peripheral blood mononuclear cells

Single-cell methods offer a high-resolution approach for characterizing cell populations. Many studies rely on single-cell transcriptomics to draw conclusions regarding cell state and behavior, with the underlying assumption that transcriptomic readouts largely parallel their protein counterparts and subsequent activity. However, the relationship between transcriptomic and proteomic measurements is imprecise, and thus datasets that probe the extent of their concordance will be useful to refine such conclusions. Additionally, novel single-cell analysis tools often lack appropriate gold standard datasets for the purposes of assessment. Integrative (combining the two data modalities) and predictive (using one modality to improve results from the other) approaches in particular, would benefit from transcriptomic and proteomic data from the same sample of cells. For these reasons, we performed single-cell RNA sequencing, mass cytometry, and flow cytometry on a split-sample of human peripheral blood mononuclear cells. We directly compare the proportions of specific cell types resolved by each technique, and further describe the extent to which protein and mRNA measurements correlate within distinct cell types.

Spatial transcriptomic brain imaging reveals the effects of immunomodulation therapy on specific regional brain cells in a mouse dementia model

Increasing evidence of brain-immune crosstalk raises expectations for the efficacy of novel immunotherapies in Alzheimer's disease (AD), but the lack of methods to examine brain tissues makes it difficult to evaluate therapeutics. Here, we investigated the changes in spatial transcriptomic signatures and brain cell types using the 10x Genomics Visium platform in immune-modulated AD models after various treatments. To proceed with an analysis suitable for barcode-based spatial transcriptomics, we first organized a workflow for segmentation of neuroanatomical regions, establishment of appropriate gene combinations, and comprehensive review of altered brain cell signatures. Ultimately, we investigated spatial transcriptomic changes following administration of immunomodulators, NK cell supplements and an anti-CD4 antibody, which ameliorated behavior impairment, and designated brain cells and regions showing probable associations with behavior changes. We provided the customized analytic pipeline into an application named STquantool. Thus, we anticipate that our approach can help researchers interpret the real action of drug candidates by simultaneously investigating the dynamics of all transcripts for the development of novel AD therapeutics.

Highly sensitive spatial transcriptomics using FISHnCHIPs of multiple co-expressed genes

High-dimensional, spatially resolved analysis of intact tissue samples promises to transform biomedical research and diagnostics, but existing spatial omics technologies are costly and labor-intensive. We present Fluorescence In Situ Hybridization of Cellular HeterogeneIty and gene expression Programs (FISHnCHIPs) for highly sensitive in situ profiling of cell types and gene expression programs. FISHnCHIPs achieves this by simultaneously imaging ~2-35 co-expressed genes (clustered into modules) that are spatially co-localized in tissues, resulting in similar spatial information as single-gene Fluorescence In Situ Hybridization (FISH), but with ~2-20-fold higher sensitivity. Using FISHnCHIPs, we image up to 53 modules from the mouse kidney and mouse brain, and demonstrate high-speed, large field-of-view profiling of a whole tissue section. FISHnCHIPs also reveals spatially restricted localizations of cancer-associated fibroblasts in a human colorectal cancer biopsy. Overall, FISHnCHIPs enables fast, robust, and scalable cell typing of tissues with normal physiology or undergoing pathogenesis.

CaSSiDI: novel single-cell "Cluster Similarity Scoring and Distinction Index" reveals critical functions for PirB and context-dependent Cebpb repression

PirB is an inhibitory cell surface receptor particularly prominent on myeloid cells. PirB curtails the phenotypes of activated macrophages during inflammation or tumorigenesis, but its functions in macrophage homeostasis are obscure. To elucidate PirB-related functions in macrophages at steady-state, we generated and compared single-cell RNA-sequencing (scRNAseq) datasets obtained from myeloid cell subsets of wild type (WT) and PirB-deficient knockout (PirB KO) mice. To facilitate this analysis, we developed a novel approach to clustering parameter optimization called "Cluster Similarity Scoring and Distinction Index" (CaSSiDI). We demonstrate that CaSSiDI is an adaptable computational framework that facilitates tandem analysis of two scRNAseq datasets by optimizing clustering parameters. We further show that CaSSiDI offers more advantages than a standard Seurat analysis because it allows direct comparison of two or more independently clustered datasets, thereby alleviating the need for batch-correction while identifying the most similar and different clusters. Using CaSSiDI, we found that PirB is a novel regulator of Cebpb expression that controls the generation of Ly6Clo patrolling monocytes and the expansion properties of peritoneal macrophages. PirB's effect on Cebpb is tissue-specific since it was not observed in splenic red pulp macrophages (RPMs). However, CaSSiDI revealed a segregation of the WT RPM population into a CD68loIrf8+ "neuronal-primed" subset and an CD68hiFtl1+ "iron-loaded" subset. Our results establish the utility of CaSSiDI for single-cell assay analyses and the determination of optimal clustering parameters. Our application of CaSSiDI in this study has revealed previously unknown roles for PirB in myeloid cell populations. In particular, we have discovered homeostatic functions for PirB that are related to Cebpb expression in distinct macrophage subsets.

Distinct ontogenetic lineages dictate cDC2 heterogeneity

Conventional dendritic cells (cDCs) include functionally and phenotypically diverse populations, such as cDC1s and cDC2s. The latter population has been variously subdivided into Notch-dependent cDC2s, KLF4-dependent cDC2s, T-bet+ cDC2As and T-bet- cDC2Bs, but it is unclear how all these subtypes are interrelated and to what degree they represent cell states or cell subsets. All cDCs are derived from bone marrow progenitors called pre-cDCs, which circulate through the blood to colonize peripheral tissues. Here, we identified distinct mouse pre-cDC2 subsets biased to give rise to cDC2As or cDC2Bs. We showed that a Siglec-H+ pre-cDC2A population in the bone marrow preferentially gave rise to Siglec-H- CD8α+ pre-cDC2As in tissues, which differentiated into T-bet+ cDC2As. In contrast, a Siglec-H- fraction of pre-cDCs in the bone marrow and periphery mostly generated T-bet- cDC2Bs, a lineage marked by the expression of LysM. Our results showed that cDC2A versus cDC2B fate specification starts in the bone marrow and suggest that cDC2 subsets are ontogenetically determined lineages, rather than cell states imposed by the peripheral tissue environment.

A comparison of marker gene selection methods for single-cell RNA sequencing data

BACKGROUND

The development of single-cell RNA sequencing (scRNA-seq) has enabled scientists to catalog and probe the transcriptional heterogeneity of individual cells in unprecedented detail. A common step in the analysis of scRNA-seq data is the selection of so-called marker genes, most commonly to enable annotation of the biological cell types present in the sample. In this paper, we benchmark 59 computational methods for selecting marker genes in scRNA-seq data.

RESULTS

We compare the performance of the methods using 14 real scRNA-seq datasets and over 170 additional simulated datasets. Methods are compared on their ability to recover simulated and expert-annotated marker genes, the predictive performance and characteristics of the gene sets they select, their memory usage and speed, and their implementation quality. In addition, various case studies are used to scrutinize the most commonly used methods, highlighting issues and inconsistencies.

CONCLUSIONS

Overall, we present a comprehensive evaluation of methods for selecting marker genes in scRNA-seq data. Our results highlight the efficacy of simple methods, especially the Wilcoxon rank-sum test, Student's t-test, and logistic regression.

MarkerMap: nonlinear marker selection for single-cell studies

Single-cell RNA-seq data allow the quantification of cell type differences across a growing set of biological contexts. However, pinpointing a small subset of genomic features explaining this variability can be ill-defined and computationally intractable. Here we introduce MarkerMap, a generative model for selecting minimal gene sets which are maximally informative of cell type origin and enable whole transcriptome reconstruction. MarkerMap provides a scalable framework for both supervised marker selection, aimed at identifying specific cell type populations, and unsupervised marker selection, aimed at gene expression imputation and reconstruction. We benchmark MarkerMap's competitive performance against previously published approaches on real single cell gene expression data sets. MarkerMap is available as a pip installable package, as a community resource aimed at developing explainable machine learning techniques for enhancing interpretability in single-cell studies.

Leukaemia exposure alters the transcriptional profile and function of BCR::ABL1 negative macrophages in the bone marrow niche

Macrophages are fundamental cells of the innate immune system that support normal haematopoiesis and play roles in both anti-cancer immunity and tumour progression. Here we use a chimeric mouse model of chronic myeloid leukaemia (CML) and human bone marrow (BM) derived macrophages to study the impact of the dysregulated BM microenvironment on bystander macrophages. Utilising single-cell RNA sequencing (scRNA-seq) of Philadelphia chromosome (Ph) negative macrophages we reveal unique subpopulations of immature macrophages residing in the CML BM microenvironment. CML exposed macrophages separate from their normal counterparts by reduced expression of the surface marker CD36, which significantly reduces clearance of apoptotic cells. We uncover aberrant production of CML-secreted factors, including the immune modulatory protein lactotransferrin (LTF), that suppresses efferocytosis, phagocytosis, and CD36 surface expression in BM macrophages, indicating that the elevated secretion of LTF is, at least partially responsible for the supressed clearance function of Ph- macrophages.

A single cell RNAseq benchmark experiment embedding "controlled" cancer heterogeneity

Single-cell RNA sequencing (scRNA-seq) has emerged as a vital tool in tumour research, enabling the exploration of molecular complexities at the individual cell level. It offers new technical possibilities for advancing tumour research with the potential to yield significant breakthroughs. However, deciphering meaningful insights from scRNA-seq data poses challenges, particularly in cell annotation and tumour subpopulation identification. Efficient algorithms are therefore needed to unravel the intricate biological processes of cancer. To address these challenges, benchmarking datasets are essential to validate bioinformatics methodologies for analysing single-cell omics in oncology. Here, we present a 10XGenomics scRNA-seq experiment, providing a controlled heterogeneous environment using lung cancer cell lines characterised by the expression of seven different driver genes (EGFR, ALK, MET, ERBB2, KRAS, BRAF, ROS1), leading to partially overlapping functional pathways. Our dataset provides a comprehensive framework for the development and validation of methodologies for analysing cancer heterogeneity by means of scRNA-seq.

cellMarkerPipe: Cell Marker Identification and Evaluation Pipeline in Single Cell Transcriptomes

Assessing marker genes from all cell clusters can be time-consuming and lack systematic strategy. Streamlining this process through a unified computational platform that automates identification and benchmarking will greatly enhance efficiency and ensure a fair evaluation. We therefore developed a novel computational platform, cellMarkerPipe (https://github.com/yao-laboratory/cellMarkerPipe), for automated cell-type specific marker gene identification from scRNA-seq data, coupled with comprehensive evaluation schema. CellMarkerPipe adaptively wraps around a collection of commonly used and state-of-the-art tools, including Seurat, COSG, SC3, SCMarker, COMET, and scGeneFit. From rigorously testing across diverse samples, we ascertain SCMarker's overall reliable performance in single marker gene selection, with COSG showing commendable speed and comparable efficacy. Furthermore, we demonstrate the pivotal role of our approach in real-world medical datasets. This general and opensource pipeline stands as a significant advancement in streamlining cell marker gene identification and evaluation, fitting broad applications in the field of cellular biology and medical research.

scMD facilitates cell type deconvolution using single-cell DNA methylation references

The proliferation of single-cell RNA-sequencing data has led to the widespread use of cellular deconvolution, aiding the extraction of cell-type-specific information from extensive bulk data. However, those advances have been mostly limited to transcriptomic data. With recent developments in single-cell DNA methylation (scDNAm), there are emerging opportunities for deconvolving bulk DNAm data, particularly for solid tissues like brain that lack cell-type references. Due to technical limitations, current scDNAm sequences represent a small proportion of the whole genome for each single cell, and those detected regions differ across cells. This makes scDNAm data ultra-high dimensional and ultra-sparse. To deal with these challenges, we introduce scMD (single cell Methylation Deconvolution), a cellular deconvolution framework to reliably estimate cell type fractions from tissue-level DNAm data. To analyze large-scale complex scDNAm data, scMD employs a statistical approach to aggregate scDNAm data at the cell cluster level, identify cell-type marker DNAm sites, and create precise cell-type signature matrixes that surpass state-of-the-art sorted-cell or RNA-derived references. Through thorough benchmarking in several datasets, we demonstrate scMD's superior performance in estimating cellular fractions from bulk DNAm data. With scMD-estimated cellular fractions, we identify cell type fractions and cell type-specific differentially methylated cytosines associated with Alzheimer's disease.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.

Intestinal microbiota-specific Th17 cells possess regulatory properties and suppress effector T cells via c-MAF and IL-10

Commensal microbes induce cytokine-producing effector tissue-resident CD4+ T cells, but the function of these T cells in mucosal homeostasis is not well understood. Here, we report that commensal-specific intestinal Th17 cells possess an anti-inflammatory phenotype marked by expression of interleukin (IL)-10 and co-inhibitory receptors. The anti-inflammatory phenotype of gut-resident commensal-specific Th17 cells was driven by the transcription factor c-MAF. IL-10-producing commensal-specific Th17 cells were heterogeneous and derived from a TCF1+ gut-resident progenitor Th17 cell population. Th17 cells acquired IL-10 expression and anti-inflammatory phenotype in the small-intestinal lamina propria. IL-10 production by CD4+ T cells and IL-10 signaling in intestinal macrophages drove IL-10 expression by commensal-specific Th17 cells. Intestinal commensal-specific Th17 cells possessed immunoregulatory functions and curbed effector T cell activity in vitro and in vivo in an IL-10-dependent and c-MAF-dependent manner. Our results suggest that tissue-resident commensal-specific Th17 cells perform regulatory functions in mucosal homeostasis.

MAGNETO: Cell type marker panel generator from single-cell transcriptomic data

Single-cell RNA sequencing experiments produce data useful to identify different cell types, including uncharacterized and rare ones. This enables us to study the specific functional roles of these cells in different microenvironments and contexts. After identifying a (novel) cell type of interest, it is essential to build succinct marker panels, composed of a few genes referring to cell surface proteins and clusters of differentiation molecules, able to discriminate the desired cells from the other cell populations. In this work, we propose a fully-automatic framework called MAGNETO, which can help construct optimal marker panels starting from a single-cell gene expression matrix and a cell type identity for each cell. MAGNETO builds effective marker panels solving a tailored bi-objective optimization problem, where the first objective regards the identification of the genes able to isolate a specific cell type, while the second conflicting objective concerns the minimization of the total number of genes included in the panel. Our results on three public datasets show that MAGNETO can identify marker panels that identify the cell populations of interest better than state-of-the-art approaches. Finally, by fine-tuning MAGNETO, our results demonstrate that it is possible to obtain marker panels with different specificity levels.

Programs, Origins, and Niches of Immunomodulatory Myeloid Cells in Gliomas

Gliomas are incurable malignancies notable for an immunosuppressive microenvironment with abundant myeloid cells whose immunomodulatory properties remain poorly defined. Here, utilizing scRNA-seq data for 183,062 myeloid cells from 85 human tumors, we discover that nearly all glioma-associated myeloid cells express at least one of four immunomodulatory activity programs: Scavenger Immunosuppressive, C1Q Immunosuppressive, CXCR4 Inflammatory, and IL1B Inflammatory. All four programs are present in IDH1 mutant and wild-type gliomas and are expressed in macrophages, monocytes, and microglia whether of blood or resident myeloid cell origins. Integrating our scRNA-seq data with mitochondrial DNA-based lineage tracing, spatial transcriptomics, and organoid explant systems that model peripheral monocyte infiltration, we show that these programs are driven by microenvironmental cues and therapies rather than myeloid cell type, origin, or mutation status. The C1Q Immunosuppressive program is driven by routinely administered dexamethasone. The Scavenger Immunosuppressive program includes ligands with established roles in T-cell suppression, is induced in hypoxic regions, and is associated with immunotherapy resistance. Both immunosuppressive programs are less prevalent in lower-grade gliomas, which are instead enriched for the CXCR4 Inflammatory program. Our study provides a framework to understand immunomodulatory myeloid cells in glioma, and a foundation to develop more effective immunotherapies.

On the use of QDE-SVM for gene feature selection and cell type classification from scRNA-seq data

Cell type identification is one of the fundamental tasks in single-cell RNA sequencing (scRNA-seq) studies. It is a key step to facilitate downstream interpretations such as differential expression, trajectory inference, etc. scRNA-seq data contains technical variations that could affect the interpretation of the cell types. Therefore, gene selection, also known as feature selection in data science, plays an important role in selecting informative genes for scRNA-seq cell type identification. Generally speaking, feature selection methods are categorized into filter-, wrapper-, and embedded-based approaches. From the existing literature, methods from filter- and embedded-based approaches are widely applied in scRNA-seq gene selection tasks. The wrapper-based method that gives promising results in other fields has yet been extensively utilized for selecting gene features from scRNA-seq data; in addition, most of the existing wrapper methods used in this field are clustering instead of classification-based. With a large number of annotated data available today, this study applied a classification-based approach as an alternative to the clustering-based wrapper method. In our work, a quantum-inspired differential evolution (QDE) wrapped with a classification method was introduced to select a subset of genes from twelve well-known scRNA-seq transcriptomic datasets to identify cell types. In particular, the QDE was combined with different machine-learning (ML) classifiers namely logistic regression, decision tree, support vector machine (SVM) with linear and radial basis function kernels, as well as extreme learning machine. The linear SVM wrapped with QDE, namely QDE-SVM, was chosen by referring to the feature selection results from the experiment. QDE-SVM showed a superior cell type classification performance among QDE wrapping with other ML classifiers as well as the recent wrapper methods (i.e., FSCAM, SSD-LAHC, MA-HS, and BSF). QDE-SVM achieved an average accuracy of 0.9559, while the other wrapper methods achieved average accuracies in the range of 0.8292 to 0.8872.

Microfluidics-free single-cell genomics reveals complex central-peripheral immune crosstalk in the mouse brain during peripheral inflammation

Inflammation is a realized detriment to brain health in a growing number of neurological diseases, but querying neuroinflammation in its cellular complexity remains a challenge. This manuscript aims to provide a reliable and accessible strategy for examining the brain's immune system. We compare the efficacy of cell isolation methods in producing ample and pure immune samples from mouse brains. Then, with the high-input single-cell genomics platform PIPseq, we generate a rich neuroimmune dataset containing microglia and many peripheral immune populations. To demonstrate this strategy's utility, we interrogate the well-established model of LPS-induced neuroinflammation with single-cell resolution. We demonstrate the activation of crosstalk between microglia and peripheral phagocytes and highlight the unique contributions of microglia and peripheral immune cells to neuroinflammation. Our approach enables the high-depth evaluation of inflammation in longstanding rodent models of neurological disease to reveal novel insight into the contributions of the immune system to brain health.

Systems level identification of a matrisome-associated macrophage polarisation state in multi-organ fibrosis

Tissue fibrosis affects multiple organs and involves a master-regulatory role of macrophages which respond to an initial inflammatory insult common in all forms of fibrosis. The recently unravelled multi-organ heterogeneity of macrophages in healthy and fibrotic human disease suggests that macrophages expressing osteopontin (SPP1) associate with lung and liver fibrosis. However, the conservation of this SPP1+ macrophage population across different tissues and its specificity to fibrotic diseases with different etiologies remain unclear. Integrating 15 single-cell RNA-sequencing datasets to profile 235,930 tissue macrophages from healthy and fibrotic heart, lung, liver, kidney, skin, and endometrium, we extended the association of SPP1+ macrophages with fibrosis to all these tissues. We also identified a subpopulation expressing matrisome-associated genes (e.g., matrix metalloproteinases and their tissue inhibitors), functionally enriched for ECM remodelling and cell metabolism, representative of a matrisome-associated macrophage (MAM) polarisation state within SPP1+ macrophages. Importantly, the MAM polarisation state follows a differentiation trajectory from SPP1+ macrophages and is associated with a core set of regulon activity. SPP1+ macrophages without the MAM polarisation state (SPP1+MAM-) show a positive association with ageing lung in mice and humans. These results suggest an advanced and conserved polarisation state of SPP1+ macrophages in fibrotic tissues resulting from prolonged inflammatory cues within each tissue microenvironment.

A Roadmap for the Human Gut Cell Atlas

The number of studies investigating the human gastrointestinal tract using various single-cell profiling methods has increased substantially in the past few years. Although this increase provides a unique opportunity for the generation of the first comprehensive Human Gut Cell Atlas (HGCA), there remains a range of major challenges ahead. Above all, the ultimate success will largely depend on a structured and coordinated approach that aligns global efforts undertaken by a large number of research groups. In this Roadmap, we discuss a comprehensive forward-thinking direction for the generation of the HGCA on behalf of the Gut Biological Network of the Human Cell Atlas. Based on the consensus opinion of experts from across the globe, we outline the main requirements for the first complete HGCA by summarizing existing data sets and highlighting anatomical regions and/or tissues with limited coverage. We provide recommendations for future studies and discuss key methodologies and the importance of integrating the healthy gut atlas with related diseases and gut organoids. Importantly, we critically overview the computational tools available and provide recommendations to overcome key challenges.

LIDER: cell embedding based deep neural network classifier for supervised cell type identification

Background

Automatic cell type identification has been an urgent task for the rapid development of single-cell RNA-seq techniques. Generally, the current approach for cell type identification is to generate cell clusters by unsupervised clustering and later assign labels to each cell cluster with manual annotation.

Methods

Here, we introduce LIDER (celL embeddIng based Deep nEural netwoRk classifier), a deep supervised learning method that combines cell embedding and deep neural network classifier for automatic cell type identification. Based on a stacked denoising autoencoder with a tailored and reconstructed loss function, LIDER identifies cell embedding and predicts cell types with a deep neural network classifier. LIDER was developed upon a stacked denoising autoencoder to learn encoder-decoder structures for identifying cell embedding.

Results

LIDER accurately identifies cell types by using stacked denoising autoencoder. Benchmarking against state-of-the-art methods across eight types of single-cell data, LIDER achieves comparable or even superior enhancement performance. Moreover, LIDER suggests comparable robust to batch effects. Our results show a potential in deep supervised learning for automatic cell type identification of single-cell RNA-seq data. The LIDER codes are available at https://github.com/ShiMGLab/LIDER.

scMD: cell type deconvolution using single-cell DNA methylation references

The proliferation of single-cell RNA sequencing data has led to the widespread use of cellular deconvolution, aiding the extraction of cell type-specific information from extensive bulk data. However, those advances have been mostly limited to transcriptomic data. With recent development in single-cell DNA methylation (scDNAm), new avenues have been opened for deconvolving bulk DNAm data, particularly for solid tissues like the brain that lack cell-type references. Due to technical limitations, current scDNAm sequences represent a small proportion of the whole genome for each single cell, and those detected regions differ across cells. This makes scDNAm data ultra-high dimensional and ultra-sparse. To deal with these challenges, we introduce scMD (single cell Methylation Deconvolution), a cellular deconvolution framework to reliably estimate cell type fractions from tissue-level DNAm data. To analyze large-scale complex scDNAm data, scMD employs a statistical approach to aggregate scDNAm data at the cell cluster level, identify cell-type marker DNAm sites, and create a precise cell-type signature matrix that surpasses state-of-the-art sorted-cell or RNA-derived references. Through thorough benchmarking in several datasets, we demonstrate scMD's superior performance in estimating cellular fractions from bulk DNAm data. With scMD-estimated cellular fractions, we identify cell type fractions and cell type-specific differentially methylated cytosines associated with Alzheimer's disease.

The spatiotemporal dynamics of spatially variable genes in developing mouse brain revealed by a novel computational scheme

To understand how brain regions form and work, it is important to explore the spatially variable genes (SVGs) enriched in specific brain regions during development. Spatial transcriptomics techniques provide opportunity to select SVGs in the high-throughput way. However, previous methods neglected the ranking order and combinatorial effect of SVGs, making them difficult to automatically select the high-priority SVGs from spatial transcriptomics data. Here, we proposed a novel computational pipeline, called SVGbit, to rank the individual and combinatorial SVGs for marker selection in various brain regions, which was tested in different kinds of public datasets for both human and mouse brains. We then generated the spatial transcriptomics and immunohistochemistry data from mouse brain at critical embryonic and neonatal stages. The results show that our ranking and clustering scheme captures the key SVGs which coincide with known anatomic regions in the developing mouse brain. More importantly, SVGbit can facilitate the identification of multiple gene combination sets in different brain regions. We identified three dynamical sub-regions which can be segregated by the staining of Sox2 and Calb2 in thalamus, and we also found that Nr4a2 expression gradually segregates the neocortex and hippocampus during the development. In summary, our work not only reveals the spatiotemporal dynamics of individual and combinatorial SVGs in developing mouse brain, but also provides a novel computational pipeline to facilitate the selection of marker genes from spatial transcriptomics data.

B-cell-specific checkpoint molecules that regulate anti-tumour immunity

The role of B cells in anti-tumour immunity is still debated and, accordingly, immunotherapies have focused on targeting T and natural killer cells to inhibit tumour growth1,2. Here, using high-throughput flow cytometry as well as bulk and single-cell RNA-sequencing and B-cell-receptor-sequencing analysis of B cells temporally during B16F10 melanoma growth, we identified a subset of B cells that expands specifically in the draining lymph node over time in tumour-bearing mice. The expanding B cell subset expresses the cell surface molecule T cell immunoglobulin and mucin domain 1 (TIM-1, encoded by Havcr1) and a unique transcriptional signature, including multiple co-inhibitory molecules such as PD-1, TIM-3, TIGIT and LAG-3. Although conditional deletion of these co-inhibitory molecules on B cells had little or no effect on tumour burden, selective deletion of Havcr1 in B cells both substantially inhibited tumour growth and enhanced effector T cell responses. Loss of TIM-1 enhanced the type 1 interferon response in B cells, which augmented B cell activation and increased antigen presentation and co-stimulation, resulting in increased expansion of tumour-specific effector T cells. Our results demonstrate that manipulation of TIM-1-expressing B cells enables engagement of the second arm of adaptive immunity to promote anti-tumour immunity and inhibit tumour growth.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding the pathologies of diseases. However, several experimental and computational challenges remain in developing and implementing transcriptomics-based deconvolution approaches, especially those using a single cell/nuclei RNA-seq reference atlas, which are becoming rapidly available across many tissues. Notably, deconvolution algorithms are frequently developed using samples from tissues with similar cell sizes. However, brain tissue or immune cell populations have cell types with substantially different cell sizes, total mRNA expression, and transcriptional activity. When existing deconvolution approaches are applied to these tissues, these systematic differences in cell sizes and transcriptomic activity confound accurate cell proportion estimates and instead may quantify total mRNA content. Furthermore, there is a lack of standard reference atlases and computational approaches to facilitate integrative analyses, including not only bulk and single cell/nuclei RNA-seq data, but also new data modalities from spatial -omic or imaging approaches. New multi-assay datasets need to be collected with orthogonal data types generated from the same tissue block and the same individual, to serve as a "gold standard" for evaluating new and existing deconvolution methods. Below, we discuss these key challenges and how they can be addressed with the acquisition of new datasets and approaches to analysis.

scMAGS: Marker gene selection from scRNA-seq data for spatial transcriptomics studies

Single-Cell RNA sequencing (scRNA-seq) has provided unprecedented opportunities for exploring gene expression and thus uncovering regulatory relationships between genes at the single-cell level. However, scRNA-seq relies on isolating cells from tissues. Therefore, the spatial context of the regulatory processes is lost. A recent technological innovation, spatial transcriptomics, allows for the measurement of gene expression while preserving spatial information. An initial step in the spatial transcriptomic analysis is to identify the cell type, which requires a careful selection of cell-specific marker genes. For this purpose, currently, scRNA-seq data is used to select a limited number of marker genes from among all genes that distinguish cell types from each other. This study proposes scMAGS (single-cell MArker Gene Selection), a novel method for marker gene selection from scRNA-seq data for spatial transcriptomics studies. scMAGS uses a filtering step in which the candidate genes are identified before the marker gene selection step. For the selection of marker genes, cluster validity indices, the Silhouette index, or the Calinski-Harabasz index (for large datasets) are utilized. Experimental results showed that, in comparison to the existing methods, scMAGS is scalable, fast, and accurate. Even for large datasets with millions of cells, scMAGS could find the required number of marker genes in a reasonable amount of time with fewer memory requirements. scMAGS is made freely available at https://github.com/doganlab/scmags and can be downloaded from the Python Package Directory (PyPI) software repository with the command pip install scmags.

Using "Galaxy-rCASC": A Public Galaxy Instance for Single-Cell RNA-Seq Data Analysis

rCASC is a modular workflow providing an integrated environment for single-cell RNA-seq (scRNA-Seq) data analysis exploiting Docker containers to achieve functional and computational reproducibility. It was initially developed as an R package usable also through a Java GUI. However, the Java frontend cannot be employed when running rCASC on a remote server, a typical setup due to the significant computational resources commonly needed to analyze scRNA-Seq data.To allow the use of rCASC through a graphical user interface on the client side and to harness the many advantages provided by the Galaxy platform, we have made rCASC available as a Galaxy set of tools, also providing a dedicated public instance of Galaxy named "Galaxy-rCASC." To integrate rCASC into Galaxy, all its functions, originally implemented as a set of Docker containers to maximize reproducibility, have been extensively reworked to become independent from the R package functions that launch them in the original implementation. Furthermore, suitable Galaxy wrappers have been developed for most functions of rCASC. We provide a detailed reference document to the use of Galaxy-rCASC with insights and explanations on the platform functionalities, parameters, and output while guiding the reader through the typical rCASC analysis workflow of a scRNA-Seq dataset.

Identifying Gene Markers Associated with Cell Subpopulations

An important point of the analysis of a single-cell RNA experiment is the identification of the key elements, i.e., genes, characterizing each cell subpopulation cluster. In this chapter, we describe the use of sparsely connected autoencoder, as a tool to convert single-cell clusters in pseudo-RNAseq experiments to be used as input for differential expression analysis, and the use of COMET, as a tool to depict cluster-specific gene markers.

LY6D marks pre-existing resistant basosquamous tumor subpopulations

Improved response to canonical therapies requires a mechanistic understanding of dynamic tumor heterogeneity by identifying discrete cellular populations with enhanced cellular plasticity. We have previously demonstrated distinct resistance mechanisms in skin basal cell carcinomas, but a comprehensive understanding of the cellular states and markers associated with these populations remains poorly understood. Here we identify a pre-existing resistant cellular population in naive basal cell carcinoma tumors marked by the surface marker LY6D. LY6D+ tumor cells are spatially localized and possess basal cell carcinoma and squamous cell carcinoma-like features. Using computational tools, organoids, and spatial tools, we show that LY6D+ basosquamous cells represent a persister population lying on a central node along the skin lineage-associated spectrum of epithelial states with local environmental and applied therapies determining the kinetics of accumulation. Surprisingly, LY6D+ basosquamous populations exist in many epithelial tumors, such as pancreatic adenocarcinomas, which have poor outcomes. Overall, our results identify the resistant LY6D+ basosquamous population as an important clinical target and suggest strategies for future therapeutic approaches to target them.

Single-Cell Transcriptomic Census of Endothelial Changes Induced by Matrix Stiffness and the Association with Atherosclerosis

Vascular endothelial cell (EC) plasticity plays a critical role in the progression of atherosclerosis by giving rise to mesenchymal phenotypes in the plaque lesion. Despite the evidence for arterial stiffening as a major contributor to atherosclerosis, the complex interplay among atherogenic stimuli in vivo has hindered attempts to determine the effects of extracellular matrix (ECM) stiffness on endothelial-mesenchymal transition (EndMT). To study the regulatory effects of ECM stiffness on EndMT, an in vitro model is developed in which human coronary artery ECs are cultured on physiological or pathological stiffness substrates. Leveraging single-cell RNA sequencing, cell clusters with mesenchymal transcriptional features are identified to be more prevalent on pathological substrates than physiological substrates. Trajectory inference analyses reveal a novel mesenchymal-to-endothelial reverse transition, which is blocked by pathological stiffness substrates, in addition to the expected EndMT trajectory. ECs pushed to a mesenchymal character by pathological stiffness substrates are enriched in transcriptional signatures of atherosclerotic ECs from human and murine plaques. This study characterizes at single-cell resolution the transcriptional programs that underpin EC plasticity in both physiological or pathological milieus, and thus serves as a valuable resource for more precisely defining EndMT and the transcriptional programs contributing to atherosclerosis.

Progenitor Hierarchy of Chronic Myelomonocytic Leukemia Identifies Inflammatory Monocytic-Biased Trajectory Linked to Worse Outcomes

Myeloblast expansion is a hallmark of disease progression and comprises CD34+ hematopoietic stem and progenitor cells (HSPC). How this compartment evolves during disease progression in chronic myeloid neoplasms is unknown. Using single-cell RNA sequencing and high-parameter flow cytometry, we show that chronic myelomonocytic leukemia (CMML) CD34+ HSPC can be classified into three differentiation trajectories: monocytic, megakaryocyte-erythroid progenitor (MEP), and normal-like. Hallmarks of monocytic-biased trajectory were enrichment of CD120b+ inflammatory granulocyte-macrophage progenitor (GMP)-like cells, activated cytokine receptor signaling, phenotypic hematopoietic stem cell (HSC) depletion, and adverse outcomes. Cytokine receptor diversity was generally an adverse feature and elevated in CD120b+ GMPs. Hypomethylating agents decreased monocytic-biased cells in CMML patients. Given the enrichment of RAS pathway mutations in monocytic-biased cells, NRAS-competitive transplants and LPS-treated xenograft models recapitulated monocytic-biased CMML, suggesting that hematopoietic stress precipitates the monocytic-biased state. Deconvolution of HSPC compartments in other myeloid neoplasms and identifying therapeutic strategies to mitigate the monocytic-biased differentiation trajectory should be explored.

SIGNIFICANCE

Our findings establish that multiple differentiation states underlie CMML disease progression. These states are negatively augmented by inflammation and positively affected by hypomethylating agents. Furthermore, we identify HSC depletion and expansion of GMP-like cells with increased cytokine receptor diversity as a feature of myeloblast expansion in inflammatory chronic myeloid neoplasms. This article is highlighted in the In This Issue feature, p. 476.

Single cell clonal analysis identifies an AID-dependent pathway of plasma cell differentiation

Germinal centers (GC) are microstructures where B cells that have been activated by antigen can improve the affinity of their B cell receptors and differentiate into memory B cells (MBCs) or antibody-secreting plasma cells. Here, we have addressed the role of activation-induced deaminase (AID), which initiates somatic hypermutation and class switch recombination, in the terminal differentiation of GC B cells. By combining single cell transcriptome and immunoglobulin clonal analysis in a mouse model that traces AID-experienced cells, we have identified a novel subset of late-prePB cells (L-prePB), which shares the strongest clonal relationships with plasmablasts (PBs). Mice lacking AID have various alterations in the size and expression profiles of transcriptional clusters. We find that AID deficiency leads to a reduced proportion of L-prePB cells and severely impairs transitions between the L-prePB and the PB subsets. Thus, AID shapes the differentiation fate of GC B cells by enabling PB generation from a prePB state.

Circulating monocytes associated with anti-PD-1 resistance in human biliary cancer induce T cell paralysis

Suppressive myeloid cells can contribute to immunotherapy resistance, but their role in response to checkpoint inhibition (CPI) in anti-PD-1 refractory cancers, such as biliary tract cancer (BTC), remains elusive. We use multiplexed single-cell transcriptomic and epitope sequencing to profile greater than 200,000 peripheral blood mononuclear cells from advanced BTC patients (n = 9) and matched healthy donors (n = 8). Following anti-PD-1 treatment, CD14+ monocytes expressing high levels of immunosuppressive cytokines and chemotactic molecules (CD14CTX) increase in the circulation of patients with BTC tumors that are CPI resistant. CD14CTX can directly suppress CD4+ T cells and induce SOCS3 expression in CD4+ T cells, rendering them functionally unresponsive. The CD14CTX gene signature associates with worse survival in patients with BTC as well as in other anti-PD-1 refractory cancers. These results demonstrate that monocytes arising after anti-PD-1 treatment can induce T cell paralysis as a distinct mode of tumor-mediated immunosuppression leading to CPI resistance.

Detection of cell markers from single cell RNA-seq with sc2marker

Background

Single-cell RNA sequencing (scRNA-seq) allows the detection of rare cell types in complex tissues. The detection of markers for rare cell types is useful for further biological analysis of, for example, flow cytometry and imaging data sets for either physical isolation or spatial characterization of these cells. However, only a few computational approaches consider the problem of selecting specific marker genes from scRNA-seq data.

Results

Here, we propose sc2marker, which is based on the maximum margin index and a database of proteins with antibodies, to select markers for flow cytometry or imaging. We evaluated the performances of sc2marker and competing methods in ranking known markers in scRNA-seq data of immune and stromal cells. The results showed that sc2marker performed better than the competing methods in accuracy, while having a competitive running time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04817-5.

Self-supervised learning of cell type specificity from immunohistochemical images

MOTIVATION

Advances in bioimaging now permit in situ proteomic characterization of cell-cell interactions in complex tissues, with important applications across a spectrum of biological problems from development to disease. These methods depend on selection of antibodies targeting proteins that are expressed specifically in particular cell types. Candidate marker proteins are often identified from single-cell transcriptomic data, with variable rates of success, in part due to divergence between expression levels of proteins and the genes that encode them. In principle, marker identification could be improved by using existing databases of immunohistochemistry for thousands of antibodies in human tissue, such as the Human Protein Atlas. However, these data lack detailed annotations of the types of cells in each image.

RESULTS

We develop a method to predict cell type specificity of protein markers from unlabeled images. We train a convolutional neural network with a self-supervised objective to generate embeddings of the images. Using non-linear dimensionality reduction, we observe that the model clusters images according to cell types and anatomical regions for which the stained proteins are specific. We then use estimates of cell type specificity derived from an independent single-cell transcriptomics dataset to train an image classifier, without requiring any human labelling of images. Our scheme demonstrates superior classification of known proteomic markers in kidney compared to selection via single-cell transcriptomics.

AVAILABILITY AND IMPLEMENTATION

Code and trained model are available at www.github.com/murphy17/HPA-SimCLR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Minimal gene set discovery in single-cell mRNA-seq datasets with ActiveSVM

Sequencing costs currently prohibit the application of single-cell mRNA-seq to many biological and clinical analyses. Targeted single-cell mRNA-sequencing reduces sequencing costs by profiling reduced gene sets that capture biological information with a minimal number of genes. Here we introduce an active learning method that identifies minimal but highly informative gene sets that enable the identification of cell types, physiological states and genetic perturbations in single-cell data using a small number of genes. Our active feature selection procedure generates minimal gene sets from single-cell data by employing an active support vector machine (ActiveSVM) classifier. We demonstrate that ActiveSVM feature selection identifies gene sets that enable ~90% cell-type classification accuracy across, for example, cell atlas and disease-characterization datasets. The discovery of small but highly informative gene sets should enable reductions in the number of measurements necessary for application of single-cell mRNA-seq to clinical tests, therapeutic discovery and genetic screens.

Gibbin mesodermal regulation patterns epithelial development

Proper ectodermal patterning during human development requires previously identified transcription factors such as GATA3 and p63, as well as positional signalling from regional mesoderm1-6. However, the mechanism by which ectoderm and mesoderm factors act to stably pattern gene expression and lineage commitment remains unclear. Here we identify the protein Gibbin, encoded by the Xia-Gibbs AT-hook DNA-binding-motif-containing 1 (AHDC1) disease gene7-9, as a key regulator of early epithelial morphogenesis. We find that enhancer- or promoter-bound Gibbin interacts with dozens of sequence-specific zinc-finger transcription factors and methyl-CpG-binding proteins to regulate the expression of mesoderm genes. The loss of Gibbin causes an increase in DNA methylation at GATA3-dependent mesodermal genes, resulting in a loss of signalling between developing dermal and epidermal cell types. Notably, Gibbin-mutant human embryonic stem-cell-derived skin organoids lack dermal maturation, resulting in p63-expressing basal cells that possess defective keratinocyte stratification. In vivo chimeric CRISPR mouse mutants reveal a spectrum of Gibbin-dependent developmental patterning defects affecting craniofacial structure, abdominal wall closure and epidermal stratification that mirror patient phenotypes. Our results indicate that the patterning phenotypes seen in Xia-Gibbs and related syndromes derive from abnormal mesoderm maturation as a result of gene-specific DNA methylation decisions.

COT: an efficient and accurate method for detecting marker genes among many subtypes

Motivation

Ideally, a molecularly distinct subtype would be composed of molecular features that are expressed uniquely in the subtype of interest but in no others-so-called marker genes (MGs). MG plays a critical role in the characterization, classification or deconvolution of tissue or cell subtypes. We and others have recognized that the test statistics used by most methods do not exactly satisfy the MG definition and often identify inaccurate MG.

Results

We report an efficient and accurate data-driven method, formulated as a Cosine-based One-sample Test (COT) in scatter space, to detect MG among many subtypes using subtype expression profiles. Fundamentally different from existing approaches, the test statistic in COT precisely matches the mathematical definition of an ideal MG. We demonstrate the performance and utility of COT on both simulated and real gene expression and proteomics data. The open source Python/R tool will allow biologists to efficiently detect MG and perform a more comprehensive and unbiased molecular characterization of tissue or cell subtypes in many biomedical contexts. Nevertheless, COT complements not replaces existing methods.

Availability and implementation

The Python COT software with a detailed user's manual and a vignette are freely available at https://github.com/MintaYLu/COT.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

The landscape of receptor-mediated precision cancer combination therapy via a single-cell perspective

Mining a large cohort of single-cell transcriptomics data, here we employ combinatorial optimization techniques to chart the landscape of optimal combination therapies in cancer. We assume that each individual therapy can target any one of 1269 genes encoding cell surface receptors, which may be targets of CAR-T, conjugated antibodies or coated nanoparticle therapies. We find that in most cancer types, personalized combinations composed of at most four targets are then sufficient for killing at least 80% of tumor cells while sparing at least 90% of nontumor cells in the tumor microenvironment. However, as more stringent and selective killing is required, the number of targets needed rises rapidly. Emerging individual targets include PTPRZ1 for brain and head and neck cancers and EGFR in multiple tumor types. In sum, this study provides a computational estimate of the identity and number of targets needed in combination to target cancers selectively and precisely.

TCR-sequencing in cancer and autoimmunity: barcodes and beyond

The T cell receptor (TCR) endows T cells with antigen specificity and is central to nearly all aspects of T cell function. Each naïve T cell has a unique TCR sequence that is stably maintained during cell division. In this way, the TCR serves as a molecular barcode that tracks processes such as migration, differentiation, and proliferation of T cells. Recent technological advances have enabled sequencing of the TCR from single cells alongside deep molecular phenotypes on an unprecedented scale. In this review, we discuss strengths and limitations of TCR sequences as molecular barcodes and their application to study immune responses following Programmed Death-1 (PD-1) blockade in cancer. Additionally, we consider applications of TCR data beyond use as a barcode.

A Sight of the Diagnostic Value of Aberrant Cell-Free DNA Methylation in Lung Cancer

Background

Lung cancer is one of the most commonly diagnosed cancer worldwide. As one of the liquid biopsy analytes, alternations in cell-free DNA (cfDNA) methylation could function as promising biomarkers for lung cancer detection.

Methods

In this study, differential methylation analysis was performed to identify candidate markers, and lasso regression with 10-fold cross-validation (CV) was used to establish the diagnostic marker panel. The performance of the binary classifier was evaluated using the receiver operating characteristic (ROC) curve and the precision-recall (PR) curve.

Results

We identified 4072 differentially methylated regions (DMRs) based on cfDNA methylation data, and then a 10-DMR marker panel was established. The panel achieved an area under the ROC curve (AUROC) of 0.922 and an area under the PR curve (AUPR) of 0.899 in a cfDNA cohort containing 29 lung cancer and 74 normal samples, showing outstanding performance. Besides, the cfDNA-derived markers also performed well in primary tissue datasets, which were more robust than the tissue-derived markers.

Conclusion

Our study suggested that the 10-DMR marker panel attained high accuracy and robustness and may function as a novel and promising target for lung cancer detection.