Population genetics and external proficiency testing for HLA disease associations

Numerous associations of HLA variants with susceptibility to diseases, namely, those with an immunopathological component, have been described to date. The strongest HLA associations were incorporated into the standard algorithms for the diagnostics. Disease-associated HLA variants are routinely detected by various techniques including DNA-based assays. For the identification of HLA markers or their combinations with the highest diagnostic value and those with frequent clinical indications (e.g., HLA-B*27, -B*57:01, -DQ2/-DQ8, -DQB1*06:02), diagnostic tests that focus on a single or limited number of specific HLA antigens/alleles, have already been developed; the use of complete typing for particular HLA loci is a relevant alternative. Importantly, external proficiency testing (EPT) became an integral part of good laboratory practice for HLA disease associations in accredited laboratories and not only supports correct "technical" identification of the associated HLA variants, but also adequate interpretation of the results to the clinicians. In the present article selected aspects of EPT for HLA disease associations related to population genetics are reviewed and discussed with the emphasis on the optimal level of HLA typing resolution, population-based differences in disease associated HLA alleles within the allelic group, distribution and linkage disequilibrium of HLA alleles in particular populations and interpretation of the presence of less common HLA variants/haplotypes. In conclusion, the laboratories that perform and interpret the tests to the clinicians, producers of the certified diagnostics and EPT providers should consider, among others, the genetic characteristics of the populations in order to optimise the diagnostic value of the tests for disease-associated HLA variants.

The novel HLA-A*29:172 allele identified in a patient indicated for hematopoietic stem cell transplantation

HLA-A*29:172 allele differs from HLA-A*29:01:01:01 by one missense single C/G nucleotide exchange in codon 77.

Cellular senescence marker p16INK4a and NFKB1 gene polymorphisms in lower gastro-intestinal acute graft versus host disease

BACKGROUND

Lower gastrointestinal (GI) graft versus host disease (GVHD) represents a severe complication in allogeneic hematopoietic stem cell transplant (HSCT) recipients with high rates of transplant-related mortality. Deregulated innate immunity reactions are the features of its pathogenesis. Cellular senescence has been considered a program of the innate immunity. We focused on lower GI GVHD from the perspective of cellular senescence.

OBJECTIVE

We analyzed the impact of p16^INK4a expression, a hallmark of cellular senescence, in intestinal biopsies of patients with lower GI GVHD symptoms and NFKB1 gene polymorphisms (rs3774937 C/T and rs3774959 A/G) on HSCT outcome.

STUDY DESIGN

Fifty-two single-center patients who presented with symptoms of lower GI GVHD were analyzed in a retrospective manner. Two SNPs located in the NFKB1 gene regions (rs3774937 C/T and rs3774959 A/G) were genotyped from the peripheral blood samples collected before the start of the conditioning. All patients underwent proctosigmoidoscopy with biopsy of the mucosa. The expression of p16^INK4a was analyzed in normal intestinal crypts and stroma.

RESULTS

Fifty-two patients (50% male) received HSCT for hematological diseases (acute leukemias in 67%) and developed lower GI symptoms. Patients with p16^INK4a expression in the intestinal stroma were in lower risk of developing histological grade 3-4 aGVHD (RR 0.18 [95% CI 0.05-0.65]; p = 0.009). The multivariate linear regression confirmed the independent effect of p16^INK4a expression on time of the lower GI aGVHD symptoms onset (Coef. 38.9 [95% CI 12.7-65.1]; p = 0.005). The NFKB1 rs3774937 CC and TT/TC genotype were present in 40 and 80% of patients with p16^INK4a expression, respectively (p = 0.04). The rs3774959 AA and GG/AG genotype were present among 43 and 82% of patients with p16^INK4a expression, respectively (p = 0.02). Expression of p16^INK4a was associated with no clinical variable but NFKB1 genotype.

CONCLUSIONS

Our results address possible new mechanisms that may lead to better understanding of HSCT-related immune complications. Cellular senescence may bring novel approaches in GVHD diagnostics and therapy.

Association of multispecific red blood cell alloimmunization with HLA-Class II variants is related to Rh phenotypes

AIM

It remains unclear, why only some patients form alloantibodies against foreign RBC antigens. Transfusion of red blood cell (RBC) products and pregnancy are the most relevant causes of immunization against RBC alloantigens. Here we investigated the relationship between RBC alloantibodies, Rh phenotype, and HLA phenotype among patients with multiple RBC alloantibodiesMETHODS: In a group of 124 multi-responders ‒ including both pregnant women and transplant recipients ‒ we analysed the distribution of HLA-Class II variants in subgroups of multi-responders to RBC alloantigens according to their Rh status.

RESULTS

As expected, the RhD-negative phenotype was overrepresented in our alloimmunized group (49.2 %) compared to in the general population. Importantly, HLA-DRB1*15 carriers were significantly overrepresented among D-negative multi-responders compared to D-positive multi-responders (Pc = 0.045). Furthermore, the linked HLA-DRB1*13, HLA-DQB1*06, and HLA-DQA1*01 variants were more frequent in individuals with the DCCee phenotype than in other RhD-positive phenotypes.

CONCLUSION

Our present findings showed that RBC multispecific alloimmunization was associated with particular HLA-Class II variants based on Rh status (Tab. 3, Ref. 22).

Genetic control of CCL24, POR, and IL23R contributes to the pathogenesis of sarcoidosis

Sarcoidosis is a genetically complex systemic inflammatory disease that affects multiple organs. We present a GWAS of a Japanese cohort (700 sarcoidosis cases and 886 controls) with replication in independent samples from Japan (931 cases and 1,042 controls) and the Czech Republic (265 cases and 264 controls). We identified three loci outside the HLA complex, CCL24, STYXL1-SRRM3, and C1orf141-IL23R, which showed genome-wide significant associations (P < 5.0 × 10⁻⁸) with sarcoidosis; CCL24 and STYXL1-SRRM3 were novel. The disease-risk alleles in CCL24 and IL23R were associated with reduced CCL24 and IL23R expression, respectively. The disease-risk allele in STYXL1-SRRM3 was associated with elevated POR expression. These results suggest that genetic control of CCL24, POR, and IL23R expression contribute to the pathogenesis of sarcoidosis. We speculate that the CCL24 risk allele might be involved in a polarized Th1 response in sarcoidosis, and that POR and IL23R risk alleles may lead to diminished host defense against sarcoidosis pathogens.

Akira Meguro et al. report a genome-wide association study for sarcoidosis—a systemic inflammatory disease—in the Japanese population. They identify 3 non-HLA loci with genome-wide significance, 2 of which have not been previously associated with sarcoidosis in any population.

NFKB1 gene single-nucleotide polymorphisms: implications for graft-versus-host disease in allogeneic hematopoietic stem cell transplantation

Graft-versus-host disease (GVHD) represents a significant cause of mortality after allogeneic hematopoietic stem cell transplantation (HSCT). NF-kB system is a master regulator of innate immunity responses. It controls the expression of various cytokines and chemokines many of which are involved in GVHD pathogenesis. Chemo(radio) therapy administered during conditioning induces DNA damage and activates DNA damage response (DDR) signaling resulting in irreversible cell cycle arrest - cellular senescence which has been described to be associated with robust pro-inflammatory secretion mostly controlled by NF-kB. The NFKB1 gene encodes the DNA-binding subunit of the NF-kB complex. Using the candidate gene approach, we analyzed possible association of two single-nucleotide polymorphisms (SNPs) rs3774937 C/T and rs3774959 A/G of the NFKB1 gene with GVHD and transplant-related mortality (TRM) occurrence in 109 recipients allografted from HLA-identical donor. Both SNPs in recipients were found to be strongly associated with acute GVHD. Nevertheless, no significant association with chronic GVHD and TRM was found. Presented pilot results contribute to pre-clinical observations and suggest that NF-kB may be an important regulator of HSCT-related inflammatory reactions such as acute GVHD. Novel pathogenic mechanisms of GVHD may arise from perspectives of DDR and cellular senescence where NF-kB plays an essential role.

Modulatory Effect of the Euro-Lupus Low-Dose Intravenous Cyclophosphamide Regimen on Circulating Immune Cells in Systemic Lupus Erythematosus

A Euro-Lupus regimen of low-dose intravenous cyclophosphamide (CFA) is commonly used to treat severe organ manifestations of systemic lupus erythematosus (SLE), particularly lupus nephritis (LN). There are no data on the distributions and dynamics of immune cell populations in patients with various treatment outcomes. The circulating immune cells of 11 female SLE patients were assessed before and after Euro-Lupus regimen (cumulative dose of 3000 mg CFA) by flow cytometry together with those of 16 healthy women. A subanalysis was performed in LN patients who achieved complete remission (CR; n = 3), partial remission (PR; n = 4), and no response (NR; n = 2). In SLE, the Euro-Lupus regimen decreased the percentage and absolute count of B cells; increased the percentage of CD8⁺ T cells, T regulatory cells, neutrophils, and monocyte subsets; and activated T and NK cells compared to healthy controls (P < 0.050). Patients with LN achieving CR had significantly lower proportions of CD27⁺ B memory cells compared to poor responders (PR/NR, P = 0.035). The post-treatment percentages and absolute numbers of B cells, T cells, NK cells, monocytes, and neutrophils showed high inter-individual variability with no association with treatment outcome. Our pilot study revealed the dynamics of changes in immune cell populations in SLE patients during a Euro-Lupus regimen, mainly the lowering of B cells. In LN patients who achieved CR, a lower proportion of CD27⁺ B memory cells was evident compared to poor responders (PR/NR). Further studies on usefulness of monitoring immune cells for treatment response prediction on larger cohorts are needed.

Cross-Disease Innate Gene Signature: Emerging Diversity and Abundance in RA Comparing to SLE and SSc

Overactivation of the innate immune system together with the impaired downstream pathway of type I interferon-responding genes is a hallmark of rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), and systemic sclerosis (SSc). To date, limited data on the cross-disease innate gene signature exists among those diseases. We compared therefore an innate gene signature of Toll-like receptors (TLRs), seven key members of the interleukin (IL)1/IL1R family, and CXCL8/IL8 in peripheral blood mononuclear cells from well-defined patients with active stages of RA (n = 36, DAS28 ≥ 3.2), SLE (n = 28, SLEDAI > 6), and SSc (n = 22, revised EUSTAR index > 2.25). Emerging diversity and abundance of the innate signature in RA patients were detected: RA was characterized by the upregulation of TLR3, TLR5, IL1RAP/IL1R3, IL18R1, and SIGIRR/IL1R8 when compared to SSc (P corr < 0.02) and of TLR2, TLR5, and SIGIRR/IL1R8 when compared to SLE (P corr < 0.02). Applying the association rule analysis, six rules (combinations and expression of genes describing disease) were identified for RA (most frequently included high TLR3 and/or IL1RAP/IL1R3) and three rules for SLE (low IL1RN and IL18R1) and SSc (low TLR5 and IL18R1). This first cross-disease study identified emerging heterogeneity in the innate signature of RA patients with many upregulated innate genes compared to that of SLE and SSc.

TP53 rs1042522 and rs8064946 variants in myocardial infarction

OBJECTIVE

This study investigated the hypothesis that the single nucleotide polymorphisms (SNPs) of TP53 gene are related to a risk of myocardial infarction.

METHODS

The coding SNP at codon 72 (rs1042522) and non-coding rs8064946 SNP were genotyped by polymerase chain reaction with sequence specific primers in 205 Czech patients with myocardial infarction and 148 Czech control subjects.

RESULTS

The distribution of both SNPs was in agreement with the Hardy-Weinberg equilibrium and was similar to other European populations. Our power analysis showed 96 % of probability to detect an odd ratio equal to 2. Neither rs1042522 nor rs8064946 were associated with the risk of myocardial infarction. The haplotypes combined of rs1042522 and rs8064946 were not associated with myocardial infarction in the present study.

CONCLUSION

The TP53 SNPs are not strongly associated with genetic predisposition to myocardial infarction (Tab. 3, Fig. 3, Ref. 23).

Serum protein pattern associated with organ damage and lupus nephritis in systemic lupus erythematosus revealed by PEA immunoassay

BACKGROUND

Systemic lupus erythematosus (SLE) is a remarkably heterogeneous autoimmune disease. Despite tremendous efforts, our knowledge of serum protein patterns in severe SLE phenotypes is still limited. We investigated the serum protein pattern of SLE, with special emphasis on irreversible organ damage and active lupus nephritis (LN) as assessed by renal Systemic Lupus Erythematosus Disease Activity Index.

METHODS

We used proximity extension immunoassay (PEA, Proseek Multiplex, Olink) to assess the serum levels of ninety-two inflammation-related proteins in Czech patients with SLE (n = 75) and age-matched healthy control subjects (n = 23). Subgroup analysis was carried out on the basis of organ damage (with/without, 42/33) and biopsy-proven LN (with/without, 27/48; active LN, n = 13; inactive LN, n = 14).

RESULTS

Of thirty deregulated proteins between SLE and the healthy controls (P_corr  < 0.05), the top upregulated proteins in SLE were sirtuin 2, interleukin 18 (IL18), and caspase 8 (P_corr  < 0.0006). Of these, sirtuin 2 and caspase 8 had not yet been reported with SLE. Elevated levels of IL8, CCL2/MCP1, CCL11, and MMP10 (P_corr  < 0.05) were detected in patients with organ damage for which the serum levels of CCL11 and MMP10 were particularly informative in organ damage prediction. Comparing patients based on LN, elevated levels of CSF1, sIL15RA, sCD40, sCX3CL1, caspase 8, sIL18R1, bNGF, and GDNF (P_corr  < 0.05) were detected in active LN. Except GDNF, all LN-associated markers showed usefulness in prediction of active renal disease.

CONCLUSIONS

This highly sensitive PEA analysis identified the serum pattern of SLE, organ damage, and active LN, with many novel candidate proteins detected. Their exact role and suitability as biomarkers in SLE deserve further investigation.

SNP Variants in Major Histocompatibility Complex Are Associated with Sarcoidosis Susceptibility-A Joint Analysis in Four European Populations

Sarcoidosis is a multiorgan inflammatory disorder with heritability estimates up to 66%. Previous studies have shown the major histocompatibility complex (MHC) region to be associated with sarcoidosis, suggesting a functional role for antigen-presenting molecules and immune mediators in the disease pathogenesis. To detect variants predisposing to sarcoidosis and to identify genetic differences between patient subgroups, we studied four genes in the MHC Class III region (LTA, TNF, AGER, BTNL2) and HLA-DRA with tag-SNPs and their relation to HLA-DRB1 alleles. We present results from a joint analysis of four study populations (Finnish, Swedish, Dutch, and Czech). Patients with sarcoidosis (n = 805) were further subdivided based on the disease activity and the presence of Löfgren’s syndrome. In a joint analysis, seven SNPs were associated with non-Löfgren sarcoidosis (NL; the strongest association with rs3177928, P = 1.79E−07, OR = 1.9) and eight with Löfgren’s syndrome [Löfgren syndrome (LS); the strongest association with rs3129843, P = 3.44E−12, OR = 3.4] when compared with healthy controls (n = 870). Five SNPs were associated with sarcoidosis disease course (the strongest association with rs3177928, P = 0.003, OR = 1.9). The high linkage disequilibrium (LD) between SNPs and an HLA-DRB1 challenged the result interpretation. When the SNPs and HLA-DRB1 alleles were analyzed together, independent association was observed for four SNPs in the HLA-DRA/BTNL2 region: rs3135365 (NL; P = 0.015), rs3177928 (NL; P < 0.001), rs6937545 (LS; P = 0.012), and rs5007259 (disease activity; P = 0.002). These SNPs act as expression quantitative trait loci (eQTL) for HLA-DRB1 and/or HLA-DRB5. In conclusion, we found novel SNPs in BTNL2 and HLA-DRA regions associating with sarcoidosis. Our finding further establishes that polymorphisms in the HLA-DRA and BTNL2 have a role in sarcoidosis susceptibility. This multi-population study demonstrates that at least a part of these associations are HLA-DRB1 independent (e.g., not due to LD) and shared across ancestral origins. The variants that were independent of HLA-DRB1 associations acted as eQTL for HLA-DRB1 and/or -DRB5, suggesting a role in regulating gene expression.

Association of HLA-DRB1 and HLA-DQB1 with red-blood-cell alloimmunization in the Czech population

BACKGROUND AND OBJECTIVES

Alloimmune antibodies against red-blood-cell (RBC) antigens induced in susceptible individuals (responders) by transfusion, pregnancy or transplantation may have serious clinical consequences. The aim of this study was to investigate association of alloimmunization against selected RBC antigens with HLA-Class II.

MATERIALS AND METHODS

A total of 230 responders (106 monoresponders and 124 multiresponders) were enrolled into the study. HLA-DRB1 and HLA-DQB1 variants were determined by PCR-SSO and their frequencies compared between the patients (patient subgroups) and 375 ethnically and regionally matched controls.

RESULTS

Development of multiple RBC antibodies was associated with HLA-DRB1*15 and HLA-DQB1*06 allelic groups in the patients, with the relationship being particularly apparent in those with anti-C+D antibodies. Furthermore, DRB1*13 and DQB1*06 were more frequent in multiresponders with anti-E+c antibodies and DRB1*03 and DQB1*02 in those with anti-E+Cw.

CONCLUSION

For the first time, we confirmed the association of HLA-DRB1*15 with RBC antibody multiresponder status and found HLA-Class II associations for three frequent RBC antibody combinations. Our data support the concept that HLA restriction plays an important role in the response to RBC alloantigens.

A novel HLA-B allele, HLA-B*35:279, identified by sequencing-based typing in a Czech patient

The identification of a novel HLA-B*35:279 allele in a Czech patient is described. This allele is identical to the B*35:03:01 variant except the G/A nucleotide exchange at position 652 of the HLA-B gene that corresponds to the amino acid substitution from valine to isoleucine in alpha 3 domain of the HLA-B antigen.

Identification of Immune-Relevant Factors Conferring Sarcoidosis Genetic Risk

RATIONALE

Genetic variation plays a significant role in the etiology of sarcoidosis. However, only a small fraction of its heritability has been explained so far.

OBJECTIVES

To define further genetic risk loci for sarcoidosis, we used the Immunochip for a candidate gene association study of immune-associated loci.

METHODS

Altogether the study population comprised over 19,000 individuals. In a two-stage design, 1,726 German sarcoidosis cases and 5,482 control subjects were genotyped for 128,705 single-nucleotide polymorphisms using the Illumina Immunochip for the screening step. The remaining 3,955 cases, 7,514 control subjects, and 684 parents of affected offspring were used for validation and replication of 44 candidate and two established risk single-nucleotide polymorphisms.

MEASUREMENTS AND MAIN RESULTS

Four novel susceptibility loci were identified with genome-wide significance in the European case-control populations, located on chromosomes 12q24.12 (rs653178; ATXN2/SH2B3), 5q33.3 (rs4921492; IL12B), 4q24 (rs223498; MANBA/NFKB1), and 2q33.2 (rs6748088; FAM117B). We further defined three independent association signals in the HLA region with genome-wide significance, peaking in the BTNL2 promoter region (rs5007259), at HLA-B (rs4143332/HLA-B*0801) and at HLA-DPB1 (rs9277542), and found another novel independent signal near IL23R (rs12069782) on chromosome 1p31.3.

CONCLUSIONS

Functional predictions and protein network analyses suggest a prominent role of the drug-targetable IL23/Th17 signaling pathway in the genetic etiology of sarcoidosis. Our findings reveal a substantial genetic overlap of sarcoidosis with diverse immune-mediated inflammatory disorders, which could be of relevance for the clinical application of modern therapeutics.

ATM gene single nucleotide polymorphisms predict regimen-related gastrointestinal toxicity in patients allografted after reduced conditioning

Polymorphisms of genes involved in innate and adaptive immunity have become an object of major interest in regard to hematopoietic stem cell transplantation (HSCT) complications. Regimen-related gastrointestinal toxicity (RR-GIT) is the dominant complication during the pre-engraftment period and has been linked to increased risk of graft-versus-host disease (GVHD) development. According to our hypothesis, functional variants of genes participating in DNA damage response (DDR) may have an impact on the extent of tissue damage caused by the conditioning regimen. In our single-center study, we analyzed 62 patients who underwent HSCT from HLA-identical donors after reduced conditioning. The patients were genotyped for 5 single nucleotide polymorphisms (SNPs, rs4585 T/G, rs189037 A/G, rs227092 T/G, rs228590 C/T, and rs664677 T/C) of the ATM gene-the essential member of the DDR pathways, using allele-specific matrix-assisted laser desorption/ionization, time-of-flight (MALDI-TOF) mass spectrometry assay. Because of almost absolute linkage disequilibrium observed among all 5 SNPs, association of 2 major ATM haplotypes (ATM1/ATM2) with RR-GIT and acute GVHD (aGVHD) was analyzed. Importantly, the univariate and multivariate analysis showed that patients homozygous for ATM2 haplotype (rs4585*T, rs189037*A, rs227092*T, rs228590*C, and rs664677*T) are more likely to suffer from high-grade RR-GIT than ATM1 homozygous patients. The association with aGVHD was not significant. To our knowledge, this is the first report showing the ATM gene variability in relation to RR-GIT in the allogeneic HSCT setting.

Sequence variation in promoter regions of genes for CC chemokine ligands (CCL)19 and 21 in Czech patients with myocardial infarction

Recruitment of inflammatory cells to the arterial wall is an important pathogenic mechanism of atherosclerosis and coronary artery disease (CAD). Functional variability in the genes encoding for chemokines that promote infiltration of atherosclerotic plaques by macrophages and lymphocytes may therefore contribute to the genetic susceptibility to CAD. We, therefore, investigated the association between myocardial infarction (MI) and polymorphisms in the promoter regions of the chemokine genes CCL19 and CCL21. Based on re-sequencing screening we selected and, using PCR-SSP, determined three polymorphisms of CCL19 gene (GenBank ID rs2233872) and CCL21 gene (GenBank ID rs11574914 and rs11574915) in 211 Czech patients with MI and 150 healthy control subjects. There was no difference in allelic frequencies of the investigated SNPs between patients and controls (p>0.05). However, the proportion of homozygotes for the minor G allele of the CCL21 promoter variant (rs11574915 GG) was lower among the MI patients (1%) in comparison with the control subjects (5%, nominal p=0.03). Though rare in the Czech population, CCL21 (rs11574915) GG genotype may confer protection from myocardial infarction. Our preliminary data have to be independently replicated.

A novel HLA-DRB1 allele, HLA-DRB1*13:116, identified by sequencing-based typing in a member of the Czech National Marrow Donor Registry

We describe the identification of a novel HLA-DRB1 allele, DRB1*13:116, in a member of the Czech National Marrow Donor Registry. The novel allele differs from the known DRB1*13:17 variant by a nucleotide exchange at position 227 (T/A) of the coding HLA-DRB1 sequence, which causes an amino acid substitution (Phe47Tyr) in the HLA-DR beta 1 chain.