MasterPATH: network analysis of functional genomics screening data.

BACKGROUND: Functional genomics employs several experimental approaches to investigate gene functions. High-throughput techniques, such as loss-of-function screening and transcriptome profiling, allow to identify lists of genes potentially involved in biological processes of interest (so called hit list). Several computational methods exist to analyze and interpret such lists, the most widespread of which aim either at investigating of significantly enriched biological processes, or at extracting significantly represented subnetworks. RESULTS: Here we propose a novel network analysis method and corresponding computational software that employs the shortest path approach and centrality measure to discover members of molecular pathways leading to the studied phenotype, based on functional genomics screening data. The method works on integrated interactomes that consist of both directed and undirected networks - HIPPIE, SIGNOR, SignaLink, TFactS, KEGG, TransmiR, miRTarBase. The method finds nodes and short simple paths with significant high centrality in subnetworks induced by the hit genes and by so-called final implementers - the genes that are involved in molecular events responsible for final phenotypic realization of the biological processes of interest. We present the application of the method to the data from miRNA loss-of-function screen and transcriptome profiling of terminal human muscle differentiation process and to the gene loss-of-function screen exploring the genes that regulates human oxidative DNA damage recognition. The analysis highlighted the possible role of several known myogenesis regulatory miRNAs (miR-1, miR-125b, miR-216a) and their targets (AR, NR3C1, ARRB1, ITSN1, VAV3, TDGF1), as well as linked two major regulatory molecules of skeletal myogenesis, MYOD and SMAD3, to their previously known muscle-related targets (TGFB1, CDC42, CTCF) and also to a number of proteins such as C-KIT that have not been previously studied in the context of muscle differentiation. The analysis also showed the role of the interaction between H3 and SETDB1 proteins for oxidative DNA damage recognition. CONCLUSION: The current work provides a systematic methodology to discover members of molecular pathways in integrated networks using functional genomics screening data. It also offers a valuable instrument to explain the appearance of a set of genes, previously not associated with the process of interest, in the hit list of each particular functional genomics screening.

A Broad Set of Chromatin Factors Influences Splicing.

Several studies propose an influence of chromatin on pre-mRNA splicing, but it is still unclear how widespread and how direct this phenomenon is. We find here that when assembled in vivo, the U2 snRNP co-purifies with a subset of chromatin-proteins, including histones and remodeling complexes like SWI/SNF. Yet, an unbiased RNAi screen revealed that the outcome of splicing is influenced by a much larger variety of chromatin factors not all associating with the spliceosome. The availability of this broad range of chromatin factors impacting splicing further unveiled their very context specific effect, resulting in either inclusion or skipping, depending on the exon under scrutiny. Finally, a direct assessment of the impact of chromatin on splicing using an in vitro co-transcriptional splicing assay with pre-mRNAs transcribed from a nucleosomal template, demonstrated that chromatin impacts nascent pre-mRNP in their competence for splicing. Altogether, our data show that numerous chromatin factors associated or not with the spliceosome can affect the outcome of splicing, possibly as a function of the local chromatin environment that by default interferes with the efficiency of splicing.

A subset of the histone H3 lysine 9 methyltransferases Suv39h1, G9a, GLP, and SETDB1 participate in a multimeric complex.

Lysine 9 of histone 3 (H3K9) can be mono-, di-, or trimethylated, inducing distinct effects on gene expression and chromatin compaction. H3K9 methylation can be mediated by several histone methyltransferases (HKMTs) that possess mono-, di-, or trimethylation activities. Here we provide evidence that a subset of each of the main H3K9 HKMTs, G9a/KMT1C, GLP/KMT1D, SETDB1/KMT1E, and Suv39h1/KMT1A, coexist in the same megacomplex. Moreover, in Suv39h or G9a null cells, the remaining HKMTs are destabilized at the protein level, indicating that the integrity of these HKMTs is interdependent. The four HKMTs are recruited to major satellite repeats, a known Suv39h1 genomic target, but also to multiple G9a target genes. Moreover, we report a functional cooperation between the four H3K9 HKMTs in the regulation of known G9a target genes. Altogether, our data identify a H3K9 methylation multimeric complex.

Post-transcriptional modulation of interleukin 8 by CNOT6L regulates skeletal muscle differentiation.

CNOT6L is a deadenylase subunit belonging to the CCR4-NOT complex, a major deadenylase complex in eukaryotes involved at multiple levels in regulation of gene expression. While CNOT6L is expressed in skeletal muscle cells, its specific functions in this tissue are still largely unknown. Our previous work highlighted the functional of CNOT6L in skeletal muscle cell differentiation. To further explore how CNOT6L regulates myogenesis, we used here gene expression analysis to identify CNOT6L mRNA targets in human myoblasts. Among these novel targets, IL-8 (interleukin 8) mRNA was the most upregulated in CNOT6L knock-down (KD) cells. Biochemical approaches and poly (A) tail length assays showed that IL-8 mRNA is a direct target of CNOT6L, and further investigations by loss- and gain-of-function assays pointed out that IL-8 is an important effector of myogenesis. Therefore, we have characterized CNOT6L-IL-8 as a new signaling axis that regulates myogenesis.

miR-98 delays skeletal muscle differentiation by down-regulating E2F5.

A genome-wide screen had previously shown that knocking down miR-98 and let-7g, two miRNAs of the let-7 family, leads to a dramatic increase in terminal myogenic differentiation. In the present paper, we report that a transcriptomic analysis of human myoblasts, where miR-98 was knocked down, revealed that approximately 240 genes were sensitive to miR-98 depletion. Among these potential targets of miR-98, we identified the transcriptional repressor E2F5 and showed that it is a direct target of miR-98. Knocking down simultaneously E2F5 and miR-98 almost fully restored normal differentiation, indicating that E2F5 is involved in the regulation of skeletal muscle differentiation. We subsequently show that E2F5 can bind to the promoters of two inhibitors of terminal muscle differentiation, ID1 (inhibitor of DNA binding 1) and HMOX1 (heme oxygenase 1), which decreases their expression in skeletal myoblasts. We conclude that miR-98 regulates muscle differentiation by altering the expression of the transcription factor E2F5 and, in turn, of multiple E2F5 targets.

Kinetic signatures of microRNA modes of action.

MicroRNAs (miRNAs) are key regulators of all important biological processes, including development, differentiation, and cancer. Although remarkable progress has been made in deciphering the mechanisms used by miRNAs to regulate translation, many contradictory findings have been published that stimulate active debate in this field. Here we contribute to this discussion in three ways. First, based on a comprehensive analysis of the existing literature, we hypothesize a model in which all proposed mechanisms of microRNA action coexist, and where the apparent mechanism that is detected in a given experiment is determined by the relative values of the intrinsic characteristics of the target mRNAs and associated biological processes. Among several coexisting miRNA mechanisms, the one that will effectively be measurable is that which acts on or changes the sensitive parameters of the translation process. Second, we have created a mathematical model that combines nine known mechanisms of miRNA action and estimated the model parameters from the literature. Third, based on the mathematical modeling, we have developed a computational tool for discriminating among different possible individual mechanisms of miRNA action based on translation kinetics data that can be experimentally measured (kinetic signatures). To confirm the discriminatory power of these kinetic signatures and to test our hypothesis, we have performed several computational experiments with the model in which we simulated the coexistence of several miRNA action mechanisms in the context of variable parameter values of the translation.

The microRNA miR-181 targets the homeobox protein Hox-A11 during mammalian myoblast differentiation.

Deciphering the mechanisms underlying skeletal muscle-cell differentiation in mammals is an important challenge. Cell differentiation involves complex pathways regulated at both transcriptional and post-transcriptional levels. Recent observations have revealed the importance of small (20-25 base pair) non-coding RNAs (microRNAs or miRNAs) that are expressed in both lower organisms and in mammals. miRNAs modulate gene expression by affecting mRNA translation or stability. In lower organisms, miRNAs are essential for cell differentiation during development; some miRNAs are involved in maintenance of the differentiated state. Here, we show that miR-181, a microRNA that is strongly upregulated during differentiation, participates in establishing the muscle phenotype. Moreover, our results suggest that miR-181 downregulates the homeobox protein Hox-A11 (a repressor of the differentiation process), thus establishing a functional link between miR-181 and the complex process of mammalian skeletal-muscle differentiation. Therefore, miRNAs can be involved in the establishment of a differentiated phenotype - even when they are not expressed in the corresponding fully differentiated tissue.

Tandem affinity purification of miRNA target mRNAs (TAP-Tar).

MicroRNAs (miRNAs) bind to Argonaute proteins, and together they form the RISC complex and regulate target mRNA translation and/or stability. Identification of mRNA targets is key to deciphering the physiological functions and mode of action of miRNAs. In mammals, miRNAs are generally poorly homologous to their target sequence, and target identification cannot be based solely on bioinformatics. Here, we describe a biochemical approach, based on tandem affinity purification, in which mRNA/miRNA complexes are sequentially pulled down, first via the Argonaute moiety and then via the miRNA. Our 'TAP-Tar' procedure allows the specific pull down of mRNA targets of miRNA. It is useful for validation of targets predicted in silico, and, potentially, for discovery of previously uncharacterized targets.

Expression and cellular localization of microRNA-29b and RAX, an activator of the RNA-dependent protein kinase (PKR), in the retina of streptozotocin-induced diabetic rats.

PURPOSE: The apoptosis of retinal neurons plays a critical role in the pathogenesis of diabetic retinopathy (DR), but the molecular mechanisms underlying this phenomenon remain unclear. The purpose of this study was to investigate the cellular localization and the expression of microRNA-29b (miR-29b) and its potential target PKR associated protein X (RAX), an activator of the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway, in the retina of normal and diabetic rats. METHODS: Retinas were obtained from normal and diabetic rats within 35 days after streptozotocin (STZ) injection. In silico analysis indicated that RAX is a potential target of miR-29b. The cellular localization of miR-29b and RAX was assessed by in situ hybridization and immunofluorescence, respectively. The expression levels of miR-29b and RAX mRNA were evaluated by quantitative reverse transcription PCR (qRT-PCR), and the expression of RAX protein was evaluated by western blot. A luciferase reporter assay and inhibition of endogenous RAX were performed to confirm whether RAX is a direct target of miR-29b as predicted by the in silico analysis. RESULTS: We found that miR-29b and RAX are localized in the retinal ganglion cells (RGCs) and the cells of the inner nuclear layer (INL) of the retinas from normal and diabetic rats. Thus, the expression of miR-29b and RAX, as assessed in the retina by quantitative RT-PCR, reflects their expression in the RGCs and the cells of the INL. We also revealed that RAX protein is upregulated (more than twofold) at 3, 6, 16, and 22 days and downregulated (70%) at 35 days, whereas miR-29b is upregulated (more than threefold) at 28 and 35 days after STZ injection. We did not confirm the computational prediction that RAX is a direct target of miR-29b. CONCLUSIONS: Our results suggest that RAX expression may be indirectly regulated by miR-29b, and the upregulation of this miRNA at the early stage of STZ-induced diabetes may have a protective effect against the apoptosis of RGCs and cells of the INL by the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway.

Physical and functional interaction between heterochromatin protein 1alpha and the RNA-binding protein heterogeneous nuclear ribonucleoprotein U.

By combining biochemical purification and mass spectrometry, we identified proteins associated with human heterochromatin protein 1alpha (HP1alpha) both in the nucleoplasm and in chromatin. Some of these are RNA-binding proteins, and among them is the protein heterogeneous nuclear ribonucleoprotein U (hnRNP U)/SAF-A, which is linked to chromatin organization and transcriptional regulation. Here, we demonstrate that hnRNP U is a bona fide HP1alpha-interacting molecule. More importantly, hnRNP U depletion reduces HP1alpha-dependent gene silencing and disturbs HP1alpha subcellular localization. Thus, our data demonstrate that hnRNP U is involved in HP1alpha function, shedding new light on the mode of action of HP1alpha and on the function of hnRNP U.

[Small non-coding RNAs: new targets and new therapeutic possibilities]. / Les petits ARN non codant, nouvelles cibles et nouvelles approches thérapeutiques.

The discovery of regulatory small non-coding RNAs represents a revolution in our understanding of gene regulation. These small non-coding RNAs are powerful tools for exploring cellular pathways and for artificially controlling gene expression. Natural small RNAs also represent potential therapeutic targets.

Identification of microRNAs in skeletal muscle associated with lung cancer cachexia.

BACKGROUND: Cachexia, highly prevalent in patients with non-small cell lung cancer (NSCLC), impairs quality of life and is associated with reduced tolerance and responsiveness to cancer therapy and decreased survival. MicroRNAs (miRNAs) are small non-coding RNAs that play a central role in post-transcriptional gene regulation. Changes in intramuscular levels of miRNAs have been implicated in muscle wasting conditions. Here, we aimed to identify miRNAs that are differentially expressed in skeletal muscle of cachectic lung cancer patients to increase our understanding of cachexia and to allow us to probe their potential as therapeutic targets. METHODS: A total of 754 unique miRNAs were profiled and analysed in vastus lateralis muscle biopsies of newly diagnosed treatment-naïve NSCLC patients with cachexia (n = 8) and age-matched and sex-matched healthy controls (n = 8). miRNA expression analysis was performed using a TaqMan MicroRNA Array. In silico network analysis was performed on all significant differentially expressed miRNAs. Differential expression of the top-ranked miRNAs was confirmed using reverse transcription-quantitative real-time PCR in an extended group (n = 48) consisting of NSCLC patients with (n = 15) and without cachexia (n = 11) and healthy controls (n = 22). Finally, these miRNAs were subjected to univariate and multivariate Cox proportional hazard analysis using overall survival and treatment-induced toxicity data obtained during the follow-up of this group of patients. RESULTS: We identified 28 significant differentially expressed miRNAs, of which five miRNAs were up-regulated and 23 were down-regulated. In silico miRNA-target prediction analysis showed 158 functional gene targets, and pathway analysis identified 22 pathways related to the degenerative or regenerative processes of muscle tissue. Subsequently, the expression of six top-ranked miRNAs was measured in muscle biopsies of the entire patient group. Five miRNAs were detectable with reverse transcription-quantitative real-time PCR analysis, and their altered expression (expressed as fold change, FC) was confirmed in muscle of cachectic NSCLC patients compared with healthy control subjects: miR-424-5p (FC = 4.5), miR-424-3p (FC = 12), miR-450a-5p (FC = 8.6), miR-144-5p (FC = 0.59), and miR-451a (FC = 0.57). In non-cachectic NSCLC patients, only miR-424-3p was significantly increased (FC = 5.6) compared with control. Although the statistical support was not sufficient to imply these miRNAs as individual predictors of overall survival or treatment-induced toxicity, when combined in multivariate analysis, miR-450-5p and miR-451a resulted in a significant stratification between short-term and long-term survival. CONCLUSIONS: We identified differentially expressed miRNAs putatively involved in lung cancer cachexia. These findings call for further studies to investigate the causality of these miRNAs in muscle atrophy and the mechanisms underlying their differential expression in lung cancer cachexia.

Basic, simple and extendable kinetic model of protein synthesis.

Protein synthesis is one of the most fundamental biological processes. Despite existence of multiple mathematical models of translation, surprisingly, there is no basic and simple chemical kinetic model of this process, derived directly from the detailed kinetic scheme. One of the reasons for this is that the translation process is characterized by indefinite number of states, because of the structure of the polysome. We bypass this difficulty by applying lumping of multiple states of translated mRNA into few dynamical variables and by introducing a variable describing the pool of translating ribosomes. The simplest model can be solved analytically. The simplest model can be extended, if necessary, to take into account various phenomena such as the limited amount of ribosomal units or regulation of translation by microRNA. The introduced model is more suitable to describe the protein synthesis in eukaryotes but it can be extended to prokaryotes. The model can be used as a building block for more complex models of cellular processes. We demonstrate the utility of the model in two examples. First, we determine the critical parameters of the synthesis of a single protein for the case when the ribosomal units are abundant. Second, we demonstrate intrinsic bi-stability in the dynamics of the ribosomal protein turnover and predict that a minimal number of ribosomes should pre-exists in a living cell to sustain its protein synthesis machinery, even in the absence of proliferation.

Mathematical modeling reveals the factors involved in the phenomena of cancer stem cells stabilization.

Cancer Stem Cells (CSC), a subset of cancer cells resembling normal stem cells with self-renewal and asymmetric division capabilities, are present at various but low proportions in many tumors and are thought to be responsible for tumor relapses following conventional cancer therapies. In vitro, most intriguingly, isolated CSCs rapidly regenerate the original population of stem and non-stem cells (non-CSCs) as shown by various investigators. This phenomenon still remains to be explained. We propose a mathematical model of cancer cell population dynamics, based on the main parameters of cell population growth, including the proliferation rates, the rates of cell death and the frequency of symmetric and asymmetric cell divisions both in CSCs and non-CSCs sub-populations, and taking into account the stabilization phenomenon. The analysis of the model allows determination of time-varying corridors of probabilities for different cell fates, given the particular dynamics of cancer cells populations; and determination of a cell-cell communication factors influencing these time-varying probabilities of cell behavior (division, transition) scenarios. Though the results of the model have to be experimentally confirmed, we can anticipate the development of several fundamental and practical applications based on the theoretical results of the model.

The retinoblastoma protein binds the promoter of the survival gene bcl-2 and regulates its transcription in epithelial cells through transcription factor AP-2.

The retinoblastoma (RB) gene product has been shown to restrict cell proliferation, promote cell differentiation, and inhibit apoptosis. Loss of RB function can induce both p53-dependent apoptosis and p53-independent apoptosis; little is known about the mechanisms of RB-regulated p53-independent apoptosis. Here we show that RB specifically activates transcription of the survival gene bcl-2 in epithelial cells but not in NIH 3T3 mesenchymal cells. This transcriptional activity is mediated by the transcription factor AP-2. By monitoring protein-DNA interactions in living cells using formaldehyde cross-linking and chromatin immunoprecipitation, we show that endogenous RB and AP-2 both bind to the same bcl-2 promoter sequence. In addition, we demonstrate that RB and AP-2 also bind to the E-cadherin gene promoter in vivo, consistent with regulation of this promoter by both AP-2 and RB in epithelial cells. This study provides evidence that RB activates bcl-2 and E-cadherin by binding directly to the respective promoter sequences and not indirectly by repressing an inhibitor. This recruitment is mediated by a transcription factor, in this case AP-2. For the first time, our results suggest a direct molecular mechanism by which RB might inhibit apoptosis independently of p53. The results are discussed in a context where RB and Bcl-2 contribute under nonpathological conditions to the maintenance of cell viability in association with a differentiated phenotype, contributing to the tumor suppressor function of RB and playing important roles in normal development.

miR-600 Acts as a Bimodal Switch that Regulates Breast Cancer Stem Cell Fate through WNT Signaling.

Breast cancer stem cells (bCSCs) have been implicated in tumor progression and therapeutic resistance; however, the molecular mechanisms that define this state are unclear. We have performed two microRNA (miRNA) gain- and loss-of-function screens to identify miRNAs that regulate the choice between bCSC self-renewal and differentiation. We find that micro-RNA (miR)-600 silencing results in bCSC expansion, while its overexpression reduces bCSC self-renewal, leading to decreased in vivo tumorigenicity. miR-600 targets stearoyl desaturase 1 (SCD1), an enzyme required to produce active, lipid-modified WNT proteins. In the absence of miR-600, WNT signaling is active and promotes self-renewal, whereas overexpression of miR-600 inhibits the production of active WNT and promotes bCSC differentiation. In a series of 120 breast tumors, we found that a low level of miR-600 is correlated with active WNT signaling and a poor prognosis. These findings highlight a miR-600-centered signaling network that governs bCSC-fate decisions and influences tumor progression.

SiRNA-mediated inhibition of vascular endothelial growth factor severely limits tumor resistance to antiangiogenic thrombospondin-1 and slows tumor vascularization and growth.

In the past few years, several laboratories have developed antiangiogenic molecules that starve tumors by targeting their vasculature and we have shown that, when produced in tumors, the antiangiogenic molecule thrombospondin-1 (TSP1) reduces the vascularization and delays tumor onset. Yet over time, tumor cells producing active TSP1 do eventually form exponentially growing tumors. These tumors are composed of cells secreting unusually high amounts of the angiogenic stimulator vascular endothelial growth factor (VEGF) that are sufficient to overcome the inhibitory TSP1. Here, we use short double-stranded RNA (siRNA) to trigger RNA interference and thereby impair the synthesis of VEGF and ask if this inability to produce VEGF prevents the development of TSP1 resistance. Systemic in vivo administration of crude anti-VEGF siRNA reduced the growth of unaltered fibrosarcoma tumor cells, and when the anti-VEGF siRNA was expressed from tumor cells themselves, such inhibition was synergistic with the inhibitory effects derived from TSP1 secretion by the tumor cells. Anti-VEGF siRNA delayed the emergence of TSP1-resistant tumors and strikingly reduced their subsequent growth rate.

UV-dependent phosphorylation of COP9/signalosome in UV-induced apoptosis.

The COP9/signalosome (CSN) multi-protein complex regulates the activity of cullin-RING ubiquitin ligases (CRLs), including the DDB2 and CSA CRL4 ligases (CRL4DDB2 and CRL4CSA), which are involved in the repair of UV-induced DNA damages. In the present study, we demonstrated that the protein kinase ATM, a key component of the DNA damage response (DDR), phosphorylates CSN1 and CSN7a, two subunits of the CSN complex, in a UV-dependent manner. The phosphorylation of CSN1 on serine 474 was detected as early as 3 h after UV-exposure, peaked at 8 h and persisted until 48 h post-UV irradiation. Such a time course suggests a role in late DDR rather than in DNA repair. Consistently, overexpression of a phosphorylation-resistant S474A CSN1 mutant reduced UV-induced apoptosis. Thus, CSN1 appears to play a role not only in DNA repair but also in UV-induced apoptosis.

The WEE1 regulators CPEB1 and miR-15b switch from inhibitor to activators at G2/M.

MicroRNAs (miRNAs) in the AGO-containing RISC complex control messenger RNA (mRNA) translation by binding to mRNA 3' untranslated region (3'UTR). The relationship between miRNAs and other regulatory factors that also bind to mRNA 3'UTR, such as CPEB1 (cytoplasmic polyadenylation element-binding protein), remains elusive. We found that both CPEB1 and miR-15b control the expression of WEE1, a key mammalian cell cycle regulator. Together, they repress WEE1 protein expression during G1 and S-phase. Interestingly, the 2 factors lose their inhibitory activity at the G2/M transition, at the time of the cell cycle when WEE1 expression is maximal, and, moreover, rather activate WEE1 translation in a synergistic manner. Our data show that translational regulation by RISC and CPEB1 is essential in cell cycle control and, most importantly, is coordinated, and can be switched from inhibition to activation during the cell cycle.

HDAC4 mediates transcriptional repression by the acute promyelocytic leukaemia-associated protein PLZF.

PLZF, the promyelocytic leukaemia zinc-finger protein, is a transcriptional repressor essential to development. In some acute leukaemias, a chromosomal translocation fusing the PLZF gene to that encoding the retinoic acid receptor RARalpha gives rise to a fusion protein, PLZF-RARalpha, thought to be responsible for constitutive repression of differentiation-associated genes in these cells. Repression by both PLZF and PLZF-RARalpha is sensitive to the histone deacetylase inhibitor TSA, and PLZF was previously shown to interact physically with HDAC1, a class I histone deacetylase. We here asked whether class II histone deacetylases, known to be generally involved in differentiation processes, participate in the repression mediated by PLZF and PLZF-RARalpha, and found that PLZF interacts with HDAC4 in both GST-pull-down and co-immunoprecipitation assays. Furthermore, HDAC4 is indeed involved in PLZF and PLZF-RARalpha-mediated repression, since an enzymatically dead mutant of HDAC4 released the repression, as did an siRNA that blocks HDAC4 expression. Taken together, our data indicate that recruitment of HDAC4 is necessary for PLZF-mediated repression in both normal and leukaemic cells.