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.

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.

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.

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.

A saturated map of common genetic variants associated with human height

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

OUP accepted manuscript

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⁻⁸, P = 1.35 × 10⁻³ and P = 1.24 × 10⁻³, respectively). This result was validated in a subsequent meta-analysis (P = 3.24 × 10⁻¹²). 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⁻⁵). 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.

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

Common and rare variant association analyses in amyotrophic lateral sclerosis identify 15 risk loci with distinct genetic architectures and neuron-specific biology

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with a lifetime risk of one in 350 people and an unmet need for disease-modifying therapies. We conducted a cross-ancestry genome-wide association study (GWAS) including 29,612 patients with ALS and 122,656 controls, which identified 15 risk loci. When combined with 8,953 individuals with whole-genome sequencing (6,538 patients, 2,415 controls) and a large cortex-derived expression quantitative trait locus (eQTL) dataset (MetaBrain), analyses revealed locus-specific genetic architectures in which we prioritized genes either through rare variants, short tandem repeats or regulatory effects. ALS-associated risk loci were shared with multiple traits within the neurodegenerative spectrum but with distinct enrichment patterns across brain regions and cell types. Of the environmental and lifestyle risk factors obtained from the literature, Mendelian randomization analyses indicated a causal role for high cholesterol levels. The combination of all ALS-associated signals reveals a role for perturbations in vesicle-mediated transport and autophagy and provides evidence for cell-autonomous disease initiation in glutamatergic neurons.

The genomics of heart failure: design and rationale of the HERMES consortium

Aims

The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure.

Methods and results

The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome‐wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow‐up following heart failure diagnosis ranged from 2 to 116 months. Forty‐nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34–90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low‐frequency variants (allele frequency 0.01–0.05) at P < 5 × 10⁻⁸ under an additive genetic model.

Conclusions

HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction.

Sugar-Sweetened Beverage Consumption May Modify Associations Between Genetic Variants in the CHREBP (Carbohydrate Responsive Element Binding Protein) Locus and HDL-C (High-Density Lipoprotein Cholesterol) and Triglyceride Concentrations

Supplemental Digital Content is available in the text.

Background:

ChREBP (carbohydrate responsive element binding protein) is a transcription factor that responds to sugar consumption. Sugar-sweetened beverage (SSB) consumption and genetic variants in the CHREBP locus have separately been linked to HDL-C (high-density lipoprotein cholesterol) and triglyceride concentrations. We hypothesized that SSB consumption would modify the association between genetic variants in the CHREBP locus and dyslipidemia.

Methods:

Data from 11 cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (N=63 599) and the UK Biobank (N=59 220) were used to quantify associations of SSB consumption, genetic variants, and their interaction on HDL-C and triglyceride concentrations using linear regression models. A total of 1606 single nucleotide polymorphisms within or near CHREBP were considered. SSB consumption was estimated from validated questionnaires, and participants were grouped by their estimated intake.

Results:

In a meta-analysis, rs71556729 was significantly associated with higher HDL-C concentrations only among the highest SSB consumers (β, 2.12 [95% CI, 1.16–3.07] mg/dL per allele; P<0.0001), but not significantly among the lowest SSB consumers (P=0.81; P_Diff <0.0001). Similar results were observed for 2 additional variants (rs35709627 and rs71556736). For triglyceride, rs55673514 was positively associated with triglyceride concentrations only among the highest SSB consumers (β, 0.06 [95% CI, 0.02–0.09] ln-mg/dL per allele, P=0.001) but not the lowest SSB consumers (P=0.84; P_Diff=0.0005).

Conclusions:

Our results identified genetic variants in the CHREBP locus that may protect against SSB-associated reductions in HDL-C and other variants that may exacerbate SSB-associated increases in triglyceride concentrations.

Registration:

URL: https://www.clinicaltrials.gov; Unique identifier: NCT00005133, NCT00005121, NCT00005487, and NCT00000479.

TCERG1L allelic variation is associated with cisplatin-induced hearing loss in childhood cancer, a PanCareLIFE study

In children with cancer, the heterogeneity in ototoxicity occurrence after similar treatment suggests a role for genetic susceptibility. Using a genome-wide association study (GWAS) approach, we identified a genetic variant in TCERG1L (rs893507) to be associated with hearing loss in 390 non-cranial irradiated, cisplatin-treated children with cancer. These results were replicated in two independent, similarly treated cohorts (n = 192 and 188, respectively) (combined cohort: P = 5.3 × 10-10, OR 3.11, 95% CI 2.2-4.5). Modulating TCERG1L expression in cultured human cells revealed significantly altered cellular responses to cisplatin-induced cytokine secretion and toxicity. These results contribute to insights into the genetic and pathophysiological basis of cisplatin-induced ototoxicity.

Association of Insulin Resistance and Type 2 Diabetes With Gut Microbial Diversity: A Microbiome-Wide Analysis From Population Studies

IMPORTANCE

Previous studies have indicated that gut microbiome may be associated with development of type 2 diabetes. However, these studies are limited by small sample size and insufficient for confounding. Furthermore, which specific taxa play a role in the development of type 2 diabetes remains unclear.

OBJECTIVE

To examine associations of gut microbiome composition with insulin resistance and type 2 diabetes in a large population-based setting controlling for various sociodemographic and lifestyle factors.

DESIGN, SETTING, AND PARTICIPANTS

This cross-sectional analysis included 2166 participants from 2 Dutch population-based prospective cohorts: the Rotterdam Study and the LifeLines-DEEP study.

EXPOSURES

The 16S ribosomal RNA method was used to measure microbiome composition in stool samples collected between January 1, 2012, and December 31, 2013. The α diversity (Shannon, richness, and Inverse Simpson indexes), β diversity (Bray-Curtis dissimilarity matrix), and taxa (from domain to genus level) were identified to reflect gut microbiome composition.

MAIN OUTCOMES AND MEASURES

Associations among α diversity, β diversity, and taxa with the Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) and with type 2 diabetes were examined. Glucose and insulin were measured to calculate the HOMA-IR. Type 2 diabetes cases were identified based on glucose levels and medical records from January 2012 to December 2013. Analyses were adjusted for technical covariates, lifestyle, sociodemographic, and medical factors. Data analysis was performed from January 1, 2018, to December 31, 2020.

RESULTS

There were 2166 participants in this study: 1418 from the Rotterdam Study (mean [SD] age, 62.4 [5.9] years; 815 [57.5%] male) and 748 from the LifeLines-DEEP study (mean [SD] age, 44.7 [13.4] years; 431 [57.6%] male); from this total, 193 type 2 diabetes cases were identified. Lower microbiome Shannon index and richness were associated with higher HOMA-IR (eg, Shannon index, -0.06; 95% CI, -0.10 to -0.02), and patients with type 2 diabetes had a lower richness than participants without diabetes (odds ratio [OR], 0.93; 95% CI, 0.88-0.99). The β diversity (Bray-Curtis dissimilarity matrix) was associated with insulin resistance (R2 = 0.004, P = .001 in the Rotterdam Study and R2 = 0.005, P = .002 in the LifeLines-DEEP study). A total of 12 groups of bacteria were associated with HOMA-IR or type 2 diabetes. Specifically, a higher abundance of Christensenellaceae (β = -0.08; 95% CI, -0.12 to -0.03: P < .001), Christensenellaceae R7 group (β = -0.07; 95% CI, -0.12 to -0.03; P < .001), Marvinbryantia (β = -0.07; 95% CI, -0.11 to -0.03; P < .001), Ruminococcaceae UCG005 (β = -0.09; 95% CI, -0.13 to -0.05; P < .001), Ruminococcaceae UCG008 (β = -0.07; 95% CI, -0.11 to -0.03; P < .001), Ruminococcaceae UCG010 (β = -0.08; 95% CI, -0.12 to -0.04; P < .001), or Ruminococcaceae NK4A214 group (β = -0.09; 95% CI, -0.13 to -0.05; P < .001) was associated with lower HOMA-IR. A higher abundance of Clostridiaceae 1 (OR, 0.51; 95% CI, 0.41-0.65; P < .001), Peptostreptococcaceae (OR, 0.56; 95% CI, 0.45-0.70; P < .001), C sensu stricto 1 (OR, 0.51; 95% CI, 0.40-0.65; P < .001), Intestinibacter (OR, 0.60; 95% CI, 0.48-0.76; P < .001), or Romboutsia (OR, 0.55; 95% CI, 0.44-0.70; P < .001) was associated with less type 2 diabetes. These bacteria are all known to produce butyrate.

CONCLUSIONS AND RELEVANCE

In this cross-sectional study, higher microbiome α diversity, along with more butyrate-producing gut bacteria, was associated with less type 2 diabetes and with lower insulin resistance among individuals without diabetes. These findings could help provide insight into the etiology, pathogenesis, and treatment of type 2 diabetes.

Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging

BACKGROUND

Biological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field.

RESULTS

Leveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels.

CONCLUSION

This study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

Prevalence of and risk factors for extended-spectrum beta-lactamase genes carriership in a population-based cohort of middle-aged and elderly

INTRODUCTION

Increasing resistance to beta-lactam antibiotics is an alarming development worldwide. Fecal carriership of TEM, SHV, CTX-M and CMY was studied in a community-dwelling population of middle-aged and elderly individuals.

PATIENTS AND METHODS

Feces was obtained from individuals of the Rotterdam Study. Carriership of the TEM, SHV, CTX-M and CMY genes was determined using real-time polymerase chain reaction (qPCR). Possible associations were investigated between carriership of these genes and several risk factors, such as the use of antimicrobial drugs, diabetes mellitus, protein pump inhibitor (PPI) use, travelling, the composition of the gut microbiota, and intake of certain foods.

RESULTS

The most prevalent gene was TEM (53.0%), followed by SHV (18.4%), CTX-M (5.4%) and CMY (3.6%). Use of penicillins with extended spectrum was associated with TEM carriership, whereas use of macrolides and lincosamides was associated with TEM and SHV carriership. Interestingly, use of PPIs was associated with a higher prevalence of carriership of TEM, SHV and CMY (TEM: odds ratio [OR] 1.34; 95% confidence interval [CI] 1.05-1.77; SHV: OR 2.17; 95%CI 1.55-2.87; CMY: OR 2.26; 95%CI 1.23-4.11). Furthermore, associations were found between the richness and composition of the gut microbiota and TEM and SHV carriership.

CONCLUSIONS

The prevalence of carriership of TEM was substantial, but the prevalence of carriership of the extended-spectrum β-lactamase gene, CTX-M and the AmpC β-lactamase gene, CMY was relatively low in this community-dwelling, population-based cohort. The composition of the microbiota might play a role in the retention of resistance genes, but future studies are necessary to further elucidate this relationship.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

Age-dependent sex differences in calcium and phosphate homeostasis

BACKGROUND

Sex differences in calcium and phosphate have been observed. We aimed to assess a relation with age.

METHODS

We used the laboratory values of serum calcium, phosphate and albumin from three different samples ( 2005, 2010 and 2014 years) using the hospital information system of Erasmus MC, Rotterdam. The samples were divided into three age groups: 1-17, 18-44 and ≥45 years. Sex differences in calcium and phosphate were analyzed using ANCOVA, adjusting for age and serum albumin. Furthermore, sex by age interactions were determined and we analyzed differences between age groups stratified by sex.

RESULTS

In all three samples there was a significant sex × age interaction for serum calcium and phosphate, whose levels were significantly higher in women compared to men above 45 years. No sex differences in the younger age groups were found. In men, serum calcium and phosphate levels were highest in the youngest age group compared to age groups of 18-44 and ≥45 years. In women, serum calcium levels were significantly higher in the age group 1-17 and the age group ≥45 years compared to the 18-44 years age group. In women, serum phosphate was different between the three different age groups with highest level in the group 1-17 years and lowest in the group 18-44 years.

CONCLUSION

There are age- dependent sex differences in serum calcium and phosphate. Furthermore, we found differences in serum calcium and phosphate between different age groups. Underlying mechanisms for these age- and sex- differences are not yet fully elucidated.

Large-scale association analyses identify host factors influencing human gut microbiome composition

To study the effect of host genetics on gut microbiome composition, the MiBioGen consortium curated and analyzed genome-wide genotypes and 16S fecal microbiome data from 18,340 individuals (24 cohorts). Microbial composition showed high variability across cohorts: only 9 of 410 genera were detected in more than 95% of samples. A genome-wide association study of host genetic variation regarding microbial taxa identified 31 loci affecting the microbiome at a genome-wide significant (P < 5 × 10-8) threshold. One locus, the lactase (LCT) gene locus, reached study-wide significance (genome-wide association study signal: P = 1.28 × 10-20), and it showed an age-dependent association with Bifidobacterium abundance. Other associations were suggestive (1.95 × 10-10 < P < 5 × 10-8) but enriched for taxa showing high heritability and for genes expressed in the intestine and brain. A phenome-wide association study and Mendelian randomization identified enrichment of microbiome trait loci in the metabolic, nutrition and environment domains and suggested the microbiome might have causal effects in ulcerative colitis and rheumatoid arthritis.

Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women

Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256,523 Europeans aged 60 years and over from 22 cohorts we identify 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n = 48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (HLA-DQA1 p = 4 × 10-17), arthritis (GDF5 p = 4 × 10-13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, haematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing.

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 (r_g ≈ 0.15-0.5). In contrast, relative carbohydrate and sugar intake have negative genetic correlations with waist circumference, waist-hip ratio, and neighborhood deprivation (|r_g| ≈ 0.1-0.3) and positive genetic correlations with physical activity (r_g ≈ 0.1 and 0.2). Relative fat intake has no consistent pattern of genetic correlations with poor health but has a negative genetic correlation with educational attainment (r_g ≈-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.

Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity

Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10^-16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity.

Genetic basis of falling risk susceptibility in the UK Biobank Study

Both extrinsic and intrinsic factors predispose older people to fall. We performed a genome-wide association analysis to investigate how much of an individual’s fall susceptibility can be attributed to genetics in 89,076 cases and 362,103 controls from the UK Biobank Study. The analysis revealed a small, but significant SNP-based heritability (2.7%) and identified three novel fall-associated loci (P_combined ≤ 5 × 10⁻⁸). Polygenic risk scores in two independent settings showed patterns of polygenic inheritance. Risk of falling had positive genetic correlations with fractures, identifying for the first time a pathway independent of bone mineral density. There were also positive genetic correlations with insomnia, neuroticism, depressive symptoms, and different medications. Negative genetic correlations were identified with muscle strength, intelligence and subjective well-being. Brain, and in particular cerebellum tissue, showed the highest gene expression enrichment for fall-associated variants. Overall, despite the highly heterogenic nature underlying fall risk, a proportion of the susceptibility can be attributed to genetics.

Katerina Trajanoska et al. report a genome-wide association study of self-reported falls in UK Biobank participants. They identify three novel fall-associated loci and find that risk of falling shows patterns of polygenic inheritance and a SNP-based heritability of 2.7%.

Skeletal maturation in relation to ethnic background in children of school age: The Generation R Study

Ethnicity is a well-established determinant of pediatric maturity, but the underlying genetic and environmental contributions to these ethnic differences are poorly comprehended. We aimed to evaluate the influence of ethnicity on skeletal age (SA), an assessment of pediatric maturation widely used in clinical settings. We included children from the Generation R Study, a multiethnic population-based pregnancy cohort, assessed at a mean age of 9.78 (±0.33) years. SA was evaluated by a trained observer on hand DXA scans using the Greulich and Pyle method. Ethnic background was defined as geographic ancestry (questionnaire-based assessment) (N = 5325) and genetic ancestry (based on admixture analysis) (N = 3413). Associations between the ethnic background and SA were investigated separately in boys and girls, using linear regression models adjusted for age, height and BMI. Based on geographic ancestry, 84% of the children were classified as European, 6% as Asian and 10% as African. Children of European background had on average younger SA than those of Asian or African descent. Asian boys had 0.46 (95% CI 0.26-0.66, p-value < 0.0001) and African boys 0.36 years (95% CI 0.20-0.53, p-value < 0.0001) older SA as compared to European boys. Similarly, Asian girls showed 0.64 (95% CI 0.51-0.77, p-value < 0.0001) and African girls 0.38 years (95% CI 0.27-0.48, p-value < 0.0001) older SA as compared to European girls. A similar pattern was observed in the analysis with genetically-defined ancestry. Furthermore, an increase in the proportion of Asian or African component was associated with older SA in both boys (log[Non-European/European]proportion = 0.10, 95% CI 0.06-0.13, p-value < 0.0001) and girls (log[Non-European/European]proportion = 0.06, 95% CI 0.04-0.08, p-value < 0.0001). In summary, children of Asian and African backgrounds have on average older SA as compared to children of European descent, partially explained by a genetic component.

Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure

Heart failure (HF) is a leading cause of morbidity and mortality worldwide. A small proportion of HF cases are attributable to monogenic cardiomyopathies and existing genome-wide association studies (GWAS) have yielded only limited insights, leaving the observed heritability of HF largely unexplained. We report results from a GWAS meta-analysis of HF comprising 47,309 cases and 930,014 controls. Twelve independent variants at 11 genomic loci are associated with HF, all of which demonstrate one or more associations with coronary artery disease (CAD), atrial fibrillation, or reduced left ventricular function, suggesting shared genetic aetiology. Functional analysis of non-CAD-associated loci implicate genes involved in cardiac development (MYOZ1, SYNPO2L), protein homoeostasis (BAG3), and cellular senescence (CDKN1A). Mendelian randomisation analysis supports causal roles for several HF risk factors, and demonstrates CAD-independent effects for atrial fibrillation, body mass index, and hypertension. These findings extend our knowledge of the pathways underlying HF and may inform new therapeutic strategies.

Influence of Dietary Approaches to Stop Hypertension-Type Diet, Known Genetic Variants and Their Interplay on Blood Pressure in Early Childhood

There is limited evidence on association between adherence to the Dietary Approaches to Stop Hypertension (DASH diet) and a lower blood pressure (BP) in children. In a population-based cohort study, among 1068 Dutch children aged 5 to 7, we evaluated the association between a DASH-type diet, 29 known genetic variants incorporated in a genetic risk score, and their interaction on BP. We calculated DASH score based on the food intake data measured through a validated 71-item food frequency questionnaire. In our sample, DASH score ranged from 9 (low adherence to the DASH diet) to 33 (median=21), and genetic score ranged from 18 (low genetic risk on high BP) to 41 (median=29). After adjustment for covariates, each 10 unit increase in DASH score was associated with a lower systolic BP of 0.7 mm Hg ( P =0.033). DASH score was negatively associated with hypertension (odds ratio=0.96 [0.92–0.99], P =0.044). Similarly, each SD increment in genetic score was associated with 0.5 mm Hg higher diastolic BP ( P =0.002). We found a positive interaction between low DASH score and high genetic score on diastolic BP adjusted for BP risk factors (β=1.52, P interaction =0.019 in additive scale and β=0.03, P interaction =0.021 in multiplicative scale). Our findings show that adherence to the DASH-type diet, as well as a low (adult-derived) genetic risk profile for BP, is associated with lower BP in children and that the genetic basis of BP phenotypes at least partly overlaps between adults and children. In addition, we found evidence of a gene-diet interaction on BP in children.

Relationship between gut microbiota and circulating metabolites in population-based cohorts

Gut microbiota has been implicated in major diseases affecting the human population and has also been linked to triglycerides and high-density lipoprotein levels in the circulation. Recent development in metabolomics allows classifying the lipoprotein particles into more details. Here, we examine the impact of gut microbiota on circulating metabolites measured by Nuclear Magnetic Resonance technology in 2309 individuals from the Rotterdam Study and the LifeLines-DEEP cohort. We assess the relationship between gut microbiota and metabolites by linear regression analysis while adjusting for age, sex, body-mass index, technical covariates, medication use, and multiple testing. We report an association of 32 microbial families and genera with very-low-density and high-density subfractions, serum lipid measures, glycolysis-related metabolites, ketone bodies, amino acids, and acute-phase reaction markers. These observations provide insights into the role of microbiota in host metabolism and support the potential of gut microbiota as a target for therapeutic and preventive interventions.

KCND3 potassium channel gene variant confers susceptibility to electrocardiographic early repolarization pattern

BACKGROUNDThe presence of an early repolarization pattern (ERP) on the surface ECG is associated with risk of ventricular fibrillation and sudden cardiac death. Family studies have shown that ERP is a highly heritable trait, but molecular genetic determinants are unknown.METHODSTo identify genetic susceptibility loci for ERP, we performed a GWAS and meta-analysis in 2,181 cases and 23,641 controls of European ancestry.RESULTSWe identified a genome-wide significant (P < 5 × 10-8) locus in the potassium voltage-gated channel subfamily D member 3 (KCND3) gene that was successfully replicated in additional 1,124 cases and 12,510 controls. A subsequent joint meta-analysis of the discovery and replication cohorts identified rs1545300 as the lead SNP at the KCND3 locus (OR 0.82 per minor T allele, P = 7.7 × 10-12) but did not reveal additional loci. Colocalization analyses indicate causal effects of KCND3 gene expression levels on ERP in both cardiac left ventricle and tibial artery.CONCLUSIONSIn this study, we identified for the first time to our knowledge a genome-wide significant association of a genetic variant with ERP. Our findings of a locus in the KCND3 gene provide insights not only into the genetic determinants but also into the pathophysiological mechanism of ERP, discovering a promising candidate for functional studies.FUNDINGThis project was funded by the German Center for Cardiovascular Research (DZHK Shared Expertise SE081 - STATS). For detailed funding information per study, see the Supplemental Acknowledgments.

Association of Birth Weight With Type 2 Diabetes and Glycemic Traits: A Mendelian Randomization Study

IMPORTANCE

Observational studies have shown associations of birth weight with type 2 diabetes (T2D) and glycemic traits, but it remains unclear whether these associations represent causal associations.

OBJECTIVE

To test the association of birth weight with T2D and glycemic traits using a mendelian randomization analysis.

DESIGN, SETTING, AND PARTICIPANTS

This mendelian randomization study used a genetic risk score for birth weight that was constructed with 7 genome-wide significant single-nucleotide polymorphisms. The associations of this score with birth weight and T2D were tested in a mendelian randomization analysis using study-level data. The association of birth weight with T2D was tested using both study-level data (7 single-nucleotide polymorphisms were used as an instrumental variable) and summary-level data from the consortia (43 single-nucleotide polymorphisms were used as an instrumental variable). Data from 180 056 participants from 49 studies were included.

MAIN OUTCOMES AND MEASURES

Type 2 diabetes and glycemic traits.

RESULTS

This mendelian randomization analysis included 49 studies with 41 155 patients with T2D and 80 008 control participants from study-level data and 34 840 patients with T2D and 114 981 control participants from summary-level data. Study-level data showed that a 1-SD decrease in birth weight due to the genetic risk score was associated with higher risk of T2D among all participants (odds ratio [OR], 2.10; 95% CI, 1.69-2.61; P = 4.03 × 10-5), among European participants (OR, 1.96; 95% CI, 1.42-2.71; P = .04), and among East Asian participants (OR, 1.39; 95% CI, 1.18-1.62; P = .04). Similar results were observed from summary-level analyses. In addition, each 1-SD lower birth weight was associated with 0.189 SD higher fasting glucose concentration (β = 0.189; SE = 0.060; P = .002), but not with fasting insulin, 2-hour glucose, or hemoglobin A1c concentration.

CONCLUSIONS AND RELEVANCE

In this study, a genetic predisposition to lower birth weight was associated with increased risk of T2D and higher fasting glucose concentration, suggesting genetic effects on retarded fetal growth and increased diabetes risk that either are independent of each other or operate through alterations of integrated biological mechanisms.

GWAS on longitudinal growth traits reveals different genetic factors influencing infant, child, and adult BMI

Early childhood growth patterns are associated with adult health, yet the genetic factors and the developmental stages involved are not fully understood. Here, we combine genome-wide association studies with modeling of longitudinal growth traits to study the genetics of infant and child growth, followed by functional, pathway, genetic correlation, risk score, and colocalization analyses to determine how developmental timings, molecular pathways, and genetic determinants of these traits overlap with those of adult health. We found a robust overlap between the genetics of child and adult body mass index (BMI), with variants associated with adult BMI acting as early as 4 to 6 years old. However, we demonstrated a completely distinct genetic makeup for peak BMI during infancy, influenced by variation at the LEPR/LEPROT locus. These findings suggest that different genetic factors control infant and child BMI. In light of the obesity epidemic, these findings are important to inform the timing and targets of prevention strategies.

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10^-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r² > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10^-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Publisher Correction: Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

B-vitamins and body composition: integrating observational and experimental evidence from the B-PROOF study

Purpose

Higher folate and vitamin-B12 have been linked to lower risk of overweight. However, whether this is a causal effect of these B-vitamins on obesity risk remains unclear and evidence in older individuals is scarce. This study aimed to assess the role of B-vitamin supplementation and levels on body composition in older individuals.

Methods

A double-blind, randomized controlled trial in 2919 participants aged ≥ 65 years with elevated homocysteine levels. The intervention comprised a 2-year supplementation with a combination of folic acid (400 µg) and vitamin B12 (500 µg), or with placebo. Serum folate, vitamin-B12, active vitamin-B12 (HoloTC), methylmalonic acid (MMA), and anthropometrics were measured at baseline and after 2 years of follow-up. Dietary intake of folate and vitamin-B12 was measured at baseline in a subsample (n = 603) using a validated food-frequency questionnaire. Fat mass index (FMI) and fat-free mass index (FFMI) were assessed with Dual Energy X-ray absorptiometry (DXA).

Results

Cross-sectional analyses showed that a 1 nmol/L higher serum folate was associated with a 0.021 kg/m² lower BMI (95% CI − 0.039; − 0.004). Higher HoloTC (per pmol/L log-transformed) was associated with a 0.955 kg/m² higher FMI (95% CI 0.262; 1.647), and higher MMA (per μgmol/L) was associated with a 1.108 kg/m² lower FMI (95% CI − 1.899; − 0.316). However, random allocation of B-vitamins did not have a significant effect on changes in BMI, FMI or FFMI during 2 years of intervention.

Conclusions

Although observational data suggested that folate and vitamin B12 status are associated with body composition, random allocation of a supplement with both B-vitamins combined versus placebo did not confirm an effect on BMI or body composition.

Electronic supplementary material

The online version of this article (10.1007/s00394-019-01985-8) contains supplementary material, which is available to authorized users.

GWAS of bone size yields twelve loci that also affect height, BMD, osteoarthritis or fractures

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10-42, β = -0.090) and confers risk of hip fracture (P = 1.0 × 10-8, OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

Femoral stress is prominently associated with fracture risk in children: The Generation R Study

Bone modeling is an important process in the growing skeleton. An inadequate bone modeling in response to mechanical loads would lead some children to develop weaker bones than others. The resulting higher stresses in the bones would render them more susceptible to fracture. We aimed to examine the association between femoral stress (FS) derived from structural parameters and BMD in relation to incident fractures in children. Bone stress was evaluated at the medial femoral neck, a skeletal site subject to large forces during normal locomotion. This study comprises 1840 children from the Generation R Study, with whole body and hip DXA scans at a mean age of 6.01 years. Hip structural analysis (HSA) was used to measure femur geometry for the FS calculation. Data on fractures occurring over the following 4 years after the DXA assessment were obtained by questionnaire. Incident fracture was observed in 7.6% of the participating children. Cox-multivariate regression analysis, described as hazard ratios (HR), showed that after adjustment for sex, ethnicity, age, weight and lean mass fraction, there was a significant increase in the risk of incident fracture for every standard deviation (SD) decrease in total body BMD (HR: 1.35, 95% CI 1.05-1.74, p-value = 0.021), femoral neck BMD (HR: 1.31, 95% CI 1.09-1.58, p-value = 0.005) and narrow neck BMD (HR: 1.39, 95% CI 1.14-1.68, p-value = 0.001). Whereas, every increment of one SD in femoral stress resulted in 1.33 increased risk of incident fractures (HR: 1.33, 95% CI 1.13-1.57, p-value = 0.001). This association remained (borderline) significant after the adjustment for DXA derived BMD measurements. Our results show that increased bone stress may underlie greater susceptibility to traumatic fractures in children (partially independent of BMD) and underscore the utility of hip DXA scans for the assessment of paediatric bone health and specifically fracture risk.

Fractures in school age children in relation to sex and ethnic background: The Generation R Study

Fracture rate in childhood is increasing and its consequences may affect health and developmental processes and cause school absence and restricted activity days. There are scarce epidemiologic studies regarding fractures in children. The aim of this study was to evaluate if pediatric fractures show disparities across sexes and ethnic groups. This study was conducted based on data from 3632 participants of the Generation R Study. Prevalent fractures were assessed using a questionnaire at a mean age of 9.7 years. Child's ethnicity was determined based on country of birth of the parents using questionnaires (geographic ancestry) or admixture analysis (genetic ancestry). Associations between fracture occurrence and sex or ethnicity were evaluated using logistic regression models adjusted for age, weight, lean mass fraction, bone mineral density (BMD) and sex/ethnicity. Fracture was reported for 525 (14.5%) children. The great majority of these children were classified as European (N = 3164), followed by African (N = 283) and Asian (N = 185) based on geographic ancestry. Similarly, the highest proportion of Europeans was observed based on genetic ancestry. Prevalence of fractures was not different between boys and girls, even after adjustment for possible confounders (OR: 1.03, 95% CI 0.84-1.27, p-value = 0.8). However, odds of prevalent fractures were two times higher in European when compared to Asian children (OR: 2.01, 95% CI 1.17-3.45, p-value = 0.01), and 1.5 times higher when compared to African children (OR: 1.50, 95% CI 1.00-2.26, p-value = 0.05). Overall, in this study, European children showed a highest risk of prevalent fractures independently of factors such as body composition and BMD, while no difference in the prevalence of fractures between boys and girls was observed.

Non-genetic and genetic predictors of a superficial first basal cell carcinoma

Background

Several observational studies have suggested differences in the risk factor profile between patients with superficial basal cell carcinomas (BCCs) and non‐superficial BCCs.

Objective

To test the reproducibility of previous study findings and to find new genetic and non‐genetic predictors for patients with a superficial first BCC.

Methods

A total of 14.628 participants of northwestern European descent aged 45 years or older from a prospective population‐based cohort study (Rotterdam Study) were linked with the Dutch Pathology Registry (PALGA) of whom 1528 were identified as BCC patients. After exclusion, 948 eligible BCC patients remained for further non‐genetic analyses and 1014 for genetic analyses. We included 11 phenotypic, environmental and tumour‐specific characteristics, and 20 candidate single nucleotide polymorphisms (SNP) as potential predictors for patients with a superficial first BCC. We performed binary logistic multivariable regression analyses.

Results

We found that patients with a superficial first BCC were significantly younger, almost two times more often female and 12–18 times more likely to have their BCC on the trunk or extremities than patients with a non‐superficial first BCC. One SNP (rs12203592), mapped to IRF4, looked promising (OR 1.83, 95% CI 1.13–2.97, P‐value <0.05), but after adjustment for multiple testing, no significant differences in genetic make‐up between superficial BCC and non‐superficial BCC patients were found.

Conclusion

We conclude that patients with a superficial BCC differ from non‐superficial BCC patients with respect to environmental factors (tumour localization as a proxy for UVR exposure) and phenotypic characteristics (age and sex), but we found no difference in genotype. As superficial BCC patients develop their first BCCs at a younger age, they could be at higher lifetime risk for subsequent skin cancers and therefore be an important group for secondary prevention.

Linked article: This article is commented on G. Argenziano et. al., p. 461 in this issue. To view this article visit https://doi.org/10.1111/jdv.15486.

The effect of antimicrobial drug use on the composition of the genitourinary microbiota in an elderly population

BACKGROUND

The urinary tract is inhabited by a diversity of microorganisms, known as the genitourinary microbiota. Here, we investigated the association between the use of antimicrobial drugs and the composition of the genitourinary microbiota.

RESULTS

Clean-catch urinary samples were collected from 27 participants of the Rotterdam Study. Bacterial DNA was extracted and the 16S ribosomal RNA gene variable regions V3 and V4 were analyzed using Illumina sequencing. 23 of the 27 participants were included in the analysis. The population consisted of 10 men and 13 women with a mean age of 75 ± 3 years. The time between the last prescription of an antimicrobial drug and sampling was determined and categorized. The use of antimicrobial drugs prior to urine sampling was associated with statistically significant differences in the beta-diversity of the genitourinary microbiota. No association was found between antimicrobial drug use and the alpha-diversity of the genitourinary microbiota. Operational Taxonomic Units (OTUs) that were lowest in participants who used antimicrobial drug belonged to Lactobacillus and Finegoldia. In contrast, an OTU belonging to the genus Parabacteroides had higher abundances. Also, an OTU belonging to the species E.coli was higher in the participants who used antimicrobial drugs.

CONCLUSION

Prior use of antimicrobial drugs is associated with a different composition of the genitourinary microbiota. Our results might indicate a persisting effect of antimicrobial drugs on the composition of the microbiota, but reverse causality cannot be ruled out. Future studies are needed to differentiate between two possibilities. Genitourinary dysbiosis could be the result of antimicrobial drug use or genitourinary dysbiosis could be a risk factor for urinary tract infections resulting in increased use of antimicrobial drugs. This may have important implications for treatment and prevention of (recurrent) UTIs.

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.

Association of maternal prenatal smoking GFI1-locus and cardio-metabolic phenotypes in 18,212 adults

BACKGROUND

DNA methylation at the GFI1-locus has been repeatedly associated with exposure to smoking from the foetal period onwards. We explored whether DNA methylation may be a mechanism that links exposure to maternal prenatal smoking with offspring's adult cardio-metabolic health.

METHODS

We meta-analysed the association between DNA methylation at GFI1-locus with maternal prenatal smoking, adult own smoking, and cardio-metabolic phenotypes in 22 population-based studies from Europe, Australia, and USA (n = 18,212). DNA methylation at the GFI1-locus was measured in whole-blood. Multivariable regression models were fitted to examine its association with exposure to prenatal and own adult smoking. DNA methylation levels were analysed in relation to body mass index (BMI), waist circumference (WC), fasting glucose (FG), high-density lipoprotein cholesterol (HDL-C), triglycerides (TG), diastolic, and systolic blood pressure (BP).

FINDINGS

Lower DNA methylation at three out of eight GFI1-CpGs was associated with exposure to maternal prenatal smoking, whereas, all eight CpGs were associated with adult own smoking. Lower DNA methylation at cg14179389, the strongest maternal prenatal smoking locus, was associated with increased WC and BP when adjusted for sex, age, and adult smoking with Bonferroni-corrected P < 0·012. In contrast, lower DNA methylation at cg09935388, the strongest adult own smoking locus, was associated with decreased BMI, WC, and BP (adjusted 1 × 10-7 < P < 0.01). Similarly, lower DNA methylation at cg12876356, cg18316974, cg09662411, and cg18146737 was associated with decreased BMI and WC (5 × 10-8 < P < 0.001). Lower DNA methylation at all the CpGs was consistently associated with higher TG levels.

INTERPRETATION

Epigenetic changes at the GFI1 were linked to smoking exposure in-utero/in-adulthood and robustly associated with cardio-metabolic risk factors. FUND: European Union's Horizon 2020 research and innovation programme under grant agreement no. 633595 DynaHEALTH.

Author Correction: Discovery and refinement of genetic loci associated with cardiometabolic risk using dense imputation maps

In the version of the article published, the surname of author Aaron Isaacs is misspelled as Issacs.

Meta-analysis of epigenome-wide association studies of cognitive abilities

Cognitive functions are important correlates of health outcomes across the life-course. Individual differences in cognitive functions are partly heritable. Epigenetic modifications, such as DNA methylation, are susceptible to both genetic and environmental factors and may provide insights into individual differences in cognitive functions. Epigenome-wide meta-analyses for blood-based DNA methylation levels at ~420,000 CpG sites were performed for seven measures of cognitive functioning using data from 11 cohorts. CpGs that passed a Bonferroni correction, adjusting for the number of CpGs and cognitive tests, were assessed for: longitudinal change; being under genetic control (methylation QTLs); and associations with brain health (structural MRI), brain methylation and Alzheimer's disease pathology. Across the seven measures of cognitive functioning (meta-analysis n range: 2557-6809), there were epigenome-wide significant (P < 1.7 × 10-8) associations for global cognitive function (cg21450381, P = 1.6 × 10-8), and phonemic verbal fluency (cg12507869, P = 2.5 × 10-9). The CpGs are located in an intergenic region on chromosome 12 and the INPP5A gene on chromosome 10, respectively. Both probes have moderate correlations (~0.4) with brain methylation in Brodmann area 20 (ventral temporal cortex). Neither probe showed evidence of longitudinal change in late-life or associations with white matter brain MRI measures in one cohort with these data. A methylation QTL analysis suggested that rs113565688 was a cis methylation QTL for cg12507869 (P = 5 × 10-5 and 4 × 10-13 in two lookup cohorts). We demonstrate a link between blood-based DNA methylation and measures of phonemic verbal fluency and global cognitive ability. Further research is warranted to understand the mechanisms linking genomic regulatory changes with cognitive function to health and disease.

Whole-Genome Linkage Scan Combined With Exome Sequencing Identifies Novel Candidate Genes for Carotid Intima-Media Thickness

Carotid intima-media thickness (cIMT) is an established heritable marker for subclinical atherosclerosis. In this study, we aim to identify rare variants with large effects driving differences in cIMT by performing genome-wide linkage analysis of individuals in the extremes of cIMT trait distribution (>90th percentile) in a large family-based study from a genetically isolated population in the Netherlands. Linked regions were subsequently explored by fine-mapping using exome sequencing. We observed significant evidence of linkage on chromosomes 2p16.3 [rs1017418, heterogeneity LOD (HLOD) = 3.35], 19q13.43 (rs3499, HLOD = 9.09), 20p13 (rs1434789, HLOD = 4.10), and 21q22.12 (rs2834949, HLOD = 3.59). Fine-mapping using exome sequencing data identified a non-coding variant (rs62165235) in PNPT1 gene under the linkage peak at chromosome 2 that is likely to have a regulatory function. The variant was associated with quantitative cIMT in the family-based study population (effect = 0.27, p-value = 0.013). Furthermore, we identified several genes under the linkage peak at chromosome 21 highly expressed in tissues relevant for atherosclerosis. To conclude, our linkage analysis identified four genomic regions significantly linked to cIMT. Further analyses are needed to demonstrate involvement of identified candidate genes in development of atherosclerosis.

DNA methylation age is associated with an altered hemostatic profile in a multiethnic meta-analysis

Many hemostatic factors are associated with age and age-related diseases; however, much remains unknown about the biological mechanisms linking aging and hemostatic factors. DNA methylation is a novel means by which to assess epigenetic aging, which is a measure of age and the aging processes as determined by altered epigenetic states. We used a meta-analysis approach to examine the association between measures of epigenetic aging and hemostatic factors, as well as a clotting time measure. For fibrinogen, we performed European and African ancestry-specific meta-analyses which were then combined via a random effects meta-analysis. For all other measures we could not estimate ancestry-specific effects and used a single fixed effects meta-analysis. We found that 1-year higher extrinsic epigenetic age as compared with chronological age was associated with higher fibrinogen (0.004 g/L/y; 95% confidence interval, 0.001-0.007; P = .01) and plasminogen activator inhibitor 1 (PAI-1; 0.13 U/mL/y; 95% confidence interval, 0.07-0.20; P = 6.6 × 10-5) concentrations, as well as lower activated partial thromboplastin time, a measure of clotting time. We replicated PAI-1 associations using an independent cohort. To further elucidate potential functional mechanisms, we associated epigenetic aging with expression levels of the PAI-1 protein encoding gene (SERPINE1) and the 3 fibrinogen subunit-encoding genes (FGA, FGG, and FGB) in both peripheral blood and aorta intima-media samples. We observed associations between accelerated epigenetic aging and transcription of FGG in both tissues. Collectively, our results indicate that accelerated epigenetic aging is associated with a procoagulation hemostatic profile, and that epigenetic aging may regulate hemostasis in part via gene transcription.

Fracture incidence and secular trends between 1989 and 2013 in a population based cohort: The Rotterdam Study

Fracture incidence needs to be evaluated over time to assess the impact of the enlarging population burden of fractures (due to increase in lifespan) and the efficacy of fracture prevention strategies. Therefore, we aimed to evaluate the association of femoral neck bone mineral density (FN-BMD) measured using dual-energy X-ray absorptiometry (DXA) at baseline with fracture risk over a long follow-up time period. Incident non-vertebral fractures were assessed in 14,613 individuals participating in the Rotterdam Study with up to 20 years of follow-up. During a mean follow-up of 10.7 ± 6.2 years, 2971 (20.3%) participants had at least one incident non-vertebral fracture. The risk for any non-vertebral fracture was 1.37 (95% Confidence Interval (CI): 1.25-1.49) and 1.42 (95%CI: 1.35-1.50) for men and women, respectively. The majority (79% in men and 75% in women) of all fractures occurred among participants a normal or osteopenic T-score. The incidence rates per 1000 person-years for the most common fractures were 5.3 [95%CI: 5.0-5.7] for hip, 4.9 [95%CI: 4.6-5.3] for wrist and 2.3 [95%CI: 2.0-2.5] for humerus. To examine the predictive ability of BMD through follow-up time we determined fracture hazard ratios (HR) per standard deviation decrease in femoral neck BMD across five year bins. No differences were observed, with a HR of 2.5 (95%CI: 2.0-3.1) after the first 5 years, and of 1.9 (95%CI: 1.1-3.3) after 20 years. To assess secular trends in fracture incidence at all skeletal sites we compared participants at an age of 70-80 years across two time periods: 1989-2001 (n = 2481, 60% women) and 2001-2013 (n = 2936, 58% women) and found no statistically significant difference (p < 0.05) between fracture incidence rates (i.e., incidence of non-vertebral fractures of 26.4 per 1000 PY [95%CI: 24.4-28.5]) between 1989 and 2001, and of 25.4 per 1000 PY [95%CI: 23.0-28.0] between 2001 and 2013. In conclusion, BMD is still predictive of future fracture over a long period of time. While no secular changes in fractures rates seem to be observed after a decade, the majority of fractures still occur above the osteoporosis threshold, emphasizing the need to improve the screening of osteopenic patients.

DNA Methylation Signatures of Depressive Symptoms in Middle-aged and Elderly Persons: Meta-analysis of Multiethnic Epigenome-wide Studies

IMPORTANCE

Depressive disorders arise from a combination of genetic and environmental risk factors. Epigenetic disruption provides a plausible mechanism through which gene-environment interactions lead to depression. Large-scale, epigenome-wide studies on depression are missing, hampering the identification of potentially modifiable biomarkers.

OBJECTIVE

To identify epigenetic mechanisms underlying depression in middle-aged and elderly persons, using DNA methylation in blood.

DESIGN, SETTING, AND PARTICIPANTS

To date, the first cross-ethnic meta-analysis of epigenome-wide association studies (EWAS) within the framework of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium was conducted. The discovery EWAS included 7948 individuals of European origin from 9 population-based cohorts. Participants who were assessed for both depressive symptoms and whole-blood DNA methylation were included in the study. Results of EWAS were pooled using sample-size weighted meta-analysis. Replication of the top epigenetic sites was performed in 3308 individuals of African American and European origin from 2 population-based cohorts.

MAIN OUTCOMES AND MEASURES

Whole-blood DNA methylation levels were assayed with Illumina-Infinium Human Methylation 450K BeadChip and depressive symptoms were assessed by questionnaire.

RESULTS

The discovery cohorts consisted of 7948 individuals (4104 [51.6%] women) with a mean (SD) age of 65.4 (5.8) years. The replication cohort consisted of 3308 individuals (2456 [74.2%] women) with a mean (SD) age of 60.3 (6.4) years. The EWAS identified methylation of 3 CpG sites to be significantly associated with increased depressive symptoms: cg04987734 (P = 1.57 × 10-08; n = 11 256; CDC42BPB gene), cg12325605 (P = 5.24 × 10-09; n = 11 256; ARHGEF3 gene), and an intergenic CpG site cg14023999 (P = 5.99 × 10-08; n = 11 256; chromosome = 15q26.1). The predicted expression of the CDC42BPB gene in the brain (basal ganglia) (effect, 0.14; P = 2.7 × 10-03) and of ARHGEF3 in fibroblasts (effect, -0.48; P = 9.8 × 10-04) was associated with major depression.

CONCLUSIONS AND RELEVANCE

This study identifies 3 methylated sites associated with depressive symptoms. All 3 findings point toward axon guidance as the common disrupted pathway in depression. The findings provide new insights into the molecular mechanisms underlying the complex pathophysiology of depression. Further research is warranted to determine the utility of these findings as biomarkers of depression and evaluate any potential role in the pathophysiology of depression and their downstream clinical effects.