The liver is a vital organ with critical functions in metabolism, protein synthesis, and immune defense. Most of the liver functions are not mature at birth and many changes happen during postnatal liver development. However, it is unclear what changes occur in liver after birth, at what developmental stages they occur, and how the developmental processes are regulated. Long non-coding RNAs (lncRNAs) are involved in organ development and cell differentiation. Here, we analyzed the transcriptome of lncRNAs in mouse liver from perinatal (day -2) to adult (day 60) by RNA-Sequencing, with an attempt to understand the role of lncRNAs in liver maturation. We found around 15,000 genes expressed, including about 2,000 lncRNAs. Most lncRNAs were expressed at a lower level than coding RNAs. Both coding RNAs and lncRNAs displayed three major ontogenic patterns: enriched at neonatal, adolescent, or adult stages. Neighboring coding and non-coding RNAs showed the trend to exhibit highly correlated ontogenic expression patterns. Gene ontology (GO) analysis revealed that some lncRNAs enriched at neonatal ages have their neighbor protein coding genes also enriched at neonatal ages and associated with cell proliferation, immune activation related processes, tissue organization pathways, and hematopoiesis; other lncRNAs enriched at adolescent ages have their neighbor protein coding genes associated with different metabolic processes. These data reveal significant functional transition during postnatal liver development and imply the potential importance of lncRNAs in liver maturation.

KAT6 histone acetyltransferases (HATs) are highly conserved in eukaryotes and have been shown to play important roles in transcriptional regulation. Here, we demonstrate that the Drosophila KAT6 Enok acetylates histone H3 Lys 23 (H3K23) in vitro and in vivo. Mutants lacking functional Enok exhibited defects in the localization of Oskar (Osk) to the posterior end of the oocyte, resulting in loss of germline formation and abdominal segments in the embryo. RNA sequencing (RNA-seq) analysis revealed that spire (spir) and maelstrom (mael), both required for the posterior localization of Osk in the oocyte, were down-regulated in enok mutants. Chromatin immunoprecipitation showed that Enok is localized to and acetylates H3K23 at the spir and mael genes. Furthermore, Gal4-driven expression of spir in the germline can largely rescue the defective Osk localization in enok mutant ovaries. Our results suggest that the Enok-mediated H3K23 acetylation (H3K23Ac) promotes the expression of spir, providing a specific mechanism linking oocyte polarization to histone modification.

Mutations that affect myelodysplasia/myeloid leukemia factor (MLF) proteins are associated with leukemia and several other cancers. However, with no strong homology to other proteins of known function, the role of MLF proteins in the cell has remained elusive. Here, we describe a proteomics approach that identifies MLF as a member of a nuclear chaperone complex containing a DnaJ protein, BCL2-associated anthanogene 2, and Hsc70. This complex associates with chromatin and regulates the expression of target genes. The MLF complex is bound to sites of nucleosome depletion and sites containing active chromatin marks (e.g., H3K4me3 and H3K4me1). Hence, MLF binding is enriched at promoters and enhancers. Additionally, the MLF-chaperone complex functions to regulate transcription factor stability, including the RUNX transcription factor involved in hematopoiesis. Although Hsc70 and other co-chaperones have been shown to play a role in nuclear translocation of a variety of proteins including transcription factors, our findings suggest that MLF and the associated co-chaperones play a direct role in modulating gene transcription.

The Swi/Snf chromatin remodeling complex functions to alter nucleosome positions by either sliding nucleosomes on DNA or the eviction of histones. The presence of histone acetylation and activator-dependent recruitment and retention of Swi/Snf is important for its efficient function. It is not understood, however, why such mechanisms are required to enhance Swi/Snf activity on nucleosomes. Snf2, the catalytic subunit of the Swi/Snf remodeling complex, has been shown to be a target of the Gcn5 acetyltransferase. Our study found that acetylation of Snf2 regulates both recruitment and release of Swi/Snf from stress-responsive genes. Also, the intramolecular interaction of the Snf2 bromodomain with the acetylated lysine residues on Snf2 negatively regulates binding and remodeling of acetylated nucleosomes by Swi/Snf. Interestingly, the presence of transcription activators mitigates the effects of the reduced affinity of acetylated Snf2 for acetylated nucleosomes. Supporting our in vitro results, we found that activator-bound genes regulating metabolic processes showed greater retention of the Swi/Snf complex even when Snf2 was acetylated. Our studies demonstrate that competing effects of (1) Swi/Snf retention by activators or high levels of histone acetylation and (2) Snf2 acetylation-mediated release regulate dynamics of Swi/Snf occupancy at target genes.

Reproductive aging is characterized by a marked decline in oocyte quality that contributes to infertility, miscarriages, and birth defects. This decline is multifactorial, and the underlying mechanisms are under active investigation. Here, we performed RNA-Seq on individual growing follicles from reproductively young and old mice to identify age-dependent functions in oocytes. This unbiased approach revealed genes involved in cellular processes known to change with age, including mitochondrial function and meiotic chromosome segregation, but also uncovered previously unappreciated categories of genes related to proteostasis and organelles required for protein metabolism. We further validated our RNA-Seq data by comparing nucleolar structure and function in oocytes from reproductively young and old mice, as this organelle is central for protein production. We examined key nucleolar markers, including upstream binding transcription factor (UBTF), an RNA polymerase I cofactor, and fibrillarin, an rRNA methyltransferase. In oocytes from mice of advanced reproductive age, UBTF was primarily expressed in giant fibrillar centers (GFCs), structures associated with high levels of rDNA transcription, and fibrillarin expression was increased ~2-fold. At the ultrastructural level, oocyte nucleoli from reproductively old mice had correspondingly more prominent fibrillar centers and dense fibrillar centers relative to young controls and more ribosomes were found in the cytoplasm. Taken together, our findings are significant because the growing oocyte is one of the most translationally active cells in the body and must accumulate high-quality maternally derived proteins to support subsequent embryo development. Thus, perturbations in protein metabolism are likely to have a profound impact on gamete health.

The KDM4 subfamily of JmjC domain-containing demethylases mediates demethylation of histone H3K36me3/me2 and H3K9me3/me2. Several studies have shown that human and yeast KDM4 proteins bind to specific gene promoters and regulate gene expression. However, the genome-wide distribution of KDM4 proteins and the mechanism of genomic-targeting remain elusive. We have previously identified Drosophila KDM4A (dKDM4A) as a histone H3K36me3 demethylase that directly interacts with HP1a. Here, we performed H3K36me3 ChIP-chip analysis in wild type and dkdm4a mutant embryos to identify genes regulated by dKDM4A demethylase activity in vivo. A subset of heterochromatic genes that show increased H3K36me3 levels in dkdm4a mutant embryos overlap with HP1a target genes. More importantly, binding to HP1a is required for dKDM4A-mediated H3K36me3 demethylation at a subset of heterochromatic genes. Collectively, these results show that HP1a functions to target the H3K36 demethylase dKDM4A to heterochromatic genes in Drosophila.

Cse4 is a variant of histone H3 that is incorporated into a single nucleosome at each centromere in budding yeast. We have discovered an E3 ubiquitin ligase, called Psh1, which controls the cellular level of Cse4 via ubiquitylation and proteolysis. The activity of Psh1 is dependent on both its RING and zinc finger domains. We demonstrate the specificity of the ubiquitylation activity of Psh1 toward Cse4 in vitro and map the sites of ubiquitylation. Mutation of key lysines prevents ubiquitylation of Cse4 by Psh1 in vitro and stabilizes Cse4 in vivo. While deletion of Psh1 stabilizes Cse4, elimination of the Cse4-specific chaperone Scm3 destabilizes Cse4, and the addition of Scm3 to the Psh1-Cse4 ubiquitylation reaction prevents Cse4 ubiquitylation, together suggesting Scm3 may protect Cse4 from ubiquitylation. Without Psh1, Cse4 overexpression is toxic and Cse4 is found at ectopic locations. Our results suggest Psh1 functions to prevent the mislocalization of Cse4.

The packaging of DNA into strings of nucleosomes is one of the features that allows eukaryotic cells to tightly regulate gene expression. The ordered disassembly of nucleosomes permits RNA polymerase II (Pol II) to access the DNA, whereas nucleosomal reassembly impedes access, thus preventing transcription and mRNA synthesis. Chromatin modifications, chromatin remodellers, histone chaperones and histone variants regulate nucleosomal dynamics during transcription. Disregulation of nucleosome dynamics results in aberrant transcription initiation, producing non-coding RNAs. Ongoing research is elucidating the molecular mechanisms that regulate chromatin structure during transcription by preventing histone exchange, thereby limiting non-coding RNA expression.

Hoxa1 has important functional roles in neural crest specification, hindbrain patterning and heart and ear development, yet the enhancers and genes that are targeted by Hoxa1 are largely unknown. In this study, we performed a comprehensive analysis of Hoxa1 target genes using genome-wide Hoxa1 binding data in mouse ES cells differentiated with retinoic acid (RA) into neural fates in combination with differential gene expression analysis in Hoxa1 gain- and loss-of-function mouse and zebrafish embryos. Our analyses reveal that Hoxa1-bound regions show epigenetic marks of enhancers, occupancy of Hox cofactors and differential expression of nearby genes, suggesting that these regions are enriched for enhancers. In support of this, 80 of them mapped to regions with known reporter activity in transgenic mouse embryos based on the Vista enhancer database. Two additional enhancers in Dok5 and Wls1 were shown to mediate neural expression in developing mouse and zebrafish. Overall, our analysis of the putative target genes indicate that Hoxa1 has input to components of major signaling pathways, including Wnt, TGF-β, Hedgehog and Hippo, and frequently does so by targeting multiple components of a pathway such as secreted inhibitors, ligands, receptors and down-stream components. We also identified genes implicated in heart and ear development, neural crest migration and neuronal patterning and differentiation, which may underlie major Hoxa1 mutant phenotypes. Finally, we found evidence for a high degree of evolutionary conservation of many binding regions and downstream targets of Hoxa1 between mouse and zebrafish. Our genome-wide analyses in ES cells suggests that we have enriched for in vivo relevant target genes and pathways associated with functional roles of Hoxa1 in mouse development.

Cse4 is the centromeric histone H3 variant in budding yeast. Psh1 is an E3 ubiquitin ligase that controls Cse4 levels through proteolysis. Here we report that Psh1 is phosphorylated by the Cka2 subunit of casein kinase 2 (CK2) to promote its E3 activity for Cse4. Deletion of CKA2 significantly stabilized Cse4. Consistent with phosphorylation promoting the activity of Psh1, Cse4 was stabilized in a Psh1 phosphodepleted mutant strain in which the major phosphorylation sites were changed to alanines. Phosphorylation of Psh1 did not control Psh1-Cse4 or Psh1-Ubc3(E2) interactions. Although Cse4 was highly stabilized in a cka2Δ strain, mislocalization of Cse4 was mild, suggesting that Cse4 misincorporation was prevented by the intact Psh1-Cse4 association. Supporting this idea, Psh1 was also stabilized in a cka2Δ strain. Collectively our data suggest that phosphorylation is crucial in Psh1-assisted control of Cse4 levels and that the Psh1-Cse4 association itself functions to prevent Cse4 misincorporation.

The chemosensory system plays an important role in orchestrating sexual behaviors in mammals. Pheromones trigger sexually dimorphic behaviors and different mouse strains exhibit differential responses to pheromone stimuli. It has been speculated that differential gene expression in the sensory organs that detect pheromones may underlie sexually-dimorphic and strain-specific responses to pheromone cues.

Maintenance of a regular chromatin structure over the coding regions of genes occurs co-transcriptionally via the 'chromatin resetting' pathway. One of the central players in this pathway is the histone methyltransferase Set2. Here we show that the loss of Set2 in yeast, Saccharomyces cerevisiae, results in transcription initiation of antisense RNAs embedded within body of protein-coding genes. These RNAs are distinct from the previously identified non-coding RNAs and cover 11% of the yeast genome. These RNA species have been named Set2-repressed antisense transcripts (SRATs) since the co-transcriptional addition of the H3K36 methyl mark by Set2 over their start sites results in their suppression. Interestingly, loss of chromatin resetting factor Set2 or the subsequent production of SRATs does not affect the abundance of the sense transcripts. This difference in transcriptional outcomes of overlapping transcripts due to a strand-independent addition of H3K36 methylation is a key regulatory feature of interleaved transcriptomes.

Hox genes modulate the properties of hematopoietic stem cells (HSCs) and reacquired Hox expression in progenitors contributes to leukemogenesis. Here, our transcriptome and DNA methylome analyses revealed that Hoxb cluster and retinoid signaling genes are predominantly enriched in LT-HSCs, and this coordinate regulation of Hoxb expression is mediated by a retinoid-dependent cis-regulatory element, distal element RARE (DERARE). Deletion of the DERARE reduced Hoxb expression, resulting in changes to many downstream signaling pathways (e.g., non-canonical Wnt signaling) and loss of HSC self-renewal and reconstitution capacity. DNA methyltransferases mediate DNA methylation on the DERARE, leading to reduced Hoxb cluster expression. Acute myeloid leukemia patients with DNMT3A mutations exhibit DERARE hypomethylation, elevated HOXB expression, and adverse outcomes. CRISPR-Cas9-mediated specific DNA methylation at DERARE attenuated HOXB expression and alleviated leukemogenesis. Collectively, these findings demonstrate pivotal roles for retinoid signaling and the DERARE in maintaining HSCs and preventing leukemogenesis by coordinate regulation of Hoxb genes.

It remains a significant challenge to define individual protein associations within networks where an individual protein can directly interact with other proteins and/or be part of large complexes, which contain functional modules. Here we demonstrate the topological scoring (TopS) algorithm for the analysis of quantitative proteomic datasets from affinity purifications. Data is analyzed in a parallel fashion where a prey protein is scored in an individual affinity purification by aggregating information from the entire dataset. Topological scores span a broad range of values indicating the enrichment of an individual protein in every bait protein purification. TopS is applied to interaction networks derived from human DNA repair proteins and yeast chromatin remodeling complexes. TopS highlights potential direct protein interactions and modules within complexes. TopS is a rapid method for the efficient and informative computational analysis of datasets, is complementary to existing analysis pipelines, and provides important insights into protein interaction networks.

Knowledge of the diverse DNA binding specificities of transcription factors is important for understanding their specific regulatory functions in animal development and evolution. We have examined the genome-wide binding properties of the mouse HOXB1 protein in embryonic stem cells differentiated into neural fates. Unexpectedly, only a small number of HOXB1 bound regions (7%) correlate with binding of the known HOX cofactors PBX and MEIS. In contrast, 22% of the HOXB1 binding peaks display co-occupancy with the transcriptional repressor REST. Analyses revealed that co-binding of HOXB1 with PBX correlates with active histone marks and high levels of expression, while co-occupancy with REST correlates with repressive histone marks and repression of the target genes. Analysis of HOXB1 bound regions uncovered enrichment of a novel 15 base pair HOXB1 binding motif HB1RE (HOXB1 response element). In vitro template binding assays showed that HOXB1, PBX1, and MEIS can bind to this motif. In vivo, this motif is sufficient for direct expression of a reporter gene and over-expression of HOXB1 selectively represses this activity. Our analyses suggest that HOXB1 has evolved an association with REST in gene regulation and the novel HB1RE motif contributes to HOXB1 function in part through a repressive role in gene expression.

The transition of hematopoiesis from the fetal liver (FL) to the bone marrow (BM) is incompletely characterized. We demonstrate that the Wiskott-Aldrich syndrome verprolin-homologous protein (WAVE) complex 2 is required for this transition, as complex degradation via deletion of its scaffold Hem-1 causes the premature exhaustion of neonatal BM hematopoietic stem cells (HSCs). This exhaustion of BM HSC is due to the failure of BM engraftment of Hem-1-/- FL HSCs, causing early death. The Hem-1-/- FL HSC engraftment defect is not due to the lack of the canonical function of the WAVE2 complex, the regulation of actin polymerization, because FL HSCs from Hem-1-/- mice exhibit no defects in chemotaxis, BM homing, or adhesion. Rather, the failure of Hem-1-/- FL HSC engraftment in the marrow is due to the loss of c-Abl survival signaling from degradation of the WAVE2 complex. However, c-Abl activity is dispensable for the engraftment of adult BM HSCs into the BM. These findings reveal a novel function of the WAVE2 complex and define a mechanism for FL HSC fitness in the embryonic BM niche.

Hoxa1 has diverse functional roles in differentiation and development. We identify and characterize properties of regions bound by HOXA1 on a genome-wide basis in differentiating mouse ES cells. HOXA1-bound regions are enriched for clusters of consensus binding motifs for HOX, PBX, and MEIS, and many display co-occupancy of PBX and MEIS. PBX and MEIS are members of the TALE family and genome-wide analysis of multiple TALE members (PBX, MEIS, TGIF, PREP1, and PREP2) shows that nearly all HOXA1 targets display occupancy of one or more TALE members. The combinatorial binding patterns of TALE proteins define distinct classes of HOXA1 targets, which may create functional diversity. Transgenic reporter assays in zebrafish confirm enhancer activities for many HOXA1-bound regions and the importance of HOX-PBX and TGIF motifs for their regulation. Proteomic analyses show that HOXA1 physically interacts on chromatin with PBX, MEIS, and PREP family members, but not with TGIF, suggesting that TGIF may have an independent input into HOXA1-bound regions. Therefore, TALE proteins appear to represent a wide repertoire of HOX cofactors, which may coregulate enhancers through distinct mechanisms. We also discover extensive auto- and cross-regulatory interactions among the Hoxa1 and TALE genes, indicating that the specificity of HOXA1 during development may be regulated though a complex cross-regulatory network of HOXA1 and TALE proteins. This study provides new insight into a regulatory network involving combinatorial interactions between HOXA1 and TALE proteins.

The clustered Hox genes, which are highly conserved across metazoans, encode homeodomain-containing transcription factors that provide a blueprint for segmental identity along the body axis. Recent studies have underscored that in addition to encoding Hox genes, the homeotic clusters contain key noncoding RNA genes that play a central role in development. In this study, we have taken advantage of genome-wide approaches to provide a detailed analysis of retinoic acid (RA)-induced transcriptional and epigenetic changes within the homeotic clusters of mouse embryonic stem cells. Although there is a general colinear response, our analyses suggest a lack of strict colinearity for several genes in the HoxA and HoxB clusters. We have identified transcribed novel noncoding RNAs (ncRNAs) and their cis-regulatory elements that function in response to RA and demonstrated that the expression of these ncRNAs from both strands represent some of the most rapidly induced transcripts in ES cells. Finally, we have provided dynamic analyses of chromatin modifications for the coding and noncoding genes expressed upon activation and suggest that active transcription can occur in the presence of chromatin modifications and machineries associated with repressed transcription state over the clusters. Overall, our data provide a resource for a better understanding of the dynamic nature of the coding and noncoding transcripts and their associated chromatin marks in the regulation of homeotic gene transcription during development.

KAT6 histone acetyltransferases (HATs) are highly conserved in eukaryotes and are involved in cell cycle regulation. However, information regarding their roles in regulating cell cycle progression is limited. Here, we report the identification of subunits of the Drosophila Enok complex and demonstrate that all subunits are important for its HAT activity. We further report a novel interaction between the Enok complex and the Elg1 proliferating cell nuclear antigen (PCNA)-unloader complex. Depletion of Enok in S2 cells resulted in a G1/S cell cycle block, and this block can be partially relieved by depleting Elg1. Furthermore, depletion of Enok reduced the chromatin-bound levels of PCNA in both S2 cells and early embryos, suggesting that the Enok complex may interact with the Elg1 complex and down-regulate its PCNA-unloading function to promote the G1/S transition. Supporting this hypothesis, depletion of Enok also partially rescued the endoreplication defects in Elg1-depleted nurse cells. Taken together, our study provides novel insights into the roles of KAT6 HATs in cell cycle regulation through modulating PCNA levels on chromatin.

The SAGA (Spt-Ada-Gcn5 acetyltransferase) complex is an important chromatin modifying complex that can both acetylate and deubiquitinate histones. Sgf29 is a novel component of the SAGA complex. Here, we report the crystal structures of the tandem Tudor domains of Saccharomyces cerevisiae and human Sgf29 and their complexes with H3K4me2 and H3K4me3 peptides, respectively, and show that Sgf29 selectively binds H3K4me2/3 marks. Our crystal structures reveal that Sgf29 harbours unique tandem Tudor domains in its C-terminus. The tandem Tudor domains in Sgf29 tightly pack against each other face-to-face with each Tudor domain harbouring a negatively charged pocket accommodating the first residue alanine and methylated K4 residue of histone H3, respectively. The H3A1 and K4me3 binding pockets and the limited binding cleft length between these two binding pockets are the structural determinants in conferring the ability of Sgf29 to selectively recognize H3K4me2/3. Our in vitro and in vivo functional assays show that Sgf29 recognizes methylated H3K4 to recruit the SAGA complex to its targets sites and mediates histone H3 acetylation, underscoring the importance of Sgf29 in gene regulation.