HLA-G(∗)01:05N null allele in Mayans (Guatemala) and Uros (Titikaka Lake, Peru): evolution and population genetics

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 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.

Next generation sequencing for HLA loci in full heritage Pima Indians of Arizona, Part II: HLA-A, -B, and -C with selected non-classical loci at 4-field resolution from whole genome sequences

While the samples and data from the Pima Indians of the Gila River Indian Community have been included in many international HLA workshops and conferences and have been the focus of numerous population reports and the source of novel alleles at the classical HLA loci, they have not been studied for the non-classical loci. In order to expand our HLA-disease association studies, we typed over 300 whole genome sequences from full Pima heritage members, controlled for first degree relationship, and employed recently developed computer algorithms to resolve HLA alleles. Both classical-HLA-A, -B, and -C- and non-classical- HLA-E, -F, -G, -J, -L, -W, -Y, -DPA2, -DPB2, -DMA, -DMB, -DOA, -DRB2, -DRB9, TAP1- loci were typed at the 4-field level of resolution. We present allele and selected haplotype frequencies, test the genotype distributions for population structure, discuss the issues that are created for tests of Hardy-Weinberg equilibrium over the four sample spaces of high resolution HLA typing, and address the implications for the evolution of non-classical pseudogenes that are no longer expressed in a phenotype subject to natural selection.

HLA-G in Mayas from Yucatan: An evolutionary approach

HLA-G allele frequencies were studied in Yucatán (Mexico) Maya Amerindians by a direct exon DNA sequencing technique. It is described that Mayas are probably one of the first populations together with Olmecs that populated Meso America and that important HLA genetic differences between Mexican and Guatemalan Mayas support that Maya languages were imposed to several neighbouring Amerindian groups. HLA-G*01:01:02, HLA-G*01:01:01 and HLA-G*01:04:01 are the most frequent alleles in this population. It is remarkable that HLA-G*01:05N allele was not found in the population in accordance with similar results found in another Amerindians. Also, protein allele HLA-G*01:04 frequency is found not to differ to those found in another far or close living Amerindians in contrast to other World populations. It seems that while high HLA-G*01:05N frequency is found in Iran and Middle East populations, probably where this allele appeared within an ancestral HLA-A*19 group of alleles haplotype and it is maintained by unknown evolutionary forces, Amerindians do not have a high frequency because a founder effect or because required natural evolutionary forces do not exist in America. Finally, we believe useful to study HLA-G evolution for its physiopathology understanding in addition to the many papers on statistics on HLA-G and in vitro models that are yearly published.

HLA-G: Function, polymorphisms and pathology

HLA-G immune modulatory genes and molecules are presently being studied by a widespread number of research groups. In the present study, we do not aim to be exhaustive since the number of manuscripts published every year is overwhelming. Instead, our aim is pointing out facts about HLA-G function, polymorphism and pathology that have been confirmed by several different researchers, together with exposing aspects that may have been overlooked or not sufficiently remarked in this productive field of study. On the other hand, we question whether performing mainly studies on HLA-G and disease associations is going to give a clear answer in the future, since 40 years of study of classical HLA molecules association with disease has still given no definite answer on this issue.

HLA-G in Amerindians: Epidemiology and Worldwide Population Comparison

Background:

HLA-G molecules are immunosuppressive and avoid fetal rejection by giving negative signals to maternal immune system from fetal trophoblast cell surface. HLA-G genes have been associated to different pathologies: Spontaneous abortions, autoimmunity, tumor progression, transplant rejection and infection. In addition, different World populations show remarkable different HLA-G allele frequencies in the allele that does not produce a full HLA-G molecule (HLA-G*05N); this allele is almost absent in studied Amerindians.

Objectives:

The aim is to study HLA-A.-B,-DRB1 and –G alleles and extended haplotypes in Amerindians for the first time. This may be useful to asses HLA-G epidemiology, association to disease and Preventive Medicine in Amerindians.

Methods:

HLA-A,-B and -DRB1 have been typed by using standard automatic protocols. HLA-G alleles have been detected by direct HLA-G exon 2, exon 3 and exon 4 DNA sequencing. Computer calculations have been done by specific standard methods.

Results:

HLA-A,-B,-DRB1 and –G extended haplotypes have been calculated in Amerindians for the first time. Also, their HLA-G frequencies have been compared with worldwide populations.

Conclusion:

Low frequencies of null HLA-G*01:05N allele are found in Amerindians. The extended haplotypes with this allele bear other typical Amerindian HLA-DRB1 alleles and its origin is discussed. HLA-G allele frequency profile is closer to that of Europeans than to that of Far East Asians. Our findings are useful to Preventive Medicine and Epidemiology associated to Fertility and HLA-G associated pathology and transplantation.

HLA-G, -A haplotypes in Amerindians (Ecuador): HLA-G*01:05N World distribution

HLA-G and HLA-A frequencies have been analysed in Amerindians from Ecuador. HLA-G allele frequencies are found to be closer to those of other Amerindians (Mayas from Guatemala and Uros from Peru) and closer to European ones than to Far East Asians groups, particularly, regarding to HLA-G*01:04 allele. HLA-G/-A haplotypes have been calculated for the first time in Amerindians. It is remarkable that HLA-G*01:05N "null" allele is found in a very low frequency (like in Amerindian Mayas and Uros) and is also found in haplotypes belonging to the HLA-A19 group of alleles (HLA-A*30, -A*31, -A*33). It was previously postulated that HLA-G*01:05N appeared in HLA-A*30/-B*13 haplotypes in Middle East Mediterraneans. It may be hypothesized that in Evolution, HLA-G*01:05N existed primarily in one of the HLA extant or extinct -A19 haplotype, whether this haplotype was placed in Middle East or other World areas, including America. However, the highest present day HLA-G*01:05N frequencies are found in Middle East Mediterraneans.

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.