Abstract WP216: Epigenome-Wide Association Study of Cerebral White Matter Hyperintensities

Background: Cerebral white matter hyperintensities (WMH) on MRI are typical markers of small vessel disease and strong radiological correlates of stroke and age-related cognitive decline. Despite a recognized high heritability, the molecular basis of WMH has not been fully characterized. In particular, the contribution of epigenetic modifications, such as DNA methylation, has not been explored.

Methods: We conducted a meta-analysis of epigenome-wide association studies of WMH burden in up to 5777 elderly participants of European and African ancestry from 9 population-based cohorts using blood-derived DNA methylation measured on the Infinium HumanMethylation450 BeadChip. Cohort-specific associations between burden of WMH and DNA methylation beta values were estimated using linear mixed-effect models and combined in a sample-size weighted fixed-effect meta-analysis. In addition, we used two different approaches to identify differentially methylated regions (DMRs), which may be more informative than individual loci. Bonferroni correction and False Discovery Rates were used to account for the multiple tests.

Results: Single-site analyses identified a CpG site within the CLDN5 gene significantly associated with WMH (P= 2.3x10 -7 ). Claudin-5 is an endothelial-specific component in the brain vasculature and regulates blood brain barrier permeability. Region-based analyses, which leverage the correlations between nearby CpG sites identified several DMRs significantly associated with WMH burden including: PRMT1, BTBD17, and IFITM10 (P=1.4x10 -10 , 2.3x10 -8 , and 3.6x10 -7 , respectively). PRMT1 encodes the Protein Arginine N-methylase 1, which methylates histones in genes involved in glioblastomagenesis. Mice lacking this gene are characterized by severe defects in oligodendrocyte maturation processes. The function of BTBD17 and IFITM10 is not characterized but both are expressed in the brain and exhibit changes in expression in response to viral infection.

Conclusions: Consistent with our previously reported genetic association studies, this genome-wide DNA methylation analysis supports a role of genes involved in glial cell function in WMH etiology. It also suggests a novel role of genes involved in viral response.

Genetic and lifestyle risk factors for MRI-defined brain infarcts in a population-based setting

OBJECTIVE

To explore genetic and lifestyle risk factors of MRI-defined brain infarcts (BI) in large population-based cohorts.

METHODS

We performed meta-analyses of genome-wide association studies (GWAS) and examined associations of vascular risk factors and their genetic risk scores (GRS) with MRI-defined BI and a subset of BI, namely, small subcortical BI (SSBI), in 18 population-based cohorts (n = 20,949) from 5 ethnicities (3,726 with BI, 2,021 with SSBI). Top loci were followed up in 7 population-based cohorts (n = 6,862; 1,483 with BI, 630 with SBBI), and we tested associations with related phenotypes including ischemic stroke and pathologically defined BI.

RESULTS

The mean prevalence was 17.7% for BI and 10.5% for SSBI, steeply rising after age 65. Two loci showed genome-wide significant association with BI: FBN2, p = 1.77 × 10-8; and LINC00539/ZDHHC20, p = 5.82 × 10-9. Both have been associated with blood pressure (BP)-related phenotypes, but did not replicate in the smaller follow-up sample or show associations with related phenotypes. Age- and sex-adjusted associations with BI and SSBI were observed for BP traits (p value for BI, p [BI] = 9.38 × 10-25; p [SSBI] = 5.23 × 10-14 for hypertension), smoking (p [BI] = 4.4 × 10-10; p [SSBI] = 1.2 × 10-4), diabetes (p [BI] = 1.7 × 10-8; p [SSBI] = 2.8 × 10-3), previous cardiovascular disease (p [BI] = 1.0 × 10-18; p [SSBI] = 2.3 × 10-7), stroke (p [BI] = 3.9 × 10-69; p [SSBI] = 3.2 × 10-24), and MRI-defined white matter hyperintensity burden (p [BI] = 1.43 × 10-157; p [SSBI] = 3.16 × 10-106), but not with body mass index or cholesterol. GRS of BP traits were associated with BI and SSBI (p ≤ 0.0022), without indication of directional pleiotropy.

CONCLUSION

In this multiethnic GWAS meta-analysis, including over 20,000 population-based participants, we identified genetic risk loci for BI requiring validation once additional large datasets become available. High BP, including genetically determined, was the most significant modifiable, causal risk factor for BI.

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.

Neandertal Introgression Sheds Light on Modern Human Endocranial Globularity

One of the features that distinguishes modern humans from our extinct relatives and ancestors is a globular shape of the braincase [1-4]. As the endocranium closely mirrors the outer shape of the brain, these differences might reflect altered neural architecture [4, 5]. However, in the absence of fossil brain tissue, the underlying neuroanatomical changes as well as their genetic bases remain elusive. To better understand the biological foundations of modern human endocranial shape, we turn to our closest extinct relatives: the Neandertals. Interbreeding between modern humans and Neandertals has resulted in introgressed fragments of Neandertal DNA in the genomes of present-day non-Africans [6, 7]. Based on shape analyses of fossil skull endocasts, we derive a measure of endocranial globularity from structural MRI scans of thousands of modern humans and study the effects of introgressed fragments of Neandertal DNA on this phenotype. We find that Neandertal alleles on chromosomes 1 and 18 are associated with reduced endocranial globularity. These alleles influence expression of two nearby genes, UBR4 and PHLPP1, which are involved in neurogenesis and myelination, respectively. Our findings show how integration of fossil skull data with archaic genomics and neuroimaging can suggest developmental mechanisms that may contribute to the unique modern human endocranial shape.

GWAS and colocalization analyses implicate carotid intima-media thickness and carotid plaque loci in cardiovascular outcomes

Carotid artery intima media thickness (cIMT) and carotid plaque are measures of subclinical atherosclerosis associated with ischemic stroke and coronary heart disease (CHD). Here, we undertake meta-analyses of genome-wide association studies (GWAS) in 71,128 individuals for cIMT, and 48,434 individuals for carotid plaque traits. We identify eight novel susceptibility loci for cIMT, one independent association at the previously-identified PINX1 locus, and one novel locus for carotid plaque. Colocalization analysis with nearby vascular expression quantitative loci (cis-eQTLs) derived from arterial wall and metabolic tissues obtained from patients with CHD identifies candidate genes at two potentially additional loci, ADAMTS9 and LOXL4. LD score regression reveals significant genetic correlations between cIMT and plaque traits, and both cIMT and plaque with CHD, any stroke subtype and ischemic stroke. Our study provides insights into genes and tissue-specific regulatory mechanisms linking atherosclerosis both to its functional genomic origins and its clinical consequences in humans.

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.

Negative effect of vitamin D on kidney function: a Mendelian randomization study

Background

The kidney plays a central role in the regulation of vitamin D metabolism. It is not clear, however, whether vitamin D influences kidney function. Previous studies have reported conflicting results, which may have been influenced by reverse causation and residual confounding. We conducted a Mendelian randomization (MR) study to obtain unconfounded estimates of the association between genetically instrumented vitamin D metabolites and estimated glomerular filtration rate (eGFR) as well as the urinary albumin:creatinine ratio (UACR).

Methods

We performed a two-sample MR study based on three single nucleotide variants associated with 25(OH)D levels: rs2282679, rs10741657 and rs12785878, related to the genes GC, CYP2R1 and DHCR7, respectively. Estimates of the allele-dependent effects on serum 25(OH)D and eGFR/UACR were obtained from summary statistics of published genome-wide association meta-analyses. Additionally, we performed a one-sample MR analysis for both 25(OH)D and 1,25(OH)2 D using individual-level data from six cohorts.

Results

The combined MR estimate supported a negative causal effect of log transformed 25(OH)D on log transformed eGFR (β = -0.013, P = 0.003). The analysis of individual-level data confirmed the main findings and also revealed a significant association of 1,25(OH)2 D on eGFR (β = -0.094, P = 0.008). These results show that a 10% increase in serum 25(OH)D levels causes a 0.3% decrease in eGFR. There was no effect of 25(OH)D on UACR (β = 0.032, P = 0.265).

Conclusion

Our study suggests that circulating vitamin D metabolite levels are negatively associated with eGFR. Further studies are needed to elucidate the underlying mechanisms.

Reassessing the Association between Circulating Vitamin D and IGFBP-3: Observational and Mendelian Randomization Estimates from Independent Sources

BACKGROUND

Circulating insulin-like growth factor binding protein 3 (IGFBP-3) has been associated with prostate cancer. Preclinical studies found that vitamin D regulates IGFBP-3 expression, although evidence from epidemiologic studies is conflicting.

METHODS

Mendelian randomization analyses (MR) were conducted to reassess associations between IGFBP-3 and prostate cancer risk and advanced prostate cancer using summary statistics from the PRACTICAL consortium (44,825 cases; 27,904 controls). Observational and MR analyses were conducted to assess the relationship between inactive vitamin D [25(OH)D] and IGFBP-3 using data from the ProtecT study (1,366 cases;1,071 controls) and summary statistics from the CHARGE consortium (n = 18,995).

RESULTS

The OR for prostate cancer per SD unit increase in circulating IGFBP-3 was 1.14 [95% confidence interval (CI), 1.02-1.28]. The OR for advanced prostate cancer per SD unit increase in IGFBP-3 was 1.22 (95% CI, 1.07-1.40). Observationally, a SD increase in 25(OH)D was associated with a 0.1SD (95% CI, 0.05-0.14) increase in IGFBP-3. MR analyses found little evidence for a causal relationship between circulating 25(OH)D and IGFBP-3 in the circulation.

CONCLUSIONS

This study provided confirmatory evidence that IGFBP-3 is a risk factor for prostate cancer risk and progression. Observationally, there was evidence that 25(OH)D is associated with IGFBP-3, but MR analyses suggested that these findings were unlikely to be causal. Findings may be limited by the nature of instrumentation of 25(OH)D and IGFBP-3 and the utility of circulating measures. 25(OH)D appears unlikely to be causally related to IGFBP-3 in the circulation, however, our findings do not preclude causal associations at the tissue level.

IMPACT

IGFBP-3 is a prostate cancer risk factor but 25(OH)D are unlikely to be causally related to IGFBP-3 in the circulation.

Publisher Correction: Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits

In the version of this article originally published, the name of author Martin H. de Borst was coded incorrectly in the XML. The error has now been corrected in the HTML version of the paper.

Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits

High blood pressure is a highly heritable and modifiable risk factor for cardiovascular disease. We report the largest genetic association study of blood pressure traits (systolic, diastolic and pulse pressure) to date in over 1 million people of European ancestry. We identify 535 novel blood pressure loci that not only offer new biological insights into blood pressure regulation but also highlight shared genetic architecture between blood pressure and lifestyle exposures. Our findings identify new biological pathways for blood pressure regulation with potential for improved cardiovascular disease prevention in the future.

Meta-analysis of exome array data identifies six novel genetic loci for lung function

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV ₁), forced vital capacity (FVC) and the ratio of FEV ₁ to FVC (FEV ₁/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 ^-7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

PR interval genome-wide association meta-analysis identifies 50 loci associated with atrial and atrioventricular electrical activity

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

An InDel in Phospholipase-C-B-1 Is Linked with Euthyroid Multinodular Goiter

BACKGROUND

Euthyroid multinodular goiter (MNG) is common, but little is known about the genetic variations conferring predisposition. Previously, a family with MNG of adolescent onset was reported in which some family members developed papillary thyroid carcinomas (PTC).

METHODS

Genome-wide linkage analysis and next-generation sequencing were conducted to identify genetic variants that may confer disease predisposition. A multipoint nonparametric LOD score of 3.01 was obtained, covering 19 cM on chromosome 20p. Haplotype analysis reduced the region of interest to 10 cM.

RESULTS

Analysis of copy number variation identified an intronic InDel (∼1000 bp) in the PLCB1 gene in all eight affected family members and carriers (an unaffected person who has inherited the genetic trait). This InDel is present in approximately 1% of "healthy" Caucasians. Next-generation sequencing of the region identified no additional disease-associated variant, suggesting a possible role of the InDel. Since PLCB1 contributes to thyrocyte growth regulation, the InDel was investigated in relevant Caucasian cohorts. It was detected in 0/70 PTC but 4/81 unrelated subjects with MNG (three females; age at thyroidectomy 27-59 years; no family history of MNG/PTC). The InDel frequency is significantly higher in MNG subjects compared to controls (χ² = 5.076; p = 0.024. PLCB1 transcript levels were significantly higher in thyroids with the InDel than without (p < 0.02).

CONCLUSIONS

The intronic PLCB1 InDel is the first variant found in familial multiple papilloid adenomata-type MNG and in a subset of patients with sporadic MNG. It may function through overexpression, and increased PLC activity has been reported in thyroid neoplasms. The potential role of the deletion as a biomarker to identify MNG patients more likely to progress to PTC merits exploration.

Multi-ethnic genome-wide association study for atrial fibrillation

Atrial fibrillation (AF) affects more than 33 million individuals worldwide1 and has a complex heritability2. We conducted the largest meta-analysis of genome-wide association studies (GWAS) for AF to date, consisting of more than half a million individuals, including 65,446 with AF. In total, we identified 97 loci significantly associated with AF, including 67 that were novel in a combined-ancestry analysis, and 3 that were novel in a European-specific analysis. We sought to identify AF-associated genes at the GWAS loci by performing RNA-sequencing and expression quantitative trait locus analyses in 101 left atrial samples, the most relevant tissue for AF. We also performed transcriptome-wide analyses that identified 57 AF-associated genes, 42 of which overlap with GWAS loci. The identified loci implicate genes enriched within cardiac developmental, electrophysiological, contractile and structural pathways. These results extend our understanding of the biological pathways underlying AF and may facilitate the development of therapeutics for AF.

Common Methods for Performing Mendelian Randomization

Mendelian randomization (MR) is a framework for assessing causal inference using cross-sectional data in combination with genetic information. This paper summarizes statistical methods commonly applied and strait forward to use for conducting MR analyses including those taking advantage of the rich dataset of SNP-trait associations that were revealed in the last decade through large-scale genome-wide association studies. Using these data, powerful MR studies are possible. However, the causal estimate may be biased in case the assumptions of MR are violated. The source and the type of this bias are described while providing a summary of the mathematical formulas that should help estimating the magnitude and direction of the potential bias depending on the specific research setting. Finally, methods for relaxing the assumptions and for conducting sensitivity analyses are discussed. Future researches in the field of MR include the assessment of non-linear causal effects, and automatic detection of invalid instruments.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression

Major depressive disorder (MDD) is a common illness accompanied by considerable morbidity, mortality, costs, and heightened risk of suicide. We conducted a genome-wide association meta-analysis based in 135,458 cases and 344,901 controls and identified 44 independent and significant loci. The genetic findings were associated with clinical features of major depression and implicated brain regions exhibiting anatomical differences in cases. Targets of antidepressant medications and genes involved in gene splicing were enriched for smaller association signal. We found important relationships of genetic risk for major depression with educational attainment, body mass, and schizophrenia: lower educational attainment and higher body mass were putatively causal, whereas major depression and schizophrenia reflected a partly shared biological etiology. All humans carry lesser or greater numbers of genetic risk factors for major depression. These findings help refine the basis of major depression and imply that a continuous measure of risk underlies the clinical phenotype.

Multiethnic Meta-Analysis Identifies RAI1 as a Possible Obstructive Sleep Apnea-related Quantitative Trait Locus in Men

Obstructive sleep apnea (OSA) is a common heritable disorder displaying marked sexual dimorphism in disease prevalence and progression. Previous genetic association studies have identified a few genetic loci associated with OSA and related quantitative traits, but they have only focused on single ethnic groups, and a large proportion of the heritability remains unexplained. The apnea-hypopnea index (AHI) is a commonly used quantitative measure characterizing OSA severity. Because OSA differs by sex, and the pathophysiology of obstructive events differ in rapid eye movement (REM) and non-REM (NREM) sleep, we hypothesized that additional genetic association signals would be identified by analyzing the NREM/REM-specific AHI and by conducting sex-specific analyses in multiethnic samples. We performed genome-wide association tests for up to 19,733 participants of African, Asian, European, and Hispanic/Latino American ancestry in 7 studies. We identified rs12936587 on chromosome 17 as a possible quantitative trait locus for NREM AHI in men (N = 6,737; P = 1.7 × 10-8) but not in women (P = 0.77). The association with NREM AHI was replicated in a physiological research study (N = 67; P = 0.047). This locus overlapping the RAI1 gene and encompassing genes PEMT1, SREBF1, and RASD1 was previously reported to be associated with coronary artery disease, lipid metabolism, and implicated in Potocki-Lupski syndrome and Smith-Magenis syndrome, which are characterized by abnormal sleep phenotypes. We also identified gene-by-sex interactions in suggestive association regions, suggesting that genetic variants for AHI appear to vary by sex, consistent with the clinical observations of strong sexual dimorphism.

Genetic Determinants of Circulating Estrogen Levels and Evidence of a Causal Effect of Estradiol on Bone Density in Men

Context

Serum estradiol (E2) and estrone (E1) levels exhibit substantial heritability.

Objective

To investigate the genetic regulation of serum E2 and E1 in men.

Design, Setting, and Participants

Genome-wide association study in 11,097 men of European origin from nine epidemiological cohorts.

Main Outcome Measures

Genetic determinants of serum E2 and E1 levels.

Results

Variants in/near CYP19A1 demonstrated the strongest evidence for association with E2, resolving to three independent signals. Two additional independent signals were found on the X chromosome; FAMily with sequence similarity 9, member B (FAM9B), rs5934505 (P = 3.4 × 10-8) and Xq27.3, rs5951794 (P = 3.1 × 10-10). E1 signals were found in CYP19A1 (rs2899472, P = 5.5 × 10-23), in Tripartite motif containing 4 (TRIM4; rs17277546, P = 5.8 × 10-14), and CYP11B1/B2 (rs10093796, P = 1.2 × 10-8). E2 signals in CYP19A1 and FAM9B were associated with bone mineral density (BMD). Mendelian randomization analysis suggested a causal effect of serum E2 on BMD in men. A 1 pg/mL genetically increased E2 was associated with a 0.048 standard deviation increase in lumbar spine BMD (P = 2.8 × 10-12). In men and women combined, CYP19A1 alleles associated with higher E2 levels were associated with lower degrees of insulin resistance.

Conclusions

Our findings confirm that CYP19A1 is an important genetic regulator of E2 and E1 levels and strengthen the causal importance of E2 for bone health in men. We also report two independent loci on the X-chromosome for E2, and one locus each in TRIM4 and CYP11B1/B2, for E1.