Genotypic and Allelic Variability in CYP19A1 among Populations of African and European Ancestry

CYP19A1 facilitates the bioconversion of estrogens from androgens. CYP19A1 intron single nucleotide polymorphisms (SNPs) may alter mRNA splicing, resulting in altered CYP19A1 activity, and potentially influencing disease susceptibility. Genetic studies of CYP19A1 SNPs have been well documented in populations of European ancestry; however, studies in populations of African ancestry are limited. In the present study, ten ‘candidate’ intronic SNPs in CYP19A1 from 125 African Americans (AA) and 277 European Americans (EA) were genotyped and their frequencies compared. Allele frequencies were also compared with HapMap and ASW 1000 Genomes populations. We observed significant differences in the minor allele frequencies between AA and EA in six of the ten SNPs including rs10459592 (p<0.0001), rs12908960 (p<0.0001), rs1902584 (p = 0.016), rs2470144 (p<0.0001), rs1961177 (p<0.0001), and rs6493497 (p = 0.003). While there were no significant differences in allele frequencies between EA and CEU in the HapMap population, a 1.2- to 19-fold difference in allele frequency for rs10459592 (p = 0.004), rs12908960 (p = 0.0006), rs1902584 (p<0.0001), rs2470144 (p = 0.0006), rs1961177 (p<0.0001), and rs6493497 (p = 0.0092) was observed between AA and the Yoruba (YRI) population. Linkage disequilibrium (LD) blocks and haplotype clusters that is unique to the EA population but not AA was also observed. In summary, we demonstrate that differences in the allele frequencies of CYP19A1 intron SNPs are not consistent between populations of African and European ancestry. Thus, investigations into whether CYP19A1 intron SNPs contribute to variations in cancer incidence, outcomes and pharmacological response seen in populations of different ancestry may prove beneficial.     

High-density single-nucleotide polymorphism maps of the human genome

Here we report a large, extensively characterized set of single-nucleotide polymorphisms (SNPs) covering the human genome. We determined the allele frequencies of 55,018 SNPs in African Americans, Asians (Japanese-Chinese), and European Americans as part of The SNP Consortium's Allele Frequency Project. A subset of 8333 SNPs was also characterized in Koreans. Because these SNPs were ascertained in the same way, the data set is particularly useful for modeling. Our results document that much genetic variation is shared among populations. For autosomes, some 44% of these SNPs have a minor allele frequency > or =10% in each population, and the average allele frequency differences between populations with different continental origins are less than 19%. However, the several percentage point allele frequency differences among the closely related Korean, Japanese, and Chinese populations suggest caution in using mixtures of well-established populations for case-control genetic studies of complex traits. We estimate that approximately 7% of these SNPs are private SNPs with minor allele frequencies <1%. A useful set of characterized SNPs with large allele frequency differences between populations (>60%) can be used for admixture studies. High-density maps of high-quality, characterized SNPs produced by this project are freely available.     

Three novel mtDNA restriction site polymorphisms allow exploration of population affinities of African Americans

To develop informative tools for the study of population affinities in African Americans, we sequenced the hypervariable segments I and II (HVS I and HVS II) of mitochondrial DNA (mtDNA) from 96 Sierra Leoneans; European Americans; rural, Gullah-speaking African Americans; urban African Americans living in Charleston, South Carolina; and Jamaicans. We identified single nucleotide polymorphisms (SNPs) exhibiting ethnic affinities, and developed restriction endonuclease tools to screen these SNPs. Here we show that three HVS restriction site polymorphisms (RSPs), EcoRV, FokI, and MfeI, exhibit appreciable differences in frequency (average delta = 0.4165) between putative African American parental populations (i.e., extant Africans living in Sierra Leone and European Americans). Estimates of European American mtDNA admixture, calculated from haplotypes composed of these three novel RSPs, show a cline of increasing admixture from Gullah-speaking African American (m = 0.0300) to urban Charleston African American (m = 0.0689) to West Coast African American (m = 0.1769) populations. This haplotype admixture in the Gullahs is the lowest recorded to date among African Americans, consistent with previous studies using autosomal markers. These RSPs may become valuable new tools in the study of ancestral affinities and admixture dynamics of African Americans.     

Population frequencies of single nucleotide polymorphisms (SNPs) in immuno-modulatory genes

Inherited polymorphisms in immuno-modulatory genes may contribute to variations in immune function and genetic susceptibility for complex diseases, including cancer. We report results from a comprehensive study to discover novel single nucleotide polymorphisms (SNPs) and to estimate allelic frequency for both novel and known coding and regulatory region SNPs in genes encoding proteins that have been implicated in the immune response to tumors. We identified 12 novel nucleotide substitution variants and one deletion variant in 17 genes analyzed (TGFBETA;R, BETA;2M, IFNGAMMA;, TNFALPHA;, TNFALPHA;R, LTALPHA;, IL-6, IL-12, IL-2, IL-1ALPHA;, IL-1BETA;, IL-1RN, IL-10, CTLA4, CD40L, FAS and FASL). We determined the frequency of these novel polymorphisms, as well as 17 previously identified polymorphisms, in a control sample of 158 individuals, approximately half of which were Caucasian (n = 74) and half of which were African American (n = 84). Significant differences in allele frequencies were observed between the two racial groups for 13/17 genes tested. These allelic variations maybe associated with alterations in immune function and thus susceptibility to a number of complex disease states such as cancer.     

Characterization of CYP1A2, CYP2C19, CYP3A4 and CYP3A5 polymorphisms in South Brazilians

Potential causes of variability in drug response include intrinsic factors such as ethnicity and genetic differences in the expression of enzymes that metabolize drugs, such as those from Cytochrome P450 (CYPs) superfamily. Pharmacogenetic studies search for genetic differences between populations since relevant alleles occur with varying frequencies among different ethnic populations. The Brazilian population is one of the most heterogeneous in the world, resulting from multiethnic admixture of Amerindians, Europeans, and Africans across centuries. Since the knowledge of CYP allele frequency distributions is relevant to pharmacogenetic strategies and these data are scarce in the Brazilian population, this study aimed to describe genotype and allele distributions of 15 single nucleotide polymorphisms (SNPs) at CYP 1A2, 2C19, 3A4, and 3A5 genes in African and European descents from South Brazil. A sample of 179 healthy individuals of European and African ancestry was genotyped by the MassARRAY SNP genotyping system. CYP3A5*3, CYP1A2*1F, CYP3A4*1B, and CYP2C19*2 were the most frequent alleles found in our sample. Significant differences in genotype and allelic distribution between African and European descents were observed for CYP3A4 and CYP3A5 genes. CYP3A4*1B was observed in higher frequency in African descents (0.379) than in European descents (0.098), and European descents showed higher frequency of CYP3A5*3 (0.810) than African descents (0.523). Our results indicate that only a few polymorphisms would have impact in pharmacogenetic testing in South Brazilians. Further studies with larger sample sizes are required also among other Brazilian regions.     

Glutathione S-transferase M1 and T1 polymorphisms in Brazilian African descendants

The glutathione S-transferase gene family has an important role in the biotransformation and detoxification of different xenobiotics and endogenous compounds. Two polymorphic genes of this family, GSTM1 and GSTT1, present null alleles that consequently do not produce the respective enzyme when the genotype is homozygous. These polymorphisms are also interesting for population dynamics studies because they have great frequency variations among different ethnic groups and have been reported worldwide. The distribution of these alleles in urban and Amerindian populations in Brazil has been described, but none of those studies reported on African-descended rural populations. The aim of this study was to analyze the genotype frequency distribution of the GSTM1 and GSTT1 null alleles in an urban sample from the Federal District (n = 91) and in four semi-isolated African-descended populations: Mocambo (n = 55), Rio das R?s (n = 117), Riacho de Sacutiaba (n = 34), and Kalunga (n = 68). The GSTM1 and GSTT1 null genotype frequencies in these populations range from 17% to 35% for GSTM1 and from 22% to 44% for GSTT1. These values are similar to those described in other African and African-descended populations. Despite this range, there is no distribution difference among the analyzed populations. Combined GSTM1 and GSTT1 null genotype frequencies range from 6% to 13% and are similar to European-derived populations, suggesting admixture with this ethnic group. This can be interpreted as a European contribution to these African-descended populations. Regarding the urban population in the Federal District, our results suggest an important African and European contribution.     

Genetic Variation in TLR Genes in Ugandan and South African Populations and Comparison with HapMap Data

Genetic epidemiological studies of complex diseases often rely on data from the International HapMap Consortium for identification of single nucleotide polymorphisms (SNPs), particularly those that tag haplotypes. However, little is known about the relevance of the African populations used to collect HapMap data for study populations conducted elsewhere in Africa. Toll-like receptor (TLR) genes play a key role in susceptibility to various infectious diseases, including tuberculosis. We conducted full-exon sequencing in samples obtained from Uganda (n = 48) and South Africa (n = 48), in four genes in the TLR pathway: TLR2, TLR4, TLR6, and TIRAP. We identified one novel TIRAP SNP (with minor allele frequency [MAF] 3.2%) and a novel TLR6 SNP (MAF 8%) in the Ugandan population, and a TLR6 SNP that is unique to the South African population (MAF 14%). These SNPs were also not present in the 1000 Genomes data. Genotype and haplotype frequencies and linkage disequilibrium patterns in Uganda and South Africa were similar to African populations in the HapMap datasets. Multidimensional scaling analysis of polymorphisms in all four genes suggested broad overlap of all of the examined African populations. Based on these data, we propose that there is enough similarity among African populations represented in the HapMap database to justify initial SNP selection for genetic epidemiological studies in Uganda and South Africa. We also discovered three novel polymorphisms that appear to be population-specific and would only be detected by sequencing efforts.     

Identification, replication, and fine-mapping of Loci associated with adult height in individuals of african ancestry

Adult height is a classic polygenic trait of high heritability (h ² ∼0.8). More than 180 single nucleotide polymorphisms (SNPs), identified mostly in populations of European descent, are associated with height. These variants convey modest effects and explain ∼10% of the variance in height. Discovery efforts in other populations, while limited, have revealed loci for height not previously implicated in individuals of European ancestry. Here, we performed a meta-analysis of genome-wide association (GWA) results for adult height in 20,427 individuals of African ancestry with replication in up to 16,436 African Americans. We found two novel height loci (Xp22-rs12393627, P = 3.4×10⁻¹² and 2p14-rs4315565, P = 1.2×10⁻⁸). As a group, height associations discovered in European-ancestry samples replicate in individuals of African ancestry (P = 1.7×10⁻⁴ for overall replication). Fine-mapping of the European height loci in African-ancestry individuals showed an enrichment of SNPs that are associated with expression of nearby genes when compared to the index European height SNPs (P<0.01). Our results highlight the utility of genetic studies in non-European populations to understand the etiology of complex human diseases and traits.     

A high-density admixture map for disease gene discovery in african americans

Admixture mapping (also known as "mapping by admixture linkage disequilibrium," or MALD) provides a way of localizing genes that cause disease, in admixed ethnic groups such as African Americans, with approximately 100 times fewer markers than are required for whole-genome haplotype scans. However, it has not been possible to perform powerful scans with admixture mapping because the method requires a dense map of validated markers known to have large frequency differences between Europeans and Africans. To create such a map, we screened through databases containing approximately 450000 single-nucleotide polymorphisms (SNPs) for which frequencies had been estimated in African and European population samples. We experimentally confirmed the frequencies of the most promising SNPs in a multiethnic panel of unrelated samples and identified 3011 as a MALD map (1.2 cM average spacing). We estimate that this map is approximately 70% informative in differentiating African versus European origins of chromosomal segments. This map provides a practical and powerful tool, which is freely available without restriction, for screening for disease genes in African American patient cohorts. The map is especially appropriate for those diseases that differ in incidence between the parental African and European populations.     

Linkage disequilibrium and allele-frequency distributions for 114 single-nucleotide polymorphisms in five populations

Single-nucleotide polymorphisms (SNPs) may be extremely important for deciphering the impact of genetic variation on complex human diseases. The ultimate value of SNPs for linkage and association mapping studies depends in part on the distribution of SNP allele frequencies and intermarker linkage disequilibrium (LD) across populations. Limited information is available about these distributions on a genomewide scale, particularly for LD. Using 114 SNPs from 33 genes, we compared these distributions in five American populations (727 individuals) of African, European, Chinese, Hispanic, and Japanese descent. The allele frequencies were highly correlated across populations but differed by >20% for at least one pair of populations in 35% of SNPs. The correlation in LD was high for some pairs of populations but not for others (e.g., Chinese American or Japanese American vs. any other population). Regardless of population, average minor-allele frequencies were significantly higher for SNPs in noncoding regions (20%-25%) than for SNPs in coding regions (12%-16%). Interestingly, we found that intermarker LD may be strongest with pairs of SNPs in which both markers are nonconservative substitutions, compared to pairs of SNPs where at least one marker is a conservative substitution. These results suggest that population differences and marker location within the gene may be important factors in the selection of SNPs for use in the study of complex disease with linkage or association mapping methods.     

Evolution of a length polymorphism in the human PER3 gene, a component of the circadian system

Period homologue 3 (PER3) is a component of the mammalian circa-dian system, although its precise role is unknown. A biallelic variable number tandem repeat (VNTR) polymorphism exists in human PER3, consisting of 4 or 5 repeats of a 54-bp sequence in a region encoding a putative phosphorylation domain. This polymorphism has previously been reported to associate with diurnal preference ("morningness" and "eveningness") and delayed sleep-phase syndrome. We have investigated the global allele frequencies of this variant in ethnically distinct indigenous populations. All populations were polymorphic, with the shorter (4-repeat) allele ranging in frequency from 0.19 (Papua New Guinea) to 0.89 (Mongolia). To investigate if allele frequency has been influenced by natural selection, the authors 1) tested for a correlation with latitude and mean annual insolation (incident sunlight energy), using classical markers to correct for historical population differentiation; and they 2) compared allele-frequency difference between European American, African American, and East Asian populations, as measured using F(ST), to an empirical null distribution of F(ST)values based on a genome-wide dataset of single nucleotide polymorphisms (SNPs) of presumed neutral loci that were previously typed by The SNP Consortium. The variation in allele frequencies between indigenous populations did not show a pattern that would indicate selective pressure on PER3resulting from day-length variation or mean annual insolation, and the allele-frequency difference between European Americans, African Americans, and East Asians was not an outlier when compared to the distribution for presumed neutral SNPs. We therefore find no evidence for differential or balancing selection in the contemporary pattern of global PER3allele frequencies.     

Ethnic differences in the linkage disequilibrium and distribution of single-nucleotide polymorphisms in 35 candidate genes for cardiovascular diseases

Single-nucleotide polymorphisms (SNPs) are commonly used to study genetics for common diseases and predict pharmacological response. The selection of likely informative SNPs in association studies depends on their allele frequencies and on the linkage disequilibrium (LD) between SNPs, both of which may show interethnic differences. Among three populations consisting of 207 Chinese, 858 French, and 395 Spanish, we compared the allele frequency distributions of 64 intragenic SNPs of 35 candidate genes for cardiovascular diseases. Twenty-eight of these SNPs from 12 genes were also examined for intragenic LD. About 20% of SNPs were restricted to Europeans, being monomorphic in Chinese, among them mostly nonsynonymous coding SNPs and noncoding SNPs. Only 1.6% of SNPs were specific in Chinese, commensurate with the detection of these SNPs almost exclusively in Caucasians. Similarly, these SNPs were more often rare (<0.1 minor allele frequency) in Chinese (44.3%) than in Europeans (31.1%). The variant allele frequencies and intermarker LDs in terms of D' and Delta(2) were highly correlated between French and Spanish populations (r = 0.98-0.99, p < 0.001). However, only moderate correlations of allele frequencies and D' were found between the Chinese and the European populations (r = 0.7 and 0.3, respectively) despite a high correlation of Delta(2) values (r = 0.8). These results suggest that ethnic considerations are important in the selection of SNPs for association studies of candidate genes, as this may affect the power of the study as well as the likelihood of asking relevant questions and getting medically meaningful answers.     

Association of FTO Gene Variants With Adiposity in African‐American Adolescents

The prevalence of obesity continues to increase significantly, with the largest rise in the African‐American adolescents. Genetic contributions to obesity are being identified with the advent of genome‐wide association studies (GWAS). Specifically, variants of the fat mass and obesity associated (FTO) gene have been associated with obesity in populations of European descent. The studies in African Americans have been inconclusive. To further evaluate the association of the FTO gene and adiposity in African Americans, we genotyped 47 single‐nucleotide polymorphisms (SNPs), including seven SNPs previously reported to be significant in the literature in a cohort consisting of 561 non‐Hispanic white and 497 African‐American individuals. Analysis of our data showed 17 SNPs to be associated with BMI Z‐score (BMI‐Z) in our study population. The strongest association was found in the African Americans. The most significant SNP was rs8057044, which was associated with BMI‐Z in the African Americans (P = 0.00054). SNP rs9939609 was found to be significant in the non‐Hispanic white population (P = 0.028). Our data confirm the association between FTO and adiposity suggesting that FTO is a childhood obesity susceptibility gene. Our data also identify a novel SNP of the FTO gene (rs8057044) that is associated with measures of adiposity in the African‐American population.     

Linkage disequilibrium analysis identifies an FGFR1 haplotype-tag SNP associated with normal variation in craniofacial shape

Mutations in FGFR1 and TWIST1 have been reported to affect the timing of calvarial suture fusion resulting in craniosynostosis and facial abnormalities. We screened nonpathologic populations for genetic polymorphisms that may associate with normal craniofacial variation. We identified 17 single-nucleotide polymorphisms (SNPs) in FGFR1, 6 of which were novel (g.8591855G-->A, g.8593685G-->A, g.8602303C-->T, g.8602475A-->G (p.Ile293Val), g.8605849C-->T, g.8607868G-->A). No SNPs were found in TWIST1. FGFR1 SNP haplotypes were reconstructed for Caucasian, Asian, Australian Aboriginal, and African American populations. All populations shared two linkage disequilibrium blocks, with one haplotype-tag SNP (htSNP) tagging each block. The htSNP g.8592931G-->C was found to have a significant negative correlation with the cephalic index for all populations (R = -0.187, p = 0.036), with larger correlations in Asians and females. This finding is a starting point in the identification of a set of SNPs that can be genotyped to determine both normal and disease craniofacial phenotypes.     

Inflammatory gene variants in the Tsimane, an indigenous Bolivian population with a high infectious load

The Tsimane of lowland Bolivia are an indigenous forager-farmer population living under conditions resembling pre-industrial European populations, with high infectious morbidity, high infection and inflammation, and shortened life expectancy. Analysis of 917 persons ages 5 to 60+ showed that allele frequencies of 9 SNPs examined in the apolipoprotein E (apoE), C-reactive protein (CRP), and interleukin-6 (IL-6) genes differed from some European, African, and north Asian-derived populations. The apoE2 allele was absent, whereas four SNPs related to CRP and IL-6 were monomorphic: CRP (rs1800947, rs3093061, and rs3093062) and IL-6 (rs1800795). No significant differences in apoE, CRP, and IL-6 variants across age were found CRP levels were higher in carriers of two CRP proinflammatory SNPs, whereas they were lower in carriers of apoE4. Taken together the evidence for (1) different allele frequencies between the Tsimane and other populations and (2) the correlations of CRP and apoE alleles with blood CRP may suggest that these variants are under selection in response to a high infection environment.     

Single nucleotide polymorphisms derived from ancestral populations show no evidence for biased diversity estimates in Drosophila melanogaster

Single nucleotide polymorphisms (SNPs) are about to become one of the most popular genetic markers for genetic model organisms. To test the usefulness of SNPs for estimating genetic diversity, we surveyed three genomic regions in two Drosophila melanogaster populations, one from Africa and one European, collected in Austria. Diversity estimates based on the full SNP set indicated higher levels of variability in the African than in the European flies. When the analysis was based on the European SNP set, European and African flies had similar levels of variability. Interestingly, this bias was not observed for diversity estimates using SNPs derived from the ancestral African population. This result suggests that diversity estimates based on SNPs from ancestral populations could provide a general strategy to avoid biased SNP diversity estimates. Finally, the potential of SNPs for nonmodel organisms is discussed.     

A/G Gln20Arg (exon 1) and G/A Val156Met (exon 5) polymorphisms of the human orosomucoid 1 gene in Mexico

The human orosomucoid 1 gene (ORM1) codes an alpha-1-acid glycoprotein that has been classified as an acute-phase reactive protein, and a major drug-binding serum component, as well as an immunomodulatory protein with genetic polymorphisms. Evaluation of ORM variation through isoelectric focusing and immunobloting has revealed a world-wide distribution of the ORM1 F and ORM1 S alleles. We evaluated and examined the genetic characteristics of two Mexican populations that have different anthropological and cultural antecedents, examining two ORM1 genotypes (exon 1 - A/G (Gln20Arg) and exon 5 G/A (Val156Met)) in 145 individuals, using nested polymerase chain reaction, sequencing, and restricted fragment length polymorphism. Mexican Mestizos had higher frequencies of the exon 1 A allele (P = 0.020) and AA genotype (P = 0.018) and lower frequency of the G allele (P = 0.020) when compared to Teenek Amerindians. When we examined exon 5 G/A (Val156Met) polymorphisms, we found significantly higher frequencies of the G allele (P = 0.0007) and the GG genotype (P = 0.0003) in the Mexican Mestizo population. The Teenek population had a significantly higher frequency of the A allele than has been reported for Chinese and African (P < 0.05) populations, and the G/A genotype was more frequently found in this Mexican population than in Chinese, African and European populations (P < 0.05).     

Polymorphisms of DNA repair genes ERCC2 and XRCC1 in populations of Russia

We have conducted a comparative study of allele frequencies of single nucleotide polymorphisms (SNPs) rs1799793 and rs13181 of the ERCC2 gene as well as rs1799782 and rs25487 of the XRCC1 gene in population samples from European regions of Russia as well as in populations of Izhemsk and Priluzsk Komi and Yakuts. Significant differences in the distribution of polymorphic variants of the ERCC2 gene were demonstrated between populations of Yakuts and populations of Russians and Komi. In case of XRCC1 gene Izhemsk Komi population exhibited dissimilar allele frequencies compared to other populations.     

Haplotypes and linkage disequilibrium at the phenylalanine hydroxylase locus, PAH, in a global representation of populations

Because defects in the phenylalanine hydroxylase gene (PAH) cause phenylketonuria (PKU), PAH was studied for normal polymorphisms and linkage disequilibrium soon after the gene was cloned. Studies in the 1980s concentrated on European populations in which PKU was common and showed that haplotype-frequency variation exists between some regions of the world. In European populations, linkage disequilibrium generally was found not to exist between RFLPs at opposite ends of the gene but was found to exist among the RFLPs clustered at each end. We have now undertaken the first global survey of normal variation and disequilibrium across the PAH gene. Four well-mapped single-nucleotide polymorphisms (SNPs) spanning approximately 75 kb, two near each end of the gene, were selected to allow linkage disequilibrium across most of the gene to be examined. These SNPs were studied as PCR-RFLP markers in samples of, on average, 50 individuals for each of 29 populations, including, for the first time, multiple populations from Africa and from the Americas. All four sites are polymorphic in all 29 populations. Although all but 5 of the 16 possible haplotypes reach frequencies >5% somewhere in the world, no haplotype was seen in all populations. Overall linkage disequilibrium is highly significant in all populations, but disequilibrium between the opposite ends is significant only in Native American populations and in one African population. This study demonstrates that the physical extent of linkage disequilibrium can differ substantially among populations from different regions of the world, because of both ancient genetic drift in the ancestor common to a large regional group of modern populations and recent genetic drift affecting individual populations.     

Molecular anatomy of the cytoplasmic domain of bovine growth hormone receptor, a quantitative trait locus

Quantitative trait loci (QTL) studies have indicated growth hormone receptor (GHR) as a candidate gene affecting cattle milk yield and composition. In order to characterize genetic variation at GHR in cattle, we studied European and East African breeds with different histories of selection, and Bos grunniens, Ovis aries, Sus scrofa, Bison bison and Rangifer tarandus as references. We sequenced most of the cytoplasmic domain (900 bp of exon 10), 89 bp of exon 8, including the putative causative mutation for the QTL effect, and 390 bp of intron 8 for comparison. In the cytoplasmic domain, seven synonymous and seven non-synonymous single nucleotide polymorphisms (SNP) were identified in cattle. Three non-synonymous SNPs were found in sheep and one synonymous SNP in yak, while other studied species were monomorphic. Three major haplotypes were observed, one unique to African breeds, one unique to European breeds and one shared. Bison and yak haplotypes are derivatives of the European haplotype lineage. Most of the exon 10 non-synonymous cattle SNPs appear at phylogenetically highly conserved sites. The polymorphisms in exon 10 cluster around a ruminant-specific tyrosine residue, suggesting that this site may act as an additional signalling domain of GHR in ruminants. Alternative explanations for the persistent polymorphism include balancing selection, hitch-hiking, pleiotropic or sexually antagonistic fitness effects or relaxed functional constraints.