Stochastic Turing patterns of trichomes in Arabidopsis leaves.

Biological patterns that emerge during the morphogenesis of multicellular organisms can display high precision at large scales, while at cellular scales, cells exhibit large fluctuations stemming from cell-cell differences in molecular copy numbers also called demographic noise. We study the conflicting interplay between high precision and demographic noise in trichome patterns on the epidermis of wild-type Arabidopsis thaliana leaves, as a two-dimensional model system. We carry out a statistical characterization of these patterns and show that their power spectra display fat tails-a signature compatible with noise-driven stochastic Turing patterns-which are absent in power spectra of patterns driven by deterministic instabilities. We then present a theoretical model that includes demographic noise stemming from birth-death processes of genetic regulators which we study analytically and by stochastic simulations. The model captures the observed experimental features of trichome patterns.

Target search by an imported conjugative DNA element for a unique integration site along a bacterial chromosome during horizontal gene transfer.

Integrative and conjugative elements (ICEs) are mobile genetic elements that can transfer by conjugation to recipient cells. Some ICEs integrate into a unique site in the genome of their hosts. We studied quantitatively the process by which an ICE searches for its unique integration site in the Bacillus subtilis chromosome. We followed the motion of both ICEBs1 and the chromosomal integration site in real time within individual cells. ICEBs1 exhibited a wide spectrum of dynamical behaviors, ranging from rapid sub-diffusive displacements crisscrossing the cell, to kinetically trapped states. The chromosomal integration site moved sub-diffusively and exhibited pronounced dynamical asymmetry between longitudinal and transversal motions, highlighting the role of chromosomal structure and the heterogeneity of the bacterial interior in the search. The successful search for and subsequent recombination into the integration site is a key step in the acquisition of integrating mobile genetic elements. Our findings provide new insights into intracellular transport processes involving large DNA molecules.

Robust stochastic Turing patterns in the development of a one-dimensional cyanobacterial organism.

Under nitrogen deprivation, the one-dimensional cyanobacterial organism Anabaena sp. PCC 7120 develops patterns of single, nitrogen-fixing cells separated by nearly regular intervals of photosynthetic vegetative cells. We study a minimal, stochastic model of developmental patterns in Anabaena that includes a nondiffusing activator, two diffusing inhibitor morphogens, demographic fluctuations in the number of morphogen molecules, and filament growth. By tracking developing filaments, we provide experimental evidence for different spatiotemporal roles of the two inhibitors during pattern maintenance and for small molecular copy numbers, justifying a stochastic approach. In the deterministic limit, the model yields Turing patterns within a region of parameter space that shrinks markedly as the inhibitor diffusivities become equal. Transient, noise-driven, stochastic Turing patterns are produced outside this region, which can then be fixed by downstream genetic commitment pathways, dramatically enhancing the robustness of pattern formation, also in the biologically relevant situation in which the inhibitors' diffusivities may be comparable.

Transcript degradation and noise of small RNA-controlled genes in a switch activated network in Escherichia coli.

Post-transcriptional regulatory processes may change transcript levels and affect cell-to-cell variability or noise. We study small-RNA downregulation to elucidate its effects on noise in the iron homeostasis network of Escherichia coli In this network, the small-RNA RyhB undergoes stoichiometric degradation with the transcripts of target genes in response to iron stress. Using single-molecule fluorescence in situ hybridization, we measured transcript numbers of the RyhB-regulated genes sodB and fumA in individual cells as a function of iron deprivation. We observed a monotonic increase of noise with iron stress but no evidence of theoretically predicted, enhanced stoichiometric fluctuations in transcript numbers, nor of bistable behavior in transcript distributions. Direct detection of RyhB in individual cells shows that its noise is much smaller than that of these two targets, when RyhB production is significant. A generalized two-state model of bursty transcription that neglects RyhB fluctuations describes quantitatively the dependence of noise and transcript distributions on iron deprivation, enabling extraction of in vivo RyhB-mediated transcript degradation rates. The transcripts' threshold-linear behavior indicates that the effective in vivo interaction strength between RyhB and its two target transcripts is comparable. Strikingly, the bacterial cell response exhibits Fur-dependent, switch-like activation instead of a graded response to iron deprivation.

Spatial fluctuations in expression of the heterocyst differentiation regulatory gene hetR in Anabaena filaments.

Under nitrogen deprivation, filaments of the cyanobacterium Anabaena undergo a process of development, resulting in a one-dimensional pattern of nitrogen-fixing heterocysts separated by about ten photosynthetic vegetative cells. Many aspects of gene expression before nitrogen deprivation and during the developmental process remain to be elucidated. Furthermore, the coupling of gene expression fluctuations between cells along a multicellular filament is unknown. We studied the statistics of fluctuations of gene expression of HetR, a transcription factor essential for heterocyst differentiation, both under steady-state growth in nitrogen-rich conditions and at different times following nitrogen deprivation, using a chromosomally-encoded translational hetR-gfp fusion. Statistical analysis of fluorescence at the individual cell level in wild-type and mutant filaments demonstrates that expression fluctuations of hetR in nearby cells are coupled, with a characteristic spatial range of circa two to three cells, setting the scale for cellular interactions along a filament. Correlations between cells predominantly arise from intercellular molecular transfer and less from cell division. Fluctuations after nitrogen step-down can build up on those under nitrogen-replete conditions. We found that under nitrogen-rich conditions, basal, steady-state expression of the HetR inhibitor PatS, cell-cell communication influenced by the septal protein SepJ and positive HetR auto-regulation are essential determinants of fluctuations in hetR expression and its distribution along filaments. A comparison between the expression of hetR-gfp under nitrogen-rich and nitrogen-poor conditions highlights the differences between the two HetR inhibitors PatS and HetN, as well as the differences in specificity between the septal proteins SepJ and FraC/FraD. Activation, inhibition and cell-cell communication lie at the heart of developmental processes. Our results show that proteins involved in these basic ingredients combine together in the presence of inevitable stochasticity in gene expression, to control the coupled fluctuations of gene expression that give rise to a one-dimensional developmental pattern in this organism.

Location of the unique integration site on an Escherichia coli chromosome by bacteriophage lambda DNA in vivo.

The search for specific sequences on long genomes is a key process in many biological contexts. How can specific target sequences be located with high efficiency, within physiologically relevant times? We addressed this question for viral integration, a fundamental mechanism of horizontal gene transfer driving prokaryotic evolution, using the infection of Escherichia coli bacteria with bacteriophage λ and following the establishment of a lysogenic state. Following the targeting process in individual live E. coli cells in real time revealed that λ DNA remains confined near the entry point of a cell following infection. The encounter between the 15-bp-long target sequence on the chromosome and the recombination site on the viral genome is facilitated by the directed motion of bacterial DNA generated during chromosome replication, in conjunction with constrained diffusion of phage DNA. Moving the native bacterial integration site to different locations on the genome and measuring the integration frequency in these strains reveals that the frequencies of the native site and a site symmetric to it relative to the origin are similar, whereas both are significantly higher than when the integration site is moved near the terminus, consistent with the replication-driven mechanism we propose. This novel search mechanism is yet another example of the exquisite coevolution of λ with its host.

Effects of post-transcriptional regulation on phenotypic noise in Escherichia coli.

Cell-to-cell variations in protein abundance, called noise, give rise to phenotypic variability between isogenic cells. Studies of noise have focused on stochasticity introduced at transcription, yet the effects of post-transcriptional regulatory processes on noise remain unknown. We study the effects of RyhB, a small-RNA of Escherichia coli produced on iron stress, on the phenotypic variability of two of its downregulated target proteins, using dual chromosomal fusions to fluorescent reporters and measurements in live individual cells. The total noise of each of the target proteins is remarkably constant over a wide range of RyhB production rates despite cells being in stress. In fact, coordinate downregulation of the two target proteins by RyhB reduces the correlation between their levels. Hence, an increase in phenotypic variability under stress is achieved by decoupling the expression of different target proteins in the same cell, rather than by an increase in the total noise of each. Extrinsic noise provides the dominant contribution to the total protein noise over the total range of RyhB production rates. Stochastic simulations reproduce qualitatively key features of our observations and show that a feed-forward loop formed by transcriptional extrinsic noise, an sRNA and its target genes exhibits strong noise filtration capabilities.

Spatio-temporal coherence of circadian clocks and temporal control of differentiation in Anabaena filaments.

Circadian clock arrays in multicellular filaments of the heterocyst-forming cyanobacterium Anabaena sp. strain PCC 7120 display remarkable spatio-temporal coherence under nitrogen-replete conditions. To shed light on the interplay between circadian clocks and the formation of developmental patterns, we followed the expression of a clock-controlled gene under nitrogen deprivation, at the level of individual cells. Our experiments showed that differentiation into heterocysts took place preferentially within a limited interval of the circadian clock cycle, that gene expression in different vegetative intervals along a developed filament was discoordinated, and that the circadian clock was active in individual heterocysts. Furthermore, Anabaena mutants lacking the kaiABC genes encoding the circadian clock core components produced heterocysts but failed in diazotrophy. Therefore, genes related to some aspect of nitrogen fixation, rather than early or mid-heterocyst differentiation genes, are likely affected by the absence of the clock. A bioinformatics analysis supports the notion that RpaA may play a role as master regulator of clock outputs in Anabaena, the temporal control of differentiation by the circadian clock and the involvement of the clock in proper diazotrophic growth. Together, these results suggest that under nitrogen-deficient conditions, the clock coherent unit in Anabaena is reduced from a full filament under nitrogen-rich conditions to the vegetative cell interval between heterocysts.IMPORTANCECircadian clocks, from unicellular organisms to animals, temporally align biological processes to day and night cycles. We study the dynamics of a circadian clock-controlled gene at the individual cell level in the multicellular filamentous cyanobacterium Anabaena, under nitrogen-stress conditions. Under these conditions, some cells along filaments differentiate to carry out atmospheric nitrogen fixation and lose their capability for oxygenic photosynthesis. We found that clock synchronization is limited to organismic units of contiguous photosynthetic cells, contrary to nitrogen-replete conditions in which clocks are synchronized over a whole filament. We provided evidence that the circadian clock regulates the process of differentiation, allowing it to occur preferentially within a limited time window during the circadian clock period. Lastly, we present evidence that the signal from the core clock to clock-regulated genes is conveyed in Anabaena as in unicellular cyanobacteria.

Stability and gene strand bias of lambda prophages and chromosome organization in Escherichia coli.

Temperate phage-mediated horizontal gene transfer is a potent driver of genetic diversity in the evolution of bacteria. Most lambdoid prophages in Escherichia coli are integrated into the chromosome with the same orientation with respect to the direction of chromosomal replication, and their location on the chromosome is far from homogeneous. To better understand these features, we studied the interplay between lysogenic and lytic states of phage lambda in both native and inverted integration orientations at the wild-type integration site as well as at other sites on the bacterial chromosome. Measurements of free phage released by spontaneous induction showed that the stability of lysogenic states is affected by location and orientation along the chromosome, with stronger effects near the origin of replication. Competition experiments and range expansions between lysogenic strains with opposite orientations and insertion loci indicated that there are no major differences in growth. Moreover, measurements of the level of transcriptional bursts of the cI gene coding for the lambda phage repressor using single-molecule fluorescence in situ hybridization resulted in similar levels of transcription for both orientations and prophage location. We postulate that the preference for a given orientation and location is a result of a balance between the maintenance of lysogeny and the ability to lyse.IMPORTANCEThe integration of genetic material of temperate bacterial viruses (phages) into the chromosomes of bacteria is a potent evolutionary force, allowing bacteria to acquire in one stroke new traits and restructure the information in their chromosomes. Puzzlingly, this genetic material is preferentially integrated in a particular orientation and at non-random sites on the bacterial chromosome. The work described here reveals that the interplay between the maintenance of the stability of the integrated phage, its ability to excise, and its localization along the chromosome plays a key role in setting chromosomal organization in Escherichia coli.

Caught in the act: the lifetime of synaptic intermediates during the search for homology on DNA.

Homologous recombination plays pivotal roles in DNA repair and in the generation of genetic diversity. To locate homologous target sequences at which strand exchange can occur within a timescale that a cell's biology demands, a single-stranded DNA-recombinase complex must search among a large number of sequences on a genome by forming synapses with chromosomal segments of DNA. A key element in the search is the time it takes for the two sequences of DNA to be compared, i.e. the synapse lifetime. Here, we visualize for the first time fluorescently tagged individual synapses formed by RecA, a prokaryotic recombinase, and measure their lifetime as a function of synapse length and differences in sequence between the participating DNAs. Surprisingly, lifetimes can be approximately 10 s long when the DNAs are fully heterologous, and much longer for partial homology, consistently with ensemble FRET measurements. Synapse lifetime increases rapidly as the length of a region of full homology at either the 3'- or 5'-ends of the invading single-stranded DNA increases above 30 bases. A few mismatches can reduce dramatically the lifetime of synapses formed with nearly homologous DNAs. These results suggest the need for facilitated homology search mechanisms to locate homology successfully within the timescales observed in vivo.

Impaired cell-cell communication in the multicellular cyanobacterium Anabaena affects carbon uptake, photosynthesis, and the cell wall.

Cell-cell communication is an essential attribute of multicellular organisms. The effects of perturbed communication were studied in septal protein mutants of the heterocyst-forming filamentous cyanobacterium Anabaena sp. PCC 7120 model organism. Strains bearing sepJ and sepJ/fraC/fraD deletions showed differences in growth, pigment absorption spectra, and spatial patterns of expression of the hetR gene encoding a heterocyst differentiation master regulator. Global changes in gene expression resulting from deletion of those genes were mapped by RNA sequencing analysis of wild-type and mutant strains, both under nitrogen-replete and nitrogen-poor conditions. The effects of sepJ and fraC/fraD deletions were non-additive, and perturbed cell-cell communication led to significant changes in global gene expression. Most significant effects, related to carbon metabolism, included increased expression of genes encoding carbon uptake systems and components of the photosynthetic apparatus, as well as decreased expression of genes encoding cell wall components related to heterocyst differentiation and to polysaccharide export.

Robust, coherent, and synchronized circadian clock-controlled oscillations along Anabaena filaments.

Circadian clocks display remarkable reliability despite significant stochasticity in biomolecular reactions. We study the dynamics of a circadian clock-controlled gene at the individual cell level in Anabaena sp. PCC 7120, a multicellular filamentous cyanobacterium. We found significant synchronization and spatial coherence along filaments, clock coupling due to cell-cell communication, and gating of the cell cycle. Furthermore, we observed low-amplitude circadian oscillatory transcription of kai genes encoding the post-transcriptional core oscillatory circuit and high-amplitude oscillations of rpaA coding for the master regulator transducing the core clock output. Transcriptional oscillations of rpaA suggest an additional level of regulation. A stochastic one-dimensional toy model of coupled clock cores and their phosphorylation states shows that demographic noise can seed stochastic oscillations outside the region where deterministic limit cycles with circadian periods occur. The model reproduces the observed spatio-temporal coherence along filaments and provides a robust description of coupled circadian clocks in a multicellular organism.

Single cell biology-a Keystone Symposia report.

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Noise⁻Seeded Developmental Pattern Formation in Filamentous Cyanobacteria.

Under nitrogen-poor conditions, multicellular cyanobacteria such as Anabaena sp. PCC 7120 undergo a process of differentiation, forming nearly regular, developmental patterns of individual nitrogen-fixing cells, called heterocysts, interspersed between intervals of vegetative cells that carry out photosynthesis. Developmental pattern formation is mediated by morphogen species that can act as activators and inhibitors, some of which can diffuse along filaments. We survey the limitations of the classical, deterministic Turing mechanism that has been often invoked to explain pattern formation in these systems, and then, focusing on a simpler system governed by birth-death processes, we illustrate pedagogically a recently proposed paradigm that provides a much more robust description of pattern formation: stochastic Turing patterns. We emphasize the essential role that cell-to-cell differences in molecular numbers-caused by inevitable fluctuations in gene expression-play, the so called demographic noise, in seeding the formation of stochastic Turing patterns over a much larger region of parameter space, compared to their deterministic counterparts.

Method for Labeling Transcripts in Individual Escherichia coli Cells for Single-molecule Fluorescence In Situ Hybridization Experiments.

A method is described for labeling individual messenger RNA (mRNA) transcripts in fixed bacteria for use in single-molecule fluorescence in situ hybridization (smFISH) experiments in E. coli. smFISH allows the measurement of cell-to-cell variability in mRNA copy number of genes of interest, as well as the subcellular location of the transcripts. The main steps involved are fixation of the bacterial cell culture, permeabilization of cell membranes, and hybridization of the target transcripts with sets of commercially available short fluorescently-labeled oligonucleotide probes. smFISH can allow the imaging of the transcripts of multiple genes in the same cell, with limitations imposed by the spectral overlap between different fluorescent markers. Following completion of the protocol illustrated below, cells can be readily imaged using a microscope coupled with a camera suitable for low-intensity fluorescence. These images, together with cell contours obtained from segmentation of phase contrast frames, or from cell membrane staining, allow the calculation of the mRNA copy number distribution of a sample of cells using open-source or custom-written software. The labeling method described here can also be applied to image transcripts with stochastic optical reconstruction microscopy (STORM).

Phenotypic noise: effects of post-transcriptional regulatory processes affecting mRNA.

The inherently stochastic nature of biomolecular processes is one of the main sources giving rise to cell-to-cell variations in protein and mRNA abundance, termed noise. Noise in isogenic populations can enhance survival under adverse conditions and stress, and has therefore played a fundamental role in evolution. On the other hand, noise may have detrimental effects and therefore cells must also display robustness to fluctuations and possess mechanisms of control in order to function properly. Noise can be introduced at every step in the cascade of intermediate events resulting in the production of functional proteins. While initial studies of noise focused on stochasticity introduced at the transcriptional level, recent years have witnessed a gradual shift of emphasis into the effects that post-transcriptional processes have on phenotypic noise. Here, we survey the insights that have been gained on the effects of processes that modify RNA transcript populations on phenotypic noise, including regulation by noncoding RNAs in prokaryotes and eukaryotes, alternative splicing and transcriptional interference.

Cancer-associated mutations in BRC domains of BRCA2 affect homologous recombination induced by Rad51.

The tumor suppressor BRCA2 protein plays a major role in the regulation of Rad51-catalyzed homologous recombination. BRCA2 interacts with monomeric Rad51 primarily via conserved BRC domains and coordinates the formation of Rad51 filaments at double-stranded DNA (dsDNA) breaks. A number of cancer-associated mutations in BRC4 and BRC2 domains have been reported. To elucidate their effects on homologous recombination, we studied Rad51 filament formation on single-stranded DNA and dsDNA substrates and Rad51-catalyzed strand exchange, in the presence of wild-type and mutated peptides of either BRC4 or BRC2. While the wild-type BRC2 and BRC4 peptides inhibited filament formation and, thus, strand exchange, the mutated forms decreased significantly these inhibitory effects. These results are consistent with a three-dimensional model for the interface between individual BRC repeats and Rad51. We suggest that mutations at sites crucial for the association between Rad51 and BRC domains impair the ability of BRCA2 to recruit Rad51 to dsDNA breaks, hampering recombinational repair.

Cyclin-dependent kinase inhibitors and JNK act as molecular switches, regulating the choice between growth arrest and apoptosis induced by galectin-8.

Galectin-8, a mammalian beta-galactoside binding lectin, functions as an extracellular matrix protein that forms high affinity interactions with integrins. Here we demonstrated that soluble galectin-8 inhibits cell cycle progression and induces growth arrest. These effects cannot be attributed to interference with cell adhesion but can be attributed to a 4-5-fold increase in the cellular content of the cyclin-dependent kinase inhibitor p21, which was already evident following a 4-h incubation of H1299 cells with galectin-8. The increase in p21 levels was preceded by a 3-5-fold increase in JNK and protein kinase B (PKB) activities. Accordingly, SP600125, the inhibitor of JNK, and wortmannin, the inhibitor of phosphatidylinositol 3-kinase, which is the upstream activator of PKB, inhibited the increase in the cellular content of p21. Furthermore, overexpression of a dominant inhibitory form of SEK1, the upstream kinase regulator of JNK, inhibited both JNK activation and p21 accumulation. When p21 expression was inhibited by cycloheximide, galectin-8 directed the cells toward apoptosis, which involves induction of poly(ADP-ribose) polymerase cleavage. Indeed, galectin-8-induced apoptosis was 2-fold higher in HTC (p21-null) cells when compared with parental HTC cells. Because overexpression of galectin-8 attenuates the rate of DNA synthesis, stable colonies that overexpress and secrete galectin-8 can be generated only in cells overexpressing a growth factor receptor, such as the insulin receptor. These results implicate galectin-8 as a modulator of cellular growth through up-regulation of p21. This process involves activation of JNK, which enhances the synthesis of p21, combined with the activation of PKB, which inhibits p21 degradation. These effects of the lectin depended upon protein-sugar interactions and were induced when galectin-8 was present as a soluble ligand or when it was overexpressed in cells.

Immune regulatory activity of CD34+ progenitor cells: evidence for a deletion-based mechanism mediated by TNF-alpha.

Previous studies suggest that cells within the CD34(+) hematopoietic stem cell compartment are endowed with immune regulatory activity. Furthermore, it is possible to expand the human regulatory cells upon short-term culture of purified CD34+ cells with an early-acting cytokine cocktail. We now show that addition of anti-CD28, anti-CD2, interleukin-2 (IL-2), anti-IL-10, or IL-12 to the bulk mixed lymphocyte reaction (MLR) cannot reverse the inhibitory activity of the CD34+ cells, ruling out anergy-based mechanisms or mechanisms involving Th1-Th2 skewing. Furthermore, phenotyping of cells present after addition of CD34+ cells to the bulk MLR ruled out potential induction of plasmacytoid dendritic precursors, known to be endowed with regulatory activity. In contrast, the inhibitory activity of CD34+ cells could be reversed by adding the caspase inhibitor BD-FMK to the bulk MLR, indicating a deletion-based mechanism. The deletion can be inhibited by anti-tumor necrosis factor-alpha (anti-TNF-alpha) and not by anti-transforming growth factor-beta (anti-TGF-beta), suggesting a potential role for TNF-alpha in the regulatory activity of CD34+ cells.