Tuning cell differentiation patterns and single cell dynamics by regulating proteins' functionalities in a toggle switch

Identifying toggle genes from transcriptome-wide scatter: A new perspective for biological regulation

The study of gene expression variability, especially for cancer and cell differentiation studies, has become important. Here, we investigate transcriptome-wide scatter of 23 cell types and conditions across different levels of biological complexity. We focused on genes that act like toggle switches between pairwise replicates of the same cell type, i.e. genes expressed in one replicate and not expressed in the other, sometimes also referred as ON/OFF genes. The proportion of these toggle genes dramatically increases from unicellular to multicellular organization, especially for development and cancer cells. A relevant portion of toggle switches are non-coding genes: in unicellular systems the most represented classes are tRNA and rRNA, while multicellular systems more frequently show lncRNA, sncRNA and pseudogenes. Notably, disease associated microRNAs (miRNAs), pseudogenes and numerous uncharacterized transcripts are present in both development and cancer cells. On top of the known intrinsic and extrinsic factors, our work indicates toggle genes as a novel collective component creating transcriptome-wide variability. This requires further investigation for elucidating both evolutionary and disease processes.

Efficiency and Robustness of Processes Driven by Nucleoid Exclusion in Escherichia coli

The internal spatial organization of prokaryotic organisms, including Escherichia coli, is essential for the proper functioning of processes such as cell division. One source of this organization in E. coli is the nucleoid, which causes the exclusion of macromolecules - e.g. protein aggregates and the chemotaxis network - from midcell. Similarly, following DNA replication, the nucleoid(s) assist in placing the Z-ring at midcell. These processes need to be efficient in optimal conditions and robust to suboptimal conditions. After reviewing recent findings on these topics, we make use of past data to study the efficiency of the spatial constraining of Z-rings, chemotaxis networks, and protein aggregates, as a function of the nucleoid(s) morphology. Also, we compare the robustness of these processes to nonoptimal temperatures. We show that Z-rings, Tsr clusters, and protein aggregates have temperature-dependent spatial distributions along the major cell axis that are consistent with the nucleoid(s) morphology and the volume-exclusion phenomenon. Surprisingly, the consequences of the changes in nucleoid size with temperature are most visible in the kurtosis of these spatial distributions, in that it has a statistically significant linear correlation with the mean nucleoid length and, in the case of Z-rings, with the distance between nucleoids prior to cell division. Interestingly, we also find a negative, statistically significant linear correlation between the efficiency of these processes at the optimal condition and their robustness to suboptimal conditions, suggesting a trade-off between these traits.

In silico analysis of division times of Escherichia coli populations as a function of the partitioning scheme of non-functional proteins

Recent evidence suggests that cells employ functionally asymmetric partitioning schemes in division to cope with aging. We explore various schemes in silico, with a stochastic model of Escherichia coli that includes gene expression, non-functional proteins generation, aggregation and polar retention, and molecule partitioning in division. The model is implemented in SGNS2, which allows stochastic, multi-delayed reactions within hierarchical, transient, interlinked compartments. After setting parameter values of non-functional proteins’ generation and effects that reproduce realistic intracellular and population dynamics, we investigate how the spatial organization of non-functional proteins affects mean division times of cell populations in lineages and, thus, mean cell numbers over time. We find that division times decrease for increasingly asymmetric partitioning. Also, increasing the clustering of non-functional proteins decreases division times. Increasing the bias in polar segregation further decreases division times, particularly if the bias favors the older pole and aggregates’ polar retention is robust. Finally, we show that the non-energy consuming retention of inherited non-functional proteins at the older pole via nucleoid occlusion is a source of functional asymmetries and, thus, is advantageous. Our results suggest that the mechanisms of intracellular organization of non-functional proteins, including clustering and polar retention, affect the vitality of E. coli populations.

Cell-to-cell diversity in protein levels of a gene driven by a tetracycline inducible promoter

BACKGROUND

Gene expression in Escherichia coli is regulated by several mechanisms. We measured in single cells the expression level of a single copy gene coding for green fluorescent protein (GFP), integrated into the genome and driven by a tetracycline inducible promoter, for varying induction strengths. Also, we measured the transcriptional activity of a tetracycline inducible promoter controlling the transcription of a RNA with 96 binding sites for MS2-GFP.

RESULTS

The distribution of GFP levels in single cells is found to change significantly as induction reaches high levels, causing the Fano factor of the cells' protein levels to increase with mean level, beyond what would be expected from a Poisson-like process of RNA transcription. In agreement, the Fano factor of the cells' number of RNA molecules target for MS2-GFP follows a similar trend. The results provide evidence that the dynamics of the promoter complex formation, namely, the variability in its duration from one transcription event to the next, explains the change in the distribution of expression levels in the cell population with induction strength.

CONCLUSIONS

The results suggest that the open complex formation of the tetracycline inducible promoter, in the regime of strong induction, affects significantly the dynamics of RNA production due to the variability of its duration from one event to the next.

Information propagation within the Genetic Network of Saccharomyces cerevisiae

Background

A gene network's capacity to process information, so as to bind past events to future actions, depends on its structure and logic. From previous and new microarray measurements in Saccharomyces cerevisiae following gene deletions and overexpressions, we identify a core gene regulatory network (GRN) of functional interactions between 328 genes and the transfer functions of each gene. Inferred connections are verified by gene enrichment.

Results

We find that this core network has a generalized clustering coefficient that is much higher than chance. The inferred Boolean transfer functions have a mean p-bias of 0.41, and thus similar amounts of activation and repression interactions. However, the distribution of p-biases differs significantly from what is expected by chance that, along with the high mean connectivity, is found to cause the core GRN of S. cerevisiae's to have an overall sensitivity similar to critical Boolean networks. In agreement, we find that the amount of information propagated between nodes in finite time series is much higher in the inferred core GRN of S. cerevisiae than what is expected by chance.

Conclusions

We suggest that S. cerevisiae is likely to have evolved a core GRN with enhanced information propagation among its genes.