cytoKernel: Robust kernel embeddings for assessing differential expression of single cell data.

High-throughput sequencing of single-cell data can be used to rigorously evlauate cell specification and enable intricate variations between groups or conditions. Many popular existing methods for differential expression target differences in aggregate measurements (mean, median, sum) and limit their approaches to detect only global differential changes. We present a robust method for differential expression of single-cell data using a kernel-based score test, cytoKernel. cytoKernel is specifically designed to assess the differential expression of single cell RNA sequencing and high-dimensional flow or mass cytometry data using the full probability distribution pattern. cytoKernel is based on kernel embeddings which employs the probability distributions of the single cell data, by calculating the pairwise divergence/distance between distributions of subjects. It can detect both patterns involving aggregate changes, as well as more elusive variations that are often overlooked due to the multimodal characteristics of single cell data. We performed extensive benchmarks across both simulated and real data sets from mass cytometry data and single-cell RNA sequencing. The cytoKernel procedure effectively controls the False Discovery Rate (FDR) and shows favourable performance compared to existing methods. The method is able to identify more differential patterns than existing approaches. We apply cytoKernel to assess gene expression and protein marker expression differences from cell subpopulations in various publicly available single-cell RNAseq and mass cytometry data sets. The methods described in this paper are implemented in the open-source R package cytoKernel, which is freely available from Bioconductor at http://bioconductor.org/packages/cytoKernel.

A partial human LCK defect causes a T cell immunodeficiency with intestinal inflammation.

Lymphocyte-specific protein tyrosine kinase (LCK) is essential for T cell antigen receptor (TCR)-mediated signal transduction. Here, we report two siblings homozygous for a novel LCK variant (c.1318C>T; P440S) characterized by T cell lymphopenia with skewed memory phenotype, infant-onset recurrent infections, failure to thrive, and protracted diarrhea. The patients' T cells show residual TCR signal transduction and proliferation following anti-CD3/CD28 and phytohemagglutinin (PHA) stimulation. We demonstrate in mouse models that complete (Lck-/-) versus partial (LckP440S/P440S) loss-of-function LCK causes disease with differing phenotypes. While both Lck-/- and LckP440S/P440S mice exhibit arrested thymic T cell development and profound T cell lymphopenia, only LckP440S/P440S mice show residual T cell proliferation, cytokine production, and intestinal inflammation. Furthermore, the intestinal disease in the LckP440S/P440S mice is prevented by CD4+ T cell depletion or regulatory T cell transfer. These findings demonstrate that P440S LCK spares sufficient T cell function to allow the maturation of some conventional T cells but not regulatory T cells-leading to intestinal inflammation.

Dysregulated Lymphocyte Antigen Receptor Signaling in Common Variable Immunodeficiency with Granulomatous Lymphocytic Interstitial Lung Disease.

PURPOSE: A subset of common variable immunodeficiency (CVID) patients either presents with or develops autoimmune and lymphoproliferative complications, such as granulomatous lymphocytic interstitial lung disease (GLILD), a major cause of morbidity and mortality in CVID. While a myriad of phenotypic lymphocyte derangements has been associated with and described in GLILD, defects in T and B cell antigen receptor (TCR/BCR) signaling in CVID and CVID with GLILD (CVID/GLILD) remain undefined, hindering discovery of biomarkers for disease monitoring, prognostic prediction, and personalized medicine approaches. METHODS: To identify perturbations of immune cell subsets and TCR/BCR signal transduction, we applied mass cytometry analysis to peripheral blood mononuclear cells (PBMCs) from healthy control participants (HC), CVID, and CVID/GLILD patients. RESULTS: Patients with CVID, regardless of GLILD status, had increased frequency of HLADR+CD4+ T cells, CD57+CD8+ T cells, and CD21lo B cells when compared to healthy controls. Within these cellular populations in CVID/GLILD patients only, engagement of T or B cell antigen receptors resulted in discordant downstream signaling responses compared to CVID. In CVID/GLILD patients, CD21lo B cells showed perturbed BCR-mediated phospholipase C gamma and extracellular signal-regulated kinase activation, while HLADR+CD4+ T cells and CD57+CD8+ T cells displayed disrupted TCR-mediated activation of kinases most proximal to the receptor. CONCLUSION: Both CVID and CVID/GLILD patients demonstrate an activated T and B cell phenotype compared to HC. However, only CVID/GLILD patients exhibit altered TCR/BCR signaling in the activated lymphocyte subsets. These findings contribute to our understanding of the mechanisms of immune dysregulation in CVID with GLILD.