The gene family with sequence similarity 13 member A (FAM13A) has recently been identified as a marker gene in insulin sensitivity and lipolysis. In this study, we first analyzed the expression patterns of this gene in different tissues of adult cattle and then constructed a phylogenetic tree based on the FAM13A amino acid sequence. This showed that subcutaneous adipose tissue had the highest expression in all tissues except lung tissue. Then we summarized the gene structure. The promoter region sequence of the gene was successfully amplified, and the -241/+54 region has been identified as the core promoter region. The core promoter region was determined by the unidirectional deletion of the 5' flanking promoter region of the FAM13A gene. Based on the bioinformatics analysis, we examined the dual luciferase activity of the vector constructed by the mutation site, and the transcription factors ACSL1 and ASCL2 were found as transcriptional regulators of FAM13A. Moreover, electrophoretic mobility shift assay (EMSA) further validated the regulatory role of ACSL1 and ASCL2 in the regulation of FAM13A. ACSL1 and ASCL2 were finally identified as activating transcription factors. Our results provide a basis for the function of the FAM13A gene in bovine adipocytes in order to improve the deposition of fat deposition in beef cattle muscle.

Access to genome-wide data provides the opportunity to address questions concerning the ability of transcription factors (TFs) to assemble in distinct macromolecular complexes. Here, we introduce the PAnDA (Protein And DNA Associations) approach to characterize DNA associations with human TFs using expression profiles, protein-protein interactions and recognition motifs. Our method predicts TF binding events with >0.80 accuracy revealing cell-specific regulatory patterns that can be exploited for future investigations. Even when the precise DNA-binding motifs of a specific TF are not available, the information derived from protein-protein networks is sufficient to perform high-confidence predictions (area under the ROC curve of 0.89). PAnDA is freely available at http://service.tartaglialab.com/new_submission/panda.

Paralogous transcription factors (TFs) are oftentimes reported to have identical DNA-binding motifs, despite the fact that they perform distinct regulatory functions. Differential genomic targeting by paralogous TFs is generally assumed to be due to interactions with protein co-factors or the chromatin environment. Using a computational-experimental framework called iMADS (integrative modeling and analysis of differential specificity), we show that, contrary to previous assumptions, paralogous TFs bind differently to genomic target sites even in vitro. We used iMADS to quantify, model, and analyze specificity differences between 11 TFs from 4 protein families. We found that paralogous TFs have diverged mainly at medium- and low-affinity sites, which are poorly captured by current motif models. We identify sequence and shape features differentially preferred by paralogous TFs, and we show that the intrinsic differences in specificity among paralogous TFs contribute to their differential in vivo binding. Thus, our study represents a step forward in deciphering the molecular mechanisms of differential specificity in TF families.

Transcriptional signatures are an indispensible source of correlative information on disease-related molecular alterations on a genome-wide level. Numerous candidate genes involved in disease and in factors of predictive, as well as of prognostic, value have been deduced from such molecular portraits, e.g. in cancer. However, mechanistic insights into the regulatory principles governing global transcriptional changes are lagging behind extensive compilations of deregulated genes. To identify regulators of transcriptome alterations, we used an integrated approach combining transcriptional profiling of colorectal cancer cell lines treated with inhibitors targeting the receptor tyrosine kinase (RTK)/RAS/mitogen-activated protein kinase pathway, computational prediction of regulatory elements in promoters of co-regulated genes, chromatin-based and functional cellular assays. We identified commonly co-regulated, proliferation-associated target genes that respond to the MAPK pathway. We recognized E2F and NFY transcription factor binding sites as prevalent motifs in those pathway-responsive genes and confirmed the predicted regulatory role of Y-box binding protein 1 (YBX1) by reporter gene, gel shift, and chromatin immunoprecipitation assays. We also validated the MAPK-dependent gene signature in colorectal cancers and provided evidence for the association of YBX1 with poor prognosis in colorectal cancer patients. This suggests that MEK/ERK-dependent, YBX1-regulated target genes are involved in executing malignant properties.

The exosome is a conserved multi-subunit ribonuclease complex that functions in 3' end processing, turnover and surveillance of nuclear and cytoplasmic RNAs. In the yeast nucleus, the 10-subunit core complex of the exosome (Exo-10) physically and functionally interacts with the Rrp6 exoribonuclease and its associated cofactor Rrp47, the helicase Mtr4 and Mpp6. Here, we show that binding of Mtr4 to Exo-10 in vitro is dependent upon both Rrp6 and Rrp47, whereas Mpp6 binds directly and independently of other cofactors. Crystallographic analyses reveal that the N-terminal domains of Rrp6 and Rrp47 form a highly intertwined structural unit. Rrp6 and Rrp47 synergize to create a composite and conserved surface groove that binds the N-terminus of Mtr4. Mutation of conserved residues within Rrp6 and Mtr4 at the structural interface disrupts their interaction and inhibits growth of strains expressing a C-terminal GFP fusion of Mtr4. These studies provide detailed structural insight into the interaction between the Rrp6-Rrp47 complex and Mtr4, revealing an important link between Mtr4 and the core exosome.

Regulatory interactions mediated by transcription factors (TFs) make up complex networks that control cellular behavior. Fully understanding these gene regulatory networks (GRNs) offers greater insight into the consequences of disease-causing perturbations than can be achieved by studying single TF binding events in isolation. Chromosomal translocations of the lysine methyltransferase 2A (KMT2A) gene produce KMT2A fusion proteins such as KMT2A-AFF1 (previously MLL-AF4), causing poor prognosis acute lymphoblastic leukemias (ALLs) that sometimes relapse as acute myeloid leukemias (AMLs). KMT2A-AFF1 drives leukemogenesis through direct binding and inducing the aberrant overexpression of key genes, such as the anti-apoptotic factor BCL2 and the proto-oncogene MYC However, studying direct binding alone does not incorporate possible network-generated regulatory outputs, including the indirect induction of gene repression. To better understand the KMT2A-AFF1-driven regulatory landscape, we integrated ChIP-seq, patient RNA-seq, and CRISPR essentiality screens to generate a model GRN. This GRN identified several key transcription factors such as RUNX1 that regulate target genes downstream of KMT2A-AFF1 using feed-forward loop (FFL) and cascade motifs. A core set of nodes are present in both ALL and AML, and CRISPR screening revealed several factors that help mediate response to the drug venetoclax. Using our GRN, we then identified a KMT2A-AFF1:RUNX1 cascade that represses CASP9, as well as KMT2A-AFF1-driven FFLs that regulate BCL2 and MYC through combinatorial TF activity. This illustrates how our GRN can be used to better connect KMT2A-AFF1 behavior to downstream pathways that contribute to leukemogenesis, and potentially predict shifts in gene expression that mediate drug response.

Core-binding factor β (Cbfβ) is a subunit of the Cbf family of heterodimeric transcription factors, which plays a critical role in skeletal development through its interaction with the Cbfα subunits, also known as Runt-related transcription factors (Runxs). However, the mechanism by which Cbfβ regulates cartilage and bone development remains unclear. Existing Cbfβ-deficient mouse models cannot specify the role of Cbfβ in skeletal cell lineage. Herein, we sought to specifically address the role of Cbfβ in cartilage and bone development by using a conditional knockout (CKO) approach. A mesenchymal-specific Cbfβ CKO mouse model was generated by using the Dermo1-Cre mouse line to specifically delete Cbfβ in mesenchymal stem cells, which give rise to osteoblasts and chondrocytes. Surprisingly, the mutant mice had under-developed larynx and tracheal cartilage, causing alveolus defects that led to death shortly after birth from suffocation. Also, the mutant mice exhibited severe skeletal deformities from defective intramembranous and endochondral ossification, owing to delayed chondrocyte maturation and impaired osteoblast differentiation. Almost all bones of the mutant mice, including the calvariae, vertebrae, tibiae, femurs, ribs, limbs and sternums were defective. Importantly, we showed that Cbfβ was expressed throughout the skeleton during both embryonic and postnatal development, which explains the multiple-skeletal defects observed in the mutant mice. Consistently, Cbfβ deficiency impaired both chondrocyte proliferation and hypertrophy zone hypertrophy during growth-plate development in the long bones of mutant mice. Notably, Cbfβ, Runx1 and Runx2 displayed different expression patterns in the growth plates of the wild-type mice, indicating that Cbfβ/Runx1 complex and Cbfβ/Runx2 complex may regulate chondrocyte proliferation and hypertrophy, respectively, in a spatial and temporal manner. Cbfβ deletion in the mesenchymal progenitors affected bone development by dramatically down-regulating Collagen X (Col X) and Osterix (Osx) but had a dispensable effect on osteoclast development. Collectively, the results demonstrate that Cbfβ mediates cartilage and bone development by interacting with Runx1 and Runx2 to regulate the expressions of Col X and Osx for chondrocyte and osteoblast development. These findings not only reveal a critical role for Cbfβ in cartilage and bone development but also facilitate the design of novel therapeutic approaches for skeletal diseases.

WRKY transcription factors have been shown to play a major role in regulating, both positively and negatively, the plant defense transcriptome. Nearly all studied WRKY factors appear to have a stereotypic binding preference to one DNA element termed the W-box. How specificity for certain promoters is accomplished therefore remains completely unknown. In this study, we tested five distinct Arabidopsis WRKY transcription factor subfamily members for their DNA binding selectivity towards variants of the W-box embedded in neighboring DNA sequences. These studies revealed for the first time differences in their binding site preferences, which are partly dependent on additional adjacent DNA sequences outside of the TTGACY-core motif. A consensus WRKY binding site derived from these studies was used for in silico analysis to identify potential target genes within the Arabidopsis genome. Furthermore, we show that even subtle amino acid substitutions within the DNA binding region of AtWRKY11 strongly impinge on its binding activity. Additionally, all five factors were found localized exclusively to the plant cell nucleus and to be capable of trans-activating expression of a reporter gene construct in vivo.

Ets transcription factors, which share the conserved Ets DNA-binding domain, number nearly 30 members in humans and are particularly involved in developmental processes. Their deregulation following changes in expression, transcriptional activity, or by chromosomal translocation plays a critical role in carcinogenesis. Ets DNA binding, selectivity, and regulation have been extensively studied; however, questions still arise regarding binding specificity outside the core GGA recognition sequence and the mode of action of Ets post-translational modifications. Here, we report the crystal structures of Etv1, Etv4, Etv5, and Fev, alone and in complex with DNA. We identify previously unrecognized features of the protein-DNA interface. Interactions with the DNA backbone account for most of the binding affinity. We describe a highly coordinated network of water molecules acting in base selection upstream of the GGAA core and the structural features that may account for discrimination against methylated cytidine residues. Unexpectedly, all proteins crystallized as disulfide-linked dimers, exhibiting a novel interface (distant to the DNA recognition helix). Homodimers of Etv1, Etv4, and Etv5 could be reduced to monomers, leading to a 40-200-fold increase in DNA binding affinity. Hence, we present the first indication of a redox-dependent regulatory mechanism that may control the activity of this subset of oncogenic Ets transcription factors.

The development of high-throughput genomic sequencing coupled with chromatin immunoprecipitation technologies allows studying the binding sites of the protein transcription factors (TF) in the genome scale. The growth of data volume on the experimentally determined binding sites raises qualitatively new problems for the analysis of gene expression regulation, prediction of transcription factors target genes, and regulatory gene networks reconstruction. Genome regulation remains an insufficiently studied though plants have complex molecular regulatory mechanisms of gene expression and response to environmental stresses. It is important to develop new software tools for the analysis of the TF binding sites location and their clustering in the plant genomes, visualization, and the following statistical estimates. This study presents application of the analysis of multiple TF binding profiles in three evolutionarily distant model plant organisms. The construction and analysis of non-random ChIP-seq binding clusters of the different TFs in mammalian embryonic stem cells were discussed earlier using similar bioinformatics approaches. Such clusters of TF binding sites may indicate the gene regulatory regions, enhancers and gene transcription regulatory hubs. It can be used for analysis of the gene promoters as well as a background for transcription networks reconstruction. We discuss the statistical estimates of the TF binding sites clusters in the model plant genomes. The distributions of the number of different TFs per binding cluster follow same power law distribution for all the genomes studied. The binding clusters in Arabidopsis thaliana genome were discussed here in detail.

Myeloid progenitor cells give rise to a variety of progenies including dendritic cells. However, the mechanism controlling the diversification of myeloid progenitors into each progeny is largely unknown. PU.1 and CCAAT/enhancing binding protein (C/EBP) family transcription factors have been characterized as key regulators for the development and function of the myeloid system. However, the roles of C/EBP transcription factors have not been fully identified because of functional redundancy among family members. Using high titer--retroviral infection, we demonstrate that a dominant-negative C/EBP completely blocked the granulocyte--macrophage commitment of human myeloid progenitors. Alternatively, Langerhans cell (LC) commitment was markedly facilitated in the absence of tumor necrosis factor (TNF)alpha, a strong inducer of LC development, whereas expression of wild-type C/EBP in myeloid progenitors promoted granulocytic differentiation, and completely inhibited TNFalpha-dependent LC development. On the other hand, expression of wild-type PU.1 in myeloid progenitors triggered LC development in the absence of TNFalpha, and its instructive effect was canceled by coexpressed C/EBP. Our findings establish reciprocal roles for C/EBP and PU.1 in LC development, and provide new insight into the molecular mechanism of LC development, which has not yet been well characterized.

We present a major update of the HOCOMOCO collection that provides DNA binding specificity patterns of 949 human transcription factors and 720 mouse orthologs. To make this release, we performed motif discovery in peak sets that originated from 14 183 ChIP-Seq experiments and reads from 2554 HT-SELEX experiments yielding more than 400 thousand candidate motifs. The candidate motifs were annotated according to their similarity to known motifs and the hierarchy of DNA-binding domains of the respective transcription factors. Next, the motifs underwent human expert curation to stratify distinct motif subtypes and remove non-informative patterns and common artifacts. Finally, the curated subset of 100 thousand motifs was supplied to the automated benchmarking to select the best-performing motifs for each transcription factor. The resulting HOCOMOCO v12 core collection contains 1443 verified position weight matrices, including distinct subtypes of DNA binding motifs for particular transcription factors. In addition to the core collection, HOCOMOCO v12 provides motif sets optimized for the recognition of binding sites in vivo and in vitro, and for annotation of regulatory sequence variants. HOCOMOCO is available at https://hocomoco12.autosome.org and https://hocomoco.autosome.org.

Background: Alveolar bone defect regeneration has long been problematic in the field of dentistry. Gingival stromal progenitor cells (GSPCs) offer a promising solution for alveolar bone regeneration. In order to optimally differentiate and proliferate progenitor cells, growth factors (GFs) are required. Platelet rich fibrin (PRF) has many GFs and can be easily manufactured. Core-binding factor subunit-α1 (CBF-α1) constitutes a well-known osteogenic differentiation transcription factor in SPCs. Sox9, as a chondrogenic transcription factor, interacts and inhibits CBF-α1, but its precise role in direct in vitro osteogenesis remains unknown. GSPCs cultured in vitro in PRF to optimally stimulate osteogenic differentiation has been largely overlooked. The aim of this study was to analyze GSPCs cultured in PRF osteogenic differentiation predicted by CBF-α1/Sox9. Methods: This study used a true experimental with post-test only control group design and random sampling. GPSCs isolated from the lower gingiva of four healthy, 250-gram, 1-month old, male Wistar rats ( Rattus Novergicus) were cultured for two weeks, passaged every 4-5 days. GSPCs in passage 3-5 were cultured in five M24 plates (N=108; n=6/group) for Day 7, Day 14, and Day 21 in three different mediums (control negative group: αModified Eagle Medium; control positive group: High Glucose-Dulbecco's Modified Eagle Medium (DMEM-HG) + osteogenic medium; Treatment group: DMEM-HG + osteogenic medium + PRF). CBF-α1 and Sox9 were examined with ICC monoclonal antibody. A one-way ANOVA continued with Tukey HSD test (p<0.05) based on Kolmogorov-Smirnov and Levene's tests (p>0.05) was performed. Results: The treatment group showed the highest CBF-α1/Sox9 ratio (16.00±3.000/14.33±2.517) on Day 7, while the lowest CBF-α1/Sox9 ratio (3.33±1.528/3.67±1.155) occurred in the control negative group on Day 21, with significant difference between the groups (p<0.05). Conclusion: GSPCs cultured in PRF had potential osteogenic differentiation ability predicted by the CBF-α1/sox9 ratio.

In pluripotent stem cells, there is increasing evidence for crosstalk between post-transcriptional and transcriptional networks, offering multifold steps at which pluripotency can be controlled. In addition to well-studied transcription factors, chromatin modifiers and miRNAs, RNA-binding proteins are emerging as fundamental players in pluripotency regulation. Here, we report a new role for the RNA-binding protein ESRP1 in the control of pluripotency. Knockdown of Esrp1 in mouse embryonic stem cells induces, other than the well-documented epithelial to mesenchymal-like state, also an increase in expression of the core transcription factors Oct4, Nanog and Sox2, thereby enhancing self-renewal of these cells. Esrp1-depleted embryonic stem cells displayed impaired early differentiation in vitro and formed larger teratomas in vivo when compared to control embryonic stem cells. We also show that ESRP1 binds to Oct4 and Sox2 mRNAs and decreases their polysomal loading. ESRP1 thus acts as a physiological regulator of the finely-tuned balance between self-renewal and commitment to a restricted developmental fate. Importantly, both mouse and human epithelial stem cells highly express ESRP1, pinpointing the importance of this RNA-binding protein in stem cell biology.

Mechanical stimulation and biological factors coordinately regulate bone development and regeneration; however, the underlying mechanisms are poorly understood. Microgravity induces bone loss, which may be partly related to the development of resistance to local cytokines, including insulin-like growth factor 1 (IGF-1). Here, we report the involvement of integrin αvβ3 in microgravity-associated bone loss. An established OSE-3T3 cell model was stably transfected with a 6OSE2 (Osteoblast-Specific Element 2)-luciferase reporter and cultured under simulated microgravity (SMG) and hypergravity (HG) conditions in the presence or absence of IGF-1, the disintegrin echistatin, the phosphoinositide 3-kinase (PI3K) inhibitor LY294002, or combinations of these agents. Activity of core-binding factor α1 (Cbfa1), an essential transcription factor for osteoblastic differentiation and osteogenesis, was reflected by luciferase activity. Different gravity conditions affected the induction of IGF-1 and subsequent effects on Cbfa1 transcription activity. SMG and HG influenced the expression and activity of integrin αvβ3 and phosphorylation level of p85. LY294002 inhibited the effects of HG or IGF-1 on Cbfa1 activity, indicating that HG and IGF-1 could increase Cbfa1 activity via PI3K signaling. Inhibition of integrin αvβ3 by echistatin attenuated the induction of IGF-1 and thus its effect on Cbfa1 activity under normal and HG conditions. Co-immunoprecipitation demonstrated that integrin β3 interacted with insulin receptor substrate 1, and that this interaction was decreased under SMG and increased under HG conditions. These results suggest that integrin αvβ3 mediates the synergetic regulation of Cbfa1 transcription activity by gravity and IGF-1 via PI3K signaling.

Circadian rhythms in transcription are generated by rhythmic abundances and DNA binding activities of transcription factors. Propagation of rhythms to transcriptional initiation involves the core promoter, its chromatin state, and the basal transcription machinery. Here, I characterize core promoters and chromatin states of genes transcribed in a circadian manner in mouse liver and in Drosophila. It is shown that the core promoter is a critical determinant of circadian mRNA expression in both species. A distinct core promoter class, strong circadian promoters (SCPs), is identified in mouse liver but not Drosophila. SCPs are defined by specific core promoter features, and are shown to drive circadian transcriptional activities with both high averages and high amplitudes. Data analysis and mathematical modeling further provided evidence for rhythmic regulation of both polymerase II recruitment and pause release at SCPs. The analysis provides a comprehensive and systematic view of core promoters and their link to circadian mRNA expression in mouse and Drosophila, and thus reveals a crucial role for the core promoter in regulated, dynamic transcription.

Position weight matrices (PWMs) have become a tool of choice for the identification of transcription factor binding sites in DNA sequences. DNA-binding proteins often show degeneracy in their binding requirement and thus the overall binding specificity of many proteins is unknown and remains an active area of research. Although existing PWMs are more reliable predictors than consensus string matching, they generally result in a high number of false positive hits. Our previous study introduced a promising approach to PWM refinement in which known motifs are used to computationally mine putative binding sites directly from aligned promoter regions using composition of similar sites. In the present study, we extended this technique originally tested on single examples of transcription factors (TFs) and showed its capability to optimize PWM performance to predict new binding sites in the fruit fly genome. We propose refined PWMs in mono- and dinucleotide versions similarly computed for a large variety of transcription factors of Drosophila melanogaster. Along with the addition of many auxiliary sites the optimization includes variation of the PWM motif length, the binding sites location on the promoters and the PWM score threshold. To assess the predictive performance of the refined PWMs we compared them to conventional TRANSFAC and JASPAR sources. The results have been verified using performed tests and literature review. Overall, the refined PWMs containing putative sites derived from real promoter content processed using optimized parameters had better general accuracy than conventional PWMs.

Metabolic syndrome is a newly identified risk factor for hepatocellular carcinoma (HCC); however, tumor-specific biomarkers still remain unclear. We performed cross-species analysis to compare gene signatures of HCC from human patients and melanocortin 4 receptor-knockout mice, which develop HCC with obesity, insulin resistance, and dyslipidemia. Unsupervised hierarchical clustering and principle component analysis of 746 differentially expressed orthologous genes classified HCC of 152 human patients and melanocortin 4 receptor-knockout mice into two distinct subgroups, one of which included mouse HCC and was causatively associated with metabolic risk factors. Nine genes commonly overexpressed in human and mouse metabolic disease-associated HCC were identified; fatty acid binding protein 4 (FABP4) was remarkably enriched in intratumoral activated hepatic stellate cells (HSCs). Subclones constitutively expressing FABP4 were established from a human HSC cell line in which expression levels of inflammatory chemokines, including IL-1A and IL-6, were up-regulated through NF-κB nuclear translocation, resulting in recruitment of macrophages. An immunohistochemical validation study of 106 additional human HCC samples indicated that FABP4-positive HSCs were distributed in tumors of 38 cases, and the FABP4-high group consisted of patients with nonviral and nonalcoholic HCC (P = 0.027) and with multiple metabolic risk factors (P < 0.001) compared with the FABP4-low group. Thus, FABP4 overexpression in HSCs may contribute to hepatocarcinogenesis in patients with metabolic risk factors by modulation of inflammatory pathways.

Despite the importance of proper cell death regulation across broad evolutionary distances, an understanding of the molecular machinery underpinning this fundamental process in plants remains largely elusive. This is despite its critical importance to development, homeostasis, and proper responses to stress. The identification of endogenous plant regulators of cell death has been hindered by the fact that many core regulators of cell death in animals are absent in plant genomes. Remarkably, numerous studies have shown that the ectopic expression of animal prosurvival genes in plants can suppress cell death imposed by many stresses. In this study, we capitalize on the ectopic expression of one of these animal prosurvival genes, an inhibitor of apoptosis from Spodoptera frugiperda (SfIAP), to identify novel cell death regulators in plants. A yeast two-hybrid assay was conducted using SfIAP as bait to screen a tomato cDNA library. This screen identified several transcription factors of the SQUAMOSA promoter-binding protein (SBP) family as potential SfIAP binding partners. We confirmed this interaction in vivo for our top two interactors, SlySBP8b and SlySBP12a, using coimmunoprecipitation. Interestingly, overexpression of SlySBP8b and SlySBP12a induced cell death in Nicotiana benthamiana leaves. Overexpression of these two transcription factors also induced the accumulation of reactive oxygen species and enhanced the growth of the necrotrophic pathogen Alternaria alternata. Fluorescence microscopy confirmed the nuclear localization of both SlySBP8b and SlySBP12a, while SlySBP12a was also localized to the ER membrane. These results suggest a prodeath role for SlySBP8b and SlySBP12a and implicate ER membrane tethering as a means of regulating SlySBP12a activity.

Emerging evidence points to an unexpected diversification of core promoter recognition complexes that serve as important regulators of cell-type specific gene transcription. Here, we report that the orphan TBP-associated factor TAF9B is selectively up-regulated upon in vitro motor neuron differentiation, and is required for the transcriptional induction of specific neuronal genes, while dispensable for global gene expression in murine ES cells. TAF9B binds to both promoters and distal enhancers of neuronal genes, partially co-localizing at binding sites of OLIG2, a key activator of motor neuron differentiation. Surprisingly, in this neuronal context TAF9B becomes preferentially associated with PCAF rather than the canonical TFIID complex. Analysis of dissected spinal column from Taf9b KO mice confirmed that TAF9B also regulates neuronal gene transcription in vivo. Our findings suggest that alternative core promoter complexes may provide a key mechanism to lock in and maintain specific transcriptional programs in terminally differentiated cell types.DOI: http://dx.doi.org/10.7554/eLife.02559.001.