👤 Bjarke Feenstra

Acute myeloid leukemia (AML) is a complex hematologic malignancy with multiple disease subgroups defined by somatic mutations and heterogeneous outcomes. Although genome-wide association studies (GWAS) have identified a small number of common genetic variants influencing AML risk, the heritable component of this disease outside of familial susceptibility remains largely undefined. Here, we perform a meta-analysis of 4 published GWAS plus 2 new GWAS, totaling 4710 AML cases and 12 938 controls. We identify a new genome-wide significant risk locus for pan-AML at 2p23.3 (rs4665765; P = 1.35 × 10-8; EFR3B, POMC, DNMT3A, and DNAJC27), which also significantly associates with patient survival (P = 6.09 × 10-3). Our analysis also identifies 3 new genome-wide significant risk loci for disease subgroups, including AML with deletions of chromosome 5 and/or 7 at 1q23.3 (rs12078864; P = 7.0 × 10-10; DUSP23) and cytogenetically complex AML at 2q33.3 (rs12988876; P = 3.28 × 10-8; PARD3B) and 2p21 (rs79918355; P = 1.60 × 10-9; EPCAM). We also investigated loci previously associated with the risk of clonal hematopoiesis (CH) or CH of indeterminate potential and identified several variants associated with the risk of AML. Our results further inform on AML etiology and demonstrate the existence of disease subgroup specific risk loci. Show less

A large number of genetic loci are associated with adult body mass index. However, the genetics of childhood body mass index are largely unknown. We performed a meta-analysis of genome-wide association studies of childhood body mass index, using sex- and age-adjusted standard deviation scores. We included 35 668 children from 20 studies in the discovery phase and 11 873 children from 13 studies in the replication phase. In total, 15 loci reached genome-wide significance (P-value < 5 × 10(-8)) in the joint discovery and replication analysis, of which 12 are previously identified loci in or close to ADCY3, GNPDA2, TMEM18, SEC16B, FAIM2, FTO, TFAP2B, TNNI3K, MC4R, GPR61, LMX1B and OLFM4 associated with adult body mass index or childhood obesity. We identified three novel loci: rs13253111 near ELP3, rs8092503 near RAB27B and rs13387838 near ADAM23. Per additional risk allele, body mass index increased 0.04 Standard Deviation Score (SDS) [Standard Error (SE) 0.007], 0.05 SDS (SE 0.008) and 0.14 SDS (SE 0.025), for rs13253111, rs8092503 and rs13387838, respectively. A genetic risk score combining all 15 SNPs showed that each additional average risk allele was associated with a 0.073 SDS (SE 0.011, P-value = 3.12 × 10(-10)) increase in childhood body mass index in a population of 1955 children. This risk score explained 2% of the variance in childhood body mass index. This study highlights the shared genetic background between childhood and adult body mass index and adds three novel loci. These loci likely represent age-related differences in strength of the associations with body mass index. Show less

The Wnt/β-catenin signaling pathway is important for multiple developmental processes and tissue maintenance in adults. Consequently, deregulated signaling is involved in a range of human diseases including cancer and developmental defects. A better understanding of the intricate regulatory mechanism and effect of physiological (active) and pathophysiological (hyperactive) WNT signaling is important for predicting treatment response and developing novel therapies. The constitutively expressed CTNNB1 (commonly and hereafter referred to as β-catenin) is degraded by a destruction complex, composed of amongst others AXIN1 and GSK3. The destruction complex is inhibited during active WNT signaling, leading to β-catenin stabilization and induction of β-catenin/TCF target genes. In this study we investigated the mechanism and effect of β-catenin stabilization during active and hyperactive WNT signaling in a combined in silico and in vitro approach. We constructed a Petri net model of Wnt/β-catenin signaling including main players from the plasma membrane (WNT ligands and receptors), cytoplasmic effectors and the downstream negative feedback target gene AXIN2. We validated that our model can be used to simulate both active (WNT stimulation) and hyperactive (GSK3 inhibition) signaling by comparing our simulation and experimental data. We used this experimentally validated model to get further insights into the effect of the negative feedback regulator AXIN2 upon WNT stimulation and observed an attenuated β-catenin stabilization. We furthermore simulated the effect of APC inactivating mutations, yielding a stabilization of β-catenin levels comparable to the Wnt-pathway activities observed in colorectal and breast cancer. Our model can be used for further investigation and viable predictions of the role of Wnt/β-catenin signaling in oncogenesis and development. Show less

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee. Show less

Dental caries is the most common chronic disease in children and a major public health concern due to its increasing incidence, serious health and social co-morbidities, and socio-demographic disparities in disease burden. We performed the first genome-wide association scan for dental caries to identify associated genetic loci and nominate candidate genes affecting tooth decay in 1305 US children ages 3-12 yrs. Affection status was defined as 1 or more primary teeth with evidence of decay based on intra-oral examination. No associations met strict criteria for genome-wide significance (p < 10E-7); however, several loci (ACTN2, MTR, and EDARADD, MPPED2, and LPO) with plausible biological roles in dental caries exhibited suggestive evidence for association. Analyses stratified by home fluoride level yielded additional suggestive loci, including TFIP11 in the low-fluoride group, and EPHA7 and ZMPSTE24 in the sufficient-fluoride group. Suggestive loci were tested but not significantly replicated in an independent sample (N = 1695, ages 2-7 yrs) after adjustment for multiple comparisons. This study reinforces the complexity of dental caries, suggesting that numerous loci, mostly having small effects, are involved in cariogenesis. Verification/replication of suggestive loci may highlight biological mechanisms and/or pathways leading to a fuller understanding of the genetic risks for dental caries. 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

Heterochromatin is important for gene regulation and chromosome structure, but the genes that are occupied by heterochromatin proteins in the mammalian genome are largely unknown. We have adapted the DamID method to systematically identify target genes of the heterochromatin proteins HP1 and SUV39H1 in human and mouse cells. Unexpectedly, we found that CBX1 (formerly HP1beta) and SUV39H1 bind to genes encoding KRAB domain containing zinc finger (KRAB-ZNF) transcriptional repressors. These genes constitute one of the largest gene families and are organized in clusters in the human genome. Preference of CBX1 for this gene family was observed in both human and mouse cells. High-resolution mapping on human chromosome 19 revealed that CBX1 coats large domains 0.1-4 Mb in size, which coincide with the position of KRAB-ZNF gene clusters. These domains show an intricate CBX1 binding pattern: While CBX1 is globally elevated throughout the domains, it is absent from the promoters and binds more strongly to the 3' ends of KRAB-ZNF genes. KRAB-ZNF domains contain large numbers of LINE elements, which may contribute to CBX1 recruitment. These results uncover a surprising link between heterochromatin and a large family of regulatory genes in mammals. We suggest a role for heterochromatin in the evolution of the KRAB-ZNF gene family. Show less