The complex global pattern of genetic variation and linkage disequilibrium at catechol-O-methyltransferase.

Genetic variation at the catechol-O-methyltransferase (COMT) gene has been significantly associated with risk for various neuropsychiatric conditions such as schizophrenia, panic disorder, bipolar disorders, anorexia nervosa and others. It has also been associated with nicotine dependence, sensitivity to pain and cognitive dysfunctions especially in schizophrenia. The non-synonymous single nucleotide polymorphism (SNP) in exon 4--Val108/158Met--is the most studied SNP at COMT and is the basis for most associations. It is not, however, the only variation in the gene; several haplotypes exist across the gene. Some studies indicate that the haplotypic combinations of alleles at the Val108/158Met SNP with those in the promoter region and in the 3'-untranslated region are responsible for the associations with disorders and not the non-synonymous SNP by itself. We have now studied DNA samples from 45 populations for 63 SNPs in a region of 172 kb across the region of 22q11.2 encompassing the COMT gene. We focused on 28 SNPs spanning the COMT-coding region and immediately flanking DNA, and found that the haplotypes are from diverse evolutionary lineages that could harbor as yet undetected variants with functional consequences. Future association studies should be based on SNPs that define the common haplotypes in the population(s) being studied.

Global patterns of variation in allele and haplotype frequencies and linkage disequilibrium across the CYP2E1 gene.

Cytochrome P450 2E1, gene symbol CYP2E1, is one of a family of enzymes with a central role in activating and detoxifying xenobiotics and endogenous compounds. Genetic variation at this gene has been reported in different human populations, and some association studies have reported increased risk for cancers and other diseases. To the best of our knowledge, multi-single-nucleotide polymorphism haplotypes and linkage disequilibrium (LD) have not been systematically studied for CYP2E1 in multiple populations. Haplotypes can greatly increase the power both to identify patterns of genetic variation relevant for gene expression as well as to detect disease-related susceptibility mutations. We present frequency and LD data and analyses for 11 polymorphisms and their haplotypes that we have studied on over 2600 individuals from 50 human population samples representing the major geographical regions of the world. The diverse patterns of haplotype variation found in the different populations we have studied show that ethnicity may be an important variable helping to explain inconsistencies that have been reported by association studies. More studies clearly are needed of the variants we have studied, especially those in the 5' region, such as the variable number of tandem repeats, as well as studies of additional polymorphisms known for this gene to establish evidence relating any systematic differences in gene expression that exist to the haplotypes at this gene.

Motoneuron-specific NR3B gene: no association with ALS and evidence for a common null allele.

OBJECTIVE: The GRIN3B gene encodes NR3B, a motoneuron-specific member of the NMDA type of ionotropic glutamate receptors. NR3B reduces the Ca(2+)-permeability as well as the overall current of the receptor response and may thereby protect motoneurons against glutamate-mediated excitotoxicity. We tested whether genetic dysfunction of GRIN3B is implicated in the pathogenesis of amyotrophic lateral sclerosis (ALS). METHODS: We searched for mutations in the GRIN3B coding region (3.1 kb) in 117 individuals with familial ALS and in 46 individuals with sporadic ALS. We genotyped the newly identified GRIN3B null allele and four "tag single nucleotide polymorphisms (SNPs)" at the GRIN3B locus in 342 individuals with sporadic ALS and in 374 matched controls. The GRIN3B null allele frequency was determined in 2,128 individuals from a worldwide panel of 42 populations. We furthermore compared the GRIN3B coding sequence in primates (human-macaque) and rodents (rat-mouse) to evaluate the molecular evolution of GRIN3B. RESULTS: Thirty-two SNPs, including 16 previously unreported SNPs, one 27-bp deletion, a polymorphic CAG repeat, and a 4-bp insertion (insCGTT), were identified. Mutational and case-control studies did not reveal variants that cause or modify disease in ALS. Intriguing is an insCGTT variant that truncates the protein at its amino terminus and results in a GRIN3B null allele. We demonstrated a global distribution of the null allele with allele frequencies ranging between 0 and 0.38, and we delineated a null allele specific haplotype of 9.89 kb. Comparative genomic analysis across four taxa demonstrated accelerated evolution of NR3B in primates. CONCLUSIONS: Our study supports the conclusions that 1) GRIN3B does not seem to be associated with familial or sporadic ALS, 2) the GRIN3B null allele is a common polymorphism, 3) the GRIN3B null allele has arisen once and early in human evolution, and 4) the GRIN3B gene belongs to a group of nervous system-related genes that have been subjected to faster evolution during evolution.

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.

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.

Genealogy reconstruction from short tandem repeat genotypes in an Amazonian population.

We have reconstructed partial genealogies in a sample of 44 SW Amazonian Rondonian Surui, in which 45 dinucleotide short tandem repeat polymorphisms had previously been typed. The genotypes of 488 pairs of individuals having an age difference of 13 or greater were compared, and parentage was excluded if a pair failed to share an allele at more than one locus. In order to test the power of this method, we computed the expected distribution of the number of exclusionary loci for such pairs of unrelated individuals, as well as that for individuals with different degrees of relatedness, and compared it to the observed distribution. We estimated that the pairs compared contained approximately 20% of individual pairs with a first-cousin relation or closer. A total of 25 pairs were identified as possible parent-child. In three instances, we could identify two or more children having a common parent; we computed a relatedness coefficient in order to establish whether the children were full or half sibs. The genealogies inferred show that instances of polygyny and polyandry (or, alternatively, serial mating), in addition to apparent monogamy, can be found among the Surui. The Surui sample can be used as a model for paleoanthropological populations, in which the determination of relatedness can provide further insights into the social structure of past populations. We estimate that, depending on the history of the populations and the degree of inbreeding, 10-20 highly informative nuclear loci should be typed in order to infer genealogies with acceptable confidence.

Short tandem repeat polymorphism evolution in humans.

Forty-five dinucleotide short tandem repeat polymorphisms were typed in ten large samples of a globally distributed set of populations. Although these markers had been selected for high heterozygosity in European populations, we found them to be sufficiently informative for linkage analysis in non-Europeans. Heterozygosity, mean number of alleles, and mean number of private alleles followed a common trend: they were highest in the African samples, were somewhat lower in Europeans and East Asians, and were lowest in Amerindians. Genetic distances also reflected this pattern, and distances modelled after the stepwise mutation model yielded trees that were less in agreement with other genetic and archaeological evidence than distances based on differentiation by drift (FST). Genetic variation in non-Africans seems to be a subset of that in Africans, supporting the replacement hypothesis for the origin of modern humans.

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

A four-site haplotype system at the dopamine D2 receptor locus (DRD2) has been studied in a global sample of 28 distinct populations. The haplotype system spans about 25 kb, encompassing the coding region of the gene. The four individual markers include three TaqI restriction site polymorphisms (RSPs) -- TaqI "A", "B", and "D" sites -- and one dinucleotide short tandem repeat polymorphism (STRP). All four of the marker systems are polymorphic in all regions of the world and in most individual populations. The haplotype system shows the highest average heterozygosity in Africa, a slightly lower average heterozygosity in Europe, and the lowest average heterozygosities in East Asia and the Americas. Across all populations, 20 of the 48 possible haplotypes reached a frequency of at least 5% in at least one population sample. However, no single population had more than six haplotypes reaching that frequency. In general, African populations had more haplotypes present in each population and more haplotypes occurring at a frequency of at least 5% in that population. Permutation tests for significance of overall disequilibrium (all sites considered simultaneously) were highly significant (P<0.001) in all 28 populations. Except for three African samples, the pairwise disequilibrium between the outermost RSP markers, TaqI "B" and "A", was highly significant with D' values greater than 0.8; in two of those exceptions the RSP marker was not polymorphic. Except for those same two African populations, the 16-repeat allele at the STRP also showed highly significant disequilibrium with the TaqI "B" site in all populations, with D' values usually greater than 0.7. Only four haplotypes account for more than 70% of all chromosomes in virtually all non-African populations, and two of those haplotypes account for more than 70% of all chromosomes in most East Asian and Amerindian populations. A new measure of the amount of overall disequilibrium shows least disequilibrium in African populations, somewhat more in European populations, and the greatest amount in East Asian and Amerindian populations. This pattern seems best explained by random genetic drift with low levels of recombination, a low mutation rate at the STRP, and essentially no recurrent mutation at the RSP sites, all in conjunction with an "Out of Africa" model for recent human evolution.

A global haplotype analysis of the myotonic dystrophy locus: implications for the evolution of modern humans and for the origin of myotonic dystrophy mutations.

Haplotypes consisting of the (CTG)n repeat, as well as several flanking markers at the myotonic dystrophy (DM) locus, were analyzed in normal individuals from 25 human populations (5 African, 2 Middle Eastern, 3 European, 6 East Asian, 3 Pacific/Australo-Melanesian, and 6 Amerindian) and in five nonhuman primate species. Non-African populations have a subset of the haplotype diversity present in Africa, as well as a shared pattern of allelic association. (CTG)18-35 alleles (large normal) were observed only in northeastern African and non-African populations and exhibit strong linkage disequilibrium with three markers flanking the (CTG)n repeat. The pattern of haplotype diversity and linkage disequilibrium observed supports a recent African-origin model of modern human evolution and suggests that the original mutation event that gave rise to DM-causing alleles arose in a population ancestral to non-Africans prior to migration of modern humans out of Africa.

Susceptibility loci for distinct components of developmental dyslexia on chromosomes 6 and 15.

Six extended dyslexic families with at least four affected individuals were genotyped with markers in three chromosomal regions: 6p23-p21.3, 15pter-qter, and 16pter-qter. Five theoretically derived phenotypes were used in the linkage analyses: (1) phonological awareness; (2) phonological decoding; (3) rapid automatized naming; (4) single-word reading; and (5) discrepancy between intelligence and reading performance, an empirically derived, commonly used phenotype. Two-point and multipoint allele-sharing analyses of chromosome 6 markers revealed significant evidence (P < 10(-6)) for linkage of the phonological awareness phenotype to five adjacent markers (D6S109, D6S461, D6S299, D6S464, and D6S306). The least compelling results were obtained with single-word reading. In contrast, with chromosome 15 markers, a LOD score of 3.15 was obtained for marker D15S143 at theta = 0.0 with single-word reading. Multipoint analyses with markers adjacent to D15S143 (D15S126, D15S132, D15S214, and D15S128) were positive, but none reached acceptable significance levels. Chromosome 15 analyses with the phonological awareness phenotype were negative. Parametric and nonparametric linkage analyses with chromosome 16 markers were negative. The most intriguing aspect of the current findings is that two very distinct reading-related phenotypes, reflecting different levels in the hierarchy of reading-related skills, each contributing to different processes, appear to be linked to two different chromosomal regions.

DRD2 haplotypes containing the TaqI A1 allele: implications for alcoholism research.

In recent years, a possible role of the dopamine D2 receptor (DRD2) locus in the etiology of alcoholism has been the focus of considerable attention. The literature now contains a mix of association studies with positive and negative conclusions. Various methodological flaws undermine the claims in many of the studies that conclude a positive association exists between alcoholism and the DRD2*A1 allele at the Taql "A" site. Although the studies with negative findings have more often come from studies using better analytic methodology, satisfactory resolution of whether or not genetic variation at the DRD2 locus plays some role in the etiology of alcoholism is unlikely to come from additional studies of the kind conducted thus far; an approach enlightened by a more thorough understanding of the population genetics of DRD2 and the phylogenetic origins of the DRD2 alleles is one alternative. If genetic variation at the DRD2 locus affects susceptibility to alcoholism, then such variation has a mutational and evolutionary history that can be traced with the aid of the various genetic polymorphisms that have been identified at the DRD2 locus. In this study, a third Taql restriction fragment-length polymorphism at DRD2, the Taql "D" site, has been converted to polymerase chain reaction-based typing and its frequencies determined in 22 populations from around the world. Haplotypes defined by the polymorphisms at the Taql "B" and "A" sites, and the short tandem repeat polymorphism in intron 2 have been constructed and the diversity of haplotypes containing the DRD2*A1 allele examined for all 22 populations. The ancestral origins of the three Taql polymorphisms have also been determined by sequencing the homologous regions in other higher primates. Because A1-containing haplotypes in populations of European, Middle Eastern, and African origin show considerable diversity within and among populations, properly designed association studies in populations descended from those areas of the world need to use haplotypes, not a single allelic system, and need to use appropriate methods to compensate for the near impossibility of genetically matching unrelated control samples.

Evolution of haplotypes at the DRD2 locus.

We present here the first evolutionary perspective on haplotypes at DRD2, the locus for the dopamine D2 receptor. The dopamine D2 receptor plays a critical role in the functioning of many neural circuits in the human brain. If functionally relevant variation at the DRD2 locus exists, understanding the evolution of haplotypes on the basis of polymorphic sites encompassing the gene should provide a powerful framework for identifying that variation. Three DRD2 polymorphisms (TaqI "A" and "B" RFLPs and the (CA)n short tandem repeat polymorphic in all the populations studied, and they display strong and significant linkage disequilibria with each other. The common haplotypes for the two TaqI RFLPs are separately derived from the ancestral haplotype but predate the spread of modern humans around the world. The knowledge of how the various haplotypes have evolved, the allele frequencies of the haplotypes in human populations, and the physical relationships of the polymorphisms to each other and to the functional parts of the gene should now allow proper design and interpretation of association studies.