Racial differences in genome-wide methylation profiling and gene expression in breast tissues from healthy women

Breast cancer is more common in European Americans (EAs) than in African Americans (AAs) but mortality from breast cancer is higher among AAs. While there are racial differences in DNA methylation and gene expression in breast tumors, little is known whether such racial differences exist in breast tissues of healthy women. Genome-wide DNA methylation and gene expression profiling was performed in histologically normal breast tissues of healthy women. Linear regression models were used to identify differentially-methylated CpG sites (CpGs) between EAs (n = 61) and AAs (n = 22). Correlations for methylation and expression were assessed. Biological functions of the differentially-methylated genes were assigned using the Ingenuity Pathway Analysis. Among 485 differentially-methylated CpGs by race, 203 were hypermethylated in EAs, and 282 were hypermethylated in AAs. Promoter-related differentially-methylated CpGs were more frequently hypermethylated in EAs (52%) than AAs (27%) while gene body and intergenic CpGs were more frequently hypermethylated in AAs. The differentially-methylated CpGs were enriched for cancer-associated genes with roles in cell death and survival, cellular development, and cell-to-cell signaling. In a separate analysis for correlation in EAs and AAs, different patterns of correlation were found between EAs and AAs. The correlated genes showed different biological networks between EAs and AAs; networks were connected by Ubiquitin C. To our knowledge, this is the first comprehensive genome-wide study to identify differences in methylation and gene expression between EAs and AAs in breast tissues from healthy women. These findings may provide further insights regarding the contribution of epigenetic differences to racial disparities in breast cancer.     

Evidence for human meiotic recombination interference obtained through construction of a short tandem repeat-polymorphism linkage map of chromosome 19

An improved linkage map for human chromosome 19 containing 35 short tandem repeat polymorphisms (STRPs) and one VNTR (D19S20) was constructed. The map included 12 new (GATA)_n tetranucleotide STRPs. Although total lengths of the male (114 cM) and female (128 cM) maps were similar, at both ends of the chromosome male recombination exceeded female recombination, while in the interior portion of the map female recombination was in excess. Cosmid clones containing the STRP sequences were identified and were positioned along the chromosome by fluorescent in situ hybridization. Four rounds of careful checking and removal of genotyping errors allowed biologically relevant conclusions to be made concerning the numbers and distributions of recombination events on chromosome 19. The average numbers of recombinations per chromosome matched closely the lengths of the genetic maps computed by using the program CRIMAP. Significant numbers of chromosomes with zero, one, two, or three recombinations were detected as products of both female and male meioses. On the basis of the total number of observed pairs of recombination events in which only a single informative marker was situated between the two recombinations, a maximal estimate for the rate of meiotic STRP “gene” conversion without recombination was calculated as 3 × 10⁻⁴/meiosis. For distances up to 30 cM between recombinations, many fewer chromosomes which had undergone exactly two recombinations were observed than were expected on the basis of the assumption of independent recombination locations. This strong new evidence for human meiotic interference will help to improve the accuracy of interpretation of clinical DNA test results involving polymorphisms flanking a genetic abnormality.     

Linkage relationships and gene order around the locus for X-linked retinoschisis.

X-linked recessive retinoschisis (RS) is a hereditary disorder with variable clinical features. The main symptoms are poor sight; radial, cystic macula degeneration; and peripheral superficial retinal detachment. The disease is quite common in Finland, where at least 300 hemizygous males have been diagnosed. We used nine polymorphic DNA markers to study the localization of RS on the short arm of the X chromosome in 31 families comprising 88 affected persons. Two-point linkage results confirmed close linkage of the RS gene to the marker loci DXS43, DXS16, DXS207, and DXS41 and also revealed close linkage to the marker loci DXS197 and DXS9. Only one recombination was observed between DXS43 and RS in 59 informative meioses, giving a maximum lod score of 13.87 at the recombination fraction .02. No recombinations were observed between the RS locus and DXS9 and DXS197 (lods between 3 and 4), but at neither locus was the number of informative meioses sufficient to provide reliable estimates of recombination fractions. The most likely gene order on the basis of multilocus analysis was Xpter-DXS85-(DXS207,DXS43)-RS-DXS41-DXS 164-Xcen. Because multilocus linkage analysis indicated that the most probable location of RS is proximal to DXS207 and DXS43 and distal to DXS41, these three flanking markers are the closest and most informative markers currently available for carrier detection.     

Genetic association analyses of PDYN polymorphisms with heroin and cocaine addiction

Genetic factors are believed to account for 30-50% of the risk for cocaine and heroin addiction. Dynorphin peptides, derived from the prodynorphin (PDYN) precursor, bind to opioid receptors, preferentially the kappa-opioid receptor, and may mediate the aversive effects of drugs of abuse. Dynorphin peptides produce place aversion in animals and produce dysphoria in humans. Cocaine and heroin have both been shown to increase expression of PDYN in brain regions relevant for drug reward and use. Polymorphisms in PDYN are therefore hypothesized to increase risk for addiction to drugs of abuse. In this study, 3 polymorphisms in PDYN (rs1022563, rs910080 and rs1997794) were genotyped in opioid-addicted [248 African Americans (AAs) and 1040 European Americans (EAs)], cocaine-addicted (1248 AAs and 336 EAs) and control individuals (674 AAs and 656 EAs). Sex-specific analyses were also performed as a previous study identified PDYN polymorphisms to be more significantly associated with female opioid addicts. We found rs1022563 to be significantly associated with opioid addiction in EAs [P = 0.03, odds ratio (OR) = 1.31; false discovery rate (FDR) corrected q-value]; however, when we performed female-specific association analyses, the OR increased from 1.31 to 1.51. Increased ORs were observed for rs910080 and rs199774 in female opioid addicts also in EAs. No statistically significant associations were observed with cocaine or opioid addiction in AAs. These data show that polymorphisms in PDYN are associated with opioid addiction in EAs and provide further evidence that these risk variants may be more relevant in females.     

Is there a high rate of mitotic recombination between the loci encoding immunoglobulin VH and CH regions in gonial cells?

Rabbit sera from approximately 6000 offspring of matings informative for recombination at the immunoglobulin (Ig) heavy chain locus were tested with allotype antisera by double diffusion in gel. Seven recombinants were found, R1K-R7K, and in every case the recombinational event had taken place in the male parent. Two single males each fathered two recombinant offspring: R2K, R4K and R5K, R6K. Four further recombinants, reported from other laboratories as well as two new, as yet not fully documented recombinants in our laboratory, also occurred in the male parent. The recombinational events either seperated the a locus which encodes the V_H region of Ig, from the loci encoding its C_H region or separated a-msm from the remaining linked loci d-e-f-g. The recombination frequency is 7/6142, approximately 0.1%. Our finding suggests that the recombinations took place in mitotic divisions during spermatogenesis. This work is the first evidence of gonial crossing-over in a mammal.     

Adiponectin multimers and metabolic syndrome traits: relative adiponectin resistance in African Americans

African Americans (AAs) tend to have lower total adiponectin levels compared to European Americans (EA); however, it is not known whether race affects adiponectin multimer distribution and their relationships to metabolic traits. We measured total adiponectin, high molecular weight (HMW), low molecular weight (LMW) (i.e., hexamer), and trimer adiponectin in 132 normoglycemic premenopausal women (75 AAs, 57 EAs), together with measures of total and abdominal fat, plasma lipids, insulin sensitivity (S(i)), and genetic admixture estimates. We found that lower total adiponectin in AAs was explained by reduced LMW, and trimer forms because levels of HMW did not differ between races. In EAs, HMW was highly correlated with multiple metabolic syndrome traits. In contrast, the LMW and trimer forms were most highly correlated with metabolic traits in AAs, including abdominal adiposity, lipids, and S(i). At similar levels of visceral adiposity, AAs exhibited significantly lower LMW adiponectin than EAs. Similarly, at comparable levels of HMW and LMW adiponectin, AAs were more insulin resistant than their EA counterparts. In conclusion, (i) serum adiponectin is lower in AAs predominantly as a result of reduced concentrations of LMW and trimers multimeric forms; (ii) LMW and trimer, not HMW, are most broadly correlated with metabolic traits in AAs. Thus, HMW adiponectin may exert less bioactivity in explaining the metabolic syndrome trait cluster in populations of predominant African genetic background.     

Defining the Proximal Border of the Huntington Disease Candidate Region by Multipoint Recombination Analyses

The candidate region for the Huntington disease (HD) gene has been narrowed down to a 2.2-Mb region between D4S10 and D4S98 on the short arm of chromosome 4. To map the HD gene within this candidate region 65 Dutch HD families were studied. In total 338 informative meioses were analyzed and 11 multiple informative crossovers were detected. Assuming a minimum number of recombinations and no double recombinations, our multiple informative crossovers are consistent with one specific genetic order for 12 loci: D4S10-(D4S81, D4S126)-D4S125-(D4S127, D4S95)-D4S43-(D4S115, D4S96, D4S111, D4S90, D4S141). This is in agreement with the known data derived from similar and other methods. The loci between brackets could not be mapped relative to each other. In our family material, two informative three-point marker recombination events were detected in the proximal HD candidate region, which are also informative for HD. Both recombination events map the HD gene distal to D4S81 and most likely distal to D4S125, narrowing down the HD candidate region to a 1.7-Mb region between D4S125 and D4S98.     

Multiple cholinergic nicotinic receptor genes affect nicotine dependence risk in African and European Americans

Several independent studies show that the chromosome 15q25.1 region, which contains the CHRNA5–CHRNA3–CHRNB4 gene cluster, harbors variants strongly associated with nicotine dependence, other smoking behaviors, lung cancer and chronic obstructive pulmonary disease. We investigated whether variants in other cholinergic nicotinic receptor subunit (CHRN) genes affect the risk of nicotine dependence in a new sample of African Americans (AAs) (N = 710). We also analyzed this AA sample together with a European American (EA) sample (N = 2062, 1608 of which have been previously studied), allowing for differing effects in the two populations. Cases are current nicotine‐dependent smokers and controls are non‐dependent smokers. Variants in or near CHRND–CHRNG, CHRNA7 and CHRNA10 show modest association with nicotine dependence risk in the AA sample. In addition, CHRNA4, CHRNB3–CHRNA6 and CHRNB1 show association in at least one population. CHRNG and CHRNA4 harbor single nucleotide polymorphisms (SNPs) that have opposite directions of effect in the two populations. In each of the population samples, these loci substantially increase the trait variation explained, although no loci meet Bonferroni‐corrected significance in the AA sample alone. The trait variation explained by three key associated SNPs in CHRNA5–CHRNA3–CHRNB4 is 1.9% in EAs and also 1.9% in AAs; this increases to 4.5% in EAs and 7.3% in AAs when we add six variants representing associations at other CHRN genes. Multiple nicotinic receptor subunit genes outside chromosome 15q25 are likely to be important in the biological processes and development of nicotine dependence, and some of these risks may be shared across diverse populations.     

Estimating the contribution of mutation, recombination and gene conversion in the generation of haplotypic diversity

Recombination occurs through both homologous crossing over and homologous gene conversion during meiosis. The contribution of recombination relative to mutation is expected to be dramatically reduced in inbreeding organisms. We report coalescent-based estimates of the recombination parameter (rho) relative to estimates of the mutation parameter (theta) for 18 genes from the highly self-fertilizing grass, wild barley, Hordeum vulgare ssp. spontaneum. Estimates of rho/theta are much greater than expected, with a mean rho/theta approximately 1.5, similar to estimates from outcrossing species. We also estimate rho with and without the contribution of gene conversion. Genotyping errors can mimic the effect of gene conversion, upwardly biasing estimates of the role of conversion. Thus we report a novel method for identifying genotyping errors in nucleotide sequence data sets. We show that there is evidence for gene conversion in many large nucleotide sequence data sets including our data that have been purged of all detectable sequencing errors and in data sets from Drosophila melanogaster, D. simulans, and Zea mays. In total, 13 of 27 loci show evidence of gene conversion. For these loci, gene conversion is estimated to contribute an average of twice as much as crossing over to total recombination.     

Haplotypes histories as pathways of recombinations

MOTIVATION: The diversity of a haplotype, represented as a string of polymorphic sites along a DNA sequence, increases exponentially with the number of sites if recombinations are taking place. Reconstructing the history of recombinations compared with that of the polymorphic sites is thus extremely difficult. However, in the human genome, because of the relatively simple pattern of haplotype diversity dominated by a few ancestral haplotypes, the complexity of the recombinational network can be reduced, thus making its reconstruction feasible. We focus on the problem of inferring the recombination pathways starting with putative ancestral haplotypes and leading to new rare recombinant haplotypes. RESULTS: We describe classes of recombinations that represent the whole set of minimal recombination pathways leading to a new haplotype. We present an O(n(2)) algorithm that outputs such representative recombination pathways. We apply it to haplotypes of the 8 kb dystrophin gene segment dys44. AVAILABILITY: A software implementing the algorithm and some other extentions has been developed on a Java platform (JDK 1.3.1). It is freely available at http://www.iro.umontreal.ca/~mabrouk/     

Unraveling the consecutive recombination events in the human IGK locus

In addition to the classical Vkappa-Jkappa, Vkappa-kappa deleting element (Kde), and intron-Kde gene rearrangements, atypical recombinations involving Jkappa recombination signal sequence (RSS) or intronRSS elements can occur in the Igkappa (IGK) locus, as observed in human B cell malignancies. In-depth analysis revealed that atypical JkappaRSS-intronRSS, Vkappa-intronRSS, and JkappaRSS-Kde recombinations not only occur in B cell malignancies, but rather reflect physiological gene rearrangements present in normal human B cells as well. Excision circle analysis and recombination substrate assays can discriminate between single-step vs multistep rearrangements. Using this combined approach, we unraveled that the atypical Vkappa-intronRSS and JkappaRSS-Kde pseudohybrid joints most probably result from ongoing recombination following an initial aberrant JkappaRSS-intronRSS signal joint formation. Based on our observations in normal and malignant human B cells, a model is presented to describe the sequential (classical and atypical) recombination events in the human IGK locus and their estimated relative frequencies (0.2-1.0 vs < 0.03). The initial JkappaRSS-intronRSS signal joint formation (except for Jkappa1RSS-intronRSS) might be a side event of an active V(D)J recombination mechanism, but the subsequent formation of Vkappa-intronRSS and JkappaRSS-Kde pseudohybrid joints can represent an alternative pathway for IGK allele inactivation and allelic exclusion, in addition to classical Ckappa deletions. Although usage of this alternative pathway is limited, it seems essential for inactivation of those IGK alleles that have undergone initial aberrant recombinations, which might otherwise hamper selection of functional Ig L chain proteins.     

Common vitamin D pathway gene variants reveal contrasting effects on serum vitamin D levels in African Americans and European Americans

Vitamin D deficiency is more common among African Americans (AAs) than among European Americans (EAs), and epidemiologic evidence links vitamin D status to many health outcomes. Two genome-wide association studies (GWAS) in European populations identified vitamin D pathway gene single-nucleotide polymorphisms (SNPs) associated with serum vitamin D [25(OH)D] levels, but a few of these SNPs have been replicated in AAs. Here, we investigated the associations of 39 SNPs in vitamin D pathway genes, including 19 GWAS-identified SNPs, with serum 25(OH)D concentrations in 652 AAs and 405 EAs. Linear and logistic regression analyses were performed adjusting for relevant environmental and biological factors. The pattern of SNP associations was distinct between AAs and EAs. In AAs, six GWAS-identified SNPs in GC, CYP2R1, and DHCR7/NADSYN1 were replicated, while nine GWAS SNPs in GC and CYP2R1 were replicated in EAs. A CYP2R1 SNP, rs12794714, exhibited the strongest signal of association in AAs. In EAs, however, a different CYP2R1 SNP, rs1993116, was the most strongly associated. Our models, which take into account genetic and environmental variables, accounted for 20 and 28 % of the variance in serum vitamin D levels in AAs and EAs, respectively.     

Identification of POMC Exonic Variants Associated with Substance Dependence and Body Mass Index

Background

Risk of substance dependence (SD) and obesity has been linked to the function of melanocortin peptides encoded by the proopiomelanocortin gene (POMC).

Methods and Results

POMC exons were Sanger sequenced in 280 African Americans (AAs) and 308 European Americans (EAs). Among them, 311 (167 AAs and 114 EAs) were affected with substance (alcohol, cocaine, opioid and/or marijuana) dependence and 277 (113 AAs and164 EAs) were screened controls. We identified 23 variants, including two common polymorphisms (rs10654394 and rs1042571) and 21 rare variants; 12 of which were novel. We used logistic regression to analyze the association between the two common variants and SD or body mass index (BMI), with sex, age, and ancestry proportion as covariates. The common variant rs1042571 in the 3′UTR was significantly associated with BMI in EAs (Overweight: P _adj = 0.005; Obese: P _adj = 0.018; Overweight+Obese: P _adj = 0.002) but not in AAs. The common variant, rs10654394, was not associated with BMI and neither common variant was associated with SD in either population. To evaluate the association between the rare variants and SD or BMI, we collapsed rare variants and tested their prevalence using Fisher’s exact test. In AAs, rare variants were nominally associated with SD overall and with specific SD traits (SD: P _FET,1df = 0.026; alcohol dependence: P _FET,1df = 0.027; cocaine dependence: P _FET,1df = 0.007; marijuana dependence: P _FET,1df = 0.050) (the P-value from cocaine dependence analysis survived Bonferroni correction). There was no such effect in EAs. Although the frequency of the rare variants did not differ significantly between the normal-weight group and the overweight or obese group in either population, certain rare exonic variants occurred only in overweight or obese subjects without SD.

Conclusion

These findings suggest that POMC exonic variants may influence risk for both SD and elevated BMI, in a population-specific manner. However, common and rare variants in this gene may exert different effects on these two phenotypes.     

Population admixture modulates risk for alcohol dependence

The admixture of different ancestral populations in America may have important implications for the risk for psychiatric disorders, as it appears to have for other medical disorders. The present study investigated the role of population admixture in risk for several psychiatric disorders in European-Americans (EAs) and African-Americans (AAs). This is a multisite study with 3,792 subjects recruited from across the United States, including 3,119 EAs and 673 AAs. These subjects included healthy controls and those with substance dependence (SD) [including alcohol dependence (AD), cocaine dependence, and opioid dependence], social phobia, affective disorders, and schizophrenia. In addition, DNA was included from 78 West Africans. The degree of admixture for each subject was estimated by analysis of a set of ancestry-informative genetic markers using the program STRUCTURE, and was compared between cases and controls. As noted previously, the degree of admixture in AAs was higher than EAs. In EAs, the degree of admixture (with African ancestry) was significantly lower in patients with SD (mainly AD) than controls (P = 0.009 for SD; P = 0.008 for AD). This finding suggests that population admixture may modulate risk for alcohol dependence. Population admixture might protect against alcohol dependence by increasing average heterozygosity and reducing the risk of deleterious recessive alleles. We cannot exclude the possibility that the results might have been influenced by selection bias due to the multisite nature of the study.     

Complex recombination events at the hypermutable minisatellite CEB1 (D2S90).

Some minisatellite structures are the site of high rates of DNA recombination in non-pathological situations, with an excess of motif insertion events and a locus-dependent sex-specific mutation bias. We previously reported the cloning of the hypermutable minisatellite locus CEB1 (D2S90), remarkable for its 13% mutation rate in the male germline (compared to approximately 0.4% in female). We have sought to analyse the mechanisms underlying the addition or deletion of motifs at this locus using the minisatellite variant repeat mapping technique. This is possible with a high precision due to the extreme sequence polymorphism seen between different motifs. No crossing-over event was observed among 38 informative neomutations. Four of the 19 informative mutant alleles with an addition of motifs are interallelic events, the others are intra-allelic. Overall, the insertion and deletion mutations are spread along the alleles, although the subset of interallelic events shows clustering towards the analysed end. The apparently complex recombination events observed can all be interpreted as a succession of elementary duplications-deletions of inter- as well as intra-chromosomal origin, suggesting a model in which sister chromatid as well as conversion-like exchanges are involved in these mutation processes.     

Mutation Induction in the Male Recombination Strains of DROSOPHILA MELANOGASTER

One group of the second chromosome lines isolated from a southern Texas population of Drosophila melanogaster, which has been known to show relatively high frequencies of male recombinations, was found to increase the frequency of sex-linked recessive lethal mutations from a control frequency of 0.18% to 1.63%. The second group, which showed a very much reduced frequency of male recombinations, was found to cause a slight increase to 0.48%, although it was not statistically significant. The first group was also tested for the recessive lethal mutation frequency in the second chromosome; the frequency increased from a control frequency of 0.28% to 2.82%. Mapping of a portion of the sex-linked lethals indicated a distribution along the entire X chromosome, although there was a tendency of clustering towards the tip of the X chromosome. One sex-linked lethal line so far tested was found to be associated with an inversion (approximate breakpoints, 14A-18A). It was suggested that the element causing male recombination might be similar to the hi mutator gene studied earlier by Ives (1950).     

Molecular Characterization and Recombination Analysis of the Complete Genome of Apple Chlorotic Leaf Spot Virus

The complete nucleotide sequence of an Indian isolate of Apple chlorotic leaf spot virus (ACLSV) was determined and found to be 7,525 nt in length. The genome organization was similar to known isolates of ACLSV, encoding three ORFs. Comparisons indicated high sequence variability among known isolates with overall nucleotide sequence identities of 80 to 84%. A striking variable region was identified among the replicase protein upstream of the RNA‐dependent RNA polymerase (aa 1510–1590), which showed a 41–43% match with the corresponding region in other isolates. Phylogenetic analysis at the nucleotide level clustered the isolates into three groups, without any relation to geographical origin. Recombination analysis showed that the isolate is a recombinant with recombination sites spread throughout the genome, especially in the polymerase gene region (nt 4700–5400). Most recombination sites were bordered by an upstream region (5′) of GC‐rich and downstream region (3′) of AU‐rich sequences of similar length. Correlation of recombination site with host type is discussed, and it was found that there were more interlineage recombinations in the apple host compared with intralineage recombinations.     

Recombination as a Point Process along Sequences

Histories of sequences in the coalescent model with recombination can be simulated using an algorithm that takes as input a sample of extant sequences. The algorithm traces the history of the sequences going back in time, encountering recombinations and coalescence (duplications) until the ancestral material is located on one sequence for homologous positions in the present sequences. Here an alternative algorithm is formulated not as going back in time and operating on sequences, but by moving spatially along the sequences, updating the history of the sequences as recombination points are encountered. This algorithm focuses on spatial aspects of the coalescent with recombination rather than on temporal aspects as is the case of familiar algorithms. Mathematical results related to spatial aspects of the coalescent with recombination are derived.     

The effect of ACACB cis-variants on gene expression and metabolic traits

Background

Acetyl Coenzyme A carboxylase β (ACACB) is the rate-limiting enzyme in fatty acid oxidation, and continuous fatty acid oxidation in Acacb knock-out mice increases insulin sensitivity. Systematic human studies have not been performed to evaluate whether ACACB variants regulate gene expression and insulin sensitivity in skeletal muscle and adipose tissues. We sought to determine whether ACACB transcribed variants were associated with ACACB gene expression and insulin sensitivity in non-diabetic African American (AA) and European American (EA) adults.

Methods

ACACB transcribed single nucleotide polymorphisms (SNPs) were genotyped in 105 EAs and 46 AAs whose body mass index (BMI), lipid profiles and ACACB gene expression in subcutaneous adipose and skeletal muscle had been measured. Allelic expression imbalance (AEI) was assessed in lymphoblast cell lines from heterozygous subjects in an additional EA sample (n = 95). Selected SNPs were further examined for association with insulin sensitivity in a cohort of 417 EAs and 153 AAs.

Results

ACACB transcribed SNP rs2075260 (A/G) was associated with adipose ACACB messenger RNA expression in EAs and AAs (p = 3.8×10⁻⁵, dominant model in meta-analysis, Stouffer method), with the (A) allele representing lower gene expression in adipose and higher insulin sensitivity in EAs (p = 0.04). In EAs, adipose ACACB expression was negatively associated with age and sex-adjusted BMI (r = −0.35, p = 0.0002).

Conclusions

Common variants within the ACACB locus appear to regulate adipose gene expression in humans. Body fat (represented by BMI) may further regulate adipose ACACB gene expression in the EA population.     

The Bf locus in the HLA region of chromosome 6: linkage and association studies.

Bf allele frequencies in a material of 172 unrelated Norwegians are given. Bf/HLA linkage relations in 49 informative matings with 178 children, and Bf/HLA association data of a material of 212 Bf-HLA haplotypes are presented. Of 171 informative meioses, there were no Bf-HLA-B recombinations, while 3 out of 158 Bf-HLA-A informative meioses showed recombination. There is significant association between the BfF and the HLA-BW35 allele. It is concluded that the Bf locus is situated on the HLA-B side of HLA-A within the HLA region, in very close proximity to HLA-B.