Physiological Reports ORIGINAL RESEARCH Physiological Reports ISSN 2051-817X Epigenetic heredity of human height Pasquale Simeone & Saverio Alberti Unit of Cancer Pathology, Department of Neuroscience and Imaging and CeSI, University "G. d'Annunzio" Foundation, Chieti Scalo, Italy Keywords Bioinformatics, DNA methylation, epigenetic regulation, genome-wide association studies, meta-analysis. Correspondence Saverio Alberti, Unit of Cancer Pathology, Center for Excellence in Research on Aging, University "G. D' Annunzio", Via Colle dell' Ara, 66100 Chieti Scalo (Chieti), Italy. Tel: (+39) 0871-541.551 Fax: (+39) 0871-541.550 E-mail: s.alberti@unich.it Funding Information We are grateful for the support of the Fondazione Cassa di Risparmio della Provincia di Chieti, the Fondazione Compagnia di San Paolo and Oncoxx Biotech. Received: 18 April 2014; Revised: 19 May 2014; Accepted: 19 May 2014 doi: 10.14814/phy2.12047 Physiol Rep, 2 (6), 2014, e12047, doi: 10.14814/phy2.12047 Abstract Genome-wide SNP analyses have identified genomic variants associated with adult human height. However, these only explain a fraction of human height variation, suggesting that significant information might have been systematically missed by SNP sequencing analysis. A candidate for such non-SNP-linked information is DNA methylation. Regulation by DNA methylation requires the presence of CpG islands in the promoter region of candidate genes. Seventy two of 87 (82.8%), height-associated genes were indeed found to contain CpG islands upstream of the transcription start site (USC CpG island searcher; validation: UCSC Genome Browser), which were shown to correlate with gene regulation. Consistent with this, DNA hypermethylation modules were detected in 42 height-associated genes, versus 1.5% of control genes (P = 8.0199e~17), as were dynamic methylation changes and gene imprinting. Epigenetic heredity thus appears to be a determinant of adult human height. Major findings in mouse models and in human genetic diseases support this model. Modulation of DNA methylation are candidate to mediate environmental influence on epigenetic traits. This may help to explain progressive height changes over multiple generations, through trans-generational heredity of progressive DNA methylation patterns. Introduction Genomic loci linked to human height have been identified by genome-wide SNP-association analyses (GWAS) (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon et al. 2008). The findings from different studies overlap substantially (Table 1, Table SI), strongly supporting identification of a significant fraction of height-determining genes. Cho et al. (2009) performed a corresponding GWAS in a cohort of Korean descent. Notably, fifteen of the genes identified (ACAN, BCAS3 also known as TBX2, EFEMP1, HHIP, HMGA1, HMGA2, LCORL, NCAPG, PLAGL1, PTCH1, SOCS2, SPAG1, UQCC also known as GDF5, ZBTB38, ZNF678) were also identified in the previous studies (Weedon and Frayling 2008; Cho et al. 2009). Correspondingly, the Caucasian height-associated genes were confirmed as such in cohorts of Japanese descent (Okada et al. 2010). These findings showed that the same genes are largely associated with height in populations of both Caucasian and non-Caucasian ancestry, and that inherited height-controlling mechanisms are conserved. Genes close to the SNP most strongly associated with body size were shown to encode extracellular matrix components, proteases, cell cycle controllers, transcription factors and signaling molecules [(Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon et al. 2008) and manuscript in prep.]. Intriguingly, although, up to 90% of the variation in adult height may be explained by genetic factors (Silventoinen et al. 2003; Weedon and Frayling 2008), © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of 2014 | Vol. 2 | Iss. 6 | e12047 the American Physiological Society and The Physiological Society. Page 1 This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited Table 1. CpG islands in human height-associated gene promoters. Upstream TSS1 Genome sequence RefSeq CpG % obsCpG/ Length assembly nt upstream Genes number Chromosome TSS (nt) islands G+C expCpG (nt) 5' (nt) 3' (nt) TSS Notes2 ACAN (Weedon et al. 2008; NM_013227 15 87147678 1 65.7 0.72 1485 87146790 87148274 888 I Cho et al. 2009) ADAMTSL3 (Gudbjartsson NM_207517 15 82113842 1 68,7 0,82 1974 82113076 82115050 766 I et al. 2008; Lettre et al. 2008; Weedon et al. 2008) ADAMTS17 (Gudbjartsson NM_139057 15 98699706 2 60,7 0,65 888 9908373 9907485 8667 * et al. 2008) 69,1 0,89 2026 98700323 98698297 617 i AGPAT6 (Lettre et al. 2008) NMJ78819 8 41554876 None ANAPC13 (Weedon et al. NM_155391 3 135687519 1 59,2 0,788 1548 135688471 135686923 952 i 2008) ANKS1 (Gudbjartsson et al. NM_015245 6 34965016 2 55 0,80 660 349633552 349634012 1664 2008) 61,7 0,78 1497 34964519 34966016 497 i ATAD5 (Gudbjartsson et al. NM_024857 17 26183149 2 62,6 0,78 1368 26175115 26176483 8031 * 2008) 56,3 0,812 980 26182826 26183806 320 i ATXN3 (Gudbjartsson et al. NM_030660 14 91642718 1 60,4 0,744 1226 91643223 91641997 516 i 2008) BCAS3 (Gudbjartsson et al. NM_001099432 17 56109954 2 55 0,70 820 56104151 56104971 5803 * 2008; Cho et al. 2009) 59,5 0,65 1058 56109366 56110424 588 i BMP2 (Gudbjartsson et al. NM_001200 20 6696745 1 66,7 0,95 1988 6695759 6696745 986 i 2008) BMP6 (Gudbjartsson et al. NM_001718 6 7672010 1 66,8 0,87 2687 7670904 7673591 1106 i 2008) CABLES1 (Gudbjartsson et al. NM_001107404 18 18969725 1 67,4 0,95 2620 18968684 18971304 1041 i 2008) CDK6 (Gudbjartsson et al. NM_001259 7 92301148 2 61,2 0,748 2768 92304527 92301759 3379 2008; Lettre et al. 2008; 62,1 0,819 1951 92301759 92299808 611 i Weedon et al. 2008) CENTA2 (Gudbjartsson et al. NM_018404 17 26272880 2 55,1 0,69 532 26263024 26263556 9856 * 2008) 65,9 0,74 1687 26272611 26274298 269 i CHCHD7 (Gudbjartsson et al. NM_024300 8 57286869 1 63,7 0,893 1970 57285833 57287803 1036 i 2008; Lettre et al. 2008) COIL (Gudbjartsson et al. NM_004645 17 52393410 1 56,3 0,772 1404 52393685 52392281 275 i 2008) CPSF2 (Gudbjartsson et al. NM_175808 13 91658090 1 58 0,8 2129 91656710 91658839 1372 i 2008) CRLF3 (Gudbjartsson et al. NM_015986 17 26175904 2 55,9 0,817 1084 26183907 26182823 8003 * 2008) 62,6 0,775 1368 26176480 26175112 576 i P. Simeone & S. Alberti Epigenetic Regulation of Human Height £ E U uj Q Q N O CO '—CO rsi dJ c o co m O Q ^ i/i c t O "° 'IT > ^ 00 ^t" o m lti 01 m co ud o oi oi o oi m ui UD ^- rsi m Ln rsi ^- ld co m CO UD LD r~- UD rsi r-- ud o o m ■^t co o m r~- o m o rsi o r-- o rsi o co ^t co co rsi Ln '— oo m oi o rsi o ■^t co o m o m o rsi o LD '— '— '— '— LD k£) rsi rsi rsi rsi rsi ^- Ln Ln co ud '— r-- o O Ol 03 Ol Ol Lfl Ln ^- Ln co m m co m co ld ld ud o co ud r-- ud co co ud o" o" o" o" o" o" m rsi r~- rsi <— LD m LD o o co o — O CO t> CO o ^ ^ Ä 3 I § Q Q " _Q S 00 o o co r--LTl ^f- o oi co UD r~- m rsi ^- o Ln rsi m r-- ^t co >- rsi m oo ^ ° o ™ Q rsi <— rsi <— £ £ -2 o cu in 5 r 5 Z rsi ■^t ^f- <— rsi rsi oi rsi m co oi Ol rsi ■^t rsi ^ <-rsi oi m oi ud rsi LD m rsi co oo >- rsi r-~ o rsi r--<- rsi <— m rsi oi r-- m oi ^ oi r-- r~- m m ud co oi m m LD o o rsi <— ^ 'IT a; +-> cu lu rsi u ^ n3 00 o cu o Ol c rsi o o o 1/1 1/1 ~n5 rsi rt3 cu cu ~a cu Gu cu ho -—' u co" o o a. to rsi CQ 2S 53 00 ifl rsi ^~ rsi ^~ 00 00 CO ^1- ^— rsi <— ^— <— m ■^t rsi co ld co Ln o m ^- m oi co rsi r~- co i-n o ■^t oi Ol CO ■^t rsi co rsi ^— m '—rsi no 00 99 o o ° — 6S J~ ra o U o .. cu rsi co "a Ö oo 2 ° c? ° £ a; S! t S J3 S 3 "§ 2 S § « % ° § ^ dJ dJ 5 % © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society. 2014 I Vol. 2 I Iss. 6 I e12047 Page 3 Epigenetic Regulation of Human Height P. Simeone & S. Alberti UD rsi rsi UD co ^ <- m rn rsi m rsi ■^t co m cn rsi cn o ^ <- m m rsi m LTi co cn rsi m rsi ^ <- m m rsi m o cn rsi cn m co lti lti o o I I ™ § 8 § S a. qj i/i rsi ^ __. ^ Si f0 "o3 (0 ^ ^ rsi ■^t rsi rsi rsi <— m co cn m o rsi m cn m o rsi m rsi m ^ o m cn cn m o rsi m <— lti 00 £: +j1 cn O .EH, V ^ O T3 H O rsi w a; a: a: a: rsi m ■^t rsi cn LTi cn m r~- rsi ^- lti ud rsi ^f- co m ud o ^t- <- <- o 01 co >- <— ud r~- <— 00^0 o rsi o o o u 01 lti 00 co rsi m cn cn co o rsi cn m cn cn r--rsi m co 0000 N CTi CTi o r--LTi rsi rsi LTi LTi rsi ^- cn o ro 10 lh lti o <— lti lti rsi lti <— O p O rsi sj rsi ^, u ^1 u ^ < o O o R CO r~- r~- rsi m cn rsi m co o rsi rsi ud co ud ^f- cn rsi i_n i_n LTi m '— rsi <— lti m ^f- ud ^f- ud ^f- co cn 00 n n m st o m cn r~- co m o rn t— 00 lti rsi o lti ud ud LTi cn rsi rsi en o o o r~- LTi i_n en ud ud ud r~- r~- m n m ^ ui in 01 \t st en rsi Lfi st LTI rs CO CO CTl O ID co m o i/i m is (t| ^f- rsi <— lh m rsi ud ud ud ud cn rsi rsi cn o o o r~- LTi i_n cTi o m m ud ud ud cn ud m ud LTi ud m LTi m m rsi <— rsi <— cn '-co r~- LTi rsi LTi LTi o rsi ud ud r-- co ud co cn ö" ö" Ö" o" ö" ö" Ö" m LTi ^f- LTi co rsi <— ^— ^— m '—rsi ^~ lti ^~ § § - i/i rsi r- j5 ■ O n3 03 ^ ¥ ^ .i, Z5 O -Q So ^ o ■- S — cr co co S o o o ^ " O 0 0 rsi rsi rsj CO Q o Ü 3 g uj rs, uj rs, CO CO 2014 I Vol. Page 4 2 I Iss. 6 I e12047 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society P. Simeone & S. Alberti Epigenetic Regulation of Human Height £ E u o 60" r- 0 n~ O S; co r--r~- co <— LTi LTi LD LT) CI m rn r~- ld ld ud m en en en oi LD ld '— r-- m rn o rsi [n. co co cn r~- r- <- o LD ^- m rsi ld ud rsi 00 <- rsi m ^ o rsi ■FT -a S £ £ la1 Z5 Oj o o rsi <— <— O a: ^ u s: q. u X tNi CQ CQ *t CQ CO [>j co" o == o 0 rS 0 rsi ^ rsi ^ 00 5- ud m rsi ■^t co o o m rsi co m ■^t '—rsi UD O 03 1 03 ^ o SO o ^ 00 ^ S n fn rsi co m P S ^-1 ID ^ Cr, U I o 5 ™ st ID U1 ^- m co o rsi rsi LD rsi LD '— O m rsi rsi ld co LTl r~- en co co co co co co UD UD UD LD LD LD LD m cn o LD rsi m m id 03 CO CO CO CO CO CO UD UD UD LD LD LD ■^t O CO UD r~- o m m r~- rsi ld '— '— LD LD <— "a "a 00 j IN O O EN O >< O >< O > rsi fa rsi fa rsi on I—■ I—■ © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 I Vol. 2 I Iss. 6 I e12047 Page 5 -Q M QJ O ID -» cn— < o Table 1. Continued. Upstream TSS1 Genome sequence RefSeq CpG % obsCpG/ Length assembly nt upstream Genes number Chromosome TSS (nt) islands G+C expCpG (nt) 5' (nt) 3' (nt) TSS Notes2 TCP11 (Gudbjartsson et al. NM_018679 6 35217165 1 67,1 0,718 1134 35217617 35216483 452 I 2008) TGS1 (Gudbjartsson et al. NM_024831 8 56848345 1 59,5 0,8 1400 56847635 56849035 710 I 2008) TMED3 (Lettre et al. 2008) NM_007364 15 77390546 1 66,8 0,71 1133 77389985 77391118 561 I TRIM25 (Gudbjartsson et al. NM_005082 17 52346408 1 63,6 0,751 1460 52347061 52345601 653 i 2008) TRIP11 (Gudbjartsson et al. NM_004239 14 91576570 1 56,4 0,768 1310 91576732 91575422 593 i 2008; Lettre et al. 2008) UQCC/GDF5 (Gudbjartsson NM_018244 20 33463247 1 56,3 0,651 774 33463669 33462895 422 i et al. 2008; Lettre et al. 2008; Weedon et al. 2008; Cho et al. 2009) WDR60 (Lettre et al. 2008) NM_01851 7 158342030 1 60,7 0,82 1440 158341361 158342801 669 i ZBTB38 (Gudbjartsson et al. NM_001080412 3 142525745 None 2008; Lettre et al. 2008; Weedon et al. 2008; Cho et al. 2009) ZFHX4 (Gudbjartsson et al. NM_024721 8 77756070 3 55 0,78 540 77748585 77749125 7485 * 2008) 56,6 0,72 1064 77752435 77743499 3635 55 0,85 545 77755726 77746271 344 i ZNF76 (Gudbjartsson et al. NM_003427 6 35335488 1 64,5 0,78 1490 35334502 35335992 986 i 2008) ZNF462 (Gudbjartsson et al. NM_021224 9 108665199 2 55,4 0,72 1475 108661904 108663379 3295 2008) 61,9 0,65 507 108664296 108664803 903 ZNF678 (Weedon et al. NM_178549 1 225817867 1 66,4 0,66 1300 225814864 108666164 3003 2008; Cho et al. 2009) 1TSS: trascription start site. CpG islands were searched for in the region upstream the TSS, within an upper limit of - 10,000 bp. 2*CpG islands that overlap the TSS;*CpG islands in the region - 10,000 < -> -4000 bp upstream the TSS. P. Simeone & S. Alberti Epigenetic Regulation of Human Height stature-associated polymorphisms have been found to only explain between 2% and 3.7% of height variation (Gudbj-artsson et al. 2008; Lettre et al. 2008; Weedon and Frayling 2008; Weedon et al. 2008). More recent analyses have increased the combined predictive power of the identified traits (Yang et al. 2010). Nevertheless, a large fraction of heritable height-associated factors has escaped detection by conventional GWAS (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon and Frayling 2008; Weedon et al. 2008; Yang et al. 2010), consistent with difficulties of previous association studies in finding variants robustly associated with height (Lettre et al. 2007; Weedon et al. 2007; Sanna et al. 2008; Mackay et al. 2009). Notably, rather large cohorts were analyzed in the human height-SNP association studies (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon et al. 2008; Cho et al. 2009), which were designed to provide adequate power of detection (approximately 500,000 SNP spaced across the genome (Barrett and Cardon 2006)). Specific cohorts reached an up to 98% power to detect variants associated with >0.5% of the height variation (Lettre et al. 2008). The significant overlaps in findings obtained by independent groups (23/87, i.e., 26.4% of the genes, were identified by more than one study) (Table 1, Table SI), supported the analytical power of the studies performed. In other words, a significant coverage of height-associated variation appeared to have been reached, and a significant fraction of the most strongly height-associated loci was actually identified. Nevertheless, a large fraction of heritable height appeared unaccounted for, indicating that a corresponding large fraction of information might have been systematically missed by SNP-based GWAS. A candidate for such non-DNA sequence-linked information is epigenetic heredity. Strategies designed to detect sequence polymorphisms would, indeed, systematically miss genetic information that is not associated with DNA sequence changes. Functionally-relevant DNA methylation patterns were thus candidates to be associated with adult stature subgroups in addition to DNA sequence variants. Functionally-relevant DNA methylation patterns may affect selective mechanisms, thus behaving as true hereditary traits. Consistent with this, a metastable epigenetic heredity of the DWARF1 locus was shown to affect plant size (Miura et al. 2009), and this phenotype was inherited through mitosis and meiosis. Notably, environmental conditions, nutrition in particular, have been shown to affect height (Silventoi-nen et al. 2000). DNA methylation patterns can keep record of the nutritional status (El-Osta et al. 2008; Guer-rero-Bosagna et al. 2008) and affect, in turn, morphometric parameters (Guerrero-Bosagna et al. 2008). Modifications of DNA methylation patterns in growth-related genes can be inherited trans-generationally (Guerrero-Bosagna et al. 2008; Hollingsworth et al. 2008; Nadeau 2009; Roth et al. 2009; Braunschweig et al. 2012), through incomplete erasure of epigenetic patterning in the germline. An example of trans-generational epigenetic heredity of complex traits is that of longevity in Caenorhabditis elegans (Greer et al. 2011) . An accumulation of epigenetic changes through generations would then provide a valuable, reversible mechanism of adaptation to progressively changing environments (Nadeau 2009; Roth et al. 2009; Verginelli et al. 2009). These findings led us to assess the functional relevance of trait-associated DNA methylation patterns. Our results indicate that inheritance of CpG island methylation patterns may indeed be involved in the control of body development. They also suggest that environmental influence on height may be mediated by modulation of epigenetic heredity. This may help to account for progressive height changes over multiple generations, through trans-generational heredity of progressive DNA methylation patterns. Material and Methods CpG island - DNA methylation analysis Most genes regulated by DNA methylation contain one or more CpG islands, most frequently in their promoter region (Lander et al. 2001; Saxonov et al. 2006; Esteller 2008; Illingworth and Bird 2009; Jin et al. 2012). The presence of CpG islands in the 87 genes most strongly associated with the variation in human height (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon et al. 2008; Cho et al. 2009) was investigated using capabilities of the USC CpG island searcher (cpgislands.usc.edu/) (Takai and Jones 2002). The human genome sequence assembly 36.3 was used (www.ncbi.nlm.nih.gov/sites/entrez). To exclude GC-rich Alu repetitive elements, CpG island limits of >500 bp, >55% in G + C content and observed CpG/expected CpG >0.65 were imposed. Validation analyses were performed at the UCSC Genome Browser (genome.ucsc.edu/cgi-bin/ hgGateway) (Illingworth and Bird 2009). Relevance for disease pathological development was investigated using OMIM (Online Mendelian Inheritance in Man) resources (omim.org). Gene imprinting Genomic imprinting is the mechanism by which monoall-elic expression is achieved in a parent-specific fashion. Several human genes are known to be imprinted (Nazor et al. 2012) . These comprise epigenetic changes, when the imprints are established in the germline, or somatic changes, when they arise during early embryonic development. Genomic imprinting defects are associated with developmental disorders, including Silver-Russell, Beckwith-Wiedemann, and Prader-Willi syndromes. Geno- © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 I Vol. 2 I Iss. 6 I e12047 Page 7 Epigenetic Regulation of Human Height P. Simeone & S. Alberti mic imprints are affected by environmental factors, and also associate with several human cancers. Gene imprints were then analyzed for paternal or maternal patterns and gene expression regulation (www.geneimprint.org). DNA hypermethylation modules Association with DNA hypermethylation modules was evaluated as described (Easwaran et al. 2012). Briefly, at a genome-wide level, hypermethylated genes are marked by Polycomb complex signatures. Most of these genes comprise developmental regulators, which contain DNA hypermethylation modules, within bivalent, that is divergently regulated, chromatin. DNA methylation dynamics The status of DNA methylation of regulatory regions of height-associated genes was assessed to provide information on DNA methylation dynamics (www.methdb.de/; 202.97.205.78/diseasemeth/) (Lv et al. 2011). Regulation of height-associated genes by DNA methylation Mechanistic relevance of DNA methylation was assessed versus a regulatory role of promoter methylation on the expression levels of height-associated genes. Signaling networks of height-associated genes To detect potential associations to cancer growth control pathways, network-based analyses were designed: SNOW (studying networks in the omics world; babelo-mics.bioinfo.cipf.es/) builds protein-protein contact networks and maps lists of genes or proteins over a reference interactome, where nodes are proteins and edges are interaction events. A SNOW interactome (minimal connected networks, MCN) of height-associated genes was built using the HPRD (human protein reference database), IntAct, BIND (biomolecular interaction network database), DIP (database of interacting proteins) and MINT (molecular interaction database). Topological parameters included the distributions of node connections degrees (i.e., of the number of edges for each node); the betweenness (i.e., a measure of centrality of all nodes and of their distribution); the distribution of clustering coefficients (i.e., of the connectivity of the neighborhood of each node) and the number of components of the network (i.e., the different groups of nodes that are generated in a network analysis). MCNs were identified adopting the option of only one bridging protein between any pair of proteins analyzed. Statistical significance was computed by comparing the obtained networks with the entire reference interactome or with networks of random composition. All networks presented below have significantly higher betweenness and lower connections degrees as compared with irrelevant networks (P < 0.0001) (Fig. 1). MetaCore™ network analysis (http://lsresearch.thomson-reuters.com/pages/solutions/ 1/metacore) builds on a curated database. MetaCore™ is based on a proprietary manually curated database of human protein-protein, pro-tein-DNA and protein-com pound interactions, metabolic and signaling pathways for human, mouse and rat, supported by proprietary ontologies and controlled vocabulary. Over 2000 multi step canonical pathway maps of human protein-protein and protein-DNA interactions are utilized as reference maps. The analyze network TF feature was exploited to identify transcriptionally regulated pathways. The shortest pathway algorithm was used to link the height-associated genes with additional database objects along a directed path, enforcing the stringent option of one-intermediate-only. Ingenuity pathways analysis (Ingenuity Systems, www. ingenuity.com) was utilized to map height-associated genes onto the Ingenuity knowledge base (focus points). Networks of focus molecules were generated by maximizing specific connectivities. Networks were ranked by score (negative log P-value by right-tailed Fisher's exact test). Each score takes into account the number of molecules in the network, the final network size, as well as the dataset size and the total number of network-forming molecules in the Ingenuity Knowledge Base. Data meta-analysis (Tables S1-S3) validated identity and relationships between network members (www. signal -ing-gateway.org/; www.wikipathways.org/index.php/Wiki-Pathways and www.genome.jp/kegg/pathway.html). Yes/ no relevance for specific signaling pathways (SMAD; c-Myc; p53; ERcc) or processes (development, cell growth; apoptosis) was verified by wet-lab approaches (www.ncbi. nlm.nih.gov/sites/entrez; and manuscript in prep.). Statistical analysis Statistical significance of differences between groups was assessed by Fisher exact test (SISA, www.quantitativeskills. com/sisa/ and GraphPad Prism 6.02). Results DNA methylation of height-associated genes Most genes regulated by DNA methylation contain one or more CpG islands, most frequently in their promoter 2014 I Vol. 2 I Iss. 6 I e12047 Page 8 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society P. Simeone & S. Alberti Epigenetic Regulation of Human Height C VEGF Cyclin a Figure 1. Pathway analysis. Graphical representation of the height-associated proteins relationships retrieved through SNOW network analysis. Proteins are represented as nodes (hubs), the biological relationships between the nodes (edges) are represented as lines. Height-associated proteins are in red; linker proteins are in white; miRNA are in gray. Major hubs are in magenta; SMAD isoforms are in blue. region (Lander et al. 2001; Saxonov et al. 2006; Esteller 2008; Illingworth and Bird 2009). The presence of CpG islands in the genes most strongly associated with the variation in human height (87 loci) (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon et al. 2008; Cho et al. 2009) was investigated (Takai and Jones 2002) (Table 1). Remarkably, 72 of 87 height-associated genes (82.8%) were found to contain at least one CpG island in the 2,000 bp upstream of the transcription start site (TSS) (99 CpG islands overall) (Table 1). Notably, in all CpG islands-associated height genes, CpG islands overlapped with the TSS, supporting an actual regulatory role in gene transcription. As gene expression regulatory regions can significantly extend upstream of the TSS (Koudritsky and Domany 2008), we extended our analysis to the 2,000-4,000 bp upstream of the TSS (Table 1), and identified CpG islands in three additional genes (11 islands overall). More extensive investigation to 10,000 bp 5' from the TSS identified four additional genes as containing a CpG island (only 15 additional islands overall), suggesting adequate coverage of our search (Table 1 and Table SI). These observed frequencies (90.8% for 10 kb analysis) were strikingly higher than those measured across the whole genome (Jiang et al. 2007; Zhu et al. 2008). Housekeeping, that is universally expressed, genes contain CpG islands at higher-than-average frequency [78.7% (Larsen et al. 1992; Zhu et al. 2008)]. Thus, we repeated our analysis after removal of housekeeping height genes [five cases (Jiang et al. 2007; Zhu et al. 2008)] (Table 2). Non-housekeeping height genes were confirmed to host CpG islands at a much higher-than-expected frequency [87.8% vs. 45% for nonhousekeeping genes (Jiang et al. 2007)] (P = l.lOe-11; Fisher exact test). It should be noted that most tissue-specific/nonhousekeeping genes possess neither CpG-islands nor TATA-boxes in their core promoters, and only 19.2% of tissue-specific genes have a TATA-less CpG-island associated core promoter (Zhu et al. 2008), further supporting the specificity of the association of height genes to CpG islands. As an internal control, searches for CpG islands were conducted in the transcribed region of height-associated genes. The regulatory force of CpG islands in transcribed regions is lower than that of islands upstream the TSS (Nguyen et al. 2001; Weber et al. 2007). Consistently, only 35 height-associated genes (40.23% of the total) were found to contain CpG islands in transcribed segments (Table SI), which did not exceed the expected frequency. This supported the relevance of a quasi-universal presence of CpG islands in the promoter of height-associated genes. Independent searches for CpG islands were conducted with the UCSC Genome Browser, using broader/less stringent settings [length of >200 bp and for G+C content © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 I Vol. 2 I Iss. 6 I e12047 Page 9 Epigenetic Regulation of Human Height P. Simeone & S. Alberti Table 2. Functional relevance of DNA methylation of human height-associated genes. Gene symbol Gene name DNA methylation Notes' ACAN Aggrecan ADAMTSL3 ADAMTS like 3 Dunn et al. (2004) * ADAMTS17 ADAM metallopeptidase with thrombospondin type 1 motif 17 Dunn et al. (2004) * AGPAT6 1-acyglycerol-3-phosphate O-acyltransferase 6 ANAPC13 Anaphase promoting complex subunit 13 ANKS1 Ankyrin repeat and sterile alpha motif domain containing 1 ATAD5 ATPase family AAA domain containing 5 ATXN3 Ataxin 3 BCAS3 Breast carcinoma amplified sequence 3 BMP2 Bone morphogenic protein 2 Wen et al. (2006) BMP6 Bone morphogenic protein 6 Taniguchi et al. (2008) * CABLE S1 cdkB and Abl enzyme substrate 1 Sakamoto et al. (2008) *• CDK6 Cyclin-dependent kinase 6 CENTA2 Centaurin alpha 2 CHCHD7 Coiled-coil-helix-coiled-coil-helix domain containing 7 COIL Coilin CPSF2 Cleavage and polyadenylation specific factor 2 • CRLF3 Cytokine receptor like factor 3 DCC Deleted in colon carcinoma Park et al. (2008) * DEF6 Differentially expressed in FDCP6 homolog DGKE Diacylglycerol kinase epsilon 64 kDa DLEU7 Deleted in lymphocytic leukemia Hammarsund et al. (2004) * DNM3 Dynamin 3 DNMT3A DNA methyl transferase 3 alpha Esteller (2008) t D0J1L DOT1-like,histone H3 methyltransferase Jones et al. (2008) Í DYM Dymeclin EFEMP1 EGF containing fibulin like extracellular matrix protein 1 Sadr-Nabavi et al. (2009) * E4F1 E4F transcription factor 1 FBLN5 Fibulin5 FUBP3 Far upstream element (FUSE)-binding protein 3 GATAD1 GATA zinc finger domain containing 1 GLT25D2 Glycosyltransferase 25 domain containing 2 GNA12 Guanine nucleotide-binding protein (G protein) alpha 12 • GPR126 G protein coupled receptor 126 GRB10 Growth factor receptor bound protein 10 Li et al. (2008) HHIP Hedgehog interacting protein Tada et al. (2008) * HIST1H1D Histone H1D Albig et al. (1993) i HMGA1 High mobility group AT-hook 1 Sgarra et al. (2003) i HMGA2 High mobility group AT-hook 2 Reeves and Beckerbauer (2001) í IHH Indian hedgehog LBH Limb bud and heart development homolog LCORL Ligand-dependent nuclear receptor corepressor-like protein LIN28B Lin 28 homolog B LYAR Ly1 antibody reactive homolog LYN v-yes-1 Yamaguchi sarcoma viral related oncogene homolog MOS v-mos Moloney murine sarcoma viral oncogene homolog Scholz et al. (2005) * MTMR11 Myotubularin-related protein NACA2 Nascent polypeptide-associated complex alpha subunit 2 NCAPG Non-SMC condensin I complex subunit G NOG Noggin NKX2-1/TTF1 NK2 homeobox/thyroid transcription factor 1 Kondo et al. (2009) * PAPPA Pregnancy-associated plasma protein A PENK Proenkephalin Comb and Goodman (1990) * PEX1 Peroxisome biogenesis factor 1 PLAGL1 Pleomorphic adenoma gene 1 Arima and Wake (2006) * 2014 I Vol. 2 I Iss. 6 I e12047 Page 10 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society. P. Simeone & S. Alberti Table 2. Continued. Epigenetic Regulation of Human Height Gene symbol Gene name DNA methylation Notes1 PNPT1 Polyribonucleotide nucleotidyltransferase PRKG2 Protein kinase, cGMP-dependent, type II PTCH1 Patched homolog 1 (hedgehog signaling) Wolf et al. (2007) PXMP3 Peroxisomal membrane protein 3 35 kDa • RAB40C Member RAS oncogene family RBBP8 Retinoblastoma-binding protein 8 Li et al. (2014) RDHE2 Epidermal dehydrogenase 2 RNF135 Ring finger protein 135 RPS20 Ribosomal protein S20 SCMH1 Sex comb on midleg homolog 1 • SCUBE3 Signal peptide CUB-domain EGF-like 3 SF3B4 Splicing factor 3b subunit 4 49 kDa SH3GL3 SH3-domain GRB2-like 3 SOCS2 Suppressor of cytokine signaling 2 Sutherland et al. (2004) SPAG1 Sperm antigen 17 SV2A Synaptic vescicle glycoprotein 2A TBX2 T-box 2 TBX4 T-box 4 TCP11 t-complex protein 11 TGS1 Trimethylguanosine synthase homolog TMED3 Transmembrane emp24 protein transport-domain containing 3 TRIM25 Tripartite motif containing 25 TRIP11 Thyroid hormone receptor interactor 11 UQCOGDF5 Ubiquinol cytochrome C reductase chaperone/growth differentiation factor 5 WDR60 WD repeat domain 60 ZBTB38 Zinc finger and BTB domain containing 38 Filion et al. (2006) t ZFHX4 Zinc finger homeobox 4 ZNF76 Zinc finger protein 76 ZNF462 Zinc finger protein 462 ZNF678 Zinc finger protein 678 1*Genes regulated by DNA methylation; *Genes regulating DNA methylation; 'Housekeeping genes. of >50% (Illingworth and Bird 2009)] (Table S2). Largely concordant results were obtained with the USC CpG Island Searcher and the UCSC Genome Browser algorithms. The highest concordance was found in the promoter region (84% of the cases), consistent with the actual identification of bona fide transcriptionally active CpG islands. Gene imprinting Gene imprinting was analyzed for paternal or maternal patterns and regulation of gene expression (www.geneim-print.org). Listed in Table S3 are the height-associated genes/gene clusters that were found associated with imprinting (Nazor et al. 2012). Twenty five height-associated genes were found to be imprinted. Additionally, imprinting of GNAS was shown to lead to either gynoge-netic or androgenetic patterns. Growth factor receptor bound protein 10 (GRB10) expression was shown an imprinted gene (Li et al. 2008). PLAGL1 lies within a DNA methylation imprinting center on human chromosome region 6q24. Imprinting of PLAGL1 was also shown by expression and phenotypic profiling of parthenogenetic fetuses (Bischoff et al. 2009). Hypomethylation patterns were demonstrated for both PLAGL1 and GRB10. Using Affymetrix GeneChip microarrays (Carletti et al. 2006) and/or semiquantitative PCR, organs expression patterns were profiled, and PLAGL1 was found expressed in a tissue-specific manner; imprinting was confirmed by pyrosequencing. A paternally imprinted regulatory role was also shown, as paternal noncoding RNAs at this locus can interact with regulators of active chromatin (Iglesias-Platas et al. 2012). DNA hypermethylation modules The functional relevance of DNA methylation of height-associated genes was first assessed by association with DNA hypermethylation modules (Easwaran et al. 2012) (Table S3), as marked by Polycomb-complex signatures © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 I Vol. 2 I Iss. 6 I e12047 Page 11 Epigenetic Regulation of Human Height within bivalent chromatin. Most of Polycomb-signed genes comprise developmental regulators, and were shown to contain DNA hypermethylation modules (Easwaran et al. 2012). As compared to control H3K4Me3 genes, which are hypermethylated in 1.5% of cases (Easwaran et al. 2012), DNA hypermethylation modules were detected in 42/87 height-associated genes (P = 8.02e~17; Fisher exact test, two tailed). DNA methylation dynamics The extent of DNA methylation of regulatory regions of height-associated genes was experimentally assessed as described by Lv et al. (2011) (www.methdb.de/; 202.97. 205.78/diseasemeth/) (Table S3). All genes where DNA methylation plays a regulatory role are expected to undergo shifts in methylation states. We challenged this prediction in our gene-set. Only five genes (ACAN, ANKS1, FBP2, NACA2, ZBTB38) were found to have no evidence of DNA methylation. The remaining genes (94.3%) were shown to undergo broad changes of DNA methylation levels across experimental conditions (Table S3), consistent with shifting between distinct, DNA meth-ylation-associated regulatory states. Regulation of height-associated genes by DNA methylation A mechanistic relevance of these findings was investigated. Al assessed genes were shown to possess a CpG island in the promoter region that overlaps the TSS (Tables SI, S2), consistent with functional relevance on gene transcription (Esteller 2008). Loss of ADAMTS proteases in nonsmall-cell lung carcinomas is caused by hypermethylation in the ADAMTS genes promoters (Dunn et al. 2004). BMPs are important regulators of cell growth, differentiation, and apoptosis. CpG island methylation in the BMP2 promoter causes loss of BMP-2 protein expression in transformed cells (Wen et al. 2006). Shut-down of the BMP6 gene by promoter methylation was observed in malignant lymphomas (Taniguchi et al. 2008). Cdk5 and Abl enzyme substrate 1 (CABLES1) is a cyclin-dependent kinase binding protein. Loss of nuclear CABLES 1 expression is due to epigenetic modifications of the CABLES1 locus. Full ablation is caused by LOH of the transcriptionally-competent allele (Sakamoto et al. 2008). The expression of the deleted in colorectal cancer (DCC) gene is frequently lost in intestinal cancers. In up to three quarters of the cases the loss of expression of DCC is due to DNA methylation (Park et al. 2008). P. Simeone & S. Albert Loss of the deleted in lymphocytic leukemia 7 (DLEU7) gene is frequently observed in chronic lymphocytic leukemia, due to hypermethylation of the DLEU7 promoter (Hammarsund et al. 2004). Decrease or loss of EGF-containing fibulin-like extracellular matrix protein 1 (EFEMP1) has been shown in human breast cancer. EFEMP1 promoter methylation is a major cause of down-regulation (Sadr-Nabavi et al. 2009). Hedgehog (Hh) signaling activation is frequently mediated by the epigenetic repression of the Hh-interacting protein (HHIP) gene, a negative regulator of Hh signaling. Consistently, HHIP transcription is rescued by deme-thylation (Tada et al. 2008). Hypermethylation of MOS is associated with the development of acute lymphoblastic leukemia (Scholz et al. 2005). CpG methylation inhibits proenkephalin (PENK) gene expression by interfering with the binding of the stimulatory transcription factor AP-2 (Comb and Goodman 1990). Differential methylation of the zinc finger protein gene pleomorphic adenoma gene 1 (PLAGL1) was found in diabetes and cancer and was shown to modulate the expression levels of this gene (Arima and Wake 2006). Patched homolog 1 (PTCH1) is the Hh receptor. Methylation of the PTCH1 promoter inhibits PTCH1 expression in breast cancers (Wolf et al. 2007). The suppressor of cytokine signaling (SOCS) 2 CpG islands are hypermethylated in ovarian cancers. Aberrant methylation of this gene correlates with transcriptional silencing also in hepatocellular carcinomas (Martinez-Chantar et al. 2008). DNA methylation suppresses the expression of thyroid transcription factor-1 (TTF-1) in 60% of undifferentiated thyroid carcinomas (Kondo et al. 2009). c-Myc is a widely expressed TF, that can bind 10-15% of all genes (Orian et al. 2003). c-Myc regulates at least seven height-associated genes (CDK6, COIL, HMGA1, LIN28B, RBBP8, RPS20, TRIM25/EFP), and its binding to genomic loci is dependent on chromatin structure and CpG methylation. Epigenetic defects and hereditary growth anomalies Distinct epigenetic defects have been linked to separate hereditary growth anomalies, providing evidence for a broad regulatory role of DNA methylation on body growth. The Beckwith-Wiedemann syndrome (130650) is caused by deregulation of imprinted genes within the llpl5 chromosomal region, i.e., KIP2, H19 and LIT1, whether alone or as interacting regulatory units (Niemitz 2014 I Vol. 2 I Iss. 6 I e12047 Page 12 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society P. Simeone & S. Alberti et al. 2004). Hypermethylation at the llpl5 telomeric imprinting control region (ICR1), are observed in about 5 to 10% of affected patients (see below for opposite epigenetic changes in Silver-Russel patients). Both H19 and LIT1, which encode untranslated RNAs, and IGF2 are either maternally imprinted genes with growth enhancing activity or paternally imprinted genes with growth suppressing activity. Key features of the syndrome are exom-phalos, macroglossia, and gigantism in the neonate, supporting a direct involvement of epigenetic regulatory mechanisms in body size and body development. Consistently, alterations include visceromegaly, adrenocortical cytomegaly and overgrowth of the external genitalia in both males and females. Affected children reach an average height of 2.5 SD above the mean at or after puberty, and their growth velocity is above the ninetieth percentile until 4-6 years of age. Prader-Willi syndrome (176270) is a complex childhood disorder that affects the nervous and hormonal systems and leads to excessive weight gain. It arises from the loss of activity of several imprinted genes on chromosome 15, that are usually expressed from the paternally inherited copy. Affected children suckle poorly for the first months of their lives; their voracious appetite and consequent obesity develop from weaning onwards. Angelman's syndrome (AS) (105830) arises from opposite imprinting defects to those described in Prader-Willi syndrome. The Angelman syndrome is characterized by mental retardation, movement disorders, altered behavior, and speech abnormalities. Most cases are caused by absence of a maternal contribution to the imprinted region on chromosome 15qll-ql3, whereas the Prader-Willi syndrome results from deletion of the same region in the paternal chromosome. Other patients with an Angelman syndrome carry mutations in the gene encoding methyl-CpG-binding protein-2 (MECP2). MECP2 (300005) is the gene mutated in the Rett syndrome, a neurological (regression of acquired skills, loss of speech, stereotypical movements, seizures, and mental retardation) and developmental disorder that occurs in females. Notably, affected girls are characterized by microcephaly and arrested development between 6 and 18 months of age. Imprinting defects are associated with other developmental disorders, for example Silver-Russell dwarfism. Up to 60% of cases of Silver-Russell syndrome are caused by hypomethylation at the ICR1 on chromosome llpl5, involving the H19 (103280) and IGF2 (147470) genes (Penaherrera et al. 2012). Affected patients demonstrate severe intrauterine growth retardation, poor postnatal growth, craniofacial alterations, and a variety of minor malformations. Epigenetic Regulation of Human Height Epigenetic regulation of body size The findings on Beckwith-Wiedemann syndrome, Prader-Willi syndrome, Angelman's syndrome, Rett syndrome and Silver-Russell syndrome suggest an epigenetic regulation of body size. Consistently, several other imprinting disorders have been shown to affect placental and fetal size (Constancia et al. 2004). The connection between imprinting and placentation is supported by the finding that erasing all genomic imprints results in the outgrowth of extra-embryonic tissues and placenta in animal models. Correspondingly, imprinting defects in the female human germline cause the appearance of hydatidiform moles, that is uncontrolled growth of placental cells. Height-associated regulators of epigenetic heredity Height-associated genes include several key epigenetic regulators of gene expression (Esteller 2008). The activity of transcription factors is modulated by associations with co-repressors, including histone deacetylase 7. BMP2 regulates the transcriptional activity of HDAC7 by inducing its export from the nucleus. The zinc finger and BTB domain containing (ZBTB38) is a methyl-DNA-binding transcriptional repressor gene that binds single methylated CpGs. Chromatin immunoprecipitation indicates that ZBTB38 recognizes the methylated alleles of the H19 and insulin-like growth factor II (IGF-2) genes and represses their transcription (Filion et al. 2006). Methylation-linked inactivation of IGF-2/H19 can be inherited in a parental-specific manner (Riccio et al. 2009). The DOTl-like and NSD1 histone methyltransferases, the two high-mobility group A (HMGA) genes, HMGA1 and HMGA2, and the histone clusters 1 and 2 are all involved in the assembly of chromatin structure. Notably, haploinsufficiency of the histone methyltransferase NSD1 causes the Sotos syndrome, that is characterized by very high stature (Kurotaki et al. 2002) and the gene can be inactivated epigenetically in tumors (Berdasco et al. 2009). PRMT5 contributes to coupling DNA methylation and chromatin organization via the methylation of histone H4R3, which then recruits DNMT3A. DNMT3A mediates additional regulatory loops in chromatin organization, as the HMGA1 gene itself is regulated by DNA methylation. Heterozygous 17qll microdeletions that encompass NF1 and RNF135 cause a subset of neurofibromatosis type 1 (Douglas et al. 2007). Notably, individuals with microdeletions are typically taller than individuals with NF1 mutations, consistently with the removal of a gene that negatively regulates human growth (RNF135 haploinsufficiency) (Douglas et al. 2007). © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 | Vol. 2 | Iss. 6 | e12047 Page 13 Epigenetic Regulation of Human Height P. Simeone & S. Albert Prenatal environmental factors and persistence of DNA methylation changes in the adult Epigenetic patterns in germ cells can persist into adulthood, and many imprinted genes continue to be expressed from a single parental copy. Environmental factors during pregnancy affect epigenetic marks, and in utero conditions have consequences on adult health and disease (Gluckman et al. 2008). This has been proven due to DNA methylation differences after exposure to prenatal famine (Tobi et al. 2009). Animal studies and human data on the imprinted IGF2 locus indicated a link between prenatal nutritional and DNA methylation. Additionally, methylation of the INSIGF, ILIO, LEP, ABCA1, GNASAS and MEG3 genes was persistently modified by exposure to prenatal famine (Tobi et al. 2009). In animal models manipulation of the maternal diet during pregnancy leads to a persistent shift in average DNA methylation levels of specific genes in offspring resulting in permanent changes in coat color or tail shape. Moreover, widespread dynamics of DNA methylation were shown in response to biotic stress (Dowen et al. 2012). Methyl donor supplementation prevents transgen-erational amplification of obesity (Waterland et al. 2008) and maternal methyl supplements were shown to increase DNA methylation in offspring (Waterland et al. 2006). Epigenetic control of height-associated genes networks Data-bank web network meta-analyses allowed us to connect height-associated genes in a function-driven web. These identified p53, c-Myc, estrogen receptor alpha (ERa), HNF4A (Guerra et al. 2012; Trerotola et al. 2013) and SMADs as major hubs (Fig. 1). The Hh pathway and regulatory clusters for programmed cell death/apoptosis were also identified as key control pathways, that are shared between human height and cancer. MYC: Metacore analysis revealed that the c-Myc regulates at least seven height-associated genes (CDK6, COIL, HMGA1, LIN28B, RBBP8, RPS20 and TRIM25/EFP). A SNOW analysis showed that c-Myc is a major hub of the height-associated protein network (19 connections, betweenness 0.366). Thirty-seven height-associated proteins were shown to be associated in the network. The binding of c-Myc to genomic loci is highly dependent on chromatin structure and DNA methylation (Guccione et al. 2006). c-Myc modulates gene expression also by increasing methylation of the 5' mRNA guanine or 'cap', which regulates the translation of individual mRNAs. Methylation of the cap is required for eukaryotic translation initiation fac-tor-4E binding and recruitment onto ribosomes for transla- tion (Cole and Cowling 2008). c-Myc regulates the cell cycle, and plays a major role in cell growth during interphase, by regulating genes required for the production of energy and metabolites. The c-Myc network widely interacts with those driven by other major hubs. c-Myc is repressed by transforming growth factor ß (TGF-ß) through the binding of SMAD3 to the MYC promoter (Frederick et al. 2004). p53 represses c-Myc through the induction of the tumor suppressor miR-145 (Sachdeva et al. 2009). c-Myc amply interacts also with the ER network: almost all of the acutely estrogen-regulated genes with roles in cell growth are c-Myc targets. Notably, estrogen-mediated activation of rRNA and protein synthesis depends on c-Myc (Musgrove et al. 2008). Equally c-Myc dependent is the estrogen-induced suppression of apoptosis caused by growth factor deprivation (Rodrik et al. 2005). TP53: A major hub of height-associated genes is TP53. Remarkably, at least 47 of the 87 height-associated genes were identified as members of p53 signaling networks (Table 1), including 36 via direct protein-protein contacts. p53 was the most important hub (25 connections, betweenness 0.446). Metacore analysis similarly identified p53 as a hub of the height-associated genes network. Input of p53 in the shortest pathway analysis indicated that p53 has a central role with a total of 125 edges pointing either in (37) or out (88) from it. Among its functions, p53 regulates the expression of target genes that modulate chromatin structure and function (Cimoli et al. 2004; Vousden and Lane 2007), cell growth, aging and apoptosis (Ambrogi et al. 2006; Biganzoli et al. 2011). p53 interacts with components of multiple different histone remodeling complexes, including CBP/EP300 (CBP/p300), GCN5, PCAF, and SETD7 modifying histones at the promoters [(Kaneshiro et al. 2007) and reference therein]. We have previously shown that p53 also controls DNA methylation levels, and that this affects genome stability (Alberti et al. 1994; Nasr et al. 2003). p53 is involved in bone remodeling, wound healing and neural tube development (Vousden and Lane 2007). It also performs an anti-teratogenic function, reducing the rate of birth defects, through the induction of apoptosis of aberrant cells (Vousden and Lane 2007). Defects of induction of apoptosis play a role also at later stages of development [(Rossi et al. 2008) and references therein]. Female p53-nu\\ mice can present neural tube-closure defects, due to a failure to induce apoptosis of progenitor cells. This is caused by abnormal regulation of mitochondrial death pathways and results in overproduction of neural tissue. p53 has a role in the regulation of both glycolysis and oxidative phosphorylation. Low levels of p53 determine a switch from oxidative phosphorylation to glycolysis. p53 enhances oxidative phosphorylation by inducing the expression of the cytochrome c oxidase 2 2014 I Vol. 2 I Iss. 6 I e12047 Page 14 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society P. Simeone & S. Alberti subunit 1 and of the ribonucleotide reductase subunit p52R2, which are involved in the maintenance of mitochondrial DNA. p53 promotes cell survival under conditions of low or basal stress. The reduction in nutrient or energy levels causes a failure to stimulate the AKT-mTOR pathway and induces the activation of AMPK both of which lead to the induction of p53. ESR1: ERa is a transcription factor which binds to estrogen response elements upstream of the target genes. META-CORE analysis shows that ERa regulates at least eight height-associated genes (BCAS3, BMP2, BMP6, DCC, GLT25D2, PENK, RBBP8, TRIM25/EPP). ERa is a major hub (16 connections, betweenness 0,376). Estrogen is candidate to be the principal hormone stimulating the pubertal growth spurt in boys as well as girls (Juul 2001). This action is mediated by both ERa and ER/?. Polymorphisms in the ER gene may influence adult height in healthy male and female subjects (Schuit et al. 2004; Dahlgren et al. 2008). Consistently, men with a disruptive mutation in the ER gene have no pubertal growth spurt and continue to grow into adulthood. ERa blockade diminishes the secretion of endogenous growth hormone, the key hormone regulator of linear growth in childhood (Juul 2001). This action is mediated by SOCS-2 (Leung et al. 2003). The ERa network widely interconnects with the p53, Hh and BMP/TGF-/J pathways. p53 regulates ER expression through transcriptional control of the ER promoter (Shirley et al. 2009). Of interest, DNA methyltransferase expression in the human endometrium is down-regulated by progesterone and estrogen (Yamagata et al. 2009). Estrogen receptor-negative human breast cancer cells inactivate the promoter of the ESR1 gene by methylation of its CpG island (Ottaviano et al. 1994). Correspondingly, DNA methylation-mediated gene silencing is the mechanism of ERa inactivation in prostatic epithelial cells (Lau et al. 2000). The SNOW-connection analysis revealed that 33 of the height-associated proteins take part in a unique network of protein-protein contacts (Fig. 1). Relevant hubs in this network are COIL (seven connections, betweenness 0.247) and LYN (seven connections, betweenness 0.146), followed by GRB10 (six connections, betweenness 0.140) and SF3B4 (six connections, betweenness 0.102). Coilin is a key constituent of Cajal bodies, which are responsible for the biogenesis of spliceosomal small nuclear ribonucleoprotein (snRNP) (Whittom et al. 2008). Cajal bodies are the initial nuclear sites for the assembly of macromolecular complexes involved in RNA processing and RNA transcription (Bogolyubov et al. 2009). Lyn is a Src tyrosine-kinase family member. Signaling molecules phosphorylated by Lyn are the PI3-kinase, STAT5 and MAP kinase. GrblO is an adaptor protein that binds to phosphorylated tyrosines, for example, in the insulin receptor (IR) in response to insulin stimulation. GrblO Epigenetic Regulation of Human Height knockout mice show embryo and placenta overgrowth and a larger size at birth [(Wang et al. 2007) and references therein]. SF3b4 is a member of the spliceosome complex U2 snRNP and contributes to the recognition of the intron's branch point. SF3b4 also binds the BMPR-IA serine/threonine kinase receptor and specifically inhibits BMP-mediated SMAD1/5/8 pathway important for osteochondral cell differentiation. Discussion Genomic loci linked to human height have been identified by genome-wide SNP-association analysis (Gudbjarts-son et al. 2008; Lettre et al. 2008; Weedon et al. 2008). Corresponding GWAS in Korean (Cho et al. 2009) and Japanese (Okada et al. 2010) cohorts identified overlapping gene-sets (Weedon and Frayling 2008), showing a conserved association with height in populations of both Caucasian and non-Caucasian ancestry. Intriguingly, although, albeit up to 90% of variation in adult height is explained by genetic factors (Silventoinen et al. 2003; Weedon and Frayling 2008), stature-associated polymorphisms were found to only explain <10% of height variation (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon and Frayling 2008; Weedon et al. 2008; Eichler et al. 2010; Lango Allen et al. 2010). More recent analyses have increased this predictive power (Yang et al. 2010). However, a large fraction of heritable height-associated factors still escapes detection by conventional GWAS (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon and Frayling 2008; Weedon et al. 2008; Eichler et al. 2010; Lango Allen et al. 2010; Yang et al. 2010; Lanktree et al. 2011), suggesting systematic loss of information by SNP-based analyses. Candidate for such non-DNA sequence-linked information is epigenetic heredity. Our findings indicate that most of the height-associated genes contain CpG islands and their transcriptional activity is regulated by DNA methylation. Distinct epigenetic defects have been linked to hereditary growth anomalies, indicating a broad regulatory role of DNA methylation on body growth. Further support is now provided by the association of DNA methylation/regulatory transcription to genes involved in normal body growth (Gudbjartsson et al. 2008; Lettre et al. 2008; Weedon and Frayling 2008; Weedon et al. 2008; Eichler et al. 2010; Lango Allen et al. 2010; Yang et al. 2010; Lanktree et al. 2011). The DNA methylome thus appears to contain an as yet untapped fraction of heritable traits associated with human height, that may be exploited by genome-wide analysis of DNA methylation patterns (Liu et al. 2013). DNA methylation patterns are faithfully propagated through successive cell division cycles (Alberti and Her-zenberg 1988; Esteller 2008). Some of these patterns can © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 | Vol. 2 | Iss. 6 | e12047 Page 15 Epigenetic Regulation of Human Height P. Simeone & S. Alberti be inherited across generations (Kaminsky et al. 2009; Nadeau 2009; Roth et al. 2009; Braunschweig et al. 2012), thus becoming true heritable traits. DNA methylation can regulate gene expression, through the inhibition/activation of gene transcription of methylated/unmethylated genes, respectively (Alberti and Herzenberg 1988; Esteller 2008), and genome stability/recombination (Alberti et al. 1994; Nasr et al. 2003). As DNA methylation patterns are affected by environmental stimuli (El-Osta et al. 2008; Guerrero-Bosagna et al. 2008), this mechanism allows for a dynamic and reversible modulation of the functional content of the genome (Esteller 2008), in the absence of DNA sequence variation. This may help explain as yet unaccounted for multi-generational stature trends, for example, decrease in stature in agricultural populations compared to their paleolithic predecessors (Verginelli et al. 2009), and the recent increase in average height (Cole 2003). Acknowledgments We thank R. Tripaldi, D. D' Ostilio and L. Apicella for support during the course of this work. Disclosures The funding agencies had no role in the research presented in the paper and the researchers were fully independent in pursuing this research. Conflict of Interest None declared. References Alberti, S., and L. A. Herzenberg. 1988. DNA methylation prevents transfection of genes for specific surface antigens. Proc. Natl Acad. Sei. USA 85:8391-8394. Alberti, S., M. Nutini, and L. A. Herzenberg. 1994. DNA methylation prevents the amplification of TROP1, a tumor associated cell surface antigen gene. Proc. Natl Acad. Sei. USA 91:5833-5837. Albig, W., B. Drabent, I. Kunz, M. Kalff-Suske, K. H. Grzeschik, and D. Doenecke. 1993. All known human HI histone genes except the H1(0) gene are clustered on chromosome 6. Genomics 16:649-654. Ambrogi, F., E. Biganzoli, P. Querzoli, S. Ferretti, P. Boracchi, S. Alberti, et al. 2006. Molecular subtyping of breast cancer from traditional tumor marker profiles using parallel clustering methods. Clin. Cancer Res. 12:781-790. Arima, T., and N. Wake. 2006. Establishment of the primary imprint of the HYMAI/PLAGL1 imprint control region during oogenesis. Cytogenet. Genome Res. 113:247-252. Barrett, I. C, and L. R. Cardon. 2006. Evaluating coverage of genome-wide association studies. Nat. Genet. 38:659-662. Berdasco, M., S. Ropero, F. Setien, M. F. Fraga, P. Lapunzina, R. Losson, et al. 2009. Epigenetic inactivation of the Sotos overgrowth syndrome gene histone methyltransferase NSD1 in human neuroblastoma and glioma. Proc. Natl Acad. Sci. USA 106:21830-21835. Biganzoli, E., D. Coradini, F. Ambrogi, I. M. Garibaldi, P. Lisboa, D. Soria, et al. 2011. p53 status identifies two subgroups of triple-negative breast cancers with distinct biological features. Ipn. I. Clin. Oncol. 41:172-179. Bischoff, S. R., S. Tsai, N. Hardison, A. A. Motsinger-Reif, B. A. Freking, D. Nonneman, et al. 2009. Characterization of conserved and nonconserved imprinted genes in swine. Biol. Reprod. 81:906-920. Bogolyubov, D., I. Stepanova, and V. Parfenov. 2009. Universal nuclear domains of somatic and germ cells: some lessons from oocyte interchromatin granule cluster and Cajal body structure and molecular composition. BioEssays 31:400-409. Braunschweig, M., V. lagannathan, A. Gutzwiller, and G. Bee. 2012. Investigations on transgenerational epigenetic response down the male line in f2 pigs. PLoS ONE 7:e30583. Carletti, E., E. Guerra, and S. Alberti. 2006. The forgotten variables of DNA array hybridization. Trends Biotechnol. 24:443-448. Cho, Y. S., M. I. Go, Y. I. Kim, I. Y. Heo, I. H. Oh, H. I. Ban, et al. 2009. A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat. Genet. 41:527-534. Cimoli, G., D. Malacarne, R. Ponassi, M. Valenti, S. Alberti, and S. Parodi. 2004. Meta-analysis of the role of p53 status in isogenic systems tested for sensitivity to cytotoxic antineoplastic drugs. Biochim. Biophys. Acta 1705:103-120. Cole, T. I. 2003. The secular trend in human physical growth: a biological view. Econ. Hum. Biol. 1:161-168. Cole, M. D., and V. H. Cowling. 2008. Transcription-independent functions of MYC: regulation of translation and DNA replication. Nat. Rev. Mol. Cell Biol. 9:810-815. Comb, M., and H. M. Goodman. 1990. CpG methylation inhibits proenkephalin gene expression and binding of the transcription factor AP-2. Nucleic Acids Res. 18:3975-3982. Constancia, M., G. Kelsey, and W. Reik. 2004. Resourceful imprinting. Nature 432:53-57. Dahlgren, A., P. Lundmark, T. Axelsson, L. Lind, and A. C. Syvanen. 2008. Association of the estrogen receptor 1 (ESR1) gene with body height in adult males from two Swedish population cohorts. PLoS ONE 3:el807. Douglas, I., D. Cilliers, K. Coleman, K. Tatton-Brown, K. Barker, B. Bernhard, et al. 2007. Mutations in RNF135, a gene within the NF1 microdeletion region, cause phenotypic abnormalities including overgrowth. Nat. Genet. 39:963-965. 2014 I Vol. 2 I Iss. 6 I e12047 Page 16 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society P. Simeone & S. Alberti Epigenetic Regulation of Human Height Dowen, R. H., M. Pelizzola, R. J. Schmitz, R. Lister, J. M. Dowen, J. R. Neiy, et al. 2012. Widespread dynamic DNA methylation in response to biotic stress. Proc. Natl Acad. Sei. USA 109:E2183-E2191. Dunn, J. R., D. Panutsopulos, M. W. Shaw, J. Heighway, R. Dormer, E. N. Salmo, et al. 2004. METH-2 silencing and promoter hypermethylation in NSCLC. Br. J. Cancer 91:1149-1154. Easwaran, LI., S. E. Johnstone, L. van Neste, J. Ohm, T. Mosbruger, Q. Wang, et al. 2012. A DNA hypermethylation module for the stem/progenitor cell signature of cancer. Genome Res. 22:837-849. Eichler, E. E., J. Flint, G. Gibson, A. Kong, S. M. Leal, J. H. Moore, et al. 2010. Missing heritability and strategies for finding the underlying causes of complex disease. Nat. Rev. Genet. 11:446-450. El-Osta, A., D. Brasacchio, D. Yao, A. Pocai, P. L. Jones, R. G. Roeder, et al. 2008. Transient high glucose causes persistent epigenetic changes and altered gene expression during subsequent normoglycemia. J. Exp. Med. 205:2409-2417. Esteller, M. 2008. Epigenetics in cancer. N. Engl. J. Med. 358:1148-1159. Filion, G. J., S. Zhenilo, S. Salozhin, D. Yamada, E. Prokhortchouk, and P. A. Defossez. 2006. A family of human zinc finger proteins that bind methylated DNA and repress transcription. Mol. Cell. Biol. 26:169-181. Frederick, J. P., N. T. Liberati, D. S. Waddell, Y. Shi, and X. F. Wang. 2004. Transforming growth factor beta-mediated transcriptional repression of c-myc is dependent on direct binding of Smad3 to a novel repressive Smad binding element. Mol. Cell. Biol. 24:2546-2559. Gluckman, P. D., M. A. Hanson, C. Cooper, and K. L. Thornburg. 2008. Effect of in utero and early-life conditions on adult health and disease. N. Engl. J. Med. 359:61-73. Greer, E. L., T. J. Maures, D. Ucar, A. G. Hauswirth, E. Mancini, J. P. Lim, et al. 2011. Transgenerational epigenetic inheritance of longevity in Caenorhabditis elegans. Nature 479:365-371. Guccione, E., F. Martinato, G. Finocchiaro, L. Luzi, L. Tizzoni, V. Dali' Olio, et al. 2006. Myc-binding-site recognition in the human genome is determined by chromatin context. Nat. Cell Biol. 8:764-770. Gudbjartsson, D. F., G. B. Walters, G. Thorleifsson, H. Stefansson, B. V. Halldorsson, P. Zusmanovich, et al. 2008. Many sequence variants affecting diversity of adult human height. Nat. Genet. 40:609-615. Guerra, E., R. Lattanzio, R. la Sorda, F. Dini, G. M. Tiboni, M. Piantelli, et al. 2012. mTropl/Epcam knockout mice develop congenital tufting enteropathy through dysregulation of intestinal E-cadherin//ä-catenin. PLoS ONE 7:e49302. Guerrero-Bosagna, C. M., P. Sabat, F. S. Valdovinos, L. E. Valladares, and S. J. Clark. 2008. Epigenetic and phenotypic changes result from a continuous pre and post natal dietary exposure to phytoestrogens in an experimental population of mice. BMC Physiol. 8:17. Hammarsund, M., M. M. Corcoran, W. Wilson, C. Zhu, S. Einhorn, O. Sangfelt, et al. 2004. Characterization of a novel B-CLL candidate gene-DLEU7-located in the 13ql4 tumor suppressor locus. FEBS Lett. 556:75-80. Hollingsworth, J. W., S. Maruoka, K. Boon, S. Garantziotis, Z. Li, J. Tomfohr, et al. 2008. In utero supplementation with methyl donors enhances allergic airway disease in mice. J. Clin. Invest. 118:3462-3469. Iglesias-Platas, I., A. Martin-Trujillo, D. Cirillo, F. Court, A. Guillaumet-Adkins, C. Camprubi, et al. 2012. Characterization of novel paternal ncRNAs at the Plagll locus, including Hymai, predicted to interact with regulators of active chromatin. PLoS ONE 7:e38907. Illingworth, R. S., and A. P. Bird. 2009. CpG islands - 'A rough guide'. FEBS Lett. 583:1713-1720. Jiang, C, L. Han, B. Su, W. H. Li, and Z. Zhao. 2007. Features and trend of loss of promoter-associated CpG islands in the human and mouse genomes. Mol. Biol. Evol. 24:1991-2000. Jin, B., J. Ernst, R. L. Tiedemann, H. Xu, S. Sureshchandra, M. Kellis, et al. 2012. Linking DNA methyltransferases to epigenetic marks and nucleosome structure genome-wide in human tumor cells. Cell Rep. 2:1411-1424. Jones, B., H. Su, A. Bhat, H. Lei, J. Bajko, S. Hevi, et al. 2008. The histone H3K79 methyltransferase DotlL is essential for mammalian development and heterochromatin structure. PLoS Genet. 4:el000190. Juul, A. 2001. The effects of oestrogens on linear bone growth. Hum. Reprod. Update 7:303-313. Kaminsky, Z. A., T. Tang, S. C. Wang, C. Ptak, G. H. Oh, A. H. Wong, et al. 2009. DNA methylation profiles in monozygotic and dizygotic twins. Nat. Genet. 41:240-245. Kaneshiro, K, S. Tsutsumi, S. Tsuji, K. Shirahige, and H. Aburatani. 2007. An integrated map of p53-binding sites and histone modification in the human ENCODE regions. Genomics 89:178-188. Kondo, T., T. Nakazawa, D. Ma, D. Niu, K. Mochizuki, T. Kawasaki, et al. 2009. Epigenetic silencing of TTF-1/ NKX2-1 through DNA hypermethylation and histone H3 modulation in thyroid carcinomas. Lab. Invest. 89: 791-799. Koudritsky, M., and E. Domany. 2008. Positional distribution of human transcription factor binding sites. Nucleic Acids Res. 36:6795-6805. Kurotaki, N., K. Imaizumi, N. Harada, M. Masuno, T. Kondoh, T. Nagai, et al. 2002. Haploinsufficiency of NSD1 causes Sotos syndrome. Nat. Genet. 30:365-366. Lander, E. S., L. M. Linton, B. Birren, C. Nusbaum, M. C. Zody, J. Baldwin, et al. 2001. Initial sequencing and analysis of the human genome. Nature 409:860-921. © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 I Vol. 2 I Iss. 6 I e12047 Page 17 Epigenetic Regulation of Human Height P. Simeone & S. Alberti Lango Allen, H., K. Estrada, G. Lettre, S. I. Berndt, M. N. Weedon, F. Rivadeneira, et al. 2010. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature 467:832-838. Lanktree, M. B., Y. Guo, M. Murtaza, J. T. Glessner, S. D. Bailey, N. C. Onland-Moret, et al. 2011. Meta-analysis of dense genecentric association studies reveals common and uncommon variants associated with height. Am. J. Hum. Genet. 88:6-18. Larsen, F., G. Gundersen, R. Lopez, and H. Prydz. 1992. CpG islands as gene markers in the human genome. Genomics 13:1095-1107. Lau, K. M., M. Laspina, J. Long, and S. M. Ho. 2000. Expression of estrogen receptor (ER)-alpha and ER-beta in normal and malignant prostatic epithelial cells: regulation by methylation and involvement in growth regulation. Cancer Res. 60:3175-3182. Lettre, G., J. L. Butler, K. G. Ardlie, and J. N. Hirschhorn. 2007. Common genetic variation in eight genes of the GH/ IGF1 axis does not contribute to adult height variation. Hum. Genet. 122:129-139. Lettre, G., A. U. Jackson, C. Gieger, F. R. Schumacher, S. I. Berndt, S. Sanna, et al. 2008. Identification of ten loci associated with height highlights new biological pathways in human growth. Nat. Genet. 40:584-591. Leung, K. C, N. Doyle, M. Ballesteros, K. Sjogren, C. K. Watts, T. H. Low, et al. 2003. Estrogen inhibits GH signaling by suppressing GH-induced JAK2 phosphorylation, an effect mediated by SOCS-2. Proc. Natl Acad. Sei. USA 100:1016-1021. Li, Y., G. Meng, and Q. N. Guo. 2008. Changes in genomic imprinting and gene expression associated with transformation in a model of human osteosarcoma. Exp. Mol. Pathol. 84:234-239. Li, L., J. Tang, B. Zhang, W. Yang, M. Liugao, R. Wang, et al. 2014. Epigenetic modification of MiR-429 promotes liver tumour-initiating cell properties by targeting Rb binding protein 4. Gut doi: 10.1136/gutjnl-2013-305715. [Epub ahead of print]. Liu, Y, M. J. Aryee, L. Padyukov, M. D. Fallin, E. Hesselberg, A. Runarsson, et al. 2013. Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat. Biotechnol. 31:142-147. Lv, J., H. Liu, J. Su, X. Wu, B. Li, X. Xiao, et al. 2011. DiseaseMeth: a human disease methylation database. Nucleic Acids Res. 40:D1030-D1035. Mackay, T. F., E. A. Stone, and J. F. Ayroles. 2009. The genetics of quantitative traits: challenges and prospects. Nat. Rev. Genet. 10:565-577. Martinez-Chantar, M. L., M. Vazquez-Chantada, U. Ariz, N. Martinez, M. Varela, Z. Luka, et al. 2008. Loss of the glycine N-methyl transferase gene leads to steatosis and hepatocellular carcinoma in mice. Hepatology 47: 1191-1199. Miura, K., M. Agetsuma, H. Kitano, A. Yoshimura, M. Matsuoka, S. E. Jacobsen, et al. 2009. A metastable DWARF1 epigenetic mutant affecting plant stature in rice. Proc. Natl Acad. Sei. USA 106:11218-11223. Musgrove, E. A., C. M. Sergio, S. Loi, C. K. Inman, L. R. Anderson, M. C. Alles, et al. 2008. Identification of functional networks of estrogen- and c-Myc-responsive genes and their relationship to response to tamoxifen therapy in breast cancer. PLoS ONE 3:e2987. Nadeau, J. H. 2009. Transgenerational genetic effects on phenotypic variation and disease risk. Hum. Mol. Genet. 18: R202-R210. Nasr, A. F., M. Nutini, B. Palombo, E. Guerra, and S. Alberti. 2003. Mutations of TP53 induce loss of DNA methylation and amplification of the TROP1 gene. Oncogene 22:1668-1677. Nazor, K. L., G. Altun, C. Lynch, H. Tran, J. V. Harness, I. Slavin, et al. 2012. Recurrent variations in DNA methylation in human pluripotent stem cells and their differentiated derivatives. Cell Stem Cell 10: 620-634. Nguyen, C, G. Liang, T. T. Nguyen, D. Tsao-Wei, S. Groshen, M. Lubbert, et al. 2001. Susceptibility of nonpromoter CpG islands to de novo methylation in normal and neoplastic cells. J. Natl Cancer Inst. 93:1465-1472. Niemitz, E. L., M. R. Debaun, J. Fallon, K. Murakami, H. Kugoh, M. Oshimura, et al. 2004. Microdeletion of LIT1 in familial Beckwith-Wiedemann syndrome. Am. J. Hum. Genet. 75:844-849. Okada, Y., Y. Kamatani, A. Takahashi, K. Matsuda, N. Hosono, H. Ohmiya, et al. 2010. A genome-wide association study in 19 633 Japanese subjects identified LHX3-QSOX2 and IGF1 as adult height loci. Hum. Mol. Genet. 19:2303-2312. Orian, A., B. van Steensel, J. Delrow, H. J. Bussemaker, L. Li, T. Sawado, et al. 2003. Genomic binding by the Drosophila Myc, Max, Mad/Mnt transcription factor network. Genes Dev. 17:1101-1114. Ottaviano, Y. L., J. P. Issa, F. F. Pari, H. S. Smith, S. B. Baylin, and N. E. Davidson. 1994. Methylation of the estrogen receptor gene CpG island marks loss of estrogen receptor expression in human breast cancer cells. Cancer Res. 54:2552-2555. Park, H. L., M. S. Kim, K. Yamashita, W. Westra, A. L. Carvalho, J. Lee, et al. 2008. DCC promoter hypermethylation in esophageal squamous cell carcinoma. Int. J. Cancer 122:2498-2502. Penaherrera, M. S., S. Weindler, M. I. van Allen, S. L. Yong, D. L. Metzger, B. McGillivray, et al. 2012. Methylation profiling in individuals with Russell-Silver syndrome. Am. J. Med. Genet. A. 152A:347-355. Reeves, R., and L. Beckerbauer. 2001. HMGI/Y proteins: flexible regulators of transcription and chromatin structure. Biochim. Biophys. Acta 1519:13-29. 2014 I Vol. 2 I Iss. 6 I e12047 Page 18 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society P. Simeone & S. Alberti Epigenetic Regulation of Human Height Riccio, A., A. Sparago, G. Verde, A. de Crescenzo, V. Citro, M. V. Cubellis, et al. 2009. Inherited and Sporadic Epimutations at the IGF2-H19 locus in Beckwith-Wiedemann syndrome and Wilms' tumor. Endocr. Dev. 14:1-9. Rodrik, V., Y. Zheng, F. Harrow, Y. Chen, and D. A. Foster. 2005. Survival signals generated by estrogen and phospholipase D in MCF-7 breast cancer cells are dependent on Myc. Mol. Cell. Biol. 25:7917-7925. Rossi, D. J., C. H. lamieson, and I. L. Weissman. 2008. Stems cells and the pathways to aging and cancer. Cell 132:681-696. Roth, T. L., F. D. Lubin, A. J. Funk, and J. D. Sweatt. 2009. Lasting epigenetic influence of early-life adversity on the BDNF gene. Biol. Psychiatry 65:760-769. Sachdeva, M., S. Zhu, F. Wu, H. Wu, V. Walia, S. Kumar, et al. 2009. p53 represses c-Myc through induction of the tumor suppressor miR-145. Proc. Natl Acad. Sci. USA 106:3207-3212. Sadr-Nabavi, A., J. Ramser, J. Volkmann, J. Naehrig, F. Wiesmann, B. Betz, et al. 2009. Decreased expression of angiogenesis antagonist EFEMP1 in sporadic breast cancer is caused by aberrant promoter methylation and points to an impact of EFEMP1 as molecular biomarker. Int. J. Cancer 124:1727-1735. Sakamoto, H., A. M. Friel, A. W. Wood, L. Guo, A. Ilic, M. V. Seiden, et al. 2008. Mechanisms of Cables 1 gene inactivation in human ovarian cancer development. Cancer Biol. Ther. 7:180-188. Sanna, S., A. U. lackson, R. Nagaraja, C. J. Wilier, W. M. Chen, L. L. Bonnycastle, et al. 2008. Common variants in the GDF5-UQCC region are associated with variation in human height. Nat. Genet. 40:198-203. Saxonov, S., P. Berg, and D. L. Brutlag. 2006. A genome-wide analysis of CpG dinucleotides in the human genome distinguishes two distinct classes of promoters. Proc. Natl Acad. Sci. USA 103:1412-1417. Scholz, C, I. Nimmrich, M. Burger, E. Becker, B. Dorken, W. D. Ludwig, et al. 2005. Distinction of acute lymphoblastic leukemia from acute myeloid leukemia through microarray-based DNA methylation analysis. Ann. Hematol. 84:236-244. Schuit, S. C, J. B. van Meurs, A. P. Bergink, M. van der Klift, Y. Fang, G. Leusink, et al. 2004. Height in pre- and postmenopausal women is influenced by estrogen receptor alpha gene polymorphisms. J. Clin. Endocrinol. Metab. 89:303-309. Sgarra, R., F. Diana, C. Bellarosa, V. Dekleva, A. Rustighi, M. Toller, et al. 2003. During apoptosis of tumor cells HMGAla protein undergoes methylation: identification of the modification site by mass spectrometry. Biochemistry 42:3575-3585. Shirley, S. H., J. E. Rundhaug, J. Tian, N. Cullinan-Ammann, I. Lambertz, C. J. Conti, et al. 2009. Transcriptional regulation of estrogen receptor-alpha by p53 in human breast cancer cells. Cancer Res. 69:3405-3414. Silventoinen, K., J. Kaprio, E. Lahelma, and M. Koskenvuo. 2000. Relative effect of genetic and environmental factors on body height: differences across birth cohorts among Finnish men and women. Am. J. Public Health 90:627-630. Silventoinen, K., S. Sammalisto, M. Perola, D. I. Boomsma, B. K. Cornes, C. Davis, et al. 2003. Heritability of adult body height: a comparative study of twin cohorts in eight countries. Twin Res. 6:399-408. Sutherland, K. D., G. J. Lindeman, D. Y. Choong, S. Wittlin, L. Brentzell, W. Phillips, et al. 2004. Differential hypermethylation of SOCS genes in ovarian and breast carcinomas. Oncogene 23:7726-7733. Tada, M., F. Kanai, Y. Tanaka, K. Tateishi, M. Ohta, Y. Asaoka, et al. 2008. Down-regulation of hedgehog-interacting protein through genetic and epigenetic alterations in human hepatocellular carcinoma. Clin. Cancer Res. 14:3768-3776. Takai, D., and P. A. lones. 2002. Comprehensive analysis of CpG islands in human chromosomes 21 and 22. Proc. Natl Acad. Sci. USA 99:3740-3745. Taniguchi, A., Y. Němoto, A. Yokoyama, N. Kotani, S. Imai, T. Shuin, et al. 2008. Promoter methylation of the bone morphogenetic protein-6 gene in association with adult T-cell leukemia. Int. J. Cancer 123:1824-1831. Tobi, E. W., L. H. Lumey, R. P. Talens, D. Kremer, H. Putter, A. D. Stein, et al. 2009. DNA Methylation differences after exposure to prenatal famine are common and timing- and sex-specific. Hum. Mol. Genet. 18: 4046-4053. Trerotola, M., P. Cantanelli, E. Guerra, R. Tripaldi, A. L. Aloisi, V. Bonasera, et al. 2013. Up-regulation of Trop-2 quantitatively stimulates human cancer growth. Oncogene 32:222-233. Verginelli, F., F. Aru, P. Battista, and R. Mariáni-Costantini. 2009. Nutrigenetics in the light of human evolution. J. Nutrigenet. Nutrigenomics 2:91-102. Vousden, K. H., and D. P. Lane. 2007. p53 in health and disease. Nat. Rev. Mol. Cell Biol. 8:275-283. Wang, L., B. Balas, C. Y. Christ-Roberts, R. Y. Kim, F. J. Ramos, C. K. Kikani, et al. 2007. Peripheral disruption of the GrblO gene enhances insulin signaling and sensitivity in vivo. Mol. Cell. Biol. 27:6497-6505. Waterland, R. A., D. C. Dolinoy, J. R. Lin, C. A. Smith, X. Shi, and K. G. Tahiliani. 2006. Maternal methyl supplements increase offspring DNA methylation at Axin Fused. Genesis 44:401-406. Waterland, R. A., M. Travisano, K. G. Tahiliani, M. T. Rached, and S. Mirza. 2008. Methyl donor supplementation prevents transgenerational amplification of obesity. Int. J. Obes. (Lond). 32:1373-1379. © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society 2014 I Vol. 2 I Iss. 6 I e12047 Page 19 Epigenetic Regulation of Human Height P. Simeone & S. Alberti Weber, M., I. Hellmann, M. B. Stadler, L. Ramos, S. Paabo, M. Rebhan, et al. 2007. Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat. Genet. 39:457-466. Weedon, M. N., and T. M. Frayling. 2008. Reaching new heights: insights into the genetics of human stature. Trends Genet. 24:595-603. Weedon, M. N., G. Lettre, R. M. Freathy, C. M. Lindgren, B. F. Voight, J. R. Perry, et al. 2007. A common variant of HMGA2 is associated with adult and childhood height in the general population. Nat. Genet. 39:1245-1250. Weedon, M. N., H. Lango, C. M. Lindgren, C. Wallace, D. M. Evans, M. Mangino, et al. 2008. Genome-wide association analysis identifies 20 loci that influence adult height. Nat. Genet. 40:575-583. Wen, X. Z., Y. Akiyama, S. B. Baylin, and Y. Yuasa. 2006. Frequent epigenetic silencing of the bone morphogenetic protein 2 gene through methylation in gastric carcinomas. Oncogene 25:2666-2673. Whittom, A. A., H. Xu, and M. D. Hebert. 2008. Coilin levels and modifications influence artificial reporter splicing. Cell. Mol. Life Sei. 65:1256-1271. Wolf, I., S. Bose, J. C. Desmond, B. T. Lin, E. A. Williamson, B. Y. Karlan, et al. 2007. Unmasking of epigenetically silenced genes reveals DNA promoter methylation and reduced expression of PTCH in breast cancer. Breast Cancer Res. Treat. 105:139-155. Yamagata, Y., H. Asada, I. Tamura, L. Lee, R. Maekawa, K. Taniguchi, et al. 2009. DNA methyltransferase expression in the human endometrium: down-regulation by progesterone and estrogen. Hum. Reprod. 24:1126-1132. Yang, ]., B. Benyamin, B. P. McEvoy, S. Gordon, A. K. Henders, D. R. Nyholt, et al. 2010. Common SNPs explain a large proportion of the heritability for human height. Nat. Genet. 42:565-569. Zhu, J., F. He, S. Hu, and J. Yu. 2008. On the nature of human housekeeping genes. Trends Genet. 24:481-484. Supporting Information Additional Supporting Information may be found in the online version of this article: Table SI. CpG islands in human height-associated genes transcribed regions. Table S2. CpG islands in human height-associated genes - UCSC Genome Browser algorithm. Table S3. Structural an functional features of DNA methylation in height-associated genes. 2014 I Vol. 2 I Iss. 6 I e12047 Page 20 © 2014 The Authors. Physiological Reports published by Wiley Periodicals, Inc. on behalf of the American Physiological Society and The Physiological Society