Mutations create variation in the population, fuel evolution and cause genetic diseases. Current knowledge about de novo mutations is incomplete and mostly indirect. Here we analyze 11,020 de novo mutations from the whole genomes of 250 families. We show that de novo mutations in the offspring of older fathers are not only more numerous but also occur more frequently in early-replicating, genic regions. Functional regions exhibit higher mutation rates due to CpG dinucleotides and show signatures of transcription-coupled repair, whereas mutation clusters with a unique signature point to a new mutational mechanism. Mutation and recombination rates independently associate with nucleotide diversity, and regional variation in human-chimpanzee divergence is only partly explained by heterogeneity in mutation rate. Finally, we provide a genome-wide mutation rate map for medical and population genetics applications. Our results provide new insights and refine long-standing hypotheses about human mutagenesis.

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.

To assess candidate genes for association with osteoarthritis (OA) and identify promising genetic factors and, secondarily, to assess the candidate gene approach in OA.

Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.

Approximately 1.5 billion people worldwide are overweight or affected by obesity, and are at risk of developing type 2 diabetes, cardiovascular disease and related metabolic and inflammatory disturbances. Although the mechanisms linking adiposity to associated clinical conditions are poorly understood, recent studies suggest that adiposity may influence DNA methylation, a key regulator of gene expression and molecular phenotype. Here we use epigenome-wide association to show that body mass index (BMI; a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci with P < 1 × 10-7, range P = 9.2 × 10-8 to 6.0 × 10-46; n = 10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find that methylation loci are enriched for functional genomic features in multiple tissues (P < 0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P < 9.0 × 10-6, range P = 5.5 × 10-6 to 6.1 × 10-35, n = 1,785 samples). The methylation loci identify genes involved in lipid and lipoprotein metabolism, substrate transport and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future development of type 2 diabetes (relative risk per 1 standard deviation increase in methylation risk score: 2.3 (2.07-2.56); P = 1.1 × 10-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type 2 diabetes and other adverse clinical consequences of obesity.

DNA repair deficiencies have been postulated to play a role in the development and progression of cardiovascular disease (CVD). The hypothesis is that DNA damage accumulating with age may induce cell death, which promotes formation of unstable plaques. Defects in DNA repair mechanisms may therefore increase the risk of CVD events. We examined whether the joints effect of common genetic variants in 5 DNA repair pathways may influence the risk of CVD events.

Hypertension is a major risk factor for a spectrum of cardiovascular diseases (CVD), including myocardial infarction, sudden death, and stroke. In the US, over 65 million people have high blood pressure and a large proportion of these individuals are prescribed antihypertensive medications. Although large long-term clinical trials conducted in the last several decades have identified a number of effective antihypertensive treatments that reduce the risk of future clinical complications, responses to therapy and protection from cardiovascular events vary among individuals.

Looking young for one's age has been a desire since time immemorial. This desire is attributable to the belief that appearance reflects health and fecundity. Indeed, perceived age predicts survival [1] and associates with molecular markers of aging such as telomere length [2]. Understanding the underlying molecular biology of perceived age is vital for identifying new aging therapies among other purposes, but studies are lacking thus far. As a first attempt, we performed genome-wide association studies (GWASs) of perceived facial age and wrinkling estimated from digital facial images by analyzing over eight million SNPs in 2,693 elderly Dutch Europeans from the Rotterdam Study. The strongest genetic associations with perceived facial age were found for multiple SNPs in the MC1R gene (p < 1 × 10(-7)). This effect was enhanced for a compound heterozygosity marker constructed from four pre-selected functional MC1R SNPs (p = 2.69 × 10(-12)), which was replicated in 599 Dutch Europeans from the Leiden Longevity Study (p = 0.042) and in 1,173 Europeans of the TwinsUK Study (p = 3 × 10(-3)). Individuals carrying the homozygote MC1R risk haplotype looked on average up to 2 years older than non-carriers. This association was independent of age, sex, skin color, and sun damage (wrinkling, pigmented spots) and persisted through different sun-exposure levels. Hence, a role for MC1R in youthful looks independent of its known melanin synthesis function is suggested. Our study uncovers the first genetic evidence explaining why some people look older for their age and provides new leads for further investigating the biological basis of how old or young people look.

Genetic studies have thus far identified a limited number of loci associated with human longevity by applying age at death or survival up to advanced ages as phenotype. As an alternative approach, one could first try to identify biomarkers of healthy ageing and the genetic variants associated with these traits and subsequently determine the association of these variants with human longevity. In the present study, we used this approach by testing whether the 35 baseline serum parameters measured in the Leiden Longevity Study (LLS) meet the proposed criteria for a biomarker of healthy ageing. The LLS consists of 421 families with long-lived siblings of European descent, who were recruited together with their offspring and the spouses of the offspring (controls). To test the four criteria for a biomarker of healthy ageing in the LLS, we determined the association of the serum parameters with chronological age, familial longevity, general practitioner-reported general health, and mortality. Out of the 35 serum parameters, we identified glucose, insulin, and triglycerides as biomarkers of healthy ageing, meeting all four criteria in the LLS. We subsequently showed that the genetic variants previously associated with these parameters are significantly enriched in the largest genome-wide association study for human longevity. In conclusion, we showed that biomarkers of healthy ageing can be used to leverage genetic studies into human longevity. We identified several genetic variants influencing the variation in glucose, insulin and triglycerides that contribute to human longevity.

Genome-wide association studies have identified numerous loci linked with complex diseases, for which the molecular mechanisms remain largely unclear. Comprehensive molecular profiling of circulating metabolites captures highly heritable traits, which can help to uncover metabolic pathophysiology underlying established disease variants. We conduct an extended genome-wide association study of genetic influences on 123 circulating metabolic traits quantified by nuclear magnetic resonance metabolomics from up to 24,925 individuals and identify eight novel loci for amino acids, pyruvate and fatty acids. The LPA locus link with cardiovascular risk exemplifies how detailed metabolic profiling may inform underlying aetiology via extensive associations with very-low-density lipoprotein and triglyceride metabolism. Genetic fine mapping and Mendelian randomization uncover wide-spread causal effects of lipoprotein(a) on overall lipoprotein metabolism and we assess potential pleiotropic consequences of genetically elevated lipoprotein(a) on diverse morbidities via electronic health-care records. Our findings strengthen the argument for safe LPA-targeted intervention to reduce cardiovascular risk.

The methylome is subject to genetic and environmental effects. Their impact may depend on sex and age, resulting in sex- and age-related physiological variation and disease susceptibility. Here we estimate the total heritability of DNA methylation levels in whole blood and estimate the variance explained by common single nucleotide polymorphisms at 411,169 sites in 2,603 individuals from twin families, to establish a catalogue of between-individual variation in DNA methylation. Heritability estimates vary across the genome (mean=19%) and interaction analyses reveal thousands of sites with sex-specific heritability as well as sites where the environmental variance increases with age. Integration with previously published data illustrates the impact of genome and environment across the lifespan at methylation sites associated with metabolic traits, smoking and ageing. These findings demonstrate that our catalogue holds valuable information on locations in the genome where methylation variation between people may reflect disease-relevant environmental exposures or genetic variation.

To identify intrinsic differences in cartilage gene expression profiles between wild-type- and Dio2-/--mice, as a mechanism to investigate factors that contribute to prolonged healthy tissue homeostasis.

Aging is associated with alterations in the content and patterns of DNA methylation virtually throughout the entire human lifespan. Reasons for these variations are not well understood. However, several lines of evidence suggest that the epigenetic instability in aging may be traced back to the alteration of the expression of DNA methyltransferases. Here, the association of the expression of DNA methyltransferases DNMT1 and DNMT3B with age has been analysed in the context of the MARK-AGE study, a large-scale cross-sectional study of the European general population. Using peripheral blood mononuclear cells, we assessed the variation of DNMT1 and DNMT3B gene expression in more than two thousand age-stratified women and men (35-75 years) recruited across eight European countries. Significant age-related changes were detected for both transcripts. The level of DNMT1 gradually dropped with aging but this was only observed up to the age of 64 years. By contrast, the expression of DNMT3B decreased linearly with increasing age and this association was particularly evident in females. We next attempted to trace the age-related changes of both transcripts to the influence of different variables that have an impact on changes of their expression in the population, including demographics, dietary and health habits, and clinical parameters. Our results indicate that age affects the expression of DNMT1 and DNMT3B as an almost independent variable in respect of all other variables evaluated.

Genome-wide association studies in Japanese populations recently identified common variants in the KCNQ1 gene to be associated with type 2 diabetes. We examined the association of these variants within KCNQ1 with type 2 diabetes in a Dutch population, investigated their effects on insulin secretion and metabolic traits and on the risk of developing complications in type 2 diabetes patients.

The PHArmacogenetic study of Statins in the Elderly at risk (PHASE) is a genome wide association study in the PROspective Study of Pravastatin in the Elderly at risk for vascular disease (PROSPER) that investigates the genetic variation responsible for the individual variation in drug response to pravastatin. Statins lower LDL-cholesterol in general by 30%, however not in all subjects. Moreover, clinical response is highly variable and adverse effects occur in a minority of patients. In this report we first describe the rationale of the PROSPER/PHASE project and second show that the PROSPER/PHASE study can be used to study pharmacogenetics in the elderly.

Lumbar disc degeneration (LDD) is an important cause of low back pain, which is a common and costly problem. LDD is characterised by disc space narrowing and osteophyte growth at the circumference of the disc. To date, the agnostic search of the genome by genome-wide association (GWA) to identify common variants associated with LDD has not been fruitful. This study is the first GWA meta-analysis of LDD.

In genome-wide association studies (GWAS) of complex traits, single SNP analysis is still the most applied approach. However, the identified SNPs have small effects and provide limited biological insight. A more appropriate approach to interpret GWAS data of complex traits is to analyze the combined effect of a SNP set grouped per pathway or gene region. We used this approach to study the joint effect on human longevity of genetic variation in two candidate pathways, the insulin/insulin-like growth factor (IGF-1) signaling (IIS) pathway and the telomere maintenance (TM) pathway. For the analyses, we used genotyped GWAS data of 403 unrelated nonagenarians from long-lived sibships collected in the Leiden Longevity Study and 1,670 younger population controls. We analyzed 1,021 SNPs in 68 IIS pathway genes and 88 SNPs in 13 TM pathway genes using four self-contained pathway tests (PLINK set-based test, Global test, GRASS and SNP ratio test). Although we observed small differences between the results of the different pathway tests, they showed consistent significant association of the IIS and TM pathway SNP sets with longevity. Analysis of gene SNP sets from these pathways indicates that the association of the IIS pathway is scattered over several genes (AKT1, AKT3, FOXO4, IGF2, INS, PIK3CA, SGK, SGK2, and YWHAG), while the association of the TM pathway seems to be mainly determined by one gene (POT1). In conclusion, this study shows that genetic variation in genes involved in the IIS and TM pathways is associated with human longevity.

Postzygotic mutations are DNA changes acquired from the zygote stage onwards throughout the lifespan. These changes lead to differences in DNA sequence among cells of an individual, potentially contributing to the etiology of complex disorders. Here we compared whole genome DNA sequence data of two monozygotic twin pairs, 40 and 100 years old, to detect somatic mosaicism. DNA samples were sequenced twice on two Illumina platforms (13X and 40X read depth) for increased specificity. Using differences in allelic ratios resulted in sets of 1,720 and 1,739 putative postzygotic mutations in the 40-year-old twin pair and 100-year-old twin pair, respectively, for subsequent enrichment analysis. This set of putative mutations was strongly (p < 4.37e-91) enriched in both twin pairs for regulatory elements. The corresponding genes were significantly enriched for genes that are alternatively spliced, and for genes involved in GTPase activity. This research shows that somatic mosaicism can be detected in monozygotic twin pairs by using allelic ratios calculated from DNA sequence data and that the mutations which are found by this approach are not randomly distributed throughout the genome.

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene-gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10-10), among which transcription factors were overrepresented (Fisher's P = 3.3 × 10-7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser.

Up to 33% and 25% of patients with end-stage hip and knee osteoarthritis (OA) are considered frail by the Groningen Frailty Indicator (GFI). This study aims to assess whether frail patients have lower functional gains after arthroplasty and to assess GFI as a tool to discriminate between good and adverse change score.