Differential methylation analysis for bisulfite sequencing using DSS

Exploring early DNA methylation alterations in type 1 diabetes: implications of glycemic control

Background

Prolonged hyperglycemia causes diabetes-related micro- and macrovascular complications, which combined represent a significant burden for individuals living with diabetes. The growing scope of evidence indicates that hyperglycemia affects the development of vascular complications through DNA methylation.

Methods

A genome-wide differential DNA methylation analysis was performed on pooled peripheral blood DNA samples from individuals with type 1 diabetes (T1D) with direct DNA sequencing. Strict selection criteria were used to ensure two age- and sex-matched groups with no clinical signs of chronic complications according to persistent mean glycated hemoglobin (HbA1c) values over 5 years: HbA1c<7% (N=10) and HbA1c>8% (N=10).

Results

Between the two groups, 8385 differentially methylated CpG sites, annotated to 1802 genes, were identified. Genes annotated to hypomethylated CpG sites were enriched in 48 signaling pathways. Further analysis of key CpG sites revealed four specific regions, two of which were hypermethylated and two hypomethylated, associated with long non-coding RNA and processed pseudogenes.

Conclusions

Prolonged hyperglycemia in individuals with T1D, who have no clinical manifestation of diabetes-related complications, is associated with multiple differentially methylated CpG sites in crucial genes and pathways known to be linked to chronic complications in T1D.

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

Motivation

Whole-genome bisulfite sequencing is a powerful tool for analyzing chromatin methylation genome-wide, but analysis of whole-genome bisulfite data is hampered by slow, inaccurate, and inflexible pipelines.

Results

We developed PCBS, 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 and 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.

Availability and Implementation

PCBS is an R package available under a GNU GPLv3 license at: https://github.com/katlande/PCBS and from CRAN: https://CRAN.R-project.org/package=PCBS. Instructions for use are available at: https://katlande.github.io/PCBS/.

Cell-free DNA methylation patterns in aging and their association with inflamm-aging

Aim: Liquid biopsies analyzing cell-free DNA (cfDNA) methylation in plasma offer a noninvasive diagnostic for diseases, with the potential of aging biomarkers underexplored. Methods: Utilizing enzymatic methyl-seq (EM-seq), this study assessed cfDNA methylation patterns in aging with blood from 35 healthy individuals. Results: It found aging signatures, including higher cfDNA levels and variations in fragment sizes, plus approximately 2000 age-related differentially methylated CpG sites. A biological age predictive model based on 48 CpG sites showed a strong correlation with chronological age, verified by two datasets. Age-specific epigenetic shifts linked to inflammation were revealed through differentially methylated regions profiling and Olink proteomics. Conclusion: These findings suggest cfDNA methylation as a potential aging biomarker and might exacerbate immunoinflammatory reactivity in older individuals.

Computational deconvolution of DNA methylation data from mixed DNA samples

In this review, we provide a comprehensive overview of the different computational tools that have been published for the deconvolution of bulk DNA methylation (DNAm) data. Here, deconvolution refers to the estimation of cell-type proportions that constitute a mixed sample. The paper reviews and compares 25 deconvolution methods (supervised, unsupervised or hybrid) developed between 2012 and 2023 and compares the strengths and limitations of each approach. Moreover, in this study, we describe the impact of the platform used for the generation of methylation data (including microarrays and sequencing), the applied data pre-processing steps and the used reference dataset on the deconvolution performance. Next to reference-based methods, we also examine methods that require only partial reference datasets or require no reference set at all. In this review, we provide guidelines for the use of specific methods dependent on the DNA methylation data type and data availability.

Comparative study on genomic and epigenomic profiles of retinoblastoma or tuberous sclerosis complex via nanopore sequencing and a joint screening framework

Recurrence and extraocular metastasis in advanced intraocular retinoblastoma (RB) are still major obstacles for successful treatment of Chinese children. Tuberous sclerosis complex (TSC) is a very rare, multisystemic genetic disorder characterized by hamartomatous growth. In this study, we aimed to compare genomic and epigenomic profiles with human RB or TSC using recently developed nanopore sequencing, and to identify disease-associated variations or genes. Peripheral blood samples were collected from either RB or RB/TSC patients plus their normal siblings, followed by nanopore sequencing and identification of disease-specific structural variations (SVs) and differentially methylated regions (DMRs) by a systematic biology strategy named as multiomics-based joint screening framework. In total, 316 RB- and 1295 TSC-unique SVs were identified, as well as 1072 RB- and 1114 TSC-associated DMRs, respectively. We eventually identified 6 key genes for RB for further functional validation. Knockdown of CDK19 with specific siRNAs significantly inhibited Y79 cellular proliferation and increased sensitivity to carboplatin, whereas downregulation of AHNAK2 promoted the cell growth as well as drug resistance. Those two genes might serve as potential diagnostic markers or therapeutic targets of RB. The systematic biology strategy combined with functional validation might be an effective approach for rare pediatric malignances with limited samples and challenging collection process.

Cold exposure impacts DNA methylation patterns in cattle sperm

DNA methylation is influenced by various exogenous factors such as nutrition, temperature, toxicants, and stress. Bulls from the Pacific Northwest region of the United States and other northern areas are exposed to extreme cold temperatures during winter. However, the effects of cold exposure on the methylation patterns of bovine sperm remain unclear. To address, DNA methylation profiles of sperm collected during late spring and winter from the same bulls were analyzed using whole genome bisulfite sequencing (WGBS). Bismark (0.22.3) were used for mapping the WGBS reads and R Bioconductor package DSS was used for differential methylation analysis. Cold exposure induced 3,163 differentially methylated cytosines (DMCs) with methylation difference ≥10% and a q-value < 0.05. We identified 438 differentially methylated regions (DMRs) with q-value < 0.05, which overlapped with 186 unique genes. We also identified eight unique differentially methylated genes (DMGs) (Pax6, Macf1, Mest, Ubqln1, Smg9, Ctnnb1, Lsm4, and Peg10) involved in embryonic development, and nine unique DMGs (Prmt6, Nipal1, C21h15orf40, Slc37a3, Fam210a, Raly, Rgs3, Lmbr1, and Gan) involved in osteogenesis. Peg10 and Mest, two paternally expressed imprinted genes, exhibited >50% higher methylation. The differential methylation patterns of six distinct DMRs: Peg10, Smg9 and Mest related to embryonic development and Lmbr1, C21h15orf40 and Prtm6 related to osteogenesis, were assessed by methylation-specific PCR (MS-PCR), which confirmed the existence of variable methylation patterns in those locations across the two seasons. In summary, cold exposure induces differential DNA methylation patterns in genes that appear to affect embryonic development and osteogenesis in the offspring. Our findings suggest the importance of replicating the results of the current study with a larger sample size and exploring the potential of these changes in affecting offspring development.

Nanopore-based DNA long-read sequencing analysis of the aged human brain

Aging disrupts cellular processes such as DNA repair and epigenetic control, leading to a gradual buildup of genomic alterations that can have detrimental effects in post-mitotic cells. Genomic alterations in regions of the genome that are rich in repetitive sequences, often termed "dark loci," are difficult to resolve using traditional sequencing approaches. New long-read technologies offer promising avenues for exploration of previously inaccessible regions of the genome. Using nanopore-based long-read whole-genome sequencing of DNA extracted from aged 18 human brains, we identify previously unreported structural variants and methylation patterns within repetitive DNA, focusing on transposable elements ("jumping genes") as crucial sources of variation, particularly in dark loci. Our analyses reveal potential somatic insertion variants and provides DNA methylation frequencies for many retrotransposon families. We further demonstrate the utility of this technology for the study of these challenging genomic regions in brains affected by Alzheimer's disease and identify significant differences in DNA methylation in pathologically normal brains versus those affected by Alzheimer's disease. Highlighting the power of this approach, we discover specific polymorphic retrotransposons with altered DNA methylation patterns. These retrotransposon loci have the potential to contribute to pathology, warranting further investigation in Alzheimer's disease research. Taken together, our study provides the first long-read DNA sequencing-based analysis of retrotransposon sequences, structural variants, and DNA methylation in the aging brain affected with Alzheimer's disease neuropathology.

Genomic structures and regulation patterns at HPV integration sites in cervical cancer

Human papillomavirus (HPV) integration has been implicated in transforming HPV infection into cancer, but its genomic consequences have been difficult to study using short-read technologies. To resolve the dysregulation associated with HPV integration, we performed long-read sequencing on 63 cervical cancer genomes. We identified six categories of integration events based on HPV-human genomic structures. Of all HPV integrants, defined as two HPV-human breakpoints bridged by an HPV sequence, 24% contained variable copies of HPV between the breakpoints, a phenomenon we termed heterologous integration. Analysis of DNA methylation within and in proximity to the HPV genome at individual integration events revealed relationships between methylation status of the integrant and its orientation and structure. Dysregulation of the human epigenome and neighboring gene expression in cis with the HPV-integrated allele was observed over megabase-ranges of the genome. By elucidating the structural, epigenetic, and allele-specific impacts of HPV integration, we provide insight into the role of integrated HPV in cervical cancer.

Multi-omics comparison of malignant and normal uveal melanocytes reveals molecular features of uveal melanoma

Uveal melanoma (UM) is a rare cancer resulting from the transformation of melanocytes in the uveal tract. Integrative analysis has identified four molecular and clinical subsets of UM. To improve our molecular understanding of UM, we performed extensive multi-omics characterization comparing two aggressive UM patient-derived xenograft models with normal choroidal melanocytes, including DNA optical mapping, specific histone modifications, and DNA topology analysis using Hi-C. Our gene expression and cytogenetic analyses suggest that genomic instability is a hallmark of UM. We also identified a recurrent deletion in the BAP1 promoter resulting in loss of expression and associated with high risk of metastases in UM patients. Hi-C revealed chromatin topology changes associated with the upregulation of PRAME, an independent prognostic biomarker in UM, and a potential therapeutic target. Our findings illustrate how multi-omics approaches can improve our understanding of tumorigenesis and reveal two distinct mechanisms of gene expression dysregulation in UM.

Herbivory induced methylation changes in the Lombardy poplar: A comparison of results obtained by epiGBS and WGBS

DNA cytosine methylation is an epigenetic mechanism involved in regulation of plant responses to biotic and abiotic stress and its ability to change can vary with the sequence context in which a cytosine appears (CpG, CHG, CHH, where H = Adenine, Thymine, Cytosine). Quantification of DNA methylation in model plant species is frequently addressed by Whole Genome Bisulfite Sequencing (WGBS), which requires a good-quality reference genome. Reduced Representation Bisulfite Sequencing (RRBS) is a cost-effective potential alternative for ecological research with limited genomic resources and large experimental designs. In this study, we provide for the first time a comprehensive comparison between the outputs of RRBS and WGBS to characterize DNA methylation changes in response to a given environmental factor. In particular, we used epiGBS (recently optimized RRBS) and WGBS to assess global and sequence-specific differential methylation after insect and artificial herbivory in clones of Populus nigra cv. 'italica'. We found that, after any of the two herbivory treatments, global methylation percentage increased in CHH, and the shift was detected as statistically significant only by epiGBS. As regards to loci-specific differential methylation induced by herbivory (cytosines in epiGBS and regions in WGBS), both techniques indicated the specificity of the response elicited by insect and artificial herbivory, together with higher frequency of hypo-methylation in CpG and hyper-methylation in CHH. Methylation changes were mainly found in gene bodies and intergenic regions when present at CpG and CHG and in transposable elements and intergenic regions at CHH context. Thus, epiGBS succeeded to characterize global, genome-wide methylation changes in response to herbivory in the Lombardy poplar. Our results support that epiGBS could be particularly useful in large experimental designs aimed to explore epigenetic changes of non-model plant species in response to multiple environmental factors.

Uropathogenic Escherichia coli infection-induced epithelial trained immunity impacts urinary tract disease outcome

Previous urinary tract infections (UTIs) can predispose one to future infections; however, the underlying mechanisms affecting recurrence are poorly understood. We previously found that UTIs in mice cause differential bladder epithelial (urothelial) remodelling, depending on disease outcome, that impacts susceptibility to recurrent UTI. Here we compared urothelial stem cell (USC) lines isolated from mice with a history of either resolved or chronic uropathogenic Escherichia coli (UPEC) infection, elucidating evidence of molecular imprinting that involved epigenetic changes, including differences in chromatin accessibility, DNA methylation and histone modification. Epigenetic marks in USCs from chronically infected mice enhanced caspase-1-mediated cell death upon UPEC infection, promoting bacterial clearance. Increased Ptgs2os2 expression also occurred, potentially contributing to sustained cyclooxygenase-2 expression, bladder inflammation and mucosal wounding-responses associated with severe recurrent cystitis. Thus, UPEC infection acts as an epi-mutagen reprogramming the urothelial epigenome, leading to urothelial-intrinsic remodelling and training of the innate response to subsequent infection.

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

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

ARIC: accurate and robust inference of cell type proportions from bulk gene expression or DNA methylation data

Quantifying cell proportions, especially for rare cell types in some scenarios, is of great value in tracking signals associated with certain phenotypes or diseases. Although some methods have been proposed to infer cell proportions from multicomponent bulk data, they are substantially less effective for estimating the proportions of rare cell types which are highly sensitive to feature outliers and collinearity. Here we proposed a new deconvolution algorithm named ARIC to estimate cell type proportions from gene expression or DNA methylation data. ARIC employs a novel two-step marker selection strategy, including collinear feature elimination based on the component-wise condition number and adaptive removal of outlier markers. This strategy can systematically obtain effective markers for weighted $\upsilon$-support vector regression to ensure a robust and precise rare proportion prediction. We showed that ARIC can accurately estimate fractions in both DNA methylation and gene expression data from different experiments. We further applied ARIC to the survival prediction of ovarian cancer and the condition monitoring of chronic kidney disease, and the results demonstrate the high accuracy and robustness as well as clinical potentials of ARIC. Taken together, ARIC is a promising tool to solve the deconvolution problem of bulk data where rare components are of vital importance.