A global survey of haplotype frequencies and linkage disequilibrium at the DRD2 locus

Distinct linkage disequilibrium (LD) runs of single nucleotide polymorphisms and microsatellite markers; implications for use of mixed marker haplotypes in LD-based mapping

It has been suggested that the haplotypic relationship between microsatellite markers and single nucleotide polymorphisms (SNPs) is of considerable importance, as microsatellite markers can potentially be incorporated into haplotypes containing SNPs to increase marker density across a region of interest. However, SNPs and microsatellite markers have different mutation rates and durations, and it is conceivable that the linkage disequilibrium (LD) patterns between the genetic markers may considerably differ. We assessed the LD patterns using 1,661 SNPs and 65 microsatellite markers along chromosome 22 and investigated whether common patterns of LD between the two genetic markers are deduced from the results. The results demonstrated that the patterns of LD among microsatellite markers varied considerably and the LD runs of SNPs and microsatellite markers showed distinct patterns. Microsatellite markers have a much higher mutation rate and the evolution of microsatellite markers is a more complex process which has distinct mutation properties from those of SNPs. We consider that these might contribute to the different LD patterns between the two genetic markers. Therefore, it would seem inadvisable to make assumptions about persistence of LD across even a relatively small genetic distance among microsatellite markers and to construct mixed marker haplotypes/LD maps employing microsatellite markers.

Overview of the pharmacogenomics of cigarette smoking

Cigarette smoking increases the risk of numerous health problems, including cancer, cardiovascular and pulmonary disorders, making smoking the leading cause of preventable death in the world. Nicotine is primarily responsible for the highly addictive properties of cigarettes. Although the majority of smokers express a desire to quit, few are successful in doing so. Twin and family studies have indicated substantial genetic contributions to smoking behaviors. One major research focus has been to elucidate the specific genes involved; this has been accomplished primarily through genome-wide linkage analyses and candidate gene association studies. Much attention has focused on genes involved in the neurotransmitter pathways for the brain reward system and genes altering nicotine metabolism. This paper reviews the current state of knowledge for genetic factors implicated in smoking behaviors, and examines how genetic variations may affect therapeutic outcomes for drugs used to assist smoking cessation.

Significant variation in haplotype block structure but conservation in tagSNP patterns among global populations

The initial belief that haplotype block boundaries and haplotypes were largely shared across populations was a foundation for constructing a haplotype map of the human genome using common SNP markers. The HapMap data document the generality of a block-like pattern of linkage disequilibrium (LD) with regions of low and high haplotype diversity but differences among the populations. Studies of many additional populations demonstrate that LD patterns can be highly variable among populations both across and within geographic regions. Because of this variation, emphasis has shifted to the generalizability of tagSNPs, those SNPs that capture the bulk of variation in a region. We have examined the LD and tagSNP patterns based upon over 2000 individual samples in 38 populations and 134 SNPs in 10 genetically independent loci for a total of 517 kb with an average density of 1 SNP/5 kb. Four different 'block' definitions and the pairwise LD tagSNP selection algorithm have been applied. Our results not only confirm large variation in block partition among populations from different regions (agreeing with previous studies including the HapMap) but also show that significant variation can occur among populations within geographic regions. None of the block-defining algorithms produces a consistent pattern within or across all geographic groups. In contrast, tagSNP transferability is much greater than the similarity of LD patterns and, although not perfect, some generalizations of transferability are possible. The analyses show an asymmetric pattern of tagSNP transferability coinciding with the subsetting of variation attributed to the spread of modern humans around the world.

Nonsynonymous SNPs: validation characteristics, derived allele frequency patterns, and suggestive evidence for natural selection

We experimentally investigated more than 1,200 entries in dbSNP that would change amino-acids (nsSNPs), using various subsets of DNA samples drawn from 18 global populations (approximately 1,000 subjects in total). First, we mined the data for any SNP features that correlated with a high validation rate. Useful predictors of valid SNPs included multiple submissions to dbSNP, having a dbSNP validation statement, and being present in a low number of ESTs. Together, these features improved validation rates by almost 10-fold. Higher-abundance SNPs (e.g., T/C variants) also validated more frequently. Second, we considered derived alleles and noted a considerably (approximately 10%) increased average derived allele frequency (DAF) in Europeans vs. Africans, plus a further increase in some other populations. This was not primarily due to an SNP ascertainment bias, nor to the effects of natural selection. Instead, it can be explained as a drift-based, progressive increase in DAF that occurs over many generations and becomes exaggerated during population bottlenecks. This observation could be used as the basis for novel DAF-based tests for comparing demographic histories. Finally, we considered individual marker patterns and identified 37 SNPs with allele frequency variance or FST values consistent with the effects of population-specific natural selection. Four particularly striking clusters of these markers were apparent, and three of these coincide with genes/regions from among only several dozen such domains previously suggested by others to carry signatures of selection.

Variability in TRBV haplotype frequency and composition in Caucasian, African American, Western African and Chinese populations

The polymorphic T-cell receptor Vbeta (TRBV) genes encode much of the variable region of the T-cell receptor beta chain. Analysis of allele frequencies of three closely linked polymorphic TRBV genes, TRBV7-3, TRBV9 and TRBV6-4, was undertaken in several populations. The frequencies of these alleles are not significantly different in populations of Caucasians, African Americans and Western Africans. However, Chinese population is extremely homogenous at all three loci. The current study identifies the existence of haplotypic relationships between alleles of these genes in the Caucasian population. The ORF allele TRBV7-3*A3 is found exclusively on chromosomes bearing TRBV9*A2 and TRBV6-4*A2 in this cohort. In contrast, TRBV7-3*A1 and the null allele TRBV7-3*A2 are associated only with TRBV9*A1 and TRBV6-4*A1. This pattern of linkage disequilibrium (LD) is altered in the African American and Western African populations. In these cohorts, there is a marked reduction in LD between alleles of TRBV7-3 and TRBV9. This study is consistent with previous population genetic studies wherein African-derived samples have a greater level of genetic diversity compared to Caucasians. These data also demonstrate that patterns of LD are not consistent across the entire TRBV locus.

Extreme population differences across Neuregulin 1 gene, with implications for association studies

Neuregulin 1 (NRG1) is one of the most exciting candidate genes for schizophrenia in recent years since its first association with the disease in an Icelandic population. Since then, many association studies have analysed allele and haplotype frequencies in distinct populations yielding varying results: some have replicated the association, although with different alleles or haplotypes being associated, whereas others have failed to replicate the association. These contradictory results might be attributed to population differences in allele and haplotype frequencies. In order to approach this issue, we have typed 13 SNPs across this large 1.4 Mb gene, including two of the SNPs originally found associated with schizophrenia in the Icelandic population, the objective being to discover if the underlying cause of the association discrepancies to date may be due to population-specific genetic variation. The analyses have been performed in a total of 1088 individuals from 39 populations, covering most of the genetic diversity in the human species. Most of the SNPs analysed displayed differing frequencies according to geographical region. These allele differences are especially relevant in two SNPs located in a large intron of the gene, as shown by the extreme F(ST) values, which reveal genetic stratification correlated to broad continental areas. This finding may be indicative of the influence of some local selective forces on this gene. Furthermore, haplotype analysis reveals a clear clustering according to geographical areas. In summary, our findings suggest that NRG1 presents extreme population differences in allele and haplotype frequencies. We have given recommendations for taking this into account in future association studies since this diversity could give rise to erroneous results.

Prevalence of dopamine and 5HT2C receptor polymorphisms in Amerindians and in an urban population from Argentina

We report on the frequency of DRD2A1, DRD3A1, DRD4/2R-10R, and 5HT2CA1 variants in the population of the city of La Plata (Argentina) and in Amerindians from Argentina, Paraguay, and Chile. In the Amerindian sample, the prevalence of DRD2A1 and DRD4/4R variants were, respectively, significantly lower and significantly higher than frequencies reported in other Native Americans. Comparison of average allele and genotype frequencies between La Plata and Amerindians showed significant differences for 5HT2CA1 and DRD4. As La Plata is a population with predominant European and Amerindian components, we used mtDNA and Y-specific markers to subdivide the La Plata sample into two strata: Amerindian La Plata and non-Amerindian La Plata. Significant variations between the two strata were detected for DRD2A1, DRD3A1, and DRD4/4R allele frequencies, and for the homozygous DRD4/4R/4R genotype. Several controversial reports suggest a possible association between a variant of DRD and/or 5HT2C receptor genes and the clinical expression of several psychiatric disorders. We suggest that ethnic variations in the prevalence of the allelic forms of these genes may be a confounding factor to be taken into consideration in studies of association between dopaminergic and serotonergic receptor genotypes and neuropsychiatric and mood disorders.

Association of multiple DRD2 polymorphisms with anorexia nervosa

To investigate whether the dopaminergic system plays a role in the etiology of anorexia nervosa (AN) via the dopamine D2 receptor, we investigated association and transmission disequilibrium at seven single-nucleotide polymorphisms (SNPs) spanning about 75 kbp of the gene DRD2. We studied 191 probands with a DSM-IV diagnosis of AN, 457 parents and affected relatives with a DSM-IV eating disorder diagnosis, and 98 unrelated, female, normal weight controls. The -141 C/- insertion/deletion (-141 Indel), previously shown to affect DRD2 transcription efficiency, and multiple exon seven polymorphisms, one of which has previously been shown to affect DRD2 transcript stability, exhibited statistically significant association with diagnosis in haplotype transmission disequilibrium and in haplotype case : control analyses. Significant linkage disequilibrium between the -141 Indel and two exon seven SNPs (939Y and 957Y) was observed over a distance of >50 kbp in the AN probands but not in the controls. Genetically transmitted variation in D2 dopamine receptor expression mediated by functional polymorphisms affecting transcription and translation efficiency may play a role in vulnerability to AN.

Linkage disequilibrium patterns vary substantially among populations

A major initiative to create a global human haplotype map has recently been launched as a tool to improve the efficiency of disease gene mapping. The 'HapMap' project will study common variants in depth in four (and to a lesser degree in up to 12) populations to catalogue haplotypes that are expected to be common to all populations. A hope of the 'HapMap' project is that much of the genome occurs in regions of limited diversity such that only a few of the SNPs in each region will capture the diversity and be relevant around the world. In order to explore the implications of studying only a limited number of populations, we have analyzed linkage disequilibrium (LD) patterns of three 175-320 kb genomic regions in 16 diverse populations with an emphasis on African and European populations. Analyses of these three genomic regions provide empiric demonstration of marked differences in frequencies of the same few haplotypes, resulting in differences in the amount of LD and very different sets of haplotype frequencies. These results highlight the distinction between the statistical concept of LD and the biological reality of haplotypes and their frequencies. The significant quantitative and qualitative variation in LD among populations, even for populations within a geographic region, emphasizes the importance of studying diverse populations in the HapMap project to assure broad applicability of the results.

Haplotype tagging efficiency in worldwide populations in CTLA4 gene

The cytotoxic T lymphocyte antigen 4 (CTLA4) acts as a potent negative regulator of T-cell response, and has been suggested as a pivotal candidate gene for autoimmune disorders such as Graves' disease, type 1 diabetes and autoimmune hypothyroidism, among others. Several single-nucleotide polymorphisms (SNPs) have been proposed as the susceptibility variants, or to be in strong linkage disequilibrium (LD) with the variant. Nevertheless, contradictory results have been found, which may be due to lack of knowledge of the genetic structure of CTLA4 and its geographic variation. We have typed 17 SNPs throughout the CTLA4 gene region in order to analyze the haplotype diversity and LD structure in a worldwide population set (1262 individuals from 44 populations) to understand the variation pattern of the region. Allele and haplotype frequency differentiation between populations is consistent with genomewide averages and points to a lack of strong population-specific selection pressures. LD is high and its pattern is not significantly different within or between continents. However, haplotype composition is significantly different between geographical groups. A continent-specific set of haplotype tagging SNPs has been designed to be used for future association studies. These are portable among populations, although their efficiency might vary depending on the population haplotype spectrum.

Genetic structure in four West African population groups

BACKGROUND

Africa contains the most genetically divergent group of continental populations and several studies have reported that African populations show a high degree of population stratification. In this regard, it is important to investigate the potential for population genetic structure or stratification in genetic epidemiology studies involving multiple African populations. The presences of genetic sub-structure, if not properly accounted for, have been reported to lead to spurious association between a putative risk allele and a disease. Within the context of the Africa America Diabetes Mellitus (AADM) Study (a genetic epidemiologic study of type 2 diabetes mellitus in West Africa), we have investigated population structure or stratification in four ethnic groups in two countries (Akan and Gaa-Adangbe from Ghana, Yoruba and Igbo from Nigeria) using data from 372 autosomal microsatellite loci typed in 493 unrelated persons (986 chromosomes).

RESULTS

There was no significant population genetic structure in the overall sample. The smallest probability is associated with an inferred cluster of 1 and little of the posterior probability is associated with a higher number of inferred clusters. The distribution of members of the sample to inferred clusters is consistent with this finding; roughly the same proportion of individuals from each group is assigned to each cluster with little variation between the ethnic groups. Analysis of molecular variance (AMOVA) showed that the between-population component of genetic variance is less than 0.1% in contrast to 99.91% for the within population component. Pair-wise genetic distances between the four ethnic groups were also very similar. Nonetheless, the small between-population genetic variance was sufficient to distinguish the two Ghanaian groups from the two Nigerian groups.

CONCLUSION

There was little evidence for significant population substructure in the four major West African ethnic groups represented in the AADM study sample. Ethnicity apparently did not introduce differential allele frequencies that may affect analysis and interpretation of linkage and association studies. These findings, although not entirely surprising given the geographical proximity of these groups, provide important insights into the genetic relationships between the ethnic groups studied and confirm previous results that showed close genetic relationship between most studied West African groups.

Evolution of the DRD2 gene haplotype and its association with alcoholism in Mexican Americans

The human D2 dopamine receptor gene (DRD2) plays a central role in the neuromodulation of appetitive behaviors and is implicated in having a possible role in susceptibility to alcoholism. We genotyped an SNP in DRD2 Exon 8 in 251 nonalcoholic, unrelated, healthy controls and 200 alcoholic Mexican Americans. The DRD2 haplotypes were analyzed using the Exon 8 genotype in combination with five other SNP genotypes, which were obtained from our previous study. The ancestral origins of the DRD2 polymorphisms have been determined by sequencing the homologous region in other higher primates. Twenty DRD2 haplotypes, defined as H1 to H20 based on their frequency from high to low, were obtained in this major minority population. The ancestral haplotype "I-B2-G-C-G-A1" and two one-step mutation haplotypes were absent in our study population. The haplotype H1, "I-B1-T-C-A-A1", with the highest frequency in the population, is a three-step mutation from the ancestral form. The first five or eight major haplotypes make up 87% or 95% of the entire population, respectively. The prevalence of the haplotype H1+ (H1/H1 and H1/Hn genotypes) is significantly higher in alcoholics and alcoholic subgroups, including early onset drinkers and benders, than in their respective control groups. The Promoter -141C allele is in linkage disequilibrium (LD) with five other loci in the nonalcoholic group, but not in the alcoholic group. All of the other five loci are in LD in both the alcoholic and control groups. The DRD2 TaqI B allele is in complete LD with the allele located in intron 6. Five SNPs, Promoter -141C, TaqI B (or Intron 6), Exon 7, Exon 8, and TaqI A, are sufficient to define the DRD2 haplotypes in Mexican Americans. Our data indicate that the DRD2 haplotypes are associated with alcoholism in Mexican Americans.

A comparison of linkage disequilibrium patterns and estimated population recombination rates across multiple populations

Large-scale studies of linkage disequilibrium (LD) have shown considerable variation in the extent and distribution of pairwise LD within and between populations. Taken at face value, these results suggest that genomewide LD maps for one population may not be generalizable to other populations. However, at least part of this diversity is due to some undesirable features of pairwise LD measures, which are well documented for the D' and r2 measures. In this report, we compare patterns of LD derived from pairwise measures with statistical estimates of population recombination rates ( rho ) along a 10-Mb stretch of chromosome 20 in four population samples, comprising East Asians, African Americans, and U.K. and U.S. individuals of western European descent. The results reveal the expected variability of D' within and between populations but show better concordance in estimates of r2 for the same markers across the population samples. Estimates of rho correlate well across populations, but there is still evidence of population-specific spikes and troughs in rho values. We conclude that it is unlikely that a single haplotype map will provide a definitive guide for association studies of many populations; rather, multiple maps will need to be constructed to provide the best-possible guides for gene mapping.

Conservation of theRB1gene in human and primates

Mutations in the RB1 gene are associated with retinoblastoma, which has served as an important model for understanding hereditary predisposition to cancer. Despite the great scrutiny that RB1 has enjoyed as the prototypical tumor suppressor gene, it has never been the object of a comprehensive survey of sequence variation in diverse human populations and primates. Therefore, we analyzed the coding (2,787 bp) and adjacent intronic and untranslated (7,313 bp) sequences of RB1 in 137 individuals from a wide range of ethnicities, including 19 Asian Indian hereditary retinoblastoma cases, and five primate species. Aside from nine apparently disease-associated mutations, 52 variants were identified. They included six singleton, coding variants that comprised five amino acid replacements and one silent site. Nucleotide diversity of the coding region (pi=0.0763+/-1.35 x 10(-4)) was 52 times lower than that of the noncoding regions (pi=3.93+/-5.26 x 10(-4)), indicative of significant sequence conservation. The occurrence of purifying selection was corroborated by phylogeny-based maximum likelihood analysis of the RB1 sequences of human and five primates, which yielded an estimated ratio of replacement to silent substitutions (omega) of 0.095 across all lineages. RB1 displayed extensive linkage disequilibrium over 174 kb, and only four unique recombination events, two in Africa and one each in Europe and Southwest Asia, were observed. Using a parsimony approach, 15 haplotypes could be inferred. Ten were found in Africa, though only 12.4% of the 274 chromosomes screened were of African origin. In non-Africans, a single haplotype accounted for from 63 to 84% of all chromosomes, most likely the consequence of natural selection and a significant bottleneck in effective population size during the colonization of the non-African continents.

The linkage disequilibrium maps of three human chromosomes across four populations reflect their demographic history and a common underlying recombination pattern

The extent and patterns of linkage disequilibrium (LD) determine the feasibility of association studies to map genes that underlie complex traits. Here we present a comparison of the patterns of LD across four major human populations (African-American, Caucasian, Chinese, and Japanese) with a high-resolution single-nucleotide polymorphism (SNP) map covering almost the entire length of chromosomes 6, 21, and 22. We constructed metric LD maps formulated such that the units measure the extent of useful LD for association mapping. LD reaches almost twice as far in chromosome 6 as in chromosomes 21 or 22, in agreement with their differences in recombination rates. By all measures used, out-of-Africa populations showed over a third more LD than African-Americans, highlighting the role of the population's demography in shaping the patterns of LD. Despite those differences, the long-range contour of the LD maps is remarkably similar across the four populations, presumably reflecting common localization of recombination hot spots. Our results have practical implications for the rational design and selection of SNPs for disease association studies.

The dopamine D2 receptor C32806T polymorphism (DRD2 Taq1A RFLP) exhibits no association with smoking behaviour in a healthy UK population

A single nucleotide polymorphism (SNP) in the Taq1A site near the DRD2 gene has been associated in several studies with smoking behaviour. We genotyped 732 current smokers (241 low, one to nine cigarettes a day, 250 mid, 10-19 cigarettes, 241 high, 20+cigarettes) and 243 never-smokers at this site (C32806T), to test for effects on smoking initiation and amount of tobacco consumed. No significant association between minor allele frequency and smoking status was detected. Multiple regression analysis including DRD2 genotype, sex, age and alcohol consumption as predictors showed that level of cigarette consumption was associated with sex (p=0.003) and age (p=0.002) but not with alcohol consumption (p=0.25) or DRD2 genotype (p=0.76).

Evaluation of microsatellite markers in association studies: a search for an immune-related susceptibility gene in sarcoidosis

Association studies using linkage disequilibrium (LD) between candidate loci and nearby markers have been proposed to identify susceptibility genes for complex diseases. We analyzed polymorphisms of microsatellites (MSs) and LD patterns of the regions in which candidate genes related to the Th1 immune response have been annotated and attempted to identify a susceptibility gene for sarcoidosis in a marker-based association study. Nineteen MSs were identified in six Th1-related genes (IFNGR1, IFNGR2, IL12RB1, IL12RB2, STAT1 and STAT4) and then eight were further characterized as useful polymorphic markers. Most of these MSs showed LD with single nucleotide polymorphisms (SNPs) on both 5' and 3' ends of these candidate genes, in which r(2) values between at least one of the MS marker alleles and the SNPs were higher than 0.1. A significant association with one MS allele near STAT4 was shown and a cluster of SNPs in LD with the MS marker was associated with sarcoidosis. These results suggest that association studies using not only SNPs but also multi-allelic MS within or near candidate loci would be useful markers to search for a disease susceptibility gene, especially in populations with unknown LD structure.

p53 polymorphisms in Russia and Belarus: correlation of the 2-1-1 haplotype frequency with longitude

Four different polymorphisms in the human p53 gene (a 16-bp duplication in intron 3, and three RFLPs: for Bsh1236I at codon 72, for MspI in intron 6 and for BamHI in the 3' flanking region) and extended haplotypes were studied in nine geographically diverse populations from Russia and Belarus. The Yakuts differed from all other populations, as they had a significantly higher frequency of the BamHI A1 allele. Most populations did not differ significantly from each other in the frequency of the Bsh1236I polymorphism. The 16-bp duplication A1 allele and MspI A2 allele frequencies were significantly higher in the Yakut and Khant populations. Linkage disequilibrium values (D') between BamHI and other polymorphic sites were not significant in many cases; for this reason we have used the 16 bp-Bsh1236I-MspI haplotype frequencies only. Of eight possible haplotypes, five were observed in the populations investigated. Haplotype 1-2-2 was the most frequent in all populations. The next most common haplotype, 1-1-2, was present at very similar frequencies among the Byelorussians and Russians from Smolensk, but was more frequent in other populations. The frequency of haplotype 2-1-1 showed a nearly continuous decrease from West to East (from 17.857% among the Byelorussians to 0.685% in the Yakuts from the Verkhoyansk) and correlated with longitude (Spearman's r = -0.8667, P = 0.0025), which may be due to natural selection and adaptation. The relationships among populations were evaluated by means of Nei's D(A) distances for the 16 bp-Bsh1236I-MspI haplotype frequencies. Based on the multidimensional scaling analysis a correlation between p53 haplotype frequencies and ethnicity is supposed.

Gametic phase estimation over large genomic regions using an adaptive window approach

The authors present ELB, an easy to programme and computationally fast algorithm for inferring gametic phase in population samples of multilocus genotypes. Phase updates are made on the basis of a window of neighbouring loci, and the window size varies according to the local level of linkage disequilibrium. Thus, ELB is particularly well suited to problems involving many loci and/or relatively large genomic regions, including those with variable recombination rate. The authors have simulated population samples of single nucleotide polymorphism genotypes with varying levels of recombination and marker density, and find that ELB provides better local estimation of gametic phase than the PHASE or HTYPER programs, while its global accuracy is broadly similar. The relative improvement in local accuracy increases both with increasing recombination and with increasing marker density. Short tandem repeat (STR, or microsatellite) simulation studies demonstrate ELB's superiority over PHASE both globally and locally. Missing data are handled by ELB; simulations show that phase recovery is virtually unaffected by up to 2 per cent of missing data, but that phase estimation is noticeably impaired beyond this amount. The authors also applied ELB to datasets obtained from random pairings of 42 human X chromosomes typed at 97 diallelic markers in a 200 kb low-recombination region. Once again, they found ELB to have consistently better local accuracy than PHASE or HTYPER, while its global accuracy was close to the best.

Comprehensive identification and characterization of diallelic insertion-deletion polymorphisms in 330 human candidate genes

Despite being the second most frequent type of polymorphism in the genome, diallelic insertion-deletion polymorphisms (indels) have received far less attention in the study of sequence variation. In this report, we describe an approach that can detect indels in the heterozygous state and can comprehensively identify indels in the target sequence. Using this approach, we identified 2393 indels in a set of 330 candidate genes, i.e. an average of seven indels per gene with about two indels per gene being common (minor allele frequency >or=0.1). We compared the population genetic characteristics of indels with substitutions in this data. Our data supported the findings that deletions occur more frequently in the human genome. 5'-UTR and coding regions of the genes showed a significantly lower diversity for indels compared with other regions, suggesting differences in effects of selection on indels and substitutions. Sequence diversity and pairwise linkage disequilibrium (LD) findings of the different populations were similar to earlier results and included a greater skew towards low-frequency variants and a faster rate of LD decay in the African-descent population compared with the non-African populations. Within populations, the allele frequency spectra and LD-decay profiles for indels were similar to substitutions. Overall, the findings suggest that, although the mechanisms giving rise to indels may be different from those causing substitutions, the evolutionary histories of indels and substitutions are similar, and that indels can play a valuable role in association studies and marker selection strategies.

Indications of linkage and association of Gilles de la Tourette syndrome in two independent family samples: 17q25 is a putative susceptibility region

Gilles de la Tourette syndrome (GTS) is characterized by multiple motor and phonic tics and high comorbidity rates with other neurobehavioral disorders. It is hypothesized that frontal-subcortical pathways and a complex genetic background are involved in the etiopathogenesis of the disorder. The genetic basis of GTS remains elusive. However, several genomic regions have been implicated. Among them, 17q25 appears to be of special interest, as suggested by various independent investigators. In the present study, we explored the possibility that 17q25 contributes to the genetic component of GTS. The initial scan of chromosome 17 performed on two large pedigrees provided a nonparametric LOD score of 2.41 near D17S928. Fine mapping with 17 additional microsatellite markers increased the peak to 2.61 (P=.002). The original families, as well as two additional pedigrees, were genotyped for 25 single-nucleotide polymorphisms (SNPs), with a focus on three genes in the indicated region that could play a role in the development of GTS, on the basis of their function and expression profile. Multiple three-marker haplotypes spanning all three genes studied provided highly significant association results (P<.001). An independent sample of 96 small families with one or two children affected with GTS was also studied. Of the 25 SNPs, 3 were associated with GTS at a statistically significant level. The transmission/disequilibrium test for a three-site haplotype moving window again provided multiple positive results. The background linkage disequilibrium (LD) of the region was studied in eight populations of European origin. A complicated pattern was revealed, with the pairwise tests producing unexpectedly high LD values at the telomeric TBCD gene. In conclusion, our findings warrant the further investigation of 17q25 as a candidate susceptibility region for GTS.

Identification and characterization of ANKK1: a novel kinase gene closely linked to DRD2 on chromosome band 11q23.1

The dopamine D2 receptor has been extensively studied in relation to alcoholism, substance abuse, and nicotine dependence. The most frequently examined polymorphism linked to this gene is the Taq1A restriction fragment length polymorphism (RFLP) (dbSNP rs1800497; g.32806C>T in GenBank AF050737.1), which has been associated with a reduction in D2 receptor density, although this is not universally accepted. The Taq1A RFLP lies 10 kB downstream of DRD2 and may therefore fall within a different coding region than the DRD2 gene or within a regulatory region. Within this downstream region, we have identified a novel kinase gene, named ankyrin repeat and kinase domain containing 1 (ANKK1), which contains a single serine/threonine kinase domain and is expressed at low levels in placenta and whole spinal cord RNA. This gene is a member of an extensive family of proteins involved in signal transduction pathways. The DRD2 Taq1A RFLP is a single nucleotide polymorphism (SNP) that causes an amino acid substitution within the 11th ankyrin repeat of ANKK1 (p.Glu713Lys), which, while unlikely to affect structural integrity, may affect substrate-binding specificity. If this is the case, then changes in ANKK1 activity may provide an alternative explanation for previously described associations between the DRD2 Taq1A RFLP and neuropsychiatric disorders such as addiction.

Understanding Human DNA Sequence Variation

Over the past century researchers have identified normal genetic variation and studied that variation in diverse human populations to determine the amounts and distributions of that variation. That information is being used to develop an understanding of the demographic histories of the different populations and the species as a whole, among other studies. With the advent of DNA-based markers in the last quarter century, these studies have accelerated. One of the challenges for the next century is to understand that variation. One component of that understanding will be population genetics. We present here examples of many of the ways these new data can be analyzed from a population perspective using results from our laboratory on multiple individual DNA-based polymorphisms, many clustered in haplotypes, studied in multiple populations representing all major geographic regions of the world. These data support an "out of Africa" hypothesis for human dispersal around the world and begin to refine the understanding of population structures and genetic relationships. We are also developing baseline information against which we can compare findings at different loci to aid in the identification of loci subject, now and in the past, to selection (directional or balancing). We do not yet have a comprehensive understanding of the extensive variation in the human genome, but some of that understanding is coming from population genetics.

The evolution and population genetics of the ALDH2 locus: random genetic drift, selection, and low levels of recombination

The catalytic deficiency of human aldehyde dehydrogenase 2 (ALDH2) is caused by a nucleotide substitution (G1510A; Glu487Lys) in exon 12 of the ALDH2 locus. This SNP, and four non-coding SNPs, including one in the promoter, span 40 kb of ALDH2; these and one downstream STRP have been tested in 37 worldwide populations. Only four major SNP-defined haplotypes account for almost all chromosomes in all populations. A fifth haplotype harbours the functional variant and is only found in East Asians. Though the SNPs showed virtually no historic recombination, LD values are quite variable because of varying haplotype frequencies, demonstrating that LD is a statistical abstraction and not a fundamental aspect of the genome, and is not a function solely of recombination. Among populations, different sets of tagging SNPs, sometimes not overlapping, can be required to identify the common haplotypes. Thus, solely because haplotype frequencies vary, there is no common minimum set of tagging SNPs globally applicable. The Fst values of the promoter region SNP and the functional SNP were about two S.D. above the mean for a reference distribution of 117 autosomal biallelic markers. These high Fst values may indicate selection has operated at these or very tightly linked sites.

Alcohol-related expectancies are associated with the D2 dopamine receptor and GABAA receptor beta3 subunit genes

Molecular genetic research has identified promising markers of alcohol dependence, including alleles of the D2 dopamine receptor (DRD2) and the GABAA receptor beta3 subunit (GABRB3) genes. Whether such genetic risk manifests itself in stronger alcohol-related outcome expectancies, or in difficulty resisting alcohol, is unknown. In the present study, A1+ (A1A1 and A1A2 genotypes) and A1- (A2A2 genotype) alleles of the DRD2 and G1+ (G1G1 and G1 non-G1 genotypes) and G1- (non-G1 non-G1 genotype) alleles of the GABRB3 gene were determined in a group of 56 medically ill patients diagnosed with alcohol dependence. Mood-related alcohol expectancy (AE) and drinking refusal self-efficacy (DRSE) were assessed using the Drinking Expectancy Profile (Manual for the Drinking Expectancy Profile, Behaviour Research and Therapy Centre, Brisbane, 1996). Patients with the DRD2 A1+ allele, compared with those with the DRD2 A1- allele, reported significantly lower DRSE in situations of social pressure. Similarly, lower DRSE was reported under social pressure by patients with the GABRB3 G1+ allele when compared to those with the GABRB3 G1- alleles. Patients with the GABRB3 G1+ allele also revealed reduced DRSE in situations characterized by negative affect than those with the GABRB3 G1- alleles. Patients carrying the GABRB3 G1+ allele showed stronger AE relating to negative affective change (for example, increased depression) than their GABRB3 G1- counterparts. Biological influence in the development of some classes of cognitions is hypothesized. The clinical implications, particularly with regard to patient-treatment matching and the development of an integrated psychological and pharmacogenetic approach, are discussed.

Haplotypic relationship between SNP and microsatellite markers at the NOS2A locus in two populations

The density of genetic markers required for successful association mapping of complex diseases depends on linkage disequilibrium (LD) between non-functional markers and functional variants. The haplotypic relationship between stable markers and potentially unstable but highly informative markers (e.g. microsatellites) indicates that LD might be maintained over considerable genetic distance in non-African populations, supporting the use of such 'mixed marker haplotypes' in LD-based mapping, and allowing inferences to be drawn about human origins. We investigated sequence variation in the proximal 2.6 kb of the inducible nitric oxide synthase (NOS2A) promoter and the relationship between SNP haplotypes and a pentanucleotide microsatellite (the 'NOS2A(-2.6) microsatellite') in Gambians and UK Caucasians. UK Caucasians exhibited a subset of sequence diversity observed in Gambians, sharing four of 11 SNPs and a similar haplotypic structure. Five SNPs were found in the sequence of interspersed repetitive DNA elements. In both populations, there was dramatic loss of LD between SNP haplotypes and microsatellite alleles across a very short physical distance, suggesting a high intrinsic mutation rate of the NOS2A(-2.6) microsatellite, the SNP haplotypes are relatively ancient, or that this was a region of frequent recombination. Understanding locus- and population-specific LD is essential when designing and interpreting genetic association studies.

Family-based association studies between 5-HT5A receptor gene and schizophrenia

BACKGROUND

Pharmacological and neurodevelopmental data support the idea that the gene, which codes for the 5-HT(5A) receptor is an important candidate gene for schizophrenia susceptibility. However, previous genetic studies focusing on this gene yielded conflicting results, potentially because of: (i) stratification biases of case-control association studies, (ii) genetic and phenotypic heterogeneity of schizophrenia, and (iii) variability in the loci analyzed (the 5-HT(5A) gene having many polymorphic sites).

METHODS

A transmission disequilibrium test was used in the present study aimed at investigating two polymorphisms in exon 1 of the 5-HT(5A) gene, the A12T silent substitution and the C43T transversion leading to a 15Pro --> Ser substitution, in 103 patients with DSM-IV diagnosis of schizophrenia, and their 206 parents.

RESULTS

We found an excess of transmission of the 12T allele from the parents to their affected children (P = 0.02), with evidence for linkage disequilibrium between the 12T-43C haplotype and schizophrenia (P = 0.002). Furthermore, patients with the 12T allele had a significantly later age at onset (P = 0.003), and the Q-TDT approach confirmed that this allele was transmitted with an older age at onset (P = 0.01).

CONCLUSIONS

These data provided convergent evidence for a significant role of the 5-HT(5A) gene in schizophrenia and more specifically in patients with later age at onset.

Possible interaction of alcohol dehydrogenase and aldehyde dehydrogenase genes with the dopamine D2 receptor gene in anxiety-depressive alcohol dependence

BACKGROUND

The role of the dopamine D2 receptor (DRD2) gene in the development of alcohol abuse or dependence is controversial. The controversy is due in part to the disparate definitions pertaining to the control groups used and to the definitions of subtypes in alcohol dependence. In the Han Chinese population, the alcohol dehydrogenase 1B*2/*2 (ADH1B*2/*2) genotype and the aldehyde dehydrogenase 2*2 (ALDH2*2) allele have been considered as protective factors against alcohol abuse or dependence. Moreover, the ADH1B and ALDH2 genes might be involved in dopamine metabolism. We hypothesized that the ADH1B and ALDH2 genes might interact with the DRD2 gene and that the association between the DRD2 gene and alcohol dependence might be affected by different ADH1B and ALDH2 genotypes. This study examined whether the DRD2 gene is associated with specific subtypes of alcohol dependence and evaluated the relationship between the DRD2 gene and alcohol-metabolizing genes in a specific subtype of alcohol dependence.

METHODS

Of the 465 Han Chinese subjects who were recruited for the study, 71 were classified with pure alcohol dependence, 113 with both alcohol dependence and anxiety-depression (ANX/DEP ALC), and 129 with anxiety-depression but without alcohol dependence (ANX/DEP). The remaining 152 subjects were supernormal controls. All subjects were interviewed with the Chinese version of the modified Schedule of Affective Disorders and Schizophrenia-Lifetime; all alcohol dependence, anxiety, and major depressive diagnoses were made according to DSM-IV criteria.

RESULTS

The DRD2 gene was not found to be associated with pure alcohol dependence or ANX/DEP, but was found to be associated with ANX/DEP ALC. Furthermore, the association between the DRD2 gene and ANX/DEP ALC was shown to be under the control of the ALDH2*1/*1 and ADH1B*1/*2 genotypes.

CONCLUSIONS

ANX/DEP ALC is a specific subtype of alcohol dependence. Because ANX/DEP ALC was associated with the DRD2 gene only under the stratification of ADH1B*1/*2 or ALDH2*1/*1, the DRD2 gene might interact with the ADH1B gene and the ALDH2 gene, respectively, in the development of ANX/DEP ALC in the Taiwan Han Chinese population.

Population stratification in the candidate gene study: fatal threat or red herring?

Advances in molecular genetics have provided behavioral scientists with a means of investigating the influence of genetic factors on human behavior. Unfortunately, recent candidate gene studies have produced inconsistent results, and a frequent scapegoat for the lack of replication across studies is the threat of population stratification. This review of the literature on population stratification suggests that the threat may be a red herring. Reliable findings will require improved specification and measurement of the behavioral phenotypes in question, a renewed focus on internal validity, and the specification and testing of genetic factors in the context of longitudinal multivariate models. In this respect, behavioral scientists are well suited to investigating genetic factors that influence psychological mechanisms.

Genetic structure and affinities among tribal populations of southern India: a study of 24 autosomal DNA markers

We describe the genetic structure and affinities of five Dravidian-speaking tribal populations inhabiting the Nilgiri hills of Tamil Nadu, in south India, using 24 autosomal DNA markers. Our goals were: (i). to examine what evolutionary forces have most significantly impacted south Indian tribal genetic variation, and (ii). to test whether the phenotypic similarities of some south Indian tribal groups to Africans represent a signature of close relationship to Africans or are due to convergence. All loci were polymorphic and average heterozygosities were substantial (range: 0.347-0.423). Genetic differentiation was high (Gst= 6.7%) and genetic distances were not significantly correlated with geographic distances. Genetic drift therefore probably played a significant role in shaping the patterns of genetic variation observed in southern Indian tribal populations. Otherwise, analyses of population relationships showed that Indian populations are closely related to one another, regardless of phenotypic characteristics, and do not show particular affinities to Africans. We conclude that the phenotypic similarities of some Indian groups to Africans do not reflect a close relationship between these groups, but are better explained by convergence.

Patterns of human genetic diversity: implications for human evolutionary history and disease

Since the completion of the human genome sequencing project, the discovery and characterization of human genetic variation is a principal focus for future research. Comparative studies across ethnically diverse human populations and across human and nonhuman primate species is important for reconstructing human evolutionary history and for understanding the genetic basis of human disease. In this review, we summarize data on patterns of human genetic diversity and the evolutionary forces (mutation, genetic drift, migration, and selection) that have shaped these patterns of variation across both human populations and the genome. African population samples typically have higher levels of genetic diversity, a complex population substructure, and low levels of linkage disequilibrium (LD) relative to non-African populations. We discuss these differences and their implications for mapping disease genes and for understanding how population and genomic diversity have been important in the evolution, differentiation, and adaptation of humans.

Dopamine receptor D2 -141C Insertion/Deletion polymorphism in a Finnish population with schizophrenia

We examined the occurrence of the -141C Ins/Del polymorphism in 93 Finnish patients with schizophrenia. In comparison with previous studies with Japanese and Caucasian populations, the incidence of this polymorphism was unexpectedly low. The findings suggest that the frequency of the -141C Ins/Del polymorphism is lower in Northern Europe compared to other Caucasian and Japanese populations.

Global survey of haplotype frequencies and linkage disequilibrium at the RET locus

We have constructed haplotypes based on normal variation at six polymorphic sites-five single nucleotide polymorphisms (SNPs) and one short tandem repeat polymorphism (STRP)-at the RET locus for samples of normal individuals from 32 populations distributed across the major continental regions of the world. The haplotyped system spans 41.6 kilobases and encompasses most of the coding region of the gene. All of the markers are polymorphic in all regions of the world and in most individual populations. Expected heterozygosities for the six-site haplotypes range from 82 to 94% for all populations studied except for two Amerindian groups from the Amazon basin at 61 and 76%. Individual populations had from four to eight haplotypes with frequencies exceeding 5%. In general, African, southwest Asian and European groups have the highest numbers of total and of commonly occurring haplotypes; the lowest numbers are observed in Amerindian populations. Overall linkage disequilibrium (LD) for the five SNP sites was very significant (P</=0.001) for all the non-African populations, but significant at that level for only one of the seven African populations. In general, the permutation-based xi coefficient that quantifies overall LD tends to increase the farther the population is from Africa, but variability of this measure of LD is often large within geographic regions. Pairwise LD measures among the SNPs also show considerable variation among populations. Association of STRP alleles with the SNP-defined background haplotypes is generally higher outside of Africa than in Africa, but is highly variable.

The loss of statistical power to distinguish populations when certain samples are ambiguous

Case-control studies are used to map loci associated with a genetic disease. The usual case-control study tests for significant differences in frequencies of alleles at marker loci. In this paper, we consider the problem of comparing two or more marker loci simultaneously and testing for significant differences in haplotype rather than allele frequencies. We consider two situations. In the first, genotypes at marker loci are resolved into haplotypes by making use of biochemical methods or by genotyping family members. In the second, genotypes at marker loci are not resolved into haplotypes, but, by assuming random mating, haplotypes can be inferred using a likelihood method such as the expectation-maximization (EM) algorithm. We assume that a causative locus has two alleles with a multiplicative effect on the penetrance of a disease, with one allele increasing the penetrance by a factor pi. We find, for small values of pi-1 and large sample sizes, asymptotic results that predict the statistical power of a test for significant differences in haplotype frequencies between cases and a random sample of the population, both when haplotypes can be resolved and when haplotypes have to be inferred. The increase in power when haplotypes can be resolved can be expressed as a ratio R, which is the increase in sample size needed to achieve the same power when haplotypes are resolved over when they are not resolved. In general, R depends on the pattern of linkage disequilibrium between the causative allele and the marker haplotypes but is independent of the frequency of the causative allele and, to a first approximation, is independent of pi. For the special situation of two di-allelic marker loci, we obtain a simple expression for R and its upper bound.

Nucleotide variation, haplotype structure, and association with end-stage renal disease of the human interleukin-1 gene cluster

A dense gene-based SNP map was constructed across a 360-kb region containing the interleukin-1 gene cluster (IL1A, IL1B, and IL1RN), focusing on IL1RN. In total, 95 polymorphisms were confirmed or identified primarily by direct sequencing. Polymorphisms were precisely mapped to completed BAC and genomic sequences spanning this region. The polymorphisms were typed in 443 case-control subjects from Caucasian and African American groups. Consecutive pair-wise marker linkage disequilibrium was not strictly correlated with distance and ranged from D'=0.0079 to 1.000 and D'=0.0521 to 1.0000 in Caucasians and African Americans, respectively. Single markers and haplotypes in IL1 cluster genes were evaluated for association with end-stage renal disease (ESRD). Eleven SNPs show some evidence of association with ESRD, with the strongest associations in two IL1A variants, one SNP, rs1516792-3, in intron 5 (p=0.0015) and a 4-bp insertion/deletion within the 3'UTR, rs16347-2 (p=0.0024), among African Americans with non-T2DM-associated ESRD.

Dopamine D2 receptor TaqIA and TaqIB polymorphisms in Parkinson's disease

In a case control study, we examined the association of DRD2 Taq1A and Taq1B polymorphisms and risk of PD, and evaluated the strength of linkage disequilibrium of the polymorphisms. The Taq1A and Taq1B polymorphisms were in strong linkage disequilibrium. There was, however, no significant association of the two polymorphisms with PD.

The application of molecular genetic approaches to the study of human evolution

The past decade of advances in molecular genetic technology has heralded a new era for all evolutionary studies, but especially the science of human evolution. Data on various kinds of DNA variation in human populations have rapidly accumulated. There is increasing recognition of the importance of this variation for medicine and developmental biology and for understanding the history of our species. Haploid markers from mitochondrial DNA and the Y chromosome have proven invaluable for generating a standard model for evolution of modern humans. Conclusions from earlier research on protein polymorphisms have been generally supported by more sophisticated DNA analysis. Co-evolution of genes with language and some slowly evolving cultural traits, together with the genetic evolution of commensals and parasites that have accompanied modern humans in their expansion from Africa to the other continents, supports and supplements the standard model of genetic evolution. The advances in our understanding of the evolutionary history of humans attests to the advantages of multidisciplinary research.

Extensive linkage disequilibrium, a common 16.7-kilobase deletion, and evidence of balancing selection in the human protocadherin alpha cluster

Regions of extensive linkage disequilibrium (LD) appear to be a common feature of the human genome. However, the mechanisms that maintain these regions are unknown. In an effort to understand whether gene density contributes to LD, we determined the degree of promoter sequence variation in a large tandem-arrayed gene family, the human protocadherin alpha cluster, on chromosome 5. These genes are expressed at synaptic junctions in the developing brain and the adult brain and may be involved in the determination of synaptic complexity. We sequenced the promoters of all 13 alpha protocadherin genes in 96 European Americans and identified polymorphisms in the promoters alpha 1, alpha 3, alpha 4, alpha 5, alpha 7, alpha 9, alpha 11, and alpha 13. In these promoters, 11 common SNPs are in extensive LD, forming two 48-kb haplotypes of equal frequency, in this population, that extend from the alpha1 through alpha 7 genes. We sequenced these promoters in East Asians and African Americans, and we estimated haplotype frequencies and calculated LD statistics for all three populations. Our results indicate that, although extensive LD is an ancient feature of the alpha cluster, it has eroded over time. SNPs 3' of alpha 7 are involved in ancestral recombination events in all populations, and overall alpha-cluster LD is reduced in African Americans. We obtained significant positive values for Tajima's D test for all alpha promoter SNPs in Europeans (D=3.03) and East Asians (D=2.64), indicating an excess of intermediate-frequency variants, which is a signature of balancing selection. We also discovered a 16.7-kb deletion that truncates the alpha 8 gene and completely removes the alpha 9 and alpha 10 genes. This deletion appears in unaffected individuals from multiple populations, suggesting that a reduction in protocadherin gene number is not obviously deleterious.

D2 dopamine receptor gene in psychiatric and neurologic disorders and its phenotypes

The D2 dopamine receptor (DRD2) has been one of the most extensively investigated gene in neuropsychiatric disorders. After the first association of the TaqI A DRD2 minor (A1) allele with severe alcoholism in 1990, a large number of international studies have followed. A meta-analysis of these studies of Caucasians showed a significantly higher DRD2 A1 allelic frequency and prevalence in alcoholics when compared to controls. Variants of the DRD2 gene have also been associated with other addictive disorders including cocaine, nicotine and opioid dependence and obesity. It is hypothesized that the DRD2 is a reinforcement or reward gene. The DRD2 gene has also been implicated in schizophrenia, posttraumatic stress disorder, movement disorders and migraine. Phenotypic differences have been associated with DRD2 variants. These include reduced D2 dopamine receptor numbers and diminished glucose metabolism in brains of subjects who carry the DRD2 A1 allele. In addition, pleiotropic effects of DRD2 variants have been observed in neurophysiologic, neuropsychologic, stress response, personality and treatment outcome characteristics. The involvement of the DRD2 gene in certain neuropsychiatric disorders opens up the potential of a targeted pharmacogenomic approach to the treatment of these disorders.

Association of seven polymorphisms of the D2 dopamine receptor gene with brain receptor-binding characteristics

Association of alleles at the Taq1 A, Taq1 B, intron 6, Taq1 D, exon 7, exon 8, and promoter-141C sites of the D2 dopamine receptor gene with D2 dopamine receptor binding characteristics in the caudate nucleus of Caucasian alcoholic and nonalcoholic subjects was determined. For the Taq1 D, exon 7, exon 8, and promoter- 141C sites there were no significant allelic differences in Bmax (number of binding sites) or Kd (binding affinity) of the D2 dopamine receptors. However, subjects having the minor alleles at the Taq1 A, Taq1 B, and intron 6 sites had significantly lower Bmax than subjects not having them. None of these three polymorphisms had any significant effect on Kd. Highly significant linkage disequilibria were observed among the Taq1 A, Taq1 B, and intron 6 polymorphic sites, but linkage disequilibria between these three sites and each of the Taq1 D, exon 7, exon 8, and promoter-141C sites were of lesser or of no significance. Taken together, these findings suggest that the Taq1 A, Taq1 B, and intron 6 polymorphisms, but not the Taq1 D, exon 7, exon 8, and promoter-141C polymorphisms, are in linkage disequilibrium with a functional allelic variant that affects D2 dopamine receptor expression.

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.

ALFRED: An allele frequency database for anthropology

The deluge of data from the human genome project (HGP) presents new opportunities for molecular anthropologists to study human variation through the promise of vast numbers of new polymorphisms (e.g., single nucleotide polymorphisms or SNPs). Collecting the resulting data into a single, easily accessible resource will be important to facilitate this research. We created a prototype Web-accessible database named ALFRED (ALelle FREquency Database, http://alfred.med.yale.edu/alfred/) to store and make publicly available allele frequency data on diverse polymorphic sites for many populations. In constructing this database, we considered many different concerns relating to the types of information needed for anthropology, population genetics, molecular genetics, and statistics, as well as issues of data integrity and ease of access to data. We also developed links to other Web-based databases as well as procedures for others to make links to the data in ALFRED. Here we present an overview of the issues considered and provisional solutions, as well as an example of data already available. It is our hope that this database will be useful for research and teaching in a wide range of fields, and that colleagues from various fields will contribute to making ALFRED an important resource for many studies as yet unforeseen.

Sequence variation and haplotype structure at the human HFE locus

The HFE locus encodes an HLA class-I-type protein important in iron regulation and segregates replacement mutations that give rise to the most common form of genetic hemochromatosis. The high frequency of one disease-associated mutation, C282Y, and the nature of this disease have led some to suggest a selective advantage for this mutation. To investigate the context in which this mutation arose and gain a better understanding of HFE genetic variation, we surveyed nucleotide variability in 11.2 kb encompassing the HFE locus and experimentally determined haplotypes. We fully resequenced 60 chromosomes of African, Asian, or European ancestry as well as one chimpanzee, revealing 41 variable sites and a nucleotide diversity of 0.08%. This indicates that linkage to the HLA region has not substantially increased the level of HFE variation. Although several haplotypes are shared between populations, one haplotype predominates in Asia but is nearly absent elsewhere, causing higher than average genetic differentiation among the three major populations. Our samples show evidence of intragenic recombination, so the scarcity of recombination events within the C282Y allele class is consistent with selection increasing the frequency of a young allele. Otherwise, the pattern of variability in this region does not clearly indicate the action of positive selection at this or linked loci.

Lipoprotein lipase and APOE/APOC-I/APOC-II gene cluster diversity in native Brazilian populations

Allele and haplotype frequencies for the T-93G, Hind III, and Pvu II variants of the lipoprotein lipase gene (LPL), and Hpa I and Ava II restriction site polymorphisms (RSP) of the APOE/C-I/C-II gene cluster were determined in 143 individuals from five Brazilian Indian tribes. These results were integrated with those previously reported for APOE. Marked interethnic variability occurs in these sites. A strong linkage disequilibrium was observed between the APOE and APOC-I loci (D' = 0.81; P < 0.00001). Linkage disequilibrium between the Hind III and Pvu II RSPs of the LPL gene was also observed (D' = 1; P < 0.001), but none of these RSPs were in linkage disequilibrium with the T-93G mutation. Considering both loci, heterozygosity was estimated as 0.45, but it was lower in the Xavante and Surui populations, in accordance with the historical and biodemographical data of these Amerindians. The results reported here may have implications for understanding interpopulation differences in lipid levels and coronary heart disease prevalences.

A global perspective on genetic variation at the ADH genes reveals unusual patterns of linkage disequilibrium and diversity

Variants of different Class I alcohol dehydrogenase (ADH) genes have been shown to be associated with an effect that is protective against alcoholism. Previous work from our laboratory has shown that the two sites showing the association are in linkage disequilibrium and has identified the ADH1B Arg47His site as causative, with the ADH1C Ile349Val site showing association only because of the disequilibrium. Here, we describe an initial study of the nature of linkage disequilibrium and genetic variation, in population samples from different regions of the world, in a larger segment of the ADH cluster (including the three Class I ADH genes and ADH7). Linkage disequilibrium across approximately 40 kb of the Class I ADH cluster is moderate to strong in all population samples that we studied. We observed nominally significant pairwise linkage disequilibrium, in some populations, between the ADH7 site and some Class I ADH sites, at moderate values and at a molecular distance as great as 100 kb. Our data indicate (1) that most ADH-alcoholism association studies have failed to consider many sites in the ADH cluster that may harbor etiologically significant alleles and (2) that the relevance of the various ADH sites will be population dependent. Some individual sites in the Class I ADH cluster show Fst values that are among the highest seen among several dozen unlinked sites that were studied in the same subset of populations. The high Fst values can be attributed to the discrepant frequencies of specific alleles in eastern Asia relative to those in other regions of the world. These alleles are part of a single haplotype that exists at high (>65%) frequency only in the eastern-Asian samples. It seems unlikely that this haplotype, which is rare or unobserved in other populations, reached such high frequency because of random genetic drift alone.

Molecular variability in Amerindians: widespread but uneven information

A review was made in relation to the molecular variability present in North, Central, and South American Indian populations. It involved results from ancient DNA, mitochondrial DNA in extant populations, HLA and other autosomal markers, X and Y chromosome variation, as well as data from parasitic viruses which could show coevolutionary changes. The questions considered were their origin, ways in which the early colonization of the continent took place, types and levels of the variability which developed, peculiarities of the Amerindian evolutionary processes, and eventual genetic heterogeneity which evolved in different geographical areas. Although much information is already available, it is highly heterogeneous in relation to populations and types of genetic systems investigated. Unfortunately, the present trend of favoring essentially applied research suggest that the situation will not basically improve in the future.