The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection.

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.

Exceptional diversity of KIR and HLA class I in Egypt.

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.

High-resolution DNA methylation screening of the major histocompatibility complex in multiple sclerosis.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods. Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls. Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen (HLA) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk. Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS.

HLA allele and haplotype frequencies of registered stem cell donors in Chile.

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.

Spherotech-EDTA combined serum treatment reduces background more effectively as compared to One Lambda Adsorb Out™ and LIFECODES Serum Cleaner in Luminex-based solid-phase immunoassays for HLA antibody detection.

Spherotech (SPT) microparticles capture non-specific binding materials present in test serum, and EDTA removes the so called" prozone effect". This study presents a novel approach of combined SPT-EDTA serum treatment prior to Luminex HLA antibody testing to remove high background, and prozone effect in a single step process, and compared the efficacy of SPT-EDTA serum pre-treatment with AdsorbOut (ADS) and Serum Cleaner (SC) to reduce background in solid phase immunoassays (SPI). A total of 21 serum samples with a history of elevated negative control (NC) values ≥500, and 20 samples with normal NC values were included to assess the potential adverse effects. A problem of high background was noted in 25% of our samples in SPI. We observed 80% effectiveness in reducing NC values <500 with SPT-EDTA serum pre-treatment compared to 72%, and 67% for ADS and SC-treated sera, respectively. Interestingly, we found a strong correlation in antibody-binding levels between SPT versus ADS; and ADS versus SC-treated sera for both phenotype and single antigen bead assays (p < 0.001). No adverse effect was noted on NC, positive control (PC) values, PC/NC ratios in the upfront use of SPT-EDTA as compared to EDTA alone. Our data revealed that combined SPT-EDTA treated sera is more effective than ADS, and SC in reducing high background in SPI. Taken together, SPT-EDTA serum treatment prior to Luminex HLA Ab testing is cost-effective, our laboratory saves nearly 30% of the annual total cost for Ab testing and improved test turnaround time by two business days.

Analysis of Cross-sectional and Longitudinal HLA and Anti-viral Responses After COVID Infection in Renal Allograft Recipients: Differences and Correlates.

BACKGROUND: Characterization of anti-HLA versus anti-severe acute respiratory syndrome coronavirus 2 (anti-SARS-CoV-2) immune globulin isotypes in organ transplant recipients after coronavirus disease 2019 (COVID-19) infection has not been reported. We aimed to determine changes in anti-HLA antibodies in renal transplant patients with COVID-19 and compare the immunoglobulin and epitope-binding pattern versus anti-SARS-CoV-2 antibodies. METHODS: This is a cross-sectional study of 46 kidney transplant recipients including 21 with longitudinal sampling. Using a semi-quantitative multiplex assay, we determined immunoglobulin (Ig) M, IgA, IgG, and IgG1-2-3-4 antibodies against Class I and Class II HLA, and 5 SARS-CoV-2 epitopes including the nucleocapsid protein and multiple regions of the spike protein. RESULTS: Fourteen of 46 (30%) patients had donor-specific anti-HLA antibodies (donor-specific antibody [DSA]), 12 (26%) had non-DSA anti-HLA antibodies and 45 (98%) had anti-SARS-CoV-2 antibodies. Most DSAs targeted HLA-DQ (71%), with a dominant IgG isotype and IgG1 subtype prevalence (93%), and/or IgG3 (64%), followed by IgG2 (36%). Comparatively, there was a higher prevalence of IgA (85% versus 14%, P = 0.0001) and IgM (87%, versus 36%, P = 0.001) in the anti-SARS-CoV-2 antibody profile, when compared to DSAs, respectively. Anti-SARS-CoV-2 antibody profile was characterized by increased prevalence of IgM and IgA, when compared to DSAs. The median calculated panel reactive antibody before COVID-19 diagnosis (24%) tended to decrease after COVID-19 diagnosis (10%) but it was not statistically significant ( P = 0.1). CONCLUSIONS: Anti-HLA antibody strength and calculated panel reactive antibody in kidney transplant recipients after COVID-19 do not significantly increase after infection. Although the IgG isotype was the dominant form in both HLA and SARS-CoV-2 antigens, the alloimmune response had a low IgA pattern, whereas anti-SARS-CoV-2 antibodies were high IgA/IgM.

Complement-Binding Donor-Specific Anti-HLA Antibodies: Biomarker for Immunologic Risk Stratification in Pediatric Kidney Transplantation Recipients.

Antibody-mediated rejection is a common cause of early kidney allograft loss but the specifics of antibody measurement, therapies and endpoints have not been universally defined. In this retrospective study, we assessed the performance of risk stratification using systematic donor-specific antibody (DSA) monitoring. Included in the study were children who underwent kidney transplantation between January 1, 2010 and March 1, 2018 at Stanford, with at least 12-months follow-up. A total of 233 patients were included with a mean follow-up time of 45 (range, 9-108) months. Median age at transplant was 12.3 years, 46.8% were female, and 76% had a deceased donor transplant. Fifty-two (22%) formed C1q-binding de novo donor-specific antibodies (C1q-dnDSA). After a standardized augmented immunosuppressive protocol was implemented, C1q-dnDSA disappeared in 31 (58.5%). Graft failure occurred in 16 patients at a median of 54 (range, 5-83) months, of whom 14 formed dnDSA. The 14 patients who lost their graft due to rejection, all had persistent C1q-dnDSA. C1q-binding status improved the individual risk assessment, with persistent; C1q binding yielding the strongest independent association of graft failure (hazard ratio, 45.5; 95% confidence interval, 11.7-177.4). C1q-dnDSA is more useful than standard dnDSA as a noninvasive biomarker for identifying patients at the highest risk of graft failure.

Current HLA testing recommendations to support HCT.

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.

A new strategy for systematically classifying HLA alleles into serological specificities.

HLA serological specificities were defined by the reactivity of HLA molecules with sets of sera and monoclonal antibodies. Many recently identified alleles defined by molecular typing lack their serotype assignment. We surveyed the literature describing the correlation of the reactivity of serologic reagents with AA residues. 20 - 25 AA residues determining epitopes (DEP) that correlated with 82 WHO serologic specificities were identified for HLA class I loci. Thirteen DEP each located in the beta-1 domains that correlated with 24 WHO serologic specificities were identified for HLA-DRB1 and -DQB1 loci. The designation of possible HLA-DPB1, -DQA1, -DPA1, and additional serological specificities that result from epitopes defined by residues located at both -DQA1 and -DQB1 subunits were also examined. HATS software was developed for automated serotype assignments to HLA alleles in one of the three hierarchical matching criteria: (1) all DEP (FULL); (2) selected DEP specific to each serological specificity (SEROTYPE); (3) one AA mismatch with one or more SEROTYPES (INCOMPLETE). Results were validated by evaluating the alleles whose serotypes do not correspond to the first field of the allele name listed in the HLA dictionary. Additional 85 and 21 DEP patterns that do not correspond to any WHO serologic specificities for common HLA class I and DRB1 alleles were identified, respectively. A comprehensive antibody identification panel would allow for accurate unacceptable antigen listing and compatibility predictions in solid organ transplantation. We propose that antibody-screening panels should include all serologic specificities identified in this study.

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

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.

Severe delayed hypersensitivity reactions to IL-1 and IL-6 inhibitors link to common HLA-DRB1*15 alleles.

OBJECTIVES: Drug reaction with eosinophilia and systemic symptoms (DRESS) is a severe, delayed hypersensitivity reaction (DHR). We observed DRESS to inhibitors of interleukin 1 (IL-1) or IL-6 in a small group of patients with Still's disease with atypical lung disease. We sought to characterise features of patients with Still's disease with DRESS compared with drug-tolerant Still's controls. We analysed human leucocyte antigen (HLA) alleles for association to inhibitor-related DHR, including in a small Kawasaki disease (KD) cohort. METHODS: In a case/control study, we collected a multicentre series of patients with Still's disease with features of inhibitor-related DRESS (n=66) and drug-tolerant Still's controls (n=65). We retrospectively analysed clinical data from all Still's subjects and typed 94/131 for HLA. European Still's-DRESS cases were ancestry matched to International Childhood Arthritis Genetics Consortium paediatric Still's cases (n=550) and compared for HLA allele frequencies. HLA association also was analysed using Still's-DRESS cases (n=64) compared with drug-tolerant Still's controls (n=30). KD subjects (n=19) were similarly studied. RESULTS: Still's-DRESS features included eosinophilia (89%), AST-ALT elevation (75%) and non-evanescent rash (95%; 88% involving face). Macrophage activation syndrome during treatment was frequent in Still's-DRESS (64%) versus drug-tolerant Still's (3%; p=1.2×10-14). We found striking enrichment for HLA-DRB1*15 haplotypes in Still's-DRESS cases versus INCHARGE Still's controls (p=7.5×10-13) and versus self-identified, ancestry-matched Still's controls (p=6.3×10-10). In the KD cohort, DRB1*15:01 was present only in those with suspected anakinra reactions. CONCLUSIONS: DRESS-type reactions occur among patients treated with IL-1/IL-6 inhibitors and strongly associate with common HLA-DRB1*15 haplotypes. Consideration of preprescription HLA typing and vigilance for serious reactions to these drugs are warranted.

Delayed Kinetics of IgG, but Not IgA, Antispike Antibodies in Transplant Recipients following SARS-CoV-2 Infection.

BACKGROUND: Kidney transplant recipients are at increased risk of severe outcomes during COVID-19. Antibodies against the virus are thought to offer protection, but a thorough characterization of anti-SARS-CoV-2 immune globulin isotypes in kidney transplant recipients following SARS-CoV-2 infection has not been reported. METHODS: We performed a cross-sectional study of 49 kidney transplant recipients and 42 immunocompetent controls at early (≤14 days) or late (>14 days) time points after documented SARS-CoV-2 infection. Using a validated semiquantitative Luminex-based multiplex assay, we determined the abundances of IgM, IgG, IgG1-4, and IgA antibodies against five distinct viral epitopes. RESULTS: Kidney transplant recipients showed lower levels of total IgG antitrimeric spike (S), S1, S2, and receptor binding domain (RBD) but not nucleocapsid (NC) at early versus late time points after SARS-CoV-2 infection. Early levels of IgG antispike protein epitopes were also lower than in immunocompetent controls. Anti-SARS-CoV-2 antibodies were predominantly IgG1 and IgG3, with modest class switching to IgG2 or IgG4 in either cohort. Later levels of IgG antispike, S1, S2, RBD, and NC did not significantly differ between cohorts. There was no significant difference in the kinetics of either IgM or IgA antispike, S1, RBD, or S2 on the basis of timing after diagnosis or transplant status. CONCLUSIONS: Kidney transplant recipients mount early anti-SARS-CoV-2 IgA and IgM responses, whereas IgG responses are delayed compared with immunocompetent individuals. These findings might explain the poor outcomes in transplant recipients with COVID-19. PODCAST: This article contains a podcast at https://www.asn-online.org/media/podcast/JASN/2021_11_23_briggsgriffin112321.mp3.

Challenges for the standardized reporting of NGS HLA genotyping: Surveying gaps between clinical and research laboratories.

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.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility.

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from 'hitchhiking' alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

Next-Generation Sequencing Identifies Extended HLA Class I and II Haplotypes Associated With Early-Onset and Late-Onset Myasthenia Gravis in Italian, Norwegian, and Swedish Populations.

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.

High-Resolution Characterization of KIR Genes in a Large North American Cohort Reveals Novel Details of Structural and Sequence Diversity.

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.

High-resolution HLA allele and haplotype frequencies in several unrelated populations determined by next generation sequencing: 17th International HLA and Immunogenetics Workshop joint report.

The primary goal of the unrelated population HLA diversity (UPHD) component of the 17th International HLA and Immunogenetics Workshop was to characterize HLA alleles at maximum allelic-resolution in worldwide populations and re-evaluate patterns of HLA diversity across populations. The UPHD project included HLA genotype and sequence data, generated by various next-generation sequencing methods, from 4,240 individuals collated from 12 different countries. Population data included well-defined large datasets from the USA and smaller samples from Europe, Australia, and Western Asia. Allele and haplotype frequencies varied across populations from distant geographical regions. HLA genetic diversity estimated at 2- and 4-field allelic resolution revealed that diversity at the majority of loci, particularly for European-descent populations, was lower at the 2-field resolution. Several common alleles with identical protein sequences differing only by intronic substitutions were found in distinct haplotypes, revealing a more detailed characterization of linkage between variants within the HLA region. The examination of coding and non-coding nucleotide variation revealed many examples in which almost complete biunivocal relations between common alleles at different loci were observed resulting in higher linkage disequilibrium. Our reference data of HLA profiles characterized at maximum resolution from many populations is useful for anthropological studies, unrelated donor searches, transplantation, and disease association studies.

Specific Class I HLA Supertypes but Not HLA Zygosity or Expression Are Associated with Outcomes following HLA-Matched Allogeneic Hematopoietic Cell Transplant: HLA Supertypes Impact Allogeneic HCT Outcomes.

Maximizing the probability of antigen presentation to T cells through diversity in HLAs can enhance immune responsiveness and translate into improved clinical outcomes, as evidenced by the association of heterozygosity and supertypes at HLA class I loci with improved survival in patients with advanced solid tumors treated with immune checkpoint inhibitors. We investigated the impact of HLA heterozygosity, supertypes, and surface expression on outcomes in adult and pediatric patients with acute myeloid leukemia (AML), myelodysplastic syndrome, acute lymphoblastic leukemia, and non-Hodgkin lymphoma who underwent 8/8 HLA-matched, T cell replete, unrelated, allogeneic hematopoietic cell transplant (HCT) from 2000 to 2015 using patient data reported to the Center for International Blood and Marrow Transplant Research. HLA class I heterozygosity and HLA expression were not associated with overall survival, relapse, transplant-related mortality (TRM), disease-free survival (DFS), and acute graft-versus-host disease following HCT. The HLA-B62 supertype was associated with decreased TRM in the entire patient cohort (hazard ratio [HR], 0.79; 95% CI, 0.69 to 0.90; P = .00053). The HLA-B27 supertype was associated with worse DFS in patients with AML (HR = 1.21; 95% CI, 1.10 to 1.32; P = .00005). These findings suggest that the survival benefit of HLA heterozygosity seen in solid tumor patients receiving immune checkpoint inhibitors does not extend to patients undergoing allogeneic HCT. Certain HLA supertypes, however, are associated with TRM and DFS, suggesting that similarities in peptide presentation between supertype members play a role in these outcomes. Beyond implications for prognosis following HCT, these findings support the further investigation of these HLA supertypes and the specific immune peptides important for transplant outcomes.

Association of Human Leukocyte Antigens Class II Variants with Susceptibility to Hidradenitis Suppurativa in a Caucasian Spanish Population.

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.