Genome-Wide Association Analyses Identify SPOCK as a Key Novel Gene Underlying Age at Menarche

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.

[1]  P. Donnelly,et al.  A new multipoint method for genome-wide association studies by imputation of genotypes , 2007, Nature Genetics.

[2]  Hui Shen,et al.  A whole-genome linkage scan suggests several genomic regions potentially containing QTLs underlying the variation of stature. , 2002, American journal of medical genetics.

[3]  A. Silman,et al.  Determinants of incident vertebral fracture in men and women: results from the European Prospective Osteoporosis Study (EPOS) , 2003, Osteoporosis International.

[4]  Chia-Hung Liu,et al.  FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization , 2006, Nucleic Acids Res..

[5]  Hui Shen,et al.  The chemokine (C-C-motif) receptor 3 (CCR3) gene is linked and associated with age at menarche in Caucasian females , 2007, Human Genetics.

[6]  Frank Dudbridge,et al.  Likelihood-Based Association Analysis for Nuclear Families and Unrelated Subjects with Missing Genotype Data , 2008, Human Heredity.

[7]  A. Tsatsoulis,et al.  Combined estrogen receptor α and estrogen receptor β genotypes influence the age of menarche , 2006 .

[8]  Hui Shen,et al.  A genomewide linkage scan for quantitative-trait loci for obesity phenotypes. , 2002, American journal of human genetics.

[9]  P. Alliel,et al.  Regulated expression of the proteoglycan SPOCK in the neuromuscular system , 2000, Mechanisms of Development.

[10]  A. Silman Risk factors for Colles' fracture in men and women: results from the European Prospective Osteoporosis Study , 2003, Osteoporosis International.

[11]  W. Willett,et al.  A genome-wide association study identifies alleles in FGFR2 associated with risk of sporadic postmenopausal breast cancer , 2007, Nature Genetics.

[12]  Peter H. Westfall,et al.  Testing Association of Statistically Inferred Haplotypes with Discrete and Continuous Traits in Samples of Unrelated Individuals , 2002, Human Heredity.

[13]  W. Rand,et al.  Recall of early menstrual history and menarcheal body size: after 30 years, how well do women remember? , 2002, American journal of epidemiology.

[14]  Judy H Cho,et al.  Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis , 2007, Nature Genetics.

[15]  S. Golub,et al.  Recollections of menarche and women's subsequent experiences with menstruation. , 1983, Women & health.

[16]  A. Duncan,et al.  Testican in human blood. , 2001, Biochemical and biophysical research communications.

[17]  J. Lakins,et al.  Differential display RT PCR of total RNA from human foreskin fibroblasts for investigation of androgen-dependent gene expression. , 1996, American journal of medical genetics.

[18]  A. Tsatsoulis,et al.  Combined estrogen receptor alpha and estrogen receptor beta genotypes influence the age of menarche. , 2006, Human reproduction.

[19]  P. Peeters,et al.  Age at menarche and breast cancer risk in nulliparous women , 2004, Breast Cancer Research and Treatment.

[20]  D. Boomsma,et al.  The Familial Clustering of Age at Menarche in Extended Twin Families , 2007, Behavior genetics.

[21]  J. Ioannidis,et al.  Association of polymorphisms of the oestrogen receptor alpha gene with the age of menarche. , 2002, Human reproduction.

[22]  Allan Birnbaum,et al.  Combining Independent Tests of Significance , 1954 .

[23]  M. Jarvelin,et al.  A Common Variant in the FTO Gene Is Associated with Body Mass Index and Predisposes to Childhood and Adult Obesity , 2007, Science.

[24]  E. Greif,et al.  The psychological impact of menarche on early adolescent females: a review of the literature. , 1982, Child development.

[25]  Lester L. Peters,et al.  Genome-wide association study identifies novel breast cancer susceptibility loci , 2007, Nature.

[26]  Ivar Heuch,et al.  Association of low age at menarche with increased all-cause mortality: a 37-year follow-up of 61,319 Norwegian women. , 2007, American journal of epidemiology.

[27]  Dong Sun Kim,et al.  Expression profiling of clonal lymphocyte cell cultures from Rett syndrome patients , 2006, BMC Medical Genetics.

[28]  A. Duncan,et al.  Testican in Human Blood1 , 2001 .

[29]  V. Henderson,et al.  Estrogen deficiency and risk of Alzheimer's disease in women. , 1994, American journal of epidemiology.

[30]  E. Koff,et al.  Flashbulb memories of menarche and adult menstrual distress. , 1987, Journal of adolescence.

[31]  H. Sato,et al.  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. , 2001, Cancer research.

[32]  J. Long,et al.  The oestrogen receptor alpha gene is linked and/or associated with age of menarche in different ethnic groups. , 2005, Journal of medical genetics.

[33]  Jing Huang,et al.  Dynamic model based algorithms for screening and genotyping over 100K SNPs on oligonucleotide microarrays , 2005, Bioinform..

[34]  L. Giudice,et al.  Human endometrial matrix metalloproteinase-2, a putative menstrual proteinase. Hormonal regulation in cultured stromal cells and messenger RNA expression during the menstrual cycle. , 1996, The Journal of clinical investigation.

[35]  S. Levy,et al.  Genome-wide association study identifies two novel loci containing FLNB and SBF2 genes underlying stature variation. , 2009, Human molecular genetics.

[36]  P. Visscher,et al.  Estimation of Variance Components for Age at Menarche in Twin Families , 2007, Behavior genetics.

[37]  S. Narod,et al.  CYP gene polymorphisms and early menarche. , 2001, Molecular genetics and metabolism.

[38]  B. Weir Genetic Data Analysis II. , 1997 .

[39]  Judy H. Cho,et al.  A Genome-Wide Association Study Identifies IL23R as an Inflammatory Bowel Disease Gene , 2006, Science.

[40]  B. Schäfer,et al.  Isolation of genes differentially expressed in human primary myoblasts and embryonal rhabdomyosarcoma , 1996, International journal of cancer.

[41]  S. Heath,et al.  Weight-adjusted genome scan analysis for mapping quantitative trait Loci for menarchal age. , 2006, The Journal of clinical endocrinology and metabolism.

[42]  R. Recker,et al.  Polymorphism in the insulin-like growth factor 1 gene is associated with age at menarche in caucasian females. , 2007, Human reproduction.

[43]  C. Edgell,et al.  Endothelial cell expression of testican mRNA. , 1997, Endothelium : journal of endothelial cell research.

[44]  J. Ioannidis,et al.  Association of polymorphisms of the oestrogen receptor α gene with the age of menarche , 2002 .

[45]  P. Jollès,et al.  Testican, a multidomain testicular proteoglycan resembling modulators of cell social behaviour. , 1993, European journal of biochemistry.

[46]  P. Visscher,et al.  A genome-wide linkage scan for age at menarche in three populations of European descent. , 2008, The Journal of clinical endocrinology and metabolism.

[47]  A. Tsatsoulis,et al.  Association of SHBG gene polymorphism with menarche. , 2005, Molecular human reproduction.

[48]  R. Brenner,et al.  Testican-1, an inhibitor of pro-MMP-2 activation, is expressed in cartilage. , 2004, Osteoarthritis and cartilage.

[49]  H. Christensen,et al.  Association of adverse childhood experiences, age of menarche, and adult reproductive behavior: Does the androgen receptor gene play a role? , 2004, American journal of medical genetics. Part B, Neuropsychiatric genetics : the official publication of the International Society of Psychiatric Genetics.

[50]  J. Long,et al.  Genomewide linkage scan for quantitative trait loci underlying variation in age at menarche. , 2006, The Journal of clinical endocrinology and metabolism.

[51]  Philippe Froguel,et al.  Genome-wide association scans identified CTNNBL1 as a novel gene for obesity. , 2008, Human molecular genetics.

[52]  R. Recker,et al.  Polymorphisms of estrogen-biosynthesis genes CYP17 and CYP19 may influence age at menarche: a genetic association study in Caucasian females. , 2006, Human molecular genetics.

[53]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[54]  Mark Daly,et al.  Haploview: analysis and visualization of LD and haplotype maps , 2005, Bioinform..

[55]  F. Hirahara,et al.  Estrogen-metabolizing gene polymorphisms, but not estrogen receptor-alpha gene polymorphisms, are associated with the onset of menarche in healthy postmenopausal Japanese women. , 2003, The Journal of clinical endocrinology and metabolism.

[56]  A. Duncan,et al.  Distribution of testican expression in human brain , 2000, Cell and Tissue Research.

[57]  R. Kaaks,et al.  Obesity , Endogenous Hormones , and Endometrial Cancer Risk : A Synthetic Review 1 , 2002 .

[58]  Terence P. Speed,et al.  Genome analysis A genotype calling algorithm for affymetrix SNP arrays , 2005 .

[59]  P. Xiao,et al.  The oestrogen receptor α gene is linked and/or associated with age of menarche in different ethnic groups , 2005, Journal of Medical Genetics.

[60]  K. Roeder,et al.  Genomic Control for Association Studies , 1999, Biometrics.

[61]  John D. Storey,et al.  Statistical significance for genomewide studies , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[62]  N M Laird,et al.  Family-based tests of association in the presence of linkage. , 2000, American journal of human genetics.

[63]  P. Donnelly,et al.  Inference of population structure using multilocus genotype data. , 2000, Genetics.