Genome-wide Methylation Profiling Research Articles

Early detection and prognosis evaluation for hepatocellular carcinoma by circulating tumour DNA methylation: A multicentre cohort study.

Early diagnosis of hepatocellular carcinoma (HCC) can significantly improve patient survival. We aimed to develop a blood-based assay to aid in the diagnosis, detection and prognostic evaluation of HCC. A three-phase multicentre study was conducted to screen, optimise and validate HCC-specific differentially methylated regions (DMRs) using next-generation sequencing and quantitative methylation-specific PCR (qMSP). Genome-wide methylation profiling was conducted to identify DMRs distinguishing HCC tumours from peritumoural tissues and healthy plasmas. The twenty most effective DMRs were verified and incorporated into a multilocus qMSP assay (HepaAiQ). The HepaAiQ model was trained to separate 293 HCC patients (Barcelona Clinic Liver Cancer (BCLC) stage 0/A, 224) from 266 controls including chronic hepatitis B (CHB) or liver cirrhosis (LC) (CHB/LC, 96), benign hepatic lesions (BHL, 23), and healthy controls (HC, 147). The model achieved an area under the curve (AUC) of 0.944 with a sensitivity of 86.0% in HCC and a specificity of 92.1% in controls. Blind validation of the HepaAiQ model in a cohort of 523 participants resulted in an AUC of 0.940 with a sensitivity of 84.4% in 205 HCC cases (BCLC stage 0/A, 167) and a specificity of 90.3% in 318 controls (CHB/LC, 100; BHL, 102; HC, 116). When evaluated in an independent test set, the HepaAiQ model exhibited a sensitivity of 70.8% in 65 HCC patients at BCLC stage 0/A and a specificity of 89.5% in 124 patients with CHB/LC. Moreover, HepaAiQ model was assessed in paired pre- and postoperative plasma samples from 103 HCC patients and correlated with 2-year patient outcomes. Patients with high postoperative HepaAiQ score showed a higher recurrence risk (Hazard ratio, 3.33, p<.001). HepaAiQ, a noninvasive qMSP assay, was developed to accurately measure HCC-specific DMRs and shows great potential for the diagnosis, detection and prognosis of HCC, benefiting at-risk populations.

Stimulatory effects of smoke solution and biogas digestate slurry application on photosynthesis, growth, and methylation profiling of solanum tuberosum

ABSTRACT Biostimulants are obtained from various sources like plants, animals, microorganisms, and industrial by-products as well as waste material. Their utilization in agriculture practices is being increased that is giving positive results. The purpose of the current study was to use plant-derived smoke (SMK) solution and biogas digestate (BGD) slurry as biostimulant to elucidate their impact on potato (Solanum tuberosum) performance. The experiment was conducted in lab as well as field conditions, and SMK and BGD solutions were prepared in varying concentrations such as SMK 1:500, SMK 1:250, BGD 50:50, and BGD 75:25. Foliar applications were performed thrice during experiments and data were collected related to photosynthesis, growth, pigments, and genome-wide methylation profiling. Net photosynthesis rate (A) and water use efficiency (WUE) were found higher in SMK- and BGD-treated lab and field grown plants. Among pigments, BGD-treated plants depicted higher levels of Chl a and Chl b while SMK-treated plants showed higher carotenoid levels. Alongside, enhancement in growth-related parameters like leaf number and dry weight was also observed in both lab- and field-treated plants. Furthermore, DNA methylation profile of SMK- and BGD-treated plants depicted variation compared to control. DNA methylation events increased in all the treatments compared to control except for SMK 1:500. These results indicate that smoke and slurry both act as efficient biostimulants which result in better performance of plants. Biostimulants also affected the genome-wide DNA methylation profile that resultantly might have changed the plant gene expression profiling and played its role in plant responsiveness to these biostimulants. However, there is need to elucidate a possible synergistic effect of SMK and BGD on plant growth along with gene expression profiling.

EpiDiP/NanoDiP: a versatile unsupervised machine learning edge computing platform for epigenomic tumour diagnostics

DNA methylation analysis based on supervised machine learning algorithms with static reference data, allowing diagnostic tumour typing with unprecedented precision, has quickly become a new standard of care. Whereas genome-wide diagnostic methylation profiling is mostly performed on microarrays, an increasing number of institutions additionally employ nanopore sequencing as a faster alternative. In addition, methylation-specific parallel sequencing can generate methylation and genomic copy number data. Given these diverse approaches to methylation profiling, to date, there is no single tool that allows (1) classification and interpretation of microarray, nanopore and parallel sequencing data, (2) direct control of nanopore sequencers, and (3) the integration of microarray-based methylation reference data. Furthermore, no software capable of entirely running in routine diagnostic laboratory environments lacking high-performance computing and network infrastructure exists. To overcome these shortcomings, we present EpiDiP/NanoDiP as an open-source DNA methylation and copy number profiling suite, which has been benchmarked against an established supervised machine learning approach using in-house routine diagnostics data obtained between 2019 and 2021. Running locally on portable, cost- and energy-saving system-on-chip as well as gpGPU-augmented edge computing devices, NanoDiP works in offline mode, ensuring data privacy. It does not require the rigid training data annotation of supervised approaches. Furthermore, NanoDiP is the core of our public, free-of-charge EpiDiP web service which enables comparative methylation data analysis against an extensive reference data collection. We envision this versatile platform as a useful resource not only for neuropathologists and surgical pathologists but also for the tumour epigenetics research community. In daily diagnostic routine, analysis of native, unfixed biopsies by NanoDiP delivers molecular tumour classification in an intraoperative time frame.

Early changes in tumor-naive cell-free methylomes and fragmentomes predict outcomes in pembrolizumab-treated solid tumors.

Early kinetics of circulating tumor DNA (ctDNA) in plasma predict response to pembrolizumab, but typically requires sequencing of matched tumor tissue or fixed gene panels. We analyzed genome-wide methylation and fragment length profiles using cell-free methylated DNA immunoprecipitation and sequencing (cfMeDIP-seq) in 204 plasma samples from 87 patients before and during treatment with pembrolizumab from a pan-cancer phase II investigator-initiated trial (INSPIRE). We trained a pan-cancer methylation signature using independent methylation array data from The Cancer Genome Atlas to quantify a cancer-specific methylation (CSM) and fragment length score (FLS) for each sample. CSM and FLS are strongly correlated with tumor-informed ctDNA levels. Early kinetics of CSM predict overall survival and progression-free survival, independently of tumor type, PD-L1, and tumor mutation burden. Early kinetics of FLS are associated with overall survival independently of CSM. Our tumor-naïve mutation-agnostic ctDNA approach integrating methylomics and fragmentomics could predict outcomes in patients treated with pembrolizumab.

Mother adversity and co-residence time impact mother–child similarity in genome-wide and gene-specific methylation profiles

BackgroundThe effects of adverse life events on physical and psychological health, with DNA methylation (DNAm) as a critical underlying mechanism, have been extensively studied. However, the epigenetic resemblance between mother and child in the context of neglectful caregiving, and whether it may be shaped by the emotional impact of maternal stressful events and the duration of co-residence (indexed by child age), remains unknown. The present study examined mother–child similarity in methylation profiles, considering the potential effect of mother adversity, mother empathy, neglect-control group, child age (an index of years of mother–child co-residence), and mother age. Using Illumina Epic arrays, we quantified DNAm in 115 mother–child saliva samples. We obtained a methylation similarity index by computing correlation coefficients between methylation profiles within dyads, for the entire epigenome, and five specific genes related to stress and empathy: NR3C1, FKPB5, OXTR, SCL6A4, and BDNF.ResultsThe methylation profiles of the mother–child familial pairs significantly correlated as compared to mother–child random pairs for the entire epigenome and NR3C1, FKBP5, OXTR and BDNF genes. Next, multiple linear regression models observed associations of mother adversity, child age, and neglect-control group on mother–child methylation similarity, only significant in mother–child familial pairs, after correcting for multiple comparisons. Higher mother adversity was associated with lower mother–child methylation similarity for the epigenome-wide analysis, for the BDNF gene, and in the neglect-control group for the OXTR gene. In turn, being an older child (longer co-residence) was associated with higher mother–child methylation similarity.ConclusionsMother adversity and co-residence time are modulating factors in the intergenerational methylation process that offer a window into development-dependent adaptations that can be affected by both hereditary and environmental factors, significantly observed only in biological dyads. A twofold implication for child well-being emerges, one is positive in that children of mothers exposed to life adversity or neglect did not necessarily inherit their methylation patterns. The other is concerning due to the influence of time spent living together, which affects similarity with the mother and potentially increases the risk of inheriting an epigenetic profile associated with future dysfunctional parenting patterns. This underscores the importance of the 'the earlier, the better' recommendation by the Child Protection System, which is not always followed.

DNA Methylation Signatures Correlate with Response to Immune Checkpoint Inhibitors in Metastatic Melanoma.

DNA methylation profiles have emerged as potential predictors of therapeutic response in various solid tumors. This study aimed to analyze the DNA methylation profiles of patients with stage IV metastatic melanoma undergoing first-line immune checkpoint inhibitor treatment and evaluate their correlation with a radiological response according to immune-related Response Evaluation Criteria in Solid Tumors (iRECIST). A total of 81 tissue samples from 71 patients with metastatic melanoma (27 female, 44 male) were included in this study. We utilized Illumina Methylation EPIC Beadchips to retrieve their genome-wide methylation profile by interrogating >850,000 CpG sites. Clustering based on the 500 most differentially methylated genes was conducted to identify distinct methylation patterns associated with immune checkpoint inhibitor response. Results were further aligned with an independent, previously published data set. The median progression-free survival was 8.5 months (range: 0-104.1 months), and the median overall survival was 30.6 months (range: 0-104.1 months). Objective responses were observed in 29 patients (40.8%). DNA methylation profiling revealed specific signatures that correlated with radiological response to immune checkpoint inhibitors. Three distinct clusters were identified based on the methylation patterns of the 500 most differentially methylated genes. Cluster 1 (12/12) and cluster 2 (12/24) exhibited a higher proportion of responders, while cluster 3 (39/45) predominantly consisted of non-responders. In the validation data set, responders also showed more frequent hypomethylation although differences in the data sets limit the interpretation. These findings suggest that DNA methylation profiling of tumor tissues might serve as a predictive biomarker for immune checkpoint inhibitor response in patients with metastatic melanoma. Further validation studies are warranted to confirm the efficiency of DNA methylation profiling as a predictive tool in the context of immunotherapy for metastatic melanoma.

Diagnostic Algorithm to Subclassify Atypical Spitzoid Tumors in Low and High Risk According to Their Methylation Status.

Current diagnostic algorithms are insufficient for the optimal clinical and therapeutic management of cutaneous spitzoid tumors, particularly atypical spitzoid tumors (AST). Therefore, it is crucial to identify new markers that allow for reliable and reproducible diagnostic assessment and can also be used as a predictive tool to anticipate the individual malignant potential of each patient, leading to tailored individual therapy. Using Reduced Representation Bisulfite Sequencing (RRBS), we studied genome-wide methylation profiles of a series of Spitz nevi (SN), spitzoid melanoma (SM), and AST. We established a diagnostic algorithm based on the methylation status of seven cg sites located in TETK4P2 (Tektin 4 Pseudogene 2), MYO1D (Myosin ID), and PMF1-BGLAP (PMF1-BGLAP Readthrough), which allows the distinction between SN and SM but is also capable of subclassifying AST according to their similarity to the methylation levels of Spitz nevi or spitzoid melanoma. Thus, our epigenetic algorithm can predict the risk level of AST and predict its potential clinical outcomes.

Genome-wide study reveals novel roles for formin-2 in axon regeneration as a microtubule dynamics regulator and therapeutic target for nerve repair

Peripheral nerves regenerate successfully; however, clinical outcome after injury is poor. We demonstrated that low-dose ionizing radiation (LDIR) promoted axon regeneration and function recovery after peripheral nerve injury (PNI). Genome-wide CpG methylation profiling identified LDIR-induced hypermethylation of the Fmn2 promoter, exhibiting injury-induced Fmn2 downregulation in dorsal root ganglia (DRGs). Constitutive knockout or neuronal Fmn2 knockdown accelerated nerve repair and function recovery. Mechanistically, increased microtubule dynamics at growth cones was observed in time-lapse imaging of Fmn2-deficient DRG neurons. Increased HDAC5 phosphorylation and rapid tubulin deacetylation were found in regenerating axons of neuronal Fmn2-knockdown mice after injury. Growth-promoting effect of neuronal Fmn2 knockdown was eliminated by pharmaceutical blockade of HDAC5 or neuronal Hdac5 knockdown, suggesting that Fmn2deletion promotes axon regeneration via microtubule post-translational modification. In silico screening of FDA-approved drugs identified metaxalone, administered either immediately or 24-h post-injury, accelerating function recovery. This work uncovers a novel axon regeneration function of Fmn2 and a small-molecule strategy for PNI.

High overall copy number variation burden by genome-wide methylation profiling holds negative prognostic value in surgically treated pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive type of cancer with an overall 5-year survival of around 10 %. New prognostic tools to stratify patients are needed. Our main aim was to evaluate the prognostic value of overall copy number variation (CNV) burden in surgically treated PDAC. DNA extracted from 108 surgical PDAC specimens was examined to collect data on the genome-wide DNA methylation status of >850,000 CpG sites in promoter, gene body, and enhancer regions (Illumina Infinium Methylation EPIC BeadChip Kit). CNV profiles were obtained and all PDACs were stratified into one of three groups: Low, moderate, or high overall CNV burden. Tumors histologically showing a dominant conventional and/or tubulopapillary pattern in 60 %–100 % and 0–59 % were categorized as Group A and Group B as per Kalimuthu. We also performed targeted next-generation sequencing (NGS) and immunohistochemistry. High overall CNV burden held independent negative prognostic value with poor survival (HR 4.01 (95%CI 1.96–8.19), p = 0.00014) and was more frequent in Group B (p = 0.0003). Most frequent chromosomal arm-level aberrations were gains of 8q (29 %) and 1q (19 %) and losses of 17p (55 %), 18q (43 %), 6q (37 %), 9p (36 %), 6p (26 %), 19p (26 %), and 8p (25 %). Most frequent mutations found were in KRAS (95 %), TP53 (62 %), CDKN2A (24 %), SMAD4 (23 %), ATM (9 %), ARID1A (7 %), RNF43 (7 %), GNAS (6 %), and KDM6A (6 %). Group A PDACs showed more frequently KRAS variants other than Gly12Val and Gly12Asp (p = 0.012). Our data indicate that overall CNV burden using genome-wide methylation profiling may be a useful prognostic tool in surgically treated PDAC. Importantly, our approach, using data from genome-wide methylation profiling for analysis of overall CNV burden, can be performed on formalin-fixed and paraffin embedded PDAC tissues. Future studies should examine the prognostic value of overall CNV burden in unresectable PDAC.

Screening and identifying of biomarkers in early colorectal cancer and adenoma based on genome-wide methylation profiles

BackgroundColorectal cancer is one of the most common malignant tumors worldwide with high morbidity and mortality. This study aimed to identify different methylation sites as new methylation markers in CRC and colorectal adenoma through tissue detection.MethodsDNA extraction and bisulfite modification as well as Infinium 450K methylation microarray detection were performed in 46 samples of sporadic colorectal cancer tissue, nine samples of colorectal adenoma, and 20 normal samples, and bioinformatic analysis was conducted involving genes enrichments of GO and KEGG. Pyrosequencing methylation detection was further performed in 68 sporadic colorectal cancer tissues, 31 samples of colorectal adenoma, and 49 normal colorectal mucosae adjacent to carcinoma to investigate the differentially methylated genes obtained from methylation microarray.ResultsThere were 65,535 differential methylation marker probes, among which 25,464 were hypermethylated markers and 40,071 were hypomethylated markers in the adenoma compared with the normal group, and 395,571 were differentially methylated markers in patients with sporadic colorectal cancer compared with the normal group, including 21,710 hypermethylated markers and 17,861 hypomethylated markers. Five hypermethylated genes including ZNF471, SND1, SPOCK1, FBLIM1, and OTX1 were detected and confirmed in 68 cases of colorectal cancer, 31 cases of adenoma, and 49 cases of normal control group.ConclusionsHypermethylated genes of ZNF471, SND1, SPOCK1, FBLIM1, and OTX1 were obtained from methylation chip detection and further confirm analysis in colorectal cancer and adenoma compared with normal tissue, which may be promising diagnostic markers of colorectal cancer and colorectal adenoma.

Analytical sensitivity of a method is critical in detection of low-level BRCA1 constitutional epimutation

Recent reports based on a substantial number of cases, warrant need for population-based research to determine implications of constitutional methylation of tumor suppressor genes such as BRCA1 occurring in healthy tissue in the prediction of cancer. However, the detection of the constitutional methylation in DNA extracted from blood has already been shown to be technologically challenging, mainly because epimutations appear to be present in blood at a very low level. The analytical sensitivity required for low-level methylation detection can be provided by NGS, but this technique is still labor and cost-intensive. We assessed if PCR-based MS-HRM and BeadChip microarray technologies, which are standardized and cost-effective technologies for methylation changes screening, provide a sufficient level of analytical sensitivity for constitutional BRCA1 methylation detection in blood samples. The study included whole blood samples from 67 healthy women, 35 with previously confirmed and 32 with no detectable BRCA1 promoter methylation for which we performed both MS-HRM based BRCA1 gene methylation screening and genome wide methylation profiling with EPIC microarray. Our results shown, that low-level BRCA1 methylation can be effectively detected in DNA extracted from blood by PCR-based MS-HRM. At the same time, EPIC microarray does not provide conclusive results to unambiguously determine the presence of BRCA1 constitutional methylation in MS-HRM epimutation positive samples. The analytical sensitivity of MS-HRM is sufficient to detect low level BRCA1 constitutional epimutation in DNA extracted from blood and BeadChip technology-based microarrays appear not to provide that level of analytical sensitivity.

Functional Insight into and Refinement of the Genomic Boundaries of the JARID2-Neurodevelopmental Disorder Episignature.

JARID2 (Jumonji, AT-rich interactive domain 2) haploinsufficiency is associated with a clinically distinct neurodevelopmental syndrome. It is characterized by intellectual disability, developmental delay, autistic features, behavior abnormalities, cognitive impairment, hypotonia, and dysmorphic features. JARID2 acts as a transcriptional repressor protein that is involved in the regulation of histone methyltransferase complexes. JARID2 plays a role in the epigenetic machinery, and the associated syndrome has an identified DNA methylation episignature derived from sequence variants and intragenic deletions involving JARID2. For this study, our aim was to determine whether patients with larger deletions spanning beyond JARID2 present a similar DNA methylation episignature and to define the critical region involved in aberrant DNA methylation in 6p22-p24 microdeletions. We examined the DNA methylation profiles of peripheral blood from 56 control subjects, 13 patients with (likely) pathogenic JARID2 variants or patients carrying copy number variants, and three patients with JARID2 VUS variants. The analysis showed a distinct and strong differentiation between patients with (likely) pathogenic variants, both sequence and copy number, and controls. Using the identified episignature, we developed a binary model to classify patients with the JARID2-neurodevelopmental syndrome. DNA methylation analysis indicated that JARID2 is the driver gene for aberrant DNA methylation observed in 6p22-p24 microdeletions. In addition, we performed analysis of functional correlation of the JARID2 genome-wide methylation profile with the DNA methylation profiles of 56 additional neurodevelopmental disorders. To conclude, we refined the critical region for the presence of the JARID2 episignature in 6p22-p24 microdeletions and provide insight into the functional changes in the epigenome observed when regulation by JARID2 is lost.

Genome-Wide DNA Methylation and Gene Expression in Patients with Indolent Systemic Mastocytosis.

Mastocytosis is a clinically heterogenous, usually acquired disease of the mast cells with a survival time that depends on the time of onset. It ranges from skin-limited to systemic disease, including indolent and more aggressive variants. The presence of the oncogenic KIT p. D816V gene somatic mutation is a crucial element in the pathogenesis. However, further epigenetic regulation may also affect the expression of genes that are relevant to the pathology. Epigenetic alterations are responsible for regulating the expression of genes that do not modify the DNA sequence. In general, it is accepted that DNA methylation inhibits the binding of transcription factors, thereby down-regulating gene expression. However, so far, little is known about the epigenetic factors leading to the clinical onset of mastocytosis. Therefore, it is essential to identify possible epigenetic predictors, indicators of disease progression, and their link to the clinical picture to establish appropriate management and a therapeutic strategy. The aim of this study was to analyze genome-wide methylation profiles to identify differentially methylated regions (DMRs) in patients with mastocytosis compared to healthy individuals, as well as the genes located in those regulatory regions. Genome-wide DNA methylation profiling was performed in peripheral blood collected from 80 adult patients with indolent systemic mastocytosis (ISM), the most prevalent subvariant of mastocytosis, and 40 healthy adult volunteers. A total of 117 DNA samples met the criteria for the bisulfide conversion step and microarray analysis. Genome-wide DNA methylation analysis was performed using a MethylationEPIC BeadChip kit. Further analysis was focused on the genomic regions rather than individual CpG sites. Co-methylated regions (CMRs) were assigned via the CoMeBack method. To identify DMRs between the groups, a linear regression model with age as the covariate on CMRs was performed using Limma. Using the available data for cases only, an association analysis was performed between methylation status and tryptase levels, as well as the context of allergy, and anaphylaxis. KEGG pathway mapping was used to identify genes differentially expressed in anaphylaxis. Based on the DNA methylation results, the expression of 18 genes was then analyzed via real-time PCR in 20 patients with mastocytosis and 20 healthy adults. A comparison of the genome-wide DNA methylation profile between the mastocytosis patients and healthy controls revealed significant differences in the methylation levels of 85 selected CMRs. Among those, the most intriguing CMRs are 31 genes located within the regulatory regions. In addition, among the 10 CMRs located in the promoter regions, 4 and 6 regions were found to be either hypo- or hypermethylated, respectively. Importantly, three oncogenes-FOXQ1, TWIST1, and ERG-were identified as differentially methylated in mastocytosis patients, for the first time. Functional annotation revealed the most important biological processes in which the differentially methylated genes were involved as transcription, multicellular development, and signal transduction. The biological process related to histone H2A monoubiquitination (GO:0035518) was found to be enriched in association with higher tryptase levels, which may be associated with more aberrant mast cells and, therefore, more atypical mast cell disease. The signal in the BAIAP2 gene was detected in the context of anaphylaxis, but no significant differential methylation was found in the context of allergy. Furthermore, increased expression of genes encoding integral membrane components (GRM2 and KRTCAP3) was found in mastocytosis patients. This study confirms that patients with mastocytosis differ significantly in terms of methylation levels in selected CMRs of genes involved in specific molecular processes. The results of gene expression profiling indicate the increased expression of genes belonging to the integral component of the membrane in mastocytosis patients (GRM2 and KRTCAP3). Further work is warranted, especially in relation to the disease subvariants, to identify links between the methylation status and the symptoms and novel therapeutic targets.

Genome-wide methylation patterns from canine nanopore assemblies.

Recent advances in long-read sequencing have enabled the creation of reference-quality genome assemblies for multiple individuals within a species. In particular, eight long-read genome assemblies have recently been published for the canine model (dogs and wolves). These assemblies were created using a range of sequencing and computational approaches, with only limited comparisons described among subsets of the assemblies. Here we present three high-quality de novo reference assemblies based upon Oxford Nanopore long-read sequencing: two Bernese Mountain Dogs (BD & OD) and a Cairn terrier (CA611). These breeds are of particular interest due to the enrichment of unresolved genetic disorders. Leveraging advancement in software technologies, we utilized published data of Labrador Retriever (Yella) to generate a new assembly, resulting in a ∼280-fold increase in continuity (N50 size of 91 kbp vs 25.75 Mbp). In conjunction with these four new assemblies, we uniformly assessed eight existing assemblies for generalized quality metrics, sequence divergence, and a detailed BUSCO assessment. We identified a set of ∼400 conserved genes during the BUSCO analysis missing in all assemblies. Genome-wide methylation profiles were generated from the nanopore sequencing, resulting in broad concordance with existing whole-genome and reduced representation bisulfite sequencing, while highlighting superior overage of mobile elements. These analyses demonstrate the ability of Nanopore sequencing to resolve sequence and epigenetic profile of canine genomes.

Abstract P390: E-value Over P-value And Its Adjustments: Improving Dna Methylation Analysis With An Application To Hypertension Studies

DNA methylation plays a crucial role in transcriptional regulation. Reduced representation bisulfite sequencing (RRBS) is a technique of increasing use for analyzing genome-wide methylation profiles. Many computational tools such as Metilene, MethylKit, BiSeq and DMRfinder have been developed to use RRBS data to detect differentially methylated regions (DMRs) that may be involved in epigenetic regulations of gene expression. For DMR detection tools, as for countless other medical literature applications, P-values and their adjustments are among the most widely reported statistics used to assess the statistical significance of biological findings. However, P-values are coming under increasing criticism related to their questionable accuracy and relatively high levels of false positive or negative indications. Here, we proposed a method to calculate E-value for DMR detection in RRBS data, which is defined as likelihood ratios falling into the null hypothesis over the entire parameter space. We also provided the corresponding R package ‘ metevalue ’ as a user-friendly interface to implement E-value calculations into various DMR detection tools. To evaluate the performance of E-value, we generated various RRBS benchmarking datasets using our simulator ‘ RRBSsim ’ with 8 samples in each experimental group. Our comprehensive benchmarking analyses showed that using E-value not only significantly improved accuracy, AUC and power, over that of P-value or adjusted P-value, but also reduced false discovery rates and Type I errors. To illustrate the utility of E-value, we applied it to identify DMRs in two real RRBS datasets. One was genome-wide DNA methylome and transcriptome of renal T lymphocytes from Dahl Salt Sensitive rats treated with high- or low-salt diets. The other was genome-wide DNA methylation and gene expression data of human arterioles from our clinical trial in which 10 subjects were placed on 2-week low-salt diet (1200-mg sodium/day). Compared to the classical adjusted P-value, the use of E-value detected biologically more relevant DMRs and improved the negative association between DNA methylation and gene expression.

Genome-wide methylation profiling identified methylated KCNA3 and OTOP2 as promising diagnostic markers for esophageal squamous cell carcinoma.

Early detection of esophageal squamous cell carcinoma (ESCC) can considerably improve the prognosis of patients. Aberrant cell-free DNA (cfDNA) methylation signatures are a promising tool for detecting ESCC. However, available markers based on cell-free DNA methylation are still inadequate. This study aimed to identify ESCC-specific cfDNA methylation markers and evaluate the diagnostic performance in the early detection of ESCC. We performed whole-genome bisulfite sequencing (WGBS) for 24 ESCC tissues and their normal adjacent tissues. Based on the WGBS data, we identified 21,469,837 eligible CpG sites (CpGs). By integrating several methylation datasets, we identified several promising ESCC-specific cell-free DNA methylation markers. Finally, we developed a dual-marker panel based on methylated KCNA3 and OTOP2, and then, we evaluated its performance in our training and validation cohorts. The ESCC diagnostic model constructed based on KCNA3 and OTOP2 had an AUC of 0.91 [95% CI: 0.85-0.95], and an optimal sensitivity and specificity of 84.91% and 94.32%, respectively, in the training cohort. In the independent validation cohort, the AUC was 0.88 [95% CI: 0.83-0.92], along with an optimal sensitivity of 81.5% and specificity of 92.9%. The model sensitivity for stage I-II ESCC was 78.4%, which was slightly lower than the sensitivity of the model (85.7%) for stage III-IV ESCC. The dual-target panel based on cfDNA showed excellent performance for detecting ESCC and might be an alternative strategy for screening ESCC.

Pan-cancer analysis of genome-wide methylation profiling discover type-specific markers targeting circulating free DNA for the detection of colorectal cancer.

Pan-cancer analysis of genome-wide methylation profiling discover type-specific markers targeting circulating free DNA for the detection of colorectal cancer.

Genome-wide methylation profile of mitochondrial DNA across bovine preimplantation development

ABSTRACT This study characterized variations in the methylation profile of mitochondrial DNA (mtDNA) during initial bovine embryo development and correlated the presence of methylation with mtDNA transcription. Bovine oocytes were obtained from abattoir ovaries and submitted to in vitro culture procedures. Oocytes and embryos were collected at various stages (immature oocyte, IM; mature oocyte, MII; zygote, ZY; 4-cells, 4C; 16-cells, 16C and blastocysts, BL). Total DNA (including mtDNA) was used for Whole Genome Enzymatic Methyl Sequencing and for quantification of mtDNA copy number. Extracted RNA was used for quantification of mitochondrial transcripts using Droplet Digital PCR. We selected ND6, CYTB, tRNA-Phe and tRNA-Gln based on their location in the mitochondrial genome, functionality and/or previous literature associating these regions with cytosine methylation. The number of mtDNA copies per oocyte/embryo was found to be similar, while methylation levels in mtDNA varied among stages. Higher total methylation levels were found mainly at 4C and 16C. In specific gene regions, higher methylation levels were also observed at 4C and 16C (ND6, CYTB and tRNA-Phe), as well as an inverse correlation with the quantity of transcripts for these regions. This is a first description of epigenetic changes occurring in mtDNA during early embryonic development. Our results indicate that methylation might regulate the mtDNA transcription at a local level, particularly around the time of embryonic genome activation.

Germline pathogenic variants in HNRNPU are associated with alterations in blood methylome

HNRNPU encodes a multifunctional RNA-binding protein that plays critical roles in regulating pre-mRNA splicing, mRNA stability, and translation. Aberrant expression and dysregulation of HNRNPU have been implicated in various human diseases, including cancers and neurological disorders. We applied a next generation sequencing based assay (EPIC-NGS) to investigate genome-wide methylation profiling for >2 M CpGs for 7 individuals with a neurodevelopmental disorder associated with HNRNPU germline pathogenic loss-of-function variants. Compared to healthy individuals, 227 HNRNPU-associated differentially methylated positions were detected. Both hyper- and hypomethylation alterations were identified but the former predominated. The identification of a methylation episignature for HNRNPU-associated neurodevelopmental disorder (NDD) implicates HNPRNPU-related chromatin alterations in the aetiopathogenesis of this disorder and suggests that episignature profiling should have clinical utility as a predictor for the pathogenicity of HNRNPU variants of uncertain significance. The detection of a methylation episignaure for HNRNPU-associated NDD is consistent with a recent report of a methylation episignature for HNRNPK-associated NDD.

E-value: a superior alternative to P-value and its adjustments in DNA methylation studies.

DNA methylation plays a crucial role in transcriptional regulation. Reduced representation bisulfite sequencing (RRBS) is a technique of increasing use for analyzing genome-wide methylation profiles. Many computational tools such as Metilene, MethylKit, BiSeq and DMRfinder have been developed to use RRBS data for the detection of the differentially methylated regions (DMRs) potentially involved in epigenetic regulations of gene expression. For DMR detection tools, as for countless other medical applications, P-values and their adjustments are among the most standard reporting statistics used to assess the statistical significance of biological findings. However, P-values are coming under increasing criticism relating to their questionable accuracy and relatively high levels of false positive or negative indications. Here, we propose a method to calculate E-values, as likelihood ratios falling into the null hypothesis over the entire parameter space, for DMR detection in RRBS data. We also provide the R package 'metevalue' as a user-friendly interface to implement E-value calculations into various DMR detection tools. To evaluate the performance of E-values, we generated various RRBS benchmarking datasets using our simulator 'RRBSsim' with eight samples in each experimental group. Our comprehensive benchmarking analyses showed that using E-values not only significantly improved accuracy, area under ROC curve and power, over that of P-values or adjusted P-values, but also reduced false discovery rates and type I errors. In applications using real RRBS data of CRL rats and a clinical trial on low-salt diet, the use of E-values detected biologically more relevant DMRs and also improved the negative association between DNA methylation and gene expression.