Rockfish: A transformer-based model for accurate 5-methylcytosine prediction from nanopore sequencing.

DNA methylation plays an important role in various biological processes, including cell differentiation, ageing, and cancer development. The most important methylation in mammals is 5-methylcytosine mostly occurring in the context of CpG dinucleotides. Sequencing methods such as whole-genome bisulfite sequencing successfully detect 5-methylcytosine DNA modifications. However, they suffer from the serious drawbacks of short read lengths and might introduce an amplification bias. Here we present Rockfish, a deep learning algorithm that significantly improves read-level 5-methylcytosine detection by using Nanopore sequencing. Rockfish is compared with other methods based on Nanopore sequencing on R9.4.1 and R10.4.1 datasets. There is an increase in the single-base accuracy and the F1 measure of up to 5 percentage points on R.9.4.1 datasets, and up to 0.82 percentage points on R10.4.1 datasets. Moreover, Rockfish shows a high correlation with whole-genome bisulfite sequencing, requires lower read depth, and achieves higher confidence in biologically important regions such as CpG-rich promoters while being computationally efficient. Its superior performance in human and mouse samples highlights its versatility for studying 5-methylcytosine methylation across varied organisms and diseases. Finally, its adaptable architecture ensures compatibility with new versions of pores and chemistry as well as modification types.

NSUN6-mediated 5-methylcytosine modification of NDRG1 mRNA promotes radioresistance in cervical cancer.

BACKGROUND: Radioresistance is the leading cause of death in advanced cervical cancer (CC). Dysregulation of RNA modification has recently emerged as a regulatory mechanism in radiation and drug resistance. We aimed to explore the biological function and clinical significance of 5-methylcytosine (m5C) in cervical cancer radiosensitivity. METHODS: The abundance of RNA modification in radiotherapy-resistant and sensitive CC specimens was quantified by liquid chromatography-tandem mass spectrometry. The essential RNA modification-related genes involved in CC radiosensitivity were screened via RNA sequencing. The effect of NSUN6 on radiosensitivity was verified in CC cell lines, cell-derived xenograft (CDX), and 3D bioprinted patient-derived organoid (PDO). The mechanisms of NSUN6 in regulating CC radiosensitivity were investigated by integrative m5C sequencing, mRNA sequencing, and RNA immunoprecipitation. RESULTS: We found a higher abundance of m5C modification in resistant CC samples, and NSUN6 was the essential m5C-regulating gene concerning radiosensitivity. NSUN6 overexpression was clinically correlated with radioresistance and poor prognosis in cervical cancer. Functionally, higher NSUN6 expression was associated with radioresistance in the 3D PDO model of cervical cancer. Moreover, silencing NSUN6 increased CC radiosensitivity in vivo and in vitro. Mechanistically, NDRG1 was one of the downstream target genes of NSUN6 identified by integrated m5C-seq, mRNA-seq, and functional validation. NSUN6 promoted the m5C modification of NDRG1 mRNA, and the m5C reader ALYREF bound explicitly to the m5C-labeled NDRG1 mRNA and enhanced NDRG1 mRNA stability. NDRG1 overexpression promoted homologous recombination-mediated DNA repair, which in turn led to radioresistance in cervical cancer. CONCLUSIONS: Aberrant m5C hypermethylation and NSUN6 overexpression drive resistance to radiotherapy in cervical cancer. Elevated NSUN6 expression promotes radioresistance in cervical cancer by activating the NSUN6/ALYREF-m5C-NDRG1 pathway. The low expression of NSUN6 in cervical cancer indicates sensitivity to radiotherapy and a better prognosis.

m6A- and m5C- modified lncRNAs orchestrate the prognosis in cutaneous melanoma and m6A- modified LINC00893 regulates cutaneous melanoma cell metastasis.

BACKGROUND: As the most important modifications on the RNA level, N6-methyladenosine (m6A-) and 5-methylcytosine (m5C-) modification could have a direct influence on the RNAs. Long non-coding RNAs (lncRNAs) could also be modified by methylcytosine modification. Compared with mRNAs, the function of lncRNAs could be more potent to some extent in biological processes like tumorigenesis. Until now, rare reports have been done associated with cutaneous melanoma. Herein, we wonder if the m6A- and m5C- modified lncRNAs could influence the immune landscape and prognosis in melanoma, and we also want to find some lncRNAs which could directly affect the malignant behaviors of melanoma. METHODS: Systematically, we explored the expression pattern of m6A- and m5C- modified lncRNAs in melanoma from datasets including UCSC Xena and NCBI GEO, and the prognostic lncRNAs were selected. Then, according to the expression pattern of lncRNAs, melanoma samples from these datasets were divided into several subtypes. Prognostic model, nomogram survival model, drug sensitivity, GO, and KEGG pathway analysis were performed. Furthermore, among several selected lncRNAs, we identified one lncRNA named LINC00893 and investigated its expression pattern and its biological function in melanoma cell lines. RESULTS: We identified 27 m6A- and m5C- related lncRNAs which were significantly associated with survival, and we made a subtype analysis of melanoma samples based on these 27 lncRNAs. Among the two subtypes, we found differences of immune cells infiltration between these two subtypes. Then, LASSO algorithm was used to screen the optimized lncRNAs combination including ZNF252P-AS1, MIAT, FAM13A-AS1, LINC-PINT, LINC00893, AGAP2-AS1, OIP5-AS1, and SEMA6A-AS1. We also found that there was a significant correlation between the different risk groups predicted based on RS model and the actual prognosis. The nomogram survival model based on independent survival prognostic factors was also constructed. Besides, sensitivity to chemotherapeutic agents, GO and KEGG analysis were performed. In different risk groups, a total of 14 drug molecules with different distributions were obtained, which included AZD6482, AZD7762, AZD8055, camptothecin, dasatinib, erlotinib, gefitinib, gemcitabine, GSK269962A, nilotinib, rapamycin, and sorafenib. A total of 55 significantly related biological processes and 17 KEGG signaling pathways were screened. At last, we noticed that LINC00893 had a relatively lower expression in melanoma tissue and cell lines compared with adjacent tissues and epidermal melanocyte, and down-regulation of LINC00893 could promote the malignant behavior of melanoma cells in A875 and MV3. In these two melanoma cell lines, down-regulation of m6A-related molecules like YTHDF3 and METTL3 could promote the expression of LINC00893. CONCLUSION: We made an analysis of m6A- and m5C- related lncRNAs in melanoma samples and a prediction of these lncRNAs' role in prognosis, tumor microenvironment, immune infiltration, and clinicopathological features. We also found that LINC00893, which is potentially regulated by m6A modification, could serve as a tumor-suppressor in melanoma and play an inhibitory role in melanoma metastasis.

Gene body DNA hydroxymethylation restricts the magnitude of transcriptional changes during aging.

DNA hydroxymethylation (5hmC), the most abundant oxidative derivative of DNA methylation, is typically enriched at enhancers and gene bodies of transcriptionally active and tissue-specific genes. Although aberrant genomic 5hmC has been implicated in age-related diseases, its functional role in aging remains unknown. Here, using mouse liver and cerebellum as model organs, we show that 5hmC accumulates in gene bodies associated with tissue-specific function and restricts the magnitude of gene expression changes with age. Mechanistically, 5hmC decreases the binding of splicing associated factors and correlates with age-related alternative splicing events. We found that various age-related contexts, such as prolonged quiescence and senescence, drive the accumulation of 5hmC with age. We provide evidence that this age-related transcriptionally restrictive function is conserved in mouse and human tissues. Our findings reveal that 5hmC regulates tissue-specific function and may play a role in longevity.

TET1 displays catalytic and non-catalytic functions in the adult mouse cortex.

TET1/2/3 dioxygenases iteratively demethylate 5-methylcytosine, beginning with the formation of 5-hydroxymethylcytosine (5hmC). The post-mitotic brain maintains higher levels of 5hmC than most peripheral tissues, and TET1 ablation studies have underscored the critical role of TET1 in brain physiology. However, deletion of Tet1 precludes the disentangling of the catalytic and non-catalytic functions of TET1. Here, we dissect these functions of TET1 by comparing adult cortex of Tet1 wildtype (Tet1 WT), a novel Tet1 catalytically dead mutant (Tet1 HxD), and Tet1 knockout (Tet1 KO) mice. Using DNA methylation array, we uncover that Tet1 HxD and KO mutations perturb the methylation status of distinct subsets of CpG sites. Gene ontology (GO) analysis on specific differential 5hmC regions indicates that TET1's catalytic activity is linked to neuronal-specific functions. RNA-Seq further shows that Tet1 mutations predominantly impact the genes that are associated with alternative splicing. Lastly, we performed High-performance Liquid Chromatography Mass-Spectrometry lipidomics on WT and mutant cortices and uncover accumulation of lysophospholipids lysophosphatidylethanolamine and lysophosphatidylcholine in Tet1 HxD cortex. In summary, we show that Tet1 HxD does not completely phenocopy Tet1 KO, providing evidence that TET1 modulates distinct cortical functions through its catalytic and non-catalytic roles.

DNA methylation enables recurrent endogenization of giant viruses in an animal relative.

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

TET2 regulates extranodal NK/T cell lymphoma progression through regulation of DNA methylation.

The biological role of Ten-11 translocation 2 (TET2) and the conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) in the development of extra-nodal natural killer/T-cell lymphoma (ENKTL) remains unclear. The level of 5mC and 5hmC was detected in 112 cases of ENKTL tissue specimens by immunohistochemical (IHC) staining. Subsequently, TET2 knockdown and the overexpression cell models were constructed in ENKTL cell lines. Biochemical analyses were used to assess proliferation, apoptosis, cell cycle and monoclonal formation in cells treated or untreated with L-Ascorbic acid sodium salt (LAASS). Dot-Blots were used to detect levels of genome 5mC and 5hmC. Additionally, the ILLUMINA 850k methylation chip was used to analyze the changes of TET2 regulatory genes. RNA-Seq was used to profile differentially expressed genes regulated by TET2. The global level of 5hmC was significantly decreased, while 5mC was highly expressed in ENKTL tissue. TET2 protein expression was negatively correlated with the ratio of 5mC/5hmC (p < 0.0001). The 5mC/5hmC status were related to the site of disease, clinical stage, PINK score and Ki-67 index, as well as the 5-year OS. TET2 knockdown prolonged the DNA synthesis period, increased the cloning ability of tumor cells, increased the level of 5mC and decreased the level of 5hmC in ENKTL cells. While overexpression of TET2 presented the opposite effect. Furthermore, treatment of ENKTL cells with LAASS significantly induced ENKTL cell apoptosis. These results suggest that TET2 plays an important role in ENKTL development via regulation of 5mC and 5hmC and may serve as a novel therapeutic target for ENKTL.

5-mC methylation study of sORFs in 3'UTR of transcription factor JUNGBRUNNEN 1-like during leaf rust pathogenesis in wheat.

BACKGROUND: JUB1, a NAC domain containing hydrogen peroxide-induced transcription factor, plays a critical role in plant immunity. Little is known about how JUB1 responds to leaf rust disease in wheat. Recent discoveries in genomics have also unveiled a multitude of sORFs often assumed to be non-functional, to argue for the necessity of including them as potential regulatory players of translation. However, whether methylation on sORFs spanning the 3'UTR of regulatory genes like JUB1 modulate gene expression, remains unclear. METHODS AND RESULTS: In this study, we identified the methylation states of two sORFs in 3'UTR of a homologous gene of JUB1 in wheat, TaJUB1-L, at cytosine residues in CpG, CHH and CHG sites at different time points of disease progression in two near-isogenic lines of wheat (HD2329), with and without Lr24 gene during leaf rust pathogenesis. Here, we report a significant demethylation of the CpG dinucleotides occurring in the sORFs of the 3'UTR in the resistant isolines after 24 h post-infection. Also, the up-regulated gene expression observed through RT-qPCR was directly proportional to the demethylation of the CpG sites in the sORFs. CONCLUSIONS: Our findings indicate that TaJUB1-L might be a positive regulator in providing tolerance during leaf rust pathogenesis and cytosine methylation at 3'UTR might act as a switch for its expression control. These results enrich the potential benefit of conventional methylation assay techniques for unraveling the unexplored enigma in epigenetics during plant-pathogen interaction in a cost-effective and confidentially conclusive manner.

An Overview of Global, Local, and Base-Resolution Methods for the Detection of 5-Hydroxymethylcytosine in Genomic DNA.

The discovery of 5-hydroxymethylcytosine (5hmC) as a common DNA modification in mammalian genomes has ushered in new areas of inquiry regarding the dynamic epigenome. The balance between 5hmC and its precursor, 5-methylcytosine (5mC), has emerged as a determinant of key processes including cell fate specification, and alterations involving these bases have been implicated in the pathogenesis of various diseases. The identification of 5hmC separately from 5mC initially posed a challenge given that legacy epigenetic sequencing technologies could not discriminate between these two most abundant modifications, a significant blind spot considering their potentially functionally opposing roles. The growing interest in 5hmC, as well as in the Ten-Eleven Translocation (TET) family enzymes that catalyze its generation and further oxidation to 5-formylcytosine (5fC) and 5-carboxycytosine (5caC), has spurred the development of versatile methods for 5hmC detection. These methods enable the quantification and localization of 5hmC in diverse biological samples and, in some cases, at the resolution of individual nucleotides. However, navigating this growing toolbox of methods for 5hmC detection can be challenging. Here, we detail existing and emerging methods for the detection, quantification, and localization of 5hmC at global, locus-specific, and base resolution levels. These methods are discussed in the context of their advantages and limitations, with the goal of providing a framework to help guide researchers in choosing the level of resolution and the associated method that could be most suitable for specific needs.

Whole-Genome Bisulfite Sequencing Protocol for the Analysis of Genome-Wide DNA Methylation and Hydroxymethylation Patterns at Single-Nucleotide Resolution.

The analysis of genome-wide epigenomic alterations including DNA methylation and hydroxymethylation has become a subject of intensive research for many biological and clinical questions. DNA methylation analysis bears the particular promise to supplement or replace biochemical and imaging-based tests for the next generation of personalized medicine. Whole-genome bisulfite sequencing (WGBS) using next-generation sequencing technologies is currently considered the gold standard for a comprehensive and quantitative analysis of DNA methylation throughout the genome. However, bisulfite conversion does not allow distinguishing between cytosine methylation and hydroxymethylation requiring an additional chemical or enzymatic step to identify hydroxymethylated cytosines. Here, we provide a detailed protocol based on a commercial kit for the preparation of sequencing libraries for the comprehensive whole-genome analysis of DNA methylation and/or hydroxymethylation. The protocol is based on the construction of sequencing libraries from limited amounts of input DNA by ligation of methylated adaptors to the fragmented DNA prior to bisulfite conversion. For analyses requiring a quantitative distinction between 5-methylcytosine and 5-hydroxymethylcytosines levels, an oxidation step is included in the same workflow to perform oxidative bisulfite sequencing (OxBs-Seq). In this case, two sequencing libraries will be generated and sequenced: a classic methylome following bisulfite conversion and analyzing modified cytosines (not distinguishing between methylated and hydroxymethylated cytosines) and a methylome analyzing only methylated cytosines, respectively. Hydroxymethylation levels are deduced from the differences between the two reactions. We also provide a step-by-step description of the data analysis using publicly available bioinformatic tools. The described protocol has been successfully applied to different human and plant samples and yields robust and reproducible results.

Substrate diversity of NSUN enzymes and links of 5-methylcytosine to mRNA translation and turnover.

Maps of the RNA modification 5-methylcytosine (m5C) often diverge markedly not only because of differences in detection methods, data depand analysis pipelines but also biological factors. We re-analysed bisulfite RNA sequencing datasets from five human cell lines and seven tissues using a coherent m5C site calling pipeline. With the resulting union list of 6,393 m5C sites, we studied site distribution, enzymology, interaction with RNA-binding proteins and molecular function. We confirmed tRNA:m5C methyltransferases NSUN2 and NSUN6 as the main mRNA m5C "writers," but further showed that the rRNA:m5C methyltransferase NSUN5 can also modify mRNA. Each enzyme recognises mRNA features that strongly resemble their canonical substrates. By analysing proximity between mRNA m5C sites and footprints of RNA-binding proteins, we identified new candidates for functional interactions, including the RNA helicases DDX3X, involved in mRNA translation, and UPF1, an mRNA decay factor. We found that lack of NSUN2 in HeLa cells affected both steady-state levels of, and UPF1-binding to, target mRNAs. Our studies emphasise the emerging diversity of m5C writers and readers and their effect on mRNA function.

Epigenetic markers of tooth eruption - DNA methylation and histone acetylation.

The present study aimed to evaluate whether epigenetic markers are expressed in the dental follicles surrounding ectopically erupting teeth. Twenty-one dental follicles were collected in 20 adolescent children through surgical exposure of ectopic teeth. The epigenetic modifications of DNA methylation and histone acetylation were evaluated by immunohistochemistry. The results showed cells positive for DNA-methyltransferase 1 (DNMT1), DNA methyltransferase 3 beta (DNMT3B), ten-eleven translocation-2 (TET2), acetyl-histone H3 (AcH3), acetyl-histone H4 (AcH4), 5-methylcytosine (5mC), and 5-hydroxymethylcytosine (5hmC) were present in all the samples. The levels of epigenetic markers representing active chromatin (5hmC, AcH3, AcH4, and TET2) were statistically significantly higher than those of markers representing inactive chromatin (5mC, DNMT3B, DNMT1). In conclusion, follicles in ectopic teeth display major epigenetic modifications. In the follicles, epigenetic markers associated with the activation of bone-related genes are more abundant than markers associated with the inactivation of bone-related genes.

Global DNA methylation is not elevated in blood samples from Machado-Joseph disease mutation carriers.

Machado-Joseph disease (MJD) is an autosomal dominant spinocerebellar ataxia (SCA) caused by a polyglutamine expansion in the ataxin-3 protein, which initiates a cascade of pathogenic events, including transcriptional dysregulation. Genotype-phenotype correlations in MJD are incomplete, suggesting an influence of additional factors, such as epigenetic modifications, underlying the MJD pathogenesis. DNA methylation is known to impact the pathophysiology of neurodegenerative disorders through gene expression regulation and increased methylation has been reported for other SCAs. In this work we aimed to analyse global methylation in MJD carriers. Global 5-mC levels were quantified in blood samples of 33 MJD mutation carriers (patients and preclinical subjects) and 33 healthy controls, matched by age, sex, and smoking status. For a subset of 16 MJD subjects, a pilot follow-up analysis with two time points was also conducted. No differences were found in median global 5-mC levels between MJD mutation carriers and controls and no correlations between methylation levels and clinical or genetic variables were detected. Also, no alterations in global 5-mC levels were observed over time. Our findings do not support an increase in global blood methylation levels associated with MJD.

Smoking-Induced DNA Hydroxymethylation Signature Is Less Pronounced than True DNA Methylation: The Population-Based KORA Fit Cohort.

Despite extensive research on 5-methylcytosine (5mC) in relation to smoking, there has been limited exploration into the interaction between smoking and 5-hydroxymethylcytosine (5hmC). In this study, total DNA methylation (5mC+5hmC), true DNA methylation (5mC) and hydroxymethylation (5hmC) levels were profiled utilizing conventional bisulphite (BS) and oxidative bisulphite (oxBS) treatment, measured with the Illumina Infinium Methylation EPIC BeadChip. An epigenome-wide association study (EWAS) of 5mC+5hmC methylation revealed a total of 38,575 differentially methylated positions (DMPs) and 2023 differentially methylated regions (DMRs) associated with current smoking, along with 82 DMPs and 76 DMRs associated with former smoking (FDR-adjusted p < 0.05). Additionally, a focused examination of 5mC identified 33 DMPs linked to current smoking and 1 DMP associated with former smoking (FDR-adjusted p < 0.05). In the 5hmC category, eight DMPs related to current smoking and two DMPs tied to former smoking were identified, each meeting a suggestive threshold (p < 1 × 10-5). The substantial number of recognized DMPs, including 5mC+5hmC (7069/38,575, 2/82), 5mC (0/33, 1/1), and 5hmC (2/8, 0/2), have not been previously reported. Our findings corroborated previously established methylation positions and revealed novel candidates linked to tobacco smoking. Moreover, the identification of hydroxymethylated CpG sites with suggestive links provides avenues for future research.

Exploration of the prognostic value of methylation regulators related to m5C in papillary thyroid carcinoma.

The incidence of papillary thyroid carcinoma (PTC) has increased significantly in recent years, and for patients with metastatic and recurrent PTC, the options for treatment currently available are insufficient. To date, the exact molecular mechanism underlying PTC is still not fully understood. 5-Methylcytosine (m5C) RNA methylation is associated with the prognosis of a variety of tumors. However, the molecular mechanisms and biomarkers associated with m5C in the diagnosis, treatment, and prognosis of this disease have not been fully elucidated. Ten m5C regulators with significantly different expression levels were included in this study. Immune infiltration analysis revealed significant negative correlations between most of these regulators and regulatory T cells. TRDMT1, NSUN5, and NSUN6 had high weights and strong correlations in the protein-protein interaction network. Using gene ontology, Kyoto Encyclopedia of Genes and Genomes, and gene set enrichment analysis, 1489 differentially expressed genes were screened from The Cancer Genome Atlas messenger RNA matrix, indicating that these differentially expressed genes were significantly enriched in various pathways and functions related to cancers. Four m5C regulators, NSUN2, NSUN4, NSUN6, and DNMT3B, were screened as prognostic markers by least absolute shrinkage and selection operator regression analysis, and NSUN2 and NSUN6 were identified as risk factors for poor prognosis. We found that the prognostic prediction model constructed using the m5C regulators NSUN2, NSUN4, NSUN6, and DNMT3B showed good prognostic prediction ability and diagnostic ability. This model was applied to predict the survival probability of patients with PTC, the prediction ability of 5-year survival was the best. The multi-factor prognostic prediction model combined with the tumor node metastasis stage and risk score grouping showed better prognostic predictive power.

Correlation between sperm DNA fragmentation and methylation in male partners of couples with idiopathic recurrent pregnancy loss.

With ∼50% recurrent pregnancy loss cases being termed idiopathic (iRPL), understanding of contribution of male factors to iRPL is still lacking. Higher prevalence of sperm DNA fragmentation index (DFI) and lower sperm 5-methylcytosine (5-mC) levels have been previously reported in male partners of iRPL couples and shed light on importance of the male gamete in maintenance of a successful pregnancy. The present study aimed to determine the serum sex steroid hormone levels, sperm DFI and 5-mC and correlation between them in male partners of fertile and iRPL couples. Further, correlation between sperm DFI and 5-mC with semen parameters and paternal age in both groups were determined. 36 male partners of fertile couples and 45 male partners of women experiencing iRPL were enrolled for this study and semen and blood samples were collected. Serum testosterone and estradiol levels were measured by ELISA; sperm DFI and global 5-mC were determined by TUNEL assay and ELISA respectively. Significantly higher serum testosterone levels were noted in the iRPL group (p = 0.028). Incidence of sperm DNA fragmentation was found to be higher in the iRPL study group but with no significance difference. No significant differences in sperm 5-mC values were noted. Upon correlation analysis within both groups, strong significant negative correlation of sperm DFI % and 5-mC % was observed in the control group (p < 0.001) but not the iRPL group (p = 0.249). Hence, we infer that with lower 5-mC levels in sperm genome, there is a higher incidence of sperm DFI in fertile men. However, this trend is not noted in men of iRPL group which could possibly be due to other underlying epigenetic alterations in genomic regions probably unsusceptible to fragmentation. On the other hand, no significant correlations of semen parameters, testosterone, estradiol and paternal age with sperm DFI and 5-mC were noted in both groups.

Obesity-driven mitochondrial dysfunction in human adipose tissue-derived mesenchymal stem/stromal cells involves epigenetic changes.

Obesity exacerbates tissue degeneration and compromises the integrity and reparative potential of mesenchymal stem/stromal cells (MSCs), but the underlying mechanisms have not been sufficiently elucidated. Mitochondria modulate the viability, plasticity, proliferative capacity, and differentiation potential of MSCs. We hypothesized that alterations in the 5-hydroxymethylcytosine (5hmC) profile of mitochondria-related genes may mediate obesity-driven dysfunction of human adipose-derived MSCs. MSCs were harvested from abdominal subcutaneous fat of obese and age/sex-matched non-obese subjects (n = 5 each). The 5hmC profile and expression of nuclear-encoded mitochondrial genes were examined by hydroxymethylated DNA immunoprecipitation sequencing (h MeDIP-seq) and mRNA-seq, respectively. MSC mitochondrial structure (electron microscopy) and function, metabolomics, proliferation, and neurogenic differentiation were evaluated in vitro, before and after epigenetic modulation. hMeDIP-seq identified 99 peaks of hyper-hydroxymethylation and 150 peaks of hypo-hydroxymethylation in nuclear-encoded mitochondrial genes from Obese- versus Non-obese-MSCs. Integrated hMeDIP-seq/mRNA-seq analysis identified a select group of overlapping (altered levels of both 5hmC and mRNA) nuclear-encoded mitochondrial genes involved in ATP production, redox activity, cell proliferation, migration, fatty acid metabolism, and neuronal development. Furthermore, Obese-MSCs exhibited decreased mitochondrial matrix density, membrane potential, and levels of fatty acid metabolites, increased superoxide production, and impaired neuronal differentiation, which improved with epigenetic modulation. Obesity elicits epigenetic changes in mitochondria-related genes in human adipose-derived MSCs, accompanied by structural and functional changes in their mitochondria and impaired fatty acid metabolism and neurogenic differentiation capacity. These observations may assist in developing novel therapies to preserve the potential of MSCs for tissue repair and regeneration in obese individuals.

Predicting gene expression state and prioritizing putative enhancers using 5hmC signal.

BACKGROUND: Like its parent base 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) is a direct epigenetic modification of cytosines in the context of CpG dinucleotides. 5hmC is the most abundant oxidized form of 5mC, generated through the action of TET dioxygenases at gene bodies of actively-transcribed genes and at active or lineage-specific enhancers. Although such enrichments are reported for 5hmC, to date, predictive models of gene expression state or putative regulatory regions for genes using 5hmC have not been developed. RESULTS: Here, by using only 5hmC enrichment in genic regions and their vicinity, we develop neural network models that predict gene expression state across 49 cell types. We show that our deep neural network models distinguish high vs low expression state utilizing only 5hmC levels and these predictive models generalize to unseen cell types. Further, in order to leverage 5hmC signal in distal enhancers for expression prediction, we employ an Activity-by-Contact model and also develop a graph convolutional neural network model with both utilizing Hi-C data and 5hmC enrichment to prioritize enhancer-promoter links. These approaches identify known and novel putative enhancers for key genes in multiple immune cell subsets. CONCLUSIONS: Our work highlights the importance of 5hmC in gene regulation through proximal and distal mechanisms and provides a framework to link it to genome function. With the recent advances in 6-letter DNA sequencing by short and long-read techniques, profiling of 5mC and 5hmC may be done routinely in the near future, hence, providing a broad range of applications for the methods developed here.

m5C methylated lncRncr3-MeCP2 interaction restricts miR124a-initiated neurogenesis.

Coordination of neuronal differentiation with expansion of the neuroepithelial/neural progenitor cell (NEPC/NPC) pool is essential in early brain development. Our in vitro and in vivo studies identify independent and opposing roles for two neural-specific and differentially expressed non-coding RNAs derived from the same locus: the evolutionarily conserved lncRNA Rncr3 and the embedded microRNA miR124a-1. Rncr3 regulates NEPC/NPC proliferation and controls the biogenesis of miR124a, which determines neuronal differentiation. Rncr3 conserved exons 2/3 are cytosine methylated and bound by methyl-CpG binding protein MeCP2, which restricts expression of miR124a embedded in exon 4 to prevent premature neuronal differentiation, and to orchestrate proper brain growth. MeCP2 directly binds cytosine-methylated Rncr3 through previously unrecognized lysine residues and suppresses miR124a processing by recruiting PTBP1 to block access of DROSHA-DGCR8. Thus, miRNA processing is controlled by lncRNA m5C methylation along with the defined m5C epitranscriptomic RNA reader protein MeCP2 to coordinate brain development.

Bioinspired Heterocoordination in Adaptable Cobalt Metal-Organic Framework for DNA Epigenetic Modification Detection.

Metal-organic frameworks (MOFs) show unique advantages in simulating the dynamics and fidelity of natural coordination. Inspired by zinc finger protein, a second linker was introduced to affect the homogeneous MOF system and thus facilitate the emergence of diverse functionalities. Under the systematic identification of 12 MOF species (i.e., metal ions, linkers) and 6 second linkers (trigger), a dissipative system consisting of Co-BDC-NO2 and o-phenylenediamine (oPD) was screened out, which can rapidly and in situ generate a high photothermal complex (η = 36.9%). Meanwhile, both the carboxylation of epigenetic modifications and metal ion (Fe3+, Ni2+, Cu2+, Zn2+, Co2+ and Mn2+) screening were utilized to improve the local coordination environment so that the adaptable Co-MOF growth on the DNA strand was realized. Thus, epigenetic modification information on DNA was converted to an amplified metal ion signal, and then oPD was further introduced to generate bimodal dissipative signals by which a simple, high-sensitivity detection strategy of 5-hydroxymethylcytosine (LOD = 0.02%) and 5-formylcytosine (LOD = 0.025‰) was developed. The strategy provides one low-cost method (< 0.01 $/sample) for quantifying global epigenetic modifications, which greatly promotes epigenetic modification-based early disease diagnosis. This work also proposes a general heterocoordination design concept for molecular recognition and signal transduction, opening a new MOF-based sensing paradigm.