The equivalence of two tests and models for HLA data with no observed double blanks

Gart and Nam (1984, Biometrics 40, 887-894) consider the analysis of an ABO-like model based on the Hardy-Weinberg law in the commonly occurring case of human leukocyte antigen (HLA) data where there are no double blanks, that is, no recessive homozygotes. They derive a score test, based on the truncated likelihood, of the hypothesis that the true recessive gene (or allele) frequency is zero. Yasuda (1968, Biometrics 24, 915-935) considers a similar codominant system wherein the true recessive gene frequency is assumed zero, but the Hardy-Weinberg law does not hold. In particular, he considers the possibility of a nonzero inbreeding coefficient. We show that the two models are equivalent; each likelihood can be shown to be a reparameterization of the other. Furthermore, the score test of the zero gene frequency in Gart and Nam is identical to the score test for a zero inbreeding coefficient given by Yasuda. The results are applied to an example wherein it appears that the Hardy-Weinberg model is appropriate. Thus, it is not possible in this population to identify homozygous individuals without error from phenotypic data alone.     

The ML estimation and testing of generalized ABO-like data with no observed double recessives

The general method of the discrepancy or heterogeneity chi-square is applied to ABO-like data in which there are no observed double blanks in either the disease or the control group. When the recessive gene frequency is assumed zero, this method leads to an approximate chi-square test identical to that suggested by Smouse and Williams (1982, Biometrics 38, 757-768). When this assumption is relaxed, there arise two cases which are determined by whether the maximum likelihood estimate of this frequency is zero or not. It is shown that the value of the simple score statistic of Gart and Nam (1984, Biometrics 40, 887-894) discriminates between the two cases. The various omnibus test statistics for comparing groups are shown to differ little in several practical examples. However, under the more general assumption the appropriate degrees of freedom is one more than the number previously suggested.     

A modular concept of HLA for comprehensive peptide binding prediction

A variety of algorithms have been successful in predicting human leukocyte antigen (HLA)-peptide binding for HLA variants for which plentiful experimental binding data exist. Although predicting binding for only the most common HLA variants may provide sufficient population coverage for vaccine design, successful prediction for as many HLA variants as possible is necessary to understand the immune response in transplantation and immunotherapy. However, the high cost of obtaining peptide binding data limits the acquisition of binding data. Therefore, a prediction algorithm, which applies the binding information from well-studied HLA variants to HLA variants, for which no peptide data exist, is necessary. To this end, a modular concept of class I HLA-peptide binding prediction was developed. Accurate predictions were made for several alleles without using experimental peptide binding data specific to those alleles. We include a comparison of module-based prediction and supertype-based prediction. The modular concept increased the number of predictable alleles from 15 to 75 of HLA-A and 12 to 36 of HLA-B proteins. Under the modular concept, binding data of certain HLA alleles can make prediction possible for numerous additional alleles. We report here a ranking of HLA alleles, which have been identified to be the most informative. Modular peptide binding prediction is freely available to researchers on the web at http://www.peptidecheck.org .     

Approximate interval estimation of the difference in binomial parameters: correction for skewness and extension to multiple tables

Recently, Beal (1987, Biometrics 43, 941-950) found Mee's modification of Anbar's approximate interval estimation for the difference in binomial parameters to be a good choice in small sample sizes. As this method can be derived from the score theory of Bartlett, it is easily corrected for skewness. Exact numerical evaluation shows that this correction is not as important for this case as for the ratio of binomial parameters (Gart and Nam, 1988, Biometrics 44, 323-338). The score theory is also used to extend this method to the stratified or multiple-table case. Thus, good approximate interval estimates for differences, ratios, and odds ratios of binomial parameters can all be derived from the same general theory.     

Testing the Hardy-Weinberg equilibrium in the HLA system

A new method of testing the Hardy-Weinberg equilibrium in the human leukocyte antigen (HLA) system is proposed and applied to real data. The derivation is based on the maximum likelihood method and closely related to standard regression theory. The test statistic has a closed representation of residual sum of squares by a projection mapping of data onto the estimated regression plane. Under the Hardy-Weinberg law the noniterative estimates for the gene frequencies are suggested by the use of the projection mapping. The test statistic and gene frequency estimates are shown to be asymptotically equivalent to the maximum likelihood method and to be more efficient than the other suggested test statistic when there are more than two identified alleles.     

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.     

Linkage analysis of human leukocyte antigen (HLA) markers in familial psoriasis: strong disequilibrium effects provide evidence for a major determinant in the HLA-B/-C region.

Although psoriasis is strongly associated with certain human leukocyte antigens (HLAs), evidence for linkage to HLA markers has been limited. The objectives of this study were (1) to provide more definitive evidence for linkage of psoriasis to HLA markers in multiplex families; (2) to compare the major HLA risk alleles in these families with those determined by previous case-control studies; and (3) to localize the gene more precisely. By applying the transmission/disequilibrium test (TDT) and parametric linkage analysis, we found evidence for linkage of psoriasis to HLA-C, -B, -DR, and -DQ, with HLA-B and -C yielding the most-significant results. Linkage was detectable by parametric methods only when marker-trait disequilibrium was considered. Case-control association tests and the TDT identified alleles belonging to the EH57.1 ancestral haplotype as the major risk alleles in our sample. Among individuals carrying recombinant ancestral haplotypes involving EH57. 1, the class I markers were retained selectively among affecteds four times more often than among unaffecteds; among the few affected individuals carrying only the class II alleles from the ancestral haplotype, all but one also carried Cw6. These data show that familial and "sporadic" psoriasis share the same risk alleles. They also illustrate that substantial parametric linkage information can be extracted by accounting for linkage disequilibrium. Finally, they strongly suggest that a major susceptibility gene resides near HLA-C.     

云南纳西族HLA—DRB1基因多态性研究及其族源分析

首次在国内采用本室改进的高分辨率的基于内含子的PCR－SBT分型方法,检测云南纳西族HLA－DRB1基因多态性。在60例纳西族个体中共检出37种HLA－DRB1等位基因,其显著特点是等位基因的类型检出较多,频率分布比较平均,除12021（17.50%）外其他的等位基因频率均低于8％,其他较常见的等位基因（＞5％）还有1404（7.50%）,1504（5.83%0,04051（5.83%0,08032（5.83%）,09012（5％）,03011（5％）。这几种中频等位基因共占可检出等位基因的35％,与12021一起共占52.49%,其中DRB1＊0305、0438、1123、1132、1310、0812为首次在我国人群中检出,并且在世界各地人群中也比较罕见。以纳西族和世界各地人群的HLA－DRB1频率进行了聚类分析。比较分析的结果显示纳西族明显属于中国南方族群,未显示出其族源来自北方的痕迹。根据遗传数据,并参照民族学、历史学研究,对其民族起源做了初步的分析。     

A statistical method for predicting classical HLA alleles from SNP data

Genetic variation at classical HLA alleles is a crucial determinant of transplant success and susceptibility to a large number of infectious and autoimmune diseases. However, large-scale studies involving classical type I and type II HLA alleles might be limited by the cost of allele-typing technologies. Although recent studies have shown that some common HLA alleles can be tagged with small numbers of markers, SNP-based tagging does not offer a complete solution to predicting HLA alleles. We have developed a new statistical methodology to use SNP variation within the region to predict alleles at key class I (HLA-A, HLA-B, and HLA-C) and class II (HLA-DRB1, HLA-DQA1, and HLA-DQB1) loci. Our results indicate that a single panel of approximately 100 SNPs typed across the region is sufficient for predicting both rare and common HLA alleles with up to 95% accuracy in both African and non-African populations. Furthermore, we show that HLA alleles can be successfully predicted by using previously genotyped SNPs that are within the MHC and that had not been chosen for their ability to predict HLA alleles, such as those included on genome-wide products. These results indicate that our methodology, combined with an extended database of reference haplotypes, will facilitate large-scale experiments, including disease-association studies and vaccine trials, in which detailed information about HLA type is valuable.     

Tracking human migrations by the analysis of the distribution of HLA alleles, lineages and haplotypes in closed and open populations

The human leucocyte antigen (HLA) system shows extensive variation in the number and function of loci and the number of alleles present at any one locus. Allele distribution has been analysed in many populations through the course of several decades, and the implementation of molecular typing has significantly increased the level of diversity revealing that many serotypes have multiple functional variants. While the degree of diversity in many populations is equivalent and may result from functional polymorphism(s) in peptide presentation, homogeneous and heterogeneous populations present contrasting numbers of alleles and lineages at the loci with high-density expression products. In spite of these differences, the homozygosity levels are comparable in almost all of them. The balanced distribution of HLA alleles is consistent with overdominant selection. The genetic distances between outbred populations correlate with their geographical locations; the formal genetic distance measurements are larger than expected between inbred populations in the same region. The latter present many unique alleles grouped in a few lineages consistent with limited founder polymorphism in which any novel allele may have been positively selected to enlarge the communal peptide-binding repertoire of a given population. On the other hand, it has been observed that some alleles are found in multiple populations with distinctive haplotypic associations suggesting that convergent evolution events may have taken place as well. It appears that the HLA system has been under strong selection, probably owing to its fundamental role in varying immune responses. Therefore, allelic diversity in HLA should be analysed in conjunction with other genetic markers to accurately track the migrations of modern humans.     

A HLA-DRB supertype chart with potential overlapping peptide binding function

HLA-DRB alleles are class II alleles that are associated with CD4+ T-cell immune response. DRB alleles are polymorphic and currently there are about 622 named in the IMGT/HLA sequence database. Each allele binds short peptides with high sensitivity and specificity. However, it has been suggested that majority of HLA alleles can be covered within few HLA supertypes, where different members of a supertype bind similar peptides showing distinct repertoires. Definition of DRB supertypes using binding data is limited to few (about 29) known alleles (< 5% of all known DRB alleles). Hence, we describe a strategy using structurally defined virtual pockets to group all known DRB alleles with regard to their overlapping peptide binding specificity.     

中国壮族和裕固族群体HLA Ⅰ类区域Alu插入多态性及其与HLA-A位点的相关性

近年来研究发现:位于HLAⅠ类基因区域的Alu插入是研究不同群体HLAⅠ类基因区域祖先单倍型和HLAⅠ类基因多样性产生、进化和重组的理想工具。文章对中国壮族和裕固族群体HLAⅠ类基因区域5个Alu插入多态性(AluMICB、AluTF、AluHJ、AluHG和AluHF)进行研究,结合HLA基因分型数据,分析壮族、裕固族、哈尼族、布朗族和傣族5个民族群体中Alu插入与HLA-A等位基因的关系。研究结果显示:(1)壮族和裕固族人群中5个Alu插入频率范围分别为1.5%~35.8%和9.2~34.8%,AluMICB、AluTF和AluHF插入频率在这两个群体中有统计学差异(P<0.05);(2)在5个研究的群体中,AluHG插入与HLA-A*02的不同亚型关联;AluHJ插入与HLA-A*2402在5个群体中都关联,但AluHJ与HLA-A*1101和HLA-A*2407只在布朗族中关联。表明不同群体HLAⅠ类基因区域内Alu插入具有各自的特征,且Alu插入与不同的HLA-A等位基因相关联。这种Alu插入及其与HLA-A的关联特征可作为研究群体中HLAⅠ类基因和单倍型系谱变化的重要遗传标记。     

Repair of UV-induced pyrimidine dimers in the individual genes Gart, Notch and white from Drosophila melanogaster cell lines.

The excision repair of UV-induced pyrimidine dimers was investigated in three genes: Gart, Notch and white in a permanent Drosophila cell line Kc, derived from wild type Drosophila melanogaster embryonic cells. In this cell line Gart and Notch are actively transcribed, whereas white is not expressed. In all three genes UV-induced pyrimidine dimers were removed with the same rate and to the same extent: 60% removal within 16 hours, up to 80-100% in 24 hours after irradiation with 10 or 15 J/m2 UV. These kinetics are similar to the time course of dimer removal measured in the genome overall. No difference in repair of the inactive white locus compared to the active Gart and Notch genes was found. Similar results were obtained using a different wild type cell line, SL2, although repair appeared to be somewhat slower in this cell line. The results are discussed with respect to the data found for gene specific repair in other eukaryotic systems.     

Transporter associated with antigen processing preselection of peptides binding to the MHC: a bioinformatic evaluation

TAP is responsible for the transit of peptides from the cytosol to the lumen of the endoplasmic reticulum. In an immunological context, this event is followed by the binding of peptides to MHC molecules before export to the cell surface and recognition by T cells. Because TAP transport precedes MHC binding, TAP preferences may make a significant contribution to epitope selection. To assess the impact of this preselection, we have developed a scoring function for TAP affinity prediction using the additive method, have used it to analyze and extend the TAP binding motif, and have evaluated how well this model acts as a preselection step in predicting MHC binding peptides. To distinguish between MHC alleles that are exclusively dependent on TAP and those exhibiting only a partial dependence on TAP, two sets of MHC binding peptides were examined: HLA-A*0201 was selected as a representative of partially TAP-dependent HLA alleles, and HLA-A*0301 represented fully TAP-dependent HLA alleles. TAP preselection has a greater impact on TAP-dependent alleles than on TAP-independent alleles. The reduction in the number of nonbinders varied from 10% (TAP-independent) to 33% (TAP-dependent), suggesting that TAP preselection is an important component in the successful in silico prediction of T cell epitopes.     

CD8 T-cell recognition of multiple epitopes within specific Gag regions is associated with maintenance of a low steady-state viremia in human immunodeficiency virus type 1-seropositive patients

The importance of HLA class I-restricted CD8 T-cell responses in the control of human immunodeficiency virus (HIV) infection is generally accepted. While several studies have shown an association of certain HLA class I alleles with slower disease progression, it is not fully established whether this effect is mediated by HIV-specific CD8 T-cell responses restricted by these alleles. In order to study the influence of the HLA class I alleles on the HIV-specific CD8 T-cell response and on viral control, we have assessed HIV-specific epitope recognition, plasma viral load, and expression of HLA class I alleles in a cohort of HIV-seropositive bar workers. Possession of the HLA class I alleles B5801, B8101, and B0702 was associated with a low median viral load and simultaneously with a broader median recognition of Gag epitopes compared to all other HLA alleles (twofold increase) (P = 0.0035). We further found an inverse linear relationship between the number of Gag epitopes recognized and the plasma viral load (R = -0.36; P = 0.0016). Particularly, recognition of multiple epitopes within two regions of Gag (amino acids [aa] 1 to 75 and aa 248 to 500) was associated with the maintenance of a low steady-state viremia, even years after acute infection.     

The spectrum of HLA-DQ and HLA-DR alleles, 2006: a listing correlating sequence and structure with function

The list of alleles in the HLA-DRB, HLA-DQA, and HLA-DQB gene loci has grown enormously since the last listing in this journal 8 years ago. Crystal structure determination of several human and mouse HLA class II alleles, representative of two gene loci in each species, enables a direct comparison of ortholog and paralog loci. A new numbering system is suggested, extending earlier suggestions by [Fremont et al. in Immunity 8:305-317, (1998)], which will bring in line all the structural features of various gene loci, regardless of animal species. This system allows for structural equivalence of residues from different gene loci. The listing also highlights all amino acid residues participating in the various functions of these molecules, from antigenic peptide binding to homodimer formation, CD4 binding, membrane anchoring, and cytoplasmic signal transduction, indicative of the variety of functions of these molecules. It is remarkable that despite the enormous number of unique alleles listed thus far (DQA = 22, DQB = 54, DRA = 2, and DRB = 409), there is invariance at many specific positions in man, but slightly less so in mouse or rat, despite their much lower number of alleles at each gene locus in the latter two species. Certain key polymorphisms (from substitutions to an eight-residue insertion in the cytoplasmic tail of certain DQB alleles) that have thus far gone unnoticed are highly suggestive of differences or diversities in function and thus call for further investigation into the properties of these specific alleles. This listing is amenable to supplementation by future additions of new alleles and the highlighting of new functions to be discovered, providing thus a unifying platform of reference in all animal species for the MHC class II allelic counterparts, aiding research in the field and furthering our understanding of the functions of these molecules.     

Immunogenetic factors in beryllium sensitization and chronic beryllium disease

Exposure to beryllium in the workplace can cause beryllium sensitization and chronic beryllium disease. Sensitization to beryllium can be detected in the laboratory using the beryllium lymphocyte proliferation test. It was shown that anti-HLA antibodies could block the beryllium-specific response in the beryllium lymphocyte proliferation test, thereby implicating HLA genes in chronic beryllium disease. A supratypic genetic marker, HLA-DPB1*E69, was found to be strongly associated with immunologic sensitization to beryllium and chronic beryllium disease in beryllium workers. However, there are 40 HLA-DPB1 gene variants that have E69 but that also have other DNA sequence variations. The purpose of the study was to evaluate the evidence for potential differential susceptibility that may be associated with the physical characteristics of HLA protein molecules for which different HLA-DPB1*E69 variants code; that is, do some HLA-DPB1*E69 variants convey higher risk of beryllium sensitization and chronic beryllium disease than others. To do this, two approaches were pursued: first, detailed analysis of the findings from the published literature was performed, and second, computational chemistry was used to seek clues concerning the physical properties of the HLA protein molecules for which these alleles code. Among the 40 HLA-DPB1 gene variants that code for E69, molecular epidemiological studies have suggested a risk hierarchy, where some variants appear to convey low to moderate risk of chronic beryllium disease (e.g., HLA-DPB1*0201, approximately 3-fold increased risk), some convey an intermediate risk (e.g., HLA-DPB1*1901, approximately 5-fold) and others convey high risk (e.g., HLA-DPB1*1701, >10-fold). Molecular modeling has been used to further investigate a potential mechanistic basis for these observations. We found a strong correlation between the hierarchical order of risk of chronic beryllium disease associated with specific alleles and the predicted surface electrostatic potential and charge of the corresponding isotypes. Therefore, when alleles were grouped by the relative negative charge on the molecules for which they code, the data suggest that those alleles associated with the most negatively charged proteins carry the greatest risk of beryllium sensitization and disease.     

Allele excess at neutrally evolving microsatellites and the implications for tests of neutrality

Skews in the observed allele-frequency spectrum are frequently viewed as an indication of non-neutral evolution. Recent surveys of microsatellite variability have used an excess of alleles as a statistical approach to infer positive selection. Using neutral coalescent simulations we demonstrate that the mean numbers of alleles expected under the stepwise-mutation model and infinite-allele model deviate from the observed numbers of alleles. The magnitude of this difference is dependent on the sample size, mutation rates (theta-values) and observed gene diversities. Moreover, we show that the number of observed alleles differs among loci with the same observed gene diversity but different mutation rates (theta-values). We propose that a reliable test statistic based on allele excess must determine the confidence interval by computer simulations conditional on the observed gene diversity and theta-values. As the latter are notoriously difficult to obtain for experimental data, we suggest that other statistics, such as lnRV, may be better suited to the identification of microsatellite loci subject to selection.     

Leveraging information across HLA alleles/supertypes improves epitope prediction.

We present a model for predicting HLA class I restricted CTL epitopes. In contrast to almost all other work in this area, we train a single model on epitopes from all HLA alleles and supertypes, yet retain the ability to make epitope predictions for specific HLA alleles. We are therefore able to leverage data across all HLA alleles and/or their supertypes, automatically learning what information should be shared and also how to combine allele-specific, supertype-specific, and global information in a principled way. We show that this leveraging can improve prediction of epitopes having HLA alleles with known supertypes, and dramatically increases our ability to predict epitopes having alleles which do not fall into any of the known supertypes. Our model, which is based on logistic regression, is simple to implement and understand, is solved by finding a single global maximum, and is more accurate (to our knowledge) than any other model.