Single-cell analysis reveals lasting immunological consequences of influenza infection and respiratory immunization in the pig lung.

The pig is a natural host for influenza viruses and integrally involved in virus evolution through interspecies transmissions between humans and swine. Swine have many physiological, anatomical, and immunological similarities to humans, and are an excellent model for human influenza. Here, we employed single cell RNA-sequencing (scRNA-seq) and flow cytometry to characterize the major leukocyte subsets in bronchoalveolar lavage (BAL), twenty-one days after H1N1pdm09 infection or respiratory immunization with an adenoviral vector vaccine expressing hemagglutinin and nucleoprotein with or without IL-1ß. Mapping scRNA-seq clusters from BAL onto those previously described in peripheral blood facilitated annotation and highlighted differences between tissue resident and circulating immune cells. ScRNA-seq data and functional assays revealed lasting impacts of immune challenge on BAL populations. First, mucosal administration of IL-1ß reduced the number of functionally active Treg cells. Second, influenza infection upregulated IFI6 in BAL cells and decreased their susceptibility to virus replication in vitro. Our data provide a reference map of porcine BAL cells and reveal lasting immunological consequences of influenza infection and respiratory immunization in a highly relevant large animal model for respiratory virus infection.

A metabolic pacer ensures smooth running of the lymphocyte activation race.

Upon lymphocyte stimulation, accumulation of intracellular NAD(H) reflects the strength of antigen receptor signals and controls the rate of cell cycle entry and proliferation (see related Research Article by Turner et al.).

Cytotoxic T lymphocytes require transcription for infiltration but not target cell lysis.

Effector cytotoxic T lymphocytes (CTLs) are critical for ridding the body of infected or cancerous cells. CTL T cell receptor (TCR) ligation not only drives the delivery and release of cytolytic granules but also initiates a new wave of transcription. In order to address whether TCR-induced transcriptomic changes impact the ability of CTLs to kill, we asked which genes are expressed immediately after CTLs encounter targets and how CTL responses change when inhibiting transcription. Our data demonstrate that while expression of cytokines/chemokines and transcriptional machinery depend on transcription, cytotoxic protein expression and cytolytic activity are relatively robust to transcription blockade, with CTLs lysing nearby target cells for several hours after actinomycin D treatment. Monitoring CTL movement among target cells after inhibiting transcription demonstrates an infiltration defect that is not rectified by provision of exogenous cytokine/chemokine gradients, indicating a cell-intrinsic transcriptional requirement for infiltration. Together, our results reveal differential molecular control of CTL functions, separating recruitment and infiltration from cytolysis.

ChAdOx1 nCoV-19 vaccination generates spike-specific CD8+ T cells in aged mice.

Effective vaccines have reduced the morbidity and mortality caused by severe acute respiratory syndrome coronavirus-2 infection; however, the elderly remain the most at risk. Understanding how vaccines generate protective immunity and how these mechanisms change with age is key for informing future vaccine design. Cytotoxic CD8+ T cells are important for killing virally infected cells, and vaccines that induce antigen-specific CD8+ T cells in addition to humoral immunity provide an extra layer of immune protection. This is particularly important in cases where antibody titers are suboptimal, as can occur in older individuals. Here, we show that in aged mice, spike epitope-specific CD8+ T cells are generated in comparable numbers to younger animals after ChAdOx1 nCoV-19 vaccination, although phenotypic differences exist. This demonstrates that ChAdOx1 nCoV-19 elicits a good CD8+ T-cell response in older bodies, but that typical age-associated features are evident on these vaccine reactive T cells.

Super-killer CTLs are generated by single gene deletion of Bach2.

Bach2 codes for a transcriptional regulator exerting major influences on T cell-mediated immune regulation. Effector CTLs derived from in vitro activation of murine CD8+ T cells showed increased proliferative and cytolytic capacity in the absence of BACH2. Before activation, BACH2-deficient splenic CD8+ T cells had a higher abundance of memory and reduced abundance of naïve cells compared to wild-type. CTLs derived from central memory T cells were more potently cytotoxic than those derived from naïve T cells, but even within separated subsets, BACH2-deficiency conferred a cytotoxic advantage. Immunofluorescence and electron microscopy revealed larger granules in BACH2-deficient compared to wild-type CTLs, and proteomic analysis showed an increase in granule content, including perforin and granzymes. Thus, the enhanced cytotoxicity observed in effector CTLs lacking BACH2 arises not only from differences in their initial differentiation state but also inherent production of enlarged cytolytic granules. These results demonstrate how a single gene deletion can produce a CTL super-killer.

Divide and Conquer: Phenotypic and Temporal Heterogeneity Within CD8+ T Cell Responses.

The advent of technologies that can characterize the phenotypes, functions and fates of individual cells has revealed extensive and often unexpected levels of diversity between cells that are nominally of the same subset. CD8+ T cells, also known as cytotoxic T lymphocytes (CTLs), are no exception. Investigations of individual CD8+ T cells both in vitro and in vivo have highlighted the heterogeneity of cellular responses at the levels of activation, differentiation and function. This review takes a broad perspective on the topic of heterogeneity, outlining different forms of variation that arise during a CD8+ T cell response. Specific attention is paid to the impact of T cell receptor (TCR) stimulation strength on heterogeneity. In particular, this review endeavors to highlight connections between variation at different cellular stages, presenting known mechanisms and key open questions about how variation between cells can arise and propagate.

Staggered starts in the race to T cell activation.

How T lymphocytes tune their responses to different strengths of stimulation is a fundamental question in immunology. Recent work using new optogenetic, single-cell genomic, and live-imaging approaches has revealed that stimulation strength controls the rate of individual cell responses within a population. Moreover, these responses have been found to use shared molecular programs, regardless of stimulation strength. However, additional data indicate that stimulation duration or cytokine feedback can impact later gene expression phenotypes of activated cells. In-depth molecular studies have suggested mechanisms by which stimulation strength might modulate the probability of T cell activation. This emerging model allows activating T cells to achieve a wide range of population responses through probabilistic control within individual cells.

Longitudinal proteomic profiling of dialysis patients with COVID-19 reveals markers of severity and predictors of death.

End-stage kidney disease (ESKD) patients are at high risk of severe COVID-19. We measured 436 circulating proteins in serial blood samples from hospitalised and non-hospitalised ESKD patients with COVID-19 (n = 256 samples from 55 patients). Comparison to 51 non-infected patients revealed 221 differentially expressed proteins, with consistent results in a separate subcohort of 46 COVID-19 patients. Two hundred and three proteins were associated with clinical severity, including IL6, markers of monocyte recruitment (e.g. CCL2, CCL7), neutrophil activation (e.g. proteinase-3), and epithelial injury (e.g. KRT19). Machine-learning identified predictors of severity including IL18BP, CTSD, GDF15, and KRT19. Survival analysis with joint models revealed 69 predictors of death. Longitudinal modelling with linear mixed models uncovered 32 proteins displaying different temporal profiles in severe versus non-severe disease, including integrins and adhesion molecules. These data implicate epithelial damage, innate immune activation, and leucocyte-endothelial interactions in the pathology of severe COVID-19 and provide a resource for identifying drug targets.

COVID-19 varies from a mild illness in some people to fatal disease in others. Patients with severe disease tend to be older and have underlying medical problems. People with kidney failure have a particularly high risk of developing severe or fatal COVID-19. Patients with severe COVID-19 have high levels of inflammation, causing damage to tissues around the body. Many drugs that target inflammation have already been developed for other diseases. Therefore, to repurpose existing drugs or design new treatments, it is important to determine which proteins drive inflammation in COVID-19. Here, Gisby, Clarke, Medjeral-Thomas et al. measured 436 proteins in the blood of patients with kidney failure and compared the levels between patients who had COVID-19 to those who did not. This revealed that patients with COVID-19 had increased levels of hundreds of proteins involved in inflammation and tissue injury. Using a combination of statistical and machine learning analyses, Gisby et al. probed the data for proteins that might predict a more severe disease progression. In total, over 200 proteins were linked to disease severity, and 69 with increased risk of death. Tracking how levels of blood proteins changed over time revealed further differences between mild and severe disease. Comparing this data with a similar study of COVID-19 in people without kidney failure showed many similarities. This suggests that the findings may apply to COVID-19 patients more generally. Identifying the proteins that are a cause of severe COVID-19 ­ rather than just correlated with it ­ is an important next step that could help to select new drugs for severe COVID-19.

ILC2-driven innate immune checkpoint mechanism antagonizes NK cell antimetastatic function in the lung.

Metastasis constitutes the primary cause of cancer-related deaths, with the lung being a commonly affected organ. We found that activation of lung-resident group 2 innate lymphoid cells (ILC2s) orchestrated suppression of natural killer (NK) cell-mediated innate antitumor immunity, leading to increased lung metastases and mortality. Using multiple models of lung metastasis, we show that interleukin (IL)-33-dependent ILC2 activation in the lung is involved centrally in promoting tumor burden. ILC2-driven innate type 2 inflammation is accompanied by profound local suppression of interferon-γ production and cytotoxic function of lung NK cells. ILC2-dependent suppression of NK cells is elaborated via an innate regulatory mechanism, which is reliant on IL-5-induced lung eosinophilia, ultimately limiting the metabolic fitness of NK cells. Therapeutic targeting of IL-33 or IL-5 reversed NK cell suppression and alleviated cancer burden. Thus, we reveal an important function of IL-33 and ILC2s in promoting tumor metastasis via their capacity to suppress innate type 1 immunity.

Stimulation strength controls the rate of initiation but not the molecular organisation of TCR-induced signalling.

Millions of naïve T cells with different TCRs may interact with a peptide-MHC ligand, but very few will activate. Remarkably, this fine control is orchestrated using a limited set of intracellular machinery. It remains unclear whether changes in stimulation strength alter the programme of signalling events leading to T cell activation. Using mass cytometry to simultaneously measure multiple signalling pathways during activation of murine CD8+ T cells, we found a programme of distal signalling events that is shared, regardless of the strength of TCR stimulation. Moreover, the relationship between transcription of early response genes Nr4a1 and Irf8 and activation of the ribosomal protein S6 is also conserved across stimuli. Instead, we found that stimulation strength dictates the rate with which cells initiate signalling through this network. These data suggest that TCR-induced signalling results in a coordinated activation program, modulated in rate but not organization by stimulation strength.

Amongst the different types of cells the body uses to protect itself, killer T cells have an unique role: they can detect and neutralize cells that have been become dangerous for the organism ­ for example, cells which are cancerous or hijacked by viruses. In a healthy organism, T cells circulate through the body in an inactivated state. When a disease emerges, receptors at the surface of the cells can detect elements coming from harmful agents; this stimulation then triggers a molecular cascade inside the T cell that leads to activation. This system is relatively simple, pairing a finite number of receptors with a limited set of internal components. At the same time, the activity of T cells is finely regulated, and their activation tightly controlled: they must kill enough cells to stop the illness without causing excess damage. How this is accomplished is still unclear. A T cell can recognize harmful agents that bind its receptors with differing strengths, but how this variability in stimulation strength affects the signaling processes within the cell is still poorly understood. To investigate this question, Ma et al. used an approach called mass cytometry and analyzed the internal processes of mouse killer T cells receiving different strengths of stimulation. This investigation revealed little change in the patterns of signaling in response to signals of different strength. Instead, what differed was the proportion of T cells that became activated, and how fast this process took place: stronger stimulations led to a larger population of killer T cells being activated more rapidly. Overall, this work sheds light on how killer T cells fine-tune their response to illness using only a simple system to control their activation.

Circulating TNF-like protein 1A (TL1A) is elevated early in rheumatoid arthritis and depends on TNF.

BACKGROUND: The tumor necrosis factor (TNF) superfamily cytokine TNF-like protein 1A (TL1A) and its receptor DR3 are essential for diverse animal models of autoimmune disease and may be pathogenic in rheumatoid arthritis (RA). However, the relationship of TL1A to disease duration, activity, and response to anti-TNF and other therapies in RA is not clear. METHODS: We measured soluble TL1A in synovial fluid (SF), serum, or plasma from RA first-degree relatives (FDRs) and in early RA and established disease. We measured the effects of anti-TNF and methotrexate (MTX) therapy on circulating TL1A from multiple independent RA treatment trials. We also determined the ability of a blocking anti-TL1A antibody to inhibit clinical disease and articular bone destruction in the murine collagen-induced arthritis (CIA) model of human RA. RESULTS: Soluble TL1A was specifically elevated in the blood and SF of patients with RA compared to patients with other diseases and was elevated early in disease and in at-risk anti-cyclic citrullinated peptide (CCP) (+) first-degree relatives (FDRs). Therapeutic TNF inhibition reduced serum TL1A in both responders and non-responders, whereas TL1A declined following MTX treatment only in responders. In murine CIA, TL1A blockade was clinically efficacious and reduced bone erosions. CONCLUSIONS: TL1A is specifically elevated in RA from early in the disease course and in at-risk FDRs. The decline in TL1A after TNF blockade suggests that TL1A levels may be a useful biomarker for TNF activity in RA. These results support the further investigation of the relationship between TL1A and TNF and TL1A blockade as a potential therapeutic strategy in RA.

PIP5 Kinases Regulate Membrane Phosphoinositide and Actin Composition for Targeted Granule Secretion by Cytotoxic Lymphocytes.

How cytotoxic T lymphocytes (CTLs) sense T cell receptor (TCR) signaling in order to specialize an area of plasma membrane for granule secretion is not understood. Here, we demonstrate that immune synapse formation led to rapid localized changes in the phosphoinositide composition of the plasma membrane, both reducing phosphoinositide-4-phosphate (PI(4)P), PI(4,5)P2, and PI(3,4,5)P3 and increasing diacylglycerol (DAG) and PI(3,4)P2 within the first 2 min of synapse formation. These changes reduced negative charge across the synapse, triggering the release of electrostatically bound PIP5 kinases that are required to replenish PI(4,5)P2. As PI(4,5)P2 decreased, actin was depleted from the membrane, allowing secretion. Forced localization of PIP5Kß across the synapse prevented actin depletion, blocking both centrosome docking and secretion. Thus, PIP5Ks act as molecular sensors of TCR activation, controlling actin recruitment across the synapse, ensuring exquisite co-ordination between TCR signaling and CTL secretion.

Correcting the Mean-Variance Dependency for Differential Variability Testing Using Single-Cell RNA Sequencing Data.

Cell-to-cell transcriptional variability in otherwise homogeneous cell populations plays an important role in tissue function and development. Single-cell RNA sequencing can characterize this variability in a transcriptome-wide manner. However, technical variation and the confounding between variability and mean expression estimates hinder meaningful comparison of expression variability between cell populations. To address this problem, we introduce an analysis approach that extends the BASiCS statistical framework to derive a residual measure of variability that is not confounded by mean expression. This includes a robust procedure for quantifying technical noise in experiments where technical spike-in molecules are not available. We illustrate how our method provides biological insight into the dynamics of cell-to-cell expression variability, highlighting a synchronization of biosynthetic machinery components in immune cells upon activation. In contrast to the uniform up-regulation of the biosynthetic machinery, CD4+ T cells show heterogeneous up-regulation of immune-related and lineage-defining genes during activation and differentiation.

Reduced monocyte and macrophage TNFSF15/TL1A expression is associated with susceptibility to inflammatory bowel disease.

Chronic inflammation in inflammatory bowel disease (IBD) results from a breakdown of intestinal immune homeostasis and compromise of the intestinal barrier. Genome-wide association studies have identified over 200 genetic loci associated with risk for IBD, but the functional mechanisms of most of these genetic variants remain unknown. Polymorphisms at the TNFSF15 locus, which encodes the TNF superfamily cytokine commonly known as TL1A, are associated with susceptibility to IBD in multiple ethnic groups. In a wide variety of murine models of inflammation including models of IBD, TNFSF15 promotes immunopathology by signaling through its receptor DR3. Such evidence has led to the hypothesis that expression of this lymphocyte costimulatory cytokine increases risk for IBD. In contrast, here we show that the IBD-risk haplotype at TNFSF15 is associated with decreased expression of the gene by peripheral blood monocytes in both healthy volunteers and IBD patients. This association persists under various stimulation conditions at both the RNA and protein levels and is maintained after macrophage differentiation. Utilizing a "recall-by-genotype" bioresource for allele-specific expression measurements in a functional fine-mapping assay, we localize the polymorphism controlling TNFSF15 expression to the regulatory region upstream of the gene. Through a T cell costimulation assay, we demonstrate that genetically regulated TNFSF15 has functional relevance. These findings indicate that genetically enhanced expression of TNFSF15 in specific cell types may confer protection against the development of IBD.

Detection and removal of barcode swapping in single-cell RNA-seq data.

Barcode swapping results in the mislabelling of sequencing reads between multiplexed samples on patterned flow-cell Illumina sequencing machines. This may compromise the validity of numerous genomic assays; however, the severity and consequences of barcode swapping remain poorly understood. We have used two statistical approaches to robustly quantify the fraction of swapped reads in two plate-based single-cell RNA-sequencing datasets. We found that approximately 2.5% of reads were mislabelled between samples on the HiSeq 4000, which is lower than previous reports. We observed no correlation between the swapped fraction of reads and the concentration of free barcode across plates. Furthermore, we have demonstrated that barcode swapping may generate complex but artefactual cell libraries in droplet-based single-cell RNA-sequencing studies. To eliminate these artefacts, we have developed an algorithm to exclude individual molecules that have swapped between samples in 10x Genomics experiments, allowing the continued use of cutting-edge sequencing machines for these assays.

T cell cytolytic capacity is independent of initial stimulation strength.

How cells respond to myriad stimuli with finite signaling machinery is central to immunology. In naive T cells, the inherent effect of ligand strength on activation pathways and endpoints has remained controversial, confounded by environmental fluctuations and intercellular variability within populations. Here we studied how ligand potency affected the activation of CD8+ T cells in vitro, through the use of genome-wide RNA, multi-dimensional protein and functional measurements in single cells. Our data revealed that strong ligands drove more efficient and uniform activation than did weak ligands, but all activated cells were fully cytolytic. Notably, activation followed the same transcriptional pathways regardless of ligand potency. Thus, stimulation strength did not intrinsically dictate the T cell-activation route or phenotype; instead, it controlled how rapidly and simultaneously the cells initiated activation, allowing limited machinery to elicit wide-ranging responses.

Cleavage of TL1A Differentially Regulates Its Effects on Innate and Adaptive Immune Cells.

TNF superfamily cytokines play major roles in the regulation of adaptive and innate immunity. The TNF superfamily cytokine TL1A (TNFSF15), through its cognate receptor DR3 (TNFRSF25), promotes T cell immunity to pathogens and directly costimulates group 2 and 3 innate lymphoid cells. Polymorphisms in the TNFSF15 gene are associated with the risk for various human diseases, including inflammatory bowel disease. Like other cytokines in the TNF superfamily, TL1A is synthesized as a type II transmembrane protein and cleaved from the plasma membrane by metalloproteinases. Membrane cleavage has been shown to alter or abrogate certain activities of other TNF family cytokines; however, the functional capabilities of membrane-bound and soluble forms TL1A are not known. Constitutive expression of TL1A in transgenic mice results in expansion of activated T cells and promotes intestinal hyperplasia and inflammation through stimulation of group 2 innate lymphoid cells. Through the generation of membrane-restricted TL1A-transgenic mice, we demonstrate that membrane TL1A promotes expression of inflammatory cytokines in the lung, dependent upon DR3 expression on T cells. Soluble TL1A alone was unable to produce this phenotype but was still able to induce intestinal type 2 inflammation independently of T cells. These data suggest differential roles for membrane and soluble TL1A on adaptive and innate immune cells and have implications for the consequences of blocking these two forms of TL1A.

Testing for differential abundance in mass cytometry data.

When comparing biological conditions using mass cytometry data, a key challenge is to identify cellular populations that change in abundance. Here, we present a computational strategy for detecting 'differentially abundant' populations by assigning cells to hyperspheres, testing for significant differences between conditions and controlling the spatial false discovery rate. Our method (http://bioconductor.org/packages/cydar) outperforms other approaches in simulations and finds novel patterns of differential abundance in real data.

Targeted genomic analysis reveals widespread autoimmune disease association with regulatory variants in the TNF superfamily cytokine signalling network.

BACKGROUND: Tumour necrosis factor (TNF) superfamily cytokines and their receptors regulate diverse immune system functions through a common set of signalling pathways. Genetic variants in and expression of individual TNF superfamily cytokines, receptors and signalling proteins have been associated with autoimmune and inflammatory diseases, but their interconnected biology has been largely unexplored. METHODS: We took a hypothesis-driven approach using available genome-wide datasets to identify genetic variants regulating gene expression in the TNF superfamily cytokine signalling network and the association of these variants with autoimmune and autoinflammatory disease. Using paired gene expression and genetic data, we identified genetic variants associated with gene expression, expression quantitative trait loci (eQTLs), in four peripheral blood cell subsets. We then examined whether eQTLs were dependent on gene expression level or the presence of active enhancer chromatin marks. Using these eQTLs as genetic markers of the TNF superfamily signalling network, we performed targeted gene set association analysis in eight autoimmune and autoinflammatory disease genome-wide association studies. RESULTS: Comparison of TNF superfamily network gene expression and regulatory variants across four leucocyte subsets revealed patterns that differed between cell types. eQTLs for genes in this network were not dependent on absolute gene expression levels and were not enriched for chromatin marks of active enhancers. By examining autoimmune disease risk variants among our eQTLs, we found that risk alleles can be associated with either increased or decreased expression of co-stimulatory TNF superfamily cytokines, receptors or downstream signalling molecules. Gene set disease association analysis revealed that eQTLs for genes in the TNF superfamily pathway were associated with six of the eight autoimmune and autoinflammatory diseases examined, demonstrating associations beyond single genome-wide significant hits. CONCLUSIONS: This systematic analysis of the influence of regulatory genetic variants in the TNF superfamily network reveals widespread and diverse roles for these cytokines in susceptibility to a number of immune-mediated diseases.