Integrative analysis of ChIP-chip and ChIP-seq dataset

The androgen receptor in mesenchymal progenitors regulates skeletal muscle mass via Igf1 expression in male mice

Androgens exert their effects primarily by binding to the androgen receptor (AR), a ligand-dependent nuclear receptor. While androgens have anabolic effects on skeletal muscle, previous studies reported that AR functions in myofibers to regulate skeletal muscle quality, rather than skeletal muscle mass. Therefore, the anabolic effects of androgens are exerted via nonmyofiber cells. In this context, the cellular and molecular mechanisms of AR in mesenchymal progenitors, which play a crucial role in maintaining skeletal muscle homeostasis, remain largely unknown. In this study, we demonstrated expression of AR in mesenchymal progenitors and found that targeted AR ablation in mesenchymal progenitors reduced limb muscle mass in mature adult, but not young or aged, male mice, although fatty infiltration of muscle was not affected. The absence of AR in mesenchymal progenitors led to remarkable perineal muscle hypotrophy, regardless of age, due to abnormal regulation of transcripts associated with cell death and extracellular matrix organization. Additionally, we revealed that AR in mesenchymal progenitors regulates the expression of insulin-like growth factor 1 (Igf1) and that IGF1 administration prevents perineal muscle atrophy in a paracrine manner. These findings indicate that the anabolic effects of androgens regulate skeletal muscle mass via, at least in part, AR signaling in mesenchymal progenitors.

NSD2 is a requisite subunit of the AR/FOXA1 neo-enhanceosome in promoting prostate tumorigenesis

Androgen receptor (AR) is a ligand-responsive transcription factor that drives terminal differentiation of the prostatic luminal epithelia. By contrast, in tumors originating from these cells, AR chromatin occupancy is extensively reprogrammed to activate malignant phenotypes, the molecular mechanisms of which remain unknown. Here, we show that tumor-specific AR enhancers are critically reliant on H3K36 dimethyltransferase activity of NSD2. NSD2 expression is abnormally induced in prostate cancer, where its inactivation impairs AR transactivation potential by disrupting over 65% of its cistrome. NSD2-dependent AR sites distinctively harbor the chimeric FOXA1:AR half-motif, which exclusively comprise tumor-specific AR enhancer circuitries defined from patient specimens. NSD2 inactivation also engenders increased dependency on the NSD1 paralog, and a dual NSD1/2 PROTAC degrader is preferentially cytotoxic in AR-dependent prostate cancer models. Altogether, we characterize NSD2 as an essential AR neo-enhanceosome subunit that enables its oncogenic activity, and position NSD1/2 as viable co-targets in advanced prostate cancer.

Inhibition of Renin Expression Is Regulated by an Epigenetic Switch From an Active to a Poised State

BACKGROUND

Renin-expressing cells are myoendocrine cells crucial for the maintenance of homeostasis. Renin is regulated by cAMP, p300 (histone acetyltransferase p300)/CBP (CREB-binding protein), and Brd4 (bromodomain-containing protein 4) proteins and associated pathways. However, the specific regulatory changes that occur following inhibition of these pathways are not clear.

METHODS

We treated As4.1 cells (tumoral cells derived from mouse juxtaglomerular cells that constitutively express renin) with 3 inhibitors that target different factors required for renin transcription: H-89-dihydrochloride, PKA (protein kinase A) inhibitor; JQ1, Brd4 bromodomain inhibitor; and A-485, p300/CBP inhibitor. We performed assay for transposase-accessible chromatin with sequencing (ATAC-seq), single-cell RNA sequencing, cleavage under targets and tagmentation (CUT&Tag), and chromatin immunoprecipitation sequencing for H3K27ac (acetylation of lysine 27 of the histone H3 protein) and p300 binding on biological replicates of treated and control As4.1 cells.

RESULTS

In response to each inhibitor, Ren1 expression was significantly reduced and reversible upon washout. Chromatin accessibility at the Ren1 locus did not markedly change but was globally reduced at distal elements. Inhibition of PKA led to significant reductions in H3K27ac and p300 binding specifically within the Ren1 super-enhancer region. Further, we identified enriched TF (transcription factor) motifs shared across each inhibitory treatment. Finally, we identified a set of 9 genes with putative roles across each of the 3 renin regulatory pathways and observed that each displayed differentially accessible chromatin, gene expression, H3K27ac, and p300 binding at their respective loci.

CONCLUSIONS

Inhibition of renin expression in cells that constitutively synthesize and release renin is regulated by an epigenetic switch from an active to poised state associated with decreased cell-cell communication and an epithelial-mesenchymal transition. This work highlights and helps define the factors necessary for renin cells to alternate between myoendocrine and contractile phenotypes.

Integration of shoot-derived polypeptide signals by root TGA transcription factors is essential for survival under fluctuating nitrogen environments

Unlike plants in the field, which experience significant temporal fluctuations in environmental conditions, plants in the laboratory are typically grown in controlled, stable environments. Therefore, signaling pathways evolved for survival in fluctuating environments often remain functionally latent in laboratory settings. Here, we show that TGA1 and TGA4 act as hub transcription factors through which the expression of genes involved in high-affinity nitrate uptake are regulated in response to shoot-derived phloem mobile polypeptides, CEP DOWNSTREAM 1 (CEPD1), CEPD2 and CEPD-like 2 (CEPDL2) as nitrogen (N) deficiency signals, and Glutaredoxin S1 (GrxS1) to GrxS8 as N sufficiency signals. CEPD1/2/CEPDL2 and GrxS1-S8 competitively bind to TGA1/4 in roots, with the former acting as transcription coactivators that enhance the uptake of nitrate, while the latter function as corepressor complexes together with TOPLESS (TPL), TPL-related 1 (TPR1) and TPR4 to limit nitrate uptake. Arabidopsis plants deficient in TGA1/4 maintain basal nitrate uptake and exhibit growth similar to wild-type plants in a stable N environment, but are impaired in regulation of nitrate acquisition in response to shoot N demand, leading to defective growth under fluctuating N environments where rhizosphere nitrate ions switch periodically between deficient and sufficient states. TGA1/4 are crucial transcription factors that enable plants to survive under fluctuating and challenging N environmental conditions.

C-terminally phosphorylated p27 activates self-renewal driver genes to program cancer stem cell expansion, mammary hyperplasia and cancer

In many cancers, a stem-like cell subpopulation mediates tumor initiation, dissemination and drug resistance. Here, we report that cancer stem cell (CSC) abundance is transcriptionally regulated by C-terminally phosphorylated p27 (p27pT157pT198). Mechanistically, this arises through p27 co-recruitment with STAT3/CBP to gene regulators of CSC self-renewal including MYC, the Notch ligand JAG1, and ANGPTL4. p27pTpT/STAT3 also recruits a SIN3A/HDAC1 complex to co-repress the Pyk2 inhibitor, PTPN12. Pyk2, in turn, activates STAT3, creating a feed-forward loop increasing stem-like properties in vitro and tumor-initiating stem cells in vivo. The p27-activated gene profile is over-represented in STAT3 activated human breast cancers. Furthermore, mammary transgenic expression of phosphomimetic, cyclin-CDK-binding defective p27 (p27CK-DD) increases mammary duct branching morphogenesis, yielding hyperplasia and microinvasive cancers that can metastasize to liver, further supporting a role for p27pTpT in CSC expansion. Thus, p27pTpT interacts with STAT3, driving transcriptional programs governing stem cell expansion or maintenance in normal and cancer tissues.

p300/CBP degradation is required to disable the active AR enhanceosome in prostate cancer

Prostate cancer is an exemplar of an enhancer-binding transcription factor-driven disease. The androgen receptor (AR) enhanceosome complex comprised of chromatin and epigenetic coregulators assembles at enhancer elements to drive disease progression. The paralog lysine acetyltransferases p300 and CBP deposit histone marks that are associated with enhancer activation. Here, we demonstrate that p300/CBP are determinant cofactors of the active AR enhanceosome in prostate cancer. Histone H2B N-terminus multisite lysine acetylation (H2BNTac), which is exclusively reliant on p300/CBP catalytic function, marked active enhancers and was notably elevated in prostate cancer lesions relative to the adjacent benign epithelia. Degradation of p300/CBP rapidly depleted acetylation marks associated with the active AR enhanceosome, which was only partially phenocopied by inhibition of their reader bromodomains. Notably, H2BNTac was effectively abrogated only upon p300/CBP degradation, which led to a stronger suppression of p300/CBP-dependent oncogenic gene programs relative to bromodomain inhibition or the inhibition of its catalytic domain. In vivo experiments using an orally active p300/CBP proteolysis targeting chimera (PROTAC) degrader (CBPD-409) showed that p300/CBP degradation potently inhibited tumor growth in preclinical models of castration-resistant prostate cancer and synergized with AR antagonists. While mouse p300/CBP orthologs were effectively degraded in host tissues, prolonged treatment with the PROTAC degrader was well tolerated with no significant signs of toxicity. Taken together, our study highlights the pivotal role of p300/CBP in maintaining the active AR enhanceosome and demonstrates how target degradation may have functionally distinct effects relative to target inhibition, thus supporting the development of p300/CBP degraders for the treatment of advanced prostate cancer.

Human ASXL1-Mutant Hematopoiesis Is Driven by a Truncated Protein Associated with Aberrant Deubiquitination of H2AK119

Mutations in additional sex combs like 1 (ASXL1) confer poor prognosis both in myeloid malignancies and in premalignant clonal hematopoiesis (CH). However, the mechanisms by which these mutations contribute to disease initiation remain unresolved, and mutation-specific targeting has remained elusive. To address this, we developed a human disease model that recapitulates the disease trajectory from ASXL1-mutant CH to lethal myeloid malignancy. We demonstrate that mutations in ASXL1 lead to the expression of a functional, truncated protein and determine that truncated ASXL1 leads to global redistribution of the repressive chromatin mark H2AK119Ub, increased transposase-accessible chromatin, and activation of both myeloid and stem cell gene-expression programs. Finally, we demonstrate that H2AK119Ub levels are tied to truncated ASXL1 expression levels and leverage this observation to demonstrate that inhibition of the PRC1 complex might be an ASXL1-mutant-specific therapeutic vulnerability in both premalignant CH and myeloid malignancy.

SIGNIFICANCE

Mutant ASXL1 is a common driver of CH and myeloid malignancy. Using primary human HSPCs, we determine that truncated ASXL1 leads to redistribution of H2AK119Ub and may affect therapeutic vulnerability to PRC1 inhibition.

NSD2 is a requisite subunit of the AR/FOXA1 neo-enhanceosome in promoting prostate tumorigenesis

The androgen receptor (AR) is a ligand-responsive transcription factor that binds at enhancers to drive terminal differentiation of the prostatic luminal epithelia. By contrast, in tumors originating from these cells, AR chromatin occupancy is extensively reprogrammed to drive hyper-proliferative, metastatic, or therapy-resistant phenotypes, the molecular mechanisms of which remain poorly understood. Here, we show that the tumor-specific enhancer circuitry of AR is critically reliant on the activity of Nuclear Receptor Binding SET Domain Protein 2 (NSD2), a histone 3 lysine 36 di-methyltransferase. NSD2 expression is abnormally gained in prostate cancer cells and its functional inhibition impairs AR trans-activation potential through partial off-loading from over 40,000 genomic sites, which is greater than 65% of the AR tumor cistrome. The NSD2-dependent AR sites distinctly harbor a chimeric AR-half motif juxtaposed to a FOXA1 element. Similar chimeric motifs of AR are absent at the NSD2-independent AR enhancers and instead contain the canonical palindromic motifs. Meta-analyses of AR cistromes from patient tumors uncovered chimeric AR motifs to exclusively participate in tumor-specific enhancer circuitries, with a minimal role in the physiological activity of AR. Accordingly, NSD2 inactivation attenuated hallmark cancer phenotypes that were fully reinstated upon exogenous NSD2 re-expression. Inactivation of NSD2 also engendered increased dependency on its paralog NSD1, which independently maintained AR and MYC hyper-transcriptional programs in cancer cells. Concordantly, a dual NSD1/2 PROTAC degrader, called LLC0150, was preferentially cytotoxic in AR-dependent prostate cancer as well as NSD2-altered hematologic malignancies. Altogether, we identify NSD2 as a novel subunit of the AR neo-enhanceosome that wires prostate cancer gene expression programs, positioning NSD1/2 as viable paralog co-targets in advanced prostate cancer.

VGLL1 cooperates with TEAD4 to control human trophectoderm lineage specification

In contrast to rodents, the mechanisms underlying human trophectoderm and early placenta specification are understudied due to ethical barriers and the scarcity of embryos. Recent reports have shown that human pluripotent stem cells (PSCs) can differentiate into trophectoderm (TE)-like cells (TELCs) and trophoblast stem cells (TSCs), offering a valuable in vitro model to study early placenta specification. Here, we demonstrate that the VGLL1 (vestigial-like family member 1), which is highly expressed during human and non-human primate TE specification in vivo but is negligibly expressed in mouse, is a critical regulator of cell fate determination and self-renewal in human TELCs and TSCs derived from naïve PSCs. Mechanistically, VGLL1 partners with the transcription factor TEAD4 (TEA domain transcription factor 4) to regulate chromatin accessibility at target gene loci through histone acetylation and acts in cooperation with GATA3 and TFAP2C. Our work is relevant to understand primate early embryogenesis and how it differs from other mammalian species.

Transcription factor-mediated direct cellular reprogramming yields cell-type specific DNA methylation signature

Direct reprogramming, inducing the conversion of one type of somatic cell into another by the forced expression of defined transcription factors, is a technology with anticipated medical applications. However, due to the many unresolved aspects of the induction mechanisms, it is essential to thoroughly analyze the epigenomic state of the generated cells. Here, we performed comparative genome-wide DNA methylation analyses of mouse embryonic fibroblasts (MEFs) and cells composing organoids formed by intestinal stem cells (ISCs) or induced ISCs (iISCs) that were directly induced from MEFs. We found that the CpG methylation state was similar between cells forming ISC organoids and iISC organoids, while they differed widely from those in MEFs. Moreover, genomic regions that were differentially methylated between ISC organoid- and iISC organoid-forming cells did not significantly affect gene expression. These results demonstrate the accuracy and safety of iISC induction, leading to the medical applications of this technology.

FoxK1 associated gene regulatory network in hepatic insulin action and its relationship to FoxO1 and insulin receptor mediated transcriptional regulation

OBJECTIVE

Insulin acts on the liver via changes in gene expression to maintain glucose and lipid homeostasis. This study aimed to the Forkhead box protein K1 (FOXK1) associated gene regulatory network as a transcriptional regulator of hepatic insulin action and to determine its role versus FoxO1 and possible actions of the insulin receptor at the DNA level.

METHODS

Genome-wide analysis of FoxK1 binding were studied by chromatin immunoprecipitation sequencing and compared to those for IR and FoxO1. These were validated by knockdown experiments and gene expression analysis.

RESULTS

Chromatin immunoprecipitation (ChIP) sequencing shows that FoxK1 binds to the proximal promoters and enhancers of over 4000 genes, and insulin enhances this interaction for about 75% of them. These include genes involved in cell cycle, senescence, steroid biosynthesis, autophagy, and metabolic regulation, including glucose metabolism and mitochondrial function and are enriched in a TGTTTAC consensus motif. Some of these genes are also bound by FoxO1. Comparing this FoxK1 ChIP-seq data to that of the insulin receptor (IR) reveals that FoxK1 may act as the transcription factor partner for some of the previously reported roles of IR in gene regulation, including for LARS1 and TIMM22, which are involved in rRNA processing and cell cycle.

CONCLUSION

These data demonstrate that FoxK1 is an important regulator of gene expression in response to insulin in liver and may act in concert with FoxO1 and IR in regulation of genes in metabolism and other important biological pathways.

Range of chromatin accessibility configurations are permissive of GABAergic fate acquisition in developing mouse brain

In recent single-cell -omics studies, both the differential activity of transcription factors regulating cell fate determination and differential genome activation have been tested for utility as descriptors of cell types. Naturally, genome accessibility and gene expression are interlinked. To understand the variability in genomic feature activation in the GABAergic neurons of different spatial origins, we have mapped accessible chromatin regions and mRNA expression in single cells derived from the developing mouse central nervous system (CNS). We first defined a reference set of open chromatin regions for scATAC-seq read quantitation across samples, allowing comparison of chromatin accessibility between brain regions and cell types directly. Second, we integrated the scATAC-seq and scRNA-seq data to form a unified resource of transcriptome and chromatin accessibility landscape for the cell types in di- and telencephalon, midbrain and anterior hindbrain of E14.5 mouse embryo. Importantly, we implemented resolution optimization at the clustering, and automatized the cell typing step. We show high level of concordance between the cell clustering based on the chromatin accessibility and the transcriptome in analyzed neuronal lineages, indicating that both genome and transcriptome features can be used for cell type definition. Hierarchical clustering by the similarity in accessible chromatin reveals that the genomic feature activation correlates with neurotransmitter phenotype, selector gene expression, cell differentiation stage and neuromere origins.

The lncRNA Sweetheart regulates compensatory cardiac hypertrophy after myocardial injury in murine males

After myocardial infarction in the adult heart the remaining, non-infarcted tissue adapts to compensate the loss of functional tissue. This adaptation requires changes in gene expression networks, which are mostly controlled by transcription regulating proteins. Long non-coding transcripts (lncRNAs) are taking part in fine-tuning such gene programs. We describe and characterize the cardiomyocyte specific lncRNA Sweetheart RNA (Swhtr), an approximately 10 kb long transcript divergently expressed from the cardiac core transcription factor coding gene Nkx2-5. We show that Swhtr is dispensable for normal heart development and function but becomes essential for the tissue adaptation process after myocardial infarction in murine males. Re-expressing Swhtr from an exogenous locus rescues the Swhtr null phenotype. Genes that depend on Swhtr after cardiac stress are significantly occupied and therefore most likely regulated by NKX2-5. The Swhtr transcript interacts with NKX2-5 and disperses upon hypoxic stress in cardiomyocytes, indicating an auxiliary role of Swhtr for NKX2-5 function in tissue adaptation after myocardial injury.

Efficient genome editing in erythroid cells unveils novel MYB target genes and regulatory functions

Targeted genome editing holds great promise in biology. However, efficient genome modification, including gene knock-in (KI), remains an unattained goal in multiple cell types and loci due to poor transfection efficiencies and low target genes expression, impeding the positive selection of recombined cells. Here, we describe a genome editing approach to achieve efficient gene targeting using hard to transfect erythroid cell lines. We demonstrate robust fluorescent protein KI efficiency in low expressed transcription factor (TF) genes (e.g., Myb or Zeb1). We further show the ability to target two independent loci in individual cells, exemplified by MYB-GFP and NuMA-Cherry double KI, allowing multicolor labeling of regulatory factors at physiological endogenous levels. Our KI tagging approach allowed us to perform genome-wide TF analysis at increased signal-to-noise ratios, and highlighted previously unidentified MYB target genes and pathways. Overall, we establish a versatile CRISPR-Cas9-based platform, offering attractive opportunities for the dissection of the erythroid differentiation process.

Loss of TAZ after YAP deletion severely impairs foregut development and worsens cholestatic hepatocellular injury

BACKGROUND

We previously showed that loss of yes-associated protein 1 (YAP) in early liver development (YAPKO) leads to an Alagille syndrome-like phenotype, with failure of intrahepatic bile duct development, severe cholestasis, and chronic hepatocyte adaptations to reduce liver injury. TAZ, a paralog of YAP, was significantly upregulated in YAPKO hepatocytes and interacted with TEA domain family member (TEAD) transcription factors, suggesting possible compensatory activity.

METHODS

We deleted both Yap1 and Wwtr1 (which encodes TAZ) during early liver development using the Foxa3 promoter to drive Cre expression, similar to YAPKO mice, resulting in YAP/TAZ double knockout (DKO) and YAPKO with TAZ heterozygosity (YAPKO TAZHET). We evaluated these mice using immunohistochemistry, serum biochemistry, bile acid profiling, and RNA sequencing.

RESULTS

DKO mice were embryonic lethal, but their livers were similar to YAPKO, suggesting an extrahepatic cause of death. Male YAPKO TAZHET mice were also embryonic lethal, with insufficient samples to determine the cause. However, YAPKO TAZHET females survived and were phenotypically similar to YAPKO mice, with increased bile acid hydrophilicity and similar global gene expression adaptations but worsened the hepatocellular injury. TAZ heterozygosity in YAPKO impacted the expression of canonical YAP targets Ctgf and Cyr61, and we found changes in pathways regulating cell division and inflammatory signaling correlating with an increase in hepatocyte cell death, cell cycling, and macrophage recruitment.

CONCLUSIONS

YAP loss (with or without TAZ loss) aborts biliary development. YAP and TAZ play a codependent critical role in foregut endoderm development outside the liver, but they are not essential for hepatocyte development. TAZ heterozygosity in YAPKO livers increased cell cycling and inflammatory signaling in the setting of chronic injury, highlighting genes that are especially sensitive to TAZ regulation.

Loss of ANCO1 Expression Regulates Chromatin Accessibility and Drives Progression of Early-Stage Triple-Negative Breast Cancer

Mutations in the gene ankyrin repeat domain containing 11 (ANKRD11/ANCO1) play a role in neurodegenerative disorders, and its loss of heterozygosity and low expression are seen in some cancers. Here, we show that low ANCO1 mRNA and protein expression levels are prognostic markers for poor clinical outcomes in breast cancer and that loss of nuclear ANCO1 protein expression predicts lower overall survival of patients with triple-negative breast cancer (TNBC). Knockdown of ANCO1 in early-stage TNBC cells led to aneuploidy, cellular senescence, and enhanced invasion in a 3D matrix. The presence of a subpopulation of ANCO1-depleted cells enabled invasion of the overall cell population in vitro and they converted more rapidly to invasive lesions in a xenograft mouse model. In ANCO1-depleted cells, ChIP-seq analysis showed a global increase in H3K27Ac signals that were enriched for AP-1, TEAD, STAT3, and NFκB motifs. ANCO1-regulated H3K27Ac peaks had a significantly higher overlap with known breast cancer enhancers compared to ANCO1-independent ones. H3K27Ac engagement was associated with transcriptional activation of genes in the PI3K-AKT, epithelial-mesenchymal transition (EMT), and senescence pathways. In conclusion, ANCO1 has hallmarks of a tumor suppressor whose loss of expression activates breast-cancer-specific enhancers and oncogenic pathways that can accelerate the early-stage progression of breast cancer.

ChIP-seq profiling of H3K4me3 and H3K27me3 in an invasive insect, Bactrocera dorsalis

Introduction: While it has been suggested that histone modifications can facilitate animal responses to rapidly changing environments, few studies have profiled whole-genome histone modification patterns in invasive species, leaving the regulatory landscape of histone modifications in invasive species unclear. Methods: Here, we screen genome-wide patterns of two important histone modifications, trimethylated Histone H3 Lysine 4 (H3K4me3) and trimethylated Histone H3 Lysine 27 (H3K27me3), in adult thorax muscles of a notorious invasive pest, the Oriental fruit fly Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), using Chromatin Immunoprecipitation with high-throughput sequencing (ChIP-seq). Results: We identified promoters featured by the occupancy of H3K4me3, H3K27me3 or bivalent histone modifications that were respectively annotated with unique genes key to muscle development and structure maintenance. In addition, we found H3K27me3 occupied the entire body of genes, where the average enrichment was almost constant. Transcriptomic analysis indicated that H3K4me3 is associated with active gene transcription, and H3K27me3 is mostly associated with transcriptional repression. Importantly, we identified genes and putative motifs modified by distinct histone modification patterns that may possibly regulate flight activity. Discussion: These findings provide the first evidence of histone modification signature in B. dorsalis, and will be useful for future studies of epigenetic signature in other invasive insect species.

The role of plastic and evolved DNA methylation in parallel adaptation of threespine stickleback (Gasterosteus aculeatus)

Repeated phenotypic patterns among populations undergoing parallel evolution in similar environments provide support for the deterministic role of natural selection. Epigenetic modifications can mediate plastic and evolved phenotypic responses to environmental change and might make important contributions to parallel adaptation. While many studies have explored the genetic basis of repeated phenotypic divergence, the role of epigenetic processes during parallel adaptation remains unclear. The parallel evolution of freshwater ecotypes of threespine stickleback fish (Gasterosteus aculeatus) following colonization of thousands of lakes and streams from the ocean is a classic example of parallel phenotypic and genotypic adaptation. To investigate epigenetic modifications during parallel adaptation of threespine stickleback, we reanalysed three independent data sets that investigated DNA methylation variation between marine and freshwater ecotypes. Although we found widespread methylation differentiation between ecotypes, there was no significant tendency for CpG sites associated with repeated methylation differentiation across studies to be parallel versus nonparallel. To next investigate the role of plastic versus evolved changes in methylation during freshwater adaptation, we explored if CpG sites exhibiting methylation plasticity during salinity change were more likely to also show evolutionary divergence in methylation between ecotypes. The directions of divergence between ecotypes were generally in the opposite direction to those observed for plasticity when ecotypes were challenged with non-native salinity conditions, suggesting that most plastic responses are likely to be maladaptive during colonization of new environments. Finally, we found a greater number of CpG sites showing evolved changes when ancestral marine ecotypes are acclimated to freshwater environments, whereas plastic changes predominate when derived freshwater ecotypes transition back to their ancestral marine environments. These findings provide evidence for an epigenetic contribution to parallel adaptation and demonstrate the contrasting roles of plastic and evolved methylation differences during adaptation to new environments.

Identification of Heme Oxygenase-1 as a Putative DNA-Binding Protein

Heme oxygenase-1 (HO-1) is a rate-limiting enzyme in degrading heme into biliverdin and iron. HO-1 can also enter the nucleus and regulate gene transcription independent of its enzymatic activity. Whether HO-1 can alter gene expression through direct binding to target DNA remains unclear. Here, we performed HO-1 CHIP-seq and then employed 3D structural modeling to reveal putative HO-1 DNA binding domains. We identified three probable DNA binding domains on HO-1. Using the Proteinarium, we identified several genes as the most highly connected nodes in the interactome among the HO-1 gene binding targets. We further demonstrated that HO-1 modulates the expression of these key genes using Hmox1 deficient cells. Finally, mutation of four conserved amino acids (E215, I211, E201, and Q27) within HO-1 DNA binding domain 1 significantly increased expression of Gtpbp3 and Eif1 genes that were identified within the top 10 binding hits normalized by gene length predicted to bind this domain. Based on these data, we conclude that HO-1 protein is a putative DNA binding protein, and regulates targeted gene expression. This provides the foundation for developing specific inhibitors or activators targeting HO-1 DNA binding domains to modulate targeted gene expression and corresponding cellular function.

Biphasic cell cycle defect causes impaired neurogenesis in down syndrome

Impaired neurogenesis in Down syndrome (DS) is characterized by reduced neurons, increased glial cells, and delayed cortical lamination. However, the underlying cause for impaired neurogenesis in DS is not clear. Using both human and mouse iPSCs, we demonstrate that DS impaired neurogenesis is due to biphasic cell cycle dysregulation during the generation of neural progenitors from iPSCs named the “neurogenic stage” of neurogenesis. Upon neural induction, DS cells showed reduced proliferation during the early phase followed by increased proliferation in the late phase of the neurogenic stage compared to control cells. While reduced proliferation in the early phase causes reduced neural progenitor pool, increased proliferation in the late phase leads to delayed post mitotic neuron generation in DS. RNAseq analysis of late-phase DS progenitor cells revealed upregulation of S phase-promoting regulators, Notch, Wnt, Interferon pathways, and REST, and downregulation of several genes of the BAF chromatin remodeling complex. NFIB and POU3F4, neurogenic genes activated by the interaction of PAX6 and the BAF complex, were downregulated in DS cells. ChIPseq analysis of late-phase neural progenitors revealed aberrant PAX6 binding with reduced promoter occupancy in DS cells. Together, these data indicate that impaired neurogenesis in DS is due to biphasic cell cycle dysregulation during the neurogenic stage of neurogenesis.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

Distinct Roles of NANOS1 and NANOS3 in the Cell Cycle and NANOS3-PUM1-FOXM1 Axis to Control G2/M Phase in a Human Primordial Germ Cell Model

Nanos RNA-binding proteins are critical factors of germline development throughout the animal kingdom and their dysfunction causes infertility. During evolution, mammalian Nanos paralogues adopted divergent roles in germ cell biology. However, the molecular basis behind this divergence, such as their target mRNAs, remains poorly understood. Our RNA-sequencing analysis in a human primordial germ cell model-TCam-2 cell line revealed distinct pools of genes involved in the cell cycle process downregulated upon NANOS1 and NANOS3 overexpression. We show that NANOS1 and NANOS3 proteins influence different stages of the cell cycle. Namely, NANOS1 is involved in the G1/S and NANOS3 in the G2/M phase transition. Many of their cell cycle targets are known infertility and cancer-germ cell genes. Moreover, NANOS3 in complex with RNA-binding protein PUM1 causes 3′UTR-mediated repression of FOXM1 mRNA encoding a transcription factor crucial for G2/M phase transition. Interestingly, while NANOS3 and PUM1 act as post-transcriptional repressors of FOXM1, FOXM1 potentially acts as a transcriptional activator of NANOS3, PUM1, and itself. Finally, by utilizing publicly available RNA-sequencing datasets, we show that the balance between FOXM1-NANOS3 and FOXM1-PUM1 expression levels is disrupted in testis cancer, suggesting a potential role in this disease.

SMYD5 catalyzes histone H3 lysine 36 trimethylation at promoters

Histone marks, carriers of epigenetic information, regulate gene expression. In mammalian cells, H3K36me3 is mainly catalyzed by SETD2 at gene body regions. Here, we find that in addition to gene body regions, H3K36me3 is enriched at promoters in primary cells. Through screening, we identify SMYD5, which is recruited to chromatin by RNA polymerase II, as a methyltransferase catalyzing H3K36me3 at promoters. The enzymatic activity of SMYD5 is dependent on its C-terminal glutamic acid-rich domain. Overexpression of full-length Smyd5, but not the C-terminal domain-truncated Smyd5, restores H3K36me3 at promoters in Smyd5 knockout cells. Furthermore, elevated Smyd5 expression contributes to tumorigenesis in liver hepatocellular carcinoma. Together, our findings identify SMYD5 as the H3K36me3 methyltransferase at promoters that regulates gene expression, providing insights into the localization and function of H3K36me3.

ALGAEFUN with MARACAS, microALGAE FUNctional enrichment tool for MicroAlgae RnA-seq and Chip-seq AnalysiS

Background

Microalgae are emerging as promising sustainable sources for biofuels, biostimulants in agriculture, soil bioremediation, feed and human nutrients. Nonetheless, the molecular mechanisms underpinning microalgae physiology and the biosynthesis of compounds of biotechnological interest are largely uncharacterized. This hinders the development of microalgae full potential as cell-factories. The recent application of omics technologies into microalgae research aims at unraveling these systems. Nevertheless, the lack of specific tools for analysing omics raw data generated from microalgae to provide biological meaningful information are hampering the impact of these technologies. The purpose of ALGAEFUN with MARACAS consists in providing researchers in microalgae with an enabling tool that will allow them to exploit transcriptomic and cistromic high-throughput sequencing data.

Results

ALGAEFUN with MARACAS consists of two different tools. First, MARACAS (MicroAlgae RnA-seq and Chip-seq AnalysiS) implements a fully automatic computational pipeline receiving as input RNA-seq (RNA sequencing) or ChIP-seq (chromatin immunoprecipitation sequencing) raw data from microalgae studies. MARACAS generates sets of differentially expressed genes or lists of genomic loci for RNA-seq and ChIP-seq analysis respectively. Second, ALGAEFUN (microALGAE FUNctional enrichment tool) is a web-based application where gene sets generated from RNA-seq analysis as well as lists of genomic loci from ChIP-seq analysis can be used as input. On the one hand, it can be used to perform Gene Ontology and biological pathways enrichment analysis over gene sets. On the other hand, using the results of ChIP-seq data analysis, it identifies a set of potential target genes and analyses the distribution of the loci over gene features. Graphical representation of the results as well as tables with gene annotations are generated and can be downloaded for further analysis.

Conclusions

ALGAEFUN with MARACAS provides an integrated environment for the microalgae research community that facilitates the process of obtaining relevant biological information from raw RNA-seq and ChIP-seq data. These applications are designed to assist researchers in the interpretation of gene lists and genomic loci based on functional enrichment analysis. ALGAEFUN with MARACAS is publicly available on https://greennetwork.us.es/AlgaeFUN/.

Pan-AMPK activator O304 prevents gene expression changes and remobilisation of histone marks in islets of diet-induced obese mice

AMP-activated protein kinase (AMPK) has an important role in cellular energy homeostasis and has emerged as a promising target for treatment of Type 2 Diabetes (T2D) due to its beneficial effects on insulin sensitivity and glucose homeostasis. O304 is a pan-AMPK activator that has been shown to improve glucose homeostasis in both mouse models of diabetes and in human T2D subjects. Here, we describe the genome-wide transcriptional profile and chromatin landscape of pancreatic islets following O304 treatment of mice fed high-fat diet (HFD). O304 largely prevented genome-wide gene expression changes associated with HFD feeding in CBA mice and these changes were associated with remodelling of active and repressive chromatin marks. In particular, the increased expression of the β-cell stress marker Aldh1a3 in islets from HFD-mice is completely abrogated following O304 treatment, which is accompanied by loss of active chromatin marks in the promoter as well as distant non-coding regions upstream of the Aldh1a3 gene. Moreover, O304 treatment restored dysfunctional glucose homeostasis as well as expression of key markers associated with β-cell function in mice with already established obesity. Our findings provide preclinical evidence that O304 is a promising therapeutic compound not only for T2D remission but also for restoration of β-cell function following remission of T2D diabetes.

Master lineage transcription factors anchor trans mega transcriptional complexes at highly accessible enhancer sites to promote long-range chromatin clustering and transcription of distal target genes

The term 'super enhancers' (SE) has been widely used to describe stretches of closely localized enhancers that are occupied collectively by large numbers of transcription factors (TFs) and co-factors, and control the transcription of highly-expressed genes. Through integrated analysis of >600 DNase-seq, ChIP-seq, GRO-seq, STARR-seq, RNA-seq, Hi-C and ChIA-PET data in five human cancer cell lines, we identified a new class of autonomous SEs (aSEs) that are excluded from classic SE calls by the widely used Rank Ordering of Super-Enhancers (ROSE) method. TF footprint analysis revealed that compared to classic SEs and regular enhancers, aSEs are tightly bound by a dense array of master lineage TFs, which serve as anchors to recruit additional TFs and co-factors in trans. In addition, aSEs are preferentially enriched for Cohesins, which likely involve in stabilizing long-distance interactions between aSEs and their distal target genes. Finally, we showed that aSEs can be reliably predicted using a single DNase-seq data or combined with Mediator and/or P300 ChIP-seq. Overall, our study demonstrates that aSEs represent a unique class of functionally important enhancer elements that distally regulate the transcription of highly expressed genes.

Integrated Analysis of Transcriptome and Histone Modifications in Granulosa Cells During Ovulation in Female Mice

The ovulatory luteinizing hormone (LH) surge induces rapid changes of gene expression and cellular functions in granulosa cells (GCs) undergoing luteinization. However, it remains unclear how the changes in genome-wide gene expression are regulated. H3K4me3 histone modifications are involved in the rapid alteration of gene expression. In this study, we investigated genome-wide changes of transcriptome and H3K4me3 status in mouse GCs undergoing luteinization. GCs were obtained from mice treated with equine chorionic gonadotropin (hCG) before, 4 hours, and 12 hours after human chorionic gonadotropin injection. RNA-sequencing identified a number of upregulated and downregulated genes, which could be classified into 8 patterns according to the time-course changes of gene expression. Many genes were transiently upregulated or downregulated at 4 hours after hCG stimulation. Gene Ontology terms associated with these genes included steroidogenesis, ovulation, cumulus-oocyte complex (COC) expansion, angiogenesis, immune system, reactive oxygen species (ROS) metabolism, inflammatory response, metabolism, and autophagy. The cellular functions of DNA repair and cell growth were newly identified as being activated during ovulation. Chromatin immunoprecipitation-sequencing revealed a genome-wide and rapid change in H3K4me3 during ovulation. Integration of transcriptome and H3K4me3 data identified many H3K4me3-associated genes that are involved in steroidogenesis, ovulation, COC expansion, angiogenesis, inflammatory response, immune system, ROS metabolism, lipid and glucose metabolism, autophagy, and regulation of cell size. The present results suggest that genome-wide changes in H3K4me3 after the LH surge are associated with rapid changes in gene expression in GCs, which enables GCs to acquire a lot of cellular functions within a short time that are required for ovulation and luteinization.

PRMT6 activates cyclin D1 expression in conjunction with the transcription factor LEF1

The establishment of cell type specific gene expression by transcription factors and their epigenetic cofactors is central for cell fate decisions. Protein arginine methyltransferase 6 (PRMT6) is an epigenetic regulator of gene expression mainly through methylating arginines at histone H3. This way it influences cellular differentiation and proliferation. PRMT6 lacks DNA-binding capability but is recruited by transcription factors to regulate gene expression. However, currently only a limited number of transcription factors have been identified, which facilitate recruitment of PRMT6 to key cell cycle related target genes. Here, we show that LEF1 contributes to the recruitment of PRMT6 to the central cell cycle regulator CCND1 (Cyclin D1). We identified LEF1 as an interaction partner of PRMT6. Knockdown of LEF1 or PRMT6 reduces CCND1 expression. This is in line with our observation that knockdown of PRMT6 increases the number of cells in G1 phase of the cell cycle and decreases proliferation. These results improve the understanding of PRMT6 activity in cell cycle regulation. We expect that these insights will foster the rational development and usage of specific PRMT6 inhibitors for cancer therapy.

The microbiota is dispensable for the early stages of peripheral regulatory T cell induction within mesenteric lymph nodes

Intestinal Foxp3+ regulatory T cell (Treg) subsets are crucial players in tolerance to microbiota-derived and food-borne antigens, and compelling evidence suggests that the intestinal microbiota modulates their generation, functional specialization, and maintenance. Selected bacterial species and microbiota-derived metabolites, such as short-chain fatty acids (SCFAs), have been reported to promote Treg homeostasis in the intestinal lamina propria. Furthermore, gut-draining mesenteric lymph nodes (mLNs) are particularly efficient sites for the generation of peripherally induced Tregs (pTregs). Despite this knowledge, the direct role of the microbiota and their metabolites in the early stages of pTreg induction within mLNs is not fully elucidated. Here, using an adoptive transfer-based pTreg induction system, we demonstrate that neither transfer of a dysbiotic microbiota nor dietary SCFA supplementation modulated the pTreg induction capacity of mLNs. Even mice housed under germ-free (GF) conditions displayed equivalent pTreg induction within mLNs. Further molecular characterization of these de novo induced pTregs from mLNs by dissection of their transcriptomes and accessible chromatin regions revealed that the microbiota indeed has a limited impact and does not contribute to the initialization of the Treg-specific epigenetic landscape. Overall, our data suggest that the microbiota is dispensable for the early stages of pTreg induction within mLNs.

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.

Deficient spermiogenesis in mice lacking Rlim

The X-linked gene Rlim plays major roles in female mouse development and reproduction, where it is crucial for the maintenance of imprinted X chromosome inactivation in extraembryonic tissues of embryos. However, while females carrying a systemic Rlim knockout (KO) die around implantation, male Rlim KO mice appear healthy and are fertile. Here, we report an important role for Rlim in testis where it is highly expressed in post-meiotic round spermatids as well as in Sertoli cells. Systemic deletion of the Rlim gene results in lower numbers of mature sperm that contains excess cytoplasm, leading to decreased sperm motility and in vitro fertilization rates. Targeting the conditional Rlim cKO specifically to the spermatogenic cell lineage largely recapitulates this phenotype. These results reveal functions of Rlim in male reproduction specifically in round spermatids during spermiogenesis.

Aging of preleukemic thymocytes drives CpG island hypermethylation in T-cell acute lymphoblastic leukemia

Cancer cells display DNA hypermethylation at specific CpG islands in comparison to their normal healthy counterparts, but the mechanism that drives this so-called CpG island methylator phenotype (CIMP) remains poorly understood. Here, we show that CpG island methylation in human T-cell acute lymphoblastic leukemia (T-ALL) mainly occurs at promoters of Polycomb Repressor Complex 2 (PRC2) target genes that are not expressed in normal or malignant T-cells and which display a reciprocal association with H3K27me3 binding. In addition, we revealed that this aberrant methylation profile reflects the epigenetic history of T-ALL and is established already in pre-leukemic, self-renewing thymocytes that precede T-ALL development. Finally, we unexpectedly uncover that this age-related CpG island hypermethylation signature in T-ALL is completely resistant to the FDA-approved hypomethylating agent Decitabine. Altogether, we here provide conceptual evidence for the involvement of a pre-leukemic phase characterized by self-renewing thymocytes in the pathogenesis of human T-ALL.

Globally altered epigenetic landscape and delayed osteogenic differentiation in H3.3-G34W-mutant giant cell tumor of bone

The neoplastic stromal cells of giant cell tumor of bone (GCTB) carry a mutation in H3F3A, leading to a mutant histone variant, H3.3-G34W, as a sole recurrent genetic alteration. We show that in patient-derived stromal cells H3.3-G34W is incorporated into the chromatin and associates with massive epigenetic alterations on the DNA methylation, chromatin accessibility and histone modification level, that can be partially recapitulated in an orthogonal cell line system by the introduction of H3.3-G34W. These epigenetic alterations affect mainly heterochromatic and bivalent regions and provide possible explanations for the genomic instability, as well as the osteolytic phenotype of GCTB. The mutation occurs in differentiating mesenchymal stem cells and associates with an impaired osteogenic differentiation. We propose that the observed epigenetic alterations reflect distinct differentiation stages of H3.3 WT and H3.3 MUT stromal cells and add to H3.3-G34W-associated changes.

The histone variant mutation H3.3-G34W occurs in the majority of giant cell tumor of bone (GCTB). By profiling patient-derived GCTB tumor cells, the authors show that this mutation associates with epigenetic alterations in heterochromatic and bivalent regions that contribute to an impaired osteogenic differentiation and the osteolytic phenotype of GCTB.

YAP/TAZ Activation Drives Uveal Melanoma Initiation and Progression

Uveal melanoma (UM), the most common ocular malignancy, is characterized by GNAQ/11 mutations. Hippo/YAP and Ras/mitogen-activated protein kinase (MAPK) emerge as two important signaling pathways downstream of G protein alpha subunits of the Q class (GαQ/11)-mediated transformation, although whether and how they contribute to UM genesis in vivo remain unclear. Here, we adapt an adeno-associated virus (AAV)-based ocular injection method to directly deliver Cre recombinase into the mouse uveal tract and demonstrate that Lats1/2 kinases suppress UM formation specifically in uveal melanocytes. We find that genetic activation of YAP, but not Kras, is sufficient to initiate UM. We show that YAP/TAZ activation induced by Lats1/2 deletion cooperates with Kras to promote UM progression via downstream transcriptional reinforcement. Furthermore, dual inhibition of YAP/TAZ and Ras/MAPK synergizes to suppress oncogenic growth of human UM cells. Our data highlight the functional significance of Lats-YAP/TAZ in UM initiation and progression in vivo and suggest combination inhibition of YAP/TAZ and Ras/MAPK as a new therapeutic strategy for UM.

Genome-Wide Transcriptional Regulation of the Long Non-coding RNA Steroid Receptor RNA Activator in Human Erythroblasts

Erythropoiesis of human hematopoietic stem cells (HSCs) maintains generation of red blood cells throughout life. However, little is known how human erythropoiesis is regulated by long non-coding RNAs (lncRNAs). By using ChIRP-seq, we report here that the lncRNA steroid receptor RNA activator (SRA) occupies chromatin, and co-localizes with CTCF, H3K4me3, and H3K27me3 genome-wide in human erythroblast cell line K562. CTCF binding sites that are also occupied by SRA are enriched for either H3K4me3 or H3K27me3. Transcriptome-wide analyses reveal that SRA facilitates expression of erythroid-associated genes, while repressing leukocyte-associated genes in both K562 and CD36-positive primary human proerythroblasts derived from HSCs. We find that SRA-regulated genes are enriched by both CTCF and SRA bindings. Further, silencing of SRA decreases expression of the erythroid-specific markers TFRC and GYPA, and down-regulates expression of globin genes in both K562 and human proerythroblast cells. Taken together, our findings establish that the lncRNA SRA occupies chromatin, and promotes transcription of erythroid genes, therefore facilitating human erythroid transcriptional program.

Distinct and temporary-restricted epigenetic mechanisms regulate human αβ and γδ T cell development

The development of TCRαβ and TCRγδ T cells comprises a step-wise process in which regulatory events control differentiation and lineage outcome. To clarify these mechanisms, we employed RNA-sequencing, ATAC-sequencing and ChIPmentation on well-defined thymocyte subsets that represent the continuum of human T cell development. The chromatin accessibility dynamics show clear stage specificity and reveal that human T cell-lineage commitment is marked by GATA3- and BCL11B-dependent closing of PU.1 sites. A temporary increase in H3K27me3 without open chromatin modifications is unique for β-selection, whereas emerging γδ T cells, which originate from common precursors of β-selected cells, show large chromatin accessibility changes due to strong T cell receptor (TCR) signaling. Furthermore, we unravel distinct chromatin landscapes between CD4⁺ and CD8⁺ αβ-lineage cells that support their effector functions and reveal gene-specific mechanisms that define mature T cells. This resource provides a framework for studying gene regulatory mechanisms that drive normal and malignant human T cell development.

Genome-wide cooperation of EMT transcription factor ZEB1 with YAP and AP-1 in breast cancer

Invasion, metastasis and therapy resistance are the major cause of cancer‐associated deaths, and the EMT‐inducing transcription factor ZEB1 is a crucial stimulator of these processes. While work on ZEB1 has mainly focused on its role as a transcriptional repressor, it can also act as a transcriptional activator. To further understand these two modes of action, we performed a genome‐wide ZEB1 binding study in triple‐negative breast cancer cells. We identified ZEB1 as a novel interactor of the AP‐1 factors FOSL1 and JUN and show that, together with the Hippo pathway effector YAP, they form a transactivation complex, predominantly activating tumour‐promoting genes, thereby synergising with its function as a repressor of epithelial genes. High expression of ZEB1, YAP, FOSL1 and JUN marks the aggressive claudin‐low subtype of breast cancer, indicating the translational relevance of our findings. Thus, our results link critical tumour‐promoting transcription factors: ZEB1, AP‐1 and Hippo pathway factors. Disturbing their molecular interaction may provide a promising treatment option for aggressive cancer types.

Transcriptional Patterning of the Ventricular Cardiac Conduction System

RATIONALE

The heartbeat is organized by the cardiac conduction system (CCS), a specialized network of cardiomyocytes. Patterning of the CCS into atrial node versus ventricular conduction system (VCS) components with distinct physiology is essential for the normal heartbeat. Distinct node versus VCS physiology has been recognized for more than a century, but the molecular basis of this regional patterning is not well understood.

OBJECTIVE

To study the genetic and genomic mechanisms underlying node versus VCS distinction and investigate rhythm consequences of failed VCS patterning.

METHODS AND RESULTS

Using mouse genetics, we found that the balance between T-box transcriptional activator, Tbx5, and T-box transcriptional repressor, Tbx3, determined the molecular and functional output of VCS myocytes. Adult VCS-specific removal of Tbx5 or overexpression of Tbx3 re-patterned the fast VCS into slow, nodal-like cells based on molecular and functional criteria. In these cases, gene expression profiling showed diminished expression of genes required for VCS-specific fast conduction but maintenance of expression of genes required for nodal slow conduction physiology. Action potentials of Tbx5-deficient VCS myocytes adopted nodal-specific characteristics, including increased action potential duration and cellular automaticity. Removal of Tbx5 in vivo precipitated inappropriate depolarizations in the atrioventricular (His)-bundle associated with lethal ventricular arrhythmias. TBX5 bound and directly activated cis-regulatory elements at fast conduction channel genes required for fast physiological characteristics of the VCS action potential, defining the identity of the adult VCS.

CONCLUSIONS

The CCS is patterned entirely as a slow, nodal ground state, with a T-box dependent, physiologically dominant, fast conduction network driven specifically in the VCS. Disruption of the fast VCS gene regulatory network allowed nodal physiology to emerge, providing a plausible molecular mechanism for some lethal ventricular arrhythmias.

Distinct features of nucleolus-associated domains in mouse embryonic stem cells

Heterochromatin in eukaryotic interphase cells frequently localizes to the nucleolar periphery (nucleolus-associated domains (NADs)) and the nuclear lamina (lamina-associated domains (LADs)). Gene expression in somatic cell NADs is generally low, but NADs have not been characterized in mammalian stem cells. Here, we generated the first genome-wide map of NADs in mouse embryonic stem cells (mESCs) via deep sequencing of chromatin associated with biochemically purified nucleoli. As we had observed in mouse embryonic fibroblasts (MEFs), the large type I subset of NADs overlaps with constitutive LADs and is enriched for features of constitutive heterochromatin, including late replication timing and low gene density and expression levels. Conversely, the type II NAD subset overlaps with loci that are not lamina-associated, but in mESCs, type II NADs are much less abundant than in MEFs. mESC NADs are also much less enriched in H3K27me3 modified regions than are NADs in MEFs. Additionally, comparision of MEF and mESC NADs revealed enrichment of developmentally regulated genes in cell-type-specific NADs. Together, these data indicate that NADs are a developmentally dynamic component of heterochromatin. These studies implicate association with the nucleolar periphery as a mechanism for developmentally regulated gene expression and will facilitate future studies of NADs during mESC differentiation.

Consequences of HSF knockdown on gene expression during the heat shock response in Tigriopus californicus

Although the existence of a cellular heat shock response is nearly universal, its relationship to organismal thermal tolerance is not completely understood. Many of the genes involved are known to be regulated by the highly conserved heat shock transcription factor-1 (HSF-1), yet the regulatory network is not fully characterized. Here, we investigated the role of HSF-1 in gene expression following thermal stress using knockdown of HSF-1 by RNA interference in the intertidal copepod Tigriopus californicus We observed some evidence for decreased transcription of heat shock protein genes following knockdown, supporting the widely acknowledged role of HSF-1 in the heat shock response. However, the majority of differentially expressed genes between the control and HSF-1 knockdown groups were upregulated, suggesting that HSF-1 normally functions to repress their expression. Differential expression observed in genes related to chitin and cuticle formation lends support to previous findings that these processes are highly regulated following heat stress. We performed a genome scan and identified a set of 396 genes associated with canonical heat shock elements. RNA-seq data did not find those genes to be more highly represented in our HSF-1 knockdown treatment, indicating that requirements for binding and interaction of HSF-1 with a given gene are not simply predicted by the presence of HSF-1 binding sites. Further study of the pathways implicated by these results and future comparisons among populations of T. californicus may help us understand the role and importance of HSF-1 in the heat shock response and, more broadly, in organismal thermal tolerance.

Optimization of Ligands Using Focused DNA-Encoded Libraries To Develop a Selective, Cell-Permeable CBX8 Chromodomain Inhibitor

Polycomb repressive complex 1 (PRC1) is critical for mediating gene expression during development. Five chromobox (CBX) homolog proteins, CBX2, CBX4, CBX6, CBX7, and CBX8, are incorporated into PRC1 complexes, where they mediate targeting to trimethylated lysine 27 of histone H3 (H3K27me3) via the N-terminal chromodomain (ChD). Individual CBX paralogs have been implicated as drug targets in cancer; however, high similarities in sequence and structure among the CBX ChDs provide a major obstacle in developing selective CBX ChD inhibitors. Here we report the selection of small, focused, DNA-encoded libraries (DELs) against multiple homologous ChDs to identify modifications to a parental ligand that confer both selectivity and potency for the ChD of CBX8. This on-DNA, medicinal chemistry approach enabled the development of SW2_110A, a selective, cell-permeable inhibitor of the CBX8 ChD. SW2_110A binds CBX8 ChD with a K_d of 800 nM, with minimal 5-fold selectivity for CBX8 ChD over all other CBX paralogs in vitro. SW2_110A specifically inhibits the association of CBX8 with chromatin in cells and inhibits the proliferation of THP1 leukemia cells driven by the MLL-AF9 translocation. In THP1 cells, SW2_110A treatment results in a significant decrease in the expression of MLL-AF9 target genes, including HOXA9, validating the previously established role for CBX8 in MLL-AF9 transcriptional activation, and defining the ChD as necessary for this function. The success of SW2_110A provides great promise for the development of highly selective and cell-permeable probes for the full CBX family. In addition, the approach taken provides a proof-of-principle demonstration of how DELs can be used iteratively for optimization of both ligand potency and selectivity.

Engagement of DNA and H3K27me3 by the CBX8 chromodomain drives chromatin association

Polycomb repressive complex 1 (PRC1) is critical for mediating gene repression during development and adult stem cell maintenance. Five CBX proteins, CBX2,4,6,7,8, form mutually exclusive PRC1 complexes and are thought to play a role in the association of PRC1 with chromatin. Specifically, the N-terminal chromodomain (CD) in the CBX proteins is thought to mediate specific targeting to methylated histones. For CBX8, however, the chromodomain has demonstrated weak affinity and specificity for methylated histones in vitro, leaving doubt as to its role in CBX8 chromatin association. Here, we investigate the function of the CBX8 CD in vitro and in vivo. We find that the CD is in fact a major driver of CBX8 chromatin association and determine that this is driven by both histone and previously unrecognized DNA binding activity. We characterize the structural basis of histone and DNA binding and determine how they integrate on multiple levels. Notably, we find that the chromatin environment is critical in determining the ultimate function of the CD in CBX8 association.

Drug-induced chromatin accessibility changes associate with sensitivity to liver tumor promotion

Liver cancer susceptibility varies amongst humans and between experimental animal models because of multiple genetic and epigenetic factors. The molecular characterization of such susceptibilities has the potential to enhance cancer risk assessment of xenobiotic exposures and disease prevention strategies. Here, using DNase I hypersensitivity mapping coupled with transcriptomic profiling, we investigate perturbations in cis-acting gene regulatory elements associated with the early stages of phenobarbital (PB)-mediated liver tumor promotion in susceptible versus resistant mouse strains (B6C3F1 versus C57BL/6J). Integrated computational analyses of strain-selective changes in liver chromatin accessibility underlying PB response reveal differential epigenetic regulation of molecular pathways associated with PB-mediated tumor promotion, including Wnt/β-catenin signaling. Complementary transcription factor motif analyses reveal mouse strain-selective gene regulatory networks and a novel role for Stat, Smad, and Fox transcription factors in the early stages of PB-mediated tumor promotion. Mapping perturbations in cis-acting gene regulatory elements provides novel insights into the molecular basis for susceptibility to xenobiotic-induced rodent liver tumor promotion and has the potential to enhance mechanism-based cancer risk assessments of xenobiotic exposures.

Two contrasting classes of nucleolus-associated domains in mouse fibroblast heterochromatin

In interphase eukaryotic cells, almost all heterochromatin is located adjacent to the nucleolus or to the nuclear lamina, thus defining nucleolus-associated domains (NADs) and lamina-associated domains (LADs), respectively. Here, we determined the first genome-scale map of murine NADs in mouse embryonic fibroblasts (MEFs) via deep sequencing of chromatin associated with purified nucleoli. We developed a Bioconductor package called NADfinder and demonstrated that it identifies NADs more accurately than other peak-calling tools, owing to its critical feature of chromosome-level local baseline correction. We detected two distinct classes of NADs. Type I NADs associate frequently with both the nucleolar periphery and the nuclear lamina, and generally display characteristics of constitutive heterochromatin, including late DNA replication, enrichment of H3K9me3, and little gene expression. In contrast, Type II NADs associate with nucleoli but do not overlap with LADs. Type II NADs tend to replicate earlier, display greater gene expression, and are more often enriched in H3K27me3 than Type I NADs. The nucleolar associations of both classes of NADs were confirmed via DNA-FISH, which also detected Type I but not Type II probes enriched at the nuclear lamina. Type II NADs are enriched in distinct gene classes, including factors important for differentiation and development. In keeping with this, we observed that a Type II NAD is developmentally regulated, and present in MEFs but not in undifferentiated embryonic stem (ES) cells.

Rlim-Dependent and -Independent Pathways for X Chromosome Inactivation in Female ESCs

During female mouse embryogenesis, two forms of X chromosome inactivation (XCI) ensure dosage compensation from sex chromosomes. Beginning at the four-cell stage, imprinted XCI (iXCI) exclusively silences the paternal X (Xp), and this pattern is maintained in extraembryonic cell types. Epiblast cells, which give rise to the embryo proper, reactivate the Xp (XCR) and undergo a random form of XCI (rXCI) around implantation. Both iXCI and rXCI depend on the long non-coding RNA Xist. The ubiquitin ligase RLIM is required for iXCI in vivo and occupies a central role in current models of rXCI. Here, we demonstrate the existence of Rlim-dependent and Rlim-independent pathways for rXCI in differentiating female ESCs. Upon uncoupling these pathways, we find more efficient Rlim-independent XCI in ESCs cultured under physiological oxygen conditions. Our results revise current models of rXCI and suggest that caution must be taken when comparing XCI studies in ESCs and mice.

The EZH2 SANT1 domain is a histone reader providing sensitivity to the modification state of the H4 tail

SANT domains are found in a number of chromatin regulators. They contain approximately 50 amino acids and have high similarity to the DNA binding domain of Myb related proteins. Though some SANT domains associate with DNA others have been found to bind unmodified histone tails. There are two SANT domains in Enhancer of Zeste 2 (EZH2), the catalytic subunit of the Polycomb Repressive Complex 2 (PRC2), of unknown function. Here we show that the first SANT domain (SANT1) of EZH2 is a histone binding domain with specificity for the histone H4 N-terminal tail. Using NMR spectroscopy, mutagenesis, and molecular modeling we structurally characterize the SANT1 domain and determine the molecular mechanism of binding to the H4 tail. Though not important for histone binding, we find that the adjacent stimulation response motif (SRM) stabilizes SANT1 and transiently samples its active form in solution. Acetylation of H4K16 (H4K16ac) or acetylation or methylation of H4K20 (H4K20ac and H4K20me3) are seen to abrogate binding of SANT1 to H4, which is consistent with these modifications being anti-correlated with H3K27me3 in-vivo. Our results provide significant insight into this important regulatory region of EZH2 and the first characterization of the molecular mechanism of SANT domain histone binding.

GWAS and colocalization analyses implicate carotid intima-media thickness and carotid plaque loci in cardiovascular outcomes

Carotid artery intima media thickness (cIMT) and carotid plaque are measures of subclinical atherosclerosis associated with ischemic stroke and coronary heart disease (CHD). Here, we undertake meta-analyses of genome-wide association studies (GWAS) in 71,128 individuals for cIMT, and 48,434 individuals for carotid plaque traits. We identify eight novel susceptibility loci for cIMT, one independent association at the previously-identified PINX1 locus, and one novel locus for carotid plaque. Colocalization analysis with nearby vascular expression quantitative loci (cis-eQTLs) derived from arterial wall and metabolic tissues obtained from patients with CHD identifies candidate genes at two potentially additional loci, ADAMTS9 and LOXL4. LD score regression reveals significant genetic correlations between cIMT and plaque traits, and both cIMT and plaque with CHD, any stroke subtype and ischemic stroke. Our study provides insights into genes and tissue-specific regulatory mechanisms linking atherosclerosis both to its functional genomic origins and its clinical consequences in humans.

Interactome mapping defines BRG1, a component of the SWI/SNF chromatin remodeling complex, as a new partner of the transcriptional regulator CTCF

The highly conserved zinc finger CCCTC-binding factor (CTCF) regulates genomic imprinting and gene expression by acting as a transcriptional activator or repressor of promoters and insulator of enhancers. The multiple functions of CTCF are accomplished by co-association with other protein partners and are dependent on genomic context and tissue specificity. Despite the critical role of CTCF in the organization of genome structure, to date, only a subset of CTCF interaction partners have been identified. Here we present a large-scale identification of CTCF-binding partners using affinity purification and high-resolution LC-MS/MS analysis. In addition to functional enrichment of specific protein families such as the ribosomal proteins and the DEAD box helicases, we identified novel high-confidence CTCF interactors that provide a still unexplored biochemical context for CTCF's multiple functions. One of the newly validated CTCF interactors is BRG1, the major ATPase subunit of the chromatin remodeling complex SWI/SNF, establishing a relationship between two master regulators of genome organization. This work significantly expands the current knowledge of the human CTCF interactome and represents an important resource to direct future studies aimed at uncovering molecular mechanisms modulating CTCF pleiotropic functions throughout the genome.

CBFβ-SMMHC Inhibition Triggers Apoptosis by Disrupting MYC Chromatin Dynamics in Acute Myeloid Leukemia

The fusion oncoprotein CBFβ-SMMHC, expressed in leukemia cases with chromosome 16 inversion, drives leukemia development and maintenance by altering the activity of the transcription factor RUNX1. Here, we demonstrate that CBFβ-SMMHC maintains cell viability by neutralizing RUNX1-mediated repression of MYC expression. Upon pharmacologic inhibition of the CBFβ-SMMHC/RUNX1 interaction, RUNX1 shows increased binding at three MYC distal enhancers, where it represses MYC expression by mediating the replacement of the SWI/SNF complex component BRG1 with the polycomb-repressive complex component RING1B, leading to apoptosis. Combining the CBFβ-SMMHC inhibitor with the BET inhibitor JQ1 eliminates inv(16) leukemia in human cells and a mouse model. Enhancer-interaction analysis indicated that the three enhancers are physically connected with the MYC promoter, and genome-editing analysis demonstrated that they are functionally implicated in deregulation of MYC expression. This study reveals a mechanism whereby CBFβ-SMMHC drives leukemia maintenance and suggests that inhibitors targeting chromatin activity may prove effective in inv(16) leukemia therapy.