Genetic diversity of the HLA-G coding region in Amerindian populations from the Brazilian Amazon: a possible role of natural selection

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.

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 Polymorphisms Associated with HIV Infection and Preeclampsia in South Africans of African Ancestry

OBJECTIVES

HLA-G, part of the major histocompatibility complex (MHC), is associated with the risk of developing preeclampsia (PE). In this study, we determined the contribution of specific HLA-G polymorphisms on the risk of developing preeclampsia in HIV-infected and uninfected South Africans of African ancestry.

METHODS

One hundred and ninety-three women of African ancestry were enrolled (74 HIV-uninfected normotensive, 60 HIV-infected normotensive, 34 HIV-uninfected, and 25 HIV-infected preeclamptics). Sanger sequencing of the untranslated region was performed to genotype six SNPs, i.e., 14 bp Ins/Del of rs66554220, rs1710, rs1063320, rs1610696, rs9380142, and rs1707).

RESULTS

For rs66554220, we have the following results: (a) based on pregnancy type-the Ins/Ins and Del/Ins genotype frequency was higher in preeclampsia (PE) compared to normotensive pregnancies (Ins/Ins vs. Del/Ins, P = 0.02∗: OR (95%CI) = 13.44 (0.7222-249.9); Del/Del vs. Del/Ins, P = 0.03∗: OR (95%CI) = 2.95 (1.10-7.920)); (b) based on HIV status-the Ins/Ins showed both genotypic and allelic association with HIV infection. HIV-infected PE has higher Ins/Ins genotypic and allelic frequencies compared to HIV-uninfected PE (Ins/Ins vs. Del/Ins, P = 0.005∗∗: OR (95%CI) = 21.32 (1.71-4.17); Ins, P = 0.005∗∗; OR (95%IC) = 21.32 (1.71-4.17)). For rs1707, we have the following results: (a) based on pregnancy type-there were CT genotypic frequencies in PE, more especially LOPE compared to normotensive pregnancies (TT vs. CT, P = 0.0092∗∗: OR (95%CI) = 5.(1.39 - 25.64)), and no allelic association was noted; (b) based on HIV status-CT was higher in HIV-infected LOPE compared to uninfected LOPE (TT vs. TC, P = 0.0006∗∗∗: OR (95%CI) = 40.00 (2.89 - 555.1)). For rs1710 and rs1063320, no significant differences in the genotype and allele frequencies were noted based on pregnancy type and HIV status. For rs9380142, we have the following results: (a) based on pregnancy type-no significant differences were noted between normotensive compared to PE pregnancies; (b) based on HIV status-AA genotypes occurred more in the HIV-infected PE group (AA vs. GG, P = 0.02∗: OR (95%CI) = 13.97 (0.73 - 269.4)), while A allelic frequency occurred more in HIV-infected PE, especially LOPE compared to uninfected groups (A vs. G, P = 0.0003∗∗∗: OR (95%CI) = 10.72 (2.380 - 48.32); P = 0.02∗: OR (95%CI) = 9.00 (1.07 - 75.74)). For rs1610696, we have the following results: (a) based on pregnancy type-genotypic and allelic frequencies of CC were higher in PE compared to normotensive pregnancies (CC vs. GG, P = 0.0003∗∗∗: OR (95%CI) = 31.87 (1.861 - 545.9); C, P = 0.0001∗∗∗: OR (95%IC) = 21.91 (2.84 - 169.0)); (b) based on HIV status-GG frequencies were higher in the HIV-infected PE more especially LOPE groups (GG vs. GC, P = 0.02∗: OR (95%CI) = 16.87 (0.81 - 352.1); GG vs. CC, P = 0.0001∗∗∗: OR (95%CI) = 159.5 (13.10 - 1942)).

CONCLUSION

Selected HLA-G 14 bp polymorphisms (Ins/Ins) and genotypic and allelic differences in rs9380142, rs1610696, and rs1707 are associated with the pathogenesis of preeclampsia in HIV-infected South African women of African ancestry. More genetic studies evaluating the association between preeclampsia and HIV infection are needed to improve diagnosis and antenatal care.

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.

The role of HLA-G in parasitic diseases

Little attention has been devoted to the role of HLA-G gene and molecule on parasitic disorders, and the available studies have focused on malaria, African and American trypanosomiasis, leishmaniosis, toxoplasmosis and echinococcosis. After reporting a brief description regarding the role of the cells of innate and adaptive immune system against parasites, we reviewed the major features of the HLA-G gene and molecule and the role of HLA-G on the major cells of immune system. Increased levels of soluble HLA-G (sHLA-G) have been observed in patients presenting toxoplasmosis and in the active phase of echinococcosis. In addition, increased sHLA-G has also been associated with increased susceptibility to malaria and increased susceptibility to develop human African trypanosomiasis (HAT). In contrast, decreased membrane-bound HLA-G has been reported in placenta of patients infected with Plasmodium falciparum and in heart and colon of patients presenting Chagas disease. The 3' untranslated region of the HLA-G gene has been the main focus of studies on malaria, HAT and Chagas disease, exhibiting distinct patterns of associations. Considering that HLA-G is an immune checkpoint molecule, inhibiting the activity of several cells of the immune system, the excessive neoexpression and the increased sHLA-G levels together with the decreased constitutive tissue expression of membrane-bound HLA-G may be detrimental to the host infected with parasite agents.

HLA-F coding and regulatory segments variability determined by massively parallel sequencing procedures in a Brazilian population sample

Human Leucocyte Antigen F (HLA-F) is a non-classical HLA class I gene distinguished from its classical counterparts by low allelic polymorphism and distinctive expression patterns. Its exact function remains unknown. It is believed that HLA-F has tolerogenic and immune modulatory properties. Currently, there is little information regarding the HLA-F allelic variation among human populations and the available studies have evaluated only a fraction of the HLA-F gene segment and/or have searched for known alleles only. Here we present a strategy to evaluate the complete HLA-F variability including its 5' upstream, coding and 3' downstream segments by using massively parallel sequencing procedures. HLA-F variability was surveyed on 196 individuals from the Brazilian Southeast. The results indicate that the HLA-F gene is indeed conserved at the protein level, where thirty coding haplotypes or coding alleles were detected, encoding only four different HLA-F full-length protein molecules. Moreover, a same protein molecule is encoded by 82.45% of all coding alleles detected in this Brazilian population sample. However, the HLA-F nucleotide and haplotype variability is much higher than our current knowledge both in Brazilians and considering the 1000 Genomes Project data. This protein conservation is probably a consequence of the key role of HLA-F in the immune system physiology.

Mitochondrial DNA variability among six South American Amerindian villages from the Pano linguistic group

Although scattered throughout a large geographic area, the members of the Pano linguistic group present strong ethnic, linguistic, and cultural homogeneity, a feature that causes them to be considered components of a same "Pano" tribe. Nevertheless, the genetic homogeneity between Pano villages has not yet been examined. To study the genetic structure of the Pano linguistic group, four major Native American mitochondrial DNA (mtDNA) founder haplogroups were analyzed in 77 Amerindians from six villages of four Pano tribes (Katukina, Kaxináwa, Marúbo, and Yaminawa) located in the Brazilian Amazon. The central position of these tribes in the continent makes them relevant for attempts to reconstruct population movements in South America. Except for a single individual that presented an African haplogroup L, all remaining individuals presented one of the four Native American haplogroups. Significant heterogeneity was observed across the six Pano villages. Although Amerindian populations are usually characterized by considerable interpopulational diversity, the high heterogeneity level observed is unexpected if the strong ethnic, linguistic, and cultural homogeneity of the Pano linguistic group is taken into account. The present findings indicate that the ethnic, linguistic, and cultural homogeneity does not imply genetic homogeneity. Even though the genetic heterogeneity uncovered may be a female-specific process, the most probable explanation for that is the joint action of isolation and genetic drift as major factors influencing the genetic structure of the Pano linguistic group.

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 14 bp deletion/insertion polymorphism and mother-to-child transmission of HIV

The human leukocyte antigen HLA-G, highly expressed at the maternal-fetal interface, has a pivotal role in mediating immune tolerance. In this study we investigated the influence of HLA-G 14 bp insertion polymorphism in human immunodeficiency virus (HIV)-1 mother-to-child HIV-1 transmission. The 14 bp insertion polymorphism was analyzed among 99 HIV-1 positive mothers and 329 infants born to HIV-positive mothers in Zambia, among whom vertical transmission status and timing had been determined. HLA-G 14 bp insertion polymorphism was detected using a custom TaqMan single nucleotide polymorphisms (SNPs) genotyping assay. Logistic regression was conducted to examine the associations between HLA-G alleles and the risk of HIV transmission. The 14 bp insertion allele was more frequent in HIV exposed-uninfected (EU) infants than in infected infants, and was associated with reduced risk of both in utero (IU) and intrapartum (IP) HIV transmission, after adjusting for maternal cluster of differentiation 4 (CD4) cell count and plasma viral load. Maternal HLA-G 14 bp insertion genotype and HLA-G concordance between mother and child were not associated with the risk of perinatal HIV transmission. The presence of the 14 bp insertion associates with protection toward IU and IP HIV infection in children from Zambia, suggesting that HLA-G could be involved in the vertical transmission of HIV.