Identification and visualization of multidimensional antigen-specific T-cell populations in polychromatic cytometry data.

An important aspect of immune monitoring for vaccine development, clinical trials, and research is the detection, measurement, and comparison of antigen-specific T-cells from subject samples under different conditions. Antigen-specific T-cells compose a very small fraction of total T-cells. Developments in cytometry technology over the past five years have enabled the measurement of single-cells in a multivariate and high-throughput manner. This growth in both dimensionality and quantity of data continues to pose a challenge for effective identification and visualization of rare cell subsets, such as antigen-specific T-cells. Dimension reduction and feature extraction play pivotal role in both identifying and visualizing cell populations of interest in large, multi-dimensional cytometry datasets. However, the automated identification and visualization of rare, high-dimensional cell subsets remains challenging. Here we demonstrate how a systematic and integrated approach combining targeted feature extraction with dimension reduction can be used to identify and visualize biological differences in rare, antigen-specific cell populations. By using OpenCyto to perform semi-automated gating and features extraction of flow cytometry data, followed by dimensionality reduction with t-SNE we are able to identify polyfunctional subpopulations of antigen-specific T-cells and visualize treatment-specific differences between them

A computational framework to emulate the human perspective in flow cytometric data analysis

Background: In recent years, intense research efforts have focused on developing methods for automated flow cytometric data analysis. However, while designing such applications, little or no attention has been paid to the human perspective that is absolutely central to the manual gating process of identifying and characterizing cell populations. In particular, the assumption of many common techniques that cell populations could be modeled reliably with pre-specified distributions may not hold true in real-life samples, which can have populations of arbitrary shapes and considerable inter-sample variation. <p/>Results: To address this, we developed a new framework flowScape for emulating certain key aspects of the human perspective in analyzing flow data, which we implemented in multiple steps. First, flowScape begins with creating a mathematically rigorous map of the high-dimensional flow data landscape based on dense and sparse regions defined by relative concentrations of events around modes. In the second step, these modal clusters are connected with a global hierarchical structure. This representation allows flowScape to perform ridgeline analysis for both traversing the landscape and isolating cell populations at different levels of resolution. Finally, we extended manual gating with a new capacity for constructing templates that can identify target populations in terms of their relative parameters, as opposed to the more commonly used absolute or physical parameters. This allows flowScape to apply such templates in batch mode for detecting the corresponding populations in a flexible, sample-specific manner. We also demonstrated different applications of our framework to flow data analysis and show its superiority over other analytical methods. <p/>Conclusions: The human perspective, built on top of intuition and experience, is a very important component of flow cytometric data analysis. By emulating some of its approaches and extending these with automation and rigor, flowScape provides a flexible and robust framework for computational cytomics

Deep generative modeling for single-cell transcriptomics.

Single-cell transcriptome measurements can reveal unexplored biological diversity, but they suffer from technical noise and bias that must be modeled to account for the resulting uncertainty in downstream analyses. Here we introduce single-cell variational inference (scVI), a ready-to-use scalable framework for the probabilistic representation and analysis of gene expression in single cells ( https://github.com/YosefLab/scVI ). scVI uses stochastic optimization and deep neural networks to aggregate information across similar cells and genes and to approximate the distributions that underlie observed expression values, while accounting for batch effects and limited sensitivity. We used scVI for a range of fundamental analysis tasks including batch correction, visualization, clustering, and differential expression, and achieved high accuracy for each task

SerpinB2 regulates stromal remodelling and local invasion in pancreatic cancer

Pancreatic cancer has a devastating prognosis, with an overall 5-year survival rate of ~8%, restricted treatment options and characteristic molecular heterogeneity. SerpinB2 expression, particularly in the stromal compartment, is associated with reduced metastasis and prolonged survival in pancreatic ductal adenocarcinoma (PDAC) and our genomic analysis revealed that SERPINB2 is frequently deleted in PDAC. We show that SerpinB2 is required by stromal cells for normal collagen remodelling in vitro, regulating fibroblast interaction and engagement with collagen in the contracting matrix. In a pancreatic cancer allograft model, co-injection of PDAC cancer cells and SerpinB2(-/-) mouse embryonic fibroblasts (MEFs) resulted in increased tumour growth, aberrant remodelling of the extracellular matrix (ECM) and increased local invasion from the primary tumour. These tumours also displayed elevated proteolytic activity of the primary biochemical target of SerpinB2-urokinase plasminogen activator (uPA). In a large cohort of patients with resected PDAC, we show that increasing uPA mRNA expression was significantly associated with poorer survival following pancreatectomy. This study establishes a novel role for SerpinB2 in the stromal compartment in PDAC invasion through regulation of stromal remodelling and highlights the SerpinB2/uPA axis for further investigation as a potential therapeutic target in pancreatic cancer

Normalizing single-cell RNA sequencing data: challenges and opportunities

Single-cell transcriptomics is becoming an important component of the molecular biologist's toolkit. A critical step when analyzing data generated using this technology is normalization. However, normalization is typically performed using methods developed for bulk RNA sequencing or even microarray data, and the suitability of these methods for single-cell transcriptomics has not been assessed. We here discuss commonly used normalization approaches and illustrate how these can produce misleading results. Finally, we present alternative approaches and provide recommendations for single-cell RNA sequencing users

COMPASS identifies T-cell subsets correlated with clinical outcomes.

Advances in flow cytometry and other single-cell technologies have enabled high-dimensional, high-throughput measurements of individual cells as well as the interrogation of cell population heterogeneity. However, in many instances, computational tools to analyze the wealth of data generated by these technologies are lacking. Here, we present a computational framework for unbiased combinatorial polyfunctionality analysis of antigen-specific T-cell subsets (COMPASS). COMPASS uses a Bayesian hierarchical framework to model all observed cell subsets and select those most likely to have antigen-specific responses. Cell-subset responses are quantified by posterior probabilities, and human subject-level responses are quantified by two summary statistics that describe the quality of an individual's polyfunctional response and can be correlated directly with clinical outcome. Using three clinical data sets of cytokine production, we demonstrate how COMPASS improves characterization of antigen-specific T cells and reveals cellular 'correlates of protection/immunity' in the RV144 HIV vaccine efficacy trial that are missed by other methods. COMPASS is available as open-source software

Gene expression profiling of tumour epithelial and stromal compartments during breast cancer progression

The progression of ductal carcinoma in situ (DCIS) to invasive ductal carcinoma (IDC) marks a critical step in the evolution of breast cancer. There is some evidence to suggest that dynamic interactions between the neoplastic cells and the tumour microenvironment play an important role. Using the whole-genome cDNA-mediated annealing, selection, extension and ligation assay (WG-DASL, Illumina), we performed gene expression profiling on 87 formalin-fixed paraffin-embedded (FFPE) samples from 17 patients consisting of matched IDC, DCIS and three types of stroma: IDC-S ( 10 mm from IDC or DCIS). Differential gene expression analysis was validated by quantitative real time-PCR, immunohistochemistry and immunofluorescence. The expression of several genes was down-regulated in stroma from cancer patients relative to normal stroma from reduction mammoplasties. In contrast, neoplastic epithelium underwent more gene expression changes during progression, including down regulation of SFRP1. In particular, we observed that molecules related to extracellular matrix (ECM) remodelling (e.g. COL11A1, COL5A2 and MMP13) were differentially expressed between DCIS and IDC. COL11A1 was overexpressed in IDC relative to DCIS and was expressed by both the epithelial and stromal compartments but was enriched in invading neoplastic epithelial cells. The contributions of both the epithelial and stromal compartments to the clinically important scenario of progression from DCIS to IDC. Gene expression profiles, we identified differential expression of genes related to ECM remodelling, and specifically the elevated expression of genes such as COL11A1, COL5A2 and MMP13 in epithelial cells of IDC. We propose that these expression changes could be involved in facilitating the transition from in situ disease to invasive cancer and may thus mark a critical point in disease development

Clinico-pathological and transcriptomic determinants of SLFN11 expression in invasive breast carcinoma

SCOPUS: ar.jinfo:eu-repo/semantics/publishe

Parity-related molecular signatures and breast cancer subtypes by estrogen receptor status

INTRODUCTION: Relationships of parity with breast cancer risk are complex. Parity is associated with decreased risk of postmenopausal hormone receptor–positive breast tumors, but may increase risk for basal-like breast cancers and early-onset tumors. Characterizing parity-related gene expression patterns in normal breast and breast tumor tissues may improve understanding of the biological mechanisms underlying this complex pattern of risk. METHODS: We developed a parity signature by analyzing microRNA microarray data from 130 reduction mammoplasty (RM) patients (54 nulliparous and 76 parous). This parity signature, together with published parity signatures, was evaluated in gene expression data from 150 paired tumors and adjacent benign breast tissues from the Polish Breast Cancer Study, both overall and by tumor estrogen receptor (ER) status. RESULTS: We identified 251 genes significantly upregulated by parity status in RM patients (parous versus nulliparous; false discovery rate = 0.008), including genes in immune, inflammation and wound response pathways. This parity signature was significantly enriched in normal and tumor tissues of parous breast cancer patients, specifically in ER-positive tumors. CONCLUSIONS: Our data corroborate epidemiologic data, suggesting that the etiology and pathogenesis of breast cancers vary by ER status, which may have implications for developing prevention strategies for these tumors