Tetraspanin 5 (TSPAN5), a Novel Gatekeeper of the Tumor Suppressor DLC1 and Myocardin-Related Transcription Factors (MRTFs), Controls HCC Growth and Senescence.

Human hepatocellular carcinoma (HCC) is among the most lethal and common cancers in the human population, and new molecular targets for therapeutic intervention are urgently needed. Deleted in liver cancer 1 (DLC1) was originally identified as a tumor suppressor gene in human HCC. DLC1 is a Rho-GTPase-activating protein (RhoGAP) which accelerates the return of RhoGTPases to an inactive state. We recently described that the restoration of DLC1 expression induces cellular senescence. However, this principle is not amenable to direct therapeutic targeting. We therefore performed gene expression profiling for HepG2 cells depleted of DLC1 to identify druggable gene targets mediating the effects of DLC1 on senescence induction. This approach revealed that versican (VCAN), tetraspanin 5 (TSPAN5) and N-cadherin (CDH2) were strongly upregulated upon DLC1 depletion in HCC cells, but only TSPAN5 affected the proliferation of HCC cells and human HCC. The depletion of TSPAN5 induced oncogene-induced senescence (OIS), mediated by the p16INK4a/pRb pathways. Mechanistically, silencing TSPAN5 reduced actin polymerization and thereby myocardin-related transcription factor A- filamin A (MRTF-A-FLNA) complex formation, resulting in decreased expression of MRTF/SRF-dependent target genes and senescence induction in vitro and in vivo. Our results identify TSPAN5 as a novel druggable target for HCC.

EBF1 binds to EBNA2 and promotes the assembly of EBNA2 chromatin complexes in B cells.

Epstein-Barr virus (EBV) infection converts resting human B cells into permanently proliferating lymphoblastoid cell lines (LCLs). The Epstein-Barr virus nuclear antigen 2 (EBNA2) plays a key role in this process. It preferentially binds to B cell enhancers and establishes a specific viral and cellular gene expression program in LCLs. The cellular DNA binding factor CBF1/CSL serves as a sequence specific chromatin anchor for EBNA2. The ubiquitous expression of this highly conserved protein raises the question whether additional cellular factors might determine EBNA2 chromatin binding selectively in B cells. Here we used CBF1 deficient B cells to identify cellular genes up or downregulated by EBNA2 as well as CBF1 independent EBNA2 chromatin binding sites. Apparently, CBF1 independent EBNA2 target genes and chromatin binding sites can be identified but are less frequent than CBF1 dependent EBNA2 functions. CBF1 independent EBNA2 binding sites are highly enriched for EBF1 binding motifs. We show that EBNA2 binds to EBF1 via its N-terminal domain. CBF1 proficient and deficient B cells require EBF1 to bind to CBF1 independent binding sites. Our results identify EBF1 as a co-factor of EBNA2 which conveys B cell specificity to EBNA2.

Calcium-mediated actin reset (CaAR) mediates acute cell adaptations.

Actin has well established functions in cellular morphogenesis. However, it is not well understood how the various actin assemblies in a cell are kept in a dynamic equilibrium, in particular when cells have to respond to acute signals. Here, we characterize a rapid and transient actin reset in response to increased intracellular calcium levels. Within seconds of calcium influx, the formin INF2 stimulates filament polymerization at the endoplasmic reticulum (ER), while cortical actin is disassembled. The reaction is then reversed within a few minutes. This Calcium-mediated actin reset (CaAR) occurs in a wide range of mammalian cell types and in response to many physiological cues. CaAR leads to transient immobilization of organelles, drives reorganization of actin during cell cortex repair, cell spreading and wound healing, and induces long-lasting changes in gene expression. Our findings suggest that CaAR acts as fundamental facilitator of cellular adaptations in response to acute signals and stress.

Periodic mRNA synthesis and degradation co-operate during cell cycle gene expression.

During the cell cycle, the levels of hundreds of mRNAs change in a periodic manner, but how this is achieved by alterations in the rates of mRNA synthesis and degradation has not been studied systematically. Here, we used metabolic RNA labeling and comparative dynamic transcriptome analysis (cDTA) to derive mRNA synthesis and degradation rates every 5 min during three cell cycle periods of the yeast Saccharomyces cerevisiae. A novel statistical model identified 479 genes that show periodic changes in mRNA synthesis and generally also periodic changes in their mRNA degradation rates. Peaks of mRNA degradation generally follow peaks of mRNA synthesis, resulting in sharp and high peaks of mRNA levels at defined times during the cell cycle. Whereas the timing of mRNA synthesis is set by upstream DNA motifs and their associated transcription factors (TFs), the synthesis rate of a periodically expressed gene is apparently set by its core promoter.

One Hand Clapping: detection of condition-specific transcription factor interactions from genome-wide gene activity data.

We present One Hand Clapping (OHC), a method for the detection of condition-specific interactions between transcription factors (TFs) from genome-wide gene activity measurements. OHC is based on a mapping between transcription factors and their target genes. Given a single case-control experiment, it uses a linear regression model to assess whether the common targets of two arbitrary TFs behave differently than expected from the genes targeted by only one of the TFs. When applied to osmotic stress data in S. cerevisiae, OHC produces consistent results across three types of expression measurements: gene expression microarray data, RNA Polymerase II ChIP-chip binding data and messenger RNA synthesis rates. Among the eight novel, condition-specific TF pairs, we validate the interaction between Gcn4p and Arr1p experimentally. We apply OHC to a large gene activity dataset in S. cerevisiae and provide a compendium of condition-specific TF interactions.

MC EMiNEM maps the interaction landscape of the Mediator.

The Mediator is a highly conserved, large multiprotein complex that is involved essentially in the regulation of eukaryotic mRNA transcription. It acts as a general transcription factor by integrating regulatory signals from gene-specific activators or repressors to the RNA Polymerase II. The internal network of interactions between Mediator subunits that conveys these signals is largely unknown. Here, we introduce MC EMiNEM, a novel method for the retrieval of functional dependencies between proteins that have pleiotropic effects on mRNA transcription. MC EMiNEM is based on Nested Effects Models (NEMs), a class of probabilistic graphical models that extends the idea of hierarchical clustering. It combines mode-hopping Monte Carlo (MC) sampling with an Expectation-Maximization (EM) algorithm for NEMs to increase sensitivity compared to existing methods. A meta-analysis of four Mediator perturbation studies in Saccharomyces cerevisiae, three of which are unpublished, provides new insight into the Mediator signaling network. In addition to the known modular organization of the Mediator subunits, MC EMiNEM reveals a hierarchical ordering of its internal information flow, which is putatively transmitted through structural changes within the complex. We identify the N-terminus of Med7 as a peripheral entity, entailing only local structural changes upon perturbation, while the C-terminus of Med7 and Med19 appear to play a central role. MC EMiNEM associates Mediator subunits to most directly affected genes, which, in conjunction with gene set enrichment analysis, allows us to construct an interaction map of Mediator subunits and transcription factors.

The spt5 C-terminal region recruits yeast 3' RNA cleavage factor I.

During transcription elongation, RNA polymerase II (Pol II) binds the general elongation factor Spt5. Spt5 contains a repetitive C-terminal region (CTR) that is required for cotranscriptional recruitment of the Paf1 complex (D. L. Lindstrom et al., Mol. Cell. Biol. 23:1368-1378, 2003; Z. Zhang, J. Fu, and D. S. Gilmour, Genes Dev. 19:1572-1580, 2005). Here we report a new role of the Spt5 CTR in the recruitment of 3' RNA-processing factors. Chromatin immunoprecipitation (ChIP) revealed that the Spt5 CTR is required for normal recruitment of pre-mRNA cleavage factor I (CFI) to the 3' ends of Saccharomyces cerevisiae genes. RNA contributes to CFI recruitment, as RNase treatment prior to ChIP further decreases CFI ChIP signals. Genome-wide ChIP profiling detected occupancy peaks of CFI subunits around 100 nucleotides downstream of the polyadenylation (pA) sites of genes. CFI recruitment to this defined region may result from simultaneous binding to the Spt5 CTR, to nascent RNA containing the pA sequence, and to the elongating Pol II isoform that is phosphorylated at serine 2 (S2) residues in its C-terminal domain (CTD). Consistent with this model, the CTR interacts with CFI in vitro but is not required for pA site recognition and transcription termination in vivo.

Measurement of genome-wide RNA synthesis and decay rates with Dynamic Transcriptome Analysis (DTA).

Standard transcriptomics measures total cellular RNA levels. Our understanding of gene regulation would be greatly improved if we could measure RNA synthesis and decay rates on a genome-wide level. To that end, the Dynamic Transcriptome Analysis (DTA) method has been developed. DTA combines metabolic RNA labeling with standard transcriptomics to measure RNA synthesis and decay rates in a precise and non-perturbing manner. Here, we present the open source R/Bioconductor software package DTA. It implements all required bioinformatics steps that allow the accurate absolute quantification and comparison of RNA turnover.

Global DNA hypomethylation prevents consolidation of differentiation programs and allows reversion to the embryonic stem cell state.

DNA methylation patterns change dynamically during mammalian development and lineage specification, yet scarce information is available about how DNA methylation affects gene expression profiles upon differentiation. Here we determine genome-wide transcription profiles during undirected differentiation of severely hypomethylated (Dnmt1⁻/⁻) embryonic stem cells (ESCs) as well as ESCs completely devoid of DNA methylation (Dnmt1⁻/⁻;Dnmt3a⁻/⁻;Dnmt3b⁻/⁻ or TKO) and assay their potential to transit in and out of the ESC state. We find that the expression of only few genes mainly associated with germ line function and the X chromosome is affected in undifferentiated TKO ESCs. Upon initial differentiation as embryoid bodies (EBs) wild type, Dnmt1⁻/⁻ and TKO cells downregulate pluripotency associated genes and upregulate lineage specific genes, but their transcription profiles progressively diverge upon prolonged EB culture. While Oct4 protein levels are completely and homogeneously suppressed, transcription of Oct4 and Nanog is not completely silenced even at late stages in both Dnmt1⁻/⁻ and TKO EBs. Despite late wild type and Dnmt1⁻/⁻ EBs showing a much higher degree of concordant expression, after EB dissociation and replating under pluripotency promoting conditions both Dnmt1⁻/⁻ and TKO cells, but not wild type cells rapidly revert to expression profiles typical of undifferentiated ESCs. Thus, while DNA methylation seems not to be critical for initial activation of differentiation programs, it is crucial for permanent restriction of developmental fate during differentiation.

A conserved GA element in TATA-less RNA polymerase II promoters.

Initiation of RNA polymerase (Pol) II transcription requires assembly of the pre-initiation complex (PIC) at the promoter. In the classical view, PIC assembly starts with binding of the TATA box-binding protein (TBP) to the TATA box. However, a TATA box occurs in only 15% of promoters in the yeast Saccharomyces cerevisiae, posing the question how most yeast promoters nucleate PIC assembly. Here we show that one third of all yeast promoters contain a novel conserved DNA element, the GA element (GAE), that generally does not co-occur with the TATA box. The distance of the GAE to the transcription start site (TSS) resembles the distance of the TATA box to the TSS. The TATA-less TMT1 core promoter contains a GAE, recruits TBP, and supports formation of a TBP-TFIIB-DNA-complex. Mutation of the promoter region surrounding the GAE abolishes transcription in vivo and in vitro. A 32-nucleotide promoter region containing the GAE can functionally substitute for the TATA box in a TATA-containing promoter. This identifies the GAE as a conserved promoter element in TATA-less promoters.

Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast.

To obtain rates of mRNA synthesis and decay in yeast, we established dynamic transcriptome analysis (DTA). DTA combines non-perturbing metabolic RNA labeling with dynamic kinetic modeling. DTA reveals that most mRNA synthesis rates are around several transcripts per cell and cell cycle, and most mRNA half-lives range around a median of 11 min. DTA can monitor the cellular response to osmotic stress with higher sensitivity and temporal resolution than standard transcriptomics. In contrast to monotonically increasing total mRNA levels, DTA reveals three phases of the stress response. During the initial shock phase, mRNA synthesis and decay rates decrease globally, resulting in mRNA storage. During the subsequent induction phase, both rates increase for a subset of genes, resulting in production and rapid removal of stress-responsive mRNAs. During the recovery phase, decay rates are largely restored, whereas synthesis rates remain altered, apparently enabling growth at high salt concentration. Stress-induced changes in mRNA synthesis rates are predicted from gene occupancy with RNA polymerase II. DTA-derived mRNA synthesis rates identified 16 stress-specific pairs/triples of cooperative transcription factors, of which seven were known. Thus, DTA realistically monitors the dynamics in mRNA metabolism that underlie gene regulatory systems.

Identification, structure, and functional requirement of the Mediator submodule Med7N/31.

Mediator is a modular multiprotein complex required for regulated transcription by RNA polymerase (Pol) II. Here, we show that the middle module of the Mediator core contains a submodule of unique structure and function that comprises the N-terminal part of subunit Med7 (Med7N) and the highly conserved subunit Med31 (Soh1). The Med7N/31 submodule shows a conserved novel fold, with two proline-rich stretches in Med7N wrapping around the right-handed four-helix bundle of Med31. In vitro, Med7N/31 is required for activated transcription and can act in trans when added exogenously. In vivo, Med7N/31 has a predominantly positive function on the expression of a specific subset of genes, including genes involved in methionine metabolism and iron transport. Comparative phenotyping and transcriptome profiling identify specific and overlapping functions of different Mediator submodules.

Identification of the rapamycin-sensitive phosphorylation sites within the Ser/Thr-rich domain of the yeast Npr1 protein kinase.

The Saccharomyces cerevisae nitrogen permease reactivator Npr1 is a hyperphosphorylated protein that belongs to a fungus-specific family of Ser/Thr protein kinases dedicated to the regulation of plasma membrane transporters. Its activity is regulated by the TOR (target of rapamycin) signalling pathway. Inhibition of the TOR proteins by treating yeast cells with the immunosuppressant drug rapamycin promotes rapid dephosphorylation of Npr1. To identify the rapamycin-sensitive phosphorylation sites in yeast Npr1, glutathione-S-transferase (GST)-tagged Npr1 was isolated from untreated or rapamycin-treated cells, and analyzed by mass spectrometry. Here, we report for the first time 22 phosphorylation sites that are clustered in two regions of the N-terminal serine-rich domain. All phosphorylation sites, except two, were found to be rapamycin-sensitive. Some phosphorylation sites are contained in motifs that closely resemble those in mammalian S6K (serines followed by a tyrosine or a phenylalanine) and 4E-BP1 (serines followed by a proline). Other sites, such as serines followed by Ala, Asn, Gln, His, Ile, Leu, or Val, appear to define new motifs. Thus, TOR controls an unusually broad array of phosphorylation sites in Npr1. In addition to phosphorylation by upstream kinases, Npr1 undergoes autophosphorylation that was mapped to three distinct serines in the N-terminal domain of which Ser257 appears to be the main autophosphorylation site. Site-directed mutagenesis confirmed the mass spectral assignments of the autophosphorylation sites and shows that Ser257 is specifically involved in forming an in vitro substrate-binding site.

Genome-associated RNA polymerase II includes the dissociable Rpb4/7 subcomplex.

Yeast RNA polymerase (Pol) II consists of a 10-subunit core enzyme and the Rpb4/7 subcomplex, which is dispensable for catalytic activity and dissociates in vitro. To investigate whether Rpb4/7 is an integral part of DNA-associated Pol II in vivo, we used chromatin immunoprecipitation coupled to high resolution tiling microarray analysis. We show that the genome-wide occupancy profiles for Rpb7 and the core subunit Rpb3 are essentially identical. Thus, the complete Pol II associates with DNA in vivo, consistent with functional roles of Rpb4/7 throughout the transcription cycle.

TOR1 and TOR2 have distinct locations in live cells.

TOR is a structurally and functionally conserved Ser/Thr kinase found in two multiprotein complexes that regulate many cellular processes to control cell growth. Although extensively studied, the localization of TOR is still ambiguous, possibly because endogenous TOR in live cells has not been examined. Here, we examined the localization of green fluorescent protein (GFP) tagged, endogenous TOR1 and TOR2 in live S. cerevisiae cells. A DNA cassette encoding three copies of green fluorescent protein (3XGFP) was inserted in the TOR1 gene (at codon D330) or the TOR2 gene (at codon N321). The TORs were tagged internally because TOR1 or TOR2 tagged at the N or C terminus was not functional. The TOR1(D330-3XGFP) strain was not hypersensitive to rapamycin, was not cold sensitive, and was not resistant to manganese toxicity caused by the loss of Pmr1, all indications that TOR1-3XGFP was expressed and functional. TOR2-3XGFP was functional, as TOR2 is an essential gene and TOR2(N321-3XGFP) haploid cells were viable. Thus, TOR1 and TOR2 retain function after the insertion of 748 amino acids in a variable region of their noncatalytic domain. The localization patterns of TOR1-3XGFP and TOR2-3XGFP were documented by imaging of live cells. TOR1-3XGFP was diffusely cytoplasmic and concentrated near the vacuolar membrane. The TOR2-3XGFP signal was cytoplasmic but predominately in dots at the plasma membrane. Thus, TOR1 and TOR2 have distinct localization patterns, consistent with the regulation of cellular processes as part of two different complexes.

Regulation of ribosome biogenesis: where is TOR?

Li et al., (2006) have shown that TOR complex 1 in yeast binds directly to the rDNA promoter and thereby activates Pol I-dependent synthesis of 35S RNA. This is an important advance in the understanding of how ribosome biogenesis is regulated in response to environmental conditions.

Mutual antagonism of target of rapamycin and calcineurin signaling.

Growth and stress are generally incompatible states. Stressed cells adapt to an insult by restraining growth, and conversely, growing cells keep stress responses at bay. This is evident in many physiological settings, including for example, the effect of stress on the immune or nervous system, but the underlying signaling mechanisms mediating such mutual antagonism are poorly understood. In eukaryotes, a central activator of cell growth is the protein kinase target of rapamycin (TOR) and its namesake signaling network. Calcineurin is a conserved, Ca(2+)/calmodulin-dependent protein phosphatase and target of the immunosuppressant FK506 (tacrolimus) that is activated in yeast during stress to promote cell survival. Here we show yeast mutants defective for TOR complex 2 (TORC2) or the essential homologous TORC2 effectors, SLM1 and SLM2, exhibited constitutive activation of calcineurin-dependent transcription and actin depolarization. Conversely, cells defective in calcineurin exhibited SLM1 hyperphosphorylation and enhanced interaction between TORC2 and SLM1. Furthermore, a mutant SLM1 protein (SLM1(DeltaC14)) lacking a sequence related to the consensus calcineurin docking site (PxIxIT) was insensitive to calcineurin, and SLM1(Delta)(C14) slm2 mutant cells were hypersensitive to oxidative stress. Thus, TORC2-SLM signaling negatively regulates calcineurin, and calcineurin negatively regulates TORC2-SLM. These findings provide a molecular basis for the mutual antagonism of growth and stress.

The expanding TOR signaling network.

Cell growth (increase in cell mass or size) is tightly coupled to nutrient availability, growth factors and the energy status of the cell. The target of rapamycin (TOR) integrates all three inputs to control cell growth. The discovery of upstream regulators of TOR (AMPK, the TSC1-TSC2 complex and Rheb) has provided new insights into the mechanism by which TOR integrates its various inputs. A recent finding in flies reveals that TOR controls not only growth of the cell in which it resides (cell-autonomous growth) but also the growth of distant cells, thereby determining organ and organism size in addition to the size of isolated cells. In yeast and mammals, the identification of two structurally and functionally distinct multiprotein TOR complexes (TORC1 and TORC2) has provided a molecular basis for the complexity of TOR signaling. Furthermore, TOR has emerged as a regulator of growth-related processes such as development, aging and the response to hypoxia. Thus, TOR is part of an intra- and inter-cellular signaling network with a remarkably broad role in eukaryotic biology.

TOR regulates ribosomal protein gene expression via PKA and the Forkhead transcription factor FHL1.

The regulation of ribosome biogenesis in response to environmental conditions is a key aspect of cell growth control. Ribosomal protein (RP) genes are regulated by the nutrient-sensitive, conserved target of rapamycin (TOR) signaling pathway. TOR controls the subcellular localization of protein kinase A (PKA) and the PKA-regulated kinase YAK1. However, the target transcription factor(s) of the TOR-PKA pathway are unknown. We show that regulation of RP gene transcription via TOR and PKA in yeast involves the Forkhead-like transcription factor FHL1 and the two cofactors IFH1 (a coactivator) and CRF1 (a corepressor). TOR, via PKA, negatively regulates YAK1 and maintains CRF1 in the cytoplasm. Upon TOR inactivation, activated YAK1 phosphorylates and activates CRF1. Phosphorylated CRF1 accumulates in the nucleus and competes with IFH1 for binding to FHL1 at RP gene promoters, and thereby inhibits transcription of RP genes. Thus, we describe a signaling mechanism linking an environmental sensor to ribosome biogenesis.

Rank Difference Analysis of Microarrays (RDAM), a novel approach to statistical analysis of microarray expression profiling data.

BACKGROUND: A key step in the analysis of microarray expression profiling data is the identification of genes that display statistically significant changes in expression signals between two biological conditions. RESULTS: We describe a new method, Rank Difference Analysis of Microarrays (RDAM), which estimates the total number of truly varying genes and assigns a p-value to each signal variation. Information on a group of differentially expressed genes includes the sensitivity and the false discovery rate. We demonstrate the feasibility and efficiency of our approach by applying it to a large synthetic expression data set and to a biological data set obtained by comparing vegetatively-growing wild type and tor2-mutant yeast strains. In both cases we observed a significant improvement of the power of analysis when our method is compared to another popular nonparametric method. CONCLUSIONS: This study provided a valuable new statistical method to analyze microarray data. We conclude that the good quality of the results obtained by RDAM is mainly due to the quasi-perfect equalization of variation distribution, which is related to the standardization procedure used and to the measurement of variation by rank difference.