Endothelial deletion of EPH receptor A4 alters single-cell profile and Tie2/Akap12 signaling to preserve blood-brain barrier integrity.

Neurobiological consequences of traumatic brain injury (TBI) result from a complex interplay of secondary injury responses and sequela that mediates chronic disability. Endothelial cells are important regulators of the cerebrovascular response to TBI. Our work demonstrates that genetic deletion of endothelial cell (EC)-specific EPH receptor A4 (EphA4) using conditional EphA4f/f/Tie2-Cre and EphA4f/f/VE-Cadherin-CreERT2 knockout (KO) mice promotes blood-brain barrier (BBB) integrity and tissue protection, which correlates with improved motor function and cerebral blood flow recovery following controlled cortical impact (CCI) injury. scRNAseq of capillary-derived KO ECs showed increased differential gene expression of BBB-related junctional and actin cytoskeletal regulators, namely, A-kinase anchor protein 12, Akap12, whose presence at Tie2 clustering domains is enhanced in KO microvessels. Transcript and protein analysis of CCI-injured whole cortical tissue or cortical-derived ECs suggests that EphA4 limits the expression of Cldn5, Akt, and Akap12 and promotes Ang2. Blocking Tie2 using sTie2-Fc attenuated protection and reversed Akap12 mRNA and protein levels cortical-derived ECs. Direct stimulation of Tie2 using Vasculotide, angiopoietin-1 memetic peptide, phenocopied the neuroprotection. Finally, we report a noteworthy rise in soluble Ang2 in the sera of individuals with acute TBI, highlighting its promising role as a vascular biomarker for early detection of BBB disruption. These findings describe a contribution of the axon guidance molecule, EphA4, in mediating TBI microvascular dysfunction through negative regulation of Tie2/Akap12 signaling.

Phosphorylation of RPT6 Controls Its Ability to Bind DNA and Regulate Gene Expression in the Hippocampus of Male Rats during Memory Formation.

Memory formation requires coordinated control of gene expression, protein synthesis, and ubiquitin-proteasome system (UPS)-mediated protein degradation. The catalytic component of the UPS, the 26S proteasome, contains a 20S catalytic core surrounded by two 19S regulatory caps, and phosphorylation of the 19S cap regulatory subunit RPT6 at serine 120 (pRPT6-S120) has been widely implicated in controlling activity-dependent increases in proteasome activity. Recently, RPT6 was also shown to act outside the proteasome where it has a transcription factor-like role in the hippocampus during memory formation. However, little is known about the proteasome-independent function of "free" RPT6 in the brain or during memory formation and whether phosphorylation of S120 is required for this transcriptional control function. Here, we used RNA-sequencing along with novel genetic approaches and biochemical, molecular, and behavioral assays to test the hypothesis that pRPT6-S120 functions independently of the proteasome to bind DNA and regulate gene expression during memory formation. RNA-sequencing following siRNA-mediated knockdown of free RPT6 revealed 46 gene targets in the dorsal hippocampus of male rats following fear conditioning, where RPT6 was involved in transcriptional activation and repression. Through CRISPR-dCas9-mediated artificial placement of RPT6 at a target gene, we found that RPT6 DNA binding alone may be important for altering gene expression following learning. Further, CRISPR-dCas13-mediated conversion of S120 to glycine on RPT6 revealed that phosphorylation at S120 is necessary for RPT6 to bind DNA and properly regulate transcription during memory formation. Together, we reveal a novel function for phosphorylation of RPT6 in controlling gene transcription during memory formation.

Monoubiquitination of histone H2B is a crucial regulator of the transcriptome during memory formation.

Posttranslational modification of histone proteins is critical for memory formation. Recently, we showed that monoubiquitination of histone H2B at lysine 120 (H2Bub) is critical for memory formation in the hippocampus. However, the transcriptome controlled by H2Bub remains unknown. Here, we found that fear conditioning in male rats increased or decreased the expression of 86 genes in the hippocampus but, surprisingly, siRNA-mediated knockdown of the H2Bub ligase, Rnf20, abolished changes in all but one of these genes. These findings suggest that monoubiquitination of histone H2B is a crucial regulator of the transcriptome during memory formation.

Impairment of DNA Methylation Maintenance Is the Main Cause of Global Demethylation in Naive Embryonic Stem Cells.

Global demethylation is part of a conserved program of epigenetic reprogramming to naive pluripotency. The transition from primed hypermethylated embryonic stem cells (ESCs) to naive hypomethylated ones (serum-to-2i) is a valuable model system for epigenetic reprogramming. We present a mathematical model, which accurately predicts global DNA demethylation kinetics. Experimentally, we show that the main drivers of global demethylation are neither active mechanisms (Aicda, Tdg, and Tet1-3) nor the reduction of de novo methylation. UHRF1 protein, the essential targeting factor for DNMT1, is reduced upon transition to 2i, and so is recruitment of the maintenance methylation machinery to replication foci. Concurrently, there is global loss of H3K9me2, which is needed for chromatin binding of UHRF1. These mechanisms synergistically enforce global DNA hypomethylation in a replication-coupled fashion. Our observations establish the molecular mechanism for global demethylation in naive ESCs, which has key parallels with those operating in primordial germ cells and early embryos.

Dynamics of RNA m5C modification during brain development.

Post-transcriptional RNA modifications have been recognized as key regulators of neuronal differentiation and synapse development in the mammalian brain. While distinct sets of 5-methylcytosine (m5C) modified mRNAs have been detected in neuronal cells and brain tissues, no study has been performed to characterize methylated mRNA profiles in the developing brain. Here, together with regular RNA-seq, we performed transcriptome-wide bisulfite sequencing to compare RNA cytosine methylation patterns in neural stem cells (NSCs), cortical neuronal cultures, and brain tissues at three postnatal stages. Among 501 m5C sites identified, approximately 6% are consistently methylated across all five conditions. Compared to m5C sites identified in NSCs, 96% of them were hypermethylated in neurons and enriched for genes involved in positive transcriptional regulation and axon extension. In addition, brains at the early postnatal stage demonstrated substantial changes in both RNA cytosine methylation and gene expression of RNA cytosine methylation readers, writers, and erasers. Furthermore, differentially methylated transcripts were significantly enriched for genes regulating synaptic plasticity. Altogether, this study provides a brain epitranscriptomic dataset as a new resource and lays the foundation for further investigations into the role of RNA cytosine methylation during brain development.

Folate regulates RNA m5C modification and translation in neural stem cells.

BACKGROUND: Folate is an essential B-group vitamin and a key methyl donor with important biological functions including DNA methylation regulation. Normal neurodevelopment and physiology are sensitive to the cellular folate levels. Either deficiency or excess of folate may lead to neurological disorders. Recently, folate has been linked to tRNA cytosine-5 methylation (m5C) and translation in mammalian mitochondria. However, the influence of folate intake on neuronal mRNA m5C modification and translation remains largely unknown. Here, we provide transcriptome-wide landscapes of m5C modification in poly(A)-enriched RNAs together with mRNA transcription and translation profiles for mouse neural stem cells (NSCs) cultured in three different concentrations of folate. RESULTS: NSCs cultured in three different concentrations of folate showed distinct mRNA methylation profiles. Despite uncovering only a few differentially expressed genes, hundreds of differentially translated genes were identified in NSCs with folate deficiency or supplementation. The differentially translated genes induced by low folate are associated with cytoplasmic translation and mitochondrial function, while the differentially translated genes induced by high folate are associated with increased neural stem cell proliferation. Interestingly, compared to total mRNAs, polysome mRNAs contained high levels of m5C. Furthermore, an integrative analysis indicated a transcript-specific relationship between RNA m5C methylation and mRNA translation efficiency. CONCLUSIONS: Altogether, our study reports a transcriptome-wide influence of folate on mRNA m5C methylation and translation in NSCs and reveals a potential link between mRNA m5C methylation and mRNA translation.

Advances in methods and software for RNA cytosine methylation analysis.

Our understanding of RNA modifications has been growing rapidly over the last decade. Epitranscriptomics has recently emerged as an exciting, new field for understanding the fundamental mechanisms underlying RNA modifications and their impact on gene expression. Among the over one hundred different kinds of RNA modifications, cytosine methylation in mRNA (5-mrC) is now recognized as an important epigenetic mark that modulates mRNA transportation, translation, and stability at the post-transcriptional level. Across plant and animal species, recent studies have revealed the roles of mRNA cytosine methylation in several fundamental biological processes. In mammals, genome-wide profiling has determined thousands of mRNA transcripts carrying the 5-mrC modification in a tissue specific manner. Here, we summarize the experimental techniques that were exploited to determine 5-mrC in mRNA and the computational procedures implemented for RNA bisulfite sequencing data analysis.

EGR2 is elevated and positively regulates inflammatory IFNγ production in lupus CD4+ T cells.

BACKGROUND: Recent studies have shown that early growth response 2 (EGR2) is highly induced in activated T cells and regulates T cell functions. In normal C57BL/6 (B6) mice, deletion of EGR2 in lymphocytes results in the development of lupus-like systemic autoimmune disease, which implies indirectly an autoimmune protective role of EGR2. Conversely, increased EGR2 gene expression is suggested to link with high risk of human lupus. In the present studies we sought to clarify the expression and inflammation regulatory role of EGR2 in murine lupus T cells directly. RESULTS: We performed RT-qPCR analysis and found a significant increase of EGR2 mRNA expression in human lupus PBMCs and in CD4+ T cells from three different murine lupus models including MRL-lpr, B6-lpr, and B6.sle123 mice at diseased stage when compared to age-matched control MRL or B6 mice. By performing intracellular flow cytometry analysis, we found that EGR2 protein expression was significantly increased in resting lupus (either MRL-lpr or B6.sle123) CD4+ T cells when compared to CD4+ T cells from their respective non-autoimmune controls. However, there was no difference of EGR2 protein expression in anti-CD3 and anti-CD28 stimulated control and lupus CD4+ T cells since there was a stronger induction of EGR2 in activated control CD4+ T cells. EGR2 expression was significantly increased in MRL-lpr mice at an age when lupus is manifested. To understand further the function of elevated EGR2 in lupus CD4+ T cells, we inhibited EGR2 with a specific siRNA in vitro in splenocytes from MRL-lpr and control MRL mice at 15 weeks-of-age. We found that EGR2 inhibition significantly reduced IFNγ production in PMA and ionomycin activated MRL-lpr lupus CD4+ T cells, but not control MRL CD4+ T cells. We also found that inhibition of EGR2 in vitro suppressed the Th1 differentiation in both MRL and MRL-lpr naïve CD4+ T cells. CONCLUSIONS: EGR2 is highly upregulated in human and murine lupus cells. Our in vitro data suggest a positive role of EGR2 in the regulation of Th1 differentiation and IFNγ production in lupus effector CD4+ T cells.

Integrative single-cell omics analyses reveal epigenetic heterogeneity in mouse embryonic stem cells.

Embryonic stem cells (ESCs) consist of a population of self-renewing cells displaying extensive phenotypic and functional heterogeneity. Research towards the understanding of the epigenetic mechanisms underlying the heterogeneity among ESCs is still in its initial stage. Key issues, such as how to identify cell-subset specifically methylated loci and how to interpret the biological meanings of methylation variations remain largely unexplored. To fill in the research gap, we implemented a computational pipeline to analyze single-cell methylome and to perform an integrative analysis with single-cell transcriptome data. According to the origins of variation in DNA methylation, we determined the genomic loci associated with allelic-specific methylation or asymmetric DNA methylation, and explored a beta mixture model to infer the genomic loci exhibiting cell-subset specific methylation (CSM). We observed that the putative CSM loci in ESCs are significantly enriched in CpG island (CGI) shelves and regions with histone marks for promoter and enhancer, and the genes hosting putative CSM loci show wide-ranging expression among ESCs. More interestingly, the putative CSM loci may be clustered into co-methylated modules enriching the binding motifs of distinct sets of transcription factors. Taken together, our study provided a novel tool to explore single-cell methylome and transcriptome to reveal the underlying transcriptional regulatory networks associated with epigenetic heterogeneity of ESCs.

The dynamics of DNA methylation fidelity during mouse embryonic stem cell self-renewal and differentiation.

The faithful transmission of DNA methylation patterns through cell divisions is essential for the daughter cells to retain a proper cell identity. To achieve a comprehensive assessment of methylation fidelity, we implemented a genome-scale hairpin bisulfite sequencing approach to generate methylation data for DNA double strands simultaneously. We show here that methylation fidelity increases globally during differentiation of mouse embryonic stem cells (mESCs), and is particularly high in the promoter regions of actively expressed genes and positively correlated with active histone modification marks and binding of transcription factors. The majority of intermediately (40%-60%) methylated CpG dinucleotides are hemi-methylated and have low methylation fidelity, particularly in the differentiating mESCs. While 5-hmC and 5-mC tend to coexist, there is no significant correlation between 5-hmC levels and methylation fidelity. Our findings may shed new light on our understanding of the origins of methylation variations and the mechanisms underlying DNA methylation transmission.

Characterization and machine learning prediction of allele-specific DNA methylation.

A large collection of Single Nucleotide Polymorphisms (SNPs) has been identified in the human genome. Currently, the epigenetic influences of SNPs on their neighboring CpG sites remain elusive. A growing body of evidence suggests that locus-specific information, including genomic features and local epigenetic state, may play important roles in the epigenetic readout of SNPs. In this study, we made use of mouse methylomes with known SNPs to develop statistical models for the prediction of SNP associated allele-specific DNA methylation (ASM). ASM has been classified into parent-of-origin dependent ASM (P-ASM) and sequence-dependent ASM (S-ASM), which comprises scattered-S-ASM (sS-ASM) and clustered-S-ASM (cS-ASM). We found that P-ASM and cS-ASM CpG sites are both enriched in CpG rich regions, promoters and exons, while sS-ASM CpG sites are enriched in simple repeat and regions with high frequent SNP occurrence. Using Lasso-grouped Logistic Regression (LGLR), we selected 21 out of 282 genomic and methylation related features that are powerful in distinguishing cS-ASM CpG sites and trained the classifiers with machine learning techniques. Based on 5-fold cross-validation, the logistic regression classifier was found to be the best for cS-ASM prediction with an ACC of 0.77, an AUC of 0.84 and an MCC of 0.54. Lastly, we applied the logistic regression classifier on human brain methylome and predicted 608 genes associated with cS-ASM. Gene ontology term enrichment analysis indicated that these cS-ASM associated genes are significantly enriched in the category coding for transcripts with alternative splicing forms. In summary, this study provided an analytical procedure for cS-ASM prediction and shed new light on the understanding of different types of ASM events.

Nonparametric Bayesian clustering to detect bipolar methylated genomic loci.

BACKGROUND: With recent development in sequencing technology, a large number of genome-wide DNA methylation studies have generated massive amounts of bisulfite sequencing data. The analysis of DNA methylation patterns helps researchers understand epigenetic regulatory mechanisms. Highly variable methylation patterns reflect stochastic fluctuations in DNA methylation, whereas well-structured methylation patterns imply deterministic methylation events. Among these methylation patterns, bipolar patterns are important as they may originate from allele-specific methylation (ASM) or cell-specific methylation (CSM). RESULTS: Utilizing nonparametric Bayesian clustering followed by hypothesis testing, we have developed a novel statistical approach to identify bipolar methylated genomic regions in bisulfite sequencing data. Simulation studies demonstrate that the proposed method achieves good performance in terms of specificity and sensitivity. We used the method to analyze data from mouse brain and human blood methylomes. The bipolar methylated segments detected are found highly consistent with the differentially methylated regions identified by using purified cell subsets. CONCLUSIONS: Bipolar DNA methylation often indicates epigenetic heterogeneity caused by ASM or CSM. With allele-specific events filtered out or appropriately taken into account, our proposed approach sheds light on the identification of cell-specific genes/pathways under strong epigenetic control in a heterogeneous cell population.

Investigating bisulfite short-read mapping failure with hairpin bisulfite sequencing data.

BACKGROUND: DNA methylation is an important epigenetic mark relevant to normal development and disease genesis. A common approach to characterizing genome-wide DNA methylation is using Next Generation Sequencing technology to sequence bisulfite treated DNA. The short sequence reads are mapped to the reference genome to determine the methylation statuses of Cs. However, despite intense effort, a much smaller proportion of the reads derived from bisulfite treated DNA (usually about 40-80%) can be mapped than regular short reads mapping (> 90%), and it is unclear what factors lead to this low mapping efficiency. RESULTS: To address this issue, we used the hairpin bisulfite sequencing technology to determine sequences of both DNA double strands simultaneously. This enabled the recovery of the original non-bisulfite-converted sequences. We used Bismark for bisulfite read mapping and Bowtie2 for recovered read mapping. We found that recovering the reads improved unique mapping efficiency by 9-10% compared to the bisulfite reads. Such improvement in mapping efficiency is related to sequence entropy. CONCLUSIONS: The hairpin recovery technique improves mapping efficiency, and sequence entropy relates to mapping efficiency.

Deciphering the heterogeneity in DNA methylation patterns during stem cell differentiation and reprogramming.

BACKGROUND: Human induced pluripotent stem cells (iPSCs) have a wide range of applications throughout the fields of basic research, disease modeling and drug screening. Epigenetic instable iPSCs with aberrant DNA methylation may divide and differentiate into cancer cells. Unfortunately, little effort has been taken to compare the epigenetic variation in iPSCs with that in differentiated cells. Here, we developed an analytical procedure to decipher the DNA methylation heterogeneity of mixed cells and further exploited it to quantitatively assess the DNA methylation variation in the methylomes of adipose-derived stem cells (ADS), mature adipocytes differentiated from ADS cells (ADS-adipose) and iPSCs reprogrammed from ADS cells (ADS-iPSCs). RESULTS: We observed that the degree of DNA methylation variation varies across distinct genomic regions with promoter and 5'UTR regions exhibiting low methylation variation and Satellite showing high methylation variation. Compared with differentiated cells, ADS-iPSCs possess globally decreased methylation variation, in particular in repetitive elements. Interestingly, DNA methylation variation decreases in promoter regions during differentiation but increases during reprogramming. Methylation variation in promoter regions is negatively correlated with gene expression. In addition, genes showing a bipolar methylation pattern, with both completely methylated and completely unmethylated reads, are related to the carbohydrate metabolic process, cellular development, cellular growth, proliferation, etc. CONCLUSIONS: This study delivers a way to detect cell-subset specific methylation genes in a mixed cell population and provides a better understanding of methylation dynamics during stem cell differentiation and reprogramming.

DMEAS: DNA methylation entropy analysis software.

SUMMARY: DMEAS is the first user-friendly tool dedicated to analyze the distribution of DNA methylation patterns for the quantification of epigenetic heterogeneity. It supports the analysis of both locus-specific and genome-wide bisulfite sequencing data. DMEAS progressively scans the mapping results of bisulfite sequencing reads to extract DNA methylation patterns for contiguous CpG dinucleotides. It determines the DNA methylation level and calculates methylation entropy for genomic segments to enable the quantitative assessment of DNA methylation variations observed in cell populations. AVAILABILITY AND IMPLEMENTATION: DMEAS program, user guide and all the testing data are freely available from http://sourceforge.net/projects/dmeas/files/

Risk of Excess Maternal Folic Acid Supplementation in Offspring.

Folate, also known as vitamin B9, facilitates the transfer of methyl groups among molecules, which is crucial for amino acid metabolism and nucleotide synthesis. Adequate maternal folate supplementation has been widely acknowledged for its pivotal role in promoting cell proliferation and preventing neural tube defects. However, in the post-fortification era, there has been a rising concern regarding an excess maternal intake of folic acid (FA), the synthetic form of folate. In this review, we focused on recent advancements in understanding the influence of excess maternal FA intake on offspring. For human studies, we summarized findings from clinical trials investigating the effects of periconceptional FA intake on neurodevelopment and molecular-level changes in offspring. For studies using mouse models, we compiled the impact of high maternal FA supplementation on gene expression and behavioral changes in offspring. In summary, excessive maternal folate intake could potentially have adverse effects on offspring. Overall, we highlighted concerns regarding elevated maternal folate status in the population, providing a comprehensive perspective on the potential adverse effects of excessive maternal FA supplementation on offspring.

Histone deacetylase 1 deficiency impairs differentiation and electrophysiological properties of cardiomyocytes derived from induced pluripotent cells.

Epigenetic and chromatin modifications play particularly important roles in embryonic and induced pluripotent stem cells (ESCs and iPSCs) allowing for the cells to both differentiate and dedifferentiate back to a pluripotent state. We analyzed how the loss of a key chromatin-modifying enzyme, histone deacetylase 1 (HDAC1), affects early and cardiovascular differentiation of both ESCs and iPSCs. We also investigated potential differences between these two cell types when differentiation is induced. Our data indicate an essential role for HDAC1 in deacetylating regulatory regions of key pluripotency-associated genes during early differentiation. Although HDAC1 functions primarily as a HDAC, its loss also affects DNA methylation in ESCs and iPSCs both during pluripotency and differentiation. We show that HDAC1 plays a crucial, nonredundant role in cardiomyocyte differentiation and maturation. Our data also elucidate important differences between ESCs and iPSCs, when levels of this enzyme are reduced, that affect their ability to differentiate into functional cardiomyocytes. As varying levels of chromatin-modifying enzymes are likely to exist in patient-derived iPSCs, understanding the molecular circuitry of these enzymes in ESCs and iPSCs is critical for their potential use in cardiovascular therapeutic applications

Methylation alterations of WT1 and homeobox genes in inflamed muscle biopsy samples from patients with untreated juvenile dermatomyositis suggest self-renewal capacity.

OBJECTIVE: To determine the effect of methylation alteration in inflamed muscles from children with juvenile dermatomyositis (DM) and other idiopathic inflammatory myopathies (IIMs). METHODS: Magnetic resonance imaging-directed diagnostic muscle biopsies yielded samples from 20 children with juvenile DM, which were used for genome-wide DNA methylation profiling, as were muscle biopsy samples from 4 healthy controls. Bisulfite treatment followed by pyrosequencing confirmed methylation status in juvenile DM and other IIMs. Immunohistochemistry defined localization and expression levels of WT1. RESULTS: Comparison of genome-wide DNA methylation profiling between juvenile DM muscle and normal control muscle revealed 27 genes with a significant methylation difference between the groups. These genes were enriched with transcription factors and/or cell cycle regulators and were unrelated to duration of untreated disease. Six homeobox genes were among them; ALX4, HOXC11, HOXD3, and HOXD4 were hypomethylated, while EMX2 and HOXB1 were hypermethylated. WT1 was significantly hypomethylated in juvenile DM (Δß = -0.41, P < 0.001). Bisulfite pyrosequencing verification in samples from 56 patients with juvenile DM confirmed the methylation alterations of these genes. Similar methylation alterations were observed in juvenile polymyositis (n = 5) and other IIMs (n = 9). Concordant with the other findings, WT1 protein was increased in juvenile DM muscle, with average positive staining of 11.6%, but was undetectable in normal muscle (P < 0.001). CONCLUSION: These results suggest that affected muscles of children with juvenile DM and IIMs have the capacity to be repaired, and that homeobox and WT1 genes are epigenetically marked to facilitate this repair process, potentially suggesting new avenues of therapeutic intervention.

Genome-wide quantitative assessment of variation in DNA methylation patterns.

Genomic DNA methylation contributes substantively to transcriptional regulations that underlie mammalian development and cellular differentiation. Much effort has been made to decipher the molecular mechanisms governing the establishment and maintenance of DNA methylation patterns. However, little is known about genome-wide variation of DNA methylation patterns. In this study, we introduced the concept of methylation entropy, a measure of the randomness of DNA methylation patterns in a cell population, and exploited it to assess the variability in DNA methylation patterns of Alu repeats and promoters. A few interesting observations were made: (i) within a cell population, methylation entropy varies among genomic loci; (ii) among cell populations, the methylation entropies of most genomic loci remain constant; (iii) compared to normal tissue controls, some tumors exhibit greater methylation entropies; (iv) Alu elements with high methylation entropy are associated with high GC content but depletion of CpG dinucleotides and (v) Alu elements in the intronic regions or far from CpG islands are associated with low methylation entropy. We further identified 12 putative allelic-specific methylated genomic loci, including four Alu elements and eight promoters. Lastly, using subcloned normal fibroblast cells, we demonstrated the highly variable methylation patterns are resulted from low fidelity of DNA methylation inheritance.