Comparative analysis of sequence variability in the upstream regulatory region of the HLA-DQB1 gene

Positive selection in coding regions and motif duplication in regulatory regions of bottlenose dolphin MHC class II genes

The vertebrate immune response is mediated through highly adaptive, quickly evolving cell surface receptors, the major histocompatibility complex (MHC). MHC molecules bind and present a diverse array of pathogenic molecules and trigger a cascade of defenses. Use of MHC variation as a marker for population health has also evolved quickly following advances in sequencing methods. We applied a combination of traditional and next generation sequencing methodology to characterize coding (peptide binding region) and regulatory (proximal promoter) sequence variation in MHC Class II DQA and DQB genes between estuarine and coastal populations of the bottlenose dolphin, Tursiops truncatus, an apex predator whose health status is indicative of anthropogenic impacts on the ecosystem. The coding regions had 10 alleles each at DQA and DQB; the promoters had 6 and 7 alleles at DQA and DQB, respectively with variation within key regulatory motifs. Positive selection was observed for the coding regions of both genes while both coding and promoter regions exhibited geographic differences in allele composition that likely indicates diversifying selection across habitats. Most notable was the discovery of a complete duplication of a 14-bp T-box motif in the DQA promoter. Four class II promoter regions (DQA, DQB, DRA, DRB) were characterized in species from four cetacean families (Delphinidae, Monodontidae, Lipotidae, and Physeteridae) and revealed substantial promoter structural diversity across this order. Peptide binding regions may not be the only source of adaptive potential within cetacean MHC for responding to pathogenic threats. These findings are the first analysis of cetacean MHC regulatory motifs, which may divulge unique immunogenetic strategies among cetaceans and reveal how MHC transcriptional control continues to evolve. The combined MHC regulatory and coding data provide new genetic context for distinct vulnerability profiles between coastal and estuarine populations, which are key concerns for health and risk management.

Association of genetic polymorphisms and autoimmune Addison's disease

Autoimmune Addison's disease (AAD) is a complex genetic disease that results from the interaction of a predisposing genetic background with as yet unknown environmental factors. The disease is marked by the appearance of circulating autoantibodies against steroid 21-hydroxylase. Mutations of the autoimmune regulator gene are responsible for the so-called autoimmune polyendocrine syndrome type I (APS I), of which AAD is a major disease component. Among genetic factors for isolated AAD and APS II, a major role is played by HLA class II genes: HLA-DRB1 0301-DQA1 0501-DQB1 0201 and DRB1 04-DQA1 0301-DQB1 0302 are positively, and RB1 0403 is negatively, associated with a genetic risk for AAD. The MHC class I chain-related gene A allele 5.1 is strongly and positively associated with AAD. Other gene polymorphisms contributing to genetic risk for AAD are MHC2TA, the gene coding for class II transactivator, the master regulator of class II expression, cytotoxic T lymphocyte antigen-4, PTPN22 and the vitamin D receptor.

Differential expression of HLA-DQ alleles in peripheral blood mononuclear cells: alleles associated with susceptibility to and protection from autoimmune type 1 diabetes

Differential expression of human leucocyte antigen (HLA) class II genes has been postulated to influence the risk of developing autoimmune disease. In this study, we investigated the relationship between the level of mRNA expression of DQA1 and DQB1 alleles in peripheral blood mononuclear cells and the influence of the alleles on susceptibility to type 1 diabetes (T1D). Transcripts from pairs of DQA1 and DQB1 alleles were quantified in 59 DQ-heterozygous individuals (29 patients with T1D and 30 healthy control subjects). Luciferase reporter gene assays were used to investigate the relative promoter activities of the alleles associated with high and low risk of disease. DQA1*0301 and the DQB1*06 group of alleles (*0601, *0602, *0603 and *0604) were generally overexpressed in comparison to other alleles. In contrast, mRNA for DQB1*0201/*0202 was generally less abundant than other DQB1 transcripts. These data correlated well with the relative promoter activities observed for the diabetes-associated alleles; the strongest promoters were those derived from DQA1*0301 and DQB1*0602, while a 700-bp fragment derived from the DQB1*0201 promoter showed the lowest activity of the DQB1 constructs. There was no simple correlation between the level of expression of specific DQ alleles and their influence on the risk of diabetes. The functional relevance of our findings and their implications for the pathogenesis of autoimmunity remain to be determined.

Sequence variations in the transcriptional regulatory region and intron1 of HLA-DQB1 gene and their linkage in southern Chinese ethnic groups

Sequence polymorphism in the transcriptional regulatory region extending to intron1 (DQRRI1) of HLA-DQB1 gene, and their haplotypic distributions were investigated in southern Chinese populations. We cloned and sequenced a 1.1-kb segment containing the transcriptional regulatory region, exon1, and partial intron1 of HLA-DQB1 gene of 37 individuals from nine different ethnic groups in southern China. A high-density map of 162 polymorphisms, including 152 single nucleotide polymorphisms (SNPs) and 10 insertion-deletion polymorphisms, was revealed. By comparing these data with SNPs deposited in dbSNP database in National Center for Biotechnology Information and polymorphisms that have been reported, 66 polymorphisms were firstly reported. A total of 16 different haplotypes were detected based on these 162 polymorphisms. The distribution of 16 alleles of DQRRI1 as well as their linkage with DQB1 exon2 alleles was also investigated based on the population study and phylogenetic analysis. Tight linkage between these two regions were discovered, as each of DQB1*02, DQB1*03, DQB1*04, DQB1*05, and DQB1*06 alleles was seen to be linked with specific DQRRI1 allele. Our study showed different pattern of transcriptional regulatory region, exon1, and intron1 of different DQB1 alleles.

Polymorphisms of the upstream regulatory region of the major histocompatibility complex DRB genes in domestic horses

Sequence information was obtained on the variation of the ELA-DRB upstream regulatory region (URR) after polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) cloning and sequencing of approximately 220 bp upstream of the first exon of horse DRB genes. The sequence of the proximal URR of equine DRB is composed of highly conserved sequence motifs, showing the presence of the W, X, Y, CAAT and TATA conserved boxes of major histocompatibility complex (MHC) class II promoters. Five different polymorphic horse DRB promoter sequences were detected in five horse breeds. The results demonstrate the existence of polymorphism in the nucleotide sequences of the ELA-DRB URR, located in the functionally important conserved consensus sequences, the X2 box, the Y box and the TATA box, while conservation were observed in X1 and CAAT boxes. The nucleotide diversity among horse URRs was intermediate between that seen within human and mouse DRB promoters, suggesting the existence of another important source of variability in ELA-DRB genes. In addition, phylogenetic comparisons, identity analysis and sequence organization suggested that the reported sequences would correspond to an expressed ELA-DRB locus. However, further information about the functional significance of these promoter polymorphisms will probably be acquired through expression studies on the different sequences.

MHC promoter polymorphism in grey wolves and domestic dogs

A functional immune system requires a tight control over major histocompatibility complex (MHC) gene transcription, as the abnormal MHC expression patterns of severe immunodeficiency and autoimmune diseases demonstrate. Although the regulation of MHC expression has been well documented in humans and mice, little is known in other species. In this study, we detail the level of polymorphism in wolf and dog MHC gene promoters. The promoter regions of the DRB, DQA and DQB locus were sequenced in 90 wolves and 90 dogs. The level of polymorphism was high in the DQB promoters, with variation found within functionally relevant regions, including binding sites for transcription factors. Clear associations between DQB promoters and exon 2 alleles were noted in wolves, indicating strong linkage disequilibrium in this region. Low levels of polymorphism were found within the DRB and DQA promoter regions. However, a variable site was identified within the T box, a TNF-alpha response element, of the DQA promoter. Furthermore, we identified a previously unrecognised 18-base-pair deletion within exon 1 of the DQB locus.

Allele-specific quantification of HLA-DQB1 gene expression by real-time reverse transcriptase-polymerase chain reaction

In addition to coding region polymorphism, allele-specific variation in the upstream regulatory region of the HLA-DQB1 gene has been detected. Reporter gene assays and transfection studies have indicated that HLA-DQB1 promoter polymorphism may be of functional significance. The aim of this study was to utilize real-time reverse transcriptase-polymerase chain reaction (RT-PCR) for allele-specific quantification of HLA-DQB1 expression and to analyze cell-specific HLA-DQB1 expression in vivo. For the allele-specific quantification of DQB1 gene products, a real-time RT-PCR set of primer pairs (n=27) and probes (n=5) targeting exon 2 variability was established. The robustness and integrity of the assay system were confirmed by using recombinant DQB1 exon 2 plasmid clones as active exogenous controls. Sensitivity and reproducibility were assessed by serial dilution and allelic mixing analyses. In application to the study of allele-specific expression of DQB1 gene products during cytokine-driven maturation of monocyte-derived dendritic cells, differential patterns of allelic expression in heterozygous individuals were observed for DQB1*0301, compared to DQB1*0501 and DQB1*0602. At maximum, 1.9-fold (*0301/*0501) and 2.5-fold (*0301/*0602) higher induction was seen for DQB*0301. In conclusion, HLA-DQB1 expression can be analyzed by real-time RT-PCR suitable for cell- and allele-specific detection of HLA-DQB1 transcripts in homo- and heterozygous combinations.

Polymorphism in the bovine BOLA-DRB3 upstream regulatory regions detected through PCR-SSCP and DNA sequencing

In the present work, we describe through polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) and DNA sequencing the polymorphism within the URR-BoLA-DRB3 in 15 cattle breeds. In total, seven PCR-SSCP defined alleles were detected. The alignment of studied sequences showed six polymorphic sites (four transitions, one transversion and one deletion) in the interconsensus regions of the BoLA-DRB3 upstream regulatory region (URR), while the consensus boxes were invariant. Five out of six detected polymorphic sites were of one nucleotide substitution in the interconsensus regions. It is expected that these mutations do not affect significantly the level of expression. In contrast, the deletion observed in the sequence between CCAAT and TATA boxes could have some effect on affinity interactions between the promoter region and the transcription factors. The URR-BoLA-DRB3 DNA analyzed sequences showed moderate level of nucleotide diversity, high level of identity among them and were grouped in the same clade in the phylogenetic tree. In addition, the phylogenetic tree, the similarity analysis and the sequence structure confirmed that the fragment analyzed in this study corresponds to the URR-BoLA-DRB3. The functional role of the observed polymorphic sites among the regulatory motifs in bovine needs to be analyzed and confirmed by means of gene expression assays.

Unbalanced amounts of HLA-DQA1 allele mRNA: DQA1*03 shows high and DQA1*0501 low amounts of mRNA in heterozygous individuals

Genes of the HLA-DR, DQ region confer strong susceptibility to type 1 diabetes mellitus (IDDM). A possible mechanism of susceptibility is a difference in the amounts of transcripts of predisposing and neutral or protective haplotypes. In this study we developed an assay to compare the amounts of mRNA of two distinct HLA-DQA1 alleles in peripheral blood lymphocytes (PBLs) of heterozygous individuals, using a quantitative RT-PCR with an internal standard covering all HLA-DQA1 specifities. We also developed an algorithm to calculate the amounts of mRNA for two distinct alleles in heterozygous individuals based on the comparison to the same internal standard. In total, 37 HLA-DQA1 heterozygous individuals were analysed, including patients with IDDM (n = 14) and healthy controls (n = 23). Intra-individually, we observed different amounts of mRNA for different HLA-DQA1 alleles in the order: HLA-DQA1*03 > *01 > *0201 > *05. This order was observed in all individuals. We also observed a variation in the ratio of these unbalanced amounts of mRNA in individuals with the same HLA-DQA1 allele combinations. In all allele combinations the average ratio was increased in patients with IDDM compared to the control samples. HLA-DQA1*03 positive and DQA1*03, *05 heterozygous patients had the highest average ratios. Nevertheless, based on limited sample numbers, these differences did not reach significance. We therefore conclude that variations between HLA-DQA1 alleles are not limited to the nucleotide sequence but are also found at the level of amounts of mRNA.

DQA1 and DQB1 promoter diversity and linkage disequilibrium with class II haplotypes in Mexican Mestizo population

The upstream sequences in the 5' flanking region of HLA class II genes, regulate their expression and contribute to the development of immunological diseases. We analyzed 105 healthy unrelated Mexican Mestizos for QAP and QBP polymorphism. DNA typing for DRB1, DQA1, DQB1, QAP1 and QBP1 was done using a standardized PCR-SSOP. Although all QAP alleles previously described were found in Mexicans, the distribution differed as compared to other populations. QAP-3.1, 4.1 and 4.2 were the most frequent alleles and were associated with DQA1*03, *0501 and *0402 respectively. The prevalent QBP alleles were 3.21, 3.1 and 4.1 found mainly associated with DQB1*0302, *0301 and *0501. Linkage disequilibria between the promoter and the corresponding DQA1 and DQB1 allele, are in general the same as described by others. A total of 61 different haplotypes were defined, only six of them with a frequency above 4%. The haplotypes DRB1*0407-QAP-3.1-DQA1*03-QBP-3.21-DQB1*0302 (HF = 14.37%) and DRB1*0802-QAP-4.2-DQA1*0401-QBP-4.1-DQB1*0402 (HF = 14.22%), which have an Amerindian ancestry, are the most frequent in Mexicans. Some rare combinations were detected such as DRB1*0405-QAP-1.3-DQA1*0101/4-QBP-5.11/5.12-DQB1*0501 and DRB1*0403-QAP-3.2-DQA1*03-QBP-3.21-DQB1*0302, probably due to ancient recombination events. This knowledge is relevant as a basis to evaluate functional implications and to explore the role of promoter diversity in disease expression.

Different binding of NF-Y transcriptional factor to DQA1 promoter variants

Polymorphism in the HLA-DQA1 promoter (QAP) sequences could influence the gene expression through a differential binding of transcriptional factors. Considering the main role played by the Y-box in the transcription, we focused on the QAP4 variants differing for a G vs A transition from the QAP Y-box consensus sequence. Electrophoretic Mobility Shift Assay using the two Y-box sequences was performed to determine whether this mutation could be reflected in an allele-specific binding of transcriptional factors. Indeed, the NF-Y specific band, recognised by supershift experiments, was clearly observed using the Y-box consensus probe but it was barely detectable with the QAP4 one. On the contrary, two other complexes were found to more strongly interact with QAP4 Y-box in comparison to the consensus sequence. The analysis of a selected panel of HLA homozygous lymphoblastoid cell lines by competitive RT-PCR and by Northern blotting revealed that the DQA1 *0401, *0501,*0601 alleles regulated by the QAP4 promoters were less expressed at the mRNA level than the DQA1* 0201 allele regulated by the QAP2.1 variant. In conclusion, these results show an evident reduction of NF-Y binding to the mutated QAP4 Y-box and a decreased mRNA accumulation of the DQA1 alleles regulated by these variants.

Highly polymorphic promoter regions of HLA DQA1 and DQB1 genes do not help to further define disease susceptibility in insulin-dependent diabetes mellitus

HLA DQA1, HLA DQB1 genes confer susceptibility to insulin-dependent (type 1) diabetes mellitus (IDDM). Since variants of their upstream regulatory regions are linked to the exons, we investigated their promoter polymorphisms (QAP and QBP) by a combination of PCR-based typing protocols in 136 IDDM patients, 167 controls and 6 families with an IDDM proband to identify possible additional susceptibility markers. Of major interest for IDDM susceptibility are the promoter "splits" of HLA DQA1*0301 (QAP3.1 and QAP3.2) and HLA DQB1*0302 (QBP3.2 and QBP3.3). QAP 3.1 (96% in patients vs 98% in controls) and QBP3.2 (100% vs 99%) were found to be the most frequent promoter variants for HLA DQA1*0301 and DQB1*0302, respectively, whereas QAP3.2 and QBP3.3 were very rare. Furthermore the promoter "splits" were equally distributed on the respective exon alleles in all groups and cosegregated in families as expected. In conclusion, HLA DQ-mediated susceptibility and protection in IDDM is not restricted to the exon but extends to the promoter region without further defining the genetic risk.

Functional significance of polymorphism among MHC class II gene promoters

The functional significance of polymorphism among MHC class II promoters in man and mouse is here reviewed, mainly in terms of the hypothesis of differential expression. The hypothesis proposes that differences between antigen-presenting cells in MHC class II expression exert a co-dominant effect on the Th1-Th2 cytokine balance, such that class II molecules of one type come to control to a greater extent the production of one group of cytokines, and those of another type the production of the alternative group. The survey deals with the influence of signal strength and antigen-presenting cell type on T-cell subset differentiation; functional differences between MHC class II molecules not obviously related to determinant selection; disease protection mediated by HLA alleles; mechanisms possibly responsible for allotypic and isotypic bias; overdominance (heterozygous advantage) in selection for expression of class II alleles; MHC class II promoter structure and function; inter-locus and inter-allele variability within human MHC class II gene upstream regulatory regions; a comparison of these polymorphisms in mouse and man; read-out of class II promoter function; and a comparison with expression of MHC class I. We conclude that the evidence that this variation is functionally active (i.e. controls expression) is increasing, but is not yet compelling. The crucial test still to come, we suggest, is whether or not the biological effects attributable to this polymorphism will line up with molecular studies on expression.

DQB1 promoter sequence variability and linkage in caucasoids

Sequence variability in the upstream regulatory regions (URR) of HLA class II genes has been described as an additional mechanism of diversity in these polymorphic genes. For HLA-DQB1, 12 URR variants have been identified previously by sequence analysis of approx. 600 bp located immediately upstream of the first exon of the DQB1 gene. To investigate the distribution of these promoter alleles and their linkage with the structural portion of the DQB1 gene, a population-based study was carried out. Sequence information was utilized to develop 25 sequence-specific oligonucleotide probes to analyze enzymatically amplified locus-specific DNA. Supplemented with one sequence-specific primer pair to differentiate QBP1-6.2 from -6.3, all known 12 QBP1 alleles could be identified. Subsequently, 215 healthy, unrelated German controls were investigated for the distribution and linkage of DQB1 and QBP1 alleles. A total of 10 out of 12 known QBP1 alleles were observed. Since there was tight linkage between the promoter region and exon 2 of DQB1, the phenotype and genotype frequencies of the promoter alleles corresponded by and large to the frequencies observed for their linked DQB1 alleles. Exceptions were mainly seen for DQ5 and DQ6 haplotypes, as single DQB1 alleles could be linked to different, however, closely related QBP1 alleles and vice versa. Interestingly, for each DQB1 allele a single DQB1/QBP1 haplotype dominated (75.9 to 96.4%) the distribution. It is concluded that promoter and coding region variability are tightly linked by linkage disequilibrium. Exceptions are restricted to DQB1 DQ5 and DQ6 haplotypes. Since functional differences between different QBP1 alleles exist, the maintenance of haplotypic integrity may be of functional importance.

Ultrastructural allelic variation in HLA-DQB1 promoter elements

Sequence variation among HLA class II promoter elements may contribute to functional differences in transcriptional regulation of different class II alleles. In addition to influencing the binding sites for nuclear transcription factors, promoter polymorphism may also alter intrinsic structural properties of the DNA strands, such as conformation and curvature, which influence the formation of stable transcription complexes. We used SSCP analysis of PCR-amplified promoter regions from the DQB1 locus to evaluate conformational polymorphism within DQ alleles. Distinct electrophoretic migration patterns of the SSCP products were detected for six DQB1 alleles; analysis of the DQB1*0302 promoter, known to be associated with type 1 diabetes, showed no SSCP differences between IDDM patients and normal controls. Using computer modeling based on a "nearest-neighbor" energy of predicted curvature theory, we examined the effect of allelic promoter region sequence polymorphism on the predicted curvature of double-stranded DNA, and found distinct allelic differences in predicted DNA curvature, both in transcriptional consensus binding sites and in regions located between binding sites. These data are consistent with a model in which intrinsic sequence variation in the promoter region results in ultrastructural differences which may influence DNA bending and interactions with multimeric DNA-protein transcription complexes.

Transcriptional regulation of MHC class II genes

MHC class II molecules play a fundamental role in the homeostasis of the immune response, functioning as receptors for antigenic peptides to be presented to regulatory T cells. Both quantitative and qualitative alterations in the expression of these molecules on the cell surface dramatically affect the onset of the immune response, and may be the basis of a wide variety of disease states, such as autoimmunity, immunodeficiencies, and cancer. Most regulation of MHC class II molecule expression is under the control of transcription mechanisms which are both cell type and development specific. In the last few years classical genetics together with molecular biology have greatly contributed to the widening of our knowledge on the regulatory mechanisms operating in the control of class II gene expression. This review deals with the latest developments in this fundamental area of immunology.

Differential expression of HLA-DQB1 alleles in a heterozygous cell line

Regulation of HLA class II transcription is complex, with both locus-specific and allele-specific polymorphisms associated with consensus regulatory region promoter elements. In order to evaluate the potential function of allele-specific elements, HLA-DQB1 transcripts and in vivo footprinting of the DQB1 promoter were studied in human cell lines inducible for HLA class II expression. In the heterozygous melanoma cell line Me9229/18, differential expression was observed for two DQB1 alleles. This variation was reflected in both steady-state and inducible transcripts, although methylation patterns of the promoter elements were similar for both alleles. This in vivo allelic difference expression, which correlates with previous studies of in vitro reporter gene transcription, indicates the potential for functionally important differences in tissue-specific and inducible expression attributable to promoter region polymorphism in HLA genes.

Polymorphism in the upstream regulatory region of DQA1 gene in the Italian population

Polymorphism in the 5'-upstream regulatory region of the DQA1 gene has been recently described. Using PCR-SSO method and SSCP analysis we have investigated this polymorphism in a group of 111 Italian blood donors which had been oligotyped for DRB1, DQA1 and DQB1 genes. Eight allelic variants were detected. Looking at the relationships among QAP sequences and DQA1 and DRB1 genes, three alternative situations were found: 1. a one-to-one relation between QAP and DQA1 alleles, independently of the other class II genes; 2. the same QAP allele in association with different DQA1-DRB1 haplotypes; 3. the same DQA1 allele with different QAP sequences according to the DRB1 specificity. No unexpected associations with DQB1 gene were found. These results must be interpreted considering that DQA1 and DRB1 genes are transcribed in opposite directions so that the promoter region of DQA1 gene lies between DQA1 and DRB1, close to the former but several hundreds kb away from the latter.