Reduced Representation Bisulphite Sequencing Research Articles

Cell-free epigenomes enhanced fragmentomics-based model for early detection of lung cancer.

Lung cancer is a leading cause of cancer mortality, highlighting the need for innovative non-invasive early detection methods. Although cell-free DNA (cfDNA) analysis shows promise, its sensitivity in early-stage lung cancer patients remains a challenge. This study aimed to integrate insights from epigenetic modifications and fragmentomic features of cfDNA using machine learning to develop a more accurate lung cancer detection model. To address this issue, a multi-centre prospective cohort study was conducted, with participants harbouring suspicious malignant lung nodules and healthy volunteers recruited from two clinical centres. Plasma cfDNA was analysed for its epigenetic and fragmentomic profiles using chromatin immunoprecipitation sequencing, reduced representation bisulphite sequencing and low-pass whole-genome sequencing. Machine learning algorithms were then employed to integrate the multi-omics data, aiding in the development of a precise lung cancer detection model. Cancer-related changes in cfDNA fragmentomics were significantly enriched in specific genes marked by cell-free epigenomes. A total of 609 genes were identified, and the corresponding cfDNA fragmentomic features were utilised to construct the ensemble model. This model achieved a sensitivity of 90.4% and a specificity of 83.1%, with an AUC of 0.94 in the independent validation set. Notably, the model demonstrated exceptional sensitivity for stage I lung cancer cases, achieving 95.1%. It also showed remarkable performance in detecting minimally invasive adenocarcinoma, with a sensitivity of 96.2%, highlighting its potential for early detection in clinical settings. With feature selection guided by multiple epigenetic sequencing approaches, the cfDNA fragmentomics-based machine learning model demonstrated outstanding performance in the independent validation cohort. These findings highlight its potential as an effective non-invasive strategy for the early detection of lung cancer. Our study elucidated the regulatory relationships between epigenetic modifications and their effects on fragmentomic features. Identifying epigenetically regulated genes provided a critical foundation for developing the cfDNA fragmentomics-based machine learning model. The model demonstrated exceptional clinical performance, highlighting its substantial potential for translational application in clinical practice.

Combined analysis of the effects of hypoxia and oxidative stress on DNA methylation and the transcriptome in HTR-8/SVneo trophoblast cells.

The alterations in DNA methylation and transcriptome in trophoblast cells under conditions of low oxygen and oxidative stress have major implications for pregnancy-related disorders. However, the exact mechanism is still not fully understood. In this study, we established models of hypoxia (H group) and oxidative stress (HR group) using HTR-8/SVneo trophoblast cells and performed combined analysis of genome-wide DNA methylation changes using reduced representation bisulphite sequencing and transcriptome expression changes using RNA sequencing. Our findings revealed that the H group exhibited a higher number of differentially methylated genes and differentially expressed genes than the HR group. In the H group, only 0.90% of all differentially expressed genes displayed simultaneous changes in DNA methylation and transcriptome expression. After the threshold was expanded, this number increased to 6.29% in the HR group. Notably, both the H group and HR group exhibited concurrent alterations in DNA methylation and transcriptome expression within Axon guidance and MAPK signalling pathway. Among the top 25 differentially methylated KEGG pathways in the promoter region, 11 pathways were commonly enriched in H group and HR group, accounting for 44.00%. Among the top 25 KEGG pathways in transcriptome with significant differences between the H group and HR group, 10 pathways were consistent, accounting for 40.00%. By integrating our previous data on DNA methylation from preeclamptic placental tissues, we identified that the ANKRD37 and PFKFB3 genes may contribute to the pathogenesis of preeclampsia through DNA methylation-mediated transcriptome expression under hypoxic conditions.

Identification of cell-specific epigenetic patterns associated with chondroitin sulfate treatment response in an endemic arthritis, Kashin-Beck disease.

To assess the alterations in cell-specific DNA methylation associated with chondroitin sulphate response using peripheral blood collected from Kashin-Beck disease (KBD) patients before initiation of chondroitin sulphate treatment. Peripheral blood samples were collected from KBD patients at baseline of chondroitin sulphate treatment. Methylation profiles were generated using reduced representation bisulphite sequencing (RRBS) from peripheral blood. Differentially methylated regions (DMRs) were identified using MethylKit, while DMR-related genes were defined as those annotated to the gene body or 2.2-kilobase upstream regions of DMRs. Selected DMR-related genes were further validated by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) to assess expression levels. Tensor composition analysis was performed to identify cell-specific differential DNA methylation from bulk tissue. This study revealed 21,060 hypermethylated and 44,472 hypomethylated DMRs, and 13,194 hypermethylated and 22,448 hypomethylated CpG islands for differential global methylation for chondroitin sulphate treatment response. A total of 12,666 DMR-related genes containing DMRs were identified in their promoter regions, such as CHL1 (false discovery rate (FDR) = 2.11 × 10-11), RIC8A (FDR = 7.05 × 10-4), and SOX12 (FDR = 1.43 × 10-3). Additionally, RIC8A and CHL1 were hypermethylated in responders, while SOX12 was hypomethylated in responders, all showing decreased gene expression. The patterns of cell-specific differential global methylation associated with chondroitin sulphate response were observed. Specifically, we found that DMRs located in TESPA1 and ATP11A exhibited differential DNA methylation between responders and non-responders in granulocytes, monocytes, and B cells. Our study identified cell-specific changes in DNA methylation associated with chondroitin sulphate response in KBD patients.

One-carbon metabolism nutrients impact the interplay between DNA methylation and gene expression in liver, enhancing protein synthesis in Atlantic salmon

ABSTRACT Supplementation of one-carbon (1C) metabolism micronutrients, which include B-vitamins and methionine, is essential for the healthy growth and development of Atlantic salmon (Salmo salar). However, the recent shift towards non-fish meal diets in salmon aquaculture has led to the need for reassessments of recommended micronutrient levels. Despite the importance of 1C metabolism in growth performance and various cellular regulations, the molecular mechanisms affected by these dietary alterations are less understood. To investigate the molecular effect of 1C nutrients, we analysed gene expression and DNA methylation using two types of omics data: RNA sequencing (RNA-seq) and reduced-representation bisulphite sequencing (RRBS). We collected liver samples at the end of a feeding trial that lasted 220 days through the smoltification stage, where fish were fed three different levels of four key 1C nutrients: methionine, vitamin B6, B9, and B12. Our results indicate that the dosage of 1C nutrients significantly impacts genetic and epigenetic regulations in the liver of Atlantic salmon, particularly in biological pathways related to protein synthesis. The interplay between DNA methylation and gene expression in these pathways may play an important role in the mechanisms underlying growth performance affected by 1C metabolism.

Oxford nanopore sequencing as an alternative to reduced representation bisulphite sequencing for the identification of CpGs of interest in livestock populations

Epigenetic marks could potentially explain a portion of the phenotypic variance related to traits of interest in animal breeding models. DNA methylation is the most studied epigenetic mark, involving the addition of a methyl group to the fifth carbon of a cytosine, which transforms it into 5-methylcytosine. This mark is commonly associated with inhibiting gene expression without altering the nucleotide sequence, particularly when located within promoter regions. While bisulphite sequencing is the gold standard for detecting methylation marks, new techniques have emerged to address bisulphite's limitations. Some of these limitations include the confounding effects of bisulphite treatment on DNA methylation and C to T polymorphisms, as well as the inability to distinguish between DNA methylation and hydroxylation. In this context, we propose employing Nanopore sequencing to identify methylation marks. We conducted sequencing on six bull semen samples using both Nanopore and reduced representation bisulphite sequencing. CpGs were filtered based on two coverage thresholds (4x and 7x). Our results indicated that Nanopore methylation data exhibited a correlation exceeding 0.95 with bisulphite sequencing data. The correlation decreased to 0.67 when analysing only CpGs with methylation frequencies ranging from 0.1 to 0.9, due to the lack of enough coverage for ONT. We also examined differentially methylated cytosines identified by each sequencing method. The overlap between the two sets of DMCs and the associated genes was limited, as RRBS library preparation predominantly amplifies gene promoter regions, while ONT covers more intergenic regions. Interestingly, both methods highlighted differentially methylated genes and positions linked to quantitative trait loci of significant economic interest in animal breeding, suggesting promising applications for ONT in the field of agri-genomics.

Identifying key DNA methylation sites and their cis-methylation quantitative loci for intramuscular fatty acid traits using genome and methylome data in Yorkshire pigs

Single-nucleotide polymorphisms (SNPs) that are located near CpG sites have the potential to influence DNA methylation levels, leading to distinct patterns of DNA methylation. These methylation quantitative trait loci (meQTLs) can influence methylation levels across the entire genome, thereby impacting phenotypic traits. In this study, we utilised reduced representation bisulphite sequencing (RRBS) to analyse DNA methylation in muscle tissue samples collected from 140 Yorkshire pigs. Additionally, we assessed SNPs using the Porcine 50K chip, examined 17 fatty acid traits, and conducted epigenome-wide association studies (EWAS) and genome-wide association studies (GWAS). Our investigation revealed numerous associations between DNA methylation and fatty acid traits, including C16:0, C18:0, C20:3n-3, and C22:6n-3. Notably, a substantial portion of these associations was specifically identified by EWAS and not GWAS. Subsequently, we performed GWAS to identify meQTLs using methylation levels of significant CpGs from the EWAS as traits. The results revealed the predominance of trans-meQTLs; however, we also discovered a cis-meQTL on SSC12. Finally, we reported a crucial CpG site (SSC12:3735523) and its cis-meQTL on SSC12, which are located in close proximity to the FASN gene. Our findings contribute to the identification of CpG sites and their cis-meQTLs associated with fatty acid traits, utilising comprehensive genome and methylome data. Moreover, our study provides novel insights into the molecular mechanisms underlying fatty acid traits in pigs.

Case series on clinical applications of liquid biopsy in pediatric solid tumors: towards improved diagnostics and disease monitoring.

Solid tumors account for about 30% of all pediatric cancers. The diagnosis is typically based on histological and molecular analysis of a primary tumor biopsy. Liquid biopsies carry several advantages over conventional tissue biopsy. However, their use for genomic analysis and response monitoring of pediatric solid tumors is still in experimental stages and mostly performed retrospectively without direct impact on patient management. In this case series we discuss six clinical cases of children with a solid tumor for whom a liquid biopsy assay was performed and demonstrate the potential of liquid biopsy for future clinical decision making. We performed quantitative real-time PCR (RT-qPCR), droplet digital PCR (ddPCR) or reduced representation bisulphite sequencing of cell-free DNA (cfRRBS) on liquid biopsies collected from six pediatric patients with a solid tumor treated between 2017 and 2023 at the Princess Máxima Center for Pediatric Oncology in the Netherlands. Results were used to aid in clinical decision making by contribution to establish a diagnosis, by prognostication and response to therapy monitoring. In three patients cfRRBS helped to establish the diagnosis of a rhabdomyosarcoma, an Ewing sarcoma and a neuroblastoma (case 1-3). In two patients, liquid biopsies were used for prognostication, by MYCN ddPCR in a patient with neuroblastoma and by RT-qPCR testing rhabdomyosarcoma-specific mRNA in bone marrow of a patient with a rhabdomyosarcoma (case 4 and 5). In case 6, mRNA testing demonstrated disease progression and assisted clinical decision making. This case series illustrates the value of liquid biopsy. We further demonstrate and recommend the use of liquid biopsies to be used in conjunction with conventional methods for the determination of metastatic status, prognostication and monitoring of treatment response in patients with pediatric solid tumors.

Buffy coat signatures of breast cancer risk in a prospective cohort study

BackgroundEpigenetic alterations are a near-universal feature of human malignancy and have been detected in malignant cells as well as in easily accessible specimens such as blood and urine. These findings offer promising applications in cancer detection, subtyping, and treatment monitoring. However, much of the current evidence is based on findings in retrospective studies and may reflect epigenetic patterns that have already been influenced by the onset of the disease.MethodsStudying breast cancer, we established genome-scale DNA methylation profiles of prospectively collected buffy coat samples (n = 702) from a case–control study nested within the EPIC-Heidelberg cohort using reduced representation bisulphite sequencing (RRBS).ResultsWe observed cancer-specific DNA methylation events in buffy coat samples. Increased DNA methylation in genomic regions associated with SURF6 and REXO1/CTB31O20.3 was linked to the length of time to diagnosis in the prospectively collected buffy coat DNA from individuals who subsequently developed breast cancer. Using machine learning methods, we piloted a DNA methylation-based classifier that predicted case–control status in a held-out validation set with 76.5% accuracy, in some cases up to 15 years before clinical diagnosis of the disease.ConclusionsTaken together, our findings suggest a model of gradual accumulation of cancer-associated DNA methylation patterns in peripheral blood, which may be detected long before clinical manifestation of cancer. Such changes may provide useful markers for risk stratification and, ultimately, personalized cancer prevention.

DNA methylation and gene expression changes in mouse mammary tissue during successive lactations: part II – the impact of lactation rank

ABSTRACT Mastitis is among the main reasons women cease breastfeeding. In farm animals, mastitis results in significant economic losses and the premature culling of some animals. Nevertheless, the effect of inflammation on the mammary gland is not completely understood. This article discusses the changes to DNA methylation in mouse mammary tissue caused by lipopolysaccharide-induced inflammation after in vivo intramammary challenges and the differences in DNA methylation between 1st and 2nd lactations. Lactation rank induces 981 differential methylations of cytosines (DMCs) in mammary tissue. Inflammation in 1st lactation compared to inflammation in 2nd lactation results in the identification of 964 DMCs. When comparing inflammation in 1st vs. 2nd lactations with previous inflammation history, 2590 DMCs were identified. Moreover, Fluidigm PCR data show changes in the expression of several genes related to mammary function, epigenetic regulation, and the immune response. We show that the epigenetic regulation of two successive physiological lactations is not the same in terms of DNA methylation and that the effect of lactation rank on DNA methylation is stronger than that of the onset of inflammation. The conditions presented here show that few DMCs are shared between comparisons, suggesting a specific epigenetic response depending on lactation rank, the presence of inflammation, and even whether the cells had previously suffered inflammation. In the long term, this information could lead to a better understanding of the epigenetic regulation of lactation in both physiological and pathological conditions. Abbreviations: RRBS, reduced representation bisulphite sequencing; RT-qPCR, real-time quantitative polymerase chain reaction; MEC, mammary epithelial cells; MaSC, mammary stem cell; TSS, transcription start site; TTS, transcription termination site; UTR, untranslated region; SINE, short interspersed nuclear element; LINE, long interspersed nuclear element; CGI, CpG island; DEG, differentially expressed gene; DMC, differentially methylated cytosine; DMR, differentially methylated region; GO term, gene ontology term; MF, molecular function; BP, biological process

Region-based epigenetic clock design improves RRBS-based age prediction.

Recent studies suggest that epigenetic rejuvenation can be achieved using drugs that mimic calorie restriction and techniques such as reprogramming-induced rejuvenation. To effectively test rejuvenation in vivo, mouse models are the safest alternative. However, we have found that the recent epigenetic clocks developed for mouse reduced-representation bisulphite sequencing (RRBS) data have significantly poor performance when applied to external datasets. We show that the sites captured and the coverage of key CpGs required for age prediction vary greatly between datasets, which likely contributes to the lack of transferability in RRBS clocks. To mitigate these coverage issues in RRBS-based age prediction, we present two novel design strategies that use average methylation over large regions rather than individual CpGs, whereby regions are defined by sliding windows (e.g. 5 kb), or density-based clustering of CpGs. We observe improved correlation and error in our regional blood clocks (RegBCs) compared to published individual-CpG-based techniques when applied to external datasets. The RegBCs are also more robust when applied to low coverage data and detect a negative age acceleration in mice undergoing calorie restriction. Our RegBCs offer a proof of principle that age prediction of RRBS datasets can be improved by accounting for multiple CpGs over a region, which negates the lack of read depth currently hindering individual-CpG-based approaches.

Coordinated alternation of DNA methylation and alternative splicing of PBRM1 affect bovine sperm structure and motility

ABSTRACT DNA methylation and gene alternative splicing drive spermatogenesis. In screening DNA methylation markers and transcripts related to sperm motility, semen from three pairs of full-sibling Holstein bulls with high and low motility was subjected to reduced representation bisulphite sequencing. A total of 948 DMRs were found in 874 genes (gDMRs). Approximately 89% of gDMR-related genes harboured alternative splicing events, including SMAD2, KIF17, and PBRM1. One DMR in exon 29 of PBRM1 with the highest 5mC ratio was found, and hypermethylation in this region was related to bull sperm motility. Furthermore, alternative splicing events at exon 29 of PBRM1 were found in bull testis, including PBRM1-complete, PBRM1-SV1 (exon 28 deletion), and PBRM1-SV2 (exons 28–29 deletion). PBRM1-SV2 exhibited significantly higher expression in adult bull testes than in newborn bull testes. In addition, PBRM1 was localized to the redundant nuclear membrane of bull sperm, which might be related to sperm motility caused by sperm tail breakage. Therefore, the hypermethylation of exon 29 may be associated with the production of PBRM1-SV2 in spermatogenesis. These findings indicated that DNA methylation alteration at specific loci could regulate gene splicing and expression and synergistically alter sperm structure and motility.

Environmental and genealogical effects on DNA methylation in a widespread apomictic dandelion lineage.

DNA methylation in plant genomes occurs in different sequences and genomic contexts that have very different properties. DNA methylation that occurs in CG (mCG) sequence context shows transgenerational stability and high epimutation rate, and can thus provide genealogical information at short time scales. However, due to meta-stability and because mCG variants may arise due to other factors than epimutation, such as environmental stress exposure, it is not clear how well mCG captures genealogical information at micro-evolutionary time scales. Here, we analysed DNA methylation variation between accessions from a geographically widespread, apomictic common dandelion (Taraxacum officinale) lineage when grown experimentally under different light conditions. Using a reduced-representation bisulphite sequencing approach, we show that the light treatment induced differentially methylated cytosines (DMCs) in all sequence contexts, with a bias towards transposable elements. Accession differences were associated mainly with DMCs in CG context. Hierarchical clustering of samples based on total mCG profiles revealed a perfect clustering of samples by accession identity, irrespective of light conditions. Using microsatellite information as a benchmark of genetic divergence within the clonal lineage, we show that genetic divergence between accessions correlates strongly with overall mCG profiles. However, our results suggest that environmental effects that do occur in CG context may produce a heritable signal that partly dilutes the genealogical signal. Our study shows that methylation information in plants can be used to reconstruct micro-evolutionary genealogy, providing a useful tool in systems that lack genetic variation such as clonal and vegetatively propagated plants.

Loss of the crumbs cell polarity complex disrupts epigenetic transcriptional control and cell cycle progression in the developing retina.

The crumbs cell polarity complex plays a crucial role in apical-basal epithelial polarity, cellular adhesion, and morphogenesis. Homozygous variants in human CRB1 result in autosomal recessive Leber congenital amaurosis (LCA) and retinitis pigmentosa (RP), with no established genotype-phenotype correlation. The associated protein complexes have key functions in developmental pathways; however, the underlying disease mechanism remains unclear. Using the oko meduzym289/m289 (crb2a-/- ) zebrafish, we performed integrative transcriptomic (RNA-seq data) and methylomic [reduced representation bisulphite sequencing (RRBS)] analysis of whole retina to identify dysregulated genes and pathways. Delayed retinal cell specification was identified in both the crb2a-/- zebrafish and CRB1 patient-derived retinal organoids, highlighting the dysfunction of cell cycle modulation and epigenetic transcriptional control. Differential DNA methylation analysis revealed novel hypermethylated pathways involving biological adhesion, Hippo, and transforming growth factor β (TGFβ) signalling. By integrating gene expression with DNA methylation using functional epigenetic modules (FEM), we identified six key modules involving cell cycle control and disturbance of TGFβ, bone morphogenetic protein (BMP), Hippo, and SMAD protein signal transduction pathways, revealing significant interactome hotspots relevant to crb2a function and confirming the epigenetic control of gene regulation in early retinal development, which points to a novel mechanism underlying CRB1-retinopathies. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

DNA methylation changes from primary cultures through senescence-bypass in Syrian hamster fetal cells initially exposed to benzo[a]pyrene

Current chemical testing strategies are limited in their ability to detect non-genotoxic carcinogens (NGTxC). Epigenetic anomalies develop during carcinogenesis regardless of whether the molecular initiating event is associated with genotoxic (GTxC) or NGTxC events; therefore, epigenetic markers may be harnessed to develop new approach methodologies that improve the detection of both types of carcinogens. This study used Syrian hamster fetal cells to establish the chronology of carcinogen-induced DNA methylation changes from primary cells until senescence-bypass as an essential carcinogenic step. Cells exposed to solvent control for 7 days were compared to naïve primary cultures, to cells exposed for 7 days to benzo[a]pyrene, and to cells at the subsequent transformation stages: normal colonies, morphologically transformed colonies, senescence, senescence-bypass, and sustained proliferation in vitro. DNA methylation changes identified by reduced representation bisulphite sequencing were minimal at day-7. Profound DNA methylation changes arose during cellular senescence and some of these early differentially methylated regions (DMRs) were preserved through the final sustained proliferation stage. A set of these DMRs (e.g., Pou4f1, Aifm3, B3galnt2, Bhlhe22, Gja8, Klf17, and L1l) were validated by pyrosequencing and their reproducibility was confirmed across multiple clones obtained from a different laboratory. These DNA methylation changes could serve as biomarkers to enhance objectivity and mechanistic understanding of cell transformation and could be used to predict senescence-bypass and chemical carcinogenicity.

Comparison between indicine and taurine cattle DNA methylation reveals epigenetic variation associated to differences in morphological adaptive traits

ABSTRACT Indicine and taurine subspecies present distinct morphological traits as a consequence of environmental adaptation and artificial selection. Although the two subspecies have been characterized and compared at genome-wide level and at specific loci, their epigenetic diversity has not yet been explored. In this work, Reduced Representation Bisulphite Sequencing (RRBS) profiling of the taurine Angus (A) and indicine Nellore (N) cattle breeds was applied to identify methylation differences between the two subspecies. Genotyping by sequencing (GBS) of the same animals was performed to detect single nucleotide polymorphisms (SNPs) at cytosines in CpG dinucleotides and remove them from the differential methylation analysis. A total of 660,845 methylated cytosines were identified within the CpG context (CpGs) across the 10 animals sequenced (5 N and 5 A). A total of 25,765 of these were differentially methylated (DMCs). Most DMCs clustered in CpG stretches nearby genes involved in cellular and anatomical structure morphogenesis. Also, sequences flanking DMC were enriched in SNPs compared to all other CpGs, either methylated or unmethylated in the two subspecies. Our data suggest a contribution of epigenetics to the regulation and divergence of anatomical morphogenesis in the two subspecies relevant for cattle evolution and sub-species differentiation and adaptation.

Cell-Free DNA as a Diagnostic and Prognostic Biomarker in Pediatric Rhabdomyosarcoma.

Total cell-free DNA (cfDNA) and tumor-derived cfDNA (ctDNA) can be used to study tumor-derived genetic aberrations. We analyzed the diagnostic and prognostic potential of cfDNA and ctDNA, obtained from pediatric patients with rhabdomyosarcoma. cfDNA was isolated from diagnostic plasma samples from 57 patients enrolled in the EpSSG RMS2005 study. To study the diagnostic potential, shallow whole genome sequencing (shWGS) and cell-free reduced representation bisulphite sequencing (cfRRBS) were performed in a subset of samples and all samples were tested using droplet digital polymerase chain reaction to detect methylated RASSF1A (RASSF1A-M). Correlation with outcome was studied by combining cfDNA RASSF1A-M detection with analysis of our rhabdomyosarcoma-specific RNA panel in paired cellular blood and bone marrow fractions and survival analysis in 56 patients. At diagnosis, ctDNA was detected in 16 of 30 and 24 of 26 patients using shallow whole genome sequencing and cfRRBS, respectively. Furthermore, 21 of 25 samples were correctly classified as embryonal by cfRRBS. RASSF1A-M was detected in 21 of 57 patients. The presence of RASSF1A-M was significantly correlated with poor outcome (the 5-year event-free survival [EFS] rate was 46.2% for 21 RASSF1A-M‒positive patients, compared with 84.9% for 36 RASSF1A-M‒negative patients [P < .001]). RASSF1A-M positivity had the highest prognostic effect among patients with metastatic disease. Patients both negative for RASSF1A-M and the rhabdomyosarcoma-specific RNA panel (28 of 56 patients) had excellent outcome (5-year EFS 92.9%), while double-positive patients (11/56) had poor outcome (5-year EFS 13.6%, P < .001). Analyzing ctDNA at diagnosis using various techniques is feasible in pediatric rhabdomyosarcoma and has potential for clinical use. Measuring RASSF1A-M in plasma at initial diagnosis correlated significantly with outcome, particularly when combined with paired analysis of blood and bone marrow using a rhabdomyosarcoma-specific RNA panel.

DNA methylation in clonal duckweed (Lemna minor L.) lineages reflects current and historical environmental exposures.

Environmentally induced DNA methylation variants may mediate gene expression responses to environmental changes. If such induced variants are transgenerationally stable, there is potential for expression responses to persist over multiple generations. Our current knowledge in plants, however, is almost exclusively based on studies conducted in sexually reproducing species where the majority of DNA methylation changes are subject to resetting in germlines, limiting the potential for transgenerational epigenetics stress memory. Asexual reproduction circumvents germlines, and may therefore be more conducive to long-term inheritance of epigenetic marks. Taking advantage of the rapid clonal reproduction of the common duckweed Lemna minor, we hypothesize that long-term, transgenerational stress memory from exposure to high temperature can be detected in DNA methylation profiles. Using a reduced representation bisulphite sequencing approach (epiGBS), we show that temperature stress induces DNA hypermethylation at many CG and CHG cytosine contexts but not CHH. Additionally, differential methylation in CHG context that was observed was still detected in a subset of cytosines, even after 3-12 generations of culturing in a common environment. This demonstrates a memory effect of stress reflected in the methylome and that persists over multiple clonal generations. Structural annotation revealed that this memory effect in CHG methylation was enriched in transposable elements. The observed epigenetic stress memory is probably caused by stable transgenerational persistence of temperature-induced DNA methylation variants across clonal generations. To the extent that such epigenetic memory has functional consequences for gene expression and phenotypes, this result suggests potential for long-term modulation of stress responses in asexual plants.

Improving Infinium MethylationEPIC data processing: re-annotation of enhancers and long noncoding RNA genes and benchmarking of normalization methods

ABSTRACT Illumina Infinium DNA Methylation (5mC) arrays are a popular technology for low-cost, high-throughput, genome-scale measurement of 5mC distribution, especially in cancer and other complex diseases. After the success of its HumanMethylation450 array (450k), Illumina released the MethylationEPIC array (850k) featuring increased coverage of enhancers. Despite the widespread use of 850k, analysis of the corresponding data remains suboptimal: it still relies mostly on Illumina’s default annotation, which underestimates enhancerss and long noncoding RNAs. Results: We have thus developed an approach, based on the ENCODE and LNCipedia databases, which greatly improves upon Illumina’s default annotation of enhancers and long noncoding transcripts. We compared the re-annotated 850k with both 450k and reduced-representation bisulphite sequencing (RRBS), another high-throughput 5mC profiling technology. We found 850k to cover at least three times as many enhancers and long noncoding RNAs as either 450k or RRBS. We further investigated the reproducibility of the three technologies, applying various normalization methods to the 850k data. Most of these methods reduced variability to a level below that of RRBS data. We then used 850k with our new annotation and normalization to profile 5mC changes in breast cancer biopsies. 850k highlighted aberrant enhancer methylation as the predominant feature, in agreement with previous reports. Our study provides an updated processing approach for 850k data, based on refined probe annotation and normalization, allowing for improved analysis of methylation at enhancers and long noncoding RNA genes. Our findings will help to further advance understanding of the DNA methylome in health and disease.

The impact of Piscirickettsia salmonis infection on genome-wide DNA methylation profile in Atlantic Salmon

Salmon rickettsial septicaemia (SRS), caused by the bacteria Piscirickettsia salmonis (P. salmonis), is responsible for significant mortality in farmed Atlantic salmon in Chile. Currently there are no effective treatments or preventive measures for this disease, although genetic selection or genome engineering to increase salmon resistance to SRS are promising strategies. The accuracy and efficiency of these strategies are usually influenced by the available biological background knowledge of the disease. The aim of this study was to investigate DNA methylation changes in response to P. salmonis infection in the head kidney and liver tissue of Atlantic salmon, and the interaction between gene expression and DNA methylation in the same tissues. The head kidney and liver methylomes of 66 juvenile salmon were profiled using reduced representation bisulphite sequencing (RRBS), and compared between P. salmonis infected animals (3 and 9 days post infection) and uninfected controls, and between SRS resistant and susceptible fish. Methylation was correlated with matching RNA-Seq data from the same animals, revealing that methylation in the first exon leads to an important repression of gene expression. Head kidney methylation showed a clear response to the infection, associated with immunological processes such as actin cytoskeleton regulation, phagocytosis, endocytosis and pathogen associated pattern receptor signaling. Our results contribute to the growing understanding of the role of methylation in regulation of gene expression and response to infectious diseases and could inform the incorporation of epigenetic markers into genomic selection for disease resistant and the design of diagnostic epigenetic markers to better manage fish health in salmon aquaculture.

Comparative analysis of the Illumina Mouse Methylation BeadChip and Reduced Representation Bisulphite Sequencing for routine DNA methylation analysis of murine samples

Comparative analysis of the Illumina Mouse Methylation BeadChip and Reduced Representation Bisulphite Sequencing for routine DNA methylation analysis of murine samples