Age at puberty and risk of asthma: A Mendelian randomisation study

Abstract

Background

Observational studies on pubertal timing and asthma, mainly performed in females, have provided conflicting results about a possible association of early puberty with higher risk of adult asthma, possibly due to residual confounding. To overcome issues of confounding, we used Mendelian randomisation (MR), i.e., genetic variants were used as instrumental variables to estimate causal effects of early puberty on post-pubertal asthma in both females and males.

Methods and findings

MR analyses were performed in UK Biobank on 243,316 women using 254 genetic variants for age at menarche, and on 192,067 men using 46 variants for age at voice breaking. Age at menarche, recorded in years, was categorised as early (<12), normal (12–14), or late (>14); age at voice breaking was recorded and analysed as early (younger than average), normal (about average age), or late (older than average). In females, we found evidence for a causal effect of pubertal timing on asthma, with an 8% increase in asthma risk for early menarche (odds ratio [OR] 1.08; 95% CI 1.04 to 1.12; p = 8.7 × 10−5) and an 8% decrease for late menarche (OR 0.92; 95% CI 0.89 to 0.97; p = 3.4 × 10−4), suggesting a continuous protective effect of increasing age at puberty. In males, we found very similar estimates of causal effects, although with wider confidence intervals (early voice breaking: OR 1.07; 95% CI 1.00 to 1.16; p = 0.06; late voice breaking: OR 0.93; 95% CI 0.87 to 0.99; p = 0.03). We detected only modest pleiotropy, and our findings showed robustness when different methods to account for pleiotropy were applied. BMI may either introduce pleiotropy or lie on the causal pathway; secondary analyses excluding variants associated with BMI yielded similar results to those of the main analyses. Our study relies on self-reported exposures and outcomes, which may have particularly affected the power of the analyses on age at voice breaking.

Conclusions

This large MR study provides evidence for a causal detrimental effect of early puberty on asthma, and does not support previous observational findings of a U-shaped relationship between pubertal timing and asthma. Common biological or psychological mechanisms associated with early puberty might explain the similarity of our results in females and males, but further research is needed to investigate this. Taken together with evidence for other detrimental effects of early puberty on health, our study emphasises the need to further investigate and address the causes of the secular shift towards earlier puberty observed worldwide.

About the Authors:

Cosetta Minelli

Roles Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Population Health and Occupational Disease, National Heart and Lung Institute, Imperial College London, London, United Kingdom

ORCID logo http://orcid.org/0000-0001-9166-3958

Diana A. van der Plaat

Roles Formal analysis, Writing – original draft, Writing – review & editing

Affiliations UMR 1152, INSERM, Paris, France, UMR 1152, Université Paris Diderot, Paris, France

ORCID logo http://orcid.org/0000-0002-0361-357X

Bénédicte Leynaert

Roles Methodology, Writing – original draft, Writing – review & editing

Affiliations UMR 1152, INSERM, Paris, France, UMR 1152, Université Paris Diderot, Paris, France

Raquel Granell

Roles Investigation, Writing – original draft, Writing – review & editing

Affiliations MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom, Population Health Sciences, University of Bristol, Bristol, United Kingdom

Andre F. S. Amaral

Roles Investigation, Writing – original draft, Writing – review & editing

Affiliation: Population Health and Occupational Disease, National Heart and Lung Institute, Imperial College London, London, United Kingdom

ORCID logo http://orcid.org/0000-0002-0369-9449

Miguel Pereira

Roles Formal analysis, Writing – review & editing

Affiliation: Population Health and Occupational Disease, National Heart and Lung Institute, Imperial College London, London, United Kingdom

Osama Mahmoud

Roles Investigation, Writing – review & editing

Affiliation: Population Health Sciences, University of Bristol, Bristol, United Kingdom

ORCID logo http://orcid.org/0000-0003-0342-6704

James Potts

Roles Formal analysis, Writing – review & editing

Affiliation: Population Health and Occupational Disease, National Heart and Lung Institute, Imperial College London, London, United Kingdom

Nuala A. Sheehan

Roles Methodology, Writing – review & editing

Affiliation: Department of Health Sciences, University of Leicester, Leicester, United Kingdom

Jack Bowden

Roles Methodology, Writing – review & editing

Affiliations MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom, Population Health Sciences, University of Bristol, Bristol, United Kingdom

John Thompson

Roles Methodology, Writing – review & editing

Affiliation: Department of Health Sciences, University of Leicester, Leicester, United Kingdom

Debbie Jarvis

Roles Conceptualization, Writing – original draft, Writing – review & editing

Affiliation: Population Health and Occupational Disease, National Heart and Lung Institute, Imperial College London, London, United Kingdom

ORCID logo http://orcid.org/0000-0002-1753-3896

George Davey Smith

Roles Conceptualization, Methodology, Writing – review & editing

Affiliations MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom, Population Health Sciences, University of Bristol, Bristol, United Kingdom

ORCID logo http://orcid.org/0000-0002-1407-8314

John Henderson

Roles Conceptualization, Writing – original draft, Writing – review & editing

Affiliations MRC Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom, Population Health Sciences, University of Bristol, Bristol, United Kingdom

Introduction

Asthma is one of the most common chronic diseases in children and adults, and its prevalence has been increasing worldwide [1]. Multiple risk factors for asthma have been proposed, but one of the few consistent observations is sex differences in disease prevalence. Understanding the reason for these differences could lead to the identification of modifiable factors that could be targeted by preventive measures [2].

Sex differences in asthma prevalence are age dependent, with a switch from male predominance in childhood to female predominance from adolescence onwards [3,4]. It has been hypothesized that the higher prevalence of post-pubertal asthma in females, with higher incidence and lower remission rates [4], might be due to a pathogenetic role of female sex hormones [3,5]. This is further supported by observational findings of an association of early menarche with increased risk of post-pubertal asthma [6]. In addition to the hypothesis of sex hormones mediating the effect of early menarche on asthma, another proposed mechanism is through obesity: Early menarche would increase the risk of adult obesity, which would in turn affect asthma, possibly through an effect of higher leptin levels on the immune system and airway inflammation, although it could be the other way round, with higher childhood BMI and leptin increasing the risk of early menarche, which would increase the risk of asthma [7,8]. Despite the plausibility that the observed association of early menarche with asthma may represent a causal effect, this cannot be concluded based on the available evidence, and, if there is a causal effect, it is unclear what its magnitude may be. A meta-analysis of 7 observational studies suggested a 37% increase in asthma risk for women with menarche before the age of 12 years, but there was substantial heterogeneity in results across studies [9]. The meta-analysis was mainly based on young study populations, and an observational study in UK Biobank including women aged 40 to 69 years suggested a smaller increase of 6% in asthma risk after adjustment for confounders [10]. The same study also showed an association in males, with an 11% increase in risk associated with early voice breaking, a marker of male puberty, and it suggested a U-shaped relationship in females and males, with an increased risk of asthma associated with both early and late puberty [10]. More recently, a prospective study in women up to the age of 60 years suggested no association after extensive adjustment for confounders [11]. The possibility of residual confounding by unmeasured factors, including intrauterine exposures as well as childhood diet, lifestyle and socioeconomic status, limits the interpretation of these findings and may explain inconsistencies across studies. There is also a potential problem of reverse causation, as there are some data suggesting that girls with childhood asthma might have earlier menarche [12]. This could generate spurious associations, particularly if study participants with pre-pubertal asthma were not excluded from the analyses due to poor recall of their childhood symptoms.

The Mendelian randomisation (MR) approach applied here uses genetic variants known to affect age at puberty as proxies (instrumental variables) to derive estimates of the effect of age at puberty on asthma [13]. As genes are randomly allocated at the time of conception, genetic associations are not affected by typical confounding factors or reverse causation, and MR can provide indirect evidence for a causal effect if its underlying assumptions hold [13,14]. Here we describe a large MR study in UK Biobank aimed at estimating the causal effect of age at puberty on the risk of post-pubertal asthma in females and males.

Methods

The UK Biobank is a prospective study of about 500,000 adults aged 40–69 years recruited in 2006–2010 across 22 centres, aimed at identifying causes of chronic disease in middle and old age [15]. About 95% of the participants are white. Data on doctor-diagnosed asthma and age at onset (AAO), age at menarche in females, and age at voice breaking in males are all based on self-reports from questionnaires. The specific questions used to define the variables used are reported in Table 1, together with the UK Biobank data-field number for access to further information on measurement procedures through the UK Biobank webpage (http://biobank.ctsu.ox.ac.uk/crystal/search.cgi). While age at menarche was recorded in years, age at voice breaking was recorded as a 3-level categorical variable (at a younger, average, or older age compared to peers). In order to increase comparability with the findings in males, we categorised age at menarche in females as early (<12 years), normal (12 to 14), or late (>14), using “normal” as the reference. This also allowed the investigation of possible non-linear associations, such as U-shaped relationships. In females, we defined post-pubertal asthma as asthma occurring after menarche; in males, this was defined as asthma occurring after the age of 15 years. Participants with pre-pubertal asthma were excluded from the analyses.

[Figure omitted. See PDF.]

Table 1. Characteristics of the study population, for females and males.

https://doi.org/10.1371/journal.pmed.1002634.t001

Choice of instruments for MR analyses

We chose MR instruments for age at menarche based on findings from the large genome-wide association (GWA) meta-analysis by Day et al. [16], which included 329,345 women of European ancestry from the ReproGen consortium (40 studies, N = 179,117), 23andMe (N = 76,831), and UK Biobank (N = 73,397). In all, 389 independent single nucleotide polymorphisms (SNPs) were identified at a genome-wide significant level (p < 5 × 10−8), explaining overall 7.4% of the population variance in age at menarche. Of the reported 389 SNPs, 372 were available (or had a proxy, linkage disequilibrium r2 > 0.7) in UK Biobank.

The instrument strength of each SNP for both early and late menarche was assessed using the F statistic, a function of the magnitude and precision of the SNP’s effect on menarche [17]. To avoid bias associated with the use of weak instruments, we excluded SNPs with a low F statistic in UK Biobank, using the common threshold of 10 [18]. This left us with 206 and 151 SNPs (254 in total) as instruments for early and late menarche, respectively; F statistic values varied from 10.1 to 382.5 for early menarche and from 10.0 to 150.0 for late menarche (S1 Table). All SNPs had good imputation quality (info score > 0.8).

In the GWA meta-analysis by Day et al. [16], 127 of the 389 age-at-menarche SNPs were also associated in the same direction with age at voice breaking at nominal significance level; of these, only 40 were strong instruments (F statistic > 10) for either early or late puberty in males. In order to increase our number of instruments, we searched the literature for additional genetic studies on age at voice breaking, using the NHGRI GWAS Catalog [19] and HuGE Navigator [20], as well as crosschecking of references. We identified 3 other GWA studies [21–23], 2 of which added new independent SNPs, the studies by Day et al. [22] and by Pickrell et al. [23] (S2 Table). After combining all SNPs and excluding weak instruments, 24 and 37 SNPs (46 in total) remained for the MR analyses of early and late voice breaking, respectively. Their F statistic values varied from 10.2 to 52.4 for early voice breaking and from 10.8 to 101.6 for late voice breaking (S2 Table). All SNPs had good imputation quality (info score > 0.9). All 46 SNPs were independent, except for 3 pairs of correlated SNPs (linkage disequilibrium r2 of 0.3, 0.6, and 0.6). As a post hoc analysis to investigate the impact of the correlation between these SNPs on our MR results, we used the method by Burgess et al., which accounts for the correlation between instruments [24].

Details on genotyping and imputation methods and quality control procedures in UK Biobank are available elsewhere [25].

MR methods

In 1-sample MR investigations (where the genetic associations with both risk factor and outcome are estimated within the same study), the standard MR approach is the 2-stage least squares (2SLS) method: the risk factor is regressed jointly on all SNPs, and the outcome is regressed on the genetically predicted values of the risk factor from the first regression [26]. As this method does not account for pleiotropy (where SNPs chosen as instruments for the risk factor affect the outcome through additional independent pathways [14]) and we cannot rule out the possibility of pleiotropic effects for some of our SNPs, we used methods developed for 2-sample MR. Here estimates of the genetic associations with risk factor and outcome are obtained separately for each SNP. If all SNPs are valid instruments (e.g., no pleiotropy) and linear models with no interactions faithfully describe the SNP–risk factor and risk factor–outcome relationships [27], the individual MR estimates of causal effect will vary only by chance, with no between-instrument heterogeneity [28]. For each SNP, the association with early or late menarche/voice breaking (G–X, expressed as log odds ratio [OR]) and the association with asthma (G–Y, log OR) were estimated using multiple logistic regression adjusted for assessment centre, first 10 ancestry principal components, and genotyping array. It should be noted that with binary variables the linearity assumption is strictly violated [29], but in practice the magnitude of the bias has been shown to be negligible [30].

An MR estimate was obtained for each SNP using a Wald estimator (ratio of G–Y over G–X), with standard error derived using the Delta method [31], and individual MR estimates were pooled using inverse-variance weighted (IVW) fixed-effect meta-analysis [32]. The pooled MR estimate is asymptotically equal to the MR estimate obtained with the 2SLS method [32], which we also performed as a sensitivity analysis. Since the IVW fixed-effect meta-analysis method assumes no pleiotropy, we performed secondary analyses using approaches based on different assumptions about pleiotropy across SNPs: (1) IVW random-effect meta-analysis, assuming random pleiotropic effects (i.e., they cancel out) [33]; (2) MR-Egger regression with penalised weights, assuming directional pleiotropic effects [34]; and (3) weighted median analysis, assuming that at least 50% of the SNPs are valid instruments (but with no assumption about direction of pleiotropic effects) [35].

Some of the genetic variants associated with pubertal timing are also associated with BMI [16], and this may complicate the interpretation of the results given that BMI could either introduce pleiotropy (i.e., affecting asthma independently from age at puberty) or lie on the causal pathway. In secondary analyses we therefore excluded SNPs associated with BMI from our instruments. We identified SNPs associated with BMI, overweight, or obesity in previous GWA studies from PhenoScanner, a database of publicly available GWA findings (http://www.phenoscanner.medschl.cam.ac.uk/phenoscanner) [36], using a significance threshold corrected for multiple testing, p < 2.0 × 10−4 (0.05/254) for menarche SNPs and p < 1.1 × 10−3 (0.05/46) for voice breaking SNPs.

To provide further support for the absence of pleiotropy, we used pre-pubertal asthma as a “negative control” outcome [37]. If our instruments affect asthma solely through age at menarche/voice breaking, we expect to find no effect on asthma developed before puberty.

There was no formal pre-specified protocol for this study; all the analyses described above were decided a priori, and we performed additional post hoc analyses suggested by the reviewers. In particular, sensitivity analyses were carried out to assess the robustness of the results to the use of different AAO cutoffs for post-pubertal asthma definition in males, and to the presence of some correlation between genetic variants. Subgroup analyses were also performed to investigate the possibility that the effect of age at puberty on asthma might be modified by whether or not one is overweight.

All analyses were performed using Stata 15 (StataCorp).

Results

The characteristics of the study population are summarised in Table 1. We included 243,316 women and 192,067 men with available genetic data and information on age at puberty, asthma and AAO; post-pubertal asthma was present in 22,295 (9.2%) women and 12,066 (6.3%) men. Asthma prevalence by category of age at puberty was 10.4% for early, 8.9% for normal, and 8.8% for late menarche in females, and 7.5% for early, 6.2% for normal, and 7.1% for late age at voice breaking in males.

MR analyses

Detailed results from all MR analyses are reported in S3 and S4 Tables, for females and males, respectively. In females, we found evidence for a detrimental causal effect of early menarche, with an OR of 1.08 (95% CI 1.04 to 1.12; p = 8.7 × 10−5), and for a protective causal effect of late menarche, with an OR of 0.92 (95% CI 0.89 to 0.97; p = 3.4 × 10−4). This suggests a continuous protective effect of increasing age at puberty, and does not support observational evidence of a U-shaped relationship. When we repeated the analysis using age at menarche as a continuous variable (332 SNPs with F > 10 as instruments), we found an OR of 1.06 (95% CI 1.03 to 1.11; p = 1.6 × 10−4) per year decrease in age at menarche. In males, we found an OR of 1.07 (95% CI 1.00 to 1.16; p = 0.064) for early and 0.93 (95% CI 0.87 to 0.99; p = 0.029) for late voice breaking. The results remained the same in sensitivity analyses allowing for the correlation between SNPs, and became slightly stronger in sensitivity analyses using higher age cutoffs for asthma AAO to define post-pubertal asthma, specifically AAO >16 and >17 years (S4 Table). These findings in males are very similar to those in females, despite the much wider confidence intervals reflecting the lower power of the analyses (fewer instruments available for age at voice breaking and possibly larger measurement error). The use of a 1-sample 2SLS approach gave similar results in both females and males (S3 and S4 Tables).

We found only modest between-instrument heterogeneity in females (I2 of 16% and 20% for early and late puberty, respectively; S3 Table) and little evidence of heterogeneity in males (I2 of 2% and 0%, respectively; S4 Table). The results of the analyses controlling for pleiotropy were consistent with those of the main analyses (Fig 1), with differences in the precision of the estimates in line with expected differences in statistical power [33,35]. MR-Egger regression in males suggested stronger effects compared to the main analysis, particularly for early puberty (Fig 1), but the confidence intervals were wide as the analysis had limited power. This low power is due to the low number of SNPs and their relatively similar instrument strengths (MR-Egger regression works best when there is a large spread of strengths [38]).

[Figure omitted. See PDF.]

Fig 1. Estimates across different MR methods of the effects of early and late puberty on asthma risk in females (menarche) and males (voice breaking).

In red: main analysis. IVW, inverse-variance weighted; MR, Mendelian randomisation.

https://doi.org/10.1371/journal.pmed.1002634.g001

Secondary analyses with exclusion of SNPs associated with BMI (14 and 12 excluded for early and late menarche; 5 and 7 excluded for early and late voice breaking) gave similar results to those of the main analyses for both females and males (S3 and S4 Tables).

As the percentages of females in the early and late menarche categories were substantially higher than the percentages of males in the corresponding categories (Table 1), we repeated the analyses using more stringent cutoffs for age at menarche. Using a cutoff of <11 instead of <12 years for early and >15 instead of >14 years for late menarche, the percentages were very similar to those in males, 4.6% and 5.8% for early and late menarche, respectively. The results were in line with those of the main analyses, with an OR of 1.06 (95% CI 1.02 to 1.09; p = 6.4 × 10−4) for early and 0.95 (95% CI 0.91 to 0.98; p = 0.002) for late menarche.

The MR analyses on pre-pubertal asthma, used as a negative control outcome, suggested no effect for either age at menarche or age at voice breaking, further reassuring against effects of our genetic instruments on asthma independent from pubertal timing.

Finally, we repeated the analyses separately by BMI group, with low/high BMI defined as below/above the median (<26/≥26 kg/m2 in females; <27/≥27 kg/m2 in males), to investigate possible effect modification due to whether or not participants are overweight. BMI data refer to the participants’ BMI at the time of recruitment into UK Biobank (Table 1). In females, these analyses (S3 Table) showed effect estimates for early and late menarche consistent with the main analyses in the high BMI group, and no evidence for effects of early or late menarche in the low BMI group, but there was no statistical evidence of interaction (p-values of 0.23 and 0.67 for tests of the difference in estimates between BMI groups for early and late menarche, respectively). In males, there was no evidence of effect modification (S4 Table), as expected given the very low statistical power.

Discussion

This MR study provides evidence for a detrimental causal effect of earlier age at puberty on the risk of developing asthma later in life. In females, we estimated an 8% increase in asthma risk for early menarche (<12 years) and an 8% decrease for late menarche (>14 years). A very similar pattern was observed in males, with a 7% increase in risk for early puberty and a 7% decrease for late puberty, although the evidence was weaker than in females due to the much lower statistical power. Previous observational studies on pubertal timing and asthma, mainly performed in females, provided conflicting results, possibly due to residual confounding, reverse causation, or inappropriate adjustment for factors that lie on the causal pathway. This study adds to current knowledge by providing evidence for a causal effect of early puberty on asthma in both females and males, and by clarifying the pattern of the relationship between pubertal timing and asthma. Previous observational analyses in UK Biobank [10] suggested that there could be a U-shaped relationship (both early and late puberty would be detrimental), but our MR study does not support this and suggests a continuous protective effect of increasing age at puberty. The discrepancy may be due to confounding of the observational findings by factors associated with both delayed puberty and increased risk of asthma.

The results of our post hoc analysis by BMI group are compatible with the possibility of a stronger effect of age at puberty on asthma in overweight females, but these results have major limitations that make it difficult to draw conclusions. Stratifying by BMI might introduce collider bias, whereby conditioning on a common effect of exposure and outcome generates a spurious association between the two [39]; BMI is indeed affected by some of the puberty SNPs and might also be affected by asthma through medication side effects or reduced physical exercise. Moreover, the power of these analyses in females is low, in particular to detect interaction, and even more so in males. Finally, in our analyses we could only use adult BMI; although childhood obesity is likely to track to obesity in adulthood, the inability to disentangle the two makes these analyses difficult to interpret. Further research using longitudinal data is needed to investigate a possible role of obesity in modifying the effect of age at puberty on asthma.

The similarity of effect estimates in females and males provides insight into possible underlying mechanisms. Previously reported associations of early menarche with asthma have been hypothesized to be explained by earlier and greater cumulative exposure to female sex hormones such as oestrogens or an unbalanced oestrogen/progesterone ratio [5,9], although experimental evidence in mice does not support this explanation [40]. A role for female sex hormones has also been suggested by a reported association between pregnancy, characterised by increased levels of progesterone and oestrogens, and asthma development [41], and by decreased incidence of asthma after natural menopause [42,43]. The similarity of our results in both sexes suggests that early menarche may affect asthma through mechanisms that go beyond female sex hormones. Other factors accompanying early puberty in both females and males could be responsible, including biological but also psychological factors, such as depression, anxiety, and psychosomatic symptoms [44]. Further investigation of these aspects is important to identify modifiable factors that mediate the effects of early puberty on asthma and other adverse health outcomes.

The effect estimates found in our large study are of modest magnitude and therefore of limited clinical relevance to an individual, but they have important implications at a population level. In Europe, age at menarche decreased by 2 to 3 months per decade between 1790 and 1980 [45] and continued to decrease at a lower rate in the following decades [46], with similar trends observed worldwide. Although the reasons for such a shift are not fully understood, there is evidence suggesting a role for changing childhood risk factors, including diet and obesity, psychological stress, and environmental exposures [47]. Childhood obesity is likely the strongest potentially preventable factor, and currently represents a global epidemic, with half of the world’s population predicted to be overweight or obese by 2030 [48]. Other environmental risk factors for early puberty warrant further investigation, including endocrine-disrupting chemicals found in many household products that are hypothesized to have transgenerational epigenetic effects [49]. Observational evidence, particularly in females, has suggested detrimental effects of early puberty on other health outcomes including lung function, BMI, cardiometabolic traits, cardiovascular disease, and cancer [10,50–55], with recent support from MR analyses [16,22,56,57]. Taken together with this evidence, our findings emphasise the importance of investigating modifiable risk factors for early puberty with the aim of reversing the secular shift towards earlier puberty observed worldwide.

Strengths and weaknesses

MR is a valuable tool to assess causality in epidemiology as it is not affected by the classical confounding factors or reverse causation typical of observational studies. However, this technique is at its most reliable and transparent in the absence of pleiotropy, i.e., in this case, where the genetic instruments do not have direct effects on asthma independent of age at puberty [14]. As the risk of pleiotropy is high in MR investigations where many SNPs of uncertain biological function are included, we thoroughly investigated its potential impact on our findings. Obesity was identified a priori as a potential source of pleiotropic effects given its genetic overlap with pubertal timing, and in fact around 7% of our instruments for age at puberty were also associated with BMI. In practice, what role BMI might play in our study is unclear, as it could induce “horizontal” pleiotropy by affecting asthma independently from puberty, but it could also lie on the causal pathway either preceding puberty (high childhood BMI → early puberty → high asthma risk) or mediating its effects (early puberty → high adult BMI → high asthma risk), which is referred to as “vertical pleiotropy” [58]. Only horizontal pleiotropy would affect the validity of our MR findings. When excluding SNPs associated with BMI, we found very similar results to those of the main analyses. Pleiotropy could act through mechanisms other than obesity. We therefore investigated the likely magnitude of pleiotropy by assessing the heterogeneity among causal estimates from the individual SNPs (no heterogeneity if all valid instruments [28]), and found evidence of only modest heterogeneity. This was supported by the consistency in the findings from 3 different methods to adjust for pleiotropy (IVW random-effect, weighted median, and MR-Egger regression) with those of the main analyses. Finally, we repeated our analyses using pre-pubertal asthma as a negative control outcome, and the null results provided further reassurance against pleiotropic mechanisms independent from pubertal timing.

We used MR methods developed to investigate and adjust for pleiotropy in 2-sample MR, where estimates of the association of each SNP with risk factor and outcome are derived from separate studies, and the 2 estimates are assumed independent. The behaviour of these methods has yet to be studied in the 1-sample context, where independence of the 2 estimates is violated. However, the consistency of our findings with those derived from the 1-sample 2SLS method provides reassurance.

Our study is based on self-reported information on doctor-diagnosed asthma, age at menarche, and age at voice breaking, which introduces measurement error and potentially recall bias. However, since it is unlikely that the error in recalling age at puberty would correlate with the reporting of asthma or asthma AAO, and vice versa, the measurement error is likely to be random and only result in reduced power. Voice breaking represents the drop in the pitch of speech due to changes in vocal chords and larynx cartilage in response to androgen exposure during late puberty [59]. Despite having been suggested as a good marker of pubertal timing in males when assessed prospectively [60], age at voice breaking is much harder to recall in adult cross-sectional studies compared with the more clear-cut milestone of menarche in females, and may be subject to greater measurement error. This is a general problem and partly explains the limited knowledge about risk factors and effects of early puberty in males, for which more evidence from longitudinal birth cohort studies is needed. Another limitation is that we based the choice of MR instruments for age at menarche and, partly, age at voice breaking on findings from the GWA meta-analysis by Day et al. [16], of which UK Biobank genetic data (first release only) represented about 20% of the total sample. However, any overestimation of the SNP–age at puberty associations resulting from this would have pulled our MR estimates towards the null, therefore leading to underestimation of the true causal effects rather than to false positive results. Finally, findings from observational studies might be biased by “cohort effects” if individuals from younger generations are at different risk of both early (or late) puberty and asthma compared with older generations. Despite the secular trends of decreasing age at puberty [45] and increasing prevalence of asthma [1], this type of confounding would not bias the results of our MR study, unless the secular trends also affect the genetic variants used as instruments for age at puberty, which is unlikely.

In conclusion, we provide evidence for a causal detrimental effect of early puberty on asthma. Further studies are needed to investigate the mechanisms underlying this effect, and the similarity of our findings in females and males suggests a possible role for common biological or psychological factors accompanying early pubertal development. Taken together with evidence for other detrimental effects of early puberty on health, our study emphasises the need to investigate and address causes of the secular shift towards earlier puberty observed worldwide.

Supporting information

S1 Table. Data for the MR analyses of age at menarche in females.

https://doi.org/10.1371/journal.pmed.1002634.s001

(XLSX)

S2 Table. Data for the MR analyses of age at voice breaking in males.

https://doi.org/10.1371/journal.pmed.1002634.s002

(XLSX)

S3 Table. Results of main and secondary MR analyses of age at menarche in females.

https://doi.org/10.1371/journal.pmed.1002634.s003

(XLSX)

S4 Table. Results of main and secondary MR analyses of age at voice breaking in males.

https://doi.org/10.1371/journal.pmed.1002634.s004

(XLSX)

Acknowledgments

This research was conducted using the UK Biobank resource under application number 19136, and we thank the participants, field workers, and data managers for their time and cooperation.

Citation: Minelli C, van der Plaat DA, Leynaert B, Granell R, Amaral AFS, Pereira M, et al. (2018) Age at puberty and risk of asthma: A Mendelian randomisation study. PLoS Med 15(8): e1002634. https://doi.org/10.1371/journal.pmed.1002634

References

1. Eder W, Ege MJ, von Mutius E. The asthma epidemic. N Engl J Med. 2006;355(21):2226–35. pmid:17124020

2. Beasley R, Semprini A, Mitchell EA. Risk factors for asthma: is prevention possible? Lancet. 2015;386(9998):1075–85. pmid:26382999

3. Almqvist C, Worm M, Leynaert B. Impact of gender on asthma in childhood and adolescence: a GA2LEN review. Allergy. 2008;63(1):47–57. pmid:17822448

4. Vink NM, Postma DS, Schouten JP, Rosmalen JG, Boezen HM. Gender differences in asthma development and remission during transition through puberty: the TRacking Adolescents’ Individual Lives Survey (TRAILS) study. J Allergy Clin Immunol. 2010;126(3):498–504.e1–6. pmid:20816186

5. Zein JG, Erzurum SC. Asthma is different in women. Curr Allergy Asthma Rep. 2015;15(6):28. pmid:26141573

6. McCleary N, Nwaru BI, Nurmatov UB, Critchley H, Sheikh A. Endogenous and exogenous sex steroid hormones in asthma and allergy in females: a systematic review and meta-analysis. J Allergy Clin Immunol. 2018;141(4):1510–3.e8. pmid:29305316

7. Macsali F, Svanes C, Bjorge L, Omenaas ER, Gomez Real F. Respiratory health in women: from menarche to menopause. Exp Rev Resp Med. 2012;6(2):187–200. pmid:22455491

8. Mumby HS, Elks CE, Li S, Sharp SJ, Khaw KT, Luben RN, et al. Mendelian randomisation study of childhood BMI and early menarche. J Obes. 2011;2011:180729. pmid:21773002

9. Lieberoth S, Gade EJ, Brok J, Backer V, Thomsen SF. Age at menarche and risk of asthma: systematic review and meta-analysis. J Asthma. 2014;51(6):559–65. pmid:24628525

10. Day FR, Elks CE, Murray A, Ong KK, Perry JR. Puberty timing associated with diabetes, cardiovascular disease and also diverse health outcomes in men and women: the UK Biobank study. Sci Rep. 2015;5:11208. pmid:26084728

11. Hansen S, Probst-Hensch N, Bettschart R, Pons M, Leynaert B, Gomez Real F, et al. Early menarche and new onset of asthma: results from the SAPALDIA cohort study. Maturitas. 2017;101:57–63. pmid:28539170

12. Drosdzol A, Skrzypulec V, Wilk K, Nowosielski K. The influence of bronchial asthma on sexual maturation of girls. J Physiol Pharmacol. 2007;58(Suppl 5(Pt 1)):155–63. pmid:18204126

13. Davey Smith G, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22. pmid:12689998

14. Sheehan NA, Didelez V, Burton PR, Tobin MD. Mendelian randomisation and causal inference in observational epidemiology. PLoS Med. 2008;5(8):e177. pmid:18752343

15. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12(3):e1001779. pmid:25826379

16. Day FR, Thompson DJ, Helgason H, Chasman DI, Finucane H, Sulem P, et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat Genet. 2017;49(6):834–41. pmid:28436984

17. Palmer TM, Lawlor DA, Harbord RM, Sheehan NA, Tobias JH, Timpson NJ, et al. Using multiple genetic variants as instrumental variables for modifiable risk factors. Stat Methods Med Res. 2012;21(3):223–42. pmid:21216802

18. Lawlor DA, Harbord RM, Sterne JA, Timpson N, Davey Smith G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med. 2008;27(8):1133–63. pmid:17886233

19. Welter D, MacArthur J, Morales J, Burdett T, Hall P, Junkins H, et al. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014;42(Database issue):D1001–6. pmid:24316577

20. Yu W, Gwinn M, Clyne M, Yesupriya A, Khoury MJ. A navigator for human genome epidemiology. Nat Genet. 2008;40(2):124–5. pmid:18227866

21. Ong KK, Elks CE, Li S, Zhao JH, Luan J, Andersen LB, et al. Genetic variation in LIN28B is associated with the timing of puberty. Nat Genet. 2009;41(6):729–33. pmid:19448623

22. Day FR, Bulik-Sullivan B, Hinds DA, Finucane HK, Murabito JM, Tung JY, et al. Shared genetic aetiology of puberty timing between sexes and with health-related outcomes. Nat Commun. 2015;6:8842. pmid:26548314

23. Pickrell JK, Berisa T, Liu JZ, Segurel L, Tung JY, Hinds DA. Detection and interpretation of shared genetic influences on 42 human traits. Nat Genet. 2016;48(7):709–17. pmid:27182965

24. Burgess S, Dudbridge F, Thompson SG. Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat Med. 2016;35(11):1880–906. pmid:26661904

25. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. Genome-wide genetic data on ~500,000 UK Biobank participants. bioRxiv. 2017 Jul 20.

26. Palmer TM, Sterne JA, Harbord RM, Lawlor DA, Sheehan NA, Meng S, et al. Instrumental variable estimation of causal risk ratios and causal odds ratios in Mendelian randomization analyses. Am J Epidemiol. 2011;173(12):1392–403. pmid:21555716

27. Didelez V, Meng S, Sheehan NA. Assumptions of IV methods for observational epidemiology. Stat Sci. 2010;25(1):22–40.

28. Greco MF, Minelli C, Sheehan NA, Thompson JR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med. 2015;34(21):2926–40. pmid:25950993

29. Vansteelandt S, Bowden J, Babanezhad M, Goetghebeur E. On instrumental variables estimation of causal odds ratios. Stat Sci. 2011;26:403–22.

30. Zhao Q, Wang J, Hemani G, Bowden J, Small DS. Statistical inference in two-sample summary data Mendelian randomization using a robust adjusted profile score. arXiv. 2018 Feb 7. arXiv: 1801.09652v2.

31. Thompson JR, Minelli C, Del Greco MF. Mendelian randomization using public data from genetic consortia. Int J Biostat. 2016;12(2). pmid:27092657

32. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37(7):658–65. pmid:24114802

33. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. 2017;36(11):1783–802. pmid:28114746

34. 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(2):512–25. pmid:26050253

35. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14. pmid:27061298

36. Staley JR, Blackshaw J, Kamat MA, Ellis S, Surendran P, Sun BB, et al. PhenoScanner: a database of human genotype-phenotype associations. Bioinformatics. 2016;32(20):3207–9. pmid:27318201

37. Gage SH, Jones HJ, Burgess S, Bowden J, Davey Smith G, Zammit S, et al. Assessing causality in associations between cannabis use and schizophrenia risk: a two-sample Mendelian randomization study. Psychol Med. 2017;47(5):971–80. pmid:27928975

38. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan NA, Thompson JR. Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic. Int J Epidemiol. 2016;45(6):1961–74. pmid:27616674

39. Cole SR, Platt RW, Schisterman EF, Chu H, Westreich D, Richardson D, et al. Illustrating bias due to conditioning on a collider. Int J Epidemiol. 2010;39(2):417–20. pmid:19926667

40. Draijer C, Hylkema MN, Boorsma CE, Klok PA, Robbe P, Timens W, et al. Sexual maturation protects against development of lung inflammation through estrogen. Am J Physiol Lung Cell Mol Physiol. 2016;310(2):L166–74. pmid:26608529

41. Jenkins MA, Dharmage SC, Flander LB, Douglass JA, Ugoni AM, Carlin JB, et al. Parity and decreased use of oral contraceptives as predictors of asthma in young women. Clin Exp Allergy. 2006;36(5):609–13. pmid:16650045

42. Vrieze A, Postma DS, Kerstjens HA. Perimenstrual asthma: a syndrome without known cause or cure. J Allergy Clin Immunol. 2003;112(2):271–82. pmid:12897732

43. Troisi RJ, Speizer FE, Willett WC, Trichopoulos D, Rosner B. Menopause, postmenopausal estrogen preparations, and the risk of adult-onset asthma. A prospective cohort study. Am J Respir Crit Care Med. 1995;152(4 Pt 1):1183–8. pmid:7551368

44. Mendle J, Turkheimer E, Emery RE. Detrimental psychological outcomes associated with early pubertal timing in adolescent girls. Dev Rev. 2007;27(2):151–71. pmid:20740062

45. Wyshak G, Frisch RE. Evidence for a secular trend in age of menarche. N Engl J Med. 1982;306(17):1033–5. pmid:7062994

46. Adams Hillard PJ. Menstruation in adolescents: what’s normal, what’s not. Ann N Y Acad Sci. 2008;1135:29–35. pmid:18574205

47. Fisher MM, Eugster EA. What is in our environment that effects puberty? Reprod Toxicol. 2014;44:7–14. pmid:23602892

48. Kelly T, Yang W, Chen CS, Reynolds K, He J. Global burden of obesity in 2005 and projections to 2030. Int J Obes (Lond). 2008;32(9):1431–7. pmid:18607383

49. Zama AM, Uzumcu M. Epigenetic effects of endocrine-disrupting chemicals on female reproduction: an ovarian perspective. Front Neuroendocrinol. 2010;31(4):420–39. pmid:20609371

50. Pierce MB, Leon DA. Age at menarche and adult BMI in the Aberdeen children of the 1950s cohort study. Am J Clin Nutr. 2005;82(4):733–9. pmid:16210700

51. Gong TT, Wu QJ, Vogtmann E, Lin B, Wang YL. Age at menarche and risk of ovarian cancer: a meta-analysis of epidemiological studies. Int J Cancer. 2013;132(12):2894–900. pmid:23175139

52. He C, Zhang C, Hunter DJ, Hankinson SE, Buck Louis GM, Hediger ML, et al. Age at menarche and risk of type 2 diabetes: results from 2 large prospective cohort studies. Am J Epidemiol. 2010;171(3):334–44. pmid:20026580

53. Hsieh CC, Trichopoulos D, Katsouyanni K, Yuasa S. Age at menarche, age at menopause, height and obesity as risk factors for breast cancer: associations and interactions in an international case-control study. Int J Cancer. 1990;46(5):796–800. pmid:2228308

54. Lakshman R, Forouhi NG, Sharp SJ, Luben R, Bingham SA, Khaw KT, et al. Early age at menarche associated with cardiovascular disease and mortality. J Clin Endocrinol Metab. 2009;94(12):4953–60. pmid:19880785

55. Macsali F, Real FG, Plana E, Sunyer J, Anto J, Dratva J, et al. Early age at menarche, lung function, and adult asthma. Am J Respir Crit Care Med. 2011;183(1):8–14. pmid:20732985

56. Gill D, Sheehan NA, Wielscher M, Shrine N, Amaral AFS, Thompson JR, et al. Age at menarche and lung function: a Mendelian randomization study. Eur J Epidemiol. 2017;32(8):701–10. pmid:28624884

57. Gill D, Brewer CF, Del Greco MF, Sivakumaran P, Bowden J, Sheehan NA, et al. Age at menarche and adult body mass index: a Mendelian randomization study. Int J Obes (Lond). 2018 Feb 26. pmid:29549348

58. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23(R1):R89–98. pmid:25064373

59. Harries ML, Walker JM, Williams DM, Hawkins S, Hughes IA. Changes in the male voice at puberty. Arch Dis Child. 1997;77(5):445–7. pmid:9487971

60. Ong KK, Bann D, Wills AK, Ward K, Adams JE, Hardy R, et al. Timing of voice breaking in males associated with growth and weight gain across the life course. J Clin Endocrinol Metab. 2012;97(8):2844–52. pmid:22654120