CytoGLMM: conditional differential analysis for flow and mass cytometry experiments.

BACKGROUND: Flow and mass cytometry are important modern immunology tools for measuring expression levels of multiple proteins on single cells. The goal is to better understand the mechanisms of responses on a single cell basis by studying differential expression of proteins. Most current data analysis tools compare expressions across many computationally discovered cell types. Our goal is to focus on just one cell type. Our narrower field of application allows us to define a more specific statistical model with easier to control statistical guarantees. RESULTS: Differential analysis of marker expressions can be difficult due to marker correlations and inter-subject heterogeneity, particularly for studies of human immunology. We address these challenges with two multiple regression strategies: a bootstrapped generalized linear model and a generalized linear mixed model. On simulated datasets, we compare the robustness towards marker correlations and heterogeneity of both strategies. For paired experiments, we find that both strategies maintain the target false discovery rate under medium correlations and that mixed models are statistically more powerful under the correct model specification. For unpaired experiments, our results indicate that much larger patient sample sizes are required to detect differences. We illustrate the CytoGLMM R package and workflow for both strategies on a pregnancy dataset. CONCLUSION: Our approach to finding differential proteins in flow and mass cytometry data reduces biases arising from marker correlations and safeguards against false discoveries induced by patient heterogeneity.

Profiling of the Human Natural Killer Cell Receptor-Ligand Repertoire.

Natural killer (NK) cells are among the first responders to viral infections. The ability of NK cells to rapidly recognize and kill virally infected cells is regulated by their expression of germline-encoded inhibitory and activating receptors. The engagement of these receptors by their cognate ligands on target cells determines whether the intercellular interaction will result in NK cell killing. This protocol details the design and optimization of two complementary mass cytometry (CyTOF) panels. One panel was designed to phenotype NK cells based on receptor expression. The other panel was designed to interrogate expression of known ligands for NK cell receptors on several immune cell subsets. Together, these two panels allow for the profiling of the human NK cell receptor-ligand repertoire. Furthermore, this protocol also details the process by which we stain samples for CyTOF. This process has been optimized for improved reproducibility and standardization. An advantage of CyTOF is its ability to measure over 40 markers in each panel, with minimal signal overlap, allowing researchers to capture the breadth of the NK cell receptor-ligand repertoire. Palladium barcoding also reduces inter-sample variation, as well as consumption of reagents, making it easier to stain samples with each panel in parallel. Limitations of this protocol include the relatively low throughput of CyTOF and the inability to recover cells after analysis. These panels were designed for the analysis of clinical samples from patients suffering from acute and chronic viral infections, including dengue virus, human immunodeficiency virus (HIV), and influenza. However, they can be utilized in any setting to investigate the human NK cell receptor-ligand repertoire. Importantly, these methods can be applied broadly to the design and execution of future CyTOF panels.

Mass Cytometry Analysis of the NK Cell Receptor-Ligand Repertoire Reveals Unique Differences between Dengue-Infected Children and Adults.

Dengue virus (DENV) is a significant cause of morbidity in many regions of the world, with children at the greatest risk of developing severe dengue. NK cells, characterized by their ability to rapidly recognize and kill virally infected cells, are activated during acute DENV infection. However, their role in viral clearance versus pathogenesis has not been fully elucidated. Our goal was to profile the NK cell receptor-ligand repertoire to provide further insight into the function of NK cells during pediatric and adult DENV infection. We used mass cytometry to phenotype isolate NK cells and PBMCs from a cohort of DENV-infected children and adults. Using unsupervised clustering, we found that pediatric DENV infection leads to a decrease in total NK cell frequency with a reduction in the percentage of CD56dimCD38bright NK cells and an increase in the percentage of CD56dimperforinbright NK cells. No such changes were observed in adults. Next, we identified markers predictive of DENV infection using a differential state test. In adults, NK cell expression of activation markers, including CD69, perforin, and Fas-L, and myeloid cell expression of activating NK cell ligands, namely Fas, were predictive of infection. In contrast, increased NK cell expression of the maturation marker CD57 and myeloid cell expression of inhibitory ligands, such as HLA class I molecules, were predictive of pediatric DENV infection. These findings suggest that acute pediatric DENV infection may result in diminished NK cell activation, which could contribute to enhanced pathogenesis and disease severity.

Cytokine profile in plasma of severe COVID-19 does not differ from ARDS and sepsis.

BACKGROUNDElevated levels of inflammatory cytokines have been associated with poor outcomes among COVID-19 patients. It is unknown, however, how these levels compare with those observed in critically ill patients with acute respiratory distress syndrome (ARDS) or sepsis due to other causes.METHODSWe used a Luminex assay to determine expression of 76 cytokines from plasma of hospitalized COVID-19 patients and banked plasma samples from ARDS and sepsis patients. Our analysis focused on detecting statistical differences in levels of 6 cytokines associated with cytokine storm (IL-1ß, IL-1RA, IL-6, IL-8, IL-18, and TNF-α) between patients with moderate COVID-19, severe COVID-19, and ARDS or sepsis.RESULTSFifteen hospitalized COVID-19 patients, 9 of whom were critically ill, were compared with critically ill patients with ARDS (n = 12) or sepsis (n = 16). There were no statistically significant differences in baseline levels of IL-1ß, IL-1RA, IL-6, IL-8, IL-18, and TNF-α between patients with COVID-19 and critically ill controls with ARDS or sepsis.CONCLUSIONLevels of inflammatory cytokines were not higher in severe COVID-19 patients than in moderate COVID-19 or critically ill patients with ARDS or sepsis in this small cohort. Broad use of immunosuppressive therapies in ARDS has failed in numerous Phase 3 studies; use of these therapies in unselected patients with COVID-19 may be unwarranted.FUNDINGFunding was received from NHLBI K23 HL125663 (AJR); The Bill and Melinda Gates Foundation OPP1113682 (AJR and CAB); Burroughs Wellcome Fund Investigators in the Pathogenesis of Infectious Diseases #1016687 NIH/NIAID U19AI057229-16; Stanford Maternal Child Health Research Institute; and Chan Zuckerberg Biohub (CAB).

A single-cell atlas of the peripheral immune response to severe COVID-19.

There is an urgent need to better understand the pathophysiology of Coronavirus disease 2019 (COVID-19), the global pandemic caused by SARS-CoV-2. Here, we apply single-cell RNA sequencing (scRNA-seq) to peripheral blood mononuclear cells (PBMCs) of 7 patients hospitalized with confirmed COVID-19 and 6 healthy controls. We identify substantial reconfiguration of peripheral immune cell phenotype in COVID-19, including a heterogeneous interferon-stimulated gene (ISG) signature, HLA class II downregulation, and a novel B cell-derived granulocyte population appearing in patients with acute respiratory failure requiring mechanical ventilation. Importantly, peripheral monocytes and lymphocytes do not express substantial amounts of pro-inflammatory cytokines, suggesting that circulating leukocytes do not significantly contribute to the potential COVID-19 cytokine storm. Collectively, we provide the most thorough cell atlas to date of the peripheral immune response to severe COVID-19.

A single-cell atlas of the peripheral immune response in patients with severe COVID-19.

There is an urgent need to better understand the pathophysiology of Coronavirus disease 2019 (COVID-19), the global pandemic caused by SARS-CoV-2, which has infected more than three million people worldwide1. Approximately 20% of patients with COVID-19 develop severe disease and 5% of patients require intensive care2. Severe disease has been associated with changes in peripheral immune activity, including increased levels of pro-inflammatory cytokines3,4 that may be produced by a subset of inflammatory monocytes5,6, lymphopenia7,8 and T cell exhaustion9,10. To elucidate pathways in peripheral immune cells that might lead to immunopathology or protective immunity in severe COVID-19, we applied single-cell RNA sequencing (scRNA-seq) to profile peripheral blood mononuclear cells (PBMCs) from seven patients hospitalized for COVID-19, four of whom had acute respiratory distress syndrome, and six healthy controls. We identify reconfiguration of peripheral immune cell phenotype in COVID-19, including a heterogeneous interferon-stimulated gene signature, HLA class II downregulation and a developing neutrophil population that appears closely related to plasmablasts appearing in patients with acute respiratory failure requiring mechanical ventilation. Importantly, we found that peripheral monocytes and lymphocytes do not express substantial amounts of pro-inflammatory cytokines. Collectively, we provide a cell atlas of the peripheral immune response to severe COVID-19.

Characterization of the Impact of Daclizumab Beta on Circulating Natural Killer Cells by Mass Cytometry.

Daclizumab beta is a humanized monoclonal antibody that binds to CD25 and selectively inhibits high-affinity IL-2 receptor signaling. As a former treatment for relapsing forms of multiple sclerosis (RMS), daclizumab beta induces robust expansion of the CD56bright subpopulation of NK cells that is correlated with the drug's therapeutic effects. As NK cells represent a heterogeneous population of lymphocytes with a range of phenotypes and functions, the goal of this study was to better understand how daclizumab beta altered the NK cell repertoire to provide further insight into the possible mechanism(s) of action in RMS. We used mass cytometry to evaluate expression patterns of NK cell markers and provide a comprehensive assessment of the NK cell repertoire in individuals with RMS treated with daclizumab beta or placebo over the course of 1 year. Treatment with daclizumab beta significantly altered the NK cell repertoire compared to placebo treatment. As previously reported, daclizumab beta significantly increased expression of CD56 on total NK cells. Within the CD56bright NK cells, treatment was associated with multiple phenotypic changes, including increased expression of NKG2A and NKp44, and diminished expression of CD244, CD57, and NKp46. These alterations occurred broadly across the CD56bright population, and were not associated with a specific subset of CD56bright NK cells. While the changes were less dramatic, CD56dim NK cells responded distinctly to daclizumab beta treatment, with higher expression of CD2 and NKG2A, and lower expression of FAS-L, HLA-DR, NTB-A, NKp30, and Perforin. Together, these data indicate that the expanded CD56bright NK cells share features of both immature and mature NK cells. These findings show that daclizumab beta treatment is associated with unique changes in NK cells that may enhance their ability to kill autoreactive T cells or to exert immunomodulatory functions.

Identification of Functionally Relevant microRNAs in the Regulation of Allergic Inflammation.

Transgenic methods to manipulate CD4 T lymphocytes in vivo via forced expression of TCR transgenes and targeted "knockout" of individual genes by Cre-lox technology are fundamental to modern immunology. However, efforts to scale up functional analysis by modifying expression of larger numbers of genes in T cells ex vivo have proven surprisingly difficult. Early RNA interference experiments achieved successful small RNA transfection by using very high concentrations of short-interfering RNA (siRNA) [1], but primary T cells are generally resistant to standard electroporation, cationic liposome-, and calcium phosphate-mediated transfection methods. Moreover, although viral vectors can successfully introduce DNA fragments of varying length, expression of these constructs in primary T cells is low efficiency and the subcloning process laborious. In this context, the relatively recent discovery of dozens of highly expressed microRNAs (miRNAs) in the immune system provides both an opportunity and a new challenge [2, 3]. How can we query the miRNAome of a cell to assign particular roles to individual miRNAs? Here, we describe an optimized technique for efficient and reproducible transfection of primary mouse CD4 T cells in vitro with synthetic miRNA mimics.

Alternative splicing of interleukin-33 and type 2 inflammation in asthma.

Type 2 inflammation occurs in a large subgroup of asthmatics, and novel cytokine-directed therapies are being developed to treat this population. In mouse models, interleukin-33 (IL-33) activates lung resident innate lymphoid type 2 cells (ILC2s) to initiate airway type 2 inflammation. In human asthma, which is chronic and difficult to model, the role of IL-33 and the target cells responsible for persistent type 2 inflammation remain undefined. Full-length IL-33 is a nuclear protein and may function as an "alarmin" during cell death, a process that is uncommon in chronic stable asthma. We demonstrate a previously unidentified mechanism of IL-33 activity that involves alternative transcript splicing, which may operate in stable asthma. In human airway epithelial cells, alternative splicing of the IL-33 transcript is consistently present, and the deletion of exons 3 and 4 (Δ exon 3,4) confers cytoplasmic localization and facilitates extracellular secretion, while retaining signaling capacity. In nonexacerbating asthmatics, the expression of Δ exon 3,4 is strongly associated with airway type 2 inflammation, whereas full-length IL-33 is not. To further define the extracellular role of IL-33 in stable asthma, we sought to determine the cellular targets of its activity. Comprehensive flow cytometry and RNA sequencing of sputum cells suggest basophils and mast cells, not ILC2s, are the cellular sources of type 2 cytokines in chronic asthma. We conclude that IL-33 isoforms activate basophils and mast cells to drive type 2 inflammation in chronic stable asthma, and novel IL-33 inhibitors will need to block all biologically active isoforms.

IFN-γ-Producing T-Helper 17.1 Cells Are Increased in Sarcoidosis and Are More Prevalent than T-Helper Type 1 Cells.

RATIONALE: Pulmonary sarcoidosis is classically defined by T-helper (Th) cell type 1 inflammation (e.g., IFN-γ production by CD4(+) effector T cells). Recently, IL-17A-secreting cells have been found in lung lavage, invoking Th17 immunity in sarcoidosis. Studies also identified IL-17A-secreting cells that expressed IFN-γ, but their abundance as a percentage of total CD4(+) cells was either low or undetermined. OBJECTIVES: Based on evidence that Th17 cells can be polarized to Th17.1 cells to produce only IFN-γ, our goal was to determine whether Th17.1 cells are a prominent source of IFN-γ in sarcoidosis. METHODS: We developed a single-cell approach to define and isolate major Th-cell subsets using combinations of chemokine receptors and fluorescence-activated cell sorting. We subsequently confirmed the accuracy of subset enrichment by measuring cytokine production. MEASUREMENTS AND MAIN RESULTS: Discrimination between Th17 and Th17.1 cells revealed very high percentages of Th17.1 cells in lung lavage in sarcoidosis compared with controls in two separate cohorts. No differences in Th17 or Th1 lavage cells were found compared with controls. Lung lavage Th17.1-cell percentages were also higher than Th1-cell percentages, and approximately 60% of Th17.1-enriched cells produced only IFN-γ. CONCLUSIONS: Combined use of surface markers and functional assays to study CD4(+) T cells in sarcoidosis revealed a marked expansion of Th17.1 cells that only produce IFN-γ. These results suggest that Th17.1 cells could be misclassified as Th1 cells and may be the predominant producer of IFN-γ in pulmonary sarcoidosis, challenging the Th1 paradigm of pathogenesis.

MicroRNA regulation of lymphocyte tolerance and autoimmunity.

Understanding the cell-intrinsic cues that permit self-reactivity in lymphocytes, and therefore autoimmunity, requires an understanding of the transcriptional and posttranscriptional regulation of gene expression in these cells. In this Review, we address seminal and recent research on microRNA (miRNA) regulation of central and peripheral tolerance. Human and mouse studies demonstrate that the PI3K pathway is a critical point of miRNA regulation of immune cell development and function that affects the development of autoimmunity. We also discuss how miRNA expression profiling in human autoimmune diseases has inspired mechanistic studies of miRNA function in the pathogenesis of multiple sclerosis, rheumatoid arthritis, systemic lupus erythematosus, type 1 diabetes, and asthma.

A microRNA upregulated in asthma airway T cells promotes TH2 cytokine production.

MicroRNAs (miRNAs) exert powerful effects on immunological function by tuning networks of target genes that orchestrate cell activity. We sought to identify miRNAs and miRNA-regulated pathways that control the type 2 helper T cell (TH2 cell) responses that drive pathogenic inflammation in asthma. Profiling miRNA expression in human airway-infiltrating T cells revealed elevated expression of the miRNA miR-19a in asthma. Modulating miR-19 activity altered TH2 cytokine production in both human and mouse T cells, and TH2 cell responses were markedly impaired in cells lacking the entire miR-17∼92 cluster. miR-19 promoted TH2 cytokine production and amplified inflammatory signaling by direct targeting of the inositol phosphatase PTEN, the signaling inhibitor SOCS1 and the deubiquitinase A20. Thus, upregulation of miR-19a in asthma may be an indicator and a cause of increased TH2 cytokine production in the airways.

Interleukin-4 production by follicular helper T cells requires the conserved Il4 enhancer hypersensitivity site V.

Follicular helper T cells (Tfh cells) are the major producers of interleukin-4 (IL-4) in secondary lymphoid organs where humoral immune responses develop. Il4 regulation in Tfh cells appears distinct from the classical T helper 2 (Th2) cell pathway, but the underlying molecular mechanisms remain largely unknown. We found that hypersensitivity site V (HS V; also known as CNS2), a 3' enhancer in the Il4 locus, is essential for IL-4 production by Tfh cells. Mice lacking HS V display marked defects in type 2 humoral immune responses, as evidenced by abrogated IgE and sharply reduced IgG1 production in vivo. In contrast, effector Th2 cells that are involved in tissue responses were far less dependent on HS V. HS V facilitated removal of repressive chromatin marks during Th2 and Tfh cell differentiation and increased accessibility of the Il4 promoter. Thus, Tfh and Th2 cells utilize distinct but overlapping molecular mechanisms to regulate Il4, a finding with important implications for understanding the molecular basis of allergic diseases.

An integrated nano-scale approach to profile miRNAs in limited clinical samples.

Profiling miRNA expression in cells that directly contribute to human disease pathogenesis is likely to aid the discovery of novel drug targets and biomarkers. However, tissue heterogeneity and the limited amount of human diseased tissue available for research purposes present fundamental difficulties that often constrain the scope and potential of such studies. We established a flow cytometry-based method for isolating pure populations of pathogenic T cells from bronchial biopsy samples of asthma patients, and optimized a high-throughput nano-scale qRT-PCR method capable of accurately measuring 96 miRNAs in as little as 100 cells. Comparison of circulating and airway T cells from healthy and asthmatic subjects revealed asthma-associated and tissue-specific miRNA expression patterns. These results establish the feasibility and utility of investigating miRNA expression in small populations of cells involved in asthma pathogenesis, and set a precedent for application of our nano-scale approach in other human diseases. The microarray data from this study (Figure 7) has been submitted to the NCBI Gene Expression Omnibus (GEO; http://ncbi.nlm.nih.gov/geo) under accession no. GSE31030.

Epigenetic instability due to defective replication of structured DNA.

The accurate propagation of histone marks during chromosomal replication is proposed to rely on the tight coupling of replication with the recycling of parental histones to the daughter strands. Here, we show in the avian cell line DT40 that REV1, a key regulator of DNA translesion synthesis at the replication fork, is required for the maintenance of repressive chromatin marks and gene silencing in the vicinity of DNA capable of forming G-quadruplex (G4) structures. We demonstrate a previously unappreciated requirement for REV1 in replication of G4 forming sequences and show that transplanting a G4 forming sequence into a silent locus leads to its derepression in REV1-deficient cells. Together, our observations support a model in which failure to maintain processive DNA replication at G4 DNA in REV1-deficient cells leads to uncoupling of DNA synthesis from histone recycling, resulting in localized loss of repressive chromatin through biased incorporation of newly synthesized histones.

In vitro and in vivo high-throughput assays for the testing of anti-Trypanosoma cruzi compounds.

BACKGROUND: The two available drugs for treatment of T. cruzi infection, nifurtimox and benznidazole (BZ), have potential toxic side effects and variable efficacy, contributing to their low rate of use. With scant economic resources available for antiparasitic drug discovery and development, inexpensive, high-throughput and in vivo assays to screen potential new drugs and existing compound libraries are essential. METHODS: In this work, we describe the development and validation of improved methods to test anti-T. cruzi compounds in vitro and in vivo using parasite lines expressing the firefly luciferase (luc) or the tandem tomato fluorescent protein (tdTomato). For in vitro assays, the change in fluorescence intensity of tdTomato-expressing lines was measured as an indicator of parasite replication daily for 4 days and this method was used to identify compounds with IC(50) lower than that of BZ. FINDINGS: This method was highly reproducible and had the added advantage of requiring relatively low numbers of parasites and no additional indicator reagents, enzymatic post-processes or laborious visual counting. In vivo, mice were infected in the footpads with fluorescent or bioluminescent parasites and the signal intensity was measured as a surrogate of parasite load at the site of infection before and after initiation of drug treatment. Importantly, the efficacy of various drugs as determined in this short-term (<2 weeks) assay mirrored that of a 40 day treatment course. CONCLUSION: These methods should make feasible broader and higher-throughput screening programs needed to identify potential new drugs for the treatment of T. cruzi infection and for their rapid validation in vivo.

Insufficient TLR activation contributes to the slow development of CD8+ T cell responses in Trypanosoma cruzi infection.

During experimental infection with Trypanosoma cruzi, mice develop a strong CD8(+) T cell response focused mainly on a few immunodominant peptides encoded in trans-sialidase family genes. Despite the potency of this response, the initial emergence and peak of parasite-specific CD8(+) T cells has been noted to be relatively slow. In this study, we further document this delayed onset of T cell responses to T. cruzi as measured by the increase in frequency of parasite-specific T cells, the effector function of these cells, T cell proliferation in general, and the recruitment of cells into the draining lymph nodes. This delay does not appear to be the result of general immunosuppressive effects of the infection, a limitation in parasite numbers, or parasite trafficking to lymph nodes or to the specific epitope. Increasing the initial infecting dose or the density of parasite epitopes on APCs can modestly speed the generation of anti-T. cruzi T cell responses. Given these characteristics of the response, we propose that T. cruzi is a stealth invader, largely avoiding recognition by components of the innate immune system until the infection is well established. This conclusion is supported by the ability to accelerate the induction of T cell responses to T. cruzi by administration of ligands for TLR2 and TLR9 at the time of infection. These studies highlight a previously unappreciated mechanism of immune evasion, the surreptitious establishment of infection, by the protozoan T. cruzi.

Timing matters: error-prone gap filling and translesion synthesis in immunoglobulin gene hypermutation.

By temporarily deferring the repair of DNA lesions encountered during replication, the bypass of DNA damage is critical to the ability of cells to withstand genomic insults. Damage bypass can be achieved either by recombinational mechanisms that are generally accurate or by a process called translesion synthesis. Translesion synthesis involves replacing the stalled replicative polymerase with one of a number of specialized DNA polymerases whose active sites are able to tolerate a distorted or damaged DNA template. While this property allows the translesion polymerases to synthesize across damaged bases, it does so with the trade-off of an increased mutation rate. The deployment of these enzymes must therefore be carefully regulated. In addition to their important role in general DNA damage tolerance and mutagenesis, the translesion polymerases play a crucial role in converting the products of activation induced deaminase-catalysed cytidine deamination to mutations during immunoglobulin gene somatic hypermutation. In this paper, we specifically consider the control of translesion synthesis in the context of the timing of lesion bypass relative to replication fork progression and arrest at sites of DNA damage. We then examine how recent observations concerning the control of translesion synthesis might help refine our view of the mechanisms of immunoglobulin gene somatic hypermutation.