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

Revealing the role of natural killer cells in ankylosing spondylitis: identifying diagnostic biomarkers and therapeutic targets

BACKGROUND

Ankylosing spondylitis (AS) is a chronic autoimmune disease that primarily affects the axial joints. Immune cells play a key role in the pathogenesis of AS. This study integrated bioinformatics methods with experimental validation to explore the role of natural killer (NK) cells in AS.

METHODS

Two microarray datasets, GSE25101 and GSE73754, were selected, and the scRNA-seq data were obtained from GSE194315 and Liu's research. Differentially expressed genes (DEGs) and functional enrichment analysis were performed respectively. Weighted gene co-expression network analysis (WGCNA) was conducted to identify key modules of co-expressed genes and genes involved in NK cell function. The diagnostic value of the identified key genes was evaluated using ROC curves, logistic regression analysis, and a nomogram. Real-time PCR (RT-PCR) was used to quantified the expression of genes. Statistical analysis was conducted using the R software package, and a p-value of less than 0.05 was considered statistically significant.

RESULTS

Pathways enrichment analysis revealed the involvement of NK cell-mediated immune pathways and regulation of the innate immune response, indicating the crucial role of innate immunity, especially NK cells, in AS pathogenesis. The construction of a co-expression network revealed that the MElightyellow module was most relevant to the NK cell-mediated immune pathway. IL2RB, CD247, PLEKHF1, EOMES, S1PR5, FGFBP2 from the MElightyellow module were identified as key genes involved in NK cell-mediated immune response and served as potential diagnostic biomarkers for AS, with moderate to high diagnostic values based on AUC values. Further analysis using scRNA-seq profiling revealed the higher expression level of IL2RB, CD247, PLEKHF1, S1PR5, FGFBP2 in NK cells compared to that in other cell types. CD247, PLEKHF1, EOMES, S1PR5, and FGFBP2 were reduced expressed in AS patients as compare to control group verified by scRNA-seq data, CD247, EOMES, FGFBP2, IL2RB and S1PR5 were reduced expressed verified by RT-PCR, and PLEKHF1, S1PR5, and FGFBP2 was upregulated after TNF-α blocker therapy.

CONCLUSION

The study revealed the potential role of NK cells and identified IL2RB, CD247, PLEKHF1, EOMES, S1PR5, and FGFBP2 as key genes associated with NK cells in the pathogenesis of AS.

Disease associations of natural killer (NK) cell KIR gene content variation in 352,783 Finns

Allelic, gene presence/absence, and gene-copy number variations in the KIR genes encoding Natural Killer (NK) cell surface receptors have been reported to be associated in case-control studies with infectious and autoimmune diseases, and relapse after stem cell transplantation. To understand more comprehensively the role of KIR gene presence/absence variation and HLA-KIR interactions in disease susceptibility, we imputed from genome SNP data the presence and absence of 10 KIR genes in the FinnGen cohort. The cohort consists of 352,783 Finns with extensive phenotypes from the national health registries. We tested associations between 762 FinnGen phenotypes and presence/absence variation based on imputation of KIR genes using 5,900 SNPs located in the KIR genomic segment. Our results provide a platform to query HLA-KIR associations in a large population cohort. We found 13 phenotype - KIR gene or KIR - HLA C combination associations with false discovery rate < 0.05. These results differ from the very high number of associations between HLA alleles and diseases reported earlier in the FinnGen cohort. Five of the 13 significant associations included malignant phenotypes, e.g., melanoma, thyroid gland neoplasm, and haematopoietic malignancy, supporting the essential role of NK cells in controlling malignancy.

Identification of necroptosis-related genes in ankylosing spondylitis by bioinformatics and experimental validation

The pathogenesis of ankylosing spondylitis (AS) remains unclear, and while recent studies have implicated necroptosis in various autoimmune diseases, an investigation of its relationship with AS has not been reported. In this study, we utilized the Gene Expression Omnibus database to compare gene expressions between AS patients and healthy controls, identifying 18 differentially expressed necroptosis-related genes (DENRGs), with 8 upregulated and 10 downregulated. Through the application of three machine learning algorithms-least absolute shrinkage and selection operation, support vector machine-recursive feature elimination and random forest-two hub genes, FASLG and TARDBP, were pinpointed. These genes demonstrated high specificity and sensitivity for AS diagnosis, as evidenced by receiver operating characteristic curve analysis. These findings were further supported by external datasets and cellular experiments, which confirmed the downregulation of FASLG and upregulation of TARDBP in AS patients. Immune cell infiltration analysis suggested that CD4+ T cells, CD8+ T cells, NK cells and neutrophils may be associated with the development of AS. Notably, in the group with high FASLG expression, there was a significant infiltration of CD8+ T cells, memory-activated CD4+ T cells and resting NK cells, with relatively less infiltration of memory-resting CD4+ T cells and neutrophils. Conversely, in the group with high TARDBP expression, there was enhanced infiltration of naïve CD4+ T cells and M0 macrophages, with a reduced presence of memory-resting CD4+ T cells. In summary, FASLG and TARDBP may contribute to AS pathogenesis by regulating the immune microenvironment and immune-related signalling pathways. These findings offer new insights into the molecular mechanisms of AS and suggest potential new targets for therapeutic strategies.

Unveiling the Significance of HLA and KIR Diversity in Underrepresented Populations

Human leukocyte antigen (HLA) molecules and their relationships with natural killer (NK) cells, specifically through their interaction with killer-cell immunoglobulin-like receptors (KIRs), exhibit robust associations with the outcomes of diverse diseases. Moreover, genetic variations in HLA and KIR immune system genes offer limitless depths of complexity. In recent years, a surge of high-powered genome-wide association studies (GWASs) utilizing single nucleotide polymorphism (SNP) arrays has occurred, significantly advancing our understanding of disease pathogenesis. Additionally, advances in HLA reference panels have enabled higher resolution and more reliable imputation, allowing for finer-grained evaluation of the association between sequence variations and disease risk. However, it is essential to note that the majority of these GWASs have focused primarily on populations of Caucasian and Asian origins, neglecting underrepresented populations in Latin America and Africa. This omission not only leads to disparities in health care access but also restricts our knowledge of novel genetic variants involved in disease pathogenesis within these overlooked populations. Since the KIR and HLA haplotypes prevalent in each population are clearly modelled by the specific environment, the aim of this review is to encourage studies investigating HLA/KIR involvement in infection and autoimmune diseases, reproduction, and transplantation in underrepresented populations.

Biology of HLA class I associated inflammatory diseases

Human leukocyte antigen (HLA) class I association is a well-established feature of common and uncommon inflammatory diseases, but it is unknown whether it impacts the pathogenesis of these disorders. The "arthritogenic peptide" hypothesis proposed initially for HLA-B27-associated ankylosing spondylitis (AS) seems the most intuitive to serve as a model for other HLA class I-associated diseases, but evidence supporting it has been scarce. Recent technological advances and the discovery of epistatic relationships between disease-associated HLA class I and endoplasmic reticulum aminopeptidase (ERAP) coding variants have led to the generation of new data and conceptual approaches to the problem requiring its re-examination. Continued success in these endeavors holds promise to resolve a Gordian Knot in human immunobiology. It may ultimately benefit patients by enabling the development of new therapies and precision tools for assessing disease risk and predicting treatment responses.

EULAR study group on ‘MHC-I-opathy’: identifying disease-overarching mechanisms across disciplines and borders

The ‘MHC-I (major histocompatibility complex class I)-opathy’ concept describes a family of inflammatory conditions with overlapping clinical manifestations and a strong genetic link to the MHC-I antigen presentation pathway. Classical MHC-I-opathies such as spondyloarthritis, Behçet’s disease, psoriasis and birdshot uveitis are widely recognised for their strong association with certain MHC-I alleles and gene variants of the antigen processing aminopeptidases ERAP1 and ERAP2 that implicates altered MHC-I peptide presentation to CD8+T cells in the pathogenesis. Progress in understanding the cause and treatment of these disorders is hampered by patient phenotypic heterogeneity and lack of systematic investigation of the MHC-I pathway.

Here, we discuss new insights into the biology of MHC-I-opathies that strongly advocate for disease-overarching and integrated molecular and clinical investigation to decipher underlying disease mechanisms. Because this requires transformative multidisciplinary collaboration, we introduce the EULAR study group on MHC-I-opathies to unite clinical expertise in rheumatology, dermatology and ophthalmology, with fundamental and translational researchers from multiple disciplines such as immunology, genomics and proteomics, alongside patient partners. We prioritise standardisation of disease phenotypes and scientific nomenclature and propose interdisciplinary genetic and translational studies to exploit emerging therapeutic strategies to understand MHC-I-mediated disease mechanisms. These collaborative efforts are required to address outstanding questions in the etiopathogenesis of MHC-I-opathies towards improving patient treatment and prognostication.

The genetic backbone of ankylosing spondylitis: how knowledge of genetic susceptibility informs our understanding and management of disease

Ankylosing spondylitis (AS) is a seronegative, chronic inflammatory arthritis with high genetic burden. A strong association with HLA-B27 has long been established, but to date its contribution to disease aetiology remains unresolved. Recent insights through genome wide studies reveal an increasing array of immunogenetic risk variants extraneous to the HLA complex in AS cohorts. These genetic traits build a complex profile of disease causality, highlighting several molecular pathways associated with the condition. This and other evidence strongly implicates T-cell-driven pathology, revolving around the T helper 17 cell subset as an important contributor to disease. This prominence of the T helper 17 cell subset has presented the opportunity for therapeutic intervention through inhibition of interleukins 17 and 23 which drive T helper 17 activity. While targeting of interleukin 17 has proven effective, this success has not been replicated with interleukin 23 inhibition in AS patients. Evidence points to significant genetic diversity between AS patients which may, in part, explain the observed refractoriness among a proportion of patients. In this review we discuss the impact of genetics on our understanding of AS and its relationship with closely linked pathologies. We further explore how genetics can be used in the development of therapeutics and as a tool to assist in the diagnosis and management of patients. This evidence indicates that genetic profiling should play a role in the clinician’s choice of therapy as part of a precision medicine strategy towards disease management.

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.

Advances of Genomic Medicine in Psoriatic Arthritis

Psoriatic arthritis (PsA) is a common type of inflammatory arthritis found in up to 40% of patients with psoriasis. Although early diagnosis is important for reducing the risk of irreversible structural damage, there are no adequate screening tools for this purpose, and there are no clear markers of predisposition to the disease. Much evidence indicates that PsA disorder is complex and heterogeneous, where genetic and environmental factors converge to trigger inflammatory events and the development of the disease. Nevertheless, the etiologic events that underlie PsA are complex and not completely understood. In this review, we describe the existing data in PsA in order to highlight the need for further research in this disease to progress in the knowledge of its pathobiology and to obtain early diagnosis tools for these patients.

OUP accepted manuscript

Genetic variants in the human leukocyte antigen and killer cell immunoglobulin-like receptor regions have been associated with many brain-related diseases, but how they shape brain structure and function remains unclear. To identify the genetic variants in HLA and KIR genes associated with human brain phenotypes, we performed a genetic association study of ∼30 000 European unrelated individuals using brain MRI phenotypes generated by the UK Biobank (UKB). We identified 15 HLA alleles in HLA class I and class II genes significantly associated with at least one brain MRI-based phenotypes (P < 5 × 10⁻⁸). These associations converged on several main haplotypes within the HLA. In particular, the human leukocyte antigen alleles within an ancestral haplotype 8.1 were associated with multiple MRI measures, including grey matter volume, cortical thickness (TH) and diffusion MRI (dMRI) metrics. These alleles have been strongly associated with schizophrenia. Additionally, associations were identified between HLA-DRB1*04∼DQA1*03:01∼DQB1*03:02 and isotropic volume fraction of diffusion MRI in multiple white matter tracts. This haplotype has been reported to be associated with Parkinson’s disease. These findings suggest shared genetic associations between brain MRI biomarkers and brain-related diseases. Additionally, we identified 169 associations between the complement component 4 (C4) gene and imaging phenotypes. We found that C4 gene copy number was associated with cortical TH and dMRI metrics. No KIR gene copy numbers were associated with image-derived phenotypes at genome-wide threshold. To address the multiple testing burden in the phenome-wide association study, we performed a multi-trait association analysis using trait-based association test that uses extended Simes procedure and identified MRI image-specific associations. This study contributes to insight into how critical immune genes affect brain-related traits as well as the development of neurological and neuropsychiatric disorders.

Mechanistic diversity in MHC class I antigen recognition

Throughout its evolution, the human immune system has developed a plethora of strategies to diversify the antigenic peptide sequences that can be targeted by the CD8+ T cell response against pathogens and aberrations of self. Here we provide a general overview of the mechanisms that lead to the diversity of antigens presented by MHC class I complexes and their recognition by CD8+ T cells, together with a more detailed analysis of recent progress in two important areas that are highly controversial: the prevalence and immunological relevance of unconventional antigen peptides; and cross-recognition of antigenic peptides by the T cell receptors of CD8+ T cells.

Cellular and molecular diversity in spondyloarthritis

The spondyloarthritides are a cluster of inflammatory rheumatic diseases characterized by different diagnostic entities with heterogeneous phenotypes. The current classification system groups spondyloarthritis patients in two main categories, axial and peripheral spondyloarthritis, providing a framework wherein the clinical picture guides the treatment. However, the heterogeneity of the clinical manifestations of the pathologies, even when residing in the same group, highlights the importance of analyzing the smallest features of each entity to understand how different cellular subsets evolve, what the underlying mechanisms are and what biological markers can be identified and validated to evaluate the stage of disease and the corresponding efficacy of treatments. In this review, we will focus mostly on axial spondyloarthritis, report current knowledge concerning the cellular populations involved in its pathophysiology, and their molecular diversity. We will discuss the implications of such a diversity, and their meaning in terms of patients' stratification.

Interaction of HLA Class II rs9272219 and TMPO rs17028450 (Arg690Cys) Variants Affects Neuromyelitis Optica Spectrum Disorder Susceptibility in an Admixed Mexican Population

Neuromyelitis Optica Spectrum Disorder (NMOSD) is a demyelinating autoimmune disease of the central nervous system, more prevalent in individuals of non-European ancestry. Few studies have analyzed genetic risk factors in NMOSD, and HLA class II gene variation has been associated NMOSD risk in various populations including Mexicans. Thymopoietin (TMPO) has not been tested as a candidate gene for NMOSD or other autoimmune disease, however, experimental evidence suggests this gene may be involved in negative selection of autoreactive T cells and autoimmunity. We thus investigated whether the missense TMPO variant rs17028450 (Arg630Cys, frequent in Latin America) is associated with NMOSD, and whether this variant shows an interaction with HLA-class II rs9272219, previously associated with NMOSD risk. A total of 119 Mexican NMOSD patients, 1208 controls and 357 Native Mexican individuals were included. The HLA rs9272219 “T” risk allele frequency ranged from 21 to 68%, while the rs17028450 “T” minor allele frequency was as high as 18% in Native Mexican groups. Both rs9272219 and rs17028450 were significantly associated with NMOSD risk under additive models (OR = 2.48; p = 8 × 10^–10 and OR = 1.59; p = 0.0075, respectively), and a significant interaction between both variants was identified with logistic regression models (p = 0.048). Individuals bearing both risk alleles had an estimated 3.9-fold increased risk of NMOSD. To our knowledge, this is the first study reporting an association of TMPO gene variation with an autoimmune disorder and the interaction of specific susceptibility gene variants, that may contribute to the genetic architecture of NMOSD in admixed Latin American populations.

Nanomaterials for the Diagnosis and Treatment of Inflammatory Arthritis

Nanomaterials have received increasing attention due to their unique chemical and physical properties for the treatment of rheumatoid arthritis (RA), the most common complex multifactorial joint-associated autoimmune inflammatory disorder. RA is characterized by an inflammation of the synovium with increased production of proinflammatory cytokines (IL-1, IL-6, IL-8, and IL-10) and by the destruction of the articular cartilage and bone, and it is associated with the development of cardiovascular disorders such as heart attack and stroke. While a number of imaging tools allow for the monitoring and diagnosis of inflammatory arthritis, and despite ongoing work to enhance their sensitivity and precision, the proper assessment of RA remains difficult particularly in the early stages of the disease. Our goal here is to describe the benefits of applying various nanomaterials as next-generation RA imaging and detection tools using contrast agents and nanosensors and as improved drug delivery systems for the effective treatment of the disease.

HLA-A29 and Birdshot Uveitis: Further Down the Rabbit Hole

HLA class I alleles constitute established risk factors for non-infectious uveitis and preemptive genotyping of HLA class I alleles is standard practice in the diagnostic work-up. The HLA-A29 serotype is indispensable to Birdshot Uveitis (BU) and renders this enigmatic eye condition a unique model to better understand how the antigen processing and presentation machinery contributes to non-infectious uveitis or chronic inflammatory conditions in general. This review will discuss salient points regarding the protein structure of HLA-A29 and how key amino acid positions impact the peptide binding preference and interaction with T cells. We discuss to what extent the risk genes ERAP1 and ERAP2 uniquely affect HLA-A29 and how the discovery of a HLA-A29-specific submotif may impact autoantigen discovery. We further provide a compelling argument to solve the long-standing question why BU only affects HLA-A29-positive individuals from Western-European ancestry by exploiting data from the 1000 Genomes Project. We combine novel insights from structural and immunopeptidomic studies and discuss the functional implications of genetic associations across the HLA class I antigen presentation pathway to refine the etiological basis of Birdshot Uveitis.