Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

Nuclear TARBP2 Drives Oncogenic Dysregulation of RNA Splicing and Decay.

Post-transcriptional regulation of RNA stability is a key step in gene expression control. We describe a regulatory program, mediated by the RNA binding protein TARBP2, that controls RNA stability in the nucleus. TARBP2 binding to pre-mRNAs results in increased intron retention, subsequently leading to targeted degradation of TARBP2-bound transcripts. This is mediated by TARBP2 recruitment of the m6A RNA methylation machinery to its target transcripts, where deposition of m6A marks influences the recruitment of splicing regulators, inhibiting efficient splicing. Interactions between TARBP2 and the nucleoprotein TPR then promote degradation of these TARBP2-bound transcripts by the nuclear exosome. Additionally, analysis of clinical gene expression datasets revealed a functional role for TARBP2 in lung cancer. Using xenograft mouse models, we find that TARBP2 affects tumor growth in the lung and that this is dependent on TARBP2-mediated destabilization of ABCA3 and FOXN3. Finally, we establish ZNF143 as an upstream regulator of TARBP2 expression.

PHF6-mediated transcriptional control of NSC via Ephrin receptors is impaired in the intellectual disability syndrome BFLS.

The plant homeodomain zinc-finger protein, PHF6, is a transcriptional regulator, and PHF6 germline mutations cause the X-linked intellectual disability (XLID) Börjeson-Forssman-Lehmann syndrome (BFLS). The mechanisms by which PHF6 regulates transcription and how its mutations cause BFLS remain poorly characterized. Here, we show genome-wide binding of PHF6 in the developing cortex in the vicinity of genes involved in central nervous system development and neurogenesis. Characterization of BFLS mice harbouring PHF6 patient mutations reveals an increase in embryonic neural stem cell (eNSC) self-renewal and a reduction of neural progenitors. We identify a panel of Ephrin receptors (EphRs) as direct transcriptional targets of PHF6. Mechanistically, we show that PHF6 regulation of EphR is impaired in BFLS mice and in conditional Phf6 knock-out mice. Knockdown of EphR-A phenocopies the PHF6 loss-of-function defects in altering eNSCs, and its forced expression rescues defects of BFLS mice-derived eNSCs. Our data indicate that PHF6 directly promotes Ephrin receptor expression to control eNSC behaviour in the developing brain, and that this pathway is impaired in BFLS.

SparsePro: An efficient fine-mapping method integrating summary statistics and functional annotations.

Identifying causal variants from genome-wide association studies (GWAS) is challenging due to widespread linkage disequilibrium (LD) and the possible existence of multiple causal variants in the same genomic locus. Functional annotations of the genome may help to prioritize variants that are biologically relevant and thus improve fine-mapping of GWAS results. Classical fine-mapping methods conducting an exhaustive search of variant-level causal configurations have a high computational cost, especially when the underlying genetic architecture and LD patterns are complex. SuSiE provided an iterative Bayesian stepwise selection algorithm for efficient fine-mapping. In this work, we build connections between SuSiE and a paired mean field variational inference algorithm through the implementation of a sparse projection, and propose effective strategies for estimating hyperparameters and summarizing posterior probabilities. Moreover, we incorporate functional annotations into fine-mapping by jointly estimating enrichment weights to derive functionally-informed priors. We evaluate the performance of SparsePro through extensive simulations using resources from the UK Biobank. Compared to state-of-the-art methods, SparsePro achieved improved power for fine-mapping with reduced computation time. We demonstrate the utility of SparsePro through fine-mapping of five functional biomarkers of clinically relevant phenotypes. In summary, we have developed an efficient fine-mapping method for integrating summary statistics and functional annotations. Our method can have wide utility in understanding the genetics of complex traits and increasing the yield of functional follow-up studies of GWAS. SparsePro software is available on GitHub at https://github.com/zhwm/SparsePro.

SharePro: an accurate and efficient genetic colocalization method accounting for multiple causal signals.

MOTIVATION: Colocalization analysis is commonly used to assess whether two or more traits share the same genetic signals identified in genome-wide association studies (GWAS), and is important for prioritizing targets for functional follow-up of GWAS results. Existing colocalization methods can have suboptimal performance when there are multiple causal variants in one genomic locus. RESULTS: We propose SharePro to extend the COLOC framework for colocalization analysis. SharePro integrates linkage disequilibrium (LD) modeling and colocalization assessment by grouping correlated variants into effect groups. With an efficient variational inference algorithm, posterior colocalization probabilities can be accurately estimated. In simulation studies, SharePro demonstrated increased power with a well-controlled false positive rate at a low computational cost. Compared to existing methods, SharePro provided stronger and more consistent colocalization evidence for known lipid-lowering drug target proteins and their corresponding lipid traits. Through an additional challenging case of the colocalization analysis of the circulating abundance of R-spondin 3 GWAS and estimated bone mineral density GWAS, we demonstrated the utility of SharePro in identifying biologically plausible colocalized signals. AVAILABILITY AND IMPLEMENTATION: SharePro for colocalization analysis is written in Python and openly available at https://github.com/zhwm/SharePro_coloc.

Pyruvate Kinase M (PKM) binds ribosomes in a poly-ADP ribosylation dependent manner to induce translational stalling.

In light of the numerous studies identifying post-transcriptional regulators on the surface of the endoplasmic reticulum (ER), we asked whether there are factors that regulate compartment specific mRNA translation in human cells. Using a proteomic survey of spatially regulated polysome interacting proteins, we identified the glycolytic enzyme Pyruvate Kinase M (PKM) as a cytosolic (i.e. ER-excluded) polysome interactor and investigated how it influences mRNA translation. We discovered that the PKM-polysome interaction is directly regulated by ADP levels-providing a link between carbohydrate metabolism and mRNA translation. By performing enhanced crosslinking immunoprecipitation-sequencing (eCLIP-seq), we found that PKM crosslinks to mRNA sequences that are immediately downstream of regions that encode lysine- and glutamate-enriched tracts. Using ribosome footprint protection sequencing, we found that PKM binding to ribosomes causes translational stalling near lysine and glutamate encoding sequences. Lastly, we observed that PKM recruitment to polysomes is dependent on poly-ADP ribosylation activity (PARylation)-and may depend on co-translational PARylation of lysine and glutamate residues of nascent polypeptide chains. Overall, our study uncovers a novel role for PKM in post-transcriptional gene regulation, linking cellular metabolism and mRNA translation.

Myf6/MRF4 is a myogenic niche regulator required for the maintenance of the muscle stem cell pool.

The function and maintenance of muscle stem cells (MuSCs) are tightly regulated by signals originating from their niche environment. Skeletal myofibers are a principle component of the MuSC niche and are in direct contact with the muscle stem cells. Here, we show that Myf6 establishes a ligand/receptor interaction between muscle stem cells and their associated muscle fibers. Our data show that Myf6 transcriptionally regulates a broad spectrum of myokines and muscle-secreted proteins in skeletal myofibers, including EGF. EGFR signaling blocks p38 MAP kinase-induced differentiation of muscle stem cells. Homozygous deletion of Myf6 causes a significant reduction in the ability of muscle to produce EGF, leading to a deregulation in EGFR signaling. Consequently, although Myf6-knockout mice are born with a normal muscle stem cell compartment, they undergo a progressive reduction in their stem cell pool during postnatal life due to spontaneous exit from quiescence. Taken together, our data uncover a novel role for Myf6 in promoting the expression of key myokines, such as EGF, in the muscle fiber which prevents muscle stem cell exhaustion by blocking their premature differentiation.

Identification of sequence variation in the oocyte-derived bone morphogenetic protein 15 (BMP15) gene (BMP15) associated with litter size in New Zealand sheep (Ovis aries) breeds.

Advances in the study of reproductive traits indicate that functional variation in fertility genes may be useful for improving sheep fertility. The aim of this study was to search for variation in the bone morphogenetic protein 15 gene (BMP15) and ascertain any association with litter size in purebred Finnish Landrace sheep (n = 148), Finnish Landrace × Texel-cross sheep (n = 45), and composite sheep (of varying breed background; n = 58) from New Zealand (NZ). A 482 bp and 312 bp fragment of exon 1 and 2, respectively, of BMP15 were analysed using polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). The additive and dominance effect of BMP15 variation on litter size were estimated using animal and sire models. Two variants (A and B) were detected in exon 1; no sequence variation was detected in exon 2. Variant A had the nucleotide sequence CTT between positions c.31 and c.33, while variant B had a deletion (c.31_33del). The observed frequency for variant A in the Finnish Landrace sheep, Finnish Landrace × Texel-cross sheep and the composite sheep, was 0.77, 0.92, and 0.68, respectively while the frequency of variant B (c.31_33del) was 0.23, 0.08, and 0.32, respectively. An association between litter size and c.31_33del (P < 0.001) was observed in composite sheep. Analysis of more sheep will be required to confirm these results. Litter size did not differ significantly between sheep breeds regardless of the presence/absence of c.31_33del. Results suggested that c.31_33del might be a genetic marker for improving fecundity in some NZ sheep.

High-resolution genome-wide expression analysis of single myofibers using SMART-Seq.

Skeletal muscle is a heterogeneous tissue. Individual myofibers that make up muscle tissue exhibit variation in their metabolic and contractile properties. Although biochemical and histological assays are available to study myofiber heterogeneity, efficient methods to analyze the whole transcriptome of individual myofibers are lacking. Here, we report on a single-myofiber RNA-sequencing (smfRNA-Seq) approach to analyze the whole transcriptome of individual myofibers by combining single-fiber isolation with Switching Mechanism at 5' end of RNA Template (SMART) technology. Using smfRNA-Seq, we first determined the genes that are expressed in the whole muscle, including in nonmyogenic cells. We also analyzed the differences in the transcriptome of myofibers from young and old mice to validate the effectiveness of this new method. Our results suggest that aging leads to significant changes in the expression of metabolic genes, such as Nos1, and structural genes, such as Myl1, in myofibers. We conclude that smfRNA-Seq is a powerful tool to study developmental, disease-related, and age-related changes in the gene expression profile of skeletal muscle.

Identification of polymorphisms in the oocyte-derived growth differentiation growth factor 9 (GDF9) gene associated with litter size in New Zealand sheep (Ovis aries) breeds.

Having the ability to control litter size is important for sheep farmers and breeders worldwide. However, making genetic gain in key livestock traits like reproductive performance needs typically a lot of time, and both the fecundity and fertility traits have a great economic importance. Attention has therefore turned to better understanding the genes that control reproductive performance. Of these genes, research has focussed on the growth differentiation growth factor 9 (GDF9) gene (GDF9). In this study, a PCR-single strand conformation polymorphism (PCR-SSCP) approach was used to investigate variation in this gene in separate groups of purebred Finnish Landrace sheep, Finnish Landrace × Texel-cross sheep and composite sheep of undefined breed background, but based on New Zealand Romney-type genetics. Three GDF9 variants (named A, B and C) were found, and upon DNA sequencing, the nucleotide substitutions c.978A>G, c.994G>A and c.1111G>A were revealed. The frequency of variant A (containing nucleotides c.978A, c.994G and c.1111G) in the Finnish Landrace, Finnish Landrace × Texel-cross and composite sheep was 0.86, 0.78 and 0.76, respectively. In these three sheep groups, the frequency of B (defined by the presence of nucleotides c.978G and c.994A) was 0.01, 0.03 and 0.23 and for C (containing c.1111A) was 0.13, 0.18 and 0.01, respectively. An animal model was used to estimate the additive effect of fertility data for Finnish Landrace × Texel-cross sheep and revealed an association between litter size and the c.1111G>A variation (p = .036), but this was not observed for the Finnish Landrace sheep (p = .27) or the composite sheep (p = .17). When all the sheep were analysed together, the presence of c.1111A was associated (p < .05) with increased litter size, when compared to ewes that had c.1111G. Litter size did not differ between sheep with and without c.994A in all three groups of sheep investigated. This study suggests that c.1111A could be a useful genetic marker for improving fecundity in New Zealand sheep breeds and that it could be introgressed into other breeds, but analysis of more sheep will be required to confirm the associations that have been observed here.

Multiparameter functional diversity of human C2H2 zinc finger proteins.

C2H2 zinc finger proteins represent the largest and most enigmatic class of human transcription factors. Their C2H2-ZF arrays are highly variable, indicating that most will have unique DNA binding motifs. However, most of the binding motifs have not been directly determined. In addition, little is known about whether or how these proteins regulate transcription. Most of the ∼700 human C2H2-ZF proteins also contain at least one KRAB, SCAN, BTB, or SET domain, suggesting that they may have common interacting partners and/or effector functions. Here, we report a multifaceted functional analysis of 131 human C2H2-ZF proteins, encompassing DNA binding sites, interacting proteins, and transcriptional response to genetic perturbation. We confirm the expected diversity in DNA binding motifs and genomic binding sites, and provide motif models for 78 previously uncharacterized C2H2-ZF proteins, most of which are unique. Surprisingly, the diversity in protein-protein interactions is nearly as high as diversity in DNA binding motifs: Most C2H2-ZF proteins interact with a unique spectrum of co-activators and co-repressors. Thus, multiparameter diversification likely underlies the evolutionary success of this large class of human proteins.

A compendium of RNA-binding motifs for decoding gene regulation.

RNA-binding proteins are key regulators of gene expression, yet only a small fraction have been functionally characterized. Here we report a systematic analysis of the RNA motifs recognized by RNA-binding proteins, encompassing 205 distinct genes from 24 diverse eukaryotes. The sequence specificities of RNA-binding proteins display deep evolutionary conservation, and the recognition preferences for a large fraction of metazoan RNA-binding proteins can thus be inferred from their RNA-binding domain sequence. The motifs that we identify in vitro correlate well with in vivo RNA-binding data. Moreover, we can associate them with distinct functional roles in diverse types of post-transcriptional regulation, enabling new insights into the functions of RNA-binding proteins both in normal physiology and in human disease. These data provide an unprecedented overview of RNA-binding proteins and their targets, and constitute an invaluable resource for determining post-transcriptional regulatory mechanisms in eukaryotes.

Systematic discovery of structural elements governing stability of mammalian messenger RNAs.

Decoding post-transcriptional regulatory programs in RNA is a critical step towards the larger goal of developing predictive dynamical models of cellular behaviour. Despite recent efforts, the vast landscape of RNA regulatory elements remains largely uncharacterized. A long-standing obstacle is the contribution of local RNA secondary structure to the definition of interaction partners in a variety of regulatory contexts, including--but not limited to--transcript stability, alternative splicing and localization. There are many documented instances where the presence of a structural regulatory element dictates alternative splicing patterns (for example, human cardiac troponin T) or affects other aspects of RNA biology. Thus, a full characterization of post-transcriptional regulatory programs requires capturing information provided by both local secondary structures and the underlying sequence. Here we present a computational framework based on context-free grammars and mutual information that systematically explores the immense space of small structural elements and reveals motifs that are significantly informative of genome-wide measurements of RNA behaviour. By applying this framework to genome-wide human mRNA stability data, we reveal eight highly significant elements with substantial structural information, for the strongest of which we show a major role in global mRNA regulation. Through biochemistry, mass spectrometry and in vivo binding studies, we identified human HNRPA2B1 (heterogeneous nuclear ribonucleoprotein A2/B1, also known as HNRNPA2B1) as the key regulator that binds this element and stabilizes a large number of its target genes. We created a global post-transcriptional regulatory map based on the identity of the discovered linear and structural cis-regulatory elements, their regulatory interactions and their target pathways. This approach could also be used to reveal the structural elements that modulate other aspects of RNA behaviour.

Expression of DRD2 Is Increased in Human Pancreatic Ductal Adenocarcinoma and Inhibitors Slow Tumor Growth in Mice.

BACKGROUND & AIMS: Incidence of and mortality from pancreatic ductal adenocarcinoma (PDAC), the most common form of pancreatic cancer, are almost equivalent, so better treatments are needed. We studied gene expression profiles of PDACs and the functions of genes with altered expression to identify new therapeutic targets. METHODS: We performed microarray analysis to analyze gene expression profiles of 195 PDAC and 41 non-tumor pancreatic tissue samples. We undertook an extensive analysis of the PDAC transcriptome by superimposing interaction networks of proteins encoded by aberrantly expressed genes over signaling pathways associated with PDAC development to identify factors that might alter regulation of these pathways during tumor progression. We performed tissue microarray analysis to verify changes in expression of candidate protein using an independent set of 152 samples (40 nontumor pancreatic tissues, 63 PDAC sections, and 49 chronic pancreatitis samples). We validated the functional relevance of the candidate molecule using RNA interference or pharmacologic inhibitors in pancreatic cancer cell lines and analyses of xenograft tumors in mice. RESULTS: In an analysis of 38,276 human genes and loci, we identified 1676 genes that were significantly up-regulated and 1166 genes that were significantly down-regulated in PDAC compared with nontumor pancreatic tissues. One gene that was up-regulated and associated with multiple signaling pathways that are dysregulated in PDAC was G protein subunit αi2, which has not been previously associated with PDAC. G protein subunit αi2 mediates the effects of dopamine receptor D2 (DRD2) on cyclic adenosine monophosphate signaling; PDAC tissues had a slight but significant increase in DRD2 messenger RNA. Levels of DRD2 protein were substantially increased in PDACs, compared with non-tumor tissues, in tissue microarray analyses. RNA interference knockdown of DRD2 or inhibition with pharmacologic antagonists (pimozide and haloperidol) reduced proliferation of pancreatic cancer cells, induced endoplasmic reticulum stress and apoptosis, and reduced cell migration. RNA interference knockdown of DRD2 in pancreatic tumor cells reduced growth of xenograft tumors in mice, and administration of the DRD2 inhibitor haloperidol to mice with orthotopic xenograft tumors reduced final tumor size and metastasis. CONCLUSIONS: In gene expression profile analysis of PDAC samples, we found the DRD2 signaling pathway to be activated. Inhibition of DRD2 in pancreatic cancer cells reduced proliferation and migration, and slowed growth of xenograft tumors in mice. DRD2 antagonists routinely used for management of schizophrenia might be tested in patients with pancreatic cancer.

Motif comparison based on similarity of binding affinity profiles.

Measuring motif similarity is essential for identifying functionally related transcription factors (TFs) and RNA-binding proteins, and for annotating de novo motifs. Here, we describe Motif Similarity Based on Affinity of Targets (MoSBAT), an approach for measuring the similarity of motifs by computing their affinity profiles across a large number of random sequences. We show that MoSBAT successfully associates de novo ChIP-seq motifs with their respective TFs, accurately identifies motifs that are obtained from the same TF in different in vitro assays, and quantitatively reflects the similarity of in vitro binding preferences for pairs of TFs. AVAILABILITY AND IMPLEMENTATION: MoSBAT is available as a webserver at mosbat.ccbr.utoronto.ca, and for download at github.com/csglab/MoSBAT. CONTACT: t.hughes@utoronto.ca or hamed.najafabadi@mcgill.caSupplementary information: Supplementary data are available at Bioinformatics online.

A structural approach reveals how neighbouring C2H2 zinc fingers influence DNA binding specificity.

Development of an accurate protein-DNA recognition code that can predict DNA specificity from protein sequence is a central problem in biology. C2H2 zinc fingers constitute by far the largest family of DNA binding domains and their binding specificity has been studied intensively. However, despite decades of research, accurate prediction of DNA specificity remains elusive. A major obstacle is thought to be the inability of current methods to account for the influence of neighbouring domains. Here we show that this problem can be addressed using a structural approach: we build structural models for all C2H2-ZF-DNA complexes with known binding motifs and find six distinct binding modes. Each mode changes the orientation of specificity residues with respect to the DNA, thereby modulating base preference. Most importantly, the structural analysis shows that residues at the domain interface strongly and predictably influence the binding mode, and hence specificity. Accounting for predicted binding mode significantly improves prediction accuracy of predicted motifs. This new insight into the fundamental behaviour of C2H2-ZFs has implications for both improving the prediction of natural zinc finger-binding sites, and for prioritizing further experiments to complete the code. It also provides a new design feature for zinc finger engineering.

Identification of C2H2-ZF binding preferences from ChIP-seq data using RCADE.

UNLABELLED: Current methods for motif discovery from chromatin immunoprecipitation followed by sequencing (ChIP-seq) data often identify non-targeted transcription factor (TF) motifs, and are even further limited when peak sequences are similar due to common ancestry rather than common binding factors. The latter aspect particularly affects a large number of proteins from the Cys2His2 zinc finger (C2H2-ZF) class of TFs, as their binding sites are often dominated by endogenous retroelements that have highly similar sequences. Here, we present recognition code-assisted discovery of regulatory elements (RCADE) for motif discovery from C2H2-ZF ChIP-seq data. RCADE combines predictions from a DNA recognition code of C2H2-ZFs with ChIP-seq data to identify models that represent the genuine DNA binding preferences of C2H2-ZF proteins. We show that RCADE is able to identify generalizable binding models even from peaks that are exclusively located within the repeat regions of the genome, where state-of-the-art motif finding approaches largely fail. AVAILABILITY AND IMPLEMENTATION: RCADE is available as a webserver and also for download at http://rcade.ccbr.utoronto.ca/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. CONTACT: t.hughes@utoronto.ca.

Global identification of conserved post-transcriptional regulatory programs in trypanosomatids.

While regulatory programs are extensively studied at the level of transcription, elements that are involved in regulation of post-transcriptional processes are largely unknown, and methods for systematic identification of these elements are in early stages. Here, using a novel computational framework, we have integrated sequence information with several functional genomics data sets to characterize conserved regulatory programs of trypanosomatids, a group of eukaryotes that almost entirely rely on post-transcriptional processes for regulation of mRNA abundance. This analysis revealed a complex network of linear and structural RNA elements that potentially govern mRNA abundance across different life stages and environmental conditions. Furthermore, we show that the conserved regulatory network that we have identified is responsive to chemical perturbation of several biological functions in trypanosomatids. We have further characterized one of the most abundant regulatory RNA elements that we discovered, an AU-rich element (ARE) that can be found in 3' untranslated region of many trypanosomatid genes. Using bioinformatics approaches as well as in vitro and in vivo experiments, we have identified three ELAV-like homologs, including the developmentally critical protein TbRBP6, which regulate abundance of a large number of trypanosomatid ARE-containing transcripts. Together, these studies lay out a roadmap for characterization of mechanisms that modulate development and metabolic pathways in trypanosomatids.

Photocatalytic removal of 2-nitrophenol using silver and sulfur co-doped TiO2under natural solar light.

To overcome the drawback of poor solar light utilization brought about by the narrow photoresponse range of TiO2, a silver and sulfur co-doped TiO2was synthesized. Using the prepared catalyst, solar photocatalytic degradation of 2-nitrophenol (2-NP) by a TiO2-based catalyst was studied for the first time. Effects of the co-doping on the structural, optical and morphological properties of the synthesized nanoparticles were investigated by different characterization methods: X-ray diffraction, N2 adsorption-desorption measurements, scanning electron microscopy, transmission electron microscopy, energy-dispersive X-ray spectroscopy, UV-visible diffuse reflectance spectroscopy and Fourier transform infrared spectroscopy. Solar experiments showed that the co-doping with silver and sulfur significantly increased the photocatalytic activity. In various initial concentrations of 2-NP more than 99% of the contaminant was decomposed by Ag-S/TiO2in less than 150 minutes, while the degradation efficiency was much less in the presence of bare TiO2. Kinetic studies suggested that solar photocatalytic degradation of 2-NP is consistent with the Langmuir-Hinshelwood model. The rate constant of the reaction and adsorption constant of the modified photocatalyst were found to be 2.4 and 4.1 times larger than that of bare TiO2, respectively.

Transcriptional reprogramming of skeletal muscle stem cells by the niche environment.

Adult stem cells are indispensable for tissue regeneration, but their function declines with age. The niche environment in which the stem cells reside plays a critical role in their function. However, quantification of the niche effect on stem cell function is lacking. Using muscle stem cells (MuSC) as a model, we show that aging leads to a significant transcriptomic shift in their subpopulations accompanied by locus-specific gain and loss of chromatin accessibility and DNA methylation. By combining in vivo MuSC transplantation and computational methods, we show that the expression of approximately half of all age-altered genes in MuSCs from aged male mice can be restored by exposure to a young niche environment. While there is a correlation between gene reversibility and epigenetic alterations, restoration of gene expression occurs primarily at the level of transcription. The stem cell niche environment therefore represents an important therapeutic target to enhance tissue regeneration in aging.