Bisulfite Genomic Sequencing Research Articles

Abstract B005: Investigating distinct methylation signatures characteristic of breast cancer subtypes in residual disease via cell free DNA methylation

Abstract Breast cancer is the most common type of cancer in females and recurrence increases over time, unlike many other cancers. Treatment for recurrent cancer is often the same as the primary with no additional biopsies taken. Current research suggests that subtype switching and tumor character changes frequently occur between primary and recurrent breast cancers. Therefore, it may be beneficial to patients to switch treatment based on these changes. Since physical biopsies are cumbersome and not always feasible, liquid biopsies open a way to monitor tumor changes less invasively and more comprehensively. It has been previously established that cell free DNA (cfDNA) shed into the bloodstream from dying cells can reflect cell type of origin via methylation pattern and rate of death via concentration of cfDNA. In this pilot study, we seek to look at the cfDNA of fifteen late stage pre-/ post-surgical breast cancer patients who also received radiation treatment. We will perform cfDNA extraction on the serum of these patients and both whole genome bisulfite sequencing (WGBS), which chemically converts unmethylated cytosine to uracil/thymine in the DNA and is the current gold-standard of methylation sequencing, and a newer method, enzymatic methylation sequencing (EM-seq), which enzymatically converts (TET2/APOBEC) unmethylated cytosine to uracil/thymine and potentially preserves more of the cfDNA. To validate any signatures found in the cfDNA of the breast cancer patients, we have begun the genomic DNA (gDNA) extraction and WGBS/ EM-seq protocols on a variety of breast cancer cell lines including: MCF10A, MCFDCIS, MCF7, T47D, BT474 MDA MB453, MDA MB436 and MDA MB231 (including in-lab brain, bone, and lung metastatic clones). Bioanalyzer traces are produced from the extracted cfDNA/gDNA and also for the final sequencing libraries. The success of the methylation conversion is evaluated after sequencing data is returned and conversion rates of the cytosine to uracil/thymine are compared to unmethylated DNA control (lambda) and methylated DNA control (pUC19). Once this sequencing data is obtained, we use an in-lab deconvolution algorithm to detect cell types of origin and intend to make the algorithm more robust for cancer cell types as well. We have currently produced breast cancer cell line methylation sequencing libraries and are in the process of producing the libraries for the patient samples. Our current data suggests that there are changes in the cfDNA general fragmentation patterns and cfDNA concentrations between pre-/post- surgery samples. Once our sequencing data is obtained for the patient samples, we will run our deconvolution algorithm. The potential of characterizing breast cancer subtype and progression signatures in cfDNA of late stage pre-/ post-surgical breast cancer patients can have significant impact on patient treatment options. Identifying these breast cancer signatures less invasively and, therefore, more frequently may allow for early and more targeted intervention to improve breast cancer patient outcomes. Citation Format: Amber Alley, Megan McNamara, Sidharth Jain, Anton Wellstein. Investigating distinct methylation signatures characteristic of breast cancer subtypes in residual disease via cell free DNA methylation [abstract]. In: Proceedings of the AACR Special Conference: Liquid Biopsy: From Discovery to Clinical Implementation; 2024 Nov 13-16; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(21_Suppl):Abstract nr B005.

Abstract PR011: Advance prostate cancer detection through epigenomic profiling of cell-free DNA

Abstract Introduction: Metastatic castration-resistant prostate cancer (mCRPC) is a heterogeneous disease which can be classified into clinically relevant subtypes based on the expression of genes, such as the androgen receptor (AR) and neuroendocrine markers. Neuroendocrine prostate cancer (NEPC), characterized by gain of stem-like and neuroendocrine features and lack of AR expression is a clinically aggressive variant. Due to the lack of adequate biomarkers, NEPC is usually detected at a very advanced stage. There is mounting evidence that molecular subtype changes seen in NEPC are enforced by widespread epigenetic alterations, in particular DNA methylation changes. In this study, we aim to devise a novel DNA methylation-based assay for molecular subtyping and disease monitoring from cell-free DNA (cfDNA). Methods: We analyzed genome wide methylation patterns in 56 prostate cancer patient-derived xenograft (PDX) and 128 mCRPC tumors using array- and sequencing-based assays. We integrated DNA methylation at promoters, gene bodies and transcription factor binding site (TFBS) to determine the landscape of methylation alterations at key lineage specific genes. Using whole genome methylation derived from tissue with matched expression data we developed a deep learning framework to predict gene expression directly from tissue or cfDNA. Using key marker genes, the model was used to discern tumor molecular phenotypes from tissue and cfDNA in three independent cohorts of mCRPC patients using whole genome bisulfite sequencing and low-pass Enzymatic Methyl-Seq (EM-seq). Results: We observed a tight association between promoter, gene body and TFBS methylation with gene expression. Inferring gene expression from methylation for lineage specific markers such as AR, KLK3, ASCL1, INSM1, SRRM4 and DLL3 we classified molecular subtypes from both tissue and cfDNA. Additionally, for AR and ASCL1, we identified core sets of TFBSs whose differential methylation allowed for accurate assay-independent molecular subtype quantification. Applying the optimized quantitative model to mCRPC patients who underwent comprehensive tissue sampling by rapid autopsy we observed accurate subtype classification from both tissue samples and cfDNA for all cases. A similar analytical performance was observed in additional clinical mCRPC cohorts with cfDNA. Conclusion: Whole-genome methylation analysis of cfDNA allows for the prediction of gene expression patterns in tumor tissues, enabling non-invasive tumor subclassification and assessment of therapeutic targets. Citation Format: Mohamed Adil, Brian Hanratty, Pallabi Mustafi, Chitvan Mittal, Helen Richards, Ilsa Coleman, Radhika Patel, Anna-Lisa Doebley, Robert Patton, Eden Cruikshank, Patricia Galipeau, Ruth Dumpit, Martine Roudier, Jin-Yih Low, Navonil Sarkar, Robert Montgomery, Eva Corey, Colm Morrissey, Peter Nelson, Gavin Ha, Michael Haffner. Advance prostate cancer detection through epigenomic profiling of cell-free DNA [abstract]. In: Proceedings of the AACR Special Conference: Liquid Biopsy: From Discovery to Clinical Implementation; 2024 Nov 13-16; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(21_Suppl):Abstract nr PR011.

Sex and tissue-specificity of piRNA regulation in adult mice following perinatal lead (Pb) exposure

ABSTRACT Lead (Pb) is a neurotoxicant with early life exposure linked to long-term health effects. Piwi-interacting RNAs (piRNAs) are small non-coding RNAs that associate with PIWIL proteins to induce DNA methylation. It remains unknown whether Pb exposure influences piRNA expression. This study evaluated how perinatal Pb exposure (32 ppm in drinking water) impacts piRNA expression in adult mice and assessed piRNA dysregulation as a potential mechanism for Pb-induced toxicity. Pb exposure effects on piRNA expression and associated gene repression in the germline (testis/ovary) and soma (liver and brain) were evaluated. Small RNA sequencing was used to determine differentially expressed piRNAs, RT-qPCR to examine piRNA target expression, and whole genome bisulfite sequencing to evaluate target DNA methylation status. Three piRNAs (mmpiR-1500602, mmpiR-0201406, and mmpiR-0200026) were significant after multiple testing correction (all downregulated in the male Pb-exposed brain in comparison to control; FDR < 0.05). Within piOxiDB, TAO Kinase 3 was identified as a downstream mRNA target for one of the three Pb-sensitive piRNA. The Pb-exposed male brain exhibited increased Taok3 expression (p < 0.05) and decreased DNA methylation (FDR < 0.01). The results demonstrate that perinatal Pb exposure stably influences longitudinal piRNA expression in a tissue- and sex-specific manner, potentially via DNA methylation-directed mechanisms.

PATH-24. NOVEL MOLECULAR SUBTYPES IDENTIFIED IN INTRACRANIAL GERM CELL TUMORS SHOWED UNIQUE WHOLE-GENOME METHYLATION PROFILE

Abstract BACKGROUND By integrating whole-genome sequencing (WGS) and RNA-seq data, we found that intracranial germ cell tumors (IGCTs) can be divided into three distinct molecular subtypes: immune-hot, MYC/E2F, and SHH, each characterized by a unique profile of genomic abnormalities. To further define these subtypes, we subsequently performed whole-genome bisulfate sequencing (WGBS) on IGCT tumors METHODS Within our prior discovery cohort (n=110), DNA from 84 samples was available for profiling methylation levels using Whole Genome Bisulfite Sequencing (WGBS). Fastq data was aligned to the human reference genome (hg19), with duplicates removed, and DNA methylation was extracted using bismark pipeline (ver0.24.0). After quality control, weight methylation level was calculated for each CpG island as the sum of methylated reads divided by coverage. To identify IGCTs molecular groups, we applied hierarchical consensus k-means clustering to the most variant methylated CpG island (n=1000) across 84 IGCTs tumor tissues. The methylation level of RNA-seq subtypes was visualized to further clarify the molecular subtypes of IGCTs RESULTS Unsupervised clustering revealed that IGCTs can be divided into three groups based on the CpG island methylation levels: hypomethylation group (n=32), moderate methylation group (n=27), and hyper methylation group (n=25). The hypomethylation and moderate methylation group are predominantly observed in germinoma (&amp;gt;90%). In contrast, hypermethylation group comprises the majority of non-germinomatous germ cell tumors. Furthermore, these subgroups derived from methylation analysis are highly consistent with the molecular subtypes previously defined through gene transcription profile analysis. We observed that the majority of MYC/E2F subtype (24/33) belonged to the hypomethylation group, whereas the majority of immune-hot subtype (14/23) fell into the moderate methylation group. The SHH subtype was the most common within the hypermethylation group (18/26) CONCLUSIONS The three molecular IGCT subtypes showed unique methylation profile. These results indicate different mechanisms of tumorigenesis among the subtypes and pave the way for the clinical application of methylation-based subtyping.

Integrative analysis of whole genome bisulfite and transcriptome sequencing reveals the effect of sodium butyrate on DNA methylation in the differentiation of bovine skeletal muscle satellite cells

Butyric acid as a short-chain fatty acid (SCFA) is one of the key microbial metabolites of ruminants. Numerous studies indicate that butyrate is crucial in muscle growth and development, and plays an important molecular regulatory role mainly by inhibiting histone deacetylation. DNA methylation, a major epigenetic modification, is involved in cell differentiation. Butyrate, in addition to its role in acetylation modifications, can alter the DNA methylation status of cells. However, the impact of butyrate on the DNA methylation of bovine skeletal muscle satellite cells (SMSCs) remains unclear. In this study, we developed a differentiation model of SMSCs and employed RNA sequencing (RNA-seq) alongside whole genome bisulfite sequencing (WGBS) to explore the effects of butyrate treatment on DNA methylation status and its relationship with gene expression. Treatment of SMSCs with sodium butyrate (NaB) at 1.0 mM for 2 days significantly inhibited the expression of DNA methyltransferases (DNMT1, DNMT2, DNMT3A) at the mRNA and protein levels while promoting the expression of demethylases (TET1, TET2, TET3) at mRNA levels. WGBS identified 4292 differentially methylated regions (DMRs), comprising 2294 hypermethylated and 1998 hypomethylated regions. These DMRs were significantly enriched in the MAPK, cAMP, and Wnt signaling pathways, all of which are implicated in myogenesis and development. Combining RNA-seq and WGBS data revealed a total of 130 overlapping genes, including MDFIC, CREBBP, DMD, LTBP2 and KLF4. These genes are predominantly involved in regulating the FoxO, MAPK, PI3K-Akt, and Wnt signaling pathways. This study provides new insights into the effects of butyrate-mediated DNA methylation on SMSC development and enhances our understanding of butyrate as an epigenetic modifier beyond its role in acetylation.

Methylome of circulating cell free DNA as a novel biomarker tool: from acute coronary syndrome to broader risk stratification of cardiovascular disease

Abstract Background Cell death during tissue injury leads to release of chromosomal DNA which is promptly degraded. However, short DNA fragments wrapped around nucleosomes survive longer in the bloodstream and can be isolated as circulating cell-free DNA (ccfDNA). Oncology field pioneered the use of cfDNA as clinical biomarker by identifying sequence modifications induced in cancer cells. Importantly, ccfDNA also retains the same epigenetic information as the tissue of origin. Despite major advances for disease diagnosis and therapy, coronary artery diseases (CAD) and associated conditions such acute myocardial infarction (AMI) or heart failure (HF) are responsible for leading number of deaths worldwide. Preventive and surveillance strategies with non-invasive testing are needed. Purpose Evaluation of ccfDNA methylation profiles as prognostic non-invasive biomarker for risk stratification of acute coronary syndromes and their potential application in preventive and surveillance strategies in other cardiovascular disorders. Methods CcfDNA isolated from citrate-plasma collected from our discovery cohort, comprising healthy controls, STEMI, NSTEMI and unstable angina (UA) patients, was analyzed applying Whole Genome Bisulfite Sequencing (WGBS) using a low-input BS-seq (PBAT) protocol. Computational analysis generated differentially methylated regions (DMRs) by tiling CpG areas and setting a threshold of minimal 25% methylation difference compared to the healthy group. To assess the validity of our findings, a selection of disease-specific DMRs were tested using targeted methylation sequencing protocol in a new set of patients. This procedure consisted in enzymatic conversion of ccfDNA methylated sites followed by probe-specific enrichment libraries and sequencing at a Novaseq 6000 platform. Results The study identified a total of 1941 DMRs and from those 486 were STEMI, 223 NSTEMI and 684 UA specific. From the identified DMRs, 74.13% from STEMI, 72.68% from NSTEMI and 61.62% from UA annotated with heart-related diseases using DisGeNET gene-disease associations. Interestingly, most of identified DMRs located within intronic regions (~60%) and also annotated as enhancers and lncRNAs. The independent validation with targeted methylation sequencing identified 566, 306 and 248 DMRs, respectively, and reduced the number of disease-specific DMRs to 171 in STEMI, 115 in NSTEMI and 82 in UA. Conclusion and Future Approaches Methylation profiles of ccfDNA showed potential as a biomarker to stratify the severity of acute coronary syndromes, a novel non-invasive detection methodology that could be extrapolated to risk stratification in other cardiovascular disease such as heart failure or coronary artery disease.

PCBS: an R package for fast and accurate analysis of bisulfite sequencing data.

Whole-genome bisulfite sequencing is a powerful tool for analyzing chromatin methylation genome-wide, but analysis of whole-genome bisulfite data is slow, inflexible, and often inaccurate. We developed PCBS (Principal Component BiSulfite), a computationally efficient R package for Whole Genome Bisulfite Sequencing analysis that demonstrates remarkable accuracy and flexibility compared to current tools. PCBS identifies differentially methylated loci, differentially methylated regions, and offers novel functionality that allows for more targeted methylation analyses. PCBS uses minimal computational resources; a complete pipeline in mouse can run on a local RStudio instance in a matter of minutes. PCBS is an R package available under a GNU GPLv3 license on GitHub: https://github.com/katlande/PCBS and CRAN: https://CRAN.R-project.org/package=PCBS. Instructions for use are available at: https://katlande.github.io/PCBS/.

PRKACB is a novel imprinted gene in marsupials

BackgroundGenomic imprinting results in parent-of-origin-specific gene expression and, among vertebrates, is found only in therian mammals: marsupials and eutherians. A differentially methylated region (DMR), in which the methylation status of CpG dinucleotides differs between the two alleles, can mark the parental identity of imprinted genes. We developed a computational pipeline that detected CpG islands (CGIs) marked by both methylated and unmethylated signals in whole genome bisulfite sequencing data. This approach identified candidate marsupial DMRs in a publicly available koala methylome. One of these candidate DMRs was associated with PRKACB, a gene encoding the protein kinase A catalytic subunit beta. Nothing is known about the imprinting status of PRKACB in eutherian mammals although mutations of this gene are associated with endocrine neoplasia and other developmental disorders.ResultsIn the tammar wallaby and brushtail possum there was parent-of-origin-specific DNA methylation in the PRKACB DMR in which the maternal allele was methylated and the paternal allele was unmethylated. There were multiple RNAs transcribed from this locus. Allele-specific expression analysis identified paternal expression of a PRKACB lncRNA and an mRNA isoform. Comparison of the PRKACB gene start site between marsupials and eutherians demonstrated that the CGI is longer in marsupials. The PRKACB gene product functions in the same signalling pathway as the guanine nucleotide-binding protein alpha subunit encoded at the GNAS locus, a known eutherian imprinted gene. In a mouse methylome Gnas had three differentially methylated CGIs, while in the koala methylome the GNAS locus had two unmethylated CGIs.ConclusionsWe conclude that PRKACB is a novel, DMR-associated marsupial imprinted gene. Imprinting of PRKACB in marsupials and GNAS in eutherians may indicate a conserved selection pressure for imprinting of the protein kinase A signalling pathway in therians with the two lineages adapting by imprinting different genes.

Genome-wide DNA methylation sequencing reveals the involvement of ferroptosis in hepatotoxicity induced by dietary exposure to food-grade titanium dioxide

BackgroundFollowing the announcement by the European Food Safety Authority that the food additive titanium dioxide (E 171) is unsafe for human consumption, and the subsequent ban by the European Commission, concerns have intensified over the potential risks E 171 poses to human vital organs. The liver is the main organ for food-grade nanoparticle metabolism. It is increasingly being found that epigenetic changes may play an important role in nanomaterial-induced hepatotoxicity. However, the profound effects of E 171 on the liver, especially at the epigenetic level, remain largely unknown.MethodsMice were exposed orally to human-relevant doses of two types of E 171 mixed in diet for 28 and/or 84 days. Conventional toxicology and global DNA methylation analyses were performed to assess E 171-induced hepatotoxicity and epigenetic changes. Whole genome bisulfite sequencing and further ferroptosis protein detection were used to reveal E 171-induced changes in liver methylation profiles and toxic mechanisms.ResultsExposed to E 171 for 28 and/or 84 days resulted in reduced global DNA methylation and hydroxymethylation in the liver of mice. E 171 exposure for 84 days elicited inflammation and damage in the mouse liver, whereas 28-day exposure did not. Whole-genome DNA methylation sequencing disclosed substantial methylation alterations at the CG and non-CG sites of the liver DNA in mice exposed to E 171 for 84 days. Mechanistic analysis of the DNA methylation alterations indicated that ferroptosis contributed to the liver toxicity induced by E 171. E 171-induced DNA methylation changes triggered NCOA4-mediated ferritinophagy, attenuated the protein levels of GPX4, FTH1, and FTL in the liver, and thereby caused ferroptosis.ConclusionsLong-term oral exposure to E 171 triggers hepatotoxicity and induces methylation changes in both CG and non-CG sites of liver DNA. These epigenetic alterations activate ferroptosis in the liver through NCOA4-mediated ferritinophagy, highlighting the role of DNA methylation and ferroptosis in the potential toxicity caused by E 171 in vivo.

DNA Methylation Carries Signatures of Sublethal Effects Under Thermal Stress in Loggerhead Sea Turtles.

To date, studies of the impacts of climate warming on individuals and populations have mostly focused on mortality and thermal tolerance. In contrast, much less is known about the consequences of sublethal effects, which are more challenging to detect, particularly in wild species with cryptic life histories. This necessitates the development of molecular tools to identify their signatures. In a split-clutch field experiment, we relocated clutches of wild, nesting loggerhead sea turtles (Caretta caretta) to an insitu hatchery. Eggs were then split into two sub-clutches and incubated under shallow or deep conditions, with those in the shallow treatment experiencing significantly higher temperatures in otherwise natural conditions. Although no difference in hatching success was observed between treatments, hatchlings from the shallow, warmer treatment had different length-mass relationships and were weaker at locomotion tests than their siblings incubated in the deep, cooler treatment. To characterise the molecular signatures of these thermal effects, we performed whole genome bisulfite sequencing on blood samples collected upon emergence. We identified 287 differentially methylated sites between hatchlings from different treatments, including on genes with neurodevelopmental, cytoskeletal, and lipid metabolism functions. Taken together, our results show that higher incubation temperatures induce sublethal effects in hatchlings, which are reflected in their DNA methylation status at identified sites. These sites could be used as biomarkers of thermal stress, especially if they are retained across life stages. Overall, this study suggests that global warming reduces hatchling fitness, which has implications for dispersal capacity and ultimately a population's adaptive potential. Conservation efforts for these endangered species and similar climate-threatened taxa will therefore benefit from strategies for monitoring and mitigating exposure to temperatures that induce sublethal effects.

Role of Epigenetic Factors in Determining the Biological Behavior and Prognosis of Hepatocellular Carcinoma.

The current study's objective is to evaluate the molecular genetic mechanisms influencing the biological behavior of hepatocellular carcinoma (HCC) by analyzing the transcriptomic and epigenetic signatures of the tumors. Transcriptomic data were downloaded from the NCBI GEO database. We investigated the expression differences between the GSE46444 (48 cirrhotic tissues versus 88 HCC tissues) and GSE63898 (168 cirrhotic tissues versus 228 HCC tissues) data sets using GEO2R. Differentially expressed genes were evaluated using GO and KEGG metabolic pathway analysis websites. Whole genome bisulfite sequencing (WGBS) and Methylated DNA Immunoprecipitation Sequencing (MeDIP-Seq) data sets (26 HCC tissues versus 26 adjacent non-tumoral tissues) were also downloaded from the NCBI SRA database. These data sets were analyzed using Bismark and QSEA, respectively. The methylation differences between the groups were assessed using functional enrichment analysis. In the GSE46444 data set, 80 genes were upregulated, and 315 genes were downregulated in the tumor tissue (HCC tissue) compared to the non-tumor cirrhotic tissue. In the GSE63898 data set, 1261 genes were upregulated, and 458 genes were downregulated in the cirrhotic tissue compared to the tumor tissues. WGBS revealed that 20 protein-coding loci were hypermethylated. while the hypomethylated regions were non-protein-coding. The methylated residues of the tumor tissue, non-tumorous cirrhotic tissue, and healthy tissue were comparable. MeDIP-Seq, conducted on tumoral and non-tumoral tissues, identified hypermethylated or hypomethylated areas as protein-coding regions. The functional enrichment analysis indicated that these genes were related to pathways including peroxisome, focal adhesion, mTOR, RAP1, Phospholipase D, Ras, and PI3K/AKT signal transduction. The investigation of transcriptomic and epigenetic mechanisms identified several genes significant in the biological behavior of HCC. These genes present potential targets for the development of targeted therapy.

Comprehensive Analysis of Methylome and Transcriptome to Identify Potential Genes Regulating Porcine Testis Development.

DNA methylation plays a critical role in regulating gene expression during testicular development. However, few studies report on candidate genes related to the DNA methylation regulation of porcine testicular development. This study examined the differentially expressed genes (DEGs) and their methylation levels in testicular tissues from pigs at 60 days of age (60 d) and 180 days of age (180 d) using RNA-Seq and whole genome bisulfite sequencing (WGBS). It was determined that DNA methylation primarily occurs in the cytosine-guanine (CG) context, and the analysis identified 106,282 differentially methylated regions (DMRs) corresponding to 12,385 differentially methylated genes (DMGs). Further integrated analysis of RNA-Seq and WGBS data revealed 1083 DMGs negatively correlated with the expression of DEGs. GO analysis showed that these genes were significantly enriched in spermatogenesis, germ cell development, and spermatid differentiation. The screening of enriched genes revealed that hyper-methylation repressed ADAM30, ADAM3A, DPY19L2, H2BC1, MAK, RPL10L, SPATA16, and YBX2, while hypo-methylation elevated CACNA1I, CADM1, CTNNB1, JAM2, and PAFAH1B3 expression. Additionally, the methylation status of the key genes ADAM3A, ADAM30, YBX2, JAM2, PAFAH1B3, and CTNNB1 was detected by bisulfite sequencing PCR (BSP). This study offers insights into the epigenetic regulation mechanisms underlying porcine testicular development.

DNA methylation analysis reveals local changes in resistant and susceptible soybean lines in response to Phytophthora sansomeana.

Phytophthora sansomeana is an emerging oomycete pathogen causing root rot in many agricultural species including soybean. However, as of now, only one potential resistance gene has been identified in soybean, and our understanding of how genetic and epigenetic regulation in soybean contributes to responses against this pathogen remains largely unknown. In this study, we performed whole genome bisulfite sequencing (WGBS) on two soybean lines, Colfax (resistant) and Williams 82 (susceptible), in response to P. sansomeana at two time points: 4 and 16 hours post-inoculation to compare their methylation changes. Our findings revealed that there were no significant changes in genome-wide CG, CHG (H = A, T, or C), and CHH methylation. However, we observed local methylation changes, specially an increase in CHH methylation around genes and transposable elements (TEs) after inoculation, which occurred earlier in the susceptible line and later in the resistant line. After inoculation, we identified differentially methylated regions (DMRs) in both Colfax and Williams 82, with a predominant presence in TEs. Notably, our data also indicated that more TEs exhibited changes in their methylomes in the susceptible line compared to the resistant line. Furthermore, we discovered 837 DMRs within or flanking 772 differentially expressed genes (DEGs) in Colfax and 166 DMRs within or flanking 138 DEGs in Williams 82. These DEGs had diverse functions, with Colfax primarily showing involvement in metabolic process, defense response, plant and pathogen interaction, anion and nucleotide binding, and catalytic activity, while Williams 82 exhibited a significant association with photosynthesis. These findings suggest distinct molecular responses to P. sansomeana infection in the resistant and susceptible soybean lines.

Females with autism spectrum disorders show stronger DNA methylation signatures than males in perinatal tissues.

Autism spectrum disorder (ASD) comprises a group of neurodevelopmental conditions currently diagnosed by behavioral assessment in childhood, with reported underdiagnosis in females. Though diagnosis in early life is linked to improved outcomes, we currently lack objective screening tools for newborns. To address this gap, we sought to identify a sex-specific DNA methylation signature for ASD using perinatal tissues that reflect dysregulation in the brain. DNA methylation was assayed from ASD and typically developing (TD) newborn blood, umbilical cord blood, placenta, and post-mortem cortex samples using whole genome bisulfite sequencing (WGBS) in a total of 511 samples. We found that methylation levels of differentially methylated regions (DMRs) differentiated samples by ASD diagnosis in females more than males across the perinatal tissues. We tested three theories for ASD sex differences in newborn blood, finding epigenetic support for an X chromosome-related female protective effect, as well as a high replication rate of DMRs (48.1%) in females across two independent cohorts. In our pan-tissue analysis, three genes (X-linked BCOR, GALNT9, OPCML) mapped to ASD DMRs replicated in all four female tissues. ASD DMRs from all tissues were enriched for neuro-related processes (females) and SFARI ASD-risk genes (females and males). Overall, we found a highly replicated methylation signature of ASD in females across perinatal tissues that reflected dysregulation in the brain and involvement of X chromosome epigenetics. This comparative study of perinatal tissues shows the promise of newborn blood DNA methylation biomarkers for early detection of females at risk for ASD and emphasizes the importance of sex-stratification in ASD studies.

Creation and validation of the first infinium DNA methylation array for the human imprintome

BackgroundDifferentially methylated imprint control regions (ICRs) regulate the monoallelic expression of imprinted genes. Their epigenetic dysregulation by environmental exposures throughout life results in the formation of common chronic diseases. Unfortunately, existing Infinium methylation arrays lack the ability to profile these regions adequately. Whole genome bisulfite sequencing (WGBS) is the unique method able to profile the ICRs. However, it is very expensive and it requires not only a high coverage, but it is also computationally intensive to assess these regions.FindingsTo address this deficiency, we developed a custom methylation array containing 22,819 probes. Among them, 10,438 are CG probes targeting unique CpG sites, with 9,757 probes successfully mapping to 1,088 out of the 1,488 candidate ICRs recently described. To assess the performance of the array, we created matched samples processed with the Human Imprintome array and WGBS, which is the current standard method for assessing the methylation of the Human Imprintome. We compared the methylation levels from the shared CpG sites, and obtained a mean R2 = 0.569. We also created matched samples processed with the Human Imprintome array and the Infinium Methylation EPIC v2 array, and obtained a mean R2 = 0.796. Furthermore, replication experiments demonstrated high reliability (ICC: 0.799–0.945).ConclusionsOur custom array will be useful for replicable and accurate assessment, mechanistic insight, and targeted investigation of ICRs. This tool should accelerate the discovery of ICRs associated with a wide range of diseases and exposures, and advance our understanding of genomic imprinting and its relevance in development and disease formation throughout the life course.

SETDB1 regulates short interspersed nuclear elements and chromatin loop organization in mouse neural precursor cells

BackgroundTransposable elements play a critical role in maintaining genome architecture during neurodevelopment. Short Interspersed Nuclear Elements (SINEs), a major subtype of transposable elements, are known to harbor binding sites for the CCCTC-binding factor (CTCF) and pivotal in orchestrating chromatin organization. However, the regulatory mechanisms controlling the activity of SINEs in the developing brain remains elusive.ResultsIn our study, we conduct a comprehensive genome-wide epigenetic analysis in mouse neural precursor cells using ATAC-seq, ChIP-seq, whole genome bisulfite sequencing, in situ Hi-C, and RNA-seq. Our findings reveal that the SET domain bifurcated histone lysine methyltransferase 1 (SETDB1)-mediated H3K9me3, in conjunction with DNA methylation, restricts chromatin accessibility on a selective subset of SINEs in neural precursor cells. Mechanistically, loss of Setdb1 increases CTCF access to these SINE elements and contributes to chromatin loop reorganization. Moreover, de novo loop formation contributes to differential gene expression, including the dysregulation of genes enriched in mitotic pathways. This leads to the disruptions of cell proliferation in the embryonic brain after genetic ablation of Setdb1 both in vitro and in vivo.ConclusionsIn summary, our study sheds light on the epigenetic regulation of SINEs in mouse neural precursor cells, suggesting their role in maintaining chromatin organization and cell proliferation during neurodevelopment.

Multi-omic characterization of allele-specific regulatory variation in hybrid pigs

Hybrid mapping is a powerful approach to efficiently identify and characterize genes regulated through mechanisms in cis. In this study, using reciprocal crosses of the phenotypically divergent Duroc and Lulai pig breeds, we perform a comprehensive multi-omic characterization of regulatory variation across the brain, liver, muscle, and placenta through four developmental stages. We produce one of the largest multi-omic datasets in pigs to date, including 16 whole genome sequenced individuals, as well as 48 whole genome bisulfite sequencing, 168 ATAC-Seq and 168 RNA-Seq samples. We develop a read count-based method to reliably assess allele-specific methylation, chromatin accessibility, and RNA expression. We show that tissue specificity was much stronger than developmental stage specificity in all of DNA methylation, chromatin accessibility, and gene expression. We identify 573 genes showing allele specific expression, including those influenced by parent-of-origin as well as allele genotype effects. We integrate methylation, chromatin accessibility, and gene expression data to show that allele specific expression can be explained in great part by allele specific methylation and/or chromatin accessibility. This study provides a comprehensive characterization of regulatory variation across multiple tissues and developmental stages in pigs.

Genetic and epigenetic differentiation in response to genomic selection for avian lay date.

Anthropogenic climate change has led to globally increasing temperatures at an unprecedented pace and, to persist, wild species have to adapt to their changing world. We, however, often fail to derive reliable predictions of species' adaptive potential. Genomic selection represents a powerful tool to investigate the adaptive potential of a species, but constitutes a 'blind process' with regard to the underlying genomic architecture of the relevant phenotypes. Here, we used great tit (Parus major) females from a genomic selection experiment for avian lay date to zoom into this blind process. We aimed to identify the genetic variants that responded to genomic selection and epigenetic variants that accompanied this response and, this way, might reflect heritable genetic variation at the epigenetic level. We applied whole genome bisulfite sequencing to blood samples of individual great tit females from the third generation of bidirectional genomic selection lines for early and late lay date. Genomic selection resulted in differences at both the genetic and epigenetic level. Genetic variants that showed signatures of selection were located within genes mostly linked to brain development and functioning, including LOC107203824 (SOX3-like). SOX3 is a transcription factor that is required for normal hypothalamo-pituitary axis development and functioning, an essential part of the reproductive axis. As for epigenetic differentiation, the early selection line showed hypomethylation relative to the late selection line. Sites with differential DNA methylation were located in genes important for various biological processes, including gonadal functioning (e.g., MSTN and PIK3CB). Overall, genomic selection for avian lay date provided insights into where within the genome the heritable genetic variation for lay date, on which selection can operate, resides and indicates that some of this variation might be reflected by epigenetic variants.

Oral pimonidazole unveils clinicopathologic and epigenetic features of hypoxic tumour aggressiveness in localized prostate cancer

BackgroundTumor hypoxia is associated with prostate cancer (PCa) treatment resistance and poor prognosis. Pimonidazole (PIMO) is an investigational hypoxia probe used in clinical trials. A better understanding of the clinical significance and molecular alterations underpinning PIMO-labeled tumor hypoxia is needed for future clinical application. Here, we investigated the clinical significance and molecular alterations underpinning PIMO-labeled tumor hypoxia in patients with localized PCa, in order to apply PIMO as a prognostic tool and to identify potential biomarkers for future clinical translation.MethodsA total of 39 patients with localized PCa were recruited and administered oral PIMO before undergoing radical prostatectomy (RadP). Immunohistochemical staining for PIMO was performed on 37 prostatectomy specimens with staining patterns evaluated and clinical association analyzed. Whole genome bisulfite sequencing was performed using laser-capture of microdissected specimen sections comparing PIMO positive and negative tumor areas. A hypoxia related methylation molecular signature was generated by integrating the differentially methylated regions with previously established RNA-seq datasets.ResultsThree PIMO staining patterns were distinguished: diffuse, focal, and comedo-like. The comedo-like staining pattern was more commonly associated with adverse pathology. PIMO-defined hypoxia intensity was positively correlated with advanced pathologic stage, tumor invasion, and cribriform and intraductal carcinoma morphology. The generated DNA methylation signature was found to be a robust hypoxia biomarker, which could risk-stratify PCa patients across multiple clinical datasets, as well as be applicable in other cancer types.ConclusionsOral PIMO unveiled clinicopathologic features of disease aggressiveness in localized PCa. The generated DNA methylation signature is a novel and robust hypoxia biomarker that has the potential for future clinical translation.

ATRT-16. MAPPING THE DISORDERED METHYLOME IN ATRT

Abstract BACKGROUND Atypical Teratoid Rhabdoid Tumor (ATRT) is an aggressive pediatric brain tumor driven by inactivation of the chromatin regulator SMARCB1, with remarkably few other mutations. SMARCB1 loss leads to epigenetic dysregulation which plays a pivotal role in ATRT development and progression, and DNA methylation profiling has distinguished molecular subgroups of this disease. However, studies that provide a static view of the epigenome may not capture epigenetic variability that drives cellular plasticity and heterogeneity. Mapping epigenetic variability is crucial to understanding the phenotypic plasticity of this tumor and its ability to evade therapeutic interventions in the clinic. METHODS We performed whole genome bisulfite sequencing (WGBS) on 22 ATRT samples, 4 normal fetal brain controls, and 1 hindbrain neural stem cell control, and applied a novel analytic pipeline (InformME) to generate DNA methylation potential energy landscapes encapsulating epigenetic stochasticity. This characterizes mean methylation as well as methylation entropy. We identified key genes and regulatory regions with the most stochastic methylation signatures. RESULTS We find that ATRT exhibits a globally destabilized epigenome, with markedly elevated methylation entropy (a measure of variability/disorder). Methylation discordance mapped to key regulatory regions including bivalent promoters, quiescent chromatin regions, and enhancers, and specific genes, with an enrichment for targets of the Polycomb repressive complex 2 and regulators of pluripotency and differentiation. CONCLUSIONS The disordered methylome of ATRT could play a critical role in the tumor’s ability to evade therapies. Mapping the targets of epigenetic instability can identify targets for intervention. Moreover, strategies to constrain epigenetic variability may have therapeutic benefit by limiting the number of epigenetic states available to cancer cells.