Certain human traits such as neurodevelopmental disorders (NDs) and congenital anomalies (CAs) are believed to be primarily genetic in origin. However, even after whole-genome sequencing (WGS), a substantial fraction of such disorders remain unexplained. We hypothesize that some cases of ND–CA are caused by aberrant DNA methylation leading to dysregulated genome function. Comparing DNA methylation profiles from 489 individuals with ND–CAs against 1534 controls, we identify epivariations as a frequent occurrence in the human genome. De novo epivariations are significantly enriched in cases, while RNAseq analysis shows that epivariations often have an impact on gene expression comparable to loss-of-function mutations. Additionally, we detect and replicate an enrichment of rare sequence mutations overlapping CTCF binding sites close to epivariations, providing a rationale for interpreting non-coding variation. We propose that epivariations contribute to the pathogenesis of some patients with unexplained ND–CAs, and as such likely have diagnostic relevance.

Despite representing an important source of genetic variation, tandem repeats (TRs) remain poorly studied due to technical difficulties. We hypothesized that TRs can operate as expression (eQTLs) and methylation (mQTLs) quantitative trait loci. To test this we analyzed the effect of variation at 4849 promoter-associated TRs, genotyped in 120 individuals, on neighboring gene expression and DNA methylation. Polymorphic promoter TRs were associated with increased variance in local gene expression and DNA methylation, suggesting functional consequences related to TR variation. We identified >100 TRs associated with expression/methylation levels of adjacent genes. These potential eQTL/mQTL TRs were enriched for overlaps with transcription factor binding and DNaseI hypersensitivity sites, providing a rationale for their effects. Moreover, we showed that most TR variants are poorly tagged by nearby single nucleotide polymorphisms (SNPs) markers, indicating that many functional TR variants are not effectively assayed by SNP-based approaches. Our study assigns biological significance to TR variations in the human genome, and suggests that a significant fraction of TR variations exert functional effects via alterations of local gene expression or epigenetics. We conclude that targeted studies that focus on genotyping TR variants are required to fully ascertain functional variation in the genome.

Genomic imprinting is a mechanism in which gene expression varies depending on parental origin. Imprinting occurs through differential epigenetic marks on the two parental alleles, with most imprinted loci marked by the presence of differentially methylated regions (DMRs). To identify sites of parental epigenetic bias, here we have profiled DNA methylation patterns in a cohort of 57 individuals with uniparental disomy (UPD) for 19 different chromosomes, defining imprinted DMRs as sites where the maternal and paternal methylation levels diverge significantly from the biparental mean. Using this approach we identified 77 DMRs, including nearly all those described in previous studies, in addition to 34 DMRs not previously reported. These include a DMR at TUBGCP5 within the recurrent 15q11.2 microdeletion region, suggesting potential parent-of-origin effects associated with this genomic disorder. We also observed a modest parental bias in DNA methylation levels at every CpG analyzed across ∼1.9 Mb of the 15q11-q13 Prader-Willi/Angelman syndrome region, demonstrating that the influence of imprinting is not limited to individual regulatory elements such as CpG islands, but can extend across entire chromosomal domains. Using RNA-seq data, we detected signatures consistent with imprinted expression associated with nine novel DMRs. Finally, using a population sample of 4,004 blood methylomes, we define patterns of epigenetic variation at DMRs, identifying rare individuals with global gain or loss of methylation across multiple imprinted loci. Our data provide a detailed map of parental epigenetic bias in the human genome, providing insights into potential parent-of-origin effects.

The extent to which the DNA relaxation activities of eukaryotic topoisomerases (topo I and topo II) are redundant during gene transcription is unclear. Although both enzymes can often substitute for each other in vivo, studies in vitro had revealed that the DNA cross-inversion mechanism of topo II relaxes chromatin more efficiently than the DNA strand-rotation mechanism of topo I. Here, we show that the inactivation of topo II in budding yeast produces an abrupt decrease of virtually all polyA+ RNA transcripts of length above ∼3 kb, irrespective of their function. This reduction is not related to transcription initiation but to the stall of RNA polymerase II (Pol II) during elongation. This reduction does not occur in topo I mutants; and it is not avoided by overproducing yeast topo I or bacterial topo I, which relaxes (−) DNA supercoils. It is rescued by catalytically active topo II or a GyrBA enzyme, which relaxes (+) DNA supercoils. These findings demonstrate that DNA relaxation activities of topo I and topo II are not interchangeable in vivo. Apparently, only topo II relaxes efficiently the (+) DNA supercoils that stall the advancement of Pol II in long genes. A mechanistic model is proposed.

While population studies have resulted in detailed maps of genetic variation in humans, to date there are few robust maps of epigenetic variation. We identified sites containing clusters of CpGs with high inter-individual epigenetic variation, termed Variably Methylated Regions (VMRs) in five purified cell types. We observed that VMRs occur preferentially at enhancers and 3’ UTRs. While the majority of VMRs have high heritability, a subset of VMRs within the genome show highly correlated variation in trans, forming co-regulated networks that have low heritability, differ between cell types and are enriched for specific transcription factor binding sites and biological pathways of functional relevance to each tissue. For example, in T cells we defined a network of 95 co-regulated VMRs enriched for genes with roles in T-cell activation; in fibroblasts a network of 34 co-regulated VMRs comprising all four HOX gene clusters enriched for control of tissue growth; and in neurons a network of 18 VMRs enriched for roles in synaptic signaling. By culturing genetically-identical fibroblasts under varying environmental conditions, we experimentally demonstrated that some VMR networks are responsive to the environment, with methylation levels at these loci changing in a coordinated fashion in trans dependent on cellular growth. Intriguingly these environmentally-responsive VMRs showed a strong enrichment for imprinted loci (p<10−80), suggesting that these are particularly sensitive to environmental conditions. Our study provides a detailed map of common epigenetic variation in the human genome, showing that both genetic and environmental causes underlie this variation.

How DNA helical tension is constrained along the linear chromosomes of eukaryotic cells is poorly understood. In this study, we induced the accumulation of DNA ( þ ) helical tension in Saccharomyces cerevisiae cells and examined how DNA transcription was affected along yeast chromosomes. The results revealed that, whereas the overwinding of DNA produced a general impairment of transcription initiation, genes situated at o100 kb from the chromosomal ends gradually escaped from the transcription stall. This novel positional effect seemed to be a simple function of the gene distance to the telomere: It occurred evenly in all 32 chromosome extremities and was independent of the atypical structure and transcription activity of subtelomeric chromatin. These results suggest that DNA helical tension dissipates at chromosomal ends and, therefore, provides a functional indication that yeast chromosome extremities are topologically open. The gradual escape from the transcription stall along the chromosomal flanks also indicates that friction restrictions to DNA twist diffusion, rather than tight topological boundaries, might suffice to confine DNA helical tension along eukaryotic chromatin.

No abstract

Transfer RNA (tRNA) activity is tightly regulated to provide a physiological protein translation, and tRNA chemical modifications control its function in a complex with ribosomes and messenger RNAs (mRNAs). In this regard, the correct hypermodification of position G37 of phenylalanine-tRNA, adjacent to the anticodon, is critical to prevent ribosome frameshifting events. Here we report that the tRNA-yW Synthesizing Protein 2 (TYW2) undergoes promoter hypermethylation-associated transcriptional silencing in human cancer, particularly in colorectal tumors. The epigenetic loss of TYW2 induces guanosine hypomodification in phenylalanine-tRNA, an increase in −1 ribosome frameshift events, and down-regulation of transcripts by mRNA decay, such as of the key cancer gene ROBO1. Importantly, TYW2 epigenetic inactivation is linked to poor overall survival in patients with early-stage colorectal cancer, a finding that could be related to the observed acquisition of enhanced migration properties and epithelial-to-mesenchymal features in the colon cancer cells that harbor TYW2 DNA methylation-associated loss. These findings provide an illustrative example of how epigenetic changes can modify the epitranscriptome and further support a role for tRNA modifications in cancer biology.