Impact of Respiratory and Cardiovascular Drug Exposure on Lung Function Trajectories.

Rationale: Research suggests that respiratory and cardiovascular drugs can ameliorate the rate of lung function decline. Objectives: To investigate the impact of respiratory and cardiovascular pharmacotherapy on lung function trajectories in the general population. Methods: Repeated spirometry was performed in the Rotterdam Study, a population-based cohort of adults aged ≥45 years. Exposure to long-acting beta2-agonists (LABA), long-acting muscarinic antagonists (LAMA), inhaled corticosteroids (ICS), cardioselective beta-blockers, calcium channel blockers, ACE-inhibitors, angiotensin-II receptor blockers, and statins was quantified from pharmacy records to account for therapy adherence. Propensity-score matching and multinomial logistic regression were performed to model medication effects on lung function trajectories, which were previously identified based on FEV1 and FVC patterns. Models were additionally stratified by genetic variation in each drug target. Measurement and Main Results: Among 3,783 individuals, 2,974 (78.6%) were classified as normal lung function decliners, 432 (11.4%) as rapid decliners, and 377 (10.0%) as improvers. Exposure to LABA (odds ratio (OR) =1.09 [95%-CI: 1.03-1.16] per 10% increase in exposure), ICS (OR=1.08 [95%-CI: 1.02-1.14]), and statins (OR=1.04 [95%-CI: 1.02-1.06]) significantly increased the odds of being an improver compared to a normal decliner. Beta1-blocker use was associated with higher odds of being a rapid decliner (OR=1.04 [95%-CI: 1.00-1.09]), which was driven by incident users. Pharmacogenetic analysis suggests that the effects of LABA, ICS, and beta1-blockers are dependent on genetic variation in their drug targets. Conclusions: Our study suggests that LABA, ICS, and statins may favorably modulate lung function trajectories in adults, while initiation of beta1-blockers was associated with rapid lung function decline.

Interindividual variation in ovarian reserve after gonadotoxic treatment in female childhood cancer survivors - a genome-wide association study: results from PanCareLIFE.

OBJECTIVE: To discover new variants associated with low ovarian reserve after gonadotoxic treatment among adult female childhood cancer survivors using a genome-wide association study approach. DESIGN: Genome-wide association study. SETTING: Not applicable. PATIENTS: A discovery cohort of adult female childhood cancer survivors from the pan-European PanCareLIFE cohort (n = 743; median age: 25.8 years), excluding those who received bilateral ovarian irradiation, bilateral oophorectomy, central nervous system or total body irradiation, or stem cell transplantation. Replication was attempted in the US-based St. Jude Lifetime Cohort (n = 391; median age: 31.3 years). EXPOSURE: Female childhood cancer survivors are at risk of therapy-related gonadal impairment. Alkylating agents are well-established risk factors, and the interindividual variability in gonadotoxicity may be explained by genetic polymorphisms. Data were collected in real-life conditions, and cyclophosphamide equivalent doses were used to quantify alkylation agent exposure. MAIN OUTCOME MEASURE: Anti-Müllerian hormone (AMH) levels served as a proxy for ovarian function, and the findings were combined in a meta-analysis. RESULTS: Three genome-wide significant (<5.0 × 10-8) and 16 genome-wide suggestive (<5.0 × 10-6) loci were associated with log-transformed AMH levels, adjusted for cyclophosphamide equivalent dose of alkylating agents, age at diagnosis, and age at study in the PanCareLIFE cohort. On the basis of the effect allele frequency (EAF) (>0.01 if not genome-wide significant), and biologic relevance, 15 single nucleotide polymorphisms were selected for replication. None of the single nucleotide polymorphisms were statistically significantly associated with AMH levels. A meta-analysis indicated that rs78861946 was associated with borderline genome-wide statistical significance (reference/effect allele: C/T; effect allele frequency: 0.04, beta (SE): -0.484 (0.091). CONCLUSION: This study found no genetic variants associated with a lower ovarian reserve after gonadotoxic treatment because the findings of this genome-wide association study were not statistically significant replicated in the replication cohort. Suggestive evidence for the potential importance of 1 variant is briefly discussed, but the lack of statistical significance calls for larger cohort sizes. Because the population of childhood cancer survivors is increasing, large-scale and systematic research is needed to identify genetic variants that could aid predictive risk models of gonadotoxicity as well as fertility preservation options for childhood cancer survivors.

Bone health index in the assessment of bone health: The Generation R Study.

Bone Health Index (BHI) has been proposed as a useful instrument for assessing bone health in children. However, its relationship with fracture risk remains unknown. We aimed to investigate whether BHI is associated with bone mineral density (BMD) and prevalent fracture odds in children from the Generation R Study. We also implemented genome-wide association study (GWAS) and polygenic score (PGS) approaches to improve our understanding of BHI and its potential. In total, 4150 children (49.4 % boys; aged 9.8 years) with genotyped data and bone assessments were included in this study. BMD was measured across the total body (less head following ISCD guidelines) using a GE-Lunar iDXA densitometer; and BHI was determined from the hand DXA scans using BoneXpert®. Fractures were self-reported collected with home questionnaires. The association of BHI with BMD and fractures was evaluated using linear models corrected for age, sex, ethnicity, height, and weight. We observed a positive correlation between BHI and BMD (ρ = 0.32, p-value<0.0001). Further, every SD decrease in BHI was associated with an 11 % increased risk of prevalent fractures (OR:1.11, 95 % CI 1.00-1.24, p-value = 0.05). Our BHI GWAS identified variants (lead SNP rs1404264-A, p-value = 2.61 × 10-14) mapping to the ING3/CPED1/WNT16 locus. Children in the extreme tails of the BMD PGS presented a difference in BHI values of -0.10 standard deviations (95% CI -0.14 to -0.07; p-value<0.0001). On top of the demonstrated epidemiological association of BHI with both BMD and fracture risk, our results reveal a partially shared biological background between BHI and BMD. These findings highlight the potential value of using BHI to screen children at risk of fracture.

Disruption of TUFT1, a Desmosome-Associated Protein, Causes Skin Fragility, Woolly Hair, and Palmoplantar Keratoderma.

Desmosomes are dynamic complex protein structures involved in cellular adhesion. Disruption of these structures by loss-of-function variants in desmosomal genes leads to a variety of skin- and heart-related phenotypes. In this study, we report TUFT1 as a desmosome-associated protein, implicated in epidermal integrity. In two siblings with mild skin fragility, woolly hair, and mild palmoplantar keratoderma but without a cardiac phenotype, we identified a homozygous splice-site variant in the TUFT1 gene, leading to aberrant mRNA splicing and loss of TUFT1 protein. Patients' skin and keratinocytes showed acantholysis, perinuclear retraction of intermediate filaments, and reduced mechanical stress resistance. Immunolabeling and transfection studies showed that TUFT1 is positioned within the desmosome and that its location is dependent on the presence of the desmoplakin carboxy-terminal tail. A Tuft1-knockout mouse model mimicked the patients' phenotypes. Altogether, this study reveals TUFT1 as a desmosome-associated protein, whose absence causes skin fragility, woolly hair, and palmoplantar keratoderma.

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.

Low-Cost High-Throughput Genotyping for Diagnosing Familial Hypercholesterolemia.

BACKGROUND: Familial hypercholesterolemia (FH) is a common but underdiagnosed genetic disorder characterized by high low-density lipoprotein cholesterol levels and premature cardiovascular disease. Current sequencing methods to diagnose FH are expensive and time-consuming. In this study, we evaluated the accuracy of a low-cost, high-throughput genotyping array for diagnosing FH. METHODS: An Illumina Global Screening Array was customized to include probes for 636 variants, previously classified as FH-causing variants. First, its theoretical coverage was assessed in all FH variant carriers diagnosed through next-generation sequencing between 2016 and 2022 in the Netherlands (n=1772). Next, the performance of the array was validated in another sample of FH variant carriers previously identified in the Dutch FH cascade screening program (n=1268). RESULTS: The theoretical coverage of the array for FH-causing variants was 91.3%. Validation of the array was assessed in a sample of 1268 carriers of whom 1015 carried a variant in LDLR, 250 in APOB, and 3 in PCSK9. The overall sensitivity was 94.7% and increased to 98.2% after excluding participants with variants not included in the array design. Copy number variation analysis yielded a 89.4% sensitivity. In 18 carriers, the array identified a total of 19 additional FH-causing variants. Subsequent DNA analysis confirmed 5 of the additionally identified variants, yielding a false-positive result in 16 subjects (1.3%). CONCLUSIONS: The FH genotyping array is a promising tool for genetically diagnosing FH at low costs and has the potential to greatly increase accessibility to genetic testing for FH. Continuous customization of the array will further improve its performance.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification.

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

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.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

Advanced Glycation End Products (AGEs) in Diet and Skin in Relation to Stool Microbiota: The Rotterdam Study.

BACKGROUND: Advanced glycation end products (AGEs) are involved in age-related diseases, but the interaction of gut microbiota with dietary AGEs (dAGEs) and tissue AGEs in the population is unknown. OBJECTIVE: Our objective was to investigate the association of dietary and tissue AGEs with gut microbiota in the population-based Rotterdam Study, using skin AGEs as a marker for tissue accumulation and stool microbiota as a surrogate for gut microbiota. DESIGN: Dietary intake of three AGEs (dAGEs), namely carboxymethyl-lysine (CML), N-(5-hydro-5-methyl-4-imidazolon-2-yl)-ornithine (MGH1), and carboxyethyl-lysine (CEL), was quantified at baseline from food frequency questionnaires. Following up after a median of 5.7 years, skin AGEs were measured using skin autofluorescence (SAF), and stool microbiota samples were sequenced (16S rRNA) to measure microbial composition (including alpha-diversity, beta-dissimilarity, and taxonomic abundances) as well as predict microbial metabolic pathways. Associations of both dAGEs and SAF with microbial measures were investigated using multiple linear regression models in 1052 and 718 participants, respectively. RESULTS: dAGEs and SAF were not associated with either the alpha-diversity or beta-dissimilarity of the stool microbiota. After multiple-testing correction, dAGEs were not associated with any of the 188 genera tested, but were nominally inversely associated with the abundance of Barnesiella, Colidextribacter, Oscillospiraceae UCG-005, and Terrisporobacter, in addition to being positively associated with Coprococcus, Dorea, and Blautia. A higher abundance of Lactobacillus was associated with a higher SAF, along with several nominally significantly associated genera. dAGEs and SAF were nominally associated with several microbial pathways, but none were statistically significant after multiple-testing correction. CONCLUSIONS: Our findings did not solidify a link between habitual dAGEs, skin AGEs, and overall stool microbiota composition. Nominally significant associations with several genera and functional pathways suggested a potential interaction between gut microbiota and AGE metabolism, but validation is required. Future studies are warranted, to investigate whether gut microbiota modifies the potential impact of dAGEs on health.

Sexual dimorphisms in serum calcium and phosphate concentrations in the Rotterdam Study.

Sex differences in serum phosphate and calcium have been reported but the exact nature and underlying regulatory mechanisms remain unclear. We aimed to compare calcium and phosphate concentrations between sexes, and explore potential covariates to elucidate underlying mechanisms of sex differences in a prospective, population-based cohort study. Pooled data of subjects > 45 years from three independent cohorts of the Rotterdam Study (RS) were used: RS-I-3 (n = 3623), RS-II-1 (n = 2394), RS-III-1 (n = 3241), with separate analyses from an additional time point of the first cohort RS-I-1 (n = 2688). Compared to men, women had significantly higher total serum calcium and phosphate concentrations which was not explained by BMI, kidney function nor smoking. Adjustment for serum estradiol diminished sex differences in serum calcium while adjustment for serum testosterone diminished sex differences in serum phosphate. Adjustment for vitamin D and alkaline phosphatase did not change the association between sex and calcium or phosphate in RS-I-1. In the sex-combined group, both serum calcium and phosphate decreased with age with a significant interaction for sex differences for serum calcium but not phosphate. In sex-stratified analyses, serum estradiol but not testosterone was inversely associated with serum calcium in both sexes. Serum estradiol was inversely associated with serum phosphate in both sexes to a similar degree, while serum testosterone was inversely associated with serum phosphate in both sexes with an apparent stronger effect in men than in women. Premenopausal women had lower serum phosphate compared to postmenopausal women. Serum testosterone was inversely associated with serum phosphate in postmenopausal women only. In conclusion, women > 45 years have higher serum calcium and phosphate concentrations compared to men of similar age, not explained by vitamin D or alkaline phosphatase concentrations. Serum estradiol but not testosterone was inversely associated with serum calcium while serum testosterone was inversely associated with serum phosphate in both sexes. Serum testosterone may in part explain sex differences in serum phosphate while estradiol could partly explain sex differences in serum calcium.

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.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk.

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Lifestyle Intervention Randomized Controlled Trial for Age-Related Macular Degeneration (AMD-Life): Study Design.

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

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.

Combining genome-wide association studies highlight novel loci involved in human facial variation.

Standard genome-wide association studies (GWASs) rely on analyzing a single trait at a time. However, many human phenotypes are complex and composed by multiple correlated traits. Here we introduce C-GWAS, a method for combining GWAS summary statistics of multiple potentially correlated traits. Extensive computer simulations demonstrated increased statistical power of C-GWAS compared to the minimal p-values of multiple single-trait GWASs (MinGWAS) and the current state-of-the-art method for combining single-trait GWASs (MTAG). Applying C-GWAS to a meta-analysis dataset of 78 single trait facial GWASs from 10,115 Europeans identified 56 study-wide suggestively significant loci with multi-trait effects on facial morphology of which 17 are novel loci. Using data from additional 13,622 European and Asian samples, 46 (82%) loci, including 9 (53%) novel loci, were replicated at nominal significance with consistent allele effects. Functional analyses further strengthen the reliability of our C-GWAS findings. Our study introduces the C-GWAS method and makes it available as computationally efficient open-source R package for widespread future use. Our work also provides insights into the genetic architecture of human facial appearance.

Gut microbiome-wide association study of depressive symptoms.

Depression is one of the most poorly understood diseases due to its elusive pathogenesis. There is an urgency to identify molecular and biological mechanisms underlying depression and the gut microbiome is a novel area of interest. Here we investigate the relation of fecal microbiome diversity and composition with depressive symptoms in 1,054 participants from the Rotterdam Study cohort and validate these findings in the Amsterdam HELIUS cohort in 1,539 subjects. We identify association of thirteen microbial taxa, including genera Eggerthella, Subdoligranulum, Coprococcus, Sellimonas, Lachnoclostridium, Hungatella, Ruminococcaceae (UCG002, UCG003 and UCG005), LachnospiraceaeUCG001, Eubacterium ventriosum and Ruminococcusgauvreauiigroup, and family Ruminococcaceae with depressive symptoms. These bacteria are known to be involved in the synthesis of glutamate, butyrate, serotonin and gamma amino butyric acid (GABA), which are key neurotransmitters for depression. Our study suggests that the gut microbiome composition may play a key role in depression.

Where is bone science taking us?

Bone science has over the last decades unraveled many important pathways in bone and mineral metabolism and the interplay between genetic factors and the environment. Some of these discoveries have led to the development of pharmacological treatments of osteoporosis and rare bone diseases. Other scientific avenues have uncovered a role for the gut microbiome in regulating bone mass, which have led to investigations on the possible therapeutic role of probiotics in the prevention of osteoporosis. Huge advances have been made in identifying the genes that cause rare bone diseases, which in some cases have led to therapeutic interventions. Advances have also been made in understanding the genetic basis of the more common polygenic bone diseases, including osteoporosis and Paget's disease of bone (PDB). Polygenic profiles are used for establishing genetic risk scores aiming at early diagnosis and intervention, but also in Mendelian randomization (MR) studies to investigate both desired and undesired effects of targets for drug design.