To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of... Show moreTo date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases. Show less
Age-related macular degeneration (AMD) has a strong genetic basis, but environmental factors such as smoking and a healthy diet can decrease the genetic fate by up to 50%. Current guidelines for... Show moreAge-related macular degeneration (AMD) has a strong genetic basis, but environmental factors such as smoking and a healthy diet can decrease the genetic fate by up to 50%. Current guidelines for clinical management include recommendations for a healthy lifestyle and antioxidant supplementation. However, many ophthalmologists do not inform their patients of this AMD-beneficial lifestyle. An important reason is the lack of trust that transition of lifestyle will be feasible in persons of advanced age and lack of methodology to measure lifestyle or its biological effects. To address these issues, we set up the lifestyle intervention study AMD-Life. It aims to investigate whether personalized risk-profiling (including genetic testing) and/or additional coaching can motivate patients to change their lifestyle. It also explores which biomarkers best reflect lifestyle change beneficial for AMD. The first year is a three-arm, self-contained open-label randomized clinical trial. A total of 150 AMD patients aged 55-85 years were randomized into three arms: (A) merely standard recommendations; (B) A conditions plus personalized risk profiling based on genetics and lifestyle, (C) B conditions plus coaching. The second year tests sustainability of lifestyle changes without active intervention. AMD-Life can provide further insight into the relevance of these interventions for the clinical management of AMD. Show less
Spek, A. van der; Karamujic-Comic, H.; Pool, R.; Bot, M.; Beekman, M.; Garmaeva, S.; ... ; BBMRI Metabolomics Consortium 2022
Telomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer and mortality. Lipid and fatty acid metabolism... Show moreTelomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer and mortality. Lipid and fatty acid metabolism have been associated with telomere shortening. We have conducted an in-depth study investigating the association of metabolic biomarkers with telomere length (LTL). We performed an association analysis of 226 metabolic biomarkers with LTL using data from 11 775 individuals from six independent population-based cohorts (BBMRI-NL consortium). Metabolic biomarkers include lipoprotein lipids and subclasses, fatty acids, amino acids, glycolysis measures and ketone bodies. LTL was measured by quantitative polymerase chain reaction or F1owFISH. Linear regression analysis was performed adjusting for age, sex, lipid-lowering medication and cohort-specific covariates (model 1) and additionally for body mass index (BMI) and smoking (model 2), followed by inverse variance-weighted meta-analyses (significance threshold P me t a = 6.5 x 10(-4)). We identified four metabolic biomarkers positively associated with LTL, including two cholesterol to lipid ratios in small VLDL (S-VLDL-C % and S-VLDL-CE %) and two omega-6 fatty acid ratios (FAw6/FA and LA/FA). After additionally adjusting for BMI and smoking, these metabolic biomarkers remained associated with LTL with similar effect estimates. In addition, cholesterol esters in very small VLDL (XS-VLDL-CE) became significantly associated with LTL (P = 3.6 x 10(-4)). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (P = 1.9 x 10(-4)). To conclude, we identified multiple metabolic biomarkers involved in lipid and fatty acid metabolism that may be involved in LTL biology. Longitudinal studies are needed to exclude reversed causation. Show less
Maas, S.C.E.; Vidaki, A.; Teumer, A.; Costeira, R.; Wilson, R.; Dongen, J. van; ... ; Kayser, M. 2021
Background Information on long-term alcohol consumption is relevant for medical and public health research, disease therapy, and other areas. Recently, DNA methylation-based inference of alcohol... Show moreBackground Information on long-term alcohol consumption is relevant for medical and public health research, disease therapy, and other areas. Recently, DNA methylation-based inference of alcohol consumption from blood was reported with high accuracy, but these results were based on employing the same dataset for model training and testing, which can lead to accuracy overestimation. Moreover, only subsets of alcohol consumption categories were used, which makes it impossible to extrapolate such models to the general population. By using data from eight population-based European cohorts (N = 4677), we internally and externally validated the previously reported biomarkers and models for epigenetic inference of alcohol consumption from blood and developed new models comprising all data from all categories. Results By employing data from six European cohorts (N = 2883), we empirically tested the reproducibility of the previously suggested biomarkers and prediction models via ten-fold internal cross-validation. In contrast to previous findings, all seven models based on 144-CpGs yielded lower mean AUCs compared to the models with less CpGs. For instance, the 144-CpG heavy versus non-drinkers model gave an AUC of 0.78 +/- 0.06, while the 5 and 23 CpG models achieved 0.83 +/- 0.05, respectively. The transportability of the models was empirically tested via external validation in three independent European cohorts (N = 1794), revealing high AUC variance between datasets within models. For instance, the 144-CpG heavy versus non-drinkers model yielded AUCs ranging from 0.60 to 0.84 between datasets. The newly developed models that considered data from all categories showed low AUCs but gave low AUC variation in the external validation. For instance, the 144-CpG heavy and at-risk versus light and non-drinkers model achieved AUCs of 0.67 +/- 0.02 in the internal cross-validation and 0.61-0.66 in the external validation datasets. Conclusions The outcomes of our internal and external validation demonstrate that the previously reported prediction models suffer from both overfitting and accuracy overestimation. Our results show that the previously proposed biomarkers are not yet sufficient for accurate and robust inference of alcohol consumption from blood. Overall, our findings imply that DNA methylation prediction biomarkers and models need to be improved considerably before epigenetic inference of alcohol consumption from blood can be considered for practical applications. Show less
Aims The HERMES (HEart failure Molecular Epidemiology for Therapeutic targets) consortium aims to identify the genomic and molecular basis of heart failure.Methods and results The consortium... Show moreAims The HERMES (HEart failure Molecular Epidemiology for Therapeutic targets) consortium aims to identify the genomic and molecular basis of heart failure.Methods and results The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome-wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow-up following heart failure diagnosis ranged from 2 to 116 months. Forty-nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34-90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of >1.10 for common variants (allele frequency > 0.05) and >1.20 for low-frequency variants (allele frequency 0.01-0.05) at P < 5 x 10(-8) under an additive genetic model.Conclusions HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction. Show less
Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100... Show moreOsteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation. Show less
Introduction: Increasing resistance to beta-lactam antibiotics is an alarming development worldwide. Fecal carriership of TEM, SHV, CTX-M and CMY was studied in a community-dwelling population of... Show moreIntroduction: Increasing resistance to beta-lactam antibiotics is an alarming development worldwide. Fecal carriership of TEM, SHV, CTX-M and CMY was studied in a community-dwelling population of middle-aged and elderly individuals.Patients and methods: Feces was obtained from individuals of the Rotterdam Study. Carriership of the TEM, SHV, CTX-M and CMY genes was determined using real-time polymerase chain reaction (qPCR). Possible associations were investigated between carriership of these genes and several risk factors, such as the use of antimicrobial drugs, diabetes mellitus, protein pump inhibitor (PPI) use, travelling, the composition of the gut microbiota, and intake of certain foods.Results: The most prevalent gene was TEM (53.0%), followed by SHV (18.4%), CTX-M (5.4%) and CMY (3.6%). Use of penicillins with extended spectrum was associated with TEM carriership, whereas use of macrolides and lincosamides was associated with TEM and SHV carriership. Interestingly, use of PPIs was associated with a higher prevalence of carriership of TEM, SHV and CMY (TEM: odds ratio [OR] 1.34; 95% confidence interval [CI] 1.05-1.77; SHV: OR 2.17; 95%CI 1.55-2.87; CMY: OR 2.26; 95%CI 1.23-4.11). Furthermore, associations were found between the richness and composition of the gut microbiota and TEM and SHV carriership.Conclusions: The prevalence of carriership of TEM was substantial, but the prevalence of carriership of the extended-spectrum beta-lactamase gene, CTX-M and the AmpC beta-lactamase gene, CMY was relatively low in this community-dwelling, population-based cohort. The composition of the microbiota might play a role in the retention of resistance genes, but future studies are necessary to further elucidate this relationship. (C) 2021 Published by Elsevier Ltd. Show less
DNA methylation quantitative trait locus (mQTL) analyses on 32,851 participants identify genetic variants associated with DNA methylation at 420,509 sites in blood, resulting in a database of >... Show moreDNA methylation quantitative trait locus (mQTL) analyses on 32,851 participants identify genetic variants associated with DNA methylation at 420,509 sites in blood, resulting in a database of >270,000 independent mQTLs.Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated. Show less
ObjectivesVitamin K is hypothesised to play a role in osteoarthritis (OA) pathogenesis through effects on vitamin K-dependent bone and cartilage proteins, and therefore may represent a modifiable... Show moreObjectivesVitamin K is hypothesised to play a role in osteoarthritis (OA) pathogenesis through effects on vitamin K-dependent bone and cartilage proteins, and therefore may represent a modifiable risk factor. A genetic variant in a vitamin K-dependent protein that is an essential inhibitor for cartilage calcification, matrix Gla protein (MGP), was associated with an increased risk for OA. Vitamin K antagonist anticoagulants (VKAs), such as warfarin and acenocoumarol, act as anticoagulants through inhibition of vitamin K-dependent blood coagulation proteins. VKAs likely also affect the functioning of other vitamin K-dependent proteins such as MGP.MethodsWe investigated the effect of acenocoumarol usage on progression and incidence of radiographic OA in 3494 participants of the Rotterdam Study cohort. We also examined the effect of MGP and VKORC1 single nucleotide variants on this association.ResultsAcenocoumarol usage was associated with an increased risk of OA incidence and progression (OR=2.50, 95% CI=1.94-3.20), both for knee (OR=2.34, 95% CI=1.67-3.22) and hip OA (OR=2.74, 95% CI=1.82-4.11). Among acenocoumarol users, carriers of the high VKORC1(BB) expression haplotype together with the MGP OA risk allele (rs1800801-T) had an increased risk of OA incidence and progression (OR=4.18, 95% CI=2.69-6.50), while this relationship was not present in non-users of that group (OR=1.01, 95% CI=0.78-1.33).ConclusionsThese findings support the importance of vitamin K and vitamin K-dependent proteins, as MGP, in the pathogenesis of OA. Additionally, these results may have direct implications for the clinical prevention of OA, supporting the consideration of direct oral anticoagulants in favour of VKAs. Show less
Background and Aims Previous small studies have appraised the gut microbiome (GM) in steatosis, but large-scale studies are lacking. We studied the association of the GM diversity and composition,... Show moreBackground and Aims Previous small studies have appraised the gut microbiome (GM) in steatosis, but large-scale studies are lacking. We studied the association of the GM diversity and composition, plasma metabolites, predicted functional metagenomics, and steatosis.Approach and Results This is a cross-sectional analysis of the prospective population-based Rotterdam Study. We used 16S ribosomal RNA gene sequencing and determined taxonomy using the SILVA reference database. Alpha diversity and beta diversity were calculated using the Shannon diversity index and Bray-Curtis dissimilarities. Differences were tested across steatosis using permutational multivariate analysis of variance. Hepatic steatosis was diagnosed by ultrasonography. We subsequently selected genera using regularized regression. The functional metagenome was predicted based on the GM using Kyoto Encyclopedia of Genes and Genomes pathways. Serum metabolomics were assessed using high-throughput proton nuclear magnetic resonance. All analyses were adjusted for age, sex, body mass index, alcohol, diet, and proton-pump inhibitors. We included 1,355 participants, of whom 472 had steatosis. Alpha diversity was lower in steatosis (P = 1.1 center dot 10(-9)), and beta diversity varied across steatosis strata (P = 0.001). Lasso selected 37 genera of which three remained significantly associated after adjustment (Coprococcus3: beta = -65; Ruminococcus Gauvreauiigroup: beta = 62; and Ruminococcus Gnavusgroup: beta = 45, Q-value = 0.037). Predicted metagenome analyses revealed that pathways of secondary bile-acid synthesis and biotin metabolism were present, and D-alanine metabolism was absent in steatosis. Metabolic profiles showed positive associations for aromatic and branched chain amino acids and glycoprotein acetyls with steatosis and R. Gnavusgroup, whereas these metabolites were inversely associated with alpha diversity and Coprococcus3.Conclusions We confirmed, on a large-scale, the lower microbial diversity and association of Coprococcus and Ruminococcus Gnavus with steatosis. We additionally showed that steatosis and alpha diversity share opposite metabolic profiles. Show less
Background Despite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the... Show moreBackground Despite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the complexity of the different hand joints and heterogeneity of OA pathology.Methods We used data from the Rotterdam Study (RSI, RSII and RSIII) to create three hand OA phenotypes based on clustering patterns of radiographic OA severity to increase power in our modest discovery genome-wide association studies in the RS (n=8700), and sought replication in an independent cohort, the Framingham Heart Study (n=1203). We used multiple approaches that leverage different levels of information and functional data to further investigate the underlying biological mechanisms and candidate genes for replicated loci. We also attempted to replicate known OA loci at other joint sites, including the hips and knees.Results We found two novel genome-wide significant loci for OA in the thumb joints. We identified WNT9A as a possible novel causal gene involved in OA pathogenesis. Furthermore, several previously identified genetic loci for OA seem to confer risk for OA across multiple joints: TGFa, RUNX2, COL27A1, ASTN2, IL11 and GDF5 loci.Conclusions We identified a robust novel genetic locus for hand OA on chromosome 1, of which WNT9A is the most likely causal gene. In addition, multiple genetic loci were identified to be associated with OA across multiple joints. Our study confirms the potential for novel insight into the genetic architecture of OA by using biologically meaningful stratified phenotypes. Show less
Araghi, S.O.; Kiefte-de Jong, J.C.; Dijk, S.C. van; Swart, K.M.A.; Ploegmakers, K.J.; Zillikens, M.C.; ... ; Velde, N. van der 2021
Background & aims: In the initial B-proof, we found inconsistent results of B vitamin supplementation. However, the debate regarding the effects of B vitamins on age-related diseases continues.... Show moreBackground & aims: In the initial B-proof, we found inconsistent results of B vitamin supplementation. However, the debate regarding the effects of B vitamins on age-related diseases continues. Therefore, our aim was to investigate the long-term effects (5-7 years follow-up) of an intervention with folic acid and vitamin-B12 supplementation on fracture and cardiovascular disease risk.Methods: Extended follow-up of the B-PROOF trial, a multi-center, double-blind randomized placebo controlled trial designed to assess the effect of 2-3 years daily supplementation with folic acid (400 mg) and vitamin-B12 (500 mg) versus placebo (n = 2,919). Primary outcome was verified self reported fracture incidence and secondary outcomes were self-reported cardiovascular endpoints, which were collected through a follow-up questionnaires Proportional hazard analyses was used for the effect of the intervention on risk of fracture(s) and logistic regression for the effect of the intervention on risk of cardiovascular disease.Results: A total of 1,298 individuals (4 4.5%) participated in the second follow-up round with median of 54 months [51-58], (n = 662 and n = 636, treatment versus placebo group). Median age at baseline was 71.0 years [68.0-76.0] for both groups. No effect was observed of the intervention on osteoporotic fracture or any fracture risk after a follow-up (HR: 0.99, 95% CI: 0.62-1.59 and HR: 0.77; 95% CI: 0.50-1.19, respectively), nor on cardiovascular or cerebrovascular disease risk (OR: 1.05; 95%CI: 0.80-1.44 and OR: 0.85; 95%CI: 0.50-1.45, respectively). Potential interaction by baseline homocysteine concentration was observed for osteoporotic-and any fracture (p = 0.10 and 0.06 respectively), which indicated a significantly lower risk of any fracture in the treatment group with higher total homocysteine concentrations (>15.1 mmol/l). No age-dependent effects were present.Conclusions: This study supports and extends previous null -findings of the B-PROOF trial and shows that supplementation of folic acid and vitamin-B12 has no effect on fracture risk, nor on cardiovascular disease in older individuals over a longer follow-up period. However, B-vitamin supplementation may be beneficial in reducing fractures in individuals with high total homocysteine concentrations, a finding which needs to be replicated. (c) 2020 Elsevier Ltd and European Society for Clinical Nutrition and Metabolism. All rights reserved. Show less
Mulder, M.; Radjabzadeh, D.; Kiefte-de Jong, J.C.; Uitterlinden, A.G.; Kraaij, R.; Stricker, B.H.; Verbon, A. 2020
Introduction Antimicrobial drugs are known to have effects on the human gut microbiota. We studied the long-term temporal relationship between several antimicrobial drug groups and the composition... Show moreIntroduction Antimicrobial drugs are known to have effects on the human gut microbiota. We studied the long-term temporal relationship between several antimicrobial drug groups and the composition of the human gut microbiota determined in feces samples.Methods Feces samples were obtained from a community-dwelling cohort of middle-aged and elderly individuals (Rotterdam Study). Bacterial DNA was isolated and sequenced using V3/V4 16 S ribosomal RNA sequencing (Illumina MiSeq). The time between the last prescription of several antimicrobial drug groups and the day of sampling was categorized into 0-12, 12-24, 24-48 and >48 months. The effects of the antimicrobial drug groups on the Shannon alpha-diversity (diversity), the Bray-Curtis beta-diversity (community structure), the Firmicutes/Bacteroidetes (F/B) ratio and individual genera were determined.Results We studied the gut microbiota of 1413 individuals (57.5% female, median age 62.6 years). The alpha-diversity was significantly lower up to 4 years after prescriptions of macrolides and lincosamides. It was also lower in the first year after the use of beta-lactams. The community structure (beta-diversity) of the microbiota was significantly different up to 4 years for macrolides and lincosamides, the first year for beta-lactams and at least the first year for quinolones. For the F/B ratio, drugs with a high anaerobic activity shifted the ratio toward Firmicutes in the first year whereas other antimicrobial drugs shifted the ratio toward Bacteroidetes.Conclusion Use of antimicrobial drugs is associated with a shift in the composition of the gut microbiota.These effects differ in strength and duration, depending on the antimicrobial drug group used. These findings should be considered when prescribing antimicrobial drugs. Show less
Meessen, J.M.T.A.; Saberi-Hosnijeh, F.; Bomer, N.; Hollander, W. den; Bom, J.G. van der; Hilten, J.A. van; ... ; Meulenbelt, I. 2020
Higher body mass index (BMI) is associated with osteoarthritis (OA) in both weight-bearing and non-weight-bearing joints, suggesting a link between OA and poor metabolic health beyond mechanical... Show moreHigher body mass index (BMI) is associated with osteoarthritis (OA) in both weight-bearing and non-weight-bearing joints, suggesting a link between OA and poor metabolic health beyond mechanical loading. This risk may be influenced by systemic factors accompanying BMI. Fluctuations in concentrations of metabolites may mark or even contribute to development of OA. This study explores the association of metabolites with radiographic knee/hip OA prevalence and progression. A H-1-NMR-metabolomics assay was performed on plasma samples of 1564 cases for prevalent OA and 2,125 controls collected from the Rotterdam Study, CHECK, GARP/NORREF and LUMC-arthroplasty cohorts. OA prevalence and 5 to 10 year progression was assessed by means of Kellgren-Lawrence (KL) score and the OARSI-atlas. End-stage knee/hip OA (TJA) was defined as indication for arthroplasty surgery. Controls did not have OA at baseline or follow-up. Principal component analysis of 227 metabolites demonstrated 23 factors, of which 19 remained interpretable after quality-control. Associations of factor scores with OA definitions were investigated with logistic regression. Fatty acids chain length (FALen), which was included in two factors which associated with TJA, was individually associated with both overall OA as well as TJA. Increased Fatty Acid chain Length is associated with OA. Show less
Lakeman, I.M.M.; Rodriguez-Girondo, M.; Lee, A.; Ruiter, R.; Stricker, B.H.; Wijnant, S.R.A.; ... ; Devilee, P. 2020
Purpose We evaluated the performance of the recently extended Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA version 5) in a Dutch prospective cohort,... Show morePurpose We evaluated the performance of the recently extended Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA version 5) in a Dutch prospective cohort, using a polygenic risk score (PRS) based on 313 breast cancer (BC)-associated variants (PRS313) and other, nongenetic risk factors. Methods Since 1989, 6522 women without BC aged 45 or older of European descent have been included in the Rotterdam Study. The PRS(313)was calculated per 1 SD in controls from the Breast Cancer Association Consortium (BCAC). Cox regression analysis was performed to estimate the association between the PRS(313)and incident BC risk. Cumulative 10-year risks were calculated with BOADICEA including different sets of variables (age, risk factors and PRS313). C-statistics were used to evaluate discriminative ability. Results In total, 320 women developed BC. The PRS(313)was significantly associated with BC (hazard ratio [HR] per SD of 1.56, 95% confidence interval [CI] [1.40-1.73]). Using 10-year risk estimates including age and the PRS313, other risk factors improved the discriminatory ability of the BOADICEA model marginally, from a C-statistic of 0.636 to 0.653. Conclusions The effect size of the PRS(313)is highly reproducible in the Dutch population. Our results validate the BOADICEA v5 model for BC risk assessment in the Dutch general population. Show less
Meddens, S.F.W.; Vlaming, R. de; Bowers, P.; Burik, C.A.P.; Linner, R.K.; Lee, C.; ... ; Lifelines Cohort Study 2020
We conducted genome-wide association studies (GWAS) of relative intake from the macronutrients fat, protein, carbohydrates, and sugar in over 235,000 individuals of European ancestries. We... Show moreWe conducted genome-wide association studies (GWAS) of relative intake from the macronutrients fat, protein, carbohydrates, and sugar in over 235,000 individuals of European ancestries. We identified 21 unique, approximately independent lead SNPs. Fourteen lead SNPs are uniquely associated with one macronutrient at genome-wide significance (P < 5 x 10(-8)), while five of the 21 lead SNPs reach suggestive significance (P < 1 x 10(-5)) for at least one other macronutrient. While the phenotypes are genetically correlated, each phenotype carries a partially unique genetic architecture. Relative protein intake exhibits the strongest relationships with poor health, including positive genetic associations with obesity, type 2 diabetes, and heart disease (r(g) approximate to 0.15-0.5). In contrast, relative carbohydrate and sugar intake have negative genetic correlations with waist circumference, waist-hip ratio, and neighborhood deprivation (|r(g)| approximate to 0.1-0.3) and positive genetic correlations with physical activity (r(g) approximate to 0.1 and 0.2). Relative fat intake has no consistent pattern of genetic correlations with poor health but has a negative genetic correlation with educational attainment (r(g) approximate to-0.1). Although our analyses do not allow us to draw causal conclusions, we find no evidence of negative health consequences associated with relative carbohydrate, sugar, or fat intake. However, our results are consistent with the hypothesis that relative protein intake plays a role in the etiology of metabolic dysfunction. Show less
Suratannon, N.; Wijck, R.T.A. van; Broer, L.; Xue, L.X.; Meurs, J.B.J. van; Barendregt, B.H.; ... ; South East Asia Primary Immunode 2020
Background: Genetic tests for primary immunodeficiency disorders (PIDs) are expensive, time-consuming, and not easily accessible in developing countries. Therefore, we studied the feasibility of a... Show moreBackground: Genetic tests for primary immunodeficiency disorders (PIDs) are expensive, time-consuming, and not easily accessible in developing countries. Therefore, we studied the feasibility of a customized single nucleotide variant (SNV) microarray that we developed to detect disease-causing variants and copy number variation (CNV) in patients with PIDs for only 40 Euros.Methods: Probes were custom-designed to genotype 9,415 variants of 277 PID-related genes, and were added to the genome-wide Illumina Global Screening Array (GSA). Data analysis of GSA was performed using Illumina GenomeStudio 2.0, Biodiscovery Nexus 10.0, and R-3.4.4 software. Validation of genotype calling was performed by comparing the GSA with whole-genome sequencing (WGS) data of 56 non-PID controls. DNA samples of 95 clinically diagnosed PID patients, of which 60 patients (63%) had a genetically established diagnosis (by Next-Generation Sequencing (NGS) PID panels or Sanger sequencing), were analyzed to test the performance of the GSA. The additional SNVs detected by GSA were validated by Sanger sequencing.Results: Genotype calling of the customized array had an accuracy rate of 99.7%. The sensitivity for detecting rare PID variants was high (87%). The single sample replication in two runs was high (94.9%). The customized GSA was able to generate a genetic diagnosis in 37 out of 95 patients (39%). These 37 patients included 29 patients in whom the genetic variants were confirmed by conventional methods (26 patients by SNV and 3 by CNV analysis), while in 8 patients a new genetic diagnosis was established (6 patients by SNV and 2 patients suspected for leukemia by CNV analysis). Twenty-eight patients could not be detected due to the limited coverage of the custom probes. However, the diagnostic yield can potentially be increased when newly updated variants are added.Conclusion: Our robust customized GSA seems to be a promising first-line rapid screening tool for PIDs at an affordable price, which opens opportunities for low-cost genetic testing in developing countries. The technique is scalable, allows numerous new genetic variants to be added, and offers the potential for genetic testing not only in PIDs, but also in many other genetic diseases. Show less
Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results... Show moreLeukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1 , PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease. Show less
BACKGROUND: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a... Show moreBACKGROUND: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a plausible explanation for the overlap, and in this study we tested whether genetic variation in the major histocompatibility complex (MHC), which is associated with risk for autoimmune diseases, is also associated with risk for depression.METHODS: We fine-mapped the classical MHC (chr6: 29.6-33.1 Mb), imputing 216 human leukocyte antigen (HLA) alleles and 4 complement component 4 (C4) haplotypes in studies from the Psychiatric Genomics Consortium Major Depressive Disorder Working Group and the UK Biobank. The total sample size was 45,149 depression cases and 86,698 controls. We tested for association between depression status and imputed MHC variants, applying both a region-wide significance threshold (3.9 x 10(-6) ) and a candidate threshold (1.6 x 10(-4) ).RESULTS: No HLA alleles or C4 haplotypes were associated with depression at the region-wide threshold. HLAB*08:01 was associated with modest protection for depression at the candidate threshold for testing in HLA genes in the meta-analysis (odds ratio = 0.98, 95% confidence interval = 0.97-0.99).CONCLUSIONS: We found no evidence that an increased risk for depression was conferred by HLA alleles, which play a major role in the genetic susceptibility to autoimmune diseases, or C4 haplotypes, which are strongly associated with schizophrenia. These results suggest that any HLA or C4 variants associated with depression either are rare or have very modest effect sizes. Show less
Insights into individual differences in gene expression and its heritability (h(2)) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52... Show moreInsights into individual differences in gene expression and its heritability (h(2)) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52,844 genes measured in whole blood in the largest twin RNA-Seq sample to date (1497 individuals including 459 monozygotic twin pairs and 150 dizygotic twin pairs) from classical twin modeling and identity-by-state-based approaches. We estimated for each gene h(total)(2), composed of cis-heritability (h(cis)(2), the variance explained by single nucleotide polymorphisms in the cis-window of the gene), and trans-heritability (h(res)(2), the residual variance explained by all other genome-wide variants). Mean h(total)(2) was 0.26, which was significantly higher than heritability estimates earlier found in a microarray-based study using largely overlapping (>60%) RNA samples (mean h(2) = 0.14, p = 6.15 x 10(-258)). Mean h(cis)(2) was 0.06 and strongly correlated with beta of the top cis expression quantitative loci (eQTL, rho = 0.76, p < 10(-308)) and with estimates from earlier RNA-Seq-based studies. Mean h(res)(2) was 0.20 and correlated with the beta of the corresponding trans-eQTL (rho = 0.04, p < 1.89 x 10(-3)) and was significantly higher for genes involved in cytokine-cytokine interactions (p = 4.22 x 10(-15)), many other immune system pathways, and genes identified in genome-wide association studies for various traits including behavioral disorders and cancer. This study provides a thorough characterization of cis- and trans-h(2) estimates of gene expression, which is of value for interpretation of GWAS and gene expression studies. Show less