Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Variability of innate immune system genes in Native American populations-relationship with history and epidemiology.

OBJECTIVES: The immune system of a host, defending him/her against invading pathogens, has two main subsystems: innate immunity and acquired immunity. There are several evidences showing that Native American populations are immunologically different from non-Native populations. Our aim was to describe the variability of innate immune system genes in Native American populations. MATERIALS AND METHODS: We investigated heterozygozities and patterns of population differentiation (FST ) of 14 polymorphisms related to the innate immune response in five Native American populations (Aché, Guarani-Kaiowá, Guarani-Ñandeva, Kaingang, and Xavante) and the results were compared with the three major world population data (YRI, CEU, and CHB) available at the 1,000 genomes database. RESULTS: Mean heterozygosities ranged between 0.241 ± 0.057 (Aché) and 0.343 ± 0.033 (Kaingang), but no significant differences were observed (Friedman test, P = 0.197). Mean heterozygosities were also not significantly different when Amerindians were pooled and compared with the 1000 genomes populations (Friedman test, P = 0.506). When the Native American populations were grouped as Amerindians, a significantly higher FST value (0.194) was observed between the Amerindian and African populations. The Ewens-Watterson neutrality test showed that these markers are not under strong selective pressure. DISCUSSION: Native American populations present similar levels of heterozygosity as those of other continents, but are different from Africans in the frequency of polymorphisms of innate immune genes. This higher differentiation is probably due to demographic processes that occurred during the out-of-Africa event.

Distribution of a novel CYP2C haplotype in Native American populations.

The CYP2C19 gene, located in the CYP2C cluster, encodes the major drug metabolism enzyme CYP2C19. This gene is highly polymorphic and no-function (CYP2C19*2 and CYP2C19*3), reduced function (CYP2C19*9) and increased function (CYP2C19*17) star alleles (haplotypes) are commonly used to predict CYP2C19 metabolic phenotypes. CYP2C19*17 and the genotype-predicted rapid (RM) and ultrarapid (UM) CYP2C19 metabolic phenotypes are absent or rare in several Native American populations. However, discordance between genotype-predicted and pharmacokinetically determined CYP2C19 phenotypes in Native American cohorts have been reported. Recently, a haplotype defined by rs2860840T and rs11188059G alleles in the CYP2C cluster has been shown to encode increased rate of metabolism of the CYP2C19 substrate escitalopram, to a similar extent as CYP2C19*17. We investigated the distribution of the CYP2C:TG haplotype and explored its potential impact on CYP2C19 metabolic activity in Native American populations. The study cohorts included individuals from the One Thousand Genomes Project AMR superpopulation (1 KG_AMR), the Human Genome Diversity Project (HGDP), and from indigenous populations living in Brazil (Kaingang and Guarani). The frequency range of the CYP2C:TG haplotype in the study cohorts, 0.469 to 0.598, is considerably higher than in all 1 KG superpopulations (range: 0.014-to 0.340). We suggest that the high frequency of the CYP2C:TG haplotype might contribute to the reported discordance between CYP2C19-predicted and pharmacokinetically verified CYP2C19 metabolic phenotypes in Native American cohorts. However, functional studies involving genotypic correlations with pharmacokinetic parameters are warranted to ascertain the importance of the CYP2C:TG haplotype.

Stability or variation? Patterns of lactase gene and its enhancer region distributions in Brazilian Amerindians.

Lactase persistence (LP) is the phenotypic trait in which lactase secretion is maintained during adulthood. LP is due to mutations in the LCT enhancer region, located 14-kb upstream of the gene. In Europeans, the -13910*T allele is associated with LP. In Africans this allele is rare while other mutations in this same region were related to LP. The LCT is highly polymorphic in human populations, but so far Brazilian Amerindians had not been investigated for these polymorphisms or for the presence of LP mutations. We describe the genetic diversity of the LCT region and the presence of LP enhancer mutations in four native Brazilian populations (Guarani-Kaiowá, Guarani-Ñandeva, Kaingang, and Xavante). Twelve polymorphisms were genotyped by PCR-based methods. The -13910*T allele varied from 0.5% in the Xavante to 7.6% in the Guarani-Ñandeva. These frequencies probably derive from European sources and they correlate with non-native admixture proportions previously estimated for these groups. But since admixture is virtually absent in the Xavante, we suggest that the presence of the LP allele could have been determined by a de novo mutation. No other mutations in the -14 kb enhancer region were found. The LCT was highly polymorphic in the present sample showing 15 haplotypes with a heterogeneous distribution among the four Amerindian populations. This diversity could be due to drift, as indicated by the neutrality test performed.

Single Nucleotide Variants as Proxies for HLA-A*31:01 in Native American Populations.

Carbamazepine triggers dermatologic hypersensitivity reactions, associated with specific human leukocyte antigens (HLAs), especially HLA-B*15:02 and HLA-A*31:01. Previous efforts to identify single nucleotide variants (SNVs) with high predictive value as HLA proxies, revealed that rs1061235 and rs17179220 fulfill these requirements for HLA-A*31:01 in some but not all populations. Herein we explored the predictive performance of rs1061235 and rs17179220 as HLA-A*31:01 tags in populations of Native American ancestry, which are largely underrepresented in pharmacogenomic studies. The study cohorts comprised the overall Admixed American superpopulation of the 1000 Genomes Project (1 KG_AMR), a subcohort of individuals with predominant Native American ancestry (1 KG_NAT), the Native American population of the Human Genome Diversity Project (HGDP), plus Kaingang (KRC) and Guarani (GRC and GKW) adults from indigenous reservation areas in Brazil. The diversity of cohorts is reflected in the range of frequencies of HLA-A*31:01 (0.02-0.65), rs1061235 (0.03-0.13) and rs17179220 (0.12-0.66), as well as in the predictive performance of these SNVs as HLA-A*31:01 proxies. NPV (negative predictive value), the metric of primary interest for pharmacogenetic-informed carbamazepine prescription was maximal (NPV = 1.0) for both SNVs in 1 KG_AMR and 1 KG_NAT, for rs17179220, but not rs1061235 (NPV = 0.91) in HGDP, and for rs17179220 in GWK, but not GRC (NPV = 0.88) or KRC (NPV = 0.80). Collectively, the data support the notion that rs1061235 and rs17179220 are not optimal proxies for HLA-A*31:01 across populations of Native American ancestry.

Geographic patterns of genome admixture in Latin American Mestizos.

The large and diverse population of Latin America is potentially a powerful resource for elucidating the genetic basis of complex traits through admixture mapping. However, no genome-wide characterization of admixture across Latin America has yet been attempted. Here, we report an analysis of admixture in thirteen Mestizo populations (i.e. in regions of mainly European and Native settlement) from seven countries in Latin America based on data for 678 autosomal and 29 X-chromosome microsatellites. We found extensive variation in Native American and European ancestry (and generally low levels of African ancestry) among populations and individuals, and evidence that admixture across Latin America has often involved predominantly European men and both Native and African women. An admixture analysis allowing for Native American population subdivision revealed a differentiation of the Native American ancestry amongst Mestizos. This observation is consistent with the genetic structure of pre-Columbian populations and with admixture having involved Natives from the area where the Mestizo examined are located. Our findings agree with available information on the demographic history of Latin America and have a number of implications for the design of association studies in population from the region.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans.

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

Contrasting patterns of nuclear and mtDNA diversity in Native American populations.

We report an integrated analysis of nuclear (autosomal, X- and Y-chromosome) short tandem repeat (STR) data and mtDNA D-loop sequences obtained in the same set of 22 Native populations from across the Americas. A north to south gradient of decreasing population diversity was observed, in agreement with a settlement of the Americas from the extreme northwest of the continent. This correlation is stronger with "least cost distances," which consider the coasts as facilitators of migration. Continent-wide estimates of population structure are highest for the Y-chromosome and lowest for the autosomes, consistent with the effective size of the different marker systems examined. Population differentiation is highest in East South America and lowest in Meso America and the Andean region. Regional analyses suggest a deviation from mutation-drift equilibrium consistent with population expansion in Meso America and the Andes and population contraction in Northwest and East South America. These data hint at an early divergence of Andean and non-Andean South Americans and at a contrasting demographic history for populations from these regions.

Dopamine receptor D4 allele distribution in Amerindians: a reflection of past behavior differences?

The DRD4 variable number of tandem repeats (VNTR) allele distribution of 172 Guarani (Kaiowá and Ñandeva subgroups) and Kaingang Brazilian Amerindians is reported. These results are integrated with those previously obtained for this ethnic group. Allele frequencies for the three populations are within the interval observed for 15 other Native American populations and show intermediate values between those observed in Amazonia and Patagonia. Significant differences in allele distribution between recent past hunter-gatherer and agriculturalist populations are observed, with an increase of the 7R allele among hunter-gatherers (P < 0.001). Analysis of molecular variance (AMOVA) and pairwise F(ST) data suggest three distinct sectors for the genetic landscape of Native South America: Andes, Center/Southeast region, and Amazonia. Common traits among hunter-gatherers such as novelty-seeking temperament, hyperactivity, and impulsivity could have been important and advantageous in new environments during America's prehistoric colonization.

Genetic variation and population structure in native Americans.

We examined genetic diversity and population structure in the American landmass using 678 autosomal microsatellite markers genotyped in 422 individuals representing 24 Native American populations sampled from North, Central, and South America. These data were analyzed jointly with similar data available in 54 other indigenous populations worldwide, including an additional five Native American groups. The Native American populations have lower genetic diversity and greater differentiation than populations from other continental regions. We observe gradients both of decreasing genetic diversity as a function of geographic distance from the Bering Strait and of decreasing genetic similarity to Siberians--signals of the southward dispersal of human populations from the northwestern tip of the Americas. We also observe evidence of: (1) a higher level of diversity and lower level of population structure in western South America compared to eastern South America, (2) a relative lack of differentiation between Mesoamerican and Andean populations, (3) a scenario in which coastal routes were easier for migrating peoples to traverse in comparison with inland routes, and (4) a partial agreement on a local scale between genetic similarity and the linguistic classification of populations. These findings offer new insights into the process of population dispersal and differentiation during the peopling of the Americas.

Demographic history and selection at HLA loci in Native Americans.

The American continent was the last to be occupied by modern humans, and native populations bear the marks of recent expansions, bottlenecks, natural selection, and population substructure. Here we investigate how this demographic history has shaped genetic variation at the strongly selected HLA loci. In order to disentangle the relative contributions of selection and demography process, we assembled a dataset with genome-wide microsatellites and HLA-A, -B, -C, and -DRB1 typing data for a set of 424 Native American individuals. We find that demographic history explains a sizeable fraction of HLA variation, both within and among populations. A striking feature of HLA variation in the Americas is the existence of alleles which are present in the continent but either absent or very rare elsewhere in the world. We show that this feature is consistent with demographic history (i.e., the combination of changes in population size associated with bottlenecks and subsequent population expansions). However, signatures of selection at HLA loci are still visible, with significant evidence selection at deeper timescales for most loci and populations, as well as population differentiation at HLA loci exceeding that seen at neutral markers.

HLA-DR and HLA-DQ alleles in patients from the south of Brazil: markers for leprosy susceptibility and resistance.

BACKGROUND: Many epidemiological studies have shown that the genetic factors of the host play a role in the variability of clinical response to infection caused by M. leprae. With the purpose of identifying genes of susceptibility, the present study investigated the possible role of HLA-DRB1 and DQA1/DQB1 alleles in susceptibility to leprosy, and whether they account for the heterogeneity in immune responses observed following infection in a Southern Brazilian population. METHODS: One hundred and sixty-nine leprosy patients and 217 healthy controls were analyzed by polymerase chain reaction amplification and reverse hybridization with sequence-specific oligonucleotide probes and sequence-specific primers (One Lambda, CA, USA). RESULTS: There was a positive association of HLA-DRB1*16 (*1601 and *1602) with leprosy per se (7.3% vs. 3.2%, P = 0.01, OR = 2.52, CI = 1.26-5.01), in accord with previous serological studies, which showed DR2 as a marker of leprosy. Although, HLA-DQA1*05 frequency (29.8% vs. 20.9%, P = 0.0424, OR = 1.61, CI = 1.09-2.39) was higher in patients, and HLA-DQA1*02 (3.0% vs. 7.5%, P = 0.0392, OR = 0.39, CI = 0.16 - 0.95) and HLA-DQA1*04 (4.0% vs. 9.1%, P = 0.0314, OR = 0.42, CI = 0.19 - 0.93) frequencies lower, P-values were not significant after the Bonferroni's correction. Furthermore, HLA-DRB1*1601 (9.0% vs. 1.8%; P = 0.0016; OR = 5.81; CI = 2.05-16.46) was associated with susceptibility to borderline leprosy compared to control group, and while HLA-DRB1*08 (11.2% vs. 1.2%; P = 0.0037; OR = 12.00; CI = 1.51 - 95.12) was associated with susceptibility to lepromatous leprosy, when compared to tuberculoid leprosy, DRB1*04 was associated to protection. CONCLUSION: These data confirm the positive association of HLA-DR2 (DRB1*16) with leprosy per se, and the protector effect of DRB1*04 against lepromatous leprosy in Brazilian patients.

Linkage disequilibrium patterns and genetic structure of Amerindian and non-Amerindian Brazilian populations revealed by long-range X-STR markers.

The extent of X-chromosome linkage disequilibrium (LD) was studied in a southern Brazilian population, and in a pool of samples from Amerindian populations. For this purpose, 11 microsatellites, located mostly in a Xq region comprising approximately 86 Mb was investigated. The lower Amerindian gene diversity associated with significant differences between the populations studied indicated population structure as the main cause for the higher LD values in the Amerindian pool. On the other hand, the LD levels of the non-Amerindian Brazilian sample, although less extensive than that of the Amerindians, were probably determined by admixture events. Our results indicated that different demographic histories have significant effects on LD levels of human populations, and provide a first approach to the X-chromosome ancestry of Amerindian and non-Amerindian Brazilian populations, being valuable for future studies involving mapping and population genetic studies.

Analysis of the CCR5 gene coding region diversity in five South American populations reveals two new non-synonymous alleles in Amerindians and high CCR5*D32 frequency in Euro-Brazilians.

The CC chemokine receptor 5 (CCR5) molecule is an important co-receptor for HIV. The effect of the CCR5*D32 allele in susceptibility to HIV infection and AIDS disease is well known. Other alleles than CCR5*D32 have not been analysed before, neither in Amerindians nor in the majority of the populations all over the world. We investigated the distribution of the CCR5 coding region alleles in South Brazil and noticed a high CCR5*D32 frequency in the Euro-Brazilian population of the Paraná State (9.3%), which is the highest thus far reported for Latin America. The D32 frequency is even higher among the Euro-Brazilian Mennonites (14.2%). This allele is uncommon in Afro-Brazilians (2.0%), rare in the Guarani Amerindians (0.4%) and absent in the Kaingang Amerindians and the Oriental-Brazilians. R223Q is common in the Oriental-Brazilians (7.7%) and R60S in the Afro-Brazilians (5.0%). A29S and L55Q present an impaired response to ß-chemokines and occurred in Afro- and Euro-Brazilians with cumulative frequencies of 4.4% and 2.7%, respectively. Two new non-synonymous alleles were found in Amerindians: C323F (g.3729G > T) in Guarani (1.4%) and Y68C (g.2964A > G) in Kaingang (10.3%). The functional characteristics of these alleles should be defined and considered in epidemiological investigations about HIV-1 infection and AIDS incidence in Amerindian populations.

Two new mutations of the human BCHE gene (IVS3-14T>C and L574fsX576).

The genetic variation of human butyrylcholinesterase is associated with the majority of prolonged cases of apnea in patients submitted to the muscle relaxant succinylcholine. The present study reports two new mutations of the BCHE gene in 346 Euro-Brazilians: IVS3-14T>C found in five heterozygotes (allele frequency: 0.72+/-0.32%) and L574fsX576 found in one heterozygote (allele frequency: 0.14+/-0.14%). These two variants were not found in 85 Guarani Amerindians. It is not expected that the IVS3-14T>C mutation may interfere in the splicing process and that the mutation found in exon 4 (L574fsX576) may disturb BChE tetramerization and activity.

Y-STR analysis in Brazilian and South Amerindian populations.

A sample of 203 Brazilian males from Rio Grande do Sul (RS), the Brazilian southernmost state, was typed for 11 Y-STR markers (DYS19, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS385, DYS437, DYS438, and DYS439). We also typed 42 individuals from two South Amerindian tribes (Kaingang and Guarani) to use the data as parental Amerindian contribution to our analyses. Gene and haplotypic diversities were estimated, with the South Amerindian samples showing smaller values for these parameters than Brazilians. To obtain a more comprehensive picture of the genetic structure of the Brazilian population as a whole, the Y-STR data from the RS sample was compared with those already published. No genetic substructuring was observed in the comparisons performed. Multidimensional scaling confirmed the proposed European source of most Y-chromosome Brazilian patrilineages.

Association between KIR genes and dust mite sensitization in a Brazilian population.

BACKGROUND: Killer cell immunoglobulin-like receptors (KIRs), found on the surface of natural killer (NK) cells, play a key role in controlling the innate response. Such response depends on a series of cellular interactions between these receptors and HLA activating/inhibiting ligands. Atopic diseases have been associated with genes that regulate cytokine production and HLA genes, which may either protect or predispose to hypersensitivity. OBJECTIVE: To verify an association study of KIR genes with sensitization to the following mites: Dermatophagoides farinae, Dermatophagoides pteronyssinus, and Blomia tropicalis. METHODS: A total of 341 children aged up to 14 years, were classified as mite-sensitive or mite-insensitive after undergoing a skin prick test for immediate allergic reactions. The presence/absence of KIR genes and their human leukocyte antigen (HLA) ligands was determined by polymerase chain reaction-sequence specific oligonucleotide (PCR-SSO) with the commercial kit LabType™ using Luminex™. RESULTS: The frequencies of KIR genes and their respective class I HLA ligands and the frequency of haplotypes were performed in sensitive and insensitive individuals, and no significant differences were found. CONCLUSION: Our results suggest no influence of KIR genes on resistance/susceptibility to sensitization to dust mites.

[Association between HLA and leukemia in a mixed Brazilian population]. / Associação entre HLA e leucemia em uma população brasileira de etnia mista.

OBJECTIVE: The main purpose of this study was to investigate the class I HLA antigens and class II HLA allele frequencies in 164 patients with leukemia: 35 patients with ALL (acute lymphoid leukemia), 50 with AML (acute myeloid leukemia) and 78 with CML (chronic myeloid leukemia). METHODS: The genotyping of class I HLA was performed by microlymphocytotoxicity and of class II by PCR-SSP (polymerase chain reaction - sequence specific of primers) (One Lambda, Canoga Park, CA, USA). RESULTS: In patients with LLA, frequencies of HLA-B45 and HLA-B56 were higher (P = 0.02; OR = 3.13; 95%IC = 0.94-10.44; P = 0.03; OR = 3.61; 95%IC = 0.47-27.64, respectively), than in controls. In patients with AML, the frequency of HLA-B7 (P = 0.01; OR = 2.41; 95%IC = 1.25-4.67) was higher than in controls. The presence of HLA-B45 (P= 0.01; OR = 3.29; 95%IC = 1.46-7.40), HLA-DRB1*04 (P = 0.002; OR = 2.17; 95%IC = 1.36-3.46) and HLA-DRB1*08 (P = 0.004; OR = 2.36; 95%IC = 1.34-4.16) was associated to increased risk of CML developing. CONCLUSION: Our results suggest that variants of HLA confer susceptibility to the same forms of leukemia, and could provide new tools for the investigation of genetics and etiology of this disease.

Association between HLA genes and dust mite sensitivity in a Brazilian population.

BACKGROUND: Type I hypersensitivity, also known as IgE-mediated allergy, is a complex, multifactorial condition whose onset and severity are influenced by both genetic and environmental factors. Mite allergens stimulate the production of humoral response (IgE), especially in children, which is closely involved in atopic asthma and rhinitis. OBJECTIVE: This study aimed to investigate the association between HLA class I (-A, -B, and -C), and HLA class II (-DRB1) genes in individuals sensitive to dust mites (Dermatophagoides farinae, Dermatophagoides pteronyssinus, or Blomia tropicalis) and mite-insensitive controls. METHODS: 396 participants were grouped as mite-sensitive and mite-insensitive according to immediate hypersensitivity as determined by skin-prick tests, and to HLA genotyping by polymerase chain reaction-sequence specific oligonucleotide (PCR-SSO). RESULTS: After chi-square heterogeneity testing no significant differences were observed in HLA-A, B, and C genes, except for the HLA-DRB1 locus, which, showed a negative association for DRB1∗04, between mite-sensitive and mite-insensitive individuals. In high resolution, DRB1∗04:11 allele was significantly different from all other results (P=0.0042, OR=0.26, and 95%CI=0.09-0.70). The analysis stratified by etiologic agent confirmed these associations. CONCLUSION: Our results suggest a possible association between HLA-DRB1 genes and hypersensitivity to dust mites.

Human T-cell lymphotropic virus type II in Guaraní Indians, Southern Brazil.

Human T-cell lymphotropic virus type II (HTLV-II) is found in many New World Indian groups on the American continent. In Brazil, HTLV-II has been found among urban residents and Indians in the Amazon region, in the North. Guaraní Indians in the South of Brazil were studied for HTLV-I/II infection. Among 52 individuals, three (5.76%) showed positive anti-HTLV-II antibodies (enzyme-linked immunosorbent assay and Western blot). This preliminary report is the first seroepidemiological study showing HTLV-II infection among Indians in the South of Brazil.