For females, menarche is a most significant physiological event. Age at menarche (AAM) is a trait with high genetic determination and is associated with major complex diseases in women. However, specific genes for AAM variation are largely unknown. To identify genetic factors underlying AAM variation, a genome-wide association study (GWAS) examining about 380,000 SNPs was conducted in 477 Caucasian women. A follow-up replication study was performed to validate our major GWAS findings using two independent Caucasian cohorts with 854 siblings and 762 unrelated subjects, respectively, and one Chinese cohort of 1,387 unrelated subjects—all females. Our GWAS identified a novel gene, SPOCK (Sparc/Osteonectin, CWCV, and Kazal-like domains proteoglycan), which had seven SNPs associated with AAM with genome-wide false discovery rate (FDR) q<0.05. Six most significant SNPs of the gene were selected for validation in three independent replication cohorts. All of the six SNPs were replicated in at least one cohort. In particular, SNPs rs13357391 and rs1859345 were replicated both within and across different ethnic groups in all three cohorts, with p values of 5.09×10？3 and 4.37×10？3, respectively, in the Chinese cohort and combined p values (obtained by Fisher's method) of 5.19×10？5 and 1.02×10？4, respectively, in all three replication cohorts. Interestingly, SPOCK can inhibit activation of MMP-2 (matrix metalloproteinase-2), a key factor promoting endometrial menstrual breakdown and onset of menstrual bleeding. Our findings, together with the functional relevance, strongly supported that the SPOCK gene underlies variation of AAM.
References
[1]
Kaaks R, Lukanova A, Kurzer MS (2002) Obesity, endogenous hormones, and endometrial cancer risk: a synthetic review. Cancer Epidemiol Biomarkers Prev 11: 1531–1543.
[2]
Peeters PH, Verbeek AL, Krol A, Matthyssen MM, de Waard F (1995) Age at menarche and breast cancer risk in nulliparous women. Breast Cancer Res Treat 33: 55–61.
[3]
Paganini-Hill A, Henderson VW (1994) Estrogen deficiency and risk of Alzheimer's disease in women. Am J Epidemiol 140: 256–261.
[4]
Silman AJ (2003) Risk factors for Colles' fracture in men and women: results from the European Prospective Osteoporosis Study. Osteoporos Int 14: 213–218.
[5]
Roy DK, O'Neill TW, Finn JD, Lunt M, Silman AJ, et al. (2003) Determinants of incident vertebral fracture in men and women: results from the European Prospective Osteoporosis Study (EPOS). Osteoporos Int 14: 19–26.
[6]
Jacobsen BK, Heuch I, Kvale G (2007) Association of low age at menarche with increased all-cause mortality: a 37-year follow-up of 61,319 Norwegian women. Am J Epidemiol 166: 1431–1437.
[7]
van den Berg SM, Boomsma DI (2007) The familial clustering of age at menarche in extended twin families. Behav Genet 37: 661–667.
[8]
Anderson CA, Duffy DL, Martin NG, Visscher PM (2007) Estimation of variance components for age at menarche in twin families. Behav Genet 37: 668–677.
[9]
Anderson CA, Zhu G, Falchi M, van den Berg SM, Treloar SA, et al. (2008) A genome-wide linkage scan for age at menarche in three populations of European descent. J Clin Endocrinol Metab 93: 3965–3970.
[10]
Long JR, Xu H, Zhao LJ, Liu PY, Shen H, et al. (2005) The oestrogen receptor alpha gene is linked and/or associated with age of menarche in different ethnic groups. J Med Genet 42: 796–800.
[11]
Stavrou I, Zois C, Ioannidis JP, Tsatsoulis A (2002) Association of polymorphisms of the oestrogen receptor alpha gene with the age of menarche. Hum Reprod 17: 1101–1105.
[12]
Stavrou I, Zois C, Chatzikyriakidou A, Georgiou I, Tsatsoulis A (2006) Combined estrogen receptor alpha and estrogen receptor beta genotypes influence the age of menarche. Hum Reprod 21: 554–557.
[13]
Xita N, Tsatsoulis A, Stavrou I, Georgiou I (2005) Association of SHBG gene polymorphism with menarche. Mol Hum Reprod 11: 459–462.
[14]
Jorm AF, Christensen H, Rodgers B, Jacomb PA, Easteal S (2004) Association of adverse childhood experiences, age of menarche, and adult reproductive behavior: does the androgen receptor gene play a role? Am J Med Genet B Neuropsychiatr Genet 125: 105–111.
[15]
Zhao J, Xiong DH, Guo Y, Yang TL, Recker RR, et al. (2007) Polymorphism in the insulin-like growth factor 1 gene is associated with age at menarche in caucasian females. Hum Reprod 22: 1789–1794.
[16]
Yang F, Xiong DH, Guo Y, Shen H, Xiao P, et al. (2007) The chemokine (C-C-motif) receptor 3 (CCR3) gene is linked and associated with age at menarche in Caucasian females. Hum Genet 121: 35–42.
[17]
Guo Y, Xiong DH, Yang TL, Guo YF, Recker RR, et al. (2006) Polymorphisms of estrogen-biosynthesis genes CYP17 and CYP19 may influence age at menarche: a genetic association study in Caucasian females. Hum Mol Genet 15: 2401–2408.
[18]
Gorai I, Tanaka K, Inada M, Morinaga H, Uchiyama Y, et al. (2003) Estrogen-metabolizing gene polymorphisms, but not estrogen receptor-alpha gene polymorphisms, are associated with the onset of menarche in healthy postmenopausal Japanese women. J Clin Endocrinol Metab 88: 799–803.
[19]
Lai J, Vesprini D, Chu W, Jernstrom H, Narod SA (2001) CYP gene polymorphisms and early menarche. Mol Genet Metab 74: 449–457.
[20]
Guo Y, Shen H, Xiao P, Xiong DH, Yang TL, et al. (2006) Genomewide linkage scan for quantitative trait loci underlying variation in age at menarche. J Clin Endocrinol Metab 91: 1009–1014.
[21]
Rothenbuhler A, Fradin D, Heath S, Lefevre H, Bouvattier C, et al. (2006) Weight-adjusted genome scan analysis for mapping quantitative trait Loci for menarchal age. J Clin Endocrinol Metab 91: 3534–3537.
[22]
Hunter DJ, Kraft P, Jacobs KB, Cox DG, Yeager M, et al. (2007) A genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancer. Nat Genet.
[23]
Easton DF, Pooley KA, Dunning AM, Pharoah PD, Thompson D, et al. (2007) Genome-wide association study identifies novel breast cancer susceptibility loci. Nature.
[24]
Rioux JD, Xavier RJ, Taylor KD, Silverberg MS, Goyette P, et al. (2007) Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis. Nat Genet 39: 596–604.
[25]
Duerr RH, Taylor KD, Brant SR, Rioux JD, Silverberg MS, et al. (2006) A genome-wide association study identifies IL23R as an inflammatory bowel disease gene. Science 314: 1461–1463.
[26]
Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, et al. (2007) A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science 316: 889–894.
[27]
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.
[28]
Devlin B, Roeder K (1999) Genomic control for association studies. Biometrics 55: 997–1004.
[29]
Fisher RA (1948) Combining independent tests of significance. American Statistician 2: 30.
[30]
Alliel PM, Perin JP, Jolles P, Bonnet FJ (1993) Testican, a multidomain testicular proteoglycan resembling modulators of cell social behaviour. Eur J Biochem 214: 347–350.
[31]
Marr HS, BaSalamah MA, Bouldin TW, Duncan AW, Edgell CJ (2000) Distribution of testican expression in human brain. Cell Tissue Res 302: 139–144.
[32]
Hausser HJ, Decking R, Brenner RE (2004) Testican-1, an inhibitor of pro-MMP-2 activation, is expressed in cartilage. Osteoarthritis Cartilage 12: 870–877.
[33]
Marr HS, BaSalamah MA, Edgell CJ (1997) Endothelial cell expression of testican mRNA. Endothelium 5: 209–219.
[34]
Genini M, Schwalbe P, Scholl FA, Schafer BW (1996) Isolation of genes differentially expressed in human primary myoblasts and embryonal rhabdomyosarcoma. Int J Cancer 66: 571–577.
[35]
Nitsche EM, Moquin A, Adams PS, Guenette RS, Lakins JN, et al. (1996) Differential display RT PCR of total RNA from human foreskin fibroblasts for investigation of androgen-dependent gene expression. Am J Med Genet 63: 231–238.
[36]
Delgado IJ, Kim DS, Thatcher KN, LaSalle JM, Van dVI (2006) Expression profiling of clonal lymphocyte cell cultures from Rett syndrome patients. BMC Med Genet 7: 61.
[37]
Cifuentes-Diaz C, Alliel PM, Charbonnier F, de la PS, Molgo J, et al. (2000) Regulated expression of the proteoglycan SPOCK in the neuromuscular system. Mech Dev 94: 277–282.
[38]
BaSalamah MA, Marr HS, Duncan AW, Edgell CJ (2001) Testican in human blood. Biochem Biophys Res Commun 283: 1083–1090.
[39]
Nakada M, Yamada A, Takino T, Miyamori H, Takahashi T, et al. (2001) Suppression of membrane-type 1 matrix metalloproteinase (MMP)-mediated MMP-2 activation and tumor invasion by testican 3 and its splicing variant gene product, N-Tes. Cancer Res 61: 8896–8902.
[40]
Irwin JC, Kirk D, Gwatkin RB, Navre M, Cannon P, et al. (1996) Human endometrial matrix metalloproteinase-2, a putative menstrual proteinase. Hormonal regulation in cultured stromal cells and messenger RNA expression during the menstrual cycle. J Clin Invest 97: 438–447.
[41]
Liu YJ, Liu XG, Wang L, Dina C, Yan H, et al. (2008) Genome-wide association scans identified CTNNBL1 as a novel gene for obesity. Hum Mol Genet 17: 1803–1813.
[42]
Lei SF, Tan LJ, Liu XG, Wang L, Yan H, et al. (2008) Genome-wide association study identifies two novel loci containing FLNB and SBF2 genes underlying stature variation. Hum Mol Genet.
[43]
Deng HW, Deng H, Liu YJ, Liu YZ, Xu FH, et al. (2002) A genomewide linkage scan for quantitative-trait loci for obesity phenotypes. Am J Hum Genet 70: 1138–1151.
[44]
Deng HW, Xu FH, Liu YZ, Shen H, Deng H, et al. (2002) A whole-genome linkage scan suggests several genomic regions potentially containing QTLs underlying the variation of stature. Am J Med Genet 113: 29–39.
[45]
Must A, Phillips SM, Naumova EN, Blum M, Harris S, et al. (2002) Recall of early menstrual history and menarcheal body size: after 30 years, how well do women remember? Am J Epidemiol 155: 672–679.
[46]
Golub S, Catalano J (1983) Recollections of menarche and women's subsequent experiences with menstruation. Women Health 8: 49–61.
[47]
Greif EB, Ulman KJ (1982) The psychological impact of menarche on early adolescent females: a review of the literature. Child Dev 53: 1413–1430.
[48]
Pillemer DB, Koff E, Rhinehart ED, Rierdan J (1987) Flashbulb memories of menarche and adult menstrual distress. J Adolesc 10: 187–199.
[49]
Di X, Matsuzaki H, Webster TA, Hubbell E, Liu G, et al. (2005) Dynamic model based algorithms for screening and genotyping over 100 K SNPs on oligonucleotide microarrays. Bioinformatics 21: 1958–1963.
[50]
Rabbee N, Speed TP (2006) A genotype calling algorithm for affymetrix SNP arrays. Bioinformatics 22: 7–12.
[51]
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.
[52]
Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, et al. (2002) Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered 53: 79–91.
[53]
Weir B (2007) Genetic data analysis III. Sunderland, MA: Sinauer Associates.
[54]
Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.
[55]
Storey JD, Tibshirani R (2003) Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 100: 9440–9445.
[56]
Lake SL, Blacker D, Laird NM (2000) Family-based tests of association in the presence of linkage. Am J Hum Genet 67: 1515–1525.
[57]
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.
[58]
Dudbridge F (2008) Likelihood-based association analysis for nuclear families and unrelated subjects with missing genotype data. Hum Hered 66: 87–98.
[59]
Yuan HY, Chiou JJ, Tseng WH, Liu CH, Liu CK, et al. (2006) FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization. Nucleic Acids Res 34: W635–W641.
[60]
Marchini J, Howie B, Myers S, McVean G, Donnelly P (2007) A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet 39: 906–913.
