Person:

De Vivo, Immaculata

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, p = 1.8×(10^{−106})), PRMT6 (rs17496332, 1p13.3, p = 1.4×(10^{−11})), GCKR (rs780093, 2p23.3, p = 2.2×(10^{−16})), ZBTB10 (rs440837, 8q21.13, p = 3.4×(10^{−9})), JMJD1C (rs7910927, 10q21.3, p = 6.1×(10^{−35})), SLCO1B1 (rs4149056, 12p12.1, p = 1.9×(10^{−08})), NR2F2 (rs8023580, 15q26.2, p = 8.3×(10^{−12})), ZNF652 (rs2411984, 17q21.32, p = 3.5×(10^{−14})), TDGF3 (rs1573036, Xq22.3, p = 4.1×(10^{−14})), LHCGR (rs10454142, 2p16.3, p = 1.3×(10^{−07}), BAIAP2L1 (rs3779195, 7q21.3, p = 2.7×(10^{−08})), and UGT2B15 (rs293428, 4q13.2, p = 5.5×(10^{−06})). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men p = 2.5×(10^{−08}), women p = 0.66, heterogeneity p = 0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ∼15.6% and ∼8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance.

Genome-wide association studies (GWAS) have successfully identified common genetic variants that contribute to breast cancer risk. Discovering additional variants has become difficult, as power to detect variants of weaker effect with present sample sizes is limited. An alternative approach is to look for variants associated with quantitative traits that in turn affect disease risk. As exposure to high circulating estradiol and testosterone, and low sex hormone-binding globulin (SHBG) levels is implicated in breast cancer etiology, we conducted GWAS analyses of plasma estradiol, testosterone, and SHBG to identify new susceptibility alleles. Cancer Genetic Markers of Susceptibility (CGEMS) data from the Nurses’ Health Study (NHS), and Sisters in Breast Cancer Screening data were used to carry out primary meta-analyses among ∼1600 postmenopausal women who were not taking postmenopausal hormones at blood draw. We observed a genome-wide significant association between SHBG levels and rs727428 (joint (\beta) = -0.126; joint P = 2.09×10–16), downstream of the SHBG gene. No genome-wide significant associations were observed with estradiol or testosterone levels. Among variants that were suggestively associated with estradiol (P<10–5), several were located at the CYP19A1 gene locus. Overall results were similar in secondary meta-analyses that included ∼900 NHS current postmenopausal hormone users. No variant associated with estradiol, testosterone, or SHBG at P<10–5 was associated with postmenopausal breast cancer risk among CGEMS participants. Our results suggest that the small magnitude of difference in hormone levels associated with common genetic variants is likely insufficient to detectably contribute to breast cancer risk.

The Nurses’ Health Study (NHS), Nurses’ Health Study II (NHSII), Health Professionals Follow Up Study (HPFS) and the Physicians Health Study (PHS) have collected detailed longitudinal data on multiple exposures and traits for approximately 310,000 study participants over the last 35 years. Over 160,000 study participants across the cohorts have donated a DNA sample and to date, 20,691 subjects have been genotyped as part of genome-wide association studies (GWAS) of twelve primary outcomes. However, these studies utilized six different GWAS arrays making it difficult to conduct analyses of secondary phenotypes or share controls across studies. To allow for secondary analyses of these data, we have created three new datasets merged by platform family and performed imputation using a common reference panel, the 1,000 Genomes Phase I release. Here, we describe the methodology behind the data merging and imputation and present imputation quality statistics and association results from two GWAS of secondary phenotypes (body mass index (BMI) and venous thromboembolism (VTE)). We observed the strongest BMI association for the FTO SNP rs55872725 (β = 0.45, p = 3.48x10-22), and using a significance level of p = 0.05, we replicated 19 out of 32 known BMI SNPs. For VTE, we observed the strongest association for the rs2040445 SNP (OR = 2.17, 95% CI: 1.79–2.63, p = 2.70x10-15), located downstream of F5 and also observed significant associations for the known ABO and F11 regions. This pooled resource can be used to maximize power in GWAS of phenotypes collected across the cohorts and for studying gene-environment interactions as well as rare phenotypes and genotypes.