Iron dysregulation and inflammatory stress erythropoiesis associates with long-term outcome of COVID-19.

Persistent symptoms following SARS-CoV-2 infection are increasingly reported, although the drivers of post-acute sequelae (PASC) of COVID-19 are unclear. Here we assessed 214 individuals infected with SARS-CoV-2, with varying disease severity, for one year from COVID-19 symptom onset to determine the early correlates of PASC. A multivariate signature detected beyond two weeks of disease, encompassing unresolving inflammation, anemia, low serum iron, altered iron-homeostasis gene expression and emerging stress erythropoiesis; differentiated those who reported PASC months later, irrespective of COVID-19 severity. A whole-blood heme-metabolism signature, enriched in hospitalized patients at month 1-3 post onset, coincided with pronounced iron-deficient reticulocytosis. Lymphopenia and low numbers of dendritic cells persisted in those with PASC, and single-cell analysis reported iron maldistribution, suggesting monocyte iron loading and increased iron demand in proliferating lymphocytes. Thus, defects in iron homeostasis, dysregulated erythropoiesis and immune dysfunction due to COVID-19 possibly contribute to inefficient oxygen transport, inflammatory disequilibrium and persisting symptomatology, and may be therapeutically tractable.

A patient-centric modeling framework captures recovery from SARS-CoV-2 infection.

The biology driving individual patient responses to severe acute respiratory syndrome coronavirus 2 infection remains ill understood. Here, we developed a patient-centric framework leveraging detailed longitudinal phenotyping data and covering a year after disease onset, from 215 infected individuals with differing disease severities. Our analyses revealed distinct 'systemic recovery' profiles, with specific progression and resolution of the inflammatory, immune cell, metabolic and clinical responses. In particular, we found a strong inter-patient and intra-patient temporal covariation of innate immune cell numbers, kynurenine metabolites and lipid metabolites, which highlighted candidate immunologic and metabolic pathways influencing the restoration of homeostasis, the risk of death and that of long COVID. Based on these data, we identified a composite signature predictive of systemic recovery, using a joint model on cellular and molecular parameters measured soon after disease onset. New predictions can be generated using the online tool http://shiny.mrc-bsu.cam.ac.uk/apps/covid-19-systemic-recovery-prediction-app , designed to test our findings prospectively.

Longitudinal analysis reveals that delayed bystander CD8+ T cell activation and early immune pathology distinguish severe COVID-19 from mild disease.

The kinetics of the immune changes in COVID-19 across severity groups have not been rigorously assessed. Using immunophenotyping, RNA sequencing, and serum cytokine analysis, we analyzed serial samples from 207 SARS-CoV2-infected individuals with a range of disease severities over 12 weeks from symptom onset. An early robust bystander CD8+ T cell immune response, without systemic inflammation, characterized asymptomatic or mild disease. Hospitalized individuals had delayed bystander responses and systemic inflammation that was already evident near symptom onset, indicating that immunopathology may be inevitable in some individuals. Viral load did not correlate with this early pathological response but did correlate with subsequent disease severity. Immune recovery is complex, with profound persistent cellular abnormalities in severe disease correlating with altered inflammatory responses, with signatures associated with increased oxidative phosphorylation replacing those driven by cytokines tumor necrosis factor (TNF) and interleukin (IL)-6. These late immunometabolic and immune defects may have clinical implications.

Genetically encoded Runx3 and CD4+ intestinal epithelial lymphocyte deficiencies link SKG mouse and human predisposition to spondyloarthropathy.

Disturbances in immune regulation, intestinal dysbiosis and inflammation characterize ankylosing spondylitis (AS), which is associated with RUNX3 loss-of-function variants. ZAP70W163C mutant (SKG) mice have reduced ZAP70 signaling, spondyloarthritis and ileitis. In small intestine, Foxp3+ regulatory T cells (Treg) and CD4+CD8αα+TCRαß+ intraepithelial lymphocytes (CD4-IEL) control inflammation. TGF-ß and retinoic acid (RA)-producing dendritic cells and MHC-class II+ intestinal epithelial cells (IEC) are required for Treg and CD4-IEL differentiation from CD4+ conventional or Treg precursors, with upregulation of Runx3 and suppression of ThPOK. We show in SKG mouse ileum, that ZAP70W163C or ZAP70 inhibition prevented CD4-IEL but not Treg differentiation, dysregulating Runx3 and ThPOK. TGF-ß/RA-mediated CD4-IEL development, T-cell IFN-γ production, MHC class-II+ IEC, tissue-resident memory T-cell and Runx3-regulated genes were reduced. In AS intestine, CD4-IEL were decreased, while in AS blood CD4+CD8+ T cells were reduced and Treg increased. Thus, genetically-encoded TCR signaling dysfunction links intestinal T-cell immunodeficiency in mouse and human spondyloarthropathy.

Epistatic interactions between killer immunoglobulin-like receptors and human leukocyte antigen ligands are associated with ankylosing spondylitis.

The killer immunoglobulin-like receptors (KIRs), found predominantly on the surface of natural killer (NK) cells and some T-cells, are a collection of highly polymorphic activating and inhibitory receptors with variable specificity for class I human leukocyte antigen (HLA) ligands. Fifteen KIR genes are inherited in haplotypes of diverse gene content across the human population, and the repertoire of independently inherited KIR and HLA alleles is known to alter risk for immune-mediated and infectious disease by shifting the threshold of lymphocyte activation. We have conducted the largest disease-association study of KIR-HLA epistasis to date, enabled by the imputation of KIR gene and HLA allele dosages from genotype data for 12,214 healthy controls and 8,107 individuals with the HLA-B*27-associated immune-mediated arthritis, ankylosing spondylitis (AS). We identified epistatic interactions between KIR genes and their ligands (at both HLA subtype and allele resolution) that increase risk of disease, replicating analyses in a semi-independent cohort of 3,497 cases and 14,844 controls. We further confirmed that the strong AS-association with a pathogenic variant in the endoplasmic reticulum aminopeptidase gene ERAP1, known to alter the HLA-B*27 presented peptidome, is not modified by carriage of the canonical HLA-B receptor KIR3DL1/S1. Overall, our data suggests that AS risk is modified by the complement of KIRs and HLA ligands inherited, beyond the influence of HLA-B*27 alone, which collectively alter the proinflammatory capacity of KIR-expressing lymphocytes to contribute to disease immunopathogenesis.

HLA and KIR Associations of Cervical Neoplasia.

Background: Cervical cancer is the fourth most common cancer in women, and we recently reported human leukocyte antigen (HLA) alleles showing strong associations with cervical neoplasia risk and protection. HLA ligands are recognized by killer immunoglobulin-like receptors (KIRs) expressed on a range of immune cell subsets, governing their proinflammatory activity. We hypothesized that the inheritance of particular HLA-KIR combinations would increase cervical neoplasia risk. Methods: Here, we used HLA and KIR dosages imputed from single-nucleotide polymorphism genotype data from 2143 cervical neoplasia cases and 13858 healthy controls of European decent. Results: The following 4 novel HLA alleles were identified in association with cervical neoplasia, owing to their linkage disequilibrium with known cervical neoplasia-associated HLA-DRB1 alleles: HLA-DRB3*9901 (odds ratio [OR], 1.24; P = 2.49 × 10-9), HLA-DRB5*0101 (OR, 1.29; P = 2.26 × 10-8), HLA-DRB5*9901 (OR, 0.77; P = 1.90 × 10-9), and HLA-DRB3*0301 (OR, 0.63; P = 4.06 × 10-5). We also found that homozygosity of HLA-C1 group alleles is a protective factor for human papillomavirus type 16 (HPV16)-related cervical neoplasia (C1/C1; OR, 0.79; P = .005). This protective association was restricted to carriers of either KIR2DL2 (OR, 0.67; P = .00045) or KIR2DS2 (OR, 0.69; P = .0006). Conclusions: Our findings suggest that HLA-C1 group alleles play a role in protecting against HPV16-related cervical neoplasia, mainly through a KIR-mediated mechanism.

Contribution of HLA and KIR Alleles to Systemic Sclerosis Susceptibility and Immunological and Clinical Disease Subtypes.

Systemic sclerosis (SSc) is an autoinflammatory, fibrotic condition of unknown aetiology. The presence of detectable autoantibodies against diverse nuclear antigens, as well as strong HLA associations with disease, suggest autoimmune involvement, however the links between endogenous and exogenous risk factors and SSc pathology remain undetermined. We have conducted a genetic analysis of HLA inheritance in two independent and meta-analysed cohorts of 1,465 SSc cases and 13,273 controls, including stratified association analyses in clinical and autoantibody positive subgroups of disease. Additionally, we have used patient genotypes to impute gene dosages across the KIR locus, encoding paired activating and inhibitory lymphocyte receptors for Class I HLA ligands, to conduct the largest analysis of KIR-HLA epistatic interactions in SSc to date. We confirm previous Class II HLA associations with SSc risk and report a new Class I association with haplotype HLA-B*44:03-HLA-C*16:01 at genome-wide significance (GWS). We further report statistically significant HLA associations with clinical and serological subtypes of disease through direct case-case comparison, and report a new association of HLA-DRB1*15:01, previously shown to bind topoisomerase-1 derived peptides, with anti-topoisomerase (ATA) positive disease. Finally, we identify genetic epistasis between KIRs and HLA class I ligands, suggesting genetic modulation of lymphocyte activation may further contribute to an individual's underlying disease risk. Taken together, these findings support future functional investigation into endogenous immunological and environmental stimuli for disrupted immune tolerance in SSc.

The impact of hypoxia on B cells in COVID-19.

BACKGROUND: Prominent early features of COVID-19 include severe, often clinically silent, hypoxia and a pronounced reduction in B cells, the latter important in defence against SARS-CoV-2. This presentation resembles the phenotype of mice with VHL-deficient B cells, in which Hypoxia-Inducible Factors are constitutively active, suggesting hypoxia might drive B cell abnormalities in COVID-19. METHODS: Detailed B cell phenotyping was undertaken by flow-cytometry on longitudinal samples from patients with COVID-19 across a range of severities (NIHR Cambridge BioResource). The impact of hypoxia on the transcriptome was assessed by single-cell and whole blood RNA sequencing analysis. The direct effect of hypoxia on B cells was determined through immunisation studies in genetically modified and hypoxia-exposed mice. FINDINGS: We demonstrate the breadth of early and persistent defects in B cell subsets in moderate/severe COVID-19, including reduced marginal zone-like, memory and transitional B cells, changes also observed in B cell VHL-deficient mice. These findings were associated with hypoxia-related transcriptional changes in COVID-19 patient B cells, and similar B cell abnormalities were seen in mice kept in hypoxic conditions. INTERPRETATION: Hypoxia may contribute to the pronounced and persistent B cell pathology observed in acute COVID-19 pneumonia. Assessment of the impact of early oxygen therapy on these immune defects should be considered, as their correction could contribute to improved outcomes. FUNDING: Evelyn Trust, Addenbrooke's Charitable Trust, UKRI/NIHR, Wellcome Trust.

Altered Repertoire Diversity and Disease-Associated Clonal Expansions Revealed by T Cell Receptor Immunosequencing in Ankylosing Spondylitis Patients.

OBJECTIVE: Ankylosing spondylitis (AS) is a common spondyloarthropathy primarily affecting the axial skeleton and strongly associated with HLA-B*27 carriage. Genetic evidence implicates both autoinflammatory processes and autoimmunity against an HLA-B*27-restricted autoantigen in immunopathology. In addition to articular symptoms, up to 70% of AS patients present with concurrent bowel inflammation, suggesting that adverse interactions between a genetically primed host immune system and the gut microbiome contribute to the disease. Accordingly, this study aimed to characterize adaptive immune responses to antigenic stimuli in AS. METHODS: The peripheral CD4 and CD8 T cell receptor (TCR) repertoire was profiled in AS patients (n = 47) and HLA-B*27-matched healthy controls (n = 38). Repertoire diversity was estimated using the Normalized Shannon Diversity Entropy (NSDE) index, and univariate and multivariate statistical analyses were performed to characterize AS-associated clonal signatures. Furthermore, T cell proliferation and cytokine production in response to immunogenic antigen exposure were investigated in vitro in peripheral blood mononuclear cells from AS patients (n = 19) and HLA-B*27-matched healthy controls (n = 14). RESULTS: Based on the NSDE measure of sample diversity across CD4 and CD8 T cell repertoires, AS patients showed increased TCR diversity compared to healthy controls (for CD4 T cells, P = 7.8 × 10-6 ; for CD8 T cells, P = 9.3 × 10-4 ), which was attributed to a significant reduction in the magnitude of peripheral T cell expansions globally. Upon in vitro stimulation, fewer T cells from AS patients than from healthy controls expressed interferon-γ (for CD8 T cells, P = 0.03) and tumor necrosis factor (for CD4 T cells, P = 0.01; for CD8 T cells, P = 0.002). In addition, the CD8 TCR signature was altered in HLA-B*27+ AS patients compared to healthy controls, with significantly expanded Epstein-Barr virus-specific clonotypes (P = 0.03) and cytomegalovirus-specific clonotypes (P = 0.02). HLA-B*27+ AS patients also showed an increased incidence of "public" CD8 TCRs, representing identical clonotypes emerging in response to common antigen encounters, including homologous clonotypes matching those previously isolated from individuals with bacterial-induced reactive arthritis. CONCLUSION: The dynamics of peripheral T cell responses in AS patients are altered, suggesting that differential antigen exposure and disrupted adaptive immunity are underlying features of the disease.

The genetics, structure and function of the M1 aminopeptidase oxytocinase subfamily and their therapeutic potential in immune-mediated disease.

The oxytocinase subfamily of M1 aminopeptidases plays an important role in processing and trimming of peptides for presentation on major histocompatibility (MHC) Class I molecules. Several large-scale genomic studies have identified association of members of this family of enzymes, most notably ERAP1 and ERAP2, with immune-mediated diseases including ankylosing spondylitis, psoriasis and birdshot chorioretinopathy. Much is now known about the genetics of these enzymes and how genetic variants alter their function, but how these variants contribute to disease remains largely unresolved. Here we discuss what is known about their structure and function and highlight some of the knowledge gaps that affect development of drugs targeting these enzymes.

Genetic Variants in ERAP1 and ERAP2 Associated With Immune-Mediated Diseases Influence Protein Expression and the Isoform Profile.

OBJECTIVE: Endoplasmic reticulum aminopeptidase 1 (ERAP-1) and ERAP-2, encoded on chromosome 5q15, trim endogenous peptides for HLA-mediated presentation to the immune system. Polymorphisms in ERAP1 and/or ERAP2 are strongly associated with several immune-mediated diseases with specific HLA backgrounds, implicating altered peptide handling and presentation as prerequisites for autoreactivity against an arthritogenic peptide. Given the thorough characterization of disease risk-associated polymorphisms that alter ERAP activity, this study aimed instead to interrogate the expression effect of chromosome 5q15 polymorphisms to determine their effect on ERAP isoform and protein expression. METHODS: RNA sequencing and genotyping across chromosome 5q15 were performed to detect genetic variants in ERAP1 and ERAP2 associated with altered total gene and isoform-specific expression. The functional implication of a putative messenger RNA splice-altering variant on ERAP-1 protein levels was validated using mass spectrometry. RESULTS: Polymorphisms associated with ankylosing spondylitis (AS) significantly influenced the transcript and protein expression of ERAP-1 and ERAP-2. Disease risk-associated polymorphisms in and around both genes were also associated with increased gene expression. Furthermore, key risk-associated ERAP1 variants were associated with altered transcript splicing, leading to allele-dependent alternate expression of 2 distinct isoforms and significant differences in the type of ERAP-1 protein produced. CONCLUSION: In accordance with studies demonstrating that polymorphisms that increase aminopeptidase activity predispose to immune disease, the increased risk also attributed to increased expression of ERAP1 and ERAP2 supports the notion of using aminopeptidase inhibition to treat AS and other ERAP-associated conditions.