Evolution and molecular interactions of major histocompatibility complex (MHC)-G, -E and -F genes

Classical HLA (Human Leukocyte Antigen) is the Major Histocompatibility Complex (MHC) in man. HLA genes and disease association has been studied at least since 1967 and no firm pathogenic mechanisms have been established yet. HLA-G immune modulation gene (and also -E and -F) are starting the same arduous way: statistics and allele association are the trending subjects with the same few results obtained by HLA classical genes, i.e., no pathogenesis may be discovered after many years of a great amount of researchers’ effort. Thus, we believe that it is necessary to follow different research methodologies: (1) to approach this problem, based on how evolution has worked maintaining together a cluster of immune-related genes (the MHC) in a relatively short chromosome area since amniotes to human at least, i.e., immune regulatory genes (MHC-G, -E and -F), adaptive immune classical class I and II genes, non-adaptive immune genes like (C2, C4 and Bf) (2); in addition to using new in vitro models which explain pathogenetics of HLA and disease associations. In fact, this evolution may be quite reliably studied during about 40 million years by analyzing the evolution of MHC-G, -E, -F, and their receptors (KIR—killer-cell immunoglobulin-like receptor, NKG2—natural killer group 2-, or TCR-T-cell receptor—among others) in the primate evolutionary lineage, where orthology of these molecules is apparently established, although cladistic studies show that MHC-G and MHC-B genes are the ancestral class I genes, and that New World apes MHC-G is paralogous and not orthologous to all other apes and man MHC-G genes. In the present review, we outline past and possible future research topics: co-evolution of adaptive MHC classical (class I and II), non-adaptive (i.e., complement) and modulation (i.e., non-classical class I) immune genes may imply that the study of full or part of MHC haplotypes involving several loci/alleles instead of single alleles is important for uncovering HLA and disease pathogenesis. It would mainly apply to starting research on HLA-G extended haplotypes and disease association and not only using single HLA-G genetic markers.

HLA-G genetic diversity and evolutive aspects in worldwide populations

HLA-G is a promiscuous immune checkpoint molecule. The HLA-G gene presents substantial nucleotide variability in its regulatory regions. However, it encodes a limited number of proteins compared to classical HLA class I genes. We characterized the HLA-G genetic variability in 4640 individuals from 88 different population samples across the globe by using a state-of-the-art method to characterize polymorphisms and haplotypes from high-coverage next-generation sequencing data. We also provide insights regarding the HLA-G genetic diversity and a resource for future studies evaluating HLA-G polymorphisms in different populations and association studies. Despite the great haplotype variability, we demonstrated that: (1) most of the HLA-G polymorphisms are in introns and regulatory sequences, and these are the sites with evidence of balancing selection, (2) linkage disequilibrium is high throughout the gene, extending up to HLA-A, (3) there are few proteins frequently observed in worldwide populations, with lack of variation in residues associated with major HLA-G biological properties (dimer formation, interaction with leukocyte receptors). These observations corroborate the role of HLA-G as an immune checkpoint molecule rather than as an antigen-presenting molecule. Understanding HLA-G variability across populations is relevant for disease association and functional studies.

HLA-G 3'UTR haplotype frequencies in highland and lowland South Native American populations

Human Leukocyte Antigen (HLA)-G participates in several biological processes, including reproduction, vascular remodeling, immune tolerance, and hypoxia response. HLA-G is a potential candidate gene for high altitude adaptation since its expression is modulated in both micro and macro environment under hypoxia and constant cellular stress. Besides the promoter region, the HLA-G 3'untranslated region (UTR) influences HLA-G expression patterns through several post-transcriptional mechanisms. Currently, the 3'UTR genetic diversity in terms of altitude adaptation of Native American populations is still unexplored, particularly at high altitude ecoregions. Here, we evaluated 288 Native Americans from 9 communities located in the Andes [highland (HL); ≥2,500 m (range = 2,838-4,433 m)] and 8 populations located in lowland (LL) regions [<2,500 m (range = 80-431 m); Amazonian tropical forest, Brazilian central plateau, and Chaco] of South America. In total, nine polymorphic sites and ten haplotypes were observed. The most frequent haplotypes (UTR-1, UTR-2, and UTR-3) accounted for ∼ 77% of haplotypes found in LL, while in the HL, the same haplotypes reach ∼ 93%. Also, a remarkable high frequency of putative ancestral UTR-5 haplotype was observed in LL (21.5%), while in HL UTR-2 reaches up to 47%. Further, UTR-2 frequency positively correlates with altitude-related variables, while a negative correlation for UTR-5 was observed. From an evolutionary perspective, we observed a tendency towards balancing selection in HL and LL populations thus suggesting that haplotypes of ancient and more derived alleles may have been co-opted for relatively recent adaptations such as those experienced by modern humans in the highland and lowland of South America. We also discuss how long-term balancing selection can be a reservoir of genetic variants that can be positively selected. Finally, our study provides some pieces of evidence that HLA-G 3'UTR haplotypes may have contributed to high altitude adaptation in the Andes.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

HLA-G Polymorphisms Are Associated with Non-segmental Vitiligo among Brazilians

(1) Background: Vitiligo is characterized by white patches on the skin caused by loss of melanocyte activity or the absence of these cells. The available treatments minimize the symptoms by retarding the process of skin depigmentation or re-pigmenting the affected regions. New studies are required for a better comprehension of the mechanisms that trigger the disease and for the development of more efficient treatments. Studies have suggested an autoimmune feature for vitiligo, based on the occurrence of other autoimmune diseases in vitiligo patients and their relatives, and on the involvement of genes related to the immune response. (2) Methods: We evaluated, by massive parallel sequencing, polymorphisms of the HLA-G gene in vitiligo patients and control samples, to verify if variants of this gene could influence the susceptibility to vitiligo. (3) Results: We detected an association with non-segmental vitiligo regarding the haplotype Distal-010101a/G*01:01:01:01/UTR-1, adjusting for population stratification by using ancestry-informative markers (AIMs). (4) Conclusions: It remains unclear whether the HLA-G variants associated with vitiligo were detected because of the high linkage disequilibrium (LD) with HLA-A*02, or if the HLA-A variants previously reported as associated with vitiligo were detected because of the high LD with HLA-G*01:01:01:01/UTR-1, or if both genes jointly contribute to vitiligo susceptibility.

MHC class I polymorphic Alu insertion (POALIN) allele and haplotype frequencies in the Arabs of the United Arab Emirates and other world populations

Polymorphic Alu insertions (POALINs) are found throughout the human genome and have been used in various studies to infer geographic origin of human populations. The main aim of this study was to determine the allele and haplotype frequencies of five POALINs, AluHF, AluHG, AluHJ, AluTF and AluMICB, within the major histocompatibility complex (MHC) class I region of 95 UAE Arabs, and correlate their frequencies to those of the HLA-A, HLA-C and HLA-B class I allele lineages. Evolutionary relationships between the POALINs of the Arabs and those previously studied in populations of African, Asian and European descent were compared. At each of the five Alu loci (AluHF, AluHG, AluHJ, AluTF and AluMICB), Alu insertion was designated as Alu(locus)*02 and absence was Alu(locus)*01. The AluHG insertion (AluHG*02) had the highest frequency (0.332), followed by AluHF*02 (0.300), AluHJ*02 (0.263), AluMICB*02 (0.111) and AluTF*02 (0.058). Of the 270 Alu-HLA haplotypes pairs in the UAE Arabs, 110 had no Alu insertion, and 54 had an Alu insertion at >50% per haplotype. An Alu insertion >75% per haplotype was found between AluMICB*02 and HLA-B*14, HLA-B*22, HLA-B*44, HLA-B*55, HLA-B*57 and HLA-B*73, and with HLA-C*01 and HLA-C*18; AluHJ*02 with HLA-A*01, HLA-A*19, HLA-A*24 and HLA-A*32; AluHG*02 with HLA-A*02 and HLA-B*18; and AluHF*02 with HLA-A*10. The genotyped allele and haplotype frequencies of the MHC POALINs in UAE Arabs were compared with the results of 30 previously published Asian, European, American and African populations. Phylogenetic and multidimensional scaling (MDS) analysis of the relative MHC POALINs allele and haplotype frequencies revealed that the UAE Arabs have a similar lineage to Caucasians and the most distant genetic relationship to the Waorani native American population of Ecuador. The structure of both the phylogenetic tree and the MDS analysis supports the Out of Africa theory of human evolution. The nature of the clusters suggests the Arabian Middle East represents a crossroads from which human populations migrated towards Asia in the east and Europe to the north-west.

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.

Adaptation of human skin color in various populations

BACKGROUND

Skin color is a well-recognized adaptive trait and has been studied extensively in humans. Understanding the genetic basis of adaptation of skin color in various populations has many implications in human evolution and medicine.

DISCUSSION

Impressive progress has been made recently to identify genes associated with skin color variation in a wide range of geographical and temporal populations. In this review, we discuss what is currently known about the genetics of skin color variation. We enumerated several cases of skin color adaptation in global modern humans and archaic hominins, and illustrated why, when, and how skin color adaptation occurred in different populations. Finally, we provided a summary of the candidate loci associated with pigmentation, which could be a valuable reference for further evolutionary and medical studies.

CONCLUSION

Previous studies generally indicated a complex genetic mechanism underlying the skin color variation, expanding our understanding of the role of population demographic history and natural selection in shaping genetic and phenotypic diversity in humans. Future work is needed to dissect the genetic architecture of skin color adaptation in numerous ethnic minority groups around the world, which remains relatively obscure compared with that of major continental groups, and to unravel the exact genetic basis of skin color adaptation.

Haplotypes of the HLA-G 3' Untranslated Region Respond to Endogenous Factors of HLA-G+ and HLA-G- Cell Lines Differentially

The immune checkpoint HLA-G prevents maternal rejection of the fetus and contributes in cancer invasion and acceptance of allografts. The 5’ and 3’ regulatory regions of the HLA-G gene are polymorphic and balancing selection probably maintains this variability. It is proposed that nucleotide variations may affect the level of HLA-G expression. To investigate this issue we aimed to analyze how haplotypes of the 3’ untranslated region (3’UTR) with highest worldwide frequencies, namely UTR-1, UTR-2, UTR-3, UTR-4, UTR-5, UTR-18 and UTR-7, impact the expression of a luciferase reporter gene in vitro. Experiments performed with the HLA-G positive cell lines JEG-3 (choricarcinoma) and FON (melanoma), and with the HLA-G negative cell lines M8 (melanoma) and U251MG (glioblastoma) showed that the HLA-G 3’UTR polymorphism influences the response to endogenous cellular factors and may vary according to the cell type. UTR-5 and UTR-7 impact the activity of luciferase the most whereas UTR-2, UTR-3, UTR-4, and UTR-18 have intermediate impact, and UTR-1 has the lowest impact. These results corroborate the previous associations between amounts of plasma sHLA-G levels and 3’UTR haplotypes in healthy individuals and reinforce that 3’UTR typing may be a predictor of the genetic predisposition of an individual to express different levels of HLA-G.

Immunogenetics of prostate cancer and benign hyperplasia--the potential use of an HLA-G variant as a tag SNP for prostate cancer risk

Human leukocyte antigen G (HLA-G) is an immunomodulatory molecule with important roles both physiologically as well as an escape mechanism of cancer cells. In this study, we evaluated the impact of eight polymorphisms at the 3' untranslated region (3'UTR) of the HLA-G gene in the development of prostate cancer (PCa) and benign prostatic hyperplasia (BPH). A total of 468 DNA samples of Brazilian men predominantly Euro-descendant with PCa (N = 187), BPH (N = 152) and healthy control individuals (N = 129) were evaluated. The HLA-G 3'UTR region was amplified by polymerase chain reaction (PCR), sequenced and genotyped to identify the 14 bp insertion/deletion (rs371194629), +3003T/C (rs1707), +3010C/G (rs1710), +3027A/C (rs17179101), +3035C/T (rs17179108), +3142G/C (rs1063320), +3187A/G (rs9380142) and +3196C/G (rs1610696) polymorphisms. Regression logistic and chi-square tests were performed to verify the influence of single nucleotide polymorphisms (SNPs) in PCa and/or BPH susceptibility, as well as in PCa progression (clinicopathological status). Our data showed the UTR-4 haplotype as a risk factor to PCa in comparison with control [odds ratio (OR) 2.35, 95% confidence interval (CI) 1.39-3.96, P adjusted = 0.003) and BPH groups (OR 1.82, 95% CI 1.15-2.86, P adjusted = 0.030). Further, the 'non-14bp Ins_ + 3142G_+3187A' haplotype (OR 1.56, 95% CI 1.10-2.20, P adjusted = 0.036), the +3003CT genotype (OR 4.44, 95% CI 1.33-4.50, P adjusted = 0.032) and the +3003C allele (OR 2.33, 95% CI 1.38-3.92, P adjusted = 0.016) also conferred susceptibility to PCa. Our data suggest an important influence of HLA-G 3'UTR polymorphisms in PCa susceptibility and support the use of the +3003 variant as a tag SNP for PCa risk.

HLA-E coding and 3' untranslated region variability determined by next-generation sequencing in two West-African population samples

HLA-E is a non-classical Human Leucocyte Antigen class I gene with immunomodulatory properties. Whereas HLA-E expression usually occurs at low levels, it is widely distributed amongst human tissues, has the ability to bind self and non-self antigens and to interact with NK cells and T lymphocytes, being important for immunosurveillance and also for fighting against infections. HLA-E is usually the most conserved locus among all class I genes. However, most of the previous studies evaluating HLA-E variability sequenced only a few exons or genotyped known polymorphisms. Here we report a strategy to evaluate HLA-E variability by next-generation sequencing (NGS) that might be used to other HLA loci and present the HLA-E haplotype diversity considering the segment encoding the entire HLA-E mRNA (including 5'UTR, introns and the 3'UTR) in two African population samples, Susu from Guinea-Conakry and Lobi from Burkina Faso. Our results indicate that (a) the HLA-E gene is indeed conserved, encoding mainly two different protein molecules; (b) Africans do present several unknown HLA-E alleles presenting synonymous mutations; (c) the HLA-E 3'UTR is quite polymorphic and (d) haplotypes in the HLA-E 3'UTR are in close association with HLA-E coding alleles. NGS has proved to be an important tool on data generation for future studies evaluating variability in non-classical MHC genes.

The Role of HLA-G Molecule and HLA-G Gene Polymorphisms in Tumors, Viral Hepatitis, and Parasitic Diseases

Considering that the non-classical HLA-G molecule has well-recognized tolerogenic properties, HLA-G expression is expected to be deleterious when present in tumor cells and in cells chronically infected by viruses, whereas HLA-G expression is expected to be advantageous in autoimmune disorders. The expression of HLA-G on tissue or peripheral blood cells, the levels of soluble HLA-G and polymorphic sites along the gene have been studied in several disorders. In this study, we revised the role of the molecule and polymorphic sites along the HLA-G gene in tumors, viral hepatitis, and parasitic disorders. Overall, several lines of evidence clearly show that the induction of HLA-G expression in tumors has been associated with worse disease outcome and disease spread. In addition, the few studies conducted on hepatitis and parasitic disorders indicate that HLA-G may contribute to disease pathogenesis. Few isolated polymorphic sites, primarily located at the coding or 3′ untranslated HLA-G region, have been evaluated in these disorders, and a complete HLA-G typing together with the study of gene regulatory elements may further help on the understanding of the influence of the genetic background on disease susceptibility.

Insights into HLA-G Genetics Provided by Worldwide Haplotype Diversity

Human leukocyte antigen G (HLA-G) belongs to the family of non-classical HLA class I genes, located within the major histocompatibility complex (MHC). HLA-G has been the target of most recent research regarding the function of class I non-classical genes. The main features that distinguish HLA-G from classical class I genes are (a) limited protein variability, (b) alternative splicing generating several membrane bound and soluble isoforms, (c) short cytoplasmic tail, (d) modulation of immune response (immune tolerance), and (e) restricted expression to certain tissues. In the present work, we describe the HLA-G gene structure and address the HLA-G variability and haplotype diversity among several populations around the world, considering each of its major segments [promoter, coding, and 3′ untranslated region (UTR)]. For this purpose, we developed a pipeline to reevaluate the 1000Genomes data and recover miscalled or missing genotypes and haplotypes. It became clear that the overall structure of the HLA-G molecule has been maintained during the evolutionary process and that most of the variation sites found in the HLA-G coding region are either coding synonymous or intronic mutations. In addition, only a few frequent and divergent extended haplotypes are found when the promoter, coding, and 3′UTRs are evaluated together. The divergence is particularly evident for the regulatory regions. The population comparisons confirmed that most of the HLA-G variability has originated before human dispersion from Africa and that the allele and haplotype frequencies have probably been shaped by strong selective pressures.

Transcriptional and posttranscriptional regulations of the HLA-G gene

HLA-G has a relevant role in immune response regulation. The overall structure of the HLA-G coding region has been maintained during the evolution process, in which most of its variable sites are synonymous mutations or coincide with introns, preserving major functional HLA-G properties. The HLA-G promoter region is different from the classical class I promoters, mainly because (i) it lacks regulatory responsive elements for IFN-γ and NF-κB, (ii) the proximal promoter region (within 200 bases from the first translated ATG) does not mediate transactivation by the principal HLA class I transactivation mechanisms, and (iii) the presence of identified alternative regulatory elements (heat shock, progesterone and hypoxia-responsive elements) and unidentified responsive elements for IL-10, glucocorticoids, and other transcription factors is evident. At least three variable sites in the 3' untranslated region have been studied that may influence HLA-G expression by modifying mRNA stability or microRNA binding sites, including the 14-base pair insertion/deletion, +3142C/G and +3187A/G polymorphisms. Other polymorphic sites have been described, but there are no functional studies on them. The HLA-G coding region polymorphisms might influence isoform production and at least two null alleles with premature stop codons have been described. We reviewed the structure of the HLA-G promoter region and its implication in transcriptional gene control, the structure of the HLA-G 3'UTR and the major actors of the posttranscriptional gene control, and, finally, the presence of regulatory elements in the coding region.

HLA-G +3142 polymorphism as a susceptibility marker in two rheumatoid arthritis populations in Brazil

In this study, we sought to investigate the genetic influence of two HLA-G 3'-untranslated region (3'-UTR) polymorphisms - 14 bp (rs66554220) and +3142C>G (rs1063320) and their compounding haplotypes in susceptibility to rheumatoid arthritis (RA) in a two-region Brazilian study comprising of 539 patients and 489 controls. All subjects were polymerase chain reaction (PCR) genotyped for the referred polymorphisms and logistic regression models controlling for sex, city and age were performed. Homozygozity for the +3142G allele was associated with an increased risk of RA [odds ratio (OR) = 1.45, 95% confidence interval (CI) = 1.075-1.959, P(Bonf)  = 0.030], whereas no association was observed for the 14 bp polymorphism. Haplotype comparisons between patients and controls showed a decreased frequency of the delC haplotype in patients (OR = 0.70, 95% CI = 0.521-0.946, P(Bonf)  = 0.040), which remained significant in the rheumatoid factor (RF)-positive group (OR = 0.66, 95% CI = 0.482-0.900, P(Bonf)  = 0.018), but not in the RF-negative group. These results corroborate the hypothesis of an involvement of HLA-G in the susceptibility of RA. The +3142G allele is associated with haplotype lineages that share high identity and are regarded as low producers. The presence of the G allele in homozygosis could be responsible for a low HLA-G expression profile that could favor the triggering of RA.

Worldwide HLA-E nucleotide and haplotype variability reveals a conserved gene for coding and 3' untranslated regions

The human leukocyte antigen-E (HLA-E) locus is a human major histocompatibility complex (MHC) gene associated with immune-modulation and suppression of the immune response by the interaction with specific natural killer (NK) and T cell receptors (TCRs). It is considered one of the most conserved genes of the human MHC; however, this low nucleotide variability seems to be a consequence of the scarce number of studies focusing on this subject. In this manuscript we assessed the nucleotide variability at the HLA-E coding and 3' untranslated regions (3'UTRs) in Brazil and in the populations from the 1000Genomes Consortium. Twenty-eight variable sites arranged into 33 haplotypes were detected and most of these haplotypes (98.2%) are encoding one of the two HLA-E molecules found worldwide, E*01:01 and E*01:03. Moreover, three worldwide spread haplotypes, associated with the coding alleles E*01:01:01, E*01:03:01 and E*01:03:02, account for 85% of all HLA-E haplotypes, suggesting that they arose early before human speciation. In addition, the low nucleotide diversity found for the HLA-E coding and 3'UTR in worldwide populations suggests that the HLA-E gene is in fact a conserved gene, which might be a consequence of its key role in the modulation of the immune system.