Extended HLA-G genetic diversity and ancestry composition in a Brazilian admixed population sample: Implications for HLA-G transcriptional control and for case-control association studies

Evidence for Epistatic Interaction between HLA-G and LILRB1 in the Pathogenesis of Nonsegmental Vitiligo

Vitiligo is the most frequent cause of depigmentation worldwide. Genetic association studies have discovered about 50 loci associated with disease, many with immunological functions. Among them is HLA-G, which modulates immunity by interacting with specific inhibitory receptors, mainly LILRB1 and LILRB2. Here we investigated the LILRB1 and LILRB2 association with vitiligo risk and evaluated the possible role of interactions between HLA-G and its receptors in this pathogenesis. We tested the association of the polymorphisms of HLA-G, LILRB1, and LILRB2 with vitiligo using logistic regression along with adjustment by ancestry. Further, methods based on the multifactor dimensionality reduction (MDR) approach (MDR v.3.0.2, GMDR v.0.9, and MB-MDR) were used to detect potential epistatic interactions between polymorphisms from the three genes. An interaction involving rs9380142 and rs2114511 polymorphisms was identified by all methods used. The polymorphism rs9380142 is an HLA-G 3'UTR variant (+3187) with a well-established role in mRNA stability. The polymorphism rs2114511 is located in the exonic region of LILRB1. Although no association involving this SNP has been reported, ChIP-Seq experiments have identified this position as an EBF1 binding site. These results highlight the role of an epistatic interaction between HLA-G and LILRB1 in vitiligo pathogenesis.

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.

Prediction of eye and hair pigmentation phenotypes using the HIrisPlex system in a Brazilian admixed population sample

Human pigmentation is a complex trait, probably involving more than 100 genes. Predicting phenotypes using SNPs present in those genes is important for forensic purpose. For this, the HIrisPlex tool was developed for eye and hair color prediction, with both models achieving high accuracy among Europeans. Its evaluation in admixed populations is important, since they present a higher frequency of intermediate phenotypes, and HIrisPlex has demonstrated limitations in such predictions; therefore, the performance of this tool may be impaired in such populations. Here, we evaluate the set of 24 markers from the HIrisPlex system in 328 individuals from Ribeirão Preto (SP) region, predicting eye and hair color and comparing the predictions with their real phenotypes. We used the HaloPlex Target Enrichment System and MiSeq Personal Sequencer platform for massively parallel sequencing. The prediction of eye and hair color was accomplished by the HIrisPlex online tool, using the default prediction settings. Ancestry was estimated using the SNPforID 34-plex to observe if and how an individual's ancestry background would affect predictions in this admixed sample. Our sample presented major European ancestry (70.5%), followed by African (21.1%) and Native American/East Asian (8.4%). HIrisPlex presented an overall sensitivity of 0.691 for hair color prediction, with sensitivities ranging from 0.547 to 0.782. The lowest sensitivity was observed for individuals with black hair, who present a reduced European contribution (48.4%). For eye color prediction, the overall sensitivity was 0.741, with sensitivities higher than 0.85 for blue and brown eyes, although it failed in predicting intermediate eye color. Such struggle in predicting this phenotype category is in accordance with what has been seen in previous studies involving HIrisPlex. Individuals with brown eye color are more admixed, with European ancestry decreasing to 62.6%; notwithstanding that, sensitivity for brown eyes was almost 100%. Overall sensitivity increases to 0.791 when a 0.7 threshold is set, though 12.5% of the individuals become undefined. When combining eye and hair prediction, hit rates between 51.3 and 68.9% were achieved. Despite the difficulties with intermediate phenotypes, we have shown that HIrisPlex results can be very helpful when interpreted with caution.

Variability at the 3' untranslated region of the HLA-G gene: a study on patients with AIDS and cytomegalovirus retinochoroiditis

Cytomegalovirus retinochoroiditis (CMV-R) is the primary cause of blindness among AIDS patients. Since HLA-G is associated with the modulation of the immune response, we hypothesized that variability at the 3′ untraslated region (3′UTR) of the gene could be implicated on the predisposition to CMV-R. We evaluated whether HLA-G 3′UTR influences CMV-R development in Brazilian AIDS patients. Peripheral blood DNA was obtained from two groups of patients: (1) AIDS exhibiting CMV-R (n = 40) and (2) AIDS without CMV-R (n = 147). HLA-G 3′UTR typing was performed using sequencing analysis. Allele, genotype and haplotype frequencies were evaluated using Fisher’s exact test accompanied by the calculation of the odds ratio and its 95% confidence interval (95% CI). The etiologic (EF) and preventive fractions were also estimated. Compared to AIDS patients without CMV-R, AIDS patients with CMV-R showed increased frequencies of the: (1) + 3001T allele, (2) the + 3001C/T genotype and (3) the UTR-17 (InsTTCCGTGACG) haplotype (EFs = 0.02–0.04). The UTR-3 (DelTCCGCGACG) haplotype was associated with protection against CMV-R development. Although the risk for developing CMR-V at the population level was relatively low (EF), the identification of HLA-G 3′UTR variation sites may help to further evaluate the role of post-transcriptional factors that may contribute to the existent immunosuppresion caused by HIV per se.

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.

The history behind the mosaic of the Americas

Focusing on literature published in 2018-2020, we review inferences about: (i) how ancient DNA is contributing to clarify the peopling of the Americas and the dispersal of its first inhabitants, (ii) how the interplay between environmental diversity and culture has influenced the genetic structure and adaptation of Andean and Amazon populations, (iii) how genetics has contributed to our understanding of the Pre-Columbian Tupi expansion in Eastern South America, (iv) the subcontinental origins and dynamics of Post-Columbian admixture in the Americas, and finally, (v) episodes of adaptive natural selection in the American continent, particularly in the high altitudes of the Andes.

Human Antigen Leucocyte (HLA)-G and HLA-E are differentially expressed in pancreatic disorders

BACKGROUND

Little information is available regarding the expression of the immunomodulatory Human Leukocyte Antigen (HLA)-G and -E molecules in pancreatic disorders.

AIM

To analyze HLA-G and -E expression in specimens of alcoholic chronic pancreatitis (ACP), idiopathic chronic pancreatitis (ICP), type 1 (T1D) and type 2 diabetes (T2D) and in histologically normal pancreas (HNP).

METHODS

HLA-G and -E expression (ACP = 30, ICP = 10, T1D = 10, T2D = 30 and HNP = 20) was evaluated by immunohistochemistry in three different areas (acini, islets and inflammatory infiltrate).

RESULTS

Acini and islets from HNP specimens exhibited higher HLA-G and -E expression compared to corresponding areas from all other patient groups. In inflammatory infiltrate, HLA-G and -E expression was observed only among the pancreatic disorders. We observed higher HLA-G and -E expression in acini from T2D compared to ACP, as well as higher HLA-G expression compared to ICP.

CONCLUSION

The decreased expression of HLA-G and -E in islets and acini together with the expression of these molecules in the inflammatory infiltrating cells were shared features among chronic inflammatory and autoimmune pancreatic disorders evaluated in this study, possibly reflecting tissue damage.

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.

HLAIb worldwide genetic diversity: New HLA-H alleles and haplotype structure description

The classical HLA class I genes (HLA Ia) were extensively studied because of their implication in clinical fields and anthropology. Less is known about worldwide genetic diversity and linkage disequilibrium for non-classical HLA class I genes (HLA Ib) and HLA pseudogenes. Notably, HLA-H, which is deleted in a fraction of the population, remains scarcely explored. The aims of this study were 1/ to get further insight into HLA-H genetic diversity and into how this variability potentially affects its expression and 2/ to define HLA Ib worldwide allelic diversity and linkage. Exome sequence data from the 1000 Genomes Project were used to define second field HLA-A, -E, -F, -G and -H typing using PolyPheMe software. Allelic and two-loci haplotype frequencies were estimated using Gene[Rate] software both at worldwide and continental levels. Eleven novel HLA-H alleles identified in exome data were validated by NGS performed on 25 genomic DNA samples from the same cohort. Phylogenetic analysis and frequency distribution of HLA-H alleles revealed three clades, each predominantly represented in Admixed American, European and East Asian populations, African populations and South Asian populations. Among these eleven novel alleles, two potentially encode complete transmembrane HLA proteins. We confirm the high LD between HLA-H and -A, and between HLA-H and -G, and show the three genes have distinct worldwide allelic distribution. Conversely, HLA-E and HLA-F both showed little LD, displayed restricted allelic diversity and practically no difference in their distribution across the planet. Our work thus reveals an unexpectedly high HLA-H genetic diversity, with alleles highly represented in Asia possibly encoding a functional HLA protein. Functional implication of these results remains to be explored, both in physiological and pathological contexts.