Behçet disease, new insights in disease associations and manifestations: a next-generation sequencing study.

Behçet disease is a multi-system disease associated with human leukocyte antigen (HLA) class I polymorphism. High-resolution next-generation sequencing (NGS) with haplotype analysis has not been performed previously for this disease. Sixty Egyptian patients diagnosed according to the International Study Group (ISG) criteria for Behçet disease and 160 healthy geographic and ethnic-matched controls were genotyped for HLA class I loci (HLA-A, B, C). For HLA class II loci (DRB1, DRB3/4/5, DQA1, DQB1, DPA1, DPB1), 40 control samples were genotyped. High-resolution HLA genotyping was performed using NGS and the results were analyzed. Clinical manifestations were oral ulcers (100%), genital ulcers (100%), eye (55%) and neurological (28%) and vascular involvement (35%). HLA-B*51:08 [odds ratio (OR) = 19·75, 95% confidence interval (CI) = 6·5-79; P < 0·0001], HLA-B*15:03 (OR = 12·15, 95% CI = 3·7-50·7; P < 0·0001), HLA-C*16:02 (OR = 6·53, 95% CI = 3-14; P < 0·0001), HLA-A*68:02 (OR = 3·14, 95% CI = 1·1-8·9; P < 0·01) were found to be associated with Behçet disease, as were HLA-DRB1*13:01 and HLA-DQB1*06:03 (OR = 3·39, 95% CI = 0·9-18·9; P = 0·04 for both). By contrast, HLA-A*03:01 (OR = 0·13, 95% CI = 0-0·8; P = 0·01) and HLA-DPB1*17:01 were found to be protective (OR = 0·27, 95% CI = 0·06-1·03; P = 0·02). We identified strong linkage disequilibrium between HLA-B*51:08 and C*16:02 and A*02:01 in a haplotype associated with Behçet disease. HLA-B*51:08 was significantly associated with legal blindness (OR = 2·98, 95% CI = 1·06-8·3; P = 0·01). In Egyptian Behçet patients, HLA-B*51:08 is the most common susceptibility allele and holds poor prognosis for eye involvement.

Validation of statistical imputation of allele-level multilocus phased genotypes from ambiguous HLA assignments.

Genetic matching for loci in the human leukocyte antigen (HLA) region between a donor and a patient in hematopoietic stem cell transplantation (HSCT) is critical to outcome; however, methods for HLA genotyping of donors in unrelated stem cell registries often yield results with allelic and phase ambiguity and/or do not query all clinically relevant loci. We present and evaluate a statistical method for in silico imputation of HLA alleles and haplotypes in large ambiguous population data from the Be The Match(®) Registry. Our method builds on haplotype frequencies estimated from registry populations and exploits patterns of linkage disequilibrium (LD) across HLA haplotypes to infer high resolution HLA assignments. We performed validation on simulated and real population data from the Registry with non-trivial ambiguity content. While real population datasets caused some predictions to deviate from expectation, validations still showed high percent recall for imputed results with average recall >76% when imputing HLA alleles from registry data. We simulated ambiguity generated by several HLA genotyping methods to evaluate the imputation performance on several levels of typing resolution. On average, imputation percent recall of allele-level HLA haplotypes was >95% for allele-level typing, >92% for intermediate resolution typing and >58% for serology (low-resolution) typing. Thus, allele-level HLA assignments can be imputed through the application of a set of statistical and population genetics inferences and with knowledge of haplotype frequencies and self-identified race and ethnicities.

16(th) IHIW: extending the number of resources and bioinformatics analysis for the investigation of HLA rare alleles.

Continuing a project presented at the 15th International HLA and Immunogenetics Workshop (IHIWS) on the rarity of HLA alleles, we sought to expand the number of data sources and bioinformatics tools available in the Allele Frequencies Net Database website (AFND, www.allelefrequencies.net). In this 16th IHIWS Rare Alleles project, HLA alleles described in the latest IMGT/HLA Database (release 3.8.0) were queried against different sources including data from registries (stem cell) and from 74 different laboratories around the world. We demonstrated that approximately 40% of the alleles officially named in the IMGT/HLA Database have been reported only once across all different sources. To facilitate the large-scale analysis of rare alleles, we have produced an online tool called the Rare Allele Detector that simplifies the detection of alleles that are considered to be 'very rare', 'rare' or 'frequent'. Tools and associated data can be accessed via the www.allelefrequencies.net website.

16(th) IHIW: global distribution of extended HLA haplotypes.

This report describes the project to identify the global distribution of extended HLA haplotypes, a component of 16th International HLA and Immunogenetics Workshop (IHIW), and summarizes the initial analyses of data collected. The project aims to investigate extended HLA haplotypes, compare their distribution among different populations, assess their frequency in hematopoietic stem cell unrelated donor registries and initiate an international family studies database and DNA repository to be made publicly available. HLA haplotypes compiled in immunogenetics laboratories during the evaluation of transplant candidates and related potential donors were analysed. Haplotypes were determined using the pedigree analysis tool publicly available from the National Marrow Donor Program (NMDP) website. Nineteen laboratories from 10 countries (11 laboratories from North America, five from Asia, two from Latin America and one from Australia) contributed data on a total of 1719 families comprised of 7474 individuals. We identified 10393 HLA haplotypes, of which 1682 haplotypes included high-resolution typing at HLA-A, B, C, DRB1 and DQB1 loci. We also present haplotypes containing MICA and other HLA loci and haplotypes containing rare alleles seen in these families. The project will be extended through the 17th IHIW, and investigators interested in joining the project may communicate with the first author.

Association between sarcoidosis and HLA polymorphisms in a Czech population from Central Europe: focus on a relationship with clinical outcome and treatment.

Background: Sarcoidosis is an immune-mediated systemic disease with unknown etiology affecting the lung predominantly. The clinical manifestation of sarcoidosis is rather diverse ranging from Löfgren's syndrome to fibrotic disease. Also, it differs among patients with distinct geographical and ethnic origins, consistent with environmental and genetic factors' role in its pathogenesis. Of those, the polymorphic genes of the HLA system have been previously implicated in sarcoidosis. Therefore, we have performed an association study in a well-defined cohort of Czech patients aiming to define how variation in HLA genes, may contribute to disease origin and development. Materials and methods: Total of the 301 Czech unrelated sarcoidosis patients were diagnosed according to international guidelines. In those, HLA typing was performed using next-generation sequencing. The allele frequencies at six HLA loci (HLA-A,-B,-C,-DRB1,-DQA1, and -DQB1) observed in the patients were compared with HLA allele distribution determined in 309 unrelated healthy Czech subjects; sub-analyses of relationships between HLA and distinct sarcoidosis clinical phenotypes were performed. Associations were assessed by two-tailed Fischer's exact test with correction for multiple comparisons. Results: We report two variants, HLA-DQB1*06:02, and HLA-DQB1*06:04, as risk factors for sarcoidosis, and three variants, HLA-DRB1*01:01, HLA-DQA1*03:01, and HLA-DQB1*03:02 as protective factors. HLA-B*08:01, HLA-C*07:01, HLA-DRB1*03:01, HLA-DQA1*05:01, and HLA-DQB1*02:01 variants associated with Löfgren's syndrome, a more benign phenotype. HLA- DRB1*03:01 and HLA-DQA1*05:01 alleles were connected with better prognosis-chest X-ray (CXR) stage 1, disease remission, and non-requirement of corticosteroid treatment. The alleles HLA-DRB1*11:01 and HLA-DQA1*05:05 are associated with more advanced disease represented by the CXR stages 2-4. HLA-DQB1*05:03 associated with sarcoidosis extrapulmonary manifestation. Conclusion: In our Czech cohort, we document some associations between sarcoidosis and HLA previously described in other populations. Further, we suggest novel susceptibility factors for sarcoidosis, such as HLA-DQB1*06:04, and characterize associations between HLA and sarcoidosis clinical phenotypes in Czech patients. Our study also extends the role of the 8.1 ancestral haplotype (HLA-A*01:01∼HLA-B*08:01∼HLA-C*07:01∼HLA-DRB1*03:01∼HLA-DQA1*05:01∼HLA-DQB1*02:01), already implicated in autoimmune diseases, as a possible predictor of better prognosis in sarcoidosis. The general translational application of our newly reported findings for personalized patient care should be validated by an independent study from another, international referral center.

Association of HLA alleles with Plasmodium falciparum severity in Malian children.

Pre-erythrocytic immunity to Plasmodium falciparum malaria is likely to be mediated by T-cell recognition of malaria epitopes presented on infected host cells via class I and II major histocompatibility complex (MHC) antigens. To test for associations of human leukocyte antigen (HLA) alleles with disease severity, we performed high-resolution typing of HLA class I and II loci and compared the distributions of alleles of HLA-A, -B, -C and -DRB1 loci in 359 Malian children of Dogon ethnicity with uncomplicated or severe malaria. We observed that alleles A*30:01 and A*33:01 had higher frequency in the group of patients with cerebral disease compared to patients with uncomplicated disease [A*30:01: gf = 0.2031 vs gf = 0.1064, odds ratio (OR) = 3.17, P = 0.004, confidence interval (CI) (1.94-5.19)] and [A*33:01: gf = 0.0781 vs gf = 0.0266, 4.21, P = 0.005, CI (1.89-9.84)], respectively. The A*30:01 and A*33:01 alleles share some sequence motifs and A*30:01 appears to have a unique peptide binding repertoire compared to other A*30 group alleles. Computer algorithms predicted malaria peptides with strong binding affinity for HLA-A*30:01 and HLA-A*33:01 but not to closely related alleles. In conclusion, we identified A*30:01 and A*33:01 as potential susceptibility factors for cerebral malaria, providing further evidence that polymorphism of MHC genes results in altered malaria susceptibility.

Genetic differentiation of Jewish populations.

The Jewish diaspora can be viewed as a natural process in population dispersion and differentiation. We extend genetic studies on the Jewish diaspora to an analysis of human leukocyte antigen (HLA) haplotype distributions in the Jewish peoples, and show the value of this information for the design of Jewish marrow donor registries. HLA data from the Hadassah Bone Marrow Registry having parental country-of-origin information comprise samples of geographically discrete regions. We analyzed the HLA allele and haplotype frequencies for each national sample using population genetic and clustering methods. Population differentiation among diaspora populations was shown on the basis of HLA haplotype frequencies, including differences within the more recently diverged European groups. A method of haplotype and population clustering showed patterns of unique haplotype affinities associated with specific Jewish populations. The evidence showed that diaspora Jewish populations can be sorted into distinct clades of which the Ashkenazi are but one. Relationships among Jewish populations are interpretable in light of the historical record. We suggest that a major contributing factor to the genetic divergence between Jewish groups may have been admixture with local host populations, while, at the same time, threads of Eastern Mediterranean ancestry remain evident.

Further evidence of an Amerindian contribution to the Polynesian gene pool on Easter Island.

Available evidence suggests a Polynesian origin of the Easter Island population. We recently found that some native Easter Islanders also carried some common American Indian (Amerindian) human leukocyte antigen (HLA) alleles, which probably were introduced before Europeans discovered the island in 1722. In this study, we report molecular genetic investigations of 21 other selected native Easter Islanders. Analysis of mitochondrial DNA and Y chromosome markers showed no traces of an Amerindian contribution. However, high-resolution genomic HLA typing showed that two individuals carried some other common Amerindian HLA alleles, different from those found in our previous investigations. The new data support our previous evidence of an Amerindian contribution to the gene pool on Easter Island.

A bioinformatics approach to ascertaining the rarity of HLA alleles.

A project of the 15th International Histocompatibility Workshop examined the rarity of human leukocyte antigen (HLA) alleles. A section was constructed in the website, www.allelefrequencies.net to contain this data from different sources. A mechanism to search the data was implemented for use by any individual.

Gene for familial psoriasis susceptibility mapped to the distal end of human chromosome 17q.

A gene involved in psoriasis susceptibility was localized to the distal region of human chromosome 17q as a result of a genome-wide linkage analysis with polymorphic microsatellites and eight multiply affected psoriasis kindreds. In the family which showed the strongest evidence for linkage, the recombination fraction between a psoriasis susceptibility locus and D17S784 was 0.04 with a maximum two-point lod score of 5.33. There was also evidence for genetic heterogeneity and although none of the linked families showed any association with HLA-Cw6, two unlinked families showed weak levels of association. This study demonstrates that in some families, psoriasis susceptibility is due to variation at a single major genetic locus other than the human lymphocyte antigen locus.

Exploring the ancestry and admixture of Mexican Oaxaca Mestizos from Southeast Mexico using next-generation sequencing of 11 HLA loci.

The Mestizos of Oaxaca resulted from the admixture of Zapotecan Natives with Spaniards and Africans. We selected 112 donors from Oaxaca and applied next-generation sequencing to characterize exon and intron variants in complete or extended HLA genes. Some alleles found, are unique to Mexican Natives and most likely will be absent in most major ethnicities, namely: Caucasians, Africans or Asians. Among these are HLA-A*68:03:01, HLA-A*68:05:01, HLA-C*03:04:01:02, HLA-C*15:09, HLA-C*3:05, HLA-C*03:06:01, HLA-B*39:05:01, HLA-B*35:14:01, HLA-B*35:12:01, HLA-B*35:43:01, HLA-B*40:05, HLA-B:40:08, HLA-B*51:02:01, HLA-B*35:24:01 and HLA-B*39:08. HLA-DQA1*05:05:01:05 and some HLA-DRB1 alleles were only present in Amerindians/Mestizos. Three haplotypes are unique to Mexican Natives, five to Middle-Eastern and Sephardi-Jews. We detected a novel HLA-DQA1*04:01:01 exon 4 variant. Any novel allele may have been positively selected to enlarge the peptide-binding repertoire, and some, like HLA-B*39:02:02 and HLA-B*39:05:01 were found with unique haplotype associations, suggesting convergent evolution events and/or allele lineage diversification. The allele frequencies were fairly evenly distributed in most HLA loci with the exception of HLA-DPB1. The application of NGS in Oaxaca is novel and will lead to better use in the clinical setting. It offers deep knowledge on the population structure, origins, migration, and discovery of new alleles and haplotypes that other techniques did not achieve.

A novel HLA-B allele, B*4093, sharing sequences with B*4006 and B*0702 in an Asian Indian donor.

A new HLA-B allele (B*4093) in a North Indian Hindu donor differing from B*4006 by four nucleotide substitutions in codon 41.1, codon 44.3, codon 45.1 and codon 50.3 has been identified. This novel allele was part of the A*0211-B*4093-Cw*1502-DRB1*15-DQB1*06 HLA haplotype.

Strategies and technical challenges in allele level Class II typing in 2578 bone marrow transplantation donor-recipient pairs.

Human leukocyte antigen typing of 2578 donor-recipient pairs whose transplantation was facilitated by the National Marrow Donor Program allowed for an in-depth analysis of the accuracy of high-volume allele level testing data. The methods employed provided allele level typing at DRB1/3/5, DQA1, DQB1, DPA1, and DPB1 using sequence-specific oligonucleotide probe hybridization (SSOPH), polymerase chain reaction (PCR) restriction fragment length polymorphism analysis, sequence specific PCR, and direct sequence-based typing (SBT). Each typing was independently tested by two laboratories in Phase 1, and in subsequent phases targeted samples were typed in duplicate by SBT to monitor typing quality. Comparison with prior transplant center typing was also evaluated. SSOPH detected discrepancies ranged from 0.6% at DPB1 to 5.1% at DQB1 in Phase 1. The majority of discrepancies, 62%, resulted from human error such as sample handling, result interpretation, or clerical errors. Alleles that are frequently discrepant have been identified in this predominantly white population.

Deciphering the fine nucleotide diversity of full HLA class I and class II genes in a well-documented population from sub-Saharan Africa.

With the aim to understand how next-generation sequencing (NGS) improves both our assessment of genetic variation within populations and our knowledge on HLA molecular evolution, we sequenced and analysed 8 HLA loci in a well-documented population from sub-Saharan Africa (Mandenka). The results of full-gene NGS-MiSeq sequencing compared with those obtained by traditional typing techniques or limited sequencing strategies showed that segregating sites located outside exon 2 are crucial to describe not only class I but also class II population diversity. A comprehensive analysis of exons 2, 3, 4 and 5 nucleotide diversity at the 8 HLA loci revealed remarkable differences among these gene regions, notably a greater variation concentrated in the antigen recognition sites of class I exons 3 and some class II exons 2, likely associated with their peptide-presentation function, a lower diversity of HLA-C exon 3, possibly related to its role as a KIR ligand, and a peculiar molecular diversity of HLA-A exon 2, revealing demographic signals. Based on full-length HLA sequences, we also propose that the most frequent DRB1 allele in the studied population, DRB1*13:04, emerged from an allelic conversion involving 3 potential alleles as donors and DRB1*11:02:01 as recipient. Finally, our analysis revealed a high occurrence of the DRB1*13:04-DQA1*05:05:01-DQB1*03:19 haplotype, possibly resulting from a selective sweep due to protection to Onchorcerca volvulus, a prevalent pathogen in West Africa. This study unveils highly relevant information on the molecular evolution of HLA genes in relation to their immune function, calling for similar analyses in other populations living in contrasting environments.

Development and characterization of desmoglein-3 specific T cells from patients with pemphigus vulgaris.

Pemphigus vulgaris (PV) is a cutaneous autoimmune disease characterized by blister formation in the suprabasilar layers of skin and mucosae and anti-desmoglein-3 (Dsg3) autoantibodies bound to the surface of lesional keratinocytes and circulating in the serum of patients. This disease can be reproduced in neonatal mice by passive transfer of patients' IgG, indicating that humoral immunity plays an important role in the pathogenesis of PV. Currently, the role of T lymphocytes in the development of PV is not clear. Here, we report that three immunoreactive segments of the ectodomain of Dsg3 specifically induced proliferation of T cells from PV patients. We found that T lymphocytes from 13 out of 14 patients responded to at least one of three Dsg3 peptides. T cells from controls and other patient groups did not respond to these Dsg3 peptides. The major T cell population stimulated by these Dsg3 peptides was CD4 positive. Dsg3-specific T cell lines and clones were developed and were shown to express a CD4 positive memory T cell phenotype. Upon stimulation, these cell lines and clones secreted a Th2-like cytokine profile. The Dsg3 responses of these T cells were restricted to HLA-DR, and not -DQ and -DP, of the major histocompatibility complex. This information will help to elucidate the cellular immune abnormalities leading to production of pathogenic IgG autoantibodies in patients with PV.

Fludarabine with pharmacokinetically guided IV busulfan is superior to fixed-dose delivery in pretransplant conditioning of AML/MDS patients.

We hypothesized that IV busulfan (Bu) dosing could be safely intensified through pharmacokinetic (PK-) dose guidance to minimize the inter-patient variability in systemic exposure (SE) associated with body-sized dosing, and that this should improve outcome of AML/MDS patients undergoing allogeneic stem cell transplantation. To test this hypothesis, we treated 218 patients (median age 50.7 years, male/female 50/50%) with fludarabine 40 mg/m2 once daily x4, each dose followed by IV Bu, randomized to 130 mg/m2 (N=107) or PK-guided to average daily SE, AUC of 6000 µM min (N=111), stratified for remission status and allo-grafting from HLA-matched donors. Toxicity and GvHD rates in the groups were similar; the risk of relapse or treatment-related mortality remained higher in the fixed-dose group throughout the 80-month observation period. Further, PK-guidance yielded safer disease control, leading to improved overall and PFS, most prominently in MDS patients and in AML patients not in remission at allogeneic stem cell transplantation. We conclude that AML/MDS patients receiving pretransplant conditioning treatment with our 4-day regimen may benefit significantly from PK-guided Bu dosing. This could be considered an alternative to fixed-dose delivery since it provides the benefit of precise dose delivery to a predetermined SE without increasing risk(s) of serious toxicity and/or GvHD.

High-throughput next-generation sequencing to genotype six classical HLA loci from 96 donors in a single MiSeq run.

Next generation sequencing (NGS) methods have been established as an efficient approach for HLA typing because unlike traditional Sanger sequencing, they provide unambiguous results at a reasonable cost. We previously developed a multi-locus index method to genotype four HLA loci (A, B, C, and DRB1) on the Illumina MiSeq platform. We have now expanded this method to include two additional loci, HLA-DPB1 and DQB1. Contiguous full-length amplicons from 5'UTR through 3'UTR regions were generated using one long-range PCR reaction per locus for each of the six loci from 96 individuals of different ethnicities. The six amplicons from each donor were pooled, enzymatically fragmented and given a donor-specific index. This approach enabled sequencing of 576 loci from 96 individuals in a single MiSeq run. Donor-specific sequence reads were demultiplexed, and allele calls were generated from FASTQ files using commercially available software. Comparison to HLA genotypes generated from Sanger sequence-based typing (SBT) identified no discordances among any of the alleles analyzed in this study. Importantly, this method was able to resolve 22 DPB1 and 20 DQB1 alleles that were ambiguous with the SBT method. Furthermore, a novel allele in each of these two loci was identified, with the DQB1*05:01:24 allele having a frequency of greater than five percent. This method was subsequently validated against a blinded panel of 22 samples from the 17th International HLA and Immunogenetics Workshop. The flexibility of the method is further highlighted by successful genotyping of eight loci comprising all classical HLA loci for a subset of the samples. We now present a high-throughput, high-resolution, scalable NGS HLA typing method to accurately and efficiently genotype all classical HLA class I and II loci.

Evidence for selection against human lung cancers bearing p53 missense mutations which occur within the HLA A*0201 peptide consensus motif.

Short peptide fragments of intracellular proteins that fit a defined sequence motif bind to the most common human major histocompatibility complex class I molecule, HLA A*0201, and mediate killing by cytotoxic T-cells [D.F. Hunt et al., Science (Washington DC), 255: 1261-1263, 1992; K. Falk et al., Nature (Lond.), 351: 290-296, 1991]. The existence of such a motif allows prediction of whether novel peptides derived from mutant oncoporteins might be presented on the surface of cancer cells bearing that HLA allele. Clinical cancer might develop only when these mutations occur outside a major histocompatibility complex binding motif or in those cells that acquire defects in antigen presentation. Here, we find that missense mutations of p53 from a variety of tumors fall within the HLA A*0201 motif less often than would be expected if the location of mutations and motifs were independent. When we analyzed the HLA subtype of lung cancer cell lines with known p53 missense mutations, we found that all of the mutant oncopeptides predicted to be presentable by HLA A*0201 came from tumors that either did not carry the A*0201 allele or had lost that allele in the process of tumorigenesis. Presentation of mutant oncogene peptides on class I major histocompatibility complex might thus represent a physiologically significant selection pressure in the development of human cancer.

Evolution of HLA class II molecules: Allelic and amino acid site variability across populations.

Analysis of the highly polymorphic beta1 domains of the HLA class II molecules encoded by the DRB1, DQB1, and DPB1 loci reveals contrasting levels of diversity at the allele and amino acid site levels. Statistics of allele frequency distributions, based on Watterson's homozygosity statistic F, reveal distinct evolutionary patterns for these loci in ethnically diverse samples (26 populations for DQB1 and DRB1 and 14 for DPB1). When examined over all populations, the DQB1 locus allelic variation exhibits striking balanced polymorphism (P < 10(-4)), DRB1 shows some evidence of balancing selection (P < 0.06), and while there is overall very little evidence for selection of DPB1 allele frequencies, there is a trend in the direction of balancing selection (P < 0.08). In contrast, at the amino acid level all three loci show strong evidence of balancing selection at some sites. Averaged over polymorphic amino acid sites, DQB1 and DPB1 show similar deviation from neutrality expectations, and both exhibit more balanced polymorphic amino acid sites than DRB1. Across ethnic groups, polymorphisms at many codons show evidence for balancing selection, yet data consistent with directional selection were observed at other codons. Both antigen-binding pocket- and non-pocket-forming amino acid sites show overall deviation from neutrality for all three loci. Only in the case of DRB1 was there a significant difference between pocket- and non-pocket-forming amino acid sites. Our findings indicate that balancing selection at the MHC occurs at the level of polymorphic amino acid residues, and that in many cases this selection is consistent across populations.