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.

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.

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.

Temporal Trends and Outcomes Among Patients Admitted for Immune-Related Adverse Events: A Single-Center Retrospective Cohort Study from 2011 to 2018.

BACKGROUND: The aim of this study was to characterize severe immune-related adverse events (irAEs) seen among hospitalized patients and to examine risk factors for irAE admissions and clinically relevant outcomes, including length of stay, immune checkpoint inhibitor (ICI) discontinuation, readmission, and death. METHODS: Patients who received ICI therapy (ipilimumab, pembrolizumab, nivolumab, atezolizumab, durvalumab, avelumab, or any ICI combination) at Massachusetts General Hospital (MGH) and were hospitalized at MGH following ICI initiation between January 1, 2011, and October 24, 2018, were identified using pharmacy and hospital admission databases. Medical records of all irAE admissions were reviewed, and specialist review with defined criteria was performed. Demographic data, relevant clinical history (malignancy type and most recent ICI regimen), and key admission characteristics, including dates of admission and discharge, immunosuppressive management, ICI discontinuation, readmission, and death, were collected. RESULTS: In total, 450 admissions were classified as irAE admissions and represent the study's cohort. Alongside the increasing use of ICIs at our institution, the number of patients admitted to MGH for irAEs has gradually increased every year from 9 in 2011 to 92 in 2018. The hospitalization rate per ICI recipient has declined over that same time period (25.0% in 2011 to 8.5% in 2018). The most common toxicities leading to hospitalization in our cohort were gastrointestinal (30.7%; n = 138), pulmonary (15.8%; n = 71), hepatic (14.2%; n = 64), endocrine (12.2%; n = 55), neurologic (8.4%; n = 38), cardiac (6.7%; n = 30), and dermatologic (4.4%; n = 20). Multivariable logistic regression revealed statistically significant increases in irAE admission risk for CTLA-4 monotherapy recipients (odds ratio [OR], 2.02; p < .001) and CTLA-4 plus PD-1 combination therapy recipients (OR, 1.88; p < .001), relative to PD-1/PD-L1 monotherapy recipients, and patients with multiple toxicity had a 5-fold increase in inpatient mortality. CONCLUSION: This study illustrates that cancer centers must be prepared to manage a wide variety of irAE types and that CTLA-4 and combination ICI regimens are more likely to cause irAE admissions, and earlier. In addition, admissions for patients with multi-organ involvement is common and those patients are at highest risk of inpatient mortality. IMPLICATIONS FOR PRACTICE: The number of patients admitted to Massachusetts General Hospital for immune-related adverse events (irAEs) has gradually increased every year and the most common admissions are for gastrointestinal (30.7%), pulmonary (15/8%), and hepatic (14.2%) events. Readmission rates are high (29% at 30 days, 49% at 180 days) and 64.2% have to permanently discontinue immune checkpoint inhibitor therapy. Importantly, multiple concurrent toxicities were seen in 21.6% (97/450) of irAE admissions and these patients have a fivefold increased risk of inpatient death.

Genetic analysis of isoform usage in the human anti-viral response reveals influenza-specific regulation of ERAP2 transcripts under balancing selection.

While genetic variants are known to be associated with overall gene abundance in stimulated immune cells, less is known about their effects on alternative isoform usage. By analyzing RNA-seq profiles of monocyte-derived dendritic cells from 243 individuals, we uncovered thousands of unannotated isoforms synthesized in response to influenza infection and type 1 interferon stimulation. We identified more than a thousand quantitative trait loci (QTLs) associated with alternate isoform usage (isoQTLs), many of which are independent of expression QTLs (eQTLs) for the same gene. Compared with eQTLs, isoQTLs are enriched for splice sites and untranslated regions, but depleted of sequences upstream of annotated transcription start sites. Both eQTLs and isoQTLs explain a significant proportion of the disease heritability attributed to common genetic variants. At the ERAP2 locus, we shed light on the function of the gene and how two frequent, highly differentiated haplotypes with intermediate frequencies could be maintained by balancing selection. At baseline and following type 1 interferon stimulation, the major haplotype is associated with low ERAP2 expression caused by nonsense-mediated decay, while the minor haplotype, known to increase Crohn's disease risk, is associated with high ERAP2 expression. In response to influenza infection, we found two uncharacterized isoforms expressed from the major haplotype, likely the result of multiple perfectly linked variants affecting the transcription and splicing at the locus. Thus, genetic variants at a single locus could modulate independent gene regulatory processes in innate immune responses and, in the case of ERAP2, may confer a historical fitness advantage in response to virus.