Insulin/IGF-1-mediated longevity is marked by reduced protein metabolism.

Mutations in the daf-2 gene of the conserved Insulin/Insulin-like Growth Factor (IGF-1) pathway double the lifespan of the nematode Caenorhabditis elegans. This phenotype is completely suppressed by deletion of Forkhead transcription factor daf-16. To uncover regulatory mechanisms coordinating this extension of life, we employed a quantitative proteomics strategy with daf-2 mutants in comparison with N2 and daf-16; daf-2 double mutants. This revealed a remarkable longevity-specific decrease in proteins involved in mRNA processing and transport, the translational machinery, and protein metabolism. Correspondingly, the daf-2 mutants display lower amounts of mRNA and 20S proteasome activity, despite maintaining total protein levels equal to that observed in wild types. Polyribosome profiling in the daf-2 and daf-16;daf-2 double mutants confirmed a daf-16-dependent reduction in overall translation, a phenotype reminiscent of Dietary Restriction-mediated longevity, which was independent of germline activity. RNA interference (RNAi)-mediated knockdown of proteins identified by our approach resulted in modified C. elegans lifespan confirming the importance of these processes in Insulin/IGF-1-mediated longevity. Together, the results demonstrate a role for the metabolism of proteins in the Insulin/IGF-1-mediated extension of life.

Posterior association networks and functional modules inferred from rich phenotypes of gene perturbations.

Combinatorial gene perturbations provide rich information for a systematic exploration of genetic interactions. Despite successful applications to bacteria and yeast, the scalability of this approach remains a major challenge for higher organisms such as humans. Here, we report a novel experimental and computational framework to efficiently address this challenge by limiting the 'search space' for important genetic interactions. We propose to integrate rich phenotypes of multiple single gene perturbations to robustly predict functional modules, which can subsequently be subjected to further experimental investigations such as combinatorial gene silencing. We present posterior association networks (PANs) to predict functional interactions between genes estimated using a Bayesian mixture modelling approach. The major advantage of this approach over conventional hypothesis tests is that prior knowledge can be incorporated to enhance predictive power. We demonstrate in a simulation study and on biological data, that integrating complementary information greatly improves prediction accuracy. To search for significant modules, we perform hierarchical clustering with multiscale bootstrap resampling. We demonstrate the power of the proposed methodologies in applications to Ewing's sarcoma and human adult stem cells using publicly available and custom generated data, respectively. In the former application, we identify a gene module including many confirmed and highly promising therapeutic targets. Genes in the module are also significantly overrepresented in signalling pathways that are known to be critical for proliferation of Ewing's sarcoma cells. In the latter application, we predict a functional network of chromatin factors controlling epidermal stem cell fate. Further examinations using ChIP-seq, ChIP-qPCR and RT-qPCR reveal that the basis of their genetic interactions may arise from transcriptional cross regulation. A Bioconductor package implementing PAN is freely available online at http://bioconductor.org/packages/release/bioc/html/PANR.html.

Trichothiodystrophy-associated MPLKIP maintains DBR1 levels for proper lariat debranching and ectodermal differentiation.

The brittle hair syndrome Trichothiodystrophy (TTD) is characterized by variable clinical features, including photosensitivity, ichthyosis, growth retardation, microcephaly, intellectual disability, hypogonadism, and anaemia. TTD-associated mutations typically cause unstable mutant proteins involved in various steps of gene expression, severely reducing steady-state mutant protein levels. However, to date, no such link to instability of gene-expression factors for TTD-associated mutations in MPLKIP/TTDN1 has been established. Here, we present seven additional TTD individuals with MPLKIP mutations from five consanguineous families, with a newly identified MPLKIP variant in one family. By mass spectrometry-based interaction proteomics, we demonstrate that MPLKIP interacts with core splicing factors and the lariat debranching protein DBR1. MPLKIP-deficient primary fibroblasts have reduced steady-state DBR1 protein levels. Using Human Skin Equivalents (HSEs), we observed impaired keratinocyte differentiation associated with compromised splicing and eventually, an imbalanced proteome affecting skin development and, interestingly, also the immune system. Our data show that MPLKIP, through its DBR1 stabilizing role, is implicated in mRNA splicing, which is of particular importance in highly differentiated tissue.

Sox2-positive dermal papilla cells specify hair follicle type in mammalian epidermis.

The dermal papilla comprises the specialised mesenchymal cells at the base of the hair follicle. Communication between dermal papilla cells and the overlying epithelium is essential for differentiation of the hair follicle lineages. We report that Sox2 is expressed in all dermal papillae at E16.5, but from E18.5 onwards expression is confined to a subset of dermal papillae. In postnatal skin, Sox2 is only expressed in the dermal papillae of guard/awl/auchene follicles, whereas CD133 is expressed both in guard/awl/auchene and in zigzag dermal papillae. Using transgenic mice that express GFP under the control of the Sox2 promoter, we isolated Sox2(+) (GFP(+)) CD133(+) cells and compared them with Sox2(-) (GFP(-)) CD133(+) dermal papilla cells. In addition to the 'core' dermal papilla gene signature, each subpopulation expressed distinct sets of genes. GFP(+) CD133(+) cells had upregulated Wnt, FGF and BMP pathways and expressed neural crest markers. In GFP(-) CD133(+) cells, the hedgehog, IGF, Notch and integrin pathways were prominent. In skin reconstitution assays, hair follicles failed to form when dermis was depleted of both GFP(+) CD133(+) and GFP(-) CD133(+) cells. In the absence of GFP(+) CD133(+) cells, awl/auchene hairs failed to form and only zigzag hairs were found. We have thus demonstrated a previously unrecognised heterogeneity in dermal papilla cells and shown that Sox2-positive cells specify particular hair follicle types.

Mapping dynamic histone acetylation patterns to gene expression in nanog-depleted murine embryonic stem cells.

Embryonic stem cells (ESC) have the potential to self-renew indefinitely and to differentiate into any of the three germ layers. The molecular mechanisms for self-renewal, maintenance of pluripotency and lineage specification are poorly understood, but recent results point to a key role for epigenetic mechanisms. In this study, we focus on quantifying the impact of histone 3 acetylation (H3K9,14ac) on gene expression in murine embryonic stem cells. We analyze genome-wide histone acetylation patterns and gene expression profiles measured over the first five days of cell differentiation triggered by silencing Nanog, a key transcription factor in ESC regulation. We explore the temporal and spatial dynamics of histone acetylation data and its correlation with gene expression using supervised and unsupervised statistical models. On a genome-wide scale, changes in acetylation are significantly correlated to changes in mRNA expression and, surprisingly, this coherence increases over time. We quantify the predictive power of histone acetylation for gene expression changes in a balanced cross-validation procedure. In an in-depth study we focus on genes central to the regulatory network of Mouse ESC, including those identified in a recent genome-wide RNAi screen and in the PluriNet, a computationally derived stem cell signature. We find that compared to the rest of the genome, ESC-specific genes show significantly more acetylation signal and a much stronger decrease in acetylation over time, which is often not reflected in a concordant expression change. These results shed light on the complexity of the relationship between histone acetylation and gene expression and are a step forward to dissect the multilayer regulatory mechanisms that determine stem cell fate.

Single-cell intracellular epitope and transcript detection reveals signal transduction dynamics.

To further our understanding of how biochemical information flows through cells upon external stimulation, we require single-cell multi-omics methods that concurrently map changes in (phospho)protein levels across signaling networks and the associated gene expression profiles. Here, we present quantification of RNA and intracellular epitopes by sequencing (QuRIE-seq), a droplet-based platform for single-cell RNA and intra- and extracellular (phospho)protein quantification through sequencing. We applied QuRIE-seq to quantify cell-state changes at both the signaling and the transcriptome level after 2-, 4-, 6-, 60-, and 180-min stimulation of the B cell receptor pathway in Burkitt lymphoma cells. Using the multi-omics factor analysis (MOFA+) framework, we delineated changes in single-cell (phospho)protein and gene expression patterns over multiple timescales and revealed the effect of an inhibitory drug (ibrutinib) on signaling and gene expression landscapes.

Human Ccr4-Not complexes contain variable deadenylase subunits.

The Ccr4-Not complex is evolutionarily conserved and important for regulation of mRNA synthesis and decay. The composition of the yeast complex has been well described. Orthologues of the yeast Ccr4-Not components have been identified in human cells including multiple subunits with mRNA deadenylase activity. In the present study, we examine the composition of the human Ccr4-Not complex in an in-depth proteomic approach using stable cell lines expressing tagged CNOT proteins. We find at least four different variants of the human complex, consisting of seven stable core proteins and mutually exclusive associated mRNA deadenylase subunits. Interestingly, human CNOT4 is in a separate approximately 200 kDa complex. Furthermore, analyses of associated proteins indicate involvement of Ccr4-Not complexes in splicing, transport and localization of RNA molecules. Taken together, human Ccr4-Not complexes are heterogeneous in composition owing to differences in their deadenylase subunits, which may reflect the multi-functionality of these complexes in cellular processes.

The corepressor NCOR1 and OCT4 facilitate early reprogramming by suppressing fibroblast gene expression.

Reprogramming somatic cells to induced pluripotent stem cells (iPSC) succeeds only in a small fraction of cells within the population. Reprogramming occurs in distinctive stages, each facing its own bottlenecks. It initiates with overexpression of transcription factors OCT4, SOX2, KLF4 and c-MYC (OSKM) in somatic cells such as mouse embryonic fibroblasts (MEFs). OSKM bind chromatin, silencing the somatic identity and starting the stepwise reactivation of the pluripotency programme. However, inefficient suppression of the somatic lineage leads to unwanted epigenetic memory from the tissue of origin, even in successfully generated iPSCs. Thus, it is essential to shed more light on chromatin regulators and processes involved in dissolving the somatic identity. Recent work characterised the role of transcriptional corepressors NCOR1 and NCOR2 (also known as NCoR and SMRT), showing that they cooperate with c-MYC to silence pluripotency genes during late reprogramming stages. NCOR1/NCOR2 were also proposed to be involved in silencing fibroblast identity, however it is unclear how this happens. Here, we shed light on the role of NCOR1 in early reprogramming. We show that siRNA-mediated ablation of NCOR1 and OCT4 results in very similar phenotypes, including transcriptomic changes and highly correlated high-content colony phenotypes. Both NCOR1 and OCT4 bind to promoters co-occupied by c-MYC in MEFs. During early reprogramming, downregulation of one group of somatic MEF-expressed genes requires both NCOR1 and OCT4, whereas another group of MEF-expressed genes is downregulated by NCOR1 but not OCT4. Our data suggest that NCOR1, assisted by OCT4 and c-MYC, facilitates transcriptional repression of genes with high expression in MEFs, which is necessary to bypass an early reprogramming block; this way, NCOR1 facilitates early reprogramming progression.

Regulation of histone H3K4 tri-methylation and PAF complex recruitment by the Ccr4-Not complex.

Efficient transcription is linked to modification of chromatin. For instance, tri-methylation of lysine 4 on histone H3 (H3K4) strongly correlates with transcriptional activity and is regulated by the Bur1/2 kinase complex. We found that the evolutionarily conserved Ccr4-Not complex is involved in establishing H3K4 tri-methylation in Saccharomyces cerevisiae. We observed synthetic lethal interactions of Ccr4-Not components with BUR1 and BUR2. Further analysis indicated that the genes encoding the Not-proteins are essential for efficient regulation of H3K4me3, but not H3K4me1/2, H3K36me2 or H3K79me2/3 levels. Moreover, regulation of H3K4me3 levels by NOT4 is independent of defects in RNA polymerase II loading. We found NOT4 to be important for ubiquitylation of histone H2B via recruitment of the PAF complex, but not for recruitment or activation of the Bur1/2 complex. These results suggest a mechanism in which the Ccr4-Not complex functions parallel to or downstream of the Bur1/2 kinase to facilitate H3K4me3 via PAF complex recruitment.

The interplay of chromatin and transcription factors during cell fate transitions in development and reprogramming.

Reprogramming to induced pluripotency through expression of OCT4, SOX2, KLF4, MYC (OSKM) factors is often considered the dedifferentiation of somatic cells. This would suggest that reprogramming represents the reversal of embryonic differentiation. Indeed, molecular events involving the activity of the pluripotency network occur in opposite directions. However, reprogramming and development substantially differ as OSKM bind to accessible regulatory elements in the genome of somatic cells due to their overexpression, including regulatory elements never bound by these factors during normal differentiation. In addition, rewiring the transcriptional network back to pluripotency involves overcoming molecular barriers that protect or stabilize the somatic identity, whereas extrinsic and intrinsic cues will drive differentiation in an energetically favorable landscape in the embryo. This review focuses on how cell fate transitions in reprogramming and development are differentially governed by interactions between transcription factors and chromatin. We also discuss how these interactions shape chromatin architecture and the transcriptional output. Major technological advances have resulted in a better understanding of both differentiation and reprogramming, which is essential to exploit reprogramming regimes for regenerative medicine.

BLNCR is a long non-coding RNA adjacent to integrin beta-1 that is rapidly lost during epidermal progenitor cell differentiation.

As our understanding of transcriptional regulation improves so does our appreciation of its complexity. Both coding and (long) non-coding RNAs provide cells with multiple levels of control and thereby flexibility to adapt gene expression to the environment. However, few long non-coding RNAs (lncRNAs) have been studied in human epidermal stem cells. Here, we characterized the expression of 26 lncRNAs in human epidermal keratinocytes, 7 of which we found to be dynamically expressed during differentiation. We performed in depth analysis of a lncRNA located proximal to the epidermal stem cell marker integrin beta-1 (ITGB1) and transcribed in the opposite direction. We dubbed this gene Beta1-adjacent long non-coding RNA, or BLNCR, and found that its expression is regulated by p63 and AP1 transcription factors. Furthermore, BLNCR expression is regulated downstream the integrin and EGF signaling pathways that are key to epidermal stem cell maintenance. Finally, we found that BLNCR expression is rapidly reduced upon induction of differentiation, preceding the down regulation of integrin beta-1 expression. These dynamics closely mirror the loss of proliferative and adhesion capacity of epidermal stem cells in colony formation assays. Together, these results suggest that loss of BLNCR expression marks the switch from a proliferative state towards terminal differentiation in human epidermal stem cells.

Combined quantification of intracellular (phospho-)proteins and transcriptomics from fixed single cells.

Environmental stimuli often lead to heterogeneous cellular responses and transcriptional output. We developed single-cell RNA and Immunodetection (RAID) to allow combined analysis of the transcriptome and intracellular (phospho-)proteins from fixed single cells. RAID successfully recapitulated differentiation-state changes at the protein and mRNA level in human keratinocytes. Furthermore, we show that differentiated keratinocytes that retain high phosphorylated FAK levels, a feature associated with stem cells, also express a selection of stem cell associated transcripts. Our data demonstrates that RAID allows investigation of heterogeneous cellular responses to environmental signals at the mRNA and phospho-proteome level.

WDR5, BRCA1, and BARD1 Co-regulate the DNA Damage Response and Modulate the Mesenchymal-to-Epithelial Transition during Early Reprogramming.

Differentiated cells are epigenetically stable, but can be reprogrammed to pluripotency by expression of the OSKM transcription factors. Despite significant effort, relatively little is known about the cellular requirements for reprogramming and how they affect the properties of induced pluripotent stem cells. We have performed high-content screening with small interfering RNAs targeting 300 chromatin-associated factors and extracted colony-level quantitative features. This revealed five morphological phenotypes in early reprogramming, including one displaying large round colonies exhibiting an early block of reprogramming. Using RNA sequencing, we identified transcriptional changes associated with these phenotypes. Furthermore, double knockdown epistasis experiments revealed that BRCA1, BARD1, and WDR5 functionally interact and are required for the DNA damage response. In addition, the mesenchymal-to-epithelial transition is affected in Brca1, Bard1, and Wdr5 knockdowns. Our data provide a resource of chromatin-associated factors in early reprogramming and underline colony morphology as an important high-dimensional readout for reprogramming quality.

Modulation of Ubc4p/Ubc5p-mediated stress responses by the RING-finger-dependent ubiquitin-protein ligase Not4p in Saccharomyces cerevisiae.

The Ccr4-Not complex consists of nine subunits and acts as a regulator of mRNA biogenesis in Saccharomyces cerevisiae. The human ortholog of yeast NOT4, CNOT4, displays UbcH5B-dependent ubiquitin-protein ligase (E3 ligase) activity in a reconstituted in vitro system. However, an in vivo role for this enzymatic activity has not been identified. Site-directed mutagenesis of the RING finger of yeast Not4p identified residues required for interaction with Ubc4p and Ubc5p, the yeast orthologs of UbcH5B. Subsequent in vitro assays with purified Ccr4-Not complexes showed Not4p-mediated E3 ligase activity, which was dependent on the interaction with Ubc4p. To investigate the in vivo relevance of this activity, we performed synthetic genetic array (SGA) analyses using not4Delta and not4L35A alleles. This indicates involvement of the RING finger of Not4p in transcription, ubiquitylation, and DNA damage responses. In addition, we found a phenotypic overlap between deletions of UBC4 and mutants encoding single-amino-acid substitutions of the RING finger of Not4p. Together, our results show that Not4p functions as an E3 ligase by modulating Ubc4p/Ubc5p-mediated stress responses in vivo.

Menin links estrogen receptor activation to histone H3K4 trimethylation.

The product of the multiple endocrine neoplasia type 1 (MEN1) tumor suppressor gene, menin, is an integral component of MLL1/MLL2 histone methyltransferase complexes specific for Lys4 of histone H3 (H3K4). We show that menin is a transcriptional coactivator of the nuclear receptors for estrogen and vitamin D. Activation of the endogenous estrogen-responsive TFF1 (pS2) gene results in promoter recruitment of menin and in elevated trimethylation of H3K4. Knockdown of menin reduces both activated TFF1 (pS2) transcription and H3K4 trimethylation. In addition, menin can directly interact with the estrogen receptor-alpha (ERalpha) in a hormone-dependent manner. The majority of disease-related MEN1 mutations prevent menin-ERalpha interaction. Importantly, ERalpha-interacting mutants are also defective in coactivator function. Our results indicate that menin is a critical link between recruitment of histone methyltransferase complexes and nuclear receptor-mediated transcription.

Single-Cell ID-seq Reveals Dynamic BMP Pathway Activation Upstream of the MAF/MAFB-Program in Epidermal Differentiation.

Epidermal homeostasis requires balanced and coordinated adult stem cell renewal and differentiation. These processes are controlled by both extracellular signaling and by cell intrinsic transcription regulatory networks, yet how these control mechanisms are integrated to achieve this is unclear. Here, we developed single-cell Immuno-Detection by sequencing (scID-seq) and simultaneously measured 69 proteins (including 34 phosphorylated epitopes) at single-cell resolution to study the activation state of signaling pathways during human epidermal differentiation. Computational pseudo-timing inference revealed dynamic activation of the JAK-STAT, WNT, and BMP pathways along the epidermal differentiation trajectory. We found that during differentiation, cells start producing BMP2-ligands and activate the canonical intracellular effectors SMAD1/5/9. Mechanistically, the BMP pathway is responsible for activating the MAF/MAFB/ZNF750 transcription factor network to drive late-stage epidermal differentiation. Our work indicates that incorporating signaling pathway activation into this transcription regulatory network enables coordination of transcription programs during epidermal differentiation.

Splicing and Chromatin Factors Jointly Regulate Epidermal Differentiation.

Epidermal homeostasis requires balanced progenitor cell proliferation and loss of differentiated cells from the epidermal surface. During this process, cells undergo major changes in their transcriptional programs to accommodate new cellular functions. We found that transcriptional and post-transcriptional mechanisms underlying these changes jointly control genes involved in cell adhesion, a key process in epidermal maintenance. Using siRNA-based perturbation screens, we identified DNA and/or RNA binding regulators of epidermal differentiation. Computational modeling and experimental validation identified functional interactions between the matrin-type 2 zinc-finger protein ZMAT2 and the epigenetic modifiers ING5, SMARCA5, BRD1, UHRF1, BPTF, and SMARCC2. ZMAT2 is an interactor of the pre-spliceosome that is required to keep cells in an undifferentiated, proliferative state. RNA immunoprecipitation and transcriptome-wide RNA splicing analysis showed that ZMAT2 associates with and regulates transcripts involved in cell adhesion in conjunction with ING5. Thus, joint control by splicing regulation, histone, and DNA modification is important to maintain epidermal cells in an undifferentiated state.

Mutant p63 Affects Epidermal Cell Identity through Rewiring the Enhancer Landscape.

Transcription factor p63 is a key regulator of epidermal keratinocyte proliferation and differentiation. Mutations in the p63 DNA-binding domain are associated with ectrodactyly, ectodermal dysplasia, and cleft lip/palate (EEC) syndrome. However, the underlying molecular mechanism of these mutations remains unclear. Here, we characterized the transcriptome and epigenome of p63 mutant keratinocytes derived from EEC patients. The transcriptome of p63 mutant keratinocytes deviated from the normal epidermal cell identity. Epigenomic analyses showed an altered enhancer landscape in p63 mutant keratinocytes contributed by loss of p63-bound active enhancers and unexpected gain of enhancers. The gained enhancers were frequently bound by deregulated transcription factors such as RUNX1. Reversing RUNX1 overexpression partially rescued deregulated gene expression and the altered enhancer landscape. Our findings identify a disease mechanism whereby mutant p63 rewires the enhancer landscape and affects epidermal cell identity, consolidating the pivotal role of p63 in controlling the enhancer landscape of epidermal keratinocytes.

Immuno-detection by sequencing enables large-scale high-dimensional phenotyping in cells.

Cell-based small molecule screening is an effective strategy leading to new medicines. Scientists in the pharmaceutical industry as well as in academia have made tremendous progress in developing both large-scale and smaller-scale screening assays. However, an accessible and universal technology for measuring large numbers of molecular and cellular phenotypes in many samples in parallel is not available. Here we present the immuno-detection by sequencing (ID-seq) technology that combines antibody-based protein detection and DNA-sequencing via DNA-tagged antibodies. We use ID-seq to simultaneously measure 70 (phospho-)proteins in primary human epidermal stem cells to screen the effects of ~300 kinase inhibitor probes to characterise the role of 225 kinases. The results show an association between decreased mTOR signalling and increased differentiation and uncover 13 kinases potentially regulating epidermal renewal through distinct mechanisms. Taken together, our work establishes ID-seq as a flexible solution for large-scale high-dimensional phenotyping in fixed cell populations.

DNA damage and replication stress induced transcription of RNR genes is dependent on the Ccr4-Not complex.

Genetic experiments have indicated a role for the Ccr4-Not complex in the response to hydroxyurea (HU) induced replication stress and ionizing radiation in yeast. This response includes transcriptional induction of the four genes constituting the ribonucleotide reductase (RNR) enzymatic complex, RNR1-4 and degradation of its inhibitor, Sml1p. The Ccr4-Not complex has originally been described as a negative regulator of RNA polymerase II (pol II) transcription, but it has also been implicated in mRNA turnover and protein ubiquitination. We investigated the mechanism of the HU sensitivity conferred by mutation of CCR4-NOT genes. We found that the ubiquitin protein ligase activity of Not4p does not play a role in HU induced Sml1p degradation. We show, however, that the HU sensitivity of ccr4-not mutant strains correlated very well with a defect in accumulation of RNR2, RNR3 and RNR4 mRNA after HU or methyl-methane sulfonate (MMS) treatment. Chromatin immunoprecipitation (ChIP) experiments show that TBP, pol II and Set1p recruitment to the activated RNR3 locus is defective in cells lacking NOT4. Moreover, RNR3-promoter activity is not induced by HU in these cells. Our experiments show that induction of RNR gene transcription is defective in ccr4-not mutant strains, providing an explanation for their sensitivity to HU.