Although the core subunits of Polycomb group (PcG) complexes are well characterized, little is known about the dynamics of these protein complexes during cellular differentiation. We used quantitative interaction proteomics and genome-wide profiling to study PcG proteins in mouse embryonic stem cells (ESCs) and neural progenitor cells (NPCs). We found that the stoichiometry and genome-wide binding of PRC1 and PRC2 were highly dynamic during neural differentiation. Intriguingly, we observed a downregulation and loss of PRC2 from chromatin marked with trimethylated histone H3 K27 (H3K27me3) during differentiation, whereas PRC1 was retained at these sites. Additionally, we found PRC1 at enhancer and promoter regions independently of PRC2 binding and H3K27me3. Finally, overexpression of NPC-specific PRC1 interactors in ESCs led to increased Ring1b binding to, and decreased expression of, NPC-enriched Ring1b-target genes. In summary, our integrative analyses uncovered dynamic PcG subcomplexes and their widespread colocalization with active chromatin marks during differentiation.

Here, we report DNA methylation and hydroxymethylation dynamics at nucleotide resolution using C/EBPα-enhanced reprogramming of B cells into induced pluripotent cells (iPSCs). We observed successive waves of hydroxymethylation at enhancers, concomitant with a decrease in DNA methylation, suggesting active demethylation. Consistent with this finding, ablation of the DNA demethylase Tet2 almost completely abolishes reprogramming. C/EBPα, Klf4, and Tfcp2l1 each interact with Tet2 and recruit the enzyme to specific DNA sites. During reprogramming, some of these sites maintain high levels of 5hmC, and enhancers and promoters of key pluripotency factors become demethylated as early as 1 day after Yamanaka factor induction. Surprisingly, methylation changes precede chromatin opening in distinct chromatin regions, including Klf4 bound sites, revealing a pioneer factor activity associated with alternation in DNA methylation. Rapid changes in hydroxymethylation similar to those in B cells were also observed during compound-accelerated reprogramming of fibroblasts into iPSCs, highlighting the generality of our observations.

The molecular mechanisms underlying specification from embryonic stem cells (ESCs) and maintenance of neural progenitor cells (NPCs) are largely unknown. Recently, we reported that the Zuotin-related factor 1 (Zrf1) is necessary for chromatin displacement of the Polycomb-repressive complex 1 (PRC1). We found that Zrf1 is required for NPC specification from ESCs and that it promotes the expression of NPC markers, including the key regulator Pax6. Moreover, Zrf1 is essential to establish and maintain Wnt ligand expression levels, which are necessary for NPC self-renewal. Reactivation of proper Wnt signaling in Zrf1-depleted NPCs restores Pax6 expression and the self-renewal capacity. ESC-derived NPCs in vitro resemble most of the characteristics of the self-renewing NPCs located in the developing embryonic cortex, which are termed radial glial cells (RGCs). Depletion of Zrf1 in vivo impairs the expression of key self-renewal regulators and Wnt ligand genes in RGCs. Thus, we demonstrate that Zrf1 plays an essential role in NPC generation and maintenance.

Experimental and bioinformatic studies of transcription initiation by RNA polymerase II (RNAP2) have revealed a mechanism of RNAP2 transcription initiation less uniform across gene promoters than initially thought. However, the general transcription factor TFIIB is presumed to be universally required for RNAP2 transcription initiation. Based on bioinformatic analysis of data and effects of TFIIB knockdown in primary and transformed cell lines on cellular functionality and global gene expression, we report that TFIIB is dispensable for transcription of many human promoters, but is essential for herpes simplex virus-1 (HSV-1) gene transcription and replication. We report a novel cell cycle TFIIB regulation and localization of the acetylated TFIIB variant on the transcriptionally silent mitotic chromatids. Taken together, these results establish a new paradigm for TFIIB functionality in human gene expression, which when downregulated has potent anti-viral effects.

Colorectal carcinoma represents a heterogeneous entity, with only a fraction of the tumours responding to available therapies, requiring a better molecular understanding of the disease in precision oncology. To address this challenge, the OncoTrack consortium recruited 106 CRC patients (stages I-IV) and developed a pre-clinical platform generating a compendium of drug sensitivity data totalling >4,000 assays testing 16 clinical drugs on patient-derived in vivo and in vitro models. This large biobank of 106 tumours, 35 organoids and 59 xenografts, with extensive omics data comparing donor tumours and derived models provides a resource for advancing our understanding of CRC. Models recapitulate many of the genetic and transcriptomic features of the donors, but defined less complex molecular sub-groups because of the loss of human stroma. Linking molecular profiles with drug sensitivity patterns identifies novel biomarkers, including a signature outperforming RAS/RAF mutations in predicting sensitivity to the EGFR inhibitor cetuximab.

Reprimo (RPRM), a member of the RPRM gene family, is a tumor-suppressor gene involved in the regulation of the p53-mediated cell cycle arrest at G2/M. RPRM has been associated with malignant tumor progression and proposed as a potential biomarker for early cancer detection. However, the expression and role of RPRM, as well as its family, are poorly understood and their physiology is as yet unstudied. In this scenario, a model system like the zebrafish could serve to dissect the role of the RPRM family members in vivo. Phylogenetic analysis reveals that RPRM and RPRML have been differentially retained by most species throughout vertebrate evolution, yet RPRM3 has been retained only in a small group of distantly related species, including zebrafish. Herein, we characterized the spatiotemporal expression of RPRM (present in zebrafish as an infraclass duplication rprma/rprmb), RPRML and RPRM3 in the zebrafish. By whole-mount in situ hybridization (WISH) and fluorescent in situ hybridization (FISH), we demonstrate that rprm (rprma/rprmb) and rprml show a similar spatiotemporal expression profile during zebrafish development. At early developmental stages rprmb is expressed in somites. After one day post-fertilization, rprm (rprma/rprmb) and rprml are expressed in the notochord, brain, blood vessels and digestive tube. On the other hand, rprm3 shows the most unique expression profile, being expressed only in the central nervous system (CNS). We assessed the expression patterns of RPRM gene transcripts in adult zebrafish and human RPRM protein product in tissue samples by RT-qPCR and immunohistochemistry (IHC) staining, respectively. Strikingly, tissue-specific expression patterns of the RPRM transcripts and protein are conserved between zebrafish and humans. We propose the zebrafish as a powerful tool to elucidate the both physiological and pathological roles of the RPRM gene family.

The prediction and annotation of the genomic regions involved in gene expression has been largely explored. Most of the energy has been devoted to the development of approaches that detect transcription start sites, leaving the identification of regulatory regions and their functional transcription factor binding sites (TFBSs) largely unexplored and with important quantitative and qualitative methodological gaps.

The capacity of tumour cells to maintain continual overgrowth potential has been linked to the commandeering of normal self-renewal pathways. Using an epithelial cancer model in Drosophila melanogaster, we carried out an overexpression screen for oncogenes capable of cooperating with the loss of the epithelial apico-basal cell polarity regulator, scribbled (scrib), and identified the cell fate regulator, Abrupt, a BTB-zinc finger protein. Abrupt overexpression alone is insufficient to transform cells, but in cooperation with scrib loss of function, Abrupt promotes the formation of massive tumours in the eye/antennal disc. The steroid hormone receptor coactivator, Taiman (a homologue of SRC3/AIB1), is known to associate with Abrupt, and Taiman overexpression also drives tumour formation in cooperation with the loss of Scrib. Expression arrays and ChIP-Seq indicates that Abrupt overexpression represses a large number of genes, including steroid hormone-response genes and multiple cell fate regulators, thereby maintaining cells within an epithelial progenitor-like state. The progenitor-like state is characterised by the failure to express the conserved Eyes absent/Dachshund regulatory complex in the eye disc, and in the antennal disc by the failure to express cell fate regulators that define the temporal elaboration of the appendage along the proximo-distal axis downstream of Distalless. Loss of scrib promotes cooperation with Abrupt through impaired Hippo signalling, which is required and sufficient for cooperative overgrowth with Abrupt, and JNK (Jun kinase) signalling, which is required for tumour cell migration/invasion but not overgrowth. These results thus identify a novel cooperating oncogene, identify mammalian family members of which are also known oncogenes, and demonstrate that epithelial tumours in Drosophila can be characterised by the maintenance of a progenitor-like state.

We have isolated a microsporidium from a laboratory stock of the European grape vine moth, Lobesia botrana Den. et Schiff. (Lepidoptera, Tortricidae). Screening of this stock showed an infection rate of more than 90%, whereas field collected larvae from three different locations in Rhineland-Palatinate (Germany) did not demonstrate any signs of infection. Light and electron microscopic investigations of infected insects showed that gross pathology, morphology, and ultrastructure of the microsporidium are similar to those described earlier for Pleistophora legeri. Comparative phylogenetic analysis of the small subunit rDNA using maximum likelihood, maximum parsimony, and neighbour joining distance methods showed that our isolate was closely related to Cystosporogenes operophterae. Based on our morphological and molecular investigations we propose to rename this species Cystosporogenes legeri nov. comb.

The trithorax group (trxG) and Polycomb group (PcG) proteins are responsible for the maintenance of stable transcriptional patterns of many developmental regulators. They bind to specific regions of DNA and direct the post-translational modifications of histones, playing a role in the dynamics of chromatin structure.

Exon skipping is an effective strategy for the treatment of many Duchenne Muscular Dystrophy (DMD) mutations. Natural exon skipping observed in several DMD cases can help in identifying novel therapeutic tools. Here, we show a DMD study case where the lack of a splicing factor (Celf2a), which results in exon skipping and dystrophin rescue, is due to a maternally inherited trans-generational epigenetic silencing. We found that the study case and his mother express a repressive long non-coding RNA, DUXAP8, whose presence correlates with silencing of the Celf2a coding region. We also demonstrate that DUXAP8 expression is lost upon cell reprogramming and that, upon induction of iPSCs into myoblasts, Celf2a expression is recovered leading to the loss of exon skipping and loss of dystrophin synthesis. Finally, CRISPR/Cas9 inactivation of the splicing factor Celf2a was proven to ameliorate the pathological state in other DMD backgrounds establishing Celf2a ablation or inactivation as a novel therapeutic approach for the treatment of Duchenne Muscular Dystrophy.

Skin aging is characterized by structural and functional changes that contribute to age-associated frailty. This probably depends on synergy between alterations in the local niche and stem cell-intrinsic changes, underscored by proinflammatory microenvironments that drive pleotropic changes. The nature of these age-associated inflammatory cues, or how they affect tissue aging, is unknown. Based on single-cell RNA sequencing of the dermal compartment of mouse skin, we show a skew towards an IL-17-expressing phenotype of T helper cells, γδ T cells and innate lymphoid cells in aged skin. Importantly, in vivo blockade of IL-17 signaling during aging reduces the proinflammatory state of the skin, delaying the appearance of age-related traits. Mechanistically, aberrant IL-17 signals through NF-κB in epidermal cells to impair homeostatic functions while promoting an inflammatory state. Our results indicate that aged skin shows signs of chronic inflammation and that increased IL-17 signaling could be targeted to prevent age-associated skin ailments.

Cell cycle progression is linked to transcriptome dynamics and variations in the response of pluripotent cells to differentiation cues, mostly through unknown determinants. Here, we characterized the cell-cycle-associated transcriptome and proteome of mouse embryonic stem cells (mESCs) in naive ground state. We found that the thymine DNA glycosylase (TDG) is a cell-cycle-regulated co-factor of the tumor suppressor p53. Furthermore, TDG and p53 co-bind ESC-specific cis-regulatory elements and thereby control transcription of p53-dependent genes during self-renewal. We determined that the dynamic expression of TDG is required to promote the cell-cycle-associated transcriptional heterogeneity. Moreover, we demonstrated that transient depletion of TDG influences cell fate decisions during the early differentiation of mESCs. Our findings reveal an unanticipated role of TDG in promoting molecular heterogeneity during the cell cycle and highlight the central role of protein dynamics for the temporal control of cell fate during development.

Wound healing is an essential homeostatic mechanism that maintains the epithelial barrier integrity after tissue damage. Although we know the overall steps in wound healing, many of the underlying molecular mechanisms remain unclear. Genetically amenable systems, such as wound healing in Drosophila imaginal discs, do not model all aspects of the repair process. However, they do allow the less understood aspects of the healing response to be explored, e.g., which signal(s) are responsible for initiating tissue remodeling? How is sealing of the epithelia achieved? Or, what inhibitory cues cancel the healing machinery upon completion? Answering these and other questions first requires the identification and functional analysis of wound specific genes. A variety of different microarray analyses of murine and humans have identified characteristic profiles of gene expression at the wound site, however, very few functional studies in healing regulation have been carried out. We developed an experimentally controlled method that is healing-permissive and that allows live imaging and biochemical analysis of cultured imaginal discs. We performed comparative genome-wide profiling between Drosophila imaginal cells actively involved in healing versus their non-engaged siblings. Sets of potential wound-specific genes were subsequently identified. Importantly, besides identifying and categorizing new genes, we functionally tested many of their gene products by genetic interference and overexpression in healing assays. This non-saturated analysis defines a relevant set of genes whose changes in expression level are functionally significant for proper tissue repair. Amongst these we identified the TCP1 chaperonin complex as a key regulator of the actin cytoskeleton essential for the wound healing response. There is promise that our newly identified wound-healing genes will guide future work in the more complex mammalian wound healing response.

Cellular reprogramming to iPSCs has uncovered unsuspected links between tumor suppressors and pluripotency factors. Using this system, it was possible to identify tumor suppressor p27 as a repressor of Sox2 during differentiation. This led to the demonstration that defects in the repression of Sox2 can contribute to tumor development. The members of the retinoblastoma family of pocket proteins, pRb, p107 and p130, are negative regulators of the cell cycle with tumor suppressor activity and with roles in differentiation. In this work we studied the relative contribution of the retinoblastoma family members to the regulation of Sox2 expression. We found that deletion of Rb or p130 leads to impaired repression of Sox2, a deffect amplified by inactivation of p53. We also identified binding of pRb and p130 to an enhancer with crucial regulatory activity on Sox2 expression. Using cellular reprogramming we tested the impact of the defective repression of Sox2 and confirmed that Rb deficiency allows the generation of iPSCs in the absence of exogenous Sox2. Finally, partial depletion of Sox2 positive cells reduced the pituitary tumor development initiated by Rb loss in vivo. In summary, our results show that Sox2 repression by pRb is a relevant mechanism of tumor suppression.

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

We address the challenge of regulatory sequence alignment with a new method, Pro-Coffee, a multiple aligner specifically designed for homologous promoter regions. Pro-Coffee uses a dinucleotide substitution matrix estimated on alignments of functional binding sites from TRANSFAC. We designed a validation framework using several thousand families of orthologous promoters. This dataset was used to evaluate the accuracy for predicting true human orthologs among their paralogs. We found that whereas other methods achieve on average 73.5% accuracy, and 77.6% when trained on that same dataset, the figure goes up to 80.4% for Pro-Coffee. We then applied a novel validation procedure based on multi-species ChIP-seq data. Trained and untrained methods were tested for their capacity to correctly align experimentally detected binding sites. Whereas the average number of correctly aligned sites for two transcription factors is 284 for default methods and 316 for trained methods, Pro-Coffee achieves 331, 16.5% above the default average. We find a high correlation between a method's performance when classifying orthologs and its ability to correctly align proven binding sites. Not only has this interesting biological consequences, it also allows us to conclude that any method that is trained on the ortholog data set will result in functionally more informative alignments.

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels--to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human-mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

The common cutworm (Agrotis segetum) and the black cutworm (Agrotis ipsilon) are serious soil pests of many vegetable and field crops all over the world. We have demonstrated the cross-infectivity of two baculoviruses, A. segetum nucleopolyhedrovirus (AgseNPV) and A. ipsilon nucleopolyhedrovirus (AgipNPV) for these two insect pests. The susceptibility of A. segetum to AgipNPV was confirmed by DNA restriction endonuclease analyses of DNA isolated from virus harvested from infected A. segetum larvae. For an initial comparison of both viruses, partial polyhedrin sequences were amplified by PCR, cloned, and sequenced. Both viruses shared a very similar polyhedrin gene sequence resulting in only three amino acid substitutions. Phylogenetic analyses clearly demonstrated that both viruses belong to NPV group II and are most closely related to a clade consisting of Spodoptera exigua NPV, Spodoptera frugiperda NPV, and Spodoptera littoralis NPV. Since AgipNPV shows high virulence for both cutworm species, it appears to be a suitable candidate as a single biological control agent of A. segetum and A. ipsilon.

Changes in the epigenetic regulation of gene expression have a central role in evolution. Here, we extensively profiled a panel of human, chimpanzee, gorilla, orangutan, and macaque lymphoblastoid cell lines (LCLs), using ChIP-seq for five histone marks, ATAC-seq and RNA-seq, further complemented with whole genome sequencing (WGS) and whole genome bisulfite sequencing (WGBS). We annotated regulatory elements (RE) and integrated chromatin contact maps to define gene regulatory architectures, creating the largest catalog of RE in primates to date. We report that epigenetic conservation and its correlation with sequence conservation in primates depends on the activity state of the regulatory element. Our gene regulatory architectures reveal the coordination of different types of components and highlight the role of promoters and intragenic enhancers (gE) in the regulation of gene expression. We observe that most regulatory changes occur in weakly active gE. Remarkably, novel human-specific gE with weak activities are enriched in human-specific nucleotide changes. These elements appear in genes with signals of positive selection and human acceleration, tissue-specific expression, and particular functional enrichments, suggesting that the regulatory evolution of these genes may have contributed to human adaptation.