Genetic diversity of the LILRB1 and LILRB2 coding regions in an admixed Brazilian population sample.

Leukocyte immunoglobulin (Ig)-like receptors (LILR) LILRB1 and LILRB2 may play a pivotal role in maintaining self-tolerance and modulating the immune response through interaction with classical and nonclassical HLA molecules. Although both diversity and natural selection patterns over HLA genes have been extensively evaluated, little information is available concerning the genetic diversity and selection signatures on the LILRB1/2 regions. Therefore, we identified the LILRB1/2 genetic diversity using next-generation sequencing in a population sample from São Paulo State, Brazil. We identified 58 LILRB1 Single Nucleotide Variants (SNVs), which gave rise to 13 haplotypes, and 41 LILRB2 SNVs arranged into 11 haplotypes. Although we may not exclude as a possible effect of population structure, we found evidence of either positive or purifying selection on LILRB1/2 coding regions. Some residues in both proteins showed to be under the effect of positive selection, suggesting that amino acid replacements in these proteins resulted in beneficial functional changes. Finally, we have revealed that allelic variation (six and five amino acid exchanges in LILRB1 and LILRB2, respectively) affects the structure and/or stability of both molecules. Nonetheless, LILRB2 has shown higher average stability, with no D1/D2 residue affecting protein structure. Overall, our findings demonstrate that LILRB1 and LILRB2 are as polymorphic as HLA class Ib genes and provide strong evidence supporting the directional selection regime hypothesis.

Hla-C genetic diversity and evolutionary insights in two samples from Brazil and Benin.

Human leukocyte antigen-C (HLA-C) is a classical HLA class I molecule that binds and presents peptides to cytotoxic T lymphocytes in the cell surface. HLA-C has a dual function because it also interacts with Killer-cell immunoglobulin-like receptors (KIR) receptors expressed in natural killer and T cells, modulating their activity. The structure and diversity of the HLA-C regulatory regions, as well as the relationship among variants along the HLA-C locus, are poorly addressed, and few population-based studies explored the HLA-C variability in the entire gene in different population samples. Here we present a molecular and bioinformatics method to evaluate the entire HLA-C diversity, including regulatory sequences. Then, we applied this method to survey the HLA-C diversity in two population samples with different demographic histories, one highly admixed from Brazil with major European contribution, and one from Benin with major African contribution. The HLA-C promoter and 3'UTR were very polymorphic with the presence of few, but highly divergent haplotypes. These segments also present conserved sequences that are shared among different primate species. Nucleotide diversity was higher in other segments rather than exons 2 and 3, particularly around exon 5 and the second half of the 3'UTR region. We detected evidence of balancing selection on the entire HLA-C locus and positive selection in the HLA-C leader peptide, for both populations. HLA-C motifs previously associated with KIR interaction and expression regulation are similar between both populations. Each allele group is associated with specific regulatory sequences, reflecting the high linkage disequilibrium along the entire HLA-C locus in both populations.

HLA-F displays highly divergent and frequent haplotype lineages associated with different mRNA expression levels.

HLA-F is one of the most conserved loci among the HLA gene family. The exact function of HLA-F is still under investigation. HLA-F might present tolerogenic features, participate in the stabilization of HLA molecules in open conformation, and also participate in the recycling of HLA molecules. Here we evaluate the variability and haplotype structure of the HLA-F distal promoter segment (from -1893 to -943) and how this segment is correlated with the coding region. Variability at the promoter segment was surveyed in 196 Brazilian samples using second-generation sequencing. The HLA-F promoter region presents two major haplotype lineages. Most of the variable sites are in perfect linkage and associated with a single promoter haplotype, here named F∗distal-C. This haplotype is associated with F∗01:01:02 alleles, while alleles from the F∗01:01:01 or F∗01:03 groups present closely related promoter sequences. F∗distal-C is quite frequent in Brazil and in worldwide populations, with frequencies ranging from 8.41% at the Iberian Population in Spain to 34.34% in Vietnam. F∗distal-C is also present in Neanderthal and Denisovan samples. In silico analyses demonstrated that F∗distal-C presents a different transcription factor binding profile compared with other HLA-F promoters. Moreover, individuals carrying this haplotype present higher HLA-F mRNA expression levels. Functional studies are required to define the exact mechanism underlying this higher HLA-F mRNA expression level associated with F∗distal-C and F∗01:01:02 alleles.

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-E regulatory and coding region variability and haplotypes in a Brazilian population sample.

The HLA-E gene is characterized by low but wide expression on different tissues. HLA-E is considered a conserved gene, being one of the least polymorphic class I HLA genes. The HLA-E molecule interacts with Natural Killer cell receptors and T lymphocytes receptors, and might activate or inhibit immune responses depending on the peptide associated with HLA-E and with which receptors HLA-E interacts to. Variable sites within the HLA-E regulatory and coding segments may influence the gene function by modifying its expression pattern or encoded molecule, thus, influencing its interaction with receptors and the peptide. Here we propose an approach to evaluate the gene structure, haplotype pattern and the complete HLA-E variability, including regulatory (promoter and 3'UTR) and coding segments (with introns), by using massively parallel sequencing. We investigated the variability of 420 samples from a very admixed population such as Brazilians by using this approach. Considering a segment of about 7kb, 63 variable sites were detected, arranged into 75 extended haplotypes. We detected 37 different promoter sequences (but few frequent ones), 27 different coding sequences (15 representing new HLA-E alleles) and 12 haplotypes at the 3'UTR segment, two of them presenting a summed frequency of 90%. Despite the number of coding alleles, they encode mainly two different full-length molecules, known as E*01:01 and E*01:03, which corresponds to about 90% of all. In addition, differently from what has been previously observed for other non classical HLA genes, the relationship among the HLA-E promoter, coding and 3'UTR haplotypes is not straightforward because the same promoter and 3'UTR haplotypes were many times associated with different HLA-E coding haplotypes. This data reinforces the presence of only two main full-length HLA-E molecules encoded by the many HLA-E alleles detected in our population sample. In addition, this data does indicate that the distal HLA-E promoter is by far the most variable segment. Further analyses involving the binding of transcription factors and non-coding RNAs, as well as the HLA-E expression in different tissues, are necessary to evaluate whether these variable sites at regulatory segments (or even at the coding sequence) may influence the gene expression profile.

HLA-G variability and haplotypes detected by massively parallel sequencing procedures in the geographicaly distinct population samples of Brazil and Cyprus.

The HLA-G molecule presents immunomodulatory properties that might inhibit immune responses when interacting with specific Natural Killer and T cell receptors, such as KIR2DL4, ILT2 and ILT4. Thus, HLA-G might influence the outcome of situations in which fine immune system modulation is required, such as autoimmune diseases, transplants, cancer and pregnancy. The majority of the studies regarding the HLA-G gene variability so far was restricted to a specific gene segment (i.e., promoter, coding or 3' untranslated region), and was performed by using Sanger sequencing and probabilistic models to infer haplotypes. Here we propose a massively parallel sequencing (NGS) with a bioinformatics strategy to evaluate the entire HLA-G regulatory and coding segments, with haplotypes inferred relying more on the straightforward haplotyping capabilities of NGS, and less on probabilistic models. Then, HLA-G variability was surveyed in two admixed population samples of distinct geographical regions and demographic backgrounds, Cyprus and Brazil. Most haplotypes (promoters, coding, 3'UTR and extended ones) were detected both in Brazil and Cyprus and were identical to the ones already described by probabilistic models, indicating that these haplotypes are quite old and may be present worldwide.

Association of HLA-G 3' untranslated region variants with type 1 diabetes mellitus.

Besides the well recognized association of HLA-DRB1 and DQB1 alleles with type 1 diabetes mellitus (T1D), linkage studies have identified a gene region close to the non-classical class I HLA-G gene as an independent susceptibility marker. HLA-G is constitutively expressed in the endocrine compartment of the human pancreas and may play a role in controlling autoimmune responses. We evaluated the genetic diversity of the 3' untranslated region (3'UTR) of HLA-G, which have been associated with HLA-G mRNA post-transcriptional regulation, in 120 Brazilian T1D patients and in 120 healthy controls. We found the +3001 T allele was observed only in T1D patients. Notably, the +3001 T allele was in linkage disequilibrium with polymorphic sites associated with low production of HLA-G mRNA or soluble HLA-G levels. Moreover, T1D patients showed a low frequency of the HLA-G 3'UTR-17 (14bpINS/+3001T/+3003T/+3010C/+3027C/+3035T/+3142G/+3187A/+3196C). The +3010 CC genotype and the UTR-3 haplotype (14bpDEL/+3001C/+3003T/+3010C/+3027C/+3035C/+3142G/+3187A/+3196C), associated with low and moderate soluble HLA-G expression, respectively, were underrepresented in patients. The decreased expression of HLA-G at the pancreas level should be detrimental in individuals genetically prone to produce less HLA-G.

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.

A comprehensive study of polymorphic sites along the HLA-G gene: implication for gene regulation and evolution.

HLA-G molecule plays an important role on immune response regulation and has been implicated on the inhibition of T and natural killer cell cytolytic function and inhibition of allogeneic T-cell proliferation. Due to its immune-modulator properties, the HLA-G gene expression has been associated with the outcome of allograft and of autoimmune, infectious, and malignant disorders. Several lines of evidence indicate that HLA-G polymorphisms at the 5'-upstream regulatory region (5' URR) and 3'-untranslated region (3' UTR) may influence the HLA-G expression levels. Because Brazilians represent one of the most heterogeneous populations in the world with the widest HLA-G coding region variability already detected among the studied populations, a high level of variability and haplotype diversity would be expected in Brazilians. On this basis, the 5' URR, coding, and 3' UTR variability were evaluated in a Brazilian series consisting of 100 healthy bone marrow donors, as well as the linkage disequilibrium pattern along the gene and the extended haplotypes encompassing several gene segment variations. The HLA-G locus seems to present six different HLA-G lineages showing functional variations mainly in nucleotides of the regulatory regions. Differences were observed at the 5' URR in positions that either coincide with or are close to transcription factor-binding sites and at the 3' UTR mainly in positions that have already been reported to influence HLA-G mRNA availability. We report several lines of evidence for balancing selection acting on the regulatory regions, which may indicate that these HLA-G lineages may be related to the differential HLA-G expression profiles.

In silico analysis of microRNAS targeting the HLA-G 3' untranslated region alleles and haplotypes.

It has been reported that microRNAs (miRNA) may have allele-specific targeting for the 3' untranslated region (3' UTR) of the HLA-G locus. In a previous study, we reported 11 3'UTR haplotypes encompassing the 14-bp insertion/deletion polymorphism and seven SNPs (+3003 T/C, +3010 C/G, +3027 C/A, +3035 C/T, +3142 C/G, +3187 A/G, and +3196 C/G), of which only the +3142 C/G SNP has been reported to influence the binding of miRNAs. Using bioinformatics analyses, we identified putative miRNA-binding sites considering the haplotypes encompassing these eight polymorphic sites, and we ranked the lowest free energies that could potentially lead to an mRNA degradation or translational repression. When a specific haplotype or a particular SNP was associated with a miRNA-binding site, we defined a free energy difference of 4 kcal/mol between alleles to classify them energetically distant. The best results were obtained for the miR-513a-5p, miR-518c*, miR-1262 and miR-92a-1*, miR-92a-2*, miR-661, miR-1224-5p, and miR-433 miRNAs, all influencing one or more of the +3003, +3010, +3027, and +3035 SNPs. The miR-2110, miR-93, miR-508-5p, miR-331-5p, miR-616, miR-513b, and miR-589* miRNAs targeted the 14-bp fragment region, and miR-148a, miR-19a*, miR-152, mir-148b, and miR-218-2 also influenced the +3142 C/G polymorphism. These results suggest that these miRNAs might play a relevant role on the HLA-G expression pattern.