Hardy-Weinberg testing for HLA class II (DRB1, DQA1, DQB1, and DPB1) loci in 26 human ethnic groups

Quantifying uncertainty of molecular mismatch introduced by mislabeled ancestry using haplotype-based HLA genotype imputation

Introduction

Modern histocompatibility algorithms depend on the comparison and analysis of high-resolution HLA protein sequences and structures, especially when considering epitope-based algorithms, which aim to model the interactions involved in antibody or T cell binding. HLA genotype imputation can be performed in the cases where only low/intermediate-resolution HLA genotype is available or if specific loci are missing, and by providing an individuals' race/ethnicity/ancestry information, imputation results can be more accurate. This study assesses the effect of imputing high-resolution genotypes on molecular mismatch scores under a variety of ancestry assumptions.

Methods

We compared molecular matching scores from "ground-truth" high-resolution genotypes against scores from genotypes which are imputed from low-resolution genotypes. Analysis was focused on a simulated patient-donor dataset and confirmed using two real-world datasets, and deviations were aggregated based on various ancestry assumptions.

Results

We observed that using multiple imputation generally results in lower error in molecular matching scores compared to single imputation, and that using the correct ancestry assumptions can reduce error introduced during imputation.

Discussion

We conclude that for epitope analysis, imputation is a valuable and low-risk strategy, as long as care is taken regarding epitope analysis context, ancestry assumptions, and (multiple) imputation strategy.

A high-resolution HLA reference panel capturing global population diversity enables multi-ancestry fine-mapping in HIV host response

Fine-mapping to plausible causal variation may be more effective in multi-ancestry cohorts, particularly in the MHC, which has population-specific structure. To enable such studies, we constructed a large (n = 21,546) HLA reference panel spanning five global populations based on whole-genome sequences. Despite population-specific long-range haplotypes, we demonstrated accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT) populations). Applying HLA imputation to genome-wide association study data for HIV-1 viral load in three populations (EUR, AA and LAT), we obviated effects of previously reported associations from population-specific HIV studies and discovered a novel association at position 156 in HLA-B. We pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide-binding groove, explaining 12.9% of trait variance.

Clinical features, biochemistry, and HLA-DRB1 status in youth-onset type 1 diabetes in Sudan

OBJECTIVE

To further understand clinical and biochemical features, and HLA-DRB1 genotypes, in new cases of diabetes in Sudanese children and adolescents.

RESEARCH DESIGN AND METHODS

Demographic characteristics, clinical information, and biochemical parameters (blood glucose, HbA1c, C-peptide, autoantibodies against glutamic acid decarboxylase 65 [GADA] and insulinoma-associated protein-2 [IA-2A], and HLA-DRB1) were assessed in 99 individuals <18 years, recently (<18 months) clinically diagnosed with T1D. HLA-DRB1 genotypes for 56 of these Arab individuals with T1D were compared to a mixed control group of 198 healthy Arab (75%) and African (25%) individuals without T1D.

RESULTS

Mean ± SD age at diagnosis was 10.1 ± 4.3 years (range 0.7-17.6 years) with mode at 9-12 years. A female preponderance was observed. Fifty-two individuals (55.3%) presented in diabetic ketoacidosis (DKA). Mean ± SD serum fasting C-peptide values were 0.22 ± 0.25 nmol/L (0.66±0.74 ng/ml). 31.3% were autoantibody negative, 53.4% were GADA positive, 27.2% were IA-2A positive, with 12.1% positive for both autoantibodies. Association analysis compared to 198 controls of similar ethnic origin revealed strong locus association with HLA-DRB1 (p < 2.4 × 10-14 ). Five HLA-DRB1 alleles exhibited significant T1D association: three alleles (DRB1*03:01, DRB1*04:02, and DRB1*04:05) were positively associated, while three (DRB1*10:01, DRB1*15:02, and DRB1*15:03) were protective. DRB1*03:01 had the strongest association (odds ratio = 5.04, p = 1.7 × 10-10 ).

CONCLUSIONS

Young Sudanese individuals with T1D generally have similar characteristics to reported European-origin T1D populations. However, they have higher rates of DKA and slightly lower autoantibody rates than reported European-origin populations, and a particularly strong association with HLA-DRB1*03:01.

HLA-DRB1* and DQB1* allele and haplotype diversity in eight tribal populations: Global affinities and genetic basis of diseases in South India

The distribution of HLA class-II DRB1* and DQB1* alleles/ haplotypes were studied in 438 individuals of 8 Dravidian tribal groups inhabiting the Western Ghats, south India. The HLA typing was performed by PCR-SSP method. In order to identify the 5-locus Ancestral Extended Haplotypes (AEH), the alleles of HLA-A, -B and -C loci were typed for DNAs with predominant 2-locus haplotypes. The analyses have revealed allele HLA-DRB1*15 as the most predominant allele (Lowest / Highest range: Urali, 14.81 / Malasar, 48.94), followed by the alleles DRB1*10 (Katunayakan, 1.85 / Paliyan, 48.21), DRB1*14 (Paliyan 4.46 / Katunayakan, 40.74), DRB1*12 (Mannan, 1.64 / Katunayakan, 20.37) and DRB1*03 (Mannan, 1.64 / Urali, 29.63). The most frequent DQB1* alleles were DQB1*02 (Paliyan 3.57 / Urali, 23.15), DQB1*05 (Katunayakan, 27.77 / Paliyan 84.82) and DQB1*06 (Malasar, 8.51 / Kuruman, 33.51). The most predominant two-locus haplotypes observed were DRB1*15-DQB1*05, DRB1*10-DQB1*05, DRB1*15-DQB1*06 and DRB1*04-DQB1*05. The present study of HLA immunogenetics of south Indian tribes have revealed the presence of globally shared two and 5-locus haplotypes. Many of these haplotypes were implicated in a number of diseases in south India. We observed the presence of ancestral extended haplotypes (AEHs), hitherto not reported in Indian populations such as, A*68-B*35-C*02-DRB1*15:01-DQB1*05:01, A*24-B*57-C*06-DRB1*04:01-DQB1*05:01 and A*24-B*35-C*02-DRB1*15:01-DQB1*05:02. The dendrogram based phylogenetic analyses have revealed the Caucasian affinity of Urali, palaeo-Mediterranean and Indo-European affinity of Malasar tribes. The presence of globally shared susceptible and protective haplotypes reiterated the mosaic immunogenetic fabric of south Indian tribes.

The immunogenetic diversity of the HLA system in Mexico correlates with underlying population genetic structure

We studied HLA class I (HLA-A, -B) and class II (HLA-DRB1, -DQB1) allele groups and alleles by PCR-SSP based typing in a total of 15,318 mixed ancestry Mexicans from all the states of the country divided into 78 sample sets, providing information regarding allelic and haplotypic frequencies and their linkage disequilibrium, as well as admixture estimates and genetic substructure. We identified the presence of 4268 unique HLA extended haplotypes across Mexico and find that the ten most frequent (HF > 1%) HLA haplotypes with significant linkage disequilibrium (Δ'≥0.1) in Mexico (accounting for 20% of the haplotypic diversity of the country) are of primarily Native American ancestry (A*02~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*08~DQB1*04, A*68~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*14~DQB1*03:01, A*24~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*04~DQB1*03:02, A*02~B*40:02~DRB1*04~DQB1*03:02, A*68~B*35~DRB1*04~DQB1*03:02, A*02~B*15:01~DRB1*04~DQB1*03:02). Admixture estimates obtained by a maximum likelihood method using HLA-A/-B/-DRB1 as genetic estimators revealed that the main genetic components in Mexico as a whole are Native American (ranging from 37.8% in the northern part of the country to 81.5% in the southeastern region) and European (ranging from 11.5% in the southeast to 62.6% in northern Mexico). African admixture ranged from 0.0 to 12.7% not following any specific pattern. We were able to detect three major immunogenetic clusters correlating with genetic diversity and differential admixture within Mexico: North, Central and Southeast, which is in accordance with previous reports using genome-wide data. Our findings provide insights into the population immunogenetic substructure of the whole country and add to the knowledge of mixed ancestry Latin American population genetics, important for disease association studies, detection of demographic signatures on population variation and improved allocation of public health resources.

Clinical features, biochemistry and HLA-DRB1 status in children and adolescents with diabetes in Dhaka, Bangladesh

AIMS

Little information is published on diabetes in young people in Bangladesh. We aimed to investigate the demographic, clinical, and biochemical features, and HLA-DRB1 alleles in new cases of diabetes affecting Bangladeshi children and adolescents <22 years of age.

METHODS

The study was conducted at Bangladesh Institute of Research and Rehabilitation of Diabetes, Endocrine and Metabolic Disorders (BIRDEM) in Dhaka. One hundred subjects aged <22 years at diagnosis were enrolled. Demographic characteristics, clinical information, biochemical parameters (blood glucose, HbA1c, C-peptide, and autoantibodies against glutamic acid decarboxylase 65 (GADA) and islet antigen-2 (IA-2A) were measured. High-resolution DNA genotyping was performed for HLA-DRB1.

RESULTS

Eighty-four subjects were clinically diagnosed as type 1 diabetes (T1D), seven as type 2 diabetes (T2D), and nine as fibrocalculous pancreatic disease (FCPD). Of the 84 with T1D, 37 (44%) were males and 47 (56%) females, with median age at diagnosis 13 years (y) (range 1.6-21.7) and peak age at onset 12-15 years. 85% of subjects were assessed within one month of diagnosis and all within eleven months. For subjects diagnosed with T1D, mean C-peptide was 0.46 ± 0.22 nmol/L (1.40 ± 0.59 ng/mL), with 9 (10.7%) IA-2A positive, 22 (26%) GADA positive, and 5 (6%) positive for both autoantibodies. Analysis of HLA-DRB1 genotypes revealed locus-level T1D association (p = 6.0E-05); DRB1*04:01 appeared predisposing (p < 3.0E-06), and DRB1*14:01 appeared protective (p = 1.7E-02).

CONCLUSIONS

Atypical forms of T1D appear to be more common in young people in Bangladesh than in European populations. This will be helpful in guiding more specific assessment at onset and potentially, expanding treatment options.

Next-generation sequencing reveals new information about HLA allele and haplotype diversity in a large European American population

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

High-resolution HLA allele and haplotype frequencies in majority and minority populations of Costa Rica and Nicaragua: Differential admixture proportions in neighboring countries

The HLA system shows the most extensive polymorphism in the human genome. Allelic and haplotypic frequencies of HLA genes vary dramatically across human populations. Due to a complex history of migration, populations in Latin America show a broad variety of admixture proportions, usually varying not only between countries, but also within countries. Knowledge of HLA allele and haplotype frequencies is essential for medical fields such as transplantation, but also serves as a means to assess genetic diversity and ancestry in human populations. Here, we have determined high-resolution HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies in a sample of 713 healthy subjects from three Mestizo populations, one population of African descent, and Amerindians of five different groups from Costa Rica and Nicaragua and compared their profiles to a large set of indigenous populations from Iberia, Sub-Saharan Africa, and the Americas. Our results show a great degree of allelic and haplotypic diversity within and across these populations, with most extended haplotypes being private. Mestizo populations show alleles and haplotypes of putative European, Amerindian, and Sub-Saharan African origin, albeit with differential proportions. Despite some degree of gene flow, Amerindians and Afro-descendants show great similarity to other Amerindian and West African populations, respectively. This is the first comprehensive study reporting high-resolution HLA diversity in Central America, and its results will shed light into the genetic history of this region while also supporting the development of medical programs for organ and stem cell transplantation.

Semen quality is affected by HLA class I alleles together with sexually transmitted diseases

BACKGROUND

The human leukocyte antigen (HLA) locus includes several genes with key roles in antigen presentation and immune response, some of them inclusively found to be associated with non-obstructive azoospermia. Still, HLA connections to other infertility phenotypes such as semen hyperviscosity (SHV), asthenozoospermia (AST), and oligozoospermia (OLI) have been often neglected.

OBJECTIVES

In this work, we aimed to evaluate the association of HLA class I and II genes with SHV, AST, and OLI phenotypes while exploring a possible role in an adaptive immune response to sexually transmitted diseases (STD).

MATERIALS AND METHODS

Whole-exome sequencing was performed in a Portuguese cohort of 71 infertility cases and 68 controls, followed by HLA typing using a specific software-HLA*PRG:LA tool. Molecular screenings of seven STD were carried out in a subset of 72 samples (30 cases and 42 controls).

RESULTS

Statistical tests uncovered three protective alleles: HLA-A*11:01, associated with all forms of male infertility (p = 0.0006); HLA-DQB1*03:02 with SHV and OLI (P_SHV   = 0.0303, P_OLI   = 0.0153); and HLA-A*29:02 with OLI (p = 0.0355), which was found to interfere in sperm number together with HPV (p = 0.0313). Five risk alleles were also identified: two linked with SHV (HLA-B*50:01, p = 0.0278; and HLA-C*06:02, p = 0.0461), another one with both SHV and OLI (HLA-DQA1*05:01, P_SHV   = 0.0444 and P_OLI =0.0265), and two with OLI (HLA-C*03:03, p = 0.0480; and HLA-DQB1*03:01, p = 0.0499). Here, HLA-C*03:03 carriers tend to be HPV infected.

CONCLUSIONS

The application of HLA*PRG:LA tool to the study of male infertility provided novel insights for an HLA correlation with semen quality, namely among SHV and OLI phenotypes. The discovery of an HLA-A*29:02/HPV crosstalk, together with former reports of HLA alleles conferring resistance-susceptibility to diverse human pathogens, raises the hypothesis of a mechanistic link between male infertility, HLA polymorphism, and host response to STD.

High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis

We investigated association between HLA class I and class II alleles and haplotypes, and KIR loci and their HLA class I ligands, with multiple sclerosis (MS) in 412 European American MS patients and 419 ethnically matched controls, using next-generation sequencing. The DRB1*15:01~DQB1*06:02 haplotype was highly predisposing (odds ratio (OR) = 3.98; 95% confidence interval (CI) = 3-5.31; p-value (p) = 2.22E-16), as was DRB1*03:01~DQB1*02:01 (OR = 1.63; CI = 1.19-2.24; p = 1.41E-03). Hardy-Weinberg (HW) analysis in MS patients revealed a significant DRB1*03:01~DQB1*02:01 homozyote excess (15 observed; 8.6 expected; p = 0.016). The OR for this genotype (5.27; CI = 1.47-28.52; p = 0.0036) suggests a recessive MS risk model. Controls displayed no HW deviations. The C*03:04~B*40:01 haplotype (OR = 0.27; CI = 0.14-0.51; p = 6.76E-06) was highly protective for MS, especially in haplotypes with A*02:01 (OR = 0.15; CI = 0.04-0.45; p = 6.51E-05). By itself, A*02:01 is moderately protective, (OR = 0.69; CI = 0.54-0.87; p = 1.46E-03), and haplotypes of A*02:01 with the HLA-B Thr80 Bw4 variant (Bw4T) more so (OR = 0.53; CI = 0.35-0.78; p = 7.55E-04). Protective associations with the Bw4 KIR ligand resulted from linkage disequilibrium (LD) with DRB1*15:01, but the Bw4T variant was protective (OR = 0.64; CI = 0.49-0.82; p = 3.37-04) independent of LD with DRB1*15:01. The Bw4I variant was not associated with MS. Overall, we find specific class I HLA polymorphisms to be protective for MS, independent of the strong predisposition conferred by DRB1*15:01.

Clinical features, biochemistry and HLA-DRB1 status in youth-onset type 1 diabetes in Pakistan

Published information on diabetes in Pakistani youth is limited. We aimed to investigate the demographic, clinical, and biochemical features, and HLA-DRB1 alleles in new cases of diabetes affecting children and adolescents <22 years of age. The study was conducted at Baqai Institute of Diabetology and Endocrinology in Karachi from June 2013-December 2015. One hundred subjects aged <22 years at diagnosis were enrolled. Demographic characteristics, clinical information, biochemical parameters (blood glucose, HbA1c, C-peptide, glutamic acid decarboxylase 65 (GAD65) and islet antigen 2 (IA-2) autoantibodies) were measured. DNA from 100 subjects and 200 controls was extracted and genotyped for HLA-DRB1 using high-resolution genotyping technology. Ninety-nine subjects were clinically diagnosed as type 1 diabetes (T1D) and one as type 2 diabetes (T2D). Of the 99 with T1D, 57 (57.6%) were males and 42 (42.4%) females, with mean age at diagnosis 11.0 ± 5.2 years (range 1.6-21.7 years) and peaks at six and fifteen years. Fifty-seven subjects were assessed within one month of diagnosis and all within eleven months. For the subjects diagnosed as T1D, mean C-peptide was 0.63 ± 0.51 nmol/L (1.91 ± 1.53 ng/mL), with 16 (16.2%) IA2 positive, 53 (53.5%) GAD-65 positive, and 10 (10.1%) positive for both autoantibodies. In T1D patients, the allele DRB1*03:01 demonstrated highly significant T1D association (p < 10-16), with no apparent risk conferred by DRB1*04:xx alleles. CONCLUSIONS: Heterogeneous forms of T1D appear more common in children and youth in Pakistan than in European populations. Individual understanding of such cases could enable improved management strategies and healthier outcomes.

Epidemiology of childhood-onset type 1 diabetes in Azerbaijan: Incidence, clinical features, biochemistry, and HLA-DRB1 status

AIMS

Determine the incidence and typology of diabetes in children in Azerbaijan.

METHODS

Clinical features, C-peptide, autoantibodies (glutamic acid decarboxylase 65 (GAD65) and islet antigen 2 (IA-2)), and HLA-DRB1 status were studied in 106 subjects <18 years of age who were recently diagnosed. 104 cases were consecutive. Incidence was determined for Baku and Absheron regions, where ascertainment is estimated to be essentially 100%.

RESULTS

104 of the 106 (98%) were diagnosed with type 1 diabetes, one with type 2 diabetes and one with atypical diabetes. Type 1 diabetes incidence in Baku City and Absheron was 7.05 per 100,000 population <15 years per year. Peak age of onset was 10 years. There was a slight male preponderance (male:female 1.17:1), and no temporal association with seasons. Almost all type 1 diabetes subjects presented with classic symptoms including a high incidence (58%) of diabetic ketoacidosis. 86% presented with low C-peptide values (<0.13 nmol/L, <0.40 ng/mL) and 74% were positive for at least one type 1 diabetes-related autoantibody.

CONCLUSIONS

Azerbaijan has a moderate type 1 diabetes incidence and clinical, biochemical and genetic features similar to that in European populations.

HLA Population Genetics in Solid Organ Transplantation

HLAs are fundamental to the adaptive immune response and play critical roles in the cellular and humoral response in solid organ transplantation. The genes encoding HLA proteins are the most polymorphic within the human genome, with thousands of different allelic variants known within the population. Application of the principles of population genetics to the HLA genes has resulted in the development of a numeric metric, the calculated panel-reactive antibody (CPRA) that predicts the likelihood of a positive crossmatch as a function of a transplant candidate's unacceptable HLA antigens. The CPRA is an indispensible measure of access to transplantation for sensitized candidates and is used as the official measure of sensitization for allocation of points in the US Kidney Allocation System and Eurotransplant. Here, we review HLA population genetics and detail the mathematical basis of the CPRA. An understanding of these principles by transplant clinicians will lay the foundation for continued innovation in the care of sensitized patients.

HLA-A, HLA-B, and HLA-DRB1 Allele and Haplotype Frequencies in Renal Transplant Candidates in a Population in Southern Brazil

BACKGROUND

Very few studies have examined the diversity of human leukocyte antigens (HLA) in the Brazilian renal transplant candidates.

METHODS

The frequencies of the HLA-A, HLA-B, and HLA-DRB1 alleles, haplotypes and phenotypes were studied in 522 patients with chronic renal failure, renal transplant candidates, registered at the Transplant Centers in north/northwestern Paraná State, southern Brazil. Patients were classified according to the ethnic group (319 whites [Caucasians], 134 mestizos [mixed race descendants of Europeans, Africans, and Amerindians; browns or "pardos"] and 69 blacks). The HLA typing was performed by the polymerase chain reaction sequence-specific oligonucleotide method (PCR-SSO), combined with Luminex technology.

RESULTS

In the analysis of the total samples, 20 HLA-A, 32 HLA-B, and 13 HLA-DRB1 allele groups were identified. The most frequent allele groups for each HLA locus were HLA-A*02 (25.4%), HLA-B*44 (10.9%), and HLA-DRB1*13 (13.9%). The most frequent haplotypes were HLA-A*01-B*08-DRB1*03 (2.3%), A*02-B*44-DRB1*07 (1.2%), and A*03-B*07-DRB1*11 (1.0%). Significant differences (P < 0.05) were observed in the HLA-A*68, B*08, and B*58 allele frequencies among ethnic groups.

CONCLUSIONS

This study provides the first data on the HLA-A, HLA-B, and HLA-DRB1 allele, phenotype and haplotype frequencies of renal transplant candidates in a population in southern Brazil.

TNFA promoter alleles--frequencies and linkage with classical HLA genes in a South German Caucasian population

The tumor necrosis factor alpha (TNFA) promoter region exhibits several polymorphisms, which have been hypothesized to influence gene expression, thereby associating positively or negatively with inflammatory conditions. Many studies have focused on single nucleotide polymorphisms (SNPs) taking not into account additive or inverse effects between different SNPs. We typed 1,021 healthy Caucasian volunteer stem cell donors for their TNFA promoter as well as their HLA-A,-C,-B,-DRB1 loci. Using statistical methods, we reconstructed TNFA promoter alleles and analyzed their frequency and linkage with HLA genes. We show that the number of TNFA promoter alleles frequent enough to be analyzed in clinical studies is limited and that a strong linkage with classical HLA genes is present, especially for the extended HLA-haplotype HLA-A*01:01/HLA-C*07:01/HLA-B*08:01/TNFA promoter allele 3/HLA*DRB1*03:01. Taking into account SNP frequency information, it is possible to quite accurately deduce TNFA promoter alleles by generic Sanger sequencing, obviating the need for elaborating allele-specific sequencing. This information may enable investigators to consider the complete TNFA regulatory region in a phase-separated manner in contrast to previous approaches examining only one or few isolated SNPs.

Analytical methods for immunogenetic population data

In this chapter, we describe analyses commonly applied to immunogenetic population data, along with software tools that are currently available to perform those analyses. Where possible, we focus on tools that have been developed specifically for the analysis of highly polymorphic immunogenetic data. These analytical methods serve both as a means to examine the appropriateness of a dataset for testing a specific hypothesis, as well as a means of testing hypotheses. Rather than treat this chapter as a protocol for analyzing any population dataset, each researcher and analyst should first consider their data, the possible analyses, and any available tools in light of the hypothesis being tested. The extent to which the data and analyses are appropriate to each other should be determined before any analyses are performed.

Human leukocyte antigen class I (A, B, C) and II (DRB1) diversity in the black and Caucasian South African population

A cross-section of black and Caucasian South Africans (N = 302) were genotyped at high resolution (class I HLA-A, -B, -C and class II HLA-DRB1). Five new class I alleles (A*30:01:02, A*30:02:02, A*68:27, B*42:06, and B*45:07) and one new confirmatory allele (A*29:11) were identified in the black population. Alleles and haplotypes showed expected differences between the black and Caucasian populations, with the black population, on average, showing a broader spectrum of allele representation (less single allele dominance). The most prevalent alleles at the four loci in the black population were A*30:01, B*58:02, C*06:02, and DRB1*13:01 and in the Caucasian population were A*02:01:01, B*07:02:01, C*07:01, and DRB1*03:01. HLA-B, and HLA-C loci showed the strongest overall linkage disequilibrium (LD) and HLA-B/HLA-C two locus haplotypes also showed the strongest LD (D'(ij)) in both population groups. Bw allotype representation was similar between the two populations; however C allotypes differed significantly (C1 higher representation in Caucasians; C2 higher representation in blacks). HLA-A Supertype family phenotypic frequencies did not differ between the two populations, but four (B08, B27, B58, and B62) HLA-B Supertype families differed significantly. However, vaccine coverage estimation came close to 100% in both population groups, with inclusion of only four Supertype families (A1, A2, B7, B58).

Three microsatellites from the T1DGC MHC data set show highly significant association with type 1 diabetes, independent of the HLA-DRB1, -DQA1 and -DQB1 genes

AIM

The aim of this study was to test the microsatellites in the Type 1 Diabetes Genetics Consortium major histocompatibility complex (MHC) data set for association with type 1 diabetes (T1D) independent of the HLA-DRB1, -DQA1 and -DQB1 genes.

METHODS

The data set was edited to contain only one affected child per family, and broad ethnic subgroups were defined. Genotypes for HLA-DRB1, -DQA1 and -DQB1 were replaced by a haplotype code spanning all three loci, with phase inferred based on common haplotypes. The final data set contained 8190 samples in 2301 families, 59 microsatellites and the DRB1-DQA1-DQB1 haplotype code. Statistical analyses consisted of conditional logistic regression and haplotype estimations and linkage disequilibrium calculations.

RESULTS

The data set was screened using a main effects test approach adjusted for DRB1-DQA1-DQB1, and significant results tested for validity. After these procedures, four markers remained significant at the Bonferroni-corrected threshold: D6S2773 (p = 0.00014), DG6S185 (p = 0.00015), DG6S398 (p = 0.00043) and D6S2998 (p = 0.00015). These results were supported by allelic tests conditioned on DRB1-DQA1-DQB1 haplotypes, except for DG6S185, which may contain artefacts.

CONCLUSIONS

We have identified three microsatellites that mark additional risk factors for T1D at highly significant levels in the MHC. Further analyses are needed to establish the relationship with other possible genetic determinants in this region.

Racial/ethnic variation in EBV-positive classical Hodgkin lymphoma in California populations

Epstein-Barr virus (EBV) is detected in the tumor cells of some but not all Hodgkin lymphoma (HL) patients, and evidence indicates that EBV-positive and -negative HL are distinct entities. Racial/ethnic variation in EBV-positive HL in international comparisons suggests etiologic roles for environmental and genetic factors, but these studies used clinical series and evaluated EBV presence by differing protocols. Therefore, we evaluated EBV presence in the tumors of a large (n = 1,032), racially and sociodemographically diverse series of California incident classical HL cases with uniform pathology re-review and EBV detection methods. Tumor EBV-positivity was associated with Hispanic and Asian/Pacific Islander (API) but not black race/ethnicity, irrespective of demographic and clinical factors. Complex race-specific associations were observed between EBV-positive HL and age, sex, histology, stage, neighborhood socioeconomic status (SES), and birth place. In Hispanics, EBV-positive HL was associated not only with young and older age, male sex, and mixed cellularity histology, but also with foreign birth and lower SES in females, suggesting immune function responses to correlates of early childhood experience and later environmental exposures, respectively, as well as of pregnancy. For APIs, a lack of association with birth place may reflect the higher SES of API than Hispanic immigrants. In blacks, EBV-positive HL was associated with later-stage disease, consistent with racial/ethnic variation in certain cytokine polymorphisms. The racial/ethnic variation in our findings suggests that EBV-positive HL results from an intricate interplay of early- and later-life environmental, hormonal, and genetic factors leading to depressed immune function and poorly controlled EBV infection.

Deviations from Hardy-Weinberg proportions for multiple alleles under viability selection

Departures of genotype frequencies from Hardy-Weinberg proportions (HWP) for a single autosomal locus due to viability selection in a random mating population have been studied only for the two-allele case. In this article, the analysis of deviations from HWP due to constant viability selection is extended to multiple alleles. The deviations for an autosomal locus with k alleles are measured by means of k fii fixation indices for homozygotes and k(k-1)/2 fij fixation indices for heterozygotes, and expressions are obtained for these indices (FIS statistics) under the multiallele viability model. Furthermore, expressions for fii and fij when the multiallele polymorphism is at stable equilibrium are also derived and it is demonstrated that the pattern of multiallele Hardy-Weinberg deviations at equilibrium is characterized by a global heterozygote excess and a deficiency of each of the homozygotes. This pattern may be useful for detecting whether a given multiallelic polymorphism is at stable equilibrium in the population due to viability selection. An analysis of Hardy-Weinberg deviations from published data for the three-allele polymorphism at the beta-globin locus in human populations from West Africa is presented for illustration.

Clinical significance of HLA DRB103-DRB104 in type 1 autoimmune hepatitis

BACKGROUND

HLA DRB1*03-DRB1*04 combines both susceptibility factors for type-1 autoimmune hepatitis.

AIMS

Determine whether this phenotype is a risk factor for autoimmune hepatitis in white North American patients, assess its associations with clinical features and treatment outcome, and determine whether alleles within this phenotype affect prognosis.

METHODS

One hundred and ninety-eight patients with type 1 autoimmune hepatitis and 102 normal adults were evaluated. HLA typing was performed by DNA-based techniques.

RESULTS

Twenty-eight patients had HLA DRB1*03-DRB1*04, and the frequency was higher than in normal subjects (14% vs 4%, OR 4.0%, 95% CI 1.4-11.8, P = 0.01). Patients with DRB1*03-DRB1*04 relapsed less frequently than patients with DRB1*03 (1.3 +/- 0.3 vs 2.1 +/- 0.2, P = 0.04), but they otherwise had outcomes similar to patients with other phenotypes. Patients with DRB1*03-DRB1*04 who had 3-4 alleles encoding lysine at position DRbeta71 within the class II molecule of the major histocompatibility complex developed cirrhosis more commonly (75% vs 9%, P = 0.05) and had a higher frequency of hepatic-related death or liver transplantation (40% vs 0%, P = 0.04) than patients with fewer alleles.

CONCLUSIONS

HLA DRB1*03-DRB1*04 is a risk factor for type-1 autoimmune hepatitis, and its impact on outcome relates to the diversity of DRB1*04 alleles that encode a critical motif.

High-resolution molecular characterization of the HLA class I and class II in the Tarahumara Amerindian population

We describe for the first time the high-resolution profiling of HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1 in a culturally and geographically distinct Mexican ethnic group, the Tarahumaras. The alleles most frequently found by reference strand-mediated conformational analysis in this population were for class I: HLA-A*240201, *020101/09, *0206, *310102, *680102; HLA-B*4002, *1501, *510201, *3501/02/03, *4005, *4801; HLA-Cw*0304, *0801, *0102, *040101; and for class II: HLA-DRB1*080201, *1402, *040701; HLA-DQB1*0402, *0301, *0302/07; HLA-DPB1*0402, *0401, *020102. In addition, a novel allele, HLA-A*0257, was found. Based on comparison of presently known HLA-DRB1 and -DQB1 allele frequencies in Amerindian groups and worldwide populations, the Tarahumaras are unexpectedly more related to the geographically and linguistically distant Aymara and Terena Amerindian groups than they are to neighbouring tribes.

Hardy-Weinberg testing of a single homozygous genotype

No proper statistical test is available for the evaluation of deviation of a single homozygous genotype from Hardy-Weinberg equilibrium (HWE) proportion. We propose a 1-d.f. chi2-test. The power of the proposed test is favorable compared to existing HWE testing procedures. The applications of this test are discussed.

HLA DRB1* and DQB1* allelic diversity among nadars: A primitive South Indian Dravidian caste group

South India is one of the oldest geophysical regions occupied mainly by Dravidian language-speaking people. In the present study, 84 unrelated healthy Nadar individuals from Tamil Nadu State were analyzed for human leukocyte antigen (HLA) DRB1* and DQB1* allele frequencies and their haplotype frequencies were analyzed by high resolution polymerase chain reaction (PCR)-typing of genomic DNA. Results revealed that Nadars possessed >0.1 allele frequencies of HLA DRB1*15011, 03, 04, 1404, 0701, 1001, and DQB1*0201, 03, 0501, and 0503. The highest frequencies of DRB1*0317 (0.048), DRB1*1317 (0.012), and DRB1*1404 (0.107) in the world were found among the Nadar caste. The most frequent two-locus haplotypes identified among the Nadars were DRB1*0701-DQB1*0201 and DRB1*1000-DQB1*0501, and DRB1*1404-DQB1*0503 and DRB1*15011-DQB1*0601. When compared with world populations, HLA DRB1*1322, DRB1*1434, and DQB1*0614 were identified exclusively among Nadars. The heterozygosity and Watterson's homozygosity tests indicated that balancing selection still operates among the South Indian caste groups. The hierarchical cluster analysis for DRB1 data from the Indian caste and population groups showed that the Nadars cluster with Kashmir Brahmins and North Indians. The present study shows that the Nadar caste has several unique alleles and haplotypes at high frequency that are rare or absent in other populations of India and the world.

Perspective: detecting adaptive molecular polymorphism: lessons from the MHC

In the 1960s, when population geneticists first began to collect data on the amount of genetic variation in natural populations, balancing selection was invoked as a possible explanation for how such high levels of molecular variation are maintained. However, the predictions of the neutral theory of molecular evolution have since become the standard by which cases of balancing selection may be inferred. Here we review the evidence for balancing selection acting on the major histocompatibility complex (MHC) of vertebrates, a genetic system that defies many of the predictions of neutrality. We apply many widely used tests of neutrality to MHC data as a benchmark for assessing the power of these tests. These tests can be categorized as detecting selection in the current generation, over the history of populations, or over the histories of species. We find that selection is not detectable in MHC datasets in every generation, population, or every evolutionary lineage. This suggests either that selection on the MHC is heterogeneous or that many of the current neutrality tests lack sufficient power to detect the selection consistently. Additionally, we identify a potential inference problem associated with several tests of neutrality. We demonstrate that the signals of selection may be generated in a relatively short period of microevolutionary time, yet these signals may take exceptionally long periods of time to be erased in the absence of selection. This is especially true for the neutrality test based on the ratio of nonsynonymous to synonymous substitutions. Inference of the nature of the selection events that create such signals should be approached with caution. However, a combination of tests on different time scales may overcome such problems.

HLA-DRB1*, -DQB1* in Piramalai Kallars and Yadhavas, two Dravidian-speaking castes of Tamil Nadu, South India

Two Dravidian-speaking castes of Tamil Nadu, Piramalai Kallars (PKs, n = 205) and Yadhavas (YDs, n = 239) and a random panel (84) were studied for HLA-DRB1* and -DQB1* polymorphisms by DNA-SSOP typing methods. XI and XII International Histocompatibility primers and non-radioactive-labelled oligo probes were employed to identify the alleles. Results revealed that PKs possessed >0.1 allele frequencies of HLA-DRB1*15011, 0301, -DQB1*0201, 0501 and 0601; YDs, HLA-DRB1*0301, 0401, 07 and -DQB1*0601; and the random panel, DRB1*15021, 0401, 07, -DQB1 0201, 0301, 0302 and 0501. The highest frequency of DRB1*1501 in the world (GF = 0.225) was found in PKs. The most frequent two-locus haplotype (>500/10,000) in all the study samples was DRB1*10-DQB1*0501, while 1501-0601 was frequent in PKs and YDs. Comparison of the HLA-DRB1* data with Eastern European and South-East Asian populations suggested migration as the prime cause of the observed diversity in DRB1* allele frequencies. Nonetheless, the heterozygocity test and Watterson's homozygosity test indicated that balancing selection still operates on HLA-DRB1* locus, in this endemic region of various infectious diseases. This and spatial autocorrelation analysis support the view that selection may be a cause of "generating" new variants and allelic diversity in different ancient settlements. The study suggested that South Indian, inbred, endogamous, sympatrically isolated castes or similar well-defined breeding isolates around the world, living under the same milieu-epidemiology, may be ideal models to test the immunogenetic basis of disease susceptibility.

HLA antigens, alleles and haplotypes among the Yup'ik Alaska natives: report of the ASHI Minority Workshops, Part II

As part of the American Society for Histocompatibility and Immunogenetics coordinated studies among minority populations, human leukocyte antigen (HLA) alleles were defined for 460 volunteer Yup'ik Eskimos from the Yukon Kuskokwim delta region of southwestern Alaska. The study group included 252 adults with no other first-degree relatives and 48 informative nuclear families. Full Yupik ancestry through both maternal and paternal grandparents was claimed by 81.1% of participants. HLA-A, -B, -Cw, -DRB1, and -DQB1 alleles were determined by SBT, SSOP, reverse SSOP, and/or RSCA according to the protocols of five participating laboratories. Polymorphism was limited with 3-6 alleles comprising > 80% of the alleles observed at each locus. Homozygosity was high, particularly at the HLA-A and -DQB1 loci, with 36.6% and 44% of individuals having a single allele defined at these respective loci. HLA-A, -B, and -DRB1 alleles were in Hardy-Weinberg equilibrium, whereas HLA-Cw and -DQB1 alleles gave significant deviation (p = 0.002; 0.005). Significant linkage disequilibrium (p < or = 0.00001) was observed in all pairwise evaluations. A new Cw*0806 allele was observed in high linkage disequilibrium with B*4801(Delta = 0.099; Delta(rel) = 1.0). Three extended haplotypes were found to have frequencies > 5%, the most prevalent being A*2402; B*4801; DRB1*0401; DQB1*0301 (0.0933). Comparison of available class I data indicate that the Yup'ik share several common alleles with other Native American populations, including: A*2402, *0206, *6801; B*1501, *2705, *3501, *4002, *4801, *5101; and Cw*0202, *0304, *0401. Comparisons of class II data also confirm a close relationship of the Yup'ik to two other Eskimo populations, Siberian and East Greenland Eskimos. DRB1*0401 and *1101, which occur in high frequency among these Eskimo populations, but not in other Native Americans, were also prevalent among the Yup'ik, with respective frequencies of 0.232 and 0.107.

Molecular variability in Amerindians: widespread but uneven information

A review was made in relation to the molecular variability present in North, Central, and South American Indian populations. It involved results from ancient DNA, mitochondrial DNA in extant populations, HLA and other autosomal markers, X and Y chromosome variation, as well as data from parasitic viruses which could show coevolutionary changes. The questions considered were their origin, ways in which the early colonization of the continent took place, types and levels of the variability which developed, peculiarities of the Amerindian evolutionary processes, and eventual genetic heterogeneity which evolved in different geographical areas. Although much information is already available, it is highly heterogeneous in relation to populations and types of genetic systems investigated. Unfortunately, the present trend of favoring essentially applied research suggest that the situation will not basically improve in the future.