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.

SARS-CoV-2 Vaccine-Elicited Immunity after B Cell Depletion in Multiple Sclerosis.

The impact of B cell deficiency on the humoral and cellular responses to SARS-CoV2 mRNA vaccination remains a challenging and significant clinical management question. We evaluated vaccine-elicited serological and cellular responses in 1) healthy individuals who were pre-exposed to SARS-CoV-2 (n = 21), 2) healthy individuals who received a homologous booster (mRNA, n = 19; or Novavax, n = 19), and 3) persons with multiple sclerosis on B cell depletion therapy (MS-αCD20) receiving mRNA homologous boosting (n = 36). Pre-exposure increased humoral and CD4 T cellular responses in immunocompetent individuals. Novavax homologous boosting induced a significantly more robust serological response than mRNA boosting. MS-α CD20 had an intact IgA mucosal response and an enhanced CD8 T cell response to mRNA boosting compared with immunocompetent individuals. This enhanced cellular response was characterized by the expansion of only effector, not memory, T cells. The enhancement of CD8 T cells in the setting of B cell depletion suggests a regulatory mechanism between B and CD8 T cell vaccine responses.

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.

Innate and Adaptive Cell-Mediated Immune Responses to a COVID-19 mRNA Vaccine in Young Children.

Background: There is little information on cell-mediated immunity (CMI) to COVID-19 mRNA vaccines in children. We studied adaptive and innate CMI in vaccinated children aged 6 to 60 months. Methods: Blood obtained from participants in a randomized placebo-controlled trial of an mRNA vaccine before and 1 month after the first dose was used for antibody measurements and CMI (flow cytometry). Results: We enrolled 29 children with a mean age of 28.5 months (SD, 15.7). Antibody studies revealed that 10 participants were infected with SARS-CoV-2 prevaccination. Ex vivo stimulation of peripheral blood mononuclear cells with SARS-CoV-2 spike peptides showed significant increases pre- to postimmunization of activated conventional CD4+ and γδ T cells, natural killer cells, monocytes, and conventional dendritic cells but not mucosa-associated innate T cells. Conventional T-cell, monocyte, and conventional dendritic cell responses in children were higher immediately after vaccination than after SARS-CoV-2 infection. The fold increase in CMI pre- to postvaccination did not differ between children previously infected with SARS-CoV-2 and those uninfected. Conclusions: Children aged 6 to 60 months who were vaccinated with a COVID-19 mRNA vaccine developed robust CMI responses, including adaptive and innate immunity.

Expansion of extrafollicular B and T cell subsets in childhood-onset systemic lupus erythematosus.

Introduction: Most childhood-onset SLE patients (cSLE) develop lupus nephritis (cLN), but only a small proportion achieve complete response to current therapies. The prognosis of children with LN and end-stage renal disease is particularly dire. Mortality rates within the first five years of renal replacement therapy may reach 22%. Thus, there is urgent need to decipher and target immune mechanisms that drive cLN. Despite the clear role of autoantibody production in SLE, targeted B cell therapies such as rituximab (anti-CD20) and belimumab (anti-BAFF) have shown only modest efficacy in cLN. While many studies have linked dysregulation of germinal center formation to SLE pathogenesis, other work supports a role for extrafollicular B cell activation in generation of pathogenic antibody secreting cells. However, whether extrafollicular B cell subsets and their T cell collaborators play a role in specific organ involvement in cLN and/or track with disease activity remains unknown. Methods: We analyzed high-dimensional mass cytometry and gene expression data from 24 treatment naïve cSLE patients at the time of diagnosis and longitudinally, applying novel computational tools to identify abnormalities associated with clinical manifestations (cLN) and disease activity (SLEDAI). Results: cSLE patients have an extrafollicular B cell expansion signature, with increased frequency of i) DN2, ii) Bnd2, iii) plasmablasts, and iv) peripheral T helper cells. Most importantly, we discovered that this extrafollicular signature correlates with disease activity in cLN, supporting extrafollicular T/B interactions as a mechanism underlying pediatric renal pathogenesis. Discussion: This study integrates established and emerging themes of extrafollicular B cell involvement in SLE by providing evidence for extrafollicular B and peripheral T helper cell expansion, along with elevated type 1 IFN activation, in a homogeneous cohort of treatment-naïve cSLE patients, a point at which they should display the most extreme state of their immune dysregulation.

SARS-CoV-2 mRNA vaccination induces an antigen-specific T cell response correlating with plasma interferon-gamma in B cell depleted patients.

Emerging evidence is encouraging and suggests that a substantial proportion of patients without antibody responses (due to anti-CD20 therapy or other etiologies) to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) vaccines develop T cell responses. However, antigen-specific T cellular responses are notoriously difficult to assess clinically, given the lack of such assays under satisfactory CAP/CLIA regulation, and the laborious nature of the flow cytometric assessment. To evaluate the ability to apply a clinically feasible assay to measure T cellular responses to SARS-CoV-2 mRNA vaccination, we compared flow cytometric and enzyme-linked immunosorbent assay (ELISA) based assays in 24 participants treated with anti-CD20 therapy. T cellular activation (CD69 + CD137+ surface expression, i.e., activation induced markers [AIM]) and intracellular interferon gamma (INFγ) production via flow cytometry was compared to plasma Interferon Gamma Release Assay (IGRA) via ELISA. Plasma INFγ production measured by IGRA correlated with the percent of INFγ-producing AIM positive T cells, supporting the use of IGRA assay as a robust assessment of T cellular response to the SARS-CoV-2 vaccine for B-cell depleted patients that is clinically feasible, time efficient, and cost effective.

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.

Enhanced Reactivity for Aromatic Bromination via Halogen Bonding with Lactic Acid Derivatives.

We report a new method for regioselective aromatic bromination using lactic acid derivatives as halogen bond acceptors with N-bromosuccinimide (NBS). Several structural analogues of lactic acid affect the efficiency of aromatic brominations, presumably via Lewis acid/base halogen-bonding interactions. Rate comparisons of aromatic brominations demonstrate the reactivity enhancement available via catalytic additives capable of halogen bonding. Computational results demonstrate that Lewis basic additives interact with NBS to increase the electropositive character of bromine prior to electrophilic transfer. An optimized procedure using catalytic mandelic acid under aqueous conditions at room temperature was developed to promote aromatic bromination on a variety of arene substrates with complete regioselectivity.

Compositional Data Analysis using Kernels in mass cytometry data.

MOTIVATION: Cell-type abundance data arising from mass cytometry experiments are compositional in nature. Classical association tests do not apply to the compositional data due to their non-Euclidean nature. Existing methods for analysis of cell type abundance data suffer from several limitations for high-dimensional mass cytometry data, especially when the sample size is small. RESULTS: We proposed a new multivariate statistical learning methodology, Compositional Data Analysis using Kernels (CODAK), based on the kernel distance covariance (KDC) framework to test the association of the cell type compositions with important predictors (categorical or continuous) such as disease status. CODAK scales well for high-dimensional data and provides satisfactory performance for small sample sizes (n < 25). We conducted simulation studies to compare the performance of the method with existing methods of analyzing cell type abundance data from mass cytometry studies. The method is also applied to a high-dimensional dataset containing different subgroups of populations including Systemic Lupus Erythematosus (SLE) patients and healthy control subjects. AVAILABILITY AND IMPLEMENTATION: CODAK is implemented using R. The codes and the data used in this manuscript are available on the web at http://github.com/GhoshLab/CODAK/. CONTACT: prudra@okstate.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics Advances online.

Specialized interferon action in COVID-19.

The impacts of interferon (IFN) signaling on COVID-19 pathology are multiple, with both protective and harmful effects being documented. We report here a multiomics investigation of systemic IFN signaling in hospitalized COVID-19 patients, defining the multiomics biosignatures associated with varying levels of 12 different type I, II, and III IFNs. The antiviral transcriptional response in circulating immune cells is strongly associated with a specific subset of IFNs, most prominently IFNA2 and IFNG. In contrast, proteomics signatures indicative of endothelial damage and platelet activation associate with high levels of IFNB1 and IFNA6. Seroconversion and time since hospitalization associate with a significant decrease in a specific subset of IFNs. Additionally, differential IFN subtype production is linked to distinct constellations of circulating myeloid and lymphoid immune cell types. Each IFN has a unique metabolic signature, with IFNG being the most associated with activation of the kynurenine pathway. IFNs also show differential relationships with clinical markers of poor prognosis and disease severity. For example, whereas IFNG has the strongest association with C-reactive protein and other immune markers of poor prognosis, IFNB1 associates with increased neutrophil to lymphocyte ratio, a marker of late severe disease. Altogether, these results reveal specialized IFN action in COVID-19, with potential diagnostic and therapeutic implications.

Peripheral immunophenotyping of AITD subjects reveals alterations in immune cells in pediatric vs adult-onset AITD.

Autoimmune thyroid disease (AITD) is caused by aberrant activation of the immune system allowing autoreactive B and T cells to target the thyroid gland leading to disease. Although AITD is more frequently diagnosed in adults, children are also affected but rarely studied. Here, we performed phenotypic and functional characterization of peripheral blood immune cells from pediatric and adult-onset AITD patients and age-matched controls using mass cytometry. Major findings indicate that unlike adult-onset AITD patients, pediatric AITD patients exhibit a decrease in anergic B cells (BND) and DN2 B cells and an increase in immature B cells compared to age-matched controls. These results indicate alterations in peripheral blood immune cells seen in pediatric-onset AITD could lead to rapid progression of disease. Hence, this study demonstrates diversity of AITD by showing differences in immune cell phenotypes and function based on age of onset, and may inform future therapies.

Tuning three-dimensional nano-assembly in the mesoscale via bis(imino)pyridine molecular functionalization.

We investigate the effect of bis(imino)pyridine (BIP) ligands in guiding self-assembly of semiconducting CdSe/ZnS quantum dots (QDs) into three-dimensional multi-layered shells with diameters spanning the entire mesoscopic range, from 200 nm to 2 µm. The assembly process is directed by guest-host interactions between the BIP ligands and a thermotropic liquid crystal (LC), with the latter's phase transition driving the process. Characterization of the shell structures, through scanning electron microscopy and dynamic light scattering, demonstrates that the average shell diameter depends on the BIP structure, and that changing one functional group in the chemical scaffold allows systematic tuning of shell sizes across the entire range. Differential scanning calorimetry confirms a relationship between shell sizes and the thermodynamic perturbation of the BIP molecules to the LC phase transition temperature, allowing analytical modeling of shell assembly energetics. This novel mechanism to controllably tune shell sizes over the entire mesoscale via one standard protocol is a significant development for research on in situ cargo/drug delivery platforms using nano-assembled structures.

T and B cell abnormalities, pneumocystis pneumonia, and chronic lymphocytic leukemia associated with an AIOLOS defect in patients.

AIOLOS/IKZF3 is a member of the IKAROS family of transcription factors. IKAROS/IKZF1 mutations have been previously associated with different forms of primary immunodeficiency. Here we describe a novel combined immunodeficiency due to an IKZF3 mutation in a family presenting with T and B cell involvement, Pneumocystis jirovecii pneumonia, and/or chronic lymphocytic leukemia. Patients carrying the AIOLOS p.N160S heterozygous variant displayed impaired humoral responses, abnormal B cell development (high percentage of CD21low B cells and negative CD23 expression), and abrogated CD40 responses. Naive T cells were increased, T cell differentiation was abnormal, and CD40L expression was dysregulated. In vitro studies demonstrated that the mutant protein failed DNA binding and pericentromeric targeting. The mutant was fully penetrant and had a dominant-negative effect over WT AIOLOS but not WT IKAROS. The human immunophenotype was recapitulated in a murine model carrying the corresponding human mutation. As demonstrated here, AIOLOS plays a key role in T and B cell development in humans, and the particular gene variant described is strongly associated with immunodeficiency and likely malignancy.

The COVIDome Explorer researcher portal.

COVID-19 pathology involves dysregulation of diverse molecular, cellular, and physiological processes. To expedite integrated and collaborative COVID-19 research, we completed multi-omics analysis of hospitalized COVID-19 patients, including matched analysis of the whole-blood transcriptome, plasma proteomics with two complementary platforms, cytokine profiling, plasma and red blood cell metabolomics, deep immune cell phenotyping by mass cytometry, and clinical data annotation. We refer to this multidimensional dataset as the COVIDome. We then created the COVIDome Explorer, an online researcher portal where the data can be analyzed and visualized in real time. We illustrate herein the use of the COVIDome dataset through a multi-omics analysis of biosignatures associated with C-reactive protein (CRP), an established marker of poor prognosis in COVID-19, revealing associations between CRP levels and damage-associated molecular patterns, depletion of protective serpins, and mitochondrial metabolism dysregulation. We expect that the COVIDome Explorer will rapidly accelerate data sharing, hypothesis testing, and discoveries worldwide.

Versatile New Reagent for Nitrosation under Mild Conditions.

Here we report a new chemical reagent for transnitrosation under mild experimental conditions. This new reagent is stable to air and moisture across a broad range of temperatures and is effective for transnitrosation in multiple solvents. Compared with traditional nitrosation methods, our reagent shows high functional group tolerance for substrates that are susceptible to oxidation or reversible transnitrosation. Several challenging nitroso compounds are accessed here for the first time, including 15N isotopologues. X-ray data confirm that two rotational isomers of the reagent are configurationally stable at room temperature, although only one isomer is effective for transnitrosation. Computational analysis describes the energetics of rotamer interconversion, including interesting geometry-dependent hybridization effects.

The COVIDome Explorer Researcher Portal.

COVID-19 pathology involves dysregulation of diverse molecular, cellular, and physiological processes. In order to expedite integrated and collaborative COVID-19 research, we completed multi-omics analysis of hospitalized COVID-19 patients including matched analysis of the whole blood transcriptome, plasma proteomics with two complementary platforms, cytokine profiling, plasma and red blood cell metabolomics, deep immune cell phenotyping by mass cytometry, and clinical data annotation. We refer to this multidimensional dataset as the COVIDome. We then created the COVIDome Explorer, an online researcher portal where the data can be analyzed and visualized in real time. We illustrate here the use of the COVIDome dataset through a multi-omics analysis of biosignatures associated with C-reactive protein (CRP), an established marker of poor prognosis in COVID-19, revealing associations between CRP levels and damage-associated molecular patterns, depletion of protective serpins, and mitochondrial metabolism dysregulation. We expect that the COVIDome Explorer will rapidly accelerate data sharing, hypothesis testing, and discoveries worldwide.

Seroconversion stages COVID19 into distinct pathophysiological states.

COVID19 is a heterogeneous medical condition involving diverse underlying pathophysiological processes including hyperinflammation, endothelial damage, thrombotic microangiopathy, and end-organ damage. Limited knowledge about the molecular mechanisms driving these processes and lack of staging biomarkers hamper the ability to stratify patients for targeted therapeutics. We report here the results of a cross-sectional multi-omics analysis of hospitalized COVID19 patients revealing that seroconversion status associates with distinct underlying pathophysiological states. Low antibody titers associate with hyperactive T cells and NK cells, high levels of IFN alpha, gamma and lambda ligands, markers of systemic complement activation, and depletion of lymphocytes, neutrophils, and platelets. Upon seroconversion, all of these processes are attenuated, observing instead increases in B cell subsets, emergency hematopoiesis, increased D-dimer, and hypoalbuminemia. We propose that seroconversion status could potentially be used as a biosignature to stratify patients for therapeutic intervention and to inform analysis of clinical trial results in heterogenous patient populations.

Seroconversion stages COVID19 into distinct pathophysiological states.

COVID19 is a heterogeneous medical condition involving a suite of underlying pathophysiological processes including hyperinflammation, endothelial damage, thrombotic microangiopathy, and end-organ damage. Limited knowledge about the molecular mechanisms driving these processes and lack of staging biomarkers hamper the ability to stratify patients for targeted therapeutics. We report here the results of a cross-sectional multi-omics analysis of hospitalized COVID19 patients revealing that seroconversion status associates with distinct underlying pathophysiological states. Seronegative COVID19 patients harbor hyperactive T cells and NK cells, high levels of IFN alpha, gamma and lambda ligands, markers of systemic complement activation, neutropenia, lymphopenia and thrombocytopenia. In seropositive patients, all of these processes are attenuated, observing instead increases in B cell subsets, emergency hematopoiesis, increased markers of platelet activation, and hypoalbuminemia. We propose that seroconversion status could potentially be used as a biosignature to stratify patients for therapeutic intervention and to inform analysis of clinical trial results in heterogenous patient populations.

Modeling broadband cloaking using 3D nano-assembled plasmonic meta-structures.

The concept of "cloaking" an object is a very attractive one, especially in the visible (VIS) and near infra-red (NIR) regions of the electromagnetic spectrum, as that would reduce the visibility of an object to the eye. One possible route to achieving this goal is by leveraging the plasmonic property of metallic nanoparticles (NPs). We model and simulate light in the VIS and NIR scattered by a core of a homogeneous medium, covered by plasmonic cloak that is a spherical shell composed of gold nanoparticles (AuNPs). To consider realistic, scalable, and robust plasmonic cloaks that are comparable, or larger, in size to the wavelength, we introduce a multiscale simulation platform. This model uses the multiple scattering theory of Foldy and Lax to model interactions of light with AuNPs combined with the method of fundamental solutions to model interactions with the core. Numerical results of our simulations for the scattering cross-sections of core-shell composite indicate significant scattering suppression of up to 50% over a substantial portion of the desired spectral range (400 - 600 nm) for cores as large as 900 nm in diameter by a suitable combination of AuNP sizes and filling fractions of AuNPs in the shell.

Minimizing Batch Effects in Mass Cytometry Data.

Cytometry by Time-Of-Flight (CyTOF) uses antibodies conjugated to isotopically pure metals to identify and quantify a large number of cellular features with single-cell resolution. A barcoding approach allows for 20 unique samples to be pooled and processed together in one tube, reducing the intra-barcode technical variability. However, with only 20 samples per barcode, multiple barcode sets (batches) are required to address questions in robustly powered study designs. A batch adjustment procedure is required to reduce variability across batches and to facilitate direct comparison of runs performed across multiple barcodes run over weeks, months, or years. We describe a method using technical replicates that are included in each run to determine and apply an appropriate adjustment per batch without manual intervention. The use of technical replicate samples (i.e., anchors or reference samples) avoids assumptions of sample homogeneity among batches, and allows direct estimation of batch effects and appropriate adjustment parameters applicable to all samples within a batch. Quantification of cell subpopulations and mean signal intensity pre- and post-adjustment using both manual gating and unsupervised clustering demonstrate substantial mitigation of batch effects in the anchor samples used for this adjustment calculation, and in a second validation set of technical replicates.