Transcriptional networks controlling the cell cycle

Expression and Functional Analysis of Immuno-Micro-RNAs mir-146a and mir-326 in Colorectal Cancer

Micro-RNAs (miRNAs) are non-coding RNAs with importance in the development of cancer. They are involved in both tumor development and immune processes in tumors. The present study aims to characterize the behavior of two miRNAs, the proinflammatory miR-326-5p and the anti-inflammatory miR-146a-5p, in colorectal cancer (CRC), to decipher the mechanisms that regulate their expression, and to study potential applications. Tissue levels of miR-326-5p and miR-146a-5p were determined by qrt-PCR (real-time quantitative reverse transcription polymerase chain reaction) in 45 patients with colorectal cancer in tumoral and normal adjacent tissue. Subsequent bioinformatic analysis was performed to characterize the transcriptional networks that control the expression of the two miRNAs. The biomarker potential of miRNAs was assessed. The expression of miR-325-5p and miR-146a-5p was decreased in tumors compared to normal tissue. The two miRNAs are regulated through a transcriptional network, which originates in the inflammatory and proliferative pathways and regulates a set of cellular functions related to immunity, proliferation, and differentiation. The miRNAs coordinate distinct modules in the network. There is good biomarker potential of miR-326 with an AUC (Area under the curve) of 0.827, 0.911 sensitivity (Sn), and 0.689 specificity (Sp), and of the combination miR-326-miR-146a, with an AUC of 0.845, Sn of 0.75, and Sp of 0.89. The miRNAs are downregulated in the tumor tissue. They are regulated by a transcriptional network in which they coordinate distinct modules. The structure of the network highlights possible therapeutic approaches. MiR-326 and the combination of the two miRNAs may serve as biomarkers in CRC.

A Group of Tumor-Suppressive micro-RNAs Changes Expression Coordinately in Colon Cancer

MicroRNAs (miRNAs) are molecules with a role in the post-transcriptional regulation of messenger RNA, being involved in a wide range of biological and pathological processes. In the present study, we aim to characterize the behavior of a few miRNAs with roles in the cell cycle and differentiation of colon cancer (CC) cells. The present work considers miRNAs as reflections of the complex cellular processes in which they are generated, their observed variations being used to characterize the molecular networks in which they are part and through which cell proliferation is achieved. Tumoral and adjacent normal tissue samples were obtained from 40 CC patients, and the expression of miR-29a, miR-146a, miR-215 and miR-449 were determined by qRT-PCR analysis. Subsequent bioinformatic analysis was performed to highlight the transcription factors (TFs) network that regulate the miRNAs and functionally characterizes this network. There was a significant decrease in the expression of all miRNAs in tumor tissue. All miRNAs were positively correlated with each other. The analysis of the TF network showed tightly connected functional modules related to the cell cycle and associated processes. The four miRNAs are downregulated in CC; they are strongly correlated, showing coherence within the cellular network that regulates them and highlighting possible approach strategies.

The Green Valley of Drosophila melanogaster Constitutive Heterochromatin: Protein-Coding Genes Involved in Cell Division Control

Constitutive heterochromatin represents a significant fraction of eukaryotic genomes (10% in Arabidopsis, 20% in humans, 30% in D. melanogaster, and up to 85% in certain nematodes) and shares similar genetic and molecular properties in animal and plant species. Studies conducted over the last few years on D. melanogaster and other organisms led to the discovery of several functions associated with constitutive heterochromatin. This made it possible to revise the concept that this ubiquitous genomic territory is incompatible with gene expression. The aim of this review is to focus the attention on a group of protein-coding genes resident in D. melanogaster constitutive of heterochromatin, which are implicated in different steps of cell division.

MYC drives autophagy to adapt to stress in Penaeus vannamei

MYC proto-oncogene (MYC), a first oncogenic nuclear transcription factor isolated from the human genome, belongs to the helix loop helix/leucine zipper protein family (bHLHzip). MYC plays an important part in the process of various physiological and biochemical of vertebrate, such as cell growth, proliferation, cycle, and autophagy. However, its molecular regulation mechanism and function in invertebrates are still unclear. In this study, a novel transcription factor MYC gene was screened, cloned, and characterized from Penaeus vannamei. The open reading frame of PvMYC was 1593bp, encode a polypeptide of 530 amino acids with molecular weight of 58.5 kDa, and a theoretical PI of 5.75. The results of tissue distribution showed that PvMYC was constitutively expressed in all detected tissues, and highest expression in hepatopancreas. The expression level of PvMYC up-regulated significantly and responded to low temperature stress by nuclear ectopic after low temperature stress. Overexpression of PvMYC in shrimp hemocytes negatively regulated the expression of Beclin-1 and reduced the conversion from LC3I to LC3II, yet p62 was decreased significantly. Meanwhile, RAPA eliminated the inhibition of autophagy caused by overexpression of PvMYC. ROS levels and autophagy flux showed the similar trend under low temperature stress after silencing PvMYC. The expression levels of Beclin-1, key ATG gene and LC3II increased significantly, while p62 decreased significantly under the same conditions. In addition, the Total hemocyte count (THC) decreased sharply, and accelerated the injury of hepatopancreas under low temperature stress after silencing PvMYC. Collectively, these results suggest that PvMYC has vital role in the cold adaptation mechanism of P. vannamei by negatively regulating autophagy.

Upregulation of ribosome biogenesis via canonical E-boxes is required for Myc-driven proliferation

The transcription factor Myc drives cell growth across animal phyla and is activated in most forms of human cancer. However, it is unclear which Myc target genes need to be regulated to induce growth and whether multiple targets act additively or if induction of each target is individually necessary. Here, we identified Myc target genes whose regulation is conserved between humans and flies and deleted Myc-binding sites (E-boxes) in the promoters of fourteen of these genes in Drosophila. E-box mutants of essential genes were homozygous viable, indicating that the E-boxes are not required for basal expression. Eight E-box mutations led to Myc-like phenotypes; the strongest mutant, ppan^Ebox-/-, also made the flies resistant to Myc-induced cell growth without affecting Myc-induced apoptosis. The ppan^Ebox-/- flies are healthy and display only a minor developmental delay, suggesting that it may be possible to treat or prevent tumorigenesis by targeting individual downstream targets of Myc.

Cell size homeostasis: Metabolic control of growth and cell division

Joint regulation of growth rate and cell division rate determines cell size. Here we discuss how animal cells achieve cell size homeostasis potentially involving multiple signaling pathways converging at metabolic regulation of growth rate and cell cycle progression. While several models have been developed to explain cell size control, comparison of the two predominant models shows that size homeostasis is dependent on the ability to adjust cellular growth rate based on cell size. Consequently, maintenance of size homeostasis requires that larger cells can grow slower than small cells in relative terms. We review recent experimental evidence showing that such size adjustment occurs primarily at or immediately before the G1/S transition of the cell cycle. We further propose that bidirectional feedback between growth rate and size results in cell size sensing and discuss potential mechanisms how this may be accomplished.

PAF1 complex component Leo1 helps recruit Drosophila Myc to promoters

The Myc oncogene is a transcription factor with a powerful grip on cellular growth and proliferation. The physical interaction of Myc with the E-box DNA motif has been extensively characterized, but it is less clear whether this sequence-specific interaction is sufficient for Myc's binding to its transcriptional targets. Here we identify the PAF1 complex, and specifically its component Leo1, as a factor that helps recruit Myc to target genes. Since the PAF1 complex is typically associated with active genes, this interaction with Leo1 contributes to Myc targeting to open promoters.

NetProphet 2.0: mapping transcription factor networks by exploiting scalable data resources

Motivation

Cells process information, in part, through transcription factor (TF) networks, which control the rates at which individual genes produce their products. A TF network map is a graph that indicates which TFs bind and directly regulate each gene. Previous work has described network mapping algorithms that rely exclusively on gene expression data and ‘integrative’ algorithms that exploit a wide range of data sources including chromatin immunoprecipitation sequencing (ChIP-seq) of many TFs, genome-wide chromatin marks, and binding specificities for many TFs determined in vitro. However, such resources are available only for a few major model systems and cannot be easily replicated for new organisms or cell types.

Results

We present NetProphet 2.0, a ‘data light’ algorithm for TF network mapping, and show that it is more accurate at identifying direct targets of TFs than other, similarly data light algorithms. In particular, it improves on the accuracy of NetProphet 1.0, which used only gene expression data, by exploiting three principles. First, combining multiple approaches to network mapping from expression data can improve accuracy relative to the constituent approaches. Second, TFs with similar DNA binding domains bind similar sets of target genes. Third, even a noisy, preliminary network map can be used to infer DNA binding specificities from promoter sequences and these inferred specificities can be used to further improve the accuracy of the network map.

Availability and implementation

Source code and comprehensive documentation are freely available at https://github.com/yiming-kang/NetProphet_2.0.

Supplementary information

Supplementary data are available at Bioinformatics online.

SlgA, encoded by the homolog of the human schizophrenia-associated gene PRODH, acts in clock neurons to regulate Drosophila aggression

Mutations in the proline dehydrogenase gene PRODH are linked to behavioral alterations in schizophrenia and as part of DiGeorge and velo-cardio-facial syndromes, but the role of PRODH in their etiology remains unclear. Here, we establish a Drosophila model to study the role of PRODH in behavioral disorders. We determine the distribution of the Drosophila PRODH homolog slgA in the brain and show that knockdown and overexpression of human PRODH and slgA in the lateral neurons ventral (LNv) lead to altered aggressive behavior. SlgA acts in an isoform-specific manner and is regulated by casein kinase II (CkII). Our data suggest that these effects are, at least partially, due to effects on mitochondrial function. We thus show that precise regulation of proline metabolism is essential to drive normal behavior and we identify Drosophila aggression as a model behavior relevant for the study of the mechanisms that are impaired in neuropsychiatric disorders.

Editors' choice: A Drosophila model to study the role of PRODH, a schizophrenia-associated gene, in behavioral disorders.

Transcriptional reprogramming in yeast using dCas9 and combinatorial gRNA strategies

BACKGROUND

Transcriptional reprogramming is a fundamental process of living cells in order to adapt to environmental and endogenous cues. In order to allow flexible and timely control over gene expression without the interference of native gene expression machinery, a large number of studies have focused on developing synthetic biology tools for orthogonal control of transcription. Most recently, the nuclease-deficient Cas9 (dCas9) has emerged as a flexible tool for controlling activation and repression of target genes, by the simple RNA-guided positioning of dCas9 in the vicinity of the target gene transcription start site.

RESULTS

In this study we compared two different systems of dCas9-mediated transcriptional reprogramming, and applied them to genes controlling two biosynthetic pathways for biobased production of isoprenoids and triacylglycerols (TAGs) in baker's yeast Saccharomyces cerevisiae. By testing 101 guide-RNA (gRNA) structures on a total of 14 different yeast promoters, we identified the best-performing combinations based on reporter assays. Though a larger number of gRNA-promoter combinations do not perturb gene expression, some gRNAs support expression perturbations up to ~threefold. The best-performing gRNAs were used for single and multiplex reprogramming strategies for redirecting flux related to isoprenoid production and optimization of TAG profiles. From these studies, we identified both constitutive and inducible multiplex reprogramming strategies enabling significant changes in isoprenoid production and increases in TAG.

CONCLUSION

Taken together, we show similar performance for a constitutive and an inducible dCas9 approach, and identify multiplex gRNA designs that can significantly perturb isoprenoid production and TAG profiles in yeast without editing the genomic context of the target genes. We also identify a large number of gRNA positions in 14 native yeast target pomoters that do not affect expression, suggesting the need for further optimization of gRNA design tools and dCas9 engineering.

Past Roadblocks and New Opportunities in Transcription Factor Network Mapping

One of the principal mechanisms by which cells differentiate and respond to changes in external signals or conditions is by changing the activity levels of transcription factors (TFs). This changes the transcription rates of target genes via the cell's TF network, which ultimately contributes to reconfiguring cellular state. Since microarrays provided our first window into global cellular state, computational biologists have eagerly attacked the problem of mapping TF networks, a key part of the cell's control circuitry. In retrospect, however, steady-state mRNA abundance levels were a poor substitute for TF activity levels and gene transcription rates. Likewise, mapping TF binding through chromatin immunoprecipitation proved less predictive of functional regulation and less amenable to systematic elucidation of complete networks than originally hoped. This review explains these roadblocks and the current, unprecedented blossoming of new experimental techniques built on second-generation sequencing, which hold out the promise of rapid progress in TF network mapping.

Systems proteomics of liver mitochondria function

Recent improvements in quantitative proteomics approaches, including Sequential Window Acquisition of all Theoretical Mass Spectra (SWATH-MS), permit reproducible large-scale protein measurements across diverse cohorts. Together with genomics, transcriptomics, and other technologies, transomic data sets can be generated that permit detailed analyses across broad molecular interaction networks. Here, we examine mitochondrial links to liver metabolism through the genome, transcriptome, proteome, and metabolome of 386 individuals in the BXD mouse reference population. Several links were validated between genetic variants toward transcripts, proteins, metabolites, and phenotypes. Among these, sequence variants in Cox7a2l alter its protein's activity, which in turn leads to downstream differences in mitochondrial supercomplex formation. This data set demonstrates that the proteome can now be quantified comprehensively, serving as a key complement to transcriptomics, genomics, and metabolomics--a combination moving us forward in complex trait analysis.

Inferring gene targets of drugs and chemical compounds from gene expression profiles

Motivation: Finding genes which are directly perturbed or targeted by drugs is of great interest and importance in drug discovery. Several network filtering methods have been created to predict the gene targets of drugs from gene expression data based on an ordinary differential equation model of the gene regulatory network (GRN). A critical step in these methods involves inferring the GRN from the expression data, which is a very challenging problem on its own. In addition, existing network filtering methods require computationally intensive parameter tuning or expression data from experiments with known genetic perturbations or both.

Results: We developed a method called DeltaNet for the identification of drug targets from gene expression data. Here, the gene target predictions were directly inferred from the data without a separate step of GRN inference. DeltaNet formulation led to solving an underdetermined linear regression problem, for which we employed least angle regression (DeltaNet-LAR) or LASSO regularization (DeltaNet-LASSO). The predictions using DeltaNet for expression data of Escherichia coli, yeast, fruit fly and human were significantly more accurate than those using network filtering methods, namely mode of action by network identification (MNI) and sparse simultaneous equation model (SSEM). Furthermore, DeltaNet using LAR did not require any parameter tuning and could provide computational speed-up over existing methods.

Conclusion: DeltaNet is a robust and numerically efficient tool for identifying gene perturbations from gene expression data. Importantly, the method requires little to no expert supervision, while providing accurate gene target predictions.

Availability and implementation: DeltaNet is available on http://www.cabsel.ethz.ch/tools/DeltaNet.

Contact:rudi.gunawan@chem.ethz.ch

Supplementary information:Supplementary data are available at Bioinformatics online.

snoRNAs are a novel class of biologically relevant Myc targets

BACKGROUND

Myc proteins are essential regulators of animal growth during normal development, and their deregulation is one of the main driving factors of human malignancies. They function as transcription factors that (in vertebrates) control many growth- and proliferation-associated genes, and in some contexts contribute to global gene regulation.

RESULTS

We combine chromatin immunoprecipitation-sequencing (ChIPseq) and RNAseq approaches in Drosophila tissue culture cells to identify a core set of less than 500 Myc target genes, whose salient function resides in the control of ribosome biogenesis. Among these genes we find the non-coding snoRNA genes as a large novel class of Myc targets. All assayed snoRNAs are affected by Myc, and many of them are subject to direct transcriptional activation by Myc, both in Drosophila and in vertebrates. The loss of snoRNAs impairs growth during normal development, whereas their overexpression increases tumor mass in a model for neuronal tumors.

CONCLUSIONS

This work shows that Myc acts as a master regulator of snoRNP biogenesis. In addition, in combination with recent observations of snoRNA involvement in human cancer, it raises the possibility that Myc's transforming effects are partially mediated by this class of non-coding transcripts.

Conservation of transcription factor binding specificities across 600 million years of bilateria evolution

Divergent morphology of species has largely been ascribed to genetic differences in the tissue-specific expression of proteins, which could be achieved by divergence in cis-regulatory elements or by altering the binding specificity of transcription factors (TFs). The relative importance of the latter has been difficult to assess, as previous systematic analyses of TF binding specificity have been performed using different methods in different species. To address this, we determined the binding specificities of 242 Drosophila TFs, and compared them to human and mouse data. This analysis revealed that TF binding specificities are highly conserved between Drosophila and mammals, and that for orthologous TFs, the similarity extends even to the level of very subtle dinucleotide binding preferences. The few human TFs with divergent specificities function in cell types not found in fruit flies, suggesting that evolution of TF specificities contributes to emergence of novel types of differentiated cells.

Spatial analysis of gene regulation reveals new insights into the molecular basis of upper thermal limits

The cellular stress response has long been the primary model for studying the molecular basis of thermal adaptation, yet the link between gene expression, RNA metabolism and physiological responses to thermal stress remains largely unexplored. We address this by comparing the transcriptional and physiological responses of three geographically distinct populations of Drosophila melanogaster from eastern Australia in response to, and recovery from, a severe heat stress with and without a prestress hardening treatment. We focus on starvin (stv), recently identified as an important thermally responsive gene. Intriguingly, stv encodes seven transcripts from alternative transcription sites and alternative splicing, yet appears to be rapidly heat inducible. First, we show genetic differences in upper thermal limits of the populations tested. We then demonstrate that the stv locus does not ubiquitously respond to thermal stress but is expressed as three distinct thermal and temporal RNA phenotypes (isoforms). The shorter transcript isoforms are rapidly upregulated under stress in all populations and show similar molecular signatures to heat-shock proteins. Multiple stress exposures seem to generate a reserve of pre-mRNAs, effectively 'priming' the cells for subsequent stress. Remarkably, we demonstrate a bypass in the splicing blockade in these isoforms, suggesting an essential role for these transcripts under heat stress. Temporal profiles for the weakly heat responsive stv isoform subset show opposing patterns in the two most divergent populations. Innate and induced transcriptome responses to hyperthermia are complex, and warrant moving beyond gene-level analyses.

Cyclin-dependent kinase 8 module expression profiling reveals requirement of mediator subunits 12 and 13 for transcription of Serpent-dependent innate immunity genes in Drosophila

The Cdk8 (cyclin-dependent kinase 8) module of Mediator integrates regulatory cues from transcription factors to RNA polymerase II. It consists of four subunits where Med12 and Med13 link Cdk8 and cyclin C (CycC) to core Mediator. Here we have investigated the contributions of the Cdk8 module subunits to transcriptional regulation using RNA interference in Drosophila cells. Genome-wide expression profiling demonstrated separation of Cdk8-CycC and Med12-Med13 profiles. However, transcriptional regulation by Cdk8-CycC was dependent on Med12-Med13. This observation also revealed that Cdk8-CycC and Med12-Med13 often have opposite transcriptional effects. Interestingly, Med12 and Med13 profiles overlapped significantly with that of the GATA factor Serpent. Accordingly, mutational analyses indicated that GATA sites are required for Med12-Med13 regulation of Serpent-dependent genes. Med12 and Med13 were also found to be required for Serpent-activated innate immunity genes in defense to bacterial infection. The results reveal a novel role for the Cdk8 module in Serpent-dependent transcription and innate immunity.

Drosophila mbm is a nucleolar myc and casein kinase 2 target required for ribosome biogenesis and cell growth of central brain neuroblasts

Proper cell growth is a prerequisite for maintaining repeated cell divisions. Cells need to translate information about intracellular nutrient availability and growth cues from energy-sensing organs into growth-promoting processes, such as sufficient supply with ribosomes for protein synthesis. Mutations in the mushroom body miniature (mbm) gene impair proliferation of neural progenitor cells (neuroblasts) in the central brain of Drosophila melanogaster. Yet the molecular function of Mbm has so far been unknown. Here we show that mbm does not affect the molecular machinery controlling asymmetric cell division of neuroblasts but instead decreases their cell size. Mbm is a nucleolar protein required for small ribosomal subunit biogenesis in neuroblasts. Accordingly, levels of protein synthesis are reduced in mbm neuroblasts. Mbm expression is transcriptionally regulated by Myc, which, among other functions, relays information from nutrient-dependent signaling pathways to ribosomal gene expression. At the posttranslational level, Mbm becomes phosphorylated by casein kinase 2 (CK2), which has an impact on localization of the protein. We conclude that Mbm is a new part of the Myc target network involved in ribosome biogenesis, which, together with CK2-mediated signals, enables neuroblasts to synthesize sufficient amounts of proteins required for proper cell growth.

Myc function in Drosophila

Drosophila contains a single MYC gene. Like its vertebrate homologs, it encodes a transcription factor that activates many targets, including prominently genes involved in ribosome biogenesis and translation. This activity makes Myc a central regulator of growth and/or proliferation of many cell types, such as imaginal disc cells, polyploid cells, stem cells, and blood cells. Importantly, not only does Myc act cell autonomously but it also affects the fate of adjacent cells and tissues. This potential of Myc is harnessed by many different signaling pathways, involving, among others, Wg, Dpp, Hpo, ecdysone, insulin, and mTOR.