Hierarchical patterns of global human Y-chromosome diversity.

We examined 43 biallelic polymorphisms on the nonrecombining portion of the Y chromosome (NRY) in 50 human populations encompassing a total of 2,858 males to study the geographic structure of Y-chromosome variation. Patterns of NRY diversity varied according to geographic region and method/level of comparison. For example, populations from Central Asia had the highest levels of heterozygosity, while African populations exhibited a higher level of mean pairwise differences among haplotypes. At the global level, 36% of the total variance of NRY haplotypes was attributable to differences among populations (i.e., Phi(ST) = 0.36). When a series of AMOVA analyses was performed on different groupings of the 50 populations, high levels of among-groups variance (Phi(CT)) were found between Africans, Native Americans, and a single group containing all 36 remaining populations. The same three population groupings formed distinct clusters in multidimensional scaling plots. A nested cladistic analysis (NCA) demonstrated that both population structure processes (recurrent gene flow restricted by isolation by distance and long-distance dispersals) and population history events (contiguous range expansions and long-distance colonizations) were instrumental in explaining this tripartite division of global NRY diversity. As in our previous analyses of smaller NRY data sets, the NCA detected a global contiguous range expansion out of Africa at the level of the total cladogram. Our new results support a general scenario in which, after an early out-of-Africa range expansion, global-scale patterns of NRY variation were mainly influenced by migrations out of Asia. Two other notable findings of the NCA were (1) Europe as a "receiver" of intercontinental signals primarily from Asia, and (2) the large number of intracontinental signals within Africa. Our AMOVA analyses also supported the hypothesis that patrilocality effects are evident at local and regional scales, rather than at intercontinental and global levels. Finally, our results underscore the importance of subdivision of the human paternal gene pool and imply that caution should be exercised when using models and experimental strategies based on the assumption of panmixia.     

DNA sequence variation in a 3.7-kb noncoding sequence 5' of the CYP1A2 gene: implications for human population history and natural selection

CYP1A2 is a cytochrome P450 gene that is involved in human physiological responses to a variety of drugs and toxins. To investigate the role of population history and natural selection in shaping genetic diversity in CYP1A2, we sequenced a 3.7-kb region 5' from CYP1A2 in a diverse collection of 113 individuals from three major continental regions of the Old World (Africa, Asia, and Europe). We also examined sequences in the 90-member National Institutes of Health DNA Polymorphism Discovery Resource (PDR). Eighteen single-nucleotide polymorphisms (SNPs) were found. Most of the high-frequency SNPs found in the Old World sample were also found in the PDR sample. However, six SNPs were detected in the Old World sample but not in the PDR sample, and two SNPs found in the PDR sample were not found in the Old World sample. Most pairs of SNPs were in complete linkage disequilibrium with one another, and there was no indication of a decline of disequilibrium with physical distance in this region. The average +/- SD nucleotide diversity in the Old World sample was 0.00043+/-0.00026. The African population had the highest level of nucleotide diversity and the lowest level of linkage disequilibrium. Two distinct haplotype clusters with broadly overlapping geographical distributions were present. Of the 17 haplotypes found in the Old World sample, 12 were found in the African sample, 8 were found in Indians, 5 were found in non-Indian Asians, and 5 were found in Europeans. Haplotypes found outside Africa were mostly a subset of those found within Africa. These patterns are all consistent with an African origin of modern humans. Seven SNPs were singletons, and the site-frequency spectrum showed a significant departure from neutral expectations, suggesting population expansion and/or natural selection. Comparison with outgroup species showed that four derived SNPs have achieved high (>0.90) frequencies in human populations, a trend consistent with the action of positive natural selection. These patterns have a number of implications for disease-association studies in CYP1A2 and other genes.     

Separating Population Structure from Population History: A Cladistic Analysis of the Geographical Distribution of Mitochondrial DNA Haplotypes in the Tiger Salamander, Ambystoma Tigrinum

Nonrandom associations of alleles or haplotypes with geographical location can arise from restricted gene flow, historical events (fragmentation, range expansion, colonization), or any mixture of these factors. In this paper, we show how a nested cladistic analysis of geographical distances can be used to test the null hypothesis of no geographical association of haplotypes, test the hypothesis that significant associations are due to restricted gene flow, and identify patterns of significant association that are due to historical events. In this last case, criteria are given to discriminate among contiguous range expansion, long-distance colonization, and population fragmentation. The ability to make these discriminations depends critically upon an adequate geographical sampling design. These points are illustrated with a worked example: mitochondrial DNA haplotypes in the salamander Ambystoma tigrinum. For this example, prior information exists about restricted gene flow and likely historical events, and the nested cladistic analyses were completely concordant with this prior information. This concordance establishes the plausibility of this nested cladistic approach, but much future work will be necessary to demonstrate robustness and to explore the power and accuracy of this procedure.     

The Geographic Distribution of Human Y Chromosome Variation

We examined variation on the nonrecombining portion of the human Y chromosome to investigate human evolution during the last 200,000 years. The Y-specific polymorphic sites included the Y Alu insertional polymorphism or ``YAP' element (DYS287), the poly(A) tail associated with the YAP element, three point mutations in close association with the YAP insertion site, an A-G polymorphic transition (DYS271), and a tetranucleotide microsatellite (DYS19). Global variation at the five bi-allelic sites (DYS271, DYS287, and the three point mutations) gave rise to five ``YAP haplotypes' in 60 populations from Africa, Europe, Asia, Australasia, and the New World (n = 1500). Combining the multi-allelic variation at the microsatellite loci (poly(A) tail and DYS19) with the YAP haplotypes resulted in a total of 27 ``combination haplotypes'. All five of the YAP haplotypes and 21 of the 27 combination haplotypes were found in African populations, which had greater haplotype diversity than did populations from other geographical locations. Only subsets of the five YAP haplotypes were found outside of Africa. Patterns of observed variation were compatible with a variety of hypotheses, including multiple human migrations and range expansions.     

东北亚地区野猪种群mtDNA遗传结构及系统地理发生

研究测定了中国东北、华北及四川西部72个野猪(Susscrofa)个体线粒体控制区全序列,并结合GenBank报道的日本野猪(S.s.leucomystax)、琉球野猪(S.s.riukiuanus)72个同源区序列,分析了东北亚地区野猪线粒体DNA的变异及系统地理格局。在东北亚地区野猪的线粒体控制区共发现42个变异位点,均为转换,共定义了34个单元型。单元型之间的系统发生分析表明,东北亚地区野猪来自同一祖先。东北亚地区野猪现生种群具有显著的种群遗传结构,其中日本野猪与分布于中国东北地区的东北野猪之间亲缘关系较近;而琉球野猪则与华北野猪间亲缘关系较近,与日本野猪和东北野猪间的关系相对较远。嵌套进化枝系统地理分析(Nestedcladephylogeographicalanalysis,NCPA)表明:东北亚地区野猪由同一祖先经过长距离的迁徙而形成现生各种群(或亚种);琉球野猪应起源于大陆野猪,其种群演化可能经历了片断化事件;华北野猪呈现南部种群遗传多样性高的特点,其种群内部曾经历了一次分布区由南向北的扩张     

Phylogeography and recent emergence of the Old World screwworm fly, Chrysomya bezziana, based on mitochondrial and nuclear gene sequences

Abstract.  A previous study had identified an African and an Asian race of the Old World screwworm fly, Chrysomya bezziana Villeneuve (Diptera: Calliphoridae), based on the 3' terminal 279 basepairs (bp) of the mitochondrial cytochrome b gene. The current study improved the phylogeographic resolution of cytochrome b for this species by characterizing more of the gene (the 3' terminal 715 bp) and by sampling more geographical populations, including Oman, Iran, Hong Kong and the Indonesian islands of Sulawesi and East Sumba. Strong support was found for recognizing an African race, but not for a monophyletic Asian race. The cladistic and genealogical relationships among the Asian populations were complex. There was sufficient genetic homogeneity throughout separate regions (mainland Asia and each Indonesian island) to suggest that there are no reproductive barriers within each region that might necessitate the production of more than one strain for control by the sterile insect technique (SIT). Primers were designed for the amplification by polymerase chain reaction of two nuclear loci, the highly conserved elongation factor-1α gene and the less conserved white gene, and the preliminary results indicated that these genes showed the same pattern of small-scale regional variation as cytochrome b . The cytochrome b haplotypes are useful markers for identifying the geographical origins of any emerging infestations of the species: the absence of Indonesian and African haplotypes in the Middle East demonstrates that the large-scale transport of livestock is not spreading Old World screwworm.     

Y‐specific microsatellites reveal an African subfamily in taurine (Bos taurus) cattle

Five cattle Y‐specific microsatellites, totalling six loci, were selected from a set of 44 markers and genotyped on 608 Bos taurus males belonging to 45 cattle populations from Europe and Africa. A total of 38 haplotypes were identified. Haplogroups (Y1 and Y2) previously defined using single nucleotide polymorphisms did not share haplotypes. Nine of the 27 Y2‐haplotypes were only present in African cattle. Network and correspondence analyses showed that this African‐specific subfamily clustered separately from the main Y2‐subfamily and the Y1 haplotypes. Within‐breed genetic variability was generally low, with most breeds (78%) showing haplotypes belonging to a single haplogroup. amova analysis showed that partitioning of genetic variation among breeds can be mainly explained by their geographical and haplogroup assignment. Between‐breed genetic variability summarized via Principal Component Analysis allowed the identification of three principal components explaining 94.2% of the available information. Projection of principal components on geographical maps illustrated that cattle populations located in mainland Europe, the three European Peninsulas and Mediterranean Africa presented similar genetic variation, whereas those breeds from Atlantic Europe and British Islands (mainly carrying Y1 haplotypes) and those from Sub‐Saharan Africa (belonging to Y2‐haplogroup) showed genetic variation of a different origin. Our study confirmed the existence of two large Y‐chromosome lineages (Y1 and Y2) in taurine cattle. However, Y‐specific microsatellites increased analytical resolution and allowed at least two different Y2‐haplotypic subfamilies to be distinguished, one of them restricted to the African continent.     

Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history

Since the 1920s, population geneticists have had measures that describe how genetic variation is distributed spatially within a species' geographical range. Modern genetic survey techniques frequently yield information on the evolutionary relationships among the alleles or haplotypes as well as information on allele frequencies and their spatial distributions. This evolutionary information is often expressed in the form of an estimated haplotype or allele tree. Traditional statistics of population structure, such as F statistics, do not make use of evolutionary genealogical information, so it is necessary to develop new statistical estimators and tests that explicitly incorporate information from the haplotype tree. One such technique is to use the haplotype tree to define a nested series of branches (clades), thereby allowing an evolutionary nested analysis of the spatial distribution of genetic variation. Such a nested analysis can be performed regarding the geographical sampling locations either as categorical or continuous variables (i.e. some measure of spatial distance). It is shown that such nested phylogeographical analyses have more power to detect geographical associations than traditional, nonhistorical analyses and, as a consequence, allow a broader range of gene-flow parameters to be estimated in a precise fashion. More importantly, such nested analyses can discriminate between phylogeographical associations due to recurrent but restricted gene flow vs. historical events operating at the population level (e.g. past fragmentation, colonization, or range expansion events). Restricted gene flow and historical events can be intertwined, and the cladistic analyses can reconstruct their temporal juxtapositions, thereby yielding great insight into both the evolutionary history and population structure of the species. Examples are given that illustrate these properties, concentrating on the detection of range expansion events.     

Global Analysis of ATM Polymorphism Reveals Significant Functional Constraint

ATM, the gene that is mutated in ataxia-telangiectasia, is associated with cerebellar degeneration, abnormal proliferation of small blood vessels, and cancer. These clinically important manifestations have stimulated interest in defining the sequence variation in the ATM gene. Therefore, we undertook a comprehensive survey of sequence variation in ATM in diverse human populations. The protein-encoding exons of the gene (9,168 bp) and the adjacent intron and untranslated sequences (14,661 bp) were analyzed in 93 individuals from seven major human populations. In addition, the coding sequence was analyzed in one chimpanzee, one gorilla, one orangutan, and one Old World monkey. In human ATM, 88 variant sites were discovered by denaturing high-performance liquid chromatography, which is 96%-100% sensitive for detection of DNA sequence variation. ATM was compared to 14 other autosomal genes for nucleotide diversity. The noncoding regions of ATM had diversity values comparable to other genes, but the coding regions had very low diversity, especially in the last 29% of the protein sequence. A test of the neutral evolution hypothesis, through use of the Hudson/Kreitman/Aguadé statistic, revealed that this region of the human ATM gene was significantly constrained relative to that of the orangutan, the Old World monkey, and the mouse, but not relative to that of the chimpanzee or the gorilla. ATM displayed extensive linkage disequilibrium, consistent with suppression of meiotic recombination at this locus. Seven haplotypes were defined. Two haplotypes accounted for 82% of all chromosomes analyzed in all major populations; two others carrying the same D126E missense polymorphism accounted for 33% of chromosomes in Africa but were never observed outside of Africa. The high frequency of this polymorphism may be due either to a population expansion within Africa or to selective pressure.     

Recombination networks as genetic markers in a human variation study of the Old World

We have analyzed human genetic diversity in 33 Old World populations including 23 populations obtained through Genographic Project studies. A set of 1,536 SNPs in five X chromosome regions were genotyped in 1,288 individuals (mostly males). We use a novel analysis employing subARG network construction with recombining chromosomal segments. Here, a subARG is constructed independently for each of five gene-free regions across the X chromosome, and the results are aggregated across them. For PCA, MDS and ancestry inference with STRUCTURE, the subARG is processed to obtain feature vectors of samples and pairwise distances between samples. The observed population structure, estimated from the five short X chromosomal segments, supports genome-wide frequency-based analyses: African populations show higher genetic diversity, and the general trend of shared variation is seen across the globe from Africa through Middle East, Europe, Central Asia, Southeast Asia, and East Asia in broad patterns. The recombinational analysis was also compared with established methods based on SNPs and haplotypes. For haplotypes, we also employed a fixed-length approach based on information-content optimization. Our recombinational analysis suggested a southern migration route out of Africa, and it also supports a single, rapid human expansion from Africa to East Asia through South Asia.     

Ancestral Asian source(s) of new world Y-chromosome founder haplotypes

Haplotypes constructed from Y-chromosome markers were used to trace the origins of Native Americans. Our sample consisted of 2,198 males from 60 global populations, including 19 Native American and 15 indigenous North Asian groups. A set of 12 biallelic polymorphisms gave rise to 14 unique Y-chromosome haplotypes that were unevenly distributed among the populations. Combining multiallelic variation at two Y-linked microsatellites (DYS19 and DXYS156Y) with the unique haplotypes results in a total of 95 combination haplotypes. Contra previous findings based on Y- chromosome data, our new results suggest the possibility of more than one Native American paternal founder haplotype. We postulate that, of the nine unique haplotypes found in Native Americans, haplotypes 1C and 1F are the best candidates for major New World founder haplotypes, whereas haplotypes 1B, 1I, and 1U may either be founder haplotypes and/or have arrived in the New World via recent admixture. Two of the other four haplotypes (YAP+ haplotypes 4 and 5) are probably present because of post-Columbian admixture, whereas haplotype 1G may have originated in the New World, and the Old World source of the final New World haplotype (1D) remains unresolved. The contrasting distribution patterns of the two major candidate founder haplotypes in Asia and the New World, as well as the results of a nested cladistic analysis, suggest the possibility of more than one paternal migration from the general region of Lake Baikal to the Americas.     

Spatial and temporal distribution of the neutral polymorphisms in the last ZFX intron: analysis of the haplotype structure and genealogy.

With 10 segregating sites (simple nucleotide polymorphisms) in the last intron (1089 bp) of the ZFX gene we have observed 11 haplotypes in 336 chromosomes representing a worldwide array of 15 human populations. Two haplotypes representing 77% of all chromosomes were distributed almost evenly among four continents. Five of the remaining haplotypes were detected in Africa and 4 others were restricted to Eurasia and the Americas. Using the information about the ancestral state of the segregating positions (inferred from human-great ape comparisons), we applied coalescent analysis to estimate the age of the polymorphisms and the resulting haplotypes. The oldest haplotype, with the ancestral alleles at all the sites, was observed at low frequency only in two groups of African origin. Its estimated age of 740 to 1100 kyr corresponded to the time to the most recent common ancestor. The two most frequent worldwide distributed haplotypes were estimated at 550 to 840 and 260 to 400 kyr, respectively, while the age of the continentally restricted polymorphisms was 120 to 180 kyr and smaller. Comparison of spatial and temporal distribution of the ZFX haplotypes suggests that modern humans diverged from the common ancestral stock in the Middle Paleolithic era. Subsequent range expansion prevented substantial gene flow among continents, separating African groups from populations that colonized Eurasia and the New World.     

Population Structure and the Colonization Route of One of the Oldest North American Invasive Insects: Stories from the Worn Road of the Hessian Fly,Mayetiola destructor (Say)

An integral part to understanding the biology of an invasive species is determining its origin, particularly in pest species. As one of the oldest known invasive species, the goals of this study were to evaluate the evidence of a westward expansion of Hessian fly into North America, from a potential singular introduction event, and the population genetic structure of current populations. Levels of genetic diversity and population structure in the Hessian fly were compared across North America, Europe, North Africa, Western Asia, and New Zealand. Furthermore, Old World populations were evaluated as possible sources of introduction. We tested diversity and population structure by examining 18 microsatellite loci with coverage across all four Hessian fly chromosomes. Neither genetic diversity nor population genetic structure provided evidence of a westward movement from a single introduction in North America. Introduced populations in North America did not show identity or assignment to any Old World population, likely indicating a multiple introduction scenario with subsequent gene flow between populations. Diversity and selection were assessed on a chromosomal level, with no differences in diversity or selection between chromosomes or between native and introduced populations.     

Phylogeny of human beta-globin haplotypes and its implications for recent human evolution

The evolutionary histories and relationships among African, Eurasian, and Pacific Island populations are investigated by using observations on five polymorphic restriction sites in the beta-globin gene cluster. We present new data on 222 chromosomes from a global sample and combine these with previously published observations on 591 chromosomes. It is shown that the data are rich in rare haplotypes and that rare variants are not helpful for standard methods of population structure analysis. Consequently, a new approach is developed. We first consider the phylogeny of beta-globin haplotypes. The roles of mutation, gene conversion, and recombination in the generation of haplotype diversity are specifically focused upon. The relationships among human populations are then inferred from the phylogenetic relationships among the haplotypes, their presence or absence, and frequencies within populations. Questions regarding whether or not a phyletic process can account for relationships among the major geographical populations and whether or not an extant human population exhibits the qualities that would be expected of an ancestral group are addressed. The results of this analysis support an African origin for modern Homo sapiens and a phyletic structuring of the major geographical regions. However, it is shown that divergence times for the various populations cannot be determined from these data.     

Spatial diversity in mitochondrial cytochrome c oxidase in house flies

DNA sequence analysis at mitochondrial gene COI was surveyed in 293 house flies, Musca domestica Linneaus (Diptera: Muscidae), in 29 populations from North, Central and South America, Europe, Asia, Africa, and the Western Pacific. Nei's gene diversity index (H(S)) was 0.47, the chance that two randomly chosen flies have different COI haplotypes. Haplotype diversity was greater in the Old World (H(S) = 0.58) than the New World (H(S) = 0.31). The hierarchical partition of the total diversity indicated substantial differentiation at all levels (G(ST) = 0.30), and highly structured populations. All pairwise estimates of gene flow between zoogeographical regions were less than 0.70 reproducing females per generation. The results are compared to those of a similar study based on the single-strand conformation polymorphism method. Probable colonization scenarios for house flies into the New World are discussed and it is concluded that house flies are a recent addition to the fauna of the Western Hemisphere.     

Geographical survey of βS-globin gene haplotypes: Evidence for an independent Asian origin of the sickle-cell mutation

The haplotypes of 152 beta S-chromosomes were characterized in six different population groups. The chromosomes of individuals from Nigeria and from the southwest of the Arabian peninsula have the haplotype - - - - + + - + previously found in west African, Jamaican, and U.S. American blacks, whereas those from the eastern oases of Saudi Arabia and from the west and the east coast of India showed a different haplotype not found in Africa (+ + - + + + + -). These data are most consistent with an independent Asian origin of the sickle-cell mutation and provide further information about the geographic distribution of beta S-haplotypes in the Old World. The distribution of the Asian beta S-haplotype corresponds to the reported geographic distribution of a mild clinical phenotype of homozygous SS disease.     

Phylogenetic status and matrilineal structure of the biting midge, Culicoides imicola, in Portugal, Rhodes and Israel

The biting midge Culicoides imicola Kieffer (Diptera: Ceratopogonidae) is the most important Old World vector of African horse sickness (AHS) and bluetongue (BT). Recent increases of BT incidence in the Mediterranean basin are attributed to its increased abundance and distribution. The phylogenetic status and genetic structure of C. imicola in this region are unknown, despite the importance of these aspects for BT epidemiology in the North American BT vector. In this study, analyses of partial mitochondrial cytochrome oxidase subunit I gene (COI) sequences were used to infer phylogenetic relationships among 50 C. imicola from Portugal, Rhodes, Israel, and South Africa and four other species of the Imicola Complex from southern Africa, and to estimate levels of matrilineal subdivision in C. imicola between Portugal and Israel. Eleven haplotypes were detected in C. imicola, and these formed one well-supported clade in maximum likelihood and Bayesian trees implying that the C. imicola samples comprise one phylogenetic species. Molecular variance was distributed mainly between Portugal and Israel, with no haplotypes shared between these countries, suggesting that female-mediated gene flow at this scale has been either limited or non-existent. Our results provide phylogenetic evidence that C. imicola in the study areas are potentially competent AHS and BT vectors. The geographical structure of the C. imicola COI haplotypes was concordant with that of BT virus serotypes in recent BT outbreaks in the Mediterranean basin, suggesting that population subdivision in its vector can impose spatial constraints on BT virus transmission.     

High-resolution analysis of human Y-chromosome variation shows a sharp discontinuity and limited gene flow between northwestern Africa and the Iberian Peninsula

In the present study we have analyzed 44 Y-chromosome biallelic polymorphisms in population samples from northwestern (NW) Africa and the Iberian Peninsula, which allowed us to place each chromosome unequivocally in a phylogenetic tree based on >150 polymorphisms. The most striking results are that contemporary NW African and Iberian populations were found to have originated from distinctly different patrilineages and that the Strait of Gibraltar seems to have acted as a strong (although not complete) barrier to gene flow. In NW African populations, an Upper Paleolithic colonization that probably had its origin in eastern Africa contributed 75% of the current gene pool. In comparison, approximately 78% of contemporary Iberian Y chromosomes originated in an Upper Paleolithic expansion from western Asia, along the northern rim of the Mediterranean basin. Smaller contributions to these gene pools (constituting 13% of Y chromosomes in NW Africa and 10% of Y chromosomes in Iberia) came from the Middle East during the Neolithic and, during subsequent gene flow, from Sub-Saharan to NW Africa. Finally, bidirectional gene flow across the Strait of Gibraltar has been detected: the genetic contribution of European Y chromosomes to the NW African gene pool is estimated at 4%, and NW African populations may have contributed 7% of Iberian Y chromosomes. The Islamic rule of Spain, which began in a.d. 711 and lasted almost 8 centuries, left only a minor contribution to the current Iberian Y-chromosome pool. The high-resolution analysis of the Y chromosome allows us to separate successive migratory components and to precisely quantify each historical layer.     

A back migration from Asia to sub-Saharan Africa is supported by high-resolution analysis of human Y-chromosome haplotypes

The variation of 77 biallelic sites located in the nonrecombining portion of the Y chromosome was examined in 608 male subjects from 22 African populations. This survey revealed a total of 37 binary haplotypes, which were combined with microsatellite polymorphism data to evaluate internal diversities and to estimate coalescence ages of the binary haplotypes. The majority of binary haplotypes showed a nonuniform distribution across the continent. Analysis of molecular variance detected a high level of interpopulation diversity (PhiST=0.342), which appears to be partially related to the geography (PhiCT=0.230). In sub-Saharan Africa, the recent spread of a set of haplotypes partially erased pre-existing diversity, but a high level of population (PhiST=0.332) and geographic (PhiCT=0.179) structuring persists. Correspondence analysis shows that three main clusters of populations can be identified: northern, eastern, and sub-Saharan Africans. Among the latter, the Khoisan, the Pygmies, and the northern Cameroonians are clearly distinct from a tight cluster formed by the Niger-Congo-speaking populations from western, central western, and southern Africa. Phylogeographic analyses suggest that a large component of the present Khoisan gene pool is eastern African in origin and that Asia was the source of a back migration to sub-Saharan Africa. Haplogroup IX Y chromosomes appear to have been involved in such a migration, the traces of which can now be observed mostly in northern Cameroon.     

Extent of linkage disequilibrium between the androgen receptor gene CAG and GGC repeats in human populations: implications for prostate cancer risk

While studies have implicated alleles at the CAG and GGC trinucleotide repeats of the androgen receptor gene with high-grade, aggressive prostate cancer disease, little is known about the normal range of variation for these two loci, which are separated by about 1.1 kb. More importantly, few data exist on the extent of linkage disequilibrium (LD) between the two loci in different human populations. Here we present data on CAG and GGC allelic variation and LD in six diverse populations. Alleles at the CAG and GGC repeat loci of the androgen receptor were typed in over 1000 chromosomes from Africa, Asia, and North America. Levels of linkage disequilibrium between the two loci were compared between populations. Haplotype variation and diversity were estimated for each population. Our results reveal that populations of African descent possess significantly shorter alleles for the two loci than non-African populations (P<0.0001). Allelic diversity for both markers was higher among African Americans than any other population, including indigenous Africans from Sierra Leone and Nigeria. Analysis of molecular variance revealed that approx. 20% of CAG and GGC repeat variance could be attributed to differences between the populations. All non-African populations possessed the same common haplotype while the three populations of African descent possessed three divergent common haplotypes. Significant LD was observed in our sample of healthy African Americans. The LD observed in the African American population may be due to several reasons; recent migration of African Americans from diverse rural communities following urbanization, recurrent gene flow from diverse West African populations, and admixture with European Americans. This study represents the largest genotyping effort to be performed on the two androgen receptor trinucleotide repeat loci in diverse human populations.