MethCP: Differentially Methylated Region Detection with Change Point Models

Whole-genome bisulfite sequencing (WGBS) provides a precise measure of methylation across the genome, yet presents a challenge in identifying differentially methylated regions (DMRs) between different conditions. Many methods have been developed, which focus primarily on the setting of two-group comparison. We develop a DMR detecting method MethCP for WGBS data, which is applicable for a wide range of experimental designs beyond the two-group comparisons, such as time-course data. MethCP identifies DMRs based on change point detection, which naturally segments the genome and provides region-level differential analysis. For simple two-group comparison, we show that our method outperforms developed methods in accurately detecting the complete DMR on a simulated data set and an Arabidopsis data set. Moreover, we show that MethCP is capable of detecting wide regions with small effect sizes, which can be common in some settings, but existing techniques are poor in detecting such DMRs. We also demonstrate the use of MethCP for time-course data on another data set after methylation throughout seed germination in Arabidopsis.

Epigenetic Footprints of CRISPR/Cas9-Mediated Genome Editing in Plants

CRISPR/Cas9 has been widely applied to various plant species accelerating the pace of plant genome editing and precision breeding in crops. Unintended effects beyond off-target nucleotide mutations are still somewhat unexplored. We investigated the degree and patterns of epigenetic changes after gene editing. We examined changes in DNA methylation in genome-edited promoters of naturally hypermethylated genes (AT1G72350 and AT1G09970) and hypomethylated genes (AT3G17320 and AT5G28770) from Arabidopsis. Transgenic plants were developed via Agrobacterium-mediated floral dip transformation. Homozygous edited lines were selected from segregated T2 plants by an in vitro digestion assay using ribonucleoprotein complex. Bisulfite sequencing comparisons were made between paired groups of edited and non-edited plants to identify changes in DNA methylation of the targeted loci. We found that directed mutagenesis via CRISPR/Cas9 resulted in no unintended morphological or epigenetic alterations. Phenotypes of wild-type, transgenic empty vector, and transgenic edited plants were similar. Epigenetic profiles revealed that methylation patterns of promoter regions flanking target sequences were identical among wild-type, transgenic empty vector, and transgenic edited plants. There was no effect of mutation type on epigenetic status. We also evaluated off-target mutagenesis effects in the edited plants. Potential off-target sites containing up to 4-bp mismatch of each target were sequenced. No off-target mutations were detected in candidate sites. Our results showed that CRISPR/Cas9 did not leave an epigenetic footprint on either the immediate gene-edited DNA and flanking DNA or introduce off-target mutations.

Differential methylation analysis for bisulfite sequencing using DSS

Bisulfite sequencing (BS-seq) technology measures DNA methylation at single nucleotide resolution. A key task in BS-seq data analysis is to identify differentially methylation (DM) under different conditions. Here we provide a tutorial for BS-seq DM analysis using Bioconductor package DSS. DSS uses a beta-binomial model to characterize the sequence counts from BS-seq, and implements rigorous statistical method for hypothesis testing. It provides flexible functionalities for a variety of DM analyses.

13-gene DNA Methylation Analysis from Oral Brushing: A Promising Non Invasive Tool in the Follow-up of Oral Cancer Patients

Background: This study aimed to evaluate the prognostic value of a non-invasive sampling procedure based on 13-gene DNA methylation analysis in the follow-up of patients previously treated for oral squamous cell carcinoma (OSCC). Methods: The study population included 49 consecutive patients treated for OSCC. Oral brushing sample collection was performed at two different times: before any cancer treatment in the tumor mass and during patient follow-up almost 6 months after OSCC treatment, within the regenerative area after OSCC resection. Each sample was considered positive or negative in relation to a predefined cut-off value. Results: Before any cancer treatment, 47/49 specimens exceeded the score and were considered as positive. Six months after OSCC resection, 16/49 specimens also had positive scores in the samples collected from the regenerative area. During the follow-up period, 7/49 patients developed locoregional relapse: 6/7 patients had a positive score in the regenerative area after OSCC resection. The presence of a positive score after oral cancer treatment was the most powerful variable related to the appearance of locoregional relapse. Conclusion: 13-gene DNA methylation analysis by oral brushing may have a clinical application as a prognostic non-invasive tool in the follow-up of patients surgically treated for OSCC.

Circulating Tumor Cell Clustering Shapes DNA Methylation to Enable Metastasis Seeding

The ability of circulating tumor cells (CTCs) to form clusters has been linked to increased metastatic potential. Yet biological features and vulnerabilities of CTC clusters remain largely unknown. Here, we profile the DNA methylation landscape of single CTCs and CTC clusters from breast cancer patients and mouse models on a genome-wide scale. We find that binding sites for stemness- and proliferation-associated transcription factors are specifically hypomethylated in CTC clusters, including binding sites for OCT4, NANOG, SOX2, and SIN3A, paralleling embryonic stem cell biology. Among 2,486 FDA-approved compounds, we identify Na⁺/K⁺ ATPase inhibitors that enable the dissociation of CTC clusters into single cells, leading to DNA methylation remodeling at critical sites and metastasis suppression. Thus, our results link CTC clustering to specific changes in DNA methylation that promote stemness and metastasis and point to cluster-targeting compounds to suppress the spread of cancer.

Whole-Genome Bisulfite Sequencing Reveals a Role for DNA Methylation in Variants from Callus Culture of Pineapple (Ananas comosus L.)

DNA methylation changes can occur in some loci during callus culture, resulting in somaclonal variations (SVs). In the present study, we applied whole genome bisulfite sequencing to analyze context-specific DNA methylation changes in the pineapple genome between the cutting seedings and 5 SV plants. In general, SV plants exhibited methylation patterns analogous to those of cutting seedlings (CK). A total of 27.98% of the genomic cytosines of CK were methylcytosines, which was higher than that of 5 SV plants. Moreover, mCG and mCHG was hypermethylated, whereas mCHH was hypomethylated among the 5 SV plants genomic when compared with the CK. Most of the variation of DNA methylation was distributed in gene bodies, thus suggesting that phenotypic differences are probably perturbed by genes methylated from callus culture. In addition, the methylated genes were highly enriched for the Gene Ontology (GO) categories of binding and catalytic activity, cell part and organelle, cellular process, abiotic stimulus, and DNA modification. These results suggest that methylation mediates these pathways in the callus culture of pineapple. The results also suggested that the callus culture induced DNA methylation may result in the SV.

Genome-wide effects of social status on DNA methylation in the brain of a cichlid fish, Astatotilapia burtoni

BACKGROUND

Successful social behavior requires real-time integration of information about the environment, internal physiology, and past experience. The molecular substrates of this integration are poorly understood, but likely modulate neural plasticity and gene regulation. In the cichlid fish species Astatotilapia burtoni, male social status can shift rapidly depending on the environment, causing fast behavioral modifications and a cascade of changes in gene transcription, the brain, and the reproductive system. These changes can be permanent but are also reversible, implying the involvement of a robust but flexible mechanism that regulates plasticity based on internal and external conditions. One candidate mechanism is DNA methylation, which has been linked to social behavior in many species, including A. burtoni. But, the extent of its effects after A. burtoni social change were previously unknown.

RESULTS

We performed the first genome-wide search for DNA methylation patterns associated with social status in the brains of male A. burtoni, identifying hundreds of Differentially Methylated genomic Regions (DMRs) in dominant versus non-dominant fish. Most DMRs were inside genes supporting neural development, synapse function, and other processes relevant to neural plasticity, and DMRs could affect gene expression in multiple ways. DMR genes were more likely to be transcription factors, have a duplicate elsewhere in the genome, have an anti-sense lncRNA, and have more splice variants than other genes. Dozens of genes had multiple DMRs that were often seemingly positioned to regulate specific splice variants.

CONCLUSIONS

Our results revealed genome-wide effects of A. burtoni social status on DNA methylation in the brain and strongly suggest a role for methylation in modulating plasticity across multiple biological levels. They also suggest many novel hypotheses to address in mechanistic follow-up studies, and will be a rich resource for identifying the relationships between behavioral, neural, and transcriptional plasticity in the context of social status.

Reassessment of Viroid RNA Cytosine Methylation Status at the Single Nucleotide Level

Composed of a few hundreds of nucleotides, viroids are infectious, circular, non-protein coding RNAs able to usurp plant cellular enzymes and molecular machineries to replicate and move in their hosts. Several secondary and tertiary RNA structural motifs have been implicated in the viroid infectious cycle, but whether modified nucleotides, such as 5C-methylcytosine (m5C), also play a role has not been deeply investigated so far. Here, the possible existence of m5C in both RNA polarity strands of potato spindle tuber viroid and avocado sunblotch viroid -which are representative members of the nucleus- and chloroplast-replicating viroids, respectively- has been assessed at single nucleotide level. We show that a standard bisulfite protocol efficiently used for identifying m5C in cellular RNAs may generate false positive results in the case of the highly structured viroid RNAs. Applying a bisulfite conversion protocol specifically adapted to RNAs with high secondary structure, no m5C was identified in both polarity strands of both viroids, indicating that this specific nucleotide modification does not likely play a role in viroid biology.

Aberrant DNA methylation as a diagnostic biomarker of diabetic embryopathy

PURPOSE

Maternal diabetes is a known teratogen that can cause a wide spectrum of birth defects, collectively referred to as diabetic embryopathy (DE). However, the pathogenic mechanisms underlying DE remain uncertain and there are no definitive tests to establish the diagnosis. Here, we explore the potential of DNA methylation as a diagnostic biomarker for DE and to inform disease pathogenesis.

METHODS

Bisulfite sequencing was used to identify gene regions with differential methylation between DE neonates and healthy infants born with or without prenatal exposure to maternal diabetes, and to investigate the role of allele-specific methylation at implicated sites.

RESULTS

We identified a methylation signature consisting of 237 differentially methylated loci that distinguished infants with DE from control infants. These loci were found proximal to genes associated with Mendelian syndromes that overlap the DE phenotype (e.g., CACNA1C, TRIO, ANKRD11) or genes known to influence embryonic development (e.g., BRAX1, RASA3). Further, we identified allele-specific methylation (ASM) at 11 of these loci, within which 61.5% of ASM single-nucleotide variants are known expression quantitative trait loci (eQTLs).

CONCLUSIONS

Our study suggests a role for aberrant DNA methylation and cis-sequence variation in the pathogenesis of DE and highlights the diagnostic potential of DNA methylation for teratogenic birth defects.

Opportunities and limitations of reduced representation bisulfite sequencing in plant ecological epigenomics

Contents Summary 738 I. Introduction 738 II. RRBS loci as genome-wide epigenetic markers 739 III. Exploiting functional annotation of RRBS loci 739 IV. Limitations of RRBS methods for nonmodel species 740 V. Maximising the impact of RRBS in plants 741 VI. Conclusions 741 Acknowledgements 741 SUMMARY: Investigating the features and implications of epigenetic mechanisms across the breadth of organisms and ecosystems is important for understanding the ecological relevance of epigenetics. Several cost-effective reduced representation bisulfite sequencing approaches (RRBS) have been recently developed and applied to different organisms that lack a well annotated reference genome. These new approaches improve the assessment of epigenetic diversity in ecological settings and may provide functional insights. We assess here the opportunities and limitations of RRBS in nonmodel plant species. Well thought out experimental designs that include complementary gene expression studies, and the improvement of genomics resources for the target group, promise to maximize the effect of future RRBS studies.

[Suppression of NR0B2 gene in Clear Cell Renal Cell Carcinoma Is Associated with Hypermethylation of Its Promoter]

Clear cell renal cell carcinoma (ccRCC) is a common urologic malignancy. Understanding of the transcriptional regulation of oncogenes and tumor suppressor genes involved is critical for the development of the treatments for renal tumors. Using ccRCC subdivision of the TCGA dataset, we identified NR0B2 encoding orphan nuclear receptor as a tumor suppressor candidate in renal tissue. In independent cohort of primary renal tumors, quantitative PCR experiments confirmed significant suppression of NR0B2 mRNA in 86% of ccRCC samples studied. In 80% of these cases, we detected the hypermethylation of the NR0B2 pro-moter region. These results suggest that NR0B2 is a tumor suppressor gene in ccRCC, and that the hypermethylation of promoter region is the main mechanism of its downregulation.

BS-virus-finder: virus integration calling using bisulfite sequencing data

Background

DNA methylation plays a key role in the regulation of gene expression and carcinogenesis. Bisulfite sequencing studies mainly focus on calling single nucleotide polymorphism, different methylation region, and find allele-specific DNA methylation. Until now, only a few software tools have focused on virus integration using bisulfite sequencing data.

Findings

We have developed a new and easy-to-use software tool, named BS-virus-finder (BSVF, RRID:SCR_015727), to detect viral integration breakpoints in whole human genomes. The tool is hosted at https://github.com/BGI-SZ/BSVF.

Conclusions

BS-virus-finder demonstrates high sensitivity and specificity. It is useful in epigenetic studies and to reveal the relationship between viral integration and DNA methylation. BS-virus-finder is the first software tool to detect virus integration loci by using bisulfite sequencing data.

DNA Methylation Analysis of Dormancy Release in Almond (Prunus dulcis) Flower Buds Using Epi-Genotyping by Sequencing

DNA methylation and histone post-translational modifications have been described as epigenetic regulation mechanisms involved in developmental transitions in plants, including seasonal changes in fruit trees. In species like almond (Prunus dulcis (Mill.) D.A: Webb), prolonged exposure to cold temperatures is required for dormancy release and flowering. Aiming to identify genomic regions with differential methylation states in response to chill accumulation, we carried out Illumina reduced-representation genome sequencing on bisulfite-treated DNA from floral buds. To do this, we analyzed almond genotypes with different chilling requirements and flowering times both before and after dormancy release for two consecutive years. The study was performed using epi-Genotyping by Sequencing (epi-GBS). A total of 7317 fragments were sequenced and the samples compared. Out of these fragments, 677 were identified as differentially methylated between the almond genotypes. Mapping these fragments using the Prunus persica (L.) Batsch v.2 genome as reference provided information about coding regions linked to early and late flowering methylation markers. Additionally, the methylation state of ten gene-coding sequences was found to be linked to the dormancy release process.

DNA Methylation Analysis of Dormancy Release in Almond (Prunus dulcis) Flower Buds Using Epi-Genotyping by Sequencing

DNA methylation and histone post-translational modifications have been described as epigenetic regulation mechanisms involved in developmental transitions in plants, including seasonal changes in fruit trees. In species like almond (Prunus dulcis (Mill.) D.A: Webb), prolonged exposure to cold temperatures is required for dormancy release and flowering. Aiming to identify genomic regions with differential methylation states in response to chill accumulation, we carried out Illumina reduced-representation genome sequencing on bisulfite-treated DNA from floral buds. To do this, we analyzed almond genotypes with different chilling requirements and flowering times both before and after dormancy release for two consecutive years. The study was performed using epi-Genotyping by Sequencing (epi-GBS). A total of 7317 fragments were sequenced and the samples compared. Out of these fragments, 677 were identified as differentially methylated between the almond genotypes. Mapping these fragments using the Prunus persica (L.) Batsch v.2 genome as reference provided information about coding regions linked to early and late flowering methylation markers. Additionally, the methylation state of ten gene-coding sequences was found to be linked to the dormancy release process.

Impact of vitamin D depletion during development on mouse sperm DNA methylation

Suboptimal environmental conditions during development can substantially alter the epigenome. Stable environmentally-induced changes to the germline epigenome, in particular, have important implications for the health of the next generation. We showed previously that developmental vitamin D depletion (DVD) resulted in loss of DNA methylation at several imprinted loci over two generations. Here, we assessed the impact of DVD on genome-wide methylation in mouse sperm in order to characterize the number, extent and distribution of methylation changes in response to DVD and to find genes that may be susceptible to this prevalent environmental perturbation. We detected 15,827 loci that were differentially methylated in DVD mouse sperm vs. controls. Most epimutations (69%) were loss of methylation, and the extent of methylation change and number of CpGs affected in a region were associated with genic location and baseline methylation state. Methylation response to DVD at validated loci was only detected in offspring that exhibited a phenotypic response to DVD (increased body and testes weight) suggesting the two types of responses are linked, though a causal relationship is unclear. Epimutations localized to regions enriched for developmental and metabolic genes and pathway analyses showed enrichment for Cadherin, Wnt, PDGF and Integrin signaling pathways. These findings show for the first time that vitamin D status during development leads to substantial DNA methylation changes across the sperm genome and that locus susceptibility is linked to genomic and epigenomic context.

Transcriptome-wide analysis suggests that temporal changes in the relative contributions of hyperplasia, hypertrophy and apoptosis underlie liver growth in pregnant mice

Maternal liver undergoes structural and metabolic changes during pregnancy to meet the demands of the developing fetus. In rodents, this involves increased liver weight, but the mechanism remains unclear. To address this, we analyzed the histology, gene expression, and DNA methylation of livers of nonpregnant and pregnant C57/BL6 mice. Gestational liver growth in pregnant mice was accompanied by increased hepatocyte area and lower cell density (days 14 and 18). Expression of cell proliferation markers was increased on days 14 and 18. A total of 115 genes were differentially expressed on day 14 and 123 genes on day 18 (79 on both days). Pathway analysis indicated that pregnancy involves progressive increase in cell proliferation and decreased apoptosis. This was confirmed using archived data from the FVB wild-type mouse liver transcriptome. Four differentially DNA methylated and two differentially DNA hydroxymethylated regions identified on days 14 and 18 by methylome-wide analysis, but were not associated with altered gene expression. Long interspersed nuclear element-1 hypomethylation on days 14 and 18 was accompanied by increased ten-eleven translocase-2 and decreased DNA methyltransferase 3a and 3b expression. These findings suggest that gestational liver growth involves increased mitosis and hypertrophy, and decreased apoptosis contingent on pregnancy stage. Such changes may involve repetitive sequence, but not gene specific, DNA methylation.

DNA methylation regulates the neonatal CD4+ T-cell response to pneumonia in mice

Pediatric acute lung injury, usually because of pneumonia, has a mortality rate of more than 20% and an incidence that rivals that of all childhood cancers combined. CD4⁺ T-cells coordinate the immune response to pneumonia but fail to function robustly among the very young, who have poor outcomes from lung infection. We hypothesized that DNA methylation represses a mature CD4⁺ T-cell transcriptional program in neonates with pneumonia. Here, we found that neonatal mice (3-4 days old) aspirated with Escherichia coli bacteria had a higher mortality rate than juvenile mice (11-14 days old). Transcriptional profiling with an unsupervised RNA-Seq approach revealed that neonates displayed an attenuated lung CD4⁺ T-cell transcriptional response to pneumonia compared with juveniles. Unlike neonates, juveniles up-regulated a robust set of canonical T-cell immune response genes. DNA methylation profiling with modified reduced representation bisulfite sequencing revealed 44,119 differentially methylated CpGs, which preferentially clustered around transcriptional start sites and CpG islands. A methylation difference-filtering algorithm detected genes with a high likelihood of differential promoter methylation regulating their expression; these 731 loci encoded important immune response and tissue-protective T-cell pathway components. Disruption of DNA methylation with the hypomethylating agent decitabine induced plasticity in the lung CD4⁺ T-cell marker phenotype. Altogether, multidimensional profiling suggested that DNA methylation within the promoters of a core set of CD4⁺ T-cell pathway genes contributes to the hyporesponsive neonatal immune response to pneumonia. These findings also suggest that DNA methylation could serve as a mechanistic target for disease-modifying therapies in pediatric lung infection and injury.

Heritable Epigenomic Changes to the Maize Methylome Resulting from Tissue Culture

DNA methylation can contribute to the maintenance of genome integrity and regulation of gene expression. In most situations, DNA methylation patterns are inherited quite stably. However, changes in DNA methylation can occur at some loci as a result of tissue culture resulting in somaclonal variation. To investigate heritable epigenetic changes as a consequence of tissue culture, a sequence-capture bisulfite sequencing approach was implemented to monitor context-specific DNA methylation patterns in ∼15 Mb of the maize genome for a population of plants that had been regenerated from tissue culture. Plants that have been regenerated from tissue culture exhibit gains and losses of DNA methylation at a subset of genomic regions. There was evidence for a high rate of homozygous changes to DNA methylation levels that occur consistently in multiple independent tissue culture lines, suggesting that some loci are either targeted or hotspots for epigenetic variation. The consistent changes inherited following tissue culture include both gains and losses of DNA methylation and can affect CG, CHG, or both contexts within a region. Only a subset of the tissue culture changes observed in callus plants are observed in the primary regenerants, but the majority of DNA methylation changes present in primary regenerants are passed onto offspring. This study provides insights into the susceptibility of some loci and potential mechanisms that could contribute to altered DNA methylation and epigenetic state that occur during tissue culture in plant species.

Methodology for Accurate Detection of Mitochondrial DNA Methylation

Quantification of DNA methylation can be achieved using bisulfite sequencing, which takes advantage of the property of sodium bisulfite to convert unmethylated cytosine into uracil, in a single-stranded DNA context. Bisulfite sequencing can be targeted (using PCR) or performed on the whole genome and provides absolute quantification of cytosine methylation at the single base-resolution. Given the distinct nature of nuclear- and mitochondrial DNA, notably in the secondary structure, adaptions of bisulfite sequencing methods for investigating cytosine methylation in mtDNA should be made. Secondary and tertiary structure of mtDNA can indeed lead to bisulfite sequencing artifacts leading to false-positives due to incomplete denaturation poor access of bisulfite to single-stranded DNA. Here, we describe a protocol using an enzymatic digestion of DNA with BamHI coupled with bioinformatic analysis pipeline to allow accurate quantification of cytosine methylation levels in mtDNA. In addition, we provide guidelines for designing the bisulfite sequencing primers specific to mtDNA, in order to avoid targeting undesirable NUclear MiTochondrial segments (NUMTs) inserted into the nuclear genome.

Quantitative Analysis of the DNA Methylation Sensitivity of Transcription Factor Complexes

Although DNA modifications play an important role in gene regulation, the underlying mechanisms remain elusive. We developed EpiSELEX-seq to probe the sensitivity of transcription factor binding to DNA modification in vitro using massively parallel sequencing. Feature-based modeling quantifies the effect of cytosine methylation (⁵mC) on binding free energy in a position-specific manner. Application to the human bZIP proteins ATF4 and C/EBPβ and three different Pbx-Hox complexes shows that ⁵mCpG can both increase and decrease affinity, depending on where the modification occurs within the protein-DNA interface. The TF paralogs tested vary in their methylation sensitivity, for which we provide a structural rationale. We show that ⁵mCpG can also enhance in vitro p53 binding and provide evidence for increased in vivo p53 occupancy at methylated binding sites, correlating with primed enhancer histone marks. Our results establish a powerful strategy for dissecting the epigenomic modulation of protein-DNA interactions and their role in gene regulation.

Whole genome DNA methylation: beyond genes silencing

The combination of DNA bisulfite treatment with high-throughput sequencing technologies has enabled investigation of genome-wide DNA methylation at near base pair level resolution, far beyond that of the kilobase-long canonical CpG islands that initially revealed the biological relevance of this covalent DNA modification. The latest high-resolution studies have revealed a role for very punctual DNA methylation in chromatin plasticity, gene regulation and splicing. Here, we aim to outline the major biological consequences of DNA methylation recently discovered. We also discuss the necessity of tuning DNA methylation resolution into an adequate scale to ease the integration of the methylome information with other chromatin features and transcription events such as gene expression, nucleosome positioning, transcription factors binding dynamic, gene splicing and genomic imprinting. Finally, our review sheds light on DNA methylation heterogeneity in cell population and the different approaches used for its assessment, including the contribution of single cell DNA analysis technology.

MethyMer: Design of combinations of specific primers for bisulfite sequencing of complete CpG islands

We present MethyMer, a Python-based tool aimed at selecting primers for amplification of complete CpG islands. These regions are difficult in terms of selecting appropriate primers because of their low-complexity, high GC content. Moreover, bisulfite treatment, in fact, leads to the reduction of the 4-letter alphabet (ATGC) to 3-letter one (ATG, except for methylated cytosines), and this also reduces region complexity and increases mispriming potential. MethyMer has a flexible scoring system, which optimizes the balance between various characteristics such as nucleotide composition, thermodynamic features (melting temperature, dimers [Formula: see text]G, etc.), the presence of CpG sites and polyN tracts, and primer specificity, which is assessed with aligning primers to the bisulfite-treated genome using bowtie (up to three mismatches are allowed). Users are able to customize desired or limit ranges of various parameters as well as penalties for non-desired values. Moreover, MethyMer allows picking up the optimal combination of PCR primer pairs to perform the amplification of a large genomic locus, e.g. CpG island or other hard-to-study region, with minimal overlap of the individual amplicons. MethyMer incorporates ENCODE genome annotation records (promoter/enhancer/insulator), The Cancer Genome Atlas (TCGA) CpG methylation data derived with Illumina Infinium 450K microarrays, and records on correlations between TCGA RNA-Seq and CpG methylation data for 20 cancer types. These databases are included in the MethyMer release. Our tool is available at https://sourceforge.net/projects/methymer/ .

Intracellular African swine fever virus DNA remains unmethylated in infected Vero cells

AIM

Sequence-specific CpG methylation of eukaryotic promoters is an important epigenetic signal for long-term gene silencing. We have now studied the methylation status of African swine fever virus (ASFV) DNA at various times after infection of Vero cells in culture.

METHODS & RESULTS

ASFV DNA was detectable throughout the infection cycle and was found unmethylated in productively infected Vero cells as documented by bisulfite sequencing of 13 viral DNA segments.

CONCLUSION

ASFV DNA does not become de novo methylated in the course of infection in selected segments spread across the entire genome. Thus DNA methylation does not interfere with ASFV genome transcription. Lack of de novo methylation has previously been observed for free intracellular viral DNA in cells permissively infected with human adenoviruses, with human papillomaviruses and others.

Evidence Suggesting Absence of Mitochondrial DNA Methylation

Methylation of nuclear genes encoding mitochondrial proteins participates in the regulation of mitochondria function. The existence of cytosine methylation in the mitochondrial genome is debated. To investigate whether mitochondrial DNA (mtDNA) is methylated, we used both targeted- and whole mitochondrial genome bisulfite sequencing in cell lines and muscle tissue from mouse and human origin. While unconverted cytosines were detected in some portion of the mitochondrial genome, their abundance was inversely associated to the sequencing depth, indicating that sequencing analysis can bias the estimation of mtDNA methylation levels. In intact mtDNA, few cytosines remained 100% unconverted. However, removal of supercoiled structures of mtDNA with the restriction enzyme BamHI prior to bisulfite sequencing decreased cytosine unconversion rate to <1.5% at all the investigated regions: D-loop, tRNA-F+12S, 16S, ND5 and CYTB, suggesting that mtDNA supercoiled structure blocks the access to bisulfite conversion. Here, we identified an artifact of mtDNA bisulfite sequencing that can lead to an overestimation of mtDNA methylation levels. Our study supports that cytosine methylation is virtually absent in mtDNA.

BAT: Bisulfite Analysis Toolkit: BAT is a toolkit to analyze DNA methylation sequencing data accurately and reproducibly. It covers standard processing and analysis steps from raw read mapping up to annotation data integration and calculation of correlating DMRs

Here, we present BAT, a modular bisulfite analysis toolkit, that facilitates the analysis of bisulfite sequencing data. It covers the essential analysis steps of read alignment, quality control, extraction of methylation information, and calling of differentially methylated regions, as well as biologically relevant downstream analyses, such as data integration with gene expression, histone modification data, or transcription factor binding site annotation.