# Pubertal development and prostate cancer risk: Mendelian randomization study in a population-based cohort

Carolina Bonilla, Sarah J Lewis, Richard M Martin, Jenny L Donovan, Freddie C Hamdy, David E Neal, Rosalind Eeles, Doug Easton, Zsofia Kote-Jarai, Ali Amin Al Olama, et al.

► To cite this version:

Carolina Bonilla, Sarah J Lewis, Richard M Martin, Jenny L Donovan, Freddie C Hamdy, et al. Pubertal development and prostate cancer risk: Mendelian randomization study in a population-based cohort. BMC Medicine, 2016, 14 (1), pp.66. 10.1186/s12916-016-0602-x. cea-04551279

HAL Id: cea-04551279

https://cea.hal.science/cea-04551279

Submitted on 18 Apr 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_1.png" /></figure>

<figure><img src="image_2.png" /></figure>

# Pubertal development and prostate cancer risk: Mendelian randomization study in a population-based cohort

Carolina Bonilla<fnref n="1" />,<fnref n="2" />, Sarah J. Lewis<fnref n="1" />,<fnref n="2" />, Richard M. Martin<fnref n="1" />,<fnref n="2" />,³*, Jenny L. Donovan<fnref n="1" />, Freddie C. Hamdy⁴, David E. Neal⁴,⁵, Rosalind Eeles⁶,⁷, Doug Easton⁸, Zsofia Kote-Jarai⁶, Ali Amin Al Olama⁸, Sara Benlloch⁸, Kenneth Muir⁹,¹⁰, Graham G. Giles¹¹,¹², Fredrik Wiklund¹³, Henrik Gronberg¹³, Christopher A. Haiman¹⁴, Johanna Schleutker¹⁵,¹⁶, Børge G. Nordestgaard¹⁷, Ruth C. Travis¹⁸, Nora Pashayan¹⁹,²⁰, Kay-Tee Khaw²¹, Janet L. Stanford²²,²³, William J. Blot²⁴, Stephen Thibodeau²⁵, Christiane Maier²⁶,²⁷, Adam S. Kibel²⁸,²⁹, Cezary Cybulski³⁰, Lisa Cannon-Albright³¹, Hermann Brenner³²,³³,³⁴, Jong Park³⁵, Radka Kaneva³⁶, Jyotsna Batra³⁷, Manuel R. Teixeira³⁸,³⁹, Hardev Pandha⁴⁰, Mark Lathrop⁴¹,⁴², George Davey Smith<fnref n="1" />,<fnref n="2" /> and The PRACTICAL consortium

**Abstract**

**Background:** Epidemiological studies have observed a positive association between an earlier age at sexual development and prostate cancer, but markers of sexual maturation in boys are imprecise and observational estimates are likely to suffer from a degree of uncontrolled confounding. To obtain causal estimates, we examined the role of pubertal development in prostate cancer using genetic polymorphisms associated with Tanner stage in adolescent boys in a Mendelian randomization (MR) approach.

**Methods:** We derived a weighted genetic risk score for pubertal development, combining 13 SNPs associated with male Tanner stage. A higher score indicated a later puberty onset. We examined the association of this score with prostate cancer risk, stage and grade in the UK-based ProtecT case-control study (n = 2,927), and used the PRACTICAL consortium (n = 43,737) as a replication sample.

**Results:** In ProtecT, the puberty genetic score was inversely associated with prostate cancer grade (odds ratio (OR) of high- vs. low-grade cancer, per tertile of the score: 0.76; 95% CI, 0.64–0.89). In an instrumental variable estimation of the causal OR, later physical development in adolescence (equivalent to a difference of one Tanner stage between pubertal boys of the same age) was associated with a 77% (95% CI, 43–91%) reduced odds of high Gleason prostate cancer. In PRACTICAL, the puberty genetic score was associated with prostate cancer stage (OR of advanced vs. localized cancer, per tertile: 0.95; 95% CI, 0.91–1.00) and prostate cancer-specific mortality (hazard ratio amongst cases, per tertile: 0.94; 95% CI, 0.90–0.98), but not with disease grade.

**Conclusions:** Older age at sexual maturation is causally linked to a reduced risk of later prostate cancer, especially aggressive disease.

**Keywords:** Boys, Mendelian randomization, Prostate cancer, Puberty, Tanner scale

<footnote>* Correspondence: richard.martin@bristol.ac.uk</footnote>

<footnote n="1">¹School of Social and Community Medicine, University of Bristol, Bristol, UK</footnote>

<footnote n="2">²MRC Integrative Epidemiology Unit at the University of Bristol, Bristol, UK</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Pubertal development and prostate cancer risk: Mendelian randomization study in a population-based cohort

Carolina Bonilla<fnref n="1" />,<fnref n="2" />, Sarah J. Lewis<fnref n="1" />,<fnref n="2" />, Richard M. Martin<fnref n="1" />,<fnref n="2" />,³*, Jenny L. Donovan<fnref n="1" />, Freddie C. Hamdy⁴, David E. Neal⁴,⁵, Rosalind Eeles⁶,⁷, Doug Easton⁸, Zsofia Kote-Jarai⁶, Ali Amin Al Olama⁸, Sara Benlloch⁸, Kenneth Muir⁹,¹⁰, Graham G. Giles¹¹,¹², Fredrik Wiklund¹³, Henrik Gronberg¹³, Christopher A. Haiman¹⁴, Johanna Schleutker¹⁵,¹⁶, Børge G. Nordestgaard¹⁷, Ruth C. Travis¹⁸, Nora Pashayan¹⁹,²⁰, Kay-Tee Khaw²¹, Janet L. Stanford²²,²³, William J. Blot²⁴, Stephen Thibodeau²⁵, Christiane Maier²⁶,²⁷, Adam S. Kibel²⁸,²⁹, Cezary Cybulski³⁰, Lisa Cannon-Albright³¹, Hermann Brenner³²,³³,³⁴, Jong Park³⁵, Radka Kaneva³⁶, Jyotsna Batra³⁷, Manuel R. Teixeira³⁸,³⁹, Hardev Pandha⁴⁰, Mark Lathrop⁴¹,⁴², George Davey Smith<fnref n="1" />,<fnref n="2" /> and The PRACTICAL consortium

**Abstract**

**Background:** Epidemiological studies have observed a positive association between an earlier age at sexual development and prostate cancer, but markers of sexual maturation in boys are imprecise and observational estimates are likely to suffer from a degree of uncontrolled confounding. To obtain causal estimates, we examined the role of pubertal development in prostate cancer using genetic polymorphisms associated with Tanner stage in adolescent boys in a Mendelian randomization (MR) approach.

**Methods:** We derived a weighted genetic risk score for pubertal development, combining 13 SNPs associated with male Tanner stage. A higher score indicated a later puberty onset. We examined the association of this score with prostate cancer risk, stage and grade in the UK-based ProtecT case-control study (n = 2,927), and used the PRACTICAL consortium (n = 43,737) as a replication sample.

**Results:** In ProtecT, the puberty genetic score was inversely associated with prostate cancer grade (odds ratio (OR) of high- vs. low-grade cancer, per tertile of the score: 0.76; 95% CI, 0.64–0.89). In an instrumental variable estimation of the causal OR, later physical development in adolescence (equivalent to a difference of one Tanner stage between pubertal boys of the same age) was associated with a 77% (95% CI, 43–91%) reduced odds of high Gleason prostate cancer. In PRACTICAL, the puberty genetic score was associated with prostate cancer stage (OR of advanced vs. localized cancer, per tertile: 0.95; 95% CI, 0.91–1.00) and prostate cancer-specific mortality (hazard ratio amongst cases, per tertile: 0.94; 95% CI, 0.90–0.98), but not with disease grade.

**Conclusions:** Older age at sexual maturation is causally linked to a reduced risk of later prostate cancer, especially aggressive disease.

**Keywords:** Boys, Mendelian randomization, Prostate cancer, Puberty, Tanner scale

<footnote>* Correspondence: richard.martin@bristol.ac.uk</footnote>

<footnote n="1">¹School of Social and Community Medicine, University of Bristol, Bristol, UK</footnote>

<footnote n="2">²MRC Integrative Epidemiology Unit at the University of Bristol, Bristol, UK</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
## Background

Prostate cancer is now the most frequently detected cancer among men in Westernized countries [1]. Prostatic intraepithelial neoplasia, a precursor of prostate cancer, has been observed among men in their 20s, suggesting that early-life exposures may play a role in the development of prostate cancer [2] and provide novel opportunities for prostate cancer prevention [3].

Circulating hormones, which rise during puberty, in particular androgens and insulin-like growth factors (IGFs), may play a role in prostate cancer initiation and progression [4, 5], although the relevance of serum androgen levels has recently been called into question [6]. Age at onset of puberty may be a risk factor for prostate adenocarcinoma in men given that exposure to high levels of hormones takes place during the critical window of prostate development in adolescence [3]. Age of menarche is a well-known risk factor for breast cancer [7], but it is yet unclear whether sexual maturation similarly influences later life cancer events in men. However, timing of puberty in boys is difficult to measure as it is not defined by a specific event as in women (menarche); thus, assessing it as a risk factor for prostate cancer in men is challenging.

We investigated whether pubertal development influences risk of prostate cancer in a population-based cohort. We used a genetic score comprised of single nucleotide polymorphisms (SNPs) associated with Tanner genital stage in adolescent boys [8, 9], as a surrogate for the onset and progression of pubertal changes, and we determined associations of this genetic score with prostate cancer risk, stage and grade. The Tanner scale is a widely used 5-point scale that rates breast development in girls, genital development in boys, and pubic hair growth in both [10]. Using a genetic score instead of directly assessed Tanner stage, in an approach known as Mendelian randomization (MR) [11], allows stronger causal inferences because genetic variants are usually unaffected by non-genetic confounding, reverse causality, or measurement error, which underlie the problematic interpretation of observational studies [11, 12].

## Methods

### Subjects

This is a case-control study nested within a multicenter randomized controlled trial of treatments for prostate-specific antigen (PSA)-detected prostate cancer: the Prostate Testing for cancer and Treatment (ProtecT) study (ISRCTN20141297) [13]. During recruitment to the ProtecT study between 2001 and 2009, over 100,000 men aged 50–69 years at 337 general practices in nine UK centres (Birmingham, Bristol, Cambridge, Cardiff, Edinburgh, Leeds, Leicester, Newcastle, Sheffield) were offered a PSA test at a community-based 'prostate check

clinic', and those with raised levels (≥3 ng/mL) were offered a diagnostic biopsy [14]. Detected tumours were all histologically confirmed and clinically staged using the TNM system [15]. In the current analysis, cancer stages T1-T2 were categorized as 'localized'; and T3-T4 as 'locally advanced', because few tumors had metastasized. Histologic material obtained at biopsy was assigned a Gleason score by specialist uropathologists following a standard proforma and, for the purposes of this study, categorized as low- (score ≤6) or high-grade (score ≥7) cancers. All men without evidence of prostate cancer were eligible for selection as controls; that is, men with a PSA <3 ng/mL or a raised PSA (≥3 ng/mL) combined with at least one negative biopsy and no subsequent prostate cancer diagnosis during the follow-up protocol. We selected one stratum-matched control for each case from those men who had provided a non-fasting blood sample at the prostate check clinic. Controls were randomly selected from the same stratum, i.e. 5-year age-band (age at PSA test) and GP/family practice, as cases.

The working dataset consisted of 2,927 individuals (1,136 cases, 1,791 controls) of European descent with available genotype and phenotype information. All men provided written informed consent prior to inclusion in the study. Trent Multicentre Research Ethics Committee (MREC) approved the ProtecT study (MREC/01/4/025) and the associated ProMPT study which collected biological material (MREC/01/4/061; see Additional file 1: Supplementary Methods for further details).

### Genetic risk score

We derived a genetic risk score for pubertal development in boys based on associations between 13 SNPs and Tanner genital stage in males between 12.6 and 15 years of age described in two recent genome-wide association studies (GWAS) of sexual maturation [8, 9]. All SNPs in the score were associated with Tanner stage in boys (independently of whether they were also associated with Tanner stage in girls or in a combined sample of boys and girls), and they had previously been associated with age at menarche [8, 16], although not always in the direction consistent with their association with Tanner genital stage [9]. Variants in the same gene were included in the score provided their linkage disequilibrium, r², was lower than 0.8.

Scores are used instead of individual genetic variants because they are likely to explain a larger proportion of trait variability and therefore represent stronger proxies for the exposure [12]. Scores were calculated by summing up the dosages of the risk alleles at all 13 SNPs in each individual, weighted by the effect size of the variant in males as reported in the discovery GWAS [8, 9], in such a way that a unit increase in the score corresponded
## Background

Prostate cancer is now the most frequently detected cancer among men in Westernized countries [1]. Prostatic intraepithelial neoplasia, a precursor of prostate cancer, has been observed among men in their 20s, suggesting that early-life exposures may play a role in the development of prostate cancer [2] and provide novel opportunities for prostate cancer prevention [3].

Circulating hormones, which rise during puberty, in particular androgens and insulin-like growth factors (IGFs), may play a role in prostate cancer initiation and progression [4, 5], although the relevance of serum androgen levels has recently been called into question [6]. Age at onset of puberty may be a risk factor for prostate adenocarcinoma in men given that exposure to high levels of hormones takes place during the critical window of prostate development in adolescence [3]. Age of menarche is a well-known risk factor for breast cancer [7], but it is yet unclear whether sexual maturation similarly influences later life cancer events in men. However, timing of puberty in boys is difficult to measure as it is not defined by a specific event as in women (menarche); thus, assessing it as a risk factor for prostate cancer in men is challenging.

We investigated whether pubertal development influences risk of prostate cancer in a population-based cohort. We used a genetic score comprised of single nucleotide polymorphisms (SNPs) associated with Tanner genital stage in adolescent boys [8, 9], as a surrogate for the onset and progression of pubertal changes, and we determined associations of this genetic score with prostate cancer risk, stage and grade. The Tanner scale is a widely used 5-point scale that rates breast development in girls, genital development in boys, and pubic hair growth in both [10]. Using a genetic score instead of directly assessed Tanner stage, in an approach known as Mendelian randomization (MR) [11], allows stronger causal inferences because genetic variants are usually unaffected by non-genetic confounding, reverse causality, or measurement error, which underlie the problematic interpretation of observational studies [11, 12].

## Methods

### Subjects

This is a case-control study nested within a multicenter randomized controlled trial of treatments for prostate-specific antigen (PSA)-detected prostate cancer: the Prostate Testing for cancer and Treatment (ProtecT) study (ISRCTN20141297) [13]. During recruitment to the ProtecT study between 2001 and 2009, over 100,000 men aged 50–69 years at 337 general practices in nine UK centres (Birmingham, Bristol, Cambridge, Cardiff, Edinburgh, Leeds, Leicester, Newcastle, Sheffield) were offered a PSA test at a community-based 'prostate check

clinic', and those with raised levels (≥3 ng/mL) were offered a diagnostic biopsy [14]. Detected tumours were all histologically confirmed and clinically staged using the TNM system [15]. In the current analysis, cancer stages T1-T2 were categorized as 'localized'; and T3-T4 as 'locally advanced', because few tumors had metastasized. Histologic material obtained at biopsy was assigned a Gleason score by specialist uropathologists following a standard proforma and, for the purposes of this study, categorized as low- (score ≤6) or high-grade (score ≥7) cancers. All men without evidence of prostate cancer were eligible for selection as controls; that is, men with a PSA <3 ng/mL or a raised PSA (≥3 ng/mL) combined with at least one negative biopsy and no subsequent prostate cancer diagnosis during the follow-up protocol. We selected one stratum-matched control for each case from those men who had provided a non-fasting blood sample at the prostate check clinic. Controls were randomly selected from the same stratum, i.e. 5-year age-band (age at PSA test) and GP/family practice, as cases.

The working dataset consisted of 2,927 individuals (1,136 cases, 1,791 controls) of European descent with available genotype and phenotype information. All men provided written informed consent prior to inclusion in the study. Trent Multicentre Research Ethics Committee (MREC) approved the ProtecT study (MREC/01/4/025) and the associated ProMPT study which collected biological material (MREC/01/4/061; see Additional file 1: Supplementary Methods for further details).

### Genetic risk score

We derived a genetic risk score for pubertal development in boys based on associations between 13 SNPs and Tanner genital stage in males between 12.6 and 15 years of age described in two recent genome-wide association studies (GWAS) of sexual maturation [8, 9]. All SNPs in the score were associated with Tanner stage in boys (independently of whether they were also associated with Tanner stage in girls or in a combined sample of boys and girls), and they had previously been associated with age at menarche [8, 16], although not always in the direction consistent with their association with Tanner genital stage [9]. Variants in the same gene were included in the score provided their linkage disequilibrium, r², was lower than 0.8.

Scores are used instead of individual genetic variants because they are likely to explain a larger proportion of trait variability and therefore represent stronger proxies for the exposure [12]. Scores were calculated by summing up the dosages of the risk alleles at all 13 SNPs in each individual, weighted by the effect size of the variant in males as reported in the discovery GWAS [8, 9], in such a way that a unit increase in the score corresponded
approximately to one risk allele. Risk alleles were those associated with a lower Tanner stage (i.e. delayed pubertal development). Polymorphisms included in the score are shown in Table 1.

## Statistical analysis

Associations of individual SNPs and the multiple SNP score for pubertal development with case/control status and other binary outcomes (localized [reference] vs. locally advanced stage and low [reference] vs. high grade) were determined using logistic regression, with adjustment for age, study center, and the 10 principal components which defined the population structure.

The genetic score was entered into the regression models as a categorical variable with three levels (tertiles). We also used this variable to test for linearity of effect and compare prostate cancer risk among men in the lowest and highest tertiles (i.e. with the earliest and latest sexual maturation, respectively). All analyses were carried out in Stata 13 (StataCorp LP, College Station, TX).

We plotted the effect of each SNP in the genetic score on Tanner stage in approximately 13- to 15-year-old boys against the corresponding effect on high-grade prostate cancer, the disease outcome with the strongest association with the score. The likelihood of bias due to overall directional pleiotropy was formally evaluated with MR-Egger regression [17]. MR-Egger regression also provides an unbiased effect estimate (see definitions in Additional file 1: Supplementary Methods).

We did not have data on Tanner stage measured in adolescence in ProtecT men, and therefore could not estimate its association with the genetic score in ProtecT, in order to run a typical instrumental variable analysis. However, we used a recently developed MR method (summarized data allele score with correlated variants) that provides an estimate of the causal effect of the exposure (i.e. Tanner stage) on the outcome of interest (i.e. prostate cancer) using information on the association of individual SNPs in the score with exposure and outcome [18]. We obtained the effect estimates

<table><caption>**Table 1** SNPs included in the pubertal development genetic risk score in the ProtecT study</caption><thead><tr><th>SNP</th><th>Nearest gene (distance)</th><th>Chr</th><th>Position<sup>a</sup></th><th>Tanner stage decreasing/ other allele</th><th>Tanner stage decreasing allele frequency (ProtecT controls)<sup>b</sup></th><th>Tanner stage decreasing allele frequency (CEU)</th><th>Hardy-Weinberg equilibrium P value</th><th>Gene function</th></tr></thead><tbody><tr><td>rs2274465</td><td>KDM4A</td><td>1</td><td>44121557</td><td>C/G</td><td>0.664</td><td>0.676</td><td>0.9</td><td>Histone demethylation</td></tr><tr><td>rs6427782</td><td>NR5A2 (+198 kb)</td><td>1</td><td>199798339</td><td>A/G</td><td>0.510</td><td>0.556</td><td>0.005</td><td>DNA binding/ steroid hormone receptor activity</td></tr><tr><td>rs6762477</td><td>RBM6</td><td>3</td><td>50093209</td><td>A/G</td><td>0.547</td><td>0.551</td><td>0.6</td><td>Regulation of alternative splicing</td></tr><tr><td>rs2153127</td><td>LIN28B (+36 kb)</td><td>6</td><td>105348544</td><td>T/C</td><td>0.530</td><td>0.515</td><td>0.8</td><td>Cell reprogramming</td></tr><tr><td>rs7759938</td><td>LIN28B (+6 kb)</td><td>6</td><td>105378954</td><td>C/T</td><td>0.318</td><td>0.373</td><td>1.0</td><td>Cell reprogramming</td></tr><tr><td>rs7821178</td><td>PEX2 (-181 kb)</td><td>8</td><td>78093837</td><td>C/A</td><td>0.665</td><td>0.658</td><td>0.1</td><td>Peroxisome biogenesis</td></tr><tr><td>rs10453225</td><td>TMEM38B (-381 kb)</td><td>9</td><td>108920220</td><td>G/T</td><td>0.681</td><td>0.700</td><td>0.2</td><td>Maintenance of intracellular calcium release</td></tr><tr><td>rs2090409</td><td>TMEM38B (-428 kb)</td><td>9</td><td>108967088</td><td>C/A</td><td>0.680</td><td>0.688</td><td>0.3</td><td>Maintenance of intracellular calcium release</td></tr><tr><td>rs10739221</td><td>TMEM38B (-522 kb)</td><td>9</td><td>109060830</td><td>C/T</td><td>0.772</td><td>0.770</td><td>0.5</td><td>Maintenance of intracellular calcium release</td></tr><tr><td>rs1324913</td><td>KLF12</td><td>13</td><td>74635588</td><td>T/G</td><td>0.338</td><td>0.312</td><td>0.6</td><td>Transcription factor/ gene expression regulation</td></tr><tr><td>rs12915845</td><td>DET1</td><td>15</td><td>89042467</td><td>C/T</td><td>0.587</td><td>0.582</td><td>0.1</td><td>Development regulation</td></tr><tr><td>rs246185</td><td>MKL2 (-35 kb)</td><td>16</td><td>14395432</td><td>C/T</td><td>0.313</td><td>0.300</td><td>0.6</td><td>Regulation of immediate early genes/muscle genes</td></tr><tr><td>rs12446632</td><td>GPRC5B (-38 kb)</td><td>16</td><td>19935389</td><td>A/G</td><td>0.146</td><td>0.127</td><td>0.8</td><td>Modulation of insulin secretion</td></tr></tbody></table>

<footnote><sup>a</sup>Position based on GRCh37.p13 assembly</footnote>

<footnote><sup>b</sup>n = 1,791</footnote>
approximately to one risk allele. Risk alleles were those associated with a lower Tanner stage (i.e. delayed pubertal development). Polymorphisms included in the score are shown in Table 1.

## Statistical analysis

Associations of individual SNPs and the multiple SNP score for pubertal development with case/control status and other binary outcomes (localized [reference] vs. locally advanced stage and low [reference] vs. high grade) were determined using logistic regression, with adjustment for age, study center, and the 10 principal components which defined the population structure.

The genetic score was entered into the regression models as a categorical variable with three levels (tertiles). We also used this variable to test for linearity of effect and compare prostate cancer risk among men in the lowest and highest tertiles (i.e. with the earliest and latest sexual maturation, respectively). All analyses were carried out in Stata 13 (StataCorp LP, College Station, TX).

We plotted the effect of each SNP in the genetic score on Tanner stage in approximately 13- to 15-year-old boys against the corresponding effect on high-grade prostate cancer, the disease outcome with the strongest association with the score. The likelihood of bias due to overall directional pleiotropy was formally evaluated with MR-Egger regression [17]. MR-Egger regression also provides an unbiased effect estimate (see definitions in Additional file 1: Supplementary Methods).

We did not have data on Tanner stage measured in adolescence in ProtecT men, and therefore could not estimate its association with the genetic score in ProtecT, in order to run a typical instrumental variable analysis. However, we used a recently developed MR method (summarized data allele score with correlated variants) that provides an estimate of the causal effect of the exposure (i.e. Tanner stage) on the outcome of interest (i.e. prostate cancer) using information on the association of individual SNPs in the score with exposure and outcome [18]. We obtained the effect estimates

<table><caption>**Table 1** SNPs included in the pubertal development genetic risk score in the ProtecT study</caption><thead><tr><th>SNP</th><th>Nearest gene (distance)</th><th>Chr</th><th>Position<sup>a</sup></th><th>Tanner stage decreasing/ other allele</th><th>Tanner stage decreasing allele frequency (ProtecT controls)<sup>b</sup></th><th>Tanner stage decreasing allele frequency (CEU)</th><th>Hardy-Weinberg equilibrium P value</th><th>Gene function</th></tr></thead><tbody><tr><td>rs2274465</td><td>KDM4A</td><td>1</td><td>44121557</td><td>C/G</td><td>0.664</td><td>0.676</td><td>0.9</td><td>Histone demethylation</td></tr><tr><td>rs6427782</td><td>NR5A2 (+198 kb)</td><td>1</td><td>199798339</td><td>A/G</td><td>0.510</td><td>0.556</td><td>0.005</td><td>DNA binding/ steroid hormone receptor activity</td></tr><tr><td>rs6762477</td><td>RBM6</td><td>3</td><td>50093209</td><td>A/G</td><td>0.547</td><td>0.551</td><td>0.6</td><td>Regulation of alternative splicing</td></tr><tr><td>rs2153127</td><td>LIN28B (+36 kb)</td><td>6</td><td>105348544</td><td>T/C</td><td>0.530</td><td>0.515</td><td>0.8</td><td>Cell reprogramming</td></tr><tr><td>rs7759938</td><td>LIN28B (+6 kb)</td><td>6</td><td>105378954</td><td>C/T</td><td>0.318</td><td>0.373</td><td>1.0</td><td>Cell reprogramming</td></tr><tr><td>rs7821178</td><td>PEX2 (-181 kb)</td><td>8</td><td>78093837</td><td>C/A</td><td>0.665</td><td>0.658</td><td>0.1</td><td>Peroxisome biogenesis</td></tr><tr><td>rs10453225</td><td>TMEM38B (-381 kb)</td><td>9</td><td>108920220</td><td>G/T</td><td>0.681</td><td>0.700</td><td>0.2</td><td>Maintenance of intracellular calcium release</td></tr><tr><td>rs2090409</td><td>TMEM38B (-428 kb)</td><td>9</td><td>108967088</td><td>C/A</td><td>0.680</td><td>0.688</td><td>0.3</td><td>Maintenance of intracellular calcium release</td></tr><tr><td>rs10739221</td><td>TMEM38B (-522 kb)</td><td>9</td><td>109060830</td><td>C/T</td><td>0.772</td><td>0.770</td><td>0.5</td><td>Maintenance of intracellular calcium release</td></tr><tr><td>rs1324913</td><td>KLF12</td><td>13</td><td>74635588</td><td>T/G</td><td>0.338</td><td>0.312</td><td>0.6</td><td>Transcription factor/ gene expression regulation</td></tr><tr><td>rs12915845</td><td>DET1</td><td>15</td><td>89042467</td><td>C/T</td><td>0.587</td><td>0.582</td><td>0.1</td><td>Development regulation</td></tr><tr><td>rs246185</td><td>MKL2 (-35 kb)</td><td>16</td><td>14395432</td><td>C/T</td><td>0.313</td><td>0.300</td><td>0.6</td><td>Regulation of immediate early genes/muscle genes</td></tr><tr><td>rs12446632</td><td>GPRC5B (-38 kb)</td><td>16</td><td>19935389</td><td>A/G</td><td>0.146</td><td>0.127</td><td>0.8</td><td>Modulation of insulin secretion</td></tr></tbody></table>

<footnote><sup>a</sup>Position based on GRCh37.p13 assembly</footnote>

<footnote><sup>b</sup>n = 1,791</footnote>
<table><caption>**Table 2** Clinical characteristics of prostate cancer cases and controls in the ProtecT study</caption><thead><tr><th></th><th>Cases</th><th>Controls</th><th>n</th><th>P value</th></tr></thead><tbody><tr><td>Total, n</td><td>1,136</td><td>1,791</td><td>2,927</td><td></td></tr><tr><td>Age, years</td><td>62.2 ± 5.1</td><td>61.6 ± 5.2</td><td>2,927</td><td>0.003</td></tr><tr><td>PSA, ng/mL</td><td>8.5 ± 15.4</td><td>1.1 ± 1.3</td><td>2,925</td><td>&lt;0.001</td></tr><tr><td>Gleason grade, % (&lt;7/≥7)</td><td>70.0/30.0</td><td>n/a</td><td>1,135</td><td></td></tr><tr><td>TNM stage, % (localized/advanced)</td><td>88.4/11.6</td><td>n/a</td><td>1,136</td><td></td></tr><tr><td>BPH, % (no/yes)<sup>a</sup></td><td>92.7/7.3</td><td>89.0/11.0</td><td>1,363</td><td>0.02</td></tr><tr><td>BMI, kg/m<sup>2</sup></td><td>27.0 ± 3.7</td><td>27.5 ± 4.1</td><td>1,973</td><td>0.01</td></tr><tr><td>Height, cm</td><td>176.4 ± 7.0</td><td>175.9 ± 6.9</td><td>2,078</td><td>0.2</td></tr><tr><td>Weight, kg</td><td>84.5 ± 13.0</td><td>85.9 ± 14.6</td><td>2,677</td><td>0.02</td></tr><tr><td>leg length, cm</td><td>76.7 ± 4.8</td><td>76.5 ± 4.6</td><td>2,055</td><td>0.3</td></tr><tr><td>Birthweight, g</td><td>3,437.2 ± 744.9</td><td>3,476.1 ± 663.4</td><td>939</td><td>0.4</td></tr><tr><td>Family history, % (no/yes)<sup>b</sup></td><td>92.7/7.3</td><td>95.0/5.0</td><td>2,908</td><td>0.01</td></tr><tr><td>Diabetes, % (no/yes)<sup>c</sup></td><td>92.8/7.2</td><td>91.1/8.9</td><td>1,895</td><td>0.2</td></tr><tr><td>IGF-I, ng/mL</td><td>156.1 ± 55.8</td><td>163.2 ± 57.0</td><td>1,756</td><td>0.01</td></tr><tr><td>IGF-II, ng/mL</td><td>862.1 ± 323.6</td><td>733.7 ± 265.3</td><td>1,720</td><td>&lt;0.001</td></tr><tr><td>IGFBP-2, ng/mL</td><td>731.5 ± 426.0</td><td>726.2 ± 444.7</td><td>1,745</td><td>0.5<sup>d</sup></td></tr><tr><td>IGFBP-3, ng/mL</td><td>4,673.6 ± 1,041.9</td><td>4,370.8 ± 1,055.7</td><td>1,711</td><td>&lt;0.001</td></tr><tr><td>IGF-I:IGFBP-3 molar ratio<sup>e</sup></td><td>0.12 ± 0.04</td><td>0.14 ± 0.06</td><td>1,711</td><td>&lt;0.001<sup>d</sup></td></tr></tbody></table>

Continuous variables: mean ± SD

Two-sided t-tests and χ² tests were used to analyze continuous and categorical variables, respectively

PSA Prostate-specific antigen, BMI Body mass index, BPH Benign prostatic hyperplasia, IGF Insulin-like growth factor, IGFBP Insulin-like growth factor binding protein

<sup>a</sup>n cases = 682, n controls = 681

<sup>b</sup>n cases = 1,131, n controls = 1,777

<sup>c</sup>n cases = 735, n controls = 1,160

<sup>d</sup>P value obtained from analysis of natural log transformed variable

<sup>e</sup>IGF-I:IGFBP-3 molar ratio = 0.13*[[IGF-I]:0.036*[[IGFBP-3]]]

The genetic score was associated with 15-year prostate cancer-specific mortality amongst men with prostate cancer (hazard ratio (HR) per tertile: 0.94; 95 % CI, 0.90–0.98; *P* = 0.01), and marginally with 15-year all-cause mortality (HR per tertile: 0.97; 95 % CI, 0.95–1.00; *P* = 0.04). The proportional hazards assumption was not met (*P* < 0.001), probably due to the fact that up to 5 years post-diagnosis there do not appear to be differences in survival between individuals with different

numbers of risk alleles (i.e. alleles associated with later pubertal timing; Additional file 1: Figure S3).

We estimated the effect of being ranked a unit lower in the Tanner stage (for the same age) on 15-year prostate cancer-specific mortality as HR 0.62 (95 % CI, 0.49–0.78; *P* < 0.001). The corresponding funnel plot and MR-Egger results, which did not uncover evidence of pleiotropy, are shown in Additional file 1: Figure S4.

<table><caption>**Table 3** Pubertal development genetic risk score and prostate cancer risk, stage and grade in the ProtecT study</caption><thead><tr><th>Trait</th><th>n</th><th>OR</th><th>95 % CI</th><th>P value</th></tr></thead><tbody><tr><td>Control/case (0/1)</td><td>2,927</td><td>0.95</td><td>0.87–1.04</td><td>0.3</td></tr><tr><td>Gleason score (0:≤6/1:≥7)<sup>a</sup></td><td>1,135</td><td>0.76</td><td>0.64–0.89</td><td>0.001</td></tr><tr><td>Stage (0:localised/1:locally advanced)</td><td>1,136</td><td>0.80</td><td>0.64–1.01</td><td>0.06</td></tr><tr><td>BPH (0:no/1:yes)</td><td>1,363</td><td>1.11</td><td>0.88–1.40</td><td>0.4</td></tr></tbody></table>

ORs indicate effects per tertile increase in the genetic score, adjusted by age, recruitment centre and 10 principal components for population structure

<sup>a</sup>0 corresponds to the reference category

BPH Benign prostatic hyperplasia, CI Confidence intervals, OR Odds ratio

<table><caption>**Table 4** Odds ratios (ORs) for high- vs low-grade prostate cancer by pubertal development genetic risk score tertiles in the ProtecT study</caption><thead><tr><th>Genetic score tertiles</th><th>OR</th><th>95 % CI</th><th>P value</th></tr></thead><tbody><tr><td>T1</td><td>Reference</td><td></td><td>0.003</td></tr><tr><td>T2</td><td>0.79</td><td>0.58–1.07</td><td></td></tr><tr><td>T3</td><td>0.57</td><td>0.41–0.79</td><td></td></tr></tbody></table>

High-grade prostate cancer = Gleason ≥7

Low-grade prostate cancer = Gleason ≤6

OR adjusted by age, recruitment centre and 10 principal components for population structure

*n* = 1,135

CI Confidence intervals
<table><caption>**Table 2** Clinical characteristics of prostate cancer cases and controls in the ProtecT study</caption><thead><tr><th></th><th>Cases</th><th>Controls</th><th>n</th><th>P value</th></tr></thead><tbody><tr><td>Total, n</td><td>1,136</td><td>1,791</td><td>2,927</td><td></td></tr><tr><td>Age, years</td><td>62.2 ± 5.1</td><td>61.6 ± 5.2</td><td>2,927</td><td>0.003</td></tr><tr><td>PSA, ng/mL</td><td>8.5 ± 15.4</td><td>1.1 ± 1.3</td><td>2,925</td><td>&lt;0.001</td></tr><tr><td>Gleason grade, % (&lt;7/≥7)</td><td>70.0/30.0</td><td>n/a</td><td>1,135</td><td></td></tr><tr><td>TNM stage, % (localized/advanced)</td><td>88.4/11.6</td><td>n/a</td><td>1,136</td><td></td></tr><tr><td>BPH, % (no/yes)<sup>a</sup></td><td>92.7/7.3</td><td>89.0/11.0</td><td>1,363</td><td>0.02</td></tr><tr><td>BMI, kg/m<sup>2</sup></td><td>27.0 ± 3.7</td><td>27.5 ± 4.1</td><td>1,973</td><td>0.01</td></tr><tr><td>Height, cm</td><td>176.4 ± 7.0</td><td>175.9 ± 6.9</td><td>2,078</td><td>0.2</td></tr><tr><td>Weight, kg</td><td>84.5 ± 13.0</td><td>85.9 ± 14.6</td><td>2,677</td><td>0.02</td></tr><tr><td>leg length, cm</td><td>76.7 ± 4.8</td><td>76.5 ± 4.6</td><td>2,055</td><td>0.3</td></tr><tr><td>Birthweight, g</td><td>3,437.2 ± 744.9</td><td>3,476.1 ± 663.4</td><td>939</td><td>0.4</td></tr><tr><td>Family history, % (no/yes)<sup>b</sup></td><td>92.7/7.3</td><td>95.0/5.0</td><td>2,908</td><td>0.01</td></tr><tr><td>Diabetes, % (no/yes)<sup>c</sup></td><td>92.8/7.2</td><td>91.1/8.9</td><td>1,895</td><td>0.2</td></tr><tr><td>IGF-I, ng/mL</td><td>156.1 ± 55.8</td><td>163.2 ± 57.0</td><td>1,756</td><td>0.01</td></tr><tr><td>IGF-II, ng/mL</td><td>862.1 ± 323.6</td><td>733.7 ± 265.3</td><td>1,720</td><td>&lt;0.001</td></tr><tr><td>IGFBP-2, ng/mL</td><td>731.5 ± 426.0</td><td>726.2 ± 444.7</td><td>1,745</td><td>0.5<sup>d</sup></td></tr><tr><td>IGFBP-3, ng/mL</td><td>4,673.6 ± 1,041.9</td><td>4,370.8 ± 1,055.7</td><td>1,711</td><td>&lt;0.001</td></tr><tr><td>IGF-I:IGFBP-3 molar ratio<sup>e</sup></td><td>0.12 ± 0.04</td><td>0.14 ± 0.06</td><td>1,711</td><td>&lt;0.001<sup>d</sup></td></tr></tbody></table>

Continuous variables: mean ± SD

Two-sided t-tests and χ² tests were used to analyze continuous and categorical variables, respectively

PSA Prostate-specific antigen, BMI Body mass index, BPH Benign prostatic hyperplasia, IGF Insulin-like growth factor, IGFBP Insulin-like growth factor binding protein

a n cases = 682, n controls = 681

b n cases = 1,131, n controls = 1,777

c n cases = 735, n controls = 1,160

d P value obtained from analysis of natural log transformed variable

e IGF-I:IGFBP-3 molar ratio = 0.13*[[IGF-I]:0.036*[[IGFBP-3]]]

The genetic score was associated with 15-year prostate cancer-specific mortality amongst men with prostate cancer (hazard ratio (HR) per tertile: 0.94; 95 % CI, 0.90–0.98; *P* = 0.01), and marginally with 15-year all-cause mortality (HR per tertile: 0.97; 95 % CI, 0.95–1.00; *P* = 0.04). The proportional hazards assumption was not met (*P* <0.001), probably due to the fact that up to 5 years post-diagnosis there do not appear to be differences in survival between individuals with different

numbers of risk alleles (i.e. alleles associated with later pubertal timing; Additional file 1: Figure S3).

We estimated the effect of being ranked a unit lower in the Tanner stage (for the same age) on 15-year prostate cancer-specific mortality as HR 0.62 (95 % CI, 0.49–0.78; *P* <0.001). The corresponding funnel plot and MR-Egger results, which did not uncover evidence of pleiotropy, are shown in Additional file 1: Figure S4.

<table><caption>**Table 3** Pubertal development genetic risk score and prostate cancer risk, stage and grade in the ProtecT study</caption><thead><tr><th>Trait</th><th>n</th><th>OR</th><th>95 % CI</th><th>P value</th></tr></thead><tbody><tr><td>Control/case (0/1)</td><td>2,927</td><td>0.95</td><td>0.87–1.04</td><td>0.3</td></tr><tr><td>Gleason score (0:≤6/1:≥7)<sup>a</sup></td><td>1,135</td><td>0.76</td><td>0.64–0.89</td><td>0.001</td></tr><tr><td>Stage (0:localised/1:locally advanced)</td><td>1,136</td><td>0.80</td><td>0.64–1.01</td><td>0.06</td></tr><tr><td>BPH (0:no/1:yes)</td><td>1,363</td><td>1.11</td><td>0.88–1.40</td><td>0.4</td></tr></tbody></table>

ORs indicate effects per tertile increase in the genetic score, adjusted by age, recruitment centre and 10 principal components for population structure

a 0 corresponds to the reference category

BPH Benign prostatic hyperplasia, CI Confidence intervals, OR Odds ratio

<table><caption>**Table 4** Odds ratios (ORs) for high- vs low-grade prostate cancer by pubertal development genetic risk score tertiles in the ProtecT study</caption><thead><tr><th>Genetic score tertiles</th><th>OR</th><th>95 % CI</th><th>P value</th></tr></thead><tbody><tr><td>T1</td><td>Reference</td><td></td><td>0.003</td></tr><tr><td>T2</td><td>0.79</td><td>0.58–1.07</td><td></td></tr><tr><td>T3</td><td>0.57</td><td>0.41–0.79</td><td></td></tr></tbody></table>

High-grade prostate cancer = Gleason ≥7

Low-grade prostate cancer = Gleason ≤6

OR adjusted by age, recruitment centre and 10 principal components for population structure

n = 1,135

CI Confidence intervals
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Effect of Tanner stage change in boys on the risk of developing high-grade prostate cancer. In ProtecT, proportional risk reduction for high-grade prostate cancer (Gleason ≥7) for each SNP plotted against each SNP's absolute effect on lowering Tanner stage. The trend line, set to intercept the axes at the origin, represents the percentage risk reduction for high-grade disease per unit decrease in Tanner stage. Excluding SNP rs6427782, which was out of Hardy-Weinberg equilibrium pre-Bonferroni correction for multiple testing, from the plot did not produce an appreciable change in the results. Tanner genital stage in boys was treated as a quantitative trait on a scale of 1–5, according to the studies where the associated SNPs were first described [8, 9]</figcaption></figure>

# Discussion

## Main findings

In a study of PSA-detected prostate cancer cases and controls we found strong evidence that a genetic score, comprised of SNPs associated with Tanner stage in approximately 13- to 15-year-old boys, was inversely associated with prostate cancer progression. A later pubertal development (expected among those with higher genetic score values) lowered the risk of developing high-grade disease, a possible clinically relevant subtype because of its stronger relationship than low-grade disease to progression.

Replication analysis using 21 prostate cancer studies across Europe, the USA and Australia in the PRACTICAL consortium, uncovered a weak inverse association between prostate cancer risk and stage and the puberty genetic score, with a reduced effect detected on disease grade. However, we found a stronger association of the score with prostate cancer-specific mortality up to 15 years after diagnosis, indicating that (on average) men whose sexual maturation was later than their peers were less likely to die due to the disease than those whose onset of puberty was earlier. This is in agreement with our findings in ProtecT, regarding the association of earlier puberty with high-grade disease, as men with more aggressive cancer tend to have a poorer prognosis [20]. It is possible that differing definitions of low- and high-grade prostate cancer across studies may have prevented the detection of an effect of the puberty score on this phenotype, with mortality being a stronger and more clear-cut marker of an aggressive disease. There

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Pubertal development genetic risk score and prostate cancer risk (top) and stage (bottom) in the PRACTICAL consortium</figcaption></figure>

were also differences between studies in method of disease detection: the cases enrolled in ProtecT were PSA-detected, whereas the majority of men in the PRACTICAL studies were clinically identified. Additionally, in PRACTICAL, there was a wide variation in the proportion of men with a positive family history of prostate cancer, which ranged from as low as 2.4 % in EPIC to 42.4 % in WUGS (conversely, the proportion of men with a positive family history in ProtecT was ~6 %).

## Mechanisms explaining the observed associations

It has been suggested that endogenous androgen and IGF-I hormones may underlie the relationship between puberty timing and prostate cancer risk [4, 5]. The concentrations of these hormones increases markedly during puberty and are likely to be especially influential on the prostate gland as it becomes fully developed at
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Effect of Tanner stage change in boys on the risk of developing high-grade prostate cancer. In ProtecT, proportional risk reduction for high-grade prostate cancer (Gleason ≥7) for each SNP plotted against each SNP's absolute effect on lowering Tanner stage. The trend line, set to intercept the axes at the origin, represents the percentage risk reduction for high-grade disease per unit decrease in Tanner stage. Excluding SNP rs6427782, which was out of Hardy-Weinberg equilibrium pre-Bonferroni correction for multiple testing, from the plot did not produce an appreciable change in the results. Tanner genital stage in boys was treated as a quantitative trait on a scale of 1–5, according to the studies where the associated SNPs were first described [8, 9]</figcaption></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 2** Pubertal development genetic risk score and prostate cancer risk (top) and stage (bottom) in the PRACTICAL consortium</figcaption></figure>

# Discussion

## Main findings

In a study of PSA-detected prostate cancer cases and controls we found strong evidence that a genetic score, comprised of SNPs associated with Tanner stage in approximately 13- to 15-year-old boys, was inversely associated with prostate cancer progression. A later pubertal development (expected among those with higher genetic score values) lowered the risk of developing high-grade disease, a possible clinically relevant subtype because of its stronger relationship than low-grade disease to progression.

Replication analysis using 21 prostate cancer studies across Europe, the USA and Australia in the PRACTICAL consortium, uncovered a weak inverse association between prostate cancer risk and stage and the puberty genetic score, with a reduced effect detected on disease grade. However, we found a stronger association of the score with prostate cancer-specific mortality up to 15 years after diagnosis, indicating that (on average) men whose sexual maturation was later than their peers were less likely to die due to the disease than those whose onset of puberty was earlier. This is in agreement with our findings in ProtecT, regarding the association of earlier puberty with high-grade disease, as men with more aggressive cancer tend to have a poorer prognosis [20]. It is possible that differing definitions of low- and high-grade prostate cancer across studies may have prevented the detection of an effect of the puberty score on this phenotype, with mortality being a stronger and more clear-cut marker of an aggressive disease. There

were also differences between studies in method of disease detection: the cases enrolled in ProtecT were PSA-detected, whereas the majority of men in the PRACTICAL studies were clinically identified. Additionally, in PRACTICAL, there was a wide variation in the proportion of men with a positive family history of prostate cancer, which ranged from as low as 2.4 % in EPIC to 42.4 % in WUGS (conversely, the proportion of men with a positive family history in ProtecT was ~6 %).

## Mechanisms explaining the observed associations

It has been suggested that endogenous androgen and IGF-I hormones may underlie the relationship between puberty timing and prostate cancer risk [4, 5]. The concentrations of these hormones increases markedly during puberty and are likely to be especially influential on the prostate gland as it becomes fully developed at
this time [21]. Because the prostate is still maturing, puberty may be an important biological window at which early life exposures could have long-term effects on the prostate [3].

Androgens play a central role in the etiology of prostate cancer, as prostate cancer is dependent on androgen receptor activation for growth and survival [22]. A delayed pubertal onset may reduce the length of time an individual is exposed to high androgen levels during a period when the prostate is particularly susceptible to carcinogenic exposures [3].

IGF-I is a potent mitogen and inhibitor of apoptosis that mediates growth during childhood and adolescence and, consequently, stimulates carcinogenesis. IGF-I levels increase from birth to a pubertal peak before declining steadily from young adulthood [23, 24]. It is possible that a younger age at the initiation of puberty is linked to an increase in prostate cancer risk by its association with higher IGF-I levels. Serum IGF-I has been positively associated with an earlier pubertal age in an observational study examining adult IGF-I levels, suggesting that higher pre-pubertal IGF-I (if reflected by higher adult IGF-I) may accelerate childhood growth and the start of puberty [24].

In our study, the puberty genetic score was not associated with circulating IGF-I, although the effect appeared to be in the direction anticipated [24], but we observed an inverse association with the IGF-I:IGFBP-3 molar ratio, suggesting that an earlier pubertal development may be influenced by higher levels of bioavailable IGF-I.

The weakly positive association of the genetic score with adult leg length and height agrees with studies that showed that earlier age at puberty was associated with shorter stature, primarily attributable to shorter leg length, in US women [25] and Swedish men [19]. Growth in leg length, an indicator of pre-pubertal living conditions, on the other hand, has been positively associated with IGF-I levels in UK children, particularly boys [26].

Our finding that a younger age at sexual maturation increases the risk of developing high-grade prostate cancer, and of dying due to the disease, strengthens the idea of the existence of trade-offs between reproductive success and health. MR analysis suggests that there may be a causal relationship between early life environments that promote an accelerated onset of puberty under conditions of uncertainty (e.g. in cases of familial stress due to low income, marital conflict or father absence), so as to favour reproduction, and a detrimental effect on health and longevity in the long term [27, 28].

## Strengths and limitations

Studies of puberty in men are problematic because its initiation is not defined as a single event in the way that

menarche is. In addition, in studies of middle aged and elderly men, an attempt to measure puberty is likely to suffer from recall bias. The genetic score represents a more accurate instrument to assess the causality of the association of pubertal development and prostate cancer risk. The association with leg length provides to some extent a validation of the genetic score in the ProtecT population.

A genetic score is unlikely to be associated with non-genetic confounders, which frequently obscure the interpretation of observational data, and this is, in fact, the case in our study with respect to a few measured confounders.

As the genetic score in our study was not associated with age, PSA, BMI, diabetes, or BPH, we believe that its association with prostate cancer does not represent an artefact of detection due to, for instance, men who are seen more frequently by a doctor having an incidental diagnosis of prostate cancer.

One important assumption in MR is that the instrument (i.e. the genetic score) should be associated with the outcome of interest (i.e. prostate cancer) only via the exposure (i.e. pubertal development). Some SNPs in the score have been associated with height (P <0.05, http://www.gwascentral.org/index) while four of them are located near genes (LIN28B and TMEM38B) recently associated with sitting height ratio and found to disproportionately affect leg length [29]. However, this could be an example of mediated pleiotropy (where a single process influences a cascade of events) [30] and as such it does not undermine our findings. Furthermore, a formal test of the assumption of no pleiotropy, implemented using MR-Egger's regression, found no evidence of a violation of this principle.

Replication of our findings as well as uncovering the potential mechanisms through which the timing of puberty might affect the progression of prostate cancer were likely hindered by differences in phenotype ascertainment in PRACTICAL studies.

## Comparison with existing literature

Few studies have examined the role that pubertal development has on the initiation and progression of prostate cancer, in contrast to the more extensive research on age at menarche and breast cancer. This research shows that an earlier age at menarche is reliably associated with greater breast cancer risk [31]. Given the difficulties in defining puberty among males, studies have used a variety of traits as proxies, mainly age at different life events such as shaving initiation [21, 32], first sexual intercourse [33–36], first ejaculation [37], peak height velocity [24], attainment of adult height [38], as well as the well-established Tanner scale [10]. Results from studies that assessed these variables with respect to
this time [21]. Because the prostate is still maturing, puberty may be an important biological window at which early life exposures could have long-term effects on the prostate [3].

Androgens play a central role in the etiology of prostate cancer, as prostate cancer is dependent on androgen receptor activation for growth and survival [22]. A delayed pubertal onset may reduce the length of time an individual is exposed to high androgen levels during a period when the prostate is particularly susceptible to carcinogenic exposures [3].

IGF-I is a potent mitogen and inhibitor of apoptosis that mediates growth during childhood and adolescence and, consequently, stimulates carcinogenesis. IGF-I levels increase from birth to a pubertal peak before declining steadily from young adulthood [23, 24]. It is possible that a younger age at the initiation of puberty is linked to an increase in prostate cancer risk by its association with higher IGF-I levels. Serum IGF-I has been positively associated with an earlier pubertal age in an observational study examining adult IGF-I levels, suggesting that higher pre-pubertal IGF-I (if reflected by higher adult IGF-I) may accelerate childhood growth and the start of puberty [24].

In our study, the puberty genetic score was not associated with circulating IGF-I, although the effect appeared to be in the direction anticipated [24], but we observed an inverse association with the IGF-I:IGFBP-3 molar ratio, suggesting that an earlier pubertal development may be influenced by higher levels of bioavailable IGF-I.

The weakly positive association of the genetic score with adult leg length and height agrees with studies that showed that earlier age at puberty was associated with shorter stature, primarily attributable to shorter leg length, in US women [25] and Swedish men [19]. Growth in leg length, an indicator of pre-pubertal living conditions, on the other hand, has been positively associated with IGF-I levels in UK children, particularly boys [26].

Our finding that a younger age at sexual maturation increases the risk of developing high-grade prostate cancer, and of dying due to the disease, strengthens the idea of the existence of trade-offs between reproductive success and health. MR analysis suggests that there may be a causal relationship between early life environments that promote an accelerated onset of puberty under conditions of uncertainty (e.g. in cases of familial stress due to low income, marital conflict or father absence), so as to favour reproduction, and a detrimental effect on health and longevity in the long term [27, 28].

## Strengths and limitations

Studies of puberty in men are problematic because its initiation is not defined as a single event in the way that

menarche is. In addition, in studies of middle aged and elderly men, an attempt to measure puberty is likely to suffer from recall bias. The genetic score represents a more accurate instrument to assess the causality of the association of pubertal development and prostate cancer risk. The association with leg length provides to some extent a validation of the genetic score in the ProtecT population.

A genetic score is unlikely to be associated with non-genetic confounders, which frequently obscure the interpretation of observational data, and this is, in fact, the case in our study with respect to a few measured confounders.

As the genetic score in our study was not associated with age, PSA, BMI, diabetes, or BPH, we believe that its association with prostate cancer does not represent an artefact of detection due to, for instance, men who are seen more frequently by a doctor having an incidental diagnosis of prostate cancer.

One important assumption in MR is that the instrument (i.e. the genetic score) should be associated with the outcome of interest (i.e. prostate cancer) only via the exposure (i.e. pubertal development). Some SNPs in the score have been associated with height (P <0.05, http://www.gwascentral.org/index) while four of them are located near genes (LIN28B and TMEM38B) recently associated with sitting height ratio and found to disproportionately affect leg length [29]. However, this could be an example of mediated pleiotropy (where a single process influences a cascade of events) [30] and as such it does not undermine our findings. Furthermore, a formal test of the assumption of no pleiotropy, implemented using MR-Egger's regression, found no evidence of a violation of this principle.

Replication of our findings as well as uncovering the potential mechanisms through which the timing of puberty might affect the progression of prostate cancer were likely hindered by differences in phenotype ascertainment in PRACTICAL studies.

## Comparison with existing literature

Few studies have examined the role that pubertal development has on the initiation and progression of prostate cancer, in contrast to the more extensive research on age at menarche and breast cancer. This research shows that an earlier age at menarche is reliably associated with greater breast cancer risk [31]. Given the difficulties in defining puberty among males, studies have used a variety of traits as proxies, mainly age at different life events such as shaving initiation [21, 32], first sexual intercourse [33–36], first ejaculation [37], peak height velocity [24], attainment of adult height [38], as well as the well-established Tanner scale [10]. Results from studies that assessed these variables with respect to
prostate cancer showed that later growth relative to peers [37], height attainment [38], age at first sexual intercourse [33–36], and being older than 16 at first ejaculation [37] were all associated in a protective direction with prostate cancer. Older age at initiation of shaving was also protective but only among African Americans [21]. As far as we know, no observational study to date has used Tanner stage to investigate the relationship between pubertal development and prostate cancer, so a comparison with our MR findings is not possible.

## Conclusions

Using an MR approach, we have found evidence that experiencing a later sexual maturation reduces prostate cancer risk, especially that of aggressive prostate cancer, as well as mortality owing to the disease. The mechanisms that underlie this relationship may involve the androgenic or IGF pathways, but additional MR studies, using specific instruments for these exposures, should be carried out to investigate this further. Although altering pubertal timing is not a viable prostate cancer prevention strategy there is public health value in identifying those individuals who are more likely to have a worse prognosis [3]. On the other hand, if pubertal development is shown to be driven by increased IGF-I levels, then dietary interventions to regulate its course could potentially be considered.

## Availability of data and materials

Summary data is provided in Additional file 2.

## Additional files

Additional file 1: Tables S1–S6. Figures S1–S4. Supplementary Methods. Definitions. Supplementary References. (DOCX 214 kb)

Additional file 2: Summary data of the association of each SNP in the pubertal development genetic risk score with exposure and outcomes in the PRACTICAL consortium. (XLSX 13 kb)

### Abbreviations

BMI: Body mass index; BPH: Benign prostatic hyperplasia; GWAS: Genome-wide association study; IGF: Insulin-like growth factor; IGFBP: Insulin-like growth factor binding protein; LD: Linkage disequilibrium; MR: Mendelian randomization; PSA: Prostate-specific antigen; SNP: Single nucleotide polymorphism.

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

RMM, SJL and GDS developed the hypothesis and secured funding. CB carried out statistical analyses. CB, SJL and RMM wrote the first draft of the paper. ML organized the genome-wide genotyping of participants in the ProtecT study. JLD, FCH and DEN are PIs of the ProtecT (Prostate Testing for Cancer and Treatment) study. The remaining authors are PIs of the studies contributing to the PRACTICAL consortium. All PIs contributed to the conception, design, and management of each study and of the consortium as a whole. All authors critically commented on and approved the final submitted version of the paper.

### Authors’ information

JLD, FCH and DEN are NIHR Senior Investigators.

### Acknowledgments

The authors thank the tremendous contribution of all members of the ProtecT study research group, and especially the following, who were involved in this research (Prasad Bollina, Sue Bonnington, Lynn Bradshaw, James Catto, Debbie Cooper, Michael Davis, Liz Down, Andrew Doble, Alan Doherty, Garrett Durkan, Emma Elliott, David Gillatt, Pippa Herbert, Peter Holding, Joanne Howson, Mandy Jones, Roger Kockelbergh, Howard Kynaston, Teresa Lennon, Norma Lyons, Hing Leung, Malcolm Mason, Hilary Moody, Philip Powell, Alan Paul, Stephen Prescott, Derek Rosario, Patricia O'Sullivan, Pauline Thompson, Sarah Tidball). We thank Gemma Marsden and Luke Marsden, who processed the blood samples at the biorepository, and Rajeev Kumar, data manager. We also would like to thank the Center National de Génotypage, Evry, France, for genotyping the ProtecT samples. We also thank Stephen Burgess for his help with some of the Mendelian randomization analyses, and Theresa Redaniel for assistance with the survival analysis.

The authors are grateful for the provision of the additional epidemiological data by the NHS R&D Directorate supported Prodigal study and the ProMPT (Prostate Mechanisms of Progression and Treatment) collaboration which is supported by the National Cancer Research Institute (NCRI) formed by the Department of Health, the Medical Research Council and Cancer Research UK (G0500966/75466).

The Collaborative Oncological Gene-environment Study (COGS), within which the PRACTICAL consortium was assembled, would not have been possible without the contributions of the following: Per Hall (COGS), Douglas F. Easton, Paul Pharoah, Kyriaki Michailidou, Manjeet K. Bolla, Qin Wang (BCAC), Andrew Berchuck (OCAC), Rosalind A. Eeles, Douglas F. Easton, Ali Amin Al Olama, Zsofia Kote-Jarai, Sara Benlloch (PRACTICAL), Georgia Chenevix-Trench, Antonis Antoniou, Lesley McGuffog, Fergus Couch, Ken Offit (CIMBA), Joe Dennis, Alison M. Dunning, Andrew Lee, Ed Dicks, Craig Luccarini, and the staff of the Centre for Genetic Epidemiology Laboratory; Javier Benitez, Anna Gonzalez-Neira, and the staff of the CNIO genotyping unit; Jacques Simard, Daniel C. Tessier, Francois Bacot, Daniel Vincent, Sylvie LaBoissière, Frederic Robidoux, and the staff of the McGill University and Génome Québec Innovation Centre; Stig E. Bojesen, Sune F. Nielsen, Maren Weischer, Børge G. Nordestgaard, and the staff of the Copenhagen DNA laboratory; and Julie M. Cunningham, Sharon A. Windebank, Christopher A. Hilker, Jeffrey Meyer, and the staff of the Mayo Clinic Genotyping Core Facility.

The funding sources had no role in the design and conduct of the study, collection, management, analysis and interpretation or preparation, review, or approval of the article.

### Funding

This work was supported by the World Cancer Research Fund (2011/419) and Cancer Research UK (C18281/A19169). The Integrative Epidemiology Unit (IEU) is supported by the MRC and the University of Bristol (G0600705, MC_UU_12013/19), and the Integrative Cancer Epidemiology Programme is supported by Cancer Research UK programme grant C18281/A19169. The National Institute for Health Research (NIHR) Bristol Nutrition Biomedical Research Unit is funded by the NIHR and is a partnership between University Hospitals Bristol NHS Foundation Trust and the University of Bristol. The ProtecT study is supported by the UK NIHR Health Technology Assessment (HTA) Programme (HTA 96/20/99; ISRCTN20141297). Funding for PRACTICAL and the iCOGS infrastructure came from: the European Community's Seventh Framework Programme under grant agreement n° 223175 (HEALTH-F2-2009-223175) (COGS), Cancer Research UK (C1287/A10118, C1287/A10710, C12292/A11174, C1281/A12014, C5047/A8384, C5047/A15007, C5047/A10692, C8197/A16565), the National Institutes of Health (CA128978), and Post-Cancer GWAS initiative (1U19 CA148537, 1U19 CA148065 and 1U19 CA148112 – the GAME-ON initiative), the Department of Defence (W81XWH-10-01-0341), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer, Komen Foundation for the Cure, the Breast Cancer Research Foundation, and the Ovarian Cancer Research Fund. We acknowledge support from the NIHR to the Biomedical Research Centre at The Institute of Cancer Research and The Royal Marsden NHS Foundation Trust.
prostate cancer showed that later growth relative to peers [37], height attainment [38], age at first sexual intercourse [33–36], and being older than 16 at first ejaculation [37] were all associated in a protective direction with prostate cancer. Older age at initiation of shaving was also protective but only among African Americans [21]. As far as we know, no observational study to date has used Tanner stage to investigate the relationship between pubertal development and prostate cancer, so a comparison with our MR findings is not possible.

## Conclusions

Using an MR approach, we have found evidence that experiencing a later sexual maturation reduces prostate cancer risk, especially that of aggressive prostate cancer, as well as mortality owing to the disease. The mechanisms that underlie this relationship may involve the androgenic or IGF pathways, but additional MR studies, using specific instruments for these exposures, should be carried out to investigate this further. Although altering pubertal timing is not a viable prostate cancer prevention strategy there is public health value in identifying those individuals who are more likely to have a worse prognosis [3]. On the other hand, if pubertal development is shown to be driven by increased IGF-I levels, then dietary interventions to regulate its course could potentially be considered.

## Availability of data and materials

Summary data is provided in Additional file 2.

## Additional files

**Additional file 1: Tables S1 – S6.** Figures S1–S4. Supplementary Methods. Definitions. Supplementary References. (DOCX 214 kb)

**Additional file 2:** Summary data of the association of each SNP in the pubertal development genetic risk score with exposure and outcomes in the PRACTICAL consortium. (XLSX 13 kb)

### Abbreviations

BMI: Body mass index; BPH: Benign prostatic hyperplasia; GWAS: Genome-wide association study; IGF: Insulin-like growth factor; IGFBP: Insulin-like growth factor binding protein; LD: Linkage disequilibrium; MR: Mendelian randomization; PSA: Prostate-specific antigen; SNP: Single nucleotide polymorphism.

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

RMM, SJL and GDS developed the hypothesis and secured funding. CB carried out statistical analyses. CB, SJL and RMM wrote the first draft of the paper. ML organized the genome-wide genotyping of participants in the ProtecT study. JLD, FCH and DEN are PIs of the ProtecT (Prostate Testing for Cancer and Treatment) study. The remaining authors are PIs of the studies contributing to the PRACTICAL consortium. All PIs contributed to the conception, design, and management of each study and of the consortium as a whole. All authors critically commented on and approved the final submitted version of the paper.

### Authors’ information

JLD, FCH and DEN are NIHR Senior Investigators.

### Acknowledgments

The authors thank the tremendous contribution of all members of the ProtecT study research group, and especially the following, who were involved in this research (Prasad Bollina, Sue Bonnington, Lynn Bradshaw, James Catto, Debbie Cooper, Michael Davis, Liz Down, Andrew Doble, Alan Doherty, Garrett Durkan, Emma Elliott, David Gillatt, Pippa Herbert, Peter Holding, Joanne Howson, Mandy Jones, Roger Kockelbergh, Howard Kynaston, Teresa Lennon, Norma Lyons, Hing Leung, Malcolm Mason, Hilary Moody, Philip Powell, Alan Paul, Stephen Prescott, Derek Rosario, Patricia O'Sullivan, Pauline Thompson, Sarah Tidball). We thank Gemma Marsden and Luke Marsden, who processed the blood samples at the biorepository, and Rajeev Kumar, data manager. We also would like to thank the Center National de Génotypage, Evry, France, for genotyping the ProtecT samples. We also thank Stephen Burgess for his help with some of the Mendelian randomization analyses, and Theresa Redaniel for assistance with the survival analysis.

The authors are grateful for the provision of the additional epidemiological data by the NHS R&D Directorate supported Prodigal study and the ProMPT (Prostate Mechanisms of Progression and Treatment) collaboration which is supported by the National Cancer Research Institute (NCRI) formed by the Department of Health, the Medical Research Council and Cancer Research UK (G0500966/75466).

The Collaborative Oncological Gene-environment Study (COGS), within which the PRACTICAL consortium was assembled, would not have been possible without the contributions of the following: Per Hall (COGS), Douglas F. Easton, Paul Pharoah, Kyriaki Michailidou, Manjeet K. Bolla, Qin Wang (BCAC), Andrew Berchuck (OCAC), Rosalind A. Eeles, Douglas F. Easton, Ali Amin Al Olama, Zsofia Kote-Jarai, Sara Benlloch (PRACTICAL), Georgia Chenevix-Trench, Antonis Antoniou, Lesley McGuffog, Fergus Couch, Ken Offit (CIMBA), Joe Dennis, Alison M. Dunning, Andrew Lee, Ed Dicks, Craig Luccarini, and the staff of the Centre for Genetic Epidemiology Laboratory; Javier Benitez, Anna Gonzalez-Neira, and the staff of the CNIO genotyping unit; Jacques Simard, Daniel C. Tessier, Francois Bacot, Daniel Vincent, Sylvie LaBoissière, Frederic Robidoux, and the staff of the McGill University and Génome Québec Innovation Centre; Stig E. Bojesen, Sune F. Nielsen, Maren Weischer, Børge G. Nordestgaard, and the staff of the Copenhagen DNA laboratory; and Julie M. Cunningham, Sharon A. Windebank, Christopher A. Hilker, Jeffrey Meyer, and the staff of the Mayo Clinic Genotyping Core Facility.

The funding sources had no role in the design and conduct of the study, collection, management, analysis and interpretation or preparation, review, or approval of the article.

### Funding

This work was supported by the World Cancer Research Fund (2011/419) and Cancer Research UK (C18281/A19169). The Integrative Epidemiology Unit (IEU) is supported by the MRC and the University of Bristol (G0600705, MC_UU_12013/19), and the Integrative Cancer Epidemiology Programme is supported by Cancer Research UK programme grant C18281/A19169. The National Institute for Health Research (NIHR) Bristol Nutrition Biomedical Research Unit is funded by the NIHR and is a partnership between University Hospitals Bristol NHS Foundation Trust and the University of Bristol. The ProtecT study is supported by the UK NIHR Health Technology Assessment (HTA) Programme (HTA 96/20/99; ISRCTN20141297). Funding for PRACTICAL and the iCOGS infrastructure came from: the European Community's Seventh Framework Programme under grant agreement n° 223175 (HEALTH-F2-2009-223175) (COGS), Cancer Research UK (C1287/A10118, C1287/A10710, C12292/A11174, C1281/A12014, C5047/A8384, C5047/A15007, C5047/A10692, C8197/A16565), the National Institutes of Health (CA128978), and Post-Cancer GWAS initiative (1U19 CA148537, 1U19 CA148065 and 1U19 CA148112 – the GAME-ON initiative), the Department of Defence (W81XWH-10-01-0341), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer, Komen Foundation for the Cure, the Breast Cancer Research Foundation, and the Ovarian Cancer Research Fund. We acknowledge support from the NIHR to the Biomedical Research Centre at The Institute of Cancer Research and The Royal Marsden NHS Foundation Trust.
## The PRACTICAL CONSORTIUM (in addition to those named in the author list)

Additional members from the consortium are: Margaret Cook¹, Angela Morga², Artitaya Lophatananon³,⁴, Cyril Fisher², Daniel Leongamornlert², Edward J. Saunders², Emma J. Sawyer², Koveela Govindasami², Magorzata Tymrakiewicz², Michelle Guy², Naomi Livni², Rosemary Wilkinson², Sara Jugurnauth-Little², Steve Hazel², Tokhir Dadaev², Melissa C. Southey⁵, Liesel M. Fitzgerald⁶, John Pedersen⁷, John Hopper⁸, Ami Karlsson⁹, Carin Cavalli-Bjoerkmann⁹, Jan-Erik Johansson⁹, Jan Adolfson⁹, Markus Aly⁹,¹⁰, Michael Broms⁹, Paer Stattin⁹, Brian E. Henderson¹¹, Fredrick Schumacher¹¹, Anssi Auvinen¹², Kimmo Taari¹³, Liisa Maaeettaen¹⁴, Paula Kujala¹⁵, Teemu Murtola¹⁶,¹⁷, Teuvo LJ Tammela¹⁷, Csilla Sipeky¹⁸, Martin Andreas Roder¹⁹, Peter Iversen¹⁹, Peter Klarskov²⁰, Sune F. Nielsen²¹,²², Maren Weischer²¹,²², Tim J. Key²³, Hans Wallinder²⁴, Sven Gustafsson²⁴, Angela Cox²⁵, Anne George²⁶, Athene Lane²⁷, Gemma Marsden²⁸, Michael Davis²⁷, Paul Brown²⁷, Paul Pharoah²⁹, Lisa B. Signorello³⁰,³¹, Wei Zheng³², Shannon K. McDonnell³³, Daniel J. Schaid³³, Liang Wang³³, Lori Tillmans³³, Shaun Riska³³, Thomas Schnoeller³⁴, Kathleen Herkommer³⁵, Manuel Luedeke³⁴, Walther Vogel³⁶, Dominika Wokozorczyk³⁷, Jan Lubiski³⁷, Wojciech Kluzniak³⁷, Katja Butterbach³⁸, Christa Stegmaier³⁹, Bernd Holleczek³⁹, Babu Charchariah⁴⁰, Hui-Yi Lim⁴¹, Hyun Park⁴⁰, James Haley⁴⁰, Julio Pow-Sang⁴⁰, Maria Rincon⁴⁰, Selina Radlein⁴⁰, Thomas Sellers⁴⁰, Chavdar Slavov⁴², Aleksandrina Vlahova⁴³, Atanaska Mitkova⁴⁴, Darina Kachakova⁴⁴, Elenko Popov⁴², Svetlana Christova⁴³, Thiomir Dikov⁴³, Vanio Mitev⁴⁴, Allison Eckert⁴⁵, Angus Collins⁴⁵, Glenn Wood⁴⁵, Greg Malone⁴⁵, Judith A. Clements⁴⁵,⁴⁶, Kris Kerr⁴⁵, Megan Turner⁴⁵, Pamela Saunders⁴⁵, Peter Heathcote⁴⁵, Gail Risbridge⁴⁵, Wayne Tiley⁴⁵, Lisa Horvath⁴⁵, Trina Yeadon⁴⁵, Srilakshmi Srinivasan⁴⁶, Leire Moya⁴⁶, Amanda Spurdle⁴⁷, Joana Santos⁴⁸, Carmen Jerónimo⁴⁸, Paula Paulo⁴⁸, Pedro Pinto⁴⁸, Rui Henrique⁴⁸, Sofia Maia⁴⁸, Agnieszka Michael⁴⁹, Andrzej Kierzek⁴⁹, and Huihai Wu⁴⁹.

¹Centre for Cancer Genetic Epidemiology, Department of Public Health and Primary Care, University of Cambridge, Strangeways Laboratory, Worts Causeway, Cambridge CB1 8RN, UK; ²The Institute of Cancer Research, Sutton, UK; ³Institute of Population Health, University of Manchester, Manchester, UK; ⁴Warwick Medical School, University of Warwick, Coventry, UK; ⁵Genetic Epidemiology Laboratory, Department of Pathology, The University of Melbourne, Grattan Street, Parkville, Victoria 3010, Australia; ⁶Cancer Epidemiology Centre, The Cancer Council Victoria, 615 St Kilda Road, Melbourne, Victoria, Australia; ⁷Tissupath Pty Ltd., Melbourne, Victoria 3122, Australia; ⁸Centre for Epidemiology and Biostatistics, Melbourne School of Population and Global Health, The University of Melbourne, Victoria, Australia; ⁹Department of Medical Epidemiology and Biostatistics, Karolinska Institute, Stockholm, Sweden; ¹⁰Department of Clinical Sciences at Danderyds Hospital, Stockholm, Sweden; ¹¹Department of Preventive Medicine, Keck School of Medicine, University of Southern California/Norris Comprehensive Cancer Center, Los Angeles, California, USA; ¹²Department of Epidemiology, School of Health Sciences, University of Tampere, Tampere, Finland; ¹³Department of Urology, Helsinki University Central Hospital and University of Helsinki, Helsinki, Finland; ¹⁴Finnish Cancer Registry, Helsinki, Finland; ¹⁵Fimlab Laboratories, Tampere University Hospital, Tampere, Finland; ¹⁶School of Medicine, University of Tampere, Tampere, Finland; ¹⁷Department of Urology, Tampere University Hospital and Medical School, University of Tampere, Finland; ¹⁸Medical Biochemistry and Genetics, Institute of Biomedicine, University of Turku, Finland; ¹⁹Copenhagen Prostate Cancer Center, Department of Urology, Rigshospitalet, Copenhagen University Hospital, Tagensevej 20, 7521, DK-2200 Copenhagen, Denmark; ²⁰Department of Urology, Herlev and Gentofte Hospital, Copenhagen University Hospital, Herlev Ringvej 75, DK-2730 Herlev, Denmark; ²¹Department of Clinical Biochemistry, Herlev and Gentofte Hospital, Copenhagen University Hospital, Herlev Ringvej 75, DK-2730 Herlev, Denmark; ²²Faculty of Health and Medical Sciences, University of Copenhagen, Denmark; ²³Cancer Epidemiology Unit, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK; ²⁴Department of Epidemiology and Biostatistics, School of Public Health, Imperial College, London, UK; ²⁵CR-UK/YCR Sheffield Cancer Research Centre, University of Sheffield, Sheffield, UK; ²⁶University of Cambridge, Department of Oncology, Box 279, Addenbrooke's Hospital, Hills Road Cambridge CB2 0QQ, UK; ²⁷School of Social and Community Medicine, University of Bristol, Canyng Hall, 39 Whatley Road, Bristol, BS8 2PS, UK; ²⁸Nuffield Department of Surgical Sciences, Faculty of Medical Science, University of Oxford, John Radcliffe Hospital, Oxford, UK; ²⁹Centre for Cancer Genetic Epidemiology, Department of Oncology, University of Cambridge, Strangeways Laboratory, Worts Causeway, Cambridge, UK; ³⁰International Epidemiology Institute, 1555

Research Blvd., Suite 550, Rockville, MD 20850, USA; ³¹Department of Epidemiology, Harvard School of Public Health, 677 Huntington Avenue, Boston, MA 02115, USA; ³²Division of Epidemiology, Department of Medicine, Vanderbilt University Medical Center, 2525 West End Avenue, Suite 800, Nashville, TN 37232, USA; ³³Mayo Clinic, Rochester, MN, USA; ³⁴Department of Urology, University Hospital Ulm, Germany; ³⁵Department of Urology, Klinikum rechts der Isar der Technischen Universitaet Muenchen, Munich, Germany; ³⁶Institute of Human Genetics, University Hospital Ulm, Germany; ³⁷International Hereditary Cancer Center, Department of Genetics and Pathology, Pomeranian Medical University, Szczecin, Poland; ³⁸Division of Clinical Epidemiology and Aging Research, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; ³⁹Saarland Cancer Registry, 66119 Saarbruecken, Germany; ⁴⁰Department of Cancer Epidemiology, Moffitt Cancer Center, 12902 Magnolia Drive, Tampa, FL 33612, USA; ⁴²Department of Urology and Alexandrovsksa University Hospital, Medical University, Sofia, Bulgaria; ⁴³Department of General and Clinical Pathology, Medical University, Sofia, Bulgaria; ⁴⁴Department of Medical Chemistry and Biochemistry, Molecular Medicine Center, Medical University, Sofia, 2 Zdrave Str., 1431 Sofia, Bulgaria; ⁴⁶Australian Prostate Cancer BioResource, Brisbane, Australia; ⁴⁶Australian Prostate Cancer Research Centre-Queensland, Institute of Health and Biomedical Innovation and School of Biomedical Science, Queensland University of Technology, Brisbane, Australia; ⁴⁷Molecular Cancer Epidemiology Laboratory, Queensland Institute of Medical Research, Brisbane, Australia; ⁴⁸Department of Genetics, Portuguese Oncology Institute, Porto, Portugal; ⁴⁹The University of Surrey, Guildford, Surrey, GU2 7XH, UK.

## Author details

¹School of Social and Community Medicine, University of Bristol, Bristol, UK.

²MRC Integrative Epidemiology Unit at the University of Bristol, Bristol, UK.

³National Institute for Health Research, Bristol Biomedical Research Unit in Nutrition, Bristol, UK. ⁴Nuffield Department of Surgery, University of Oxford, Oxford, UK. ⁵Surgical Oncology (Uro-Oncology: S4), University of Cambridge, Box 279, Addenbrooke's Hospital, Hills Road, Cambridge, UK. ⁶The Institute of Cancer Research, 15 Cotswold Road, Sutton SM2 5NG Surrey, UK. ⁷The Royal Marsden NHS Foundation Trust, Fulham and Sutton, London and Surrey, UK.

⁸Centre for Cancer Genetic Epidemiology, Department of Public Health and Primary Care, University of Cambridge, Strangeways Research Laboratory, Worts Causeway, Cambridge, UK. ⁹University of Warwick, Coventry, UK.

¹⁰Institute of Population Health, The University of Manchester, Manchester M13 9PL, UK. ¹¹The Cancer Council Victoria, 615 St Kilda Road, Melbourne, Victoria 3004, Australia. ¹²Centre for Epidemiology and Biostatistics, Melbourne School of Population and Global Health, The University of Melbourne, Melbourne, Victoria 3010, Australia. ¹³Department of Medical Epidemiology and Biostatistics, Karolinska Institute, Stockholm, Sweden.

¹⁴Department of Preventive Medicine, Keck School of Medicine, University of Southern California/Norris Comprehensive Cancer Center, Los Angeles, CA, USA. ¹⁵Department of Medical Biochemistry and Genetics, University of Turku, Turku, Finland. ¹⁶Institute of Biomedical Technology/BioMediTech, University of Tampere and FimLab Laboratories, Tampere, Finland.

¹⁷Department of Clinical Biochemistry, Herlev Hospital, Copenhagen University Hospital, Herlev Ringvej 75, Herlev DK-2730, Denmark. ¹⁸Cancer Epidemiology Unit, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK. ¹⁹Centre for Cancer Genetic Epidemiology, Department of Oncology, University of Cambridge, Strangeways Research Laboratory, Worts Causeway, Cambridge, UK. ²⁰Department of Applied Health Research, University College London, 1-19 Torrington Place, London WC1E 7HB, UK.

²¹Cambridge Institute of Public Health, University of Cambridge, Forvie Site, Robinson Way, Cambridge CB2 0SR, UK. ²²Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA, USA. ²³Department of Epidemiology, School of Public Health, University of Washington, Seattle, WA, USA. ²⁴International Epidemiology Institute, 1455 Research Blvd., Suite 550, Rockville, MD 20850, USA. ²⁵Mayo Clinic, Rochester, MN, USA. ²⁶Department of Urology, University Hospital Ulm, Ulm, Germany. ²⁷Institute of Human Genetics, University Hospital Ulm, Ulm, Germany. ²⁸Brigham and Women's Hospital/Dana-Farber Cancer Institute, 45 Francis Street - ASB II-3, Boston, MA 02115, USA. ²⁹Washington University, St Louis, MO, USA. ³⁰International Hereditary Cancer Center, Department of Genetics and Pathology, Pomeranian Medical University, Szczecin, Poland. ³¹Division of Genetic Epidemiology, Department of Medicine, University of Utah School of Medicine, Salt Lake City, UT, USA. ³²Division of Clinical Epidemiology and
## The PRACTICAL CONSORTIUM (in addition to those named in the author list)

Additional members from the consortium are: Margaret Cook¹, Angela Morga², Artitaya Lophatananon³,⁴, Cyril Fisher², Daniel Leongamornlert², Edward J. Saunders², Emma J. Sawyer², Koveela Govindasami², Magorzata Tymrakiewicz², Michelle Guy², Naomi Livni², Rosemary Wilkinson², Sara Jugurnauth-Little², Steve Hazel², Tokhir Dadaev², Melissa C. Southey⁵, Liesel M. Fitzgerald⁶, John Pedersen⁷, John Hopper⁸, Ami Karlsson⁹, Carin Cavalli-Bjoerkmann⁹, Jan-Erik Johansson⁹, Jan Adolfson⁹, Markus Aly⁹,¹⁰, Michael Broms⁹, Paer Stattin⁹, Brian E. Henderson¹¹, Fredrick Schumacher¹¹, Anssi Auvinen¹², Kimmo Taari¹³, Liisa Maaeettaenen¹⁴, Paula Kujala¹⁵, Teemu Murtola¹⁶,¹⁷, Teuvo LJ Tammela¹⁷, Csilla Sipeky¹⁸, Martin Andreas Roder¹⁹, Peter Iversen¹⁹, Peter Klarskov²⁰, Sune F. Nielsen²¹,²², Maren Weischer²¹,²², Tim J. Key²³, Hans Wallinder²⁴, Sven Gustafsson²⁴, Angela Cox²⁵, Anne George²⁶, Athene Lane²⁷, Gemma Marsden²⁸, Michael Davis²⁷, Paul Brown²⁷, Paul Pharoah²⁹, Lisa B. Signorello³⁰,³¹, Wei Zheng³², Shannon K. McDonnell³³, Daniel J. Schaid³³, Liang Wang³³, Lori Tillmans³³, Shaun Riska³³, Thomas Schnoeller³⁴, Kathleen Herkommer³⁵, Manuel Luedeke³⁴, Walther Vogel³⁶, Dominika Wokozorczyk³⁷, Jan Lubiski³⁷, Wojciech Kluzniak³⁷, Katja Butterbach³⁸, Christa Stegmaier³⁹, Bernd Holleczek³⁹, Babu Chatariah⁴⁰, Hui-Yi Lim⁴¹, Hyun Park⁴⁰, James Haley⁴⁰, Julio Pow-Sang⁴⁰, Maria Rincon⁴⁰, Selina Radlein⁴⁰, Thomas Sellers⁴⁰, Chavdar Slavov⁴², Aleksandrina Vlahova⁴³, Atanaska Mitkova⁴⁴, Darina Kachakova⁴⁴, Elenko Popov⁴², Svetlana Christova⁴³, Thiomir Dikov⁴³, Vanio Mitev⁴⁴, Allison Eckert⁴⁵, Angus Collins⁴⁵, Glenn Wood⁴⁵, Greg Malone⁴⁵, Judith A. Clements⁴⁵,⁴⁶, Kris Kerr⁴⁵, Megan Turner⁴⁵, Pamela Saunders⁴⁵, Peter Heathcote⁴⁵, Gail Risbridge⁴⁵, Wayne Tiley⁴⁵, Lisa Horvath⁴⁵, Trina Yeadon⁴⁵, Srilakshmi Srinivasan⁴⁶, Leire Moya⁴⁶, Amanda Spurdle⁴⁷, Joana Santos⁴⁸, Carmen Jerónimo⁴⁸, Paula Paulo⁴⁸, Pedro Pinto⁴⁸, Rui Henrique⁴⁸, Sofia Maia⁴⁸, Agnieszka Michael⁴⁹, Andrzej Kierzek⁴⁹, and Huihai Wu⁴⁹.

¹Centre for Cancer Genetic Epidemiology, Department of Public Health and Primary Care, University of Cambridge, Strangeways Laboratory, Worts Causeway, Cambridge CB1 8RN, UK; ²The Institute of Cancer Research, Sutton, UK; ³Institute of Population Health, University of Manchester, Manchester, UK; ⁴Warwick Medical School, University of Warwick, Coventry, UK; ⁵Genetic Epidemiology Laboratory, Department of Pathology, The University of Melbourne, Grattan Street, Parkville, Victoria 3010, Australia; ⁶Cancer Epidemiology Centre, The Cancer Council Victoria, 615 St Kilda Road, Melbourne, Victoria, Australia; ⁷Tissupath Pty Ltd., Melbourne, Victoria 3122, Australia; ⁸Centre for Epidemiology and Biostatistics, Melbourne School of Population and Global Health, The University of Melbourne, Victoria, Australia; ⁹Department of Medical Epidemiology and Biostatistics, Karolinska Institute, Stockholm, Sweden; ¹⁰Department of Clinical Sciences at Danderyds Hospital, Stockholm, Sweden; ¹¹Department of Preventive Medicine, Keck School of Medicine, University of Southern California/Norris Comprehensive Cancer Center, Los Angeles, California, USA; ¹²Department of Epidemiology, School of Health Sciences, University of Tampere, Tampere, Finland; ¹³Department of Urology, Helsinki University Central Hospital and University of Helsinki, Helsinki, Finland; ¹⁴Finnish Cancer Registry, Helsinki, Finland; ¹⁵Fimlab Laboratories, Tampere University Hospital, Tampere, Finland; ¹⁶School of Medicine, University of Tampere, Tampere, Finland; ¹⁷Department of Urology, Tampere University Hospital and Medical School, University of Tampere, Finland; ¹⁸Medical Biochemistry and Genetics, Institute of Biomedicine, University of Turku, Finland; ¹⁹Copenhagen Prostate Cancer Center, Department of Urology, Rigshospitalet, Copenhagen University Hospital, Tagensevej 20, 7521, DK-2200 Copenhagen, Denmark; ²⁰Department of Urology, Herlev and Gentofte Hospital, Copenhagen University Hospital, Herlev Ringvej 75, DK-2730 Herlev, Denmark; ²¹Department of Clinical Biochemistry, Herlev and Gentofte Hospital, Copenhagen University Hospital, Herlev Ringvej 75, DK-2730 Herlev, Denmark; ²²Faculty of Health and Medical Sciences, University of Copenhagen, Denmark; ²³Cancer Epidemiology Unit, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK; ²⁴Department of Epidemiology and Biostatistics, School of Public Health, Imperial College, London, UK; ²⁵CR-UK/YCR Sheffield Cancer Research Centre, University of Sheffield, Sheffield, UK; ²⁶University of Cambridge, Department of Oncology, Box 279, Addenbrooke's Hospital, Hills Road Cambridge CB2 0QQ, UK; ²⁷School of Social and Community Medicine, University of Bristol, Canyng Hall, 39 Whatley Road, Bristol, BS8 2PS, UK; ²⁸Nuffield Department of Surgical Sciences, Faculty of Medical Science, University of Oxford, John Radcliffe Hospital, Oxford, UK; ²⁹Centre for Cancer Genetic Epidemiology, Department of Oncology, University of Cambridge, Strangeways Laboratory, Worts Causeway, Cambridge, UK; ³⁰International Epidemiology Institute, 1555

Research Blvd., Suite 550, Rockville, MD 20850, USA; ³¹Department of Epidemiology, Harvard School of Public Health, 677 Huntington Avenue, Boston, MA 02115, USA; ³²Division of Epidemiology, Department of Medicine, Vanderbilt University Medical Center, 2525 West End Avenue, Suite 800, Nashville, TN 37232, USA; ³³Mayo Clinic, Rochester, MN, USA; ³⁴Department of Urology, University Hospital Ulm, Germany; ³⁵Department of Urology, Klinikum rechts der Isar der Technischen Universitaet Muenchen, Munich, Germany; ³⁶Institute of Human Genetics, University Hospital Ulm, Germany; ³⁷International Hereditary Cancer Center, Department of Genetics and Pathology, Pomeranian Medical University, Szczecin, Poland; ³⁸Division of Clinical Epidemiology and Aging Research, German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany; ³⁹Saarland Cancer Registry, 66119 Saarbruecken, Germany; ⁴⁰Department of Cancer Epidemiology, Moffitt Cancer Center, 12902 Magnolia Drive, Tampa, FL 33612, USA; ⁴²Department of Urology and Alexandrovsksa University Hospital, Medical University, Sofia, Bulgaria; ⁴³Department of General and Clinical Pathology, Medical University, Sofia, Bulgaria; ⁴⁴Department of Medical Chemistry and Biochemistry, Molecular Medicine Center, Medical University, Sofia, 2 Zdrave Str., 1431 Sofia, Bulgaria; ⁴⁶Australian Prostate Cancer BioResource, Brisbane, Australia; ⁴⁶Australian Prostate Cancer Research Centre-Queensland, Institute of Health and Biomedical Innovation and School of Biomedical Science, Queensland University of Technology, Brisbane, Australia; ⁴⁷Molecular Cancer Epidemiology Laboratory, Queensland Institute of Medical Research, Brisbane, Australia; ⁴⁸Department of Genetics, Portuguese Oncology Institute, Porto, Portugal; ⁴⁹The University of Surrey, Guildford, Surrey, GU2 7XH, UK.

## Author details

¹School of Social and Community Medicine, University of Bristol, Bristol, UK.
²MRC Integrative Epidemiology Unit at the University of Bristol, Bristol, UK.
³National Institute for Health Research, Bristol Biomedical Research Unit in Nutrition, Bristol, UK.
⁴Nuffield Department of Surgery, University of Oxford, Oxford, UK.
⁵Surgical Oncology (Uro-Oncology: S4), University of Cambridge, Box 279, Addenbrooke's Hospital, Hills Road, Cambridge, UK.
⁶The Institute of Cancer Research, 15 Cotswold Road, Sutton SM2 5NG Surrey, UK.
⁷The Royal Marsden NHS Foundation Trust, Fulham and Sutton, London and Surrey, UK.
⁸Centre for Cancer Genetic Epidemiology, Department of Public Health and Primary Care, University of Cambridge, Strangeways Research Laboratory, Worts Causeway, Cambridge, UK.
⁹University of Warwick, Coventry, UK.
¹⁰Institute of Population Health, The University of Manchester, Manchester M13 9PL, UK.
¹¹The Cancer Council Victoria, 615 St Kilda Road, Melbourne, Victoria 3004, Australia.
¹²Centre for Epidemiology and Biostatistics, Melbourne School of Population and Global Health, The University of Melbourne, Melbourne, Victoria 3010, Australia.
¹³Department of Medical Epidemiology and Biostatistics, Karolinska Institute, Stockholm, Sweden.
¹⁴Department of Preventive Medicine, Keck School of Medicine, University of Southern California/Norris Comprehensive Cancer Center, Los Angeles, CA, USA.
¹⁵Department of Medical Biochemistry and Genetics, University of Turku, Turku, Finland.
¹⁶Institute of Biomedical Technology/BioMediTech, University of Tampere and FimLab Laboratories, Tampere, Finland.
¹⁷Department of Clinical Biochemistry, Herlev Hospital, Copenhagen University Hospital, Herlev Ringvej 75, Herlev DK-2730, Denmark.
¹⁸Cancer Epidemiology Unit, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, UK.
¹⁹Centre for Cancer Genetic Epidemiology, Department of Oncology, University of Cambridge, Strangeways Research Laboratory, Worts Causeway, Cambridge, UK.
²⁰Department of Applied Health Research, University College London, 1-19 Torrington Place, London WC1E 7HB, UK.
²¹Cambridge Institute of Public Health, University of Cambridge, Forvie Site, Robinson Way, Cambridge CB2 0SR, UK.
²²Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA, USA.
²³Department of Epidemiology, School of Public Health, University of Washington, Seattle, WA, USA.
²⁴International Epidemiology Institute, 1455 Research Blvd., Suite 550, Rockville, MD 20850, USA.
²⁵Mayo Clinic, Rochester, MN, USA.
²⁶Department of Urology, University Hospital Ulm, Ulm, Germany.
²⁷Institute of Human Genetics, University Hospital Ulm, Ulm, Germany.
²⁸Brigham and Women's Hospital/Dana-Farber Cancer Institute, 45 Francis Street - ASB II-3, Boston, MA 02115, USA.
²⁹Washington University, St Louis, MO, USA.
³⁰International Hereditary Cancer Center, Department of Genetics and Pathology, Pomeranian Medical University, Szczecin, Poland.
³¹Division of Genetic Epidemiology, Department of Medicine, University of Utah School of Medicine, Salt Lake City, UT, USA.
³²Division of Clinical Epidemiology and
Aging Research, German Cancer Research Center (DKFZ), Heidelberg, Germany. ³³Division of Preventive Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany. ³⁴German Cancer Consortium (DKTK), German Cancer Research Center (DKFZ), Heidelberg, Germany. ³⁵Division of Cancer Prevention and Control, H. Lee Moffitt Cancer Center, 12902 Magnolia Dr., Tampa, FL, USA. ³⁶Molecular Medicine Center and Department of Medical Chemistry and Biochemistry, Medical University-Sofia, 2 Zdrave St., Sofia 1431, Bulgaria. ³⁷Australian Prostate Cancer Research Centre - Qld, Institute of Health and Biomedical Innovation and School of Biomedical Sciences, Queensland University of Technology, Brisbane, Australia.

³⁸Department of Genetics, Portuguese Oncology Institute, Porto, Portugal.

³⁹Biomedical Sciences Institute (ICBAS), Porto University, Porto, Portugal.

⁴⁰The University of Surrey, Guildford, Surrey GU2 7XH, UK. ⁴¹Commissariat à l'Energie Atomique, Center National de Génotypage, Evry, France. ⁴²McGill University-Génome Québec Innovation Centre, Montreal, Canada.

Received: 20 November 2015 Accepted: 16 March 2016

Published online: 04 April 2016

## References

1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61:69–90.

2. Sakr WA, Haas GP, Cassin BF, Pontes JE, Crissman JD. The frequency of carcinoma and intraepithelial neoplasia of the prostate in young male patients. J Urol. 1993;150:379–85.

3. Sutcliffe S, Colditz GA. Prostate cancer: is it time to expand the research focus to early-life exposures? Nat Rev Cancer. 2013;13:208–518.

4. Salonia A, Abdollah F, Capitanio U, Suardi N, Gallina A, Castagna G, Clementi MC, Briganti A, Rigatti P, Montorsi F. Circulating sex steroids and prostate cancer: introducing the time-dependency theory. World J Urol. 2013;31:267–73.

5. Rowlands M-A, Gunnell D, Harris R, Vatten LJ, Holly JMP, Martin RM. Circulating insulin-like growth factor peptides and prostate cancer risk: a systematic review and meta-analysis. Int J Cancer. 2009;124:2416–29.

6. Isbam H, Pinthus JH, Marks LS, Montorsi F, Morales A, Morgentaler A, Schulman C. Testosterone and prostate cancer: revisiting old paradigms. Eur Urol. 2009;56:48–56.

7. García-Closas M, Brinton LA, Lissowska J, Chatterjee N, Peplonska B, Anderson WF, Szeszenia-Dabrowska N, Bardin-Mikolajczak A, Zatonski W, Blair A, Kalaylioglu Z, Rymkiewicz G, Mazepa-Sikora D, Kordek R, Lukaszek S, Sherman ME. Established breast cancer risk factors by clinically important tumour characteristics. Br J Cancer. 2006;95:123–9.

8. Perry JRB, Day F, Elks CE, Sulem P, Thompson DJ, Ferreira T, He C, Chasman DL, Esko T, Thorleifsson G, Albrecht E, Ang WQ, Corre T, Cousminer DL, Feenstra B, Franceschini N, Ganna A, Johnson AD, Kjellqvist S, Lunetta KL, McMahon G, Nolte IM, Paternoster L, Porcu E, Smith A V, Stolk L, Teumer A, Tšemikova N, Tikkanen E, Ulivi S, et al. Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche. Nature. 2014;514:92–7.

9. Cousminer DL, Stergiakouli E, Berry DJ, Ang W, Groen-Blokhuis MM, Körner A, Siitonen N, Ntalla I, Marinelli M, Perry JRB, Kettunen J, Jansen R, Surakka I, Timpson NJ, Ring S, Mcmahon G, Power C, Wang C, Kähönen M, Viikari J, Lehtimäki T, Middeldorp CM, Hulshoff Pol HE, Neef M, Weise S, Pahkala K, Niinikoski H, Zeggini E, Panoutsopoulou K, Bustamante M, et al. Genome-wide association study of sexual maturation in males and females highlights a role for body mass and menarche loci in male puberty. Hum Mol Genet. 2014;23:4452–64.

10. Tanner JM, Whitehouse RH. Clinical longitudinal standards for height, weight, height velocity, weight velocity, and stages of puberty. Arch Dis Child. 1976;51:170–9.

11. Davey Smith G, Ebrahim S. "Mendelian randomization": can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32:1–22.

12. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, Davey Smith G, Sterne JAC. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21:223–42.

13. Lane JA, Hamdy FC, Martin RM, Turner EL, Neal DE, Donovan JL. Latest results from the UK trials evaluating prostate cancer screening and treatment: the CAP and ProtecT studies. Eur J Cancer. 2010;46:3095–101.

14. Lane JA, Donovan JL, Davis M, Walsh E, Dedman D, Down L, Turner EL, Mason MD, Metcalfe C, Peters TJ, Neal DE, Hamdy FC. Active monitoring, radical prostatectomy, or radiotherapy for localised prostate cancer: study design and diagnostic and baseline results of the ProtecT randomised phase 3 trial. Lancet Oncol. 2014;15:1109–18.

15. Ohori M, Wheeler TM, Scardino PT. The New American Joint Committee on Cancer and International Union Against Cancer TNM classification of prostate cancer. Cancer. 1994;74:104–14.

16. Elks CE, Perry JRB, Sulem P, Chasman DL, Franceschini N, He C, Lunetta KL, Visser J a, Byrne EM, Cousminer DL, Gudbjartsson DF, Esko T, Feenstra B, Hottenga J-J, Koller DL, Kutalik Z, Lin P, Mangino M, Marongiu M, McArdle PF, Smith A V, Stolk L, van Wingerden SH, Zhao JH, Albrecht E, Corre T, Ingelsson E, Hayward C, Magnusson PKE, Smith EN, et al. Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies. Nat Genet. 2010;42:1077–85.

17. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.

18. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30:543–52.

19. Lorentzon M, Norjavaara E, Kindblom JM. Pubertal timing predicts leg length and childhood body mass index predicts sitting height in young adult men. J Pediatr. 2011;158:452–7.

20. Albertsen PC, Fryback DG, Storer BE, Kolon TF, Fine J. Long-term survival among men with conservatively treated localized prostate cancer. JAMA. 1995;274:626–31.

21. Habel LA, Van Den Eeden SK, Friedman GD. Body size, age at shaving initiation, and prostate cancer in a large, multiracial cohort. Prostate. 2000;43:136–43.

22. Grossmann M, Cheung AS, Zajac JD. Androgens and prostate cancer; pathogenesis and deprivation therapy. Best Pract Res Clin Endocrinol Metab. 2013;27:603–16.

23. Bidlingmaier M, Friedrich N, Emery RT, Spranger J, Wolthers OD, Roswall J, Körner A, Obermayer-Pietsch B, Hübener C, Dahlgren J, Frystyk J, Pfeiffer AFH, Doering A, Bielohuby M, Wallaschofski H, Arafat AM. Reference intervals for insulin-like growth factor-1 (IGF-1) from birth to senescence: results from a multicenter study using a new automated chemiluminescence IGF-1 immunoassay conforming to recent international recommendations. J Clin Endocrinol Metab. 2014;99:1712–21.

24. Sandhu J, Davey Smith G, Holly J, Cole TJ, Ben-Shlomo Y. Timing of puberty determines serum insulin-like growth factor-I in late adulthood. J Clin Endocrinol Metab. 2006;91:3150–7.

25. McIntyre MH. Adult stature, body proportions and age at menarche in the United States National Health and Nutrition Survey (NHANES) III. Ann Hum Biol. 2011;38:716–20.

26. Rogers I, Metcalfe C, Gunnell D, Emmett P, Dunger D, Holly J. Insulin-like growth factor-I and growth in height, leg length, and trunk length between ages 5 and 10 years. J Clin Endocrinol Metab. 2006;91:2514–9.

27. Belsky J, Steinberg L, Draper P. Childhood experience, interpersonal development, and reproductive strategy: and evolutionary theory of socialization. Child Dev. 1991;62:647–70.

28. Thomas F, Elguero E, Brodeur J, Roche B, Missé D, Raymond M. Malignancies and high birth weight in human: which cancers could result from antagonistic pleiotropy? J Evol Med. 2012;1:1–5.

29. Chan Y, Salem RM, Hsu Y-HH, McMahon G, Pers TH, Vedantam S, Esko T, Guo MH, Lim ET, Franke L, Smith GD, Strachan DP, Hirschhorn JN. Genome-wide analysis of body proportion classifies height-associated variants by mechanism of action and implicates genes important for skeletal development. Am J Hum Genet. 2015;96:695–708.

30. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23:R89–98.

31. Velie EM, Nechuta S, Osuch JR. Lifetime reproductive and anthropometric risk factors for breast cancer in postmenopausal women. Breast Dis. 2006;24:17–35.

32. Barba M, Terrenato I, Schünemann HJ, Fuhrman B, Sperati F, Teter B, Gallucci M, D'Amato A, Muti P. Indicators of sexual and somatic development and adolescent body size in relation to prostate cancer risk: results from a case-control study. Urology. 2008;72:183–7.

33. Andersson SO, Baron J, Bergström R, Lindgren C, Wolk A, Adami HO. Lifestyle factors and prostate cancer risk: a case-control study in Sweden. Cancer Epidemiol Biomarkers Prev. 1996;5:509–13.

34. Honda GD, Bernstein L, Ross RK, Greenland S, Gerkins V, Henderson BE. Vasectomy, cigarette smoking, and age at first sexual intercourse as risk factors for prostate cancer in middle-aged men. Br J Cancer. 1988;57:326–31.

35. Schuman LM, Mandel J, Blackard C, Bauer H, Scarlett J, McHugh R. Epidemiologic study of prostatic cancer: preliminary report. Cancer Treat Rep. 1977;61:181–6.
Aging Research, German Cancer Research Center (DKFZ), Heidelberg, Germany. ³³Division of Preventive Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany. ³⁴German Cancer Consortium (DKTK), German Cancer Research Center (DKFZ), Heidelberg, Germany. ³⁵Division of Cancer Prevention and Control, H. Lee Moffitt Cancer Center, 12902 Magnolia Dr., Tampa, FL, USA. ³⁶Molecular Medicine Center and Department of Medical Chemistry and Biochemistry, Medical University-Sofia, 2 Zdrave St., Sofia 1431, Bulgaria. ³⁷Australian Prostate Cancer Research Centre - Qld, Institute of Health and Biomedical Innovation and School of Biomedical Sciences, Queensland University of Technology, Brisbane, Australia.

³⁸Department of Genetics, Portuguese Oncology Institute, Porto, Portugal.

³⁹Biomedical Sciences Institute (ICBAS), Porto University, Porto, Portugal.

⁴⁰The University of Surrey, Guildford, Surrey GU2 7XH, UK. ⁴¹Commissariat à l'Energie Atomique, Center National de Génotypage, Evry, France. ⁴²McGill University-Génome Québec Innovation Centre, Montreal, Canada.

Received: 20 November 2015 Accepted: 16 March 2016

Published online: 04 April 2016

## References

1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61:69–90.

2. Sakr WA, Haas GP, Cassin BF, Pontes JE, Crissman JD. The frequency of carcinoma and intraepithelial neoplasia of the prostate in young male patients. J Urol. 1993;150:379–85.

3. Sutcliffe S, Colditz GA. Prostate cancer: is it time to expand the research focus to early-life exposures? Nat Rev Cancer. 2013;13:208–518.

4. Salonia A, Abdollah F, Capitanio U, Suardi N, Gallina A, Castagna G, Clementi MC, Briganti A, Rigatti P, Montorsi F. Circulating sex steroids and prostate cancer: introducing the time-dependency theory. World J Urol. 2013;31:267–73.

5. Rowlands M-A, Gunnell D, Harris R, Vatten LJ, Holly JMP, Martin RM. Circulating insulin-like growth factor peptides and prostate cancer risk: a systematic review and meta-analysis. Int J Cancer. 2009;124:2416–29.

6. Isbam H, Pinthus JH, Marks LS, Montorsi F, Morales A, Morgentaler A, Schulman C. Testosterone and prostate cancer: revisiting old paradigms. Eur Urol. 2009;56:48–56.

7. García-Closas M, Brinton LA, Lissowska J, Chatterjee N, Peplonska B, Anderson WF, Szeszenia-Dabrowska N, Bardin-Mikolajczak A, Zatonski W, Blair A, Kalaylioglu Z, Rymkiewicz G, Mazepa-Sikora D, Kordek R, Lukaszek S, Sherman ME. Established breast cancer risk factors by clinically important tumour characteristics. Br J Cancer. 2006;95:123–9.

8. Perry JRB, Day F, Elks CE, Sulem P, Thompson DJ, Ferreira T, He C, Chasman DL, Esko T, Thorleifsson G, Albrecht E, Ang WQ, Corre T, Cousminer DL, Feenstra B, Franceschini N, Ganna A, Johnson AD, Kjellqvist S, Lunetta KL, McMahon G, Nolte IM, Paternoster L, Porcu E, Smith A V, Stolk L, Teumer A, Tšemikova N, Tikkanen E, Ulivi S, et al. Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche. Nature. 2014;514:92–7.

9. Cousminer DL, Stergiakouli E, Berry DJ, Ang W, Groen-Blokhuis MM, Körner A, Siitonen N, Ntalla I, Marinelli M, Perry JRB, Kettunen J, Jansen R, Surakka I, Timpson NJ, Ring S, Mcmahon G, Power C, Wang C, Kähönen M, Viikari J, Lehtimäki T, Middeldorp CM, Hulshoff Pol HE, Neef M, Weise S, Pahkala K, Niinikoski H, Zeggini E, Panoutsopoulou K, Bustamante M, et al. Genome-wide association study of sexual maturation in males and females highlights a role for body mass and menarche loci in male puberty. Hum Mol Genet. 2014;23:4452–64.

10. Tanner JM, Whitehouse RH. Clinical longitudinal standards for height, weight, height velocity, weight velocity, and stages of puberty. Arch Dis Child. 1976;51:170–9.

11. Davey Smith G, Ebrahim S. "Mendelian randomization": can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32:1–22.

12. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, Davey Smith G, Sterne JAC. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21:223–42.

13. Lane JA, Hamdy FC, Martin RM, Turner EL, Neal DE, Donovan JL. Latest results from the UK trials evaluating prostate cancer screening and treatment: the CAP and ProtecT studies. Eur J Cancer. 2010;46:3095–101.

14. Lane JA, Donovan JL, Davis M, Walsh E, Dedman D, Down L, Turner EL, Mason MD, Metcalfe C, Peters TJ, Neal DE, Hamdy FC. Active monitoring, radical prostatectomy, or radiotherapy for localised prostate cancer: study design and diagnostic and baseline results of the ProtecT randomised phase 3 trial. Lancet Oncol. 2014;15:1109–18.

15. Ohori M, Wheeler TM, Scardino PT. The New American Joint Committee on Cancer and International Union Against Cancer TNM classification of prostate cancer. Cancer. 1994;74:104–14.

16. Elks CE, Perry JRB, Sulem P, Chasman DL, Franceschini N, He C, Lunetta KL, Visser J a, Byrne EM, Cousminer DL, Gudbjartsson DF, Esko T, Feenstra B, Hottenga J-J, Koller DL, Kutalik Z, Lin P, Mangino M, Marongiu M, McArdle PF, Smith A V, Stolk L, van Wingerden SH, Zhao JH, Albrecht E, Corre T, Ingelsson E, Hayward C, Magnusson PKE, Smith EN, et al. Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies. Nat Genet. 2010;42:1077–85.

17. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.

18. Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur J Epidemiol. 2015;30:543–52.

19. Lorentzon M, Norjavaara E, Kindblom JM. Pubertal timing predicts leg length and childhood body mass index predicts sitting height in young adult men. J Pediatr. 2011;158:452–7.

20. Albertsen PC, Fryback DG, Storer BE, Kolon TF, Fine J. Long-term survival among men with conservatively treated localized prostate cancer. JAMA. 1995;274:626–31.

21. Habel LA, Van Den Eeden SK, Friedman GD. Body size, age at shaving initiation, and prostate cancer in a large, multiracial cohort. Prostate. 2000;43:136–43.

22. Grossmann M, Cheung AS, Zajac JD. Androgens and prostate cancer; pathogenesis and deprivation therapy. Best Pract Res Clin Endocrinol Metab. 2013;27:603–16.

23. Bidlingmaier M, Friedrich N, Emery RT, Spranger J, Wolthers OD, Roswall J, Körner A, Obermayer-Pietsch B, Hübener C, Dahlgren J, Frystyk J, Pfeiffer AFH, Doering A, Bielohuby M, Wallaschofski H, Arafat AM. Reference intervals for insulin-like growth factor-1 (IGF-1) from birth to senescence: results from a multicenter study using a new automated chemiluminescence IGF-1 immunoassay conforming to recent international recommendations. J Clin Endocrinol Metab. 2014;99:1712–21.

24. Sandhu J, Davey Smith G, Holly J, Cole TJ, Ben-Shlomo Y. Timing of puberty determines serum insulin-like growth factor-I in late adulthood. J Clin Endocrinol Metab. 2006;91:3150–7.

25. McIntyre MH. Adult stature, body proportions and age at menarche in the United States National Health and Nutrition Survey (NHANES) III. Ann Hum Biol. 2011;38:716–20.

26. Rogers I, Metcalfe C, Gunnell D, Emmett P, Dunger D, Holly J. Insulin-like growth factor-I and growth in height, leg length, and trunk length between ages 5 and 10 years. J Clin Endocrinol Metab. 2006;91:2514–9.

27. Belsky J, Steinberg L, Draper P. Childhood experience, interpersonal development, and reproductive strategy: and evolutionary theory of socialization. Child Dev. 1991;62:647–70.

28. Thomas F, Elguero E, Brodeur J, Roche B, Missé D, Raymond M. Malignancies and high birth weight in human: which cancers could result from antagonistic pleiotropy? J Evol Med. 2012;1:1–5.

29. Chan Y, Salem RM, Hsu Y-HH, McMahon G, Pers TH, Vedantam S, Esko T, Guo MH, Lim ET, Franke L, Smith GD, Strachan DP, Hirschhorn JN. Genome-wide analysis of body proportion classifies height-associated variants by mechanism of action and implicates genes important for skeletal development. Am J Hum Genet. 2015;96:695–708.

30. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23:R89–98.

31. Velie EM, Nechuta S, Osuch JR. Lifetime reproductive and anthropometric risk factors for breast cancer in postmenopausal women. Breast Dis. 2006;24:17–35.

32. Barba M, Terrenato I, Schünemann HJ, Fuhrman B, Sperati F, Teter B, Gallucci M, D'Amato A, Muti P. Indicators of sexual and somatic development and adolescent body size in relation to prostate cancer risk: results from a case-control study. Urology. 2008;72:183–7.

33. Andersson SO, Baron J, Bergström R, Lindgren C, Wolk A, Adami HO. Lifestyle factors and prostate cancer risk: a case-control study in Sweden. Cancer Epidemiol Biomarkers Prev. 1996;5:509–13.

34. Honda GD, Bernstein L, Ross RK, Greenland S, Gerkins V, Henderson BE. Vasectomy, cigarette smoking, and age at first sexual intercourse as risk factors for prostate cancer in middle-aged men. Br J Cancer. 1988;57:326–31.

35. Schuman LM, Mandel J, Blackard C, Bauer H, Scarlett J, McHugh R. Epidemiologic study of prostatic cancer: preliminary report. Cancer Treat Rep. 1977;61:181–6.
36. Dennis LK, Dawson DV. Meta-analysis of measures of sexual activity and prostate cancer. Epidemiology. 2002;13:72–9.

37. Giles GG, Severi G, English DR, McCredie MRE, MacInnis R, Boyle P, Hopper JL. Early growth, adult body size and prostate cancer risk. Int J Cancer. 2003;103:241–5.

38. Hayes RB, de Jong FH, Raatgever J, Bogdanovicz J, Schroeder FH, van der Maas P, Oishi K, Yoshida O. Physical characteristics and factors related to sexual development and behaviour and the risk for prostatic cancer. Eur J Cancer Prev. 1992;1:239–45.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central
36. Dennis LK, Dawson DV. Meta-analysis of measures of sexual activity and prostate cancer. Epidemiology. 2002;13:72–9.

37. Giles GG, Severi G, English DR, McCredie MRE, MacInnis R, Boyle P, Hopper JL. Early growth, adult body size and prostate cancer risk. Int J Cancer. 2003;103:241–5.

38. Hayes RB, de Jong FH, Raatgever J, Bogdanovicz J, Schroeder FH, van der Maas P, Oishi K, Yoshida O. Physical characteristics and factors related to sexual development and behaviour and the risk for prostatic cancer. Eur J Cancer Prev. 1992;1:239–45.

Submit your next manuscript to BioMed Central
and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central