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.

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.

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, 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.

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.

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.

Application of Algebraic Topology to Homologous Recombination of DNA.

Brouwer's fixed point theorem, a fundamental theorem in algebraic topology proved more than a hundred years ago, states that given any continuous map from a closed, simply connected set into itself, there is a point that is mapped unto itself. Here we point out the connection between a one-dimensional application of Brouwer's fixed point theorem and a mechanism proposed to explain how extension of single-stranded DNA substrates by recombinases of the RecA superfamily facilitates significantly the search for homologous sequences on long chromosomes.

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.

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).

Sub-cellular mRNA localization modulates the regulation of gene expression by small RNAs in bacteria.

Small non-coding RNAs can exert significant regulatory activity on gene expression in bacteria. In recent years, substantial progress has been made in understanding bacterial gene expression by sRNAs. However, recent findings that demonstrate that families of mRNAs show non-trivial sub-cellular distributions raise the question of how localization may affect the regulatory activity of sRNAs. Here we address this question within a simple mathematical model. We show that the non-uniform spatial distributions of mRNA can alter the threshold-linear response that characterizes sRNAs that act stoichiometrically, and modulate the hierarchy among targets co-regulated by the same sRNA. We also identify conditions where the sub-cellular organization of cofactors in the sRNA pathway can induce spatial heterogeneity on sRNA targets. Our results suggest that under certain conditions, interpretation and modeling of natural and synthetic gene regulatory circuits need to take into account the spatial organization of the transcripts of participating genes.

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.

The multicellular nature of filamentous heterocyst-forming cyanobacteria.

Cyanobacteria carry out oxygenic photosynthesis, play a key role in the cycling of carbon and nitrogen in the biosphere, and have had a large impact on the evolution of life and the Earth itself. Many cyanobacterial strains exhibit a multicellular lifestyle, growing as filaments that can be hundreds of cells long and endowed with intercellular communication. Furthermore, under depletion of combined nitrogen, filament growth requires the activity of two interdependent cell types: vegetative cells that fix CO2 and heterocysts that fix N2. Intercellular molecular transfer is essential for signaling involved in the regulation of heterocyst differentiation and for reciprocal nutrition of heterocysts and vegetative cells. Here we review various aspects of multicellularity in cyanobacterial filaments and their differentiation, including filament architecture with emphasis on the structures used for intercellular communication; we survey theoretical models that have been put forward to understand heterocyst patterning and discuss the factors that need to be considered for these models to reflect the biological entity; and finally, since cell division in filamentous cyanobacteria has the peculiarity of producing linked instead of independent cells, we review distinct aspects of cell division in these organisms.

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.

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.

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.

Damped oscillations in the adaptive response of the iron homeostasis network of E. coli.

Living organisms often have to adapt to sudden environmental changes and reach homeostasis. To achieve adaptation, cells deploy motifs such as feedback in their genetic networks, endowing the cellular response with desirable properties. We studied the iron homeostasis network of E. coli, which employs feedback loops to regulate iron usage and uptake, while maintaining intracellular iron at non-toxic levels. Using fluorescence reporters for iron-dependent promoters in bulk and microfluidics-based, single-cell experiments, we show that E. coli cells exhibit damped oscillations in gene expression, following sudden reductions in external iron levels. The oscillations, lasting for several generations, are independent of position along the cell cycle. Experiments with mutants in network components demonstrate the involvement of iron uptake in the oscillations. Our findings suggest that the response is driven by intracellular iron oscillations large enough to induce nearly full network activation/deactivation. We propose a mathematical model based on a negative feedback loop closed by rapid iron uptake, and including iron usage and storage, which captures the main features of the observed behaviour. Taken together, our results shed light on the control of iron metabolism in bacteria and suggest that the oscillations represent a compromise between the requirements of stability and speed of response.

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.

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.