Identification of cell cycle-regulated genes periodically expressed in U2OS cells and their regulation by FOXM1 and E2F transcription factors

Towards understanding the Regulation of Histone H1 Somatic Subtypes with OMICs

Histone H1 is involved in the regulation of chromatin higher-order structure and compaction. In humans, histone H1 is a multigene family with seven subtypes differentially expressed in somatic cells. Which are the regulatory mechanisms that determine the variability of the H1 complement is a long-standing biological question regarding histone H1. We have used a new approach based on the integration of OMICs data to address this issue. We have examined the 3D-chromatin structure, the binding of transcription factors (TFs), and the expression of somatic H1 genes in human cell lines, using data from public repositories, such as ENCODE. Analysis of Hi-C, ChIP-seq, and RNA-seq data, have revealed that transcriptional control has a greater impact on H1 regulation than previously thought. Somatic H1 genes located in topologically associated domains (TADs) show higher expression than in boundary regions. H1 genes are targeted by a variable number of transcription factors including cell cycle-related TFs, and tissue-specific TFs, suggesting a fine-tuned, subtype-specific transcriptional control. We describe, for the first time, that all H1 somatic subtypes are under transcriptional co-regulation. The replication-independent subtypes, which are encoded in different chromosomes isolated from other histone genes, are also co-regulated with the rest of the somatic H1 genes, indicating that transcriptional co-regulation extends beyond the histone cluster. Transcriptional control and transcriptional co-regulation explain, at least in part, the variability of H1 complement, the fluctuations of H1 subtypes during development, and also the compensatory effects observed, in model systems, after perturbation of one or more H1 subtypes.

In silico APC/C substrate discovery reveals cell cycle-dependent degradation of UHRF1 and other chromatin regulators

The anaphase-promoting complex/cyclosome (APC/C) is an E3 ubiquitin ligase and critical regulator of cell cycle progression. Despite its vital role, it has remained challenging to globally map APC/C substrates. By combining orthogonal features of known substrates, we predicted APC/C substrates in silico. This analysis identified many known substrates and suggested numerous candidates. Unexpectedly, chromatin regulatory proteins are enriched among putative substrates, and we show experimentally that several chromatin proteins bind APC/C, oscillate during the cell cycle, and are degraded following APC/C activation, consistent with being direct APC/C substrates. Additional analysis revealed detailed mechanisms of ubiquitylation for UHRF1, a key chromatin regulator involved in histone ubiquitylation and DNA methylation maintenance. Disrupting UHRF1 degradation at mitotic exit accelerates G1-phase cell cycle progression and perturbs global DNA methylation patterning in the genome. We conclude that APC/C coordinates crosstalk between cell cycle and chromatin regulatory proteins. This has potential consequences in normal cell physiology, where the chromatin environment changes depending on proliferative state, as well as in disease.

The forkhead box M1 (FOXM1) expression and antitumor effect of FOXM1 inhibition in malignant rhabdoid tumor

PURPOSE

Malignant rhabdoid tumor (MRT) is a rare, highly aggressive sarcoma with an uncertain cell of origin. Despite the existing standard of intensive multimodal therapy, the prognosis of patients with MRT is very poor. Novel antitumor agents are needed for MRT patients. Forkhead box transcription factor 1 (FOXM1) is overexpressed and is correlated with the pathogenesis in several human malignancies. In this study, we identified the clinicopathological and prognostic values of the expression of FOXM1 and its roles in the progression of MRT.

METHODS

We investigated the FOXM1 expression levels and their clinical significance in 23 MRT specimens using immunohistochemistry and performed clinicopathologic and prognostic analyses. We also demonstrated correlations between the downregulation of FOXM1 and oncological characteristics using small interfering RNA (siRNA) and FOXM1 inhibitor in MRT cell lines.

RESULTS

Histopathological analyses revealed that primary renal MRTs showed significantly low FOXM1 protein expression levels (p = 0.032); however, there were no significant differences in other clinicopathological characteristics or the survival rate. FOXM1 siRNA and FOXM1 inhibitor (thiostrepton) successfully downregulated the mRNA and protein expression of FOXM1 in vitro and the downregulation of FOXM1 inhibited cell proliferation, drug resistance to chemotherapeutic agents, migration, invasion, and caused the cell cycle arrest and apoptosis of MRT cell lines. A cDNA microarray analysis showed that FOXM1 regulated FANCD2 and NBS1, which are key genes for DNA damage repair.

CONCLUSION

This study demonstrates that FOXM1 may serve as a promising therapeutic target for MRT.

Pluripotency-associated genes reposition during early embryonic developmental stages in pigs

We examined the allelic expression and positioning of two pluripotency‐associated genes, OCT4 and SOX2, and two housekeeping genes, ACTB and TUBA, in 4‐ and 8‐cell porcine embryos utilizing RNA and DNA fluorescence in situ hybridization (FISH) in single blastomeres. The proportion of blastomeres expressing SOX2 bi‐allelically increased from 45% at the 4‐cell stage to 60% at the 8‐cell stage. Moreover, in 8‐cell embryos, SOX2 was expressed bi‐allelically in significantly more blastomeres than was the case for OCT4, and this was associated with a tendency for SOX2 alleles to move toward the nuclear interior during 4‐ to 8‐cell transition. However, the radial location of OCT4 alleles did not change significantly during this transition. The locations of active and inactive alleles based on DNA and RNA FISH signals were also calculated. Inactive OCT4 alleles were located in very close proximity to the nuclear membrane, whereas active OCT4 alleles were more centrally disposed in the nucleus. Nevertheless, the nuclear location of active and inactive SOX2 alleles did not change in either 4‐ or 8‐cell blastomeres. Our RNA and DNA FISH data provide novel information on the allelic expression patterns and positioning of pluripotency‐associated genes, OCT4 and SOX2, during embryonic genome activation in pigs.

Mapping metabolic oscillations during cell cycle progression

Proliferating cells must synthesize a wide variety of macromolecules while progressing through the cell cycle, but the coordination between cell cycle progression and cellular metabolism is still poorly understood. To identify metabolic processes that oscillate over the cell cycle, we performed comprehensive, non-targeted liquid chromatography-high resolution mass spectrometry (LC-HRMS) based metabolomics of HeLa cells isolated in the G₁ and SG₂M cell cycle phases, capturing thousands of diverse metabolite ions. When accounting for increased total metabolite abundance due to cell growth throughout the cell cycle, 18% of the observed LC-HRMS peaks were at least twofold different between the stages, consistent with broad metabolic remodeling throughout the cell cycle. While most amino acids, phospholipids, and total ribonucleotides were constant across cell cycle phases, consistent with the view that total macromolecule synthesis does not vary across the cell cycle, certain metabolites were oscillating. For example, ribonucleotides were highly phosphorylated in SG₂M, indicating an increase in energy charge, and several phosphatidylinositols were more abundant in G₁, possibly indicating altered membrane lipid signaling. Within carbohydrate metabolism, pentose phosphates and methylglyoxal metabolites were associated with the cycle. Interestingly, hundreds of yet uncharacterized metabolites similarly oscillated between cell cycle phases, suggesting previously unknown metabolic activities that may be synchronized with cell cycle progression, providing an important resource for future studies.

FUBP1 and FUBP2 enforce distinct epigenetic setpoints for MYC expression in primary single murine cells

Physiologically, MYC levels must be precisely set to faithfully amplify the transcriptome, but in cancer MYC is quantitatively misregulated. Here, we study the variation of MYC amongst single primary cells (B-cells and murine embryonic fibroblasts, MEFs) for the repercussions of variable cellular MYC-levels and setpoints. Because FUBPs have been proposed to be molecular “cruise controls” that constrain MYC expression, their role in determining basal or activated MYC-levels was also examined. Growing cells remember low and high-MYC setpoints through multiple cell divisions and are limited by the same expression ceiling even after modest MYC-activation. High MYC MEFs are enriched for mRNAs regulating inflammation and immunity. After strong stimulation, many cells break through the ceiling and intensify MYC expression. Lacking FUBPs, unstimulated MEFs express levels otherwise attained only with stimulation and sponsor MYC chromatin changes, revealed by chromatin marks. Thus, the FUBPs enforce epigenetic setpoints that restrict MYC expression.

Ying Zheng et al. characterize MYC gene and protein expression in single mammalian cells in response to various external signals. They find that individual cells show either high or low basal MYC expression setpoints, and that adherence to these setpoints as well as the magnitude of the response of MYC to stimulation, is controlled by FUBP1 and FUBP2.

A database of integrated molecular and phytochemical interactions of the foxm1 pathway for lung cancer

The FoxM1 pathway is an oncogenic signaling pathway involved in essential mechanisms including control cell-cycle progression, apoptosis and cell growth which are the common hallmarks of various cancers. Although its biological functions in the tumor development and progression are known, the mechanism by which it participates in those processes is not understood. The present work reveals images of the oncogenic FoxM1 pathway controlling the cell cycle process with alternative treatment options via phytochemical substances in the lung cancer study. The downstream significant protein modules of the FoxM1 pathway were extracted by the Molecular Complex Detection (MCODE) and the maximal clique (Mclique) algorithms. Furthermore, the effects of post-transcriptional modification by microRNA, transcription factor binding and the phytochemical compounds are observed through their interactions with the lung cancer protein modules. We provided two case studies to demonstrate the usefulness of our database. Our results suggested that the combination of various phytochemicals is effective in the treatment of lung cancer. The ultimate goal of the present work is to partly support the discovery of plant-derived compounds in combination treatment of classical chemotherapeutic agents to increase the efficacy of lung cancer method probably with minor side effects. Furthermore, a web-based system displaying results of the present work is set up for investigators posing queries at http://sit.mfu.ac.th/lcgdb/index_FoxM1.php.Communicated by Ramaswamy H. Sarma.

Complex Cartography: Regulation of E2F Transcription Factors by Cyclin F and Ubiquitin

The E2F family of transcriptional regulators sits at the center of cell cycle gene expression and plays vital roles in normal and cancer cell cycles. Whereas control of E2Fs by the retinoblastoma family of proteins is well established, much less is known about their regulation by ubiquitin pathways. Recent studies placed the Skp1-Cul1-F-box-protein (SCF) family of E3 ubiquitin ligases with the F-box protein Cyclin F at the center of E2F regulation, demonstrating temporal proteolysis of both activator and atypical repressor E2Fs. Importantly, these E2F members, in particular activator E2F1 and repressors E2F7 and E2F8, form a feedback circuit at the crossroads of cell cycle and cell death. Moreover, Cyclin F functions in a reciprocal circuit with the cell cycle E3 ligase anaphase-promoting complex/cyclosome (APC/C), which also controls E2F7 and E2F8. This review focuses on the complex contours of feedback within this circuit, highlighting the deep crosstalk between E2F, SCF-Cyclin F, and APC/C in regulating the oscillator underlying human cell cycles.

Head-to-head comparisons of Toxoplasma gondii and its near relative Hammondia hammondi reveal dramatic differences in the host response and effectors with species-specific functions

Toxoplasma gondii and Hammondia hammondi are closely-related coccidian intracellular parasites that differ in their ability to cause disease in animal and (likely) humans. The role of the host response in these phenotypic differences is not known and to address this we performed a transcriptomic analysis of a monocyte cell line (THP-1) infected with these two parasite species. The pathways altered by infection were shared between species ~95% the time, but the magnitude of the host response to H. hammondi was significantly higher compared to T. gondii. Accompanying this divergent host response was an equally divergent impact on the cell cycle of the host cell. In contrast to T. gondii, H. hammondi infection induces cell cycle arrest via pathways linked to DNA-damage responses and cellular senescence and robust secretion of multiple chemokines that are known to be a part of the senescence associated secretory phenotype (SASP). Remarkably, prior T. gondii infection or treatment with T. gondii-conditioned media suppressed responses to H. hammondi infection, and promoted the replication of H. hammondi in recipient cells. Suppression of inflammatory responses to H. hammondi was found to be mediated by the T. gondii effector IST, and this finding was consistent with reduced functionality of the H. hammondi IST ortholog compared to its T. gondii counterpart. Taken together our data suggest that T. gondii manipulation of the host cell is capable of suppressing previously unknown stress and/or DNA-damage induced responses that occur during infection with H. hammondi, and that one important impact of this T. gondii mediated suppression is to promote parasite replication.

S-phase Enriched Non-coding RNAs Regulate Gene Expression and Cell Cycle Progression

Many proteins that are needed for progression through S-phase are produced from transcripts that peak in the S-phase, linking temporal expression of those proteins to the time that they are required in cell cycle. Here, we explore the potential roles of long non-coding RNAs in cell cycle progression. We use a sensitive click-chemistry approach to isolate nascent RNAs in a human cell line, and we identify more than 900 long non-coding RNAs (lncRNAs) whose synthesis peaks during the S-phase. More than 200 of these are long intergenic non-coding RNAs (lincRNAs) with S-phase-specific expression. We characterize three of these lincRNAs by knockdown and find that all three lincRNAs are required for appropriate S-phase progression. We infer that non-coding RNAs are key regulatory effectors during the cell cycle, acting on distinct regulatory networks, and herein, we provide a large catalog of candidate cell-cycle regulatory RNAs.

Personalized Network Modeling of the Pan-Cancer Patient and Cell Line Interactome

PURPOSE

Personalized network inference on diverse clinical and in vitro model systems across cancer types can be used to delineate specific regulatory mechanisms, uncover drug targets and pathways, and develop individualized predictive models in cancer.

METHODS

We developed TransPRECISE (personalized cancer-specific integrated network estimation model), a multiscale Bayesian network modeling framework, to analyze the pan-cancer patient and cell line interactome to identify differential and conserved intrapathway activities, to globally assess cell lines as representative models for patients, and to develop drug sensitivity prediction models. We assessed pan-cancer pathway activities for a large cohort of patient samples (> 7,700) from the Cancer Proteome Atlas across ≥ 30 tumor types, a set of 640 cancer cell lines from the MD Anderson Cell Lines Project spanning 16 lineages, and ≥ 250 cell lines' response to > 400 drugs.

RESULTS

TransPRECISE captured differential and conserved proteomic network topologies and pathway circuitry between multiple patient and cell line lineages: ovarian and kidney cancers shared high levels of connectivity in the hormone receptor and receptor tyrosine kinase pathways, respectively, between the two model systems. Our tumor stratification approach found distinct clinical subtypes of the patients represented by different sets of cell lines: patients with head and neck tumors were classified into two different subtypes that are represented by head and neck and esophagus cell lines and had different prognostic patterns (456 v 654 days of median overall survival; P = .02). High predictive accuracy was observed for drug sensitivities in cell lines across multiple drugs (median area under the receiver operating characteristic curve > 0.8) using Bayesian additive regression tree models with TransPRECISE pathway scores.

CONCLUSION

Our study provides a generalizable analytic framework to assess the translational potential of preclinical model systems and to guide pathway-based personalized medical decision making, integrating genomic and molecular data across model systems.

Enigmatic MELK: The controversy surrounding its complex role in cancer

The Ser/Thr protein kinase MELK (maternal embryonic leucine zipper kinase) has been considered an attractive therapeutic target for managing cancer since 2005. Studies using expression analysis have indicated that MELK expression is higher in numerous cancer cells and tissues than in their normal, nonneoplastic counterparts. Further, RNAi-mediated MELK depletion impairs proliferation of multiple cancers, including triple-negative breast cancer (TNBC), and these growth defects can be rescued with exogenous WT MELK, but not kinase-dead MELK complementation. Pharmacological MELK inhibition with OTS167 (alternatively called OTSSP167) and NVS-MELK8a, among other small molecules, also impairs cancer cell growth. These collective results led to MELK being classified as essential for cancer proliferation. More recently, in 2017, the proliferation of TNBC and other cancer cell lines was reported to be unaffected by genetic CRISPR/Cas9-mediated MELK deletion, calling into question the essentiality of this kinase in cancer. To date, the requirement of MELK in cancer remains controversial, and mechanisms underlying the disparate growth effects observed with RNAi, pharmacological inhibition, and CRISPR remain unclear. Our objective with this review is to highlight the evidence on both sides of this controversy, to provide commentary on the purported requirement of MELK in cancer, and to emphasize the need for continued elucidation of the functions of MELK.

FOXM1 Deubiquitination by USP21 Regulates Cell Cycle Progression and Paclitaxel Sensitivity in Basal-like Breast Cancer

The transcription factor FOXM1 contributes to cell cycle progression and is significantly upregulated in basal-like breast cancer (BLBC). Despite its importance in normal and cancer cell cycles, we lack a complete understanding of mechanisms that regulate FOXM1. We identified USP21 in an RNAi-based screen for deubiquitinases that control FOXM1 abundance. USP21 increases the stability of FOXM1, and USP21 binds and deubiquitinates FOXM1 in vivo and in vitro, indicating a direct enzyme-substrate relationship. Depleting USP21 downregulates the FOXM1 transcriptional network and causes a signifi-cant delay in cell cycle progression. Significantly, USP21 depletion sensitized BLBC cell lines and mouse xenograft tumors to paclitaxel, an anti-mitotic, frontline therapy in BLBC treatment. USP21 is the most frequently amplified deubiquitinase in BLBC patient tumors, and its amplification co-occurs with the upregulation of FOXM1 protein. Altogether, these data suggest a role for USP21 in the proliferation and potentially treatment of FOXM1-high, USP21-high BLBC.

The cell cycle transcription factor FOXM1 is activated in basal-like breast cancer (BLBC) and associated with therapeutic resistance and poor patient outcomes. Arceci et al. show USP21 antagonizes FOXM1 degradation, thereby promoting proliferation and paclitaxel resistance. USP21 is catalytically active and recurrently overexpressed in BLBC, representing a potential therapeutic target.

Enhanced identification of significant regulators of gene expression

Background

Diseases like cancer will lead to changes in gene expression, and it is relevant to identify key regulatory genes that can be linked directly to these changes. This can be done by computing a Regulatory Impact Factor (RIF) score for relevant regulators. However, this computation is based on estimating correlated patterns of gene expression, often Pearson correlation, and an assumption about a set of specific regulators, normally transcription factors. This study explores alternative measures of correlation, using the Fisher and Sobolev metrics, and an extended set of regulators, including epigenetic regulators and long non-coding RNAs (lncRNAs). Data on prostate cancer have been used to explore the effect of these modifications.

Results

A tool for computation of RIF scores with alternative correlation measures and extended sets of regulators was developed and tested on gene expression data for prostate cancer. The study showed that the Fisher and Sobolev metrics lead to improved identification of well-documented regulators of gene expression in prostate cancer, and the sets of identified key regulators showed improved overlap with previously defined gene sets of relevance to cancer. The extended set of regulators lead to identification of several interesting candidates for further studies, including lncRNAs. Several key processes were identified as important, including spindle assembly and the epithelial-mesenchymal transition (EMT).

Conclusions

The study has shown that using alternative metrics of correlation can improve the performance of tools based on correlation of gene expression in genomic data. The Fisher and Sobolev metrics should be considered also in other correlation-based applications.

Cell cycle transcriptomics of Capsaspora provides insights into the evolution of cyclin-CDK machinery

Progression through the cell cycle in eukaryotes is regulated on multiple levels. The main driver of the cell cycle progression is the periodic activity of cyclin-dependent kinase (CDK) complexes. In parallel, transcription during the cell cycle is regulated by a transcriptional program that ensures the just-in-time gene expression. Many core cell cycle regulators are widely conserved in eukaryotes, among them cyclins and CDKs; however, periodic transcriptional programs are divergent between distantly related species. In addition, many otherwise conserved cell cycle regulators have been lost and independently evolved in yeast, a widely used model organism for cell cycle research. For a better understanding of the evolution of the cell cycle regulation in opisthokonts, we investigated the transcriptional program during the cell cycle of the filasterean Capsaspora owczarzaki, a unicellular species closely related to animals. We developed a protocol for cell cycle synchronization in Capsaspora cultures and assessed gene expression over time across the entire cell cycle. We identified a set of 801 periodic genes that grouped into five clusters of expression over time. Comparison with datasets from other eukaryotes revealed that the periodic transcriptional program of Capsaspora is most similar to that of animal cells. We found that orthologues of cyclin A, B and E are expressed at the same cell cycle stages as in human cells and in the same temporal order. However, in contrast to human cells where these cyclins interact with multiple CDKs, Capsaspora cyclins likely interact with a single ancestral CDK1-3. Thus, the Capsaspora cyclin-CDK system could represent an intermediate state in the evolution of animal-like cyclin-CDK regulation. Overall, our results demonstrate that Capsaspora could be a useful unicellular model system for animal cell cycle regulation.

Mass spectrometry-based selectivity profiling identifies a highly selective inhibitor of the kinase MELK that delays mitotic entry in cancer cells

The maternal embryonic leucine zipper kinase (MELK) has been implicated in the regulation of cancer cell proliferation. RNAi-mediated MELK depletion impairs growth and causes G2/M arrest in numerous cancers, but the mechanisms underlying these effects are poorly understood. Furthermore, the MELK inhibitor OTSSP167 has recently been shown to have poor selectivity for MELK, complicating the use of this inhibitor as a tool compound to investigate MELK function. Here, using a cell-based proteomics technique called multiplexed kinase inhibitor beads/mass spectrometry (MIB/MS), we profiled the selectivity of two additional MELK inhibitors, NVS-MELK8a (8a) and HTH-01-091. Our results revealed that 8a is a highly selective MELK inhibitor, which we further used for functional studies. Resazurin and crystal violet assays indicated that 8a decreases triple-negative breast cancer cell viability, and immunoblotting revealed that impaired growth is due to perturbation of cell cycle progression rather than induction of apoptosis. Using double-thymidine synchronization and immunoblotting, we observed that MELK inhibition delays mitotic entry, which was associated with delayed activation of Aurora A, Aurora B, and cyclin-dependent kinase 1 (CDK1). Following this delay, cells entered and completed mitosis. Using live-cell microscopy of cells harboring fluorescent proliferating cell nuclear antigen, we confirmed that 8a significantly and dose-dependently lengthens G2 phase. Collectively, our results provide a rationale for using 8a as a tool compound for functional studies of MELK and indicate that MELK inhibition delays mitotic entry, likely via transient G2/M checkpoint activation.

Improved recovery of cell-cycle gene expression in Saccharomyces cerevisiae from regulatory interactions in multiple omics data

BACKGROUND

Gene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task.

RESULTS

We employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression in Saccharomyces cerevisiae. Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators (STE12-TEC1 and RAP1-HAP1-MSN4).

CONCLUSION

Our findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

Analysis of five deep-sequenced trio-genomes of the Peninsular Malaysia Orang Asli and North Borneo populations

Background

Recent advances in genomic technologies have facilitated genome-wide investigation of human genetic variations. However, most efforts have focused on the major populations, yet trio genomes of indigenous populations from Southeast Asia have been under-investigated.

Results

We analyzed the whole-genome deep sequencing data (~ 30×) of five native trios from Peninsular Malaysia and North Borneo, and characterized the genomic variants, including single nucleotide variants (SNVs), small insertions and deletions (indels) and copy number variants (CNVs). We discovered approximately 6.9 million SNVs, 1.2 million indels, and 9000 CNVs in the 15 samples, of which 2.7% SNVs, 2.3% indels and 22% CNVs were novel, implying the insufficient coverage of population diversity in existing databases. We identified a higher proportion of novel variants in the Orang Asli (OA) samples, i.e., the indigenous people from Peninsular Malaysia, than that of the North Bornean (NB) samples, likely due to more complex demographic history and long-time isolation of the OA groups. We used the pedigree information to identify de novo variants and estimated the autosomal mutation rates to be 0.81 × 10^− 8 – 1.33 × 10^− 8, 1.0 × 10^− 9 – 2.9 × 10^− 9, and ~ 0.001 per site per generation for SNVs, indels, and CNVs, respectively. The trio-genomes also allowed for haplotype phasing with high accuracy, which serves as references to the future genomic studies of OA and NB populations. In addition, high-frequency inherited CNVs specific to OA or NB were identified. One example is a 50-kb duplication in DEFA1B detected only in the Negrito trios, implying plausible effects on host defense against the exposure of diverse microbial in tropical rainforest environment of these hunter-gatherers. The CNVs shared between OA and NB groups were much fewer than those specific to each group. Nevertheless, we identified a 142-kb duplication in AMY1A in all the 15 samples, and this gene is associated with the high-starch diet. Moreover, novel insertions shared with archaic hominids were identified in our samples.

Conclusion

Our study presents a full catalogue of the genome variants of the native Malaysian populations, which is a complement of the genome diversity in Southeast Asians. It implies specific population history of the native inhabitants, and demonstrated the necessity of more genome sequencing efforts on the multi-ethnic native groups of Malaysia and Southeast Asia.

A novel FoxM1-PSMB4 axis contributes to proliferation and progression of cervical cancer

The abnormally high activity of the proteasome system is closely related to the occurrence and development of various tumors. PSMB4 is a non-catalytic subunit for the proteasome assembly. Although the reports from genetic screening have demonstrated it's a driver gene for cell growth in several types of solid tumor, its expression pattern and regulatory mechanisms in malignant diseases are still elusive. Here, we found that PSMB4 is overexpressed in cervical cancer tissues. And knockdown of PSMB4 significantly inhibited cervical cancer cell proliferation. The mechanistic study revealed that FoxM1, a master regulator of cell division, binds directly to the promoter region of PSMB4 and regulates the PSMB4 expression in the mRNA level. In addition, the data analysis from TCGA showed a positive correlation between FxoM1 and PSMB4 in cervical cancer. Furthermore, the loss of functional and rescue experiments confirmed that PSMB4 is required for FoxM1-driven cervical cancer cell proliferation. Collectively, our study explains the phenomenon of dysregulated expression of PSMB4 in cervical cancer tissues and verifies its driver effect on cancer cell proliferation. More importantly, it highlights a FoxM1-PSMB4 axis could be a potential target for the treatment of cervical cancer.

Preparation for DNA replication: the key to a successful S phase

Successful genome duplication is required for cell proliferation and demands extraordinary precision and accuracy. The mechanisms by which cells enter, progress through, and exit S phase are intense areas of focus in the cell cycle and genome stability fields. Key molecular events in the G1 phase of the cell division cycle, especially origin licensing, are essential for pre-establishing conditions for efficient DNA replication during the subsequent S phase. If G1 events are poorly regulated or disordered, then DNA replication can be compromised leading to genome instability, a hallmark of tumorigenesis. Upon entry into S phase, coordinated origin firing and replication progression ensure complete, timely, and precise chromosome replication. Both G1 and S phase progressions are controlled by master cell cycle protein kinases and ubiquitin ligases that govern the activity and abundance of DNA replication factors. In this short review, we describe current understanding and recent developments related to G1 progression and S phase entrance and exit with a particular focus on origin licensing regulation in vertebrates.

Assembly of a parts list of the human mitotic cell cycle machinery

The set of proteins required for mitotic division remains poorly characterized. Here, an extensive series of correlation analyses of human and mouse transcriptomics data were performed to identify genes strongly and reproducibly associated with cells undergoing S/G2-M phases of the cell cycle. In so doing, 701 cell cycle-associated genes were defined and while it was shown that many are only expressed during these phases, the expression of others is also driven by alternative promoters. Of this list, 496 genes have known cell cycle functions, whereas 205 were assigned as putative cell cycle genes, 53 of which are functionally uncharacterized. Among these, 27 were screened for subcellular localization revealing many to be nuclear localized and at least three to be novel centrosomal proteins. Furthermore, 10 others inhibited cell proliferation upon siRNA knockdown. This study presents the first comprehensive list of human cell cycle proteins, identifying many new candidate proteins.

Tyrosine Threonine Kinase Inhibition Eliminates Lung Cancers by Augmenting Apoptosis and Polyploidy

The spindle assembly checkpoint maintains genomic integrity. A key component is tyrosine threonine kinase (TTK, also known as Mps1). TTK antagonism is hypothesized to cause genomic instability and cell death. Interrogating The Cancer Genome Atlas revealed high TTK expression in lung adenocarcinomas and squamous cell cancers versus the normal lung (P < 0.001). This correlated with an unfavorable prognosis in examined lung adenocarcinoma cases (P = 0.007). TTK expression profiles in lung tumors were independently assessed by RNA in situ hybridization. CFI-402257 is a highly selective TTK inhibitor. Its potent antineoplastic effects are reported here against a panel of well-characterized murine and human lung cancer cell lines. Significant antitumorigenic activity followed independent treatments of athymic mice bearing human lung cancer xenografts (6.5 mg/kg, P < 0.05; 8.5 mg/kg, P < 0.01) and immunocompetent mice with syngeneic lung cancers (P < 0.001). CFI-402257 antineoplastic mechanisms were explored. CFI-402257 triggered aneuploidy and apoptotic death of lung cancer cells without changing centrosome number. Reverse phase protein arrays (RPPA) of vehicle versus CFI-402257-treated lung cancers were examined using more than 300 critical growth-regulatory proteins. RPPA bioinformatic analyses discovered CFI-402257 enhanced MAPK signaling, implicating MAPK antagonism in augmenting TTK inhibitory effects. This was independently confirmed using genetic and pharmacologic repression of MAPK that promoted CFI-402257 anticancer actions. TTK antagonism exerted marked antineoplastic effects against lung cancers and MAPK inhibition cooperated. Future work should determine whether CFI-402257 treatment alone or with a MAPK inhibitor is active in the lung cancer clinic.

Vitamin B12 , folate, and the methionine remethylation cycle-biochemistry, pathways, and regulation

Vitamin B₁₂ (cobalamin, Cbl) is a nutrient essential to human health. Due to its complex structure and dual cofactor forms, Cbl undergoes a complicated series of absorptive and processing steps before serving as cofactor for the enzymes methylmalonyl-CoA mutase and methionine synthase. Methylmalonyl-CoA mutase is required for the catabolism of certain (branched-chain) amino acids into an anaplerotic substrate in the mitochondrion, and dysfunction of the enzyme itself or in production of its cofactor adenosyl-Cbl result in an inability to successfully undergo protein catabolism with concomitant mitochondrial energy disruption. Methionine synthase catalyzes the methyl-Cbl dependent (re)methylation of homocysteine to methionine within the methionine cycle; a reaction required to produce this essential amino acid and generate S-adenosylmethionine, the most important cellular methyl-donor. Disruption of methionine synthase has wide-ranging implications for all methylation-dependent reactions, including epigenetic modification, but also for the intracellular folate pathway, since methionine synthase uses 5-methyltetrahydrofolate as a one-carbon donor. Folate-bound one-carbon units are also required for deoxythymidine monophosphate and de novo purine synthesis; therefore, the flow of single carbon units to each of these pathways must be regulated based on cellular needs. This review provides an overview on Cbl metabolism with a brief description of absorption and intracellular metabolic pathways. It also provides a description of folate-mediated one-carbon metabolism and its intersection with Cbl at the methionine cycle. Finally, a summary of recent advances in understanding of how both pathways are regulated is presented.

Unbiased Boolean analysis of public gene expression data for cell cycle gene identification

Cell proliferation is essential for the development and maintenance of all organisms and is dysregulated in cancer. Using synchronized cells progressing through the cell cycle, pioneering microarray studies defined cell cycle genes based on cyclic variation in their expression. However, the concordance of the small number of synchronized cell studies has been limited, leading to discrepancies in definition of the transcriptionally regulated set of cell cycle genes within and between species. Here we present an informatics approach based on Boolean logic to identify cell cycle genes. This approach used the vast array of publicly available gene expression data sets to query similarity to CCNB1, which encodes the cyclin subunit of the Cdk1-cyclin B complex that triggers the G2-to-M transition. In addition to highlighting conservation of cell cycle genes across large evolutionary distances, this approach identified contexts where well-studied genes known to act during the cell cycle are expressed and potentially acting in nondivision contexts. An accessible web platform enables a detailed exploration of the cell cycle gene lists generated using the Boolean logic approach. The methods employed are straightforward to extend to processes other than the cell cycle.

CINSARC signature as a prognostic marker for clinical outcome in sarcomas and beyond

Prognostication is a key issue for sarcoma patients' care as it triggers the therapeutic approach including chemotherapy, which is still not standard for localized patients. Current prognostic evaluation, based on the FNCLCC grading system, has recently been improved by the CINSARC signature outperforming histology-based grading system by identifying high-risk patients in every grade, even in those considered as low. CINSARC is an expression-based signature related to mitosis and chromosome integrity with prognostic value in a wide range of cancers additional to sarcoma. First developed with frozen material, CINSARC is now coupled with NanoString technology allowing evaluation from FFPE blocks used in clinical practice. Consequently, CINSARC is currently evaluated in clinical trials with a dual objective of demonstrating the benefit of chemotherapy in sarcoma patients and testing its response prediction. Considering its overarching value in oncology, its development is welcome in any cancers where the prognostication needs to be improved.

Mapping Metabolic Events in the Cancer Cell Cycle Reveals Arginine Catabolism in the Committed SG2M Phase

Alterations in cell-cycle regulation and cellular metabolism are associated with cancer transformation, and enzymes active in the committed cell-cycle phase may represent vulnerabilities of cancer cells. Here, we map metabolic events in the G1 and SG2M phases by combining cell sorting with mass spectrometry-based isotope tracing, revealing hundreds of cell-cycle-associated metabolites. In particular, arginine uptake and ornithine synthesis are active during SG2M in transformed but not in normal cells, with the mitochondrial arginase 2 (ARG2) enzyme as a potential mechanism. While cancer cells exclusively use ARG2, normal epithelial cells synthesize ornithine via ornithine aminotransferase (OAT). Knockdown of ARG2 markedly reduces cancer cell growth and causes G2M arrest, while not inducing compensation via OAT. In human tumors, ARG2 is highly expressed in specific tumor types, including basal-like breast tumors. This study sheds light on the interplay between metabolism and cell cycle and identifies ARG2 as a potential metabolic target.

Cell Cycle-Regulated Transcription of CENP-A by the MBF Complex Ensures Optimal Level of CENP-A for Centromere Formation

The centromere plays an essential role in chromosome segregation. In most eukaryotes, centromeres are epigenetically defined by the conserved histone H3 variant CENP-A. Proper centromere assembly is dependent upon the tight regulation of CENP-A level. Cell cycle regulation of CENP-A transcription appears to be a universal feature across eukaryotes, but the molecular mechanism underlying the temporal control of CENP-A transcription and how such regulation contributes to centromere function remains elusive. CENP-A in fission yeast has been shown to be transcribed before S phase. Using various synchronization methods, we confirmed that CENP-A transcription occurs at G1, leading to an almost twofold increase of the protein during S phase. Through a genetic screen, we identified the MBF (MluI box-binding factors) complex as a key regulator of temporal control of CENP-A transcription. The periodic transcription of CENP-A is lost in MBF mutants, resulting in CENP-A mislocalization and chromosome segregation defects. We identified the MCB (MluI cell cycle box) motif in the CENP-A promoter, and further showed that the MBF complex binds to the motif to restrict CENP-A transcription to G1. Mutations of the MCB motif cause constitutive CENP-A expression and deleterious effects on cell survival. Using promoters driving transcription to different cell cycle stages, we found that timing of CENP-A transcription is dispensable for its centromeric localization. Our data instead indicate that cell cycle-regulated CENP-A transcription is a key step to ensure that a proper amount of CENP-A is generated across generations. This study provides mechanistic insights into the regulation of cell cycle-dependent CENP-A transcription, as well as its importance on centromere function.

Developmental Heterogeneity of Microglia and Brain Myeloid Cells Revealed by Deep Single-Cell RNA Sequencing

Microglia are increasingly recognized for their major contributions during brain development and neurodegenerative disease. It is currently unknown whether these functions are carried out by subsets of microglia during different stages of development and adulthood or within specific brain regions. Here, we performed deep single-cell RNA sequencing (scRNA-seq) of microglia and related myeloid cells sorted from various regions of embryonic, early postnatal, and adult mouse brains. We found that the majority of adult microglia expressing homeostatic genes are remarkably similar in transcriptomes, regardless of brain region. By contrast, early postnatal microglia are more heterogeneous. We discovered a proliferative-region-associated microglia (PAM) subset, mainly found in developing white matter, that shares a characteristic gene signature with degenerative disease-associated microglia (DAM). Such PAM have amoeboid morphology, are metabolically active, and phagocytose newly formed oligodendrocytes. This scRNA-seq atlas will be a valuable resource for dissecting innate immune functions in health and disease.

FOXM1 modulates 5-fluorouracil sensitivity in cholangiocarcinoma through thymidylate synthase (TYMS): implications of FOXM1-TYMS axis uncoupling in 5-FU resistance

Fluorouracil (5-FU) is the first-line chemotherapeutic drug for cholangiocarcinoma (CCA), but its efficacy has been compromised by the development of resistance. Development of 5-FU resistance is associated with elevated expression of its cellular target, thymidylate synthase (TYMS). E2F1 transcription factor has previously been shown to modulate the expression of FOXM1 and TYMS. Immunohistochemical (IHC) analysis revealed a strong correlated upregulation of FOXM1 (78%) and TYMS (48%) expression at the protein levels in CCA tissues. In agreement, RT-qPCR and western blot analyses of four human CCA cell lines at the baseline level and in response to high doses of 5-FU revealed good correlations between FOXM1 and TYMS expression in the CCA cell lines tested, except for the highly 5-FU-resistant HuCCA cells. Consistently, siRNA-mediated knockdown of FOXM1 reduced the clonogenicity and TYMS expression in the relatively sensitive KKU-D131 but not in the highly resistant HuCCA cells. Interestingly, silencing of TYMS sensitized both KKU-D131 and HuCCA to 5-FU treatment, suggesting that resistance to very high levels of 5-FU is due to the inability of the genotoxic sensor FOXM1 to modulate TYMS expression. Consistently, ChIP analysis revealed that FOXM1 binds efficiently to the TYMS promoter and modulates TYMS expression at the promoter level upon 5-FU treatment in KKU-D131 but not in HuCCA cells. In addition, E2F1 expression did not correlate with either FOXM1 or TYMS expression and E2F1 depletion has no effects on the clonogenicity and TYMS expression in the CCA cells. In conclusion, our data show that FOXM1 regulates TYMS expression to modulate 5-FU resistance in CCA and that severe 5-FU resistance can be caused by the uncoupling of the regulation of TYMS by FOXM1. Our findings suggest that the FOXM1–TYMS axis can be a novel diagnostic, predictive and prognostic marker as well as a therapeutic target for CCA.

Layers of regulation of cell-cycle gene expression in the budding yeast Saccharomyces cerevisiae

In the budding yeast Saccharomyces cerevisiae, transcription factors (TFs) regulate the periodic expression of many genes during the cell cycle, including gene products required for progression through cell-cycle events. Experimental evidence coupled with quantitative models suggests that a network of interconnected TFs is capable of regulating periodic genes over the cell cycle. Importantly, these dynamical models were built on transcriptomics data and assumed that TF protein levels and activity are directly correlated with mRNA abundance. To ask whether TF transcripts match protein expression levels as cells progress through the cell cycle, we applied a multiplexed targeted mass spectrometry approach (parallel reaction monitoring) to synchronized populations of cells. We found that protein expression of many TFs and cell-cycle regulators closely followed their respective mRNA transcript dynamics in cycling wild-type cells. Discordant mRNA/protein expression dynamics was also observed for a subset of cell-cycle TFs and for proteins targeted for degradation by E3 ubiquitin ligase complexes such as SCF (Skp1/Cul1/F-box) and APC/C (anaphase-promoting complex/cyclosome). We further profiled mutant cells lacking B-type cyclin/CDK activity ( clb1-6) where oscillations in ubiquitin ligase activity, cyclin/CDKs, and cell-cycle progression are halted. We found that a number of proteins were no longer periodically degraded in clb1-6 mutants compared with wild type, highlighting the importance of posttranscriptional regulation. Finally, the TF complexes responsible for activating G1/S transcription (SBF and MBF) were more constitutively expressed at the protein level than at periodic mRNA expression levels in both wild-type and mutant cells. This comprehensive investigation of cell-cycle regulators reveals that multiple layers of regulation (transcription, protein stability, and proteasome targeting) affect protein expression dynamics during the cell cycle.

Cell cycle regulation of folate-mediated one-carbon metabolism

Folate-mediated one-carbon metabolism (FOCM) comprises a network of interconnected folate-dependent metabolic pathways responsible for serine and glycine interconversion, de novo purine synthesis, de novo thymidylate synthesis and homocysteine remethylation to methionine. These pathways are compartmentalized in the cytosol, nucleus and mitochondria. Individual enzymes within the FOCM network compete for folate cofactors because intracellular folate concentrations are limiting. Although there are feedback mechanisms that regulate the partitioning of folate cofactors among the folate-dependent pathways, less recognized is the impact of cell cycle regulation on FOCM. This review summarizes the evidence for temporal regulation of expression, activity and cellular localization of enzymes and pathways in the FOCM network in mammalian cells through the cell cycle. This article is categorized under: Biological Mechanisms > Metabolism Physiology > Mammalian Physiology in Health and Disease.

Who guards the guardian? Mechanisms that restrain APC/C during the cell cycle

The cell cycle is principally controlled by Cyclin Dependent Kinases (CDKs), whose oscillating activities are determined by binding to Cyclin coactivators. Cyclins exhibit dynamic changes in abundance as cells pass through the cell cycle. The sequential, timed accumulation and degradation of Cyclins, as well as many other proteins, imposes order on the cell cycle and contributes to genome maintenance. The destruction of many cell cycle regulated proteins, including Cyclins A and B, is controlled by a large, multi-subunit E3 ubiquitin ligase termed the Anaphase Promoting Complex/Cyclosome (APC/C). APC/C activity is tightly regulated during the cell cycle. Its activation state increases dramatically in mid-mitosis and it remains active until the end of G1 phase. Following its mandatory inactivation at the G1/S boundary, APC/C activity remains low until the subsequent mitosis. Due to its role in guarding against the inappropriate or untimely accumulation of Cyclins, the APC/C is a core component of the cell cycle oscillator. In addition to the regulation of Cyclins, APC/C controls the degradation of many other substrates. Therefore, it is vital that the activity of APC/C itself be tightly guarded. The APC/C is most well studied for its role and regulation during mitosis. However, the APC/C also plays a similarly important and conserved role in the maintenance of G1 phase. Here we review the diverse mechanisms counteracting APC/C activity throughout the cell cycle and the importance of their coordinated actions on cell growth, proliferation, and disease.

GABPA inhibits invasion/metastasis in papillary thyroid carcinoma by regulating DICER1 expression

The ETS family transcription factor GABPA is suggested as an oncogenic element, which is further supported by the recent reporting of it as the sole ETS member to activate the mutant TERT promoter in thyroid carcinomas (TC). However, it remains unclear how GABPA contributes to TC pathogenesis. The present study is designed to address this issue. TERT expression was significantly diminished in TERT promoter-mutated TC cells upon GABPA inhibition. Surprisingly, GABPA depletion led to robustly increased cellular invasion independently of TERT promoter mutations and TERT expression. DICER1, a component of the microRNA machinery, was identified as a downstream effector of GABPA. GABPA facilitated Dicer1 transcription while its depletion reduced Dicer1 expression. The mutation of the GABPA binding site in the DICER1 promoter led to diminished basal levels of DICER1 promoter activity and abolishment of GABPA-stimulated promoter activity as well. The forced DICER1 expression abrogated the invasiveness of GABPA-depleted TC cells. Consistently, the analyses of 93 patients with papillary thyroid carcinoma (PTC) revealed a positive correlation between GABPA and DICER1 expression. GABPA expression was negatively associated with TERT expression and promoter mutations, in contrast to published observations in cancer cell lines. Lower GABPA expression was associated with distant metastasis and shorter overall/disease-free survival in PTC patients. Similar results were obtained for PTC cases in the TCGA dataset. In addition, a positive correlation between GABPA and DICER1 expression was seen in multiple types of malignancies. Taken together, despite its stimulatory effect on the mutant TERT promoter and telomerase activation, GABPA may itself act as a tumor suppressor rather than an oncogenic factor to inhibit invasion/metastasis in TCs and be a useful predictor for patient outcomes.

Cezanne/OTUD7B is a cell cycle-regulated deubiquitinase that antagonizes the degradation of APC/C substrates

The anaphase-promoting complex/cyclosome (APC/C) is an E3 ubiquitin ligase and key regulator of cell cycle progression. Since APC/C promotes the degradation of mitotic cyclins, it controls cell cycle-dependent oscillations in cyclin-dependent kinase (CDK) activity. Both CDKs and APC/C control a large number of substrates and are regulated by analogous mechanisms, including cofactor-dependent activation. However, whereas substrate dephosphorylation is known to counteract CDK, it remains largely unknown whether deubiquitinating enzymes (DUBs) antagonize APC/C substrate ubiquitination during mitosis. Here, we demonstrate that Cezanne/OTUD7B is a cell cycle-regulated DUB that opposes the ubiquitination of APC/C targets. Cezanne is remarkably specific for K11-linked ubiquitin chains, which are formed by APC/C in mitosis. Accordingly, Cezanne binds established APC/C substrates and reverses their APC/C-mediated ubiquitination. Cezanne depletion accelerates APC/C substrate degradation and causes errors in mitotic progression and formation of micronuclei. These data highlight the importance of tempered APC/C substrate destruction in maintaining chromosome stability. Furthermore, Cezanne is recurrently amplified and overexpressed in numerous malignancies, suggesting a potential role in genome maintenance and cancer cell proliferation.

The Temporal Regulation of S Phase Proteins During G1

Successful DNA replication requires intimate coordination with cell-cycle progression. Prior to DNA replication initiation in S phase, a series of essential preparatory events in G₁ phase ensures timely, complete, and precise genome duplication. Among the essential molecular processes are regulated transcriptional upregulation of genes that encode replication proteins, appropriate post-transcriptional control of replication factor abundance and activity, and assembly of DNA-loaded protein complexes to license replication origins. In this chapter we describe these critical G₁ events necessary for DNA replication and their regulation in the context of both cell-cycle entry and cell-cycle progression.

Identification of Drivers of Aneuploidy in Breast Tumors

Although aneuploidy is found in the majority of tumors, the degree of aneuploidy varies widely. It is unclear how cancer cells become aneuploid or how highly aneuploid tumors are different from those of more normal ploidy. We developed a simple computational method that measures the degree of aneuploidy or structural rearrangements of large chromosome regions of 522 human breast tumors from The Cancer Genome Atlas (TCGA). Highly aneuploid tumors overexpress activators of mitotic transcription and the genes encoding proteins that segregate chromosomes. Overexpression of three mitotic transcriptional regulators, E2F1, MYBL2, and FOXM1, is sufficient to increase the rate of lagging anaphase chromosomes in a non-transformed vertebrate tissue, demonstrating that this event can initiate aneuploidy. Highly aneuploid human breast tumors are also enriched in TP53 mutations. TP53 mutations co-associate with the overexpression of mitotic transcriptional activators, suggesting that these events work together to provide fitness to breast tumors.

Computational analysis of the oscillatory behavior at the translation level induced by mRNA levels oscillations due to finite intracellular resources

Recent studies have demonstrated how the competition for the finite pool of available gene expression factors has important effect on fundamental gene expression aspects. In this study, based on a whole-cell model simulation of translation in S. cerevisiae, we evaluate for the first time the expected effect of mRNA levels fluctuations on translation due to the finite pool of ribosomes. We show that fluctuations of a single gene or a group of genes mRNA levels induce periodic behavior in all S. cerevisiae translation factors and aspects: the ribosomal densities and the translation rates of all S. cerevisiae mRNAs oscillate. We numerically measure the oscillation amplitudes demonstrating that fluctuations of endogenous and heterologous genes can cause a significant fluctuation of up to 50% in the steady-state translation rates of the rest of the genes. Furthermore, we demonstrate by synonymous mutations that oscillating the levels of mRNAs that experience high ribosomal occupancy (e.g. ribosomal “traffic jam”) induces the largest impact on the translation of the S. cerevisiae genome.

The results reported here should provide novel insights and principles related to the design of synthetic gene expression circuits and related to the evolutionary constraints shaping gene expression of endogenous genes.

Each cell contains a limited number of macromolecules and factors that participate in the gene expression process. These expression resources are shared between the different molecules that encode the genetic code, resulting in non-trivial couplings and competitions between the different gene expression stages. Such competitions should be considered when analyzing the cellular economy of the cell, the genome evolution, and the design of synthetic expression circuits. Here we study the effect of couplings and competitions for ribosomes by performing a whole-cell simulation of translation of S. cerevisiae, with parameters estimated from experimental data. We demonstrate that by periodically changing the mRNA levels of a single gene (endogenous or heterologous) or a set of genes, the translation of all S. cerevisiae genes are affected in a periodic manner. We numerically estimate the exact impact of the mRNA levels periodicity on the translation process dynamics, as well as on the dynamics of the free ribosomal pool and the way it is affected by parameters such as the codon composition of the oscillating gene, its initiation rate and mRNA levels. Furthermore, we show that the codon compositions of synthetically highly expressed heterologous genes that are expected to oscillate must be carefully considered. For example, synonymous mutations resulting in “traffic jams” of ribosomes along the fluctuated mRNAs may cause significant fluctuations of up to 50% in the steady-state translation rates of all genes.

Separase Inhibitor Sepin-1 Inhibits Foxm1 Expression and Breast Cancer Cell Growth

Sepin-1, a potent non-competitive inhibitor of separase, inhibits cancer cell growth, but the mechanisms of Sepin-1-mediated growth inhibition are not fully understood. Here we report that Sepin-1 hinders growth of breast cancer cells, cell migration, and wound healing. Inhibition of cell growth induced by Sepin-1 in vitro doesn't appear to be through apoptosis but rather due to growth inhibition. Following Sepin-1 treatment caspases 3 and 7 are not activated and Poly (ADP-ribose) polymerase (Parp) is not cleaved. The expression of Forkhead box protein M1 (FoxM1), a transcription factor, and its target genes in the cell cycle, including Plk1, Cdk1, Aurora A, and Lamin B1, are reduced in a Sepin-1-dependent manner. Expressions of Raf kinase family members A-Raf, B-Raf, and C-Raf also are inhibited following treatment with Sepin-1. Raf is an intermediator in the Raf-Mek-Erk signaling pathway that phosphorylates FoxM1. Activated FoxM1 can promote its own transcription via a positive feedback loop. Sepin-1-induced downregulation of Raf and FoxM1 may inhibit expression of cell cycle-driving genes, resulting in inhibition of cell growth.

Treatment with docetaxel in combination with Aneustat leads to potent inhibition of metastasis in a patient-derived xenograft model of advanced prostate cancer

Background:

Docetaxel used for first-line treatment of advanced prostate cancer (PCa) is only marginally effective. We previously showed, using the LTL-313H subrenal capsule patient-derived metastatic PCa xenograft model, that docetaxel combined with Aneustat (OMN54), a multivalent plant-derived therapeutic, led to marked synergistic tumour growth inhibition. Here, we investigated the effect of docetaxel+Aneustat on metastasis.

Methods:

C4-2 cells were incubated with docetaxel, Aneustat and docetaxel+Aneustat to assess effects on cell migration. The LTL-313H model, similarly treated, was analysed for effects on lung micro-metastasis and kidney invasion. The LTL-313H gene expression profile was compared with profiles of PCa patients (obtained from Oncomine) and subjected to IPA to determine involvement of cancer driver genes.

Results:

Docetaxel+Aneustat markedly inhibited C4-2 cell migration and LTL-313H lung micro-metastasis/kidney invasion. Oncomine analysis indicated that treatment with docetaxel+Aneustat was associated with improved patient outcome. The drug combination markedly downregulated expression of cancer driver genes such as FOXM1 (and FOXM1-target genes). FOXM1 overexpression reduced the anti-metastatic activity of docetaxel+Aneustat.

Conclusions:

Docetaxel+Aneustat can inhibit PCa tissue invasion and metastasis. This activity appears to be based on reduced expression of cancer driver genes such as FOXM1. Use of docetaxel+Aneustat may provide a new, more effective regimen for therapy of metastatic PCa.

Transcriptome-wide analysis of the Trypanosoma cruzi proliferative cycle identifies the periodically expressed mRNAs and their multiple levels of control

Trypanosoma cruzi is the protozoan parasite causing American trypanosomiasis or Chagas disease, a neglected parasitosis with important human health impact in Latin America. The efficacy of current therapy is limited, and its toxicity is high. Since parasite proliferation is a fundamental target for rational drug design, we sought to progress into its understanding by applying a genome-wide approach. Treating a TcI linage strain with hydroxyurea, we isolated epimastigotes in late G1, S and G2/M cell cycle stages at 70% purity. The sequencing of each phase identified 305 stage-specific transcripts (1.5-fold change, p≤0.01), coding for conserved cell cycle regulated proteins and numerous proteins whose cell cycle dependence has not been recognized before. Comparisons with the parasite T. brucei and the human host reveal important differences. The meta-analysis of T. cruzi transcriptomic and ribonomic data indicates that cell cycle regulated mRNAs are subject to sub-cellular compartmentalization. Compositional and structural biases of these genes- including CAI, GC content, UTR length, and polycistron position- may contribute to their regulation. To discover nucleotide motifs responsible for the co-regulation of cell cycle regulated genes, we looked for overrepresented motifs at their UTRs and found a variant of the cell cycle sequence motif at the 3' UTR of most of the S and G2 stage genes. We additionally identified hairpin structures at the 5' UTRs of a high proportion of the transcripts, suggesting that periodic gene expression might also rely on translation initiation in T. cruzi. In summary, we report a comprehensive list of T. cruzi cell cycle regulated genes, including many previously unstudied proteins, we show evidence favoring a multi-step control of their expression, and we identify mRNA motifs that may mediate their regulation. Our results provide novel information of the T. cruzi proliferative proteins and the integrated levels of their gene expression control.

Cross platform analysis of transcriptomic data identifies ageing has distinct and opposite effects on tendon in males and females

The development of tendinopathy is influenced by a variety of factors including age, gender, sex hormones and diabetes status. Cross platform comparative analysis of transcriptomic data elucidated the connections between these entities in the context of ageing. Tissue-engineered tendons differentiated from bone marrow derived mesenchymal stem cells from young (20-24 years) and old (54-70 years) donors were assayed using ribonucleic acid sequencing (RNA-seq). Extension of the experiment to microarray and RNA-seq data from tendon identified gender specific gene expression changes highlighting disparity with existing literature and published pathways. Separation of RNA-seq data by sex revealed underlying negative binomial distributions which increased statistical power. Sex specific de novo transcriptome assemblies generated fewer larger transcripts that contained miRNAs, lincRNAs and snoRNAs. The results identify that in old males decreased expression of CRABP2 leads to cell proliferation, whereas in old females it leads to cellular senescence. In conjunction with existing literature the results explain gender disparity in the development and types of degenerative diseases as well as highlighting a wide range of considerations for the analysis of transcriptomic data. Wider implications are that degenerative diseases may need to be treated differently in males and females because alternative mechanisms may be involved.

Transcriptomic Characterization of the Human Cell Cycle in Individual Unsynchronized Cells

The highly fine-tuned dynamics of cell cycle gene expression have been intensely studied for several decades. However, some previous observations may be difficult to fully decouple from artifacts induced by traditional cell synchronization procedures. In addition, bulk cell measurements may have disguised intricate details. Here, we address this by sorting and transcriptomic sequencing of single cells progressing through the cell cycle without prior synchronization. Genes and pathways with known cell cycle roles are confirmed, associated regulatory sequence motifs are determined, and we also establish ties between other biological processes and the unsynchronized cell cycle. Importantly, we find the G1 phase to be surprisingly heterogeneous, with transcriptionally distinct early and late time points. We additionally note that mRNAs accumulate to reach maximum total levels at mitosis and find that stable transcripts show reduced cell-to-cell variability, consistent with the transcriptional burst model of gene expression. Our study provides the first detailed transcriptional profiling of an unsynchronized human cell cycle.

Accelerated parallel algorithm for gene network reverse engineering

Background

The Algorithm for the Reconstruction of Accurate Cellular Networks (ARACNE) represents one of the most effective tools to reconstruct gene regulatory networks from large-scale molecular profile datasets. However, previous implementations require intensive computing resources and, in some cases, restrict the number of samples that can be used. These issues can be addressed elegantly in a GPU computing framework, where repeated mathematical computation can be done efficiently, but requires extensive redesign to apply parallel computing techniques to the original serial algorithm, involving detailed optimization efforts based on a deep understanding of both hardware and software architecture.

Result

Here, we present an accelerated parallel implementation of ARACNE (GPU-ARACNE). By taking advantage of multi-level parallelism and the Compute Unified Device Architecture (CUDA) parallel kernel-call library, GPU-ARACNE successfully parallelizes a serial algorithm and simplifies the user experience from multi-step operations to one step. Using public datasets on comparable hardware configurations, we showed that GPU-ARACNE is faster than previous implementations and is able to reconstruct equally valid gene regulatory networks.

Conclusion

Given that previous versions of ARACNE are extremely resource demanding, either in computational time or in hardware investment, GPU-ARACNE is remarkably valuable for researchers who need to build complex regulatory networks from large expression datasets, but with limited budget on computational resources. In addition, our GPU-centered optimization of adaptive partitioning for Mutual Information (MI) estimation provides lessons that are applicable to other domains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0458-5) contains supplementary material, which is available to authorized users.

A map of protein dynamics during cell-cycle progression and cell-cycle exit

The cell-cycle field has identified the core regulators that drive the cell cycle, but we do not have a clear map of the dynamics of these regulators during cell-cycle progression versus cell-cycle exit. Here we use single-cell time-lapse microscopy of Cyclin-Dependent Kinase 2 (CDK2) activity followed by endpoint immunofluorescence and computational cell synchronization to determine the temporal dynamics of key cell-cycle proteins in asynchronously cycling human cells. We identify several unexpected patterns for core cell-cycle proteins in actively proliferating (CDK2-increasing) versus spontaneously quiescent (CDK2-low) cells, including Cyclin D1, the levels of which we find to be higher in spontaneously quiescent versus proliferating cells. We also identify proteins with concentrations that steadily increase or decrease the longer cells are in quiescence, suggesting the existence of a continuum of quiescence depths. Our single-cell measurements thus provide a rich resource for the field by characterizing protein dynamics during proliferation versus quiescence.

Cancer cell line specific co-factors modulate the FOXM1 cistrome

ChIP-seq has been commonly applied to identify genomic occupation of transcription factors (TFs) in a context-specific manner. It is generally assumed that a TF should have similar binding patterns in cells from the same or closely related tissues. Surprisingly, this assumption has not been carefully examined. To this end, we systematically compared the genomic binding of the cell cycle regulator FOXM1 in eight cell lines from seven different human tissues at binding signal, peaks and target genes levels. We found that FOXM1 binding in ER-positive breast cancer cell line MCF-7 are distinct comparing to those in not only other non-breast cell lines, but also MDA-MB-231, ER-negative breast cancer cell line. However, binding sites in MDA-MB-231 and non-breast cell lines were highly consistent. The recruitment of estrogen receptor alpha (ERα) caused the unique FOXM1 binding patterns in MCF-7. Moreover, the activity of FOXM1 in MCF-7 reflects the regulatory functions of ERα, while in MDA-MB-231 and non-breast cell lines, FOXM1 activities regulate cell proliferation. Our results suggest that tissue similarity, in some specific contexts, does not hold precedence over TF-cofactors interactions in determining transcriptional states and that the genomic binding of a TF can be dramatically affected by a particular co-factor under certain conditions.

Cell cycle transcription control: DREAM/MuvB and RB-E2F complexes

The precise timing of cell cycle gene expression is critical for the control of cell proliferation; de-regulation of this timing promotes the formation of cancer and leads to defects during differentiation and development. Entry into and progression through S phase requires expression of genes coding for proteins that function in DNA replication. Expression of a distinct set of genes is essential to pass through mitosis and cytokinesis. Expression of these groups of cell cycle-dependent genes is regulated by the RB pocket protein family, the E2F transcription factor family, and MuvB complexes together with B-MYB and FOXM1. Distinct combinations of these transcription factors promote the transcription of the two major groups of cell cycle genes that are maximally expressed either in S phase (G1/S) or in mitosis (G2/M). In this review, we discuss recent work that has started to uncover the molecular mechanisms controlling the precisely timed expression of these genes at specific cell cycle phases, as well as the repression of the genes when a cell exits the cell cycle.

Regulation of Poly(A) Tail and Translation during the Somatic Cell Cycle

Poly(A) tails are critical for mRNA stability and translation. However, recent studies have challenged this view, showing that poly(A) tail length and translation efficiency are decoupled in non-embryonic cells. Using TAIL-seq and ribosome profiling, we investigate poly(A) tail dynamics and translational control in the somatic cell cycle. We find dramatic changes in poly(A) tail lengths of cell-cycle regulatory genes like CDK1, TOP2A, and FBXO5, explaining their translational repression in M phase. We also find that poly(A) tail length is coupled to translation when the poly(A) tail is <20 nucleotides. However, as most genes have >20 nucleotide poly(A) tails, their translation is regulated mainly via poly(A) tail length-independent mechanisms during the cell cycle. Specifically, we find that terminal oligopyrimidine (TOP) tract-containing transcripts escape global translational suppression in M phase and are actively translated. Our quantitative and comprehensive data provide a revised view of translational control in the somatic cell cycle.

Functional interrelationship between TFII-I and E2F transcription factors at specific cell cycle gene loci

Transcription factor TFII-I is a multifunctional protein implicated in the regulation of cell cycle and stress-response genes. Previous studies have shown that a subset of TFII-I associated genomic sites contained DNA-binding motifs for E2F family transcription factors. We analyzed the co-association of TFII-I and E2Fs in more detail using bioinformatics, chromatin immunoprecipitation, and co-immunoprecipitation experiments. The data show that TFII-I interacts with E2F transcription factors. Furthermore, TFII-I, E2F4, and E2F6 interact with DNA-regulatory elements of several genes implicated in the regulation of the cell cycle, including DNMT1, HDAC1, CDKN1C, and CDC27. Inhibition of TFII-I expression led to a decrease in gene expression and in the association of E2F4 and E2F6 with these gene loci in human erythroleukemia K562 cells. Finally, TFII-I deficiency reduced the proliferation of K562 cells and increased the sensitivity toward doxorubicin toxicity. The results uncover novel interactions between TFII-I and E2Fs and suggest that TFII-I mediates E2F function at specific cell cycle genes.

A positive feedback loop of long noncoding RNA CCAT2 and FOXM1 promotes hepatocellular carcinoma growth

Hepatocellular carcinoma (HCC) is one of the most common malignancies around the world. Long noncoding RNAs (lncRNAs) are greater than 200 nucleotides without protein-coding potential and play critical roles in tumorigenesis, cell differentiation, and cancer metastasis. Colon cancer-associated transcript 2 (CCAT2), a newly identified lncRNA, was shown to be dysregulated in cancers. However, the functional role of CCAT2 in HCC remains questionable. In the present study, we found a significant upregulation of CCAT2 in HCC tissues as compared to non-tumor tissues. Functional assays showed that CCAT2 promotes cell growth in vivo and in vitro. In addition, we found a positive feedback loop between CCAT2 and FOXM1. CCAT2 upregulates FOXM1 expression through interaction with, and suppression of, miR-34a, and FOXM1 activates CCAT2 transcription. We evaluated the therapeutic potential of ultrasound-targeted microbubble destruction (UTMD)-mediated siRNA delivery to specifically target CCAT2. UTMD-mediated siCCAT2 delivery significantly suppressed tumor growth in vivo. Thus, CCAT2-FOXM1 may be a novel target for the treatment of HCC.