High-throughput phenotyping reveals expansive genetic and structural underpinnings of immune variation.

By developing a high-density murine immunophenotyping platform compatible with high-throughput genetic screening, we have established profound contributions of genetics and structure to immune variation (http://www.immunophenotype.org). Specifically, high-throughput phenotyping of 530 unique mouse gene knockouts identified 140 monogenic 'hits', of which most had no previous immunologic association. Furthermore, hits were collectively enriched in genes for which humans show poor tolerance to loss of function. The immunophenotyping platform also exposed dense correlation networks linking immune parameters with each other and with specific physiologic traits. Such linkages limit freedom of movement for individual immune parameters, thereby imposing genetically regulated 'immunologic structures', the integrity of which was associated with immunocompetence. Hence, we provide an expanded genetic resource and structural perspective for understanding and monitoring immune variation in health and disease.

Epithelia Use Butyrophilin-like Molecules to Shape Organ-Specific γδ T Cell Compartments.

Many body surfaces harbor organ-specific γδ T cell compartments that contribute to tissue integrity. Thus, murine dendritic epidermal T cells (DETCs) uniquely expressing T cell receptor (TCR)-Vγ5 chains protect from cutaneous carcinogens. The DETC repertoire is shaped by Skint1, a butyrophilin-like (Btnl) gene expressed specifically by thymic epithelial cells and suprabasal keratinocytes. However, the generality of this mechanism has remained opaque, since neither Skint1 nor DETCs are evolutionarily conserved. Here, Btnl1 expressed by murine enterocytes is shown to shape the local TCR-Vγ7(+) γδ compartment. Uninfluenced by microbial or food antigens, this activity evokes the developmental selection of TCRαß(+) repertoires. Indeed, Btnl1 and Btnl6 jointly induce TCR-dependent responses specifically in intestinal Vγ7(+) cells. Likewise, human gut epithelial cells express BTNL3 and BTNL8 that jointly induce selective TCR-dependent responses of human colonic Vγ4(+) cells. Hence, a conserved mechanism emerges whereby epithelia use organ-specific BTNL/Btnl genes to shape local T cell compartments.

Neutralization potency of monoclonal antibodies recognizing dominant and subdominant epitopes on SARS-CoV-2 Spike is impacted by the B.1.1.7 variant.

Interaction of the SARS-CoV-2 Spike receptor binding domain (RBD) with the receptor ACE2 on host cells is essential for viral entry. RBD is the dominant target for neutralizing antibodies, and several neutralizing epitopes on RBD have been molecularly characterized. Analysis of circulating SARS-CoV-2 variants has revealed mutations arising in the RBD, N-terminal domain (NTD) and S2 subunits of Spike. To understand how these mutations affect Spike antigenicity, we isolated and characterized >100 monoclonal antibodies targeting epitopes on RBD, NTD, and S2 from SARS-CoV-2-infected individuals. Approximately 45% showed neutralizing activity, of which ∼20% were NTD specific. NTD-specific antibodies formed two distinct groups: the first was highly potent against infectious virus, whereas the second was less potent and displayed glycan-dependant neutralization activity. Mutations present in B.1.1.7 Spike frequently conferred neutralization resistance to NTD-specific antibodies. This work demonstrates that neutralizing antibodies targeting subdominant epitopes should be considered when investigating antigenic drift in emerging variants.

Prolidase Deficiency Causes Spontaneous T Cell Activation and Lupus-like Autoimmunity.

Prolidase deficiency (PD) is a multisystem disorder caused by mutations in the PEPD gene, which encodes a ubiquitously expressed metallopeptidase essential for the hydrolysis of dipeptides containing C-terminal proline or hydroxyproline. PD typically presents in childhood with developmental delay, skin ulcers, recurrent infections, and, in some patients, autoimmune features that can mimic systemic lupus erythematosus. The basis for the autoimmune association is uncertain, but might be due to self-antigen exposure with tissue damage, or indirectly driven by chronic infection and microbial burden. In this study, we address the question of causation and show that Pepd-null mice have increased antinuclear autoantibodies and raised serum IgA, accompanied by kidney immune complex deposition, consistent with a systemic lupus erythematosus-like disease. These features are associated with an accumulation of CD4 and CD8 effector T cells in the spleen and liver. Pepd deficiency leads to spontaneous T cell activation and proliferation into the effector subset, which is cell intrinsic and independent of Ag receptor specificity or antigenic stimulation. However, an increase in KLRG1+ effector CD8 cells is not observed in mixed chimeras, in which the autoimmune phenotype is also absent. Our findings link autoimmune susceptibility in PD to spontaneous T cell dysfunction, likely to be acting in combination with immune activators that lie outside the hemopoietic system but result from the abnormal metabolism or loss of nonenzymatic prolidase function. This knowledge provides insight into the role of prolidase in the maintenance of self-tolerance and highlights the importance of treatment to control T cell activation.

Global patterns of antigen receptor repertoire disruption across adaptive immune compartments in COVID-19.

Whereas pathogen-specific T and B cells are a primary focus of interest during infectious disease, we have used COVID-19 to ask whether their emergence comes at a cost of broader B cell and T cell repertoire disruption. We applied a genomic DNA-based approach to concurrently study the immunoglobulin-heavy (IGH) and T cell receptor (TCR) ß and δ chain loci of 95 individuals. Our approach detected anticipated repertoire focusing for the IGH repertoire, including expansions of clusters of related sequences temporally aligned with SARS-CoV-2-specific seroconversion, and enrichment of some shared SARS-CoV-2-associated sequences. No significant age-related or disease severity-related deficiencies were noted for the IGH repertoire. By contrast, whereas focusing occurred at the TCRß and TCRδ loci, including some TCRß sequence-sharing, disruptive repertoire narrowing was almost entirely limited to many patients aged older than 50 y. By temporarily reducing T cell diversity and by risking expansions of nonbeneficial T cells, these traits may constitute an age-related risk factor for COVID-19, including a vulnerability to new variants for which T cells may provide key protection.

The CD46-Jagged1 interaction is critical for human TH1 immunity.

CD46 is a complement regulator with important roles related to the immune response. CD46 functions as a pathogen receptor and is a potent costimulator for the induction of interferon-γ (IFN-γ)-secreting effector T helper type 1 (T(H)1) cells and their subsequent switch into interleukin 10 (IL-10)-producing regulatory T cells. Here we identified the Notch family member Jagged1 as a physiological ligand for CD46. Furthermore, we found that CD46 regulated the expression of Notch receptors and ligands during T cell activation and that disturbance of the CD46-Notch crosstalk impeded induction of IFN-γ and switching to IL-10. Notably, CD4(+) T cells from CD46-deficient patients and patients with hypomorphic mutations in the gene encoding Jagged1 (Alagille syndrome) failed to mount appropriate T(H)1 responses in vitro and in vivo, which suggested that CD46-Jagged1 crosstalk is responsible for the recurrent infections in subpopulations of these patients.

COVID-19: Using high-throughput flow cytometry to dissect clinical heterogeneity.

Here we consider how high-content flow cytometric methodology at appropriate scale and throughput rapidly provided meaningful biological data in our recent studies of COVID-19, which we discuss in the context of other similar investigations. In our work, high-throughput flow cytometry was instrumental to identify a consensus immune signature in COVID-19 patients, and to investigate the impact of SARS-CoV-2 exposure on patients with either solid or hematological cancers. We provide here some examples of our 'holistic' approach, in which flow cytometry data generated by lymphocyte and myelomonocyte panels were integrated with other analytical metrics, including SARS-CoV-2-specific serum antibody titers, plasma cytokine/chemokine levels, and in-depth clinical annotation. We report how selective differences between T cell subsets were revealed by a newly described flow cytometric TDS assay to distinguish actively cycling T cells in the peripheral blood. By such approaches, our and others' high-content flow cytometry studies collectively identified overt abnormalities and subtle but critical changes that discriminate the immuno-signature of COVID-19 patients from those of healthy donors and patients with non-COVID respiratory infections. Thereby, these studies offered several meaningful biomarkers of COVID-19 severity that have the potential to improve the management of patients and of hospital resources. In sum, flow cytometry provides an important means for rapidly obtaining data that can guide clinical decision-making without requiring highly expensive, sophisticated equipment, and/or "-omics" capabilities. We consider how this approach might be further developed.

Safety and immunogenicity of one versus two doses of the COVID-19 vaccine BNT162b2 for patients with cancer: interim analysis of a prospective observational study.

BACKGROUND: The efficacy and safety profiles of vaccines against SARS-CoV-2 in patients with cancer is unknown. We aimed to assess the safety and immunogenicity of the BNT162b2 (Pfizer-BioNTech) vaccine in patients with cancer. METHODS: For this prospective observational study, we recruited patients with cancer and healthy controls (mostly health-care workers) from three London hospitals between Dec 8, 2020, and Feb 18, 2021. Participants who were vaccinated between Dec 8 and Dec 29, 2020, received two 30 µg doses of BNT162b2 administered intramuscularly 21 days apart; patients vaccinated after this date received only one 30 µg dose with a planned follow-up boost at 12 weeks. Blood samples were taken before vaccination and at 3 weeks and 5 weeks after the first vaccination. Where possible, serial nasopharyngeal real-time RT-PCR (rRT-PCR) swab tests were done every 10 days or in cases of symptomatic COVID-19. The coprimary endpoints were seroconversion to SARS-CoV-2 spike (S) protein in patients with cancer following the first vaccination with the BNT162b2 vaccine and the effect of vaccine boosting after 21 days on seroconversion. All participants with available data were included in the safety and immunogenicity analyses. Ongoing follow-up is underway for further blood sampling after the delayed (12-week) vaccine boost. This study is registered with the NHS Health Research Authority and Health and Care Research Wales (REC ID 20/HRA/2031). FINDINGS: 151 patients with cancer (95 patients with solid cancer and 56 patients with haematological cancer) and 54 healthy controls were enrolled. For this interim data analysis of the safety and immunogenicity of vaccinated patients with cancer, samples and data obtained up to March 19, 2021, were analysed. After exclusion of 17 patients who had been exposed to SARS-CoV-2 (detected by either antibody seroconversion or a positive rRT-PCR COVID-19 swab test) from the immunogenicity analysis, the proportion of positive anti-S IgG titres at approximately 21 days following a single vaccine inoculum across the three cohorts were 32 (94%; 95% CI 81-98) of 34 healthy controls; 21 (38%; 26-51) of 56 patients with solid cancer, and eight (18%; 10-32) of 44 patients with haematological cancer. 16 healthy controls, 25 patients with solid cancer, and six patients with haematological cancer received a second dose on day 21. Of the patients with available blood samples 2 weeks following a 21-day vaccine boost, and excluding 17 participants with evidence of previous natural SARS-CoV-2 exposure, 18 (95%; 95% CI 75-99) of 19 patients with solid cancer, 12 (100%; 76-100) of 12 healthy controls, and three (60%; 23-88) of five patients with haematological cancers were seropositive, compared with ten (30%; 17-47) of 33, 18 (86%; 65-95) of 21, and four (11%; 4-25) of 36, respectively, who did not receive a boost. The vaccine was well tolerated; no toxicities were reported in 75 (54%) of 140 patients with cancer following the first dose of BNT162b2, and in 22 (71%) of 31 patients with cancer following the second dose. Similarly, no toxicities were reported in 15 (38%) of 40 healthy controls after the first dose and in five (31%) of 16 after the second dose. Injection-site pain within 7 days following the first dose was the most commonly reported local reaction (23 [35%] of 65 patients with cancer; 12 [48%] of 25 healthy controls). No vaccine-related deaths were reported. INTERPRETATION: In patients with cancer, one dose of the BNT162b2 vaccine yields poor efficacy. Immunogenicity increased significantly in patients with solid cancer within 2 weeks of a vaccine boost at day 21 after the first dose. These data support prioritisation of patients with cancer for an early (day 21) second dose of the BNT162b2 vaccine. FUNDING: King's College London, Cancer Research UK, Wellcome Trust, Rosetrees Trust, and Francis Crick Institute.

Heteromeric interactions regulate butyrophilin (BTN) and BTN-like molecules governing γδ T cell biology.

The long-held view that gamma delta (γδ) T cells in mice and humans are fundamentally dissimilar, as are γδ cells in blood and peripheral tissues, has been challenged by emerging evidence of the cells' regulation by butyrophilin (BTN) and butyrophilin-like (BTNL) molecules. Thus, murine Btnl1 and the related gene, Skint1, mediate T cell receptor (TCR)-dependent selection of murine intraepithelial γδ T cell repertoires in gut and skin, respectively; BTNL3 and BTNL8 are TCR-dependent regulators of human gut γδ cells; and BTN3A1 is essential for TCR-dependent activation of human peripheral blood Vγ9Vδ2+ T cells. However, some observations concerning BTN/Btnl molecules continue to question the extent of mechanistic conservation. In particular, murine and human gut γδ cell regulation depends on pairings of Btnl1 and Btnl6 and BTNL3 and BTNL8, respectively, whereas blood γδ cells are reported to be regulated by BTN3A1 independent of other BTNs. Addressing this paradox, we show that BTN3A2 regulates the subcellular localization of BTN3A1, including functionally important associations with the endoplasmic reticulum (ER), and is specifically required for optimal BTN3A1-mediated activation of Vγ9Vδ2+ T cells. Evidence that BTNL3/BTNL8 and Btnl1/Btnl6 likewise associate with the ER reinforces the prospect of broadly conserved mechanisms underpinning the selection and activation of γδ cells in mice and humans, and in blood and extralymphoid sites.

flowLearn: fast and precise identification and quality checking of cell populations in flow cytometry.

Motivation: Identification of cell populations in flow cytometry is a critical part of the analysis and lays the groundwork for many applications and research discovery. The current paradigm of manual analysis is time consuming and subjective. A common goal of users is to replace manual analysis with automated methods that replicate their results. Supervised tools provide the best performance in such a use case, however they require fine parameterization to obtain the best results. Hence, there is a strong need for methods that are fast to setup, accurate and interpretable. Results: flowLearn is a semi-supervised approach for the quality-checked identification of cell populations. Using a very small number of manually gated samples, through density alignments it is able to predict gates on other samples with high accuracy and speed. On two state-of-the-art datasets, our tool achieves median(F1)-measures exceeding 0.99 for 31%, and 0.90 for 80% of all analyzed populations. Furthermore, users can directly interpret and adjust automated gates on new sample files to iteratively improve the initial training. Availability and implementation: FlowLearn is available as an R package on https://github.com/mlux86/flowLearn. Evaluation data is publicly available online. Details can be found in the Supplementary Material. Supplementary information: Supplementary data are available at Bioinformatics online.

High throughput automated analysis of big flow cytometry data.

The rapid expansion of flow cytometry applications has outpaced the functionality of traditional manual analysis tools used to interpret flow cytometry data. Scientists are faced with the daunting prospect of manually identifying interesting cell populations in 50-dimensional datasets, equalling the complexity previously only reached in mass cytometry. Data can no longer be analyzed or interpreted fully by manual approaches. While automated gating has been the focus of intense efforts, there are many significant additional steps to the analytical pipeline (e.g., cleaning the raw files, event outlier detection, extracting immunophenotypes). We review the components of a customized automated analysis pipeline that can be generally applied to large scale flow cytometry data. We demonstrate these methodologies on data collected by the International Mouse Phenotyping Consortium (IMPC).

Giant taxon-character matrices: the future of morphological systematics.

Simões et al. () argued that large matrices are linked to the construction of "problematic" characters, and that those characters negatively affect tree topology. In their re-evaluation of two squamate datasets, however, Simões et al. () simply eliminated what they termed "problematic" characters, rather than recode them. This practice ignores potential sources of phylogenetic information and, if it were to be more widely followed, would inhibit the advancement of the field of systematics. Here, we defend the necessity and inevitability of large morphological (phenomic) datasets and discuss best practices for morphological data collection in contemporary phylogenetics.

Gro/TLE enables embryonic stem cell differentiation by repressing pluripotent gene expression.

Gro/TLE proteins (TLE1-4) are a family of transcriptional corepressors acting downstream of multiple signalling pathways. Several TLEs are expressed in a dynamic manner throughout embryonic development and at high levels in embryonic stem cells (ESCs). Here we find that Gro/TLE is not required in ESC for sustaining pluripotency and suppressing differentiation genes, but rather is important for the shutting down of the pluripotency network in differentiation. Consistent with this view, we found that one of the Gro/TLE family, TLE4 is expressed heterogeneously in ESCs in a population that corresponds to a Nanog low subset of ESC culture. TLE4 expression is also increased in response to LIF withdrawal and Fgf/Mek/Erk stimulation. To explore the role of Gro/TLE in more detail we generated an allelic series of knockout ESCs of two TLE genes expressed most dynamically in early differentiation, TLE3 and TLE4. Genetic reduction in TLE dose resulted in an increase in the expression of pluripotency markers and inhibition of ESC differentiation towards both epiblast and endoderm lineages. Overexpression of a drug inducible TLE4 could both rescue TLE3/TLE4 compound phenotypes and induce early expression of endoderm (Hhex-Venus) and neural (Sox1-GFP) reporter genes. Taken together, our results suggest that TLE activity is essential for early differentiation where it acts to suppress the pluripotency network, allowing for the initiation of lineage specific gene expression programs.

Innate TCRß-chain engagement drives human T cells toward distinct memory-like effector phenotypes with immunotherapeutic potentials.

Clonotypic αß T cell responses to cargoes presented by major histocompatibility complex (MHC), MR1, or CD1 proteins underpin adaptive immunity. Those responses are mostly mediated by complementarity-determining region 3 motifs created by quasi-random T cell receptor (TCR) gene rearrangements, with diversity being highest for TCRγδ. Nonetheless, TCRγδ also displays nonclonotypic innate responsiveness following engagement of germline-encoded Vγ-specific residues by butyrophilin (BTN) or BTN-like (BTNL) proteins that uniquely mediate γδ T cell subset selection. We now report that nonclonotypic TCR engagement likewise induces distinct phenotypes in TCRαß+ cells. Specifically, antibodies to germline-encoded human TCRVß motifs consistently activated naïve or memory T cells toward core states distinct from those induced by anti-CD3 or superantigens and from others commonly reported. Those states combined selective proliferation and effector function with activation-induced inhibitory receptors and memory differentiation. Thus, nonclonotypic TCRVß targeting broadens our perspectives on human T cell response modes and might offer ways to induce clinically beneficial phenotypes in defined T cell subsets.

Impact of the B.1.1.7 variant on neutralizing monoclonal antibodies recognizing diverse epitopes on SARS-CoV-2 Spike.

The interaction of the SARS-CoV-2 Spike receptor binding domain (RBD) with the ACE2 receptor on host cells is essential for viral entry. RBD is the dominant target for neutralizing antibodies and several neutralizing epitopes on RBD have been molecularly characterized. Analysis of circulating SARS-CoV-2 variants has revealed mutations arising in the RBD, the N-terminal domain (NTD) and S2 subunits of Spike. To fully understand how these mutations affect the antigenicity of Spike, we have isolated and characterized neutralizing antibodies targeting epitopes beyond the already identified RBD epitopes. Using recombinant Spike as a sorting bait, we isolated >100 Spike-reactive monoclonal antibodies from SARS-CoV-2 infected individuals. ≈45% showed neutralizing activity of which ≈20% were NTD-specific. None of the S2-specific antibodies showed neutralizing activity. Competition ELISA revealed that NTD-specific mAbs formed two distinct groups: the first group was highly potent against infectious virus, whereas the second was less potent and displayed glycan-dependant neutralization activity. Importantly, mutations present in B.1.1.7 Spike frequently conferred resistance to neutralization by the NTD-specific neutralizing antibodies. This work demonstrates that neutralizing antibodies targeting subdominant epitopes need to be considered when investigating antigenic drift in emerging variants.

Acute Immune Signatures and Their Legacies in Severe Acute Respiratory Syndrome Coronavirus-2 Infected Cancer Patients.

Given the immune system's importance for cancer surveillance and treatment, we have investigated how it may be affected by SARS-CoV-2 infection of cancer patients. Across some heterogeneity in tumor type, stage, and treatment, virus-exposed solid cancer patients display a dominant impact of SARS-CoV-2, apparent from the resemblance of their immune signatures to those for COVID-19+ non-cancer patients. This is not the case for hematological malignancies, with virus-exposed patients collectively displaying heterogeneous humoral responses, an exhausted T cell phenotype and a high prevalence of prolonged virus shedding. Furthermore, while recovered solid cancer patients' immunophenotypes resemble those of non-virus-exposed cancer patients, recovered hematological cancer patients display distinct, lingering immunological legacies. Thus, while solid cancer patients, including those with advanced disease, seem no more at risk of SARS-CoV-2-associated immune dysregulation than the general population, hematological cancer patients show complex immunological consequences of SARS-CoV-2 exposure that might usefully inform their care.

Butyrophilin-like proteins display combinatorial diversity in selecting and maintaining signature intraepithelial γδ T cell compartments.

Butyrophilin-like (Btnl) genes are emerging as major epithelial determinants of tissue-associated γδ T cell compartments. Thus, the development of signature, murine TCRγδ+ intraepithelial lymphocytes (IEL) in gut and skin depends on Btnl family members, Btnl1 and Skint1, respectively. In seeking mechanisms underlying these profound effects, we now show that normal gut and skin γδ IEL development additionally requires Btnl6 and Skint2, respectively, and furthermore that different Btnl heteromers can seemingly shape different intestinal γδ+ IEL repertoires. This formal genetic evidence for the importance of Btnl heteromers also applied to the steady-state, since sustained Btnl expression is required to maintain the signature TCR.Vγ7+ IEL phenotype, including specific responsiveness to Btnl proteins. In sum, Btnl proteins are required to select and to maintain the phenotypes of tissue-protective γδ IEL compartments, with combinatorially diverse heteromers having differential impacts on different IEL subsets.

A dynamic COVID-19 immune signature includes associations with poor prognosis.

Improved understanding and management of COVID-19, a potentially life-threatening disease, could greatly reduce the threat posed by its etiologic agent, SARS-CoV-2. Toward this end, we have identified a core peripheral blood immune signature across 63 hospital-treated patients with COVID-19 who were otherwise highly heterogeneous. The signature includes discrete changes in B and myelomonocytic cell composition, profoundly altered T cell phenotypes, selective cytokine/chemokine upregulation and SARS-CoV-2-specific antibodies. Some signature traits identify links with other settings of immunoprotection and immunopathology; others, including basophil and plasmacytoid dendritic cell depletion, correlate strongly with disease severity; while a third set of traits, including a triad of IP-10, interleukin-10 and interleukin-6, anticipate subsequent clinical progression. Hence, contingent upon independent validation in other COVID-19 cohorts, individual traits within this signature may collectively and individually guide treatment options; offer insights into COVID-19 pathogenesis; and aid early, risk-based patient stratification that is particularly beneficial in phasic diseases such as COVID-19.

Author Correction: A dynamic COVID-19 immune signature includes associations with poor prognosis.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Author Correction: A dynamic COVID-19 immune signature includes associations with poor prognosis.

A Correction to this paper has been published: https://doi.org/10.1038/s41591-020-01186-5.