Recent advances in next generation sequencing (NGS) technology now provide the opportunity to rapidly interrogate the methylation status of the genome. However, there are challenges in handling and interpretation of the methylation sequence data because of its large volume and the consequences of bisulphite modification. We sequenced reduced representation human genomes on the Illumina platform and efficiently mapped and visualized the data with different pipelines and software packages. We examined three pipelines for aligning bisulphite converted sequencing reads and compared their performance. We also comment on pre-processing and quality control of Illumina data. This comparison highlights differences in methods for NGS data processing and provides guidance to advance sequence-based methylation data analysis for molecular biologists.

The sex drive hypothesis predicts that stronger selection on male traits has resulted in masculinization of the genome. Here we test whether such masculinizing effects can be detected at the level of the transcriptome and methylome in the adult zebrafish brain. Although methylation is globally similar, we identified 914 specific differentially methylated CpGs (DMCs) between males and females (435 were hypermethylated and 479 were hypomethylated in males compared to females). These DMCs were prevalent in gene body, intergenic regions and CpG island shores. We also discovered 15 distinct CpG clusters with striking sex-specific DNA methylation differences. In contrast, at transcriptome level, more female-biased genes than male-biased genes were expressed, giving little support for the male sex drive hypothesis. Our study provides genome-wide methylome and transcriptome assessment and sheds light on sex-specific epigenetic patterns and in zebrafish for the first time.

Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a complex disease with variable severity. Patients experience frequent relapses where symptoms increase in severity, leaving them with a marked reduction in quality of life. Previous work has investigated molecular differences between ME/CFS patients and healthy controls, but not the dynamic changes specific to each individual patient. We applied precision medicine here to map genomic changes in two selected ME/CFS patients through a period that contained a relapse recovery cycle. DNA was isolated from two patients and a healthy age/gender matched control at regular intervals and captured the patient relapse in each case. Reduced representation DNA methylation sequencing profiles were obtained spanning the relapse recovery cycle. Both patients showed a significantly larger methylome variability (10–20-fold) through the period of sampling compared with the control. During the relapse, changes in the methylome profiles of the two patients were detected in regulatory-active regions of the genome that were associated, respectively, with 157 and 127 downstream genes, indicating disturbed metabolic, immune and inflammatory functions. Severe health relapses in the ME/CFS patients resulted in functionally important changes in their DNA methylomes that, while differing between the two patients, led to very similar compromised physiology. DNA methylation as a signature of disease variability in ongoing ME/CFS may have practical applications for strategies to decrease relapse frequency.

Formalin fixed paraffin embedded (FFPE) tumor samples are a major source of DNA from patients in cancer research. However, FFPE is a challenging material to work with due to macromolecular fragmentation and nucleic acid crosslinking. FFPE tissue particularly possesses challenges for methylation analysis and for preparing sequencing-based libraries relying on bisulfite conversion. Successful bisulfite conversion is a key requirement for sequencing-based methylation analysis. Here we describe a complete and streamlined workflow for preparing next generation sequencing libraries for methylation analysis from FFPE tissues. This includes, counting cells from FFPE blocks and extracting DNA from FFPE slides, testing bisulfite conversion efficiency with a polymerase chain reaction (PCR) based test, preparing reduced representation bisulfite sequencing libraries and massively parallel sequencing. The main features and advantages of this protocol are: We provide a complete workflow and an integrated protocol for performing DNA methylation analysis at the genome-scale and we believe this will facilitate clinical epigenetic research that involves the use of FFPE tissue.

Melanoma is the most aggressive type of skin cancer, with increasing incidence worldwide. Advances in targeted therapy and immunotherapy have improved the survival of melanoma patients experiencing recurrent disease, but unfortunately treatment resistance frequently reduces patient survival. Resistance to targeted therapy is associated with transcriptomic changes and has also been shown to be accompanied by increased expression of programmed death ligand 1 (PD-L1), a potent inhibitor of immune response. Intrinsic upregulation of PD-L1 is associated with genome-wide DNA hypomethylation and widespread alterations in gene expression in melanoma cell lines. However, an in-depth analysis of the transcriptomic landscape of melanoma cells with intrinsically upregulated PD-L1 expression is lacking. To determine the transcriptomic landscape of intrinsically upregulated PD-L1 expression in melanoma, we investigated transcriptomes in melanomas with constitutive versus inducible PD-L1 expression (referred to as PD-L1_CON and PD-L1_IND). RNA-Seq analysis was performed on seven PD-L1_CON melanoma cell lines and ten melanoma cell lines with low inducible PD-L1_IND expression. We observed that PD-L1_CON melanoma cells had a reprogrammed transcriptome with a characteristic pattern of dedifferentiated gene expression, together with active interferon (IFN) and tumour necrosis factor (TNF) signalling pathways. Furthermore, we identified key transcription factors that were also differentially expressed in PD-L1_CON versus PD-L1_IND melanoma cell lines. Overall, our studies describe transcriptomic reprogramming of melanomas with PD-L1_CON expression.

DNA methylation is a heritable epigenetic mark that is fundamental to mammalian development. Aberrant DNA methylation is an epigenetic hallmark of cancer cells. Cell lines are a critical in vitro model and very widely used to unravel mechanisms of cancer cell biology. However, limited data are available to assess whether DNA methylation patterns in tissues are retained when cell lines are established. Here, we provide the first genome-scale sequencing-based methylation map of metastatic melanoma tumour tissues and their derivative cell lines. We show that DNA methylation profiles are globally conserved in vitro compared to the tumour tissue of origin. However, we identify sites that are consistently hypermethylated in cell lines compared to their tumour tissue of origin. The genes associated with these common differentially methylated regions are involved in cell metabolism, cell cycle and apoptosis and are also strongly enriched for the H3K27me3 histone mark and PRC2 complex-related genes. Our data indicate that although global methylation patterns are similar between tissues and cell lines, there are site-specific epigenomic differences that could potentially impact gene expression. Our work provides a valuable resource for identifying false positives due to cell culture and for better interpretation of cancer epigenetics studies in the future.