Experimental studies suggest that matrix metalloproteinases (MMPs) are involved in the degradation of amyloid beta (Abeta) protein which plays a key role in the pathogenesis of Alzheimer's disease. Whether MMPs are associated with changes in beta amyloid levels in humans remains unclear. We related common haplotypes within the gene encoding MMP-3 with plasma levels of Abeta(1-40) and Abeta(1-42) in 1621 non-demented participants of the population-based Rotterdam Study. In non-demented persons, plasma Abeta concentration reflect levels of Abeta in the brain. DNA was genotyped for polymorphisms 1187 (5A6A, rs3025058), 2092A>G (rs522616), 9775T>A (rs563096) and 6658T>C (rs3025066) and haplotypes reconstructed (coding from 1187 (5A6A), 2092A>G, 9775T>A and 6658T>C: haplotype 1=5A-A-T-T, haplotype 2=6A-G-T-T, haplotype 3=6A-A-T-T, haplotype 4=6A-A-A-T and haplotype 5=5A-A-T-C). Then the associations of these haplotypes with plasma Abeta(1-40) and Abeta(1-42) levels were assessed using the program Haplo.Stats. Compared with haplotype 1, haplotype 4 was associated with significantly lower levels of plasma Abeta(1-40) (beta=-8.04, 95% CI (-14.79; -1.28), p=0.02) after adjusting for age and sex. Haplotype 2 was associated with significantly higher levels of plasma Abeta(1-42) (beta=3.70, 95% CI (1.75; 5.65), p=0.0002). Our observations suggest that variation in the gene encoding MMP-3 is associated with changes in amyloid beta levels in humans. Factors modulating secretion or activity of MMP-3 may have the potential to influence the amount of Abeta concentration and deposition in the brain.

An inverse association between birth weight and the risk of developing type 2 diabetes (T2D) in adulthood has been reported. This association may be explained by common genetic variants related to insulin secretion and resistance, since insulin is the most important growth factor in fetal life. The objective of this study was to examine whether T2D gene polymorphism TCF7L2 rs7903146 is associated with growth patterns from fetal life until infancy.

To identify genetic loci influencing central obesity and fat distribution, we performed a meta-analysis of 16 genome-wide association studies (GWAS, N = 38,580) informative for adult waist circumference (WC) and waist-hip ratio (WHR). We selected 26 SNPs for follow-up, for which the evidence of association with measures of central adiposity (WC and/or WHR) was strong and disproportionate to that for overall adiposity or height. Follow-up studies in a maximum of 70,689 individuals identified two loci strongly associated with measures of central adiposity; these map near TFAP2B (WC, P = 1.9x10(-11)) and MSRA (WC, P = 8.9x10(-9)). A third locus, near LYPLAL1, was associated with WHR in women only (P = 2.6x10(-8)). The variants near TFAP2B appear to influence central adiposity through an effect on overall obesity/fat-mass, whereas LYPLAL1 displays a strong female-only association with fat distribution. By focusing on anthropometric measures of central obesity and fat distribution, we have identified three loci implicated in the regulation of human adiposity.

Age-related macular degeneration (AMD) is the leading cause of irreversible visual loss in the developed countries and is caused by both environmental and genetic factors. A recent study (Tuo et al., PNAS) reported an association between AMD and a single nucleotide polymorphism (SNP) (rs3793784) in the ERCC6 (NM_000124) gene. The risk allele also increased ERCC6 expression. ERCC6 is involved in DNA repair and mutations in ERCC6 cause Cockayne syndrome (CS). Amongst others, photosensitivity and pigmentary retinopathy are hallmarks of CS.

The objective of this study was to examine the relationship between common genetic variation of the ESR2 gene and osteoarthritis.

Migraine is a common episodic neurological disorder, typically presenting with recurrent attacks of severe headache and autonomic dysfunction. Apart from rare monogenic subtypes, no genetic or molecular markers for migraine have been convincingly established. We identified the minor allele of rs1835740 on chromosome 8q22.1 to be associated with migraine (P = 5.38 × 10⁻⁹, odds ratio = 1.23, 95% CI 1.150-1.324) in a genome-wide association study of 2,731 migraine cases ascertained from three European headache clinics and 10,747 population-matched controls. The association was replicated in 3,202 cases and 40,062 controls for an overall meta-analysis P value of 1.69 × 10⁻¹¹ (odds ratio = 1.18, 95% CI 1.127-1.244). rs1835740 is located between MTDH (astrocyte elevated gene 1, also known as AEG-1) and PGCP (encoding plasma glutamate carboxypeptidase). In an expression quantitative trait study in lymphoblastoid cell lines, transcript levels of the MTDH were found to have a significant correlation to rs1835740 (P = 3.96 × 10⁻⁵, permuted threshold for genome-wide significance 7.7 × 10⁻⁵. To our knowledge, our data establish rs1835740 as the first genetic risk factor for migraine.

Major depression (MD) is one of the most prevalent psychiatric disorders and a leading cause of loss in work productivity. A combination of genetic and environmental risk factors probably contributes to MD. We present data from a genome-wide association study revealing a neuron-specific neutral amino acid transporter (SLC6A15) as a susceptibility gene for MD. Risk allele carrier status in humans and chronic stress in mice were associated with a downregulation of the expression of this gene in the hippocampus, a brain region implicated in the pathophysiology of MD. The same polymorphisms also showed associations with alterations in hippocampal volume and neuronal integrity. Thus, decreased SLC6A15 expression, due to genetic or environmental factors, might alter neuronal circuits related to the susceptibility for MD. Our convergent data from human genetics, expression studies, brain imaging, and animal models suggest a pathophysiological mechanism for MD that may be accessible to drug targeting.

Migraine is a common neurological disorder with a genetically complex background. This paper describes a meta-analysis of genome-wide association (GWA) studies on migraine, performed by the Dutch-Icelandic migraine genetics (DICE) consortium, which brings together six population-based European migraine cohorts with a total sample size of 10,980 individuals (2446 cases and 8534 controls). A total of 32 SNPs showed marginal evidence for association at a P-value<10(-5). The best result was obtained for SNP rs9908234, which had a P-value of 8.00 × 10(-8). This top SNP is located in the nerve growth factor receptor (NGFR) gene. However, this SNP did not replicate in three cohorts from the Netherlands and Australia. Of the other 31 SNPs, 18 SNPs were tested in two replication cohorts, but none replicated. In addition, we explored previously identified candidate genes in the meta-analysis data set. This revealed a modest gene-based significant association between migraine and the metadherin (MTDH) gene, previously identified in the first clinic-based GWA study (GWAS) for migraine (Bonferroni-corrected gene-based P-value=0.026). This finding is consistent with the involvement of the glutamate pathway in migraine. Additional research is necessary to further confirm the involvement of glutamate.

Information about growth of kidney structures in early life is limited. In a population-based prospective cohort study, from foetal life onwards, we constructed reference curves for kidney growth from the third trimester of pregnancy until early childhood, using data from 1,158 healthy children. Kidney size, defined as length, width, depth and volume, was measured in the third trimester of pregnancy and at the postnatal ages of 6 months and 24 months. Analyses were based on more than 2,500 kidney measurements. In the third trimester of pregnancy and at 6 months of age all kidney measurements were larger in boys than in girls. At 24 months of age, these gender differences were only significant for left kidney structures and right kidney length. Both groups showed trends towards smaller left kidney measurements than right kidney measurements at all ages. Gender-specific reference curves based on post-conceptional and postnatal ages were constructed for left and right kidney length, width, depth and volume. We concluded that kidney size is influenced by age and gender. Left kidney size tended to be smaller than right kidney size, except for kidney length. The reference curves can be used for assessing kidney structures by ultrasound in foetal life and early childhood.

SIRT1 has pleiotropic metabolic functions. We investigated whether SIRT1 genetic variation is associated with obesity.

Systemic sclerosis (SSc) is an orphan, complex, inflammatory disease affecting the immune system and connective tissue. SSc stands out as a severely incapacitating and life-threatening inflammatory rheumatic disease, with a largely unknown pathogenesis. We have designed a two-stage genome-wide association study of SSc using case-control samples from France, Italy, Germany, and Northern Europe. The initial genome-wide scan was conducted in a French post quality-control sample of 564 cases and 1,776 controls, using almost 500 K SNPs. Two SNPs from the MHC region, together with the 6 loci outside MHC having at least one SNP with a P<10(-5) were selected for follow-up analysis. These markers were genotyped in a post-QC replication sample of 1,682 SSc cases and 3,926 controls. The three top SNPs are in strong linkage disequilibrium and located on 6p21, in the HLA-DQB1 gene: rs9275224, P = 9.18×10(-8), OR = 0.69, 95% CI [0.60-0.79]; rs6457617, P = 1.14×10(-7) and rs9275245, P = 1.39×10(-7). Within the MHC region, the next most associated SNP (rs3130573, P = 1.86×10(-5), OR = 1.36 [1.18-1.56]) is located in the PSORS1C1 gene. Outside the MHC region, our GWAS analysis revealed 7 top SNPs (P<10(-5)) that spanned 6 independent genomic regions. Follow-up of the 17 top SNPs in an independent sample of 1,682 SSc and 3,926 controls showed associations at PSORS1C1 (overall P = 5.70×10(-10), OR:1.25), TNIP1 (P = 4.68×10(-9), OR:1.31), and RHOB loci (P = 3.17×10(-6), OR:1.21). Because of its biological relevance, and previous reports of genetic association at this locus with connective tissue disorders, we investigated TNIP1 expression. A markedly reduced expression of the TNIP1 gene and also its protein product were observed both in lesional skin tissue and in cultured dermal fibroblasts from SSc patients. Furthermore, TNIP1 showed in vitro inhibitory effects on inflammatory cytokine-induced collagen production. The genetic signal of association with TNIP1 variants, together with tissular and cellular investigations, suggests that this pathway has a critical role in regulating autoimmunity and SSc pathogenesis.

Sudden cardiac death (SCD) continues to be one of the leading causes of mortality worldwide, with an annual incidence estimated at 250,000-300,000 in the United States and with the vast majority occurring in the setting of coronary disease. We performed a genome-wide association meta-analysis in 1,283 SCD cases and >20,000 control individuals of European ancestry from 5 studies, with follow-up genotyping in up to 3,119 SCD cases and 11,146 controls from 11 European ancestry studies, and identify the BAZ2B locus as associated with SCD (P = 1.8×10(-10)). The risk allele, while ancestral, has a frequency of ~1.4%, suggesting strong negative selection and increases risk for SCD by 1.92-fold per allele (95% CI 1.57-2.34). We also tested the role of 49 SNPs previously implicated in modulating electrocardiographic traits (QRS, QT, and RR intervals). Consistent with epidemiological studies showing increased risk of SCD with prolonged QRS/QT intervals, the interval-prolonging alleles are in aggregate associated with increased risk for SCD (P = 0.006).

Runs of homozygosity (ROH) are extended tracts of adjacent homozygous single nucleotide polymorphisms (SNPs) that are more common in unrelated individuals than previously thought. It has been proposed that estimating ROH on a genome-wide level, by making use of the genome-wide single nucleotide polymorphism (SNP) data, will enable to indentify recessive variants underlying complex traits. Here, we examined ROH larger than 1.5 Mb individually and in combination for association with survival in 5974 participants of the Rotterdam Study. In addition, we assessed the role of overall homozygosity, expressed as a percentage of the autosomal genome that is in ROH longer than 1.5 Mb, on survival during a mean follow-up period of 12 years. None of these measures of homozygosity was associated with survival to old age.

In the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological investigations.

Obesity is highly heritable. Genetic variants showing robust associations with obesity traits have been identified through genome-wide association studies. We investigated whether a composite score representing healthy diet modifies associations of these variants with obesity traits. Totally, 32 body mass index (BMI)- and 14 waist-hip ratio (WHR)-associated single nucleotide polymorphisms were genotyped, and genetic risk scores (GRS) were calculated in 18 cohorts of European ancestry (n = 68 317). Diet score was calculated based on self-reported intakes of whole grains, fish, fruits, vegetables, nuts/seeds (favorable) and red/processed meats, sweets, sugar-sweetened beverages and fried potatoes (unfavorable). Multivariable adjusted, linear regression within each cohort followed by inverse variance-weighted, fixed-effects meta-analysis was used to characterize: (a) associations of each GRS with BMI and BMI-adjusted WHR and (b) diet score modification of genetic associations with BMI and BMI-adjusted WHR. Nominally significant interactions (P = 0.006-0.04) were observed between the diet score and WHR-GRS (but not BMI-GRS), two WHR loci (GRB14 rs10195252; LYPLAL1 rs4846567) and two BMI loci (LRRN6C rs10968576; MTIF3 rs4771122), for the respective BMI-adjusted WHR or BMI outcomes. Although the magnitudes of these select interactions were small, our data indicated that associations between genetic predisposition and obesity traits were stronger with a healthier diet. Our findings generate interesting hypotheses; however, experimental and functional studies are needed to determine their clinical relevance.

Genome-wide association studies (GWAS) have been successful in identifying loci associated with a wide range of complex human traits and diseases. Up to now, the majority of GWAS have focused on European populations. However, the inclusion of other ethnic groups as well as admixed populations in GWAS studies is rapidly rising following the pressing need to extrapolate findings to non-European populations and to increase statistical power. In this paper, we describe the methodological steps surrounding genetic data generation, quality control, study design and analytical procedures needed to run GWAS in the multiethnic and highly admixed Generation R Study, a large prospective birth cohort in Rotterdam, the Netherlands. Furthermore, we highlight a number of practical considerations and alternatives pertinent to the quality control and analysis of admixed GWAS data.

Variants associated with blood lipid levels may be population-specific. To identify low-frequency variants associated with this phenotype, population-specific reference panels may be used. Here we impute nine large Dutch biobanks (~35,000 samples) with the population-specific reference panel created by the Genome of The Netherlands Project and perform association testing with blood lipid levels. We report the discovery of five novel associations at four loci (P value <6.61 × 10(-4)), including a rare missense variant in ABCA6 (rs77542162, p.Cys1359Arg, frequency 0.034), which is predicted to be deleterious. The frequency of this ABCA6 variant is 3.65-fold increased in the Dutch and its effect (βLDL-C=0.135, βTC=0.140) is estimated to be very similar to those observed for single variants in well-known lipid genes, such as LDLR.

Several genes implicated in high-density lipoprotein (HDL) metabolism have been reported to be associated with age-related macular degeneration (AMD). Furthermore, HDL transport the two carotenoids, lutein and zeaxanthin, which are highly suspected to play a key-role in the protection against AMD. The objective is to confirm the associations of HDL-related loci with AMD and to assess their associations with plasma lutein and zeaxanthin concentrations.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: β = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Chronic obstructive pulmonary disease (COPD) independently associates with an increased risk of coronary artery disease (CAD), but it has not been fully investigated whether this co-morbidity involves shared pathophysiological mechanisms. To identify potential common pathways across the two diseases, we tested all recently published single nucleotide polymorphisms (SNPs) associated with human lung function (spirometry) for association with carotid intima-media thickness (cIMT) in 3,378 subjects with multiple CAD risk factors, and for association with CAD in a case-control study of 5,775 CAD cases and 7,265 controls. SNPs rs2865531, located in the CFDP1 gene, and rs9978142, located in the KCNE2 gene, were significantly associated with CAD. In addition, SNP rs9978142 and SNP rs3995090 located in the HTR4 gene, were associated with average and maximal cIMT measures. Genetic risk scores combining the most robustly spirometry-associated SNPs from the literature were modestly associated with CAD, (odds ratio (OR) (95% confidence interval (CI95) = 1.06 (1.03, 1.09); P-value = 1.5 × 10(-4), per allele). In conclusion, our study suggests that some genetic loci implicated in determining human lung function also influence cIMT and susceptibility to CAD. The present results should help elucidate the molecular underpinnings of the co-morbidity observed across COPD and CAD.