- Research article
- Open Access
Characterization of the HLA-DRβ1 third hypervariable region amino acid sequence according to charge and parental inheritance in systemic sclerosis
Arthritis Research & Therapy volume 19, Article number: 46 (2017)
Specific HLA class II alleles are associated with systemic sclerosis (SSc) risk, clinical characteristics, and autoantibodies. HLA nomenclature initially developed with antibodies as typing reagents defining DRB1 allele groups. However, alleles from different DRB1 allele groups encode the same third hypervariable region (3rd HVR) sequence, the primary T-cell recognition site, and 3rd HVR charge differences can affect interactions with T cells. We considered 3rd HVR sequences (amino acids 67–74) irrespective of the allele group and analyzed parental inheritance considered according to the 3rd HVR charge, comparing SSc patients with controls.
In total, 306 families (121 SSc and 185 controls) were HLA genotyped and parental HLA-haplotype origin was determined. Analysis was conducted according to DRβ1 3rd HVR sequence, charge, and parental inheritance.
The distribution of 3rd HVR sequences differed in SSc patients versus controls (p = 0.007), primarily due to an increase of specific DRB1*11 alleles, in accord with previous observations. The 3rd HVR sequences were next analyzed according to charge and parental inheritance. Paternal transmission of DRB1 alleles encoding a +2 charge 3rd HVR was significantly reduced in SSc patients compared with maternal transmission (p = 0.0003, corrected for analysis of four charge categories p = 0.001). To a lesser extent, paternal transmission was increased when charge was 0 (p = 0.021, corrected for multiple comparisons p = 0.084). In contrast, paternal versus maternal inheritance was similar in controls.
SSc patients differed from controls when DRB1 alleles were categorized according to 3rd HVR sequences. Skewed parental inheritance was observed in SSc patients but not in controls when the DRβ1 3rd HVR was considered according to charge. These observations suggest that epigenetic modulation of HLA merits investigation in SSc.
HLA class II genes contribute the largest portion of genetic risk for many autoimmune diseases including systemic sclerosis (SSc) . The HLA DRB1 allele group DRB1*11 has been described consistently in association with SSc in Caucasian adults [1–3]. HLA nomenclature evolved utilizing serological typing reagents resulting in a “language of HLA”, with initial groupings based primarily on antibody recognition of the HLA molecule. However, the β1 chain of the HLA-DR molecule encoded by HLA-DRB1 is characterized by three regions of amino acid sequence hypervariability and the third hypervariable region (3rd HVR) is thought to be the most important site for T-cell recognition .
The 3rd HVR comprises amino acids 67–74 on the alpha helix of the HLA β1 chain. Allelic variation of the DRB1 DNA sequence can result in 3rd HVR amino acid sequences that differ for overall 3rd HVR charge and can impact the steric interaction between HLA molecules and T cells. The importance of charge in this region has been well demonstrated in rheumatoid arthritis, where different risk-associated DRβ1 molecules carry 3rd HVR sequences with a positive charge, whereas a negatively charged 3rd HVR sequence is associated with protection from the disease [5, 6].
In the current study we investigated the DRβ1 3rd HVR, charge, and parental transmission in SSc patients and healthy controls. We sought to examine inheritance according to parental origin because biased HLA transmission has been described in some other autoimmune diseases, including multiple sclerosis, and in some but not all studies of type 1 diabetes [7–9]. The primary purpose of this study was to test the hypothesis that SSc risk is modulated according to whether the mother or the father transmitted an HLA-DRB1 allele and is impacted by DRβ1 3rd HVR charge.
A total of 306 unrelated families were studied, 121 in which the proband had SSc and 185 healthy control families. SSc patients were recruited primarily from the Seattle, WA, USA area, with some patients from Alaska, Montana, Oregon and other states and some identified through the SSc family registry based in Houston, TX, USA. Healthy controls were recruited primarily from the Seattle, WA area. The median age at the time of sample collection for HLA genotyping for SSc patients was 43 years (range 18–62) and for controls was 31 years (range 4–73). SSc patients were female and controls were predominantly female (99 female, 22 male). All patients and controls were Caucasian. Study subjects were included providing maternal versus paternal transmission of HLA-DRB1 could be determined by HLA genotyping studies of the patient or control and family members. Eight SSc patients and one control were HLA genotyped but not included in the analysis because parental haplotype transmission could not be determined. All subjects provided written informed consent. The study was approved by the Institutional Review Board of Fred Hutchinson Cancer Research Center.
HLA genotyping was conducted from peripheral blood, buccal swabs, or mouthwash specimens. Genomic DNA was extracted from whole blood or peripheral blood mononuclear cells using the Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA), from mouthwash specimens using the High Pure PCR Template Preparation Kit (Roche Diagnostics, Indianapolis, IN, USA), or from buccal swabs using the BuccalAmp DNA Extraction Kit (Epicentre Biotechnologies, WI, USA). All subjects were genotyped for HLA-DRB1, as well as for DQA1 and DQB1 loci. DNA-based typing was conducted with the Luminex-based PCR-sequence-specific oligonucleotide probe technique (Luminex; One Lambda, Canoga Park, CA, USA) with alleles assigned using HLA Fusion™ 3.0 standard analysis software, or Dynal strip detection with sequence-specific oligonucleotide probes (Dynal RELITM SSO, UK) followed by identification of specific alleles by sequencing (Applied Biosystems, Foster city, CA, USA). Amino acid sequences were determined based on the specific alleles and the International Immunogenetics Information System for HLA (IGMT/HLA; www.ebi.ac.uk/ipd/imgt/hla/allele.html).
HLA-DRβ1 3rd HVR classification
All DRB1 alleles were first categorized according to the encoded 3rd HVR amino acid sequence, from positions 67 through 74 of the HLA-DRβ1 chain (Table 1). Among all possible DRβ1 3rd HVR sequences, 17 were present in at least one study subject with eight sequences present in at least 5% of either patients or controls. Next, 3rd HVR sequences were grouped according to their overall charge: +2, +1, 0, −1, or −2.
In our patient and control populations no subject had a DRB1 allele encoding a 3rd HVR sequence with a −1 charge. Therefore, our study compares four charge groups between SSc and controls.
After testing the normality of the distribution of continuous variables with the Shapiro–Wilk test, comparisons were made with the Wilcoxon rank-sum test as appropriate. Categorical variables were compared utilizing chi-square tests or, when indicated, Fisher exact tests. The primary outcome analyzed was maternal versus paternal inheritance of the DRβ1 3rd HVR sequence according to charge. In addition, we compared the overall distribution of 3rd HVR sequences and 3rd HVR charge in SSc patients and controls. Analysis was conducted when 5% or more of patients or controls were represented in a category and correction made for multiple comparisons by Bonferroni adjustment of the threshold for significance. Logistic regression was carried out adjusting for sex because some controls were male, but this revealed no differences from the unadjusted associations.
Because a parent of origin effect has not been studied previously in SSc, a precise estimate of statistical power is precluded. However, our available sample size of 121 cases and 185 controls is estimated to provide more than 90% power to detect a two-fold difference in parental origin from an expected 50/50 in controls to 25/75 in cases, assuming α = 0.05 and a two-tailed test.
Among SSc patients, 64% had diffuse SSc and 36% limited SSc (77 and 44 respectively). Antibodies to topoisomerase I (ATA) were positive in 34% of patients (37 of 109 tested) and to centromere (ACA) in 19% of patients (17 of 90 tested). The median age of SSc onset was 37 years, range 16–58.
HLA-DRB1 alleles were classified based on the amino acid sequence of the encoded 3rd HVR (amino acids 67–74) and the overall charge of each unique 3rd HVR sequence was determined (Table 1). Comparing SSc patients with healthy controls, the overall distribution of 3rd HVR sequences was significantly different (p = 0.007) (Table 2). The sequence “FLEDRRAA” (sequence category 5 in Table 1) was significantly enriched in SSc patients compared with controls (p = 0.004, p value corrected for multiple comparisons p c = 0.032), primarily due to enrichment of DRB1*11:04 among SSc patients, as reported previously . The sequence “FLEDRRAL” (sequence category 11 in Table 1) was also enriched in SSc patients (p = 0.006, p c = 0.048) (Table 2), due to some increase of DRB1*08:01 and DRB1*08:02 in SSc patients. These observations are concordant with a previous report of an increase in the motif “FLEDR” in a French population . In contrast to the French study, however, we did not see enrichment of this sequence encoded by DRB5, which instead was decreased in our population (data not shown).
Overall, the 3rd HVR charge on SSc haplotypes was similar to that of controls; of the 242 SSc and 370 control haplotypes respectively, 11% and 12% had a +2 charge, 26% and 31% had a +1 charge, 50% and 44% had a charge of 0, and 13% and 13% had a −2 charge. As noted previously, no 3rd HVR sequence represented in our populations carried a −1 charge. Thirty-six percent of SSc patients and 34% of controls carried the same charge on both parental haplotypes. When considering both haplotypes together, 19%, 32%, and 49% of the SSc patients had overall negative, neutral, or positive charge respectively, compared with 20%, 22%, and 58% of healthy controls (Fig. 1).
Analysis was conducted considering maternal or paternal DRB1 allele inheritance according to the charge of the DRβ1 3rd HVR. The overall distribution of inheritance indicated significant skewing among SSc patients (p < 0.001). Among SSc patients the father transmitted a 3rd HVR sequence with a +2 charge significantly less often than the mother (p = 0.0003). Almost 90% of the time, transmission was maternal rather than paternal (Fig. 1). This skewed parental inheritance remained significant (p = 0.001) when corrected for comparing four different categories of charge. The skewed inheritance was observed among patients with both diffuse and limited SSc and was not different for those with younger versus older age of SSc onset. Skewed inheritance for overall distribution, and for each charge category, was not observed in the control population in which inheritance was not different from 50–50, as expected. To a lesser extent, 3rd HVR sequences with a 0 charge were inherited more often from the father than the mother among SSc patients (p = 0.021, p c = 0.084). Parental transmissions of 3rd HVR sequences with +1 or −2 charge were not significantly skewed (Fig. 1).
While our study was designed to test a hypothesis regarding charge of the DRβ1 3rd HVR sequence and parental inheritance, linkage disequilibrium of DRB1 is very strong with DQA1 and DQB1, and secondary analysis was conducted for these loci. DQβ1 has a similarly located hypervariable region from amino acids 70 through 77. DQB1 alleles represented in our study populations had an overall net charge of +3, +1, 0, −1, or −2 (Table 3). Not surprisingly, the +3 charge for DQβ1 was also significantly decreased for paternal versus maternal transmission, because DRB1*03:01 (+2 DRβ1 3rd HVR) haplotypes in our populations with few exceptions carried DQB1*02:01, which has a +3 charge. Patterns of linkage disequilibrium did not permit clear separation of the role of DQB1 from DRB1. However, the skewing for DQβ1 with a +3 charge was less pronounced than that described for DRβ1; eight were transmitted paternally versus 31 transmitted maternally. Additionally, DRB1*07:01 haplotypes can carry either DQB1*02:02 (+3) or DQB1*03:03 (0), and there was no suggestion of skewed paternal inheritance for these two haplotypes or an increase in DRB1*07:01–DQB1*02:02 haplotypes in the SSc population overall (Additional file 1: Table S1).
Polymorphisms on DQα1 are distributed somewhat differently from DRβ1 and DQβ1, but considered for amino acids 47–56 the net charge was +4, +2, or +1 for alleles represented in our study population and did not differ due to parental inheritance (Table 4). This result was also consistent with linkage disequilibrium in the HLA class II region, but results for DQA1 could be distinguished from DRB1 because multiple different DRB1 alleles are on haplotypes with the similar DQA1 alleles and some DQA1 alleles other than DQA1*05:01 also encode a +2 charge.
As expected because of linkage disequilibrium, the DRB1*03:01–DQA1*05:01–DQB1*02:01 haplotype was similarly skewed with three inherited paternally versus 21 inherited maternally (p c = 0.001), although again skewing was somewhat less pronounced than when considered for DRβ1. (Haplotypes that were common in our study population are presented in Additional file 1: Table S1.)
In the current study we investigated patients with SSc and healthy controls, analyzing the DRβ1 3rd HVR according to charge irrespective of specific DRB1 alleles and considered parental inheritance. The importance of the DRβ1 3rd HVR and of charge in this region is well established in another autoimmune disease, rheumatoid arthritis, for which underlying HLA disease susceptibility is believed to be due to amino acid motifs of the DRβ1 3rd HVR that carry a similar charge [5, 6]. Our primary interest was to ask whether SSc patients differed from healthy controls if considered for parental origin of the DRβ1 3rd HVR categorized according to charge. In contrast to healthy individuals among whom parental inheritance was 50–50 as expected, we found significant skewing of parental inheritance in SSc patients. The most striking finding was that the father transmitted a DRB1 allele encoding a +2 charge 3rd HVR significantly less often than the mother in SSc patients. With respect to DRB1 allele frequencies, our study population was similar to prior reports of similar populations [1–3], with the patient population primarily enriched for DRB1*11:04 compared with controls.
Skewed parental inheritance of HLA-DRB1 alleles in SSc strongly suggests a parent-of-origin effect. Parent-of-origin effects are well described for a number of genes and are especially well recognized as important in early development . Skewing from the expected random (50–50) inactivation of maternal versus paternal X-chromosomes has been reported previously in women with SSc, as well as in women with rheumatoid arthritis and some other diseases [12, 13]. A few studies have evaluated parental inheritance of HLA in other autoimmune diseases. In type 1 diabetes, parent-of-origin skewing in HLA inheritance has been described in some, but not all, studies [8, 9]. In patients with multiple sclerosis, HLA transmission was distorted by the parent of origin and by gender of the affected offspring . A differential methylation signal with a peak at HLA-DRB1 was observed in CD4+ T cells of multiple sclerosis patients compared with controls in another study, suggesting a potential explanation for this observation .
A number of phenomena can underlie a parent-of-origin effect. The best characterized parent-of-origin effect is that associated with genomic imprinting . Embryonic lethality is one explanation for a parent-of-origin effect and, although skewing was not observed in our healthy population, it is possible that an extreme phenotype could result in loss of embryos destined to develop SSc in later life. The rate of decay for small molecules or epigenetic marks in sperm and ova affecting survival to gametogenesis has also been hypothesized to play a role in distorted haplotype transmission . Additionally, the molecules CCCTC-binding factor (CTCF) and Class II transactivator (CIITA), two units of a transcription complex involved in HLA class II gene transcription, are thought to be subject to epigenetic modifications [15, 16]. The current results are of further interest in light of a recent study that implicated fetal programming in SSc  as well as increasing evidence for a role of epigenetic regulation in SSc pathogenesis .
There are a number of limitations to our study. SSc patients in our study were all female and it would be of additional interest to know whether males and/or children with SSc differ. Some of our controls were male and younger than the patient population. However, neither would be anticipated to affect the overall study observations because controls had the expected random (50–50) distribution regardless of gender or age. SSc is also a rare disorder offsetting any potential bias due to development of SSc later in life by a younger control. Another limitation is that all of our study subjects were Caucasian because too few patients from other racial/ethnic backgrounds were available for analysis. Also, the number of SSc patients was modest and a larger study would be needed to evaluate potential differences according to clinical and autoantibody characteristics. Another limitation is that alleles carrying a +2 charge other than DRB1*03:01 were uncommon in our population so the skewed inheritance observed could be specific to genes on the DRB1*03:01 haplotype, including DQB1*0201, rather than the DRβ1 3rd HVR +2 charge. It should also be added that some amino acids are subject to posttranslational modification that results in a change of charge  and, while this has largely been examined on autoantigens, HLA molecules are themselves also presented as self-peptides by other HLA molecules .
Skewed parental HLA inheritance was observed among SSc patients compared with healthy controls. Paternal transmission was significantly reduced compared with maternal transmission when the HLA-DRB1 allele encoded a 3rd HVR with a +2 charge and to a lesser extent increased when charge was 0, in SSc patients but not in controls. To our knowledge, skewed parental inheritance of HLA-DRB1 alleles has not been reported previously in SSc. Provided our results are replicated, future investigations into the underlying mechanism(s) for a parent-of-origin effect of HLA inheritance in SSc may lead to new insight into the role of HLA molecules in SSc pathogenesis.
Human leukocyte antigen (or allele)
Mayes MD, Bossini-Castillo L, Gorlova O, Martin JE, Zhou X, Chen WV, et al. Immunochip analysis identifies multiple susceptibility loci for systemic sclerosis. Am J Hum Genet. 2014;94:47–61.
Arnett FC, Gourh P, Shete S, Ahn CW, Honey RE, Agarwal SK, et al. Major histocompatibility complex (MHC) class II alleles, haplotypes and epitopes which confer susceptibility or protection in systemic sclerosis: analysis in 1300 Caucasian, African-American and Hispanic cases and 1000 controls. Ann Rheum Dis. 2010;69:822–7.
Loubière LS, Lambert NC, Madeleine MM, Porter AJ, Mullarkey ME, Pang JM, et al. HLA allelic variants encoding DR11 in diffuse and limited systemic sclerosis in Caucasian women. Rheumatol Oxf Engl. 2005;44:318–22.
Reinsmoen NL, Bach FH. Structural Model for T-cell recognition of HLA class II-associated alloepitopes. Human Immunol. 1990;27(1):51–72.
Gregersen PK, Silver J, Winchester RJ. The shared epitope hypothesis. An approach to understanding the molecular genetics of susceptibility to rheumatoid arthritis. Arthritis Rheum. 1987;30(11):1205–13.
Reviron D, Perdriger A, Toussirot E, Wendling D, Balandraud N, Guid S, Semana G, Tiberghien P, Mercier P, Roudier J. Influence of shared epitope-negative HLA-DRB1 alleles on genetic susceptibility to rheumatoid arthritis. Arthritis Rheum. 2001;44(3):535–40.
Chao MJ, Herrera BM, Ramagopalan SV, Deluca G, Handunetthi L, Orton SM, Lincoln MR, Sadovnick AD, Ebers GC. Parent-of-origin effects at the major histocompatibility complex in multiple sclerosis. Hum Mol Gen. 2010;19(18):3679–89.
Sasaki T, Nemoto M, Yamasaki K, Tajima N. Preferential transmission of maternal allele with DQA1*0301-DQB1*0302 haplotype to affected offspring in families with type 1 diabetes. J Hum Genet. 1999;44:318–22.
Bronson PG, Ramsay PP, Thomson G, Barcellos LF, and the Type 1 Diabetes Genetics Consortium. Analysis of maternal-offspring HLA compatibility, parent-of-origin and non-inherited maternal effects for the classical HLA loci in type 1 diabetes. Diabetes Obs Metabol. 2009;11 Suppl 1:74–83.
Azzouz D, Rak JM, Fajardy I, Allanore Y, Tiev KP, et al. Comparing HLA shared epitopes in French Caucasian patients with Scleroderma. PLoS One. 2012;7(5):e36870.
Lawson HA, Cheverud JM, Wolf JB. Genomic imprinting and parent-of-origin effects on complex traits. Nat Rev. 2013;14:609–17.
Uz E, Loubiere LS, Gadi VK, Ozbalkan Z, Stewart J, Nelson JL, Ozcelik T. Skewed X-chromosome inactivation in scleroderma. Clin Rev Allergy Immunol. 2008;34(3):352–5.
Chabchoub G, Uz E, Maalej A, Mustafa CA, Rebai A, Mnif M, Bahloul Z, Farid NR, Ozcelik T, Ayadi H. Analysis of skewed X-chromosome inactivation in females with rheumatoid arthritis and autoimmune thyroid disease. Arthritis Res Therapy. 2009;11(4):R106.
Mc G, Benton M, Lea RA, Boyle M, Tajouri L, Macartney-Coxson D, Scott RJ, Lechner-Scott J. Methylation differences at the HLA-DRB1 locus in CD4+ T-cells are associated with multiple sclerosis. Mult Scler J. 2014;20(8):1033–41.
Phillips JE, Corces VG. CTCF: master weaver of the genome. Cell. 2011;137(7):1194–211.
Wright KL, Ting JPY. Epigenetic regulation of MHC-II and CIITA genes. Trends Immunol. 2006;27(9):405–12.
Donzelli G, Carnesecchi G, Amador C, di Tommaso M, Filippi L, Caporali R, et al. Fetal programming and systemic sclerosis. Am J Obstet Gynecol. 2015;213(6):839.e1–8–8.
Altorok N, Almeshal N, Wang Y, Kahaleh B. Epigenetics, the holy grail in the pathogenesis of systemic sclerosis. Rheumatology. 2015;54(10):1759–70.
Zavala-Cerna M, Martinez-Garcia E, Torres-Bugarin O, Rubio-Jurado B, Riebeling C, Nava A. The clinical significance of posttranslational modification of autoantigens. Clin Rev Allergy Immunol. 2014;47:73–90.
Mohme M, Hotz C, Stevanovic S, Binder T, Lee JH, Okoniewski M, et al. HLA-DR15-derived self-peptides are involved in increased autologous T cell proliferation in multiple sclerosis. Brain. 2013;136:1783–98.
The authors are grateful to The Scleroderma Family Registry and to all of the patients and family members as well as the healthy controls and their family members for their participation in this study and to Judy Allen MPH and Samantha Bell MPH for study coordination. They thank Alex Forsyth for additional review of HLA alleles and amino acid sequences.
This work was supported by NIH grant RO1 AI-41721 and a grant from The Scleroderma Foundation.
Availability of data and materials
The datasets used and analyzed in the current study are not part of a public database but are available from the corresponding author upon reasonable request.
CAG and CTL conducted experiments. CAG, HSG, and JLN analyzed and interpreted the data. MDM and DEF contributed patients and clinical assessments and contributed to the summary of results. CAG and JLN primarily wrote the manuscript. JLN conceived the study. All authors approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All subjects provided written informed consent. The study was approved by the Institutional Review Board of Fred Hutchinson Cancer Research Center.
About this article
Cite this article
Gentil, C.A., Gammill, H.S., Luu, C.T. et al. Characterization of the HLA-DRβ1 third hypervariable region amino acid sequence according to charge and parental inheritance in systemic sclerosis. Arthritis Res Ther 19, 46 (2017) doi:10.1186/s13075-017-1253-9
- Systemic sclerosis
- Human leukocyte antigen
- Skewed parental inheritance