Functional epistasis on a common MHC haplotype associated with multiple sclerosis

HLA-transgenic mouse models to study autoimmune central nervous system diseases

It is known that certain human leukocyte antigen (HLA) genes are associated with autoimmune central nervous system (CNS) diseases, such as multiple sclerosis (MS), but their exact role in disease susceptibility and etiopathogenesis remains unclear. The best studied HLA-associated autoimmune CNS disease is MS, and thus will be the primary focus of this review. Other HLA-associated autoimmune CNS diseases, such as autoimmune encephalitis and neuromyelitis optica will be discussed. The lack of animal models to accurately capture the complex human autoimmune response remains a major challenge. HLA transgenic (tg) mice provide researchers with powerful tools to investigate the underlying mechanisms promoting susceptibility and progression of HLA-associated autoimmune CNS diseases, as well as for elucidating the myelin epitopes potentially targeted by T cells in autoimmune disease patients. We will discuss the potential role(s) of autoimmune disease-associated HLA alleles in autoimmune CNS diseases and highlight information provided by studies using HLA tg mice to investigate the underlying pathological mechanisms and opportunities to use these models for development of novel therapies.

Validation and functional follow-up of cervical cancer risk variants at the HLA locus

Cervical cancer is the fourth most common cancer in females. Genome-wide association studies (GWASs) have proposed cervical cancer susceptibility variants at the HLA locus on chromosome 6p21. To corroborate these findings and investigate their functional impact in cervical tissues and cell lines, we genotyped nine variants from cervical cancer GWASs (rs17190106, rs535777, rs1056429, rs2763979, rs143954678, rs113937848, rs3117027, rs3130214, and rs9477610) in a German hospital-based series of 1122 invasive cervical cancers, 1408 dysplasias, and 1196 healthy controls. rs17190106, rs1056429 and rs143954678/rs113937848 associated with cervical malignancies overall, while rs17190106 and rs535777 associated specifically with invasive cancer (OR = 0.69, 95% CI = 0.55-0.86, p = 0.001) or adenocarcinomas (OR = 1.63, 95%CI = 1.17-2.27, p = 0.004), respectively. We tested these and one previously genotyped GWAS variant, rs9272117, for potential eQTL effects on 36 gene transcripts at the HLA locus in 280 cervical epithelial tissues. The strongest eQTL pairs were rs9272117 and HLA-DRB6 (p = 1.9x10E-5), rs1056429 and HLA-DRB5 (p = 2.5x10E-4), and rs535777 and HLA-DRB1 (p = 2.7x10E-4). We also identified transcripts that were specifically upregulated (DDX39B, HCP5, HLA-B, LTB, NFKBIL1) or downregulated (HLA-C, HLA-DPB2) in HPV+ or HPV16+ samples. In comparison, treating cervical epithelial cells with proinflammatory cytokine γ-IFN led to a dose-dependent induction of HCP5, HLA-B, HLA-C, HLA-DQB1, HLA-DRB1, HLA-DRB6, and repression of HSPA1L. Taken together, these results identify relevant genes from both the MHC class I and II regions that are inflammation-responsive in cervical epithelium and associate with HPV (HCP5, HLA-B, HLA-C) and/or with genomic cervical cancer risk variants (HLA-DRB1, HLA-DRB6). They may thus constitute important contributors to the immune escape of precancerous cells after HPV-infection.

Learning epistatic polygenic phenotypes with Boolean interactions

Detecting epistatic drivers of human phenotypes is a considerable challenge. Traditional approaches use regression to sequentially test multiplicative interaction terms involving pairs of genetic variants. For higher-order interactions and genome-wide large-scale data, this strategy is computationally intractable. Moreover, multiplicative terms used in regression modeling may not capture the form of biological interactions. Building on the Predictability, Computability, Stability (PCS) framework, we introduce the epiTree pipeline to extract higher-order interactions from genomic data using tree-based models. The epiTree pipeline first selects a set of variants derived from tissue-specific estimates of gene expression. Next, it uses iterative random forests (iRF) to search training data for candidate Boolean interactions (pairwise and higher-order). We derive significance tests for interactions, based on a stabilized likelihood ratio test, by simulating Boolean tree-structured null (no epistasis) and alternative (epistasis) distributions on hold-out test data. Finally, our pipeline computes PCS epistasis p-values that probabilisticly quantify improvement in prediction accuracy via bootstrap sampling on the test set. We validate the epiTree pipeline in two case studies using data from the UK Biobank: predicting red hair and multiple sclerosis (MS). In the case of predicting red hair, epiTree recovers known epistatic interactions surrounding MC1R and novel interactions, representing non-linearities not captured by logistic regression models. In the case of predicting MS, a more complex phenotype than red hair, epiTree rankings prioritize novel interactions surrounding HLA-DRB1, a variant previously associated with MS in several populations. Taken together, these results highlight the potential for epiTree rankings to help reduce the design space for follow up experiments.

HLA Haplotypes and Relapse After Hematopoietic Cell Transplantation

PURPOSE

Recurrence of blood malignancy is the major cause of hematopoietic cell transplant failure. HLA class II molecules play a fundamental role in antitumor responses but the role of class II haplotypes is not known.

METHODS

HLA-DR, -DQ, -DM, and -DO allele variation was determined in 1,629 related haploidentical transplants to study the clinical significance of individual molecules and haplotypes.

RESULTS

Outcome correlated with patient and donor variation for HLA-DRβ residue 86 (Gly/Val), HLA-DQ (G1/G2) heterodimers, and donor HLA-DM (DM11,11/nonDM11,11) molecules, and depended on patient-donor mismatching. Risks of relapse were lower for DRβ-86 GlyGly patients when the donor was GlyVal (hazard ratio [HR], 0.46 [95% CI, 0.30 to 0.68]; P < .001); GlyVal patients benefited from HLA-DRB1-matched donors, whereas no donor was superior to another for ValVal patients. G1G2 patients with G1G2-mismatched donors had lower relapse. Transplantation from donors with DMα residue 184 ArgHis was associated with higher risk of relapse (HR, 1.60 [95% CI, 1.09 to 2.36]; P = .02) relative to ArgArg. Relapse and mortality risks differed across HLA-DR-DQ-DM haplotypes.

CONCLUSION

HLA class II haplotypes may be functional constituents of the transplantation barrier, and their consideration in patients and donors may improve the success of transplantation.

Elevated genetic risk for multiple sclerosis emerged in steppe pastoralist populations

Multiple sclerosis (MS) is a neuro-inflammatory and neurodegenerative disease that is most prevalent in Northern Europe. Although it is known that inherited risk for MS is located within or in close proximity to immune-related genes, it is unknown when, where and how this genetic risk originated1. Here, by using a large ancient genome dataset from the Mesolithic period to the Bronze Age2, along with new Medieval and post-Medieval genomes, we show that the genetic risk for MS rose among pastoralists from the Pontic steppe and was brought into Europe by the Yamnaya-related migration approximately 5,000 years ago. We further show that these MS-associated immunogenetic variants underwent positive selection both within the steppe population and later in Europe, probably driven by pathogenic challenges coinciding with changes in diet, lifestyle and population density. This study highlights the critical importance of the Neolithic period and Bronze Age as determinants of modern immune responses and their subsequent effect on the risk of developing MS in a changing environment.

cLD: Rare-variant linkage disequilibrium between genomic regions identifies novel genomic interactions

Linkage disequilibrium (LD) is a fundamental concept in genetics; critical for studying genetic associations and molecular evolution. However, LD measurements are only reliable for common genetic variants, leaving low-frequency variants unanalyzed. In this work, we introduce cumulative LD (cLD), a stable statistic that captures the rare-variant LD between genetic regions, which reflects more biological interactions between variants, in addition to lack of recombination. We derived the theoretical variance of cLD using delta methods to demonstrate its higher stability than LD for rare variants. This property is also verified by bootstrapped simulations using real data. In application, we find cLD reveals an increased genetic association between genes in 3D chromatin interactions, a phenomenon recently reported negatively by calculating standard LD between common variants. Additionally, we show that cLD is higher between gene pairs reported in interaction databases, identifies unreported protein-protein interactions, and reveals interacting genes distinguishing case/control samples in association studies.

Global analysis of suppressor mutations that rescue human genetic defects

BACKGROUND

Genetic suppression occurs when the deleterious effects of a primary "query" mutation, such as a disease-causing mutation, are rescued by a suppressor mutation elsewhere in the genome.

METHODS

To capture existing knowledge on suppression relationships between human genes, we examined 2,400 published papers for potential interactions identified through either genetic modification of cultured human cells or through association studies in patients.

RESULTS

The resulting network encompassed 476 unique suppression interactions covering a wide spectrum of diseases and biological functions. The interactions frequently linked genes that operate in the same biological process. Suppressors were strongly enriched for genes with a role in stress response or signaling, suggesting that deleterious mutations can often be buffered by modulating signaling cascades or immune responses. Suppressor mutations tended to be deleterious when they occurred in absence of the query mutation, in apparent contrast with their protective role in the presence of the query. We formulated and quantified mechanisms of genetic suppression that could explain 71% of interactions and provided mechanistic insight into disease pathology. Finally, we used these observations to predict suppressor genes in the human genome.

CONCLUSIONS

The global suppression network allowed us to define principles of genetic suppression that were conserved across diseases, model systems, and species. The emerging frequency of suppression interactions among human genes and range of underlying mechanisms, together with the prevalence of suppression in model organisms, suggest that compensatory mutations may exist for most genetic diseases.

An overall view of the most common experimental models for multiple sclerosis

Multiple sclerosis (MS) is a complex chronic disease with an unknown etiology. It is considered an inflammatory demyelinating and neurodegenerative disorder of the central nervous system (CNS) characterized, in most cases, by an unpredictable onset of relapse and remission phases. The disease generally starts in subjects under 40; it has a higher incidence in women and is described as a multifactorial disorder due to the interaction between genetic and environmental risk factors. Unfortunately, there is currently no definitive cure for MS. Still, therapies can modify the disease's natural history, reducing the relapse rate and slowing the progression of the disease or managing symptoms. The limited access to human CNS tissue slows down. It limits the progression of research on MS. This limit has been partially overcome over the years by developing various experimental models to study this disease. Animal models of autoimmune demyelination, such as experimental autoimmune encephalomyelitis (EAE) and viral and toxin or transgenic MS models, represent the most significant part of MS research approaches. These models have now been complemented by ex vivo studies, using organotypic brain slice cultures and in vitro, through induced Pluripotent Stem cells (iPSCs). We will discuss which clinical features of the disorders might be reproduced and investigated in vivo, ex vivo, and in vitro in models commonly used in MS research to understand the processes behind the neuropathological events occurring in the CNS of MS patients. The primary purpose of this review is to give the reader a global view of the main paradigms used in MS research, spacing from the classical animal models to transgenic mice and 2D and 3D cultures.

Transcriptome-wide gene-gene interaction associations elucidate pathways and functional enrichment of complex traits

It remains unknown to what extent gene-gene interactions contribute to complex traits. Here, we introduce a new approach using predicted gene expression to perform exhaustive transcriptome-wide interaction studies (TWISs) for multiple traits across all pairs of genes expressed in several tissue types. Using imputed transcriptomes, we simultaneously reduce the computational challenge and improve interpretability and statistical power. We discover (in the UK Biobank) and replicate (in independent cohorts) several interaction associations, and find several hub genes with numerous interactions. We also demonstrate that TWIS can identify novel associated genes because genes with many or strong interactions have smaller single-locus model effect sizes. Finally, we develop a method to test gene set enrichment of TWIS associations (E-TWIS), finding numerous pathways and networks enriched in interaction associations. Epistasis is may be widespread, and our procedure represents a tractable framework for beginning to explore gene interactions and identify novel genomic targets.

I-Ag7 β56/57 polymorphisms regulate non-cognate negative selection to CD4+ T cell orchestrators of type 1 diabetes

Genetic susceptibility to type 1 diabetes is associated with homozygous expression of major histocompatibility complex class II alleles that carry specific beta chain polymorphisms. Why heterozygous expression of these major histocompatibility complex class II alleles does not confer a similar predisposition is unresolved. Using a nonobese diabetic mouse model, here we show that heterozygous expression of the type 1 diabetes-protective allele I-Ag7 β56P/57D induces negative selection to the I-Ag7-restricted T cell repertoire, including beta-islet-specific CD4+ T cells. Surprisingly, negative selection occurs despite I-Ag7 β56P/57D having a reduced ability to present beta-islet antigens to CD4+ T cells. Peripheral manifestations of non-cognate negative selection include a near complete loss of beta-islet-specific CXCR6+ CD4+ T cells, an inability to cross-prime islet-specific glucose-6-phosphatase catalytic subunit-related protein and insulin-specific CD8+ T cells and disease arrest at the insulitis stage. These data reveal that negative selection on non-cognate self-antigens in the thymus can promote T cell tolerance and protection from autoimmunity.

Genetics of Primary Biliary Cholangitis

Primary biliary cholangitis (PBC) is a rare disease of the liver characterized by an autoimmune attack on the small bile ducts. PBC is a complex trait, meaning that a large list of genetic factors interacts with environmental agents to determine its onset. Genome-wide association studies have had a huge impact in fostering research in PBC, but many steps need still to be done compared with other autoimmune diseases of similar prevalence. This review presents the state-of-the-art regarding the genetic architecture of PBC and provides some thoughtful reflections about possible future lines of research, which can be helpful to fill the missing heritability gap in PBC.

Missing Causality and Heritability of Autoimmune Hepatitis

BACKGROUND

Autoimmune hepatitis has an unknown cause and genetic associations that are not disease-specific or always present. Clarification of its missing causality and heritability could improve prevention and management strategies.

AIMS

Describe the key epigenetic and genetic mechanisms that could account for missing causality and heritability in autoimmune hepatitis; indicate the prospects of these mechanisms as pivotal factors; and encourage investigations of their pathogenic role and therapeutic potential.

METHODS

English abstracts were identified in PubMed using multiple key search phases. Several hundred abstracts and 210 full-length articles were reviewed.

RESULTS

Environmental induction of epigenetic changes is the prime candidate for explaining the missing causality of autoimmune hepatitis. Environmental factors (diet, toxic exposures) can alter chromatin structure and the production of micro-ribonucleic acids that affect gene expression. Epistatic interaction between unsuspected genes is the prime candidate for explaining the missing heritability. The non-additive, interactive effects of multiple genes could enhance their impact on the propensity and phenotype of autoimmune hepatitis. Transgenerational inheritance of acquired epigenetic marks constitutes another mechanism of transmitting parental adaptations that could affect susceptibility. Management strategies could range from lifestyle adjustments and nutritional supplements to precision editing of the epigenetic landscape.

CONCLUSIONS

Autoimmune hepatitis has a missing causality that might be explained by epigenetic changes induced by environmental factors and a missing heritability that might reflect epistatic gene interactions or transgenerational transmission of acquired epigenetic marks. These unassessed or under-evaluated areas warrant investigation.

Concurrent outcomes from multiple approaches of epistasis analysis for human body mass index associated loci provide insights into obesity biology

Genome wide association studies (GWAS) have focused on elucidating the genetic architecture of complex traits by assessing single variant effects in additive genetic models, albeit explaining a fraction of the trait heritability. Epistasis has recently emerged as one of the intrinsic mechanisms that could explain part of this missing heritability. We conducted epistasis analysis for genome-wide body mass index (BMI) associated SNPs in Alzheimer’s Disease Neuroimaging Initiative (ADNI) and followed up top significant interacting SNPs for replication in the UK Biobank imputed genotype dataset. We report two pairwise epistatic interactions, between rs2177596 (RHBDD1) and rs17759796 (MAPK1), rs1121980 (FTO) and rs6567160 (MC4R), obtained from a consensus of nine different epistatic approaches. Gene interaction maps and tissue expression profiles constructed for these interacting loci highlights co-expression, co-localisation, physical interaction, genetic interaction, and shared pathways emphasising the neuronal influence in obesity and implicating concerted expression of associated genes in liver, pancreas, and adipose tissues insinuating to metabolic abnormalities characterized by obesity. Detecting epistasis could thus be a promising approach to understand the effect of simultaneously interacting multiple genetic loci in disease aetiology, beyond single locus effects.

Genomics of Ancient Pathogens: First Advances and Prospects

Paleogenomics is one of the urgent and promising areas of interdisciplinary research in the today's world science. New genomic methods of ancient DNA (aDNA) analysis, such as next generation sequencing (NGS) technologies, make it possible not only to obtain detailed genetic information about historical and prehistoric human populations, but also to study individual microbial and viral pathogens and microbiomes from different ancient and historical objects. Studies of aDNA of pathogens by reconstructing their genomes have so far yielded complete sequences of the ancient pathogens that played significant role in the history of the world: Yersinia pestis (plague), Variola virus (smallpox), Vibrio cholerae (cholera), HBV (hepatitis B virus), as well as the equally important endemic human infectious agents: Mycobacterium tuberculosis (tuberculosis), Mycobacterium leprae (leprosy), and Treponema pallidum (syphilis). Genomic data from these pathogens complemented the information previously obtained by paleopathologists and allowed not only to identify pathogens from the past pandemics, but also to recognize the pathogen lineages that are now extinct, to refine chronology of the pathogen appearance in human populations, and to reconstruct evolutionary history of the pathogens that are still relevant to public health today. In this review, we describe state-of-the-art genomic research of the origins and evolution of many ancient pathogens and viruses and examine mechanisms of the emergence and spread of the ancient infections in the mankind history.

Experimental Autoimmune Encephalomyelitis in the Mouse

This article details the materials and methods required for both active induction and adoptive transfer of experimental autoimmune encephalomyelitis (EAE) in the SJL mouse strain using intact proteins or peptides from the two major myelin proteins: proteolipid protein (PLP) and myelin basic protein (MBP). Additionally, active induction of EAE in the C57BL/6 strain using myelin oligodendrocyte glycoprotein (MOG) peptide is also discussed. Detailed materials and methods required for the purification of both PLP and MBP are described, and a protocol for isolating CNS-infiltrating lymphocytes in EAE mice is included. Modifications of the specified protocols may be necessary for efficient induction of active or adoptive EAE in other mouse strains. © 2021 Wiley Periodicals LLC. Basic Protocol: Active induction of EAE with PLP, MBP, and MOG protein or peptide Alternate Protocol: Adoptive induction of EAE with PLP-, MBP-, or MOG-specific lymphocytes Support Protocol 1: Purification of proteolipid protein Support Protocol 2: Purification of myelin basic protein Support Protocol 3: Isolation of CNS-infiltrating lymphocytes.

Human leukocyte antigen class I association with occult hepatitis B virus infection in the Shaanxi Han group: Analysis at the haplotype level

BACKGROUND

Human leucocyte antigen (HLA) class I genes and haplotypes correlate with hepatitis B virus (HBV) infection. Occult HBV infection (OBI), a special type of chronic HBV infection, is defined as HBV surface antigen negative patients with or those without serologic markers by the means of HBV DNA detection in human plasma or in liver tissue by a diagnostic test. So far, the associations of HLA I haplotypes with OBI have not been reported previously in China.

METHODS

A case-control study between 107 OBI subjects and 280 healthy controls from blood donors in the Blood Center of the Shaanxi Province was conducted in the present association analysis. The HLA-A, -B and -C loci of case-control subjects were detected and genotyped by polymerase chain reaction-sequence based typing. The HLA-A, -B and -C haplotypic frequencies were calculated by the maximum likelihood method.

RESULTS

The HLA-A*33:03-C*07:01G (p_c  = 0.039, odds ratio [OR] = 8.996, 95% confidence interval [CI] = 1.825-44.338), B*44:03-C*07:01G (p_c  = 0.0069, OR = 12.000, 95% CI = 2.507-57.436) and A*33:03-B*44:03-C*07:01G (p_c  = 0.04, OR = 7.094, 95% CI = 1.387-36.288) haplotypes showed a a significant positive association with OBI. Independent effects demonstrated that HLA-B*44:03 and HLA-C*07:01G gave the main contribution to risk, whereas HLA-A*33:03 was associated only by linkage disequilibrium.

CONCLUSIONS

This present study is the first to demonstrate that HLA I haplotypes are associated with OBI in the Shaanxi Han population. The present results suggest that HLA-B*44:03-C*07:01G might be a potential risk factor for OBI. Comparisons of the frequencies of HLA I haplotypes at high resolution were made between OBI from northern Chinese Han and controls. Frequencies of HLA-A-B, A-C, B-C, A-B-C between the OBI group (gray column) and the control group (blank column) were compared. Haplotypes with p < 0.05 in either group were shown. p values for multiple comparisons (p_c ) were corrected by Bonferroni correction (*p_c  < 0.05).

The genetic architecture of primary biliary cholangitis

Primary biliary cholangitis (PBC) is a rare autoimmune disease of the liver affecting the small bile ducts. From a genetic point of view, PBC is a complex trait and several genetic and environmental factors have been called in action to explain its etiopathogenesis. Similarly to other complex traits, PBC has benefited from the introduction of genome-wide association studies (GWAS), which identified many variants predisposing or protecting toward the development of the disease. While a progressive endeavour toward the characterization of candidate loci and downstream pathways is currently ongoing, there is still a relatively large portion of heritability of PBC to be revealed. In addition, genetic variation behind progression of the disease and therapeutic response are mostly to be investigated yet. This review outlines the state-of-the-art regarding the genetic architecture of PBC and provides some hints for future investigations, focusing on the study of gene-gene interactions, the application of whole-genome sequencing techniques, and the investigation of X chromosome that can be helpful to cover the missing heritability gap in PBC.

Toward an Information Theory of Quantitative Genetics

Quantitative genetics has evolved dramatically in the past century, and the proliferation of genetic data, in quantity as well as type, enables the characterization of complex interactions and mechanisms beyond the scope of its theoretical foundations. In this article, we argue that revisiting the framework for analysis is important and we begin to lay the foundations of an alternative formulation of quantitative genetics based on information theory. Information theory can provide sensitive and unbiased measures of statistical dependencies among variables, and it provides a natural mathematical language for an alternative view of quantitative genetics. In the previous work, we examined the information content of discrete functions and applied this approach and methods to the analysis of genetic data. In this article, we present a framework built around a set of relationships that both unifies the information measures for the discrete functions and uses them to express key quantitative genetic relationships. Information theory measures of variable interdependency are used to identify significant interactions, and a general approach is described for inferring functional relationships in genotype and phenotype data. We present information-based measures of the genetic quantities: penetrance, heritability, and degrees of statistical epistasis. Our scope here includes the consideration of both two- and three-variable dependencies and independently segregating variants, which captures additive effects, genetic interactions, and two-phenotype pleiotropy. This formalism and the theoretical approach naturally apply to higher multivariable interactions and complex dependencies, and can be adapted to account for population structure, linkage, and nonrandomly segregating markers. This article thus focuses on presenting the initial groundwork for a full formulation of quantitative genetics based on information theory.

Multiple sclerosis: doubling down on MHC

Human leukocyte antigen (HLA)-encoded surface molecules present antigenic peptides to T lymphocytes and play a key role in adaptive immune responses. Besides their physiological role of defending the host against infectious pathogens, specific alleles serve as genetic risk factors for autoimmune diseases. For multiple sclerosis (MS), an autoimmune disease that affects the brain and spinal cord, an association with the HLA-DR15 haplotype was described in the early 1970s. This short opinion piece discusses the difficulties of disentangling the details of this association and recent observations about the functional involvement of not only one, but also the second gene of the HLA-DR15 haplotype. This information is not only important for understanding the pathomechanism of MS, but also for antigen-specific therapies.

Association of HLA-DR2-Related Haplotype (HLA-DRB5*01-DRB1*1501-DQB1*0602) in Patients with Multiple Sclerosis in Khuzestan Province

OBJECTIVE

Multiple sclerosis (MS) is a partially heritable autoimmune disease. HLA-DR2 is the largest identified genetic risk factor for MS. The largest identified genetic risk factor is haplotype from the MHC class II HLA-DR2, which increases the disease risk. The HLA-DR2 distribution in MS patients has been confirmed, but contradictory outcomes have been found. Moreover, the HLA-DR2 effect on ethnicity and gender is unclear. There are no data regarding the HLA-DR2 (HLA-DRB1*1501-DRB5*01-DQB1*0602) association with MS in Khuzestan Province, Iran. This study aimed to investigate the association of HLA-DR2 with MS regarding both sex and ethnicity in this province.

MATERIALS & METHODS

A total of 399 individuals were recruited. HLA typing was conducted using the polymerase chain reaction amplification with sequence-specific primers technology. The HLA-DR2 association with MS was analyzed, and also its probable association with gender, ethnicity, the expanded disability status scale (EDSS), and MS clinical course was examined using the Chi-square test.

RESULTS

HLA-DRB5*01 - -DQB1*0602 - as the most common HLA haplotype was found in both patient and control groups. In contrast, the DRB5*01 + -DRB1*1501 + -DQB1*0602 - frequency was very low in the groups. It was observed that haplotypes had no association with MS susceptibility. Most of the haplotypes showed no association with ethnicity, sex, EDSS, and MS course except for the HLA-DRB5*01 + -DRB1*1501 + -DQB1*0602 - haplotype that was positively associated with EDSS steps 5 to 10 (p=0.014) and non-RRMS (p=0.023).

CONCLUSION

There was no association between HLA-DR2 and MS susceptibility. However, the higher HLA-DRB5*01 + -DRB1*1501 + -DQB1*0602 - frequency may play a role in MS development. Also, HLA-DR2 did not increase significantly concerning clinical course, ethnicity, sex, and EDSS. This study further supports the importance of replication studies as susceptible loci that might differ in various ethnicities. Therefore, it is concluded that the association between HLA-DR2 and MS is more allelic than haplotypic in Khuzestan.

Genomic, proteomic, and systems biology approaches in biomarker discovery for multiple sclerosis

Multiple sclerosis (MS) is a neuroinflammatory disorder characterized by autoimmune-mediated inflammatory lesions in CNS leading to myelin damage and axonal loss. MS is a heterogenous disease with variable and unpredictable disease course. Due to its complex nature, MS is difficult to diagnose and responses to specific treatments may vary between individuals. Therefore, there is an indisputable need for biomarkers for early diagnosis, prediction of disease exacerbations, monitoring the progression of disease, and for measuring responses to therapy. Genomic and proteomic studies have sought to understand the molecular basis of MS and find biomarker candidates. Advances in next-generation sequencing and mass-spectrometry techniques have yielded an unprecedented amount of genomic and proteomic data; yet, translation of the results into the clinic has been underwhelming. This has prompted the development of novel data science techniques for exploring these large datasets to identify biologically relevant relationships and ultimately point towards useful biomarkers. Herein we discuss optimization of omics study designs, advances in the generation of omics data, and systems biology approaches aimed at improving biomarker discovery and translation to the clinic for MS.

HLA-DR15 Molecules Jointly Shape an Autoreactive T Cell Repertoire in Multiple Sclerosis

The HLA-DR15 haplotype is the strongest genetic risk factor for multiple sclerosis (MS), but our understanding of how it contributes to MS is limited. Because autoreactive CD4⁺ T cells and B cells as antigen-presenting cells are involved in MS pathogenesis, we characterized the immunopeptidomes of the two HLA-DR15 allomorphs DR2a and DR2b of human primary B cells and monocytes, thymus, and MS brain tissue. Self-peptides from HLA-DR molecules, particularly from DR2a and DR2b themselves, are abundant on B cells and thymic antigen-presenting cells. Furthermore, we identified autoreactive CD4⁺ T cell clones that can cross-react with HLA-DR-derived self-peptides (HLA-DR-SPs), peptides from MS-associated foreign agents (Epstein-Barr virus and Akkermansia muciniphila), and autoantigens presented by DR2a and DR2b. Thus, both HLA-DR15 allomorphs jointly shape an autoreactive T cell repertoire by serving as antigen-presenting structures and epitope sources and by presenting the same foreign peptides and autoantigens to autoreactive CD4⁺ T cells in MS.

•

HLA-DR15 present abundant HLA-DR-derived self-peptides on B cells

•

Autoreactive T cells in MS recognize HLA-DR-derived self-peptides/DR15 complexes

•

Foreign peptides/DR15 complexes trigger potential autoreactive T cells in MS

•

HLA-DR15 shape an autoreactive T cell repertoire by cross-reactivity/restriction

Gene Expression in Spontaneous Experimental Autoimmune Encephalomyelitis Is Linked to Human Multiple Sclerosis Risk Genes

Recent genome-wide association studies have identified over 230 genetic risk loci for multiple sclerosis. Current experimental autoimmune encephalomyelitis (EAE) models requiring active induction of disease may not be optimally suited for the characterization of the function of these genes. We have thus used gene expression profiling to study whether spontaneous opticospinal EAE (OSE) or MOG-induced EAE mirrors the genetic contribution to the pathogenesis of multiple sclerosis more faithfully. To this end, we compared gene expression in OSE and MOG EAE models and analyzed the relationship of both models to human multiple sclerosis risk genes and T helper cell biology. We observed stronger gene expression changes and an involvement of more pathways of the adaptive immune system in OSE than MOG EAE. Furthermore, we demonstrated a more extensive enrichment of human MS risk genes among transcripts differentially expressed in OSE than was the case for MOG EAE. Transcripts differentially expressed only in diseased OSE mice but not in MOG EAE were significantly enriched for T helper cell-specific transcripts. These transcripts are part of immune-regulatory pathways. The activation of the adaptive immune system and the enrichment of both human multiple sclerosis risk genes and T helper cell-specific transcripts were also observed in OSE mice showing only mild disease signs. These expression changes may, therefore, be indicative of processes at disease onset. In summary, more human multiple sclerosis risk genes were differentially expressed in OSE than was observed for MOG EAE, especially in T_H1 cells. When studying the functional role of multiple sclerosis risk genes and pathways during disease onset and their interactions with the environment, spontaneous OSE may thus show advantages over MOG-induced EAE.

HLA transgenic mice: application in reproducing idiosyncratic drug toxicity

Various types of transgenic mice carrying either class I or II human leukocyte antigen (HLA) molecules are readily available, and reports describing their use in a variety of studies have been published for more than 30 years. Examples of their use include the discovery of HLA-specific antigens against viral infection as well as the reproduction of HLA-mediated autoimmune diseases for the development of therapeutic strategies. Recently, HLA transgenic mice have been used to reproduce HLA-mediated idiosyncratic drug toxicity (IDT), a rare and unpredictable adverse drug reaction that can result in death. For example, abacavir-induced IDT has successfully been reproduced in HLA-B*57:01 transgenic mice. Several reports using HLA transgenic mice for IDT have proven the utility of this concept for the evaluation of IDT using various HLA allele combinations and drugs. It has become apparent that such models may be a valuable tool to investigate the mechanisms underlying HLA-mediated IDT. This review summarizes the latest findings in the area of HLA transgenic mouse models and discusses the current challenges that must be overcome to maximize the potential of this unique animal model.

Multiple Sclerosis: Lipids, Lymphocytes, and Vitamin D

Multiple sclerosis (MS) is an inflammatory demyelinating disease of the central nervous system. We review the two core MS features, myelin instability, fragmentation, and remyelination failure, and dominance of pathogenic CD4+ Th17 cells over protective CD4+ Treg cells. To better understand myelin pathology, we describe myelin biosynthesis, structure, and function, then highlight stearoyl-CoA desaturase (SCD) in nervonic acid biosynthesis and nervonic acid's contribution to myelin stability. Noting that vitamin D deficiency decreases SCD in the periphery, we propose it also decreases SCD in oligodendrocytes, disrupting the nervonic acid supply and causing myelin instability and fragmentation. To better understand the distorted Th17/Treg cell balance, we summarize Th17 cell contributions to MS pathogenesis, then highlight how 1,25-dihydroxyvitamin D3 signaling from microglia to CD4+ T cells restores Treg cell dominance. This signaling rapidly increases flux through the methionine cycle, removing homocysteine, replenishing S-adenosyl-methionine, and improving epigenetic marking. Noting that DNA hypomethylation and inappropriate DRB1*1501 expression were observed in MS patient CD4+ T cells, we propose that vitamin D deficiency thwarts epigenetic downregulation of DRB1*1501 and Th17 cell signature genes, and upregulation of Treg cell signature genes, causing dysregulation within the CD4+ T cell compartment. We explain how obesity reduces vitamin D status, and how estrogen and vitamin D collaborate to promote Treg cell dominance in females. Finally, we discuss the implications of this new knowledge concerning myelin and the Th17/Treg cell balance, and advocate for efforts to address the global epidemics of obesity and vitamin D deficiency in the expectation of reducing the impact of MS.

Discriminative T cell recognition of cross-reactive islet-antigens is associated with HLA-DQ8 transdimer-mediated autoimmune diabetes

Human leukocyte antigen (HLA)-DQ8 transdimer (HLA-DQA1*0501/DQB1*0302) confers exceptionally high risk in autoimmune diabetes. However, little is known about HLA-DQ8 transdimer-restricted CD4 T cell recognition, an event crucial for triggering HLA-DQ8 transdimer-specific anti-islet immunity. Here, we report a high degree of epitope overlap and T cell promiscuity between susceptible HLA-DQ8 and HLA-DQ8 transdimer. Despite preservation of putative residues for T cell receptor (TCR) contact, stronger disease-associated responses to cross-reactive, immunodominant islet epitopes are elicited by HLA-DQ8 transdimer. Mutagenesis at the α chain of HLA-DQ8 transdimer in complex with the disease-relevant GAD65_250-266 peptide and in silico analysis reveal the DQ α52 residue located within the N-terminal edge of the peptide-binding cleft for the enhanced T cell reactivity, altering avidity and biophysical affinity between TCR and HLA-peptide complexes. Accordingly, a structurally promiscuous but nondegenerate TCR-HLA-peptide interface is pivotal for HLA-DQ8 transdimer-mediated autoimmune diabetes.

Complete nucleotide sequence characterization of DRB5 alleles reveals a homogeneous allele group that is distinct from other DRB genes

Next Generation Sequencing allows for testing and typing of entire genes of the HLA region. A better and comprehensive sequence assessment can be achieved by the inclusion of full gene sequences of all the common alleles at a given locus. The common alleles of DRB5 are under-characterized with the full exon-intron sequence of two alleles available. In the present study the DRB5 genes from 18 subjects alleles were cloned and sequenced; haplotype analysis showed that 17 of them had a single copy of DRB5 and one consanguineous subject was homozygous at all HLA loci. Methodological approaches including robust and efficient long-range PCR amplification, molecular cloning, nucleotide sequencing and de novo sequence assembly were combined to characterize DRB5 alleles. DRB5 sequences covering from 5'UTR to the end of intron 5 were obtained for DRB5*01:01, 01:02 and 02:02; partial coverage including a segment spanning exon 2 to exon 6 was obtained for DRB5*01:03, 01:08N and 02:03. Phylogenetic analysis of the generated sequences showed that the DRB5 alleles group together and have distinctive differences with other DRB loci. Novel intron variants of DRB5*01:01:01, 01:02 and 02:02 were identified. The newly characterized DRB5 intron variants of each DRB5 allele were found in subjects harboring distinct associations with alleles of DRB1, B and/or ethnicity. The new information provided by this study provides reference sequences for HLA typing methodologies. Extending sequence coverage may lead to identify the disease susceptibility factors of DRB5 containing haplotypes while the unexpected intron variations may shed light on understanding of the evolution of the DRB region.

Progress of immunotherapy of anti-α-synuclein in Parkinson's disease

Many neurodegenerative diseases are characterized by progressive loss of neurons and abnormal protein accumulation, including amyloid (A)β and tau in Alzheimer's disease and Lewy bodies and α-synuclein (α-syn) in Parkinson's disease (PD). Recent evidence suggests that adaptive immunity plays an important role in PD, and that anti-α-syn antibodies can be used as therapy in neurodegenerative diseases; monoclonal antibodies were shown to inhibit α-syn propagation and aggregation in PD models and patients. In this review, we summarize the different pathological states of α-syn, including gene mutations, truncation, phosphorylation, and the high molecular weight form, and describe the specific antibodies that recognize the α-syn monomer or oligomer, some of which have been tested in clinic trials. We also discuss future research directions and potential targets in PD therapy.

PR3-ANCA-associated vasculitis is associated with a specific motif in the peptide-binding cleft of HLA-DP molecules

Objectives

This study aimed to characterize the association between HLA alleles and ANCA-associated vasculitis (AAV) in a genetically homogeneous population, and to analyse the contribution of specific HLA molecule amino acid sequences to the risk of AAV.

Methods

We included 187 Danish patients with AAV and 1070 healthy controls. All were HLA typed at two-field resolution. The association of HLA alleles to PR3- or MPO-AAV was analysed. The contribution of the dominant molecular motifs of the HLA-DPB1 molecule to the risk of AAV was investigated by association studies that included specific amino acid sequences of the hypervariable regions in exon 2.

Results

Ninety-four percent of patients with PR3-AAV were carriers of HLA-DPB1*04:01 while all patients with PR3-AAV were carriers of an HLA-DPB1*04 allele, and 85% were homozygous. This was significantly more than in the control group (P < 0.0001). The association was even stronger when HLA-DPB1*04:02 and -DPB1*23:01 were included. HLA-DPB1*04:01, -DPB1*04:02 and -DPB1*23:01 share amino acids in positions 8–9, 69, 76 and 84–87 within the hypervariable regions, but only positions 69 and 84–87 contributed significantly to the disease risk. HLA-DRB1*15 was associated with an increased risk of developing PR3-AAV, while HLA-DRB1*04, -DRB1*07 and -DQB1*03 were associated with a reduced risk of kidney involvement in PR3-AAV. MPO-AAV was only weakly associated with HLA class I alleles.

Conclusion

PR3-AAV is strongly associated with the HLA-DPB1 alleles HLA-DPB1*04:01, -DPB1*04:02 and -DPB1*23:01, which share amino acid sequences crucial for the peptide-binding groove.

Experimental Autoimmune Encephalomyelitis (EAE) as Animal Models of Multiple Sclerosis (MS)

Multiple sclerosis (MS) is a multifocal demyelinating disease of the central nervous system (CNS) leading to the progressive destruction of the myelin sheath surrounding axons. It can present with variable clinical and pathological manifestations, which might reflect the involvement of distinct pathogenic processes. Although the mechanisms leading to the development of the disease are not fully understood, numerous evidences indicate that MS is an autoimmune disease, the initiation and progression of which are dependent on an autoimmune response against myelin antigens. In addition, genetic susceptibility and environmental triggers likely contribute to the initiation of the disease. At this time, there is no cure for MS, but several disease-modifying therapies (DMTs) are available to control and slow down disease progression. A good number of these DMTs were identified and tested using animal models of MS referred to as experimental autoimmune encephalomyelitis (EAE). In this review, we will recapitulate the characteristics of EAE models and discuss how they help shed light on MS pathogenesis and help test new treatments for MS patients.

CTLA-4, PD-1 and TIM-3 expression predominantly downregulated in MS patients

Multiple sclerosis (MS), is an autoimmune disorder of central nervous system (CNS) characterized by inflammation and demyelination. Self-tolerance impairment is considered to be induced by a combination of inherited susceptibility and environmental agents. In this work, we demonstrate that a reduction in the comparative expression of well-known inhibitory receptors (i.e., CTLA-4, PD-1, and TIM-3) is importantly linked with MS patients compared to healthy controls. The relative expression of interested genes was performed on peripheral blood mononuclear cells (PBMCs), by using quantitative real time-PCR (qRT-PCR). Our data highlighted the role of inhibitory receptors in the maintenance of immune homeostasis in autoimmune disease.

Gene buddies: linked balanced polymorphisms reinforce each other even in the absence of epistasis

The fates of genetic polymorphisms maintained by balancing selection depend on evolutionary dynamics at linked sites. While coevolution across linked, epigenetically-interacting loci has been extensively explored, such supergenes may be relatively rare. However, genes harboring adaptive variation can occur in close physical proximity while generating independent effects on fitness. Here, I present a model in which two linked loci without epistasis are both under balancing selection for unrelated reasons. Using forward-time simulations, I show that recombination rate strongly influences the retention of adaptive polymorphism, especially for intermediate selection coefficients. A locus is more likely to retain adaptive variation if it is closely linked to another locus under balancing selection, even if the two loci have no interaction. Thus, two linked polymorphisms can both be retained indefinitely even when they would both be lost to drift if unlinked. While these results may be intuitive, they have important implications for genetic architecture: clusters of mutually reinforcing genes may underlie phenotypic variation in natural populations, and such genes cannot be assumed to be functionally associated. Future studies that measure selection coefficients and recombination rates among closely linked genes will be fruitful for characterizing the extent of this phenomenon.

Goodpasture's autoimmune disease - A collagen IV disorder

Goodpasture's (GP) disease is an autoimmune disorder characterized by the deposition of pathogenic autoantibodies in basement membranes of kidney and lung eliciting rapidly progressive glomerulonephritis and pulmonary hemorrhage. The principal autoantigen is the α345 network of collagen IV, which expression is restricted to target tissues. Recent discoveries include a key role of chloride and bromide for network assembly, a novel posttranslational modification of the antigen, a sulfilimine bond that crosslinks the antigen, and the mechanistic role of HLA in genetic susceptibility and resistance to GP disease. These advances provide further insights into molecular mechanisms of initiation and progression of GP disease and serve as a basis for developing of novel diagnostic tools and therapies for treatment of Goodpasture's disease.

Ancient DNA study reveals HLA susceptibility locus for leprosy in medieval Europeans

Leprosy, a chronic infectious disease caused by Mycobacterium leprae (M. leprae), was very common in Europe till the 16th century. Here, we perform an ancient DNA study on medieval skeletons from Denmark that show lesions specific for lepromatous leprosy (LL). First, we test the remains for M. leprae DNA to confirm the infection status of the individuals and to assess the bacterial diversity. We assemble 10 complete M. leprae genomes that all differ from each other. Second, we evaluate whether the human leukocyte antigen allele DRB1*15:01, a strong LL susceptibility factor in modern populations, also predisposed medieval Europeans to the disease. The comparison of genotype data from 69 M. leprae DNA-positive LL cases with those from contemporary and medieval controls reveals a statistically significant association in both instances. In addition, we observe that DRB1*15:01 co-occurs with DQB1*06:02 on a haplotype that is a strong risk factor for inflammatory diseases today.

Leprosy, caused by infection with Mycobacterium leprae, was common in Europe in the Middle Ages. Here, Krause-Kyora et al. analyze ancient DNA from a medieval Danish leprosarium to assemble 10 complete bacterial genomes and perform association analysis of the DRB1*15:01 allele with risk of leprosy infection.

Deconstruction of HLA-DRB1*04:01:01 and HLA-DRB1*15:01:01 class II haplotypes using next-generation sequencing in European-Americans with multiple sclerosis

BACKGROUND

The association between HLA-DRB1*15:01 with multiple sclerosis (MS) susceptibility is well established, but the contribution of the tightly associated HLA-DRB5*01:01 allele has not yet been completely ascertained. Similarly, the effects of HLA-DRB1*04:01 alleles and haplotypes, defined at the full-gene resolution level with MS risk remains to be elucidated.

OBJECTIVES

To characterize the molecular architecture of class II HLA-DR15 and HLA-DR4 haplotypes associated with MS.

METHODS

Next-generation sequencing was used to determine HLA-DQB1, HLA-DQA1, and HLA-DRB1/4/5 alleles in 1403 unrelated European-American patients and 1425 healthy unrelated controls. Effect sizes of HLA alleles and haplotypes on MS risk were measured by odds ratio (OR) with 95% confidence intervals.

RESULTS

HLA-DRB1*15:01:01:01SG (OR = 3.20, p < 2.2E-16), HLA-DRB5*01:01:01 (OR = 2.96, p < 2.2E-16), and HLA-DRB5*01:01:01v1_STR1 (OR = 8.18, p = 4.3E-05) alleles all occurred at significantly higher frequencies in MS patients compared to controls. The most significant predis-posing haplotypes were HLA-DQB1*06:02:01~ HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01 and HLA-DQB1*06:02:01~HLA-DQA1*01:02:01:01SG~HLA-DRB1*15:01:01:01SG~HLA-DRB5*01:01:01v1_STR1 (OR = 3.19, p < 2.2E-16; OR = 9.30, p = 9.7E-05, respectively). Analyses of the HLA-DRB1*04 cohort in the absence of HLA-DRB1*15:01 haplotypes revealed that the HLA-DQB1*03:01:01:01~HLA-DQA1*03:03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was protective (OR = 0.64, p = 0.028), whereas the HLA-DQB1*03:02:01~HLA-DQA1*03:01:01~HLA-DRB1*04:01:01:01SG~HLA-DRB4*01:03:01:01 haplotype was associated with MS susceptibility (OR = 1.66, p = 4.9E-03).

CONCLUSION

HLA-DR15 haplotypes, including genomic variants of HLA-DRB5, and HLA-DR4 haplotypes affect MS risk.

HLA variation and disease

Fifty years since the first description of an association between HLA and human disease, HLA molecules have proven to be central to physiology, protective immunity and deleterious, disease-causing autoimmune reactivity. Technological advances have enabled pivotal progress in the determination of the molecular mechanisms that underpin the association between HLA genetics and functional outcome. Here, we review our current understanding of HLA molecules as the fundamental platform for immune surveillance and responsiveness in health and disease. We evaluate the scope for personalized antigen-specific disease prevention, whereby harnessing HLA-ligand interactions for clinical benefit is becoming a realistic prospect.

A genomic perspective on HLA evolution

Several decades of research have convincingly shown that classical human leukocyte antigen (HLA) loci bear signatures of natural selection. Despite this conclusion, many questions remain regarding the type of selective regime acting on these loci, the time frame at which selection acts, and the functional connections between genetic variability and natural selection. In this review, we argue that genomic datasets, in particular those generated by next-generation sequencing (NGS) at the population scale, are transforming our understanding of HLA evolution. We show that genomewide data can be used to perform robust and powerful tests for selection, capable of identifying both positive and balancing selection at HLA genes. Importantly, these tests have shown that natural selection can be identified at both recent and ancient timescales. We discuss how findings from genomewide association studies impact the evolutionary study of HLA genes, and how genomic data can be used to survey adaptive change involving interaction at multiple loci. We discuss the methodological developments which are necessary to correctly interpret genomic analyses involving the HLA region. These developments include adapting the NGS analysis framework so as to deal with the highly polymorphic HLA data, as well as developing tools and theory to search for signatures of selection, quantify differentiation, and measure admixture within the HLA region. Finally, we show that high throughput analysis of molecular phenotypes for HLA genes-namely transcription levels-is now a feasible approach and can add another dimension to the study of genetic variation.

Autoimmune aspects of psoriasis: Heritability and autoantigens

Chronic immune-mediated disorders (IMDs) constitute a major health burden. Understanding IMD pathogenesis is facing two major constraints: Missing heritability explaining familial clustering, and missing autoantigens. Pinpointing IMD risk genes and autoimmune targets, however, is of fundamental importance for developing novel causal therapies. The strongest association of all IMDs is seen with human leukocyte antigen (HLA) alleles. Using psoriasis as an IMD model this article reviews the pathogenic role HLA molecules may have within the polygenic predisposition of IMDs. It concludes that disease-associated HLA alleles account for both missing heritability and autoimmune mechanisms by facilitating tissue-specific autoimmune responses through autoantigen presentation.

Dominant protection from HLA-linked autoimmunity by antigen-specific regulatory T cells

Susceptibility and protection against human autoimmune diseases, including type I diabetes, multiple sclerosis, and Goodpasture disease, is associated with particular human leukocyte antigen (HLA) alleles. However, the mechanisms underpinning such HLA-mediated effects on self-tolerance remain unclear. Here we investigate the molecular mechanism of Goodpasture disease, an HLA-linked autoimmune renal disorder characterized by an immunodominant CD4⁺ T-cell self-epitope derived from the α3 chain of type IV collagen (α3_135-145). While HLA-DR15 confers a markedly increased disease risk, the protective HLA-DR1 allele is dominantly protective in trans with HLA-DR15 (ref. 2). We show that autoreactive α3_135-145-specific T cells expand in patients with Goodpasture disease and, in α3_135-145-immunized HLA-DR15 transgenic mice, α3_135-145-specific T cells infiltrate the kidney and mice develop Goodpasture disease. HLA-DR15 and HLA-DR1 exhibit distinct peptide repertoires and binding preferences and present the α3_135-145 epitope in different binding registers. HLA-DR15-α3_135-145 tetramer⁺ T cells in HLA-DR15 transgenic mice exhibit a conventional T-cell phenotype (T_conv) that secretes pro-inflammatory cytokines. In contrast, HLA-DR1-α3_135-145 tetramer⁺ T cells in HLA-DR1 and HLA-DR15/DR1 transgenic mice are predominantly CD4⁺Foxp3⁺ regulatory T cells (T_reg cells) expressing tolerogenic cytokines. HLA-DR1-induced T_reg cells confer resistance to disease in HLA-DR15/DR1 transgenic mice. HLA-DR15⁺ and HLA-DR1⁺ healthy human donors display altered α3_135-145-specific T-cell antigen receptor usage, HLA-DR15-α3_135-145 tetramer⁺ Foxp3^- T_conv and HLA-DR1-α3_135-145 tetramer⁺ Foxp3⁺CD25^hiCD127^lo T_reg dominant phenotypes. Moreover, patients with Goodpasture disease display a clonally expanded α3_135-145-specific CD4⁺ T-cell repertoire. Accordingly, we provide a mechanistic basis for the dominantly protective effect of HLA in autoimmune disease, whereby HLA polymorphism shapes the relative abundance of self-epitope specific T_reg cells that leads to protection or causation of autoimmunity.

Immunodominant T-Cell Epitopes in the VP1 Capsid Protein of Rhinovirus Species A and C

Rhinovirus (RV) species A and C are the most frequent cause of respiratory viral illness worldwide, and RV-C has been linked to more severe exacerbations of asthma in young children. Little is known about the immune responses to the different RV species, although studies comparing IgG1 antibody titers found impaired antibody responses to RV-C. Therefore, the aim of this study was to assess whether T-cell immunity to RV-C is similarly impaired. We measured T-cell proliferation to overlapping synthetic peptides covering the entire VP1 capsid protein of an RV-A and RV-C genotype for 20 healthy adult donors. Human leukocyte antigen (HLA) was typed in all the donors in order to investigate possible associations between the HLA type and RV peptide recognition. Total and specific IgG1 antibody titers to the VP1 proteins of both RV-A and RV-C were also measured to examine associations between the antibody and T-cell responses. We identified T-cell epitopes that are specific to and representative of each RV-A and RV-C species. These epitopes stimulated CD4⁺-specific T-cell proliferation, with similar magnitudes of response for both RV species. All the donors, independent of their HLA-DR or -DQ type, were able to recognize the immunodominant RV-A and -C regions of VP1. Furthermore, the presence or absence of specific antibody titers was not related to changes in T-cell recognition. Our results indicate a dissociation between the antibody and T-cell responses to rhinoviruses. The species-representative T-cell epitopes identified in this study are valuable tools for future studies investigating T-cell responses to the different RV species.

IMPORTANCE

Rhinoviruses (RVs) are mostly associated with the common cold and asthma exacerbations, although their contributions to most upper and lower respiratory tract diseases have increasingly been reported. Species C (RV-C) has been associated with more frequent and severe asthma exacerbations in young children and, along with RV-A, is the most clinically relevant species. Little is known about how our immune system responds to rhinoviruses, and there are limited tools to study specific adaptive immunity against each rhinovirus species. In this study, we identified immunodominant T-cell epitopes of the VP1 proteins of RV-A and RV-C, which are representative of each species. The study found that T-cell responses to RV-A and RV-C were of similar magnitudes, in contrast with previous findings showing RV-C-specific antibody responses were low. These findings will provide the basis for future studies on the immune response to rhinoviruses and can help elucidate the mechanisms of severity of rhinovirus-induced infections.

Atlas of human diseases influenced by genetic variants with extreme allele frequency differences

Genetic variants with extreme allele frequency differences (EAFD) may underlie some human health disparities across populations. To identify EAFD loci, we systematically analyzed and characterized 81 million genomic variants from 2504 unrelated individuals of 26 world populations (phase III of the 1000 Genomes Project). Our analyses revealed a total of 434 genes, 15 pathways, and 18 diseases and traits influenced by EAFD variants from five continental populations. They included known EAFD genes, such as LCT (lactose tolerance), SLC24A5 (skin pigmentation), and EDAR (hair morphology). We found many novel EAFD genes, including TBC1D2B (autophagy mediator), TRIM40 (gastrointestinal inflammatory regulator), KRT71, KRT75, KRT83, and KRTAP10-1 (hair and epithelial keratin synthesis), PIK3R3 (insulin receptor interaction), DARS (neurological disorders), and NACA2 (skin inflammatory response). Our results also showed four complex diseases significantly associated with EAFD loci, including asthma (adjusted enrichment P = 4 × 10^-8), type I diabetes (P = 6 × 10^-9), alcohol consumption (P = 0.0002), and attention deficit/hyperactivity disorder (P = 0.003). This study provides a comprehensive atlas of genes, pathways, and human diseases significantly influenced by EAFD variants.

Pathway Analysis Incorporating Protein-Protein Interaction Networks Identified Candidate Pathways for the Seven Common Diseases

Pathway analysis has become popular as a secondary analysis strategy for genome-wide association studies (GWAS). Most of the current pathway analysis methods aggregate signals from the main effects of single nucleotide polymorphisms (SNPs) in genes within a pathway without considering the effects of gene-gene interactions. However, gene-gene interactions can also have critical effects on complex diseases. Protein-protein interaction (PPI) networks have been used to define gene pairs for the gene-gene interaction tests. Incorporating the PPI information to define gene pairs for interaction tests within pathways can increase the power for pathway-based association tests. We propose a pathway association test, which aggregates the interaction signals in PPI networks within a pathway, for GWAS with case-control samples. Gene size is properly considered in the test so that genes do not contribute more to the test statistic simply due to their size. Simulation studies were performed to verify that the method is a valid test and can have more power than other pathway association tests in the presence of gene-gene interactions within a pathway under different scenarios. We applied the test to the Wellcome Trust Case Control Consortium GWAS datasets for seven common diseases. The most significant pathway is the chaperones modulate interferon signaling pathway for Crohn’s disease (p-value = 0.0003). The pathway modulates interferon gamma, which induces the JAK/STAT pathway that is involved in Crohn’s disease. Several other pathways that have functional implications for the seven diseases were also identified. The proposed test based on gene-gene interaction signals in PPI networks can be used as a complementary tool to the current existing pathway analysis methods focusing on main effects of genes. An efficient software implementing the method is freely available at http://puppi.sourceforge.net.

The search for the target antigens of multiple sclerosis, part 1: autoreactive CD4+ T lymphocytes as pathogenic effectors and therapeutic targets

Identification of the target antigens of pathogenic antibodies and T cells is of fundamental importance for understanding the pathogenesis of multiple sclerosis, and for the development of personalised treatments for the disease. Myelin-specific CD4+ T cells emerged long ago as a key player in animal models of multiple sclerosis. Taking a forward-translational approach, autoreactive CD4+ T cells have been studied extensively in patients with multiple sclerosis, and there is evidence, but as yet no direct proof, that autoreactive CD4+ T cells are a key player in the pathogenesis of the disorder. Several therapies that selectively target myelin-specific CD4+ T cells have been investigated in clinical trials up to phase 3. So far, however, none of these (mostly underpowered) therapeutic trials have provided definitive evidence of clinical efficacy. One major obstacle to personalised, highly selective immunotherapy is the absence of standardised and reliable assays to assess antigen-specific human T-cell responses. Such assays would be essential for stratification of patients with multiple sclerosis according to their individual target antigens.

Class II HLA interactions modulate genetic risk for multiple sclerosis

Association studies have greatly refined the understanding of how variation within the human leukocyte antigen (HLA) genes influences risk of multiple sclerosis. However, the extent to which major effects are modulated by interactions is poorly characterized. We analyzed high-density SNP data on 17,465 cases and 30,385 controls from 11 cohorts of European ancestry, in combination with imputation of classical HLA alleles, to build a high-resolution map of HLA genetic risk and assess the evidence for interactions involving classical HLA alleles. Among new and previously identified class II risk alleles (HLA-DRB1*15:01, HLA-DRB1*13:03, HLA-DRB1*03:01, HLA-DRB1*08:01 and HLA-DQB1*03:02) and class I protective alleles (HLA-A*02:01, HLA-B*44:02, HLA-B*38:01 and HLA-B*55:01), we find evidence for two interactions involving pairs of class II alleles: HLA-DQA1*01:01-HLA-DRB1*15:01 and HLA-DQB1*03:01-HLA-DQB1*03:02. We find no evidence for interactions between classical HLA alleles and non-HLA risk-associated variants and estimate a minimal effect of polygenic epistasis in modulating major risk alleles.

Multiple sclerosis associated genetic variants of CD226 impair regulatory T cell function

Recent association studies have linked numerous genetic variants with an increased risk for multiple sclerosis, although their functional relevance remains largely unknown. Here we investigated phenotypical and functional consequences of a genetic variant in the CD226 gene that, among other autoimmune diseases, predisposes to multiple sclerosis. Phenotypically, effector and regulatory CD4(+) memory T cells of healthy individuals carrying the predisposing CD226 genetic variant showed, in comparison to carriers of the protective variant, reduced surface expression of CD226 and an impaired induction of CD226 after stimulation. This haplotype-dependent reduction in CD226 expression on memory T cells was abrogated in patients with multiple sclerosis, as CD226 expression was comparable to healthy risk haplotype carriers irrespective of genetic variant. Functionally, FOXP3-positive regulatory T cells from healthy carriers of the genetic protective variant showed superior suppressive capacity, which was again abrogated in multiple sclerosis patients. Mimicking the phenotype of human CD226 genetic risk variant carriers, regulatory T cells derived from Cd226-deficient mice showed similarly reduced inhibitory activity, eventually resulting in an exacerbated disease course of experimental autoimmune encephalomyelitis, the animal model of multiple sclerosis. Therefore, by combining human and mouse analyses we show that CD226 exhibits an important role in the activation of regulatory T cells, with its genetically imposed dysregulation impairing regulatory T cell function.

Puzzling role of genetic risk factors in human longevity: "risk alleles" as pro-longevity variants

Complex diseases are major contributors to human mortality in old age. Paradoxically, many genetic variants that have been associated with increased risks of such diseases are found in genomes of long-lived people, and do not seem to compromise longevity. Here we argue that trade-off-like and conditional effects of genes can play central role in this phenomenon and in determining longevity. Such effects may occur as result of: (i) antagonistic influence of gene on the development of different health disorders; (ii) change in the effect of gene on vulnerability to death with age (especially, from “bad” to “good”); (iii) gene–gene interaction; and (iv) gene–environment interaction, among other factors. A review of current knowledge provides many examples of genetic factors that may increase the risk of one disease but reduce chances of developing another serious health condition, or improve survival from it. Factors that may increase risk of a major disease but attenuate manifestation of physical senescence are also discussed. Overall, available evidence suggests that the influence of a genetic variant on longevity may be negative, neutral or positive, depending on a delicate balance of the detrimental and beneficial effects of such variant on multiple health and aging related traits. This balance may change with age, internal and external environments, and depend on genetic surrounding. We conclude that trade-off-like and conditional genetic effects are very common and may result in situations when a disease “risk allele” can also be a pro-longevity variant, depending on context. We emphasize importance of considering such effects in both aging research and disease prevention.