Extended HLA-DPB1 polymorphism: an RNA approach for HLA-DPB1 typing

HLA class II peptide-binding-region analysis reveals funneling of polymorphism in action

BACKGROUND

HLA-class II proteins hold important roles in key physiological processes. The purpose of this study was to compile all class II alleles reported in human population and investigate patterns in pocket variants and their combinations, focusing on the peptide-binding region (PBR).

METHODS

For this purpose, all protein sequences of DPA1, DQA1, DPB1, DQB1 and DRB1 were selected and filtered, in order to have full PBR sequences. Proportional representation was used for pocket variants while population data were also used.

RESULTS

All pocket variants and PBR sequences were retrieved and analyzed based on the preference of amino acids and their properties in all pocket positions. The observed number of pocket variants combinations was much lower than the possible inferred, suggesting that PBR formation is under strict funneling. Also, although class II proteins are very polymorphic, in the majority of the reported alleles in all populations, a significantly less polymorphic pocket core was found.

CONCLUSIONS

Pocket variability of five HLA class II proteins was studied revealing favorable properties of each protein. The actual PBR sequences of HLA class II proteins appear to be governed by restrictions that lead to the establishment of only a fraction of the possible combinations and the polymorphism recorded is the result of intense funneling based on function.

Genome-wide minor histocompatibility matching as related to the risk of graft-versus-host disease

The risk of acute graft-versus-host disease (GVHD) is higher after allogeneic hematopoietic cell transplantation (HCT) from unrelated donors as compared with related donors. This difference has been explained by increased recipient mismatching for major histocompatibility antigens or minor histocompatibility antigens. In the current study, we used genome-wide arrays to enumerate single nucleotide polymorphisms (SNPs) that produce graft-versus-host (GVH) amino acid coding differences between recipients and donors. We then tested the hypothesis that higher degrees of genome-wide recipient GVH mismatching correlate with higher risks of GVHD after allogeneic HCT. In HLA-genotypically matched sibling recipients, the average recipient mismatching of coding SNPs was 9.35%. Each 1% increase in genome-wide recipient mismatching was associated with an estimated 20% increase in the hazard of grades III-IV GVHD (hazard ratio [HR], 1.20; 95% confidence interval [CI], 1.05-1.37; P = .007) and an estimated 22% increase in the hazard of stage 2-4 acute gut GVHD (HR, 1.22; 95% CI, 1.02-1.45; P = .03). In HLA-A, B, C, DRB1, DQA1, DQB1, DPA1, DPB1-phenotypically matched unrelated recipients, the average recipient mismatching of coding SNPs was 17.3%. The estimated risks of GVHD-related outcomes in HLA-phenotypically matched unrelated recipients were low, relative to the large difference in genome-wide mismatching between the 2 groups. In contrast, the risks of GVHD-related outcomes were higher in HLA-DP GVH-mismatched unrelated recipients than in HLA-matched sibling recipients. Taken together, these results suggest that the increased GVHD risk after unrelated HCT is predominantly an effect of HLA-mismatching.

The influence of positive selection and trans-species evolution on DPB diversity in the golden snub-nosed monkeys (Rhinopithecus roxellana)

Genetic variation plays a significant role in the adaptive potential of the endangered species. The variation at major histocompatibility complex (MHC) genes can offer valuable information on selective pressure related to natural selection and environmental adaptation, particularly the ability of a host to continuously resist evolving parasites. Thus, the genetic polymorphism on exon 2 of the MHC DPB1 gene in the golden snub-nosed monkeys (Rhinopithecus roxellana) was specifically analyzed. The results show that the 6 Rhro-DPB1 alleles identified from 87 individuals exhibit positive selection and trans-species polymorphism. The results also imply that although the populations of the species have experienced dramatic reduction and severe habitat fragmentation in recent Chinese history, balancing selection still maintains relatively consistent, with moderate DPB1 polymorphism. Thus, the study provides valuable information and evidence in developing effective strategies and tactics for genetic health and population size expansion of the species. It also offers strong genetic background for further studies on other primate species, particularly those in Rhinopithecus-a further endeavor that would result in fully understanding the MHC genetic information of the Asian colobines.

Effects of transmembrane region variability on cell surface expression and allorecognition of HLA-DP3

The functional relevance of polymorphisms outside the peptide binding groove of HLA molecules is poorly understood. Here we have addressed this issue by studying HLA-DP3, a common antigen relevant for functional matching algorithms of unrelated hematopoietic stem cell transplantation (HSCT) encoded by two transmembrane (TM) region variants, DPB1(*)03:01 and DPB1(*)104:01. The two HLA-DP3 variants were found at a overall allelic frequency of 10.4% in 201 volunteer stem cell donors, at a ratio of 4.2:1. No significant differences were observed in cell surface expression levels of the two variants on B lymphoblastoid cell lines (BLCL), primary B cells or monocytes. Three different alloreactive T cell lines or clones showed similar levels of activation marker CD107a and/or CD137 upregulation in response to HLA-DP3 encoded by DPB1(*)03:01 and DPB1(*)104:01, either endogenously on BLCL or after lentiveral-vector mediated transfer into the same cellular background. These data provide, for the first time, direct evidence for a limited functional role of a TM region polymorphism on expression and allorecognition of HLA-DP3 and are compatible with the notion that the two variants can be considered as a single functional entity for unrelated stem cell donor selection.

Molecular characterization, balancing selection, and genomic organization of the tree shrew (Tupaia belangeri) MHC class I gene

The major histocompatibility complex (MHC) class I genes play a pivotal role in the adaptive immune response among vertebrates. Accordingly, in numerous mammals the genomic structure and molecular characterization of MHC class I genes have been thoroughly investigated. To date, however, little is known about these genes in tree shrews, despite the increasingly popularity of its usage as an animal model. To address this deficiency, we analyzed the structure and characteristic of the tree shrew MHC class I genes (Tube-MHC I) and performed a comparative gene analysis of the tree shrew and other mammal species. We found that the full-length cDNA sequence of the tree shrew MHC class I is 1074bp in length. The deduced peptide is composed of 357 amino acids containing a leader peptide, an α1 and α2 domain, an α3 domain, a transmembrane domain and a cytoplasmic domain. Among these peptides, the cysteines, CD8(+) interaction and N-glycosylation sites are all well conserved. Furthermore, the genomic sequence of the tree shrew MHC class I gene was identified to be 3180bp in length, containing 8 exons and 7 introns. In 21 MHC class I sequences, we conducted an extensive study of nucleotide substitutions. The results indicated that in the peptide binding region (PBR) the rate of non-synonymous substitutions (dN) to synonymous substitutions (dS) was greater than 1, suggesting balancing selection at the PBR. These findings provide valuable contributions in furthering our understanding of the structure, molecular polymorphism, and function of the MHC class I genes in tree shrews, further improving their utility as an animal model in biomedical research.

Allelic diversity of the MHC class II DRB genes in brown bears (Ursus arctos) and a comparison of DRB sequences within the family Ursidae

The allelic diversity of the DRB locus in major histocompatibility complex (MHC) genes was analyzed in the brown bear (Ursus arctos) from the Hokkaido Island of Japan, Siberia, and Kodiak of Alaska. Nineteen alleles of the DRB exon 2 were identified from a total of 38 individuals of U. arctos and were highly polymorphic. Comparisons of non-synonymous and synonymous substitutions in the antigen-binding sites of deduced amino acid sequences indicated evidence for balancing selection on the bear DRB locus. The phylogenetic analysis of the DRB alleles among three genera (Ursus, Tremarctos, and Ailuropoda) in the family Ursidae revealed that DRB allelic lineages were not separated according to species. This strongly shows trans-species persistence of DRB alleles within the Ursidae.

Strong association of the HLA-DP6 supertype with childhood leukaemia is due to a single allele, DPB1*0601

We previously reported that susceptibility to childhood B cell precursor ALL (BCP ALL) is associated with HLA-DPB1 alleles having glutamic acid (E) rather than lysine (K) in the P4 antigenic peptide-binding pocket. Clustering approximately 90% of DPB1 alleles into DPB69E (DP2, 6, 8) and DPB69K (DP1, 3, 4) supertypes revealed that DP2 and DP8 are associated with BCP ALL, but DP6 is also associated with non-BCP leukaemia. Here, we report that only one of seven alleles with the DP6 supertype (DPB1(*)0601) is associated with childhood leukaemia (leukaemia vs controls: odds ratio, 95% confidence interval [OR, CI]: 4.6, 2.0-10.4; corrected P=0.019), but not with childhood solid tumours or lymphomas. DPB1(*)0601 is also significantly associated with leukaemia subtypes, including BCP ALL, Pro-B ALL, T-ALL and AML. DPB1(*)0601 is significantly over-transmitted (76.9%) from parents to children with BCP ALL (OR; CI: 4.7; 1.01-22.2). Sequencing the coding region of DPB1(*)0601 revealed an exon 1-4 haplotype [T-DEAV-KIL-RVI] shared with DPB1(*)0301 and 0901, but no evidence of germline mutations in childhood leukaemia. These results suggest that the DPbeta0601 molecule may be functionally involved in childhood leukaemia. Analysis of peptide binding and T-cell activation by DPbeta0601-peptide complexes should help determine its role in childhood leukaemia causation.

Pinpointing a selective sweep to the chimpanzee MHC class I region by comparative genomics

Chimpanzees experienced a reduction of the allelic repertoire at the major histocompatibility complex (MHC) class I A and B loci, which may have been caused by a retrovirus belonging to the simian immunodeficiency virus (SIV) family. Extended MHC haplotypes were defined in a pedigreed chimpanzee colony. Comparison of genetic variation at microsatellite markers mapping inside and outside the Mhc region was carried out in humans and chimpanzees to investigate the genomic extent of the repertoire reduction. Multilocus demographic analyses underscored that chimpanzees indeed experienced a selective sweep that mainly targeted the chromosomal segment carrying the Mhc class I region. Probably due to genetic linkage, the sweep also affected other polymorphic loci, mapping in the close vicinity of the Mhc class I region genes. Nevertheless, although the allelic repertoire at particular Mhc class I and II loci appears to be limited, naturally occurring recombination events allowed the establishment of haplotype diversity after the sweep. However, recombination did not have sufficient time to erase the signal of the selective sweep.

Frequency and targeted detection of HLA-DPB1 T cell epitope disparities relevant in unrelated hematopoietic stem cell transplantation

The majority of unrelated donor (UD) hematopoietic stem cell (HSC) transplants are performed across HLA-DP mismatches, which, if involving disparity in a host-versus-graft (HVG) direction for an alloreactive T cell epitope (TCE), have been shown by our group to be associated with poor clinical outcome in 2 cohorts of patients transplanted for hematopoietic malignancies and beta-thalassemia, respectively. Using site-directed mutagenesis of DPB1*0901, we show here that the TCE is abrogated by the presence of amino acids LFQG in positions 8-11 of the DP beta-chain. Based on this and on alloreactive T cell responsiveness, we have determined the presence or absence of the TCE for 72 DPB1 alleles reported in the ethnic groups representative of the worldwide UD registries, and predict that 67%-87% (mean 77%) of UD recipient pairs will not present a DPB1 TCE disparity in the HVG direction. We developed and validated in 112 healthy Italian blood donors an innovative approach of DPB1 epitope-specific typing (EST), based on 2 PCR reactions. Our data show that DPB1 TCE disparities may hamper the clinical success of a considerable number of transplants when DPB1 matching is not included into the donor selection criteria, and that a donor without DPB1 TCE disparities in the HVG direction can be found for the majority of patients. Moreover, the study describes the first protocol of targeted epitope-specific DPB1 donor-recipient matching for unrelated HSC transplantation. This protocol will facilitate large-scale retrospective clinical studies warranted to more precisely determine the clinical relevance of DPB1 TCE disparities in different transplant conditions.

Sequence-based typing of HLA-DQA1: comprehensive approach showed molecular heterogeneity

Within the human leukocyte antigen-DQA1 workshop project the level of molecular heterogeneity of the DQA1 gene was investigated. An improved sequence-based typing protocol was used, enabling analysis of the complete coding sequence, comprising exons 1-4. The participating laboratories implemented the amplification and sequencing primers in their own sequence-based typing approach. The method proved to be sufficiently robust to handle the differences in protocols. All reference samples used for validation were correctly typed for DQA1 by all participating laboratories. Three different populations with a total of 736 individuals were investigated: a population of Korean origin (n= 467), a British Caucasian (n= 114), and a Dutch Caucasian (n= 155) population. Sixteen of the known 28 DQA1 alleles were detected and six new alleles were identified. All novel alleles showed a nucleotide substitution outside exon 2. Comparison of the calculated allele frequencies revealed major differences between the Korean and the Caucasian populations but also between Dutch and British Caucasians. A tight association between DQA1 and DRB1/DQB1 alleles was observed in all three populations.

Comparative genetics of MHC polymorphisms in different primate species: duplications and deletions

Gene products of the major histocompatibility complex (MHC) play a crucial role in the activation of adaptive (antigen-dependent) immune responses. In this paper similarities and dissimilarities among the MHCs of different primate species and their functional implications are reviewed. The human HLA system represents the most thoroughly investigated MHC of any contemporary living primate species, and so it will serve as a reference.