Skeletal muscle tissue demonstrates global hypermethylation with age. However, methylome changes across the time-course of differentiation in aged human muscle derived cells, and larger coverage arrays in aged muscle tissue have not been undertaken. Using 850K DNA methylation arrays we compared the methylomes of young (27 ± 4.4 years) and aged (83 ± 4 years) human skeletal muscle and that of young/aged heterogenous muscle-derived human primary cells (HDMCs) over several time points of differentiation (0, 72 h, 7, 10 days). Aged muscle tissue was hypermethylated compared with young tissue, enriched for; pathways-in-cancer (including; focal adhesion, MAPK signaling, PI3K-Akt-mTOR signaling, p53 signaling, Jak-STAT signaling, TGF-beta and notch signaling), rap1-signaling, axon-guidance and hippo-signalling. Aged cells also demonstrated a hypermethylated profile in pathways; axon-guidance, adherens-junction and calcium-signaling, particularly at later timepoints of myotube formation, corresponding with reduced morphological differentiation and reductions in MyoD/Myogenin gene expression compared with young cells. While young cells showed little alterations in DNA methylation during differentiation, aged cells demonstrated extensive and significantly altered DNA methylation, particularly at 7 days of differentiation and most notably in focal adhesion and PI3K-AKT signalling pathways. While the methylomes were vastly different between muscle tissue and HDMCs, we identified a small number of CpG sites showing a hypermethylated state with age, in both muscle tissue and cells on genes KIF15, DYRK2, FHL2, MRPS33, ABCA17P. Most notably, differential methylation analysis of chromosomal regions identified three locations containing enrichment of 6–8 CpGs in the HOX family of genes altered with age. With HOXD10, HOXD9, HOXD8, HOXA3, HOXC9, HOXB1, HOXB3, HOXC-AS2 and HOXC10 all hypermethylated in aged tissue. In aged cells the same HOX genes (and additionally HOXC-AS3) displayed the most variable methylation at 7 days of differentiation versus young cells, with HOXD8, HOXC9, HOXB1 and HOXC-AS3 hypermethylated and HOXC10 and HOXC-AS2 hypomethylated. We also determined that there was an inverse relationship between DNA methylation and gene expression for HOXB1, HOXA3 and HOXC-AS3. Finally, increased physical activity in young adults was associated with oppositely regulating HOXB1 and HOXA3 methylation compared with age. Overall, we demonstrate that a considerable number of HOX genes are differentially epigenetically regulated in aged human skeletal muscle and HDMCs and increased physical activity may help prevent age-related epigenetic changes in these HOX genes.

Whole-exome sequencing is an attractive alternative to microarray analysis because of the low cost and potential ability to detect copy number variations (CNV) of various sizes (from 1–2 exons to several Mb). Previous comparison of the most popular CNV calling tools showed a high portion of false-positive calls. Moreover, due to a lack of a gold standard CNV set, the results are limited and incomparable. Here, we aimed to perform a comprehensive analysis of tools capable of germline CNV calling available at the moment using a single CNV standard and reference sample set. Compiling variants from previous studies with Bayesian estimation approach, we constructed an internal standard for NA12878 sample (pilot National Institute of Standards and Technology Reference Material) including 110,050 CNV or non-CNV exons. The standard was used to evaluate the performance of 16 germline CNV calling tools on the NA12878 sample and 10 correlated exomes as a reference set with respect to length distribution, concordance, and efficiency. Each algorithm had a certain range of detected lengths and showed low concordance with other tools. Most tools are focused on detection of a limited number of CNVs one to seven exons long with a false-positive rate below 50%. EXCAVATOR2, exomeCopy, and FishingCNV focused on detection of a wide range of variations but showed low precision. Upon unified comparison, the tools were not equivalent. The analysis performed allows choosing algorithms or ensembles of algorithms most suitable for a specific goal, e.g. population studies or medical genetics.

Nuclear noncoding RNAs (ncRNAs) are key regulators of gene expression and chromatin organization. The progress in studying nuclear ncRNAs depends on the ability to identify the genome-wide spectrum of contacts of ncRNAs with chromatin. To address this question, a panel of RNA–DNA proximity ligation techniques has been developed. However, neither of these techniques examines proteins involved in RNA–chromatin interactions. Here, we introduce RedChIP, a technique combining RNA–DNA proximity ligation and chromatin immunoprecipitation for identifying RNA–chromatin interactions mediated by a particular protein. Using antibodies against architectural protein CTCF and the EZH2 subunit of the Polycomb repressive complex 2, we identify a spectrum of cis- and trans-acting ncRNAs enriched at Polycomb- and CTCF-binding sites in human cells, which may be involved in Polycomb-mediated gene repression and CTCF-dependent chromatin looping. By providing a protein-centric view of RNA–DNA interactions, RedChIP represents an important tool for studies of nuclear ncRNAs.