Broodstock nutritional programming differentially affects the hepatic transcriptome and genome-wide DNA methylome of farmed gilthead sea bream (Sparus aurata) depending on genetic background.

BACKGROUND: Broodstock nutritional programming improves the offspring utilization of plant-based diets in gilthead sea bream through changes in hepatic metabolism. Attention was initially focused on fatty acid desaturases, but it can involve a wide range of processes that remain largely unexplored. How all this can be driven by a different genetic background is hardly underlined, and the present study aimed to assess how broodstock nutrition affects differentially the transcriptome and genome-wide DNA methylome of reference and genetically selected fish within the PROGENSA® selection program. RESULTS: After the stimulus phase with a low fish oil diet, two offspring subsets of each genetic background received a control or a FUTURE-based diet. This highlighted a different hepatic transcriptome (RNA-seq) and genome-wide DNA methylation (MBD-seq) pattern depending on the genetic background. The number of differentially expressed transcripts following the challenge phase varied from 323 in reference fish to 2,009 in genetically selected fish. The number of discriminant transcripts, and associated enriched functions, were also markedly higher in selected fish. Moreover, correlation analysis depicted a hyper-methylated and down-regulated gene expression state in selected fish with the FUTURE diet, whereas the opposite pattern appeared in reference fish. After filtering for highly represented functions in selected fish, 115 epigenetic markers were retrieved in this group. Among them, lipid metabolism genes (23) were the most reactive following ordering by fold-change in expression, rendering a final list of 10 top markers with a key role on hepatic lipogenesis and fatty acid metabolism (cd36, pitpna, cidea, fasn, g6pd, lipt1, scd1a, acsbg2, acsl14, acsbg2). CONCLUSIONS: Gene expression profiles and methylation signatures were dependent on genetic background in our experimental model. Such assumption affected the magnitude, but also the type and direction of change. Thus, the resulting epigenetic clock of reference fish might depict an older phenotype with a lower methylation for the epigenetically responsive genes with a negative methylation-expression pattern. Therefore, epigenetic markers will be specific of each genetic lineage, serving the broodstock programming in our selected fish to prevent and mitigate later in life the risk of hepatic steatosis through changes in hepatic lipogenesis and fatty acid metabolism.

Epigenetic analysis in a murine genetic model of Gulf War illness.

Of the nearly 1 million military personnel who participated in the 1990-1991 Gulf War, between 25% and 35% became ill with what now is referred to as Gulf War Illness (GWI) by the Department of Defense. Symptoms varied from gastrointestinal distress to lethargy, memory loss, inability to concentrate, depression, respiratory, and reproductive problems. The symptoms have persisted for 30 years in those afflicted but the basis of the illness remains largely unknown. Nerve agents and other chemical exposures in the war zone have been implicated but the long-term effects of these acute exposures have left few if any identifiable signatures. The major aim of this study is to elucidate the possible genomic basis for the persistence of symptoms, especially of the neurological and behavioral effects. To address this, we performed a whole genome epigenetic analysis of the proposed cause of GWI, viz., exposure to organophosphate neurotoxicants combined with high circulating glucocorticoids in two inbred mouse strains, C57BL/6J and DBA/2J. The animals received corticosterone in their drinking water for 7 days followed by injection of diisopropylfluorophosphate, a nerve agent surrogate. Six weeks after DFP injection, the animals were euthanized and medial prefrontal cortex harvested for genome-wide DNA methylation analysis using high-throughput sequencing. We observed 67 differentially methylated genes, notably among them, Ttll7, Akr1c14, Slc44a4, and Rusc2, all related to different symptoms of GWI. Our results support proof of principle of genetic differences in the chronic effects of GWI-related exposures and may reveal why the disease has persisted in many of the now aging Gulf War veterans.

Transgenerational Transcriptomic and DNA Methylome Profiling of Mouse Fetal Testicular Germline and Somatic Cells after Exposure of Pregnant Mothers to Tributyltin, a Potent Obesogen.

Obesogens such as tributyltin (TBT) are xenobiotic compounds that promote obesity, in part by distorting the normal balance of lipid metabolism. The obesogenic effects of TBT can be observed in directly exposed (F1 and F2 generations) and also subsequent generations (F3 and beyond) that were never exposed. To address the effects of TBT exposure on germ cells, we exposed pregnant transgenic OG2 mouse dams (F0), which specifically express EGFP in germline cells, to an environmentally relevant dose of TBT or DMSO throughout gestation through drinking water. When fed with a high-fat diet, F3 male offspring of TBT-exposed F0 dams (TBT-F3) accumulated much more body fat than did DMSO-F3 males. TBT-F3 males also lost more body fluid and lean compositions than did DMSO-F3 males. Expression of genes involved in transcriptional regulation or mesenchymal differentiation was up-regulated in somatic cells of TBT-F1 (but not TBT-F3) E18.5 fetal testes, and promoter-associated CpG islands were hyper-methylated in TBT-F1 somatic cells. Global mRNA expression of protein-coding genes in F1 or F3 fetal testicular cells was unaffected by F0 exposure to TBT; however, expression of a subset of endogenous retroviruses was significantly affected in F1 and F3. We infer that TBT may directly target testicular somatic cells in F1 testes to irreversibly affect epigenetic suppression of endogenous retroviruses in both germline and somatic cells.

Computational challenges in detection of cancer using cell-free DNA methylation.

Cell-free DNA(cfDNA) methylation profiling is considered promising and potentially reliable for liquid biopsy to study progress of diseases and develop reliable and consistent diagnostic and prognostic biomarkers. There are several different mechanisms responsible for the release of cfDNA in blood plasma, and henceforth it can provide information regarding dynamic changes in the human body. Due to the fragmented nature, low concentration of cfDNA, and high background noise, there are several challenges in its analysis for regular use in diagnosis of cancer. Such challenges in the analysis of the methylation profile of cfDNA are further aggravated due to heterogeneity, biomarker sensitivity, platform biases, and batch effects. This review delineates the origin of cfDNA methylation, its profiling, and associated computational problems in analysis for diagnosis. Here we also contemplate upon the multi-marker approach to handle the scenario of cancer heterogeneity and explore the utility of markers for 5hmC based cfDNA methylation pattern. Further, we provide a critical overview of deconvolution and machine learning methods for cfDNA methylation analysis. Our review of current methods reveals the potential for further improvement in analysis strategies for detecting early cancer using cfDNA methylation.

Cell-free DNA methylome profiling by MBD-seq with ultra-low input.

Methylation signatures in cell-free DNA (cfDNA) have shown great sensitivity and specificity in the characterization of tumour status and classification of tumour types, as well as the response to therapy and recurrence. Currently, most cfDNA methylation studies are based on bisulphite conversion, especially targeted bisulphite sequencing, while enrichment-based methods such as cfMeDIP-seq are beginning to show potential. Here, we report an enrichment-based ultra-low input cfDNA methylation profiling method using methyl-CpG binding proteins capture, termed cfMBD-seq. We optimized the conditions for cfMBD capture by adjusting the amount of MethylCap protein along with using methylated filler DNA. Our data show high correlation between low input cfMBD-seq and standard MBD-seq (>1000 ng input). When compared to cfMEDIP-seq, cfMBD-seq demonstrates higher sequencing data quality with more sequenced reads passed filter and less duplicate rate. cfMBD-seq also outperforms cfMeDIP-seq in the enrichment of CpG islands. This new bisulphite-free ultra-low input methylation profiling technology has great potential in non-invasive and cost-effective cancer detection and classification.

Identification of a novel differentially methylated region adjacent to ATG16L2 in lung cancer cells using methyl-CpG binding domain protein-enriched genome sequencing.

Lung cancer is the most common cancer worldwide. Epigenetic modifications like DNA methylation play fundamental roles in the dynamic process of lung cancer. The objective of this study was to use methyl-CpG binding domain protein-enriched genome sequencing (MBD-Seq) to identify novel and high-confidence DNA methylation in lung tumor. We first compared the whole-genome DNA methylation of three lung cancer cell lines, including A549, H1299, and SK-MES-1, against BEAS-2B, a lung/bronchial normal epithelial cell line. We then used pyrosequencing and OneStep qMethyl kit methods to verify the results in the cell line specimens. MBD-Seq identified 279, 8046, and 22 887 differentially methylated regions (DMRs), respectively, with 120 common DMRs among three comparison groups. Three DMRs were consistent with the MBD-Seq results by both pyrosequencing and OneStep qMethyl validations. Furthermore, OneStep qMethyl kit was also performed for functional validation of these three potential DMRs in sputum DNA from clinical participants. We successfully identified one new DMR adjacent to ATG16L2. The novel DMR might have an important function in lung carcinogenesis. Further validation of the finding in clinical specimens of lung cancer patients and functional analysis of this novel DMR in the development of lung cancer through transcriptional silencing of ATG16L2 are warranted.

The Profiling of DNA Methylation and Its Regulation on Divergent Tenderness in Angus Beef Cattle.

Beef is an essential food source in the world. Beef quality, especially tenderness, has a significant impact on consumer satisfaction and industry profit. Many types of research to date have focused on the exploration of physiological and developmental mechanisms of beef tenderness. Still, the role and impact of DNA methylation status on beef tenderness have yet to be elucidated. In this study, we exhaustively analyzed the DNA methylation status in divergent tenderness observed in Angus beef. We characterized the methylation profiles related to beef tenderness and explored methylation distributions on the whole genome. As a result, differentially methylated regions (DMRs) associated with tenderness and toughness of beef were identified. Importantly, we annotated these DMRs on the bovine genome and explored bio-pathways of underlying genes and methylation biomarkers in beef quality. Specifically, we observed that the ATP binding cassette subfamily and myosin-related genes were highly methylated gene sets, and generation of neurons, regulation of GTPase activity, ion transport and anion transport, etc., were the significant pathways related with beef tenderness. Moreover, we explored the relationship between DNA methylation and gene expression in DMRs. Some methylated genes were identified as candidate biomarkers for beef tenderness. These results provide not only novel epigenetic information associated with beef quality but offer more significant insights into meat science, which will further help us explore the mechanism of muscle biology.

MBD-seq - realities of a misunderstood method for high-quality methylome-wide association studies.

The majority of methylome-wide association studies (MWAS) have been performed using commercially available array-based technologies such as the Infinium Human Methylation 450K and the Infinium MethylationEPIC arrays (Illumina). While these arrays offer a convenient and relatively robust assessment of the probed sites they only allow interrogation of 2-4% of all CpG sites in the human genome. Methyl-binding domain sequencing (MBD-seq) is an alternative approach for MWAS that provides near-complete coverage of the methylome at similar costs as the array-based technologies. However, despite publication of multiple positive evaluations, the use of MBD-seq for MWAS is often fiercely criticized. Here we discuss key features of the method and debunk misconceptions using empirical data. We conclude that MBD-seq represents an excellent approach for large-scale MWAS and that increased utilization is likely to result in more discoveries, advance biological knowledge, and expedite the clinical translation of methylome-wide research findings.

Independent Methylome-Wide Association Studies of Schizophrenia Detect Consistent Case-Control Differences.

Methylome-wide association studies (MWASs) are promising complements to sequence variation studies. We used existing sequencing-based methylation data, which assayed the majority of all 28 million CpGs in the human genome, to perform an MWAS for schizophrenia in blood, while controlling for cell-type heterogeneity with a recently generated platform-specific reference panel. Next, we compared the MWAS results with findings from 3 existing large-scale array-based schizophrenia methylation studies in blood that assayed up to ~450 000 CpGs. Our MWAS identified 22 highly significant loci (P < 5 × 10-8) and 852 suggestively significant loci (P < 1 × 10-5). The top finding (P = 5.62 × 10-11, q = 0.001) was located in MFN2, which encodes mitofusin-2 that regulates Ca2+ transfer from the endoplasmic reticulum to mitochondria in cooperation with DISC1. The second-most significant site (P = 1.38 × 10-9, q = 0.013) was located in ALDH1A2, which encodes an enzyme for astrocyte-derived retinoic acid-a key neuronal morphogen with relevance for schizophrenia. Although the most significant MWAS findings were not assayed on the arrays, we observed significant enrichment of overlapping findings with 2 of the 3 array datasets (P = 0.0315, 0.0045, 0.1946). Overrepresentation analysis of Gene Ontology terms for the genes in the significant overlaps suggested high similarity in the biological functions detected by the different datasets. Top terms were related to immune and/or stress responses, cell adhesion and motility, and a broad range of processes essential for neurodevelopment.

Demethylation and derepression of genomic retroelements in the skeletal muscles of aged mice.

Changes in DNA methylation influence the aging process and contribute to aging phenotypes, but few studies have been conducted on DNA methylation changes in conjunction with skeletal muscle aging. We explored the DNA methylation changes in a variety of retroelement families throughout aging (at 2, 20, and 28 months of age) in murine skeletal muscles by methyl-binding domain sequencing (MBD-seq). The two following contrasting patterns were observed among the members of each repeat family in superaged mice: (a) hypermethylation in weakly methylated retroelement copies and (b) hypomethylation in copies with relatively stronger methylation levels, representing a pattern of "regression toward the mean" within a single retroelement family. Interestingly, these patterns depended on the sizes of the copies. While the majority of the elements showed a slight increase in methylation, the larger copies (>5 kb) displayed evident demethylation. All these changes were not observed in T cells. RNA sequencing revealed a global derepression of retroelements during the late phase of aging (between 20 and 28 months of age), which temporally coincided with retroelement demethylation. Following this methylation drift trend of "regression toward the mean," aging tended to progressively lose the preexisting methylation differences and local patterns in the genomic regions that had been elaborately established during the early period of development.

Genome-wide DNA methylation changes in transformed foci induced by nongenotoxic carcinogens.

In vitro cell transformation assays (CTA) have been proposed as a method to identify possible nongenotoxic carcinogens. However, the current protocols do not provide information on the mechanism of action of the test articles. In this study, we combined an in vitro Bhas 42 CTA and sequencing-based DNA methylation profiling analysis to elucidate the carcinogenic mechanism associated with nongenotoxic carcinogens. Three nongenotoxic carcinogens were evaluated: cadmium chloride, methyl carbamate, and lithocholic acid. Methylation profiles were generated for the two nongenotoxic carcinogens (cadmium chloride and lithocholic acid) that were positive in Bhas 42 CTA. Methyl carbamate did not exhibit any promoter activity. Approximately 9.8% of all differentially methylated regions (DMRs) identified in cadmium chloride-induced transformed foci overlapped with DMRs in lithocholic acid-induced transformed foci. Interestingly, overlapping DMRs showed more hypermethylation than individual DMRs. In addition, the DMRs in CpG island elements common to both nongenotoxic carcinogens showed considerably more bias toward hypermethylated DMRs than those unique to either cadmium chloride or lithocholic acid. Pathway enrichment analysis revealed that genes harboring hypermethylated DMRs were significantly enriched in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways including pathways in cancer, basal cell carcinoma, and Wnt signaling. The genes harboring hypomethylated DMRs were significantly related to mRNA surveillance pathway, RNA transport, and autophagy. Taken together, our preliminary results on genome-wide methylation analysis of cell clones from nongenotoxic carcinogen-induced foci could be exploited for CTAs improvement, but further research will be required to standardize and assess the specificity and sensitivity of this combined approach. Environ. Mol. Mutagen. 2019. © 2019 Wiley Periodicals, Inc.

Transcription factors Tp73, Cebpd, Pax6, and Spi1 rather than DNA methylation regulate chronic transcriptomics changes after experimental traumatic brain injury.

Traumatic brain injury (TBI) induces a wide variety of cellular and molecular changes that can continue for days to weeks to months, leading to functional impairments. Currently, there are no pharmacotherapies in clinical use that favorably modify the post-TBI outcome, due in part to limited understanding of the mechanisms of TBI-induced pathologies. Our system biology analysis tested the hypothesis that chronic transcriptomics changes induced by TBI are controlled by altered DNA-methylation in gene promoter areas or by transcription factors. We performed genome-wide methyl binding domain (MBD)-sequencing (seq) and RNA-seq in perilesional, thalamic, and hippocampal tissue sampled at 3 months after TBI induced by lateral fluid percussion in adult male Sprague-Dawley rats. We investigated the regulated molecular networks and mechanisms underlying the chronic regulation, particularly DNA methylation and transcription factors. Finally, we identified compounds that modulate the transcriptomics changes and could be repurposed to improve recovery. Unexpectedly, DNA methylation was not a major regulator of chronic post-TBI transcriptomics changes. On the other hand, the transcription factors Cebpd, Pax6, Spi1, and Tp73 were upregulated at 3 months after TBI (False discovery rate < 0.05), which was validated using digital droplet polymerase chain reaction. Transcription regulatory network analysis revealed that these transcription factors regulate apoptosis, inflammation, and microglia, which are well-known contributors to secondary damage after TBI. Library of Integrated Network-based Cellular Signatures (LINCS) analysis identified 118 pharmacotherapies that regulate the expression of Cebpd, Pax6, Spi1, and Tp73. Of these, the antidepressant and/or antipsychotic compounds trimipramine, rolipramine, fluspirilene, and chlorpromazine, as well as the anti-cancer therapies pimasertib, tamoxifen, and vorinostat were strong regulators of the identified transcription factors, suggesting their potential to modulate the regulated transcriptomics networks to improve post-TBI recovery.

Methyl-CpG-Binding Domain Sequencing: MBD-seq.

Detailed biological knowledge about the potential importance of the methylome is typically lacking for common diseases. Therefore, methylome-wide association studies (MWAS) are critical to detect disease relevant methylation sites. Methyl-CpG-binding domain sequencing (MBD-seq) offers potential advantages compared to antibody-based enrichment, but performance depends critically on using an optimal protocol. Using an optimized protocol, MBD-seq can approximate the sensitivity/specificity obtained with whole-genome bisulfite sequencing, but at a fraction of the costs and time to complete the project. Thus, MBD-seq offers a comprehensive first pass at the CpG methylome and is economically feasible with the samples sizes required for MWAS.

Developmental dynamics of the epigenome: A longitudinal study of three toddlers.

Epigenetic regulation plays an important role in development, at the embryonic stages and later during the lifespan. Some epigenetic marks are highly conserved throughout the lifespan whereas others are closely associated with specific age periods and/or particular environmental factors. Little is known about the dynamics of epigenetic regulation during childhood, especially during the period of rapid early development. Our study was aimed to determine whether the developmental program at the early stages of human development is accompanied by significant changes in the systems of genome regulation, specifically, by genome-wide changes in DNA methylation. Using a sequencing approach (MBD-seq) we investigated genome-wide DNA methylation patterns in the T-lymphocytes of three healthy toddlers at two timepoints within the second year of life. Pairwise comparison of the methylation patterns across the individuals and time points was conducted to determine common longitudinal changes in the DNA methylation patterns. Despite relatively high interindividual variability in their epigenetic profiles and the dynamics of these profiles during the second year of life, all children showed consistent changes in the DNA methylation patterns of genes involved in the control of the immune system and genes related to the development of the CNS. Thereby, we provide evidence that early development might be accompanied by epigenetic changes in specific functional groups of genes; many such epigenetic changes appear to be related to the rapid development of the CNS.

Genome-wide analysis of DNA methylation in hypothalamus and ovary of Capra hircus.

BACKGROUND: DNA methylation is a frequently studied epigenetic modification due to its role in regulating gene expression and hence in biological processes and in determining phenotypic plasticity in organisms. Rudimentary DNA methylation patterns for some livestock species are publically available: among these, goat methylome deserves to be further explored. RESULTS: Genome-wide DNA methylation maps of the hypothalamus and ovary from Saanen goats were generated using Methyl-CpG binding domain protein sequencing (MBD-seq). Analysis of DNA methylation patterns indicate that the majority of methylation peaks found within genes are located gene body regions, for both organs. Analysis of the distribution of methylated sites per chromosome showed that chromosome X had the lowest number of methylation peaks. The X chromosome has one of the highest percentages of methylated CpG islands in both organs, and approximately 50% of the CpG islands in the goat epigenome are methylated in hypothalamus and ovary. Organ-specific Differentially Methylated Genes (DMGs) were correlated with the expression levels. CONCLUSIONS: The comparison between transcriptome and methylome in hypothalamus and ovary showed that a higher level of methylation is not accompanied by a higher gene suppression. The genome-wide DNA methylation map for two goat organs produced here is a valuable starting point for studying the involvement of epigenetic modifications in regulating goat reproduction performance.

Comparative analysis of MBD-seq and MeDIP-seq and estimation of gene expression changes in a rodent model of schizophrenia.

We conducted a comparative study of multiplexed affinity enrichment sequence methodologies (MBD-seq and MeDIP-seq) in a rodent model of schizophrenia, induced by in utero methylazoxymethanol acetate (MAM) exposure. We also examined related gene expression changes using a pooled sample approach. MBD-seq and MeDIP-seq identified 769 and 1771 differentially methylated regions (DMRs) between F2 offspring of MAM-exposed rats and saline control rats, respectively. The assays showed good concordance, with ~56% of MBD-seq-detected DMRs being identified by or proximal to MeDIP-seq DMRs. There was no significant overlap between DMRs and differentially expressed genes, suggesting that DNA methylation regulatory effects may act upon more distal genes, or are too subtle to detect using our approach. Methylation and gene expression gene ontology enrichment analyses identified biological processes important to schizophrenia pathophysiology, including neuron differentiation, prepulse inhibition, amphetamine response, and glutamatergic synaptic transmission regulation, reinforcing the utility of the MAM rodent model for schizophrenia research.

A Whole Methylome CpG-SNP Association Study of Psychosis in Blood and Brain Tissue.

Mutated CpG sites (CpG-SNPs) are potential hotspots for human diseases because in addition to the sequence variation they may show individual differences in DNA methylation. We performed methylome-wide association studies (MWAS) to test whether methylation differences at those sites were associated with schizophrenia. We assayed all common CpG-SNPs with methyl-CpG binding domain protein-enriched genome sequencing (MBD-seq) using DNA extracted from 1408 blood samples and 66 postmortem brain samples (BA10) of schizophrenia cases and controls. Seven CpG-SNPs passed our FDR threshold of 0.1 in the blood MWAS. Of the CpG-SNPs methylated in brain, 94% were also methylated in blood. This significantly exceeded the 46.2% overlap expected by chance (P-value < 1.0×10(-8)) and justified replicating findings from blood in brain tissue. CpG-SNP rs3796293 in IL1RAP replicated (P-value = .003) with the same direction of effects. This site was further validated through targeted bisulfite pyrosequencing in 736 independent case-control blood samples (P-value < 9.5×10(-4)). Our top result in the brain MWAS (P-value = 8.8×10(-7)) was CpG-SNP rs16872141 located in the potential promoter of ENC1. Overall, our results suggested that CpG-SNP methylation may reflect effects of environmental insults and can provide biomarkers in blood that could potentially improve disease management.

Paediatric leukaemia DNA methylation profiling using MBD enrichment and SOLiD sequencing on archival bone marrow smears.

BACKGROUND: Acute Lymphoblastic Leukaemia (ALL) is the most common cancer in children. Over the past four decades, research has advanced the treatment of this cancer from a less than 60% chance of survival to over 85% today. The causal molecular mechanisms remain unclear. Here, we performed sequencing-based genomic DNA methylation profiling of eight paediatric ALL patients using archived bone marrow smear microscope slides. FINDINGS: SOLiD™ sequencing data was collected from Methyl-Binding Domain (MBD) enriched fractions of genomic DNA. The primary tumour and remission bone marrow sample was analysed from eight patients. Four patients relapsed and the relapsed tumour was analysed. Input and MBD-enriched DNA from each sample was sequenced, aligned to the hg19 reference genome and analysed for enrichment peaks using MACS (Model-based Analysis for ChIP-Seq) and HOMER (Hypergeometric Optimization of Motif EnRichment). In total, 3.67 gigabases (Gb) were sequenced, 2.74 Gb were aligned to the reference genome (average 74.66% alignment efficiency). This dataset enables the interrogation of differential DNA methylation associated with paediatric ALL. Preliminary results reveal concordant regions of enrichment indicative of a DNA methylation signature. CONCLUSION: Our dataset represents one of the first SOLiD™MBD-Seq studies performed on paediatric ALL and is the first to utilise archival bone marrow smears. Differential DNA methylation between cancer and equivalent disease-free tissue can be identified and correlated with existing and published genomic studies. Given the rarity of paediatric haematopoietic malignancies, relative to adult counterparts, our demonstration of the utility of archived bone marrow smear samples to high-throughput methylation sequencing approaches offers tremendous potential to explore the role of DNA methylation in the aetiology of cancer.

Mining for viral fragments in methylation enriched sequencing data.

Most next generation sequencing experiments generate more data than is usable for the experimental set up. For example, methyl-CpG binding domain (MBD) affinity purification based sequencing is often used for DNA-methylation profiling, but up to 30% of the sequenced fragments cannot be mapped uniquely to the reference genome. Here we present and evaluate a methodology for the identification of viruses in these otherwise unused paired-end MBD-seq data. Viral detection is accomplished by mapping non-reference alignable reads to a comprehensive set of viral genomes. As viruses play an important role in epigenetics and cancer development, 92 (pre)malignant and benign samples, originating from two different collections of cervical samples and related cell lines, were used in this study. These samples include primary carcinomas (n = 22), low- and high-grade cervical intraepithelial neoplasia (CIN1 and CIN2/3 - n = 2/n = 30) and normal tissue (n = 20), as well as control samples (n = 17). Viruses that were detected include phages, adenoviruses, herpesviridae and HPV. HPV, which causes virtually all cervical cancers, was identified in 95% of the carcinomas, 100% of the CIN2/3 samples, both CIN1 samples and in 55% of the normal samples. Comparing the amount of mapped fragments on HPV for each HPV-infected sample yielded a significant difference between normal samples and carcinomas or CIN2/3 samples (adjusted p-values resp. <10(-5), <10(-5)), reflecting different viral loads and/or methylation degrees in non-normal samples. Fragments originating from different HPV types could be distinguished and were independently validated by PCR-based assays in 71% of the detections. In conclusion, although limited by the a priori knowledge of viral reference genome sequences, the proposed methodology can provide a first confined but substantial insight into the presence, concentration and types of methylated viral sequences in MBD-seq data at low additional cost.