Systems Immunology: Learning the Rules of the Immune System.

Given the many cell types and molecular components of the human immune system, along with vast variations across individuals, how should we go about developing causal and predictive explanations of immunity? A central strategy in human studies is to leverage natural variation to find relationships among variables, including DNA variants, epigenetic states, immune phenotypes, clinical descriptors, and others. Here, we focus on how natural variation is used to find patterns, infer principles, and develop predictive models for two areas: (a) immune cell activation-how single-cell profiling boosts our ability to discover immune cell types and states-and (b) antigen presentation and recognition-how models can be generated to predict presentation of antigens on MHC molecules and their detection by T cell receptors. These are two examples of a shift in how we find the drivers and targets of immunity, especially in the human system in the context of health and disease.

Intra- and Inter-cellular Rewiring of the Human Colon during Ulcerative Colitis.

Genome-wide association studies (GWAS) have revealed risk alleles for ulcerative colitis (UC). To understand their cell type specificities and pathways of action, we generate an atlas of 366,650 cells from the colon mucosa of 18 UC patients and 12 healthy individuals, revealing 51 epithelial, stromal, and immune cell subsets, including BEST4+ enterocytes, microfold-like cells, and IL13RA2+IL11+ inflammatory fibroblasts, which we associate with resistance to anti-TNF treatment. Inflammatory fibroblasts, inflammatory monocytes, microfold-like cells, and T cells that co-express CD8 and IL-17 expand with disease, forming intercellular interaction hubs. Many UC risk genes are cell type specific and co-regulated within relatively few gene modules, suggesting convergence onto limited sets of cell types and pathways. Using this observation, we nominate and infer functions for specific risk genes across GWAS loci. Our work provides a framework for interrogating complex human diseases and mapping risk variants to cell types and pathways.

Defining T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma.

Treatment of cancer has been revolutionized by immune checkpoint blockade therapies. Despite the high rate of response in advanced melanoma, the majority of patients succumb to disease. To identify factors associated with success or failure of checkpoint therapy, we profiled transcriptomes of 16,291 individual immune cells from 48 tumor samples of melanoma patients treated with checkpoint inhibitors. Two distinct states of CD8+ T cells were defined by clustering and associated with patient tumor regression or progression. A single transcription factor, TCF7, was visualized within CD8+ T cells in fixed tumor samples and predicted positive clinical outcome in an independent cohort of checkpoint-treated patients. We delineated the epigenetic landscape and clonality of these T cell states and demonstrated enhanced antitumor immunity by targeting novel combinations of factors in exhausted cells. Our study of immune cell transcriptomes from tumors demonstrates a strategy for identifying predictors, mechanisms, and targets for enhancing checkpoint immunotherapy.

The activation trajectory of plasmacytoid dendritic cells in vivo during a viral infection.

Plasmacytoid dendritic cells (pDCs) are a major source of type I interferon (IFN-I). What other functions pDCs exert in vivo during viral infections is controversial, and more studies are needed to understand their orchestration. In the present study, we characterize in depth and link pDC activation states in animals infected by mouse cytomegalovirus by combining Ifnb1 reporter mice with flow cytometry, single-cell RNA sequencing, confocal microscopy and a cognate CD4 T cell activation assay. We show that IFN-I production and T cell activation were performed by the same pDC, but these occurred sequentially in time and in different micro-anatomical locations. In addition, we show that pDC commitment to IFN-I production was marked early on by their downregulation of leukemia inhibitory factor receptor and was promoted by cell-intrinsic tumor necrosis factor signaling. We propose a new model for how individual pDCs are endowed to exert different functions in vivo during a viral infection, in a manner tightly orchestrated in time and space.

Plasmacytoid dendritic cells: Welcome back to the DC fold.

The presumed common origin of plasmacytoid and conventional dendritic cells has been the contentious subject of recent debate. In this issue of Immunity, Feng et al. employed an inducible cell barcoding system to track clonal relationships and uncovered a surprising close developmental relationship between cDC1s and pDCs.

A γδ T cell-IL-3 axis controls allergic responses through sensory neurons.

In naive individuals, sensory neurons directly detect and respond to allergens, leading to both the sensation of itch and the activation of local innate immune cells, which initiate the allergic immune response1,2. In the setting of chronic allergic inflammation, immune factors prime sensory neurons, causing pathologic itch3-7. Although these bidirectional neuroimmune circuits drive responses to allergens, whether immune cells regulate the set-point for neuronal activation by allergens in the naive state is unknown. Here we describe a γδ T cell-IL-3 signalling axis that controls the allergen responsiveness of cutaneous sensory neurons. We define a poorly characterized epidermal γδ T cell subset8, termed GD3 cells, that produces its hallmark cytokine IL-3 to promote allergic itch and the initiation of the allergic immune response. Mechanistically, IL-3 acts on Il3ra-expressing sensory neurons in a JAK2-dependent manner to lower their threshold for allergen activation without independently eliciting itch. This γδ T cell-IL-3 signalling axis further acts by means of STAT5 to promote neuropeptide production and the initiation of allergic immunity. These results reveal an endogenous immune rheostat that sits upstream of and governs sensory neuronal responses to allergens on first exposure. This pathway may explain individual differences in allergic susceptibility and opens new therapeutic avenues for treating allergic diseases.

Aryl Hydrocarbon Receptor Controls Monocyte Differentiation into Dendritic Cells versus Macrophages.

After entering tissues, monocytes differentiate into cells that share functional features with either macrophages or dendritic cells (DCs). How monocyte fate is directed toward monocyte-derived macrophages (mo-Macs) or monocyte-derived DCs (mo-DCs) and which transcription factors control these differentiation pathways remains unknown. Using an in vitro culture model yielding human mo-DCs and mo-Macs closely resembling those found in vivo in ascites, we show that IRF4 and MAFB were critical regulators of monocyte differentiation into mo-DCs and mo-Macs, respectively. Activation of the aryl hydrocarbon receptor (AHR) promoted mo-DC differentiation through the induction of BLIMP-1, while impairing differentiation into mo-Macs. AhR deficiency also impaired the in vivo differentiation of mouse mo-DCs. Finally, AHR activation correlated with mo-DC infiltration in leprosy lesions. These results establish that mo-DCs and mo-Macs are controlled by distinct transcription factors and show that AHR acts as a molecular switch for monocyte fate specification in response to micro-environmental factors.

Identification of regulators of the innate immune response to cytosolic DNA and retroviral infection by an integrative approach.

The innate immune system senses viral DNA that enters mammalian cells, or in aberrant situations self-DNA, and triggers type I interferon production. Here we present an integrative approach that combines quantitative proteomics, genomics and small molecule perturbations to identify genes involved in this pathway. We silenced 809 candidate genes, measured the response to dsDNA and connected resulting hits with the known signaling network. We identified ABCF1 as a critical protein that associates with dsDNA and the DNA-sensing components HMGB2 and IFI204. We also found that CDC37 regulates the stability of the signaling molecule TBK1 and that chemical inhibition of the CDC37-HSP90 interaction and several other pathway regulators potently modulates the innate immune response to DNA and retroviral infection.

Targeting the CBM complex causes Treg cells to prime tumours for immune checkpoint therapy.

Solid tumours are infiltrated by effector T cells with the potential to control or reject them, as well as by regulatory T (Treg) cells that restrict the function of effector T cells and thereby promote tumour growth1. The anti-tumour activity of effector T cells can be therapeutically unleashed, and is now being exploited for the treatment of some forms of human cancer. However, weak tumour-associated inflammatory responses and the immune-suppressive function of Treg cells remain major hurdles to broader effectiveness of tumour immunotherapy2. Here we show that, after disruption of the CARMA1-BCL10-MALT1 (CBM) signalosome complex, most tumour-infiltrating Treg cells produce IFNγ, resulting in stunted tumour growth. Notably, genetic deletion of both or even just one allele of CARMA1 (also known as Card11) in only a fraction of Treg cells-which avoided systemic autoimmunity-was sufficient to produce this anti-tumour effect, showing that it is not the mere loss of suppressive function but the gain of effector activity by Treg cells that initiates tumour control. The production of IFNγ by Treg cells was accompanied by activation of macrophages and upregulation of class I molecules of the major histocompatibility complex on tumour cells. However, tumour cells also upregulated the expression of PD-L1, which indicates activation of adaptive immune resistance3. Consequently, blockade of PD-1 together with CARMA1 deletion caused rejection of tumours that otherwise do not respond to anti-PD-1 monotherapy. This effect was reproduced by pharmacological inhibition of the CBM protein MALT1. Our results demonstrate that partial disruption of the CBM complex and induction of IFNγ secretion in the preferentially self-reactive Treg cell pool does not cause systemic autoimmunity but is sufficient to prime the tumour environment for successful immune checkpoint therapy.

Decoding human fetal liver haematopoiesis.

Definitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells and multipotent progenitors (HSC/MPPs) but remains poorly defined in humans. Here, using single-cell transcriptome profiling of approximately 140,000 liver and 74,000 skin, kidney and yolk sac cells, we identify the repertoire of human blood and immune cells during development. We infer differentiation trajectories from HSC/MPPs and evaluate the influence of the tissue microenvironment on blood and immune cell development. We reveal physiological erythropoiesis in fetal skin and the presence of mast cells, natural killer and innate lymphoid cell precursors in the yolk sac. We demonstrate a shift in the haemopoietic composition of fetal liver during gestation away from being predominantly erythroid, accompanied by a parallel change in differentiation potential of HSC/MPPs, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a reference for harnessing the therapeutic potential of HSC/MPPs.

IL-2 family cytokines IL-9 and IL-21 differentially regulate innate and adaptive type 2 immunity in asthma.

BACKGROUND: Asthma is often accompanied by type 2 immunity rich in IL-4, IL-5, and IL-13 cytokines produced by TH2 lymphocytes or type 2 innate lymphoid cells (ILC2s). IL-2 family cytokines play a key role in the differentiation, homeostasis, and effector function of innate and adaptive lymphocytes. OBJECTIVE: IL-9 and IL-21 boost activation and proliferation of TH2 and ILC2s, but the relative importance and potential synergism between these γ common chain cytokines are currently unknown. METHODS: Using newly generated antibodies, we inhibited IL-9 and IL-21 alone or in combination in various murine models of asthma. In a translational approach using segmental allergen challenge, we recently described elevated IL-9 levels in human subjects with allergic asthma compared with nonasthmatic controls. Here, we also measured IL-21 in both groups. RESULTS: IL-9 played a central role in controlling innate IL-33-induced lung inflammation by promoting proliferation and activation of ILC2s in an IL-21-independent manner. Conversely, chronic house dust mite-induced airway inflammation, mainly driven by adaptive immunity, was solely dependent on IL-21, which controlled TH2 activation, eosinophilia, total serum IgE, and formation of tertiary lymphoid structures. In a model of innate on adaptive immunity driven by papain allergen, a clear synergy was found between both pathways, as combined anti-IL-9 or anti-IL-21 blockade was superior in reducing key asthma features. In human bronchoalveolar lavage samples we measured elevated IL-21 protein within the allergic asthmatic group compared with the allergic control group. We also found increased IL21R transcripts and predicted IL-21 ligand activity in various disease-associated cell subsets. CONCLUSIONS: IL-9 and IL-21 play important and nonredundant roles in allergic asthma by boosting ILC2s and TH2 cells, revealing a dual IL-9 and IL-21 targeting strategy as a new and testable approach.

Gene-based association study of rare variants in children of diverse ancestries implicates TNFRSF21 in the development of allergic asthma.

BACKGROUND: Most genetic studies of asthma and allergy have focused on common variation in individuals primarily of European ancestry. Studying the role of rare variation in quantitative phenotypes and in asthma phenotypes in populations of diverse ancestries can provide additional, important insights into the development of these traits. OBJECTIVE: We sought to examine the contribution of rare variants to different asthma- or allergy-associated quantitative traits in children with diverse ancestries and explore their role in asthma phenotypes. METHODS: We examined whole-genome sequencing data from children participants in longitudinal studies of asthma (n = 1035; parent-identified as 67% Black and 25% Hispanic) to identify rare variants (minor allele frequency < 0.01). We assigned variants to genes and tested for associations using an omnibus variant-set test between each of 24,902 genes and 8 asthma-associated quantitative traits. On combining our results with external data on predicted gene expression in humans and mouse knockout studies, we identified 3 candidate genes. A burden of rare variants in each gene and in a combined 3-gene score was tested for its associations with clinical phenotypes of asthma. Finally, published single-cell gene expression data in lower airway mucosal cells after allergen challenge were used to assess transcriptional responses to allergen. RESULTS: Rare variants in USF1 were significantly associated with blood neutrophil count (P = 2.18 × 10-7); rare variants in TNFRSF21 with total IgE (P = 6.47 × 10-6) and PIK3R6 with eosinophil count (P = 4.10 × 10-5) reached suggestive significance. These 3 findings were supported by independent data from human and mouse studies. A burden of rare variants in TNFRSF21 and in a 3-gene score was associated with allergy-related phenotypes in cohorts of children with mild and severe asthma. Furthermore, TNFRSF21 was significantly upregulated in bronchial basal epithelial cells from adults with allergic asthma but not in adults with allergies (but not asthma) after allergen challenge. CONCLUSIONS: We report novel associations between rare variants in genes and allergic and inflammatory phenotypes in children with diverse ancestries, highlighting TNFRSF21 as contributing to the development of allergic asthma.

Corrigendum: Landscape of X chromosome inactivation across human tissues.

This corrects the article DOI: 10.1038/nature24265.

Altered ratio of dendritic cell subsets in skin-draining lymph nodes promotes Th2-driven contact hypersensitivity.

Plasmacytoid dendritic cells (pDCs) specialize in the production of type I IFN (IFN-I). pDCs can be depleted in vivo by injecting diphtheria toxin (DT) in a mouse in which pDCs express a diphtheria toxin receptor (DTR) transgene driven by the human CLEC4C promoter. This promoter is enriched for binding sites for TCF4, a transcription factor that promotes pDC differentiation and expression of pDC markers, including CLEC4C. Here, we found that injection of DT in CLEC4C-DTR+ mice markedly augmented Th2-dependent skin inflammation in a model of contact hypersensitivity (CHS) induced by the hapten fluorescein isothiocyanate. Unexpectedly, this biased Th2 response was independent of reduced IFN-I accompanying pDC depletion. In fact, DT treatment altered the representation of conventional dendritic cells (cDCs) in the skin-draining lymph nodes during the sensitization phase of CHS; there were fewer Th1-priming CD326+ CD103+ cDC1 and more Th2-priming CD11b+ cDC2. Single-cell RNA-sequencing of CLEC4C-DTR+ cDCs revealed that CD326+ DCs, like pDCs, expressed DTR and were depleted together with pDCs by DT treatment. Since CD326+ DCs did not express Tcf4, DTR expression might be driven by yet-undefined transcription factors activating the CLEC4C promoter. These results demonstrate that altered DC representation in the skin-draining lymph nodes during sensitization to allergens can cause Th2-driven CHS.

Cumulus provides cloud-based data analysis for large-scale single-cell and single-nucleus RNA-seq.

Massively parallel single-cell and single-nucleus RNA sequencing has opened the way to systematic tissue atlases in health and disease, but as the scale of data generation is growing, so is the need for computational pipelines for scaled analysis. Here we developed Cumulus-a cloud-based framework for analyzing large-scale single-cell and single-nucleus RNA sequencing datasets. Cumulus combines the power of cloud computing with improvements in algorithm and implementation to achieve high scalability, low cost, user-friendliness and integrated support for a comprehensive set of features. We benchmark Cumulus on the Human Cell Atlas Census of Immune Cells dataset of bone marrow cells and show that it substantially improves efficiency over conventional frameworks, while maintaining or improving the quality of results, enabling large-scale studies.

Landscape of X chromosome inactivation across human tissues.

X chromosome inactivation (XCI) silences transcription from one of the two X chromosomes in female mammalian cells to balance expression dosage between XX females and XY males. XCI is, however, incomplete in humans: up to one-third of X-chromosomal genes are expressed from both the active and inactive X chromosomes (Xa and Xi, respectively) in female cells, with the degree of 'escape' from inactivation varying between genes and individuals. The extent to which XCI is shared between cells and tissues remains poorly characterized, as does the degree to which incomplete XCI manifests as detectable sex differences in gene expression and phenotypic traits. Here we describe a systematic survey of XCI, integrating over 5,500 transcriptomes from 449 individuals spanning 29 tissues from GTEx (v6p release) and 940 single-cell transcriptomes, combined with genomic sequence data. We show that XCI at 683 X-chromosomal genes is generally uniform across human tissues, but identify examples of heterogeneity between tissues, individuals and cells. We show that incomplete XCI affects at least 23% of X-chromosomal genes, identify seven genes that escape XCI with support from multiple lines of evidence and demonstrate that escape from XCI results in sex biases in gene expression, establishing incomplete XCI as a mechanism that is likely to introduce phenotypic diversity. Overall, this updated catalogue of XCI across human tissues helps to increase our understanding of the extent and impact of the incompleteness in the maintenance of XCI.

Alveolar, Endothelial, and Organ Injury Marker Dynamics in Severe COVID-19.

Rationale: Alveolar and endothelial injury may be differentially associated with coronavirus disease (COVID-19) severity over time. Objectives: To describe alveolar and endothelial injury dynamics and associations with COVID-19 severity, cardiorenovascular injury, and outcomes. Methods: This single-center observational study enrolled patients with COVID-19 requiring respiratory support at emergency department presentation. More than 40 markers of alveolar (including receptor for advanced glycation endproducts [RAGE]), endothelial (including angiopoietin-2), and cardiorenovascular injury (including renin, kidney injury molecule-1, and troponin-I) were serially compared between invasively and spontaneously ventilated patients using mixed-effects repeated-measures models. Ventilatory ratios were calculated for intubated patients. Associations of biomarkers with modified World Health Organization scale at Day 28 were determined with multivariable proportional-odds regression. Measurements and Main Results: Of 225 patients, 74 (33%) received invasive ventilation at Day 0. RAGE was 1.80-fold higher in invasive ventilation patients at Day 0 (95% confidence interval [CI], 1.50-2.17) versus spontaneous ventilation, but decreased over time in all patients. Changes in alveolar markers did not correlate with changes in endothelial, cardiac, or renal injury markers. In contrast, endothelial markers were similar to lower at Day 0 for invasive ventilation versus spontaneous ventilation, but then increased over time only among intubated patients. In intubated patients, angiopoietin-2 was similar (fold difference, 1.02; 95% CI, 0.89-1.17) to nonintubated patients at Day 0 but 1.80-fold higher (95% CI, 1.56-2.06) at Day 3; cardiorenovascular injury markers showed similar patterns. Endothelial markers were not consistently associated with ventilatory ratios. Endothelial markers were more often significantly associated with 28-day outcomes than alveolar markers. Conclusions: Alveolar injury markers increase early. Endothelial injury markers increase later and are associated with cardiorenovascular injury and 28-day outcome. Alveolar and endothelial injury likely contribute at different times to disease progression in severe COVID-19.

Single-cell immunophenotyping of the fetal immune response to maternal SARS-CoV-2 infection in late gestation.

BACKGROUND: During the COVID-19 pandemic, thousands of pregnant women have been infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The implications of maternal SARS-CoV-2 infection on fetal and childhood well-being need to be characterized. We aimed to characterize the fetal immune response to maternal SARS-CoV-2 infection. METHODS: We performed single-cell RNA-sequencing and T cell receptor sequencing on cord blood mononuclear cells (CBMCs) from newborns of mothers infected with SARS-CoV-2 in the third trimester (cases) or without SARS-CoV-2 infection (controls). RESULTS: We identified widespread gene expression changes in CBMCs from cases, including upregulation of interferon-stimulated genes and major histocompatibility complex genes in CD14+ monocytes, transcriptional changes suggestive of activation of plasmacytoid dendritic cells, and activation and exhaustion of natural killer cells. Lastly, we observed fetal T cell clonal expansion in cases compared to controls. CONCLUSIONS: As none of the infants were infected with SARS-CoV-2, our results suggest that maternal SARS-CoV-2 infection might modulate the fetal immune system in the absence of vertical transmission. IMPACT: The implications of maternal SARS-CoV-2 infection in the absence of vertical transmission on fetal and childhood well-being are poorly understood. Maternal SARS-CoV-2 infection might modulate the fetal immune system in the absence of vertical transmission. This study raises important questions about the untoward effects of maternal SARS-CoV-2 on the fetus, even in the absence of vertical transmission.