Search

Background: After introducing IL-1/IL-6 inhibitors, some patients with Still and Still-like disease developed unusual, often fatal, pulmonary disease. This complication was associated with scoring as DReSS (drug reaction with eosinophilia and systemic symptoms) implicating these inhibitors, although DReSS can be difficult to recognize in the setting of systemic inflammatory disease., Objective: To facilitate recognition of IL-1/IL-6 inhibitor-DReSS in systemic inflammatory illnesses (Still/Still-like) by looking at timing and reaction-associated features. We evaluated outcomes of stopping or not stopping IL-1/IL-6 inhibitors after DReSS reaction began., Methods: In an international study collaborating primarily with pediatric specialists, we characterized features of 89 drug-reaction cases versus 773 drug-exposed controls and compared outcomes of 52 cases stopping IL-1/IL-6 inhibitors with 37 cases not stopping these drugs., Results: Before the reaction began, drug-reaction cases and controls were clinically comparable, except for younger disease-onset age for reaction cases with preexisting cardiothoracic comorbidities. After the reaction began, increased rates of pulmonary complications and macrophage activation syndrome differentiated drug-reaction cases from drug-tolerant controls (P = 4.7 × 10 -35 and P = 1.1 × 10 -24 , respectively). The initial DReSS feature was typically reported 2 to 8 weeks after initiating IL-1/IL-6 inhibition. In drug-reaction cases stopping versus not stopping IL-1/IL-6-inhibitor treatment, reaction-related features were indistinguishable, including pulmonary complication rates (75% [39 of 52] vs 76% [28 of 37]). Those stopping subsequently required fewer medications for treatment of systemic inflammation, had decreased rates of macrophage activation syndrome, and improved survival (P = .005, multivariate regression). Resolution of pulmonary complications occurred in 67% (26 of 39) of drug-reaction cases who stopped and in none who continued inhibitors., Conclusions: In systemic inflammatory illnesses, recognition of IL-1/IL-6-inhibitor-associated reactions followed by avoidance of IL-1/IL-6 inhibitors significantly improved outcomes., (Copyright © 2024 American Academy of Allergy, Asthma & Immunology. All rights reserved.)

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Genetically determined variation of killer cell immunoglobulin like receptors (KIR) and their HLA class I ligands affects multiple aspects of human health. Their extreme diversity is generated through complex interplay of natural selection for pathogen resistance and reproductive health, combined with demographic structure and dispersal. Despite significant importance to multiple health conditions of differential effect across populations, the nature and extent of immunogenetic diversity is under-studied for many geographic regions. Here, we describe the first high-resolution analysis of KIR and HLA class I combinatorial diversity in Northern Africa. Analysis of 125 healthy unrelated individuals from Cairo in Egypt yielded 186 KIR alleles arranged in 146 distinct centromeric and 79 distinct telomeric haplotypes. The most frequent haplotypes observed were KIR-A, encoding two inhibitory receptors specific for HLA-C, two that are specific for HLA-A and -B, and no activating receptors. Together with 141 alleles of HLA class I, 75 of which encode a KIR ligand, we identified a mean of six distinct interacting pairs of inhibitory KIR and HLA allotypes per individual. We additionally characterize 16 KIR alleles newly identified in the study population. Our findings place Egyptians as one of the most highly diverse populations worldwide, with important implications for transplant matching and studies of immune-mediated diseases., (© 2023 The Authors. HLA: Immune Response Genetics published by John Wiley & Sons Ltd.)

Patients in need of hematopoietic stem cell transplantation often rely on unrelated stem cell donors matched in certain human leukocyte antigen (HLA) genes. Donor search is complicated by the extensive allelic variability of the HLA system. Therefore, large registries of potential donors are maintained in many countries worldwide. Population-specific HLA characteristics determine the registry benefits for patients and also the need for further regional donor recruitment. In this work, we analyzed HLA allele and haplotype frequencies of donors of DKMS Chile, the first Chilean donor registry, with self-assessed "non-Indigenous" ( n =92,788) and "Mapuche" ( n =1,993) ancestry. We identified HLA alleles that were distinctly more abundant in the Chilean subpopulations than in worldwide reference populations, four of them particularly characteristic for the Mapuche subpopulation, namely B*39:09g, B*35:09, DRB1*04:07g, and DRB1*16:02g. Both population subsamples carried haplotypes of both Native American and European origin at high frequencies, reflecting Chile's complex history of admixture and immigration. Matching probability analysis revealed limited benefits for Chilean patients (both non-Indigenous and Mapuche) from donor registries of non-Chilean donors, thus indicating a need for ongoing significant donor recruitment efforts in Chile., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Solloch, Giani, Pattillo Garnham, Sauter, Bernas, Lange, Barriga, Fernández-Viña and Schmidt.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Sikorova, Osoegawa, Kocourkova, Strnad, Petrkova, Fernández-Viña, Doubkova and Petrek.)

Supporting allogeneic hematopoietic cell transplantation (HCT) is an integral function of the clinical HLA laboratory, which provides HLA testing for recipients and donors. However, the timing, scope, and methods of the HLA tests vary significantly in the field. This summary provides a comprehensive and practical HLA testing approach to maximize the efficiency of the donor search process, minimize donor-specific HLA antibody (DSA) associated risks, enable optimal donor selections, and support HCT multidisciplinary teams. This is not a comprehensive donor selection guide, but pertinent donor selection considerations and publicly available online selection tools are highlighted. In the absence of healthy HLA identical siblings, younger 8/8 (HLA-A, -B, -C, -DRB1) HLA-matched unrelated donors remain the most favorable choice for HCT. Emerging practices in preparative regimens and graft versus host disease (GvHD) prophylaxis as well as building evidence of the importance of other HLA (e.g., HLA-DPB1 allele and functional matching) and non-HLA (e.g., age, CMV, and KIR) donor attributes urge the transplant centers and the HLA laboratories to construct a comprehensive approach for the routine histocompatibility testing., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022. Published by Elsevier Inc.)

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Next generation sequencing (NGS) is being applied for HLA typing in research and clinical settings. NGS HLA typing has made it feasible to sequence exons, introns and untranslated regions simultaneously, with significantly reduced labor and reagent cost per sample, rapid turnaround time, and improved HLA genotype accuracy. NGS technologies bring challenges for cost-effective computation, data processing and exchange of NGS-based HLA data. To address these challenges, guidelines and specifications such as Genotype List (GL) String, Minimum Information for Reporting Immunogenomic NGS Genotyping (MIRING), and Histoimmunogenetics Markup Language (HML) were proposed to streamline and standardize reporting of HLA genotypes. As part of the 17th International HLA and Immunogenetics Workshop (IHIW), we implemented standards and systems for HLA genotype reporting that included GL String, MIRING and HML, and found that misunderstanding or misinterpretations of these standards led to inconsistencies in the reporting of NGS HLA genotyping results. This may be due in part to a historical lack of centralized data reporting standards in the histocompatibility and immunogenetics community. We have worked with software and database developers, clinicians and scientists to address these issues in a collaborative fashion as part of the Data Standard Hackathons (DaSH) for NGS. Here we report several categories of challenges to the consistent exchange of NGS HLA genotyping data we have observed. We hope to address these challenges in future DaSH for NGS efforts., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

Genetic susceptibility to myasthenia gravis (MG) associates with specific HLA alleles and haplotypes at the class I and II regions in various populations. Previous studies have only examined alleles at a limited number of HLA loci that defined only broad serotypes or alleles defined at the protein sequence level. Consequently, genetic variants in noncoding and untranslated HLA gene segments have not been fully explored but could also be important determinants for MG. To gain further insight into the role of HLA in MG, we applied next-generation sequencing to analyze sequence variation at eleven HLA genes in early-onset (EO) and late-onset (LO) non-thymomatous MG patients positive for the acetylcholine receptor (AChR) antibodies and ethnically matched controls from Italy, Norway, and Sweden. For all three populations, alleles and haplotype blocks present on the ancestral haplotype AH8.1 were associated with risk in AChR-EOMG patients. HLA-B*08:01:01:01 was the dominant risk allele in Italians (OR = 3.28, P = 1.83E-05), Norwegians (OR = 3.52, P = 4.41E-16), and in Swedes HLA-B*08:01 was the primary risk allele (OR = 4.24, P <2.2E-16). Protective alleles and haplotype blocks were identified on the HLA-DRB7 , and HLA-DRB13.1 class II haplotypes in Italians and Norwegians, whereas in Swedes HLA-DRB7 exhibited the main protective effect. For AChR-LOMG patients, the HLA-DRB15.1 haplotype and associated alleles were significantly associated with susceptibility in all groups. The HLA-DR13-HLA-DR-HLA-DQ haplotype was associated with protection in all AChR-LOMG groups. This study has confirmed and extended previous findings that the immunogenetic predisposition profiles for EOMG and LOMG are distinct. In addition, the results are consistent with a role for non-coding HLA genetic variants in the pathogenesis of MG., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Creary, Gangavarapu, Caillier, Cavalcante, Frangiamore, Lie, Bengtsson, Harbo, Brauner, Hollenbach, Oksenberg, Bernasconi, Maniaol, Hammarström, Mantegazza and Fernández-Viña.)

The KIR (killer-cell immunoglobulin-like receptor ) region is characterized by structural variation and high sequence similarity among genes, imposing technical difficulties for analysis. We undertook the most comprehensive study to date of KIR genetic diversity in a large population sample, applying next-generation sequencing in 2,130 United States European-descendant individuals. Data were analyzed using our custom bioinformatics pipeline specifically designed to address technical obstacles in determining KIR genotypes. Precise gene copy number determination allowed us to identify a set of uncommon gene-content KIR haplotypes accounting for 5.2% of structural variation. In this cohort, KIR2DL4 is the framework gene that most varies in copy number (6.5% of all individuals). We identified phased high-resolution alleles in large multi-locus insertions and also likely founder haplotypes from which they were deleted. Additionally, we observed 250 alleles at 5-digit resolution, of which 90 have frequencies ≥1%. We found sequence patterns that were consistent with the presence of novel alleles in 398 (18.7%) individuals and contextualized multiple orphan dbSNPs within the KIR complex. We also identified a novel KIR2DL1 variant, Pro151Arg, and demonstrated by molecular dynamics that this substitution is predicted to affect interaction with HLA-C. No previous studies have fully explored the full range of structural and sequence variation of KIR as we present here. We demonstrate that pairing high-throughput sequencing with state-of-art computational tools in a large cohort permits exploration of all aspects of KIR variation including determination of population-level haplotype diversity, improving understanding of the KIR system, and providing an important reference for future studies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Amorim, Augusto, Nemat-Gorgani, Montero-Martin, Marin, Shams, Dandekar, Caillier, Parham, Fernández-Viña, Oksenberg, Norman and Hollenbach.)

Hidradenitis suppurativa (HS) is a chronic inflammatory cutaneous disease of the hair follicle typically presenting recurrent, painful, and inflamed lesions on the inverse areas of the body. Although its pathogenesis remains unknown, the immune system appears to play a potential role. To date, two previous studies have not found any association between the Human Leukocyte Antigen system (HLA) and HS. In this study we analyzed the HLA-A , -B, -C ; and DRB1 , - DQA1, and -DQB1 allele distribution in 106 HS patients and 262 healthy controls from a Caucasian population in Cantabria (northern Spain). HLA-A *29 and B *50 were significantly more common in HS patients and A *30 and B *37 in controls, but these associations disappeared after statistical correction. DRB1 *07, DQA1 *02, and DQB1 *02 were significantly more common in controls ( p 0.026, p 0.0012, and p 0.0005, respectively) and the HLA allele DQB1 *03:01 was significantly more common in HS patients ( p 0.00007) after the Bonferroni correction. The DRB1*07~DQA1*02~DQB1*02 haplotype was significantly more common in controls ( p < 0.0005). This is the first study showing an association between HLA-class II and HS. Our results suggest that HLA-II alleles ( DRB1 *07, DQA1 *02, DQB1 *02, and DQB1 *03:01) and the DRB1 *07~ DQA1 *02~ DQB1 *02 haplotype could influence resistance or susceptibility to HS.

The human leukocyte antigen (HLA) genes are extremely polymorphic and are useful molecular markers to make inferences about human population history. However, the accuracy of the estimation of genetic diversity at HLA loci very much depends on the technology used to characterize HLA alleles; high-resolution genotyping of long-range HLA gene products improves the assessment of HLA population diversity as well as other population parameters compared to lower resolution typing methods. In this study we examined allelic and haplotype HLA diversity in a large healthy European American population sourced from the UCSF-DNA bank. A high-resolution next-generation sequencing method was applied to define non-ambiguous 3- and 4-field alleles at the HLA-A, HLA-C, HLA-B, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1 loci in samples provided by 2248 unrelated individuals. A number of population parameters were examined including balancing selection and various measurements of linkage disequilibrium were calculated. There were no detectable deviations from Hardy-Weinberg proportions at HLA-A, HLA-DRB1, HLA-DQA1 and HLA-DQB1. For the remaining loci moderate and significant deviations were detected at HLA-C, HLA-B, HLA-DRB3/4/5, HLA-DPA1 and HLA-DPB1 loci mostly from population substructures. Unique 4-field associations were observed among alleles at 2 loci and haplotypes extending large intervals that were not apparent in results obtained using testing methodologies with limited sequence coverage and phasing. The high diversity at HLA-DPA1 results from detection of intron variants of otherwise well conserved protein sequences. It may be speculated that divergence in exon sequences may be negatively selected. Our data provides a valuable reference source for future population studies that may allow for precise fine mapping of coding and non-coding sequences determining disease susceptibility and allo-immunogenicity., (Copyright © 2019 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

The highly polymorphic classical human leukocyte antigen (HLA) genes display strong linkage disequilibrium (LD) that results in conserved multi-locus haplotypes. For unrelated individuals in defined populations, HLA haplotype frequencies can be estimated using the expectation-maximization (EM) method. Haplotypes can also be constructed using HLA allele segregation from nuclear families. It is straightforward to identify many HLA genotyping inconsistencies by visually reviewing HLA allele segregation in family members. It is also possible to identify potential crossover events when two or more children are available in a nuclear family. This process of visual inspection can be unwieldy, and we developed the "HaplObserve" program to standardize the process and automatically build haplotypes using family-based HLA allele segregation. HaplObserve facilitates systematically building haplotypes, and reporting potential crossover events. HLA Haplotype Validator (HLAHapV) is a program originally developed to impute chromosomal phase from genotype data using reference haplotype data. We updated and adapted HLAHapV to systematically compare observed and estimated haplotypes. We also used HLAHapV to identify haplotypes when uninformative HLA genotypes are present in families. Finally, we developed "pould", an R package that calculates haplotype frequencies, and estimates standard measures of global (locus-level) LD from both observed and estimated haplotypes., (Copyright © 2019 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

The 17th International HLA and Immunogenetics Workshop (IHIW) conducted a project entitled "The Study of Haplotypes in Families by NGS HLA". We investigated the HLA haplotypes of 1017 subjects in 263 nuclear families sourced from five US clinical immunogenetics laboratories, primarily as part of the evaluation of related donor candidates for hematopoietic stem cell and solid organ transplantation. The parents in these families belonged to five broad groups - African (72 parents), Asian (115), European (210), Hispanic (118) and "Other" (11). High-resolution HLA genotypes were generated for each subject using next-generation sequencing (NGS) HLA typing systems. We identified the HLA haplotypes in each family using HaplObserve, software that builds haplotypes in families by reviewing HLA allele segregation from parents to children. We calculated haplotype frequencies within each broad group, by treating the parents in each family as unrelated individuals. We also calculated standard measures of global linkage disequilibrium (LD) and conditional asymmetric LD for each ethnic group, and used untruncated and two-field allele names to investigate LD patterns. Finally we demonstrated the utility of consensus DNA sequences in identifying novel variants, confirming them using HLA allele segregation at the DNA sequence level., (Copyright © 2019 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2019 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

Next-generation sequencing (NGS) at the HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1 and -DRB3/4/5 loci was performed on 282 healthy unrelated individuals from different major regions of Spain. High-resolution HLA genotypes defined by full sequencing of class I loci and extended coverage of class II loci were obtained to determine allele frequencies and also to estimate extended haplotype frequencies. HLA alleles were typed at the highest resolution level (4-field level, 4FL); with exception of a minor deviation in HLA-DPA1, no statistically significant deviations from expected Hardy Weinberg Equilibrium (HWE) proportions were observed for all other HLA loci. This study provides new 4FL-allele and -haplotype frequencies estimated for the first time in the Spanish population. Furthermore, our results describe extended haplotypes (including the less frequently typed HLA-DPA1 and HLA-DQA1 loci) and show distinctive haplotype associations found at 4FL-allele definition in this Spanish population study. The distinctive allelic and haplotypic diversity found at the 4FL reveals the high level of heterozygosity and specific haplotypic associations displayed that were not apparent at 2-field level (2FL). Overall, these results may contribute as a useful reference source for future population studies, for HLA-disease association studies as a healthy control group dataset and for improving donor recruitment strategies of bone marrow registries. HLA genotyping data of this Spanish population cohort was also included in the 17th International Histocompatibility and Immunogenetics Workshop (IHIW) as part of the study of HLA diversity in unrelated worldwide populations using NGS., (Copyright © 2019 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2019. Published by Elsevier Inc.)

The human leukocyte antigen (HLA) genes are highly variable and are known to play an important role in disease outcomes, including infectious diseases. Prior knowledge of HLA polymorphisms in a population usually forms the basis for an effective case-control study design. As a prelude to future disease association analyses, we report HLA class I and II diversity in 334 unrelated donors from a Dengue vaccine efficacy trial conducted in Thailand. Long-range PCR amplification of six HLA loci was performed on DNA extracted from saliva samples. HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1 were genotyped using a next-generation sequencing method presented at the 17th International HLA and Immunogenetics Workshop. In total, we identified 201 HLA alleles, including 35 HLA-A, 57 HLA-B, 28 HLA-C, 24 HLA-DPB1, 21 HLA-DQB1 and 36 HLA-DRB1 alleles. Very common HLA alleles with frequencies greater than 10 percent were A∗11:01:01, A∗33:03:01, A∗24:02:01, B∗46:01:01, C∗07:02:01, C∗01:02:01, C∗08:01:01, DPB1∗05:01:01, DPB1∗13:01:01, DPB1∗04:01:01, DPB1∗02:01:02, DQB1∗03:01:01, DQB1∗05:02:01, DQB1∗03:03:02, DRB1∗12:02:01, DRB1∗09:01:02, and DRB1∗15:02:01. A novel HLA allele, B∗15:450, had a non-synonymous substitution and occurred in more than one donor. Population-based full-length NGS HLA typing is more conclusive and provides a sound foundation for exploring disease association in a given population., (Copyright © 2018. Published by Elsevier Inc.)

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., (© 2017 The Authors HLA: Immune Response Genetics Published by John Wiley & Sons Ltd.)

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., (© 2017 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

A single mismatch in highly expressed HLA-A, -B, -C, and -DRB1 loci (HEL) is associated with worse outcomes in hematopoietic stem cell transplantation, while less is known about the cumulative impact of mismatches in the lesser expressed HLA loci DRB3/4/5, DQ, and DP (LEL). We studied whether accumulation of LEL mismatches is associated with deleterious effects in 3853 unrelated donor transplants stratified according to number of matches in the HEL. In the 8/8 matched HEL group, LEL mismatches were not associated with any adverse outcome. Mismatches at HLA-DRB1 were associated with occurrence of multiple LEL mismatches. In the 7/8 HEL group, patients with 3 or more LEL mismatches scored in the graft-versus-host vector had a significantly higher risk of mortality (1.45 and 1.43) and transplant-related mortality (1.68 and 1.54) than the subgroups with 0 or 1 LEL mismatches. No single LEL locus had a more pronounced effect on clinical outcome. Three or more LEL mismatches are associated with lower survival after 7/8 HEL matched transplantation. Prospective evaluation of matching for HLA-DRB3/4/5, -DQ, and -DP loci is warranted to reduce posttransplant risks in donor-recipient pairs matched for 7/8 HEL.

DRB1 null alleles are extremely rare and always sporadic, suggesting their biological selective disadvantage.

Most archaeological and linguistic evidence suggest a Polynesian origin of the population of Easter Island (Rapanui), and this view has been supported by the identification of Polynesian mitochondrial DNA (mtDNA) polymorphisms in prehistoric skeletal remains. However, some evidence of an early South American contact also exists (the sweet potato, bottle gourd etc.), but genetic studies have so far failed to show an early Amerindian contribution to the gene pool on Easter Island. To address this issue, we analyzed mtDNA and Y chromosome markers and performed high-resolution human leukocyte antigen (HLA) genotyping of DNA harvested from previously collected sera of 48 reputedly nonadmixed native Easter Islanders. All individuals carried mtDNA types and HLA alleles previously found in Polynesia, and most men carried Y chromosome markers of Polynesian origin, providing further evidence of a Polynesian origin of the population of Easter Island. A few individuals carried HLA alleles and/or Y chromosome markers of European origin. More interestingly, some individuals carried the HLA alleles A*0212 and B*3905, which are of typical Amerindian origin. The genealogy of some of the individuals carrying these non-Polynesian HLA alleles and their haplotypic backgrounds suggest an introduction into Easter Island in the early 1800s, or earlier. Thus, there may have been an early European and Amerindian contribution to the Polynesian gene pool of Easter Island.

The allelic and haplotypic diversity of the HLA-A, HLA-B, and HLA-C loci was investigated in 852 subjects from five sub-Saharan populations from Kenya (Nandi and Luo), Mali (Dogon), Uganda, and Zambia. Distributions of genotypes at all loci and in all populations fit Hardy-Weinberg equilibrium expectations. There was not a single allele predominant at any of the loci in these populations, with the exception of A*3002 [allele frequency (AF) = 0.233] in Zambians and Cw*1601 (AF = 0.283) in Malians. This distribution was consistent with balancing selection for all class I loci in all populations, which was evidenced by the homozygosity F statistic that was less than that expected under neutrality. Only in the A locus in Zambians and the C locus in Malians, the AF distribution was very close to neutrality expectations. There were six instances in which there were significant deviations of allele distributions from neutrality in the direction of balancing selection. All allelic lineages from each of the class I loci were found in all the African populations. Several alleles of these loci have intermediate frequencies (AF = 0.020-0.150) and seem to appear only in the African populations. Most of these alleles are widely distributed in the African continent and their origin may predate the separation of linguistic groups. In contrast to native American and other populations, the African populations do not seem to show extensive allelic diversification within lineages, with the exception of the groups of alleles A*02, A*30, B*57, and B*58. The alleles of human leukocyte antigen (HLA)-B are in strong linkage disequilibrium (LD) with alleles of the C locus, and the sets of B/C haplotypes are found in several populations. The associations between A alleles with C-blocks are weaker, and only a few A/B/C haplotypes (A*0201-B*4501-Cw*1601; A*2301-B*1503-Cw*0202; A*7401-B* 1503-Cw*0202; A*2902-B*4201-Cw*1701; A*3001-B*4201-Cw*1701; and A*3601-B*5301-Cw*0401) are found in multiple populations with intermediate frequencies [haplotype frequency (HF) = 0.010-0.100]. The strength of the LD associations between alleles of HLA-A and HLA-B loci and those of HLA-B and HLA-C loci was on average of the same or higher magnitude as those observed in other non-African populations for the same pairs of loci. Comparison of the genetic distances measured by the distribution of alleles at the HLA class I loci in the sub-Saharan populations included in this and other studies indicate that the Luo population from western Kenya has the closest distance with virtually all sub-Saharan population so far studied for HLA-A, a finding consistent with the putative origin of modern humans in East Africa. In all African populations, the genetic distances between each other are greater than those observed between European populations. The remarkable current allelic and haplotypic diversity in the HLA system as well as their variable distribution in different sub-Saharan populations is probably the result of evolutionary forces and environments that have acted on each individual population or in their ancestors. In this regard, the genetic diversity of the HLA system in African populations poses practical challenges for the design of T-cell vaccines and for the transplantation medical community to find HLA-matched unrelated donors for patients in need of an allogeneic transplant.

The HLA system is the most polymorphic of all human genetic systems. The frequency of HLA class I alleles and their linkage disequilibrium patterns differ significantly among human populations as shown in studies using serologic methods. Many DNA-defined alleles with identical serotypes may have variable frequencies in different populations. We typed HLA-A, B, and C loci at the allele level by PCR-based methods in 1,296 unrelated subjects from five major outbred groups living in the U.S.A (African, AFAM; Caucasians, CAU; Asian, ORI; Hispanic, HIS, and North American Natives, NAI). We detected 46, 100 and 32 HLA-A, B, and C alleles, respectively. ORI and HIS presented more alleles at each of these loci. There was lack of correlation between the levels of heterozygosity and the number of alleles detected in each population. In AFAM, heterozygosity (>90%) is maximized at all class I loci. HLA-A had the lowest heterozygosity in all populations but CAU. Tight LD was observed between HLA-B and C alleles. AFAM had weaker or nonexistent associations between alleles of HLA-A and B than other populations. Analysis of the genetic distances between these and other populations showed a close relationship between specific US populations and a population from their original continents. ORI exhibited the largest genetic distance with all the other U.S. groups and were closer to NAI. Evidence of admixture with CAU was observed for AFAM and HIS. HIS also had significant frequencies of AFAM and Mexican Indian alleles. Differences in both LD and heterozygosity levels suggest distinct evolutionary histories of the HLA loci in the geographical regions from where the U.S. populations originated.

We have studied the HLA alleles of 60 unrelated healthy Terena and 10 Terena families. They are members of an isolated Brazilian tribe located in Mato Grosso do Sul (South Central Brazil). Six novel alleles were found in this population: HLA-A*0219 (gf = 0.02), A*0222 (gf = 0.15), HLA-B* 3520 (gf = 0.01), B*3521 (gf = 0.03), B*3912 (gf = 0.03) and B*4803 (gf = 0.16). Five of the six novel alleles differ from their putative progenitors by amino acid replacements in residues that contribute to the pockets of the peptide-binding site. Many of the variants defined by molecular methods were not identified correctly by serological typing. We calculated heterozygosity values (H) for HLA-A, -B, -C, DRB1, DQB1 and DPB . The highest values were observed at the HLA-B locus, followed by HLA-A, -DRB1 and DQB1. Residue positions 9, 24, 45, 62, 67, 95, 114, 116, 156, and 163 of HLA class I showed heterozygosity values greater than 0.50. Nine of them contribute to the peptide-binding specificity pockets and one to the T cell receptor binding site. If HLA antigens are useful for defense against pathogenic agents, heterozygosity would offer an advantage by allowing binding of a larger repertoire of peptides to the class I molecules. Individuals that are heterozygous at these positions would probably have a wider repertoire of peptide presentation to T cells. The observed results including the presence of novel alleles in the class I HLA loci suggest a functionally significant, more rapid evolution of class I compared to class II loci in this South American isolated population.

The human leukocyte antigen (HLA) complex includes the most polymorphic genes in humans. More than 600 allelic variants have been described in different populations. The HLA-B locus has contributed the largest number of alleles. Although Native American populations display a restricted number of HLA-alleles, many novel HLA class I alleles have been identified in indigenous communities of Central and South America. We have studied 248 unrelated individuals from three tribes of North-East Argentina and one from South-West Brazil, as well as 80 related individuals from the Brazilian tribe. In the course of this work, we found 8 new B-locus alleles and 2 novel A-locus alleles in these populations. Here we report the nucleotide sequences of A*0219, A*0222, B*3519, B*3520, B*3521, B*3912, B*4009 and B*4803 and we show their relationship with similar alleles. The new alleles B*35092 and B*3518 have been described by us in a previous paper. The possible mechanisms that may have produced these alleles over evolutionary time are discussed.

DNA typing systems for alleles of the HLA class I loci A, B, C at intermediate (IR) and high resolution (HR) levels were developed. The approaches combine locus-specific amplification of genomic DNA by the polymerase chain reaction (PCR) and hybridization with sequence-specific oligonucleotide probes (SSOP). The SSOP were designed to match nucleotide sequences at all polymorphic sites of exons 2 and 3 at these loci. Alleles and genotypes for these loci are assigned by their unique hybridization patterns. Some genotypes with particular allele combinations resulted in the same hybridization pattern. These genotype ambiguities were resolved by performing additional group-specific amplifications with appropriate GSA primers and hybridization with informative SSOP. At intermediate resolution level, many groups of alleles of HLA-A and B with the same serologic equivalence resulted in the same hybridization pattern. Both HR and IR typing approaches required the design, validation and testing of locus- and group-specific primers and sequence-specific oligonucleotide probes (SSOP). Single locus-specific amplification and hybridization with sets of 67 SSOP for HLA-A, 99 for HLA-B and 57 for HLA-C allowed us to identify unequivocally the majority of A, B, C alleles at HR level. To resolve ambiguous genotypes at HLA-B, we performed 4 GSA with 5' primers at codon 45-46 and hybridization with selected sets of SSOP. About 22,415 high resolution typing results were obtained (4,953A, 6,621B, 10,841C). In these samples, 63 HLA-A, 170 HLA-B and 40 HLA-C alleles were observed. In the course of these studies, more than 30 new alleles have been identified. In IR testing, sets of 39 SSOP for HLA-A typing and 59 SSOP for HLA-B typing allowed us to obtain maximal resolution of the majority of common genotypes. To achieve IR level, the majority of SSOP selected were those that span codons encoding amino acid residues located in the alpha helical segments of the class I molecules. A total of 50,522 samples were typed for HLA-A and B at IR level. Approximately 2.0% of them carried ambiguous genotypes associated with alternative switches of Bw4/w6 related sequences. All these ambiguities were resolved by Bw4/Bw6 GSA by 2 primer pairs (77N-IALR-83/3B.1; 77S-DLRG-83/3B.1) and hybridization with 9 selected SSOP Testing was performed by individual hybridization of replicate dot blot membranes with the SSOP of the corresponding set. The approach was robust and cost-effective in large-scale HLA class I molecular typing. The resolution provided by HLA-A, B IR reached serologic split-level or higher. The description of primers and probes for HLA class I typing may be utilized as starting elements for development of second generation methods with a more rapid turn around time.

To assess the influence of HLA on AIDS-free survival, human immunodeficiency virus load, and CD4 cell counts, 91 Caucasian and 48 African-American seroprevalent men were typed for HLA classes I and II and TAP alleles. HLA associations with these markers were assessed by assigning sum integer scores based on 7 class I allele-TAP variants (+1) and 13 class I-class II-TAP combinations (-1) with different AIDS-free survival times found in a prior study. Subjects in both racial groups and combined with positive sum scores were less likely to have CD4 cell decline (P=.0004), to have increased virus burden (P=.014), and to develop AIDS (P=.034) in the follow-up period than were Caucasians and African Americans with scores of 0 or -1. These results confirm the reported associations of specific major histocompatibility complex genes with AIDS-free survival time in Caucasians and specifically extend them to African Americans and to two established markers of disease progression.

Previous studies have shown that susceptibility to Pauciarticular Juvenile Arthritis is associated with HLA-A*0201. Recently, autoantibodies against the protein of the DEK oncogene have been found in sera of patients with this disease. If T cells are involved in the pathogenesis of joint lesions, it is possible that they target autoantigens presented by HLA-A*0201. Therefore, we investigated whether DEK-derived peptides can bind efficiently to HLA-A*0201. Nonameric peptides selected considering anchor positions 2 and 9, and preferred amino acids at other positions, were incubated either with the human TAP-deficient cell line 174CEM.T2 (T2) or with the homozygous B cell line JESTHOM (A*0201, B*2705, Cw1), previously depleted of endogenous peptides. Binding was measured as the increase of detection of fully assembled, HLA-A*0201 molecules by flow cytometry with the anti-HLA-A2 monoclonal antibody (mAb) BB7.2. Three out of ten selected DEK-derived peptides showed binding to HLA-A*0201, which was peptide concentration-dependent (1 microM to 100 microM). DEK155-163 (AMLKSICEV), which also has two preferred amino acid residues at positions 6 and 8, yielded the highest binding. DEK163-171 (VLDLERSGV) and DEK72-80 (SLQREPFTI), which also has one preferred amino acid residue at position 8, also were able to bind to HLA-A*0201. Furthermore, peptide-induced, fully assembled, HLA-A*0201 molecules were immunoprecipitated with the BB7.2 mAb from metabolically-labeled T2 cells incubated with DEK72-80, DEK155-163, and DEK163-171. A faint band was observed in the immunoprecipitates of cells incubated with DEK65-73 (it carries a preferred amino acid residue at position 6), suggesting that this peptide interacts weakly with HLA-A*0201. These results indicate that several nonameric peptides derived from the DEK protein can bind to HLA-A 0201 and suggest that the complexes formed may be able to stimulate CD8+ T cells in patients with Pauciarticular Juvenile Arthritis.

MICA is a new polymorphic gene in the HLA region expressed in epithelial cell lines and gastrointestinal epithelium. Little is yet known about the MICA protein, and the pattern of its expression by freshly isolated cells has not been established. In the present experiments, we used antibodies raised in rabbits against alpha1 and alpha2 domain-peptides to study the expression of MICA. By western blot and immunoprecipitation, we detected a band of 62 000 Mr in various cell lines (THP-1, U937, HeLa, A431, Raji, MOLT-4, and HUV-EC-C) and in freshly isolated keratinocytes, endothelial cells, and monocytes but not in CD4+ and CD8+ T cells, and CD19+ cells (B lymphocytes). It was not possible to up-regulate the expression of MICA in different cells by stimulation with gamma-interferon, but the expression of MICA was induced in phytohemagglutinin-stimulated T cells. We confirmed that MICA is expressed at the cell surface by flow cytometry. Results of immunoprecipitation studies of beta2-microglobulin (beta2m)- or MICA-depleted, metabolically labeled HeLa cells indicated that MICA was not associated with beta2m. Although the function of MICA is still unknown, its restricted pattern of tissue expression, the fact that it is expressed on the cell surface, and its polymorphic nature suggest that this new molecule, encoded close to HLA class I, may play a role in the interaction between epithelial cells and cells of the immune system.

Native American populations have a limited HLA polymorphism compared with other ethnic groups. In spite of this, many novel HLA-B locus alleles, not observed in other populations, have been identified in South American tribes, and rapid evolution of this locus has been suggested. We have studied unrelated subjects of the Toba (TOB n = 116), Wichi (WIC n = 46) and Pilaga (PIL n = 14) tribes from northeastern Argentina to investigate the extent of the HLA polymorphism and obtain clues of selective forces that may have acted in these populations. In these tribes the number of HLA alleles is small at all loci except HLA-B, which presents 22 alleles. Seven novel alleles were characterized including 5 of HLA-B (B*35092, B*3518, B*3519, B*4009, B*4803) 1 at HLA-A (A*0219) and 1 at DRB1 (DRB1*0417). All these variants may have arisen by gene conversion events. Some of the novel variants represent the most frequent alleles of these populations (B*4803 in TOB and PIL; B*3519 in WIC) or are the most frequent subtypes in their lineages. HLA-A, B, DRB1,DQA1 and DQB1, but not DPB1, display relatively similar gene frequencies. This results in high heterozygosity in all the tribes for all the loci studied except HLA-DPB1. The larger polymorphism and the generation and maintenance of novel alleles at the HLA-B locus suggests a more specialized response of this locus to evolutionary forces. These effects may be related to the nature of the polymorphism, to the number of founder alleles and to the functional characteristics of the individual alleles.