Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

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 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.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci.

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Fat metabolism is associated with telomere length in six population-based studies.

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 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 FlowFISH. 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 Pmeta = 6.5 × 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 × 10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (P = 1.9 × 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.

Whole-exome sequencing of 14 389 individuals from the ESP and CHARGE consortia identifies novel rare variation associated with hemostatic factors.

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors.

Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues.

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

Genetic and clinical determinants of abdominal aortic diameter: genome-wide association studies, exome array data and Mendelian randomization study.

Progressive dilation of the infrarenal aortic diameter is a consequence of the ageing process and is considered the main determinant of abdominal aortic aneurysm (AAA). We aimed to investigate the genetic and clinical determinants of abdominal aortic diameter (AAD). We conducted a meta-analysis of genome-wide association studies in 10 cohorts (n = 13 542) imputed to the 1000 Genome Project reference panel including 12 815 subjects in the discovery phase and 727 subjects [Partners Biobank cohort 1 (PBIO)] as replication. Maximum anterior-posterior diameter of the infrarenal aorta was used as AAD. We also included exome array data (n = 14 480) from seven epidemiologic studies. Single-variant and gene-based associations were done using SeqMeta package. A Mendelian randomization analysis was applied to investigate the causal effect of a number of clinical risk factors on AAD. In genome-wide association study (GWAS) on AAD, rs74448815 in the intronic region of LDLRAD4 reached genome-wide significance (beta = -0.02, SE = 0.004, P-value = 2.10 × 10-8). The association replicated in the PBIO1 cohort (P-value = 8.19 × 10-4). In exome-array single-variant analysis (P-value threshold = 9 × 10-7), the lowest P-value was found for rs239259 located in SLC22A20 (beta = 0.007, P-value = 1.2 × 10-5). In the gene-based analysis (P-value threshold = 1.85 × 10-6), PCSK5 showed an association with AAD (P-value = 8.03 × 10-7). Furthermore, in Mendelian randomization analyses, we found evidence for genetic association of pulse pressure (beta = -0.003, P-value = 0.02), triglycerides (beta = -0.16, P-value = 0.008) and height (beta = 0.03, P-value < 0.0001), known risk factors for AAA, consistent with a causal association with AAD. Our findings point to new biology as well as highlighting gene regions in mechanisms that have previously been implicated in the genetics of other vascular diseases.

Advanced glycation end products measured by skin autofluorescence and subclinical cardiovascular disease: the Rotterdam Study.

BACKGROUND: Advanced glycation end products (AGEs) have been linked to cardiovascular disease (CVD), especially coronary heart disease (CHD), but their role in CVD pathogenesis remains unclear. Therefore, we investigated cross-sectional associations of skin AGEs with subclinical atherosclerosis, arterial stiffness, and hypertension after confirming their relation with CHD. METHODS: In the population-based Rotterdam Study, skin AGEs were measured as skin autofluorescence (SAF). Prevalent MI was obtained from digital medical records. Carotid plaques, carotid intima-media thickness (IMT), coronary artery calcification (CAC), pulse wave velocity (PWV), and hypertension were assessed. Associations of SAF with endophenotypes were investigated in logistic and linear regression models adjusting for common cardiovascular risk factors. Effect modification by sex, diabetes mellitus, and chronic kidney disease (CKD) was tested. RESULTS: 3001 participants were included (mean age 73 (SD 9) years, 57% women). One unit higher SAF was associated with the presence of carotid plaques (OR 1.2 (0.92, 1.57)), a higher max IMT (0.08 SD (0.01, 0.15)), higher CAC (OR 2.2 (1.39, 3.48)), and PWV (0.09 SD (0.01, 0.16)), but not with hypertension (OR 0.99 (0.81, 1.21)). The associations with endophenotypes were more pronounced in men and participants with diabetes or CKD with significant interactions. CONCLUSIONS: Previously documented associations between SAF and CVD, also found in our study, may be explained by the endophenotypes atherosclerosis and arterial stiffness, especially in men and individuals with diabetes or CKD, but not by hypertension. Longitudinal studies are needed to replicate these findings and to test if SAF is an independent risk factor or biomarker of CVD. TRIAL REGISTRATION: The Rotterdam Study has been entered into the Netherlands National Trial Register (NTR; www.trialregister.nl ) and the WHO International Clinical Trials Registry Platform (ICTRP; www.who.int/ictrp/network/primary/en/ ) under shared catalogue number NTR6831.

The gut microbiome and child mental health: A population-based study.

The link between the gut microbiome and the brain has gained increasing scientific and public interest for its potential to explain psychiatric risk. While differences in gut microbiome composition have been associated with several mental health problems, evidence to date has been largely based on animal models and human studies with modest sample sizes. In this cross-sectional study in 1,784 ten-year-old children from the multi-ethnic, population-based Generation R Study, we aimed to characterize associations of the gut microbiome with child mental health problems. Gut microbiome was assessed from stool samples using 16S rRNA sequencing. We focused on overall psychiatric symptoms as well as with specific domains of emotional and behavioral problems, assessed via the maternally rated Child Behavior Checklist. While we observed lower gut microbiome diversity in relation to higher overall and specific mental health problems, associations were not significant. Likewise, we did not identify any taxonomic feature associated with mental health problems after multiple testing correction, although suggestive findings indicated depletion of genera previously associated with psychiatric disorders, including Hungatella, Anaerotruncus and Oscillospiraceae. The identified compositional abundance differences were found to be similar across all mental health problems. Finally, we did not find significant enrichment for specific microbial functions in relation to mental health problems. In conclusion, based on the largest sample examined to date, we do not find clear evidence of associations between gut microbiome diversity, taxonomies or functions and mental health problems in the general pediatric population. In future, the use of longitudinal designs with repeated measurements of microbiome and psychiatric outcomes will be critical to identify whether and when associations between the gut microbiome and mental health emerge across development and into adulthood.

Coronary Artery Calcium Score and Polygenic Risk Score for the Prediction of Coronary Heart Disease Events.

Importance: Coronary artery calcium score and polygenic risk score have each separately been proposed as novel markers to identify risk of coronary heart disease (CHD), but no prior studies have directly compared these markers in the same cohorts. Objective: To evaluate change in CHD risk prediction when a coronary artery calcium score, a polygenic risk score, or both are added to a traditional risk factor-based model. Design, Setting, and Participants: Two observational population-based studies involving individuals aged 45 years through 79 years of European ancestry and free of clinical CHD at baseline: the Multi-Ethnic Study of Atherosclerosis (MESA) study involved 1991 participants at 6 US centers and the Rotterdam Study (RS) involved 1217 in Rotterdam, the Netherlands. Exposure: Traditional risk factors were used to calculate CHD risk (eg, pooled cohort equations [PCEs]), computed tomography for the coronary artery calcium score, and genotyped samples for a validated polygenic risk score. Main Outcomes and Measures: Model discrimination, calibration, and net reclassification improvement (at the recommended risk threshold of 7.5%) for prediction of incident CHD events were assessed. Results: The median age was 61 years in MESA and 67 years in RS. Both log (coronary artery calcium+1) and polygenic risk score were significantly associated with 10-year risk of incident CHD (hazards ratio per SD, 2.60; 95% CI, 2.08-3.26 and 1.43; 95% CI, 1.20-1.71, respectively), in MESA. The C statistic for the coronary artery calcium score was 0.76 (95% CI, 0.71-0.79) and for the polygenic risk score, 0.69 (95% CI, 0.63-0.71). The change in the C statistic when each was added to the PCEs was 0.09 (95% CI, 0.06-0.13) for the coronary artery calcium score, 0.02 (95% CI, 0.00-0.04) for the polygenic risk score, and 0.10 (95% CI, 0.07-0.14) for both. Overall categorical net reclassification improvement was significant when the coronary artery calcium score (0.19; 95% CI, 0.06-0.28) but was not significant when the polygenic risk score (0.04; 95% CI, -0.05 to 0.10) was added to the PCEs. Calibration of the PCEs and models with coronary artery calcium and/or polygenic risk scores was adequate (all χ2<20). Subgroup analysis stratified by the median age demonstrated similar findings. Similar findings were observed for 10-year risk in RS and in longer-term follow-up in MESA (median, 16.0 years). Conclusions and Relevance: In 2 cohorts of middle-aged to older adults from the US and the Netherlands, the coronary artery calcium score had better discrimination than the polygenic risk score for risk prediction of CHD. In addition, the coronary artery calcium score but not the polygenic risk score significantly improved risk discrimination and risk reclassification for CHD when added to traditional risk factors.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Microbiomics, Metabolomics, Predicted Metagenomics, and Hepatic Steatosis in a Population-Based Study of 1,355 Adults.

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, 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∙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: ß = -65; Ruminococcus Gauvreauiigroup: ß = 62; and Ruminococcus Gnavusgroup: ß = 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.

Association between an AMH promoter polymorphism and serum AMH levels in PCOS patients.

STUDY QUESTION: Do polymorphisms in the anti-Müllerian hormone (AMH) promoter have an effect on AMH levels in patients with polycystic ovary syndrome (PCOS)? SUMMARY ANSWER: We have identified a novel AMH promoter polymorphism rs10406324 that is associated with lower serum AMH levels and is suggested to play a role in the mechanism of regulation of AMH gene expression in women. WHAT IS KNOWN ALREADY: Follicle number is positively correlated with serum AMH levels, reflected by elevated AMH levels in women with PCOS. In addition, it is suggested that AMH production per follicle is higher in women with PCOS than in normo-ovulatory women, implying an altered regulation of AMH in PCOS. STUDY DESIGN, SIZE, DURATION: A discovery cohort of 655 PCOS women of Northern European ancestry and both an internal and external validation PCOS cohort (n = 458 and n = 321, respectively) were included in this study. Summary-level data of an AMH genome-wide association study meta-analysis including 7049 normo-ovulatory women was included as a control cohort. A genetic approach was taken through association analysis and in silico analysis of the associated variants in the AMH promoter. In vitro analysis was performed to investigate the functional mechanisms. PARTICIPANTS/MATERIALS, SETTING, METHODS: All common two-allelic single-nucleotide polymorphisms (SNPs) in the region Chr19:2 245 353-2 250 827 bp (Build 37) were selected for the analysis. Linear regression analyses were performed to determine the association between SNPs in the AMH promoter region and serum AMH levels. For the in silico analysis, the webtools 'HaploReg' v4.1 for ENCODE prediction weight matrices and 'atSNP' were used. In vitro analysis was performed using KK1 cells, a mouse granulosa cell line and COV434 cells, a human granulosa tumor cell line. Cells were transfected with the reference or the variant human AMH promoter reporter construct together with several transcription factors (TFs). Dual-Glo® Luciferase Assay was performed to measure the luciferase activity. MAIN RESULTS AND THE ROLE OF CHANCE: Polymorphism rs10406324 was significantly associated with serum AMH levels in all three PCOS cohorts. Carriers of the minor allele G had significantly lower log-transformed serum AMH levels compared to non-carriers (P = 8.58 × 10-8, P = 1.35 × 10-3 and P = 1.24 × 10-3, respectively). This result was validated in a subsequent meta-analysis (P = 3.24 × 10-12). Interestingly, rs10406324 was not associated with follicle count, nor with other clinical traits. Also, in normo-ovulatory women, the minor allele of this variant was associated with lower serum AMH levels (P = 1.04 × 10-5). These findings suggest that polymorphism rs10406324 plays a role in the regulation of AMH expression, irrespective of clinical background. In silico analysis suggested a decreased binding affinity of the TFs steroidogenenic factor 1, estrogen-related receptor alpha and glucocorticoid receptor to the minor allele G variant, however in vitro analysis did not show a difference in promoter activity between the A and G allele. LIMITATIONS, REASONS FOR CAUTION: Functional analyses were performed in a mouse and a human granulosa cell line using an AMH promoter reporter construct. This may have limited assessment of the impact of the polymorphism on higher order chromatin structures. Human granulosa cells generated from induced pluripotent stem cells, combined with gene editing, may provide a method to elucidate the exact mechanism behind the decrease in serum AMH levels in carriers of the -210 G allele. We acknowledge that the lack of follicle number in the external validation and the control cohort is a limitation of the paper. Although we observed that the association between rs10406324 and AMH levels was independent of follicle number in our discovery and internal validation PCOS cohorts, we cannot fully rule out that the observed effects on serum AMH levels are, in part, caused by differences in follicle number. WIDER IMPLICATIONS OF THE FINDINGS: These results suggest that variations in serum AMH levels are not only caused by differences in follicle number but also by genetic factors. Therefore, the genetic context should be taken into consideration when assessing serum AMH levels in women. This may have clinical consequences when serum AMH levels are used as a marker for the polycystic ovarian morphology phenotype. STUDY FUNDING/COMPETING INTEREST(S): No external funding was used. J.S.E.L. has received consultancy fees from the following companies: Ferring, Roche Diagnostics and Ansh Labs and has received travel reimbursement from Ferring. J.A.V. has received royalties from AMH assays, paid to the institute/lab with no personal financial gain. The other authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

Serum Phosphate, BMI, and Body Composition of Middle-Aged and Older Adults: A Cross-Sectional Association Analysis and Bidirectional Mendelian Randomization Study.

BACKGROUND: Observational studies have reported associations between serum phosphate and BMI in specific clinical settings, but the nature of this relation in the general population is unclear. OBJECTIVES: The aim of this study was twofold: to investigate the association between serum phosphate and BMI and body composition, as well as to explore evidence of causality through a bidirectional one-sample Mendelian randomization (MR) in the population-based Rotterdam Study (RS). METHODS: Observational associations between phosphate (mg/dL) and BMI, lean mass, and fat percentage (fat%), estimated by DXA, were analyzed using multivariable regression models in 9202 participants aged 45-100 y from 3 RS cohorts. The role of serum leptin was examined in a subgroup of 1089 participants. For MR analyses, allele scores with 6 single-nucleotide polymorphisms (SNPs) for phosphate and 905 SNPs for BMI were constructed in 7983 participants. RESULTS: Phosphate was inversely associated with BMI in the total population (ß: -0.89; 95% CI: -1.17, -0.62), and stronger in women (ß: -1.92; 95% CI: -2.20, -1.65) than in men (ß: -0.37; 95% CI: -0.68, -0.06) (P-interaction < 0.05). Adjustment for leptin did not change results in men. In women, adjustment for leptin attenuated the association, but it was not abolished (ß: -0.94; 95% CI: -1.45, -0.42). Phosphate was inversely associated with fat%, but not with lean mass, in both sexes. MR analyses suggested a causal effect of BMI on serum phosphate (ß: -0.01; 95% CI: -0.02, 0.00) but not vice versa. CONCLUSIONS: Serum phosphate was inversely associated with BMI and fat% in a population-based study of middle-aged and older adults, with a stronger effect in women than in men. Adjusting for leptin attenuated this relation in women only. MR results suggest a causal effect of BMI on phosphate but not vice versa. An underlying sex dimorphism in phosphate homeostasis should be further explored.

Role of the Microbiome in Regulating Bone Metabolism and Susceptibility to Osteoporosis.

The human microbiota functions at the interface between diet, medication-use, lifestyle, host immune development and health. It is therefore closely aligned with many of the recognised modifiable factors that influence bone mass accrual in the young, and bone maintenance and skeletal decline in older populations. While understanding of the relationship between micro-organisms and bone health is still in its infancy, two decades of broader microbiome research and discovery supports a role of the human gut microbiome in the regulation of bone metabolism and pathogenesis of osteoporosis as well as its prevention and treatment. Pre-clinical research has demonstrated biological interactions between the microbiome and bone metabolism. Furthermore, observational studies and randomized clinical trials have indicated that therapeutic manipulation of the microbiota by oral administration of probiotics may influence bone turnover and prevent bone loss in humans. In this paper, we summarize the content, discussion and conclusions of a workshop held by the Osteoporosis and Bone Research Academy of the Royal Osteoporosis Society in October, 2020. We provide a detailed review of the literature examining the relationship between the microbiota and bone health in animal models and in humans, as well as formulating the agenda for key research priorities required to advance this field. We also underscore the potential pitfalls in this research field that should be avoided and provide methodological recommendations to facilitate bridging the gap from promising concept to a potential cause and intervention target for osteoporosis.

Genomic analysis of diet composition finds novel loci and associations with health and lifestyle.

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 identified 21 unique, approximately independent lead SNPs. Fourteen lead SNPs are uniquely associated with one macronutrient at genome-wide significance (P < 5 × 10-8), while five of the 21 lead SNPs reach suggestive significance (P < 1 × 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 (rg ≈ 0.15-0.5). In contrast, relative carbohydrate and sugar intake have negative genetic correlations with waist circumference, waist-hip ratio, and neighborhood deprivation (|rg| ≈ 0.1-0.3) and positive genetic correlations with physical activity (rg ≈ 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 (rg ≈-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.

Correction: Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

[This corrects the article DOI: 10.1371/journal.pgen.1007813.].

Vitamin K antagonist anticoagulant usage is associated with increased incidence and progression of osteoarthritis.

OBJECTIVES: Vitamin 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. METHODS: We 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. RESULTS: Acenocoumarol 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). CONCLUSIONS: These 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.

Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene.

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 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.

Mendelian randomization study on vitamin D levels and osteoarthritis risk: a concise report.

OBJECTIVE: The role of vitamin D in OA is unclear and previous epidemiological studies have provided inconsistent results. We conducted a two-sample Mendelian randomization (MR) study to investigate the causal relationship between genetically determined serum vitamin D levels and hip/knee OA. METHODS: Six single-nucleotide polymorphisms (SNPs) associated with vitamin D levels in the Study of Underlying Genetic Determinants of Vitamin D and Highly Related Traits Consortium were selected as instrumental variables. Summary statistics of the SNPs effects on OA were derived from the Iceland and UK Biobank, comprising 23 877 knee OA cases, 17 151 hip OA cases and >562 000 controls. The control samples match the OA cases in age, sex and county of origin. RESULTS: The MR analyses showed no causal association between genetically determined vitamin D levels and knee OA [odds ratio (OR) 1.03 (95% CI 0.84, 1.26)] or hip OA [OR 1.06 (95% CI 0.83, 1.35)]. CONCLUSION: Genetic variations associated with low vitamin D serum levels are not associated with increased risk of hip or knee OA in community-dwelling older adults, suggesting that vitamin D levels are not causally linked to OA. It is therefore unlikely that vitamin D supplementation protects against hip or knee OA.