👤 G Willemsen

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success. Show less

CLN3 disease is a disorder of lysosomal homeostasis predominantly affecting the retina and the brain. The severity of the underlying mutations in We present three patients with an identical Contrasting our initial expectations, patients with a start codon variant in Show less

The nuclear receptor subfamily 1 group H member 4 (NR1H4 or farnesoid X receptor [FXR]) regulates bile acid synthesis, transport, and catabolism. FXR also regulates postprandial lipid and glucose metabolism. We performed quantitative proteomic analyses of liver tissues from mice to evaluate these functions and investigate whether FXR regulates amino acid metabolism. To study the role of FXR in mouse liver, we used mice with a disruption of Nr1h4 (FXR-knockout mice) and compared them with floxed control mice. Mice were gavaged with the FXR agonist obeticholic acid or vehicle for 11 days. Proteome analyses, as well as targeted metabolomics and chromatin immunoprecipitation, were performed on the livers of these mice. Primary rat hepatocytes were used to validate the role of FXR in amino acid catabolism by gene expression and metabolomics studies. Finally, control mice and mice with liver-specific disruption of Nr1h4 (liver FXR-knockout mice) were re-fed with a high-protein diet after 6 hours fasting and gavaged a In livers of control mice and primary rat hepatocytes, activation of FXR with obeticholic acid increased expression of proteins that regulate amino acid degradation, ureagenesis, and glutamine synthesis. We found FXR to bind to regulatory sites of genes encoding these proteins in control livers. Liver tissues from FXR-knockout mice had reduced expression of urea cycle proteins, and accumulated precursors of ureagenesis, compared with control mice. In liver FXR-knockout mice on a high-protein diet, the plasma concentration of newly formed urea was significantly decreased compared with controls. In addition, liver FXR-knockout mice had reduced hepatic expression of enzymes that regulate ammonium detoxification compared with controls. In contrast, obeticholic acid increased expression of genes encoding enzymes involved in ureagenesis compared with vehicle in C57Bl/6 mice. In livers of mice, FXR regulates amino acid catabolism and detoxification of ammonium via ureagenesis and glutamine synthesis. Failure of the urea cycle and hyperammonemia are common in patients with acute and chronic liver diseases; compounds that activate FXR might promote ammonium clearance in these patients. Show less

The Koolen-de Vries syndrome (KdVS; OMIM #610443), also known as the 17q21.31 microdeletion syndrome, is a clinically heterogeneous disorder characterised by (neonatal) hypotonia, developmental delay, moderate intellectual disability, and characteristic facial dysmorphism. Expressive language development is particularly impaired compared with receptive language or motor skills. Other frequently reported features include social and friendly behaviour, epilepsy, musculoskeletal anomalies, congenital heart defects, urogenital malformations, and ectodermal anomalies. The syndrome is caused by a truncating variant in the KAT8 regulatory NSL complex unit 1 (KANSL1) gene or by a 17q21.31 microdeletion encompassing KANSL1. Herein we describe a novel cohort of 45 individuals with KdVS of whom 33 have a 17q21.31 microdeletion and 12 a single-nucleotide variant (SNV) in KANSL1 (19 males, 26 females; age range 7 months to 50 years). We provide guidance about the potential pitfalls in the laboratory testing and emphasise the challenges of KANSL1 variant calling and DNA copy number analysis in the complex 17q21.31 region. Moreover, we present detailed phenotypic information, including neuropsychological features, that contribute to the broad phenotypic spectrum of the syndrome. Comparison of the phenotype of both the microdeletion and SNV patients does not show differences of clinical importance, stressing that haploinsufficiency of KANSL1 is sufficient to cause the full KdVS phenotype. Show less

Metformin is used as a first-line oral treatment for type 2 diabetes (T2D). However, the underlying mechanism is not fully understood. Here, we aimed to comprehensively investigate the pleiotropic effects of metformin. We analyzed both metabolomic and genomic data of the population-based KORA cohort. To evaluate the effect of metformin treatment on metabolite concentrations, we quantified 131 metabolites in fasting serum samples and used multivariable linear regression models in three independent cross-sectional studies (n = 151 patients with T2D treated with metformin [mt-T2D]). Additionally, we used linear mixed-effect models to study the longitudinal KORA samples (n = 912) and performed mediation analyses to investigate the effects of metformin intake on blood lipid profiles. We combined genotyping data with the identified metformin-associated metabolites in KORA individuals (n = 1,809) and explored the underlying pathways. We found significantly lower (P < 5.0E-06) concentrations of three metabolites (acyl-alkyl phosphatidylcholines [PCs]) when comparing mt-T2D with four control groups who were not using glucose-lowering oral medication. These findings were controlled for conventional risk factors of T2D and replicated in two independent studies. Furthermore, we observed that the levels of these metabolites decreased significantly in patients after they started metformin treatment during 7 years' follow-up. The reduction of these metabolites was also associated with a lowered blood level of LDL cholesterol (LDL-C). Variations of these three metabolites were significantly associated with 17 genes (including FADS1 and FADS2) and controlled by AMPK, a metformin target. Our results indicate that metformin intake activates AMPK and consequently suppresses FADS, which leads to reduced levels of the three acyl-alkyl PCs and LDL-C. Our findings suggest potential beneficial effects of metformin in the prevention of cardiovascular disease. Show less

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10); and (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD. Show less

Concentrations of liver enzymes in plasma are widely used as indicators of liver disease. We carried out a genome-wide association study in 61,089 individuals, identifying 42 loci associated with concentrations of liver enzymes in plasma, of which 32 are new associations (P = 10(-8) to P = 10(-190)). We used functional genomic approaches including metabonomic profiling and gene expression analyses to identify probable candidate genes at these regions. We identified 69 candidate genes, including genes involved in biliary transport (ATP8B1 and ABCB11), glucose, carbohydrate and lipid metabolism (FADS1, FADS2, GCKR, JMJD1C, HNF1A, MLXIPL, PNPLA3, PPP1R3B, SLC2A2 and TRIB1), glycoprotein biosynthesis and cell surface glycobiology (ABO, ASGR1, FUT2, GPLD1 and ST3GAL4), inflammation and immunity (CD276, CDH6, GCKR, HNF1A, HPR, ITGA1, RORA and STAT4) and glutathione metabolism (GSTT1, GSTT2 and GGT), as well as several genes of uncertain or unknown function (including ABHD12, EFHD1, EFNA1, EPHA2, MICAL3 and ZNF827). Our results provide new insight into genetic mechanisms and pathways influencing markers of liver function. Show less

C-reactive protein (CRP) is a heritable marker of chronic inflammation that is strongly associated with cardiovascular disease. We sought to identify genetic variants that are associated with CRP levels. We performed a genome-wide association analysis of CRP in 66 185 participants from 15 population-based studies. We sought replication for the genome-wide significant and suggestive loci in a replication panel comprising 16 540 individuals from 10 independent studies. We found 18 genome-wide significant loci, and we provided evidence of replication for 8 of them. Our results confirm 7 previously known loci and introduce 11 novel loci that are implicated in pathways related to the metabolic syndrome (APOC1, HNF1A, LEPR, GCKR, HNF4A, and PTPN2) or the immune system (CRP, IL6R, NLRP3, IL1F10, and IRF1) or that reside in regions previously not known to play a role in chronic inflammation (PPP1R3B, SALL1, PABPC4, ASCL1, RORA, and BCL7B). We found a significant interaction of body mass index with LEPR (P<2.9×10(-6)). A weighted genetic risk score that was developed to summarize the effect of risk alleles was strongly associated with CRP levels and explained ≈5% of the trait variance; however, there was no evidence for these genetic variants explaining the association of CRP with coronary heart disease. We identified 18 loci that were associated with CRP levels. Our study highlights immune response and metabolic regulatory pathways involved in the regulation of chronic inflammation. Show less

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation. Show less

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes. Show less

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing. Show less

Anomalies in neuropeptides and neuroactive amino acids have been postulated to play a role in neurodegeneration in a variety of diseases including the inherited neuronal ceroid lipofuscinoses (NCLs, Batten disease). These are often indicated by concentration changes in cerebrospinal fluid (CSF). Here we compare CSF neuropeptide concentrations in patients with the classical juvenile CLN3 form of NCL and the classical late infantile CLN2 form with neuropeptide and neuroactive amino acid concentrations in CSF from sheep with the late infantile variant CLN6 form. A marked disease related increase in CSF concentrations of neuron specific enolase and tau protein was noted in the juvenile CLN3 patients but this was not observed in an advanced CLN2 patient nor CLN6 affected sheep. No changes were noted in S-100b, GFAP or MBP in patients or of S-100b, GFAP or IGF-1 in affected sheep. There were no disease related changes in CSF concentrations of the neuroactive amino acids, aspartate, glutamate, serine, glutamine, glycine, taurine and GABA in these sheep. The changes observed in the CLN3 patients may be progressive markers of neurodegeneration, or of underlying metabolic changes perhaps associated with CLN3 specific changes in neuroactive amino acids, as have been postulated. The lack of changes in the CLN2 and CLN6 subjects indicate that these changes are not shared by the CLN2 or CLN6 forms and changes in CSF concentrations of these compounds are unreliable as biomarkers of neurodegeneration in the NCLs in general. Show less

Recent genome-wide association (GWA) studies of lipids have been conducted in samples ascertained for other phenotypes, particularly diabetes. Here we report the first GWA analysis of loci affecting total cholesterol (TC), low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglycerides sampled randomly from 16 population-based cohorts and genotyped using mainly the Illumina HumanHap300-Duo platform. Our study included a total of 17,797-22,562 persons, aged 18-104 years and from geographic regions spanning from the Nordic countries to Southern Europe. We established 22 loci associated with serum lipid levels at a genome-wide significance level (P < 5 x 10(-8)), including 16 loci that were identified by previous GWA studies. The six newly identified loci in our cohort samples are ABCG5 (TC, P = 1.5 x 10(-11); LDL, P = 2.6 x 10(-10)), TMEM57 (TC, P = 5.4 x 10(-10)), CTCF-PRMT8 region (HDL, P = 8.3 x 10(-16)), DNAH11 (LDL, P = 6.1 x 10(-9)), FADS3-FADS2 (TC, P = 1.5 x 10(-10); LDL, P = 4.4 x 10(-13)) and MADD-FOLH1 region (HDL, P = 6 x 10(-11)). For three loci, effect sizes differed significantly by sex. Genetic risk scores based on lipid loci explain up to 4.8% of variation in lipids and were also associated with increased intima media thickness (P = 0.001) and coronary heart disease incidence (P = 0.04). The genetic risk score improves the screening of high-risk groups of dyslipidemia over classical risk factors. Show less