Age at menarche and endometrial cancer risk: a dose-response meta-analysis of prospective studies.

Evidence between age at menarche and endometrial cancer risk have been controversial. Therefore, we conducted a meta-analysis of prospective studies to analyze the aforementioned association. Relevant studies were identified by searching PubMed and EMBASE databases until the end of June 2015. A random-effects model was used to estimate summary relative risks (RRs) and 95% confidence intervals (CIs) for associations between menarcheal age and endometrial cancer risk. Our meta-analysis included eight prospective studies involving 4553 subjects with endometrial cancer. The summarized RRs of endometrial cancer for menarcheal age were 0.68 (95%CI = 0.58-0.81, I(2) = 41.9%, P = 0.099, n = 8) when comparing women with oldest category of menarcheal age with women with youngest category of menarcheal age. Notably, there was an 4% reduction in risk for per 2 years delay in menarcheal age (summarized RR = 0.96; 95%CI = 0.94-0.98, I(2) = 45.7%, P = 0.101, n = 6). Additionally, significant inverse associations were consistent within all stratified analyses. There was no evidence of publication bias or significant heterogeneity between subgroups detected by meta-regression analyses. Our findings support the hypothesis that late menarcheal age is inversely associated with endometrial cancer risk. Further larger prospective or pooled studies are warranted to fully adjust for potential confounders and distinguish whether the associations differ by histological subtypes of endometrial cancer.

Endometrial cancer is the second most common gynecologic malignancy worldwide, with approximately 280,000 new cases in 20121. Recently, the development of endometrial cancer may be partly attributed to both endogenous and exogenous of estrogen and progesterone2. For example, estrogens are the primary stimulants of endometrial proliferation3, which is a prerequisite for carcinogenesis4. Apparently, unchecked proliferation can lead to malignant transformation5, which follows that estrogens are a cause in the development of at least some endometrial cancer. In contrast, the surge of luteal progesterone secretion which begins just prior to ovulation each month in premenopausal women may play a role in arrest endometrial proliferation and promoting secretory differentiation of the endometrium, which can mitigate against the occurrence of endometrial cancer267. In 1986, Pettersson et al.8 first suggested that menstruation span was an important risk factor for endometrial cancer. On the basis of this finding, they addressed that the number of ovulatory cycles might be important in the etiology of this disease9. Therefore, reproductive characteristics, such as age at menarche, age at menopause, and nulliparity, that are closely associated with endogenous hormone changes may affect the risk of endometrial cancer2.

Menarche was considered to be the milestone of ovulation initiation as well as the initiation of hormone changes in the childhood and adolescent period10. Our previous study found a statistically significant inverse association between later menarcheal age and ovarian cancer risk11. Furthermore, some studies suggested that increasing rate of overweight/obesity in childhood and adolescent period was inverse associated with the menarcheal age121314. In addition, early menarcheal age increased the risk of all-cause or cardiovascular disease mortality in several studies during the past decade15161718. Although the increased risk of endometrial cancer associated with early menarcheal age has been attributed to a longer lifetime exposure to endogenous estrogen and progesterone deficiency associated with anovulatory cycles192021, the evidence from observational studies are conflicting, which might be partly attributed to different study designs and limited sample sizes of individual studies. Notably, minority of these studies provided the detail information of the risk estimation of per year advance/delay in menarcheal age by dose-response analysis. Therefore, we carried out this meta-analysis of all prospective studies published up to May 2014 to systematically and quantitatively evaluate this issue.

Results

Literature search

We identified 4607 potentially relevant articles from our search of the MEDLINE (PubMed) and EMBASE databases. Of these, 4593 articles were excluded after the first screening based on the abstract or title, leaving 14 articles for full-text review. On this review, two articles2223 were excluded because of duplicate reports from the same study population and four articles24252627 were excluded because they did not report sufficient data of risk estimates. Overall, 8 prospective studies were included in this study (Fig. 1).

Figure 1

Selection of studies for inclusion in meta-analysis.

Study characteristics

Table 1 summarizes the characteristics of the 8 included studies. These included articles which represent 4,553 cases and 949,945 non-cases, were published between 1988 and 2013, and consist of 7 cohort studies9212829303132 and one nested case-control study33. Of the 8 prospective studies, five were conducted in the United States2829303133, one each in China9, Norway32, and the European Prospective Investigation into Cancer and Nutrition (EPIC)21. Sample sizes of these included studies ranged from 15,52833 to 301,60121, and the number of endometrial cases varied from 3933 to 145028. The median duration of follow-up was 11 years.

Table 1

Characteristics of included studies of menarcheal age and endometrial cancer risk.

First author, publication year [reference], Country, Study design	Cases/subject (age), duration of follow up	Source of cases	Menarcheal age categories (exposure assessment)	RR/HR (95%CI)	Matched/Adjusted factors
Yang et al.28, 2013, USA, CS	1450/114,409 (50–69y), 9.4y	Cancer registry	<13	1.00 (ref)	Age, OC use, menopausal hormone therapy use, parity, BMI, age at menopause, race, and smoking status
			13–14	0.91 (0.82–1.02)
			≥15 (Self-administered questionnaire)	0.80 (0.65–0.97)
Dossus et al.21, 2010, Europe, CS	1017/301,601 (35–70y), 8.7y	Cancer registry	<12	1.00 (ref)	Age, center, BMI, physical activity, alcohol, diabetes, smoking status and education
			12	0.97 (0.80–1.19)
			13	0.78 (0.63–0.95)
			14	0.80 (0.65–0.98)
			≥15 (Self-administered questionnaire)	0.72 (0.58–0.90)
Karageorgi et al.29, 2010, USA, CS	778/121,700 (30–55y), 28y	Medical records	≤11	1.00 (ref)	Age, parity, age at first birth, age at last birth, OC duration, postmenopausal hormone use duration, BMI, smoking, diabetes, family history of endometrial cancer and age at menopause
			12	0.99 (0.81–1.20)
			13	0.80 (0.66–0.98)
			14	0.90 (0.70–1.16)
			≥15 (Self-administered questionnaire)	0.76 (0.55–1.04)
Setiawan et al.31, 2007, USA, CS	321/46,933 (45–75y), 7.3y	Cancer registry	≤12	1.00 (ref)	Race/ethnicity, BMI, age at natural menopause, parity, hormone therapy use, OC use, smoking, diabetes, hypertension, family history of endometrial cancer
			13–14	0.86 (0.68–1.10)
			≥15 (Self-administered questionnaire)	0.85 (0.60–1.22)
Wernli et al.9, 2006, China, CS	206/267,400 (≥30y), 10y	Medical records	≤13	1.00 (ref)	Age and parity number
			14	0.48 (0.28–0.82)
			15	0.73 (0.47–1.16)
			16	0.77 (0.49–1.20)
			≥17 (Self-administered questionnaire)	0.40 (0.25–0.63)
Hisada et al.33, 2001†, USA, NC-CS	39/194/15,528 (N/A), 27y	Cancer registry	≤11	1.00 (ref)	Age and race
			12	0.89 (0.18–4.30)
			13	1.44 (0.34–6.09)
			≥14 (Trained interviewer)	0.67 (0.10–4.33)
Olson et al.30, 1999, USA, CS	322/24,848 (55–69y), 10y	Cancer registry	≤11	1.00 (ref)	Age
			11–12	0.7 (0.5–1.2)
			12–13	0.6 (0.4–1.0)
			>14 (Self-administered questionnaire)	0.4 (0.2–0.8)
Kvale et al.32, 1988, Norway, CS	420/62,079 (27–69y), 20y	Medical records	≤12	1.00 (ref)	Age at start of follow-up, urban/rural place of residence, and parity
			13	N/A
			14	N/A
			15	N/A
			16	N/A
			≥17 (Self-administered questionnaire)	0.58 (0.39–0.88)

RR: relative risk; CI: confidence interval; CS: cohort study; NC-CS: nested case-control study; N/A: not available; BMI: body mass index; OC: oral contraceptive.

†Odds ratio and 95% CI calculated from published data using EpiCalc 2000.

Late versus early age at menarche

The summary RR of endometrial cancer for the oldest versus the youngest categories of menarcheal age was 0.68 (95% CI: 0.58–0.81) with moderate heterogeneity (I2 = 41.9%, P = 0.099) (Table 2, Fig. 2). There was no indication of publication bias with Egger’s test (P-bias = 0.189) and no asymmetry in the funnel plots when inspected visually.

Table 2

Summary risk estimates of the association between menarcheal age and endometrial cancer risk, oldest versus youngest category.

	No. of	Summary RR	I2	Ph*	Ph**
	studies	(95% CI)	Value (%)
Overall	8	0.68 (0.58–0.81)	41.9	0.099	—
Subgroup analyses
Duration of follow-up					0.085
<10y	3	0.78 (0.68–0.89)	0	0.674
≥10y	5	0.56 (0.42–0.74)	37.2	0.173
Number of cases					0.281
<400	4	0.74 (0.65–0.84)	62.4	0.046
≥400	4	0.55 (0.33–0.90)	0	0.563
Exposure Assessment					0.992
Trained interviewer	1	1.01 (0.82–1.24)	N/A	N/A
Self-administered questionnaire	7	0.68 (0.57–0.81)	50.2	0.061
Study population					0.232
Non-American	3	0.58 (0.41–0.81)	62.2	0.071
American	5	0.77 (0.67–0.90)	0	0.421
Adjustment for confounders or important risk factors
Body mass index					0.057
Yes	4	0.77 (0.68–0.88)	0	0.849
No	4	0.48 (0.36–0.63)	0	0.612
Parity					0.710
Yes	5	0.69 (0.55–0.86)	57.0	0.054
No	3	0.64 (0.45–0.91)	20.3	0.285
Oral contraceptive use					0.085
Yes	3	0.80 (0.69–0.93)	0	0.900
No	5	0.56 (0.42–0.74)	41.8	0.143
Exogenous hormones use					0.085
Yes	3	0.80 (0.69–0.93)	0	0.900
No	5	0.56 (0.42–0.74)	41.8	0.143
Menopause status					0.085
Yes	3	0.80 (0.69–0.93)	0	0.900
No	5	0.56 (0.42–0.74)	41.8	0.143
Smoking status					0.057
Yes	4	0.77 (0.68–0.88)	0	0.849
No	4	0.48 (0.36–0.63)	0	0.612

RR: relative risk; CI: confidence interval; N/A: not available.

*P value for heterogeneity within each subgroup.

**P value for heterogeneity between subgroups with meta-regression analysis.

Figure 2

Forest plot (random effects model) of menarcheal age and endometrial cancer risk in prospective studies.

Squares indicate study-specific relative risks (size of the square reflects the study-specific statistical weight); horizontal lines indicate 95% CIs; diamond indicates the summary relative risk estimate with its 95% CI. CI: confidence interval; RR: relative risk.

Dose-response analysis

Six studies were included in the dose-response analysis92128293033. There was an 4% reduction in risk for per 2 years delay in menarcheal age (summarized RR = 0.96; 95% CI = 0.94–0.98) with moderate but non-significant heterogeneity (I2 = 45.7%, P = 0.101) (Table 3, Fig. 3). There was no indication of a publication bias by using Egger’s test (P-bias = 0.152) and no asymmetry was observed in funnel plots when inspected visually. There was no evidence of a nonlinear association between the menarcheal age and endometrial cancer risk (P for nonlinearity = 0.16).

Table 3

Summary risk estimates of the association between menarcheal age and endometrial cancer risk, dose-response analysis of per 2 year delay.

	No. of	Summary RR	I2	Ph*	Ph**
	studies	(95% CI)	Value (%)
Overall	6	0.96 (0.94–0.98)	45.7	0.101	—
Subgroup analyses
Duration of follow-up					0.303
<10y	2	0.97 (0.95–0.99)	44.7	0.179
≥10y	4	0.95 (0.91–0.98)	33.7	0.210
Number of cases					0.067
<400	3	0.92 (0.89–0.96)	0	0.868
≥400	3	0.97 (0.96–0.99)	0	0.388
Exposure Assessment					0.680
Trained interviewer	1	1.01 (0.82–1.24)	N/A	N/A
Self-administered questionnaire	6	0.96 (0.94–0.98)	55.8	0.066
Study population					0.259
Non-American	2	0.95 (0.91–0.98)	56.3	0.130
American	4	0.98 (0.96–0.99)	7.5	0.356
Adjustment for confounders or important risk factors
Body mass index					0.067
Yes	3	0.97 (0.96–0.99)	0	0.388
No	3	0.92 (0.89–0.96)	0	0.868
Parity					0.652
Yes	3	0.96 (0.94–0.99)	66.5	0.051
No	3	0.96 (0.93–0.98)	0	0.469
Oral contraceptive use					0.080
Yes	2	0.98 (0.97–0.99)	0	0.560
No	4	0.95 (0.92–0.97)	12.2	0.332
Exogenous hormones use					0.080
Yes	2	0.98 (0.97–0.99)	0	0.560
No	4	0.95 (0.92–0.97)	12.2	0.332
Menopause status					0.080
Yes	2	0.98 (0.97–0.99)	0	0.560
No	4	0.95 (0.92–0.97)	12.2	0.332
Smoking status					0.067
Yes	3	0.97 (0.96–0.99)	0	0.388
No	3	0.92 (0.89–0.96)	0	0.868

RR: relative risk; CI: confidence interval; N/A: not available.

*P value for heterogeneity within each subgroup.

**P value for heterogeneity between subgroups with meta-regression analysis.

Figure 3

Dose-response analysis (random effects model) between per 2 year delay in menarcheal age and risk of endometrial cancer.

Squares indicate study-specific relative risks (size of the square reflects the study-specific statistical weight); horizontal lines indicate 95% CIs; diamond indicates the summary relative risk estimate with its 95% CI. CI: confidence interval; RR: relative risk.

Subgroup analyses

In subgroup analyses, possible differences between risk estimates by various study characteristics were examined. The finding of decreased endometrial cancer risk with later menarcheal age was consistently observed in almost all of the stratified analyses, although not all strata showed statistical significance. Furthermore, there was no evidence of significant heterogeneity between subgroups detected by meta-regression analyses (Table 2 and Table 3). When stratified by whether considering or adjustment for potential confounders, we did not find a significant difference between estimates adjusted and those not adjusted for specific factors (Table 2). Likewise, similar results were also observed in the dose-response analysis (Table 3).

Sensitivity analyses

In a sensitivity analysis, we sequentially removed one study at a time and re-analyzed the data. The 8 study-specific RRs of the oldest versus youngest menarcheal age ranged from a low of 0.65 (95% CI: 0.53–0.79, I2 = 40.8%, P = 0.119) after omission of the study by Yang et al.28 to a high of 0.74 (95% CI: 0.66–0.83, I2 = 0%, P = 0.461) after omission of the study by Wernli et al.9. Similar analyses were also carried out in the dose-response analysis, with study specific RRs ranged from a low of 0.95 (95% CI: 0.93–0.97, I2 = 15.7%, P = 0.314) after omission of the study by Yang et al.28 to a high of 0.97 (95% CI: 0.96–0.99, I2 = 10.7%, P = 0.345) after omission of the study by Wernli et al.9.

Discussion

To the best of our knowledge, this is the first comprehensive and quantitative assessment to explore the relationship between menarcheal age and risk of endometrial cancer. Overall, menarcheal age was found to be associated with an estimated 32% reduction in endometrial cancer risk when the oldest menarcheal age category compared to the youngest. Additionally, the results of dose-response analyses suggested an 4% reduction in endometrial cancer risk for per 2 years delay in menarcheal age (summarized RR = 0.96, 95% CI: 0.94–0.98). These significant results were consistent within subgroup analyses and across sensitivity analyses.

The association between hormonal risk factors and endometrial cancer risk could be partly attributed to the “unopposed estrogen” hypothesis. This hypothesis proposes that exposure to estrogen, in the absence of sufficient progestin, leads to increase mitotic activity, DNA replication, and somatic mutations of endometrial cells that may result in malignant transformations23453435. In addition, as the beginning of the menstruation span, later menarcheal age might decrease the risk of endometrial cancer by decreasing a woman’s lifetime number of ovulations which was characterized by exposure to endogenous estrogen and progesterone deficiency. Several epidemiologic studies8920212335 have reported that longer year of menstruation span significantly increased the risk of endometrial cancer, which partly supported the findings in this meta-analysis, as well as the hypothesis which suggested that endometrial cancer risk might related to the lifetime number of ovulations. A recent study36, the Epidemiology of Endometrial Cancer Consortium (E2C2), has pooled of 15 observational studies (10 prospective and 5 case-control studies), provided evidence that older menarcheal age was significant inverse associated with endometrioid endometrial carcinoma (≥15 versus. <11 menarcheal age, OR: 0.63, 95% CI: 0.54–0.72, P trend = 0.04), which partly supported the findings of this meta-analysis. However, none of the included studies of this meta-analysis demonstrated the information of different histological type which limited us to confirm their finding. On the other hand, Type I endometrial cancers are mostly endometrioid adenocarcinomas, which seem to develop from abnormal glandular proliferations (i.e., endometrial hyperplasia) driven by hormonal mechanisms. In contrast, Type II endometrial cancers often display serous or clear cell histology and arise from atrophic endometrium in a less hormonally dependent manner. In addition, subtypes of these cancers are characterized by distinctive molecular alterations, and endometrioid carcinomas are more clearly linked to elevated levels of sex-steroid hormones and expression of hormone receptors283738. However, only Yang et al.28 divided the main analyses by the subtype of endometrial cancer. Given this, further studies are warranted to focus on this issue. In our previous meta-analysis of menarcheal age and ovarian cancer risk11, we found that the categories used for the age at menarche varied from studies. Although the results of meta-regression of our study found no evidence for heterogeneity of different study populations, each population was unique and the category of menarcheal age might different from other studies. For example, Wernli et al.9 reported that the youngest and oldest category of menarcheal age in a cohort study with 267,400 women which was carried out in China was less than 13 years and over than 17 years, respectively. By comparison, Olson et al. reported that the youngest and oldest category of menarcheal age in a cohort study with 24,848 women which was carried out in Iowa was less than 11 years and over than 14 years, respectively. Since this, the interpretation of the results when using the oldest category of menarcheal age compared to the youngest should be cautiously. However, the results of dose-response analysis which could reduce the impact from the different category of menarcheal age of these included studies were in accordance with the results of oldest versus youngest category analysis which not only demonstrated the findings of this meta-analysis were robust but suggested that this aforementioned association was independent of the category methods. On the other hand, several genome-wide association studies (GWAS) have identified that genetic factors may associated with menarcheal age39404142 which also might attributed to the different category of menarcheal age. Therefore, given six prospective studies were included in the dose-response analysis of this study, more epidemiologic studies should provide the risk estimates which modeling a one or two-year increment in menarcheal age continuously, as well as focus on whether genetic polymorphisms may modify the association between menarcheal age and endometrial cancer risk.

A strength of this study is that it included a relatively larger sample sizes than any other individual studies with a total of 4553 endometrial cancer cases and approximately 949,945 non-cases. Thus, we had statistical power to detect moderate and weak associations and carry out subgroup analyses by potential confounding factors and study characteristics. Another important strength of this meta-analysis was that we included all the prospective design studies and the results are unlikely to be explained by the bias of traditional retrospective studies.

This study, however, also has limitations. First, since all included prospective studies were based on the observational design, it is possible that these aforementioned observed significant inverse associations could be due to the unmeasured or residual confounding. Many factors known to affect endometrial cancer risk, especially BMI1243 is either associated with menarcheal age or as an important risk factor for endometrial cancer. However, since we could not get the access of raw data of each included studies and none of them have clarified the detail information of potential confounders in their published articles, thus we carried out the stratified analyses by these potential confounders which were adjusted in the multivariable models of their primary analyses. On the other hand, cohort studies more likely to present the age at menarche or other important risk factors of endometrial cancer by the distribution of exposure instead of by case status, which might partly results in limited studies included in this meta-analysis. Although numerous subgroup analyses were carries out and the results generally showed inverse associations, the limitation by the numbers of included prospective studies still should be considered. Second, because of the majority of the prospective study included the participants later than the focused exposure happened, thus females included in this study generally reported their exposure information long after menarche which may result in potential information bias. However, given that menarche is the milestone of puberty initiation and Bean et al.44 previously reported that age at menarche has reasonably good recall accuracy. Therefore, this kind of information bias could be limited in our study. Third, to our knowledge, the category of menarcheal age differed between studies and may have contributed to the heterogeneity in results. However, few of the included studies reported how they categorized the menarcheal age, and thus, we hardly considered this point in the subgroup analysis and ruled out the heterogeneity thoroughly. Furthermore, publication bias can also be a problem in meta-analyses of published studies but we found no statistical evidence of publication bias in this meta-analysis. Finally, as the starting point of a woman’s menstrual history, menarcheal age only focused on the early period of menstrual history, future studies should provide more evidence of the reproductive factors of other periods.

In summary, this meta-analysis suggested that later menarcheal age was inversely associated with the risk of endometrial cancer. Further studies, especially large consortium or pooled studies would be of interest to fully adjust for potential confounders and distinguish whether other risk factors which closely associated with both menstruation span and hormone change (e.g., age at menopause and multiparity) might have the similar inverse association with endometrial cancer risk and whether these associations might differ by the histological subtypes of endometrial cancer.

Methods

We performed a comprehensive literature search to the end of June 2015 using the MEDLINE (PubMed) and EMBASE databases limited to the studies of humans by using the following search key words and Medical Subject Headings terms: (menarche OR reproductive factors OR reproduction) AND (endometrial OR endometrium OR corpus uteri OR uterine corpus) AND (cancer OR neoplasm OR carcinoma OR tumor). Additionally, the reference lists of retrieved articles were also scrutinized for additional studies by manual search. A similar search strategy was utilized for our previous meta-analysis of menarcheal age and ovarian cancer11. This study was planned, conducted, and reported in adherence to standards of quality for reporting meta-analyses45.

Study selection criteria

Published studies were included by the following selection criteria if they 1) used a prospective study design, including cohort, case-cohort, and nested case-control studies; 2) evaluated the association between menarcheal age and incident endometrial cancer risk; and 3) presented relative risk (RR) or hazard ratio (HR) estimates with 95% confidence intervals (CI), standard errors (SE) or data necessary to calculate these.

Published manuscripts were excluded by the following exclusion criteria if they 1) had a retrospective design; 2) the estimates were presented without SE or other information that allowed calculation of SE; 3) reported exclusively on endometrial cancer mortality. When multiple publications from the same study were available, we used the publication with the largest number of cases and most applicable information.

All above study selection and exclusion procedures were carried out by two independent investigators (T-TG and X-XM).

Data abstraction and quality assessment

For each eligible study, two investigators (T-TG and X-XM) independently performed the eligibility evaluation and data abstraction. Discrepancies were settled by consensus or by involving a third reviewer (Y-LW) for adjudication. Data abstracted from each study were: author list, year of publication, country where the study was performed, study design, study sample size (number of cases and cohort size), range of follow-up for cohort studies, exposure and outcome assessment and menarcheal age categories, study-specific adjusted RRs or HRs with their 95% CIs for the oldest versus youngest category of menarcheal age (we extracted the RRs that reflected the greatest degree of control for potential confounders for use in the main analyses), and factors matched between cases and controls in nested case-control study and potential confounders adjusted in the data analysis.

Since all included studies are prospective study and similar to our previous studies1146, we did not assign quality scores47 to assess the methodological quality of all the included studies which lacks demonstrated validity, and sometimes results may not be associated with quality48, but investigated whether specific study characteristics, such as study design and adjustment for confounders, which are indicators of study quality, influenced the results in subgroup analyses.

Statistical analysis

The study-specific adjusted RRs were used as the common measure of association across studies. For the meta-analysis, we assumed that estimates of risk, rate or hazard ratios from prospective studies were all valid estimates of the RR and we therefore report all results as the RR for simplicity. For study of Yang et al.28 that reported results separately of type I and type II endometrial cancer, but not combined, we pooled the results using a fixed-effect model to obtain an overall combined estimate before combining with the rest of the studies.

The possible heterogeneity in results across studies was examined by using the I2 statistics49. The I2 statistic represents the proportion of total variation contributed by between-study variation49. We used a random-effects model to calculate summary RR of endometrial cancer with 95% CIs, considering within- and between-study variation50. Heterogeneity between subgroups was evaluated by meta-regression46.

For the dose-response analysis, we used the generalized least-squares trend estimation method developed by Greenland et al.51 and Orsini et al.52 to compute study-specific slopes (linear trends) and 95% CIs from the natural logs of the RRs and CIs across categories of the menarcheal age. The aforementioned method requires the following information: 1) the distribution of cases and person-years or non-cases and the RRs with the variance estimates for at least three quantitative exposure categories; 2) the median or mean level of these exposures in each category (if reported by ranges, mean level were calculated by averaging the lower and upper bound; if the lowest category was open ended, the lowest boundary was considered to be zero; if the highest category was open ended, the open-ended interval length was assumed to be the same as the adjacent interval). Since these criteria, six studies92128293033 were included in the dose-response analysis of menarcheal age and endometrial cancer risk. Furthermore, a potential nonlinear dose-response relationship between the menarcheal age and endometrial cancer risk was modeled by using restricted cubic splines with 3 knots at fixed percentiles (10%, 50%, and 90%) of the distribution of exposure. We calculated an overall P value by testing that these two regression coefficients were simultaneously equal to zero. We calculated a P value for nonlinearity by testing that the coefficient of the second spline was equal to zero. The details of this method has been published elsewhere5354. The dose-response results are presented for a two years delay in menarcheal age.

To investigate possible sources of heterogeneity of main results, we carried out stratified analyses by the following study features: duration of follow-up (<10 versus ≥10 years), number of cases (<400 versus ≥400), method of exposure assessment (self-administered questionnaire versus trained interviewers), study population (American versus Non-American), and factors matched between cases and controls in nested case-control study and potential confounders adjusted in the data analyses (body mass index, parity, oral contraceptive use, exogenous hormones use, menopause status, and smoking status). Finally, sensitivity analysis was executed by deleting each study in turn to reflect the influence of individual data set on the overall estimate.

Small study bias, such as publication bias, was evaluated by Egger’s linear regression55 and funnel plots. A P-value less than 0.05 for Egger’s test was considered representative of significant statistical publication bias. Statistical analyses were performed with Stata (version 11.2; StataCorp, College Station, TX). P-values were two sided with a significance level of 0.05.

Additional Information

How to cite this article: Gong, T.-T. et al. Age at menarche and endometrial cancer risk: a dose-response meta-analysis of prospective studies. Sci. Rep. 5, 14051; doi: 10.1038/srep14051 (2015).