Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Genetic variations, reproductive aging, and breast cancer risk in African American and European American women: The Women's Circle of Health Study

  • Marie V. Coignet,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, NY, United States of America

    ⨯
  • Gary Robert Zirpoli,

    Roles Data curation, Resources, Software, Validation, Visualization, Writing – review & editing

    Affiliations Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, NY, United States of America, Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States of America

    ⨯
  • Michelle R. Roberts,

    Roles Data curation, Resources, Software, Validation, Visualization, Writing – review & editing

    Affiliations Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, NY, United States of America, Channing Division of Network Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, United States of America

    ⨯
  • Thaer Khoury,

    Roles Data curation, Resources, Validation, Visualization, Writing – review & editing

    Affiliation Department of Pathology, Roswell Park Cancer Institute, Buffalo, NY, United States of America

    ⨯
  • Elisa V. Bandera,

    Roles Data curation, Funding acquisition, Methodology, Resources, Visualization, Writing – review & editing

    Affiliation Cancer Prevention and Control Program, Rutgers Cancer Institute of New Jersey, New Brunswick, NJ, United States of America

    ⨯
  • Qianqian Zhu,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – review & editing

    Affiliation Department of Biostatistics, Roswell Park Cancer Institute, Buffalo, NY, United States of America

    ⨯
  • Song Yao

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    song.yao@roswellpark.org

    Affiliation Department of Cancer Prevention and Control, Roswell Park Cancer Institute, Buffalo, NY, United States of America

    ⨯

Abstract

Reproductive aging phenotypes, including age at menarche (AM) and age at natural menopause (ANM), are well-established risk factors for breast cancer. In recent years, many genetic variants have been identified in association with AM and ANM in genome-wide association studies among European populations. Using data from the Women’s Circle of Health Study (WCHS) of 1,307 European-American (EA) and 1,365 African-American (AA) breast cancer cases and controls, we aimed to replicate 53 earlier GWAS variants for AM and ANM in AA and EA groups and to perform analyses on total and net reproductive lifespan (TRLS; NRLS). Breast cancer risk was also examined in relation to a polygenic risk score (PRS) for each of the reproductive aging phenotypes. We replicated a number of variants in EA women, including rs7759938 in LIN28B for AM and rs16991615 in MCM8 for ANM; whereas in the AA group, only one SNP (rs2947411 in TMEM18) for AM was directionally consistent and nominally significant. In analysis of TRLS and NRLS, several SNPs were significant, including rs466639 in RXRG that was associated with both phenotypes in both AA and EA groups. None of the PRS was associated with breast cancer risk. Given the paucity of data available among AA populations, our study contributes to the literature of genetics of reproductive aging in AA women and highlights the importance of cross population replication of GWAS variants.

Citation: Coignet MV, Zirpoli GR, Roberts MR, Khoury T, Bandera EV, Zhu Q, et al. (2017) Genetic variations, reproductive aging, and breast cancer risk in African American and European American women: The Women's Circle of Health Study. PLoS ONE 12(10): e0187205. https://doi.org/10.1371/journal.pone.0187205

Editor: Joseph Devaney, GeneDx, UNITED STATES

Received: March 5, 2017; Accepted: October 16, 2017; Published: October 26, 2017

Copyright: © 2017 Coignet et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Per the restrictions from our Institutional Review Board (IRB) and the consent agreed by the study participants, we cannot grant open access to the genotypic or phenotypic data without restriction. However, on a case-by-case basis, de-identified data from this work may be requested from the senior author Dr. Song Yao, should human subject protection and patient privacy be secured. Request can be made to Dr. Song Yao at song.yao@roswellpark.org.

Funding: This work was supported in part by grants from the US National Institutes of Health (R01 CA100598, P01 CA151135, and R03 CA192205 and P30 CA016056) and US Army Medical Research and Material Command (DAMD-17-01-1-0334). The genotyping work was performed by the Roswell Park Cancer Institute Genomics Shared Resource, which is a National Cancer Institute Cancer Center Supported Grant (CCSG) shared resource (P30 CA016056). The New Jersey State Cancer Registry is supported by the National Program of Cancer Registries of the Centers for Disease Control and Prevention under cooperative agreement 1US58DP003931-01 awarded to the New Jersey Department of Health. The collection of New Jersey cancer incidence data is also supported by the Surveillance, Epidemiology, and End Results program of the National Cancer Institute under contract N01-PC-2010-0027 and the State of New Jersey. The funding agents played no role in the design, collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Menarche and menopause are two fundamental physiological events in a woman’s life, respectively marking the beginning and the end of reproductive age. At those turning points, estrogen milieu inside a woman’s body undergoes drastic changes, and it is the lifetime exposure to sex hormones that is long speculated to have a profound impact on women’s health, extended from reproduction to hormone-related morbidities including breast cancer. Early age at menarche (AM) and late age at natural menopause (ANM) are well-established risk factors for breast cancer. The risk was estimated to increase by a factor of 1.050 for each year earlier at menarche and by 1.029 for each year older at menopause[1].

As two complex numeric human traits, the timing of AM and ANM varies greatly among individuals. They are under the influence of an intricate set of social, environmental and genetic factors. It was estimated that genetic inheritance explains over half of the variation (heritability for AM at 53–74% and for ANM at 44–65%)[2–9], which provides a strong rationale for genetic association studies to identify genes and variants determining their timing. To date, at least 17 genome-wide association studies (GWAS) have been conducted on AM and ANM[10–27], which have discovered a large number of common (minor allele frequency >5%) variants reliably associated with each of them. Results from these studies suggest that the genetic architecture of AM and ANM is polygenic, involving hundreds of variants, each with small effect, in multiple genes and biological pathways. It is worth noting that AM and ANM are typically studied separated and few studies have examined both in the same population; thus, it is largely unclear to what extent those two reproductive aging phenotypes share genetic determinants.

Racial/ethnic minority populations, particularly African Americans (AA), are under-represented in the GWAS of reproductive aging phenotypes in previous studies. Among the published GWAS, only one on AM[22] and one on ANM[23] were conducted among AA women. There are at least two reasons that more research efforts are warranted in genomic studies in minority populations. First, given the racial/ethnic diversity in genetic background and possibly in pathogenesis of cancer, it is not only prudent but imperative to replicate variants initially identified from populations of European ancestry (EA) in racial minorities, before the findings can be generalized across populations[28]. Second, AA women have distinct reproductive profiles, compared to EA women; for example, it has been shown that compared to EAs, AA women experience an earlier AM[29–31] as well as an earlier ANM[32, 33]. It is thus possible that the genetic architecture underlying reproductive aging phenotypes is distinct in AA populations. Further, reproductive factors have been postulated as important contributors to breast cancer health disparities between AA and EA women, including an earlier age at diagnosis and more aggressive pathological presentation in the former. Two recent studies show evidence that genetic variants associated with AM are also associated with breast cancer risk[27, 34]. It yet remains to be explored whether those variants contribute to the health disparities.

In the Women’s Circle of Health Study (WCHS), which enrolled a similar number of AA and EA women, with and without breast cancer, we aimed to replicate earlier GWAS findings in both AA and EA women, as well as to examine them with breast cancer risk in the two groups.

Methods

Study population

The WCHS was designed to study genetic and non-genetic factors in association with early and aggressive breast cancer risk in EA and AA women. A more detailed description of the study population has been published elsewhere[35–38]. In brief, case ascertainment in early years of the study was conducted through hospital-based referrals at targeted hospitals with large AA patient populations in 4 boroughs of New York City, and later shifted to population-based in 7 counties of New Jersey through the NJ State Cancer Registry. Case eligibility was determined by being a self-identified AA or EA woman, English speaking, between 20–75 years of age at diagnosis of primary, histologically confirmed breast cancer, with no history of cancer except for non-melanoma skin cancer. Since AA women are more likely to be diagnosed with breast cancer at an earlier age, younger women were oversampled in the enrollment. As a result, the population was younger relative to many other studies of breast cancer. Controls with no prior cancer history and matched to the cases by self-reported race, 5 year age categories and geographic residence were initially identified through random digit dialing (RDD), and later through community recruitment, especially to boost enrollment of AA controls. In a methodological analysis in the WCHS, a combination of community controls and RDD controls were found to be better representative of the general population than RDD controls alone[39]. Written informed consent was obtained from all participants and an in-person interview was administered to collect a wide range of epidemiological data, anthropometric measures and biospecimens. Blood samples were initially collected and later switched to saliva as a source of genomic DNA. The informed consent included permission to request pathology data and tumor tissue blocks and slides from attending hospitals. This study was approved by the Institutional Review Boards at Roswell Park Cancer Institute (RPCI), Rutgers Cancer Institute of New Jersey (CINJ), Icahn School of Medicine at Mount Sinai, and the participating hospitals in NYC.

SNP selection and genotyping

Our goal was to replicate single nucleotide polymorphisms (SNPs) that were previously identified in GWAS of AM and ANM in EA populations. Therefore, we selected 53 SNPs from 5 GWAS published by June 2011[13, 16, 24–26] when a custom genotyping chip was designed. These include 41 AM SNPs and 12 ANM SNP. Detailed information of the selected SNPs is summarized in S1 Table. Genotyping was attempted for 2,762 participants who had been enrolled into the WCHS and had genomic DNA samples available at that time, using the Illumina GoldenGate assay (Illumina Inc., San Diego, CA) performed by the Genomics Shared Resource at Roswell Park Cancer Institute. The average per-sample successful genotyping rate was 95.9%; 38 samples had calls less than 90% and were excluded. Quality control samples consisted of 5% blind duplicates and two sets of in-house trio samples. Concordance was >99.9% between duplicate pairs and there were no Mendelian errors in trio control samples. No SNPs selected for this analysis had a successful call rate <90%, or violated Hardy Weinberg Equilibrium when checked separately in AA and EA controls. SNPs with significant associations identified in the analyses were manually checked post hoc for the clustering plots, to ensure call robustness.

Genetic ancestry estimation

A panel of ancestry informative markers (AIMs) developed by the Black Women’s Health Study was included as part of the customized genotyping array and the data were used to ascertain genetic ancestry and to control for population admixture[40]. Estimates of the proportion of European and African ancestry were obtained using the STRUCTURE program[41]. Individuals (n = 11 in AAs and n = 1 in EAs) with a greater than 85% estimated ancestry that was discordant with the self-reported ancestry were excluded from analyses. The estimated percent European genetic ancestry was included as a numeric confounder in multivariate models.

Statistical analysis

The final dataset included 2,672 women, consisting of 658 EA breast cancer cases, 649 EA controls, 621 AA cases, and 744 AA controls. All analyses were performed separately for AA and EA women. T-tests or chi-square tests were used to compare descriptive characteristics between breast cancer cases and controls in univariate analyses. The two primary phenotypes, AM and ANM, were obtained via self-reported menstrual and pregnancy history. Women were defined as postmenopausal if they reported that they had ceased menstruation naturally for at least 12 consecutive months and those who had hysterectomy or oophorectomy were excluded from analysis of ANM, TRLS, and NRLS (n = 356). In addition to AM and ANM, we also computed another two phenotypes, total reproductive lifespan (TRLS) and and net reproductive lifespan (NRLS), using the numeric variables of AM and ANM. We considered TRLS to be the difference between ANM and AM, and NRLS to be equal to TRLS with time spent on all pregnancies and breastfeeding subtracted. Except for AM, the other three phenotypes were analyzed only among postmenopausal women. For AM analyses, 1365 AA women and 1307 EA women contributed in these analyses. For ANM/TRLS/NRLS analyses, 381 AA women and 465 EA women were included in the analyses.

Each single SNP was related to one of the four reproductive aging phenotypes, AM, ANM, TRLS and NRLS, using generalized linear modeling (GLM) with adjustment chosen a priori for birth year, estimated proportion of European ancestry, smoking (not for AM), and number of pregnancies (not for AM or NRLS) ensuring a parsimonious model. Additive genetic models were tested, with the common allele in the EA group designated for both AA and EA groups as the reference allele. Genotype analyses were first conducted in the case and control groups separately, which produced estimates of similar magnitude in the same direction. Thus, all analyses were conducted with case-control combined and the case-control status was included as a covariate in the GLMs. For the purpose of our study, replication of previous GWAS hits was defined by nominal significance of p≤0.05 and concordance in direction of effect for the same phenotype.

To test whether SNPs associated with reproductive aging had an impact on breast cancer risk, a polygenic risk score (PRS) was calculated for each of the four reproductive phenotypes, combining the effects of all SNPs associated with the phenotype at a nominal p-value <0.10. PRS for each phenotype and for each individual was calculated as a sum of the number of risk alleles, weighted by the regression coefficients from GLM. To align SNPs with opposite directions in their associations with a phenotype, the reference allele of those with a negative regression coefficient was flipped, so the regression coefficient became positive and could be tallied with other SNPs. As a result, a higher PRS indicated an older AM or ANM, or a longer TRLS and NRLS. Each PRS was then categorized into quartiles based on the distribution among controls, and related to breast cancer case-control status with a priori adjustment for age, estimated proportion of European ancestry, and family history of breast cancer, ensuring a parsimonious model. All analyses were performed using SAS 9.4 (SAS Institute, Cary, NC).

Results

Descriptive characteristics of the study population

Table 1 summarizes selected descriptive characteristics of the study population by self-reported race. Because of the oversampling of younger women, more than half of the study population was premenopausal. Compared to EA women, AA women were slightly younger at enrollment, had a significantly higher BMI and lower education attainment, and were less likely to be a never smoker and to have family history of breast cancer. For reproductive aging phenotypes, AA women had a similar AM, but younger ANM and shorter TRLS and NRLS than EA women. Further comparisons of the descriptive characteristics between breast cancer cases and controls within each racial group are shown in S2 Table. In both AA and EA groups, there were no significant differences in AM or ANM between cases and controls; whereas cases had longer TRLS and NRLS than controls (p≤0.016 except for NRLS in AAs).

thumbnail
Download:
Table 1. Descriptive characteristics of the Women’s Circle of Health Study.

https://doi.org/10.1371/journal.pone.0187205.t001

Age at menarche (AM)

In EA women, 31 out of 37 SNPs selected for association with AM were in the same direction as in the previous GWAS (p = 3.7e-6, binomial sign test), including 4 with a nominal significance of p≤0.05 and another 5 with a suggestive significance of p≤0.10 (Table 2). The top SNP, rs7759938 in LIN28B, was one of the most consistent variants in previous GWAS for AM. In our study, rs7759938 was associated with a 0.2 year of delay in menarche per copy of the variant G allele (p = 0.004). In addition, rs4843747 in 16q24.2–3, which was associated with later ANM in a previous GWAS, was shown to be significantly associated with a 0.2-year delay in AM per copy of the variant A allele (p = 0.01) in EA women in our study. In AA women, 20 out of 37 SNPs selected for association with AM were in the same direction as in the previous GWAS (p = 0.37, binomial sign test); yet only rs2947411 in TMEM18 reached a nominal significance (p = 0.03) in our study. rs757647 in KDM3B, which was associated with earlier AM in EA women in our study, was also significant in AA women but in an opposite direction.

thumbnail
Download:
Table 2. Top-ranked SNPs associated with age at menarche (p≤0.10).

https://doi.org/10.1371/journal.pone.0187205.t002

Age at natural menopause (ANM)

In EA women, 2 out of 16 SNPs selected for association with ANM were in the same direction as in the previous GWAS and reached a nominal significance of p≤0.05 (Table 3). The top SNP, rs16991615 in MCM8 was also one of the most significant SNPs identified in GWAS for ANM, with each copy of the variant A allele associated with 1.3-year delay in the onset of menopause. rs236114, another SNP in MCM8, was also replicated with ANM in EA women in our study with a suggestive level of significance (p = 0.08). In addition, two SNPs previously associated with AM, rs7821178 in PXMP3 and rs6589964 in BSX, were associated with ANM in EA women in our study. In AA women, none of the selected SNPs for association with ANM were in the same direction and reached a nominal significance. Instead, five SNPs previously associated with AM were significant with a nominal p≤0.10, including rs1859345 in SPOCK.

thumbnail
Download:
Table 3. Top-ranked SNPs associated with age at menopause (p≤0.10).

https://doi.org/10.1371/journal.pone.0187205.t003

Total reproductive lifespan (TRLS) and net reproductive lifespan (NRLS)

Tables 4 and 5 present the top ranked SNPs associated with TRLS and NRLS, respectively. It is noteworthy that a majority of those SNPs were associated with either AM or ANM in our study as described above, and most of them were associated with both TRLS and NRLS. The top SNP associated with TRLS was rs16991615 in MCM8 in EA women and rs1859345 in SPOCK in AA women, both of which were also associated with ANM and NRLS. Their effect sizes were similar across the three phenotypes (1.1 to 1.3 years per copy of the variant A allele for rs16991615 and -1.6 to -1.7 years per copy of the variant G allele). The top SNP associated with NRLS in EA was rs1079866 in INHBA, which was also associated with AM in our analysis. Per copy of the variant G allele was associated with a 0.9 year shorter in NRLS, whereas it was associated with only a 0.2-year delay in AM. Another noteworthy SNP is rs466639 in RXRG, which was associated with both TRLS and NRLS in both AA and EA women. Its effect size appeared to be slightly larger for NRLS (1.1 to 1.4 years per copy of the variant A allele) than for TRLS (0.8 to 1.0 year), the latter of which was similar to its effect size for ANM in AAs (0.8 year per copy).

thumbnail
Download:
Table 4. Top-ranked SNPs associated with total reproductive lifespan (p≤0.10).

https://doi.org/10.1371/journal.pone.0187205.t004

thumbnail
Download:
Table 5. Top-ranked SNPs associated with net reproductive lifespan (p≤0.10).

https://doi.org/10.1371/journal.pone.0187205.t005

Polygenic risk scores (PRS) for reproductive aging phenotypes and breast cancer risk

The correlations between the PRS for each reproductive aging phenotype and the phenotype itself were highly significant as expected (p≤0.0002); yet the strength of the correlations was weak (correlation coefficient 0.16–0.21) (Table 6). None of the PRS was in significant association with breast cancer risk in either AA or EA women, with the exception that the third highest quartile of the PRS for NRSL was associated with an increased risk compared to the lowest quartile (odds ratio 2.19, 95% confident interval 1.14–4.23). Nevertheless, the trend test was not statistically significant (p = 0.10).

thumbnail
Download:
Table 6. Associations of polygenic scores of reproductive aging and breast cancer risk.

https://doi.org/10.1371/journal.pone.0187205.t006

Discussion

In a breast cancer case-control study with similar number of AA and EA women, we replicated several SNPs identified in previous GWAS of AM and ANM in the EA group, including the top GWAS signals rs7759938 in LIN28B and rs16991615 in MCM8. Nevertheless, the replication of the same set of SNPs in AA women was not as successful. Only one out of the 37 selected SNPs for AM was significant in the same direction, and none of the 16 selected SNPs for ANM was significant. In analyses of TRLS and NRLS, a number of SNPs associated with AN or ANM in previous GWAS were associated with the two reproductive lifespan phenotypes, and most of the associations were consistent between the two phenotypes with similar effect sizes. We found little evidence for a strong impact of those reproductive aging-associated SNPs on breast cancer risk. Our analysis suggests that the overlap of genetic determinants of reproductive aging phenotypes between AA and EA populations may be limited, and highlights the importance of evaluating GWAS signals initially identified from EA populations in AA populations, before they can be generalized across populations.

In the past decade, GWAS has been a powerhouse in discovering common genetic variants reliably associated with a wide range of complex phenotypes and common diseases. To date, the only GWAS on AM and ANM among AA women did not identify any variants at a genome-wide significant level (5e-8)[22, 23], which is a striking difference from the remarkable success of similar studies conducted among populations of European descent. An apparent reason is the scarcity of large epidemiological populations of African descent. Another reason may be related to their distinct genetic architecture. Compared to Europeans and East Asians, African populations have more unique genetic variations and shorter linkage disequilibrium (LD) blocks across the genome, features that require GWAS in African populations to have larger sample sizes. Due to these limitations, replication studies provide a feasible alternative approach to evaluate previous EA GWAS hits in people of African ancestry. Yet, in studies among AA women in the Population Architecture Using Genomics and Epidemiology (PAGE) Study, the generalization of the previously identified GWAS variants, including rs7759938 in LIN28B for AM and rs16991615 in MCM8 for ANM, was mostly unsuccessful[10, 21]. This is consistent with our data, where only one out of 37 SNPs for AM and none out of 16 SNPs for ANM was replicated in AA women. It should be noted that rs7759938 was associated with AM in AA women in a replication analysis nested in a previous GWAS[22], but it was not significant in either our study or in PAGE. The lack of generalization of GWAS variants for AM and ANM in AAs is in contrast to similar studies in East Asians and Hispanics, where a majority of the variants evaluated were replicated[42–44]. Considering the wide difference in allele frequency and LD between European and African populations, it is possible that genes or loci associated with AM and ANM are the same between the two populations, but index SNPs that are mostly likely to be identified in association studies may be different. Fine-mapping analyses of the entire regions of the previously identified loci, but not only the index SNPs, in the two GWAS among AA women provide some support for this hypothesis [22, 23]. However, the fact that none of the variants in AAs met the commonly accepted threshold for genome-wide significance calls for definitive evidence from future studies based on larger sample sizes.

Although most of the identified associations with AM and ANM in our study did not overlap between AA and EA women, a phenomenon we have repeatedly observed in previous studies of breast cancer risk[35–37, 45, 46], we found that rs466639 in RXRG was associated with TRLS and NRLS in both AA and EA women with a similar effect size. This SNP was initially identified in association with AM in Europeans, with an effect size of 4.2 months earlier in AM per copy of the variant T allele[25]. We replicated the association of this SNP with AM in EA women. Although it was not associated with AM in AA women, this SNP was associated with a delay in ANM in AA women but not in EA women. It is thus possible that its association with reproductive lifespan was driven by the effect on AM in EAs and the effect on ANM in AAs. This SNP is intronic in RXRG, which encodes retinoid X receptor gamma, a nuclear receptor that dimerizes with the receptors for retinoic acid, thyroid hormone and vitamin D and increases their DNA binding and transcriptional activity[47]. However, the biological impact of the SNP on gene expression and the mechanism of RXRG in regulating reproductive aging, remain unclear. In a recent study, rs466639 was associated with involution of terminal duct lobular unit (TDLU), the milk-producing structures of the breast[48]. The variant T allele associated with earlier AM was also associated with a higher count of TDLU, which is a recognized risk factor for breast cancer. In another study of AM SNPs with pubertal growth in male and female adolescents, rs466639 was not associated with changes in weight, height or BMI; nevertheless, four other AM SNPs were, suggesting shared genetic influences on onset of menstrual cycle and pubertal body growth[49].

To our knowledge, the literature on the overlapping of genetic underpinning between AM and ANM has been rather limited. In an early study in 1998 by Sneider et al., no correlation was found between age at menopause and age at menarche, suggesting distinct genetic mechanisms for the two phenotypes.[6] In another study by Spencer et al., two genes, ARHGAP42 and PHACTR1, which were associated with both AM and ANM; nevertheless, the associated SNPs within each gene were different between the two phenotypes. Although the directions of the associations were similar, those SNPs were not in strong LD.[21] Therefore, current data appear not to support large overlapping of genetic determinants for AM and ANM. Although the reasons are unknown, it could be due to the fact that the occurrence of ANM has a much longer time window susceptible to impact of various environmental factors in comparison to AM, and thus the identification of genetic determinants could be more difficult for the former. It may be interesting to note that far more GWAS have been done for AM than ANM. Clearly, future genetic association studies are warranted, especially for ANM.

Given the well-established associations of early AM and late ANM with breast cancer risk, it is a reasonable next step to examine whether GWAS SNPs for reproductive aging are associated with breast cancer risk, an approach known as Mendelian randomization analysis that is less likely to be biased by confounding factors and can be used to infer causality. In the first Mendelian randomization analysis for breast cancer risk, He and colleagues examined 19 AM SNPs and 17 ANM SNPs, individually and as PRS for early AM to combine the effects of multiple markers, in the ReproGen Consortium[11]. Two AM SNPs and one ANM SNP were associated with breast cancer in the expected direction. In addition, the 4th and 5th highest quintiles of PRS for AM were associated with increased breast cancer risk, with a borderline significant trend; yet no association was found with PRS for ANM. In a recent large study of PRS based on 375 autosomal AM variants, increasing AM was associated with lower risk of breast cancer, particularly estrogen receptor (ER) positive cancer, as well as ovarian cancer and endometrial cancer, which was independent of the effects of those variants on obesity[50]. The lack of associations of PRS for reproductive aging phenotypes with breast cancer risk in either AA or EA women in our study could be due to the limited number of GWAS variants included in our study. As a result, the proportion of phenotypic variations captured by those SNPs is likely low, and their effects on breast cancer risk might be too modest to be detected with our sample size.

We noted a few limitations in our study. Because genotyping was performed in 2011, which predated a number of large-scale GWAS for AM and ANM, our analysis included only variants identified from earlier studies and missed a substantial proportion of variants known to date. The sample size was also limited, especially for analysis of phenotypes among postmenopausal women. In addition, none of the SNPs for the four phenotypes passed the significance threshold (α = 0.05) after controlling for multiple comparisons via FDR correction. Therefore, our results should be interpreted with caution of small sample size. This is particularly relevant to our speculation of different genetic underpinning for AM and ANM based on the marked lower replication success rate in AAs than in EAs. Due to limited statistical power, we evaluated the concordance of the direction of associations between our study and previous GWAS literature, as an additional criterion for replication success, which is not stringent statistically speaking. However, our study is among the first to include both AA and EA women from a well-characterized study for breast cancer health disparities. Given the challenges in conducting fully powered GWAS among AA populations and the paucity of data available, our study contributes to the literature for genetics of reproductive aging phenotypes in AA women and highlights the importance of cross-race/ethnicity replication of GWAS variants. Another strength is the inclusion of reproductive lifespans in our analysis, which may have a higher impact on the risk of breast cancer than either AM or ANM alone.

In conclusion, in the WCHS, we replicated several GWAS variants for AM and ANM in EA women, but few were replicated in AA women. We identified rs466639 in RXRG associated with TRLS and NRLS in both AA and EA groups. However, none of the PRS for reproductive aging phenotypes were in significant association with breast cancer risk, which could be due to the limited number of GWAS variants included in our study, small sample size, and possibly low power. Future studies are warranted to investigate whether differences in reproductive aging-related genetic variants explain health disparities in breast cancer.

Supporting information

S1 Table. Summary of detailed information on selected SNPs.

https://doi.org/10.1371/journal.pone.0187205.s001

(PDF)

S2 Table. Descriptive characteristics of the Women's Circle of Health Study.

https://doi.org/10.1371/journal.pone.0187205.s002

(PDF)

References

  1. 1. Menarche, menopause, and breast cancer risk: individual participant meta-analysis, including 118 964 women with breast cancer from 117 epidemiological studies. The Lancet Oncology. 2012;13(11):1141–51. Epub 2012/10/23. pmid:23084519; PubMed Central PMCID: PMC3488186.
  2. 2. Morris DH, Jones ME, Schoemaker MJ, Ashworth A, Swerdlow AJ. Determinants of age at menarche in the UK: analyses from the Breakthrough Generations Study. British journal of cancer. 2010;103(11):1760–4. Epub 2010/11/04. pmid:21045834; PubMed Central PMCID: PMCPMC2994234.
  3. 3. Towne B, Czerwinski SA, Demerath EW, Blangero J, Roche AF, Siervogel RM. Heritability of age at menarche in girls from the Fels Longitudinal Study. American journal of physical anthropology. 2005;128(1):210–9. Epub 2005/03/22. pmid:15779076.
  4. 4. de Bruin JP, Bovenhuis H, van Noord PA, Pearson PL, van Arendonk JA, te Velde ER, et al. The role of genetic factors in age at natural menopause. Human reproduction (Oxford, England). 2001;16(9):2014–8. Epub 2001/08/31. pmid:11527915.
  5. 5. Murabito JM, Yang Q, Fox C, Wilson PW, Cupples LA. Heritability of age at natural menopause in the Framingham Heart Study. The Journal of clinical endocrinology and metabolism. 2005;90(6):3427–30. Epub 2005/03/17. pmid:15769979.
  6. 6. Snieder H, MacGregor AJ, Spector TD. Genes control the cessation of a woman's reproductive life: a twin study of hysterectomy and age at menopause. The Journal of clinical endocrinology and metabolism. 1998;83(6):1875–80. Epub 1998/06/17. pmid:9626112.
  7. 7. van Asselt KM, Kok HS, Pearson PL, Dubas JS, Peeters PH, Te Velde ER, et al. Heritability of menopausal age in mothers and daughters. Fertility and sterility. 2004;82(5):1348–51. Epub 2004/11/10. pmid:15533358.
  8. 8. Treloar SA, Do KA, Martin NG. Genetic influences on the age at menopause. Lancet (London, England). 1998;352(9134):1084–5. Epub 1998/11/03. pmid:9798581.
  9. 9. Cramer DW, Xu H, Harlow BL. Family history as a predictor of early menopause. Fertility and sterility. 1995;64(4):740–5. Epub 1995/10/01. pmid:7672145.
  10. 10. Carty CL, Spencer KL, Setiawan VW, Fernandez-Rhodes L, Malinowski J, Buyske S, et al. Replication of genetic loci for ages at menarche and menopause in the multi-ethnic Population Architecture using Genomics and Epidemiology (PAGE) study. Human reproduction (Oxford, England). 2013;28(6):1695–706. Epub 2013/03/20. pmid:23508249; PubMed Central PMCID: PMCPMC3657124.
  11. 11. He C, Chasman DI, Dreyfus J, Hwang SJ, Ruiter R, Sanna S, et al. Reproductive aging-associated common genetic variants and the risk of breast cancer. Breast cancer research: BCR. 2012;14(2):R54. Epub 2012/03/22. pmid:22433456; PubMed Central PMCID: PMCPMC3446388.
  12. 12. He C, Murabito JM. Genome-wide association studies of age at menarche and age at natural menopause. Molecular and cellular endocrinology. 2014;382(1):767–79. Epub 2012/05/23. pmid:22613007.
  13. 13. Liu YZ, Guo YF, Wang L, Tan LJ, Liu XG, Pei YF, et al. Genome-wide association analyses identify SPOCK as a key novel gene underlying age at menarche. PLoS genetics. 2009;5(3):e1000420. Epub 2009/03/14. pmid:19282985; PubMed Central PMCID: PMCPMC2652107.
  14. 14. Mbarek H, Steinberg S, Nyholt DR, Gordon SD, Miller MB, McRae AF, et al. Identification of Common Genetic Variants Influencing Spontaneous Dizygotic Twinning and Female Fertility. American journal of human genetics. 2016;98(5):898–908. Epub 2016/05/03. pmid:27132594; PubMed Central PMCID: PMCPMC4863559.
  15. 15. Stolk L, Perry JR, Chasman DI, He C, Mangino M, Sulem P, et al. Meta-analyses identify 13 loci associated with age at menopause and highlight DNA repair and immune pathways. Nature genetics. 2012;44(3):260–8. Epub 2012/01/24. pmid:22267201; PubMed Central PMCID: PMCPMC3288642.
  16. 16. Stolk L, Zhai G, van Meurs JB, Verbiest MM, Visser JA, Estrada K, et al. Loci at chromosomes 13, 19 and 20 influence age at natural menopause. Nature genetics. 2009;41(6):645–7. Epub 2009/05/19. pmid:19448619; PubMed Central PMCID: PMCPMC3000545.
  17. 17. Perry JR, Day F, Elks CE, Sulem P, Thompson DJ, Ferreira T, et al. Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche. Nature. 2014;514(7520):92–7. Epub 2014/09/19. pmid:25231870; PubMed Central PMCID: PMCPMC4185210.
  18. 18. 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. Nature genetics. 2009;41(6):729–33. Epub 2009/05/19. pmid:19448623; PubMed Central PMCID: PMCPMC3000552.
  19. 19. Pyun JA, Kim S, Cho NH, Koh I, Lee JY, Shin C, et al. Genome-wide association studies and epistasis analyses of candidate genes related to age at menarche and age at natural menopause in a Korean population. Menopause (New York, NY). 2014;21(5):522–9. Epub 2013/09/21. pmid:24045676.
  20. 20. Rahmani M, Earp MA, Ramezani Tehrani F, Ataee M, Wu J, Treml M, et al. Shared genetic factors for age at natural menopause in Iranian and European women. Human reproduction (Oxford, England). 2013;28(7):1987–94. Epub 2013/04/18. pmid:23592221.
  21. 21. Spencer KL, Malinowski J, Carty CL, Franceschini N, Fernandez-Rhodes L, Young A, et al. Genetic variation and reproductive timing: African American women from the Population Architecture using Genomics and Epidemiology (PAGE) Study. PloS one. 2013;8(2):e55258. Epub 2013/02/21. pmid:23424626; PubMed Central PMCID: PMCPMC3570525.
  22. 22. Demerath EW, Liu CT, Franceschini N, Chen G, Palmer JR, Smith EN, et al. Genome-wide association study of age at menarche in African-American women. Human molecular genetics. 2013;22(16):3329–46. Epub 2013/04/20. pmid:23599027; PubMed Central PMCID: PMCPMC3723312.
  23. 23. Chen CT, Liu CT, Chen GK, Andrews JS, Arnold AM, Dreyfus J, et al. Meta-analysis of loci associated with age at natural menopause in African-American women. Human molecular genetics. 2014;23(12):3327–42. Epub 2014/02/05. pmid:24493794; PubMed Central PMCID: PMCPMC4030781.
  24. 24. Sulem P, Gudbjartsson DF, Rafnar T, Holm H, Olafsdottir EJ, Olafsdottir GH, et al. Genome-wide association study identifies sequence variants on 6q21 associated with age at menarche. Nature genetics. 2009;41(6):734–8. Epub 2009/05/19. pmid:19448622.
  25. 25. Elks CE, Perry JR, Sulem P, Chasman DI, Franceschini N, He C, et al. Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies. Nat Genet. 2010;42(12):1077–85. Epub 2010/11/26. pmid:21102462; PubMed Central PMCID: PMCPMC3140055.
  26. 26. He C, Kraft P, Chen C, Buring JE, Pare G, Hankinson SE, et al. Genome-wide association studies identify loci associated with age at menarche and age at natural menopause. Nature genetics. 2009;41(6):724–8. Epub 2009/05/19. pmid:19448621; PubMed Central PMCID: PMCPMC2888798.
  27. 27. 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. Nature genetics. 2017;49(6):834–41. pmid:28436984.
  28. 28. Spratt DE, Chan T, Waldron L, Speers C, Feng FY, Ogunwobi OO, et al. Racial/Ethnic Disparities in Genomic Sequencing. JAMA oncology. 2016;2(8):1070–4. Epub 2016/07/02. pmid:27366979; PubMed Central PMCID: PMCPMC5123755.
  29. 29. Morabia A, Costanza MC. International variability in ages at menarche, first livebirth, and menopause. World Health Organization Collaborative Study of Neoplasia and Steroid Contraceptives. American journal of epidemiology. 1998;148(12):1195–205. Epub 1998/12/29. pmid:9867266.
  30. 30. Freedman DS, Khan LK, Serdula MK, Dietz WH, Srinivasan SR, Berenson GS. Relation of age at menarche to race, time period, and anthropometric dimensions: the Bogalusa Heart Study. Pediatrics. 2002;110(4):e43. Epub 2002/10/03. pmid:12359816.
  31. 31. Karapanou O, Papadimitriou A. Determinants of menarche. Reproductive biology and endocrinology: RB&E. 2010;8:115. Epub 2010/10/06. pmid:20920296; PubMed Central PMCID: PMCPMC2958977.
  32. 32. Bromberger JT, Matthews KA, Kuller LH, Wing RR, Meilahn EN, Plantinga P. Prospective study of the determinants of age at menopause. American journal of epidemiology. 1997;145(2):124–33. Epub 1997/01/15. pmid:9006309.
  33. 33. MacMahon B, Worcester J. Age at menopause. United States—1960-1962. Vital and health statistics Series 11, Data from the national health survey. 1966;(19):1–20. Epub 1966/10/01. pmid:5298639.
  34. 34. Johnson N, Dudbridge F, Orr N, Gibson L, Jones ME, Schoemaker MJ, et al. Genetic variation at CYP3A is associated with age at menarche and breast cancer risk: a case-control study. Breast cancer research: BCR. 2014;16(3):R51. Epub 2014/06/03. pmid:24887515; PubMed Central PMCID: PMCPMC4522594.
  35. 35. Quan L, Hong CC, Zirpoli G, Roberts MR, Khoury T, Sucheston-Campbell LE, et al. Variants of estrogen-related genes and breast cancer risk in European and African American women. Endocr Relat Cancer. 2014;21(6):853–64. Epub 2014/09/18. pmid:25228414; PubMed Central PMCID: PMC4214251.
  36. 36. Yao S, Graham K, Shen J, Campbell LE, Singh P, Zirpoli G, et al. Genetic variants in microRNAs and breast cancer risk in African American and European American women. Breast cancer research and treatment. 2013;141(3):447–59. Epub 2013/09/26. pmid:24062209.
  37. 37. Yao S, Zirpoli G, Bovbjerg DH, Jandorf L, Hong CC, Zhao H, et al. Variants in the vitamin D pathway, serum levels of vitamin D, and estrogen receptor negative breast cancer among African-American women: a case-control study. Breast cancer research: BCR. 2012;14(2):R58. Epub 2012/04/07. bcr3162 [pii] pmid:22480149
  38. 38. Ambrosone CB, Ciupak GL, Bandera EV, Jandorf L, Bovbjerg DH, Zirpoli G, et al. Conducting Molecular Epidemiological Research in the Age of HIPAA: A Multi-Institutional Case-Control Study of Breast Cancer in African-American and European-American Women. J Oncol. 2009;2009:871250. Epub 2009/10/30. pmid:19865486.
  39. 39. Bandera EV, Chandran U, Zirpoli G, McCann SE, Ciupak G, Ambrosone CB. Rethinking sources of representative controls for the conduct of case-control studies in minority populations. BMC medical research methodology. 2013;13:71. Epub 2013/06/01. pmid:23721229; PubMed Central PMCID: PMC3681587.
  40. 40. Ruiz-Narvaez EA, Rosenberg L, Wise LA, Reich D, Palmer JR. Validation of a small set of ancestral informative markers for control of population admixture in African Americans. American journal of epidemiology. 2011;173(5):587–92. Epub 2011/01/26. pmid:21262910; PubMed Central PMCID: PMCPMC3105437.
  41. 41. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59. pmid:10835412; PubMed Central PMCID: PMCPMC1461096.
  42. 42. Chen CT, Fernandez-Rhodes L, Brzyski RG, Carlson CS, Chen Z, Heiss G, et al. Replication of loci influencing ages at menarche and menopause in Hispanic women: the Women's Health Initiative SHARe Study. Human molecular genetics. 2012;21(6):1419–32. Epub 2011/12/02. pmid:22131368; PubMed Central PMCID: PMCPMC3284121.
  43. 43. Delahanty RJ, Beeghly-Fadiel A, Long JR, Gao YT, Lu W, Xiang YB, et al. Evaluation of GWAS-identified genetic variants for age at menarche among Chinese women. Human reproduction (Oxford, England). 2013;28(4):1135–43. Epub 2013/02/15. pmid:23406970; PubMed Central PMCID: PMCPMC3600840.
  44. 44. Shen C, Delahanty RJ, Gao YT, Lu W, Xiang YB, Zheng Y, et al. Evaluating GWAS-identified SNPs for age at natural menopause among chinese women. PloS one. 2013;8(3):e58766. Epub 2013/03/29. pmid:23536822; PubMed Central PMCID: PMCPMC3607593.
  45. 45. Cheng TY, Shankar J, Zirpoli G, Roberts MR, Hong CC, Bandera EV, et al. Genetic variants in the mTOR pathway and interaction with body size and weight gain on breast cancer risk in African-American and European American women. Cancer causes & control: CCC. 2016;27(8):965–76. Epub 2016/06/18. pmid:27314662; PubMed Central PMCID: PMCPMC4958603.
  46. 46. Gong Z, Yao S, Zirpoli G, David Cheng TY, Roberts M, Khoury T, et al. Genetic variants in one-carbon metabolism genes and breast cancer risk in European American and African American women. International journal of cancer. 2015;137(3):666–77. Epub 2015/01/20. pmid:25598430; PubMed Central PMCID: PMCPMC4437886.
  47. 47. Giner XC, Cotnoir-White D, Mader S, Levesque D. Selective ligand activity at Nur/retinoid X receptor complexes revealed by dimer-specific bioluminescence resonance energy transfer-based sensors. FASEB journal: official publication of the Federation of American Societies for Experimental Biology. 2015;29(10):4256–67. Epub 2015/07/08. pmid:26148973.
  48. 48. Oh H, Bodelon C, Palakal M, Chatterjee N, Sherman ME, Linville L, et al. Ages at menarche- and menopause-related genetic variants in relation to terminal duct lobular unit involution in normal breast tissue. Breast cancer research and treatment. 2016;158(2):341–50. Epub 2016/06/28. pmid:27342457; PubMed Central PMCID: PMCPMC5144736.
  49. 49. Tu W, Wagner EK, Eckert GJ, Yu Z, Hannon T, Pratt JH, et al. Associations between menarche-related genetic variants and pubertal growth in male and female adolescents. The Journal of adolescent health: official publication of the Society for Adolescent Medicine. 2015;56(1):66–72. Epub 2014/10/08. pmid:25287989; PubMed Central PMCID: PMCPMC4275397.
  50. 50. Day F, Thompson D, Helgason H, Chasman D, Finucane H, Sulem P, et al. Genomic analyses for age at menarche identify 389 independent signals and indicate BMI-independent effects of puberty timing on cancer susceptibility. bioRxiv. 2016.