Metabolic diversity in human populations and correlation with genetic and ancestral geographic distances

DNA polymorphic markers and self-defined ethnicity groupings are used to group individuals with shared ancient geographic ancestry. Here we studied whether ancestral relationships between individuals could be identified from metabolic screening data reported by the California newborn screening (NBS) program. NBS data includes 41 blood metabolites measured by tandem mass spectrometry from singleton babies in 17 parent-reported ethnicity groupings. Ethnicity-associated differences identified for 71% of NBS metabolites (29 of 41, Cohen's d > 0.5) showed larger differences in blood levels of acylcarnitines than of amino acids (P < 1e-4). A metabolic distance measure, developed to compare ethnic groupings based on metabolic differences, showed low positive correlation with genetic and ancient geographic distances between the groups' ancestral world populations. Several outlier group pairs were identified with larger genetic and smaller metabolic distances (Black versus White) or with smaller genetic and larger metabolic distances (Chinese versus Japanese) indicating the influence of genetic and of environmental factors on metabolism. Using machine learning, comparison of metabolic profiles between all pairs of ethnic groupings distinguished individuals with larger genetic distance (Black versus Chinese, AUC = 0.96), while genetically more similar individuals could not be separated metabolically (Hispanic versus Native American, AUC = 0.51). Additionally, we identified metabolites informative for inferring metabolic ancestry in individuals from genetically similar populations, which included biomarkers for inborn metabolic disorders (C10:1, C12:1, C3, C5OH, Leucine-Isoleucine). This work sheds new light on metabolic differences in healthy newborns in diverse populations, which could have implications for improving genetic disease screening.

A multipurpose panel of microhaplotypes for use with STR markers in casework

A small panel of highly informative loci that can be genotyped on the same equipment as the standard CODIS short tandem repeat (STR) markers has strong potential for application in forensic casework. Single nucleotide polymorphisms (SNPs) can be typed by a couple of methods on capillary electrophoresis (CE) machines and on sequencers, but the amount of information relative to the laboratory effort has hindered use of SNPs in actual casework. Insertion-deletion markers (InDels) suffer from similar problems. Microhaplotypes (MHs) are much more informative per locus but have similar technical difficulties unless they are typed by massively parallel sequencing (MPS). As forensic labs are acquiring sequencing machines, MHs become more likely to be used in casework, especially if multiplexed with STRs. Here we present the details of a multipurpose panel of 24 MHs with the highest effective number of alleles (Ae) from previous work. An augmented STR panel of 24 loci (20 CODIS markers plus four commonly typed STRs) is also considered. The Ae and ancestry informativeness (In) distributions of these two datasets are compared. The MH panel is shown to have better individualization and population distinction than the augmented CODIS STRs. We note that the 24 MHs should be better for mixture analyses than the STRs. Finally, we suggest that a commercial kit including both the standard CODIS markers and this set of 24 MH would greatly improve the discrimination power over that of current commercial assays.

North Asian population relationships in a global context

Population genetic studies of North Asian ethnic groups have focused on genetic variation of sex chromosomes and mitochondria. Studies of the extensive variation available from autosomal variation have appeared infrequently. We focus on relationships among population samples using new North Asia microhaplotype data. We combined genotypes from our laboratory on 58 microhaplotypes, distributed across 18 autosomes, on 3945 individuals from 75 populations with corresponding data extracted for 26 populations from the Thousand Genomes consortium and for 22 populations from the GenomeAsia 100 K project. A total of 7107 individuals in 122 total populations are analyzed using STRUCTURE, Principal Component Analysis, and phylogenetic tree analyses. North Asia populations sampled in Mongolia include: Buryats, Mongolians, Altai Kazakhs, and Tsaatans. Available Siberians include samples of Yakut, Khanty, and Komi Zyriane. Analyses of all 122 populations confirm many known relationships and show that most populations from North Asia form a cluster distinct from all other groups. Refinement of analyses on smaller subsets of populations reinforces the distinctiveness of North Asia and shows that the North Asia cluster identifies a region that is ancestral to Native Americans.

New Insight into the human genetic diversity in North African populations by genotyping of SNPs in DRD3 , CSMD1 and NRG1 genes

Background

The single nucleotide polymorphisms (SNPs) of the dopamine D3 receptor (DRD3), the CUB and sushi multiple domains 1 (CSMD1) and the neuregulin 1 (NRG1) genes were used to study the genetic diversity and affinity among North African populations and to examine their genetic relationships in worldwide populations.

Methods

The rs3773678, rs3732783 and rs6280 SNPs of the DRD3 gene located on chromosome 3, the rs10108270 SNP of the CSMD1 gene and the rs383632, rs385396 and rs1462906 SNPs of the NRG1 gene located on chromosome 8 were analysed in 366 individuals from seven North African populations (Libya, Kairouan, Mehdia, Sousse, Kesra, Smar and Kerkennah).

Results

The low values of F_ST indicated that only 0.27%–1.65% of the genetic variability was due to the differences between the populations. The Kairouan population has the lowest average heterozygosity among the North African populations. Haplotypes composed of the ancestral alleles ACC and ACAT were more frequent in the Kairouan population than in other North African populations. The PCA and the haplotypic analysis showed that the genetic structure of populations in North Africa was closer to that of Europeans, Admixed Americans, South Asians and East Asians. However, analysis of the rs3732783 and rs6280 SNPs revealed that the CT microhaplotype was specific to the North African population.

Conclusions

The Kairouan population exhibited a relatively low rate of genetic variability. The North African population has undergone significant gene flow but also evolutionary forces that have made it genetically distinct from other populations.

The population genetics characteristics of a 90 locus panel of microhaplotypes

Single-nucleotide polymorphisms (SNPs) and small genomic regions with multiple SNPs (microhaplotypes, MHs) are rapidly emerging as novel forensic investigative tools to assist in individual identification, kinship analyses, ancestry inference, and deconvolution of DNA mixtures. Here, we analyzed information for 90 microhaplotype loci in 4009 individuals from 79 world populations in 6 major biogeographic regions. The study included multiplex microhaplotype sequencing (mMHseq) data analyzed for 524 individuals from 16 populations and genotype data for 3485 individuals from 63 populations curated from public repositories. Analyses of the 79 populations revealed excellent characteristics for this 90-plex MH panel for various forensic applications achieving an overall average effective number of allele values (A_e) of 4.55 (range 1.04–19.27) for individualization and mixture deconvolution. Population-specific random match probabilities ranged from a low of 10^–115 to a maximum of 10^–66. Mean informativeness (I_n) for ancestry inference was 0.355 (range 0.117–0.883). 65 novel SNPs were detected in 39 of the MHs using mMHseq. Of the 3018 different microhaplotype alleles identified, 1337 occurred at frequencies > 5% in at least one of the populations studied. The 90-plex MH panel enables effective differentiation of population groupings for major biogeographic regions as well as delineation of distinct subgroupings within regions. Open-source, web-based software is available to support validation of this technology for forensic case work analysis and to tailor MH analysis for specific geographical regions.

STAT3 polymorphisms in North Africa and its implication in breast cancer

BACKGROUND

Only a few studies have investigated the association of single nucleotide polymorphisms in STAT3 gene with the susceptibility to cancer and response to chemotherapy. Our aim was to determine the allele frequencies of rs3869550, rs957971, and rs7211777 at the STAT3 gene in North African populations and compare them to 1000 genomes populations, and to investigate their relation with cancer.

METHODS

The targeted SNPs have been analyzed in six Tunisian populations and a sample of Libyans using TaqMan® Assay. The results were compared to 1000 Genomes Project population samples. Targeting of the regions encompassing the three SNPs by micro-ARN was assessed using miR databases.

RESULTS

The analysis of the 3 SNPs showed that North African populations were close to South Asians. As expected, African populations presented a significant frequency of the ancestral CCG haplotype in contrast to other populations where the fully derived TGA haplotype was more frequent. The presence and diversity of rare haplotypes at STAT3 in North African populations could have been generated by recombination between the two major haplotypes. A screening of the micro-RNA databases showed that the STAT3 region with the mutated allele of rs7211777 (G>A) could be targeted by miR hsa-miR-3606-5p, which also targets genes involved in breast cancer.

Genetic diversity of the North African population revealed by the typing of SNPs in the DRD2/ANKK1 genomic region

The dopamine - related genes, like dopamine D2 receptor (DRD2) gene and ankyrin repeat and kinase domain containing 1 (ANKK1) gene are implicated in neurological functions. Some polymorphisms of the DRD2/ANKK1 locus (TaqIA, TaqIB, TaqID) have been used to study genetic diversity and the evolution of human populations. The present investigation aims to assess the genetic diversity in seven North African populations in order to explore their genetic structure and to compare them to others worldwide populations studied for the same locus. Nine single nucleotide polymorphisms (SNPs) from the DRD2/ANKK1 locus (rs1800497 TaqIA, rs2242592, rs1124492, rs6277, rs6275, rs1079727, rs2002453, rs2234690 and rs1079597 TaqIB) were typed in 366 individuals from seven North African populations: six from Tunisia (Sousse, Smar, Kesra, Kairouan, Mehdia and Kerkennah) and one from Libya. The allelic frequencies of rs2002453 and rs2234690 were higher in the Smar population than in the other North African populations. More, the Smar population showed the lowest average heterozygosity (0.313). The principal component analysis (PCA) showed that the Smar population was clearly separated from others. Furthermore, linkage disequilibrium analysis shown a high linkage disequilibrium in the North African population and essentially in Smar population. Comparison with other world populations has shown that the heterozygosity of North African population was very close to that of the African and European populations. The PCA and the haplotypic analysis suggested the presence of an important Eurasian genetic component for the North African population. These results suggested that the Smar population was isolated from the others North Africans ones by its peculiar genetic structure because of isolation, endogamy and genetic drift. On the other hand, the North African population is characterized by a multi ancestral gene pool from Eurasia and sub-Saharan Africa due to human migration since prehistoric times.

A Distinctive Pattern of Diversity for the TAS2R38 Gene in North Africa

The TAS2R38 gene is involved in bitter taste perception. This study documents the distinctive diversity patterns in northern Africa of functional single-nucleotide polymorphisms (SNPs) rs713598 and rs1726866 at the TAS2R38 locus and places those patterns in the context of global TAS2R38 diversity. Data previously genotyped with TaqMan assay were analyzed for rs713598 and rs1726866 for 375 unrelated subjects (305 Tunisians from seven locations: Mahdia, Sousse, Kesra, Nebeur, Kairouan, Smar, and Kerkennah; plus 70 Libyans). Data were analyzed to present haplotypes and genotypes before comparison with data from worldwide populations. This study provides information about TAS2R38 diversity in a part of the world that is relatively understudied. Considering the two SNPs rs713598 and rs1726866, the CA nucleotide haplotype leading to the PV amino acid haplotype is extremely rare almost everywhere, but it is relatively frequent (between 6% and 15%) in northern Africa, where it coexists with the globally common amino acid haplotypes PA, AA, and AV. Given its higher frequency in North Africa, the authors propose the CA nucleotide haplotype as a biogeographic marker for forensic purposes.

The distinctive geographic patterns of common pigmentation variants at the OCA2 gene

Oculocutaneous Albinism type 2 (OCA2) is a gene of great interest because of genetic variation affecting normal pigmentation variation in humans. The diverse geographic patterns for variant frequencies at OCA2 have been evident but have not been systematically investigated, especially outside of Europe. Here we examine population genetic variation in and near the OCA2 gene from a worldwide perspective. The very different patterns of genetic variation found across world regions suggest strong selection effects may have been at work over time. For example, analyses involving the variants that affect pigmentation of the iris argue that the derived allele of the rs1800407 single nucleotide polymorphism, which produces a hypomorphic protein, may have contributed to the previously demonstrated positive selection in Europe for the enhancer variant responsible for light eye color. More study is needed on the relationships of the genetic variation at OCA2 to variation in pigmentation in areas beyond Europe.

Recent Selection on a Class I ADH Locus Distinguishes Southwest Asian Populations Including Ashkenazi Jews

The derived human alcohol dehydrogenase (ADH)1B*48His allele of the ADH1B Arg48His polymorphism (rs1229984) has been identified as one component of an East Asian specific core haplotype that underwent recent positive selection. Our study has been extended to Southwest Asia and additional markers in East Asia. F_st values (Sewall Wright’s fixation index) and long-range haplotype analyses identify a strong signature of selection not only in East Asian but also in Southwest Asian populations. However, except for the ADH2B*48His allele, different core haplotypes occur in Southwest Asia compared to East Asia and the extended haplotypes also differ. Thus, the ADH1B*48His allele, as part of a core haplotype of 10 kb, has undergone recent rapid increases in frequency independently in the two regions after divergence of the respective populations. Emergence of agriculture may be the common factor underlying the evident selection.

Selecting microhaplotypes optimized for different purposes

Massively parallel sequencing is transforming forensic work by allowing various useful forensic markers, such as STRPs and SNPs, to be multiplexed providing information on ancestry, individual and familial identification, phenotypes for eye/hair/skin pigmentation, and the deconvolution of mixtures. Microhaplotypes also become feasible with massively parallel sequencing, these are DNA segments (smaller than 300 nucleotides) that are selected to contain multiple SNPs unambiguously defining three or more haplotype alleles occurring at common frequencies. The physical extent of a microhaplotype can thus be covered by a single sequence read making these loci phase-known codominant genetic systems. Such microhaplotypes supply significantly more information than a single SNP can. Our efforts to develop useful sets of microhaplotypes have already identified 182 such loci that we have studied on a large number of human populations from around the world. We present various analyses on 83 populations in our ongoing study for a subset of the best microhaplotypes currently available illustrating their characteristics and potential utility for ancestry, identification, and mixture deconvolution.

Ancestry inference of 96 population samples using microhaplotypes

Microhaplotypes have become a new type of forensic marker with a great ability to identify and deconvolute mixtures because massively parallel sequencing (MPS) allows the alleles (haplotypes) of the multi-SNP loci to be determined directly for an individual. As originally defined, a microhaplotype locus is a short segment of DNA with two or more SNPs defining three or more haplotypes. The length is short enough, less than about 300 bp, that the read length of current MPS technology can produce a phase-known sequence of each chromosome of an individual. As part of the discovery phase of our studies, data on 130 microhaplotype loci with estimates of haplotype frequency data on 83 populations have been published. To provide a better picture of global allele frequency variation, we have now tested 13 more populations for 65 of the microhaplotype loci from among those with higher levels of inter-population gene frequency variation, including 8 loci not previously published. These loci provide clear distinctions among 6 biogeographic regions and provide some information distinguishing up to 10 clusters of populations.

ALFRED: An Allele Frequency Database for Microevolutionary Studies

Many kinds of microevolutionary studies require data on multiple polymorphisms in multiple populations. Increasingly, and especially for human populations, multiple research groups collect relevant data and those data are dispersed widely in the literature. ALFRED has been designed to hold data from many sources and make them available over the web. Data are assembled from multiple sources, curated, and entered into the database. Multiple links to other resources are also established by the curators. A variety of search options are available and additional geographic based interfaces are being developed. The database can serve the human anthropologic genetic community by identifying what loci are already typed on many populations thereby helping to focus efforts on a common set of markers. The database can also serve as a model for databases handling similar DNA polymorphism data for other species.

Genetic variation in Tunisia in the context of human diversity worldwide

OBJECTIVES

North Africa has a complex demographic history of migrations from within Africa, Europe, and the Middle East. However, population genetic studies, especially for autosomal genetic markers, are few relative to other world regions. We examined autosomal markers for eight Tunisian and Libyan populations in order to place them in a global context.

MATERIALS AND METHODS

Data were collected by TaqMan on 399 autosomal single nucleotide polymorphisms on 331 individuals from Tunisia and Libya. These data were combined with data on the same SNPs previously typed on 2585 individuals from 57 populations from around the world. Where meaningful, close by SNPs were combined into multiallelic haplotypes. Data were evaluated by clustering, principal components, and population tree analyses. For a subset of 102 SNPs, data from the literature on seven additional North African populations were included in analyses.

RESULTS

Average heterozygosity of the North African populations is high relative to our global samples, consistent with a complex demographic history. The Tunisian and Libyan samples form a discrete cluster in the global and regional views and can be separated from sub-Sahara, Middle East, and Europe. Within Tunisia the Nebeur and Smar are outlier groups. Across North Africa, pervasive East-West geographical patterns were not found.

DISCUSSION

Known historical migrations and invasions did not displace or homogenize the genetic variation in the region but rather enriched it. Even a small region like Tunisia contains considerable genetic diversity. Future studies across North Africa have the potential to increase our understanding of the historical demographic factors influencing the region. Am J Phys Anthropol 161:62-71, 2016. © 2016 The Authors American Journal of Physical Anthropology Published by Wiley Periodicals, Inc.

Vibratory Urticaria Associated with a Missense Variant in ADGRE2

Patients with autosomal dominant vibratory urticaria have localized hives and systemic manifestations in response to dermal vibration, with coincident degranulation of mast cells and increased histamine levels in serum. We identified a previously unknown missense substitution in ADGRE2 (also known as EMR2), which was predicted to result in the replacement of cysteine with tyrosine at amino acid position 492 (p.C492Y), as the only nonsynonymous variant cosegregating with vibratory urticaria in two large kindreds. The ADGRE2 receptor undergoes autocatalytic cleavage, producing an extracellular subunit that noncovalently binds a transmembrane subunit. We showed that the variant probably destabilizes an autoinhibitory subunit interaction, sensitizing mast cells to IgE-independent vibration-induced degranulation. (Funded by the National Institutes of Health.).

Cross-disorder genome-wide analyses suggest a complex genetic relationship between Tourette's syndrome and OCD

OBJECTIVE

Obsessive-compulsive disorder (OCD) and Tourette's syndrome are highly heritable neurodevelopmental disorders that are thought to share genetic risk factors. However, the identification of definitive susceptibility genes for these etiologically complex disorders remains elusive. The authors report a combined genome-wide association study (GWAS) of Tourette's syndrome and OCD.

METHOD

The authors conducted a GWAS in 2,723 cases (1,310 with OCD, 834 with Tourette's syndrome, 579 with OCD plus Tourette's syndrome/chronic tics), 5,667 ancestry-matched controls, and 290 OCD parent-child trios. GWAS summary statistics were examined for enrichment of functional variants associated with gene expression levels in brain regions. Polygenic score analyses were conducted to investigate the genetic architecture within and across the two disorders.

RESULTS

Although no individual single-nucleotide polymorphisms (SNPs) achieved genome-wide significance, the GWAS signals were enriched for SNPs strongly associated with variations in brain gene expression levels (expression quantitative loci, or eQTLs), suggesting the presence of true functional variants that contribute to risk of these disorders. Polygenic score analyses identified a significant polygenic component for OCD (p=2×10(-4)), predicting 3.2% of the phenotypic variance in an independent data set. In contrast, Tourette's syndrome had a smaller, nonsignificant polygenic component, predicting only 0.6% of the phenotypic variance (p=0.06). No significant polygenic signal was detected across the two disorders, although the sample is likely underpowered to detect a modest shared signal. Furthermore, the OCD polygenic signal was significantly attenuated when cases with both OCD and co-occurring Tourette's syndrome/chronic tics were included in the analysis (p=0.01).

CONCLUSIONS

Previous work has shown that Tourette's syndrome and OCD have some degree of shared genetic variation. However, the data from this study suggest that there are also distinct components to the genetic architectures of these two disorders. Furthermore, OCD with co-occurring Tourette's syndrome/chronic tics may have different underlying genetic susceptibility compared with OCD alone.

Copy number variation in obsessive-compulsive disorder and tourette syndrome: a cross-disorder study

OBJECTIVE

Obsessive-compulsive disorder (OCD) and Tourette syndrome (TS) are heritable neurodevelopmental disorders with a partially shared genetic etiology. This study represents the first genome-wide investigation of large (>500 kb), rare (<1%) copy number variants (CNVs) in OCD and the largest genome-wide CNV analysis in TS to date.

METHOD

The primary analyses used a cross-disorder design for 2,699 case patients (1,613 ascertained for OCD, 1,086 ascertained for TS) and 1,789 controls. Parental data facilitated a de novo analysis in 348 OCD trios.

RESULTS

Although no global CNV burden was detected in the cross-disorder analysis or in secondary, disease-specific analyses, there was a 3.3-fold increased burden of large deletions previously associated with other neurodevelopmental disorders (p = .09). Half of these neurodevelopmental deletions were located in a single locus, 16p13.11 (5 case patient deletions: 0 control deletions, p = .08 in the current study, p = .025 compared to published controls). Three 16p13.11 deletions were confirmed de novo, providing further support for the etiological significance of this region. The overall OCD de novo rate was 1.4%, which is intermediate between published rates in controls (0.7%) and in individuals with autism or schizophrenia (2-4%).

CONCLUSION

Several converging lines of evidence implicate 16p13.11 deletions in OCD, with weaker evidence for a role in TS. The trend toward increased overall neurodevelopmental CNV burden in TS and OCD suggests that deletions previously associated with other neurodevelopmental disorders may also contribute to these phenotypes.

Current sequencing technology makes microhaplotypes a powerful new type of genetic marker for forensics

SNPs that are molecularly very close (<10kb) will generally have extremely low recombination rates, much less than 10(-4). Multiple haplotypes will often exist because of the history of the origins of the variants at the different sites, rare recombinants, and the vagaries of random genetic drift and/or selection. Such multiallelic haplotype loci are potentially important in forensic work for individual identification, for defining ancestry, and for identifying familial relationships. The new DNA sequencing capabilities currently available make possible continuous runs of a few hundred base pairs so that we can now determine the allelic combination of multiple SNPs on each chromosome of an individual, i.e., the phase, for multiple SNPs within a small segment of DNA. Therefore, we have begun to identify regions, encompassing two to four SNPs with an extent of <200bp that define multiallelic haplotype loci. We have identified candidate regions and have collected pilot data on many candidate microhaplotype loci. Here we present 31 microhaplotype loci that have at least three alleles, have high heterozygosity, are globally informative, and are statistically independent at the population level. This study of microhaplotype loci (microhaps) provides proof of principle that such markers exist and validates their usefulness for ancestry inference, lineage-clan-family inference, and individual identification. The true value of microhaplotypes will come with sequencing methods that can establish alleles unambiguously, including disentangling of mixtures, because a single sequencing run on a single strand of DNA will encompass all of the SNPs.

Worldwide population variation and haplotype analysis at the serotonin transporter gene SLC6A4 and implications for association studies

BACKGROUND

Variation at the serotonin transporter gene, SLC6A4, has been associated with a variety of neuropsychiatric disorders and could be involved in other health-related phenotypes.

METHODS

To determine the extent of variation at SLC6A4, we genotyped 23 markers on approximately 2500 individuals from 47 global populations, including the promoter variable number tandem repeat (VNTR) and 2 single nucleotide polymorphisms (SNPs) immediately flanking its variable region (rs25531 and rs25532), the intron 2 VNTR, and 19 additional SNPs.

RESULTS

We observed several rare alleles at the promoter VNTR (some novel) and population-specific distributions of the reported functional SNPs rs25531, rs25532, and rs6355, as well as two alleles at the intron 2 VNTR. Alleles of interest at the VNTRs occurred on specific haplotype backgrounds. The repeat-number variants at the promoter VNTR and the intron 2 VNTR, as well as the putative functional SNPs, showed ethnic variation in frequencies. The more common alleles at the VNTR polymorphisms show wide geographic distributions, whereas rare alleles at both show more restricted distributions. The derived alleles at the two functional SNPs in the promoter VNTR show restricted distributions and occur primarily on different repeat number alleles.

CONCLUSIONS

Our findings illustrate significant variation worldwide at SLC6A4 and that the functionally implicated alleles at the SNPs rs25531, rs25532, and rs6355 occur on limited haplotypes and vary significantly in global distribution. Association studies at SLC6A4 cannot a priori extrapolate across populations and should account for the multiple polymorphisms with possible functional variation across this locus, rather than focusing solely on one or two polymorphisms as commonly seen.

An historical perspective on "The world-wide distribution of allele frequencies at the human dopamine D4 receptor locus"

Human population genetics is a completely different science today compared to two decades ago, at least at the empiric level. Our paper [Chang (Hum Genet 98:91-101, 1996a)] demonstrated that three different alleles were common when one considered many populations although other low frequency alleles occurred. Because previous work had been largely done on European subjects, our findings involved 36 distinct populations and showed that East Asian populations had nearly lost the 7-repeat allele, and that Native American populations had the highest frequencies of that allele globally, was a significant early empiric demonstration of the potential magnitude of population variation at important genes. There are thousands of loci tested on many of the same populations and the gene frequency pattern seen for the DRD4 7-repeat allele is seen at other loci, arguing that this pattern commonly reflects the pattern of divergence of populations and accumulated random genetic drift.

Partitioning the heritability of Tourette syndrome and obsessive compulsive disorder reveals differences in genetic architecture

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.58 (se = 0.09, p = 5.64e-12) for TS, and 0.37 (se = 0.07, p = 1.5e-07) for OCD. In addition, we conducted multiple genomic partitioning analyses to identify genomic elements that concentrate this heritability. We examined genomic architectures of TS and OCD by chromosome, MAF bin, and functional annotations. In addition, we assessed heritability for early onset and adult onset OCD. Among other notable results, we found that SNPs with a minor allele frequency of less than 5% accounted for 21% of the TS heritability and 0% of the OCD heritability. Additionally, we identified a significant contribution to TS and OCD heritability by variants significantly associated with gene expression in two regions of the brain (parietal cortex and cerebellum) for which we had available expression quantitative trait loci (eQTLs). Finally we analyzed the genetic correlation between TS and OCD, revealing a genetic correlation of 0.41 (se = 0.15, p = 0.002). These results are very close to previous heritability estimates for TS and OCD based on twin and family studies, suggesting that very little, if any, heritability is truly missing (i.e., unassayed) from TS and OCD GWAS studies of common variation. The results also indicate that there is some genetic overlap between these two phenotypically-related neuropsychiatric disorders, but suggest that the two disorders have distinct genetic architectures.

Gene variants associated with antisocial behaviour: a latent variable approach

OBJECTIVE

The aim of this study was to determine if a latent variable approach might be useful in identifying shared variance across genetic risk alleles that is associated with antisocial behaviour at age 15 years.

METHODS

Using a conventional latent variable approach, we derived an antisocial phenotype in 328 adolescents utilizing data from a 15-year follow-up of a randomized trial of a prenatal and infancy nurse-home visitation programme in Elmira, New York. We then investigated, via a novel latent variable approach, 450 informative genetic polymorphisms in 71 genes previously associated with antisocial behaviour, drug use, affiliative behaviours and stress response in 241 consenting individuals for whom DNA was available. Haplotype and Pathway analyses were also performed.

RESULTS

Eight single-nucleotide polymorphisms (SNPs) from eight genes contributed to the latent genetic variable that in turn accounted for 16.0% of the variance within the latent antisocial phenotype. The number of risk alleles was linearly related to the latent antisocial variable scores. Haplotypes that included the putative risk alleles for all eight genes were also associated with higher latent antisocial variable scores. In addition, 33 SNPs from 63 of the remaining genes were also significant when added to the final model. Many of these genes interact on a molecular level, forming molecular networks. The results support a role for genes related to dopamine, norepinephrine, serotonin, glutamate, opioid and cholinergic signalling as well as stress response pathways in mediating susceptibility to antisocial behaviour.

CONCLUSIONS

This preliminary study supports use of relevant behavioural indicators and latent variable approaches to study the potential 'co-action' of gene variants associated with antisocial behaviour. It also underscores the cumulative relevance of common genetic variants for understanding the aetiology of complex behaviour. If replicated in future studies, this approach may allow the identification of a 'shared' variance across genetic risk alleles associated with complex neuropsychiatric dimensional phenotypes using relatively small numbers of well-characterized research participants.

Haplotype structure and positive selection at TLR1

Toll-like receptor 1, when dimerized with Toll-like receptor 2, is a cell surface receptor that, upon recognition of bacterial lipoproteins, activates the innate immune system. Variants in TLR1 associate with the risk of a variety of medical conditions and diseases, including sepsis, leprosy, tuberculosis, and others. The foremost of these is rs5743618 c.2079T>G(p.(Ile602Ser)), the derived allele of which is associated with reduced risk of sepsis, leprosy, and other diseases. Interestingly, 602Ser, which shows signatures of selection, inhibits TLR1 surface trafficking and subsequent activation of NFκB upon recognition of a ligand. This suggests that reduced TLR1 activity may be beneficial for human health. To better understand TLR1 variation and its link to human health, we have typed all 7 high-frequency missense variants (>5% in at least one population) along with 17 other variants in and around TLR1 in 2548 individuals from 56 populations from around the globe. We have also found additional signatures of selection on missense variants not associated with rs5743618, suggesting that there may be multiple functional alleles under positive selection in this gene.

Introducing the Forensic Research/Reference on Genetics knowledge base, FROG-kb

BACKGROUND

Online tools and databases based on multi-allelic short tandem repeat polymorphisms (STRPs) are actively used in forensic teaching, research, and investigations. The Fst value of each CODIS marker tends to be low across the populations of the world and most populations typically have all the common STRP alleles present diminishing the ability of these systems to discriminate ethnicity. Recently, considerable research is being conducted on single nucleotide polymorphisms (SNPs) to be considered for human identification and description. However, online tools and databases that can be used for forensic research and investigation are limited.

METHODS

The back end DBMS (Database Management System) for FROG-kb is Oracle version 10. The front end is implemented with specific code using technologies such as Java, Java Servlet, JSP, JQuery, and GoogleCharts.

RESULTS

We present an open access web application, FROG-kb (Forensic Research/Reference on Genetics-knowledge base, http://frog.med.yale.edu), that is useful for teaching and research relevant to forensics and can serve as a tool facilitating forensic practice. The underlying data for FROG-kb are provided by the already extensively used and referenced ALlele FREquency Database, ALFRED (http://alfred.med.yale.edu). In addition to displaying data in an organized manner, computational tools that use the underlying allele frequencies with user-provided data are implemented in FROG-kb. These tools are organized by the different published SNP/marker panels available. This web tool currently has implemented general functions possible for two types of SNP panels, individual identification and ancestry inference, and a prediction function specific to a phenotype informative panel for eye color.

CONCLUSION

The current online version of FROG-kb already provides new and useful functionality. We expect FROG-kb to grow and expand in capabilities and welcome input from the forensic community in identifying datasets and functionalities that will be most helpful and useful. Thus, the structure and functionality of FROG-kb will be revised in an ongoing process of improvement. This paper describes the state as of early June 2012.

Variant in the glucokinase regulatory protein (GCKR) gene is associated with fatty liver in obese children and adolescents

Recently, the single nucleotide polymorphism (SNP) identified as rs1260326, in the glucokinase regulatory protein (GCKR), was associated with hypertriglyceridemia in adults. Because accumulation of triglycerides in hepatocytes represents the hallmark of steatosis, we aimed to investigate whether this variant might be associated with fatty liver (hepatic fat content, HFF%). Moreover, because recently rs738409 in the PNPLA3 and rs2854116 in the APOC3 were associated with fatty liver, we explored how the GCKR SNP and these two variants jointly influence hepatosteatosis. We studied 455 obese children and adolescents (181 Caucasians, 139 African Americans, and 135 Hispanics). All underwent an oral glucose tolerance test and fasting lipoprotein subclasses measurement by proton nuclear magnetic resonance. A subset of 142 children underwent a fast gradient magnetic resonance imaging to measure the HFF%. The rs1260326 was associated with elevated triglycerides (Caucasians P = 0.00014; African Americans P = 0.00417), large very low-density lipoprotein (VLDL) (Caucasians P = 0.001; African Americans, P = 0.03), and with fatty liver (Caucasians P = 0.034; African Americans P = 0.00002; and Hispanics P = 0.016). The PNPLA3, but not the APOC3 rs2854116 SNP, was associated with fatty liver but not with triglyceride levels. There was a joint effect between the PNPLA3 and GCKR SNPs, explaining 32% of HFF% variance in Caucasians (P = 0.00161), 39.0% in African Americans (P = 0.00000496), and 15% in Hispanics (P = 0.00342).

CONCLUSION

The rs1260326 in GCKR is associated with hepatic fat accumulation along with large VLDL and triglyceride levels. GCKR and PNPLA3 act together to convey susceptibility to fatty liver in obese youths.

ALFRED: an allele frequency resource for research and teaching

ALFRED (http://alfred.med.yale.edu) is a free, web accessible, curated compilation of allele frequency data on DNA sequence polymorphisms in anthropologically defined human populations. Currently, ALFRED has allele frequency tables on over 663 400 polymorphic sites; 170 of them have frequency tables for more than 100 different population samples. In ALFRED, a population may have multiple samples with each ‘sample’ consisting of many individuals on which an allele frequency is based. There are 3566 population samples from 710 different populations with allele frequency tables on at least one polymorphism. Fifty of those population samples have allele frequency data for over 650 000 polymorphisms. Records also have active links to relevant resources (dbSNP, PharmGKB, OMIM, Ethnologue, etc.). The flexible search options and data display and download capabilities available through the web interface allow easy access to the large quantity of high-quality data in ALFRED.

Single nucleotide polymorphisms and haplotypes in Native American populations

Autosomal DNA polymorphisms can provide new information and understanding of both the origins of and relationships among modern Native American populations. At the same time that autosomal markers can be highly informative, they are also susceptible to ascertainment biases in the selection of the markers to use. Identifying markers that can be used for ancestry inference among Native American populations can be considered separate from identifying markers to further the quest for history. In the current study, we are using data on nine Native American populations to compare the results based on a large haplotype-based dataset with relatively small independent sets of single nucleotide polymorphisms. We are interested in what types of limited datasets an individual laboratory might be able to collect are best for addressing two different questions of interest. First, how well can we differentiate the Native American populations and/or infer ancestry by assigning an individual to her population(s) of origin? Second, how well can we infer the historical/evolutionary relationships among Native American populations and their Eurasian origins? We conclude that only a large comprehensive dataset involving multiple autosomal markers on multiple populations will be able to answer both questions; different small sets of markers are able to answer only one or the other of these questions. Using our largest dataset, we see a general increasing distance from Old World populations from North to South in the New World except for an unexplained close relationship between our Maya and Quechua samples.

Diversification of the ADH1B gene during expansion of modern humans

A variant allele, ADH1B*48His, also known as ADH1B*2, at the human Alcohol Dehydrogenase 1B gene (ADH1B) is strongly associated with alcoholism in some populations and has an unusual geographic distribution. Strong evidence implies selection has increased the frequency of this allele in some East Asian populations but does not fully explain its geographic pattern. We have studied haplotypes of 10 single nucleotide polymorphisms (SNPs) and two short tandem repeat polymorphisms (STRPs) in the ADH1B region in 2,206 individuals from a worldwide set of populations. These SNPs and STRPs define nine common haplogroups most of which have distinct geographic patterns. The haplogroups H5 and H6, both with the derived ADH1B*48His allele, appear restricted to the Middle East and East Asia, respectively. The positively selected H7 is derived from H6 by a new regulatory region variant defining SNP rs3811801 restricted to East Asia. Age estimates of the haplogroups based on the STRPs also agree with the time of the migration events estimated by other studies. H7 is estimated to have expanded recently, around 2,800 years ago, and ancient DNA samples from North China confirm its presence about that time. The dating of the H7 expansion may help understand the selective force on the ADH1B gene.

Analyses of a set of 128 ancestry informative single-nucleotide polymorphisms in a global set of 119 population samples

Background

Using DNA to determine an individual's ancestry from among human populations is generally interesting and useful for many purposes, including admixture mapping, controlling for population structure in disease or trait association studies and forensic ancestry inference. However, to estimate ancestry, including possible admixture within an individual, as well as heterogeneity within a group of individuals, allele frequencies are necessary for what are believed to be the contributing populations. For this purpose, panels of ancestry informative markers (AIMs) have been developed.

Results

We are presenting our work on one such panel, composed of 128 ancestry informative single-nucleotide polymorphisms (AISNPs) already proposed in the literature. Compared to previous studies of these AISNPs, we have studied three times the number of individuals (4,871) in three times as many population samples (119). We have validated this panel for many ancestry assignment and admixture studies, especially those that were the rationale for the original selection of the 128 SNPs: African Americans and Mexican Americans. At the same time, the limitations of the panel for distinguishing ancestry and quantifying admixture among Eurasian populations are noted.

Conclusion

We demonstrate the simultaneous importance of the specific set of population samples and their relative sample sizes in the use of the structure program to determine which groups cluster together and consequently influence the ability of a marker panel to infer ancestry. We demonstrate the strengths and weaknesses of this particular panel of AISNPs in a global context.

Aggressive behavior, related conduct problems, and variation in genes affecting dopamine turnover

A number of dopamine-related genes have been implicated in the etiology of violent behavior and conduct problems. Of these genes, the ones that code for the enzymes that influence the turnover of dopamine (DA) have received the most attention. In this study, we investigated 12 genetic polymorphisms in four genes involved with DA functioning (COMT, MAOA and MAOB, and DbetaH) in 179 incarcerated male Russian adolescents and two groups of matched controls: boys without criminal records referred to by their teachers as (a) "troubled-behavior-free" boys, n=182; and (b) "troubled-behavior" boys, n=60. The participants were classified as (1) being incarcerated or not, (2) having the DSM-IV diagnosis of conduct disorder (CD) or not, and (3) having committed violent or nonviolent crimes (for the incarcerated individuals only). The findings indicate that, although no single genetic variant in any of the four genes differentiated individuals in the investigated groups, various linear combinations (i.e., haplotypes) and nonlinear combinations (i.e., interactions between variants within and across genes) of genetic variants resulted in informative and robust classifications for two of the three groupings. These combinations of genetic variants differentiated individuals in incarceration vs. nonincarcerated and CD vs. no-CD groups; no informative combinations were established consistently for the grouping by crime within the incarcerated individuals. This study underscores the importance of considering multiple rather than single markers within candidate genes and their additive and interactive combinations, both with themselves and with nongenetic indicators, while attempting to understand the genetic background of such complex behaviors as serious conduct problems.

SNPs for a universal individual identification panel

An efficient method to uniquely identify every individual would have value in quality control and sample tracking of large collections of cell lines or DNA as is now often the case with whole genome association studies. Such a method would also be useful in forensics. SNPs represent the best markers for such purposes. We have developed a globally applicable resource of 92 SNPs for individual identification (IISNPs) with extremely low probabilities of any two unrelated individuals from anywhere in the world having identical genotypes. The SNPs were identified by screening over 500 likely/candidate SNPs on samples of 44 populations representing the major regions of the world. All 92 IISNPs have an average heterozygosity [0.4 and the F(st) values are all\0.06 on our 44 populations making these a universally applicable panel irrespective of ethnicity or ancestry. No significant linkage disequilibrium (LD) occurs for all unique pairings of 86 of the 92 IISNPs (median LD = 0.011) in all of the 44 populations. The remaining 6 IISNPs show strong LD in most of the 44 populations for a small subset (7) of the unique pairings in which they occur due to close linkage. 45 of the 86 SNPs are spread across the 22 human autosomes and show very loose or no genetic linkage with each other. These 45 IISNPs constitute an excellent panel for individual identification including paternity testing with associated probabilities of individual genotypes less than 10(-15), smaller than achieved with the current panels of forensic markers. This panel also improves on an interim panel of 40 IISNPs previously identified using 40 population samples. The unlinked status of the subset of 45 SNPs we have identified also makes them useful for situations involving close biological relationships. Comparisons with random sets of SNPs illustrate the greater discriminating power, efficiency, and more universal applicability of this IISNP panel to populations around the world. The full set of 86 IISNPs that do not show LD can be used to provide even smaller genotype match probabilities in the range of 10(-31)-10(-35) based on the 44 population samples studied.

The distribution and most recent common ancestor of the 17q21 inversion in humans

The polymorphic inversion on 17q21, sometimes called the microtubular associated protein tau (MAPT) inversion, is an approximately 900 kb inversion found primarily in Europeans and Southwest Asians. We have identified 21 SNPs that act as markers of the inverted, i.e., H2, haplotype. The inversion is found at the highest frequencies in Southwest Asia and Southern Europe (frequencies of approximately 30%); elsewhere in Europe, frequencies vary from < 5%, in Finns, to 28%, in Orcadians. The H2 inversion haplotype also occurs at low frequencies in Africa, Central Asia, East Asia, and the Americas, though the East Asian and Amerindian alleles may be due to recent gene flow from Europe. Molecular evolution analyses indicate that the H2 haplotype originally arose in Africa or Southwest Asia. Though the H2 inversion has many fixed differences across the approximately 900 kb, short tandem repeat polymorphism data indicate a very recent date for the most recent common ancestor, with dates ranging from 13,600 to 108,400 years, depending on assumptions and estimation methods. This estimate range is much more recent than the 3 million year age estimated by Stefansson et al. in 2005.

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.

Origin and dispersal of atypical aldehyde dehydrogenase ALDH2487Lys

The East Asian respond with a marked facial flushing and mild to moderate symptoms of intoxication after drinking the amounts of alcohol that has no detectable effect on European. The alcohol sensitivity in Orientals is due to a delayed oxidation of acetaldehyde by an atypical aldehyde dehydrogenase ALDH2487Lys, which is resulted from a structural mutation in gene ALDH2. The atypical ALDH2487Lys allele has been associated with various phenotypic statuses, such as protective against alcohol dependence and the risk of alcohol-related digestive tract cancers. Here, we have examined this SNP, adjacent four non-coding SNPs, and one downstream STRP on ALDH2 gene, in total of 1072 unrelated healthy individuals from 14 Chinese populations and 130 Indian individuals. Five major haplotypes based on five SNPs across the ALDH2 gene 40 kb were found in all East Asian populations. The frequencies of the ancestral haplotype GCCTG and the East Asian special haplotype GCCTA containing the atypical ALDH2487Lys allele were 44.8% and 14.9%, respectively. The frequency of the atypical ALDH2487Lys allele or the East Asian specific haplotype GCCTA is high in Yunnan, South coastal, east coastal of China, and decreased gradually toward inland China, West, Northwest and North China. Combined with demographic history in East Asian, our results showed that the presence of ALDH2487Lys allele in peripheral regions of China might be the results of historical migration events from China to these regions. The origin of ALDH2487Lys could be possibly traced back to ancient Pai-Yuei tribe in South China.

Haplotype evolution of SLITRK1, a candidate gene for Gilles de la Tourette syndrome

Gilles de la Tourette syndrome (GTS) is a complex disorder with a clear genetic component but no clearly identified genes with variation of etiologic relevance. Various candidate regions and genes show some evidence of affecting risk, though clearly not all patients/families can be explained by any one of them. Resequencing one candidate gene, SLITRK1, has identified four new variants. Including them, we have typed over 2,300 normal individuals from 44 populations for 11 SNPs spanning the gene. The unusual global pattern seen is that one non-ancestral haplotype is the single most common haplotype worldwide. Other haplotypes appear to result from accumulation of mutations with no evidence of historical recombination. Although there is no evidence of selection, the haplotype frequency variation seen around the world will need to be considered in any future association studies of this locus with GTS or any other neuropsychiatric disorder.

Ethnic related selection for an ADH Class I variant within East Asia

Background

The alcohol dehydrogenases (ADH) are widely studied enzymes and the evolution of the mammalian gene cluster encoding these enzymes is also well studied. Previous studies have shown that the ADH1B*47His allele at one of the seven genes in humans is associated with a decrease in the risk of alcoholism and the core molecular region with this allele has been selected for in some East Asian populations. As the frequency of ADH1B*47His is highest in East Asia, and very low in most of the rest of the world, we have undertaken more detailed investigation in this geographic region.

Methodology/Principal Findings

Here we report new data on 30 SNPs in the ADH7 and Class I ADH region in samples of 24 populations from China and Laos. These populations cover a wide geographic region and diverse ethnicities. Combined with our previously published East Asian data for these SNPs in 8 populations, we have typed populations from all of the 6 major linguistic phyla (Altaic including Korean-Japanese and inland Altaic, Sino-Tibetan, Hmong-Mien, Austro-Asiatic, Daic, and Austronesian). The ADH1B genotyping data are strongly related to ethnicity. Only some eastern ethnic phyla or subphyla (Korean-Japanese, Han Chinese, Hmong-Mien, Daic, and Austronesian) have a high frequency of ADH1B*47His. ADH1B haplotype data clustered the populations into linguistic subphyla, and divided the subphyla into eastern and western parts. In the Hmong-Mien and Altaic populations, the extended haplotype homozygosity (EHH) and relative EHH (REHH) tests for the ADH1B core were consistent with selection for the haplotype with derived SNP alleles. In the other ethnic phyla, the core showed only a weak signal of selection at best.

Conclusions/Significance

The selection distribution is more significantly correlated with the frequency of the derived ADH1B regulatory region polymorphism than the derived amino-acid altering allele ADH1B*47His. Thus, the real focus of selection may be the regulatory region. The obvious ethnicity-related distributions of ADH1B diversities suggest the existence of some culture-related selective forces that have acted on the ADH1B region.

Evidence of positive selection on a class I ADH locus

The alcohol dehydrogenase (ADH) family of enzymes catalyzes the reversible oxidation of alcohol to acetaldehyde. Seven ADH genes exist in a segment of ~370 kb on 4q21. Products of the three class I ADH genes that share 95% sequence identity are believed to play the major role in the first step of ethanol metabolism. Because the common belief that selection has operated at the ADH1B*47His allele in East Asian populations lacks direct biological or statistical evidence, we used genomic data to test the hypothesis. Data consisted of 54 single-nucleotide polymorphisms (SNPs) across the ADH clusters in a global sampling of 42 populations. Both the F(st) statistic and the long-range haplotype (LRH) test provided positive evidence of selection in several East Asian populations. The ADH1B Arg47His functional polymorphism has the highest F(st) of the 54 SNPs in the ADH cluster, and it is significantly above the mean F(st) of 382 presumably neutral sites tested on the same 42 population samples. The LRH test that uses cores including that site and extending on both sides also gives significant evidence of positive selection in some East Asian populations for a specific haplotype carrying the ADH1B*47His allele. Interestingly, this haplotype is present at a high frequency in only some East Asian populations, whereas the specific allele also exists in other East Asian populations and in the Near East and Europe but does not show evidence of selection with use of the LRH test. Although the ADH1B*47His allele conveys a well-confirmed protection against alcoholism, that modern phenotypic manifestation does not easily translate into a positive selective force, and the nature of that selective force, in the past and/or currently, remains speculative.

Candidate SNPs for a universal individual identification panel

Single nucleotide polymorphisms (SNPs) are likely in the near future to have a fundamental role both in human identification and description. However, because allele frequencies can vary greatly among populations, a critical issue is the population genetics underlying calculation of the probabilities of unrelated individuals having identical multi-locus genotypes. Here we report on progress in identifying SNPs that show little allele frequency variation among a worldwide sample of 40 populations, i.e., have a low F(st), while remaining highly informative. Such markers have match probabilities that are nearly uniform irrespective of population and become candidates for a universally applicable individual identification panel applicable in forensics and paternity testing. They are also immediately useful for efficient sample identification/tagging in large biomedical, association, and epidemiologic studies. Using our previously described strategy for both identifying and characterizing such SNPs (Kidd et al. in Forensic Sci Int 164:20-32, 2006), we have now screened a total of 432 SNPs likely a priori to have high heterozygosity and low allele frequency variation and from these have selected the markers with the lowest F(st) in our set of 40 populations to produce a panel of 40 low F(st), high heterozygosity SNPs. Collectively these SNPs give average match probabilities of less than 10(-16) in most of the 40 populations and less than 10(-14) in all but one small isolated population; the range is 2.02 x 10(-17) to 1.29 x 10(-13). These 40 SNPs constitute excellent candidates for the global forensic community to consider for a universally applicable SNP panel for human identification. The relative ease with which these markers could be identified also provides a cautionary lesson for investigations of possible balancing selection.

ALFRED: an allele frequency database for microevolutionary studies

Many kinds of microevolutionary studies require data on multiple polymorphisms in multiple populations. Increasingly, and especially for human populations, multiple research groups collect relevant data and those data are dispersed widely in the literature. ALFRED has been designed to hold data from many sources and make them available over the web. Data are assembled from multiple sources, curated, and entered into the database. Multiple links to other resources are also established by the curators. A variety of search options are available and additional geographic based interfaces are being developed. The database can serve the human anthropologic genetic community by identifying what loci are already typed on many populations thereby helping to focus efforts on a common set of markers. The database can also serve as a model for databases handling similar DNA polymorphism data for other species.

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.

Intra- and interpopulation genotype reconstruction from tagging SNPs

The optimal method to be used for tSNP selection, the applicability of a reference LD map to unassayed populations, and the scalability of these methods to genome-wide analysis, all remain subjects of debate. We propose novel, scalable matrix algorithms that address these issues and we evaluate them on genotypic data from 38 populations and four genomic regions (248 SNPs typed for approximately 2000 individuals). We also evaluate these algorithms on a second data set consisting of genotypes available from the HapMap database (1336 SNPs for four populations) over the same genomic regions. Furthermore, we test these methods in the setting of a real association study using a publicly available family data set. The algorithms we use for tSNP selection and unassayed SNP reconstruction do not require haplotype inference and they are, in principle, scalable even to genome-wide analysis. Moreover, they are greedy variants of recently developed matrix algorithms with provable performance guarantees. Using a small set of carefully selected tSNPs, we achieve very good reconstruction accuracy of "untyped" genotypes for most of the populations studied. Additionally, we demonstrate in a quantitative manner that the chosen tSNPs exhibit substantial transferability, both within and across different geographic regions. Finally, we show that reconstruction can be applied to retrieve significant SNP associations with disease, with important genotyping savings.

Conservative evolution in duplicated genes of the primate Class I ADH cluster

Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.

Considerable haplotype diversity within the 23kb encompassing the ADH7 gene

BACKGROUND

Of the seven known human alcohol dehydrogenase (ADH) genes, the non-liver expressed ADH7 gene codes for the enzyme with the highest maximal activity for ethanol. Previous study from our laboratory has suggested that ADH7 has an epistatic role for protection against alcoholism based on a single ADH7 SNP.

METHODS

We have now studied seven SNPs, additional populations for the SNP previously examined, and six more new SNPs, across 23 kb of ADH7 in 38 population samples originating from different geographical regions of the world.

RESULTS

The overall linkage disequilibrium is moderate to strong across this region even though considerable 7-SNP haplotype diversity is observed. This uncommonly high haplotype diversity is explained by high LD within each "half," the three upstream SNPs and the four downstream SNPs, but near randomization between the "halves." This division significantly simplified the haplotype pattern: only four major haplotypes account for almost all chromosomes in all populations in each "half."

CONCLUSIONS

The low linkage disequilibrium between these two "halves" suggests multiple recombination(s) have occurred in this region, specifically, within intron 7. The absence of strong LD between the functional variation in ADH1B that is strongly associated with alcoholism and any of the variation in ADH7 supports the genetic independence of ADH7 in association studies. Thus, the previously observed epistatic effect of ADH7 cannot be explained by its linkage disequilibrium with a causative factor in ADH1B.

HAPLOT: a graphical comparison of haplotype blocks, tagSNP sets and SNP variation for multiple populations

Understanding of human variation relevant to association studies can benefit from population comparison, especially comparing populations in the same geographical region. Variations in linkage disequilibrium patterns, in tagSNP sets, and in SNP heterozygosities among populations can be used to infer the evolutionary pattern. We present here a win32 system based Perl/Tk application for visual comparisons of these variations in different populations.

AVAILABILITY

The application package is available at http://info.med.yale.edu/genetics/kkidd/programs.html

CONTACT

sheng.gu@yale.edu.

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.