Association of HLA-A and Non-Classical HLA Class I Alleles

Federico Carlini, Virginia Ferreira, Stéphane Buhler, Audrey Tous, Jean-François Eliaou, Céline René, Jacques Chiaroni, Christophe Picard, Julie Di Cristofaro

► To cite this version:

Federico Carlini, Virginia Ferreira, Stéphane Buhler, Audrey Tous, Jean-François Eliaou, et al. As-sociation of HLA-A and Non-Classical HLA Class I Alleles. PLoS ONE, 2016, 11 (10), pp.e0163570. 10.1371/journal.pone.0163570. hal-01447663

HAL Id: hal-01447663

https://hal.science/hal-01447663v1

Submitted on 3 Oct 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

CrossMark

click for updates

RESEARCH ARTICLE

# Association of HLA-A and Non-Classical HLA Class I Alleles

Federico Carlini¹,², Virginia Ferreira², Stéphane Buhler³,⁴, Audrey Tous¹,², Jean-François Eliaou⁵,⁶, Céline René⁵,⁶, Jacques Chiaroni¹,², Christophe Picard¹,², Julie Di Cristofaro¹,²*

1 Etablissement Français du Sang Alpes Méditerranée, Marseille, France, 2 Aix Marseille Univ, CNRS, EFS, ADES, "Biologie des Groupes Sanguins", Marseille, France, 3 Laboratory of Anthropology, Genetics and Peopling history (AGP), Department of Genetics and Evolution-Anthropology Unit, University of Geneva, Geneva, Switzerland, 4 Transplantation Immunology Unit and National Reference Laboratory for Histocompatibility, Department of Genetic and Laboratory Medicine, Geneva University Hospitals, Geneva, Switzerland, 5 Department of Immunology, CHRU de Montpellier, University Hospital Saint-Eloi, Montpellier, France, 6 Faculté de Médecine, University of Montpellier 1, Montpellier, France

* julie.dicristofaro@efs.sante.fr

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Carlini F, Ferreira V, Buhler S, Tous A, Eliaou J-F, René C, et al. (2016) Association of HLA-A and Non-Classical HLA Class I Alleles. PLoS ONE 11(10): e0163570. doi:10.1371/journal.pone.0163570

**Editor:** Antony Nicodemus Antoniou, University of East London, UNITED KINGDOM

**Received:** March 2, 2016

**Accepted:** September 11, 2016

**Published:** October 4, 2016

**Copyright:** © 2016 Carlini et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** The author(s) received no specific funding for this work.

**Competing Interests:** The authors have declared that no competing interests exist.

**Abbreviations:** HLA, Human Leukocyte Antigen; CTL, Cytotoxic T-Lymphocyte; NK, Natural Killer cells; MHC, Major Histocompatibility Complex; LD, Linkage Disequilibrium; ILT, Immunoglobulin-Like Transcript; KIR, Killer-cell Immunoglobulin-like

## Abstract

The HLA-A locus is surrounded by HLA class Ib genes: HLA-E, HLA-H, HLA-G and HLA-F. HLA class Ib molecules are involved in immuno-modulation with a central role for HLA-G and HLA-E, an emerging role for HLA-F and a yet unknown function for HLA-H. Thus, the principal objective of this study was to describe the main allelic associations between HLA-A and HLA-H, -G, -F and -E. Therefore, HLA-A, -E, -G, -H and -F coding polymorphisms, as well as HLA-G UnTranslated Region haplotypes (referred to as HLA-G UTRs), were explored in 191 voluntary blood donors. Allelic frequencies, Global Linkage Disequilibrium (GLD), Linkage Disequilibrium (LD) for specific pairs of alleles and two-loci haplotype frequencies were estimated. We showed that HLA-A, HLA-H, HLA-F, HLA-G and HLA-G UTRs were all in highly significant pairwise GLD, in contrast to HLA-E. Moreover, HLA-A displayed restricted associations with HLA-G UTR and HLA-H. We also confirmed several associations that were previously found to have a negative impact on transplantation outcome. In summary, our results suggest complex functional and clinical implications of the HLA-A genetic region.

## Introduction

A complete Human Leukocyte Antigen (HLA) match is associated with better long term survival in transplant patients, however the effects of HLA-A, -B and -DRB1 matching are unequal and also depend upon the organ being transplanted. HLA-A mismatches were found to have less influence on kidney allografts than HLA-B mismatches [1]. This difference could reflect the higher number of HLA-B alleles compared to HLA-A and/or the cytotoxic T-cell allo-repertoire for HLA-A and -B antigens, as lower frequencies of Cytotoxic T-Lymphocyte precursor (CTLp) are found for HLA-A mismatches than for HLA-B mismatches [2, 3]. Few amino acid
Receptor; ER, Endoplasmic Reticulum; OC, Open
Conformer; HSCT, Hematopoietic Stem Cell
Transplantation; β2m, β2-microglobulin; UTR,
UnTranslated Region; URR, Upstream Regulatory
Region; LTX, Lung Transplantation; GVHD, Graft
Versus Host Disease; CLAD, Chronic Lung Allograft
Dysfunction; gDNA, genomic DNA.

polymorphisms at functionally important positions of the antigen recognition site of the HLA-A molecule have been shown to have a significant influence on clinical outcome [4]. In liver transplant patients however, only HLA-A mismatching was associated to a negative effect on patient survival and a higher recurrence of the hepatitis C virus, whereas no HLA-B and -DR mismatches had no effect [5]. In Lung Transplantation (LTx), the risk of Bronchiolitis Obliterans Syndrome only increased in the presence of two HLA-A mismatches, not for 0 or 1 mismatch or for HLA-B or -DR mismatches [6]. These data support the hypothesis that HLA-A mismatching leads to a particular immunological situation in organ transplantation compared to other loci. The HLA-A molecule displays specific features compared to HLA-B that could account for its particular role: HLA-A and HLA-B alleles carry an unpaired cysteine at different positions of the cytoplasmic tail domain (respectively at positions 339 and 325), which is reported to be involved in recycling, targeting for degradation and influencing recognition by NK receptors as well as in the formation of fully folded MHC class I dimers in exosomes [7]. Other genes, in association/linkage disequilibrium with HLA-A, could also induce either alloreactivity or tolerance. Indeed, the HLA-A locus on chromosome 6 is surrounded by HLA class Ib (non-classical) genes: HLA-E (at the centromeric end), HLA-G, -H and HLA-F (at the telomeric end).

HLA class Ib molecules, HLA-G, HLA-E, HLA-H and HLA-F, are described for their involvement in immune effector cells interactions. They display specific features that differentiate them from HLA class Ia (classical) molecules: a highly conserved peptide binding groove with restricted coding polymorphism and specific expression patterns [8, 9]. Moreover, HLA-G lacks the canonical MHC class I cytoplasmic tail as it is truncated due to the deletion of 6 amino acids. Yet, HLA-G is capable of initiating cellular signaling events upon cross-linking, as HLA-G protein associates with lipid rafts and thus can induce cell proliferation [10, 11].

The immuno-modulating and proliferative properties of HLA class Ib are of great interest in the clinical field, particularly in transplantation. Notably, higher HLA-G expression has been correlated to a better clinical outcome after heart, kidney and Hematopoietic Stem Cell Transplantation (HSCT) (see [12] and [13] for reviews). However, the exploration of HLA-G protein expression is currently limited because micro-environmental expression analysis requires invasive techniques. A further complication is that HLA-G is expressed as many different protein isoforms, generated by alternative splicing, dimer formation, and association with β₂ microglobulin. Unfortunately, these isoforms are not accurately discriminated by commercial antibodies whereas, most importantly, they are reported to display differential functionality. Primary HLA-G mRNA splicing generates both membrane-bound (HLA-G1 to 4) and soluble (HLA-G5 to 7) isoforms [14] and membrane-bound isoforms can also be shed via proteolytic cleavage by matrix metalloproteinase-2 [15]. The HLA-G1 and -G5 isoforms display a similar structure to HLA class Ia molecules (comprising α1, α2 and α3 domains), whereas HLA-G2 and -G6 are only composed of α1 and α3 domains, -G3 and -G7 contain an α1 domain, and -G4 is composed of α1 and α2 domains [14]. These truncated isoforms appear to be differently processed as HLA-G2, -G3 and -G4 are sensitive to endoglycosidase H, suggesting a non-involvement of the Golgi apparatus, in contrast to HLA-G1 and classical HLA class I [16]. Both HLA-G1 and -G5 isoforms reduce NK cell cytotoxicity with an additive effect [17]. Two studies concluded functional activity for HLA-G2, -G4 and -G6 [16], whereas conflicting results were published concerning HLA-G3 [16, 18]. Besides these differentially spliced isoforms, HLA-G forms disulfide-linked dimers via Cys-42 of the HLA-G α1 domain. Most of the dimers are reported to be expressed on the cell surface, but soluble HLA-G is also able to dimerize [19]. HLA-G dimerization, which may depend on cell-surface density, has an impact on HLA-G receptor affinity [20]. Hence, possibly due to the many HLA-G isoforms and the lack of appropriate serological tools, conflicting results have been reported regarding the
association between HLA-G expression and clinical outcome [21–24]. As a consequence, many studies focus on HLA-G genetic polymorphisms to predict clinical outcome.

The many HLA-G regulatory region polymorphisms (5'URR and 3'UTR) are in strong Linkage Disequilibrium (LD) with each other, forming haplotypes referred to as HLA-G UTRs in this manuscript [25]. These HLA-G UTRs are well conserved in different physiological populations, as well as in different clinical cohorts [22–24]. Several clinical studies have described favorable or deleterious effects of HLA-G alleles and HLA-G UTRs. For example, in LTx, HLA-G*01:06~UTR2 was associated with a worse evolution of cystic fibrosis and HLA-G*01:04~UTR3 impaired long term survival [24]. HLA-G*01:06 was associated with pregnancy complications in French [26] and Singaporeans women [27], and an increase of HLA-G*01:04 frequency was observed in couples with recurrent miscarriages [28, 29]. Conversely, HLA-G*01:04 has been shown to be protective in acute renal rejection and end stage renal disease (chronic inflammatory disease) [30].

Another non-classical HLA class I molecule, HLA-E, regulates NK cell function and is the only known ligand for the C-type lectin receptor CD94 combined with different NKG2 sub-units expressed on NK and CD8+ αβ T cells [31, 32]. HLA-E mRNA is expressed in most tissues [33] and its cell surface expression is controlled by leader peptides of HLA class I molecules, except HLA-F, the highest affinity being reported for HLA-G. Similarly to HLA-G, HLA-E's intrinsic peptide repertoire is constrained by a larger number of primary anchor sites compared to most class Ia molecules, which further supports a specific immunological role for this molecule [34]. HLA-E can also bind peptide ligands from stress proteins and viruses [9, 35]. Although HLA-E possesses less anchor sites than HLA-G, the HLA-E peptide repertoire appears to be more restricted than that of HLA-G. Equal frequencies of the HLA-E coding alleles (HLA-E*01:01 and HLA-E*01:03, R107G, respectively associated to low and high protein expression) suggest that balancing selection may have maintained heterozygosis of high- and low-expressing alleles [36–39]. In HSCT, HLA-E*01:03 homozygosity was associated with a lower risk of Graft Versus Host Disease (GVHD), decreased mortality and a higher disease-free survival rate, whereas HLA-E*01:01 homozygosity was associated with an increased risk of bacterial infection (See [13] for a review). In LTx, HLA-E*01:03 was correlated to Chronic Lung Allograft Dysfunction (CLAD) occurrence and HLA-E homozygosity was associated with a worse survival compared to heterozygosity (Di Cristofaro et al, submitted).

Much less is known about the other HLA class Ib molecules. HLA-F mRNA is expressed in most cell types and the protein is localized in the Endoplasmic Reticulum (ER) and Golgi apparatus. Unlike other HLA class I molecules, the function of HLA-F seems to be independent of peptide loading in the ER [31, 40], as it is exclusively expressed in an open conformer form (OC, deprived of binding peptide and/or β2 microglobulin) since no peptide could be eluted from both membrane-bound and intracellular HLA-F captured proteins [41]. HLA-F is expressed at the surface of the extravillous trophoblast and is upregulated upon activation in monocytes, NK, B and T cells, but not T-Reg CD4+ CD25+ cells [33, 42]. HLA-F seems to possess different functions: interaction with ILT-2 and ILT-4 and with KIR receptors [9, 43]; interaction with HLA-E in the OC form, possibly modifying the interaction with its receptors [44], and an implication in an original pathway of antigen (Ag) cross-presentation by HLA class I molecules [45]. HLA-F is characterized by 4 alleles defined at a 4 digit resolution, with F*01:01 representing 90% of allelic diversity.

Finally, HLA-H has been defined as a pseudogene because of a deletion in exon 4, causing a premature stop codon before cysteine 164, impairing the Ag-presenting function [46]. Of note, HLA-H may be mistaken for the HFE gene (associated with haemochromatosis) [9, 47, 48], but, although the HFE protein is similar to the proteins of the MHC class I family, it does not seem to inhibit NK cell activity [48]. Increased HLA-H nucleotide variation could be the
consequence of balancing selection acting on nearby loci or of relaxed functional pressures on this pseudogene [49]. No data is currently available concerning HLA-F or HLA-H alleles and clinical outcome.

Few data exist concerning the associations between HLA class Ib alleles; all located within a region of 724 kb on chromosome 6 (Table 1) [50–52]. So far, the studies on HLA-E reported no LD with other HLA class Ib loci [35], except with HLA-G in isolated populations [53]. However, some HLA-E and HLA-A alleles have been described as being fully associated, e.g. HLA-A*01 and HLA-E*01:01, or HLA-A*23/24/29 and HLA-E*01:03 [36]. HLA-A and HLA-G also display LD [54–56]: HLA-A*23/24 is associated to HLA-G*01:04 and this haplotype displays a deletion of more than 50 kb, comprising HLA-H [51, 57–59]. No data is available concerning HLA-F allelic associations with other HLA class I genes.

Many interactions are described between HLA class Ib molecules and immunological effectors cells, particularly NK cells. Thus, the principal objective of this study is to describe the main allelic associations between the HLA class Ib loci and HLA-A in a healthy population and to address their putative involvement in the synergy of these proteins.

## Materials and Methods

### Clinical samples

Ethylene Diamine Tetra Acetate-anticoagulated (EDTA) peripheral blood samples and sera were collected after written informed consent from 191 voluntary blood donors from south-eastern France. The donations were collected in accordance with the French blood donation regulations and ethics and with the French Public Health Code (art L1221-1). The samples were approved by the French Committee for the Protection of Persons in Biomedical Research (CCPPRB) and the entire French collection was also declared to and approved by the French Ministry of Higher Education and Research. Blood samples were anonymized according to the French Blood Center (Etablissement Français du Sang) procedure. Genomic DNA (gDNA) was extracted from a 200-µl total blood sample using the QIAmp Blood DNA Mini kit (Qiagen, Courtaboeuf, France) according to manufacturer's instructions.

### HLA-A genotyping

Luminex™ technology (HLA-A-One Lambda LABType® SSO, InGen, Chilly Mazarin, France) was used to determine HLA-A alleles at a low resolution level, i. e. 2-digit or first field level typing, using the manufacturer's kit.

### HLA-G, 5'URR and 3'UTR genotyping

A home-made primer extension method described in [22, 53] was used to simultaneously analyze 16 SNPs of the HLA-G gene: 8 SNPs (codons 13C/T (rs17875397); 31 A/T (rs41551813);

<table><caption>Table 1. HLA-E, -H, -G, -F and HFE relative location and distance from HLA-A on chromosome 6.</caption><thead><tr><td>Locus</td><td>NCBI Gene ID</td><td>Start (bp from pter)</td><td>End (bp from pter)</td><td>Size (bp)</td><td>Distance from HLA-A (Kbp)</td><td>Location relative to HLA-A</td></tr></thead><tbody><tr><td>HLA-E</td><td>3133</td><td>30489406</td><td>30494205</td><td>4799</td><td>547</td><td>Centromeric</td></tr><tr><td>HLA-A</td><td>3105</td><td>29942470</td><td>29945884</td><td>3414</td><td>0</td><td>-</td></tr><tr><td>HLA-H</td><td>3136</td><td>29887760</td><td>29891080</td><td>3320</td><td>55</td><td>Telomeric</td></tr><tr><td>HLA-G</td><td>3135</td><td>29826979</td><td>29831122</td><td>4143</td><td>115</td><td>Telomeric</td></tr><tr><td>HLA-F</td><td>3134</td><td>29887760</td><td>29891080</td><td>3320</td><td>55</td><td>Telomeric</td></tr><tr><td>HFE</td><td>3077</td><td>26087281</td><td>26096117</td><td>8836</td><td>3855</td><td>Telomeric</td></tr></tbody></table>

doi:10.1371/journal.pone.0163570.t001
54 A/G (rs41545515); 110 C/A (rs12722477); 130 C/T (rs41557518); 159 T/C (rs55916353); 219 C/T (rs45530733) and 258 C/T (rs12722482)) defining HLA-G*01:01 to G*01:09 alleles, 4 SNPs in the 5' URR region (-725 C/G/T (rs1233334), -716 G/T (rs2249863), -201 G/A (rs1233333) and -56 C/T (rs17875397)) and 4 SNPs in the 3'UTR region (exon 8 ins/del (rs66554220); 3142 C/G (rs1063320); 3187 G/A (rs9380142) and 3196 C/G (rs1610696)) defining HLA-G UTR-1 to HLA-G UTR-8. HLA-G genotypes were analyzed using GeneMapper® 4.0 with specific detection parameters.

## HLA-E and HLA-F genotyping

Next Generation Sequencing was used to genotype HLA-E and HLA-F alleles. A fragment of both genes was amplified by PCR with primers designed using the « Primer 3 » program (http://bioinfo.ut.ee/primer3-0.4.0/primer3/). The HLA-E primers flanked the genomic sequence from position 257 to 3072 [2] allowing the typing of E*01:01:01:01/E*01:01:01:02 (non resolved because the ambiguities are at positions -188 and -187, outside the amplified region), E*01:01:01:03, E*01:03:01:01, E*01:03:01:02, E*01:03:02:01, E*01:03:02:02, E*01:03:04 and E*01:06 alleles. The forward and reverse primer sequences were 5'TCATCTCTGTGGG CTACGTG3' and 5'TCAGACCCCCAGAATCTCAC3', respectively. The HLA-F primers flanked the genomic sequence from position -266 to 3249 [2] allowing the genotyping of all of the HLA-F alleles described to date. The forward and reverse primer sequences were 5'GTGGCTCTCAAGGCTCAG3' and 5'TACGATGCTTTGGTTGTTGC3', respectively. Both PCR amplifications were performed using the Long Range PCR Kit (Qiagen, Courta-boeuf, France) with 200 ng gDNA in a final volume of 25 µL containing 1x PCR buffer, 0.5 mM of dNTPs, 2 units of Long Range PCR Enzyme Mix and 0.4 µM of each primer. Amplification was carried out as follows: 1 cycle at 95°C for 3 min; 35 cycles at 95°C for 15 sec, 65°C for 30 sec, and 68°C for 5 min. Following PCR amplification, PCR fragments (HLA-E: 2815 bp, HLA-F: 3515 bp) were sequenced using a Next Generation Sequencing platform according to manufacturer's instructions (MiSeq, Illumina, The Netherlands).

## HLA-H genotyping

Direct sequencing was used to determine HLA-H alleles. A fragment of the HLA-H gene was amplified by PCR with primers designed using the « Primer 3 » program. The primers flanked the genomic sequence from position 284 to 488 [2] allowing the typing of H*01:01:01:01/H*01:01:01:03 (non resolved because the ambiguities are at positions 1252, 1261, 1415, 1427, 1429, 1430 and 1881, outside the amplified region), H*01:01:01:02, H*01:02, H*02:01:01:01, H*02:01:01:02, H*02:02, H*02:03, H*02:04 and H*del alleles. The forward and reverse primer sequences were 5' ATCTCCGTCGGCTACGTG 3' and 5' CGAGTCTCTGGGTCCGAGAT 3', respectively. PCR amplification was performed on 50 ng of gDNA in a final volume of 25 µL containing 1x PCR buffer, 1.5 mM MgCl2, 0.2 mM of dNTPs, 0.1 unit of Taq DNA-polymerase (Invitrogen, Cergy Pontoise) and 0.16 µM of each primer. Amplification was carried out as follows: 1 cycle at 95°C for 5 min; 30 cycles at 95°C for 30 sec, 57°C for 30 sec, and 72°C for 1 min; and 1 cycle at 72°C for 7 min. PCR amplification was checked by agarose electrophoresis and PCR fragments (length: 293 bp) were sequenced using a Big Dye Terminator V1.1 kit (Life Technologies) with each PCR primer according to the manufacturer's protocol. Sequences were analyzed using the Codon Code Aligner program (Codon Code Corporation, Massachusetts).

## HLA allelic assignment

The HLA-A, G, H, F and E allelic assignments were based on the HLA sequences listed in the official IMGT/HLA database [2].
HLA-G and HLA-H SNPs were automatically converted from output files (.txt) exported from GeneMapper 4.0 (HLA-G) or from Codon Code Aligner (HLA-H) into alleles and HLA-G UTRs using an in-house computer program readable by the 'Phenotype' application of the Gene[Rate] computer tool package (http://hla-net.eu/tools/) [60-62]. The HLA-G interpretation table used by Gene[Rate] is described in [53] and the HLA-H interpretation table is described in S1 Table.

HLA-E and HLA-F genotypes were assigned by the NGSengine program (GenDx, The Netherlands).

For allelic and two or more loci frequency estimations, all HLA-G and HLA-A putative homozygotes were considered either true homozygotes or heterozygotes for both the observed allele and an undefined or undetectable ('blank') allele. Due to the deletion encompassing the HLA-H locus, all HLA-H homozygous samples were encoded either as true homozygotes, as heterozygotes for both the observed allele and a 'blank' allele or as heterozygotes for both the observed allele and a deleted ('del') allele. For HLA-E and HLA-F, this procedure was not applied because typing was performed from NGS data, which allows the identification of any unknown allele (clonal amplification).

## Statistical analysis

Missing data at a locus led to the exclusion of the concerned sample from further analyses at the given locus. No multiple imputations were used.

Allelic and two or more loci haplotype frequencies were estimated using an EM algorithm implemented in the Gene[Rate] computer tools [62]. Deviations from Hardy-Weinberg equilibrium (HWE) were tested using a nested likelihood model [63]. Two-loci Global Linkage Dis-equilibrium (GLD) was assessed by a likelihood-ratio test on the frequency estimations [60, 61]. LD for specific pairs of alleles was provided as a list of standardized residuals [64] for each observed haplotype and a value greater than |2| was considered to be a significant deviation [60, 61].

## Results

### HLA-A, -E, -H, -F and -G allelic frequencies

HLA-A, -E, -H, -F and -G allele frequencies are given in Table 2. The null hypothesis of Hardy-Weinberg equilibrium was never rejected. Sixteen 2-digit HLA-A alleles were found, and their frequencies were concordant with previously published data in France [61]), HLA-A\*02, HLA-A\*03, HLA-A\*24 and HLA-A\*01 being the most frequent alleles. Six 4-digit HLA-G alleles and 8 HLA-G UTRs were observed and their frequencies were similar to previously published data [22, 53]. Putative undefined or undetectable blank 'alleles' represented at the most 1.6% of HLA-G UTRs according to the EM estimation (Table 2). Seven 8-digit HLA-E alleles were found; their frequencies were concordant with published data showing an equal distribution between HLA-E\*01:01 and HLA-E\*01:03, with HLA-E\*01:03:02:01 being the most frequent allele in Western Europe [53, 65]. The HLA-E\*01:03:01:02 and HLA-E\*01:03:04 alleles were not observed, while HLA-E\*01:06 and HLA-E\*01:03:02:02 displayed very low frequencies (<0.5%). Eight 6-digit HLA-H alleles were found, HLA-H\*01:01:01 and HLA-H\*02:04 being the most frequent. A new HLA-H allele (H\*02:04new) was identified at a frequency of 8.7% in this study, defined by a silent G>A substitution at position 368 compared to the IMGT gDNA reference for H\*02:04 (which was observed at a frequency of 15.9%). The HLA-H\*del was found at a frequency of 13.9%. HLA-H\*01:01:01:02, H\*01:01:01:03 and H\*02:01:01:02 were not observed and are further referred to as H\*01:01:01 and H\*02:01:01, respectively. Finally, four
<table><caption>Table 2. HLA-A, HLA-G UTRs, -G, -E, -H and -F allelic frequencies in 191 voluntary blood donors from southeastern France.</caption><thead><tr><th>Allele</th><th>Fq (%)</th><th>Allele</th><th>Fq (%)</th></tr></thead><tbody><tr><td>A*02</td><td>31.6</td><td>G*01:01</td><td>78.3</td></tr><tr><td>A*03</td><td>15.0</td><td>G*01:04</td><td>11.5</td></tr><tr><td>A*24</td><td>10.8</td><td>G*01:06</td><td>3.9</td></tr><tr><td>A*01</td><td>10.3</td><td>G*01:03</td><td>3.9</td></tr><tr><td>A*11</td><td>7.1</td><td>G*01:05N</td><td>2.1</td></tr><tr><td>A*29</td><td>6.6</td><td>G*01:08</td><td>0.3</td></tr><tr><td>A*30</td><td>3.7</td><td>E*01:01:01:01</td><td>47.5</td></tr><tr><td>A*31</td><td>2.6</td><td>E*01:03:02:01</td><td>42.1</td></tr><tr><td>A*26</td><td>2.6</td><td>E*01:03:01:01</td><td>5.0</td></tr><tr><td>A*23</td><td>2.4</td><td>E*01:01:01:03</td><td>2.5</td></tr><tr><td>A*68</td><td>2.1</td><td>E*01:03:05</td><td>2.0</td></tr><tr><td>A*33</td><td>1.8</td><td>E*01:06</td><td>0.5</td></tr><tr><td>A*32</td><td>1.8</td><td>E*01:03:02:02</td><td>0.5</td></tr><tr><td>A*25</td><td>1.3</td><td>H*01:01:01</td><td>33.9</td></tr><tr><td>A*80</td><td>0.3</td><td>H*02:04</td><td>15.9</td></tr><tr><td>HLA-G UTR1</td><td>31.3</td><td>H*02:01:01</td><td>15.1</td></tr><tr><td>HLA-G UTR2</td><td>21.8</td><td>H*del</td><td>13.9</td></tr><tr><td>HLA-G UTR4</td><td>14.0</td><td>H*02:04new</td><td>8.7</td></tr><tr><td>HLA-G UTR3</td><td>11.7</td><td>H*02:02</td><td>6.0</td></tr><tr><td>HLA-G UTR6</td><td>7.5</td><td>H*01:02</td><td>3.8</td></tr><tr><td>HLA-G UTR7</td><td>7.1</td><td>H*02:03</td><td>2.7</td></tr><tr><td>HLA-G UTR5</td><td>4.8</td><td>F*01:01:01</td><td>46.2</td></tr><tr><td>Blank</td><td>1.6</td><td>F*01:01:03</td><td>21.8</td></tr><tr><td>HLA-G UTR8</td><td>0.3</td><td>F*01:03:01</td><td>18.6</td></tr><tr><td></td><td></td><td>F*01:01:02</td><td>13.5</td></tr></tbody></table>

doi:10.1371/journal.pone.0163570.t002

6-digit HLA-F alleles were found, F*01:01:01 being the most frequent while HLA-F*01:02 and F*01:04 were not observed.

## Two loci linkage disequilibrium and haplotypes

Analysis of two-loci Global Linkage Disequilibrium (GLD) (Table 3) showed that HLA-A, HLA-H, HLA-F, HLA-G and HLA-G UTRs were all in highly significant pairwise GLD (p<0.001). In contrast, the HLA-E locus was not in significant GLD with any of them. However, when defined at a 4-digit resolution level, HLA-E displayed a significant GLD with HLA-F but the p-value of 1% was not significant after adjusting for multiple testing (Table 3). Haplotypes were estimated between each pair of loci and their combined frequencies are given in Table 4. Strikingly, a strong LD was found between HLA-A, HLA-G UTRs and HLA-H, as all haplotypes with a frequency above 3% and concerning these pairs of loci were in significant LD. These loci were also found to be highly associated with HLA-F, as most haplotypes with a frequency above 3% were also in LD.

The detailed description of the haplotypes observed with a frequency greater than 3% (or 2% for those involving HLA-A; Table 5) reveals significant allelic associations between HLA-A, HLA-G UTRs and HLA-H, respectively, with very high standardized residual values. HLA-F also displays significant allelic associations with HLA-A, HLA-G UTRs and HLA-H, but to a lower extent. In contrast, HLA-G alleles display few significant associations with HLA-A,
<table><caption>Table 3. Analysis of two-loci global linkage disequilibrium (GLD) between HLA-E (defined at 8 (8d) or 4-digits (4d)) and HLA-A, -H, -G, HLA-G UTRs and -F.</caption><thead><tr><th>HLA-E 8d</th><th>HLA-E 4d</th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td>0.81</td><td>0.39</td><td>HLA-A</td><td></td><td></td><td></td><td></td></tr><tr><td>0.34</td><td>0.08</td><td>&lt;0.001</td><td>HLA-H</td><td></td><td></td><td></td></tr><tr><td>0.99</td><td>0.85</td><td>&lt;0.001</td><td>&lt;0.001</td><td>HLA-G</td><td></td><td></td></tr><tr><td>0.62</td><td>0.25</td><td>&lt;0.001</td><td>&lt;0.001</td><td>&lt;0.001</td><td>HLA-G UTRs</td><td></td></tr><tr><td>0.09</td><td>0.01 *</td><td>&lt;0.001</td><td>&lt;0.001</td><td>&lt;0.001</td><td>&lt;0.001</td><td>HLA-F</td></tr></tbody></table>

doi:10.1371/journal.pone.0163570.t003

HLA-H, HLA-F and HLA-G UTRs, with the exception of HLA-G*01:03, *01:04 and *01:06. Finally, HLA-E reveals some significant associations, but the standardized residuals are very close to the threshold value for significance.

The tight association found between HLA-A, HLA-G UTRs and HLA-H is further supported by the three-loci haplotype frequency estimation shown in Table 6. HLA-G UTR3 is associated with HLA-H*del and HLA-A*23/24 alleles, HLA-G UTR6 is associated with HLA-H*02:02 and HLA-A*29, and HLA-G UTR7 is associated with HLA-H*02:04new and HLA-A*11, whereas HLA-G UTR2 displays no exclusive association.

Discussion

This study was devoted to the analysis of allelic associations between the HLA class Ia locus HLA-A, and the HLA class Ib loci HLA-E, -H, -G, HLA-G UTRs and -F which are all physically located within a region of 724 kb on chromosome 6. HLA-A seems to be involved in transplantation outcome in a different way from HLA-B or HLA-DRB1 [1, 3]. This singular feature could be due to the HLA-A molecule *per se* and/or be a consequence of yet unknown patterns of linkage disequilibrium between HLA-A and nearby class Ib genes.

Currently, interactions between HLA class Ib molecules are barely known, and even less so in the clinical setting. Genetic phylogeny is one approach that might shed some light on this subject, as it has been suggested that 1) HLA-A and -H share a more recent ancestor than

<table><caption>Table 4. Number (N) of haplotypes observed at frequencies (fq) greater or equal to 3% and 1% and their combined frequencies (Fc) in 191 voluntary blood donors from southeastern France.</caption><thead><tr><th>Pairs of loci</th><th>N of haplotypes with a <i>fq</i>&ge;1% (and Fc)</th><th>N of haplotypes with a <i>fq</i>&ge;3% (and Fc)</th><th>N of haplotypes with a <i>fq</i>&ge;3% in LD</th></tr></thead><tbody><tr><td>HLA-A~HLA-G UTR</td><td>18 (94.9)</td><td>6 (74.8)</td><td>6</td></tr><tr><td>HLA-A~HLA-H</td><td>19 (94.9)</td><td>6 (73.3)</td><td>6</td></tr><tr><td>HLA-A~HLA-F</td><td>22 (97.1)</td><td>9 (76.1)</td><td>6</td></tr><tr><td>HLA-A~HLA-G</td><td>19 (97.4)</td><td>7 (77.1)</td><td>2</td></tr><tr><td>HLA-A~HLA-E</td><td>27 (96.4)</td><td>9 (71.4)</td><td>3</td></tr><tr><td>HLA-H~HLA-G UTR</td><td>13 (93.8)</td><td>8 (85.3)</td><td>8</td></tr><tr><td>HLA-F~HLA-G UTR</td><td>13 (98.3)</td><td>9 (89.2)</td><td>6</td></tr><tr><td>HLA-E~HLA-G UTR</td><td>18 (96.5)</td><td>10 (82.4)</td><td>4</td></tr><tr><td>HLA-G~HLA-G UTR</td><td>10 (98.5)</td><td>8 (95.3)</td><td>4</td></tr><tr><td>HLA-H~HLA-F</td><td>18 (97.2)</td><td>8 (79.1)</td><td>5</td></tr><tr><td>HLA-H~HLA-G</td><td>12 (98.3)</td><td>9 (91.4)</td><td>3</td></tr><tr><td>HLA-E~HLA-H</td><td>16 (95.1)</td><td>8 (80.1)</td><td>2</td></tr><tr><td>HLA-E~HLA-F</td><td>12 (98.1)</td><td>7 (90.3)</td><td>2</td></tr><tr><td>HLA-E~HLA-G</td><td>10 (94.5)</td><td>6 (86.3)</td><td>0</td></tr></tbody></table>

doi:10.1371/journal.pone.0163570.t004

<footnote>* non-significant after adjusting for multiple testing.</footnote>
<table><caption>Table 5. (Continued)</caption><tr><td>A*2-E*01:03:02:01</td><td>15.2</td><td>12.1</td><td>1.18</td><td>E*01:01:01:01-UTR2</td><td>16.9</td><td>11.0</td><td>2.34</td><td>E*01:03:02:01-H*01:01:01</td><td>17.3</td><td>13.0</td><td>1.59</td><td>F*01:01:01-E*01:03:02:01</td><td>25.8</td><td>19.8</td><td>1.49</td></tr><tr><td>A*2-E*01:01:01:01</td><td>11.1</td><td>13.5</td><td>-0.87</td><td>E*01:03:02:01-UTR1</td><td>14.1</td><td>12.2</td><td>0.69</td><td>E*01:01:01:01-H*02:01:01</td><td>13.5</td><td>7.4</td><td>3.03</td><td>F*01:01:03-E*01:01:01</td><td>17.0</td><td>10.3</td><td>2.47</td></tr><tr><td>A*3-E*01:03:02:01</td><td>10.4</td><td>6.6</td><td>2.05</td><td>E*01:01:01:01-UTR1</td><td>12.3</td><td>13.4</td><td>-0.41</td><td>E*01:01:01:01-H*01:01:01</td><td>10.6</td><td>14.6</td><td>-1.40</td><td>F*01:01:01-E*01:01:01</td><td>16.4</td><td>21.9</td><td>-1.31</td></tr><tr><td>A*1-E*01:01:01:01</td><td>8.9</td><td>4.5</td><td>2.87</td><td>E*01:03:02:01-UTR4</td><td>11.2</td><td>6.1</td><td>2.76</td><td>E*01:03:02:01-H*02:04</td><td>10.3</td><td>8.5</td><td>0.83</td><td>F*01:03:01-E*01:03:02:01</td><td>13.7</td><td>8.0</td><td>2.42</td></tr><tr><td>A*29-E*01:03:02:01</td><td>8.5</td><td>3.6</td><td>3.57</td><td>E*01:03:02:01-UTR6</td><td>7.7</td><td>3.9</td><td>2.60</td><td>E*01:01:01:01-H*02:04</td><td>9.4</td><td>9.6</td><td>-0.09</td><td>F*01:01:02-E*01:01:01</td><td>9.8</td><td>6.4</td><td>1.64</td></tr><tr><td>A*3-E*01:01:01:01</td><td>5.1</td><td>7.4</td><td>-1.15</td><td>E*01:01:01:01-UTR3</td><td>5.4</td><td>5.3</td><td>0.09</td><td>E*01:03:02:01-H*02:02</td><td>8.0</td><td>3.5</td><td>3.37</td><td>F*01:03:01-E*01:01:01</td><td>4.2</td><td>8.8</td><td>-1.84</td></tr><tr><td>A*11-E*01:01:01:01</td><td>4.6</td><td>3.3</td><td>0.97</td><td>E*01:01:01:01-UTR7</td><td>4.4</td><td>3.1</td><td>1.01</td><td>E*01:01:01:01-H*del</td><td>7.3</td><td>5.3</td><td>1.17</td><td>F*01:01:01-E*01:03:01:01</td><td>3.4</td><td>2.1</td><td>1.14</td></tr><tr><td>A*24-E*01:01:01:01</td><td>4.5</td><td>4.3</td><td>0.18</td><td>E*01:03:02:01-UTR2</td><td>4.0</td><td>10.1</td><td>-2.53</td><td>E*01:01:01-H*02:04new</td><td>4.3</td><td>3.5</td><td>0.60</td><td></td><td></td><td></td><td></td></tr><tr><td>A*33-E*01:01:01:01</td><td>3.0</td><td>1.4</td><td>1.85</td><td>E*01:01:01:01-UTR5</td><td>3.3</td><td>3.1</td><td>0.17</td><td></td><td></td><td></td><td></td><td>Loci</td><td></td><td></td><td></td></tr><tr><td>A*24-E*01:03:02:01</td><td>2.2</td><td>3.8</td><td>-1.13</td><td>E*01:01:01:01-UTR4</td><td>3.1</td><td>6.7</td><td>-1.87</td><td></td><td></td><td></td><td></td><td>F*01:01:01-E*0103</td><td>32.4</td><td>23.4</td><td>2.05</td></tr><tr><td>A*23-E*01:01:01:01</td><td>2.0</td><td>1.4</td><td>0.68</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>F*01:01:03-E*0101</td><td>18.6</td><td>10.6</td><td>2.90</td></tr><tr><td>A*31-E*01:03:02:01</td><td>2.0</td><td>1.9</td><td>0.09</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>F*01:01:01-E*0101</td><td>13.7</td><td>22.5</td><td>-2.03</td></tr><tr><td>A*32-E*01:01:01:01</td><td>2.0</td><td>1.0</td><td>1.52</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>F*01:03:01-E*0103</td><td>13.2</td><td>9.4</td><td>1.48</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>F*01:01:02-E*0101</td><td>11.0</td><td>6.6</td><td>2.10</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>F*01:03:01-E*0101</td><td>5.4</td><td>9.1</td><td>-1.47</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>F*01:01:03-E*0103</td><td>3.2</td><td>11.0</td><td>-2.79</td></tr></table>

doi:10.1371/journal.pone.0163570.t005
<table><caption>Table 6. HLA-A~HLA-H~HLA-G UTRs haplotypes with a frequency greater than 1%.</caption><thead><tr><td>Haplotypes</td><td>Frequency</td></tr></thead><tbody><tr><td>A*01~H*02:01:01~HLA-G UTR2</td><td>10.21</td></tr><tr><td>A*02~H*01:01:01~HLA-G UTR 1</td><td>27.53</td></tr><tr><td>A*02~H*01:01:01~HLA-G UTR 5</td><td>1.39</td></tr><tr><td>A*02~H*01:01:01~HLA-G UTR 2</td><td>1.12</td></tr><tr><td>A*03~H*02:04~HLA-G UTR 4</td><td>11.6</td></tr><tr><td>A*03~H*02:04new~HLA-G UTR 4</td><td>1.67</td></tr><tr><td>A*11~H*02:04new~HLA-G UTR 7</td><td>6.92</td></tr><tr><td>A*23~H*null~HLA-G UTR 3</td><td>1.11</td></tr><tr><td>A*23~H*del~HLA-G UTR 3</td><td>1.11</td></tr><tr><td>A*24~H*del~HLA-G UTR3</td><td>8.78</td></tr><tr><td>A*25~H*01:02~HLA-G UTR2</td><td>1.11</td></tr><tr><td>A*26~H*01:02~HLA-G UTR2</td><td>2.22</td></tr><tr><td>A*29~H*02:02~HLA-G UTR6</td><td>5.83</td></tr><tr><td>A*30~H*01:01:01~HLA-G UTR1</td><td>1.67</td></tr><tr><td>A*30~H*02:01:01~HLA-G UTR2</td><td>1.39</td></tr><tr><td>A*31~H*02:04~HLA-G UTR5</td><td>1.39</td></tr><tr><td>A*33~H*02:04~HLA-G UTR5</td><td>1.94</td></tr><tr><td>A*68~H*02:01:01~HLA-G UTR2</td><td>1.65</td></tr></tbody></table>

doi:10.1371/journal.pone.0163570.t006

HLA-A with either HLA-B or C and that 2) HLA-E is more closely related to HLA-H than to -F or -G [46]. Moreover, HLA-E alleles are thought to pre-date most of the HLA-A and HLA-B polymorphisms [36]. However, an evolutionary perspective on class Ib polymorphisms would probably not be very informative regarding clinical implications.

Another possible approach to tackle this question was adopted here by looking at patterns of association previously described to have a clinical impact. Furthermore we have added to the interest of analyzing HLA-G UTRs in clinical studies as some of them display very specific associations, such as HLA-G UTR-4 associated with HLA-A*03, HLA-H*02:04 and HLA-F*01:03; HLA-G UTR-6 associated with HLA-A*29, HLA-H*02:02 and HLA-F*01:01; and HLA-G UTR-7 associated with HLA-A*11, HLA-H*02:04new and HLA-F*01:01. In turn, HLA-G UTR-2 certainly needs further investigation as it was observed in association with several HLA-A, -H and -F alleles, probably consisting of at least 2 clusters with respect to HLA-H and HLA-G.

Of interest, several studies have described a deleterious association of HLA-G*01:06~UTR2 and HLA-G*01:04~UTR3, in lung transplantation and pregnancy [24, 26–29]. In this study we confirmed the exclusive association of HLA-G*01:04 with HLA-G UTR3, HLA-A*23/*24 and the deletion encompassing HLA-H. HLA-E and HLA-F did not display an exclusive association with HLA-G*01:04, as this allele was associated with both F*01:01:03 (6%) and F*01:03:01 (3.3%), and with E*01:01:01 (5.4%) and E*01:03 (combined frequencies: 5.7%). Thus, we did not observe an exclusive association between HLA A*23/*24 and E*01:01 [36] could not be confirmed here. Whether the worse LTx progression reported for the HLA-G*01:04 allele is due to HLA-A*23/24 or to the deletion (encompassing HCG4P5, HLA-U, HLA-K, HCG4B, DDX39BP, MCCD1P1, HLA-T, HLA-H, HCG4P7, MICF and HCG3-2) remains to be explored. Noteworthy, this deletion, always associated to HLA-A*23 and A*24, was also found in HLA-G*01:01 carriers (4.4%). Interestingly, the HLA-H locus, translated into a truncated protein, has a signal sequence similar to that of HLA-A, except that the second position (Val>Leu)
<table><caption>**Table 7.** HLA-H, HLA-A and HLA-G signal sequence alignment. Variable amino-acid positions among loci are highlighted in bold.</caption><thead><tr><td>Locus</td><td>Signal sequence</td></tr></thead><tbody><tr><td>HLA-A</td><td>MAVMAPRTLLLLLSGALALTQTWA</td></tr><tr><td>HLA-G</td><td>MVVMAPRTLFLLLSGALTLTETWA</td></tr><tr><td>HLA-H</td><td>MVLMAPRTLLLLLSGALALTQTWA</td></tr></tbody></table>

doi:10.1371/journal.pone.0163570.t007

compared to other HLA class I proteins (Table 7) [2]. Interestingly, the amino-acid in the 10th position is identical to HLA-A (Leu) and different from HLA-G (Phe). This characteristic of the HLA-G peptide signal is put forward to explain a higher affinity of the HLA-G-derived nonamer-HLA-E complex with CD94/NKG2C receptor complex [66]. Whether the HLA-H signal sequence is loaded by HLA-E and preferentially protects cells from NK lysis remains to be explored, but the absence of HLA-H could potentially have an impact on the interaction of HLA-E with its receptors.

Concerning HLA-G*01:06, it displayed a restricted association pattern as it was associated with HLA-G UTR-2, HLA-A*01, HLA-H*02:01:01 and HLA-F*01:01:03, and in positive LD with HLA-E*01:01:01 and negative LD with HLA-E*01:03:02. Thus, the negative clinical effect of HLA-G*01:06 could be the combined result of lower immune-tolerance due to reduced HLA-G expression expected from the association with HLA-G UTR-2 [23] and to reduced HLA-E expression expected from the association with HLA-E*01:01 [35]; and yet unknown features of HLA-H*02:01:01 and HLA-F*01:01:03. Confirming such a hypothesis demands a much larger cohort, because of the low allelic frequency of HLA-G*01:06 (3.9%), and would also greatly benefit from functional studies on HLA-H and HLA-F.

Looking now at the peculiar pattern of linkage disequilibrium revealed by HLA-E, its large physical distance from other HLA class Ib loci (>660Kb), as well as balancing selective forces acting on this locus could be put forward to explain the absence of GLD. However, the functional and clinical implications of HLA-E alleles not being associated to other HLA class Ib loci remain to be formally investigated.

In an HLA class Ia matching transplant context, however, the significant association observed between E*01:01 and A*01 on the one hand, and between E*01:03 and A*29, on the other hand, as shown in [36], implies that no HLA-E heterozygosity should be expected for most HLA-A*01 and A*29 carriers.

Finally, besides protein modifications and differential levels of expression that could impact their respective interaction and affinity with their specific receptors, non-classical HLA polymorphisms could also influence antibody production. Non-native forms of HLA-Ib, caused by alternative splicing or loss of β-2 microglobulin, may expose immunogenic cryptic epitopes and trigger an immune response [67]. Naturally occurring anti-HLA-E antibodies were shown to mimic anti-HLA-Ia antibodies and such cross-reactive binding was reported to affect the screening of anti-HLA-Ia antibodies in renal and liver transplant recipients [68]. Anti-HLA-F IgG, associated with the inflammatory status of systemic lupus erythematosus, may represent a biomarker of primary importance in transplantation, although their exact function has not yet been fully elucidated [69].

In conclusion, this study offers a unique view of non-classical HLA class I allelic associations that may be helpful for future functional studies in deciphering HLA class Ib protein interactions. Since HLA-A typing of recipients and donors is systematically and routinely performed in immunogenetics laboratories, these results could also potentially improve transplant matching strategies.
**Supporting Information**

S1 Table. Allelic conversion table for HLA-H SNPs.
(DOC)

**Author Contributions**

**Conceptualization:** JFE CR CP JDC.

**Data curation:** JFE CP.

**Formal analysis:** FC VF AT CR.

**Funding acquisition:** JFE JC.

**Investigation:** SB JDC.

**Methodology:** JFE CR CP JDC.

**Project administration:** JFE JC CP.

**Resources:** JC JFE.

**Software:** VF SB CR JDC.

**Supervision:** JDC.

**Validation:** SB CR JDC.

**Writing – original draft:** SB JDC.

**Writing – review & editing:** FC VF SB AT JFE CR JC CP JDC.

**References**

1. Bucin D, Ekberg H, Lindholm A, Persson H, Tufvesson G, Albrechtsen D. Enhanced renal graft survival in HLA-A mismatched transplants. Transplantation proceedings. 1992; 24(6):2473-4. PMID: 1465837.

2. Robinson J, Halliwell JA, McWilliam H, Lopez R, Parham P, Marsh SGE. The IMGT/HLA database 10.1093/nar/gks949. Nucleic Acids Research. 2013; 41(D1):D1222-D7. doi: 10.1093/nar/gks949 PMID: 23080122

3. Roelen DL, Stobbe I, Young NT, van Bree SP, Doxiadis II, Oudshoorn M, et al. Permissible and immunogenic HLA-A mismatches: cytotoxic T-cell precursor frequencies reflect graft survival data. Hum Immunol. 2001; 62(7):661-7. doi: 10.1016/s0198-8859(01)00263-4 PMID: 11423171.

4. Kamoun M, Holmes JH, Israni AK, Kearns JD, Teal V, Yang WP, et al. HLA-A amino acid polymorphism and delayed kidney allograft function. Proc Natl Acad Sci U S A. 2008; 105(48):18883-8. doi: 10.1073/pnas.0810308105 PMID: 19033208; PubMed Central PMCID: PMC2596272.

5. Balan V, Ruppert K, Demetris AJ, Ledneva T, Duquesnoy RJ, Detre KM, et al. Long-term outcome of human leukocyte antigen mismatching in liver transplantation: results of the National Institute of Diabetes and Digestive and Kidney Diseases Liver Transplantation Database. Hepatology. 2008; 48(3):878-88. doi: 10.1002/hep.22435 PMID: 18752327.

6. Hayes D Jr., Black SM, Tobias JD, Kopp BT, Kirkby SE, Mansour HM, et al. Influence of human leukocyte antigen mismatching on bronchiolitis obliterans syndrome in lung transplantation. J Heart Lung Transplant. 2015. doi: 10.1016/j.healun.2015.08.022 PMID: 26496787.

7. Makhadiyeva D, Lam L, Moatari M, Vallance J, Zheng Y, Campbell EC, et al. MHC class I dimer formation by alteration of the cellular redox environment and induction of apoptosis. Immunology. 2012; 135(2):133-9. doi: 10.1111/j.1365-2567.2011.03518.x PMID: 22044191; PubMed Central PMCID: PMC3277715.

8. Kochan G, Escors D, Breckpot K, Guerrero-Setas D. Role of non-classical MHC class I molecules in cancer immunosuppression. Oncoimmunology. 2013; 2(11):e26491. doi: 10.4161/onci.26491 PMID: 24482746; PubMed Central PMCID: PMC3894240.
9. Foroni I, Couto AR, Bettencourt BF, Santos M, Lima M, Bruges-Armas J. HLA-E, HLA-F and HLA-G—The Non-Classical Side of the MHC Cluster. 2014. In: HLA and Associated Important Diseases [Internet]. InTech. Distinguished Prof. Xi Yongzhi (Ed.). Available from: http://www.intechopen.com/books/hla-and-associated-important-diseases/hla-e-hla-f-and-hla-g-the-non-classical-side-of-the-mhc-cluster.

10. Geraghty DE, Koller BH, Orr HT. A human major histocompatibility complex class I gene that encodes a protein with a shortened cytoplasmic segment. Proc Natl Acad Sci U S A. 1987; 84(24):9145-9. doi: 10.1073/pnas.84.24.9145 PMID: 3480534; PubMed Central PMCID: PMC299709.

11. Comiskey M, Domino KE, Warner CM. HLA-G is found in lipid rafts and can act as a signaling molecule. Hum Immunol. 2007; 68(1):1-11. doi: 10.1016/j.humimm.2006.09.005 PMID: 17207707; PubMed Central PMCID: PMC1839863.

12. Rebmann V, da Silva Nardi F, Wagner B, Horn PA. HLA-G as a Tolerogenic Molecule in Transplantation and Pregnancy. J Immunol Res. 2014; 2014:297073. Epub 2014/08/22. doi: 10.1155/2014/297073 PMID: 25143957; PubMed Central PMCID: PMC4131093.

13. Pabon MA, Navarro CE, Osorio JC, Gomez N, Moreno JP, Donado AF, et al. Impact of human leukocyte antigen molecules e, f, and g on the outcome of transplantation. Transplantation proceedings. 2014; 46(9):2957-65. doi: 10.1016/j.transproceed.2014.07.010 PMID: 25420801.

14. Donadi EA, Castelli EC, Arnaiz-Villena A, Roger M, Rey D, Moreau P. Implications of the polymorphism of HLA-G on its function, regulation, evolution and disease association. Cell Mol Life Sci. 2011; 68(3):369-95. doi: 10.1007/s00018-010-0580-7 PMID: 21107637; PubMed Central PMCID: PMC3021195.

15. Rizzo R, Trentini A, Bortolotti D, Manfrinato MC, Rotola A, Castellazzi M, et al. Matrix metalloproteinase-2 (MMP-2) generates soluble HLA-G1 by cell surface proteolytic shedding. Mol Cell Biochem. 2013; 381(1-2):243-55. PMID: 23737137. doi: 10.1007/s11010-013-1708-5

16. Riteau B, Rouas-Freiss N, Menier C, Paul P, Dausset J, Carosella ED. HLA-G2, -G3, and -G4 isoforms expressed as nonmature cell surface glycoproteins inhibit NK and antigen-specific CTL cytolysis. J Immunol. 2001; 166(8):5018-26. Epub 2001/04/06. doi: 10.4049/jimmunol.166.8.5018 PMID: 11290782.

17. Zhang WQ, Xu DP, Liu D, Li YY, Ruan YY, Lin A, et al. HLA-G1 and HLA-G5 isoforms have an additive effect on NK cytolysis. Hum Immunol. 2014; 75(2):182-9. doi: 10.1016/j.humimm.2013.11.001 PMID: 24269702.

18. Zhao L, Teklemariam T, Hantash BM. Mutated HLA-G3 localizes to the cell surface but does not inhibit cytotoxicity of natural killer cells. Cellular immunology. 2014; 287(1):23-6. Epub 2013/12/21. doi: 10.1016/j.cellimm.2013.11.005 PMID: 24355712.

19. Boyson JE, Erskine R, Whitman MC, Chiu M, Lau JM, Koopman LA, et al. Disulfide bond-mediated dimerization of HLA-G on the cell surface. Proc Natl Acad Sci U S A. 2002; 99(25):16180-5. doi: 10.1073/pnas.212643199 PMID: 12454284; PubMed Central PMCID: PMC138585.

20. Kuroki K, Maenaka K. Immune modulation of HLA-G dimer in maternal-fetal interface. Eur J Immunol. 2007; 37(7):1727-9. doi: 10.1002/eji.200737515 PMID: 17587197.

21. Rebmann V, van der Ven K, Passler M, Pfeiffer K, Krebs D, Grosse-Wilde H. Association of soluble HLA-G plasma levels with HLA-G alleles. Tissue Antigens. 2001; 57(1):15-21. Epub 2001/02/13. doi: 10.1034/j.1399-0039.2001.057001015.x PMID: 11169254.

22. Di Cristofaro J, El Moujally D, Agnel A, Mazieres S, Cortey M, Basire A, et al. HLA-G haplotype structure shows good conservation between different populations and good correlation with high, normal and low soluble HLA-G expression. Hum Immunol. 2013; 74(2):203-6. PMID: 23127550. doi: 10.1016/j.humimm.2012.10.027

23. Carlini F, Traore K, Cherouat N, Roubertoux P, Buhler S, Cortey M, et al. HLA-G UTR haplotype conservation in the Malian population: association with soluble HLA-G. PLoS One. 2013; 8(12):e82517. doi: 10.1371/journal.pone.0082517 PMID: 24376542; PubMed Central PMCID: PMC3871591.

24. Di Cristofaro J, Reynaud-Gaubert M, Carlini F, Roubertoux P, Loundou A, Basire A, et al. HLA-G*01:04 UTR3 Recipient Correlates With Lower Survival and Higher Frequency of Chronic Rejection After Lung Transplantation. American journal of transplantation: official journal of the American Society of Transplantation and the American Society of Transplant Surgeons. 2015. doi: 10.1111/ajt.13305 PMID: 25989360.

25. Castelli EC, Mendes-Junior CT, Veiga-Castelli LC, Roger M, Moreau P, Donadi EA. A comprehensive study of polymorphic sites along the HLA-G gene: implication for gene regulation and evolution. Mol Biol Evol. 2011; 28(11):3069-86. doi: 10.1093/molbev/msr138 PMID: 21622995

26. Moreau P, Contu L, Alba F, Lai S, Simoes R, Orru S, et al. HLA-G gene polymorphism in human placentas: possible association of G*0106 allele. Biol Reprod. 2008; 79(3):459-67. Epub 2008/05/30. doi: 10.1095/biolreprod.108.068874 PMID: 18509163.
27. Tan CY, Ho JF, Chong YS, Loganath A, Chan YH, Ravichandran J, et al. Paternal contribution of HLA-G*0106 significantly increases risk for pre-eclampsia in multigravid pregnancies. Mol Hum Reprod. 2008; 14(5):317-24. doi: 10.1093/molehr/gan013 PMID: 18353802.

28. Aldrich CL, Stephenson MD, Karrison T, Odem RR, Branch DW, Scott JR, et al. HLA-G genotypes and pregnancy outcome in couples with unexplained recurrent miscarriage. Mol Hum Reprod. 2001; 7 (12):1167-72. Epub 2001/11/24. doi: 10.1093/molehr/7.12.1167 PMID: 11719594.

29. Vargas RG, Sarturi PR, Mattar SB, Bompeixe EP, Silva Jdos S, Pirri A, et al. Association of HLA-G alleles and 3' UTR 14 bp haplotypes with recurrent miscarriage in Brazilian couples. Hum Immunol. 2011; 72(6):479-85. doi: 10.1016/j.humimm.2011.02.011 PMID: 21354454.

30. Misra MK, Pandey SK, Kapoor R, Sharma RK, Kapoor R, Prakash S, et al. HLA-G gene expression influenced at allelic level in association with end stage renal disease and acute allograft rejection. Hum Immunol. 2014; 75(8):833-9. doi: 10.1016/j.humimm.2014.06.005 PMID: 24929145.

31. Allan DS, Lepin EJ, Braud VM, O'Callaghan CA, McMichael AJ. Tetrameric complexes of HLA-E, HLA-F, and HLA-G. J Immunol Methods. 2002; 268(1):43-50. doi: 10.1016/s0022-1759(02)00199-0 PMID: 12213342.

32. Pratheek BM, Nayak TK, Sahoo SS, Mohanty PK, Chattopadhyay S, Chakraborty NG, et al. Mammalian non-classical major histocompatibility complex I and its receptors: Important contexts of gene, evolution, and immunity. Indian journal of human genetics. 2014; 20(2):129-41. doi: 10.4103/0971-6866. 142855 PMID: 25400340; PubMed Central PMCID: PMC4228563.

33. Heinrichs H, Orr HT. HLA non-A,B,C class I genes: their structure and expression. Immunol Res. 1990; 9(4):265-74. doi: 10.1007/bf02935526 PMID: 2089074.

34. Walpole NG, Kjer-Nielsen L, Kostenko L, McCluskey J, Brooks AG, Rossjohn J, et al. The structure and stability of the monomorphic HLA-G are influenced by the nature of the bound peptide. Journal of molecular biology. 2010; 397(2):467-80. doi: 10.1016/j.jmb.2010.01.052 PMID: 20122941; PubMed Central PMCID: PMC2905647.

35. Kraemer T, Blasczyk R, Bade-Doeding C. HLA-E: a novel player for histocompatibility. J Immunol Res. 2014; 2014:352160. doi: 10.1155/2014/352160 PMID: 25401109; PubMed Central PMCID: PMC4221882.

36. Geraghty DE, Stockschleader M, Ishitani A, Hansen JA. Polymorphism at the HLA-E locus predates most HLA-A and -B polymorphism. Hum Immunol. 1992; 33(3):174-84. doi: 10.1016/0198-8859(92) 90069-y PMID: 1618657.

37. Grimsley C, Ober C. Population genetic studies of HLA-E: evidence for selection. Hum Immunol. 1997; 52(1):33-40. doi: 10.1016/S0198-8859(96)00241-8 PMID: 9021407.

38. Veiga-Castelli LC, Castelli EC, Mendes CT Jr., da Silva WA Jr., Faucher MC, Beauchemin K, et al. Non-classical HLA-E gene variability in Brazilians: a nearly invariable locus surrounded by the most variable genes in the human genome. Tissue Antigens. 2012; 79(1):15-24. doi: 10.1111/j.1399-0039. 2011.01801.x PMID: 22050290.

39. Tamouza R, Busson M, Fortier C, Diagne I, Diallo D, Sloma I, et al. HLA-E*0101 allele in homozygous state favors severe bacterial infections in sickle cell anemia. Hum Immunol. 2007; 68(10):849-53. doi: 10.1016/j.humimm.2007.08.260 PMID: 17961774.

40. Boyle LH, Gillingham AK, Munro S, Trowsdale J. Selective export of HLA-F by its cytoplasmic tail. J Immunol. 2006; 176(11):6464-72. doi: 10.4049/jimmunol.176.11.6464 PMID: 16709803.

41. Goodridge JP, Burian A, Lee N, Geraghty DE. HLA-F complex without peptide binds to MHC class I protein in the open conformer form. J Immunol. 2010; 184(11):6199-208. doi: 10.4049/jimmunol. 1000078 PMID: 20483783; PubMed Central PMCID: PMC3777411.

42. Lee N, Ishitani A, Geraghty DE. HLA-F is a surface marker on activated lymphocytes. Eur J Immunol. 2010; 40(8):2308-18. doi: 10.1002/eji.201040348 PMID: 20865824; PubMed Central PMCID: PMC3867582.

43. Lepin EJ, Bastin JM, Allan DS, Roncador G, Braud VM, Mason DY, et al. Functional characterization of HLA-F and binding of HLA-F tetramers to ILT2 and ILT4 receptors. Eur J Immunol. 2000; 30 (12):3552-61. doi: 10.1002/1521-4141(200012)30:12&#60;3552::AID-IMMU3552&#60;3.0.CO;2-L PMID: 11169396.

44. Goodridge JP, Burian A, Lee N, Geraghty DE. HLA-F and MHC class I open conformers are ligands for NK cell Ig-like receptors. J Immunol. 2013; 191(7):3553-62. doi: 10.4049/jimmunol.1300081 PMID: 24018270; PubMed Central PMCID: PMC3780715.

45. Goodridge JP, Lee N, Burian A, Pyo CW, Tykodi SS, Warren EH, et al. HLA-F and MHC-I open conformers cooperate in a MHC-I antigen cross-presentation pathway. J Immunol. 2013; 191(4):1567-77. doi: 10.4049/jimmunol.1300080 PMID: 23851683; PubMed Central PMCID: PMC3732835.
46. Zemmour J, Koller BH, Ennis PD, Geraghty DE, Lawlor DA, Orr HT, et al. HLA-AR, an inactivated anti-gen-presenting locus related to HLA-A. Implications for the evolution of the MHC. J Immunol. 1990; 144(9):3619-29. PMID: 2329283.

47. Feder JN, G nirke A, Thomas W, Tsuchihashi Z, Ruddy DA, Basava A, et al. A novel MHC class I-like gene is mutated in patients with hereditary haemochromatosis. Nature genetics. 1996; 13(4):399-408. doi: 10.1038/ng0896-399 PMID: 8696333.

48. Pascolo S, Ginhoux F, Laham N, Walter S, Schoor O, Probst J, et al. The non-classical HLA class I molecule HFE does not influence the NK-like activity contained in fresh human PBMCs and does not interact with NK cells. Int Immunol. 2005; 17(2):117-22. doi: 10.1093/intimm/dxh191 PMID: 15623549.

49. Grimsley C, Mather KA, Ober C. HLA-H: a pseudogene with increased variation due to balancing selection at neighboring loci. Mol Biol Evol. 1998; 15(12):1581-8. doi: 10.1093/oxfordjournals.molbev. a025886 PMID: 9866194.

50. Schmidt CM, Orr HT. A physical linkage map of HLA-A, -G, -7.5p, and -F. Hum Immunol. 31. United States 1991. p. 180-5. doi: 10.1016/0198-8859(91)90024-4 PMID: 1890019

51. Shukla H, Gillespie GA, Srivastava R, Collins F, Chorney MJ. A class I jumping clone places the HLA-G gene approximately 100 kilobases from HLA-H within the HLA-A subregion of the human MHC. Genomics. 1991; 10(4):905-14. doi: 10.1016/0888-7543(91)90178-h PMID: 1916822.

52. Morandi F, Pistoia V. Interactions between HLA-G and HLA-E in Physiological and Pathological Conditions. Frontiers in immunology. 2014; 5:394. doi: 10.3389/fimmu.2014.00394 PMID: 25202308; PubMed Central PMCID: PMC4141331.

53. Di Cristofaro J, Buhler S, Frassati C, Basire A, Galicher V, Baier C, et al. Linkage disequilibrium between HLA-G*0104 and HLA-E*0103 alleles in Tswa Pygmies. Tissue Antigens. 2011; 77(3):193-200. doi: 10.1111/j.1399-0039.2010.01599.x PMID: 21299523

54. Ober C, Rosinsky B, Grimsley C, van der Ven K, Robertson A, Runge A. Population genetic studies of HLA-G: allele frequencies and linkage disequilibrium with HLA-A1. J Reprod Immunol. 32. Ireland 1996. p. 111-23. doi: 10.1016/s0165-0378(96)01000-5 PMID: 9023816

55. Shao W, Joo Yoo Y, Tang J, Zhang K, Allen S, Kaslow RA. Tight linkage disequilibrium between HLA-G and HLA-A alleles in native africans http://dx.doi.org/10.1016/j.humimm.2006.08.180. Human Immunology. 2006; 67, Supplement 1(0):S118. doi: 10.1016/j.humimm.2006.08.180

56. Kolte AM, Steffensen R, Nielsen HS, Hviid TV, Christiansen OB. Study of the structure and impact of human leukocyte antigen (HLA)-G-A, HLA-G-B, and HLA-G-DRB1 haplotypes in families with recurrent miscarriage. Hum Immunol. 2010; 71(5):482-8. doi: 10.1016/j.humimm.2010.02.001 PMID: 20149831.

57. el Kahloun A, Vernet C, Jouanolle AM, Boretto J, Mauvieux V, Le Gall JY, et al. A continuous restriction map from HLA-E to HLA-F. Structural comparison between different HLA-A haplotypes. Immunogenetics. 1992; 35(3):183-9. Epub 1992/01/01. doi: 10.1007/bf00185112 PMID: 1537609.

58. Geraghty DE, Koller BH, Hansen JA, Orr HT. The HLA class I gene family includes at least six genes and twelve pseudogenes and gene fragments. J Immunol. 1992; 149(6):1934-46. Epub 1992/09/15. PMID: 1517563.

59. Geraghty DE, Pei J, Lipsky B, Hansen JA, Taillon-Miller P, Bronson SK, et al. Cloning and physical mapping of the HLA class I region spanning the HLA-E-to-HLA-F interval by using yeast artificial chromosomes. Proc Natl Acad Sci U S A. 1992; 89(7):2669-73. Epub 1992/04/01. doi: 10.1073/pnas.89.7. 2669 PMID: 1557372; PubMed Central PMCID: PMCPMC48723.

60. Nunes JM. Using uniformat and gene[rate] to analyse data with ambiguities in population genetics. Figshare. 2014; http://dx.doi.org/10.6084/m9 figshare.984299.

61. Nunes JM, Buhler S, Roessli D, Sanchez-Mazas A, collaboration HL-n. The HLA-net GENE[RATE] pipeline for effective HLA data analysis and its application to 145 population samples from Europe and neighbouring areas. Tissue Antigens. 2014; 83(5):307-23. doi: 10.1111/tan.12356 PMID: 24738646.

62. Nunes JM. Using uniformat and gene[rate] to Analyze Data with Ambiguities in Population Genetics. Evol Bioinform Online. 2015; 11(Suppl 2):19-26. doi: 10.4137/EBO.S32415 PMID: 26917942; PubMed Central PMCID: PMC4762493.

63. Nunes JM, Riccio ME, Buhler S, Di D, Currat M, Ries F, et al. Analysis of the HLA population data (AHPD) submitted to the 15th International Histocompatibility/Immunogenetics Workshop by using the Gene[rate] computer tools accommodating ambiguous data (AHPD project report). Tissue Antigens. 2010; 76(1):18-30. doi: 10.1111/j.1399-0039.2010.01469.x PMID: 20331842

64. Agresti A. An Introduction to Categorical Data Analysis. 2nd edn. John Wiley & Sons NY, editor 2007.
65. Maier S, Grzeschik M, Weiss EH, Ulbrecht M. Implications of HLA-E allele expression and different HLA-E ligand diversity for the regulation of NK cells. Hum Immunol. 2000; 61(11):1059-65. doi: 10.1016/s0198-8859(00)00190-7 PMID: 11137208.

66. Llano M, Lee N, Navarro F, Garcia P, Albar JP, Geraghty DE, et al. HLA-E-bound peptides influence recognition by inhibitory and triggering CD94/NKG2 receptors: preferential response to an HLA-G-derived nonamer. Eur J Immunol. 1998; 28(9):2854-63. doi: 10.1002/(SICI)1521-4141(199809)28:09&#60;2854::AID-IMMU2854&#62;3.0.CO;2-W PMID: 9754572.

67. Jucaud V RM, Terasaki PI Immunobiology of HLA Class-Ib Molecules in Transplantation. SOJ Immunol. 2015; 3(4):1-15. http://dx.doi.org/10.15226/2372-0948/3/4/00137.

68. Ravindranath MH, Terasaki PI, Pham T, Jucaud V, Kawakita S. Therapeutic preparations of IVIg contain naturally occurring anti-HLA-E antibodies that react with HLA-Ia (HLA-A/-B/-Cw) alleles. Blood. 2013; 121(11):2013-28. doi: 10.1182/blood-2012-08-447771 PMID: 23305735.

69. Jucaud V, Ravindranath MH, Terasaki PI, Morales-Buenrostro LE, Hiepe F, Rose T, et al. Serum antibodies to human leucocyte antigen (HLA)-E, HLA-F and HLA-G in patients with systemic lupus erythematosus (SLE) during disease flares: Clinical relevance of HLA-F autoantibodies. Clinical and experimental immunology. 2016; 183(3):326-40. doi: 10.1111/cei.12724 PMID: 26440212; PubMed Central PMCID: PMC4750595.