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.

Dynamic regulation of chromatin accessibility during melanocyte stem cell activation.

Melanocyte stem cells (McSCs) of the hair follicle are necessary for hair pigmentation and can serve as melanoma cells of origin when harboring cancer-driving mutations. McSCs can be released from quiescence, activated, and undergo differentiation into pigment-producing melanocytes during the hair cycle or due to environmental stimuli, such as ultraviolet-B (UVB) exposure. However, our current understanding of the mechanisms regulating McSC stemness, activation, and differentiation remains limited. Here, to capture the differing possible states in which murine McSCs can exist, we sorted melanocyte nuclei from quiescent (telogen) skin, skin actively producing hair shafts (anagen), and skin exposed to UVB. With these sorted nuclei, we then utilized single-nucleus assay for transposase-accessible chromatin with high-throughput sequencing (snATAC-seq) and characterized three melanocyte lineages: quiescent McSCs (qMcSCs), activated McSCs (aMcSCs), and differentiated melanocytes (dMCs) that co-exist in all three skin conditions. Furthermore, we successfully identified differentially accessible genes and enriched transcription factor binding motifs for each melanocyte lineage. Our findings reveal potential gene regulators that determine these melanocyte cell states and provide new insights into how aMcSC chromatin states are regulated differently under divergent intrinsic and extrinsic cues. We also provide a publicly available online tool with a user-friendly interface to explore this comprehensive dataset, which will provide a resource for further studies on McSC regulation upon natural or UVB-mediated stem cell activation.

Remodeling of gene regulatory networks underlying thermogenic stimuli-induced adipose beiging.

Beige adipocytes are induced by cold temperatures or ß3-adrenergic receptor (Adrb3) agonists. They create heat through glucose and fatty acid (FA) oxidation, conferring metabolic benefits. The distinct and shared mechanisms by which these treatments induce beiging are unknown. Here, we perform single-nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq) on adipose tissue from mice exposed to cold or an Adrb3 agonist to identify cellular and chromatin accessibility dynamics during beiging. Both stimuli induce chromatin remodeling that influence vascularization and inflammation in adipose. Beige adipocytes from cold-exposed mice have increased accessibility at genes regulating glycolytic processes, whereas Adrb3 activation increases cAMP responses. While both thermogenic stimuli increase accessibility at genes regulating thermogenesis, lipogenesis, and beige adipocyte development, the kinetics and magnitudes of the changes are distinct for the stimuli. Accessibility changes at lipogenic genes are linked to functional changes in lipid composition of adipose. Both stimuli tend to decrease the proportion of palmitic acids, a saturated FA in adipose. However, Adrb3 activation increases the proportion of monounsaturated FAs, whereas cold increases the proportion of polyunsaturated FAs. These findings reveal common and distinct mechanisms of cold and Adrb3 induced beige adipocyte biogenesis, and identify unique functional consequences of manipulating these pathways in vivo.

Chemical, Molecular, and Single-nucleus Analysis Reveal Chondroitin Sulfate Proteoglycan Aberrancy in Fibrolamellar Carcinoma.

Fibrolamellar carcinoma (FLC) is an aggressive liver cancer with no effective therapeutic options. The extracellular environment of FLC tumors is poorly characterized and may contribute to cancer growth and/or metastasis. To bridge this knowledge gap, we assessed pathways relevant to proteoglycans, a major component of the extracellular matrix. We first analyzed gene expression data from FLC and nonmalignant liver tissue (n = 27) to identify changes in glycosaminoglycan (GAG) biosynthesis pathways and found that genes associated with production of chondroitin sulfate, but not other GAGs, are significantly increased by 8-fold. We then implemented a novel LC/MS-MS based method to quantify the abundance of different types of GAGs in patient tumors (n = 16) and found that chondroitin sulfate is significantly more abundant in FLC tumors by 6-fold. Upon further analysis of GAG-associated proteins, we found that versican (VCAN) expression is significantly upregulated at the mRNA and protein levels, the latter of which was validated by IHC. Finally, we performed single-cell assay for transposase-accessible chromatin sequencing on FLC tumors (n = 3), which revealed for the first time the different cell types in FLC tumors and also showed that VCAN is likely produced not only from FLC tumor epithelial cells but also activated stellate cells. Our results reveal a pathologic aberrancy in chondroitin (but not heparan) sulfate proteoglycans in FLC and highlight a potential role for activated stellate cells. Significance: This study leverages a multi-disciplinary approach, including state-of-the-art chemical analyses and cutting-edge single-cell genomic technologies, to identify for the first time a marked chondroitin sulfate aberrancy in FLC that could open novel therapeutic avenues in the future.

Single-cell chromatin accessibility and lipid profiling reveals SCD1-dependent metabolic shift in adipocytes induced by bariatric surgery.

Obesity promotes type 2 diabetes and cardiometabolic pathologies. Vertical sleeve gastrectomy (VSG) is used to treat obesity resulting in long-term weight loss and health improvements that precede weight loss; however, the mechanisms underlying the immediate benefits remain incompletely understood. Because adipose plays a crucial role in energy homeostasis and utilization, we hypothesized that VSG exerts its influences, in part, by modulating adipose functional states. We applied single-cell ATAC sequencing and lipid profiling to inguinal and epididymal adipose depots from mice that received sham surgery or VSG. We observed depot-specific cellular composition and chromatin accessibility patterns that were altered by VSG. Specifically, accessibility at Scd1, a fatty acid desaturase, was substantially reduced after VSG in mature adipocytes of inguinal but not epididymal depots. This was accompanied by reduced accumulation of SCD1-produced unsaturated fatty acids. Given these findings and reports that reductions in Scd1 attenuate obesity and insulin resistance our results suggest VSG exerts its beneficial effects through an inguinal depot-specific reduction of SCD1 activity.

methyl-ATAC-seq measures DNA methylation at accessible chromatin.

Chromatin features are characterized by genome-wide assays for nucleosome location, protein binding sites, three-dimensional interactions, and modifications to histones and DNA. For example, assay for transposase accessible chromatin sequencing (ATAC-seq) identifies nucleosome-depleted (open) chromatin, which harbors potentially active gene regulatory sequences; and bisulfite sequencing (BS-seq) quantifies DNA methylation. When two distinct chromatin features like these are assayed separately in populations of cells, it is impossible to determine, with certainty, where the features are coincident in the genome by simply overlaying data sets. Here, we describe methyl-ATAC-seq (mATAC-seq), which implements modifications to ATAC-seq, including subjecting the output to BS-seq. Merging these assays into a single protocol identifies the locations of open chromatin and reveals, unambiguously, the DNA methylation state of the underlying DNA. Such combinatorial methods eliminate the need to perform assays independently and infer where features are coincident.

Imprinted DNA methylation reconstituted at a non-imprinted locus.

BACKGROUND: In mammals, tight regulation of cytosine methylation is required for embryonic development and cellular differentiation. The trans-acting DNA methyltransferases that catalyze this modification have been identified and characterized; however, these proteins lack sequence specificity, leaving the mechanism of targeting unknown. A cis-acting regulator within the Rasgrf1 imprinting control region (ICR) is necessary for establishment and maintenance of local imprinted methylation. Here, we investigate whether 3-kb of sequence from the Rasgrf1 ICR is sufficient to direct appropriate imprinted methylation and target gene expression patterns when ectopically inserted at the Wnt1 locus. RESULTS: The Rasgrf1 ICR at Wnt1 lacked somatic methylation when maternally transmitted and was fully methylated upon paternal transmission, consistent with its behavior at the Rasgrf1 locus. It was unmethylated in the female germline and was enriched for methylation in the male germline, though not to the levels seen at the endogenous Rasgrf1 allele. Wnt1 expression was not imprinted by the ectopic ICR, likely due to additional sequences being required for this function. CONCLUSIONS: We have identified sequences that are sufficient for partial establishment and full maintenance of the imprinted DNA methylation patterns. Because full somatic methylation can occur without full gametic methylation, we infer that somatic methylation of the Rasgrf1 ICR is not simply a consequence of maintained gametic methylation.

Analysis of Combinatorial Epigenomic States.

Hundreds of distinct chemical modifications to DNA and histone amino acids have been described. Regulation exerted by these so-called epigenetic marks is vital to normal development, stability of cell identity through mitosis, and nongenetic transmission of traits between generations through meiosis. Loss of this regulation contributes to many diseases. Evidence indicates epigenetic marks function in combinations, whereby a given modification has distinct effects on local genome control, depending on which additional modifications are locally present. This review summarizes emerging methods for assessing combinatorial epigenomic states, as well as challenges and opportunities for their refinement.

Long non-coding RNA regulation of reproduction and development.

Noncoding RNAs (ncRNAs) have long been known to play vital roles in eukaryotic gene regulation. Studies conducted over a decade ago revealed that maturation of spliced, polyadenylated coding mRNA occurs by reactions involving small nuclear RNAs and small nucleolar RNAs; mRNA translation depends on activities mediated by transfer RNAs and ribosomal RNAs, subject to negative regulation by micro RNAs; transcriptional competence of sex chromosomes and some imprinted genes is regulated in cis by ncRNAs that vary by species; and both small-interfering RNAs and piwi-interacting RNAs bound to Argonaute-family proteins regulate post-translational modifications on chromatin and local gene expression states. More recently, gene-regulating noncoding RNAs have been identified, such as long intergenic and long noncoding RNAs (collectively referred to as lncRNAs)--a class totaling more than 100,000 transcripts in humans, which include some of the previously mentioned RNAs that regulate dosage compensation and imprinted gene expression. Here, we provide an overview of lncRNA activities, and then review the role of lncRNAs in processes vital to reproduction, such as germ cell specification, sex determination and gonadogenesis, sex hormone responses, meiosis, gametogenesis, placentation, non-genetic inheritance, and pathologies affecting reproductive tissues. Results from many species are presented to illustrate the evolutionarily conserved processes lncRNAs are involved in.

Role of PRDM16 and its PR domain in the epigenetic regulation of myogenic and adipogenic genes during transdifferentiation of C2C12 cells.

The positive regulatory domain containing 16 (PRDM16) is commonly regarded as a "switch" controlling the transdifferentiation of myoblasts to brown adipocytes. The N-positive regulatory (PR) domain, which is highly homologous to SET domain, is a characteristic structure for the PRDM family. Many SET domain containing proteins and several PRDM members have been found to possess histone methyltransferase activity, yet the role of PRDM16 and its PR domain in the epigenetic regulation of myogenic and adipogenic genes during myoblasts/adipocytes transdifferentiation remains unexplored. In this study, we transfected C2C12 myoblasts to stably express PRDM16 and observed the repression of myogenic genes and activation of adipogenic genes at both proliferation and differentiation stages. Ectopic PRDM16-induced reprogramming of myogenic and adipogenic genes was associated with the hypermethylation on some CpG sites in the enhancer or promoter of MyoD and myogenin, but the methylation status of PPARγ promoter was not affected. C2C12 cells expressing truncated PRDM16 lacking PR domain (ΔPR-PRDM16) demonstrated attenuation of both adipogenic and myogenic potentials, indicated by PPARγ inactivation and decreased triglyceride deposition, as well as a downregulation of MyoD, MyHC and MCK genes, as compared with C2C12 cells expressing intact PRDM16. Furthermore, C2C12 cells expressing ΔPR-PRDM16 exhibited significant differences in histone modifications on the promoters of MyoD and PPARγ genes. Taken together, PRDM16-induced C2C12 transdifferentiation is associated with alterations in CpG methylation of myogenic factors, and PR domain affects both myogenesis and adipogenesis with modified histone methylation marks on MyoD and PPARγ promoters.

Single molecule and single cell epigenomics.

Dynamically regulated changes in chromatin states are vital for normal development and can produce disease when they go awry. Accordingly, much effort has been devoted to characterizing these states under normal and pathological conditions. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is the most widely used method to characterize where in the genome transcription factors, modified histones, modified nucleotides and chromatin binding proteins are found; bisulfite sequencing (BS-seq) and its variants are commonly used to characterize the locations of DNA modifications. Though very powerful, these methods are not without limitations. Notably, they are best at characterizing one chromatin feature at a time, yet chromatin features arise and function in combination. Investigators commonly superimpose separate ChIP-seq or BS-seq datasets, and then infer where chromatin features are found together. While these inferences might be correct, they can be misleading when the chromatin source has distinct cell types, or when a given cell type exhibits any cell to cell variation in chromatin state. These ambiguities can be eliminated by robust methods that directly characterize the existence and genomic locations of combinations of chromatin features in very small inputs of cells or ideally, single cells. Here we review single molecule epigenomic methods under development to overcome these limitations, the technical challenges associated with single molecule methods and their potential application to single cells.

Possible chemical initiators of cognitive dysfunction in phenylketonuria, Parkinson's disease and Alzheimer's disease.

Though a great deal is known of the pathophysiology of phenylketonuria (PKU), Parkinson's disease (PD) and Alzheimer's disease (AD) very little is known regarding possible chemical species responsible for initiating the cascade of events that ultimately cause cognitive dysfunction. Can these be viewed as inborn errors in metabolism, occurring at various stages in the life cycle, analogous to adult onset diabetes? One major deficiency in understanding such conditions is the paucity of information regarding the total metabolic pathway for various amino acids that may be implicated in their causation. For example in PKU, its etiology was reported in 1934 and dietary restriction of phenylalanine proved effective for individuals with unsatisfactory metabolism of phenylalanine. Yet, current phenylalanine metabolism does not take into account fully the multiple biochemical pathways operating whose role is preventing burdensome accumulations of intermediates that can contribute to morbidity and toxicity. The same may apply for metabolism of tyrosine in PD and methionine in AD. Especially important, are the presence of labile and reactive chemical species which may be causative agents in protein alteration, misfolding and the creation of prions in neurodegenerative diseases, thereby preventing normal protein catabolism and excretion. Though genetic or epigenetic factors must be responsible, the question remains how are these translated into the chemical structures responsible for disease initiation? The purpose of this presentation is to explore potential labile metabolites in those biochemical pathways, which may be contributing factors. Finally it is worth noting, that drug development has been increasingly designed based upon targeting genetic deficiencies. The effectiveness of this approach for the treatment of these neurodegenerative illnesses will be determined in the future.

Single-molecule analysis of combinatorial epigenomic states in normal and tumor cells.

Proper placement of epigenetic marks on DNA and histones is fundamental to normal development, and perturbations contribute to a variety of disease states. Combinations of marks act together to control gene expression; therefore, detecting their colocalization is important, but because of technical challenges, such measurements are rarely reported. Instead, measurements of epigenetic marks are typically performed one at a time in a population of cells, and their colocalization is inferred by association. Here, we describe a single-molecule analytical approach that can perform direct detection of multiple epigenetic marks simultaneously and use it to identify mechanisms coordinating placement of three gene silencing marks, trimethylated histone H3 lysine 9, lysine 27 (H3K9me3, H3K27me3), and cytosine methylation (mC), in the normal and cancer genome. We show that H3K9me3 and mC are present together on individual chromatin fragments in mouse embryonic stem cells and that half of the H3K9me3 marks require mC for their placement. In contrast, mC and H3K27me3 coincidence is rare, and in fact, mC antagonizes H3K27me3 in both embryonic stem cells and primary mouse fibroblasts, indicating this antagonism is shared among primary cells. However, upon immortalization or tumorigenic transformation of mouse fibroblasts, mC is required for complete H3K27me3 placement. Importantly, in human promyelocytic cells, H3K27me3 is also dependent on mC. Because aberrant placement of gene silencing marks at tumor suppressor genes contributes to tumor progression, the improper dependency of H3K27me3 by mC in immortalized cells is likely to be fundamental to cancer. Our platform can enable other studies involving coordination of epigenetic marks and leverage efforts to discover disease biomarkers and epigenome-modifying drugs.

Coordinate regulation of DNA methylation and H3K27me3 in mouse embryonic stem cells.

Chromatin is separated into functional domains distinguished by combinatorial patterns of post-translational histone modifications and DNA methylation. Recent studies examining multiple histone modifications have found numerous chromatin states with distinct profiles of chromatin marks and functional enrichments. There are data showing coordinate regulation between DNAme and H3K27me3, which are both involved in the establishment and maintenance of epigenetic gene silencing, but the data are conflicting. Multiple studies have presented evidence to support the theory that PRC2 and DNAme cooperate to achieve silencing, or alternatively that H3K27me3 and DNAme act antagonistically. Here we examine the effect loss of either PRC2 or DNA methyltransferase activity has on the placement of the reciprocal mark in mouse ES cells. We find that DNAme is acting globally to antagonize the placement of H3K27me3, in accordance with recently published results. At least 471,011 domains in the mouse genome acquire H3K27me3 when DNAme is diminished. Of these 466,563 have been shown to be fully methylated in wildtype ES cells, indicating the effects of DNAme on H3K27me3 are direct. In a reciprocal experiment, we examine the effect loss of PRC2 has on the placement of DNAme. In contrast to the global antagonism DNAme has on the placement of H3K27me3, loss of H3K27me3 has a modest effect on DNAme, with only 4% of genes undergoing changes in DNAme, including 861 showing increases and 552 showing losses of overall DNAme. We anticipate that integrating genomic datasets where the effect of loss of a particular epigenetic mark has on the placement of other marks will help elucidate the rules governing epigenetic regulation and what role coordinate regulation of epigenetic marks plays in development and disease.

Microfluidic extraction, stretching and analysis of human chromosomal DNA from single cells.

We describe a microfluidic device for the extraction, purification and stretching of human chromosomal DNA from single cells. A two-dimensional array of micropillars in a microfluidic polydimethylsiloxane channel was designed to capture a single human cell. Megabase-long DNA strands released from the cell upon lysis are trapped in the micropillar array and stretched under optimal hydrodynamic flow conditions. Intact chromosomal DNA is entangled in the array, while other cellular components are washed from the channel. To demonstrate the entrapment principle, a single chromosome was hybridized to whole chromosome paints, and imaged by fluorescence microscopy. DNA extracted from a single cell and small cell populations (less than 100) was released from the device by restriction endonuclease digestion under continuous flow and collected for off-chip analysis. Quantification of the extracted material reveals that the microdevice efficiently extracts essentially all chromosomal DNA. The device described represents a novel platform to perform a variety of analyses on chromosomal DNA at the single cell level.

Breed-dependent transcriptional regulation of 5'-untranslated GR (NR3C1) exon 1 mRNA variants in the liver of newborn piglets.

Glucocorticoids are vital for life and regulate an array of physiological functions by binding to the ubiquitously expressed glucocorticoid receptor (GR, also known as NR3C1). Previous studies demonstrate striking breed differences in plasma cortisol levels in pigs. However, investigation into the breed-dependent GR transcriptional regulation is hampered by lacking porcine GR promoter information. In this study, we sequenced 5.3 kb upstream of the translation start codon of the porcine GR gene, and identified seven alternative 5'-untranslated exons 1-4, 1-5, 1-6, 1-7, 1-8, 1-9,10 and 1-11. Among all these mRNA variants, exons 1-4 and 1-5, as well as the total GR were expressed significantly (P<0.05) higher in the liver of newborn piglets of Large White (LW) compared with Erhualian, a Chinese indigenous breed. Overall level of CpG methylation in the region flanking exons 1-4 and 1-5 did not show breed difference. However, nuclear content of Sp1, p-CREB and GR in the liver was significantly (P<0.05) higher in LW piglets, associated with enhanced binding of p-CREB, and higher level of histone H3 acetylation in 1-4 and 1-5 promoters. In contrast, GR binding to promoters of exons 1-4 and 1-5 was significantly diminished in LW piglets, implicating the presence of negative GREs. These results indicate that the difference in the hepatic expression of GR transcript variants between two breeds of pigs is determined, at least partly, by the disparity in the binding of transcription factors and the enrichment of histone H3 acetylation to the promoters.

Real-time analysis and selection of methylated DNA by fluorescence-activated single molecule sorting in a nanofluidic channel.

Epigenetic modifications, such as DNA and histone methylation, are responsible for regulatory pathways that affect disease. Current epigenetic analyses use bisulfite conversion to identify DNA methylation and chromatin immunoprecipitation to collect molecules bearing a specific histone modification. In this work, we present a proof-of-principle demonstration for a new method using a nanofluidic device that combines real-time detection and automated sorting of individual molecules based on their epigenetic state. This device evaluates the fluorescence from labeled epigenetic modifications to actuate sorting. This technology has demonstrated up to 98% accuracy in molecule sorting and has achieved postsorting sample recovery on femtogram quantities of genetic material. We have applied it to sort methylated DNA molecules using simultaneous, multicolor fluorescence to identify methyl binding domain protein-1 (MBD1) bound to full-duplex DNA. The functionality enabled by this nanofluidic platform now provides a workflow for color-multiplexed detection, sorting, and recovery of single molecules toward subsequent DNA sequencing.

Sequences sufficient for programming imprinted germline DNA methylation defined.

Epigenetic marks are fundamental to normal development, but little is known about signals that dictate their placement. Insights have been provided by studies of imprinted loci in mammals, where monoallelic expression is epigenetically controlled. Imprinted expression is regulated by DNA methylation programmed during gametogenesis in a sex-specific manner and maintained after fertilization. At Rasgrf1 in mouse, paternal-specific DNA methylation on a differential methylation domain (DMD) requires downstream tandem repeats. The DMD and repeats constitute a binary switch regulating paternal-specific expression. Here, we define sequences sufficient for imprinted methylation using two transgenic mouse lines: One carries the entire Rasgrf1 cluster (RC); the second carries only the DMD and repeats (DR) from Rasgrf1. The RC transgene recapitulated all aspects of imprinting seen at the endogenous locus. DR underwent proper DNA methylation establishment in sperm and erasure in oocytes, indicating the DMD and repeats are sufficient to program imprinted DNA methylation in germlines. Both transgenes produce a DMD-spanning pit-RNA, previously shown to be necessary for imprinted DNA methylation at the endogenous locus. We show that when pit-RNA expression is controlled by the repeats, it regulates DNA methylation in cis only and not in trans. Interestingly, pedigree history dictated whether established DR methylation patterns were maintained after fertilization. When DR was paternally transmitted followed by maternal transmission, the unmethylated state that was properly established in the female germlines could not be maintained. This provides a model for transgenerational epigenetic inheritance in mice.