5' and 3' untranslated regions contribute to the differential expression of specific HLA-A alleles

Nucleotide alterations in the HLA-C class I gene can cause aberrant splicing and marked changes in RNA levels in a polymorphic context-dependent manner

Polymorphisms of HLA genes, which play a crucial role in presenting peptides with diverse sequences in their peptide-binding pockets, are also thought to affect HLA gene expression, as many studies have reported associations between HLA gene polymorphisms and their expression levels. In this study, we devised an ectopic expression assay for the HLA class I genes in the context of the entire gene, and used the assay to show that the HLA-C*03:03:01 and C*04:01:01 polymorphic differences observed in association studies indeed cause different levels of RNA expression. Subsequently, we investigated the C*03:23N null allele, which was previously noted for its reduced expression, attributed to an alternate exon 3 3' splice site generated by G/A polymorphism at position 781 within the exon 3. We conducted a thorough analysis of the splicing patterns of C*03:23N, and revealed multiple aberrant splicing, including the exon 3 alternative splicing, which overshadowed its canonical counterpart. After confirming a significant reduction in RNA levels caused by the G781A alteration in our ectopic assay, we probed the function of the G-rich sequence preceding the canonical exon 3 3' splice site. Substituting the G-rich sequence with a typical pyrimidine-rich 3' splice site sequence on C*03:23N resulted in a marked elevation in RNA levels, likely due to the enhanced preference for the canonical exon 3 3' splice site over the alternate site. However, the same substitution led to a reduction in RNA levels for C*03:03:01. These findings suggested the dual roles of the G-rich sequence in RNA expression, and furthermore, underscore the importance of studying polymorphism effects within the framework of the entire gene, extending beyond conventional mini-gene reporter assays.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

Regulation of HLA class I expression by non-coding gene variations

The Human Leukocyte Antigen (HLA) is a critical genetic system for different outcomes after solid organ and hematopoietic cell transplantation. Its polymorphism is usually determined by molecular technologies at the DNA level. A potential role of HLA allelic expression remains under investigation in the context of the allogenic immune response between donors and recipients. In this study, we quantified the allelic expression of all three HLA class I loci (HLA-A, B and C) by RNA sequencing and conducted an analysis of expression quantitative traits loci (eQTL) to investigate whether HLA expression regulation could be associated with non-coding gene variations. HLA-B alleles exhibited the highest expression levels followed by HLA-C and HLA-A alleles. The max fold expression variation was observed for HLA-C alleles. The expression of HLA class I loci of distinct individuals demonstrated a coordinated and paired expression of both alleles of the same locus. Expression of conserved HLA-A~B~C haplotypes differed in distinct PBMC's suggesting an individual regulated expression of both HLA class I alleles and haplotypes. Cytokines TNFα /IFNβ, which induced a very similar upregulation of HLA class I RNA and cell surface expression across alleles did not modify the individually coordinated expression at the three HLA class I loci. By identifying cis eQTLs for the HLA class I genes, we show that the non-coding eQTLs explain 29%, 13%, and 31% of the respective HLA-A, B, C expression variance in unstimulated cells, and 9%, 23%, and 50% of the variance in cytokine-stimulated cells. The eQTLs have significantly higher effect sizes in stimulated cells compared to unstimulated cells for HLA-B and HLA-C genes expression. Our data also suggest that the identified eQTLs are independent from the coding variation which defines HLA alleles and thus may be influential on intra-allele expression variability although they might not represent the causal eQTLs.

HLA-A alleles influencing NK cell function impact AML relapse following allogeneic hematopoietic cell transplantation

HLA-B allotypes exhibiting the Bw4 epitope trigger variable inhibitory signaling of KIR3DL1 receptor types, where strong inhibitory HLA-B and KIR3DL1 allele combinations are associated with increased risk for relapse of acute myelogenous leukemia (AML) following allogeneic hematopoietic cell transplantation (HCT). Several HLA-A allotypes also exhibit the Bw4 epitope. Studies with natural killer (NK) cell clones have demonstrated NK inhibition via KIR3DL1 by HLA-A Bw4+ allotypes, but did not delineate strengths of inhibition or hierarchies of NK education. Using primary NK cells from healthy donors, we demonstrate that HLA-A*23, HLA-A*24, and HLA-A*32 proteins are expressed at different densities and exhibit different capacities to educate and inhibit KIR3DL1-expressing NK cells in vitro. Among the HLA-A Bw4+ allotypes, HLA-A*24 and HLA-A*32 demonstrate the strongest inhibitory capacity. To determine if HLA-A allotypes with strong inhibitory capacity have similar negative impact in allogeneic HCT as HLA-B Bw4+ allotypes, we performed a retrospective analysis of 1729 patients with AML who received an allogeneic HCT from a 9/10 or 10/10 HLA allele-matched unrelated donor. Examination of the donor-recipient pairs whose Bw4 epitope was exclusively contributed from HLA-A*24 and A*32 allotypes revealed that patients with HLA-A*24 who received an allograft from a KIR3DL1+ donor experienced a higher risk of disease relapse (hazard ratio, 1.65; 95% confidence interval, 1.17-2.32; P = .004) when compared with patients without a Bw4 epitope. These findings indicate that despite weak affinity interactions with KIR3DL1, common HLA-A allotypes with the Bw4 epitope can interact with KIR3DL1+ donor NK cells with clinically meaningful impact and provide additional insight to donor NK alloreactivity in HLA-matched HCT.

Targeting the coding sequence: opposing roles in regulating classical and non-classical MHC class I molecules by miR-16 and miR-744

Background

To control gene expression, microRNAs (miRNAs) are of key importance and their deregulation is associated with the development and progression of various cancer types. In this context, a discordant messenger RNA/protein expression pointing to extensive post-transcriptional regulation of major histocompatibility complex (MHC) class I molecules was already shown. However, only a very limited number of miRNAs targeting these molecules have yet been identified. Despite an increasing evidence of coding sequence (CDS)-located miRNA binding sites, there exists so far, no detailed study of the interaction of miRNAs with the CDS of MHC class I molecules.

Methods

Using an MS2-tethering approach in combination with small RNA sequencing, a number of putative miRNAs binding to the CDS of human leukocyte antigen (HLA)-G were identified. These candidate miRNAs were extensively screened for their effects in the HLA-G-positive JEG3 cell line. Due to the high sequence similarity between HLA-G and classical MHC class I molecules, the impact of HLA-G candidate miRNAs on HLA class I surface expression was also analyzed. The Cancer Genome Atlas data were used to correlate candidate miRNAs and HLA class I gene expression.

Results

Transfection of candidate miRNAs revealed that miR-744 significantly downregulates HLA-G protein levels. In contrast, overexpression of the candidate miRNAs miR-15, miR-16, and miR-424 sharing the same seed sequence resulted in an unexpected upregulation of HLA-G. Comparable results were obtained for classical MHC class I members after transfection of miRNA mimics into HEK293T cells. Analyses of The Cancer Genome Atlas data sets for miRNA and MHC class I expression further validated the results.

Conclusions

Our data expand the knowledge about MHC class I regulation and showed for the first time an miRNA-dependent control of MHC class I antigens mediated by the CDS. CDS-located miRNA binding sites could improve the general use of miRNA-based therapeutic approaches as these sites are highly independent of structural variations (e.g. mutations) in the gene body. Surprisingly, miR-16 family members promoted MHC class I expression potentially in a gene activation-like mechanism.

Reduced PCR-generated errors from a hybrid capture-based NGS assay for HLA typing

Next generation sequencing (NGS) assays are state of the art for HLA genotyping. To sequence on an Illumina sequencer, the DNA of interest must be enriched, fragmented, and bookended with known oligonucleotide sequences, a process known as library construction. Many HLA genotyping assays enrich the target loci by long-range PCR (LR-PCR), prior to fragmentation. This PCR step has been reported to introduce errors in the DNA to be sequenced, including inaccurate replication of repeated sequences, and the in vitro recombination of alleles encoded on separate chromosomes. An alternative library construction method involves fragmentation of genomic DNA, followed by hybrid-capture (HC) enrichment of target HLA loci. This HC-based method involves PCR, but with far fewer cycles. Consequently, the HC method had significantly fewer PCR-induced errors, including more faithful replication of repeated sequences, and the near elimination of recombinant sequences. These improvements likely produce more accurate NGS sequencing data of HLA loci.

ERAP2 Increases the Abundance of a Peptide Submotif Highly Selective for the Birdshot Uveitis-Associated HLA-A29

Birdshot Uveitis (BU) is a blinding inflammatory eye condition that only affects HLA-A29-positive individuals. Genetic association studies linked ERAP2 with BU, an aminopeptidase which trims peptides before their presentation by HLA class I at the cell surface, which suggests that ERAP2-dependent peptide presentation by HLA-A29 drives the pathogenesis of BU. However, it remains poorly understood whether the effects of ERAP2 on the HLA-A29 peptidome are distinct from its effect on other HLA allotypes. To address this, we focused on the effects of ERAP2 on the immunopeptidome in patient-derived antigen presenting cells. Using complementary HLA-A29-based and pan-class I immunopurifications, isotope-labeled naturally processed and presented HLA-bound peptides were sequenced by mass spectrometry. We show that the effects of ERAP2 on the N-terminus of ligands of HLA-A29 are shared across endogenous HLA allotypes, but discover and replicate that one peptide motif generated in the presence of ERAP2 is specifically bound by HLA-A29. This motif can be found in the amino acid sequence of putative autoantigens. We further show evidence for internal sequence specificity for ERAP2 imprinted in the immunopeptidome. These results reveal that ERAP2 can generate an HLA-A29-specific antigen repertoire, which supports that antigen presentation is a key disease pathway in BU.

HLA RNA Sequencing With Unique Molecular Identifiers Reveals High Allele-Specific Variability in mRNA Expression

The HLA gene complex is the most important single genetic factor in susceptibility to most diseases with autoimmune or autoinflammatory origin and in transplantation matching. Most studies have focused on the vast allelic variation in these genes; only a few studies have explored differences in the expression levels of HLA alleles. In this study, we quantified mRNA expression levels of HLA class I and II genes from peripheral blood samples of 50 healthy individuals. The gene- and allele-specific mRNA expression was assessed using unique molecular identifiers, which enabled PCR bias removal and calculation of the number of original mRNA transcripts. We identified differences in mRNA expression between different HLA genes and alleles. Our results suggest that HLA alleles are differentially expressed and these differences in expression levels are quantifiable using RNA sequencing technology. Our method provides novel insights into HLA research, and it can be applied to quantify expression differences of HLA alleles in various tissues and to evaluate the role of this type of variation in transplantation matching and susceptibility to autoimmune diseases.

HLA-A29 and Birdshot Uveitis: Further Down the Rabbit Hole

HLA class I alleles constitute established risk factors for non-infectious uveitis and preemptive genotyping of HLA class I alleles is standard practice in the diagnostic work-up. The HLA-A29 serotype is indispensable to Birdshot Uveitis (BU) and renders this enigmatic eye condition a unique model to better understand how the antigen processing and presentation machinery contributes to non-infectious uveitis or chronic inflammatory conditions in general. This review will discuss salient points regarding the protein structure of HLA-A29 and how key amino acid positions impact the peptide binding preference and interaction with T cells. We discuss to what extent the risk genes ERAP1 and ERAP2 uniquely affect HLA-A29 and how the discovery of a HLA-A29-specific submotif may impact autoantigen discovery. We further provide a compelling argument to solve the long-standing question why BU only affects HLA-A29-positive individuals from Western-European ancestry by exploiting data from the 1000 Genomes Project. We combine novel insights from structural and immunopeptidomic studies and discuss the functional implications of genetic associations across the HLA class I antigen presentation pathway to refine the etiological basis of Birdshot Uveitis.

Clinical relevance of donor-specific antibodies directed at HLA-C: A long road to acceptance

In solid organ transplantation (SOT), the clinical relevance of donor-specific antibodies (DSA) directed at anti-HLA-A, -B, -DR and -DQ antigens is largely recognized while it is still a matter of debate for DSA directed at HLA-C. In this review, we summarize the peculiarities of HLA-C among class I HLA antigens as well as their immunogenicity, which underlie the clinical relevance of HLA-C locus and anti-HLA-C DSA in SOT. Many factors, both intrinsic and extrinsic to the HLA-C gene and HLA-C protein, explain its lower expression in comparison with HLA-A and -B. This lower expression can explain the apparent lower immunogenicity of HLA-C leading to a lower prevalence and strength of anti-HLA-C antibodies. Nevertheless, HLA-C antigens are truly immunogenic and preformed anti-HLA-C DSA are clinically relevant. Indeed, anti-HLA-C DSA are able to bind donor cells and to activate the complement pathway both ex vivo and in vivo. In line with this, numerous clinical studies now show that preformed DSA directed at native HLA-C molecules induce poorer graft outcomes. We then plead for the inclusion of HLA-C in all transplant allocation systems and we propose a strategy to cope with anti-HLA-C DSA in SOT. Beyond SOT, anti-HLA-C antibodies generate a growing interest in the allo-HCT, transfusion and obstetrics fields, while new concepts such as the role of the "missing-self" in solid organ rejection places HLA-C as an inescapable actor in transplant tolerance.

Review of Genetic Variation as a Predictive Biomarker for Chronic Graft-Versus-Host-Disease After Allogeneic Stem Cell Transplantation

Chronic graft-versus-host disease (cGvHD) is one of the major complications of allogeneic stem cell transplantation (HSCT). cGvHD is an autoimmune-like disorder affecting multiple organs and involves a dermatological rash, tissue inflammation and fibrosis. The incidence of cGvHD has been reported to be as high as 30% to 60% and there are currently no reliable tools for predicting the occurrence of cGvHD. There is therefore an important unmet clinical need for predictive biomarkers. The present review summarizes the state of the art for genetic variation as a predictive biomarker for cGvHD. We discuss three different modes of action for genetic variation in transplantation: genetic associations, genetic matching, and pharmacogenetics. The results indicate that currently, there are no genetic polymorphisms or genetic tools that can be reliably used as validated biomarkers for predicting cGvHD. A number of recommendations for future studies can be drawn. The majority of studies to date have been under-powered and included too few patients and genetic markers. Like in all complex multifactorial diseases, large collaborative genome-level studies are now needed to achieve reliable and unbiased results. Some of the candidate genes, in particular, CTLA4, HSPE, IL1R1, CCR6, FGFR1OP, and IL10, and some non-HLA variants in the HLA gene region have been replicated to be associated with cGvHD risk in independent studies. These associations should now be confirmed in large well-characterized cohorts with fine mapping. Some patients develop cGvHD despite very extensive immunosuppression and other treatments, indicating that the current therapeutic regimens may not always be effective enough. Hence, more studies on pharmacogenetics are also required. Moreover, all of these studies should be adjusted for diagnostic and clinical features of cGvHD. We conclude that future studies should focus on modern genome-level tools, such as machine learning, polygenic risk scores and genome-wide association study-transcription meta-analyses, instead of focusing on just single variants. The risk of cGvHD may be related to the summary level of immunogenetic differences, or whole genome histocompatibility between each donor-recipient pair. As the number of genome-wide analyses in HSCT is increasing, we are approaching an era where there will be sufficient data to incorporate these approaches in the near future.

HLA-G whole gene amplification reveals linkage disequilibrium between the HLA-G 3'UTR and coding sequence

Polymorphic sites in the HLA-G gene may influence expression and function of the protein. Knowledge of the association between high-resolution HLA-G alleles and 3-prime untranslated (3'UTR) haplotypes is useful for studies on the role of HLA-G in transplantation, pregnancy, and cancer. We developed a next generation sequencing (NGS)-based typing assay enabling full phasing over the whole HLA-G gene sequence with inclusion of the 3'UTR region. DNA from 171 mother-child pairs (342 samples) was studied for: (a) HLA-G allele information by the NGSgo-AmpX HLA-G assay, (b) 3'UTR haplotype information by an in-house developed sequence-based typing method of a 699/713 base pair region in the 3'UTR, and (c) the full phase HLA-G gene sequence, by combining primers from both assays. The mother to child inheritance allowed internal verification of newly identified alleles and of association between coding and UTR regions. The NGSgo workflow compatible with Illumina platforms was employed. Data was interpreted using NGSengine software. In 99.4% of all alleles analyzed, the extended typing was consistent with the separate allele and 3'UTR typing methods. After repeated analysis of four samples that showed discrepancy, consistency reached 100%. A high-linkage disequilibrium between IPD-IMGT/HLA Database-defined HLA-G alleles and the extended 3'UTR region was identified (D' = 0.994, P < .0001). Strong associations were found particularly between HLA-G*01:04 and UTR-3, between HLA-G*01:01:03 and UTR-7, and between HLA-G*01:03:01 and UTR-5 (for all: r = 1). Six novel HLA-G alleles and three novel 3'UTR haplotype variants were identified, of which three and one, respectively, were verified in the offspring.

Dissecting Genetic Control of HLA-DPB1 Expression and Its Relation to Structural Mismatch Models in Hematopoietic Stem Cell Transplantation

HLA expression levels have been suggested to be genetically controlled by single nucleotide polymorphisms (SNP) in the untranslated regions (UTR), and expression variants have been associated with the outcome of chronic viral infection and hematopoietic stem cell transplantation (HSCT). In particular, the 3′UTR rs9277534-G/A SNP in HLA-DPB1 has been associated with graft-versus-host-disease after HSCT (Expression model); however its relevance in different immune cells and its mode of action have not been systematically addressed. In addition, there is a strong though not complete overlap between the rs9277534-G/A SNP and structural HLA-DPB1 T cell epitope (TCE) groups which have also been associated with HSCT outcome (TCE Structural model). Here we confirm and extend previous findings of significantly higher HLA-DPB1 expression in B cell lines, unstimulated primary B cells, and monocytes homozygous for rs9277534-G compared to those homozygous for rs9277534-A. However, these differences were abrogated by interferon-γ stimulation or differentiation into dendritic cells. We identify at least seven 3′UTR rs9277534-G/A haplotypes differing by a total of 37 SNP, also characterized by linkage to length variants of a short tandem repeat (STR) in intron 2 and TCE group assignment. 3′UTR mapping did not show any significant differences in post-transcriptional regulation assessed by luciferase assays between two representative rs9277534-G/A haplotypes for any of eight overlapping fragments. Moreover, no evidence for alternative splicing associated with the intron 2 STR was obtained by RT-PCR. In an exemplary cohort of 379 HLA-DPB1 mismatched donor-recipient pairs, risk prediction by the Expression model and the Structural TCE model was 36.7% concordant, with the majority of discordances due to non-applicability of the Expression model. HLA-DPB1 from different TCE groups expressed in the absence of the 3′UTR at similar levels by transfected HeLa cells elicited significantly different mean alloreactive CD4+ T-cell responses, as assessed by CD137 upregulation assays in 178 independent cultures. Taken together, our data provide new insights into the cell type-specific and mechanistic basis of the association between the rs9277534-G/A SNP and HLA-DPB1 expression, and show that, despite partial overlap between both models in HSCT risk-prediction, differential alloreactivity determined by the TCE structural model occurs independently from HLA-DPB1 differential expression.

Platelets from donors with consistently low HLA-B8, -B12, or -B35 expression do not undergo antibody-mediated internalization

HLA-B8, -B12, or -B35 expression on platelets varies significantly between donors and in certain donors is consistently low or undetectable. Antibody-mediated internalization of platelets correlates with antigen expression and is absent in platelets with low HLA expression.

Quantification of classical HLA class I mRNA by allele-specific, real-time polymerase chain reaction for most Han individuals

Recent studies have shown that expression levels of different alleles at the same HLA class I locus can vary dramatically, which might have a broad influence on human disease. However, precise quantification of the relative expression level of each HLA allele is challenging, because distinguishing different alleles on the same locus is difficult. Here, we developed a series of allele-specific, real-time polymerase chain reaction assays for quantifying HLA class I allele mRNA in most Han individuals. The alleles of almost all heterozygous genotypes with a frequency higher than 0.5% in our population (78 alleles on HLA-A locus, 124 alleles on HLA-B locus, and 74 alleles on HLA-C locus) were specifically amplified. The specificity of the amplification was strictly validated by setting the corresponding negative control for each allele of each genotype. The amplification efficiency of each reaction was determined, and the slopes of the reactions were compared. This study provides a tool for detecting the comprehensive expression profile of HLA class I alleles and will be useful not only for the investigation of the molecular mechanism underlying HLA allele expression regulation but also for exploration of immunological mechanisms involving HLA expression in the fields of tumour immune evasion, viral infection, auto-immune disorders, and graft vs host disease after haematopoietic stem cell transplantation.

HLA genotyping by next-generation sequencing of complementary DNA

Background

Genotyping of the human leucocyte antigen (HLA) is indispensable for various medical treatments. However, unambiguous genotyping is technically challenging due to high polymorphism of the corresponding genomic region. Next-generation sequencing is changing the landscape of genotyping. In addition to high throughput of data, its additional advantage is that DNA templates are derived from single molecules, which is a strong merit for the phasing problem. Although most currently developed technologies use genomic DNA, use of cDNA could enable genotyping with reduced costs in data production and analysis. We thus developed an HLA genotyping system based on next-generation sequencing of cDNA.

Methods

Each HLA gene was divided into 3 or 4 target regions subjected to PCR amplification and subsequent sequencing with Ion Torrent PGM. The sequence data were then subjected to an automated analysis. The principle of the analysis was to construct candidate sequences generated from all possible combinations of variable bases and arrange them in decreasing order of the number of reads. Upon collecting candidate sequences from all target regions, 2 haplotypes were usually assigned. Cases not assigned 2 haplotypes were forwarded to 4 additional processes: selection of candidate sequences applying more stringent criteria, removal of artificial haplotypes, selection of candidate sequences with a relaxed threshold for sequence matching, and countermeasure for incomplete sequences in the HLA database.

Results

The genotyping system was evaluated using 30 samples; the overall accuracy was 97.0% at the field 3 level and 98.3% at the G group level. With one sample, genotyping of DPB1 was not completed due to short read size. We then developed a method for complete sequencing of individual molecules of the DPB1 gene, using the molecular barcode technology.

Conclusion

The performance of the automatic genotyping system was comparable to that of systems developed in previous studies. Thus, next-generation sequencing of cDNA is a viable option for HLA genotyping.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4300-7) contains supplementary material, which is available to authorized users.

Killer Immunoglobulin-Like Receptor Allele Determination Using Next-Generation Sequencing Technology

The impact of natural killer (NK) cell alloreactivity on hematopoietic stem cell transplantation (HSCT) outcome is still debated due to the complexity of graft parameters, HLA class I environment, the nature of killer cell immunoglobulin-like receptor (KIR)/KIR ligand genetic combinations studied, and KIR⁺ NK cell repertoire size. KIR genes are known to be polymorphic in terms of gene content, copy number variation, and number of alleles. These allelic polymorphisms may impact both the phenotype and function of KIR⁺ NK cells. We, therefore, speculate that polymorphisms may alter donor KIR⁺ NK cell phenotype/function thus modulating post-HSCT KIR⁺ NK cell alloreactivity. To investigate KIR allele polymorphisms of all KIR genes, we developed a next-generation sequencing (NGS) technology on a MiSeq platform. To ensure the reliability and specificity of our method, genomic DNA from well-characterized cell lines were used; high-resolution KIR typing results obtained were then compared to those previously reported. Two different bioinformatic pipelines were used allowing the attribution of sequencing reads to specific KIR genes and the assignment of KIR alleles for each KIR gene. Our results demonstrated successful long-range KIR gene amplifications of all reference samples using intergenic KIR primers. The alignment of reads to the human genome reference (hg19) using BiRD pipeline or visualization of data using Profiler software demonstrated that all KIR genes were completely sequenced with a sufficient read depth (mean 317× for all loci) and a high percentage of mapping (mean 93% for all loci). Comparison of high-resolution KIR typing obtained to those published data using exome capture resulted in a reported concordance rate of 95% for centromeric and telomeric KIR genes. Overall, our results suggest that NGS can be used to investigate the broad KIR allelic polymorphism. Hence, these data improve our knowledge, not only on KIR⁺ NK cell alloreactivity in HSCT but also on the role of KIR⁺ NK cell populations in control of viral infections and diseases.

Limited HLA sequence variation outside of antigen recognition domain exons of 360 10 of 10 matched unrelated hematopoietic stem cell transplant donor-recipient pairs

Traditional DNA-based typing focuses primarily on interrogating the exons of human leukocyte antigen (HLA) genes that form the antigen recognition domain (ARD). The relevance of mismatching donor and recipient for HLA variation outside the ARD on hematopoietic stem cell transplantation (HSCT) outcomes is unknown. This study was designed to evaluate the frequency of variation outside the ARD in 10 of 10 (HLA-A, -B, -C, -DRB1, -DQB1) matched unrelated donor transplant pairs (n = 360). Next-generation DNA sequencing was used to characterize both HLA exons and introns for HLA-A, -B, -C alleles; exons 2, 3 and the intervening intron for HLA-DRB1 and exons only for HLA-DQA1 and -DQB1. Over 97% of alleles at each locus were matched for their nucleotide sequence outside of the ARD exons. Of the 4320 allele comparisons overall, only 17 allele pairs were mismatched for non-ARD exons, 41 for noncoding regions and 9 for ARD exons. The observed variation between donor and recipient usually involved a single nucleotide difference (88% of mismatches); 88% of the non-ARD exon variants impacted the amino acid sequence. The impact of amino acid sequence variation caused by substitutions in exons outside ARD regions in D-R pairs will be difficult to assess in HSCT outcome studies because these mismatches do not occur very frequently.

Allorecognition of HLA-C Mismatches by CD8+ T Cells in Hematopoietic Stem Cell Transplantation Is a Complex Interplay between Mismatched Peptide-Binding Region Residues, HLA-C Expression, and HLA-DPB1 Disparities

HLA-C locus mismatches (MMs) are the most frequent class I disparities in unrelated hematopoietic stem cell transplantation (HSCT) and have a detrimental impact on clinical outcome. Recently, a few retrospective clinical studies have reported some variability in the immunogenicity of HLA-C incompatibilities. To get better insight into presumably permissive HLA-C MMs, we have developed a one-way in vitro mixed lymphocyte reaction (MLR) assay allowing to quantify activated CD56^-CD137⁺CD8⁺ lymphocytes in HLA-C incompatible combinations. T cell-mediated alloresponses were correlated with genetic markers such as HLA-C mRNA expression and the number of amino acid (aa) MMs in the α1/α2 domains (peptide-binding region). Because of the high rate of HLA-DPB1 incompatibilities in HLA-A-, B-, C-, DRB1-, and DQB1-matched unrelated HSCT patient/donor pairs, the impact of HLA-DPB1 mismatching, a potential bystander of CD4⁺ T cell activation, was also considered. Heterogeneous alloresponses were measured in 63 HLA-C-mismatched pairs with a positive assay in 52% of the combinations (2.3-18.6% activated CTLs), representing 24 different HLA-A~B~DRB1~DQB1 haplotypes. There was no correlation between measured alloresponses and mRNA expression of the mismatched HLA-C alleles. The HLA-C*03:03/03:04 MM did not induce any positive alloresponse in five MLRs. We also identified HLA-C*02:02 and HLA-C*06:02 as mismatched alleles with lower immunogenicity, and HLA-C*14:02 as a more immunogenic MM. A difference of at least 10 aa residues known to impact peptide/T cell receptor (TCR) binding and a bystander HLA-DPB1 incompatibility had a significant impact on CTL alloreactivity (p = 0.021). The same HLA-C MM, when recognized by two different responders with the same HLA haplotypes, was recognized differently, emphasizing the role of the T-cell repertoire of responding cells. In conclusion, mismatched HLA-C alleles differing by 10 or more aas in the peptide/TCR-binding region, when occurring together with HLA-DPB1 incompatibilities, should be considered as high-risk MMs in unrelated HSCT.

Allele-Independent Turnover of Human Leukocyte Antigen (HLA) Class Ia Molecules

Major histocompatibility complex class I (MHCI) glycoproteins present cytosolic peptides to CD8+ T cells and regulate NK cell activity. Their heavy chains (HC) are expressed from up to three MHC gene loci (human leukocyte antigen [HLA]-A, -B, and -C in humans), whose extensive polymorphism maps predominantly to the antigen-binding groove, diversifying the bound peptide repertoire. Codominant expression of MHCI alleles is thus functionally critical, but how it is regulated is not fully understood. Here, we have examined the effect of polymorphism on the turnover rates of MHCI molecules in cell lines with functional MHCI peptide loading pathways and in monocyte-derived dendritic cells (MoDCs). Proteins were labeled biosynthetically with heavy water (2H2O), folded MHCI molecules immunoprecipitated, and tryptic digests analysed by mass spectrometry. MHCI-derived peptides were assigned to specific alleles and isotypes, and turnover rates quantified by 2H incorporation, after correcting for cell growth. MHCI turnover half-lives ranged from undetectable to a few hours, depending on cell type, activation state, donor, and MHCI isotype. However, in all settings, the turnover half-lives of alleles of the same isotype were similar. Thus, MHCI protein turnover rates appear to be allele-independent in normal human cells. We propose that this is an important feature enabling the normal function and codominant expression of MHCI alleles.

Expression of classical HLA class I molecules: regulation and clinical impacts: Julia Bodmer Award Review 2015

Human leukocyte antigen (HLA) class I genes are ubiquitously expressed, but in a tissue specific-manner. Their expression is primarily regulated at the transcriptional level and can be modulated both positively and negatively by different stimuli. Advances in sequencing technologies led to the identification of new regulatory variants located in the untranslated regions (UTRs), which could influence the expression. After a brief description of the mechanisms underlying the transcriptional regulation of HLA class I genes expression, we will review how the expression levels of HLA class I genes could affect biological and pathological processes. Then, we will discuss on the differential expression of HLA class I genes according to the locus, allele and UTR polymorphisms and its clinical impact. This interesting field of study led to a new dimension of HLA typing, going beyond a qualitative aspect.