DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Recent developments in whole genome bisulfite sequencing (WGBS) technology have enabled genome-wide measurements of DNA methylation at single base pair resolution. Many experiments have been conducted to compare DNA methylation profiles under different biological contexts, with the goal of identifying differentially methylated regions (DMRs). Due to the high cost of WGBS experiments, many studies are still conducted without biological replicates. Methods and tools available for analyzing such data are very limited. We develop a statistical method, DSS-single, for detecting DMRs from WGBS data without replicates. We characterize the count data using a rigorous model that accounts for the spatial correlation of methylation levels, sequence depth and biological variation. We demonstrate that using information from neighboring CG sites, biological variation can be estimated accurately even without replicates. DMR detection is then carried out via a Wald test procedure. Simulations demonstrate that DSS-single has greater sensitivity and accuracy than existing methods, and an analysis of H1 versus IMR90 cell lines suggests that it also yields the most biologically meaningful results. DSS-single is implemented in the Bioconductor package DSS.

Chemical modifications on DNA molecules, such as 5-methylcytosine and 5-hydroxymethylcytosine, play important roles in the mammalian brain. A novel DNA adenine modification, N(6)-methyladenine (6mA), has recently been found in mammalian cells. However, the presence and function(s) of 6mA in the mammalian brain remain unclear. Here we demonstrate 6mA dynamics in the mouse brain in response to environmental stress. We find that overall 6mA levels are significantly elevated upon stress. Genome-wide 6mA and transcriptome profiling reveal an inverse association between 6mA dynamic changes and a set of upregulated neuronal genes or downregulated LINE transposon expression. Genes bearing stress-induced 6mA changes significantly overlap with loci associated with neuropsychiatric disorders. These results suggest an epigenetic role for 6mA in the mammalian brain as well as its potential involvement in neuropsychiatric disorders.

Cytoplasmic filamentous rods and rings (RR) structures were identified using human autoantibodies as probes. In the present study, the formation of these conserved structures in mammalian cells and functions linked to these structures were examined. Distinct cytoplasmic rods (∼3–10 µm in length) and rings (∼2–5 µm in diameter) in HEp-2 cells were initially observed in immunofluorescence using human autoantibodies. Co-localization studies revealed that, although RR had filament-like features, they were not enriched in actin, tubulin, or vimentin, and not associated with centrosomes or other known cytoplasmic structures. Further independent studies revealed that two key enzymes in the nucleotide synthetic pathway cytidine triphosphate synthase 1 (CTPS1) and inosine monophosphate dehydrogenase 2 (IMPDH2) were highly enriched in RR. CTPS1 enzyme inhibitors 6-diazo-5-oxo-L-norleucine and Acivicin as well as the IMPDH2 inhibitor Ribavirin exhibited dose-dependent induction of RR in >95% of cells in all cancer cell lines tested as well as mouse primary cells. RR formation by lower concentration of Ribavirin was enhanced in IMPDH2-knockdown HeLa cells whereas it was inhibited in GFP-IMPDH2 overexpressed HeLa cells. Interestingly, RR were detected readily in untreated mouse embryonic stem cells (>95%); upon retinoic acid differentiation, RR disassembled in these cells but reformed when treated with Acivicin. RR formation represented response to disturbances in the CTP or GTP synthetic pathways in cancer cell lines and mouse primary cells and RR are the convergence physical structures in these pathways. The availability of specific markers for these conserved structures and the ability to induce formation in vitro will allow further investigations in structure and function of RR in many biological systems in health and diseases.