Immunogenetics: Untangling the immune basis of disease susceptibility
Why do people respond differently to the same infectious agent? A possible answer is that differences in the immune responses of individuals, caused by genetic variability between them, are responsible. Indeed, many studies have looked for associations between genes involved in immunity and disease outcome (Buniello et al., 2019), and it has been found that a gene complex called the human leukocyte antigen (HLA) system has a central role (Matzaraki et al., 2017).
The HLA genes encode the major histocompatibility complex (MHC). MHC proteins are found on the surface of cells and each one presents peptides from proteins within the cell (either proteins native to the cell or from foreign entities like bacteria) to immune cells. Thus, MHC proteins allow immune cells to recognize if a given cell is dangerous or not. For example, immune cells called CD8+ T lymphocytes recognize cells that have been infected with viruses because receptors on these T cells bind specific MHC proteins loaded with viral peptides on the surface of the infected cells. Likewise, other types of immune cells – such as natural killer cells, macrophages and dendritic cells – have receptors that bind to other parts of the MHC protein displayed by infected cells (Augusto and Petzl-Erler, 2015; Hudson and Allen, 2016).
The HLA is extremely diverse among individuals. It has been assumed that the association between some diseases and certain HLA alleles is evidence for a central role of CD8+ T cells in that disease. However, given that the HLA can also interact with receptors on other immune cells, disentangling the contributions of the different arms of the immune system remains a challenge. Now, in eLife, Becca Asquith of Imperial College and colleagues – including Bisrat Debebe as first author, researchers from various institutes in the UK and US, and the IAVI Protocol C Investigators – report the results of a new approach to disentangling these contributions (Debebe et al., 2020).
Debebe et al. reasoned that it may be possible to predict which immune cells are responsible for fighting a specific disease since receptors from particular immune cells interact with different regions of the MHC protein. For example, susceptibility to a disease may be associated with individuals carrying MHC proteins that are similar in the region that presents protein fragments to the T-cell receptor (which can be thought of as being similar in ‘T-cell receptor space’). If this is the case, then CD8+ T cells are likely involved in fighting the disease. The same reasoning can be made for associations with the other receptors, such as the killer immunoglobulin-like receptor (KIR) receptors in natural killer cells and leukocyte immunoglobulin-like receptor (LILR) in myeloid cells (Figure 1).
The challenge in this method is defining MHC similarity. Debebe et al. did so by taking an HLA known to associate with a specific disease outcome and predicting its ability to bind to different immune cell receptors. Next, they compared this predicted binding ability to that of all other HLAs, and found the HLAs with the highest binding similarity. Finally, they checked if these HLAs are also associated with disease susceptibility.
Debebe et al. applied this method to three infections (human T-cell leukemia virus type one or HTLV-1; hepatitis C virus; and HIV-1) and one autoimmune disease (Crohn’s disease). They used well-characterized cohorts where both the HLA and the outcome of the disease were known for each individual.
They studied three known HLA-outcome associations for HTLV-1 and found that other HLAs that were similar in T-cell receptor space were also associated with HTLV-1. However, other HLAs that were similar to these three HLAs in LILR space or KIR space (see Figure 1) were not associated with HTLV-1. Thus, they concluded that HTLV-1 is most likely controlled by CD8+ T cells. They also studied three HLA-outcome associations for hepatitis C, reaching the same conclusion: disease outcome is mostly controlled by CD8+ T cells. For HIV-1, Debebe et al. found that natural killer cells and CD8+ T cells both have a role, which is consistent with previous results (Bashirova et al., 2011; McBrien et al., 2018). They also studied two known HLA-outcome associations for Crohn’s disease in a large sample (2650 individuals), but other HLAs that were similar in T-cell, LILR or KIR space did not have significant effect on Crohn's disease outcome.
The approach suggested by Debebe et al. generates insights into which immune cells may be involved in fighting a disease. There are two issues that need to be explored further. First, while we know a great deal about T-cell receptors, we know less about KIR, and even less about LILR binding (Hirayasu and Arase, 2015). This knowledge (or lack thereof) is part of the definition of similarity used by Debebe et al., which can and should be updated as more data becomes available (Gwozdowicz et al., 2019). The second issue pertains to the fact that it may be possible to define similarity in other ways, leading to different results that provide further insights into how genetic variability affects the immune response.
References
-
HLA/KIR restraint of HIV: surviving the fittestAnnual Review of Immunology 29:295–317.https://doi.org/10.1146/annurev-immunol-031210-101332
-
KIR specificity and avidity of standard and unusual C1, C2, Bw4, Bw6 and A3/11 amino acid motifs at entire HLA:KIR interface between NK and target cells, the functional and evolutionary classification of HLA class I moleculesInternational Journal of Immunogenetics 46:217–231.https://doi.org/10.1111/iji.12433
-
Leukocyte Ig-like receptors - A model for MHC class I disease associationsFrontiers in Immunology 7:281.https://doi.org/10.3389/fimmu.2016.00281
-
Mechanisms of CD8+ T cell-mediated suppression of HIV/SIV replicationEuropean Journal of Immunology 48:898–914.https://doi.org/10.1002/eji.201747172
Article and author information
Author details
Publication history
Copyright
© 2020, Ribeiro and Graca
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,828
- views
-
- 292
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Genetics and Genomics
Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed the effects of covariate stratification and interaction on body mass index (BMI) PGS (PGSBMI) across four cohorts of European (N = 491,111) and African (N = 21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R2 being nearly double between best- and worst-performing quintiles for certain covariates. Twenty-eight covariates had significant PGSBMI–covariate interaction effects, modifying PGSBMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R2 differences among strata and interaction effects – across all covariates, their main effects on BMI were correlated with their maximum R2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGSBMI individuals have highest R2 and increase in PGS effect. Using quantile regression, we show the effect of PGSBMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGSBMI performance and effects, we investigated ways to increase model performance taking into account nonlinear effects. Machine learning models (neural networks) increased relative model R2 (mean 23%) across datasets. Finally, creating PGSBMI directly from GxAge genome-wide association studies effects increased relative R2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGSBMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.
-
- Computational and Systems Biology
- Neuroscience
Audiovisual information reaches the brain via both sustained and transient input channels, representing signals’ intensity over time or changes thereof, respectively. To date, it is unclear to what extent transient and sustained input channels contribute to the combined percept obtained through multisensory integration. Based on the results of two novel psychophysical experiments, here we demonstrate the importance of the transient (instead of the sustained) channel for the integration of audiovisual signals. To account for the present results, we developed a biologically inspired, general-purpose model for multisensory integration, the multisensory correlation detectors, which combines correlated input from unimodal transient channels. Besides accounting for the results of our psychophysical experiments, this model could quantitatively replicate several recent findings in multisensory research, as tested against a large collection of published datasets. In particular, the model could simultaneously account for the perceived timing of audiovisual events, multisensory facilitation in detection tasks, causality judgments, and optimal integration. This study demonstrates that several phenomena in multisensory research that were previously considered unrelated, all stem from the integration of correlated input from unimodal transient channels.