KIR2DL4 genetic diversity in a Brazilian population sample: implications for transcription regulation and protein diversity in samples with different ancestry backgrounds.

KIR2DL4 is an important immune modulator expressed in natural killer cells; HLA-G is its main ligand. We have characterized the KIR2DL4 genetic diversity by considering the promoter, all exons, and all introns in a highly admixed Brazilian population sample and by using massively parallel sequencing. We introduce a molecular method to amplify and to sequence the complete KIR2DL4 gene. To avoid the mapping bias and genotype errors commonly observed in gene families, we have developed and validated a bioinformatic pipeline designed to minimize these errors and applied it to survey the variability of 220 individuals from the State of São Paulo, southeastern Brazil. We have also compared the KIR2DL4 genetic diversity in the Brazilian cohort with the diversity previously reported by the 1000Genomes consortium. KIR2DL4 presents high linkage disequilibrium throughout the gene, with coding sequences associated with specific promoters. There are few but divergent promoter haplotypes. We have also detected many new KIR2DL4 sequences, all bearing nucleotide exchanges in introns and encoding previously described proteins. Exons 3 and 4, which encode the external domains, are the most variable. The ancestry background influences the KIR2DL4 allele frequencies and must be considered for association studies regarding KIR2DL4.

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

HLA-E coding and 3' untranslated region variability determined by next-generation sequencing in two West-African population samples.

HLA-E is a non-classical Human Leucocyte Antigen class I gene with immunomodulatory properties. Whereas HLA-E expression usually occurs at low levels, it is widely distributed amongst human tissues, has the ability to bind self and non-self antigens and to interact with NK cells and T lymphocytes, being important for immunosurveillance and also for fighting against infections. HLA-E is usually the most conserved locus among all class I genes. However, most of the previous studies evaluating HLA-E variability sequenced only a few exons or genotyped known polymorphisms. Here we report a strategy to evaluate HLA-E variability by next-generation sequencing (NGS) that might be used to other HLA loci and present the HLA-E haplotype diversity considering the segment encoding the entire HLA-E mRNA (including 5'UTR, introns and the 3'UTR) in two African population samples, Susu from Guinea-Conakry and Lobi from Burkina Faso. Our results indicate that (a) the HLA-E gene is indeed conserved, encoding mainly two different protein molecules; (b) Africans do present several unknown HLA-E alleles presenting synonymous mutations; (c) the HLA-E 3'UTR is quite polymorphic and (d) haplotypes in the HLA-E 3'UTR are in close association with HLA-E coding alleles. NGS has proved to be an important tool on data generation for future studies evaluating variability in non-classical MHC genes.

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.

Insights on the HLA-G evolutionary history provided by a nearby Alu insertion.

The AluyHG element belongs to the AluYb8 subfamily. It is a polymorphic insertion, located approximately 20 kb from the HLA-G 3'-untranslated region (3'-UTR), which has been used for evolution studies because it exhibits identity for descendants and it is still polymorphic in the human genome. To understand the evolutionary mechanisms acting on HLA-G, we evaluated the presence or absence of the AluyHG element, associating this variable site with others observed at HLA-G coding, 3'-UTR, or both regions in four distinct populations (Brazilian, French, Congolese, and Senegalese). The results were compared with the 1000Genomes Consortium data. The worldwide AluyHG frequencies showed an increment, starting lower in Africa and increasing following distance and time of human dispersion out of Africa. The same haplotype pattern was observed in all populations, indicating that most of the HLA-G haplotypes already detected were originated earlier in Africa, before Homo sapiens dispersion. The AluyHG insertion was associated with the G*01:01:01:01/UTR-1 haplotype, with rare recombinants. Despite its high frequency in worldwide populations, the G*01:01:01:01/UTR-1 haplotype should be very recent. The low frequency of recombinants indicates that the rate of recombination at the HLA-G gene is very low.