Detection of the signature of natural selection in humans: evidence from the Duffy blood group locus

The Duffy blood group locus, which encodes a chemokine receptor, is characterized by three alleles-FY*A, FY*B, and FY*O. The frequency of the FY*O allele, which corresponds to the absence of Fy antigen on red blood cells, is at or near fixation in most sub-Saharan African populations but is very rare outside Africa. The FST value for the FY*O allele is the highest observed for any allele in humans, providing strong evidence for the action of natural selection at this locus. Homozygosity for the FY*O allele confers complete resistance to vivax malaria, suggesting that this allele has been the target of selection by Plasmodium vivax or some other infectious agent. To characterize the signature of directional selection at this locus, we surveyed DNA sequence variation, both in a 1.9-kb region centered on the FY*O mutation site and in a 1-kb region 5-6 kb away from it, in 17 Italians and in a total of 24 individuals from five sub-Saharan African populations. The level of variation across both regions is two- to threefold lower in the Africans than in the Italians. As a result, the pooled African sample shows a significant departure from the neutral expectation for the number of segregating sites, whereas the Italian sample does not. The FY*O allele occurs on two major haplotypes in three of the five African populations. This finding could be due to recombination, recurrent mutation, population structure, and/or mutation accumulation and drift. Although we are unable to distinguish among these alternative hypotheses, it is likely that the two major haplotypes originated prior to selection on the FY*O mutation.     

Natural Selection Mediated Association of the Duffy (FY) Gene Polymorphisms with Plasmodium vivax Malaria in India

The Duffy (Fy) antigens act as receptors for chemokines as well as for Plasmodium vivax to invade human RBCs. A recent study has correlated the occurrence of the FY*A allele of Duffy gene with decreased susceptibility to vivax malaria, but no epidemiological correlation between the distribution of FY*A allele and incidences of vivax malaria has been established so far. Furthermore, if such correlations exist, whether natural selection has mediated the association, is an important question. Since India is highly endemic to P. vivax malaria with variable eco-climatic and varying vivax malaria epidemiology across different regions, such a question could well be answered in Indians. For this, we have genotyped the FY gene at the −33^rd and the 125^th nucleotide positions in 250 Indians sampled from six different zonal plus one tribal population covering the whole of India and studied possible correlations with eco-climatic and vivax malaria incidences. No FY*O allele was found, however, both the FY*A and FY*B alleles forming FY*A/FY*A, FY*A/FY*B and FY*B/FY*B genotypes were widely distributed among Indians. Five out of seven population samples significantly deviated from the Hardy-Weinberg equilibrium expectation, and two alleles (FY*A and FY*B) and the homozygote genotype, FY*B/FY*B were clinally distributed over the population coordinates. Furthermore, vivax malaria incidences over the past five years were significantly negatively and positively associated with the frequencies of the FY*A and FY*B alleles, respectively. The Northern Indians were highly differentiated from the other zonal population samples at the FY gene, as evidenced from the reconstructed Neighbor-Joining phylogenetic tree. The results specify the role of natural selection in the distribution of FY gene polymorphism in India. Furthermore, the hypotheses on the part of the FY*A allele in conferring protection to vivax malaria could be validated following population genetic studies in a vivax malaria epidemiological setting, such as India.     

Duffy antigen/receptor for chemokines (DARC): genotypes in Ashkenazi and non-Ashkenazi Jews in Israel.

Duffy genotypes were studied in Ashkenazi and non-Ashkenazi groups in Israel. The prevalence of the genotypes for the known polymorphic FY*A, FY*B, and FY*B(GATA-) alleles was similar in the two groups. The recently described FY*B(G298A) and FY*B(C265T) alleles were also found to be polymorphic. FY*B(G298A) was significantly more prevalent in the non-Ashkenazi group than in the Ashkenazi group (in 20% and 10% of FY*B, respectively). FY*B(C265T), which markedly diminishes the expression of Fy(b) antigen, was found in 3.5% of FY*B alleles, but only together with FY*B(G298A), consistent with previous suggestions that FY*B(C265T) occurred in the FY*B(G298A) allele. No difference in Duffy genotype distribution was found between schizophrenic and control groups. Duffy antigens are receptors for chemokines and bind Plasmodium vivax. Study of Duffy genotypes in additional populations might help in elucidating the possible functional significance of Duffy allele polymorphism.     

Cytochrome P4501A1 polymorphisms in the Amerindian and Mestizo populations of Mexico

Several polymorphisms in the CYP1A1 locus have been identified and their genotypes appear to exhibit population frequencies that depend on ethnicity. We studied two CYP1A1 polymorphic sites (position 4889 and 6235) in a group of 212 unrelated healthy individuals belonging to three different Mexican populations (106 Mexican Mestizos, 52 Teenek and 54 Mayos). Comparison among Mexican populations showed increased frequency of the *Ile allele (A on position 4889) in Mexican Mestizos when compared to Amerindians (p < 0.05). The analysis of position 6235 showed increased frequencies of *m2 (C in this position) allele in Teenek when compared to Mestizos and Mayos (p < 0.05) and of *m2/*m2 genotype when compared to Mestizos (p < 0.05). Amerindian populations (from Mexico and South America) presented the lowest frequencies of *Ile (position 4889) and *m1 (position 6235) alleles, however these frequencies vary according to the ethnic group studied. Mexican Amerindian groups together with other South Amerindian populations showed the highest frequencies for *Val at position 4889 and the *m2 allele at position 6235. The present study corroborates the high frequencies of*Val and *m2 alleles in the Amerindian populations and detects some differences between Mexican populations that correlate with linguistic differences. Our data could be helpful in understanding the distribution of these polymorphisms and in clarifying their roles as genetic and evolution markers in Amerindian populations.     

Cystic fibrosis gene variability in two southern Brazilian Amerindian populations: analysis of the deltaF508 mutation and the KM19 and XV2C haplotypes

The frequencies of the deltaF508 deletion, the most common cystic fibrosis mutation in Europeans and European-derived populations, and the XV2C and KM19 restriction fragment length polymorphisms that are tightly linked to the CFTR locus vary among populations. To determine the distribution of these extragenic markers and of the deltaF508 mutation, we analyzed 326 chromosomes of individuals from two South American Indian populations, the Guarani and the Kaingang. The allele and haplotype frequencies differed greatly between the two populations as well as among Amerindians and normal European Brazilians and European Brazilian cystic fibrosis patients. The absence of the deltaF508 mutation and the B haplotype are in agreement with the hypothesis that the deltaF508 mutation occurred after the divergence of these two populations. This finding is useful for populations containing a large Amerindian component and helps us to understand the origins of the deltaF508 deletion, the most common cystic fibrosis mutation in Europeans and European-derived populations, as well as the different incidences of cystic fibrosis in continental groups.     

Genotyping of the Duffy Blood Group among Plasmodium knowlesi-Infected Patients in Malaysia

The Duffy blood group is of major interest in clinical medicine as it plays an important role in Plasmodium knowlesi and Plasmodium vivax infection. In the present study, the distribution of Duffy blood group genotypes and allelic frequencies among P. knowlesi infected patients as well as healthy individuals in Peninsular Malaysia were determined. The blood group of 60 healthy blood donors and 51 P. knowlesi malaria patients were genotyped using allele specific polymerase chain reaction (ASP-PCR). The data was analyzed using Fisher''s exact test in order to assess the significance of the variables. Our results show a high proportion of the FY*A/FY*A genotype (>85% for both groups) and a high frequency of the FY*A allele (>90% for both groups). The FY*A/FY*A genotype was the most predominant genotype in both infected and healthy blood samples. The genotype frequency did not differ significantly between the donor blood and the malaria patient groups. Also, there was no significant correlation between susceptibility to P. knowlesi infection with any Duffy blood genotype.     

Analysis of HLA class II haplotypes in the Cayapa Indians of Ecuador: a novel DRB1 allele reveals evidence for convergent evolution and balancing selection at position 86.

PCR amplification, oligonucleotide probe typing, and sequencing were used to analyze the HLA class II loci (DRB1, DQA1, DQB1, and DPB1) of an isolated South Amerindian tribe. Here we report HLA class II variation, including the identification of a new DRB1 allele, several novel DR/DQ haplotypes, and an unusual distribution of DPB1 alleles, among the Cayapa Indians (N = 100) of Ecuador. A general reduction of HLA class II allelic variation in the Cayapa is consistent with a population bottle-neck during the colonization of the Americas. The new Cayapa DRB1 allele, DRB1*08042, which arose by a G-->T point mutation in the parental DRB1*0802, contains a novel Val codon (GTT) at position 86. The generation of DRB1*08042 (Val-86) from DRB1*0802 (Gly-86) in the Cayapa, by a different mechanism than the (GT-->TG) change in the creation of DRB1*08041 (Val-86) from DRB1*0802 in Africa, implicates selection in the convergent evolution of position 86 DR beta variants. The DRB1*08042 allele has not been found in > 1,800 Amerindian haplotypes and thus presumably arose after the Cayapa separated from other South American Amerindians. Selection pressure for increased haplotype diversity can be inferred in the generation and maintenance of three new DRB1*08042 haplotypes and several novel DR/DQ haplotypes in this population. The DPB1 allelic distribution in the Cayapa is also extraordinary, with two alleles, DPB1*1401, a very rare allele in North American Amerindian populations, and DPB1*0402, the most common Amerindian DPB1 allele, constituting 89% of the Cayapa DPB1.(ABSTRACT TRUNCATED AT 250 WORDS)     

Von Willebrand factor gene polymorphisms in three Brazilian ethnic groups

Von Willebrand factor gene polymorphisms were studied in three Brazilian ethnic groups: Euro-Brazilians, Afro-Brazilians, and Amerindians. Six polymorphic sites were analyzed: RsaI (exon 18), NlaIV (exon 20), HphI, KpnI, D1472H, and V1565L (all four in exon 28). The allele frequencies were significantly different between Euro- and Afro-Brazilians for RsaI, HphI, D1472H, and V1565L, while in Amerindians NlaIV and HphI showed significant differences among tribes. This is the first report of these allele frequencies in Amerindians. Eighteen haplotypes were observed, and they showed significant differences between Euro- and Afro-Brazilians, among Amerindian tribes, and among the three ethnic groups. These results furnish important background data for evolutionary and anthropological investigations; in addition, they will be useful for establishing the origin and molecular characterization of the different forms of von Willebrand disease, as well as for detecting carriers and offering genetic counseling to those with this condition.     

Heterogeneity of the Y chromosome in Afro-Brazilian populations

Sixteen biallelic markers (SRY10831a, SRY10831b, SRY4064, SRY2627, 92R7, P2, P3, M34, M9, M3, M2, YAP, M60, M89, M213, M216) located in the nonrecombinant region of the Y chromosome were analyzed in 209 individuals belonging to six Brazilian populations: four Afro-Brazilian populations, one population of white European descendants, and one population of Japanese descendants. The results showed that most of the Y chromosomes of the Afro-Brazilians were from sub-Saharan Africa and that the proportion of Y chromosomes of European origin was greater than that of Y chromosomes of Amerindian origin. No typical African or Amerindian haplogroup was detected among Japanese individuals, and only one white individual showed a typical African haplogroup. Haplogroup P-92R7, which is highly frequent in the Portuguese and Italian populations, was the most frequent among whites (54%), and haplogroup K-M9, which shows wide geographic distribution and is absent in Africa, was the most frequent among Japanese individuals (65.6%). The two semi-isolated Afro-Brazilian populations showed the highest and the lowest genetic diversity, respectively. These differences probably reflect the effect of greater or smaller gene flow between a small isolated group and other populations. These findings show that the process of admixture does not occur homogeneously, with a tendency toward preferential marriages within the ethnic group and a clear direction in unions between European men and Amerindian or African women in the past. The results agree with historical and social data about the formation of the Brazilian population and reveal some of the factors that contribute to its heterogeneity.     

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.     

Apolipoprotein E polymorphism in the Indian and Mestizo populations of Mexico

Apolipoprotein E (APOE) genotypes were determined in 75 Mazatecan Indians and 83 Mexican mestizos. APOE allele and genotype frequencies in Mazatecans and mestizos were similar, with high frequencies of the APOE*3 allele (0.900 and 0.915, respectively) and the E3/3 genotype (0.813 and 0.831, respectively) and an absence in both samples of the APOE*2 allele. Our data are similar to those previously described for Mexican-American and Mayan populations, which show the highest frequency worldwide of the APOE*3 allele and the E3/3 genotype. Mazatecans and mestizos also show a decreased frequency of the APOE*4 allele when compared to other Amerindian groups. The absence of the APOE*2 allele has also been reported in other Amerindian groups such as Mayans and Cayapa, whereas in Caucasians the average frequency of this allele is about 8%. Our data are in agreement with previous reports showing absence of the APOE*2 allele in Native American groups. These findings suggest that the APOE*2 allele was absent in humans from northern Asia who settled in the Arctic and populated the American continent.     

Blood groups in Native Americans: a look beyond ABO and Rh

The study presents comparisons between blood group frequencies beyond ABO and Rh blood systems in Native American populations and previously published data from Brazilian blood donors. The frequencies of Diego (c.2561C>T, rs2285644), Kell (c.578C>T, rs8176058), Duffy (c.125A>G, rs12075, c.1−67T>C, rs2814778) and Kidd (c.838A>G, rs1058396) variants in Kaingang (n=72) and Guarani (n=234) populations from Brazil (1990-2000) were obtained and compared with data from these populations sampled during the 1960s and with individuals of different Brazilian regions. Data showed high frequencies of DI*01 and FY*01 alleles: 11.8% and 57.6% in Kaingang and 6.8% and 75.7% in Guarani groups, respectively. The main results indicated: (1) reduction in genetic distance over time of Kaingang and Guarani in relation to other Brazilian populations is suggestive of ongoing admixture; (2) significant differences in some frequencies of blood group markers (especially Diego, Kidd and Duffy) in relation to Native Americans and individuals from different geographical regions of Brazil. Our study shows that the frequency of red blood cell polymorphisms in two Native American groups is very different from that of blood donors, when we evaluated blood groups different from ABO and Rh systems, suggesting that a better ethnic characterization of blood unit receptors is necessary.     

Complex Signatures of Natural Selection at the Duffy Blood Group Locus

The Duffy blood group locus (FY) has long been considered a likely target of natural selection, because of the extreme pattern of geographic differentiation of its three major alleles (FY*B, FY*A, and FY*O). In the present study, we resequenced the FY region in samples of Hausa from Cameroon (fixed for FY*O), Han Chinese (fixed for FY*A), Italians, and Pakistanis. Our goals were to characterize the signature of directional selection on FY*O in sub-Saharan Africa and to understand the extent to which natural selection has also played a role in the extreme geographic differentiation of the other derived allele at this locus, FY*A. The data from the FY region are compared with the patterns of variation observed at 10 unlinked, putatively neutral loci from the same populations, as well as with theoretical expectations from the neutral-equilibrium model. The FY region in the Hausa shows evidence of directional selection in two independent properties of the data (i.e., level of sequence variation and frequency spectrum), observations that are consistent with the FY*O mutation being the target. The Italian and Chinese FY data show patterns of variation that are very unusual, particularly with regard to frequency spectrum and linkage disequilibrium, but do not fit the predictions of any simple model of selection. These patterns may represent a more complex and previously unrecognized signature of positive selection.     

BF and C3 genetic polymorphisms in Kaingang Indians from southern Brazil.

The distribution of C3 and BF variants was determined in a sample of 239 Kaingang Indians. The corresponding gene frequencies were as follows: BF*S = 0.9393, BF*F = 0.0356, BF*S05 = 0.0251, C3*S = 0.9769, C3*F = 0.0231. The presence of the BF*S05 allele, which has previously been found only in a Brazilian population, suggests that this allele originated in Amerindians. The comparatively low degree of polymorphism with high frequencies of BF*S and C3*S is in accordance with the relatedness of the Kaingang with other Amerindians, Eskimos and Asian populations.     

Genetic profiling of genes from the oxidative stress pathway among North and South Indians

The case-control association study design has been extensively used for elucidating the genetic basis of complex traits. Considerable variation in frequencies of various gene polymorphisms has been reported across different populations and ethnic groups. Thus before beginning such studies, one must know the gene variants that exist in the population. Such information is not available for the ethnically distinct Indian population, which, on the basis of the languages spoken, can be further subdivided into Indo-Europeans (North Indians) and Dravidians (South Indians). In this study we provide information on allele and genotype frequencies, pairwise linkage disequilibrium, and predominant haplotypes in two populations (North India, n=96; South India, n=96) for several of the commonly investigated polymorphisms in the oxidative stress pathway genes. Of the 33 polymorphisms in 19 genes tested, significant differences in allele and genotype frequencies between the two populations were observed for SOD3 Ala58Thr, UCP1-3826 C/T, NOS3-786 T/C, and TNFA-308 G/A polymorphisms.     

Beta1 adrenergic receptor polymorphisms Arg389Gly and Ser49Gly in the Amerindian and Mestizo populations of Mexico

The beta1 adrenergic receptor genotypes (Ser49Gly and Arg389Gly) were determined in 190 individuals from 3 Mexican populations. Mestizos and Teenek present the highest frequencies for the *Arg allele and the lowest frequencies for the *Gly allele (Arg389Gly) compared to European, Asian, and African populations. Mayos present the highest frequency for the *Gly allele. The knowledge of the distribution of these alleles could help define the significance of these polymorphisms as genetic susceptibility markers in Amerindian populations.     

Interaction between ABO and haptoglobin systems in a Bengalee population.

Blood samples from 2,232 individuals of a Bengalee Caste Hindu population were investigated in an attempt to confirm the association between the ABO and haptoglobin (HP) systems previously found in populations of European origin. Indians differ from Europeans in having lower HP*1 and higher ABO*B frequencies. In spite of this, as in previous studies, a weak HP/ABO association was found with a significantly lower HP*1 allele frequency in blood group O versus other ABO groups.     

Electrophoretic and isoelectric focusing studies in Brazilian Indians: data on four systems

Three-hundred ninety-nine individuals living in seven populations of two Brazilian Indian tribes (Macushi and I?ana River Indians) were tested for the phosphoglucomutase 1 (PGM1), properdin factor B (BF), haptoglobin (HP), and alpha-1-antitrypsin (PI) systems. We observed significant internal heterogeneity in the two tribes for the PGM1 alleles and in the Macushi for the HP markers. Frequencies in three of the four systems (the exception being BF) also show clear differences in the Macushi and I?ana River Indians. Compared with other ethnic groups, South American Indians generally present high frequencies of PGM1*1B, BF*S, HP*1S, and PI*M3. On the other hand, PGM1*1A, PI*M1, and PI*M2 are reduced, and HP*1F is absent or rare. This is the first report about HP subtypes among American Indians.     

Association of Duffy Blood Group Gene Polymorphisms with IL8 Gene in Chronic Periodontitis

The antigens of the Duffy blood group system (DARC) act as a receptor for the interleukin IL-8. IL-8 plays an important role in the pathogenesis of chronic periodontitis due to its chemotactic properties on neutrophils. The aim of this study was to investigate a possible association of Duffy blood group gene polymorphisms with the -353T>A, -845T>C and -738T>A SNPs of the IL8 gene in chronic periodontitis. One hundred and twenty-four individuals with chronic periodontitis and 187 controls were enrolled. DNA was extracted using the salting-out method. The Duffy genotypes and IL8 gene promoter polymorphisms were investigated by PCR-RFLP. Statistical analyses were conducted using the Chi square test with Yates correction or Fisher''s Exact Test, and the possibility of associations were evaluated by odds ratio with a 95% confidence interval. When analyzed separately, for the Duffy blood group system, differences in the genotype and allele frequencies were not observed between all the groups analyzed; and, in nonsmokers, the -845C allele (3.6% vs. 0.4%), -845TC genotype (7.3% vs. 0.7%) and the CTA haplotype (3.6% vs. 0.4%) were positively associated with chronic periodontitis. For the first time to our knowledge, the polymorphisms of erythroid DARC plus IL8 -353T>A SNPs were associated with chronic periodontitis in Brazilian individuals. In Afro-Brazilians patients, the FY*02N.01 with IL8 -353A SNP was associated with protection to chronic periodontitis.     

Duffy Antigen Receptor for Chemokine (DARC) Polymorphisms and Its Involvement in Acquisition of Inhibitory Anti-Duffy Binding Protein II (DBPII) Immunity

The Plasmodium vivax Duffy binding protein (PvDBP) and its erythrocytic receptor, the Duffy antigen receptor for chemokines (DARC), are involved in the major P. vivax erythrocyte invasion pathway. An open cohort study to analyze DARC genotypes and their relationship to PvDBP immune responses was carried out in 620 volunteers in an agricultural settlement of the Brazilian Amazon. Three cross-sectional surveys were conducted at 6-month intervals, comprising 395, 410, and 407 subjects, respectively. The incidence rates of P. vivax infection was 2.32 malaria episodes per 100 person-months under survey (95% confidence interval [CI] of 1.92-2.80/100 person-month) and, of P. falciparum, 0.04 per 100 person-months (95% CI of 0.007–0.14/100 person-month). The distribution of DARC genotypes was consistent with the heterogeneous ethnic origins of the Amazon population, with a predominance of non-silent DARC alleles: FY*A > FY*B. The 12-month follow-up study demonstrated no association between DARC genotypes and total IgG antibodies as measured by ELISA targeting PvDBP (region II, DBPII or regions II–IV, DBPII-IV). The naturally acquired DBPII specific binding inhibitory antibodies (BIAbs) tended to be more frequent in heterozygous individuals carrying a DARC-silent allele (FY*B^ES). These results provide evidence that DARC polymorphisms may influence the naturally acquired inhibitory anti-Duffy binding protein II immunity.