Genome-wide association study of Tourette's syndrome.

Tourette's syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10(-8)); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10(-6)). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10(-7) for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

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.

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.

The genetic architecture of selection at the human dopamine receptor D4 (DRD4) gene locus.

Associations of the seven-repeat (7R) allele of the human dopamine receptor D4 (DRD4) gene with both the personality trait of novelty seeking and attention deficit/hyperactivity disorder have been reported. Recently, on the basis of the unusual DNA sequence organization of the DRD4 7R 48-bp tandem repeat (VNTR), we proposed that the 7R allele originated as a rare mutational event that increased to high frequency by positive selection. We now have resequenced the entire DRD4 locus from 103 individuals homozygous for 2R, 4R, or 7R variants of the VNTR, a method developed to directly estimate haplotype diversity. DNA from individuals of African, European, Asian, North and South American, and Pacific Island ancestry were used. 4R/4R homozygotes exhibit little linkage disequilibrium (LD) over the region examined, with more polymorphisms observed in DNA samples from African individuals. In contrast, the evidence for strong LD surrounding the 7R allele is dramatic, with all 7R/7R individuals (including those from Africa) exhibiting the same alleles at most polymorphic sites. By intra-allelic comparison at 18 high-heterozygosity sites spanning the locus, we estimate that the 7R allele arose prior to the upper Paleolithic era (approximately 40000-50000 years ago). Further, the pattern of recombination at these polymorphic sites is the pattern expected for selection acting at the 7R VNTR itself, rather than at an adjacent site. We propose a model for selection at the DRD4 locus consistent with these observed LD patterns and with the known biochemical and physiological differences between receptor variants.

COMT haplotypes suggest P2 promoter region relevance for schizophrenia.

A recent study found, in a large sample of Ashkenazi Jews, a highly significant association between schizophrenia and a particular haplotype of three polymorphic sites in the catechol-O-methyl transferase, COMT, gene: an IVS 1 SNP (dbSNP rs737865), the exon 4 functional SNP (Val158Met, dbSNP rs165688), and a downstream SNP (dbSNP rs165599). Subsequently, this haplotype was shown to be associated with lower levels of COMT cDNA derived from normal cortical brain tissue, most likely due to cis-acting element(s). As a first step toward evaluating whether this haplotype may be relevant to schizophrenia in populations other than Ashkenazi Jews, we have studied this haplotype in 38 populations representing all major regions of the world. Adding to our previous data on four polymorphic sites in the COMT gene, including the Val158Met polymorphism, we have typed the IVS 1 rs737865 and 3' rs615599 sites and also included a novel IVS 1 indel polymorphism, yielding seven-site haplotype frequencies for normal individuals in the 38 globally distributed populations, including a sample of Ashkenazi Jews. We report that the schizophrenia-associated haplotype is significantly heterogeneous in populations worldwide. The three-site, schizophrenia-associated haplotype frequencies range from 0% in South America to 37.1% in Southwest Asia, despite the fact that schizophrenia occurs at roughly equal frequency around the world. Assuming that the published associations found between the exon 4 Val158Met SNP and schizophrenia are due to linkage disequilibrium, these new haplotype data support the hypothesis of a relevant cis variant linked to the rs737865 site, possibly just upstream in the P2 promoter driving transcription of the predominant form of COMT in the brain. The previously described HindIII restriction site polymorphism, located within the P2 promoter, varies within all populations and may provide essential information in future studies of schizophrenia.

Possible epistatic role of ADH7 in the protection against alcoholism.

In recent studies of the role of the alcohol dehydrogenase genes (ADH) in alcoholism the ADH1B Arg47His polymorphism appears to affect risk via a protective effect associated with the ADH1B*47His. Here we present evidence for an additional effect from outside the Class I ADH genes, presumably from functional variation at the ADH7 gene. The protective effect is restricted to one of two haplotypes identical at ADH1B but differing at an intronic SNP at ADH7 suggesting epistasis or strong linkage disequilibrium (LD).

ALFRED: the ALelle FREquency Database. Update.

Elaboration of ALFRED (http://alfred.med.yale.edu) is being continued in two directions. One of which is developing tools for efficiently annotating the entries and checking the integrity of the data already in the database while the other is to increase the quantity and accessibility of data. Information contained in ALFRED such as, polymorphic sites, number of populations and frequency tables (one sample typed for one site) has significantly increased.

Worldwide genetic analysis of the CFTR region.

Mutations at the cystic fibrosis transmembrane conductance regulator gene (CFTR) cause cystic fibrosis, the most prevalent severe genetic disorder in individuals of European descent. We have analyzed normal allele and haplotype variation at four short tandem repeat polymorphisms (STRPs) and two single-nucleotide polymorphisms (SNPs) in CFTR in 18 worldwide population samples, comprising a total of 1,944 chromosomes. The rooted phylogeny of the SNP haplotypes was established by typing ape samples. STRP variation within SNP haplotype backgrounds was highest in most ancestral haplotypes-although, when STRP allele sizes were taken into account, differences among haplotypes became smaller. Haplotype background determines STRP diversity to a greater extent than populations do, which indicates that haplotype backgrounds are older than populations. Heterogeneity among STRPs can be understood as the outcome of differences in mutation rate and pattern. STRP sites had higher heterozygosities in Africans, although, when whole haplotypes were considered, no significant differences remained. Linkage disequilibrium (LD) shows a complex pattern not easily related to physical distance. The analysis of the fraction of possible different haplotypes not found may circumvent some of the methodological difficulties of LD measure. LD analysis showed a positive correlation with locus polymorphism, which could partly explain the unusual pattern of similar LD between Africans and non-Africans. The low values found in non-Africans may imply that the size of the modern human population that emerged "Out of Africa" may be larger than what previous LD studies suggested.

ALFRED: an allele frequency database for diverse populations and DNA polymorphisms--an update.

ALFRED (the ALelle FREquency Database) is designed to store and disseminate frequencies of alleles at human polymorphic sites for multiple populations, primarily for the population genetics and molecular anthropology communities. Currently ALFRED has information on over 180 polymorphic sites for more than 70 populations. Since our initial release of the database we have focussed on increasing the quantity and quality of data, making reciprocal links between ALFRED and other related databases, and providing useful tools to make the data more comprehensible to the end user. ALFRED is accessible from the Kidd Lab home page (http://info.med.yale. edu/genetics/kkidd/) or from ALFRED directly (http://alfred.med.yale. edu/alfred/index.asp).

Y chromosome sequence variation and the history of human populations.

Binary polymorphisms associated with the non-recombining region of the human Y chromosome (NRY) preserve the paternal genetic legacy of our species that has persisted to the present, permitting inference of human evolution, population affinity and demographic history. We used denaturing high-performance liquid chromatography (DHPLC; ref. 2) to identify 160 of the 166 bi-allelic and 1 tri-allelic site that formed a parsimonious genealogy of 116 haplotypes, several of which display distinct population affinities based on the analysis of 1062 globally representative individuals. A minority of contemporary East Africans and Khoisan represent the descendants of the most ancestral patrilineages of anatomically modern humans that left Africa between 35,000 and 89,000 years ago.

Short tandem-repeat polymorphism/alu haplotype variation at the PLAT locus: implications for modern human origins.

Two dinucleotide short tandem-repeat polymorphisms (STRPs) and a polymorphic Alu element spanning a 22-kb region of the PLAT locus on chromosome 8p12-q11.2 were typed in 1,287-1,420 individuals originating from 30 geographically diverse human populations, as well as in 29 great apes. These data were analyzed as haplotypes consisting of each of the dinucleotide repeats and the flanking Alu insertion/deletion polymorphism. The global pattern of STRP/Alu haplotype variation and linkage disequilibrium (LD) is informative for the reconstruction of human evolutionary history. Sub-Saharan African populations have high levels of haplotype diversity within and between populations, relative to non-Africans, and have highly divergent patterns of LD. Non-African populations have both a subset of the haplotype diversity present in Africa and a distinct pattern of LD. The pattern of haplotype variation and LD observed at the PLAT locus suggests a recent common ancestry of non-African populations, from a small population originating in eastern Africa. These data indicate that, throughout much of modern human history, sub-Saharan Africa has maintained both a large effective population size and a high level of population substructure. Additionally, Papua New Guinean and Micronesian populations have rare haplotypes observed otherwise only in African populations, suggesting ancient gene flow from Africa into Papua New Guinea, as well as gene flow between Melanesian and Micronesian populations.

ALFRED: a Web-accessible allele frequency database.

We present a Web-accessible database (ALFRED) that allows public access to gene frequency data for a diverse set of population samples and genetic systems. The data in ALFRED are modeled based on the experience and needs of a single laboratory, but with the expectation that the database will meet the needs of a much broader scientific community that needs population-specific gene frequency estimates. Our database currently contains data on more than 40 populations representing most major regions of the world and data on more than 150 genetic systems including SNPs, STRPs, and insertion-deletion polymorphisms. While data are not available for all population-genetic system combinations, over 2000 allele frequency tables already exist. In this paper, we enumerate the broad needs in the scientific domain, describe their significance, and describe how we have designed the database to meet those needs. We compare our database with dbSNP, the NCBI database that has a broader but overlapping purpose.

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

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

ALFRED: an allele frequency database for diverse populations and DNA polymorphisms.

We have developed a publicly accessible database (ALFRED, the ALlele FREquency Database) that catalogues allele frequency data for a wide range of population samples and DNA polymorphisms. This database is web-accessible through our laboratory (Kidd Lab) Web site: http://info.med.yale.edu/genetics/kkidd. ALFRED currently contains data on 60 populations and 156 genetic systems including single nucleotide polymorphisms (SNPs), short tandem repeat polymorphisms (STRPs), variable number of tandem repeats (VNTRs) and insertion-deletion polymorphisms. While data are not available for all population-DNA polymorphism combinations, over 2000 allele frequency tables have been entered. Our database is designed (i) to address our specific research requirements as well as broader scientific objectives; (ii) to allow researchers and interested educators to easily navigate and retrieve data of interest to them; and (iii) to integrate links to other related public databases such as dbSNP, GenBank and PubMed.

Linkage disequilibrium at the ADH2 and ADH3 loci and risk of alcoholism.

Two of the three class I alcohol dehydrogenase (ADH) genes (ADH2 and ADH3) encode known functional variants that act on alcohol with different efficiencies. Variants at both these genes have been implicated in alcoholism in some populations because allele frequencies differ between alcoholics and controls. Specifically, controls have higher frequencies of the variants with higher Vmax (ADH2*2 and ADH3*1). In samples both of alcoholics and of controls from three Taiwanese populations (Chinese, Ami, and Atayal) we found significant pairwise disequilibrium for all comparisons of the two functional polymorphisms and a third, presumably neutral, intronic polymorphism in ADH2. The class I ADH genes all lie within 80 kb on chromosome 4; thus, variants are not inherited independently, and haplotypes must be analyzed when evaluating the risk of alcoholism. In the Taiwanese Chinese we found that, only among those chromosomes containing the ADH3*1 variant (high Vmax), the proportions of chromosomes with ADH2*1 (low Vmax) and those with ADH2*2 (high Vmax) are significantly different between alcoholics and controls (P<10-5). The proportions of chromosomes with ADH3*1 and those with ADH3*2 are not significantly different between alcoholics and controls, on a constant ADH2 background (with ADH2*1, P=.83; with ADH2*2, P=.53). Thus, the observed differences in the frequency of the functional polymorphism at ADH3, between alcoholics and controls, can be accounted for by the disequilibrium with ADH2 in this population.

Population studies of polymorphisms of the serotonin transporter protein gene.

The range of allele frequency variation in humans for any locus that may have functionally important genetic variation needs to be documented. Therefore, we tested two polymorphisms at the serotonin transporter protein locus (SLC6A4) in samples from seven specific populations from five continental regions. We studied the promoter polymorphism which is reported to have functional significance and to be associated with anxiety- and depression-related phenotypes [Lesch et al., 1996: Science 274:1527-1531], and the intron 2 VNTR polymorphism [Lesch et al., 1994: J Neural Transm 95:157-162]. Allele frequencies for both systems show significant global variation, and consequently so do haplotype frequencies. Linkage disequilibrium varied among the populations, being absent in some and highly significant in others. These differences further document a large potential for population stratification in association studies of either of these SLC6A4 polymorphisms.

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.

Allele frequencies in a worldwide survey of a CA repeat in the first intron of the CFTR gene.

A dinucleotide CA repeat within intron 1 of the CFTR gene has recently been identified. We have determined the allele frequencies of this polymorphism in samples from 18 populations covering all major geographical areas, with a total of 1,816 chromosomes. When considering allele distributions, African populations presented a wider range of alleles than other populations and also presented higher expected heterozygosities. Analysis of molecular variance showed that 8.04% of the genetic variance in this locus could be attributed to differences among populations. We concluded that the polymorphism in the CA repeat in intron 1 of the CFTR gene is highly informative in populations from all geographical regions of the world. Thus, it can be applied to family studies of unknown mutations causing cystic fibrosis (CF) and can provide valuable information for genetic counseling. Moreover, its analysis should be included in the haplotypic analysis of known CF mutations.