Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, the architecture of a promoter, defined as the number, strength and regulatory role of the operators that control the promoter, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect noise in gene expression in a systematic rather than case-by-case fashion. In this article, we make such a systematic investigation, based on a simple microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcription product from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can discriminate between different kinetic models of gene regulation.Comment: 35 pages, 6 figures, Submitte

The reductive glycine pathway allows autotrophic growth of Desulfovibrio desulfuricans

Supplementary informationis available for this paper athttps://doi.org/10.1038/s41467-020-18906-7Six CO2 fixation pathways are known to operate in photoautotrophic and chemoautotrophic microorganisms. Here, we describe chemolithoautotrophic growth of the sulphate-reducing bacterium Desulfovibrio desulfuricans (strain G11) with hydrogen and sulphate as energy substrates. Genomic, transcriptomic, proteomic and metabolomic analyses reveal that D. desulfuricans assimilates CO2 via the reductive glycine pathway, a seventh CO2 fixation pathway. In this pathway, CO2 is first reduced to formate, which is reduced and condensed with a second CO2 to generate glycine. Glycine is further reduced in D. desulfuricans by glycine reductase to acetyl-P, and then to acetyl-CoA, which is condensed with another CO2 to form pyruvate. Ammonia is involved in the operation of the pathway, which is reflected in the dependence of the autotrophic growth rate on the ammonia concentration. Our study demonstrates microbial autotrophic growth fully supported by this highly ATP-efficient CO2 fixation pathway.We acknowledge Änne-Michaelis and William Newell for assistance with the LC-MS forthe metabolomics experiments and Daniel Amador-Noguez for access to the LC-MS usedfor13C intracellular metabolomic analysis. We thank Ines Cardoso Pereira and John vander Oost for critically reading the manuscript. This research was funded by the Neth-erlands Organisation for Scientific Research (NWO) through SIAM Gravitation Grant024.002.002 and the Innovation Program Microbiology (WUR), NJC acknowledgesfunding from NWO through a Rubicon Grant (019.163LW.035) and a Veni Grant(VI.Veni.192.156).info:eu-repo/semantics/publishedVersio

Evolutionarily Conserved Herpesviral Protein Interaction Networks

Herpesviruses constitute a family of large DNA viruses widely spread in vertebrates and causing a variety of different diseases. They possess dsDNA genomes ranging from 120 to 240 kbp encoding between 70 to 170 open reading frames. We previously reported the protein interaction networks of two herpesviruses, varicella-zoster virus (VZV) and Kaposi's sarcoma-associated herpesvirus (KSHV). In this study, we systematically tested three additional herpesvirus species, herpes simplex virus 1 (HSV-1), murine cytomegalovirus and Epstein-Barr virus, for protein interactions in order to be able to perform a comparative analysis of all three herpesvirus subfamilies. We identified 735 interactions by genome-wide yeast-two-hybrid screens (Y2H), and, together with the interactomes of VZV and KSHV, included a total of 1,007 intraviral protein interactions in the analysis. Whereas a large number of interactions have not been reported previously, we were able to identify a core set of highly conserved protein interactions, like the interaction between HSV-1 UL33 with the nuclear egress proteins UL31/UL34. Interactions were conserved between orthologous proteins despite generally low sequence similarity, suggesting that function may be more conserved than sequence. By combining interactomes of different species we were able to systematically address the low coverage of the Y2H system and to extract biologically relevant interactions which were not evident from single species

Genome-Scale Analysis of Translation Elongation with a Ribosome Flow Model

We describe the first large scale analysis of gene translation that is based on a model that takes into account the physical and dynamical nature of this process. The Ribosomal Flow Model (RFM) predicts fundamental features of the translation process, including translation rates, protein abundance levels, ribosomal densities and the relation between all these variables, better than alternative (‘non-physical’) approaches. In addition, we show that the RFM can be used for accurate inference of various other quantities including genes' initiation rates and translation costs. These quantities could not be inferred by previous predictors. We find that increasing the number of available ribosomes (or equivalently the initiation rate) increases the genomic translation rate and the mean ribosome density only up to a certain point, beyond which both saturate. Strikingly, assuming that the translation system is tuned to work at the pre-saturation point maximizes the predictive power of the model with respect to experimental data. This result suggests that in all organisms that were analyzed (from bacteria to Human), the global initiation rate is optimized to attain the pre-saturation point. The fact that similar results were not observed for heterologous genes indicates that this feature is under selection. Remarkably, the gap between the performance of the RFM and alternative predictors is strikingly large in the case of heterologous genes, testifying to the model's promising biotechnological value in predicting the abundance of heterologous proteins before expressing them in the desired host

Emerging Roles of PAR-1 and PAFR in Melanoma Metastasis

Melanoma growth, angiogenesis and metastatic progression are strongly promoted by the inflammatory tumor microenvironment due to high levels of cytokine and chemokine secretion by the recruited inflammatory and stromal cells. In addition, platelets and molecular components of procoagulant pathways have been recently emerging as critical players of tumor growth and metastasis. In particular, thrombin, through the activity of its receptor protease-activated receptor-1 (PAR-1), regulates tumor cell adhesion to platelets and endothelial cells, stimulates tumor angiogenesis, and promotes tumor growth and metastasis. Notably, in many tumor types including melanoma, PAR-1 expression directly correlates with their metastatic phenotype and is directly responsible for the expression of interleukin-8, matrix metalloproteinase-2 (MMP-2), vascular endothelial growth factor, platelet-derived growth factor, and integrins. Another proinflammatory receptor–ligand pair, platelet-activating factor (PAF) and its receptor (PAFR), have been shown to act as important modulators of tumor cell adhesion to endothelial cells, angiogenesis, tumor growth and metastasis. PAF is a bioactive lipid produced by a variety of cells from membrane glycerophospholipids in the same reaction that releases arachidonic acid, and can be secreted by platelets, inflammatory cells, keratinocytes and endothelial cells. We have demonstrated that in metastatic melanoma cells, PAF stimulates the phosphorylation of cyclic adenosine monophosphate response element-binding protein (CREB) and activating transcription factor 1 (ATF-1), which results in overexpression of MMP-2 and membrane type 1-MMP (membrane type 1-MMP). Since only metastatic melanoma cells overexpress CREB/ATF-1, we propose that metastatic melanoma cells are better equipped than their non-metastatic counterparts to respond to PAF within the tumor microenvironment. The evidence supporting the hypothesis that the two G-protein coupled receptors, PAR-1 and PAFR, contribute to the acquisition of the metastatic phenotype of melanoma is presented and discussed

An Information Theoretic, Microfluidic-Based Single Cell Analysis Permits Identification of Subpopulations among Putatively Homogeneous Stem Cells

An incomplete understanding of the nature of heterogeneity within stem cell populations remains a major impediment to the development of clinically effective cell-based therapies. Transcriptional events within a single cell are inherently stochastic and can produce tremendous variability, even among genetically identical cells. It remains unclear how mammalian cellular systems overcome this intrinsic noisiness of gene expression to produce consequential variations in function, and what impact this has on the biologic and clinical relevance of highly ‘purified’ cell subgroups. To address these questions, we have developed a novel method combining microfluidic-based single cell analysis and information theory to characterize and predict transcriptional programs across hundreds of individual cells. Using this technique, we demonstrate that multiple subpopulations exist within a well-studied and putatively homogeneous stem cell population, murine long-term hematopoietic stem cells (LT-HSCs). These subgroups are defined by nonrandom patterns that are distinguishable from noise and are consistent with known functional properties of these cells. We anticipate that this analytic framework can also be applied to other cell types to elucidate the relationship between transcriptional and phenotypic variation

Sequestration of Highly Expressed mRNAs in Cytoplasmic Granules, P-Bodies, and Stress Granules Enhances Cell Viability

Transcriptome analyses indicate that a core 10%–15% of the yeast genome is modulated by a variety of different stresses. However, not all the induced genes undergo translation, and null mutants of many induced genes do not show elevated sensitivity to the particular stress. Elucidation of the RNA lifecycle reveals accumulation of non-translating mRNAs in cytoplasmic granules, P-bodies, and stress granules for future regulation. P-bodies contain enzymes for mRNA degradation; under stress conditions mRNAs may be transferred to stress granules for storage and return to translation. Protein degradation by the ubiquitin-proteasome system is elevated by stress; and here we analyzed the steady state levels, decay, and subcellular localization of the mRNA of the gene encoding the F-box protein, UFO1, that is induced by stress. Using the MS2L mRNA reporter system UFO1 mRNA was observed in granules that colocalized with P-bodies and stress granules. These P-bodies stored diverse mRNAs. Granules of two mRNAs transported prior to translation, ASH1-MS2L and OXA1-MS2L, docked with P-bodies. HSP12 mRNA that gave rise to highly elevated protein levels was not observed in granules under these stress conditions. ecd3, pat1 double mutants that are defective in P-body formation were sensitive to mRNAs expressed ectopically from strong promoters. These highly expressed mRNAs showed elevated translation compared with wild-type cells, and the viability of the mutants was strongly reduced. ecd3, pat1 mutants also exhibited increased sensitivity to different stresses. Our interpretation is that sequestration of highly expressed mRNAs in P-bodies is essential for viability. Storage of mRNAs for future regulation may contribute to the discrepancy between the steady state levels of many stress-induced mRNAs and their proteins. Sorting of mRNAs for future translation or decay by individual cells could generate potentially different phenotypes in a genetically identical population and enhance its ability to withstand stress