Six-locus high resolution HLA haplotype frequencies derived from mixed-resolution DNA typing for the entire US donor registry

SNP-based HLA allele tagging, imputation and association with antiepileptic drug-induced cutaneous reactions in Hong Kong Han Chinese

Human leukocyte antigen (HLA) genes control the regulation of the human immune system and are involved in immune-related diseases. Population surveys on relationships between single nucleotide polymorphisms (SNP) and HLA alleles are essential to conduct genetic association between HLA variants and diseases. Samples were obtained from our in-house database for epilepsy genetics and pharmacogenetics research. Using 184 epilepsy patients with both genome-wide SNP array and HLA-A/B candidate gene sequencing data, we sought tagging SNPs that completely represent sixHLA risk alleles; in addition, a Hong Kong population-specific reference panel was constructed for SNP-based HLA imputation. The performance of our new panel was compared to a recent Han Chinese panel. Finally, genetic associations of HLA variants with mild skin rash were performed on the combined sample of 408 patients. Common SNPs rs2571375 and rs144295468 were found to successfully tag HLA risk alleles A*31:01 and B*13:01, respectively. HLA-B*15:02 can be predicted by rs144012689 with >95% sensitivity and specificity. The imputation reference panel for the Hong Kong population had comparable performance to the Han Chinese panel due to the large sample size for common HLA alleles, though it retained discordance for imputing rare alleles. No significant genetic associations were found between HLA genetic variants and mild skin rash induced by aromatic antiepileptic drugs. This study provides new information on the genetic structure of HLA regions in the Hong Kong population by identifying tagging SNPs and serving as a reference panel. Moreover, our comprehensive genetic analyses revealed no significant association between HLA alleles and mild skin rash in Hong Kong Han Chinese.

East Meets West-Impact of Ethnicity on Donor Match Rates in the Ezer Mizion Bone Marrow Donor Registry

HLA haplotype frequencies in a volunteer bone marrow donor registry should reflect the frequencies of potential transplant recipients served by that registry, a challenge in a country with diverse subethnicities of immigrants from Eastern and Western cultures, such as Israel. We evaluated the likelihood of finding suitable donors for hypothetical patients drawn from defined subethnicities in the Ezer Mizion Bone Marrow Donor Registry (EM BMDR) from donors both within and outside the registry now and during the coming decade. On average, bioinformatics modeling predicts that, given current donor recruitment trends, 6/6 high-resolution HLA match rates for Israelis, which currently stand at 40% to 55% for most subethnicities, will rise by up to 1% per year over the next decade. Subethnicities with historically lower rates of interethnic admixture are less likely to find matches outside of their designated group but will benefit from expansion of the registry, whereas ethnically directed drives will enhance matching rates for currently underrepresented subethnicities. Donor searches for the same cohort using a large extramural registry was of only slight benefit for most of the 19 EM BMDR subethnicities evaluated, confirming that local donor registries that reflect the ethnic diversity of the community being served are best equipped to serve the needs of their respective communities. Contemporary trends of an increasingly multiethnic admixture in Israel may impact the effect of ethnic profiling in assessing future match rates for EM BMDR.

Using high-resolution human leukocyte antigen typing of 11,423 randomized unrelated individuals to determine allelic varieties, deduce probable human leukocyte antigen haplotypes, and observe linkage disequilibria between human leukocyte antigen-B and-C and human leukocyte antigen-DRB1 and-DQB1 alleles in the Taiwanese Chinese population

OBJECTIVE

We report here the human leukocyte antigen (HLA) allelic variety and haplotype composition in a cohort of the Taiwanese Chinese population and their patterns of linkage disequilibria on HLA-B: HLA-C alleles and HLA-DRB1: HLA-DQB1 alleles at a high-resolution level.

MATERIALS AND METHODS

Peripheral whole blood from 11,423 Taiwanese Chinese unrelated individuals was collected in acid citrate dextrose. Genomic DNA was extracted using the QIAamp DNA Blood Mini Kit. The DNA material was subjected to HLA genotyping for HLA-A,-B,-C,-DRB1, and-DQB1 loci using a commercial polymerase chain reaction-sequence-based typing (PCR-SBT) kit, the SeCore® A/B/C/DRB1/DQB1 Locus Sequencing kit. High-resolution allelic sequencing was performed as previously described.

RESULTS

The number of individual HLA-B alleles detected was greater than the number of alleles recognized in the both the HLA-A and-DRB1 loci. Several novel alleles were discovered as a result of employing the SBT method and the high number of donors tested. In addition, we observed a genetic polymorphic feature of association between HLA-A and-B, HLA-B and-C, and HLA-DRB1 and-DQB1 alleles. Further, the homozygous haplotype frequencies of HLA-A and-B; HLA-A,-C, and-B; HLA-A,-C,-B, and-DRB1; and HLA-A,-C,-B,-DRB1, and-DQB1 in Taiwanese Chinese population are presented.

CONCLUSION

As increasing number of HLA alleles are being discovered, periodic HLA profile investigation in a given population is essential to recognize the HLA complexity in that population. Population study can also provide an up-to-date strategic plan for future needs in terms of compatibility measurement for HLA matching between transplant donors and patients.

Computational Approaches to Facilitate Epitope-Based HLA Matching in Solid Organ Transplantation

Epitope-based HLA matching has been emerged over the last few years as an improved method for HLA matching in solid organ transplantation. The epitope-based matching concept has been incorporated in both the PIRCHE-II and the HLAMatchmaker algorithm to find the most suitable donor for a recipient. For these algorithms, high-resolution HLA genotype data of both donor and recipient is required. Since high-resolution HLA genotype data is often not available, we developed a computational method which allows epitope-based HLA matching from serological split level HLA typing relying on HLA haplotype frequencies. To validate this method, we simulated a donor-recipient population for which PIRCHE-II and eplet values were calculated when using both high-resolution HLA genotype data and serological split level HLA typing. The majority of the serological split level HLA-determined ln(PIRCHE-II)/ln(eplet) values did not or only slightly deviate from the reference group of high-resolution HLA-determined ln(PIRCHE-II)/ln(eplet) values. This deviation was slightly increased when HLA-C or HLA-DQ was omitted from the input and was substantially decreased when using two-field resolution HLA genotype data of the recipient and serological split level HLA typing of the donor. Thus, our data suggest that our computational approach is a powerful tool to estimate PIRCHE-II/eplet values when high-resolution HLA genotype data is not available.

High resolution HLA haplotyping by imputation for a British population bioresource

This study aimed to establish the occurrence and frequency of HLA alleles and haplotypes for a healthy British Caucasian population bioresource from Oxfordshire. We present the results of imputation from HLA SNP genotyping data using SNP2HLA for 5553 individuals from Oxford Biobank, defining one- and two-field alleles together with amino acid polymorphisms. We show that this achieves a high level of accuracy with validation using sequence-specific primer amplification PCR. We define six- and eight-locus HLA haplotypes for this population by Bayesian methods implemented using PHASE. We determine patterns of linkage disequilibrium and recombination for these individuals involving classical HLA loci and show how analysis within a haplotype block structure may be more tractable for imputed data. Our findings contribute to knowledge of HLA diversity in healthy populations and further validate future large-scale use of HLA imputation as an informative approach in population bioresources.

Determination of HLA-A, -B, and -DRB1 Allele and Haplotype Frequencies in the Croatian Population Based on a Family Study

In the present study, HLA allele and haplotype frequencies among Croatian families were investigated to evaluate valuable information about HLA genotypes and to compare them with data from the Croatian Bone Marrow Donor Registry (CBMDR). A total of 350 families which have been typed for the purpose of HSCT were included. All individuals were tested using PCR-SSO or PCR-SSP methods for HLA-A, -B, and -DRB1 alleles. The HLA-A-B-DRB1 haplotypes were determined by segregation and directly counted. A total of 30 HLA-A, 54 HLA-B, and 38 HLA-DRB1 alleles and 716 different HLA-A-B-DRB1 haplotypes were identified. Of these, the three most frequent alleles at HLA-A, -B, and -DRB1 loci, respectively, were A*02:01 (30.39%), A*11:01 (13.37%), A*24:02 (10.91%); B*51:01 (12.48%), B*18:01 (8.35%), B*08:01 (8.06%); DRB1*03:01 (11.20%), DRB1*01:01 (9.84%), DRB1*16:01 (9.63%). The following HLA alleles were detected only once: A*02:09, A*24:03, A*24:04, A*24:07; B*07:04, B*15:07, B*15:08, B*39:04, B*39:10, B*39:24, B*40:04; DRB1*08:03, DRB1*11:06, DRB1*13:32, DRB1*14:05. Five most frequent haplotypes were: A*01:01-B*08:01-DRB1*03:01 (5.34%), A*02:01-B*18:01-DRB1*11:04 (1.57%), A*02:01-B*27:02-DRB1*16:01 (1.50%), A*02:01-B*27:05-DRB1*01:01 (1.42%) and A*02:01-B*44:02:01G-DRB1*16:01 (1.28%). The haplotype frequencies based on the family study were compared with the frequencies from CBMDR, and similar results were obtained for all except for the HLA-A*26:01-B*38:01-DRB1*04:02 haplotype. A significantly higher frequency (P = 0.0017) of this haplotype was observed among family individuals. Nine haplotypes were unique and data about their frequencies do not exist in current databases. The data obtained in this study could be useful for anthropology, transplantation and disease association studies.

Features of Autoimmune Hepatitis in Patients With Drug-induced Liver Injury

BACKGROUND & AIMS

Drug-induced liver injury (DILI) has features similar to those of other liver diseases including autoimmune hepatitis (AIH). We aimed to characterize the clinical and autoimmune features of liver injury caused by nitrofurantoin, minocycline, methyldopa, or hydralazine.

METHODS

We analyzed data from 88 cases of DILI attributed to nitrofurantoin, minocycline, methyldopa, or hydralazine included in the Drug-Induced Liver Injury Network prospective study from 2004 through 2014. Sera were collected from patients at baseline and follow-up examination and tested for levels of immunoglobulin G (IgG), antibodies to nuclear antigen (ANA), smooth muscle (SMA), and soluble liver antigen (SLA). An autoimmune score was derived on the basis of increases in levels of IgG, ANA, SMA, and SLA (assigned values of 0, 1+, or 2+). AIH-associated HLA-DRB1*03:01 and HLA-DRB1*04:01 allele frequencies were compared with those of the general population (controls).

RESULTS

Of the 88 cases, 80 were women (91%), 74% had hepatocellular injury, and 25% had severe injury. At the onset of DILI, 39% of cases had increased levels of IgG, 72% had increased levels of ANA, 60% had increased levels of SMA, and none had increases in SLA. A phenotype of autoimmunity (autoimmune score ≥2) was observed in 82% of cases attributed to nitrofurantoin and 73% of cases attributed to minocycline (73%) but only 55% of cases attributed to methyldopa and 43% of cases attributed to hydralazine (P = .16 for nitrofurantoin and minocycline vs methyldopa and hydralazine). We observed a decrease in numbers of serum samples positive for ANA (P = .01) or SMA (P < .001) and in autoimmune scores (P < .001) between DILI onset and follow-up. Similar percentages of patients with DILI had HLA-DRB1*03:01 (15%) and HLA-DRB1*04:01 (9%) as controls (12% and 9%, respectively).

CONCLUSIONS

In analysis of data from the DILIN prospective study, we found that most cases of DILI attributed to nitrofurantoin or minocycline and about half of cases that were due to methyldopa and hydralazine have a phenotype of autoimmunity similar to AIH. These features decrease with recovery of the injury and are not associated with the typical HLA alleles found in patients with idiopathic AIH.

HLA-DRB1*15:01 and HLA-DRB3*02:02 in PLA2R-Related Membranous Nephropathy

Idiopathic membranous nephropathy (MN) is associated with HLA; however, the HLA allele involved remains unknown. To identify the HLA risk alleles associated with phospholipase A2 receptor (PLA2R)-related MN in the Chinese population, we sequenced the entire MHC region in DNA samples from 99 patients with PLA2R-related MN, 50 patients with PLA2R-unrelated MN, and 100 healthy subjects. Two HLA risk alleles, HLA-DRB1*15:01 and HLA-DRB3*02:02, independently and strongly associated with an increased risk of PLA2R-related MN. After adjusting for HLA-DRB1*15:01 and HLA-DRB3*02:02, no other alleles showed significant association with PLA2R-related MN. A replication study in an independent cohort of 293 participants with PLA2R-related MN and 285 healthy controls validated these findings. In a joint analysis, a multivariate logistic regression model confirmed that HLA-DRB1*15:01 (odds ratio [OR], 24.9; 95% confidence interval [95% CI], 15.3 to 42.6; P=2.3×10^-35) and HLA-DRB3*02:02 (OR, 17.7; 95% CI, 11.0 to 30.3; P=8.0×10^-29) independently and strongly associated with PLA2R-related MN. As many as 98.7% of patients with PLA2R-related MN, compared with 43.9% of control subjects, carried at least one HLA risk allele. Subjects with either risk allele had higher odds of developing PLA2R-related MN than those without a risk allele (OR, 98.9; 95% CI, 44.4 to 281.7; P=2.5×10^-23). These HLA risk alleles also associated with the age at disease onset in patients with PLA2R-related MN. In conclusion, our findings provide clear evidence that the HLA-DRB1*15:01 and HLA-DRB3*02:02 alleles independently and strongly associate with PLA2R-related MN in the Chinese population.

HLA-DRB3*01:01 is a predictor of immunization against human platelet antigen-1a but not of the severity of fetal and neonatal alloimmune thrombocytopenia

BACKGROUND

Most cases of fetal and neonatal alloimmune thrombocytopenia (FNAIT) are caused by maternal alloantibodies against human platelet antigen-1a (HPA-1a). Alloimmunization mainly occurs in HPA-1a-negative mothers who are carriers of the HLA-DRB3*01:01 allele. Recently, it has been reported that the combined presence of HLA-DRB3*01:01 and HLA-DRB4*01:01P was associated with severity of FNAIT. We tested this hypothesis by analyzing a large cohort of cases and controls.

STUDY DESIGN AND METHODS

In total, 101 mothers with a history of FNAIT caused by anti-HPA-1a were investigated. HLA-DRB1, -DRB3, -DRB4, and -DRB5 genotypes were determined by Luminex technology. Haplotype frequencies were compared between cases and 100 controls. The platelet (PLT) counts of neonates and the incidence of intracranial hemorrhage (ICH) were compared between subgroups defined by genotype.

RESULTS

Of the HPA-1a-immunized mothers, 98% (99/101) carried at least one copy of HLA-DRB3*01:01. Carriage of HLA-DRB3*01:01 was significantly associated with immune response to HPA-1a (odds ratio, 92.3; 95% confidence interval, 26.9-317.1; p = 1.34 × 10^-12 ). No association between HLA-DRB3*01:01 and HLA-DRB4*01:01P alone or in combination with the PLT count of the newborns or the incidence of ICH was detected.

CONCLUSION

In contrast to HLA-DRB4*01:01P, the inheritance of HLA-DRB3*01:01 is strongly associated with the propensity for mounting a humoral immune response against fetal HPA-1a antigen. Neither a homozygous nor a compound heterozygous gene dose predicts the severity of the disease. Testing for the presence of HLA-DRB3*01:01 might be very useful in counseling women at risk of FNAIT due to anti-HPA-1a.

Fungal peptides from pneumonitis hypersensitivity etiologic agents are able to induce specific cellular immune response

PURPOSE

Hypersensitivity pneumonitis (HP) is an immunoallergic disease due to chronic exposure to high quantities of different microorganisms such as Mycobacterium immunogenum (Mi), a mycobacterium, and Lichtheimia corymbifera (Lc), a filamentous fungus. It has recently been demonstrated that the protein DLDH (dihydrolipoyl dehydrogenase), is common to these microorganisms. This study aimed to investigate the immune potential of overlapping peptide pools covering the MiDLDH and LcDLDH.

EXPERIMENTAL DESIGN

A selection of 34 peptides, from the MiDLDH and LcDLDH, able to interact with Major Histocompatibility Complex (MHC) 1 and MHC 2, was obtained using three different epitope prediction websites. By means of ELISPOT assays, we compared the frequency of Interferon gamma (IFNγ) secreting peripheral blood mononuclear cells (PBMC) after stimulation with overlapping peptide pools. Tests were performed using cells from 35 healthy blood donors.

RESULTS

One peptide pool containing five peptides from MiDLDH and able to interact with MHC 2 induced a marked IFNγ specific immune response (Pool F, p<0.001, Wilcoxon signed-rank test).

CONCLUSION

This study demonstrated that peptides from microorganisms involved in HP were able to induce a high IFNγ specific immune response after stimulation of PBMCs from healthy blood donors which could be useful to develop an effective prevention strategy.

HLA-A, -B, -C, and -DRB1 genotyping and haplotype frequencies for a Hong Kong Chinese population of 7595 individuals

HLA-A, -B, -C and -DRB1 gene and haplotype frequencies have been calculated from 7595 southern Chinese unrelated donors in a Hong Kong Bone Marrow Donor Registry. This is the first large-scale paper to report the distribution of A-C-B-DRB1 alleles in Hong Kong Chinese. This information is important for phylogenetic, comparative studies and estimating the optimal and cost-effective donor size and likelihood of obtaining appropriately matched donors for Chinese patients awaiting haematopoietic stem cell transplantation. The allele and haplotype data are available in the Allele Frequencies Net Database under the population name ''Hong Kong Chinese BMDR'' and the identifier (AFND003357).

HLA-A, -B and -DRB1 genotyping and haplotype frequencies of 3892 cord blood units in the Hong Kong Chinese Cord Blood Registry

HLA-A, -B and -DRB1 gene and haplotype frequencies have been calculated from 3892 southern Chinese unrelated cord blood units in a Hong Kong Cord Blood Registry. This is the first large-scale paper to report the distribution of A-B-DRB1 alleles in Hong Kong Chinese Cord Blood Units. This information is important for estimating the optimal and economically cost-effective donor size and likelihood of obtaining appropriately matched cord blood units for Chinese patients awaiting haematopoietic stem cell transplantation. The data are available in the Allele Frequencies Net Database under the population name ''Hong Kong Chinese Cord Blood Registry'' and the identifier (AFND003358).

Maternal HLA genotyping is not useful for predicting severity of fetal and neonatal alloimmune thrombocytopenia

Lack of reliable laboratory parameters is the main challenge in the management of fetal and neonatal alloimmune thrombocytopenia (FNAIT). Despite the long-known association between the HLA-DRB3*01:01 allele and human platelet antigen 1a (HPA-1a) alloimmunisation, maternal human leucocyte antigen (HLA) typing has been of little clinical value. Recently, other DRB3 allele variants have been suggested to predict the severity of FNAIT. In this nationwide population-based retrospective cohort study, we performed extensive HLA typing of 96 women, accounting for 87% of our cohort of 110 families with confirmed or possible HPA-1a-immunisation. The HLA type was compared with anti-HPA-1a levels, severity of neonatal disease and responsiveness to maternally administrated intravenous gammaglobulin (IVIG). HLA haplotypes were constructed to investigate further HLA associations. Despite significantly lower anti-HPA-1a levels in DRB3*01:01-negative women, the carrier status of this particular allele could not be used to confirm or rule out FNAIT in the absence of detectable antibodies. In the haplotype analysis, the DRB3*01:01 allele was the actual factor associated with FNAIT. No other HLA allele was shown to be of additional value as a predictor of severe FNAIT or non-responsiveness to IVIG treatment. Thus, HLA genotyping was not found useful in differentiating high- and low-risk pregnancies or in guiding antenatal treatment in affected families.

The effect of HLA allele and haplotype polymorphisms on donor matching in hematopoietic stem cell transplantation - Croatian experience

The knowledge of HLA characteristics of a patient's population helps to predict the probability of finding a MUD. The study included 170 transplanted patients for whom a search for a MUD in BMDW was performed and a sample of 4000 volunteer unrelated donors from the Croatian Bone Marrow Donor Registry (CBMDR). Patients and their MUDs were typed for HLA-A, -B, -C, -DRB1, and -DQB1 loci using PCR-SSO and PCR-SSP methods while donors were typed for HLA-A, -B, -C, and -DRB1 loci using the PCR-SSO method. A comparison of allele frequencies at tested HLA loci between patients and donors from CBMDR did not reveal significant differences. The majority of patients (117, 68.8%) had a 10/10 MUD, 45 (26.5%) patients had a 9/10 MUD and eight (4.7%) patients had an 8/10 MUD. The highest number of mismatches (MM) was present at HLA-DRB1 (19; 31.1%). The presence of DRB1*11 and DRB1*04 allelic groups among patients caused allelic MMs at HLA-DRB1 in most cases. The presence of an infrequent HLA-B∼C haplotype resulted in the HLA-C MM at antigen level in the majority of cases. The present study clarified HLA factors that cause difficulties in searching for a 10/10 MUD for Croatian patients.

High-resolution HLA A∼B∼DRB1 haplotype frequencies from the Ezer Mizion Bone Marrow Donor Registry in Israel

We have investigated HLA population alleles and haplotype frequencies for the ethnicities that comprise the contemporary population of Israel, using a large data set from the Ezer Mizion Bone Barrow Donor Registry. We genotyped 275,699 individuals at the HLA-A, -B and -DRB1 loci using HLA genotyping methods. HLA A∼B∼DRB1 haplotype frequencies were estimated from 19 sub-ethnic Jewish populations and other non-Jewish minorities using the maximum likelihood model, which accommodates typing ambiguities. We present overall and sub-ethnicity specific HLA diversity results of the registry, which will help guide a data-driven strategy for future registry expansion.

Racial and ethnic disparities in outcomes after heart transplantation: A systematic review of contributing factors and future directions to close the outcomes gap

The demographics of patients undergoing heart transplantation in the United States have shifted over the last 10 years, with an increasing number of racial and ethnic minorities undergoing heart transplant. Multiple studies have shown that survival of African American patients after heart transplantation is lower compared with other ethnic groups. We review the data supporting the presence of this outcome disparity and examine the multiple mechanisms that contribute. With an increasingly diverse population in the United States, knowledge of these disparities, their mechanisms, and ways to improve outcomes is essential.

Unrelated donor search prognostic score to support early HLA consultation and clinical decisions

A simple scoring system that can provide a quick search prognosis at the onset of an adult unrelated donor (URD) search could be a useful tool for transplant physicians. We aimed to determine whether patient human leukocyte Ag genotype frequency (GF) could be used as a surrogate measure of whether or not a potential 10/10 and/or 9/10 URD in the Be The Match Registry (BTMR) can be identified for the patient. GF was assigned on a training data set of 2410 patients that searched the BTMR using the reported ethnic group. A proportional odds model was used to correlate GF with defined search productivity categories as follows: 'Good' (>2 10/10), 'Fair' (1-2 10/10 or No 10/10 and >2 9/10) or 'Poor' (No 10/10 and <3 9/10). A second cohort (n=2411) was used to calculate the concordance by the ethnic group in all three categories. In addition, we validated against an independent cohort (n=1344) resolved as having a 10/10 or 9/10 matched URD. Across the ethnic groups, >90% of cases with 'Good' GF prognosis, 20-26% 'Fair' and <10% 'Poor' had a 10/10 URD. Although not a replacement for an actual URD search, GF offers a quick way for transplant physicians to get an indication of the likely search outcome.

A Rare HLA-DRB1*14:22-DQB1*04:01 Haplotype in a Kidney Donor: Implication in the Interpretation of Donor-Specific Antibody in Kidney Transplantation-A Case Report

BACKGROUND

The presence of donor-specific human leukocyte antigen antibodies (HLA-DSA) can be determined by performing Luminex assay with single-antigen beads. The single-antigen beads' panels cover the most frequent HLA alleles of the HLA-A, B, C, DRB1, DRB3/4/5, DQB, and DP loci, although the HLA typing for deceased donors often includes only HLA-A, B, and DR. Therefore, the information of haplotypic association between DRB1 and DQB1 is essential for the analysis of HLA-DSA, especially when HLA-DQ antibodies are identified in the patient's serum.

CASE REPORT

We report the finding of a rare HLA-DRB1*14:22-DQB1*04:01 haplotype in a Thai potential kidney donor. HLA class I and class II high-resolution typing were performed by a method of polymerase chain reaction with the use of sequence-specific primers. The HLA-A*24:02-C*04:06-B*13:01-DRB3*02:02-DRB1*14:22-DQA1*05:05-DQB1*04:01 haplotype in the kidney donor was confirmed by segregation analysis in the kidney donor's family. This rare haplotype was also identified in her father and the 2 of her offspring.

CONCLUSIONS

To our knowledge, this is the 1st report of the rare HLA-DRB1*14:22-DQB1*04:01 haplotype in Thai individuals. The information of the rare HLA-DR-DQ haplotypic association provides a caution for HLA laboratory personnel when analyzing HLA-DSA in a patient with HLA-DQ antibodies and the HLA-DQ typing of a deceased donor is unavailable.

Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease

The human major histocompatibility complex (MHC) region has been shown to be associated with numerous diseases. However, it remains a challenge to pinpoint the causal variants for these associations because of the extreme complexity of the region. We thus sequenced the entire 5-Mb MHC region in 20,635 individuals of Han Chinese ancestry (10,689 controls and 9,946 patients with psoriasis) and constructed a Han-MHC database that includes both variants and HLA gene typing results of high accuracy. We further identified multiple independent new susceptibility loci in HLA-C, HLA-B, HLA-DPB1 and BTNL2 and an intergenic variant, rs118179173, associated with psoriasis and confirmed the well-established risk allele HLA-C*06:02. We anticipate that our Han-MHC reference panel built by deep sequencing of a large number of samples will serve as a useful tool for investigating the role of the MHC region in a variety of diseases and thus advance understanding of the pathogenesis of these disorders.

Identification of a novel HLA-C*03 variant, C*03:303, by sequence-based typing in a Dai Chinese individual

HLA-C*03:303 has one nonsynonymous nucleotide change from HLA-C*03:04:01:01 in codon 194, exon 4 (GTC>GCC).

Charting improvements in US registry HLA typing ambiguity using a typing resolution score

Unrelated stem cell registries have been collecting HLA typing of volunteer bone marrow donors for over 25years. Donor selection for hematopoietic stem cell transplantation is based primarily on matching the alleles of donors and patients at five polymorphic HLA loci. As HLA typing technologies have continually advanced since the beginnings of stem cell transplantation, registries have accrued typings of varied HLA typing ambiguity. We present a new typing resolution score (TRS), based on the likelihood of self-match, that allows the systematic comparison of HLA typings across different methods, data sets and populations. We apply the TRS to chart improvement in HLA typing within the Be The Match Registry of the United States from the initiation of DNA-based HLA typing to the current state of high-resolution typing using next-generation sequencing technologies. In addition, we present a publicly available online tool for evaluation of any given HLA typing. This TRS objectively evaluates HLA typing methods and can help define standards for acceptable recruitment HLA typing.

Diversity in exon 5 of HLA-C(∗)04:01:01G is significant in anthropological studies

Polymorphisms in Human Leukocyte Antigen (HLA) class I genes are generally considered to be relevant only if they reside in exons 2 and 3 or if they affect the expression of the allele. HLA-C(∗)04:82 differs from the common HLA-C(∗)04:01:01 by having a 9 nucleotide, or 3 amino acid duplication, in exon 5. Having observed HLA-C(∗)04:82 in a New Zealand Maori stem cell patient, we have attempted to examine the prevalence of this allele in different ethnicities. Although our studies are in a limited number of patients and donors, they have revealed that, in the Pacific region, HLA-C(∗)04:82 appears to be the most common allele of the HLA-C(∗)04:01:01G group of alleles, notably in Filippinos and in Maori/Polynesians. In these populations this allele has characteristic HLA-ABCDRB1 haplotypes. Thus, our studies have shown that polymorphisms outside of the clinically important exons can be considered to be relevant in anthropological studies.

Completion of HLA protein sequences by automated homology-based nearest-neighbor extrapolation of HLA database sequences

The IMGT/HLA database contains every publicly available HLA sequence. However, most of these HLA protein sequences are restricted to the alpha-1/alpha-2 domain for HLA class-I and alpha-1/beta-1 domain for HLA class-II. Nevertheless, also polymorphism outside these domains may play a role in alloreactivity after transplantation. Several computational algorithms that aim for predicting alloreactivity after transplantation, HLAMatchmaker and the PIRCHE algorithm, require a major or the whole part of the HLA protein sequence as input for their prediction. In this study we describe an automated homology-based nearest-neighbor method to extrapolate incomplete HLA protein sequences. To get insight in the reliability of our extrapolation method, we performed a 10-fold cross-validation. The majority of the amino acid positions of the individual HLA class-I and -II proteins were correctly predicted. Eplets as defined by HLAMatchmaker were located both at correctly predicted as well as at lower reliably predicted amino acid positions. Moreover, six out of seven completely sequenced HLA alleles with previously unknown exons sequences were in agreement with the extrapolated amino acid sequences. In conclusion, incomplete HLA sequences can be completed by a homology-based nearest-neighbor principle. Since this method is automated, future submitted incomplete HLA sequences can easily be extrapolated.

Immunological evasion of immediate-early varicella zoster virus proteins

The varicella zoster virus (VZV) causes the childhood disease commonly known as chickenpox and can later in life reactivate as herpes zoster. The adaptive immune system is known to play an important role in suppressing VZV reactivation. A central aspect of this system is the presentation of VZV-derived peptides by the major histocompatibility complex (MHC) proteins. Here, we investigate if key VZV proteins have evolved their amino acid sequence to avoid presentation by MHC based on predictive models of MHC-peptide affinity. This study shows that the immediate-early proteins of all characterized VZV strains are profoundly depleted for high-affinity MHC-I-restricted epitopes. The same depletion can be found in its closest animal analog, the simian varicella virus. Further orthology analysis towards other herpes viruses suggests that the protein amino acid frequency is one of the primary drivers of targeted epitope depletion.

Genetic editing of HLA expression in hematopoietic stem cells to broaden their human application

Mismatch of human leukocyte antigens (HLA) adversely impacts the outcome of patients after allogeneic hematopoietic stem-cell transplantation (alloHSCT). This translates into the clinical requirement to timely identify suitable HLA-matched donors which in turn curtails the chances of recipients, especially those from a racial minority, to successfully undergo alloHSCT. We thus sought to broaden the existing pool of registered unrelated donors based on analysis that eliminating the expression of the HLA-A increases the chance for finding a donor matched at HLA-B, -C, and -DRB1 regardless of a patient’s race. Elimination of HLA-A expression in HSC was achieved using artificial zinc finger nucleases designed to target HLA-A alleles. Significantly, these engineered HSCs maintain their ability to engraft and reconstitute hematopoiesis in immunocompromised mice. This introduced loss of HLA-A expression decreases the need to recruit large number of donors to match with potential recipients and has particular importance for patients whose HLA repertoire is under-represented in the current donor pool. Furthermore, the genetic engineering of stem cells provides a translational approach to HLA-match a limited number of third-party donors with a wide number of recipients.

Alleles of HLA-DRB1*04 Associated with Pulmonary Tuberculosis in Amazon Brazilian Population

Immunogenetic host factors are associated with susceptibility or protection to tuberculosis (TB). Strong associations of HLA class II genes with TB are reported. We analyzed the HLA-DRB1*04 alleles to identify subtypes associated with pulmonary TB and their interaction with risk factors such as alcohol, smoking, and gender in 316 pulmonary TB patients and 306 healthy individuals from the Brazilian Amazon. The HLA-DRB1*04 was prevalent in patients with pulmonary TB (p<0.0001; OR = 2.94; 95% CI = 2.12 to 4.08). Direct nucleotide sequencing of DRB1 exon 2 identified nine subtypes of HLA-DRB1*04. The subtype HLA-DRB1*04:11:01 (p = 0.0019; OR = 2.23; 95% CI = 1.34 to 3.70) was associated with susceptibility to pulmonary TB while DRB1*04:07:01 (p<0.0001; OR = 0.02; 95% CI = 0.001 to 0.33) to protection. Notably, the interaction between alcohol and HLA-DRB1*04:11:01 increased the risk for developing pulmonary TB (p = 0.0001; OR = 51.3; 95% CI = 6.81 to 386). Multibacillary pulmonary TB, the clinical presentation of disease transmission, was strongly associated with interaction to alcohol (p = 0.0026; OR = 11.1; 95% CI = 3.99 to 30.9), HLA-DRB1*04:11:01 (p = 0.0442; OR = 2.01; 95% CI = 1.03 to 3.93) and DRB1*04:92 (p = 0.0112; OR = 8.62; 95% CI = 1.63 to 45.5). These results show that HLA-DRB1*04 are associated with pulmonary TB. Interestingly, three subtypes, DRB1*04:07:01, DRB1*04:11:01 and DRB1*04:92 of the HLA-DRB1*04 could be potential immunogenetic markers that may help to explain mechanisms involved in disease development.

Determining performance characteristics of an NGS-based HLA typing method for clinical applications

This study presents performance specifications of an in-house developed human leukocyte antigen (HLA) typing assay using next-generation sequencing (NGS) on the Illumina MiSeq platform. A total of 253 samples, previously characterized for HLA-A, -B, -C, -DRB1 and -DQB1 were included in this study, which were typed at high-resolution using a combination of Sanger sequencing, sequence-specific primer (SSP) and sequence-specific oligonucleotide probe (SSOP) technologies and recorded at the two-field level. Samples were selected with alleles that cover a high percentage of HLA specificities in each of five different race/ethnic groups: European, African-American, Asian Pacific Islander, Hispanic and Native American. Sequencing data were analyzed by two software programs, Omixon's target and GenDx's NGSengine. A number of metrics including allele balance, sensitivity, specificity, precision, accuracy and remaining ambiguity were assessed. Data analyzed by the two software systems are shown independently. The majority of alleles were identical in the exonic sequences (third field) with both programs for HLA-A, -B, -C and -DQB1 in 97.7% of allele determinations. Among the remaining discrepant genotype calls at least one of the analysis programs agreed with the reference typing. Upon additional manual analysis 100% of the 2530 alleles were concordant with the reference HLA genotypes; the remaining ambiguities did not exceed 0.8%. The results demonstrate the feasibility and significant benefit of HLA typing by NGS as this technology is highly accurate, eliminates virtually all ambiguities, provides complete sequencing information for the length of the HLA gene and forms the basis for utilizing a single methodology for HLA typing in the immunogenetics labs.

Identification of epitopes on HLA-DRB alleles reacting with antibodies in sera from women sensitized during pregnancy

This report describes a HLAMatchmaker-based antibody analysis of post-pregnancy sera with antibodies against child-specific HLA-DR epitopes. These sera were reactive in IgG-binding assays with single allele bead (SAB) panels on a Luminex platform. The antibody specificity analysis focused on DRB epitopes that have been recorded in the International HLA Epitope Registry (http://www.epregistry.com.br) as experimentally verified with informative antibodies but we also considered other eplets that predict potential epitopes. The SAB panel has in several instances two or more alleles corresponding to the same serologically defined DR antigen and we selected six sera were with different reactivity patterns with DR1, DR4, DR13 and/or DR52 alleles. We demonstrate here how amino acid differences between these alleles can provide useful information in the determination of new epitope specificities of antibodies in these sera. Eight newly antibody-verified epitopes were identified including three that correspond to eplets paired with self-residue configurations. Epitope specificity information appears to be useful in the prediction of mismatch acceptability of non-SAB alleles within serological DR antigen groups.

Concordance of next generation sequence-based and sequence specific oligonucleotide probe-based HLA-DRB1 genotyping

Next generation sequencing (NGS) of clonally amplified DNA, using Roche 454 technology, was used to genotype HLA-DRB1, DRB3, DRB4, and DRB5 loci (exon 2 only) from a set of 993 samples from newborns with maternally-reported African American ancestry. DRB1 exon 2 was genotyped previously on the same sample set using sequence-specific oligonucleotide probe (SSOP) technology. Comparison of the genotype calls from both methods indicated concordance of 92.3%. Some discordance was expected due to the higher resolution of NGS data, compared to SSOP data. This resulted from selection of the incorrect allele from the ambiguity string produced by SSOP genotyping. Of 76 discordant genotypes, only three were due to resolution of ambiguity with the NGS method. The low percent of changes due to the increased resolution of the NGS method instills confidence in the overall value of previous data genotyped with moderate resolution methods, i.e., the vast majority of alleles present in a population are those that are detectable at moderate resolution. The remaining 73 discordant genotypes resulted from preventable errors in sample handling, data interpretation, and data entry. These results underscore the potential for error that can result from factors such as low quality genomic DNA, manual data entry, and interpretation of marginal genotyping results. Optimization of genomic DNA quality, automation of genotyping steps wherever possible, and use of the highest resolution technology available can lead to dramatic improvements in HLA genotype data quality. NGS-based methodology generated data of superior quality and accuracy compared to the SSOP system.

HLA Haplotype Validator for quality assessments of HLA typing

HLA alleles are observed in specific haplotypes, due to Linkage Disequilibrium (LD) between particular alleles. Haplotype frequencies for alleles in strong LD have been established for specific ethnic groups and racial categories. Application of high-resolution HLA typing using Next Generation Sequencing (NGS) is becoming a common practice in research and clinical laboratory settings. HLA typing errors using NGS occasionally occur due to allelic sequence imbalance or misalignment. Manual inspection of HLA genotypes is labor intensive and requires an in-depth knowledge of HLA alleles and haplotypes. We developed the "HLA Haplotype Validator (HLAHapV)" software, which inspects an HLA genotype for both the presence of common and well-documented alleles and observed haplotypes. The software also reports warnings when rare alleles, or alleles that do not belong to recognized haplotypes, are found. The software validates observable haplotypes in genotype data, providing increased confidence regarding the accuracy of the HLA typing, and thus reducing the effort involved in correcting potential HLA typing errors. The HLAHapV software is a powerful tool for quality control of HLA genotypes prior to the application of downstream analyses. We demonstrate the use of the HLAHapV software for identifying unusual haplotypes, which can lead to finding potential HLA typing errors.

Narcolepsy in African Americans

STUDY OBJECTIVES

Although narcolepsy affects 0.02-0.05% of individuals in various ethnic groups, clinical presentation in different ethnicities has never been fully characterized. Our goal was to study phenotypic expression across ethnicities in the United States.

DESIGN/SETTING

Cases of narcolepsy from 1992 to 2013 were identified from searches of the Stanford Center for Narcolepsy Research database. International Classification of Sleep Disorders, Third Edition diagnosis criteria for type 1 and type 2 narcolepsy were used for inclusion, but subjects were separated as with and without cataplexy for the purpose of data presentation. Information extracted included demographics, ethnicity and clinical data, HLA-DQB1*06:02, polysomnography (PSG), multiple sleep latency test (MSLT) data, and cerebrospinal fluid (CSF) hypocretin-1 level.

PATIENTS

182 African-Americans, 839 Caucasians, 35 Asians, and 41 Latinos with narcolepsy.

RESULTS

Sex ratio, PSG, and MSLT findings did not differ across ethnicities. Epworth Sleepiness Scale (ESS) score was higher and age of onset of sleepiness earlier in African Americans compared with other ethnicities. HLA-DQB1*06:02 positivity was higher in African Americans (91.0%) versus others (76.6% in Caucasians, 80.0% in Asians, and 65.0% in Latinos). CSF hypocretin-1 level, obtained in 222 patients, was more frequently low (≤ 110 pg/ml) in African Americans (93.9%) versus Caucasians (61.5%), Asians (85.7%) and Latinos (75.0%). In subjects with low CSF hypocretin-1, African Americans (28.3%) were 4.5 fold more likely to be without cataplexy when compared with Caucasians (8.1%).

CONCLUSIONS

Narcolepsy in African Americans is characterized by earlier symptom onset, higher Epworth Sleepiness Scale score, higher HLA-DQB1*06:02 positivity, and low cerebrospinal fluid hypocretin-1 level in the absence of cataplexy. In African Americans, more subjects without cataplexy have type 1 narcolepsy.

FluKB: A Knowledge-Based System for Influenza Vaccine Target Discovery and Analysis of the Immunological Properties of Influenza Viruses

FluKB is a knowledge-based system focusing on data and analytical tools for influenza vaccine discovery. The main goal of FluKB is to provide access to curated influenza sequence and epitope data and enhance the analysis of influenza sequence diversity and the analysis of targets of immune responses. FluKB consists of more than 400,000 influenza protein sequences, known epitope data (357 verified T-cell epitopes, 685 HLA binders, and 16 naturally processed MHC ligands), and a collection of 28 influenza antibodies and their structurally defined B-cell epitopes. FluKB was built using a modular framework allowing the implementation of analytical workflows and includes standard search tools, such as keyword search and sequence similarity queries, as well as advanced tools for the analysis of sequence variability. The advanced analytical tools for vaccine discovery include visual mapping of T- and B-cell vaccine targets and assessment of neutralizing antibody coverage. FluKB supports the discovery of vaccine targets and the analysis of viral diversity and its implications for vaccine discovery as well as potential T-cell breadth and antibody cross neutralization involving multiple strains. FluKB is representation of a new generation of databases that integrates data, analytical tools, and analytical workflows that enable comprehensive analysis and automatic generation of analysis reports.

Asymmetric linkage disequilibrium: Tools for assessing multiallelic LD

Standard measures of linkage disequilibrium (LD) provide an incomplete description of the correlation between two loci. Recently, Thomson and Single (2014) described a new asymmetric pair of LD measures (ALD) that give a more complete description of LD. The ALD measures are symmetric and equivalent to the correlation coefficient r when both loci are bi-allelic. When the numbers of alleles at the two loci differ, the ALD measures capture this asymmetry and provide additional detail about the LD structure. In disease association studies the ALD measures are useful for identifying additional disease genes in a genetic region, by conditioning on known effects. In evolutionary genetic studies ALD measures provide insight into selection acting on individual amino acids of specific genes, or other loci in high LD (see Thomson and Single (2014) for these examples). Here we describe new software for computing and visualizing ALD. We demonstrate the utility of this software using haplotype frequency data from the National Marrow Donor Program (NMDP). This enhances our understanding of LD patterns in the NMDP data by quantifying the degree to which LD is asymmetric and also quantifies this effect for individual alleles.

Immunogenetics of type 1 diabetes: A comprehensive review

Type 1 diabetes (T1D) results from the autoimmune destruction of insulin-producing beta cells in the pancreas. Prevention of T1D will require the ability to detect and modulate the autoimmune process before the clinical onset of disease. Genetic screening is a logical first step in identification of future patients to test prevention strategies. Susceptibility to T1D includes a strong genetic component, with the strongest risk attributable to genes that encode the classical Human Leukocyte Antigens (HLA). Other genetic loci, both immune and non-immune genes, contribute to T1D risk; however, the results of decades of small and large genetic linkage and association studies show clearly that the HLA genes confer the most disease risk and protection and can be used as part of a prediction strategy for T1D. Current predictive genetic models, based on HLA and other susceptibility loci, are effective in identifying the highest-risk individuals in populations of European descent. These models generally include screening for the HLA haplotypes "DR3" and "DR4." However, genetic variation among racial and ethnic groups reduces the predictive value of current models that are based on low resolution HLA genotyping. Not all DR3 and DR4 haplotypes are high T1D risk; some versions, rare in Europeans but high frequency in other populations, are even T1D protective. More information is needed to create predictive models for non-European populations. Comparative studies among different populations are needed to complete the knowledge base for the genetics of T1D risk to enable the eventual development of screening and intervention strategies applicable to all individuals, tailored to their individual genetic background. This review summarizes the current understanding of the genetic basis of T1D susceptibility, focusing on genes of the immune system, with particular emphasis on the HLA genes.

Distribution of HLA haplotypes across Japanese Archipelago: similarity, difference and admixture

The human leukocyte antigen (HLA) region is the most polymorphic region in the human genome. The polymorphic nature of the HLA region is thought to have been shaped from balancing selection. The complex migration events during the Out-of-Africa expansion have influenced geographic patterns of HLA allele frequencies and diversities across present-day human populations. Differences in the HLA allele frequency may contribute geographic differences in the susceptibility to many diseases, such as infectious, autoimmune and metabolic diseases. Here we briefly reviewed characteristics of frequency distribution of HLA alleles and haplotypes in Japanese population. A large part of HLA alleles and haplotypes that are common in Japanese are shared with neighboring Asian populations. The differentiations in HLA alleles and haplotypes across Japanese regional populations may provide clues to model for peopling of Japanese Archipelago and for design of genetic association studies. Finally, we introduce recent topics that new HLA alleles derived from ancient admixtures with Neanderthals and Denisovans are thought to have played an important role in the adaptation of modern humans to local pathogens during Out-of-Africa expansion.

Proceedings: human leukocyte antigen haplo-homozygous induced pluripotent stem cell haplobank modeled after the california population: evaluating matching in a multiethnic and admixed population

The development of a California-based induced pluripotent stem cell (iPSC) bank based on human leukocyte antigen (HLA) haplotype matching represents a significant challenge and a valuable opportunity for the advancement of regenerative medicine. However, previously published models of iPSC banks have neither addressed the admixed nature of populations like that of California nor evaluated the benefit to the population as a whole. We developed a new model for evaluating an iPSC haplobank based on demographic and immunogenetic characteristics reflecting California. The model evaluates haplolines or cell lines from donors homozygous for a single HLA-A, HLA-B, HLA-DRB1 haplotype. We generated estimates of the percentage of the population matched under various combinations of haplolines derived from six ancestries (black/African American, American Indian, Asian/Pacific Islander, Hispanic, and white/not Hispanic) and data available from the U.S. Census Bureau, the California Institute for Regenerative Medicine, and the National Marrow Donor Program. The model included both cis (haplotype-level) and trans (genotype-level) matching between a modeled iPSC haplobank and the recipient population following resampling simulations. We showed that serving a majority (>50%) of a simulated California population through cis matching would require the creation, redundant storage, and maintenance of almost 207 different haplolines representing the top 60 most frequent haplotypes from each ancestry group. Allowances for trans matching reduced the haplobank to fewer than 141 haplolines found among the top 40 most frequent haplotypes. Finally, we showed that a model optimized, custom haplobank was able to serve a majority of the California population with fewer than 80 haplolines.

Transposon-mediated death of an ancestral A-23-like allele: evolution of TCR-positioning motifs in the HLA-A lineage

HLA-A alleles are characterized by tandem arginine and histidine/arginine motifs (i.e., R65 and H151R motifs) present on the α1- and α2-helix, respectively. In crystallographic structures, α/β T-cell receptors (TCR) contact both motifs and appear to be geometrically positioned for alloreactivity. Herein, bioinformatics of "dual-motif" MHC A-like alleles were investigated across phylogeny. While A-like alleles with the R65 motif are widespread, the H151R motif has segregated out of most species. Surprisingly, an uncharacterized orf in tarsiers (Loc-103275158) encodes R151 within a truncated A-23-like gene, which is in frame with short footprints of Tc5 and Tigger transposons (TE); the extant tarsier A-23 allele is totally missing exon-3 and part of exon-4; together, suggesting TE-mediated inactivation of an intact/ancestral A-23 allele. Since the only other (non-human) dual-motif A-like alleles are in gorilla, chimpanzee, and the Florida manatee, we speculate that dual-motif A alleles first emerged in the Afrotherian lineage and reappeared during the evolution of higher primates.

Power laws for heavy-tailed distributions: modeling allele and haplotype diversity for the national marrow donor program

Measures of allele and haplotype diversity, which are fundamental properties in population genetics, often follow heavy tailed distributions. These measures are of particular interest in the field of hematopoietic stem cell transplant (HSCT). Donor/Recipient suitability for HSCT is determined by Human Leukocyte Antigen (HLA) similarity. Match predictions rely upon a precise description of HLA diversity, yet classical estimates are inaccurate given the heavy-tailed nature of the distribution. This directly affects HSCT matching and diversity measures in broader fields such as species richness. We, therefore, have developed a power-law based estimator to measure allele and haplotype diversity that accommodates heavy tails using the concepts of regular variation and occupancy distributions. Application of our estimator to 6.59 million donors in the Be The Match Registry revealed that haplotypes follow a heavy tail distribution across all ethnicities: for example, 44.65% of the European American haplotypes are represented by only 1 individual. Indeed, our discovery rate of all U.S. European American haplotypes is estimated at 23.45% based upon sampling 3.97% of the population, leaving a large number of unobserved haplotypes. Population coverage, however, is much higher at 99.4% given that 90% of European Americans carry one of the 4.5% most frequent haplotypes. Alleles were found to be less diverse suggesting the current registry represents most alleles in the population. Thus, for HSCT registries, haplotype discovery will remain high with continued recruitment to a very deep level of sampling, but population coverage will not. Finally, we compared the convergence of our power-law versus classical diversity estimators such as Capture recapture, Chao, ACE and Jackknife methods. When fit to the haplotype data, our estimator displayed favorable properties in terms of convergence (with respect to sampling depth) and accuracy (with respect to diversity estimates). This suggests that power-law based estimators offer a valid alternative to classical diversity estimators and may have broad applicability in the field of population genetics.

The distribution of haplotypes and species tend to be heavy tailed. The heavy tail is expected from theoretical considerations and is observed in most populations. Accurate measures of diversity are difficult to achieve given that a limited number of common haplotypes represent the majority of the population, whereas the major contributor to haplotype diversity comes from unique haplotypes that are “rare” and present in only a fraction of the population. A major issue for unrelated HSCT donor registries is estimating population coverage with respect to servicing the public need. We here use a power-law methodology that accommodates heavy-tails to estimate both the population coverage by ethnicity in the US and the genetic diversity of alleles and haplotypes. For the European American population, which has the deepest sampling amongst ethnicities, we show that registry population coverage is better than 99%, but the diversity of this sample only represents 40% of the unique haplotypes expected to be found in the population. Population coverage for other ethnicities was poorer and ranged down to 92% as was the case for Native Americans that had the worst coverage. We further show that the formalism developed here produces better estimates of the population properties than existing methods.

HLA class II diversity in HIV-1 uninfected individuals from the placebo arm of the RV144 Thai vaccine efficacy trial

The RV144 HIV vaccine trial in Thailand elicited antibody responses to the envelope of HIV-1, which correlated significantly with the risk of HIV-1 acquisition. Human leukocyte antigen (HLA) class II molecules are essential in antigen presentation to CD4 T cells for activation of B cells to produce antibodies. We genotyped the classical HLA-DRB1, DQB1, and DPB1 genes in 450 individuals from the placebo arm of the RV144 study to determine the background allele and haplotype frequencies of these genes in this cohort. High-resolution 4 and 6-digit class II HLA typing data was generated using sequencing-based methods. The observed diversity for the HLA loci was 33 HLA-DRB1, 15 HLA-DQB1, and 26 HLA-DPB1 alleles. Common alleles with frequencies greater than 10% were DRB1*07:01, DRB1*09:01, DRB1*12:02, DRB1*15:02, DQB1*02:01/02, DQB1*03:01, DQB1*03:03, DQB1*05:01, DQB1*05:02, DPB1*04:01:01, DPB1*05:01:01, and DPB1*13:01:01. We identified 28 rare alleles with frequencies of less than 1% in the Thai individuals. Ambiguity for HLA-DPB1*28:01 in exon 2 was resolved to DPB1*296:01 by next-generation sequencing of all exons. Multi-locus haplotypes including HLA class I and II loci were reported in this study. This is the first comprehensive report of allele and haplotype frequencies of all three HLA class II genes from a Thai population. A high-resolution genotyping method such as next-generation sequencing avoids missing rare alleles and resolves ambiguous calls. The HLA class II genotyping data generated in this study will be beneficial not only for future disease association/vaccine efficacy studies related to the RV144 study, but also for similar studies in other diseases in the Thai population, as well as population genetics and transplantation studies.