HLA-DR4-associated haplotypes are genotypically diverse within HLA

Biochemical diversity among products of class II HLA genes has been observed in individuals who appear to be HLA-D and DR-identical by cellular and serologic typing. We used techniques of restriction enzyme fragment analysis by Southern blotting to analyze this diversity at the level of cellular DNA. A panel of 17 HLA-DR4 homozygous cell lines (HCL) were investigated by using cDNA probes homologous to DQ beta, DQ alpha, and DR beta genes. Each probe was hybridized to cellular DNA digested with a series of different restriction endonucleases. Polymorphisms were observed with the use of the enzymes Pst I, Hind III, and Bam HI: Hybridization of cellular DNA digested with Hind III and Pst I with the DQ beta probe revealed specific polymorphisms, as did hybridization of the Pst I digest with the DQ alpha cDNA probe and the Bam HI digest with the DR beta probe. The observed differences fall into two categories: first, considerable diversity was seen between HLA-DR4 HCL that represent different HLA-D-defined haplotypes; second, diversity was also observed among HCL of the same DR4-associated HLA-D cluster. In contrast to the DQ cDNA probes, hybridization with the DR beta probe revealed relatively limited polymorphism by using a panel of different restriction endonucleases. Thus, although there is a general pattern of polymorphic restriction enzyme fragments homologous to DQ probes within an HLA-D cluster, the pattern seen for any particular cell line was not sufficiently distinct to assign an HLA-D or DR specificity.     

Molecular analysis distinguishes two HLA-DR3-bearing major histocompatibility complex extended haplotypes

We detected restriction fragment length polymorphisms that distinguish the extended haplotype HLA-B8,DR3,SCO1 from HLA-B18,DR3,F1C30 at the DR beta and DQ beta loci with five of seven restriction endonucleases used. One set of restriction fragments was always found on HLA-B8,DR3,SCO1 and associated with DRw52a, while the other was present on HLA-B18, DR3,F1C30 and correlated with DRw52b (the gene encoding the subtype of DRw52 associated with the BO1 or LB-Q1 antigen). Furthermore, using a full-length DQ beta gene probe, we found division in the DQw2 haplotype, in which DQw2a always associated with HLA-B8, DR3,SCO1, while DQw2b always occurred with HLA-B18,DR3,F1C3O. Our evidence thus indicates that serologically defined HLA-DR3, HLA-DRw52, and HLA-DQw2 are each produced by two structurally very different sets of genes, one set occurring in HLA-B8, DR3,SCO1, and the other in HLA-B18,DR3,F1C30.Abbreviations used in this paper BSA bovine serum albumin - MHC major histocompatibility complex - EDTA ethylenediaminetetraacetic acid - SDS sodium dodecyl sulfate     

Structural similarity of the HLA-DQ region in DQ3 and DQ4 haplotypes and structural diversity of the HLA-DQ region in HLA-DR7 haplotypes.

Genomic DNA obtained from a B lymphoblastoid cell line was digested with appropriate restriction endonuclease and hybridized with several probes specific for genes encoding HLA-DQ. Southern hybridization with a DQA1 3'untranslated (UT) region probe showed DQ2-type hybridization pattern in DR7DQ3 haplotype. On the contrary, DQB1 3'UT probe showed DQ3-type pattern in the same haplotype. Gene cloning and DNA sequencing analysis revealed a repetitive sequence, (TG)19, between DQA1 and DQB1 gene in the DR7DQ3 haplotype. These results suggest that a recombination event has occurred near this potential Z-DNA structure in the haplotype, DR7DQ3. The 3'UT region probes of DQA1 and DQB1 genes failed to detect restriction fragment length polymorphism (RFLP) differences between DR4DQ3 and DR4DQ4 haplotypes in this experiment, suggesting that the gene structure between DQA1 and DQB1 is conserved in these haplotypes.     

Polymorphisms within the HLA-DRw6 haplotype

We have analyzed the HLA class II gene products from HLA-DRw6 homozygous cells. Epstein-Barr virus-transformed B-cell lines were internally labeled with [³⁵S]-methionine. An NP-40 lysate of the cells was subjected to immunoprecipitation, first with a DRw52-like-specific monoclonal antibody and subsequently with a DR-specific framework antibody. The DR region-encoded gene products were analyzed by one-dimensional gel isoelectric focusing and two-dimensional gel electrophoresis. It is shown that DRw6 homozygous cell lines contain at least two nonallelic DR chains, one carrying a DRw52 determinant and one DRw52-negative population. Both chains appear to be polymorphic between the cellularly defined subtypes of DRw6. The determinant responsible for the differential mixed lymphocyte culture reactivity of Dw18 and Dw19 cells resides on the DRw52-positive population, whereas the Dw6-Dw9 differences are attributed to determinants on both populations of DR light chains. The Dw16-derived DRw52⁺ chain much resembles the Dw18 DRw52⁺ light chain whereas there is a clear-cut difference between these two subtypes in the DRw52^– population. We conclude that, for DRw6 homozygous cells, the cellularly recognized D determinants are probably located on DR-encoded molecules, both DRw52⁺ and DRw52^–, and that charge shift of these chains is at least partly responsible for differential recognition of these cells in mixed lymphocyte cultures.Abbreviations used in this paper: MLC mixed lymphocyte culture - 1D-IEF one-dimensional isoelectric focusing - 2D two-dimensional - moab monoclonal antibody     

Coeliac disease patients carry conserved HLA-DR3-DQ2 haplotypes revealed by association of TNF alleles

Certain HLA-DQ alleles are known to contribute to predisposition to coeliac disease (CD). The existence of additional independent risk-modifying loci in the HLA complex is still being debated. The DR3-DQ2 haplotype has been studied most, but the evidence is conflicting. The discrepancies may stem from the absence of such an effect, insufficient statistical power to detect an effect (i.e. small studies) and/or incomplete control of linkage disequilibrium (LD) to the neighbouring DQ-loci, known to elicit a strong effect. In the present study, we aimed to undertake a statistically high-powered family-based analysis, fully controlling effects of LD between the major DQ-risk haplotypes and neighbouring candidate loci. We investigated five markers on DR3-DQ2, DR5-DQ7 and DR7-DQ2 haplotypes in 327 Norwegian and Swedish families. Our primary finding was that TNF-308A (TNF2) was significantly associated on the DR3-DQ2 haplotype [stratum specific odds ratio (OR)=2.40 (1.25–4.48), Pc=0.009, where P_c=Pn and n=number of tests performed]. Furthermore, we confirmed earlier indications that LD between TNF2 and DQA1*05-DQB1*02 on the DR3 haplotype is more strongly maintained in family-based cases than family-based controls. In conclusion, we confirmed in this study, the largest of its kind, that additional CD risk factors independent of DQ2 alleles do exist on the DR3 haplotype.     

Polymorphisms distinguishing different mouse species and t haplotypes.

Three anonymous chromosome 17 DNA markers, D17Tu36, D17Tu43, and D17Le66B, differentiate between house mouse species and/or between t chromosomes. The D17Tu36 probe, which maps near the Fu locus and to the In(17)4 on t chromosomes, identifies at least 15 haplotypes, each haplotype characterized by a particular combination of DNA fragments obtained after digestion with the Taq I restriction endonuclease. Ten of these haplotypes occur in Mus domesticus, while the remaining five occur in M. musculus. In each of these two species, one haplotype is borne by t chromosomes while the other haplotypes are present on non-t chromosomes. The D17Tu43 probe, which maps near the D17Leh122 locus and to the In(17)3 on t chromosomes, also identifies at least 15 haplotypes in Taq I DNA digests, of which nine occur in M. domesticus and six in M. musculus. One of the nine M. domesticus haplotypes is borne by t chromosomes, the other haplotypes are borne by non-t chromosomes; two of the six M. musculus haplotypes are borne by t chromosomes and the remaining four by non-t chromosomes. Some of the D17Tu43 haplotypes are widely distributed in a given species, while others appear to be population-specific. Exceptions to species-specificity are found only in a few mice captured near the M. domesticus-M. musculus hybrid zone or in t chromosomes that appear to be of hybrid origin. The D17Leh66B probe, which maps to the In(17)2, distinguishes three haplotypes of M. domesticus-derived t chromosomes and one haplotype of M. musculus-derived t chromosomes. Because of these characteristics, the three markers are well suited for the study of mouse population genetics in general and of t chromosome population genetics in particular. A preliminary survey of wild M. domesticus and M. musculus populations has not uncovered any evidence of widespread introgression of genes from one species to the other; possible minor introgressions were found only in the vicinity of the hybrid zone. Typing of inbred strains has revealed the contribution of only M. domesticus DNA to the chromosome 17 of the laboratory mouse.     

Dissection of HLA class II haplotypes in HLA-DR4 homozygous individuals

In order to identify better markers for HLA-DR4-associated autoimmune disorders, we have studied the complexity of the HLA class II region in DR4-positive cells at the DNA level and compared the DNA polymorphism with that defined by serology, mixed lymphocyte culture (MLC) reactivity, and protein chemistry. At the DNA level, HLA-DR4 can be characterized by a homogeneous pattern of bands hybridizing to HLA class II cDNA probes. Besides, subtypes can be defined within DR4 using HLA-DR , -DQ , and -DQ cDNA probes in Southern blot analysis. Three subtypes are found using the DR cDNA probe. One of these subtypes correlates with the cellularly defined Dw15 specificity, another with the serologically defined LB4 and LB14 specificities. None of the restriction fragment length polymorphism (RFLP) patterns coincide with the MLC-defined DR4 subtypes Dw4, DW10, Dw13, and Dw14 separately. Variation of two fragments hybridizing to the DQ cDNA probe obtained after either Pvu II or Taq I digestion yields three subtypes. Pvu II- and Eco RI-digested DR4 DNA give rise to three DQ detectable subtypes. Correlation between these subtypes, isoelectric point variation of DQ molecules, and the DQ-related allelic system TA10/2B3 are demonstrated. Some of the patterns obtained with DQ and DQ cDNA probes display heterozygosity in the DQ region, as demonstrated by family segregation. No correlation was observed between DQ and the cellularly defined Dw determinants. A new polymorphism has been obtained with the DQ probe, probably due to DX polymorphism. DR RFLP divides the LB 14 supertypic specificity into two new subtypes. A combination of the four different techniques applied to a panel of 16 DR4 homozygous cell lines reveals at least nine different haplotypes in DR4. These newly defined haplotypes may be of help in further studies concerning the relationship of micropolymorphism with several diseases.     

Susceptibility locus for IgA deficiency and common variable immunodeficiency in the HLA-DR3, -B8, -A1 haplotypes.

BACKGROUND: A common genetic basis for IgA deficiency (IgAD) and common variable immunodeficiency (CVID) is suggested by their occurrence in members of the same family and the similarity of the underlying B cell differentiation defects. An association between IgAD/CVID and HLA alleles DR3, B8, and A1 has also been documented. In a search for the gene(s) in the major histocompatibility complex (MHC) that predispose to IgAD/CVID, we analyzed the extended MHC haplotypes present in a large family with 8 affected members. MATERIALS AND METHODS: We examined the CVID proband, 72 immediate relatives, and 21 spouses, and determined their serum immunoglobulin concentrations. The MHC haplotype analysis of individual family members employed 21 allelic DNA and protein markers, including seven newly available microsatellite markers. RESULTS: Forty-one (56%) of the 73 relatives by common descent were heterozygous and nine (12%) were homozygous for a fragment or the entire extended MHC haplotype designated haplotype 1 that included HLA- DR3, -C4A-0, -B8, and -A1. The remarkable prevalence of haplotype 1 was due in part to marital introduction into the family of 11 different copies of the haplotype, eight sharing 20 identical genotype markers between HLA-DR3 and HLA-B8, and three that contained fragments of haplotype 1. CONCLUSION: Crossover events within the MHC indicated a susceptibility locus for IgAD/CVID between the class III markers D821/D823 and HLA-B8, a region populated by 21 genes that include tumor necrosis factor alpha and lymphotoxins alpha and beta. Inheritance of at least this fragment of haplotype 1 appears to be necessary for the development of IgAD/CVID in this family.     

The architecture of long-range haplotypes shared within and across populations

Homologous long segments along the genomes of close or remote relatives that are identical by descent (IBD) from a common ancestor provide clues for recent events in human genetics. We set out to extensively map such IBD segments in large cohorts and investigate their distribution within and across different populations. We report analysis of several data sets, demonstrating that IBD is more common than expected by na?ve models of population genetics. We show that the frequency of IBD pairs is population dependent and can be used to cluster individuals into populations, detect a homogeneous subpopulation within a larger cohort, and infer bottleneck events in such a subpopulation. Specifically, we show that Ashkenazi Jewish individuals are all connected through transitive remote family ties evident by sharing of 50 cM IBD to a publicly available data set of less than 400 individuals. We further expose regions where long-range haplotypes are shared significantly more often than elsewhere in the genome, observed across multiple populations, and enriched for common long structural variation. These are inconsistent with recent relatedness and suggest ancient common ancestry, with limited recombination between haplotypes.     

Discovery of multiple IGS haplotypes within genotypes of Puccinia striiformis

In a search for specific molecular markers for population analysis of Puccinia striiformis f. sp. tritici, the ribosomal DNA (rDNA) intergenic spacer (IGS) 1 region (rDNA-IGS1, between the 28S and the 5S rDNA genes) was amplified, cloned, and sequenced. It was found to exhibit multiple bands and length polymorphism. Surprisingly, single isolates were found to possess between three to five different IGS1 haplotypes. Bands were cloned and sequenced, and two highly variable regions (α and β) were found between conserved regions, with repeat units interspersed in both types of regions. There were 14 different repeat units, and these were sometimes grouped further into four combinations of repeat units, with a few individual nucleotides (A or C) inserted between the repeats. Among three geographically dispersed isolates, the variable region α was divided into eight types, and the variable region β was divided into two types based on repeat units. Most of the 14 repeat units were shared by the variable and the conserved regions. Among the three isolates, there were a total of 12 IGS1 haplotypes, but some of these were shared between isolates such that there were only eight unique haplotypes. The occurrence of multiple haplotypes within single isolates may be useful for analyzing the population structure, tracking the origin of annual epidemics and providing insights into evolutionary biology of this pathogen.     

The evolutionary origin of the HLA-DR3 haplotype

The human HLA-DR3 haplotype consists of two functional genes (DRB1*03 and DRB3*01) and one pseudogene (DRB2), arranged in the order DRB1... DRB2... DRB3 on the chromosome. To shed light on the origin of the haplotype, we sequenced 1480 nucleotides of the HLA-DRB2 gene and aong stretches of two other genes, Gogo-DRB2 from a gorilla, Sylvia and Patr-DRB2 from a chimpanzee, Hugo. All three sequences (HLA-DRB2, Gogo-DRB2, Patr-DRB2) are pseudogenes. The HLA-DRB2 and Gogo-DRB2 pseudogenes lack exon 2 and contain a twenty-nucleotide deletion in exon 3, which destroys the correct translational reading frame and obliterates the highly conserved cysteine residue at position 173. The Patr-DRB2 pseudogene lacks exons 1 and 2; it does not contain the twenty-nucleotide deletion, but does contain a characteristic duplication of that part of exon 6 which codes for the last four amino acid residues of the cytoplasmic region. When the nucleotide sequences of these three genes are compared to those of all other known DRB genes, the HLA-DRB2 is seen as most closely related to Gogo-DRB2, indicating orthologous relationship between the two sequences. The Patr-DRB2 gene is more distantly related to these two DRB2 genes and whether it is orthologous to them is uncertain. The three genes are in turn most closely related to HLA-DRBVI (the pseudogene of the DR2 haplotype) and Patr-DRB6 (another pseudogene of the Hugo haplotype), followed by HLA-DRB4 (the functional but nonpolymorphic gene of the DR4 haplotype). These relationships suggest that these six genes evolved from a common ancestor which existed before the separation of the human, gorilla, and chimpanzee lineages. The DRB2 and DRB6 have apparently been pseudogenes for at least six million years (myr). In the human and the gorilla haplotype, the DRB2 pseudogene is flanked on each side by what appear to be related genes. Apparently, the DR3 haplotype has existed in its present form for more than six myr.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M86691–94.     

Polymorphisms and haplotypes in methylenetetrahydrofolate reductase gene and head and neck squamous cell carcinoma risk

Functional polymorphisms in genes encoding enzymes involved in folate metabolism might modulate head and neck carcinoma risk because folate participates in DNA methylation and synthesis. We therefore conducted a case–control study of 853 individuals (322 head and neck cancer cases and 531 non-cancer controls) to investigate associations among MTHFR C677T and MTHFR A1298C polymorphisms and head and neck squamous cell carcinoma risk. Interactions between these two polymorphisms and risk factors and clinical histopathological parameters were also evaluated. The polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) technique was used to genotype the polymorphisms and Chi-square test and multiple logistic regression were used for statistical analyses. The variables age ≥49 years, male gender, tobacco habits and alcohol consumption, MTHFR 1298 AC or CC genotypes, combined genotypes with two or more polymorphic alleles and 677T and 1298C polymorphic alleles were associated with increased risk for this disease (P < 0.05). Furthermore, we found that 1298 AC or CC genotypes were associated with age ≥49 years, tobacco and alcohol habits (P < 0.05). Regarding clinical histopathological parameters, the A1298C polymorphism was more frequent in patients with oral cavity as primary site (P < 0.05). MTHFR polymorphisms may contribute for increase risk for head and neck carcinoma and the variables age ≥49 years, male gender, tobacco and alcohol habits were associated with MTHFR 1298AC or CC genotypes, confirming that individuals with these variables and MTHFR A1298C polymorphism has higher risk for this disease.     

Ten HLA-DR4 alleles defined by sequence polymorphisms within the DRB1 first domain

We have studied DRB1 sequence polymorphisms associated with DR4 subtypes using DR4-specific DNA amplification and sequence-specific oligonucleotide probe (SSOP) hybridization. The 5 amplification primer was designed to hybridize with a unique sequence in the first hypervariable region (HVR) of the DRB1 second ex-on of all known DR4 alleles; the 3 primer was designed to hybridize with an intron sequence common to all DRB1 alleles. The specificity of the amplification step was demonstrated by double-blind testing of 105 selected DNA samples. Prospective SSOP typing of DR4 alleles was performed in 104 unrelated individuals known to be DR4-positive, including 13 who were DR4-homozygous. A DRB1 subtype corresponding with the previously defined DR4-associated specificities Dw4, Dw10, Dw13.1, Dw13.2, Dw14.1, Dw14.2, Dw15, and DwKT2 could be assigned for each of the 117 DR4 haplotypes tested. In most cases, DR4-homozygous, DRB1-heterozygous individuals could be genotyped with the panel of probes. In the course of our analysis, we identified two new DR4-related alleles, DRB1*04.CB (DRB1*0410)¹ and DRB1*04.EC (DRB1*, 0411)² which were recognized by their novel hybridization patterns. The DRB1 second exon sequence of DRB1*04.CB, is identical to DRB1*0405 except at codon 86 where GTG encodes valine instead of GGT encoding glycine. DRB1*04.EC is identical to DRB1*04.CB except at codon 74 where GAG encodes glutamic acid instead of GCG encoding alanine. Our results provide further evidence that SSOP hybridization is the most effective approach available for subtyping DR4 haplotypes and identifying unrecognized variants. A similar approach should be equally informative for subtyping other DR alleles.     

Polymorphisms and haplotypes in MyD88 are associated with the development of sarcoidosis: a candidate-gene association study

Sarcoidosis is considered as a disorder of protracted immune response to an as yet unidentified causative agent that leads to granuloma formation. Material from M. tuberculosis and P. acne has been repeatedly detected in the sarcoidosis lesions, implying the involvement of the Toll-like receptor2 (TLR2) gene that responds to these intracellular pathogens. Since TLR2 association studies have produced controversial results, we sought to investigate whether the downstream signalling molecule MyD88 could be linked to disease susceptibility. We analyzed a total of 93 cases with sarcoidosis and of 89 controls for the most common MyD88 SNPs: ?938C>A (rs4988453) and 1944C>G (rs4988457). There is evidence that the genotype distributions of both variants are associated with the development of sarcoidosis (p = 0.038 for ?938C>A and p = 0.026 for 1944C>G). In particular, ?938A and 1944G carriers were associated with risk of sarcoidosis [OR = 2.48 (1.23–5.02) and OR = 0.33 (0.14–0.76)], respectively, indicating dominance of the mutant alleles; however, the adjustment of the effect size for age and sex diminished the significance. The haplotype analysis showed association for the ?938A/1944G haplotype (p < 0.001). Since genetic association studies have linked MyD88 to Hodgkin’s lymphoma it is tempting to speculate that MyD88 may contribute to the granuloma formation that characterizes sarcoidosis.     

Quantitative Analysis of Single Nucleotide Polymorphisms within Copy Number Variation

Background

Single nucleotide polymorphisms (SNPs) have been used extensively in genetics and epidemiology studies. Traditionally, SNPs that did not pass the Hardy-Weinberg equilibrium (HWE) test were excluded from these analyses. Many investigators have addressed possible causes for departure from HWE, including genotyping errors, population admixture and segmental duplication. Recent large-scale surveys have revealed abundant structural variations in the human genome, including copy number variations (CNVs). This suggests that a significant number of SNPs must be within these regions, which may cause deviation from HWE.

Results

We performed a Bayesian analysis on the potential effect of copy number variation, segmental duplication and genotyping errors on the behavior of SNPs. Our results suggest that copy number variation is a major factor of HWE violation for SNPs with a small minor allele frequency, when the sample size is large and the genotyping error rate is 0∼1%.

Conclusions

Our study provides the posterior probability that a SNP falls in a CNV or a segmental duplication, given the observed allele frequency of the SNP, sample size and the significance level of HWE testing.     

Armenian Y chromosome haplotypes reveal strong regional structure within a single ethno-national group

Armenia has been little-studied genetically, even though it is situated in an important area with respect to theories of ancient Middle Eastern population expansion and the spread of Indo-European languages. We screened 734 Armenian males for 11 biallelic and 6 microsatellite Y chromosome markers, segregated them according to paternal grandparental region of birth within or close to Armenia, and compared them with data from other population samples. We found significant regional stratification, on a level greater than that found in some comparisons between different ethno-national identities. A diasporan Armenian sub-sample (collected in London) was not sufficient to describe this stratified haplotype distribution adequately, warning against the use of such samples as surrogates for the non-diasporan population in future studies. The haplotype distribution and pattern of genetic distances suggest a high degree of genetic isolation in the mountainous southern and eastern regions, while in the northern, central and western regions there has been greater admixture with populations from neighbouring Middle Eastern countries. Georgia, to the north of Armenia, also appears genetically more distinct, suggesting that in the past Trans-Caucasia may have acted as a genetic barrier. A Bayesian full-likelihood analysis of the Armenian sample yields a mean estimate for the start of population growth of 4.8 thousand years ago (95% credible interval: 2.0-11.1), consistent with the onset of Neolithic farming. The more isolated southern and eastern regions have high frequencies of a microsatellite defined cluster within haplogroup 1 that is centred on a modal haplotype one step removed from the Atlantic Modal Haplotype, the centre of a cluster found at high frequencies in England, Friesland and Atlantic populations, and which may represent a remnant paternal signal of a Paleolithic migration event.