Balancing immunity and tolerance: genetic footprint of natural selection in the transcriptional regulatory region of HLA-G

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.

Role of the HLA-G immune checkpoint molecule in pregnancy

The non-classical HLA class I molecule HLA-G is expressed in trophoblasts where it contributes to maternal-fetal tolerance. HLA-G has been implicated in the control of trophoblast invasion, uterine vascular remodeling, and maintenance of a local immunosuppressive state. Understanding HLA-G biology at the maternal-fetal interface is therefore a critical issue in reproduction. In this regard, we review here: (i) the effects of HLA-G on decidual leucocytes and stromal cells, (ii) the contribution of trogocytosis in HLA-G expression on decidual cells, (iii) its interaction with the ILT2, ILT4 and KIR2DL4 receptors, (iv) the link between HLA-G polymorphism and pregnancy disorders, and (v) the expression of newly-described HLA-G isoforms at the maternal-fetal interface.

Fine scale human genetic structure in three regions of Cameroon reveals episodic diversifying selection

Inferences from genetic association studies rely largely on the definition and description of the underlying populations that highlight their genetic similarities and differences. The clustering of human populations into subgroups (population structure) can significantly confound disease associations. This study investigated the fine-scale genetic structure within Cameroon that may underlie disparities observed with Cameroonian ethnicities in malaria genome-wide association studies in sub-Saharan Africa. Genotype data of 1073 individuals from three regions and three ethnic groups in Cameroon were analyzed using measures of genetic proximity to ascertain fine-scale genetic structure. Model-based clustering revealed distinct ancestral proportions among the Bantu, Semi-Bantu and Foulbe ethnic groups, while haplotype-based coancestry estimation revealed possible longstanding and ongoing sympatric differentiation among individuals of the Foulbe ethnic group, and their Bantu and Semi-Bantu counterparts. A genome scan found strong selection signatures in the HLA gene region, confirming longstanding knowledge of natural selection on this genomic region in African populations following immense disease pressure. Signatures of selection were also observed in the HBB gene cluster, a genomic region known to be under strong balancing selection in sub-Saharan Africa due to its co-evolution with malaria. This study further supports the role of evolution in shaping genomes of Cameroonian populations and reveals fine-scale hierarchical structure among and within Cameroonian ethnicities that may impact genetic association studies in the country.

Review article: the genetics of the human leucocyte antigen region in inflammatory bowel disease

BACKGROUND

The human leucocyte antigen (HLA) complex, located at chromosome 6p21.3 is a highly polymorphic region containing the classical class I and II HLA genes. The region is highly associated with inflammatory bowel disease (IBD), largely through genome-wide association studies (GWAS).

AIMS

To review the role of HLA in immune function, summarise data on risk/protective HLA genotypes for IBD, discuss the role of HLA in IBD pathogenesis, treatment and examine limitations that might be addressed by future research.

METHODS

An organised search strategy was used to collate articles describing HLA genes in IBD, including Crohn's disease and ulcerative colitis.

RESULTS

All classical HLA genes with variation (including HLA-A, B, C, DRB1, DQA1, DQB1, DPA1 and DPB1) harbour IBD-associated genotypes. The most implicated gene is HLA-DRB1, with HLA-DRB1*03:01 the most associated risk allele in both Crohn's disease and ulcerative colitis. Elucidating precise disease associations is challenging due to high linkage disequilibrium between HLA genotypes. The mechanisms by which risk alleles cause disease are multifactorial, with the best evidence indicating structural and electrostatic alteration impacting antigen binding and downstream signalling. Adverse medication events have been associated with HLA genotypes including with thiopurines (pancreatitis) and anti-TNF agents (antibody formation).

CONCLUSIONS

The HLA complex is associated with multiple risk/protective alleles for IBD. Future research utilising long-read technology, ascertainment of zygosity and integration in disease modelling will improve the functional understanding and clinical translation of genetic findings.

HLA-A promoter, coding, and 3'UTR sequences in a Brazilian cohort, and their evolutionary aspects

HLA-A is the second most polymorphic locus of the human leucocyte antigen (HLA) complex encoding a key molecule for antigen presentation and NK cell modulation. Many studies have evaluated HLA-A variability in worldwide populations, focusing mainly on exons, but the regulatory segments have been poorly characterized. HLA-A variability is particularly high in the segment encoding the peptide-binding groove (exons 2 and 3), which is related to the antigen presentation function and the balancing selection in these segments. Here we evaluate the genetic diversity of the HLA-A gene considering a continuous segment encompassing the extended promoter (1.5 kb upstream of the first translated ATG), all exons and introns, and the entire 3' untranslated region, by using massively parallel sequencing. To achieve this goal, we used a freely available bioinformatics workflow that optimizes read mapping for HLA genes and defines complete sequences using either the phase among variable sites directly observed in sequencing data and probabilistic models. The HLA-A variability detected in a highly admixed population sample from Brazil shows that the HLA-A regulatory segments present few, but divergent sequences. The regulatory segments are in close association with the coding alleles. Both exons and introns are highly variable. Moreover, patterns of molecular diversity suggest that the promoter, in addition to the coding region, might be under the same selective pressure, but a different scenario arises when it comes to exon 4 and the 3'UTR segment.

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.

Population Genetics and Natural Selection in Rheumatic Disease

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Natural selection is an important influence on human genetic variation. Because immune and inflammatory function genes are enriched for signals of positive selection, the prevalence of rheumatic disease-risk alleles seen in different populations is partially the result of differing selective pressures (eg, due to pathogens). This review summarizes the genetic regions associated with susceptibility to different rheumatic diseases and concomitant evidence for natural selection, including known agents of selection exerting selective pressure in these regions.

The genetic diversity within the 1.4 kb HLA-G 5' upstream regulatory region moderately impacts on cellular microenvironment responses

The HLA-G 5'URR extending 1.4 kb from the ATG presents a unique set of regulatory elements among HLA genes. Several variable sites have been described that coincide with or are close to these elements, thus HLA-G 5'URR polymorphism might influence the HLA-G expression level. We cloned the ten most frequent HLA-G 5'URR haplotypes to evaluate their activity on a luciferase reporter gene in HLA-G⁺ cell lines (JEG-3/choriocarcinoma and FON⁺/melanoma). We also investigated associations between the plasma HLA-G (sHLA-G) levels and the HLA-G 5'URR variability in 157 healthy individuals. Cell lines were transfected with pGL3-Basic vector constructions containing HLA-G 5'URR sequences. The G010101a (in JEG-3) and G010101b (in FON⁺) haplotypes exhibited higher promoter activity, whereas the G010101d (in JEG-3) and G010102a (in FON⁺) haplotypes exhibited lower promoter activity. In the presence of HLA-G inducers (interferon-β and progesterone) or repressors (cyclopamine) HLA-G promoter activity was modulated, but certain haplotypes exhibited differential responses. No strict association was observed between plasma sHLA-G levels and the 5'URR haplotypes or genotypes; however, the G010101b haplotype was underrepresented among HLA-G-negative plasmas. Therefore, the HLA-G 5'URR polymorphism may have an impact on the modulation of HLA-G gene expression, but alone provides a limited predictive value for sHLA-G levels in vivo.

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.

Soluble human leukocyte antigen -G during pregnancy and infancy in Benin: Mother/child resemblance and association with the risk of malaria infection and low birth weight

Human leukocyte antigen (HLA) G is a tolerogenic molecule involved in the maternal-fetal immune tolerance phenomenon. Its expression during some infectious diseases leading to immune evasion has been established. A first study conducted in Benin has shown that the production of soluble HLA-G (sHLA-G) during the first months of life is strongly correlated with the maternal level at delivery and associated with low birth weight and malaria. However sHLA-G measurements during pregnancy were not available for mothers and furthermore, to date the evolution of sHLA-G in pregnancy is not documented in African populations. To extend these previous findings, between January 2010 and June 2013, 400 pregnant women of a malaria preventive trial and their newborns were followed up in Benin until the age of 2 years. Soluble HLA-G was measured 3 times during pregnancy and repeatedly during the 2 years follow-up to explore how sHLA-G evolved and the factors associated. During pregnancy, plasma levels of sHLA-G remained stable and increased significantly at delivery (p<0.001). Multigravid women seemed to have the highest levels (p = 0.039). In infants, the level was highest in cord blood and decreased before stabilizing after 18 months (p<0.001). For children, a high level of sHLA-G was associated with malaria infection during the follow-up (p = 0.02) and low birth weight (p = 0.06). The mean level of sHLA-G during infancy was strongly correlated with the mother’s level during pregnancy (<0.001), and not only at delivery. Moreover, mothers with placental malaria infection had a higher probability of giving birth to a child with a high level of sHLA-g (p = 0.006). High sHLA-G levels during pregnancy might be associated with immune tolerance related to placental malaria. Further studies are needed but this study provides a first insight concerning the potential role of sHLA-G as a biomarker of weakness for newborns and infants.

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.

Human Leukocyte Antigen-G: A Promising Prognostic Marker of Disease Progression to Improve the Control of Human African Trypanosomiasis

BACKGROUND

Human African trypanosomiasis (HAT) caused by Trypanosoma brucei gambiense can be diagnosed in the early hemolymphatic stage (stage 1 [S1]) or meningoencephalitic stage (stage 2 [S2]). Importantly, individuals harbouring high and specific antibody responses to Tbg antigens but negative parasitology are also diagnosed in the field (seropositive [SERO]). Whereas some develop the disease in the months following their initial diagnosis (SERO/HAT), others remain parasitologically negative for long periods (SERO) and are apparently able to control infection. Human leucocyte antigen (HLA)-G, an immunosuppressive molecule, could play a critical role in this variability of progression between infection and disease.

METHODS

Soluble HLA-G (sHLA-G) was measured in plasma for patients in the SERO (n = 65), SERO/HAT (n = 14), or HAT (n = 268) group and in cerebrospinal fluid for patients in S1 (n = 55), early S2 (n = 93), or late S2 (n = 110). Associations between these different statuses and the soluble level or genetic polymorphisms of HLA-G were explored.

RESULTS

Plasma sHLA-G levels were significantly higher in HAT (P = 6 × 10^-7) and SERO/HAT (P = .007) than SERO patients. No difference was observed between the SERO/HAT and HAT groups. Within the HAT group, specific haplotypes (HG010102 and HG0103) displayed increased frequencies in S1 (P = .013) and late S2 (P = .036), respectively.

CONCLUSIONS

These results strongly suggest the involvement of HLA-G in HAT disease progression. Importantly, high plasma sHLA-G levels in SERO patients could be predictive of subsequent disease development and could represent a serological marker to help guide therapeutic decision making. Further studies are necessary to assess the predictive nature of HLA-G and to estimate both sensitivity and specificity.

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.

HLA-G coding region and 3'untranslated region (3'UTR) in two Chinese Han populations

In this study, exons 2-4 and 3'untranslated region (3'UTR) of human leukocyte antigen (HLA)-G gene were investigated for 201 and 104 healthy unrelated Han samples recruited from Hunan Province, southern China and central Inner Mongolia Autonomous Region, northern China, respectively, using sequence-based typing and cloning methods. Totally 12 HLA-G alleles in the coding region, 9 variable sites in 3'UTR, 8 3'UTR haplotypes and 15 HLA-G extended haplotypes (EHs) incorporating the coding region and 3'UTR were observed. Very strong linkage disequilibrium (LD) was observed between HLA-A and HLA-G, and between HLA-G coding region and 3'UTR in each population (all global P=0.0000). Seven HLA-A-G haplotypes showed significant LD in both populations. Three HLA-G alleles in the coding region, 4 polymorphic sites in the 3'UTR, 3 3'UTR haplotypes and 4 HLA-G EHs differed significantly in their distributions between the 2 Chinese Han populations (all P≤0.0001). There was evidence for balancing selection acting on HLA-G 3'UTR positions +3010, +3142 and +3187 in the two populations. The NJ dendrograms demonstrated the existence of two basic HLA-G lineages and indicated that, HLA-G*01:01:01, the most common HLA-G allele, formed a separate lineage from other alleles. Our results shed new lights into HLA-G genetics among Chinese Han populations. The findings reported here are of importance for future studies related to post-transcriptional regulation of HLA-G allelic expression and the potential role of HLA-G in disease association in populations of Chinese ancestry.

Evolution of the levels of human leukocyte antigen G (HLA-G) in Beninese infant during the first year of life in a malaria endemic area: using latent class analysis

Background

HLA-G, a non-classical HLA class I antigen, is of crucial interest during pregnancy by inhibiting maternal immune response. Its role during infections is discussed, and it has been described that high levels of soluble HLA-G during childhood increase the risk of malaria. To explore more precisely interactions between soluble HLA-G and malaria, latent class analysis was used to test whether distinct sub-populations of children, each with distinctive soluble HLA-G evolutions may suggest the existence of groups presenting variable malaria susceptibility.

Method

A study was conducted in Benin from 2010 to 2013 and 165 children were followed from birth to 12 months. Evolution of soluble HLA-G was studied by the latent class method.

Results

Three groups of children were identified: one with consistently low levels of soluble HLA-G during follow-up, a second with very high levels and a last intermediate group. In all groups, low birth weight, high number of malaria infections and high exposure to malaria transmission were associated with high level of soluble HLA-G. Placental malaria was not. Presence of soluble HLA-G in cord blood increased the probability of belonging to the highest trajectory.

Conclusion

These results, together with previous ones, confirm the important role of HLA-G in the individual susceptibility to malaria. Assaying soluble HLA-G at birth could be a good indicator of newborns more fragile and at risk of infections during childhood.

Low variability at the HLA-E promoter region in the Brazilian population

Little information has been reported regarding the regulatory region of the HLA-E gene. Only a proximal segment (300 bp immediately before exon 1) has been described at IMGT/HLA database. We aimed to evaluate the genetic diversity of the promoter region by PCR amplification and DNA sequencing. We observed a pattern of sequencing interruption in the same position in all samples, suggesting the presence of a secondary structure hampering the DNA polymerase sliding during the sequencing process, which was confirmed after bioinformatics analysis. Considering the promoter region evaluated (nucleotide -440 to -1), only three variation sites were found, including one new variation site at position -104, and the concomitant deletions already described at positions -187 and -186. We concluded that the promoter region of the HLA-E gene presents few and rare variable sites in this population sample; however, the double-stranded branch formation hampered the evaluation of the distant promoter by conventional sequencing.

GABBR1 and SLC6A1, Two Genes Involved in Modulation of GABA Synaptic Transmission, Influence Risk for Alcoholism: Results from Three Ethnically Diverse Populations

BACKGROUND

Animal and human studies indicate that GABBR1, encoding the GABAB1 receptor subunit, and SLC6A1, encoding the neuronal gamma-aminobutyric acid (GABA) transporter GAT1, play a role in addiction by modulating synaptic GABA. Therefore, variants in these genes might predict risk/resilience for alcoholism.

METHODS

This study included 3 populations that differed by ethnicity and alcoholism phenotype: African American (AA) men: 401 treatment-seeking inpatients with single/comorbid diagnoses of alcohol and drug dependence, 193 controls; Finnish Caucasian men: 159 incarcerated alcoholics, half with comorbid antisocial personality disorder, 181 controls; and a community sample of Plains Indian (PI) men and women: 239 alcoholics, 178 controls. Seven GABBR1 tag single nucleotide polymorphisms were genotyped in the AA and Finnish samples; rs29220 was genotyped in the PI for replication. Also, a uniquely African, functional SLC6A1 insertion promoter polymorphism (IND) was genotyped in the AAs.

RESULTS

We found a significant and congruent association between GABBR1 rs29220 and alcoholism in all 3 populations. The major genotype (heterozygotes in AAs, Finns) and the major allele in PIs were significantly more common in alcoholics. Moreover, SLC6A1 IND was more abundant in controls, that is, the major genotype predicted alcoholism. An analysis of combined GABBR1 rs29220 and SLC6A1 IND genotypes showed that rs29220 heterozygotes, irrespective of their IND status, had an increased risk for alcoholism, whereas carriers of the IND allele and either rs29220 homozygote were more resilient.

CONCLUSIONS

Our results show that with both GABBR1 and SLC6A1, the minor genotypes/alleles were protective against risk for alcoholism. Finally, GABBR1 rs29220 might predict treatment response/adverse effects for baclofen, a GABAB receptor agonist.

Genetics of autoimmune diseases: insights from population genetics

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Autoimmune diseases (ADs) are a family of complex heterogeneous disorders with similar underlying mechanisms characterized by immune responses against self. Collectively, ADs are common, exhibit gender and ethnic disparities, and increasing incidence. As natural selection is an important influence on human genetic variation, and immune function genes are enriched for signals of positive selection, it is thought that the prevalence of AD risk alleles seen in different population is partially the result of differing selective pressures (for example, due to pathogens). With the advent of high-throughput technologies, new analytical methodologies and large-scale projects, evidence for the role of natural selection in contributing to the heritable component of ADs keeps growing. This review summarizes the genetic regions associated with susceptibility to different ADs and concomitant evidence for selection, including known agents of selection exerting selective pressure in these regions. Examples of specific adaptive variants with phenotypic effects are included as an evidence of natural selection increasing AD susceptibility. Many of the complexities of gene effects in different ADs can be explained by population genetics phenomena. Integrating AD susceptibility studies with population genetics to investigate how natural selection has contributed to genetic variation that influences disease risk will help to identify functional variants and elucidate biological mechanisms. As such, the study of population genetics in human population holds untapped potential for elucidating the genetic causes of human disease and more rapidly focusing to personalized medicine.