HLA-E regulatory and coding region variability and haplotypes in a Brazilian population sample

Is HLA-E with its receptors an immune checkpoint or an antigenic determinant in allo-HCT?

Hematopoietic cell transplantation (HCT) represents a potentially curative therapeutic approach for various hematologic and non-hematologic malignancies. Human leukocyte antigen (HLA) matching is still the central selection criterion for HCT donors. Nevertheless, post-transplant complications, in particular graft-versus-host disease (GvHD), relapse of disease and infectious complications, represent a major challenge and contribute significantly to morbidity and mortality. Recently, non-classical HLA class I molecules, especially HLA-E, have gained increasing attention in the context of allogeneic HCT. This review aims to summarize the latest findings on the immunomodulatory role of HLA-E, which serves as a ligand for receptors of the innate and adaptive immune system. In particular, we aim to elucidate how (i) polymorphisms within HLA-E, (ii) the NKG2A/C axis and (iii) the repertoire of peptides presented by HLA-E jointly influence the functionality of immune effector cells. Understanding this intricate network of interactions is crucial as it significantly affects NK and T cell responses and thus clinical outcomes after HCT.

Establishment of Rapid Detection Methods for rs76971248 Related to Leukemia

Background

The HLA-E gene is a member of the HLA-I gene family. Its genetic polymorphism is regarded as associated with numerous diseases. Establishing a rapid and accurate detection method of disease-related SNP sites in HLA-E is particularly important.

Methods

Blood samples from 226 healthy blood donors and 228 leukemia patients were collected, and DNA was extracted. Three typing methods based on PCR-sequence-based typing, TaqMan genotyping, and high-resolution melting curve were established to identify rs76971248 (G>T). The Chi-square test was used for statistical analysis by SPSS.

Results

Three methods based on PCR-SBT, TaqMan genotyping, and HRM were all able to identify rs76971248. The software for analyzing the results of HLA-E sequencing was easy to use, and the results were accurate. The frequency of rs76971248 in different types of leukemia patients was significantly lower than that in healthy blood donors (p < 0.05). And the frequency of the G/G genotype in leukemia patients was significantly higher than that in healthy blood donors (p < 0.05).

Conclusions

For the screening of known SNP sites in large-scale populations, among the three methods, the TaqMan genotyping method had the advantage of shortest time consumption, simplest operation, and greatest specificity, which was the most appropriate method for this experiment. The analysis software for HLA-E gene sequencing needed to be further optimized. rs76971248 had a protective effect against leukemia. And the G/G genotype was a risk factor for leukemia.

Correlation of Human Leukocyte Antigen-E Genomic Polymorphism with Leukemia and Functional Study of Human Leukocyte Antigen-E Different Type Promoters

Human leukocyte antigen (HLA)-E is one of the least polymorphic nonclassical major histocompatibility complex (MHC) I genes; its nucleotide variability can affect immune response. In this study, we assess the correlation between HLA-E polymorphism and leukemia and further study the transcriptional activity of promoter variation at nucleotide position-26. A total of 142 healthy blood donors and 111 leukemia patients were collected. The genomic sequence of HLA-E was amplified by high-fidelity reaction system and identified by Sanger and cloning sequencing. The dual luciferase reporter gene assay was used to detect the transcription activity of promoter variation at nucleotide position-26. In the HLA-E genomic sequence results, a total of 16 alleles and 32 genotypes were detected; the HLA-E*01:01:01:06 allele had a significantly lower frequency in leukemia patients than in healthy participants (p = 0.026 < 0.05). And the HLA-E*01:03:02:01, *01:03:02:01 genotype showed the greatest difference in frequency between the two groups of participants (p = 0.028 < 0.05). Eight HLA-E alleles were first reported worldwide in Chinese individuals. The results of the dual luciferase reporter gene experiment showed that the transcription activity of the mutant-type promoter (HLA-E*01:01:01:06 with "T" allele at nucleotide position-26) was significantly lower compared with the wild-type promoter (HLA-E*01:01:01:01 with "G" allele at nucleotide position-26) (p = 0.0242 < 0.05). HLA-E*01:01:01:06 allele has a protective effect against leukemia through decreasing transcription activity by "T" variation at nucleotide position-26.

HLA Class I Molecules as Immune Checkpoints for NK Cell Alloreactivity and Anti-Viral Immunity in Kidney Transplantation

Natural killer (NK) cells are innate lymphocytes that can kill diseased- or virally-infected cells, mediate antibody dependent cytotoxicity and produce type I immune-associated cytokines upon activation. NK cells also contribute to the allo-immune response upon kidney transplantation either by promoting allograft rejection through lysis of cells of the transplanted organ or by promoting alloreactive T cells. In addition, they protect against viral infections upon transplantation which may be especially relevant in patients receiving high dose immune suppression. NK cell activation is tightly regulated through the integrated balance of signaling via inhibitory- and activating receptors. HLA class I molecules are critical regulators of NK cell activation through the interaction with inhibitory- as well as activating NK cell receptors, hence, HLA molecules act as critical immune checkpoints for NK cells. In the current review, we evaluate how NK cell alloreactivity and anti-viral immunity are regulated by NK cell receptors belonging to the KIR family and interacting with classical HLA class I molecules, or by NKG2A/C and LILRB1/KIR2DL4 engaging non-classical HLA-E or -G. In addition, we provide an overview of the methods to determine genetic variation in these receptors and their HLA ligands.

Peripheral spectrum neurological disorder after arbovirus infection is associated with HLA-F variants among Northeastern Brazilians

INTRODUCTION

Non-classical class I human leukocyte antigens (HLA) molecules are known to modulate the function of cytotoxic cells (NK and T CD8+) during viral infection by interacting with inhibitory/activating receptors. However, little is known about the HLA-E/-F genetic variability on arbovirus infections.

METHODS

We evaluated by massive parallel sequencing the full HLA-E/-F genetic diversity among patients infected during the arbovirus (ZIKV, DENV, and CHIKV) outbreak leading to a broad range of neurological complications in the Brazilian State of Pernambuco. In parallel, healthy blood donors from the same area were also studied. Plink and R software were used for genetic association study. To limit the false-positive results and enhance the reliability of the results, we adopted P-values <0.01 as significant levels.

RESULTS

Compared to controls, the HLA-F alleles: -1610 C (rs17875375), +1383 G (rs17178385), and +3537 A (rs17875384), all in complete linkage disequilibrium with each other (r² = 1), were overrepresented in patients presenting peripheral spectrum disorders (PSD). The HLA-F*Distal-D haplotype that harbored the -1610 C allele exhibited a trend increase in PSD group. No associations were found for HLA-E.

CONCLUSIONS

Our findings showed that the HLA-F genetic background seems to be more important than HLA-E on the susceptibility to PSD complications.

Human leukocyte antigen (HLA)-F and -G gene polymorphisms and haplotypes are associated with malaria susceptibility in the Beninese Toffin children

BACKGROUND

Little attention has been devoted to the role of the immunoregulatory HLA-E/-F/-G genes in malaria. We evaluated the entire HLA-E/-F/-G variability in Beninese children highly exposed to Plasmodium falciparum (P.f.) malaria.

METHODS

154 unrelated children were followed-up for six months and evaluated for the presence and number of malaria episodes. HLA-E/-F/-G genes were genotyped using massively parallel sequencing. Anti P.f. antibodies were evaluated using ELISA.

RESULTS

Children carrying the G allele at HLA-F (-1499,rs183540921) showed increased P.f. asymptomatic/symptomatic ratio, suggesting that these children experienced more asymptomatic P.f. episodes than symptomatic one. Children carrying HLA-G-UTR-03 haplotype exhibited increased risk for symptomatic P.f. episodes and showed lower IgG2 response against P.f. GLURP-R2 when compared to the non-carriers. No associations were observed for the HLA-E gene.

CONCLUSION

HLA-F associations may be related to the differential expression profiles of the encoded immunomodulatory molecules, and the regulatory sites at the HLA-G 3'UTR may be associated to posttranscriptional regulation of HLA-G and to host humoral response against P.f.

HLA-E gene polymorphisms in chronic hepatitis C: Impact on HLA-E liver expression and disease severity

Hepatitis C virus usually produces chronic infection and liver damage. Considering that: i) the human leukocyte antigen-E (HLA-E) molecule may modulate the immune response, and ii) little is known about the role of HLA-E gene variability on chronic hepatitis C, we studied the impact of HLA-E polymorphisms on the magnitude of HLA-E liver expression and severity of hepatitis C. HLA-E variability was evaluated in terms of: i) single nucleotide polymorphism (SNP) alleles and genotypes along the gene (beginning of the promoter region, coding region and 3'UTR), and ii) ensemble of SNPs that defines the coding region alleles, considered individually or as genotypes. The comparisons of the HLA-E variation sites between patients and controls revealed no significant results. The HLA-E + 424 T > C (rs1059510), +756 G > A (rs1264457) and + 3777 G > A (rs1059655) variation sites and the HLA-E*01:01:01:01 and HLA-E*01:03:02:01 alleles, considered at single or double doses, were associated with the magnitude of HLA-E liver expression in Kupfer cell, steatosis, inflammatory activity and liver fibrosis. Although these associations were lost after corrections for multiple comparisons, these variable sites may propitiate biological clues for the understanding of the mechanisms associated with hepatitis C severity.

Hla-C genetic diversity and evolutionary insights in two samples from Brazil and Benin

Human leukocyte antigen-C (HLA-C) is a classical HLA class I molecule that binds and presents peptides to cytotoxic T lymphocytes in the cell surface. HLA-C has a dual function because it also interacts with Killer-cell immunoglobulin-like receptors (KIR) receptors expressed in natural killer and T cells, modulating their activity. The structure and diversity of the HLA-C regulatory regions, as well as the relationship among variants along the HLA-C locus, are poorly addressed, and few population-based studies explored the HLA-C variability in the entire gene in different population samples. Here we present a molecular and bioinformatics method to evaluate the entire HLA-C diversity, including regulatory sequences. Then, we applied this method to survey the HLA-C diversity in two population samples with different demographic histories, one highly admixed from Brazil with major European contribution, and one from Benin with major African contribution. The HLA-C promoter and 3'UTR were very polymorphic with the presence of few, but highly divergent haplotypes. These segments also present conserved sequences that are shared among different primate species. Nucleotide diversity was higher in other segments rather than exons 2 and 3, particularly around exon 5 and the second half of the 3'UTR region. We detected evidence of balancing selection on the entire HLA-C locus and positive selection in the HLA-C leader peptide, for both populations. HLA-C motifs previously associated with KIR interaction and expression regulation are similar between both populations. Each allele group is associated with specific regulatory sequences, reflecting the high linkage disequilibrium along the entire HLA-C locus in both populations.

Single molecule real-time DNA sequencing of the full HLA-E gene for 212 reference cell lines

We have developed a genotyping assay that produces fully phased, unambiguous HLA-E genotyping using Pacific Biosciences' single molecule real-time DNA sequencing. In total 212 cell lines were genotyped, including the panel of 107 established at the 10th International Histocompatibility Workshop. Our results matched the previously known HLA-E genotype in 94 (44.3%) cell lines, in all cases either improving or equalling previous genotyping resolution. Three (1.4%) cells had discrepant HLA-E genotyping data and 115 (54.2%) had no previous HLA-E data. The HLA-E genotypes for four (1.9%) cell lines resulted in a change of zygosity by identifying two distinct haplotypes. We discovered eight novel HLA-E alleles, extended the known reference sequence of seven and confirmed the existence of a further 10.

HLAIb worldwide genetic diversity: New HLA-H alleles and haplotype structure description

The classical HLA class I genes (HLA Ia) were extensively studied because of their implication in clinical fields and anthropology. Less is known about worldwide genetic diversity and linkage disequilibrium for non-classical HLA class I genes (HLA Ib) and HLA pseudogenes. Notably, HLA-H, which is deleted in a fraction of the population, remains scarcely explored. The aims of this study were 1/ to get further insight into HLA-H genetic diversity and into how this variability potentially affects its expression and 2/ to define HLA Ib worldwide allelic diversity and linkage. Exome sequence data from the 1000 Genomes Project were used to define second field HLA-A, -E, -F, -G and -H typing using PolyPheMe software. Allelic and two-loci haplotype frequencies were estimated using Gene[Rate] software both at worldwide and continental levels. Eleven novel HLA-H alleles identified in exome data were validated by NGS performed on 25 genomic DNA samples from the same cohort. Phylogenetic analysis and frequency distribution of HLA-H alleles revealed three clades, each predominantly represented in Admixed American, European and East Asian populations, African populations and South Asian populations. Among these eleven novel alleles, two potentially encode complete transmembrane HLA proteins. We confirm the high LD between HLA-H and -A, and between HLA-H and -G, and show the three genes have distinct worldwide allelic distribution. Conversely, HLA-E and HLA-F both showed little LD, displayed restricted allelic diversity and practically no difference in their distribution across the planet. Our work thus reveals an unexpectedly high HLA-H genetic diversity, with alleles highly represented in Asia possibly encoding a functional HLA protein. Functional implication of these results remains to be explored, both in physiological and pathological contexts.

Peptides as cancer vaccines

Cancer vaccines based on synthetic peptides are a safe, well-tolerated immunotherapy able to specifically stimulate tumor-reactive T cells. However, their clinical efficacy does not approach that achieved with other immunotherapies such as immune checkpoint blockade. Nevertheless, major advances have been made in selecting tumor antigens to target, identifying epitopes binding to classical and non-classical HLA molecules, and incorporating these into optimal sized peptides for formulation into a vaccine. Limited potency of currently used adjuvants and the immunosuppressive tumor microenvironment are now understood to be major impediments to vaccine efficacy that need to be overcome. Rationally designed combination therapies are now being tested and should ultimately enable peptide vaccination to be added to immuno-oncology treatment options.

HLA-G, -E and -F regulatory and coding region variability and haplotypes in the Beninese Toffin population sample

HLA-G/E/F genes exhibit immunomodulatory properties and are expressed in placenta. Little attention has been devoted to the study of these genes in sub-Saharan African populations, which are yet the most diverse. To fill this gap, we evaluated the complete gene variability, approximately 5.1 kb for HLA-G (n = 149), 7.7 kb for HLA-E (n = 150) and 6.2 kb for HLA-F (n = 152) in the remote Beninese Toffin population, using massive parallel sequencing. Overall, 96, 37 and 68 variable sites were detected along the entire HLA-G, -E and -F, respectively, arranged into region-specific haplotypes; i.e., promoter haplotypes (16, 19, and 15 respectively), coding haplotypes (19, 15, and 29 respectively), 3' untranslated region (3'UTR) haplotypes (12, 7 and 2, respectively) and extended haplotypes (33, 31 and 32 respectively). All promoter/coding/3'UTR haplotypes followed the patterns already described in worldwide populations. HLA-E was the most conserved, exhibiting mainly two full-length encoded-molecules (E*01:01 and E*01:03), followed by HLA-F, three full-length proteins (F*01:01, F*01:02 and F*01:03) and HLA-G, four proteins: three full-length (G*01:01, G*01:03 and G*01:04) and one truncated (G*01:05N). Although HLA-G/E/F alleles in the Toffin population were the most frequently observed worldwide, the frequencies of the coding haplotypes were closely similar to those described for other African populations (Guinea-Conakry and Burkina-Faso), when compared to non-African ones (Brazilian), indicating that variable sites along these genes were present in Africa before human dispersion.

Hla-mapper: An application to optimize the mapping of HLA sequences produced by massively parallel sequencing procedures

A challenging task when more than one HLA gene is evaluated together by second-generation sequencing is to achieve a reliable read mapping. The polymorphic and repetitive nature of HLA genes might bias the read mapping process, usually underestimating variability at very polymorphic segments, or overestimating variability at some segments. To overcome this issue we developed hla-mapper, which takes into account HLA sequences derived from the IPD-IMGT/HLA database and unpublished HLA sequences to apply a scoring system. This comprehends the evaluation of each read pair, addressing them to the most likely HLA gene they were derived from. Hla-mapper provides a reliable map of HLA sequences, allowing accurate downstream analysis such as variant calling, haplotype inference, and allele typing. Moreover, hla-mapper supports whole genome, exome, and targeted sequencing data. To assess the software performance in comparison with traditional mapping algorithms, we used three different simulated datasets to compare the results obtained with hla-mapper, BWA MEM, and Bowtie2. Overall, hla-mapper presented a superior performance, mainly for the classical HLA class I genes, minimizing wrong mapping and cross-mapping that are typically observed when using BWA MEM or Bowtie2 with a single reference genome.