Association between prenatal alcohol exposure and children's facial shape: a prospective population-based cohort study.

STUDY QUESTION: Is there an association between low-to-moderate levels of prenatal alcohol exposure (PAE) and children's facial shape? SUMMARY ANSWER: PAE before and during pregnancy, even at low level (<12 g of alcohol per week), was found associated with the facial shape of children, and these associations were found attenuated as children grow older. WHAT IS KNOWN ALREADY: High levels of PAE during pregnancy can have significant adverse associations with a child's health development resulting in recognizably abnormal facial development. STUDY DESIGN, SIZE, DURATION: This study was based on the Generation R Study, a prospective cohort from fetal life onwards with maternal and offspring data. We analyzed children 3-dimensional (3D) facial images taken at ages 9 (n = 3149) and 13 years (n = 2477) together with the data of maternal alcohol consumption. PARTICIPANTS/MATERIALS, SETTING, METHODS: We defined six levels of PAE based on the frequency and dose of alcohol consumption and defined three tiers based on the timing of alcohol exposure of the unborn child. For the image analysis, we used 3D graph convolutional networks for non-linear dimensionality reduction, which compressed the high-dimensional images into 200 traits representing facial morphology. These 200 traits were used for statistical analysis to search for associations with PAE. Finally, we generated heatmaps to display the facial phenotypes associated with PAE. MAIN RESULTS AND THE ROLE OF CHANCE: The results of the linear regression in the 9-year-old children survived correction for multiple testing with false discovery rate (FDR). In Tier 1 where we examined PAE only before pregnancy (exposed N = 278, unexposed N = 760), we found three traits survived FDR correction. The lowest FDR-P is 1.7e-05 (beta = 0.021, SE = 0.0040) in Trait #29; In Tier 2b where we examine any PAE during first trimester (exposed N = 756; unexposed N = 760), we found eight traits survived FDR correction. The lowest FDR-P is 9.0e-03 (beta = -0.013, SE = 0.0033) in Trait #139. Moreover, more statistically significant facial traits were found in higher levels of PAE. No FDR-significant results were found in the 13-year-old children. We map these significant traits back to the face, and found the most common detected facial phenotypes included turned-up nose tip, shortened nose, turned-out chin, and turned-in lower-eyelid-related regions. LIMITATIONS, REASONS FOR CAUTION: We had no data for alcohol consumption more than three months prior to pregnancy and thus do not know if maternal drinking had chronic effects. The self-reported questionnaire might not reflect accurate alcohol measurements because mothers may have denied their alcohol consumption. WIDER IMPLICATIONS OF THE FINDINGS: Our results imply that facial morphology, such as quantified by the approach we proposed here, can be used as a biomarker in further investigations. Furthermore, our study suggests that for women who are pregnant or want to become pregnant soon, should quit alcohol consumption several months before conception and completely during pregnancy to avoid adverse health outcomes in the offspring. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by Erasmus Medical Centre, Rotterdam, the Erasmus University Rotterdam, and the Netherlands Organization for Health Research. V.W.V.J. reports receipt of funding from the Netherlands Organization for Health Research (ZonMw 90700303). W.J.N. is a founder, a scientific lead, and a shareholder of Quantib BV. TRIAL REGISTRATION NUMBER: N/A.

Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation.

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10-8) and mtDNA replication (p = 1.2 × 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10-4).

The association of BMI and physical activity on acetabular dysplasia in children.

BACKGROUND: Acetabular dysplasia is an important pre-disposing factor for osteoarthritis of the hip. However, it is not completely known how acetabular dysplasia develops during childhood. OBJECTIVE: To study the prevalence of acetabular dysplasia and its association with body mass index (BMI) and physical activity in 9 year old children. DESIGN: The population for this cross-sectional study was drawn from the ongoing prospective cohort study: Generation R. 9,778 mothers with a delivery date from March 2002 until January 2006 were enrolled. In a random subgroup of these children Dual-energy X-ray absorptiometry (DXA) scanning was performed at age 9. EXPOSURES: BMI, standardized for the Dutch population and categorized in four groups based on extended international Obesity Task Force cut-offs: underweight, normal, overweight and obesity. Physical activity was based on time spent on playing outdoors, playing sports and walking/cycling to school. MAIN OUTCOMES AND MEASURES: The degree of acetabular dysplasia was determined with the centre-edge angle (CEA) and acetabular depth-width ratio (ADR) in DXA images of the hip. RESULTS: 1,188 DXA images of children's hips were available for analysis. The median age of the children was 9.86 years. Prevalence of dysplasia and mild dysplasia was respectively 6.3%; 25.6% with CEA and 4.8%; 25.0% with ADR. BMI was negatively associated with mild dysplasia (OR 0.80 CI 0.71-0.90). Obese children showed less mild dysplasia compared to normal children (OR 0.48 CI 0.24-0.97) in unadjusted analysis. Physical activity represented by walking to school showed a statistically significant negative association with mild dysplasia (OR 0.87 CI 0.76-0.99). After adjustment for age, ethnicity, sex, first born, breech presentation, birthweight, gestational age and Caesarean section, the patterns of association with dysplasia remained for both BMI and physical activity. CONCLUSIONS: In this study, being overweight and light physical activity were negatively associated with the development of (mild) acetabular dysplasia at the age of 9 years.

Maternal and fetal genetic contribution to gestational weight gain.

BACKGROUND: Clinical recommendations to limit gestational weight gain (GWG) imply high GWG is causally related to adverse outcomes in mother or offspring, but GWG is the sum of several inter-related complex phenotypes (maternal fat deposition and vascular expansion, placenta, amniotic fluid and fetal growth). Understanding the genetic contribution to GWG could help clarify the potential effect of its different components on maternal and offspring health. Here we explore the genetic contribution to total, early and late GWG. PARTICIPANTS AND METHODS: A genome-wide association study was used to identify maternal and fetal variants contributing to GWG in up to 10 543 mothers and 16 317 offspring of European origin, with replication in 10 660 mothers and 7561 offspring. Additional analyses determined the proportion of variability in GWG from maternal and fetal common genetic variants and the overlap of established genome-wide significant variants for phenotypes relevant to GWG (for example, maternal body mass index (BMI) and glucose, birth weight). RESULTS: Approximately 20% of the variability in GWG was tagged by common maternal genetic variants, and the fetal genome made a surprisingly minor contribution to explain variation in GWG. Variants near the pregnancy-specific beta-1 glycoprotein 5 (PSG5) gene reached genome-wide significance (P=1.71 × 10-8) for total GWG in the offspring genome, but did not replicate. Some established variants associated with increased BMI, fasting glucose and type 2 diabetes were associated with lower early, and higher later GWG. Maternal variants related to higher systolic blood pressure were related to lower late GWG. Established maternal and fetal birth weight variants were largely unrelated to GWG. CONCLUSIONS: We found a modest contribution of maternal common variants to GWG and some overlap of maternal BMI, glucose and type 2 diabetes variants with GWG. These findings suggest that associations between GWG and later offspring/maternal outcomes may be due to the relationship of maternal BMI and diabetes with GWG.

Dietary acid load, trabecular bone integrity, and mineral density in an ageing population: the Rotterdam study.

We studied the relation between a diet that is high in acid-forming nutrients (e.g. proteins) and low in base-forming nutrients (e.g. potassium) and bone structure. We showed a negative relation, which was more prominent if proteins were of animal rather than of vegetable origin and if intake of dietary fibre was high. INTRODUCTION: Studies on dietary acid load (DAL) and fractures have shown inconsistent results. Associations between DAL, bone mineral density (BMD) and trabecular bone integrity might play a role in these inconsistencies and might be influenced by renal function and dietary fibre intake. Therefore, our aim was to study (1) associations of DAL with BMD and with the trabecular bone score (TBS) and (2) the potential influence of renal function and dietary fibre in these associations. METHODS: Dutch individuals aged 45 years and over (n = 4672) participating in the prospective cohort of the Rotterdam study were included. Based on food frequency questionnaires, three indices of DAL were calculated: the net endogenous acid production (NEAP) and the ratios of vegetable or animal protein and potassium (VegPro/K and AnPro/K). Data on lumbar spinal TBS and BMD were derived from dual-energy X-ray absorptiometry measurements. RESULTS: Independent of confounders, NEAP and AnPro/K, but not VegPro/K, were associated with low TBS (standardized ß (95%) = -0.04 (-0.07, -0.01) and -0.08 (-0.11, -0.04)) but not with BMD. Associations of AnPro/K and VegPro/K with TBS were non-linear and differently shaped. Unfavourable associations between NEAP, BMD and TBS were mainly present in subgroups with high fibre intake. CONCLUSIONS: High NEAP was associated with low TBS. Associations of AnPro/K and VegPro/K and TBS were non-linear and differently shaped. No significant associations of DAL with BMD were observed, nor was there any significant interaction between DAL and renal function. Mainly in participants with high intake of dietary fibre, DAL might be detrimental to bone.

Somatic, positive and negative domains of the Center for Epidemiological Studies Depression (CES-D) scale: a meta-analysis of genome-wide association studies.

BACKGROUND: Major depressive disorder (MDD) is moderately heritable, however genome-wide association studies (GWAS) for MDD, as well as for related continuous outcomes, have not shown consistent results. Attempts to elucidate the genetic basis of MDD may be hindered by heterogeneity in diagnosis. The Center for Epidemiological Studies Depression (CES-D) scale provides a widely used tool for measuring depressive symptoms clustered in four different domains which can be combined together into a total score but also can be analysed as separate symptom domains. METHOD: We performed a meta-analysis of GWAS of the CES-D symptom clusters. We recruited 12 cohorts with the 20- or 10-item CES-D scale (32 528 persons). RESULTS: One single nucleotide polymorphism (SNP), rs713224, located near the brain-expressed melatonin receptor (MTNR1A) gene, was associated with the somatic complaints domain of depression symptoms, with borderline genome-wide significance (p discovery = 3.82 × 10-8). The SNP was analysed in an additional five cohorts comprising the replication sample (6813 persons). However, the association was not consistent among the replication sample (p discovery+replication = 1.10 × 10-6) with evidence of heterogeneity. CONCLUSIONS: Despite the effort to harmonize the phenotypes across cohorts and participants, our study is still underpowered to detect consistent association for depression, even by means of symptom classification. On the contrary, the SNP-based heritability and co-heritability estimation results suggest that a very minor part of the variation could be captured by GWAS, explaining the reason of sparse findings.

Does fetal smoke exposure affect childhood bone mass? The Generation R Study.

UNLABELLED: We assessed the intrauterine influence of maternal smoking on childhood bone mass by comparing parental prenatal and postnatal smoking habits. We observed higher bone mass in children exposed to maternal smoking, explained by higher body weight. Maternal smoking or related lifestyle factors may affect childhood weight gain rather than skeletal growth. INTRODUCTION: Maternal smoking during pregnancy may adversely affect bone health in later life. By comparing the associations of maternal and paternal smoking and of prenatal and postnatal exposure with childhood bone measures, we aimed to explore whether the suggested association could be explained by fetal programming or reflects confounding by familial factors. METHODS: In 5565 mothers, fathers and children participating in a population-based prospective cohort study, parental smoking habits during pregnancy and current household smoking habits were assessed by postal questionnaires. Total body bone mineral content (BMC), bone area (BA) and bone mineral density (BMD) were measured by dual-energy X-ray absorptiometry (DXA) at the median age of 6.0 years (IQR 0.37). RESULTS: In confounder-adjusted models, maternal smoking during pregnancy was associated with a higher BMC of 11.6 g (95 % confidence interval (CI) 5.6, 17.5), a larger BA of 9.7 cm(2) (95 % CI 3.0, 16.4), a higher BMD of 6.7 g/cm(2) (95 % CI 2.4, 11.0) and a higher BMC of 5.4 g (95 % CI 1.3, 9.6) adjusted for BA of the child. Current weight turned out to mediate these associations. Among mothers who did not smoke, paternal smoking did not show evident associations with childhood bone measures. Also, household smoking practices during childhood were not associated with childhood bone measures. CONCLUSIONS: Our results do not support the hypothesis of fetal smoke exposure affecting childhood bone mass via intrauterine mechanisms. Maternal smoking or related lifestyle factors may affect childhood weight gain rather than skeletal growth.

Infant dietary patterns and bone mass in childhood: the Generation R Study.

UNLABELLED: Early life nutrition affects peak bone mass attainment. In this prospective cohort study, children with high adherence to a "dairy and whole grains" pattern in infancy had higher bone mineral density at the age of 6 years. Although the observed effects are small, our study provides insight into mechanisms linking early nutrition to bone acquisition in childhood. INTRODUCTION: Nutrition in early life may affect peak bone mass attainment. Previous studies on childhood nutrition and skeletal health mainly focused on individual nutrients, which does not consider the cumulative effects of nutrients. We investigated the associations between dietary patterns in infancy and childhood bone health. METHODS: This study included 2850 children participating in a population-based prospective cohort study. Dietary information was obtained from a food frequency questionnaire at the age of 13 months. Using principal component analysis, three major dietary patterns were extracted, explaining in total 30% of the variation in dietary intake. At the age of 6 years, a total body dual-energy X-ray absorptiometry (DXA) scan was performed, and bone mineral density (BMD), bone mineral content (BMC), area-adjusted BMC (aBMC), and bone area (BA) were analyzed. RESULTS: Higher adherence score to a "dairy and whole grains" pattern was positively associated with BMD and aBMC, but not with BMC and BA. Accordingly, children in the highest quartile of the "dairy and whole grains" pattern had higher BMD (difference 3.98 mg/cm(2), 95% confidence interval (CI) 0.36 to 7.61) and aBMC (difference 4.96 g, 95% CI 1.27 to 8.64) than children in the lowest quartile. Stratification for vitamin D supplementation showed that the positive associations between the "dairy and whole grains" pattern and bone outcomes were only observed in children who did not receive vitamin D supplementation. A "potatoes, rice, and vegetables" and a "refined grains and confectionery" pattern were not consistently associated with bone outcomes. CONCLUSIONS: An infant dietary pattern characterized by high intakes of dairy and cheese, whole grains, and eggs is positively associated with bone development in childhood. Further research is needed to investigate the consequences for bone health in later life.

Association of polymorphisms in the beta-2 adrenergic receptor gene with fracture risk and bone mineral density.

UNLABELLED: Signaling through the beta-2 adrenergic receptor (B2AR) on the osteoblast influences bone remodeling in rodents. In the B2AR gene, three polymorphisms influence receptor function. We show that these polymorphisms are not associated with fracture risk or bone mineral density in the UCP, Rotterdam Study, and GEFOS cohorts. INTRODUCTION: Signaling through the beta-2 adrenergic receptor (B2AR) on the osteoblast influences bone remodeling in rodents. In the B2AR gene, three polymorphisms are known to influence receptor function in vitro and in vivo (rs1042713, rs1042714, and rs1800888). We examined the role of these polymorphisms in the B2AR gene on human bone metabolism. METHODS: We performed nested case-control studies to determine the association of these polymorphisms with fracture risk in the Utrecht Cardiovascular Pharmacogenetics (UCP) cohort and in three cohorts of the Rotterdam Study. We also determined the association of these polymorphisms with bone mineral density (BMD) in the GEFOS Consortium. UCP contains drug-dispensing histories from community pharmacies linked to national registrations of hospital discharges in the Netherlands. The Rotterdam Study is a prospective cohort study investigating demographics and risk factors of chronic diseases. GEFOS is a large international collaboration studying the genetics of osteoporosis. Fractures were defined by ICD-9 codes 800-829 in the UCP cohort (158 cases and 2617 unmatched controls) and by regular X-ray examinations, general practitioner, and hospital records in the Rotterdam Study (2209 cases and 8559 unmatched controls). BMD was measured at the femoral neck and lumbar spine using dual-energy X-ray absorptiometry in GEFOS (N = 32,961). RESULTS: Meta-analysis of the two nested case-control studies showed pooled odds ratios of 0.98 (0.91-1.05, p = 0.52), 1.04 (0.97-1.12, p = 0.28), and 1.16 (0.83-1.62, p = 0.38) for the associations between rs1042713, rs1042714, and rs1800888 per minor allele and fractures, respectively. There were no significant associations of the polymorphisms and BMD in GEFOS. CONCLUSION: In conclusion, polymorphisms in the beta-2 adrenergic receptor gene are not associated with fracture risk or BMD.

Vertebral Scheuermann's disease in Europe: prevalence, geographic variation and radiological correlates in men and women aged 50 and over.

UNLABELLED: In 27 centres across Europe, the prevalence of deforming spinal Scheuermann's disease in age-stratified population-based samples of over 10,000 men and women aged 50+ averaged 8% in each sex, but was highly variable between centres. Low DXA BMD was un-associated with Scheuermann's, helping the differential diagnosis from osteoporosis. INTRODUCTION: This study aims to assess the prevalence of Scheuermann's disease of the spine across Europe in men and women over 50 years of age, to quantitate its association with bone mineral density (BMD) and to assess its role as a confounder for the radiographic diagnosis of osteoporotic fracture. METHODS: In 27 centres participating in the population-based European Vertebral Osteoporosis Study (EVOS), standardised lateral radiographs of the lumbar and of the thoracic spine from T4 to L4 were assessed in all those of adequate quality. The presence of Scheuermann's disease, a confounder for prevalent fracture in later life, was defined by the presence of at least one Schmorl's node or irregular endplate together with kyphosis (sagittal Cobb angle >40° between T4 and T12) or a wedged-shaped vertebral body. Alternatively, the (rare) Edgren-Vaino sign was taken as diagnostic. The 6-point-per-vertebral-body (13 vertebrae) method was used to assess osteoporotic vertebral shape and fracture caseness. DXA BMD of the L2-L4 and femoral neck regions was measured in subsets. We also assessed the presence of Scheuermann's by alternative published algorithms when these used the radiographic signs we assessed. RESULTS: Vertebral radiographic images from 4486 men and 5655 women passed all quality checks. Prevalence of Scheuermann's varied considerably between centres, and based on random effect modelling, the overall European prevalence using our method was 8% with no significant difference between sexes. The highest prevalences were seen in Germany, Sweden, the UK and France and low prevalences were seen in Hungary, Poland and Slovakia. Centre-level prevalences in men and women were highly correlated. Scheuermann's was not associated with BMD of the spine or hip. CONCLUSIONS: Since most of the variation in population impact of Scheuermann's was unaccounted for by the radiological and anthropometric data, the search for new genetic and environmental determinants of this disease is encouraged.

Genetic taste blindness to bitter and body composition in childhood: a Mendelian randomization design.

BACKGROUND: The ability to taste 6-n-propylthiouracil (PROP) may be associated with body composition, but previous findings from observational studies are conflicting and cannot be interpreted causally. The aim of this study was to estimate the causal association between PROP taster status and body composition in a population-based cohort study. METHODS: The study was embedded in a population-based prospective birth cohort study. The TAS2R38 genotype (rs713598) was used as an instrumental variable (IV) to obtain unbiased effect estimates of the relation between PROP taster status and body weight (n=3778). Adiposity measures included body mass index (BMI) and fat mass measured by dual- energy X-ray absorptiometry scan at the child's age of 6 years. Associations were investigated using both ordinary linear regression (OLS) and two-stage least squares regression (2SLS). RESULTS: Non-taster girls had higher BMI standard deviation scores (SDS) and higher body fat as compared with taster girls (results from linear regression BMI SDS: -0.09, P=0.023, body fat mass (%): -0.49, P=0.028). The TAS2R38 genotype predicted PROP phenotype (F=240), indicating a strong IV. The 2SLS effect estimates were imprecise but similar to the observational estimates (-0.08 for BMI SDS and -0.46 for body fat mass %) and were not significantly different from the OLS results (Hausman test: P>0.10). For boys there were no differences observed between tasters and non-tasters. CONCLUSIONS: Our findings suggest a causal relation between PROP taster status and body weight among 6-year-old girls; Mendelian randomization was consistent with conventional estimates. In contrast, body weight among boys appeared to be independent of the PROP taster status. Further research should focus on possible underlying pathways, such as dietary behavior.

Drug-gene interactions and the search for missing heritability: a cross-sectional pharmacogenomics study of the QT interval.

Variability in response to drug use is common and heritable, suggesting that genome-wide pharmacogenomics studies may help explain the 'missing heritability' of complex traits. Here, we describe four independent analyses in 33 781 participants of European ancestry from 10 cohorts that were designed to identify genetic variants modifying the effects of drugs on QT interval duration (QT). Each analysis cross-sectionally examined four therapeutic classes: thiazide diuretics (prevalence of use=13.0%), tri/tetracyclic antidepressants (2.6%), sulfonylurea hypoglycemic agents (2.9%) and QT-prolonging drugs as classified by the University of Arizona Center for Education and Research on Therapeutics (4.4%). Drug-gene interactions were estimated using covariable-adjusted linear regression and results were combined with fixed-effects meta-analysis. Although drug-single-nucleotide polymorphism (SNP) interactions were biologically plausible and variables were well-measured, findings from the four cross-sectional meta-analyses were null (Pinteraction>5.0 × 10(-8)). Simulations suggested that additional efforts, including longitudinal modeling to increase statistical power, are likely needed to identify potentially important pharmacogenomic effects.

Prediction model for knee osteoarthritis incidence, including clinical, genetic and biochemical risk factors.

OBJECTIVE: To develop and validate a prognostic model for incident knee osteoarthritis (KOA) in a general population and determine the value of different risk factor groups to prediction. METHODS: The prognostic model was developed in 2628 individuals from the Rotterdam Study-I (RS-I). Univariate and multivariate analyses were performed for questionnaire/easily obtainable variables, imaging variables, genetic and biochemical markers. The extended multivariate model was tested on discrimination (receiver operating characteristic curve and area under the curve (AUC)) in two other population-based cohorts: Rotterdam Study-II and Chingford Study. RESULTS: In RS-I, there was moderate predictive value for incident KOA based on the genetic score alone in subjects aged <65 years (AUC 0.65), while it was only 0.55 for subjects aged ≥65 years. The AUC for gender, age and body mass index (BMI) in prediction for KOA was 0.66. Addition of the questionnaire variables, genetic score or biochemical marker urinary C-terminal cross-linked telopeptide of type II collagen to the model did not change the AUC. However, when adding the knee baseline KL score to the model the AUC increased to 0.79. Applying external validation, similar results were observed in the Rotterdam Study-II and the Chingford Study. CONCLUSIONS: Easy obtainable 'Questionnaire' variables, genetic markers, OA at other joint sites and biochemical markers add only modestly to the prediction of KOA incidence using age, gender and BMI in an elderly population. Doubtful minor radiographic degenerative features in the knee, however, are a very strong predictor of future KOA. This is an important finding, as many radiologists do not report minor degenerative changes in the knee.

Dissecting the relationship between high-sensitivity serum C-reactive protein and increased fracture risk: the Rotterdam Study.

UNLABELLED: Serum high-sensitivity C-reactive protein (CRP) is an inflammatory biomarker. We investigated the relationship between CRP and bone health in the Rotterdam Study. Serum high-sensitivity CRP was associated with fracture risk and lower femoral neck bending strength. Mendelian randomization analyses did not yield evidence for this relationship being causal. INTRODUCTION: Inflammatory diseases are associated with bone pathology, reflected in a higher fracture risk. Serum high-sensitivity CRP is an inflammatory biomarker. We investigated the relationship between CRP and bone mineral density (BMD), hip bone geometry, and incident fractures in the Rotterdam Study, a prospective population-based cohort. METHODS: At baseline, serum high-sensitivity CRP was measured. A weighted genetic risk score was compiled for CRP based on published studies (29 polymorphisms; Illumina HumanHap550 Beadchip genotyping and HapMap imputation). Regression models were reported per standard deviation increase in CRP adjusted for sex, age, and BMI. Complete data was available for 6,386 participants, of whom 1,561 persons sustained a fracture (mean follow-up, 11.6 years). RESULTS: CRP was associated with a risk for any type of fracture [hazard ratio (HR) = 1.06; 95 % confidence interval (CI), 1.02-1.11], hip fractures (HR = 1.09; 1.02-1.17) and vertebral fractures [odds ratio (OR) = 1.34; 1.14-1.58]. An inverse relationship between CRP levels and section modulus (-0.011 cm(3); -0.020 to -0.003 cm(3)) was observed. The combined genetic risk score of CRP single nucleotide polymorphisms (SNPs) was associated with serum CRP levels (p = 9 × 10(-56)), but not with fracture risk (HR = 1.00; 0.99-1.00; p = 0.23). CONCLUSIONS: Serum high-sensitivity CRP is associated with fracture risk and lower bending strength. Mendelian randomization analyses did not yield evidence for this relationship being causal. Future studies might reveal what factors truly underlie the relationship between CRP and fracture risk.

Lifestyle factors and metabolomic aging biomarkers: Meta-analysis of cross-sectional and longitudinal associations in three prospective cohorts.

Biological age uses biophysiological information to capture a person's age-related risk of adverse outcomes. MetaboAge and MetaboHealth are metabolomics-based biomarkers of biological age trained on chronological age and mortality risk, respectively. Lifestyle factors contribute to the extent chronological and biological age differ. The association of lifestyle factors with MetaboAge and MetaboHealth, potential sex differences in these associations, and MetaboAge's and MetaboHealth's sensitivity to lifestyle changes have not been studied yet. Linear regression analyses and mixed-effect models were used to examine the cross-sectional and longitudinal associations of scaled lifestyle factors with scaled MetaboAge and MetaboHealth in 24,332 middle-aged participants from the Doetinchem Cohort Study, Rotterdam Study, and UK Biobank. Random-effect meta-analyses were performed across cohorts. Repeated metabolomics measurements had a ten-year interval in the Doetinchem Cohort Study and a five-year interval in the UK Biobank. In the first study incorporating longitudinal information on MetaboAge and MetaboHealth, we demonstrate associations between current smoking, sleeping ≥8 hours/day, higher BMI, and larger waist circumference were associated with higher MetaboHealth, the latter two also with higher MetaboAge. Furthermore, adhering to the dietary and physical activity guidelines were inversely associated with MetaboHealth. Lastly, we observed sex differences in the associations between alcohol use and MetaboHealth.

The contribution of hip geometry to the prediction of hip osteoarthritis.

OBJECTIVE: To determine how well measures of hip geometry can predict radiological incident hip osteoarthritis (HOA) compared to well known clinical risk factors. DESIGN: The study population is part of the Rotterdam Study, a prospective population-based cohort. Baseline pelvic radiographs were used to measure hip geometry by two methods: Statistical Shape Models (SSM) and predefined geometry parameters (PGPs). Incident HOA (Kellgren and Lawrence (KL) ≥ 2) was assessed in 688 participants after 6.5 years without radiographic HOA at baseline. The ability to predict HOA was quantified using the area under the Receiver Operating Characteristics (ROC) curve (AUC). RESULTS: Comparison of the two methods showed that both contain information that is not captured by the other method. At 6.5 years follow-up 132 hips had incident HOA. Five PGPs (Wiberg angle, Neck Width (NW), Pelvic Width (PW), Hip Axis Length (HAL) and Triangular Index (TI)) and two SSM (modes 5 and 9) were significant predictors of HOA (P = 0.007). Hip geometry added 7% to the prediction obtained by clinical risk factors (AUC = 0.67 (geometry), 0.66 (gender, age, Body Mass Index (BMI)) and combining both: AUC = 0.73, respectively). Mode 12 (associated with position of the femoral head in acetabulum) and Wiberg angle were predictors of HOA in participants without radiological signs at baseline (KL = 0). Although the strength of the prediction decreased for all variables at a longer follow-up, the contribution of hip geometry was still significant (P = 0.01). CONCLUSIONS: Hip geometry has a moderate ability to predict HOA in participants with and without initial signs of osteoarthritis (OA), similar to and largely independent of the predictive value of clinical risk factors.

Genome-wide association analysis of coffee drinking suggests association with CYP1A1/CYP1A2 and NRCAM.

Coffee consumption is a model for addictive behavior. We performed a meta-analysis of genome-wide association studies (GWASs) on coffee intake from 8 Caucasian cohorts (N=18 176) and sought replication of our top findings in a further 7929 individuals. We also performed a gene expression analysis treating different cell lines with caffeine. Genome-wide significant association was observed for two single-nucleotide polymorphisms (SNPs) in the 15q24 region. The two SNPs rs2470893 and rs2472297 (P-values=1.6 × 10(-11) and 2.7 × 10(-11)), which were also in strong linkage disequilibrium (r(2)=0.7) with each other, lie in the 23-kb long commonly shared 5' flanking region between CYP1A1 and CYP1A2 genes. CYP1A1 was found to be downregulated in lymphoblastoid cell lines treated with caffeine. CYP1A1 is known to metabolize polycyclic aromatic hydrocarbons, which are important constituents of coffee, whereas CYP1A2 is involved in the primary metabolism of caffeine. Significant evidence of association was also detected at rs382140 (P-value=3.9 × 10(-09)) near NRCAM-a gene implicated in vulnerability to addiction, and at another independent hit rs6495122 (P-value=7.1 × 10(-09))-an SNP associated with blood pressure-in the 15q24 region near the gene ULK3, in the meta-analysis of discovery and replication cohorts. Our results from GWASs and expression analysis also strongly implicate CAB39L in coffee drinking. Pathway analysis of differentially expressed genes revealed significantly enriched ubiquitin proteasome (P-value=2.2 × 10(-05)) and Parkinson's disease pathways (P-value=3.6 × 10(-05)).

Gremlin 1, frizzled-related protein, and Dkk-1 are key regulators of human articular cartilage homeostasis.

OBJECTIVE: The development of osteoarthritis (OA) may be caused by activation of hypertrophic differentiation of articular chondrocytes. Healthy articular cartilage is highly resistant to hypertrophic differentiation, in contrast to other hyaline cartilage subtypes, such as growth plate cartilage. The purpose of this study was to elucidate the molecular mechanism responsible for the difference in the propensity of human articular cartilage and growth plate cartilage to undergo hypertrophic differentiation. METHODS: Whole-genome gene-expression microarray analysis of healthy human growth plate and articular cartilage derived from the same adolescent donors was performed. Candidate genes, which were enriched in the articular cartilage, were validated at the messenger RNA (mRNA) and protein levels and examined for their potential to inhibit hypertrophic differentiation in two models. In addition, we studied a possible genetic association with OA. RESULTS: Pathway analysis demonstrated decreased Wnt signaling in articular cartilage as compared to growth plate cartilage. This was at least partly due to increased expression of the bone morphogenetic protein and Wnt antagonists Gremlin 1, Frizzled-related protein (FRP), and Dkk-1 at the mRNA and protein levels in articular cartilage. Supplementation of these proteins diminished terminal hypertrophic differentiation without affecting chondrogenesis in long-bone explant cultures and chondrogenically differentiating human mesenchymal stem cells. Additionally, we found that single-nucleotide polymorphism rs12593365, which is located in a genomic control region of GREM1, was significantly associated with a 20% reduced risk of radiographic hip OA in 2 population-based cohorts. CONCLUSION: Taken together, our study identified Gremlin 1, FRP, and Dkk-1 as natural brakes on hypertrophic differentiation in articular cartilage. As hypertrophic differentiation of articular cartilage may contribute to the development of OA, our findings may open new avenues for therapeutic intervention.

The effect of thiazide and loop diuretics on urinary levels of free deoxypyridinoline: an osteoclastic bone-resorption marker.

WHAT IS KNOWN AND OBJECTIVE: Diuretics can cause changes in calcium levels due to renal effects. Moreover, calcium levels can also vary as a result of changes in intestinal absorption and in the activity of osteoclastic cells. A marker of osteoclastic bone-resorption activity is the level of urinary free deoxypyridinoline (FDP). Deoxypyridinoline (DP) acts as a cross-link between adjacent collagen chains to provide structural rigidity. Our aim was to investigate the association between use of thiazides and loop diuretics and urinary levels FDP. METHODS: In this follow-up study, data were obtained from the Rotterdam Study, a large population-based prospective cohort study. For a subset of 658 participants, urinary levels of FDP were measured at baseline. Linear regression analysis was performed to assess the association between the use of thiazides and loop diuretics and the urinary levels of FDP. RESULTS: In women, current use of loop diuretics for less than 42 days was associated with an increased level of urinary FDP (+3·43 nmol deoxypyridinoline per mmol urinary creatinine; 95% CI 1·85; 5·02) compared with no use. However, use for a period of more than 42 days was not associated with an increased level of FDP, nor was past use of loop diuretics. For thiazide diuretics, no statistically significant associations were found. WHAT IS NEW AND CONCLUSION: In women, short-term use of loop diuretics is associated with an increased level of FDP, reflecting increased bone resorption by osteoclasts. As the difference disappears with longer term use, the clinical significance is unclear and the value of FDP as a biomarker in this setting is not established. The molecular mechanism for the observed differences in bone fracture rates with use of diuretics remains unclear.

Novel Genetic Determinants of Dental Maturation in Children.

Dental occlusion requires harmonious development of teeth, jaws, and other elements of the craniofacial complex, which are regulated by environmental and genetic factors. We performed the first genome-wide association study (GWAS) on dental development (DD) using the Demirjian radiographic method. Radiographic assessments from participants of the Generation R Study (primary study population, N1 = 2,793; mean age of 9.8 y) were correlated with ~30 million genetic variants while adjusting for age, sex, and genomic principal components (proxy for population stratification). Variants associated with DD at genome-wide significant level (P < 5 × 10-8) mapped to 16q12.2 (IRX5) (lead variant rs3922616, B = 0.16; P = 2.2 × 10-8). We used Fisher's combined probability tests weighted by sample size to perform a meta-analysis (N = 14,805) combining radiographic DD at a mean age of 9.8 y from Generation R with data from a previous GWAS (N2 = 12,012) on number of teeth (NT) in infants used as proxy of DD at a mean age of 9.8 y (including the ALSPAC and NFBC1966). This GWAS meta-analysis revealed 3 novel loci mapping to 7p15.3 (IGF2BP3: P = 3.2 × 10-8), 14q13.3 (PAX9: P = 1.9 × 10-8), and 16q12.2 (IRX5: P = 1.2 × 10-9) and validated 8 previously reported NT loci. A polygenic allele score constructed from these 11 loci was associated with radiographic DD in an independent Generation R set of children (N = 703; B = 0.05, P = 0.004). Furthermore, profiling of the identified genes across an atlas of murine and human stem cells observed expression in the cells involved in the formation of bone and/or dental tissues (>0.3 frequency per kilobase of transcript per million mapped reads), likely reflecting functional specialization. Our findings provide biological insight into the polygenic architecture of the pediatric dental maturation process.