Deep learning-driven imaging of cell division and cell growth across an entire eukaryotic life cycle.

The life cycle of biomedical and agriculturally relevant eukaryotic microorganisms involves complex transitions between proliferative and non-proliferative states such as dormancy, mating, meiosis, and cell division. New drugs, pesticides, and vaccines can be created by targeting specific life cycle stages of parasites and pathogens. However, defining the structure of a microbial life cycle often relies on partial observations that are theoretically assembled in an ideal life cycle path. To create a more quantitative approach to studying complete eukaryotic life cycles, we generated a deep learning-driven imaging framework to track microorganisms across sexually reproducing generations. Our approach combines microfluidic culturing, life cycle stage-specific segmentation of microscopy images using convolutional neural networks, and a novel cell tracking algorithm, FIEST, based on enhancing the overlap of single cell masks in consecutive images through deep learning video frame interpolation. As proof of principle, we used this approach to quantitatively image and compare cell growth and cell cycle regulation across the sexual life cycle of Saccharomyces cerevisiae. We developed a fluorescent reporter system based on a fluorescently labeled Whi5 protein, the yeast analog of mammalian Rb, and a new High-Cdk1 activity sensor, LiCHI, designed to report during DNA replication, mitosis, meiotic homologous recombination, meiosis I, and meiosis II. We found that cell growth preceded the exit from non-proliferative states such as mitotic G1, pre-meiotic G1, and the G0 spore state during germination. A decrease in the total cell concentration of Whi5 characterized the exit from non-proliferative states, which is consistent with a Whi5 dilution model. The nuclear accumulation of Whi5 was developmentally regulated, being at its highest during meiotic exit and spore formation. The temporal coordination of cell division and growth was not significantly different across three sexually reproducing generations. Our framework could be used to quantitatively characterize other single-cell eukaryotic life cycles that remain incompletely described. An off-the-shelf user interface Yeastvision provides free access to our image processing and single-cell tracking algorithms.

Diploid-dominant life cycles characterize the early evolution of Fungi.

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

Competition for DNA binding between paralogous transcription factors determines their genomic occupancy and regulatory functions.

Most eukaryotic transcription factors (TFs) are part of large protein families, with members of the same family (i.e., paralogous TFs) recognizing similar DNA-binding motifs but performing different regulatory functions. Many TF paralogs are coexpressed in the cell and thus can compete for target sites across the genome. However, this competition is rarely taken into account when studying the in vivo binding patterns of eukaryotic TFs. Here, we show that direct competition for DNA binding between TF paralogs is a major determinant of their genomic binding patterns. Using yeast proteins Cbf1 and Pho4 as our model system, we designed a high-throughput quantitative assay to capture the genomic binding profiles of competing TFs in a cell-free system. Our data show that Cbf1 and Pho4 greatly influence each other's occupancy by competing for their common putative genomic binding sites. The competition is different at different genomic sites, as dictated by the TFs' expression levels and their divergence in DNA-binding specificity and affinity. Analyses of ChIP-seq data show that the biophysical rules that dictate the competitive TF binding patterns in vitro are also followed in vivo, in the complex cellular environment. Furthermore, the Cbf1-Pho4 competition for genomic sites, as characterized in vitro using our new assay, plays a critical role in the specific activation of their target genes in the cell. Overall, our study highlights the importance of direct TF-TF competition for genomic binding and gene regulation by TF paralogs, and proposes an approach for studying this competition in a quantitative and high-throughput manner.

Chytrid fungi.

Medina and Buchler provide an introduction to chytrid fungi, an early diverging fungal lineage exhibiting characteristics found in both animals and fungi.

Genetic transformation of Spizellomyces punctatus, a resource for studying chytrid biology and evolutionary cell biology.

Chytrids are early-diverging fungi that share features with animals that have been lost in most other fungi. They hold promise as a system to study fungal and animal evolution, but we lack genetic tools for hypothesis testing. Here, we generated transgenic lines of the chytrid Spizellomyces punctatus, and used fluorescence microscopy to explore chytrid cell biology and development during its life cycle. We show that the chytrid undergoes multiple rounds of synchronous nuclear division, followed by cellularization, to create and release many daughter 'zoospores'. The zoospores, akin to animal cells, crawl using actin-mediated cell migration. After forming a cell wall, polymerized actin reorganizes into fungal-like cortical patches and cables that extend into hyphal-like structures. Actin perinuclear shells form each cell cycle and polygonal territories emerge during cellularization. This work makes Spizellomyces a genetically tractable model for comparative cell biology and understanding the evolution of fungi and early eukaryotes.

Recognition of Histone Crotonylation by Taf14 Links Metabolic State to Gene Expression.

Metabolic signaling to chromatin often underlies how adaptive transcriptional responses are controlled. While intermediary metabolites serve as co-factors for histone-modifying enzymes during metabolic flux, how these modifications contribute to transcriptional responses is poorly understood. Here, we utilize the highly synchronized yeast metabolic cycle (YMC) and find that fatty acid ß-oxidation genes are periodically expressed coincident with the ß-oxidation byproduct histone crotonylation. Specifically, we found that H3K9 crotonylation peaks when H3K9 acetylation declines and energy resources become limited. During this metabolic state, pro-growth gene expression is dampened; however, mutation of the Taf14 YEATS domain, a H3K9 crotonylation reader, results in de-repression of these genes. Conversely, exogenous addition of crotonic acid results in increased histone crotonylation, constitutive repression of pro-growth genes, and disrupted YMC oscillations. Together, our findings expose an unexpected link between metabolic flux and transcription and demonstrate that histone crotonylation and Taf14 participate in the repression of energy-demanding gene expression.

Epigenetic switching as a strategy for quick adaptation while attenuating biochemical noise.

Epigenetic switches are bistable, molecular systems built from self-reinforcing feedback loops that can spontaneously switch between heritable phenotypes in the absence of DNA mutation. It has been hypothesized that epigenetic switches first evolved as a mechanism of bet-hedging and adaptation, but the evolutionary trajectories and conditions by which an epigenetic switch can outcompete adaptation through genetic mutation remain unknown. Here, we used computer simulations to evolve a mechanistic, biophysical model of a self-activating genetic circuit, which can both adapt genetically through mutation and exhibit epigenetic switching. We evolved these genetic circuits under a fluctuating environment that alternatively selected for low and high protein expression levels. In all tested conditions, the population first evolved by genetic mutation towards a region of genotypes where genetic adaptation can occur faster after each environmental transition. Once in this region, the self-activating genetic circuit can exhibit epigenetic switching, which starts competing with genetic adaptation. We show a trade-off between either minimizing the adaptation time or increasing the robustness of the phenotype to biochemical noise. Epigenetic switching was superior in a fast fluctuating environment because it adapted faster than genetic mutation after an environmental transition, while still attenuating the effect of biochemical noise on the phenotype. Conversely, genetic adaptation was favored in a slowly fluctuating environment because it maximized the phenotypic robustness to biochemical noise during the constant environment between transitions, even if this resulted in slower adaptation. This simple trade-off predicts the conditions and trajectories under which an epigenetic switch evolved to outcompete genetic adaptation, shedding light on possible mechanisms by which bet-hedging strategies might emerge and persist in natural populations.

Evolutionary innovation, fungal cell biology, and the lateral gene transfer of a viral KilA-N domain.

Fungi are found in diverse ecological niches as primary decomposers, mutualists, or parasites of plants and animals. Although animals and fungi share a common ancestor, fungi dramatically diversified their life cycle, cell biology, and metabolism as they evolved and colonized new niches. This review focuses on a family of fungal transcription factors (Swi4/Mbp1, APSES, Xbp1, Bqt4) derived from the lateral gene transfer of a KilA-N domain commonly found in prokaryotic and eukaryotic DNA viruses. These virus-derived fungal regulators play central roles in cell cycle, morphogenesis, sexual differentiation, and quiescence. We consider the possible origins of KilA-N and how this viral DNA binding domain came to be intimately associated with fungal processes.

Exact and efficient hybrid Monte Carlo algorithm for accelerated Bayesian inference of gene expression models from snapshots of single-cell transcripts.

Single cells exhibit a significant amount of variability in transcript levels, which arises from slow, stochastic transitions between gene expression states. Elucidating the nature of these states and understanding how transition rates are affected by different regulatory mechanisms require state-of-the-art methods to infer underlying models of gene expression from single cell data. A Bayesian approach to statistical inference is the most suitable method for model selection and uncertainty quantification of kinetic parameters using small data sets. However, this approach is impractical because current algorithms are too slow to handle typical models of gene expression. To solve this problem, we first show that time-dependent mRNA distributions of discrete-state models of gene expression are dynamic Poisson mixtures, whose mixing kernels are characterized by a piecewise deterministic Markov process. We combined this analytical result with a kinetic Monte Carlo algorithm to create a hybrid numerical method that accelerates the calculation of time-dependent mRNA distributions by 1000-fold compared to current methods. We then integrated the hybrid algorithm into an existing Monte Carlo sampler to estimate the Bayesian posterior distribution of many different, competing models in a reasonable amount of time. We demonstrate that kinetic parameters can be reasonably constrained for modestly sampled data sets if the model is known a priori. If there are many competing models, Bayesian evidence can rigorously quantify the likelihood of a model relative to other models from the data. We demonstrate that Bayesian evidence selects the true model and outperforms approximate metrics typically used for model selection.

Enhancer Histone Acetylation Modulates Transcriptional Bursting Dynamics of Neuronal Activity-Inducible Genes.

Neuronal activity-inducible gene transcription correlates with rapid and transient increases in histone acetylation at promoters and enhancers of activity-regulated genes. Exactly how histone acetylation modulates transcription of these genes has remained unknown. We used single-cell in situ transcriptional analysis to show that Fos and Npas4 are transcribed in stochastic bursts in mouse neurons and that membrane depolarization increases mRNA expression by increasing burst frequency. We then expressed dCas9-p300 or dCas9-HDAC8 fusion proteins to mimic or block activity-induced histone acetylation locally at enhancers. Adding histone acetylation increased Fos transcription by prolonging burst duration and resulted in higher Fos protein levels and an elevation of resting membrane potential. Inhibiting histone acetylation reduced Fos transcription by reducing burst frequency and impaired experience-dependent Fos protein induction in the hippocampus in vivo. Thus, activity-inducible histone acetylation tunes the transcriptional dynamics of experience-regulated genes to affect selective changes in neuronal gene expression and cellular function.

Daily humidity oscillation regulates the circadian clock to influence plant physiology.

Early circadian studies in plants by de Mairan and de Candolle alluded to a regulation of circadian clocks by humidity. However, this regulation has not been described in detail, nor has its influence on physiology been demonstrated. Here we report that, under constant light, circadian humidity oscillation can entrain the plant circadian clock to a period of 24 h probably through the induction of clock genes such as CIRCADIAN CLOCK ASSOCIATED 1. Under simulated natural light and humidity cycles, humidity oscillation increases the amplitude of the circadian clock and further improves plant fitness-related traits. In addition, humidity oscillation enhances effector-triggered immunity at night possibly to counter increased pathogen virulence under high humidity. These results indicate that the humidity oscillation regulates specific circadian outputs besides those co-regulated with the light-dark cycle.

Efficient analysis of stochastic gene dynamics in the non-adiabatic regime using piecewise deterministic Markov processes.

Single-cell experiments show that gene expression is stochastic and bursty, a feature that can emerge from slow switching between promoter states with different activities. In addition to slow chromatin and/or DNA looping dynamics, one source of long-lived promoter states is the slow binding and unbinding kinetics of transcription factors to promoters, i.e. the non-adiabatic binding regime. Here, we introduce a simple analytical framework, known as a piecewise deterministic Markov process (PDMP), that accurately describes the stochastic dynamics of gene expression in the non-adiabatic regime. We illustrate the utility of the PDMP on a non-trivial dynamical system by analysing the properties of a titration-based oscillator in the non-adiabatic limit. We first show how to transform the underlying chemical master equation into a PDMP where the slow transitions between promoter states are stochastic, but whose rates depend upon the faster deterministic dynamics of the transcription factors regulated by these promoters. We show that the PDMP accurately describes the observed periods of stochastic cycles in activator and repressor-based titration oscillators. We then generalize our PDMP analysis to more complicated versions of titration-based oscillators to explain how multiple binding sites lengthen the period and improve coherence. Last, we show how noise-induced oscillation previously observed in a titration-based oscillator arises from non-adiabatic and discrete binding events at the promoter site.

The evolution of a G1/S transcriptional network in yeasts.

The G1-to-S cell cycle transition is promoted by the periodic expression of a large set of genes. In Saccharomyces cerevisiae G1/S gene expression is regulated by two transcription factor (TF) complexes, the MBF and SBF, which bind to specific DNA sequences, the MCB and SCB, respectively. Despite extensive research little is known regarding the evolution of the G1/S transcription regulation including the co-evolution of the DNA binding domains with their respective DNA binding sequences. We have recently examined the co-evolution of the G1/S TF specificity through the systematic generation and examination of chimeric Mbp1/Swi4 TFs containing different orthologue DNA binding domains in S. cerevisiae (Hendler et al. in PLoS Genet 13:e1006778. doi: 10.1371/journal.pgen.1006778 , 2017). Here, we review the co-evolution of G1/S transcriptional network and discuss the evolutionary dynamics and specificity of the MBF-MCB and SBF-SCB interactions in different fungal species.

BayFish: Bayesian inference of transcription dynamics from population snapshots of single-molecule RNA FISH in single cells.

Single-molecule RNA fluorescence in situ hybridization (smFISH) provides unparalleled resolution in the measurement of the abundance and localization of nascent and mature RNA transcripts in fixed, single cells. We developed a computational pipeline (BayFish) to infer the kinetic parameters of gene expression from smFISH data at multiple time points after gene induction. Given an underlying model of gene expression, BayFish uses a Monte Carlo method to estimate the Bayesian posterior probability of the model parameters and quantify the parameter uncertainty given the observed smFISH data. We tested BayFish on synthetic data and smFISH measurements of the neuronal activity-inducible gene Npas4 in primary neurons.

Gene duplication and co-evolution of G1/S transcription factor specificity in fungi are essential for optimizing cell fitness.

Transcriptional regulatory networks play a central role in optimizing cell survival. How DNA binding domains and cis-regulatory DNA binding sequences have co-evolved to allow the expansion of transcriptional networks and how this contributes to cellular fitness remains unclear. Here we experimentally explore how the complex G1/S transcriptional network evolved in the budding yeast Saccharomyces cerevisiae by examining different chimeric transcription factor (TF) complexes. Over 200 G1/S genes are regulated by either one of the two TF complexes, SBF and MBF, which bind to specific DNA binding sequences, SCB and MCB, respectively. The difference in size and complexity of the G1/S transcriptional network across yeast species makes it well suited to investigate how TF paralogs (SBF and MBF) and DNA binding sequences (SCB and MCB) co-evolved after gene duplication to rewire and expand the network of G1/S target genes. Our data suggests that whilst SBF is the likely ancestral regulatory complex, the ancestral DNA binding element is more MCB-like. G1/S network expansion took place by both cis- and trans- co-evolutionary changes in closely related but distinct regulatory sequences. Replacement of the endogenous SBF DNA-binding domain (DBD) with that from more distantly related fungi leads to a contraction of the SBF-regulated G1/S network in budding yeast, which also correlates with increased defects in cell growth, cell size, and proliferation.

Conservation and divergence of C-terminal domain structure in the retinoblastoma protein family.

The retinoblastoma protein (Rb) and the homologous pocket proteins p107 and p130 negatively regulate cell proliferation by binding and inhibiting members of the E2F transcription factor family. The structural features that distinguish Rb from other pocket proteins have been unclear but are critical for understanding their functional diversity and determining why Rb has unique tumor suppressor activities. We describe here important differences in how the Rb and p107 C-terminal domains (CTDs) associate with the coiled-coil and marked-box domains (CMs) of E2Fs. We find that although CTD-CM binding is conserved across protein families, Rb and p107 CTDs show clear preferences for different E2Fs. A crystal structure of the p107 CTD bound to E2F5 and its dimer partner DP1 reveals the molecular basis for pocket protein-E2F binding specificity and how cyclin-dependent kinases differentially regulate pocket proteins through CTD phosphorylation. Our structural and biochemical data together with phylogenetic analyses of Rb and E2F proteins support the conclusion that Rb evolved specific structural motifs that confer its unique capacity to bind with high affinity those E2Fs that are the most potent activators of the cell cycle.

Long-term growth data of Escherichia coli at a single-cell level.

Long-term, single-cell measurement of bacterial growth is extremely valuable information, particularly in the study of homeostatic aspects such as cell-size and growth rate control. Such measurement has recently become possible due to the development of microfluidic technology. Here we present data from single-cell measurements of Escherichia coli growth over 70 generations obtained for three different growth conditions. The data were recorded every minute, and contain time course data of cell length and fluorescent intensity of constitutively expressed yellow fluorescent protein.

Punctuated evolution and transitional hybrid network in an ancestral cell cycle of fungi.

Although cell cycle control is an ancient, conserved, and essential process, some core animal and fungal cell cycle regulators share no more sequence identity than non-homologous proteins. Here, we show that evolution along the fungal lineage was punctuated by the early acquisition and entrainment of the SBF transcription factor through horizontal gene transfer. Cell cycle evolution in the fungal ancestor then proceeded through a hybrid network containing both SBF and its ancestral animal counterpart E2F, which is still maintained in many basal fungi. We hypothesize that a virally-derived SBF may have initially hijacked cell cycle control by activating transcription via the cis-regulatory elements targeted by the ancestral cell cycle regulator E2F, much like extant viral oncogenes. Consistent with this hypothesis, we show that SBF can regulate promoters with E2F binding sites in budding yeast.

Cell cycle Start is coupled to entry into the yeast metabolic cycle across diverse strains and growth rates.

Cells have evolved oscillators with different frequencies to coordinate periodic processes. Here we studied the interaction of two oscillators, the cell division cycle (CDC) and the yeast metabolic cycle (YMC), in budding yeast. Previous work suggested that the CDC and YMC interact to separate high oxygen consumption (HOC) from DNA replication to prevent genetic damage. To test this hypothesis, we grew diverse strains in chemostat and measured DNA replication and oxygen consumption with high temporal resolution at different growth rates. Our data showed that HOC is not strictly separated from DNA replication; rather, cell cycle Start is coupled with the initiation of HOC and catabolism of storage carbohydrates. The logic of this YMC-CDC coupling may be to ensure that DNA replication and cell division occur only when sufficient cellular energy reserves have accumulated. Our results also uncovered a quantitative relationship between CDC period and YMC period across different strains. More generally, our approach shows how studies in genetically diverse strains efficiently identify robust phenotypes and steer the experimentalist away from strain-specific idiosyncrasies.

The Physarum polycephalum Genome Reveals Extensive Use of Prokaryotic Two-Component and Metazoan-Type Tyrosine Kinase Signaling.

Physarum polycephalum is a well-studied microbial eukaryote with unique experimental attributes relative to other experimental model organisms. It has a sophisticated life cycle with several distinct stages including amoebal, flagellated, and plasmodial cells. It is unusual in switching between open and closed mitosis according to specific life-cycle stages. Here we present the analysis of the genome of this enigmatic and important model organism and compare it with closely related species. The genome is littered with simple and complex repeats and the coding regions are frequently interrupted by introns with a mean size of 100 bases. Complemented with extensive transcriptome data, we define approximately 31,000 gene loci, providing unexpected insights into early eukaryote evolution. We describe extensive use of histidine kinase-based two-component systems and tyrosine kinase signaling, the presence of bacterial and plant type photoreceptors (phytochromes, cryptochrome, and phototropin) and of plant-type pentatricopeptide repeat proteins, as well as metabolic pathways, and a cell cycle control system typically found in more complex eukaryotes. Our analysis characterizes P. polycephalum as a prototypical eukaryote with features attributed to the last common ancestor of Amorphea, that is, the Amoebozoa and Opisthokonts. Specifically, the presence of tyrosine kinases in Acanthamoeba and Physarum as representatives of two distantly related subdivisions of Amoebozoa argues against the later emergence of tyrosine kinase signaling in the opisthokont lineage and also against the acquisition by horizontal gene transfer.