A meta-analysis of 120 246 individuals identifies 18 new loci for fibrinogen concentration.

Genome-wide association studies have previously identified 23 genetic loci associated with circulating fibrinogen concentration. These studies used HapMap imputation and did not examine the X-chromosome. 1000 Genomes imputation provides better coverage of uncommon variants, and includes indels. We conducted a genome-wide association analysis of 34 studies imputed to the 1000 Genomes Project reference panel and including ∼120 000 participants of European ancestry (95 806 participants with data on the X-chromosome). Approximately 10.7 million single-nucleotide polymorphisms and 1.2 million indels were examined. We identified 41 genome-wide significant fibrinogen loci; of which, 18 were newly identified. There were no genome-wide significant signals on the X-chromosome. The lead variants of five significant loci were indels. We further identified six additional independent signals, including three rare variants, at two previously characterized loci: FGB and IRF1. Together the 41 loci explain 3% of the variance in plasma fibrinogen concentration.

The Genetic Overlap Between Hair and Eye Color.

We identified the genetic variants for eye color by Genome-Wide Association Study (GWAS) in a Dutch Caucasian family-based population sample and examined the genetic correlation between hair and eye color using data from unrelated participants from the Netherlands Twin Register. With the Genome-wide Complex Trait Analysis software package, we found strong genetic correlations between various combinations of hair and eye colors. The strongest positive correlations were found for blue eyes with blond hair (0.87) and brown eyes with dark hair (0.71), whereas blue eyes with dark hair and brown eyes with blond hair showed the strongest negative correlations (-0.64 and -0.94, respectively). Red hair with green/hazel eyes showed the weakest correlation (-0.14). All analyses were corrected for age and sex, and we explored the effects of correcting for principal components (PCs) that represent ancestry and describe the genetic stratification of the Netherlands. When including the first three PCs as covariates, the genetic correlations between the phenotypes disappeared. This is not unexpected since hair and eye colors strongly indicate the ancestry of an individual. This makes it difficult to separate the effects of population stratification and the true genetic effects of variants on these particular phenotypes.

DeepSAGE reveals genetic variants associated with alternative polyadenylation and expression of coding and non-coding transcripts.

Many disease-associated variants affect gene expression levels (expression quantitative trait loci, eQTLs) and expression profiling using next generation sequencing (NGS) technology is a powerful way to detect these eQTLs. We analyzed 94 total blood samples from healthy volunteers with DeepSAGE to gain specific insight into how genetic variants affect the expression of genes and lengths of 3'-untranslated regions (3'-UTRs). We detected previously unknown cis-eQTL effects for GWAS hits in disease- and physiology-associated traits. Apart from cis-eQTLs that are typically easily identifiable using microarrays or RNA-sequencing, DeepSAGE also revealed many cis-eQTLs for antisense and other non-coding transcripts, often in genomic regions containing retrotransposon-derived elements. We also identified and confirmed SNPs that affect the usage of alternative polyadenylation sites, thereby potentially influencing the stability of messenger RNAs (mRNA). We then combined the power of RNA-sequencing with DeepSAGE by performing a meta-analysis of three datasets, leading to the identification of many more cis-eQTLs. Our results indicate that DeepSAGE data is useful for eQTL mapping of known and unknown transcripts, and for identifying SNPs that affect alternative polyadenylation. Because of the inherent differences between DeepSAGE and RNA-sequencing, our complementary, integrative approach leads to greater insight into the molecular consequences of many disease-associated variants.

Multiethnic meta-analysis of genome-wide association studies in >100 000 subjects identifies 23 fibrinogen-associated Loci but no strong evidence of a causal association between circulating fibrinogen and cardiovascular disease.

BACKGROUND: Estimates of the heritability of plasma fibrinogen concentration, an established predictor of cardiovascular disease, range from 34% to 50%. Genetic variants so far identified by genome-wide association studies explain only a small proportion (<2%) of its variation. METHODS AND RESULTS: We conducted a meta-analysis of 28 genome-wide association studies including >90 000 subjects of European ancestry, the first genome-wide association meta-analysis of fibrinogen levels in 7 studies in blacks totaling 8289 samples, and a genome-wide association study in Hispanics totaling 1366 samples. Evaluation for association of single-nucleotide polymorphisms with clinical outcomes included a total of 40 695 cases and 85 582 controls for coronary artery disease, 4752 cases and 24 030 controls for stroke, and 3208 cases and 46 167 controls for venous thromboembolism. Overall, we identified 24 genome-wide significant (P<5×10(-8)) independent signals in 23 loci, including 15 novel associations, together accounting for 3.7% of plasma fibrinogen variation. Gene-set enrichment analysis highlighted key roles in fibrinogen regulation for the 3 structural fibrinogen genes and pathways related to inflammation, adipocytokines, and thyrotrophin-releasing hormone signaling. Whereas lead single-nucleotide polymorphisms in a few loci were significantly associated with coronary artery disease, the combined effect of all 24 fibrinogen-associated lead single-nucleotide polymorphisms was not significant for coronary artery disease, stroke, or venous thromboembolism. CONCLUSIONS: We identify 23 robustly associated fibrinogen loci, 15 of which are new. Clinical outcome analysis of these loci does not support a causal relationship between circulating levels of fibrinogen and coronary artery disease, stroke, or venous thromboembolism.

Loci affecting gamma-glutamyl transferase in adults and adolescents show age × SNP interaction and cardiometabolic disease associations.

Serum gamma-glutamyl transferase (GGT) activity is a marker of liver disease which is also prospectively associated with the risk of all-cause mortality, cardiovascular disease, type 2 diabetes and cancers. We have discovered novel loci affecting GGT in a genome-wide association study (rs1497406 in an intergenic region of chromosome 1, P = 3.9 × 10(-8); rs944002 in C14orf73 on chromosome 14, P = 4.7 × 10(-13); rs340005 in RORA on chromosome 15, P = 2.4 × 10(-8)), and a highly significant heterogeneity between adult and adolescent results at the GGT1 locus on chromosome 22 (maximum P(HET) = 5.6 × 10(-12) at rs6519520). Pathway analysis of significant and suggestive single-nucleotide polymorphism associations showed significant overlap between genes affecting GGT and those affecting common metabolic and inflammatory diseases, and identified the hepatic nuclear factor (HNF) family as controllers of a network of genes affecting GGT. Our results reinforce the disease associations of GGT and demonstrate that control by the GGT1 locus varies with age.

Effect of polygenic risk scores on depression in childhood trauma.

BACKGROUND: Research on gene × environment interaction in major depressive disorder (MDD) has thus far primarily focused on candidate genes, although genetic effects are known to be polygenic. AIMS: To test whether the effect of polygenic risk scores on MDD is moderated by childhood trauma. METHOD: The study sample consisted of 1645 participants with a DSM-IV diagnosis of MDD and 340 screened controls from The Netherlands. Chronic or remitted episodes (severe MDD) were present in 956 participants. The occurrence of childhood trauma was assessed with the Childhood Trauma Interview and the polygenic risk scores were based on genome-wide meta-analysis results from the Psychiatric Genomics Consortium. RESULTS: The polygenic risk scores and childhood trauma independently affected MDD risk, and evidence was found for interaction as departure from both multiplicativity and additivity, indicating that the effect of polygenic risk scores on depression is increased in the presence of childhood trauma. The interaction effects were similar in predicting all MDD risk and severe MDD risk, and explained a proportion of variation in MDD risk comparable to the polygenic risk scores themselves. CONCLUSIONS: The interaction effect found between polygenic risk scores and childhood trauma implies that (1) studies on direct genetic effect on MDD gain power by focusing on individuals exposed to childhood trauma, and that (2) individuals with both high polygenic risk scores and exposure to childhood trauma are particularly at risk for developing MDD.

Item-Level Genome-Wide Association Study of the Alcohol Use Disorders Identification Test in Three Population-Based Cohorts.

OBJECTIVE: Genome-wide association studies (GWASs) of the Alcohol Use Disorders Identification Test (AUDIT), a 10-item screen for alcohol use disorder (AUD), have elucidated novel loci for alcohol consumption and misuse. However, these studies also revealed that GWASs can be influenced by numerous biases (e.g., measurement error, selection bias), which may have led to inconsistent genetic correlations between alcohol involvement and AUD, as well as paradoxically negative genetic correlations between alcohol involvement and psychiatric disorders and/or medical conditions. The authors used genomic structural equation modeling to elucidate the genetics of alcohol consumption and problematic consequences of alcohol use as measured by AUDIT. METHODS: To explore these unexpected differences in genetic correlations, the authors conducted the first item-level and the largest GWAS of AUDIT items (N=160,824) and applied a multivariate framework to mitigate previous biases. RESULTS: The authors identified novel patterns of similarity (and dissimilarity) among the AUDIT items and found evidence of a correlated two-factor structure at the genetic level ("consumption" and "problems," rg=0.80). Moreover, by applying empirically derived weights to each of the AUDIT items, the authors constructed an aggregate measure of alcohol consumption that was strongly associated with alcohol dependence (rg=0.67), moderately associated with several other psychiatric disorders, and no longer positively associated with health and positive socioeconomic outcomes. Lastly, by conducting polygenic analyses in three independent cohorts that differed in their ascertainment and prevalence of AUD, the authors identified novel genetic associations between alcohol consumption, alcohol misuse, and health. CONCLUSIONS: This work further emphasizes the value of AUDIT for both clinical and genetic studies of AUD and the importance of using multivariate methods to study genetic associations that are more closely related to AUD.

Genome-wide analysis of DNA methylation in buccal cells: a study of monozygotic twins and mQTLs.

BACKGROUND: DNA methylation arrays are widely used in epigenome-wide association studies and methylation quantitative trait locus (mQTL) studies. Here, we performed the first genome-wide analysis of monozygotic (MZ) twin correlations and mQTLs on data obtained with the Illumina MethylationEPIC BeadChip (EPIC array) and compared the performance of the EPIC array to the Illumina HumanMethylation450 BeadChip (HM450 array) for buccal-derived DNA. RESULTS: Good-quality EPIC data were obtained for 102 buccal-derived DNA samples from 49 MZ twin pairs (mean age = 7.5 years, range = 1-10). Differences between MZ twins in the cellular content of buccal swabs were a major driver for differences in their DNA methylation profiles, highlighting the importance to adjust for cellular composition in DNA methylation studies of buccal-derived DNA. After adjusting for cellular composition, the genome-wide mean correlation (r) between MZ twins was 0.21 for the EPIC array, and cis mQTL analysis in 84 twins identified 1,296,323 significant associations (FDR 5%), encompassing 33,749 methylation sites and 616,029 genetic variants. MZ twin correlations were slightly larger (p < 2.2 × 10-16) for novel EPIC probes (N = 383,066, mean r = 0.22) compared to probes that are also present on HM450 (N = 406,822, mean r = 0.20). In line with this observation, a larger percentage of novel EPIC probes was associated with genetic variants (novel EPIC probes with significant mQTL 4.7%, HM450 probes with mQTL 3.9%, p < 2.2 × 10-16). Methylation sites with a large MZ correlation and sites associated with mQTLs were most strongly enriched in epithelial cell DNase I hypersensitive sites (DHSs), enhancers, and histone mark H3K4me3. CONCLUSIONS: We conclude that the contribution of familial factors to individual differences in DNA methylation and the effect of mQTLs are larger for novel EPIC probes, especially those within regulatory elements connected to active regions specific to the investigated tissue.

A method to customize population-specific arrays for genome-wide association testing.

As an example of optimizing population-specific genotyping assays using a whole-genome sequence reference set, we detail the approach that followed to design the Axiom-NL array which is characterized by an improved imputation backbone based on the Genome of the Netherlands (GoNL) reference sequence and, compared with earlier arrays, a more comprehensive inclusion of SNPs on chromosomes X, Y, and the mitochondria. Common variants on the array were selected to be compatible with the Illumina Psych Array and the Affymetrix UK Biobank Axiom array. About 3.5% of the array (23 977 markers) represents SNPs from the GWAS catalog, including SNPs at FTO, APOE, Ion-channels, killer-cell immunoglobulin-like receptors, and HLA. Around 26 000 markers associated with common psychiatric disorders are included, as well as 6705 markers suggested to be associated with fertility and twinning. The platform can thus be used for risk profiling, detection of new variants, as well as ancestry determination. Results of coverage tests in 249 unrelated subjects with GoNL-based sequence data show that after imputation with 1000G as a reference, the median concordance between original and imputed genotypes is above 98%. The median imputation quality R2 for MAF thresholds of 0.001, 0.01, 0.05, and >0.05 are 0.05, 0.28, 0.80, 0.99, respectively, for the 1000G imputed SNPs, with a similar quality for the autosomes and X chromosome, showing a good genome-wide coverage for association studies after imputation.

Genetic Overlap Between Schizophrenia and Developmental Psychopathology: Longitudinal and Multivariate Polygenic Risk Prediction of Common Psychiatric Traits During Development.

Background: Several nonpsychotic psychiatric disorders in childhood and adolescence can precede the onset of schizophrenia, but the etiology of this relationship remains unclear. We investigated to what extent the association between schizophrenia and psychiatric disorders in childhood is explained by correlated genetic risk factors. Methods: Polygenic risk scores (PRS), reflecting an individual's genetic risk for schizophrenia, were constructed for 2588 children from the Netherlands Twin Register (NTR) and 6127 from the Avon Longitudinal Study of Parents And Children (ALSPAC). The associations between schizophrenia PRS and measures of anxiety, depression, attention deficit hyperactivity disorder (ADHD), and oppositional defiant disorder/conduct disorder (ODD/CD) were estimated at age 7, 10, 12/13, and 15 years in the 2 cohorts. Results were then meta-analyzed, and a meta-regression analysis was performed to test differences in effects sizes over, age and disorders. Results: Schizophrenia PRS were associated with childhood and adolescent psychopathology. Meta-regression analysis showed differences in the associations over disorders, with the strongest association with childhood and adolescent depression and a weaker association for ODD/CD at age 7. The associations increased with age and this increase was steepest for ADHD and ODD/CD. Genetic correlations varied between 0.10 and 0.25. Conclusion: By optimally using longitudinal data across diagnoses in a multivariate meta-analysis this study sheds light on the development of childhood disorders into severe adult psychiatric disorders. The results are consistent with a common genetic etiology of schizophrenia and developmental psychopathology as well as with a stronger shared genetic etiology between schizophrenia and adolescent onset psychopathology.

The low single nucleotide polymorphism heritability of plasma and saliva cortisol levels.

Cortisol is an important stress hormone affected by a variety of biological and environmental factors, such as the circadian rhythm, exercise and psychological stress. Cortisol is mostly measured using blood or saliva samples. A number of genetic variants have been found to contribute to cortisol levels with these methods. While the effects of several specific single genetic variants is known, the joint genome-wide contribution to cortisol levels is unclear. Our aim was to estimate the amount of cortisol variance explained by common single nucleotide polymorphisms, i.e. the SNP heritability, using a variety of cortisol measures, cohorts and analysis approaches. We analyzed morning plasma (n=5705) and saliva levels (n=1717), as well as diurnal saliva levels (n=1541), in the Rotterdam Study using genomic restricted maximum likelihood estimation. Additionally, linkage disequilibrium score regression was fitted on the results of genome-wide association studies (GWAS) performed by the CORNET consortium on morning plasma cortisol (n=12,597) and saliva cortisol (n=7703). No significant SNP heritability was detected for any cortisol measure, sample or analysis approach. Point estimates ranged from 0% to 9%. Morning plasma cortisol in the CORNET cohorts, the sample with the most power, had a 6% [95%CI: 0-13%] SNP heritability. The results consistently suggest a low SNP heritability of these acute and short-term measures of cortisol. The low SNP heritability may reflect the substantial environmental and, in particular, situational component of these cortisol measures. Future GWAS will require very large sample sizes. Alternatively, more long-term cortisol measures such as hair cortisol samples are needed to discover further genetic pathways regulating cortisol concentrations.

Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk.

The timing of puberty is a highly polygenic childhood trait that is epidemiologically associated with various adult diseases. Using 1000 Genomes Project-imputed genotype data in up to ∼370,000 women, we identify 389 independent signals (P < 5 × 10-8) for age at menarche, a milestone in female pubertal development. In Icelandic data, these signals explain ∼7.4% of the population variance in age at menarche, corresponding to ∼25% of the estimated heritability. We implicate ∼250 genes via coding variation or associated expression, demonstrating significant enrichment in neural tissues. Rare variants near the imprinted genes MKRN3 and DLK1 were identified, exhibiting large effects when paternally inherited. Mendelian randomization analyses suggest causal inverse associations, independent of body mass index (BMI), between puberty timing and risks for breast and endometrial cancers in women and prostate cancer in men. In aggregate, our findings highlight the complexity of the genetic regulation of puberty timing and support causal links with cancer susceptibility.

Comparison of HapMap and 1000 Genomes Reference Panels in a Large-Scale Genome-Wide Association Study.

An increasing number of genome-wide association (GWA) studies are now using the higher resolution 1000 Genomes Project reference panel (1000G) for imputation, with the expectation that 1000G imputation will lead to the discovery of additional associated loci when compared to HapMap imputation. In order to assess the improvement of 1000G over HapMap imputation in identifying associated loci, we compared the results of GWA studies of circulating fibrinogen based on the two reference panels. Using both HapMap and 1000G imputation we performed a meta-analysis of 22 studies comprising the same 91,953 individuals. We identified six additional signals using 1000G imputation, while 29 loci were associated using both HapMap and 1000G imputation. One locus identified using HapMap imputation was not significant using 1000G imputation. The genome-wide significance threshold of 5×10-8 is based on the number of independent statistical tests using HapMap imputation, and 1000G imputation may lead to further independent tests that should be corrected for. When using a stricter Bonferroni correction for the 1000G GWA study (P-value < 2.5×10-8), the number of loci significant only using HapMap imputation increased to 4 while the number of loci significant only using 1000G decreased to 5. In conclusion, 1000G imputation enabled the identification of 20% more loci than HapMap imputation, although the advantage of 1000G imputation became less clear when a stricter Bonferroni correction was used. More generally, our results provide insights that are applicable to the implementation of other dense reference panels that are under development.

Identification of context-dependent expression quantitative trait loci in whole blood.

Genetic risk factors often localize to noncoding regions of the genome with unknown effects on disease etiology. Expression quantitative trait loci (eQTLs) help to explain the regulatory mechanisms underlying these genetic associations. Knowledge of the context that determines the nature and strength of eQTLs may help identify cell types relevant to pathophysiology and the regulatory networks underlying disease. Here we generated peripheral blood RNA-seq data from 2,116 unrelated individuals and systematically identified context-dependent eQTLs using a hypothesis-free strategy that does not require previous knowledge of the identity of the modifiers. Of the 23,060 significant cis-regulated genes (false discovery rate (FDR) ≤ 0.05), 2,743 (12%) showed context-dependent eQTL effects. The majority of these effects were influenced by cell type composition. A set of 145 cis-eQTLs depended on type I interferon signaling. Others were modulated by specific transcription factors binding to the eQTL SNPs.

Disease variants alter transcription factor levels and methylation of their binding sites.

Most disease-associated genetic variants are noncoding, making it challenging to design experiments to understand their functional consequences. Identification of expression quantitative trait loci (eQTLs) has been a powerful approach to infer the downstream effects of disease-associated variants, but most of these variants remain unexplained. The analysis of DNA methylation, a key component of the epigenome, offers highly complementary data on the regulatory potential of genomic regions. Here we show that disease-associated variants have widespread effects on DNA methylation in trans that likely reflect differential occupancy of trans binding sites by cis-regulated transcription factors. Using multiple omics data sets from 3,841 Dutch individuals, we identified 1,907 established trait-associated SNPs that affect the methylation levels of 10,141 different CpG sites in trans (false discovery rate (FDR) < 0.05). These included SNPs that affect both the expression of a nearby transcription factor (such as NFKB1, CTCF and NKX2-3) and methylation of its respective binding site across the genome. Trans methylation QTLs effectively expose the downstream effects of disease-associated variants.

Genome-wide Association for Major Depression Through Age at Onset Stratification: Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium.

BACKGROUND: Major depressive disorder (MDD) is a disabling mood disorder, and despite a known heritable component, a large meta-analysis of genome-wide association studies revealed no replicable genetic risk variants. Given prior evidence of heterogeneity by age at onset in MDD, we tested whether genome-wide significant risk variants for MDD could be identified in cases subdivided by age at onset. METHODS: Discovery case-control genome-wide association studies were performed where cases were stratified using increasing/decreasing age-at-onset cutoffs; significant single nucleotide polymorphisms were tested in nine independent replication samples, giving a total sample of 22,158 cases and 133,749 control subjects for subsetting. Polygenic score analysis was used to examine whether differences in shared genetic risk exists between earlier and adult-onset MDD with commonly comorbid disorders of schizophrenia, bipolar disorder, Alzheimer's disease, and coronary artery disease. RESULTS: We identified one replicated genome-wide significant locus associated with adult-onset (>27 years) MDD (rs7647854, odds ratio: 1.16, 95% confidence interval: 1.11-1.21, p = 5.2 × 10-11). Using polygenic score analyses, we show that earlier-onset MDD is genetically more similar to schizophrenia and bipolar disorder than adult-onset MDD. CONCLUSIONS: We demonstrate that using additional phenotype data previously collected by genetic studies to tackle phenotypic heterogeneity in MDD can successfully lead to the discovery of genetic risk factor despite reduced sample size. Furthermore, our results suggest that the genetic susceptibility to MDD differs between adult- and earlier-onset MDD, with earlier-onset cases having a greater genetic overlap with schizophrenia and bipolar disorder.

Blood lipids influence DNA methylation in circulating cells.

BACKGROUND: Cells can be primed by external stimuli to obtain a long-term epigenetic memory. We hypothesize that long-term exposure to elevated blood lipids can prime circulating immune cells through changes in DNA methylation, a process that may contribute to the development of atherosclerosis. To interrogate the causal relationship between triglyceride, low-density lipoprotein (LDL) cholesterol, and high-density lipoprotein (HDL) cholesterol levels and genome-wide DNA methylation while excluding confounding and pleiotropy, we perform a stepwise Mendelian randomization analysis in whole blood of 3296 individuals. RESULTS: This analysis shows that differential methylation is the consequence of inter-individual variation in blood lipid levels and not vice versa. Specifically, we observe an effect of triglycerides on DNA methylation at three CpGs, of LDL cholesterol at one CpG, and of HDL cholesterol at two CpGs using multivariable Mendelian randomization. Using RNA-seq data available for a large subset of individuals (N = 2044), DNA methylation of these six CpGs is associated with the expression of CPT1A and SREBF1 (for triglycerides), DHCR24 (for LDL cholesterol) and ABCG1 (for HDL cholesterol), which are all key regulators of lipid metabolism. CONCLUSIONS: Our analysis suggests a role for epigenetic priming in end-product feedback control of lipid metabolism and highlights Mendelian randomization as an effective tool to infer causal relationships in integrative genomics data.

Polygenic risk scores for schizophrenia and bipolar disorder predict creativity.

We tested whether polygenic risk scores for schizophrenia and bipolar disorder would predict creativity. Higher scores were associated with artistic society membership or creative profession in both Icelandic (P = 5.2 × 10(-6) and 3.8 × 10(-6) for schizophrenia and bipolar disorder scores, respectively) and replication cohorts (P = 0.0021 and 0.00086). This could not be accounted for by increased relatedness between creative individuals and those with psychoses, indicating that creativity and psychosis share genetic roots.

A Dutch family with 'familial cortical tremor with epilepsy'. Clinical characteristics and exclusion of linkage to chromosome 8q23.3-q24.1.

OBJECTIVES: To describe the clinical characteristics of a large Dutch family with cortical tremor with epilepsy (FCTE) and to test for genetic linkage of FCTE to chromosome 8q23.3-q24.1. BACKGROUND: FCTE is an idiopathic generalised epilepsy of adult onset with autosomal dominant inheritance. It is characterised by kinesiogenic tremor and myoclonus of the limbs, generalised seizures, and electrophysiological findings consistent with cortical reflex myoclonus. Genetic analysis has been performed in five Japanese families. In all families, linkage was shown to chromosome 8q23.3-q24.1. METHODS: Clinical and electrophysiological data of a four-generation family, suspected of autosomal dominant inherited FCTE, were collected and linkage analysis was performed. Results Clinical and electrophysiological findings were consistent with a diagnosis of FCTE. Of 41 relatives examined, 13 subjects were considered to be definitely affected, three were probably affected and ten were unaffected. In 15 relatives, the diagnosis could not be established. Linkage to chromosome 8q23.3-q24.1 was excluded. CONCLUSIONS: In this family with autosomal dominant FCTE, specific clinical and electrophysiological features were identified. Exclusion of linkage to chromosome 8q23.3-q24.1 indicates that genetic heterogeneity exists for FCTE.

Global similarity with local differences in linkage disequilibrium between the Dutch and HapMap-CEU populations.

The HapMap project has facilitated the selection of tagging single nucleotide polymorphisms (tagSNPs) for genome-wide association studies (GWAS) under the assumption that linkage disequilibrium (LD) in the HapMap populations is similar to the populations under investigation. Earlier reports support this assumption, although in most of these studies only a few loci were evaluated. We compared pair-wise LD and LD block structure across autosomes between the Dutch population and the CEU-HapMap reference panel. The impact of sampling distribution on the estimation of LD blocks was studied by bootstrapping. A high Pearson correlation (genome-wide; 0.93) between pair-wise r(2) for the Dutch and the CEU populations was found, indicating that tagSNPs from the CEU-HapMap panel capture common variation in the Dutch population. However, some genomic regions exhibited, significantly lower correlation than the genome-wide estimate. This might decrease the validity of HapMap tagSNPs in these regions and the power of GWAS. The LD block structure differed considerably between the Dutch and CEU-HapMap populations. This was not explained by demographic differences between the CEU and Dutch samples, as testing for population stratification was not significant. We also found that sampling variation had a large effect on the estimation of LD blocks, as shown by the bootstrapping analysis. Thus, in small samples, most of the observed differences in LD blocks between populations are most likely the result of sampling variation. This poor concordance in LD block structure suggests that large samples are required for robust estimations of local LD block structure in populations.