Cell cycle control during early embryogenesis

Pyridaben inhibits cell cycle progression and delays early embryonic development in zebrafish (Danio rerio)

Pyridaben is a broad-spectrum, contact-killing acaricide that can be used to control a variety of harmful food and plant mites. Pyridaben displays cardiotoxicity and liver toxicity toward fish, but the effects on fish embryonic development have not been characterized. We exposed early zebrafish embryos to 20, 30, and 40 μg/L concentrations of pyridaben. The exposure caused developmental abnormalities, including delayed embryonic shield formation, yolk sac resorption, decreases in body length, reduced pigmentation, and delays in hatching. Pyridaben caused a significant increase in the transcription level of the endoderm marker foxa2, but the transcription levels of the ectoderm development marker foxb1a and the mesoderm development marker snaila were not significantly altered. The transcription levels of the genes SOX17 in early embryos were significantly reduced. After exposure to pyridaben, catalase (CAT) activity and glutathione (GSH) content were increased, and cyclin D1, that is involved in early embryonic development, was abnormally expressed. This study shows that pyridaben causes anomalous development in zebrafish embryos by interfering with the cell cycle order of early embryonic development and inducing excessive oxidative stress. Colivelin, an agonist of the STAT3 signaling pathway, acted as a salvage drug to restore the cell cycle order during embryonic development following exposure to pyridaben. Thus, the toxic effects may be caused by pyridaben's regulation of the STAT3 signaling pathway.

Phosphorylation of the nuclear poly(A) binding protein (PABPN1) during mitosis protects mRNA from hyperadenylation and maintains transcriptome dynamics

Polyadenylation controls mRNA biogenesis, nucleo-cytoplasmic export, translation and decay. These processes are interdependent and coordinately regulated by poly(A)-binding proteins (PABPs), yet how PABPs are themselves regulated is not fully understood. Here, we report the discovery that human nuclear PABPN1 is phosphorylated by mitotic kinases at four specific sites during mitosis, a time when nucleoplasm and cytoplasm mix. To understand the functional consequences of phosphorylation, we generated a panel of stable cell lines inducibly over-expressing PABPN1 with point mutations at these sites. Phospho-inhibitory mutations decreased cell proliferation, highlighting the importance of PABPN1 phosphorylation in cycling cells. Dynamic regulation of poly(A) tail length and RNA stability have emerged as important modes of gene regulation. We therefore employed long-read sequencing to determine how PABPN1 phospho-site mutants affected poly(A) tails lengths and TimeLapse-seq to monitor mRNA synthesis and decay. Widespread poly(A) tail lengthening was observed for phospho-inhibitory PABPN1 mutants. In contrast, expression of phospho-mimetic PABPN1 resulted in shorter poly(A) tails with increased non-A nucleotides, in addition to increased transcription and reduced stability of a distinct cohort of mRNAs. Taken together, PABPN1 phosphorylation remodels poly(A) tails and increases mRNA turnover, supporting the model that enhanced transcriptome dynamics reset gene expression programs across the cell cycle.

ERK synchronizes embryonic cleavages in Drosophila

Extracellular-signal-regulated kinase (ERK) signaling controls development and homeostasis and is genetically deregulated in human diseases, including neurocognitive disorders and cancers. Although the list of ERK functions is vast and steadily growing, the full spectrum of processes controlled by any specific ERK activation event remains unknown. Here, we show how ERK functions can be systematically identified using targeted perturbations and global readouts of ERK activation. Our experimental model is the Drosophila embryo, where ERK signaling at the embryonic poles has thus far only been associated with the transcriptional patterning of the future larva. Through a combination of live imaging and phosphoproteomics, we demonstrated that ERK activation at the poles is also critical for maintaining the speed and synchrony of embryonic cleavages. The presented approach to interrogating phosphorylation networks identifies a hidden function of a well-studied signaling event and sets the stage for similar studies in other organisms.

Numerous Serine/Threonine Kinases Affect Blood Cell Homeostasis in Drosophila melanogaster

Blood cells in Drosophila serve primarily innate immune responses. Various stressors influence blood cell homeostasis regarding both numbers and the proportion of blood cell types. The principle molecular mechanisms governing hematopoiesis are conserved amongst species and involve major signaling pathways like Notch, Toll, JNK, JAK/Stat or RTK. Albeit signaling pathways generally rely on the activity of protein kinases, their specific contribution to hematopoiesis remains understudied. Here, we assess the role of Serine/Threonine kinases with the potential to phosphorylate the transcription factor Su(H) in crystal cell homeostasis. Su(H) is central to Notch signal transduction, and its inhibition by phosphorylation impedes crystal cell formation. Overall, nearly twenty percent of all Drosophila Serine/Threonine kinases were studied in two assays, global and hemocyte-specific overexpression and downregulation, respectively. Unexpectedly, the majority of kinases influenced crystal cell numbers, albeit only a few were related to hematopoiesis so far. Four kinases appeared essential for crystal cell formation, whereas most kinases restrained crystal cell development. This group comprises all kinase classes, indicative of the complex regulatory network underlying blood cell homeostasis. The rather indiscriminative response we observed opens the possibility that blood cells measure their overall phospho-status as a proxy for stress-signals, and activate an adaptive immune response accordingly.

Mapping single-cell developmental potential in health and disease with interpretable deep learning

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cell fate in developmental systems. However, identifying the molecular hallmarks of potency - the capacity of a cell to differentiate into other cell types - has remained challenging. Here, we introduce CytoTRACE 2, an interpretable deep learning framework for characterizing potency and differentiation states on an absolute scale from scRNA-seq data. Across 31 human and mouse scRNA-seq datasets encompassing 28 tissue types, CytoTRACE 2 outperformed existing methods for recovering experimentally determined potency levels and differentiation states covering the entire range of cellular ontogeny. Moreover, it reconstructed the temporal hierarchy of mouse embryogenesis across 62 timepoints; identified pan-tissue expression programs that discriminate major potency levels; and facilitated discovery of cellular phenotypes in cancer linked to survival and immunotherapy resistance. Our results illuminate a fundamental feature of cell biology and provide a broadly applicable platform for delineating single-cell differentiation landscapes in health and disease.

Collective effects of cell cleavage dynamics

A conserved process of early embryonic development in metazoans is the reductive cell divisions following oocyte fertilization, termed cell cleavages. Cell cleavage cycles usually start synchronously, lengthen differentially between the embryonic cells becoming asynchronous, and cease before major morphogenetic events, such as germ layer formation and gastrulation. Despite exhibiting species-specific characteristics, the regulation of cell cleavage dynamics comes down to common controllers acting mostly at the single cell/nucleus level, such as nucleus-to-cytoplasmic ratio and zygotic genome activation. Remarkably, recent work has linked cell cleavage dynamics to the emergence of collective behavior during embryogenesis, including pattern formation and changes in embryo-scale mechanics, raising the question how single-cell controllers coordinate embryo-scale processes. In this review, we summarize studies across species where an association between cell cleavages and collective behavior was made, discuss the underlying mechanisms, and propose that cell-to-cell variability in cell cleavage dynamics can serve as a mechanism of long-range coordination in developing embryos.

Drug-resilient Cancer Cell Phenotype Is Acquired via Polyploidization Associated with Early Stress Response Coupled to HIF2α Transcriptional Regulation

Therapeutic resistance and recurrence remain core challenges in cancer therapy. How therapy resistance arises is currently not fully understood with tumors surviving via multiple alternative routes. Here, we demonstrate that a subset of cancer cells survives therapeutic stress by entering a transient state characterized by whole-genome doubling. At the onset of the polyploidization program, we identified an upregulation of key transcriptional regulators, including the early stress-response protein AP-1 and normoxic stabilization of HIF2α. We found altered chromatin accessibility, ablated expression of retinoblastoma protein (RB1), and enrichment of AP-1 motif accessibility. We demonstrate that AP-1 and HIF2α regulate a therapy resilient and survivor phenotype in cancer cells. Consistent with this, genetic or pharmacologic targeting of AP-1 and HIF2α reduced the number of surviving cells following chemotherapy treatment. The role of AP-1 and HIF2α in stress response by polyploidy suggests a novel avenue for tackling chemotherapy-induced resistance in cancer.

SIGNIFICANCE

In response to cisplatin treatment, some surviving cancer cells undergo whole-genome duplications without mitosis, which represents a mechanism of drug resistance. This study presents mechanistic data to implicate AP-1 and HIF2α signaling in the formation of this surviving cell phenotype. The results open a new avenue for targeting drug-resistant cells.

A yeast cell cycle pulse generator model shows consistency with multiple oscillatory and checkpoint mutant datasets

Modeling biological systems holds great promise for speeding up the rate of discovery in systems biology by predicting experimental outcomes and suggesting targeted interventions. However, this process is dogged by an identifiability issue, in which network models and their parameters are not sufficiently constrained by coarse and noisy data to ensure unique solutions. In this work, we evaluated the capability of a simplified yeast cell-cycle network model to reproduce multiple observed transcriptomic behaviors under genomic mutations. We matched time-series data from both cycling and checkpoint arrested cells to model predictions using an asynchronous multi-level Boolean approach. We showed that this single network model, despite its simplicity, is capable of exhibiting dynamical behavior similar to the datasets in most cases, and we demonstrated the drop in severity of the identifiability issue that results from matching multiple datasets.

A cell-based model for size control in the multiple fission alga Chlamydomonas reinhardtii

Understanding how population-size homeostasis emerges from stochastic individual cell behaviors remains a challenge in biology.1,2,3,4,5,6,7 The unicellular green alga Chlamydomonas reinhardtii (Chlamydomonas) proliferates using a multiple fission cell cycle, where a prolonged G1 phase is followed by n rounds of alternating division cycles (S/M) to produce 2n daughters. A "Commitment" sizer in mid-G1 phase ensures sufficient cell growth before completing the cell cycle. A mitotic sizer couples mother-cell size to division number (n) such that daughter size distributions are uniform regardless of mother size distributions. Although daughter size distributions were highly robust to altered growth conditions, ∼40% of daughter cells fell outside of the 2-fold range expected from a "perfect" multiple fission sizer.7,8 A simple intuitive power law model with stochastic noise failed to reproduce individual division behaviors of tracked single cells. Through additional iterative modeling, we identified an alternative modified threshold (MT) model, where cells need to cross a threshold greater than 2-fold their median starting size to become division-competent (i.e., Committed), after which their behaviors followed a power law model. The Commitment versus mitotic size threshold uncoupling in the MT model was likely a key pre-adaptation in the evolution of volvocine algal multicellularity. A similar experimental approach was used in size mutants mat3/rbr and dp1 that are, respectively, missing repressor or activator subunits of the retinoblastoma tumor suppressor complex (RBC). Both mutants showed altered relationships between Commitment and mitotic sizer, suggesting that RBC functions to decouple the two sizers.

Temporal gene expression changes and affected pathways in neurodevelopment of a mouse model of Smith-Lemli-Opitz syndrome

Smith-Lemli-Opitz syndrome is an autosomal recessive disorder that arises from mutations in the gene DHCR7, which encodes the terminal enzyme of cholesterol biosynthesis, leading to decreased production of cholesterol and accumulation of the cholesterol precursor, 7-dehydrocholesterol, and its oxysterol metabolites. The disorder displays a wide range of neurodevelopmental defects, intellectual disability, and behavioral problems. However, an in-depth study on the temporal changes of gene expression in the developing brains of SLOS mice has not been done before. In this work, we carried out the transcriptomic analysis of whole brains from WT and Dhcr7-KO mice at four-time points through postnatal day 0. First, we observed the expected downregulation of the Dhcr7 gene in the Dhcr7-KO mouse model, as well as gene expression changes of several other genes involved in cholesterol biosynthesis throughout all time points. Pathway and GO term enrichment analyses revealed affected signaling pathways and biological processes that were shared amongst time points and unique to individual time points. Specifically, the pathways important for embryonic development, including Hippo, Wnt, and TGF-β signaling pathways are the most significantly affected at the earliest time point, E12.5. Additionally, neurogenesis-related GO terms were enriched in earlier time points, consistent with the timing of development. Conversely, pathways related to synaptogenesis, which occurs later in development compared to neurogenesis, are significantly affected at the later time points, E16.5 and PND0, including the cholinergic, glutamatergic, and GABAergic synapses. The impact of these transcriptomic changes and enriched pathways is discussed in the context of known biological phenotypes of SLOS.

High-throughput sequencing reveals hub genes for human early embryonic development arrest in vitro fertilization: a pilot study

Many clinical studies have shown that embryos of in vitro fertilization (IVF) are often prone to developmental arrest, which leads to recurrent failure of IVF treatment. Early embryonic arrest has always been an urgent clinical problem in assisted reproduction centers. However, the molecular mechanisms underlying early embryonic development arrest remain largely unknown. The objective of this study is to investigate potential candidate hub genes and key signaling pathways involved in early stages of embryonic development. RNA-seq analysis was performed on normal and arrest embryos to study the changes of gene expression during early embryonic development. A total of 520 genes exhibiting differential expression were identified, with 174 genes being upregulated and 346 genes being downregulated. Upregulated genes show enrichment in biosynthesis, cellular proliferation and differentiation, and epigenetic regulation. While downregulated genes exhibit enrichment in transcriptional activity, epigenetic regulation, cell cycle progression, cellular proliferation and ubiquitination. The STRING (search tool for the retravel of interacting genes/proteins) database was utilized to analyze protein-protein interactions among these genes, aiming to enhance comprehension of the potential role of these differentially expressed genes (DEGs). A total of 22 hub genes (highly connected genes) were identified among the DEGs using Cytoscape software. Of these, ERBB2 and VEGFA were upregulated, while the remaining 20 genes (CCNB1, CCNA2, DICER1, NOTCH1, UBE2B, UBE2N, PRMT5, UBE2D1, MAPK3, SOX9, UBE2C, UB2D2, EGF, ACTB, UBA52, SHH, KRAS, UBE2E1, ADAM17 and BRCA2) were downregulated. These hub genes are associated with crucial biological processes such as ubiquitination, cellular senescence, cell proliferation and differentiation, and cell cycle. Among these hub genes, CCNA2 and CCNB1 may be involved in controlling cell cycle, which are critical process in early embryonic development.

Bensulide-induced oxidative stress causes developmental defects of cardiovascular system and liver in zebrafish (Danio rerio)

Bensulide is an organophosphate herbicide commonly used in agricultural crops; however, no studies have reported on its toxic effects in the embryonic development of vertebrates, particularly gene expression level and cellular response. Therefore, to identify developmental toxicity, zebrafish eggs 8 h post-fertilization (hpf) were exposed to bensulide concentrations of up to 3 mg/L. The results indicated that exposure to 3 mg/L bensulide inhibited the hatching of all eggs and decreased the size of the body, eyes, and inner ear. There were demonstrated effects observed in the cardiovascular system and liver caused by bensulide in fli1:eGFP and L-fabp:dsRed transgenic zebrafish models, respectively. Following exposure to 3 mg/L bensulide, normal heart development, including cardiac looping, was disrupted and the heart rate of 96 hpf zebrafish larvae decreased to 16.37%. Development of the liver, the main detoxification organ, was also inhibited by bensulide, and after exposure to 3 mg/L bensulide its size reduced to 41.98%. Additionally, exposure to bensulide resulted in inhibition of antioxidant enzyme expression and an increase in ROS levels by up to 238.29%. Collectively, we identified various biological responses associated with the toxicity of bensulide, which led to various organ malformations and cytotoxic effects in zebrafish.

Vasa, a regulator of localized mRNA translation on the spindle

Localized mRNA translation is a biological process that allows mRNA to be translated on-site, which is proposed to provide fine control in protein regulation, both spatially and temporally within a cell. We recently reported that Vasa, an RNA-helicase, is a promising factor that appears to regulate this process on the spindle during the embryonic development of the sea urchin, yet the detailed roles and functional mechanisms of Vasa in this process are still largely unknown. In this review article, to elucidate these remaining questions, we first summarize the prior knowledge and our recent findings in the area of Vasa research and further discuss how Vasa may function in localized mRNA translation, contributing to efficient protein regulation during rapid embryogenesis and cancer cell regulation.

The nuclear-to-cytoplasmic ratio drives cellularization in the close animal relative Sphaeroforma arctica

The ratio of nuclear content to cytoplasmic volume (N/C ratio) is a key regulator driving the maternal-to-zygotic transition in most animal embryos. Altering this ratio often impacts zygotic genome activation and deregulates the timing and outcome of embryogenesis.^1,2,3 Despite being ubiquitous across animals, little is known about when the N/C ratio evolved to control multicellular development. Such capacity either originated with the emergence of animal multicellularity or was co-opted from the mechanisms present in unicellular organisms.⁴ An effective strategy to tackle this question is to investigate the close relatives of animals exhibiting life cycles with transient multicellular stages.⁵ Among these are ichthyosporeans, a lineage of protists undergoing coenocytic development followed by cellularization and cell release.^6,7,8 During cellularization, a transient multicellular stage resembling animal epithelia is generated, offering a unique opportunity to examine whether the N/C ratio regulates multicellular development. Here, we use time-lapse microscopy to characterize how the N/C ratio affects the life cycle of the best-studied ichthyosporean model, Sphaeroforma arctica. We uncover that the last stages of cellularization coincide with a significant increase in the N/C ratio. Increasing the N/C ratio by reducing the coenocytic volume accelerates cellularization, whereas decreasing the N/C ratio by lowering the nuclear content halts it. Moreover, centrifugation and pharmacological inhibitor experiments suggest that the N/C ratio is locally sensed at the cortex and relies on phosphatase activity. Altogether, our results show that the N/C ratio drives cellularization in S. arctica, suggesting that its capacity to control multicellular development predates animal emergence.

The Nuclear-to-Cytoplasmic Ratio: Coupling DNA Content to Cell Size, Cell Cycle, and Biosynthetic Capacity

Though cell size varies between different cells and across species, the nuclear-to-cytoplasmic (N/C) ratio is largely maintained across species and within cell types. A cell maintains a relatively constant N/C ratio by coupling DNA content, nuclear size, and cell size. We explore how cells couple cell division and growth to DNA content. In some cases, cells use DNA as a molecular yardstick to control the availability of cell cycle regulators. In other cases, DNA sets a limit for biosynthetic capacity. Developmentally programmed variations in the N/C ratio for a given cell type suggest that a specific N/C ratio is required to respond to given physiological demands. Recent observations connecting decreased N/C ratios with cellular senescence indicate that maintaining the proper N/C ratio is essential for proper cellular functioning. Together, these findings suggest a causative, not simply correlative, role for the N/C ratio in regulating cell growth and cell cycle progression.

Manipulating the nature of embryonic mitotic waves

Early embryogenesis is characterized by rapid and synchronous cleavage divisions, which are often controlled by wave-like patterns of Cdk1 activity. Two mechanisms have been proposed for mitotic waves: sweep and trigger waves.1,2 The two mechanisms give rise to different wave speeds, dependencies on physical and molecular parameters, and spatial profiles of Cdk1 activity: upward sweeping gradients versus traveling wavefronts. Both mechanisms hinge on the transient bistability governing the cell cycle and are differentiated by the speed of the cell-cycle progression: sweep and trigger waves arise for rapid and slow drives, respectively. Here, using quantitative imaging of Cdk1 activity and theory, we illustrate that sweep waves are the dominant mechanism in Drosophila embryos and test two fundamental predictions on the transition from sweep to trigger waves. We demonstrate that sweep waves can be turned into trigger waves if the cell cycle is slowed down genetically or if significant delays in the cell-cycle progression are introduced across the embryo by altering nuclear density. Our genetic experiments demonstrate that Polo kinase is a major rate-limiting regulator of the blastoderm divisions, and genetic perturbations reducing its activity can induce the transition from sweep to trigger waves. Furthermore, we show that changes in temperature cause an essentially uniform slowdown of interphase and mitosis. That results in sweep waves being observed across a wide temperature range despite the cell-cycle durations being significantly different. Collectively, our combination of theory and experiments elucidates the nature of mitotic waves in Drosophila embryogenesis, their control mechanisms, and their mutual transitions.

Quantitative models for building and growing fated small cell networks

Small cell clusters exhibit numerous phenomena typically associated with complex systems, such as division of labour and programmed cell death. A conserved class of such clusters occurs during oogenesis in the form of germline cysts that give rise to oocytes. Germline cysts form through cell divisions with incomplete cytokinesis, leaving cells intimately connected through intercellular bridges that facilitate cyst generation, cell fate determination and collective growth dynamics. Using the well-characterized Drosophila melanogaster female germline cyst as a foundation, we present mathematical models rooted in the dynamics of cell cycle proteins and their interactions to explain the generation of germline cell lineage trees (CLTs) and highlight the diversity of observed CLT sizes and topologies across species. We analyse competing models of symmetry breaking in CLTs to rationalize the observed dynamics and robustness of oocyte fate specification, and highlight remaining gaps in knowledge. We also explore how CLT topology affects cell cycle dynamics and synchronization and highlight mechanisms of intercellular coupling that underlie the observed collective growth patterns during oogenesis. Throughout, we point to similarities across organisms that warrant further investigation and comment on the extent to which experimental and theoretical findings made in model systems extend to other species.