methyl-ATAC-seq measures DNA methylation at accessible chromatin

DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins in liver and skeletal muscle

Each tissue has a dominant set of functional proteins required to mediate tissue-specific functions. Epigenetic modifications, transcription, and translational efficiency control tissue-dominant protein production. However, the coordination of these regulatory mechanisms to achieve such tissue-specific protein production remains unclear. Here, we analyzed the DNA methylome, transcriptome, and proteome in mouse liver and skeletal muscle. We found that DNA hypomethylation at promoter regions is globally associated with liver-dominant or skeletal muscle-dominant functional protein production within each tissue, as well as with genes encoding proteins involved in ubiquitous functions in both tissues. Thus, genes encoding liver-dominant proteins, such as those involved in glycolysis or gluconeogenesis, the urea cycle, complement and coagulation systems, enzymes of tryptophan metabolism, and cytochrome P450-related metabolism, were hypomethylated in the liver, whereas those encoding-skeletal muscle-dominant proteins, such as those involved in sarcomere organization, were hypomethylated in the skeletal muscle. Thus, DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins.

Ado-Mediated Depletion of Taurine Impairs Mitochondrial Respiratory Capacity and Alters the Chromatin Landscape of Inguinal Adipose Tissue

Non-shivering thermogenesis (NST) has strong potential to combat obesity; however, a safe molecular approach to activate this process has not yet been identified. The sulfur amino acid taurine has the ability to safely activate NST and confer protection against obesity and metabolic disease in both mice and humans, but the mechanism of this action is unknown. In this study, we discover that a suite of taurine biosynthetic enzymes, especially that of cysteamine dioxygenase (ADO), significantly increases in response to β3 adrenergic signaling in inguinal adipose tissue (IWAT) in order to increase intracellular concentrations of taurine. We further show that ADO is critical for thermogenic mitochondrial respiratory function as its ablation in adipocytes significantly reduces taurine levels, which leads to declines in mitochondrial oxygen consumption rates. Finally, we demonstrate via assay for transposase-accessible chromatin with sequencing (ATAC-seq) that taurine supplementation in beige adipocytes has the ability to remodel the chromatin landscape to increase the chromatin accessibility and transcription of genes, such as glucose-6-phosphate isomerase 1 (Gpi1), which are critical for NST. Taken together, our studies highlight a potential mechanism for taurine in the activation of NST that can be leveraged toward the treatment of obesity and metabolic disease.

Integrative epigenome profiling of 47XXY provides insights into whole genomic DNA hypermethylation and active chromatin accessibility

Klinefelter syndrome (KS, 47XXY) is a disorder characterized by sex chromosomal aneuploidy, which may lead to changes in epigenetic regulations of gene expression. To define epigenetic architectures in 47XXY, we annotated DNA methylation in euploid males (46XY) and females (46XX), and 47XXY individuals using whole genome bisulfite sequencing (WGBS) and integrated chromatin accessbilty, and detected abnormal hypermethylation in 47XXY. Furthermore, we detected altered chromatin accessibility in 47XXY, in particular in chromosome X, using Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq) in cultured amniotic cells. Our results construct the whole genome-wide DNA methylation map in 47XXY, and provide new insights into the early epigenomic dysregulation resulting from an extra chromosome X in 47XXY.

Cross-protection interactions in insect pests: Implications for pest management in a changing climate

Agricultural insect pests display an exceptional ability to adapt quickly to natural and anthropogenic stressors. Emerging evidence suggests that frequent and varied sources of stress play an important role in driving protective physiological responses; therefore, intensively managed agroecosystems combined with climatic shifts might be an ideal crucible for stress adaptation. Cross-protection, where responses to one stressor offers protection against another type of stressor, has been well documented in many insect species, yet the molecular and epigenetic underpinnings that drive overlapping protective responses in insect pests remain unclear. In this perspective, we discuss cross-protection mechanisms and provide an argument for its potential role in increasing tolerance to a wide range of natural and anthropogenic stressors in agricultural insect pests. By drawing from existing literature on single and multiple stressor studies, we outline the processes that facilitate cross-protective interactions, including epigenetic modifications, which are understudied in insect stress responses. Finally, we discuss the implications of cross-protection for insect pest management, focusing on the consequences of cross-protection between insecticides and elevated temperatures associated with climate change. Given the multiple ways that insect pests are intensively managed in agroecosystems, we suggest that examining the role of multiple stressors can be important in understanding the wide adaptability of agricultural insect pests. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

iDNA-ABF: multi-scale deep biological language learning model for the interpretable prediction of DNA methylations

In this study, we propose iDNA-ABF, a multi-scale deep biological language learning model that enables the interpretable prediction of DNA methylations based on genomic sequences only. Benchmarking comparisons show that our iDNA-ABF outperforms state-of-the-art methods for different methylation predictions. Importantly, we show the power of deep language learning in capturing both sequential and functional semantics information from background genomes. Moreover, by integrating the interpretable analysis mechanism, we well explain what the model learns, helping us build the mapping from the discovery of important sequential determinants to the in-depth analysis of their biological functions.

Spatiotemporal specificity of correlated DNA methylation and gene expression pairs across different human tissues and stages of brain development

DNA methylation (DNAm) that occurs on promoter regions is primarily considered to repress gene expression. Previous studies indicated that DNAm could also show positive correlations with gene expression. Both DNAm and gene expression profiles are known to be tissue- and development-specific. This study aims to investigate how DNAm and gene expression are coordinated across different human tissues and developmental stages, as well as the biological significance of such correlations. By analyzing 2,239 samples with both DNAm and gene expression data in the same human subjects obtained from six published datasets, we evaluated the correlations between gene and CpG pairs (GCPs) at cis-regions and compared significantly correlated GCPs (cGCPs) across different tissues and brains at different age groups. A total of 37,363 cGCPs was identified in the six datasets; approximately 38% of the cGCPs were positively correlated. The majority (>90%) of cGCPs was tissue- or development-specific. We also observed that the correlation direction can be opposite in different tissues and ages. Further analysis highlights the importance of cGCPs for their cellular functions and potential roles in complex traits and human diseases. For instance, the early developmental brain possessed a highly unique set of cGCPs that were associated with neurogenesis and psychiatric disorders. By assessing the epigenetic factors involved in cGCPs, we discovered novel regulatory mechanisms of positive cGCPs distinct from negative cGCPs, which were related to multiple factors, such as H3K27me3, CTCF, and JARD2. The catalogue of cGCPs compiled can be used to guide functional interpretation of genetic and epigenetic studies.

Methionine- and Choline-Deficient Diet Identifies an Essential Role for DNA Methylation in Plasmacytoid Dendritic Cell Biology

Diet plays an important role in lifestyle disorders associated with the disturbed immune system. During the study of methionine- and choline-deficient diet-induced nonalcoholic fatty liver disease, we observed a specific decrease in the plasmacytoid dendritic cell (pDC) fraction from murine spleens. While delineating the role for individual components, we identified that l-methionine supplementation correlates with representation of the pDC fraction. S-adenosylmethionine (SAM) is a key methyl donor, and we demonstrate that supplementation of methionine-deficient medium with SAM but not homocysteine reverses the defect in pDC development. l-Methionine has been implicated in maintenance of methylation status in the cell. Based on our observed effect of SAM and zebularine on DC subset development, we sought to clarify the role of DNA methylation in pDC biology. Whole-genome bisulfite sequencing analysis from the splenic DC subsets identified that pDCs display differentially hypermethylated regions in comparison with classical DC (cDC) subsets, whereas cDC1 and cDC2 exhibited comparable methylated regions, serving as a control in our study. We validated differentially methylated regions in the sorted pDC, CD8α⁺ cDC1, and CD4⁺ cDC2 subsets from spleens as well as FL-BMDC cultures. Upon analysis of genes linked with differentially methylated regions, we identified that differential DNA methylation is associated with the MAPK pathway such that its inhibition guides DC development toward the pDC subtype. Overall, our study identifies an important role for methionine in pDC biology.

Comprehensive understanding of Tn5 insertion preference improves transcription regulatory element identification

Tn5 transposase, which can efficiently tagment the genome, has been widely adopted as a molecular tool in next-generation sequencing, from short-read sequencing to more complex methods such as assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we systematically map Tn5 insertion characteristics across several model organisms, finding critical parameters that affect its insertion. On naked genomic DNA, we found that Tn5 insertion is not uniformly distributed or random. To uncover drivers of these biases, we used a machine learning framework, which revealed that DNA shape cooperatively works with DNA motif to affect Tn5 insertion preference. These intrinsic insertion preferences can be modeled using nucleotide dependence information from DNA sequences, and we developed a computational pipeline to correct for these biases in ATAC-seq data. Using our pipeline, we show that bias correction improves the overall performance of ATAC-seq peak detection, recovering many potential false-negative peaks. Furthermore, we found that these peaks are bound by transcription factors, underscoring the biological relevance of capturing this additional information. These findings highlight the benefits of an improved understanding and precise correction of Tn5 insertion preference.

Acute depletion of CTCF rewires genome-wide chromatin accessibility

BACKGROUND

The transcription factor CTCF appears indispensable in defining topologically associated domain boundaries and maintaining chromatin loop structures within these domains, supported by numerous functional studies. However, acute depletion of CTCF globally reduces chromatin interactions but does not significantly alter transcription.

RESULTS

Here, we systematically integrate multi-omics data including ATAC-seq, RNA-seq, WGBS, Hi-C, Cut&Run, and CRISPR-Cas9 survival dropout screens, and time-solved deep proteomic and phosphoproteomic analyses in cells carrying auxin-induced degron at endogenous CTCF locus. Acute CTCF protein degradation markedly rewires genome-wide chromatin accessibility. Increased accessible chromatin regions are frequently located adjacent to CTCF-binding sites at promoter regions and insulator sites associated with enhanced transcription of nearby genes. In addition, we use CTCF-associated multi-omics data to establish a combinatorial data analysis pipeline to discover CTCF co-regulatory partners. We successfully identify 40 candidates, including multiple established partners. Interestingly, many CTCF co-regulators that have alterations of their respective downstream gene expression do not show changes of their own expression levels across the multi-omics measurements upon acute CTCF loss, highlighting the strength of our system to discover hidden co-regulatory partners associated with CTCF-mediated transcription.

CONCLUSIONS

This study highlights that CTCF loss rewires genome-wide chromatin accessibility, which plays a critical role in transcriptional regulation.

Chromatin insulation dynamics in glioblastoma: challenges and future perspectives of precision oncology

Glioblastoma (GBM) is the most aggressive primary brain tumor, having a poor prognosis and a median overall survival of less than two years. Over the last decade, numerous findings regarding the distinct molecular and genetic profiles of GBM have led to the emergence of several therapeutic approaches. Unfortunately, none of them has proven to be effective against GBM progression and recurrence. Epigenetic mechanisms underlying GBM tumor biology, including histone modifications, DNA methylation, and chromatin architecture, have become an attractive target for novel drug discovery strategies. Alterations on chromatin insulator elements (IEs) might lead to aberrant chromatin remodeling via DNA loop formation, causing oncogene reactivation in several types of cancer, including GBM. Importantly, it is shown that mutations affecting the isocitrate dehydrogenase (IDH) 1 and 2 genes, one of the most frequent genetic alterations in gliomas, lead to genome-wide DNA hypermethylation and the consequent IE dysfunction. The relevance of IEs has also been observed in a small population of cancer stem cells known as glioma stem cells (GSCs), which are thought to participate in GBM tumor initiation and drug resistance. Recent studies revealed that epigenomic alterations, specifically chromatin insulation and DNA loop formation, play a crucial role in establishing and maintaining the GSC transcriptional program. This review focuses on the relevance of IEs in GBM biology and their implementation as a potential theranostic target to stratify GBM patients and develop novel therapeutic approaches. We will also discuss the state-of-the-art emerging technologies using big data analysis and how they will settle the bases on future diagnosis and treatment strategies in GBM patients.

Effect of pulmonary surfactant on the prevention of neonatal respiratory distress syndrome in premature infants

OBJECTIVE

To investigate the effect of pulmonary surfactant on the prevention of neonatal respiratory distress syndrome (NRDS) in premature infants.

METHODS

A total of 216 premature infants in our hospital were selected. According to the guardian's choice, premature infants were grouped into treatment with pulmonary surfactant (PS) combined with nasal continuous positive airway pressure (NCPAP; the observation group) and treatment with NCPAP alone (the control group). The observation group was further divided into the young gestational age group (group A, 40 cases) and old gestational age group (group B, 74 cases), while the control group was divided into the young gestational age group (group C, 35 cases) and old gestational age group (group D, 67 cases). The incidence of NRDS, the first and second use of PS, the number of cases with mechanical ventilation, use of oxygen time, length of hospital stay and complications were observed.

RESULTS

The incidence of NRDS in the observation group was lower than that in the control group, but the first use rate of PS was higher than in the control group (P<0.05). The incidence of NRDS in group A was lower than that in the other three groups (P<0.001), while the mechanical ventilation rate in group C was higher than that in the other three groups (P<0.001). Use of oxygen time and length of hospital stay in groups A and C were higher than those in groups B and D (P<0.001). The mortality rate of group A was lower than that of group C (P<0.05). There was no difference in complications among patients with different gestational ages (P>0.05).

CONCLUSION

Early use of PS can prevent the incidence of NRDS in preterm infants, significantly reduce the incidence of NRDS in neonates with young gestational age. It can also improve the survival rate of patients with young gestational age without increasing the incidence of complications.

Stress-induced transcriptional memory accelerates promoter-proximal pause release and decelerates termination over mitotic divisions

Heat shock instantly reprograms transcription. Whether gene and enhancer transcription fully recover from stress and whether stress establishes a memory by provoking transcription regulation that persists through mitosis remained unknown. Here, we measured nascent transcription and chromatin accessibility in unconditioned cells and in the daughters of stress-exposed cells. Tracking transcription genome-wide at nucleotide-resolution revealed that cells precisely restored RNA polymerase II (Pol II) distribution at gene bodies and enhancers upon recovery from stress. However, a single heat exposure in embryonic fibroblasts primed a faster gene induction in their daughter cells by increasing promoter-proximal Pol II pausing and by accelerating the pause release. In K562 erythroleukemia cells, repeated stress refined basal and heat-induced transcription over mitotic division and decelerated termination-coupled pre-mRNA processing. The slower termination retained transcripts on the chromatin and reduced recycling of Pol II. These results demonstrate that heat-induced transcriptional memory acts through promoter-proximal pause release and pre-mRNA processing at transcription termination.

Characteristics of Fetal Extrachromosomal Circular DNA in Maternal Plasma: Methylation Status and Clearance

BACKGROUND

Although the characterization of cell-free extrachromosomal circular DNA (eccDNA) has gained much research interest, the methylation status of these molecules is yet to be elucidated. We set out to compare the methylation densities of plasma eccDNA of maternal and fetal origins, and between small and large molecules. The clearance of fetal eccDNA from maternal circulation was also investigated.

METHODS

We developed a sequencing protocol for eccDNA methylation analysis using tagmentation and enzymatic conversion approaches. A restriction enzyme-based approach was applied to verify the tagmentation results. The efficiency of cell-free fetal eccDNA clearance was investigated by fetal eccDNA fraction evaluations at various postpartum time points.

RESULTS

The methylation densities of fetal eccDNA (median: 56.3%; range: 40.5-67.6%) were lower than the maternal eccDNA (median: 66.7%; range: 56.5-75.7%) (P = 0.02, paired t-test). In addition, eccDNA molecules from the smaller peak cluster (180-230 bp) were of lower methylation levels than those from the larger peak cluster (300-450 bp). Both of these findings were confirmed using the restriction enzyme approach. We also observed comparable methylation densities between linear and eccDNA of both maternal and fetal origins. The average half-lives of fetal linear and eccDNA in the maternal blood were 30.2 and 29.7 min, respectively.

CONCLUSIONS

We found that fetal eccDNA in plasma was relatively hypomethylated compared to the maternal eccDNA. The methylation densities of eccDNA were positively correlated with their sizes. In addition, fetal eccDNA was found to be rapidly cleared from the maternal blood after delivery, similar to fetal linear DNA.

Direct cell reprogramming: approaches, mechanisms and progress

The reprogramming of somatic cells with defined factors, which converts cells from one lineage into cells of another, has greatly reshaped our traditional views on cell identity and cell fate determination. Direct reprogramming (also known as transdifferentiation) refers to cell fate conversion without transitioning through an intermediary pluripotent state. Given that the number of cell types that can be generated by direct reprogramming is rapidly increasing, it has become a promising strategy to produce functional cells for therapeutic purposes. This Review discusses the evolution of direct reprogramming from a transcription factor-based method to a small-molecule-driven approach, the recent progress in enhancing reprogrammed cell maturation, and the challenges associated with in vivo direct reprogramming for translational applications. It also describes our current understanding of the molecular mechanisms underlying direct reprogramming, including the role of transcription factors, epigenetic modifications, non-coding RNAs, and the function of metabolic reprogramming, and highlights novel insights gained from single-cell omics studies.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

Deciphering the multifaceted roles of TET proteins in T-cell lineage specification and malignant transformation

TET proteins are DNA demethylases that can oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC) and other oxidized mC bases (oxi-mCs). Importantly, TET proteins govern cell fate decisions during development of various cell types by activating a cell-specific gene expression program. In this review, we focus on the role of TET proteins in T-cell lineage specification. We explore the multifaceted roles of TET proteins in regulating gene expression in the contexts of T-cell development, lineage specification, function, and disease. Finally, we discuss the future directions and experimental strategies required to decipher the precise mechanisms employed by TET proteins to fine-tune gene expression and safeguard cell identity.

ATAC-seq with unique molecular identifiers improves quantification and footprinting

ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) provides an efficient way to analyze nucleosome-free regions and has been applied widely to identify transcription factor footprints. Both applications rely on the accurate quantification of insertion events of the hyperactive transposase Tn5. However, due to the presence of the PCR amplification, it is impossible to accurately distinguish independently generated identical Tn5 insertion events from PCR duplicates using the standard ATAC-seq technique. Removing PCR duplicates based on mapping coordinates introduces increasing bias towards highly accessible chromatin regions. To overcome this limitation, we establish a UMI-ATAC-seq technique by incorporating unique molecular identifiers (UMIs) into standard ATAC-seq procedures. UMI-ATAC-seq can rescue about 20% of reads that are mistaken as PCR duplicates in standard ATAC-seq in our study. We demonstrate that UMI-ATAC-seq could more accurately quantify chromatin accessibility and significantly improve the sensitivity of identifying transcription factor footprints. An analytic pipeline is developed to facilitate the application of UMI-ATAC-seq, and it is available at https://github.com/tzhu-bio/UMI-ATAC-seq .

Msuite: A High-Performance and Versatile DNA Methylation Data-Analysis Toolkit

DNA methylation is a pervasive and important epigenetic regulator in mammalian genome. For DNA methylome profiling, emerging bisulfite-free methods have demonstrated desirable superiority over the conventional bisulfite-treatment-based approaches, although current analysis software could not make full use of their advantages. In this work, we present Msuite, an easy-to-use, all-in-one data-analysis toolkit. Msuite implements a unique 4-letter analysis mode specifically optimized for emerging protocols; it also integrates quality controls, methylation call, and data visualizations. Msuite demonstrates substantial performance improvements over current state-of-the-art tools as well as fruitful functionalities, thus holding the potential to serve as an optimal toolkit to facilitate DNA methylome studies. Source codes and testing datasets for Msuite are freely available at https://github.com/hellosunking/Msuite/.

•

Msuite provides a unique 4-letter analysis mode for emerging bisulfite-free protocols

•

Msuite outperforms current tools in terms of higher accuracy and lower resource usage

•

Msuite has integrated quality control and fruitful data-visualization utilities

•

Msuite provides an all-in-one solution for DNA methylation data analysis

DNA methylation is an essential epigenetic modification responsible for many biological regulation pathways. Despite the fact that various high-throughput methods have been developed for base-resolution DNA methylome profiling, DNA methylation data analysis remains a complex and challenging task. Here, we present Msuite, which has integrated quality control, read alignment, methylation call, and fruitful data-visualization functionalities, aiming to offer an all-in-one package for most of the current DNA methylation profiling assays. Msuite also provides dedicated support for emerging bisulfite-free protocols and outperforms the current tools in terms of higher accuracy and lower computational resource requirement. Hence, Msuite could serve as the optimal toolkit for DNA methylation data analysis as well as facilitating the popularization of emerging bisulfite-free protocols.

The DNA methylation landscape of hematological malignancies: an update

The rapid advances in high-throughput sequencing technologies have made it more evident that epigenetic modifications orchestrate a plethora of complex biological processes. During the last decade, we have gained significant knowledge about a wide range of epigenetic changes that crucially contribute to some of the most aggressive forms of leukemia, lymphoma, and myelodysplastic syndromes. DNA methylation is a key epigenetic player in the abnormal initiation, development, and progression of these malignancies, often acting in synergy with other epigenetic alterations. It also contributes to the acquisition of drug resistance. In this review, we summarize the role of DNA methylation in hematological malignancies described in the current literature. We discuss in detail the dual role of DNA methylation in normal and aberrant hematopoiesis, as well as the involvement of this type of epigenetic change in other aspects of the disease. Finally, we present a comprehensive overview of the main clinical implications, including a discussion of the therapeutic strategies that regulate or reverse aberrant DNA methylation patterns in hematological malignancies, including their combination with (chemo)immunotherapy.

An enhanced chemopreventive effect of methyl donor S-adenosylmethionine in combination with 25-hydroxyvitamin D in blocking mammary tumor growth and metastasis

Therapeutic targeting of metastatic breast cancer still remains a challenge as the tumor cells are highly heterogenous and exploit multiple pathways for their growth and metastatic spread that cannot always be targeted by a single-agent monotherapy regimen. Therefore, a rational approach through simultaneous targeting of several pathways may provide a better anti-cancer therapeutic effect. We tested this hypothesis using a combination of two nutraceutical agents S-adenosylmethionine (SAM) and Vitamin D (Vit. D) prohormone [25-hydroxyvitamin D; '25(OH)D'] that are individually known to exert distinct changes in the expression of genes involved in tumor growth and metastasis. Our results show that both SAM and 25(OH)D monotherapy significantly reduced proliferation and clonogenic survival of a panel of breast cancer cell lines in vitro and inhibited tumor growth, lung metastasis, and breast tumor cell colonization to the skeleton in vivo. However, these effects were significantly more pronounced in the combination setting. RNA-Sequencing revealed that the transcriptomic footprint on key cancer-related signaling pathways is broader in the combination setting than any of the monotherapies. Furthermore, comparison of the differentially expressed genes from our transcriptome analyses with publicly available cancer-related dataset demonstrated that the combination treatment upregulates genes from immune-related pathways that are otherwise downregulated in bone metastasis in vivo. Since SAM and Vit. D are both approved nutraceuticals with known safety profiles, this combination treatment may serve as a novel strategy to reduce breast cancer-associated morbidity and mortality.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

ATACing DNA Methylation during Differentiation

Distal regulatory elements control gene expression during differentiation. In this issue of Molecular Cell, Barnett et al. (2020) develop a new technology, called ATAC-Me, and discover that removal of DNA methylation is not a pre-requisite for the creation of accessible chromatin at active gene regulatory elements during cellular differentiation.

ATAC-Me Captures Prolonged DNA Methylation of Dynamic Chromatin Accessibility Loci during Cell Fate Transitions

DNA methylation of enhancers is dynamic, cell-type specific, and vital for cell fate progression. However, current models inadequately define its role within the hierarchy of gene regulation. Analysis of independent datasets shows an unanticipated overlap between DNA methylation and chromatin accessibility at enhancers of steady-state stem cells, suggesting that these two opposing features might exist concurrently. To define their temporal relationship, we developed ATAC-Me, which probes accessibility and methylation from single DNA library preparations. We identified waves of accessibility occurring rapidly across thousands of myeloid enhancers in a monocyte-to-macrophage cell fate model. Prolonged methylation states were observed at a majority of these sites, while transcription of nearby genes tracked closely with accessibility. ATAC-Me uncovers a significant disconnect between chromatin accessibility, DNA methylation status, and gene activity. This unexpected observation highlights the value of ATAC-Me in constructing precise molecular timelines for understanding the role of DNA methylation in gene regulation.

EpiMethylTag: simultaneous detection of ATAC-seq or ChIP-seq signals with DNA methylation

Activation of regulatory elements is thought to be inversely correlated with DNA methylation levels. However, it is difficult to determine whether DNA methylation is compatible with chromatin accessibility or transcription factor (TF) binding if assays are performed separately. We developed a fast, low-input, low sequencing depth method, EpiMethylTag, that combines ATAC-seq or ChIP-seq (M-ATAC or M-ChIP) with bisulfite conversion, to simultaneously examine accessibility/TF binding and methylation on the same DNA. Here we demonstrate that EpiMethylTag can be used to study the functional interplay between chromatin accessibility and TF binding (CTCF and KLF4) at methylated sites.