Natural Killer cells demonstrate distinct eQTL and transcriptome-wide disease associations, highlighting their role in autoimmunity

Functional genomics implicates natural killer cells in the pathogenesis of ankylosing spondylitis

Objective

Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased functional genomics approach.

Methods

We integrated genome-wide association study (GWAS) data with epigenomic and transcriptomic datasets of human immune cells. To quantify enrichment of cell type-specific open chromatin or gene expression in AS risk loci, we used three published methods that have successfully identified relevant cell types in other diseases. We performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes.

Results

Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA-seq data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Using the human Space-Time Gut Cell Atlas, we also found significant upregulation of AS-associated genes predominantly in NK cells. Co-localization analysis revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two under-studied loci, ENTR1 (aka SDCCAG3) and B3GNT2.

Conclusion

Our findings suggest that NK cells may play a crucial role in AS development and highlight four putative target genes for functional follow-up in NK cells.

The central role of natural killer cells in mediating acute myocarditis after mRNA COVID-19 vaccination

BACKGROUND

Vaccine-related acute myocarditis is recognized as a rare and specific vaccine complication following mRNA-based COVID-19 vaccinations. The precise mechanisms remain unclear. We hypothesized that natural killer (NK) cells play a central role in its pathogenesis.

METHODS

Samples from 60 adolescents with vaccine-related myocarditis were analyzed, including pro-inflammatory cytokines, cardiac troponin T, genotyping, and immunophenotyping of the corresponding activation subsets of NK cells, monocytes, and T cells. Results were compared with samples from 10 vaccinated individuals without myocarditis and 10 healthy controls.

FINDINGS

Phenotypically, high levels of serum cytokines pivotal for NK cells, including interleukin-1β (IL-1β), interferon α2 (IFN-α2), IL-12, and IFN-γ, were observed in post-vaccination patients with myocarditis, who also had high percentage of CD57+ NK cells in blood, which in turn correlated positively with elevated levels of cardiac troponin T. Abundance of the CD57+ NK subset was particularly prominent in males and in those after the second dose of vaccination. Genotypically, killer cell immunoglobulin-like receptor (KIR) KIR2DL5B(-)/KIR2DS3(+)/KIR2DS5(-)/KIR2DS4del(+) was a risk haplotype, in addition to single-nucleotide polymorphisms related to the NK cell-specific expression quantitative trait loci DNAM-1 and FuT11, which also correlated with cardiac troponin T levels in post-vaccination patients with myocarditis.

CONCLUSION

Collectively, these data suggest that NK cell activation by mRNA COVID-19 vaccine contributed to the pathogenesis of acute myocarditis in genetically and epidemiologically vulnerable subjects.

FUNDING

This work was funded by the Hong Kong Collaborative Research Fund (CRF) 2020/21 and the CRF Coronavirus and Novel Infectious Diseases Research Exercises (reference no. C7149-20G).

A common NFKB1 variant detected through antibody analysis in UK Biobank predicts risk of infection and allergy

Infectious agents contribute significantly to the global burden of diseases through both acute infection and their chronic sequelae. We leveraged the UK Biobank to identify genetic loci that influence humoral immune response to multiple infections. From 45 genome-wide association studies in 9,611 participants from UK Biobank, we identified NFKB1 as a locus associated with quantitative antibody responses to multiple pathogens, including those from the herpes, retro-, and polyoma-virus families. An insertion-deletion variant thought to affect NFKB1 expression (rs28362491), was mapped as the likely causal variant and could play a key role in regulation of the immune response. Using 121 infection- and inflammation-related traits in 487,297 UK Biobank participants, we show that the deletion allele was associated with an increased risk of infection from diverse pathogens but had a protective effect against allergic disease. We propose that altered expression of NFKB1, as a result of the deletion, modulates hematopoietic pathways and likely impacts cell survival, antibody production, and inflammation. Taken together, we show that disruptions to the tightly regulated immune processes may tip the balance between exacerbated immune responses and allergy, or increased risk of infection and impaired resolution of inflammation.

Identifying genetic variants that influence the abundance of cell states in single-cell data

To understand genetic mechanisms driving disease, it is essential but difficult to map how risk alleles affect the composition of cells present in the body. Single-cell profiling quantifies granular information about tissues, but variant-associated cell states may reflect diverse combinations of the profiled cell features that are challenging to predefine. We introduce GeNA (Genotype-Neighborhood Associations), a statistical tool to identify cell state abundance quantitative trait loci (csaQTLs) in high-dimensional single-cell datasets. Instead of testing associations to predefined cell states, GeNA flexibly identifies the cell states whose abundance is most associated with genetic variants. In a genome-wide survey of scRNA-seq peripheral blood profiling from 969 individuals,1GeNA identifies five independent loci associated with shifts in the relative abundance of immune cell states. For example, rs3003-T (p=1.96×10-11) associates with increased abundance of NK cells expressing TNF-α response programs. This csaQTL colocalizes with increased risk for psoriasis, an autoimmune disease that responds to anti-TNF treatments. Flexibly characterizing csaQTLs for granular cell states may help illuminate how genetic background alters cellular composition to confer disease risk.

eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs

The eQTL Catalogue is an open database of uniformly processed human molecular quantitative trait loci (QTLs). We are continuously updating the resource to further increase its utility for interpreting genetic associations with complex traits. Over the past two years, we have increased the number of uniformly processed studies from 21 to 31 and added X chromosome QTLs for 19 compatible studies. We have also implemented Leafcutter to directly identify splice-junction usage QTLs in all RNA sequencing datasets. Finally, to improve the interpretability of transcript-level QTLs, we have developed static QTL coverage plots that visualise the association between the genotype and average RNA sequencing read coverage in the region for all 1.7 million fine mapped associations. To illustrate the utility of these updates to the eQTL Catalogue, we performed colocalisation analysis between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. Although most GWAS loci colocalised both with eQTLs and transcript-level QTLs, we found that visual inspection could sometimes be used to distinguish primary splicing QTLs from those that appear to be secondary consequences of large-effect gene expression QTLs. While these visually confirmed primary splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases.

Systematic visualisation of molecular QTLs reveals variant mechanisms at GWAS loci

Splicing quantitative trait loci (QTLs) have been implicated as a common mechanism underlying complex trait associations. However, utilising splicing QTLs in target discovery and prioritisation has been challenging due to extensive data normalisation which often renders the direction of the genetic effect as well as its magnitude difficult to interpret. This is further complicated by the fact that strong expression QTLs often manifest as weak splicing QTLs and vice versa, making it difficult to uniquely identify the underlying molecular mechanism at each locus. We find that these ambiguities can be mitigated by visualising the association between the genotype and average RNA sequencing read coverage in the region. Here, we generate these QTL coverage plots for 1.7 million molecular QTL associations in the eQTL Catalogue identified with five quantification methods. We illustrate the utility of these QTL coverage plots by performing colocalisation between vitamin D levels in the UK Biobank and all molecular QTLs in the eQTL Catalogue. We find that while visually confirmed splicing QTLs explain just 6/53 of the colocalising signals, they are significantly less pleiotropic than eQTLs and identify a prioritised causal gene in 4/6 cases. All our association summary statistics and QTL coverage plots are freely available at https://www.ebi.ac.uk/eqtl/.

Genetic susceptibility to autoimmunity-Current status and challenges

Autoimmune diseases (ADs) often arise from a combination of genetic and environmental triggers that disrupt the immune system's capability to properly tolerate body self-antigens. Familial studies provided the earliest insights into the risk loci of such diseases, while genome-wide association studies (GWAS) significantly broadened the horizons. A drug targeting a prominent pathological pathway can be applied to multiple indications sharing overlapping mechanisms. Advances in genomic technologies used in genetic studies provide critical insights into future research on gene-environment interactions in autoimmunity. This Review summarizes the history and recent advances in the understanding of genetic susceptibility to ADs and related immune disorders, including coronavirus disease 2019 (COVID-19), and their indications for the development of diagnostic or prognostic markers for translational applications.