Pseudomonas aeruginosa PA14 produces R-bodies, extendable protein polymers with roles in host colonization and virulence

R-bodies are long, extendable protein polymers formed in the cytoplasm of some bacteria; they are best known for their role in killing of paramecia by bacterial endosymbionts. Pseudomonas aeruginosa PA14, an opportunistic pathogen of diverse hosts, contains genes (referred to as the reb cluster) with potential to confer production of R-bodies and that have been implicated in virulence. Here, we show that products of the PA14 reb cluster associate with R-bodies and control stochastic expression of R-body structural genes. PA14 expresses reb genes during colonization of plant and nematode hosts, and R-body production is required for full virulence in nematodes. Analyses of nematode ribosome content and immune response indicate that P. aeruginosa R-bodies act via a mechanism involving ribosome cleavage and translational inhibition. Our observations provide insight into the biology of R-body production and its consequences during P. aeruginosa infection.

Noise suppression in stochastic genetic circuits using PID controllers

Inside individual cells, protein population counts are subject to molecular noise due to low copy numbers and the inherent probabilistic nature of biochemical processes. We investigate the effectiveness of proportional, integral and derivative (PID) based feedback controllers to suppress protein count fluctuations originating from two noise sources: bursty expression of the protein, and external disturbance in protein synthesis. Designs of biochemical reactions that function as PID controllers are discussed, with particular focus on individual controllers separately, and the corresponding closed-loop system is analyzed for stochastic controller realizations. Our results show that proportional controllers are effective in buffering protein copy number fluctuations from both noise sources, but this noise suppression comes at the cost of reduced static sensitivity of the output to the input signal. In contrast, integral feedback has no effect on the protein noise level from stochastic expression, but significantly minimizes the impact of external disturbances, particularly when the disturbance comes at low frequencies. Counter-intuitively, integral feedback is found to amplify external disturbances at intermediate frequencies. Next, we discuss the design of a coupled feedforward-feedback biochemical circuit that approximately functions as a derivate controller. Analysis using both analytical methods and Monte Carlo simulations reveals that this derivative controller effectively buffers output fluctuations from bursty stochastic expression, while maintaining the static input-output sensitivity of the open-loop system. In summary, this study provides a systematic stochastic analysis of biochemical controllers, and paves the way for their synthetic design and implementation to minimize deleterious fluctuations in gene product levels.

How clustered protocadherin binding specificity is tuned for neuronal self/non-self-recognition

The stochastic expression of fewer than 60 clustered protocadherin (cPcdh) isoforms provides diverse identities to individual vertebrate neurons and a molecular basis for self/non-self-discrimination. cPcdhs form chains mediated by alternating cis and trans interactions between apposed membranes, which has been suggested to signal self-recognition. Such a mechanism requires that cPcdh cis dimers form promiscuously to generate diverse recognition units, and that trans interactions have precise specificity so that isoform mismatches terminate chain growth. However, the extent to which cPcdh interactions fulfill these requirements has not been definitively demonstrated. Here we report biophysical experiments showing that cPcdh cis interactions are promiscuous, but with preferences favoring formation of heterologous cis dimers. Trans-homophilic interactions are remarkably precise, with no evidence for heterophilic interactions between different isoforms. A new C-type cPcdh crystal structure and mutagenesis data help to explain these observations. Overall, the interaction characteristics we report for cPcdhs help explain their function in neuronal self/non-self-discrimination

Stochastic Expression of Invasion Genes in Plasmodium falciparum Schizonts

Genetically identical cells are known to exhibit differential phenotypes in the same environmental conditions. These phenotypic variants are linked to transcriptional stochasticity and have been shown to contribute towards adaptive flexibility of a wide range of unicellular organisms. Here, we investigated transcriptional heterogeneity and stochastic gene expression in P. falciparum by performing single cell RNA sequencing on blood stage schizonts. Our data reveals significant transcriptional variations in the schizonts stage with a distinct group of highly variable invasion gene transcripts being identified. Moreover, our data reflected several diversification processes including putative developmental “checkpoint”; transcriptomically distinct parasite sub-populations and transcriptional switches in variable gene families (var, rifin, phist). Most of these features of transcriptional variability were preserved in isogenic parasite cell populations (albeit with a lesser amplitude) suggesting a role of epigenetic factors in cell-to-cell transcriptional variations in human malaria parasites. Lastly, we applied quantitative RT-PCR and RNA-FISH approach and confirmed stochastic expression of merozoites surface proteins encoding msp1, msp3, msp7, rhoptry protein encoding rhopH2 and erythrocyte binding antigen encoding eba181, some of which represent key candidates for invasion blocking vaccines.

Optimum Threshold Minimizes Noise in Timing of Intracellular Events

ABSTRACTHow the noisy expression of regulatory proteins affects timing of intracellular events is an intriguing fundamental problem that influences diverse cellular processes. Here we use the bacteriophage λ to study event timing in individual cells where cell lysis is the result of expression and accumulation of a single protein (holin) in the Escherchia coli cell membrane up to a critical threshold level. Site-directed mutagenesis of the holin gene was used to generate phage variants that vary in their timing of lysis from 30 to 190 min. Observation of the lysis times of single cells reveals an intriguing finding – the noise in lysis timing first decreases with increasing lysis time to reach a minimum, and then sharply increases at longer longer lysis times. A mathematical model with stochastic expression of holin together with dilution from cell growth was sufficient to explain the non-monotonic noise profile, and identify holin accumulation thresholds that generate precision in lysis timing.

Stochastic Expression of Sae-Dependent Virulence Genes during Staphylococcus aureus Biofilm Development Is Dependent on SaeS

ABSTRACT The intricate process of biofilm formation in the human pathogen Staphylococcus aureus involves distinct stages during which a complex mixture of matrix molecules is produced and modified throughout the developmental cycle. Early in biofilm development, a subpopulation of cells detaches from its substrate in an event termed “exodus” that is mediated by SaePQRS-dependent stochastic expression of a secreted staphylococcal nuclease, which degrades extracellular DNA within the matrix, causing the release of cells and subsequently allowing for the formation of metabolically heterogenous microcolonies. Since the SaePQRS regulatory system is involved in the transcriptional control of multiple S. aureus virulence factors, the expression of several additional virulence genes was examined within a developing biofilm by introducing fluorescent gene reporter plasmids into wild-type S. aureus and isogenic regulatory mutants and growing these strains in a microfluidic system that supplies the bacteria with a constant flow of media while simultaneously imaging developing biofilms in 5-min intervals. This study demonstrated that multiple virulence genes, including nuc, were expressed stochastically within a specialized subpopulation of cells in nascent biofilms. We demonstrated that virulence genes regulated by SaePQRS were stochastically expressed in nearly all strains examined whereas Agr-regulated genes were expressed more homogenously within maturing microcolonies. The commonly used Newman strain contains a variant of SaeS (SaeSP) that confers constitutive kinase activity to the protein and caused this strain to lack the stochastic expression pattern observed in other strain backgrounds. Importantly, repair of the SaeSP allele resulting in reversion to the well-conserved SaeSL allele found in other strains restored stochastic expression in this strain. IMPORTANCE Staphylococcus aureus is an important human pathogen capable of colonizing diverse tissue types and inducing severe disease in both immunocompromised and otherwise healthy individuals. Biofilm infections caused by this bacterial species are of particular concern because of their persistence, even in the face of intensive therapeutic intervention. The results of the current study demonstrate the stochastic nature of Sae-mediated virulence gene expression in S. aureus and indicate that this regulatory system may function as a “bistable switch” in a manner similar to that seen with regulators controlling competence gene expression in Bacillus subtilis and persister cell formation in Escherichia coli. The results of this study provide a new perspective on the complex mechanisms utilized by S. aureus during the establishment of infections.

Synchronization of stochastic expressions drives the clustering of functionally related genes

Functionally related genes tend to be chromosomally clustered in eukaryotic genomes even after the exclusion of tandem duplicates, but the biological significance of this widespread phenomenon is unclear. We propose that stochastic expression fluctuations of neighboring genes resulting from chromatin dynamics are more or less synchronized such that their expression ratio is more stable than that for unlinked genes. Consequently, chromosomal clustering could be advantageous when the expression ratio of the clustered genes needs to stay constant, for example, because of the accumulation of toxic compounds when this ratio is altered. Evidence from manipulative experiments on the yeast GAL cluster, comprising three chromosomally adjacent genes encoding enzymes catalyzing consecutive reactions in galactose catabolism, unequivocally supports this hypothesis and elucidates how disorder in one biological phenomenon—gene expression noise—could prompt the emergence of order in another—genome organization.

Power-Law Scaling in the Internal Variability of Cumulus Cloud Size Distributions due to Subsampling and Spatial Organization

In this study, the spatial structure of cumulus cloud populations is investigated using three-dimensional snapshots from large-domain LES experiments. The aim is to understand and quantify the internal variability in cloud size distributions due to subsampling effects and spatial organization. A set of idealized shallow cumulus cases is selected with varying degrees of spatial organization, including a slowly organizing marine precipitating case and five more quickly organizing diurnal cases over land. A subdomain analysis is applied, yielding cloud number distributions at sample sizes ranging from severely undersampled to nearly complete. A strong power-law scaling is found in the relation between cloud number variability and subdomain size, reflecting an inverse linear relation. Scaling subdomain size by cloud size yields a data collapse across time points and cases, highlighting the role played by cloud spacing in controlling the stochastic variability. Spatial organization acts on top of this baseline model by increasing the maximum cloud size and by enhancing the variability in the number of smallest clouds. This reflects that the smaller clouds start to live on top of larger-scale thermodynamic structures, such as cold pools, which favor or inhibit their formation. Compositing all continental cumulus cases suggests the existence of a prototype diurnal time dependence in the spatial organization. A simple stochastic expression for cloud number variability is proposed that is formulated in terms of two dimensionless groups, which allows objective estimation of the degree of spatial organization in simulated and observed cumulus cloud populations.

Intron-targeted mutagenesis reveals roles for Dscam1 RNA pairing-mediated splicing bias in neuronal wiring

Drosophila melanogaster Down syndrome cell adhesion molecule (Dscam1) can potentially generate 38,016 different isoforms through stochastic, yet highly biased, alternative splicing. Genetic studies demonstrated that stochastic expression of multiple Dscam1 isoforms provides each neuron with a unique identity for self/non-self-discrimination. However, due to technical obstacles, the functional significance of the highly specific bias in isoform expression remains entirely unknown. Here, we provide conclusive evidence that Dscam1 splicing bias is required for precise mushroom body (MB) axonal wiring in flies in a variable exon-specific manner. We showed that targeted deletion of the intronic docking site perturbed base pairing-mediated regulation of inclusion of variable exons. Unexpectedly, we generated mutant flies with normal overall Dscam1 protein levels and an identical number but global changes in exon 4 and exon 9 isoform bias (DscamΔ4D−/− and DscamΔ9D−/−), respectively. DscamΔ9D−/− mutant exhibited remarkable mushroom body defects, which were correlated with the extent of the disrupted isoform bias. By contrast, the DscamΔ4D−/− animals exhibited a much less severe defective phenotype than DscamΔ9D−/− animals, suggestive of a variable domain-specific requirement for isoform bias. Importantly, mosaic analysis revealed that changes in isoform bias caused axonal defects but did not influence the self-avoidance of axonal branches. We concluded that, in contrast to the Dscam1 isoform number that provides the molecular basis for neurite self-avoidance, isoform bias may play a non-repulsive role in mushroom body axonal wiring.

Rare Stochastic Expression of O6-Methylguanine- DNA Methyltransferase (MGMT) in MGMT-Negative Melanoma Cells Determines Immediate Emergence of Drug-Resistant Populations upon Treatment with Temozolomide In Vitro and In Vivo

The chemotherapeutic agent temozolomide (TMZ) kills tumor cells preferentially via alkylation of the O6-position of guanine. However, cells that express the DNA repair enzyme O6-methylguanine-DNA methyltransferase (MGMT), or harbor deficient DNA mismatch repair (MMR) function, are profoundly resistant to this drug. TMZ is in clinical use for melanoma, but objective response rates are low, even when TMZ is combined with O6-benzylguanine (O6BG), a potent MGMT inhibitor. We used in vitro and in vivo models of melanoma to characterize the early events leading to cellular TMZ resistance. Melanoma cell lines were exposed to a single treatment with TMZ, at physiologically relevant concentrations, in the absence or presence of O6BG. Surviving clones and mass cultures were analyzed by Western blot, colony formation assays, and DNA methylation studies. Mice with melanoma xenografts received TMZ treatment, and tumor tissue was analyzed by immunohistochemistry. We found that MGMT-negative melanoma cell cultures, before any drug treatment, already harbored a small fraction of MGMT-positive cells, which survived TMZ treatment and promptly became the dominant cell type within the surviving population. The MGMT-negative status in individual cells was not stable, as clonal selection of MGMT-negative cells again resulted in a mixed population harboring MGMT-positive, TMZ-resistant cells. Blocking the survival advantage of MGMT via the addition of O6BG still resulted in surviving clones, although at much lower frequency and independent of MGMT, and the resistance mechanism of these clones was based on a common lack of expression of MSH6, a key MMR enzyme. TMZ treatment of mice implanted with MGMT-negative melanoma cells resulted in effective tumor growth delay, but eventually tumor growth resumed, with tumor tissue having become MGMT positive. Altogether, these data reveal stochastic expression of MGMT as a pre-existing, key determinant of TMZ resistance in melanoma cell lines. Although MGMT activity can effectively be eliminated by pharmacologic intervention with O6BG, additional layers of TMZ resistance, although considerably rarer, are present as well and minimize the cytotoxic impact of TMZ/O6BG combination treatment. Our results provide rational explanations regarding clinical observations, where the TMZ/O6BG regimen has yielded mostly disappointing outcomes in melanoma patients.