Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings.

Mathematical methods combined with measurements of single-cell dynamics provide a means to reconstruct intracellular processes that are only partly or indirectly accessible experimentally. To obtain reliable reconstructions, the pooling of measurements from several cells of a clonal population is mandatory. However, cell-to-cell variability originating from diverse sources poses computational challenges for such process reconstruction. We introduce a scalable Bayesian inference framework that properly accounts for population heterogeneity. The method allows inference of inaccessible molecular states and kinetic parameters; computation of Bayes factors for model selection; and dissection of intrinsic, extrinsic and technical noise. We show how additional single-cell readouts such as morphological features can be included in the analysis. We use the method to reconstruct the expression dynamics of a gene under an inducible promoter in yeast from time-lapse microscopy data.

Dynamics of multiple lin gene expression in Sphingomonas paucimobilis B90A in response to different hexachlorocyclohexane isomers.

Sphingomonas paucimobilis B90A is able to degrade the alpha-, beta-, gamma-, and delta-isomers of hexachlorocyclohexane (HCH). It contains the genes linA, linB, linC, linD, linE, and linR, which have been implicated in HCH degradation. In this study, dynamic expression of the lin genes was measured in chemostat-grown S. paucimobilis B90A by RNA dot blot hybridization and real-time reverse transcriptase PCR upon exposure to a pulse of different HCH isomers. Irrespective of the addition of HCH, linA, linB, and linC were all expressed constitutively. In contrast, linD and linE were induced with alpha-HCH (2 mg/liter) and gamma-HCH (7 mg/liter). A sharp increase in mRNA levels for linD and linE was observed from 10 to 45 min after the addition of alpha- or gamma-HCH. Induction of linD and linE was not detectable upon the addition of 0.7 mg of gamma-HCH per liter, although the compound was degraded by the cells. The addition of beta-HCH (5 mg/liter) or delta-HCH (20 mg/liter) did not lead to linE and linD induction, despite the fact that 50% of the compounds were degraded. This suggests that degradation of beta- and delta-HCH proceeds by a different pathway than that of alpha- and gamma-HCH.

Expression of SCG10 and stathmin proteins in the rat olfactory system during development and axonal regeneration.

The membrane-associated protein SCG10 is expressed specifically by neuronal cells. Recent experiments have suggested that it promotes neurite outgrowth by increasing microtubule dynamics in growth cones. SCG10 is related to the ubiquitous but neuron-enriched cytosolic protein stathmin. To better understand the role played by SCG10 and stathmin in vivo, we have analyzed the expression and localization of these proteins in both the olfactory epithelium and the olfactory bulb in developing and adult rats, as well as in adult bulbectomized rats. The olfactory epithelium is exceptional in that olfactory receptor neurons constantly regenerate and reinnervate the olfactory bulb throughout animal life-span. SCG10 and stathmin expression in the olfactory receptor neurons was found to be regulated during embryonic and postnatal development and to correlate with neuronal maturation. Whereas SCG10 expression was restricted to immature olfactory receptor neurons (GAP-43-positive, olfactory marker protein-negative), stathmin was also expressed by the basal cells. In the olfactory bulb of postnatal and adult rats, a moderate to strong SCG10 immunoreactivity was present in the olfactory nerve layer, whereas no labeling was detected in the glomerular layer. Olfactory glomeruli also showed no apparent immunoreactivity for several cytoskeletal proteins such as tubulin and microtubule-associated proteins. In unilaterally bulbectomized rats, SCG10 and stathmin were seen to be up-regulated in the regenerating olfactory epithelium at postsurgery stages corresponding to olfactory axon regeneration. Our data strongly suggest that, in vivo, both SCG10 and stathmin may play a role in axonal outgrowth during ontogenesis as well as during axonal regeneration.

Spatial and temporal dynamics of jasmonate synthesis and accumulation in Arabidopsis in response to wounding.

A new metabolite profiling approach combined with an ultrarapid sample preparation procedure was used to study the temporal and spatial dynamics of the wound-induced accumulation of jasmonic acid (JA) and its oxygenated derivatives in Arabidopsis thaliana. In addition to well known jasmonates, including hydroxyjasmonates (HOJAs), jasmonoyl-isoleucine (JA-Ile), and its 12-hydroxy derivative (12-HOJA-Ile), a new wound-induced dicarboxyjasmonate, 12-carboxyjasmonoyl-l-isoleucine (12-HOOCJA-Ile) was discovered. HOJAs and 12-HOOCJA-Ile were enriched in the midveins of wounded leaves, strongly differentiating them from the other jasmonate metabolites studied. The polarity of these oxylipins at physiological pH correlated with their appearance in midveins. When the time points of accumulation of different jasmonates were determined, JA levels were found to increase within 2-5 min of wounding. Remarkably, these changes occurred throughout the plant and were not restricted to wounded leaves. The speed of the stimulus leading to JA accumulation in leaves distal to a wound is at least 3 cm/min. The data give new insights into the spatial and temporal accumulation of jasmonates and have implications in the understanding of long-distance wound signaling in plants.

Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis.

Wounding in multicellular eukaryotes results in marked changes in gene expression that contribute to tissue defense and repair. Using a cDNA microarray technique, we analyzed the timing, dynamics, and regulation of the expression of 150 genes in mechanically wounded leaves of Arabidopsis. Temporal accumulation of a group of transcripts was correlated with the appearance of oxylipin signals of the jasmonate family. Analysis of the coronatine-insensitive coi1-1 Arabidopsis mutant that is also insensitive to jasmonate allowed us to identify a large number of COI1-dependent and COI1-independent wound-inducible genes. Water stress was found to contribute to the regulation of an unexpectedly large fraction of these genes. Comparing the results of mechanical wounding with damage by feeding larvae of the cabbage butterfly (Pieris rapae) resulted in very different transcript profiles. One gene was specifically induced by insect feeding but not by wounding; moreover, there was a relative lack of water stress-induced gene expression during insect feeding. These results help reveal a feeding strategy of P. rapae that may minimize the activation of a subset of water stress-inducible, defense-related genes.

Effect of roscovitine on intracellular calcium dynamics: differential enantioselective responses.

Cyclin-dependent kinases (CDKs) inhibitors have emerged as interesting therapeutic candidates. Of these, (S)-roscovitine has been proposed as potential neuroprotective molecule for stroke while (R)-roscovitine is currently entering phase II clinical trials against cancers and phase I clinical tests against glomerulonephritis. In addition, (R)-roscovitine has been suggested as potential antihypertensive and anti-inflammatory drug. Dysfunction of intracellular calcium balance is a common denominator of these diseases, and the two roscovitine enantiomers (S and R) are known to modulate calcium voltage channel activity differentially. Here, we provide a detailed description of short- and long-term responses of roscovitine on intracellular calcium handling in renal epithelial cells. Short-term exposure to (S)-roscovitine induced a cytosolic calcium peak, which was abolished after stores depletion with cyclopiazonic acid (CPA). Instead, (R)-roscovitine caused a calcium peak followed by a small calcium plateau. Cytosolic calcium response was prevented after stores depletion. Bafilomycin, a selective vacuolar H(+)-ATPase inhibitor, abolished the small calcium plateau. Long-term exposure to (R)-roscovitine significantly reduced the basal calcium level compared to control and (S)-roscovitine treated cells. However, both enantiomers increased calcium accumulation in the endoplasmic reticulum (ER). Consistently, cells treated with (R)-roscovitine showed a significant increase in SERCA activity, whereas (S)-roscovitine incubation resulted in a reduced PMCA expression. We also found a tonic decreased ability to release calcium from the ER, likely via IP3 signaling, under treatment with (S)- or (R)-roscovitine. Together our data revealed that (S)-roscovitine and (R)-roscovitine exert distinct enantiospecific effects on intracellular calcium signaling in renal epithelial cells. This distinct pharmacological profile can be relevant for roscovitine clinical use.

Integrating receptor interactions and dynamics and Interference with endocrine disruptors in the mechanisms of action of PPAR nuclear receptors

Abstract Peroxisome Proliferator-Activated Receptors (PPARs) form a family of three nuclear receptors regulating important cellular and metabolic functions. PPARs control gene expression by directly binding to target promoters as heterodimers with the Retinoid X Receptor (RXR), and their transcriptional activity is enhanced upon activation by natural or pharmacological ligands. The binding of PPAR/RXR heterodimers on target promoters allows the anchoring of a series of coactivators and corepressors involved in promoter remodeling and the recruitment of the transcription machinery. The transcriptional output finally depends on a complex interplay between (i) the respective expression levels of PPARs, RXRs and of other nuclear receptors competing for DNA binding and RXR recruitment, (ii) the availability and the nature of PPAR and RXR ligands, (iii) the expression levels and the nature of the different coactivators and corepressors and (iv) the sequence and the epigenetic status of the promoter. Understanding how all these factors and signals integrate and fine-tune transcription remains a challenge but is necessary to understand the specificity of the physiological functions regulated by PPARs. The work presented herein focuses on the molecular mechanisms of PPAR action and aims at understanding how the interactions and mobility of the receptor modulate transcription in the physiological context of a living cell: Such observations in vivo rely on the use of engineered fluorescent protein chimeras and require the development and the application of complementary imaging techniques such as Fluorescence Recovery After Photobleaching (FRAP), Fluorescence Resonance Energy Transfer (FRET) and Fluorescence Correlation Spectroscopy (FCS). Using such techniques, PPARs are shown to reside solely in the nucleus where they are constitutively associated with RXR but transcriptional activation by ligand binding -does not promote the formation of sub-nuclear structures as observed with other nuclear receptors. In addition, the engagement of unliganded PPARs in large complexes of cofactors in living cells provides a molecular basis for their ligand-independent activity. Ligand binding reduces receptor diffusion by promoting the recruitment of coactivators which further enlarge the size of PPAR complexes to acquire full transcriptional competence. Using these molecular approaches, we deciphered the molecular mechanisms through which phthalates, a class of pollutants from the plastic industry, interfere with PPARγ signaling. Mono-ethyl-hexyl-phthalate (MEHP) binding induces the recruitment of a specific subset of cofactors and translates into the expression of a specific subset of target genes, the transcriptional output being strongly conditioned by the differentiation status of the cell. This selective PPARγ modulation induces limited adipogenic effects in cellular models while exposure to phthalates in animal models leads to protective effects on glucose tolerance and diet-induced obesity. These results demonstrate that phthalates influence lipid and carbohydrate metabolism through complex mechanisms which most likely involve PPARγ but also probably PPARα and PPARß, Altogether, the molecular and physiological demonstration of the interference of pollutants with PPAR action outlines an important role of chemical exposure in metabolic regulations. Résumé Les PPARs (Peroxisome Proliferator-Activated Receptors) forment une famille de récepteurs nucléaires qui régulent des fonctions cellulaires et métaboliques importantes. Les PPARs contrôlent l'expression des gènes en se liant directement à leurs promoteurs sous forme d'hétérodimères avec les récepteurs RXR (Retinoid X Receptor), et leur activité transcriptionnelle est stimulée par la liaison de ligands naturels ou pharmacologiques. L'association des hétérodimères PPAR/RXR avec les promoteurs des gènes cibles permet le recrutement de coactivateurs et de corépresseurs qui vont permettre le remodelage de la chromatine et le recrutement de la machinerie transcriptionnelle. Les actions transcriptionnelles du récepteur dépendent toutefois d'interactions complexes qui sont régulées par (i) le niveau d'expression des PPARs, des RXRs et d'autres récepteurs nucléaires entrant en compétition pour la liaison à l'ADN et l'association avec RXR, (ii) la disponibilité et la nature de ligands de PPAR et de RXR, (iii) les niveaux d'expression et la nature des différents coactivateurs et corépresseurs et (iv) la séquence et le marquage épigénétique des promoteurs. La compréhension des mécanismes qui permettent d'intégrer ces aspects pour assurer une régulation fine de l'activité transcriptionnelle est un défi qu'il est nécessaire de relever pour comprendre la spécificité des fonctions physiologiques régulées par les PPARs. Ce travail concerne l'étude des mécanismes d'action moléculaire des PPARs et vise à mieux comprendre comment les interactions du récepteur avec d'autres protéines ainsi que la mobilité de ce dernier régulent son activité transcriptionnelle dans le contexte physiologique des cellules vivantes. De telles observations reposent sur l'emploi de protéines fusionnées à des protéines fluorescentes ainsi que sur le développement et l'utilisation de techniques d'imagerie complémentaires telles que le FRAP (Fluorescence Recovery After Photobleaching), le FRET (Fluorescence Resonance Energy Transfer) ou la FCS (Fluorescence Corrélation Spectroscopy). En appliquant ces méthodes, nous avons pu montrer que les PPARs résident toujours dans le noyau où ils sont associés de manière constitutive à RXR, mais que l'ajout de ligand n'induit pas la formation de structures sub-nucléaires comme cela a pu être décrit pour d'autres récepteurs nucléaires. De plus, les PPARs sont engagés dans de larges complexes protéiques de cofacteurs en absence de ligand, ce qui procure une explication moléculaire à leur activité ligand-indépendante. La liaison du ligand réduit la vitesse de diffusion du récepteur en induisant le recrutement de coactivateurs qui augmente encore plus la taille des complexes afin d'acquérir un potentiel d'activation maximal. En utilisant ces approches moléculaires, nous avons pu caractériser les mécanismes permettant aux phtalates, une classe de polluants provenant de l'industrie plastique, d'interférer avec PPARγ. La liaison du mono-ethyl-hexyl-phtalate (NERF) à PPARγ induit un recrutement sélectif de cofacteurs, se traduisant par l'induction spécifique d'un sous-ensemble de gènes qui varie en fonction du niveau de différentiation cellulaire. La modulation sélective de PPARγ par le MEHP provoque une adipogenèse modérée dans des modèles cellulaires alors que l'exposition de modèles animaux aux phtalates induit des effets bénéfiques sur la tolérance au glucose et sur le développement de l'obésité. Toutefois, les phtalates ont une action complexe sur le métabolisme glucido-lipidique en faisant intervenir PPARγ mais aussi probablement PPARα et PPARß. Cette démonstration moléculaire et physiologique de l'interférence des polluants avec les récepteurs nucléaires PPAR souligne un rôle important de l'exposition à de tels composés dans les régulations métaboliques.

The evolution of gene expression levels in mammalian organs.

Changes in gene expression are thought to underlie many of the phenotypic differences between species. However, large-scale analyses of gene expression evolution were until recently prevented by technological limitations. Here we report the sequencing of polyadenylated RNA from six organs across ten species that represent all major mammalian lineages (placentals, marsupials and monotremes) and birds (the evolutionary outgroup), with the goal of understanding the dynamics of mammalian transcriptome evolution. We show that the rate of gene expression evolution varies among organs, lineages and chromosomes, owing to differences in selective pressures: transcriptome change was slow in nervous tissues and rapid in testes, slower in rodents than in apes and monotremes, and rapid for the X chromosome right after its formation. Although gene expression evolution in mammals was strongly shaped by purifying selection, we identify numerous potentially selectively driven expression switches, which occurred at different rates across lineages and tissues and which probably contributed to the specific organ biology of various mammals.

SOST expression is restricted to the great arteries during embryonic and neonatal cardiovascular development.

Spatial-temporal regulation of bone morphogenetic protein (BMP) and Wnt activity is essential for normal cardiovascular development, and altered activity of these growth factors causes maldevelopment of the cardiac outflow tract and great arteries. In the present study, we show that SOST, a Dan family member reported to antagonize BMP and Wnt activity, is expressed within the medial vessel wall of the great arteries containing smooth muscle cells. The ascending aorta, aortic arch, brachiocephalic artery, common carotids, and pulmonary trunk were all associated with SOST expressing smooth muscle cells, while the heart itself, including the valves, and more distal arteries, that is, pulmonary arteries, subclavian arteries, and descending aorta, were negative. SOST was expressed from embryonic day 15.5 up to the neonatal period. SOST expression, however, did not correspond with inhibition of Smad-dependent BMP activity or beta-catenin-dependent Wnt activity in the great arteries. Activity of both signaling pathways was already down-regulated before induction of SOST expression.

Evolutionary dynamics of coding and non-coding transcriptomes.

Gene expression changes may underlie much of phenotypic evolution. The development of high-throughput RNA sequencing protocols has opened the door to unprecedented large-scale and cross-species transcriptome comparisons by allowing accurate and sensitive assessments of transcript sequences and expression levels. Here, we review the initial wave of the new generation of comparative transcriptomic studies in mammals and vertebrate outgroup species in the context of earlier work. Together with various large-scale genomic and epigenomic data, these studies have unveiled commonalities and differences in the dynamics of gene expression evolution for various types of coding and non-coding genes across mammalian lineages, organs, developmental stages, chromosomes and sexes. They have also provided intriguing new clues to the regulatory basis and phenotypic implications of evolutionary gene expression changes.

The adaptive dynamics of niche constructing traits in spatially subdivided populations: evolving posthumous extended phenotypes.

Niche construction, by which organisms modify the environment in which they live, has been proposed to affect the evolution of many phenotypic traits. But what about the evolution of a niche constructing trait itself, whose expression changes the pattern of natural selection to which the trait is exposed in subsequent generations? This article provides an inclusive fitness analysis of selection on niche constructing phenotypes, which can affect their environment from local to global scales in arbitrarily spatially subdivided populations. The model shows that phenotypic effects of genes extending far beyond the life span of the actor can be affected by natural selection, provided they modify the fitness of those individuals living in the future that are likely to have inherited the niche construction lineage of the actor. Present benefits of behaviors are thus traded off against future indirect costs. The future costs will generally result from a complicated interplay of phenotypic effects, population demography and environmental dynamics. To illustrate these points, I derive the adaptive dynamics of a trait involved in the consumption of an abiotic resource, where resource abundance in future generations feeds back to the evolutionary dynamics of the trait.

Temporal quantification of MAPK induced expression in single yeast cells.

The quantification of gene expression at the single cell level uncovers novel regulatory mechanisms obscured in measurements performed at the population level. Two methods based on microscopy and flow cytometry are presented to demonstrate how such data can be acquired. The expression of a fluorescent reporter induced upon activation of the high osmolarity glycerol MAPK pathway in yeast is used as an example. The specific advantages of each method are highlighted. Flow cytometry measures a large number of cells (10,000) and provides a direct measure of the dynamics of protein expression independent of the slow maturation kinetics of the fluorescent protein. Imaging of living cells by microscopy is by contrast limited to the measurement of the matured form of the reporter in fewer cells. However, the data sets generated by this technique can be extremely rich thanks to the combinations of multiple reporters and to the spatial and temporal information obtained from individual cells. The combination of these two measurement methods can deliver new insights on the regulation of protein expression by signaling pathways.

Evolutionary dynamics of interlinked public goods traits: an experimental study of siderophore production in Pseudomonas aeruginosa.

Public goods cooperation is common in microbes, and there is much interest in understanding how such traits evolve. Research in recent years has identified several important factors that shape the evolutionary dynamics of such systems, yet few studies have investigated scenarios involving interactions between multiple public goods. Here, we offer general predictions about the evolutionary trajectories of two public goods traits having positive, negative or neutral regulatory influence on one another's expression, and we report on a test of some of our predictions in the context of Pseudomonas aeruginosa's production of two interlinked iron-scavenging siderophores. First, we confirmed that both pyoverdine and pyochelin siderophores do operate as public goods under appropriate environmental conditions. We then tracked their production in lines experimentally evolved under different iron-limitation regimes known to favour different siderophore expression profiles. Under strong iron limitation, where pyoverdine represses pyochelin, we saw a decline in pyoverdine and a concomitant increase in pyochelin - consistent with expansion of pyoverdine-defective cheats derepressed for pyochelin. Under moderate iron limitation, pyochelin declined - again consistent with an expected cheat invasion scenario - but there was no concomitant shift in pyoverdine because cross-suppression between the traits is unidirectional only. Alternating exposure to strong and moderate iron limitation caused qualitatively similar though lesser shifts compared to the constant-environment regimes. Our results confirm that the regulatory interconnections between public goods traits can significantly modulate the course of evolution, yet also suggest how we can start to predict the impacts such complexities will have on phenotypic divergence and community stability.

Blood-Borne Circadian Signal Stimulates Daily Oscillations in Actin Dynamics and SRF Activity.

In peripheral tissues circadian gene expression can be driven either by local oscillators or by cyclic systemic cues controlled by the master clock in the brain's suprachiasmatic nucleus. In the latter case, systemic signals can activate immediate early transcription factors (IETFs) and thereby control rhythmic transcription. In order to identify IETFs induced by diurnal blood-borne signals, we developed an unbiased experimental strategy, dubbed Synthetic TAndem Repeat PROMoter (STAR-PROM) screening. This technique relies on the observation that most transcription factor binding sites exist at a relatively high frequency in random DNA sequences. Using STAR-PROM we identified serum response factor (SRF) as an IETF responding to oscillating signaling proteins present in human and rodent sera. Our data suggest that in mouse liver SRF is regulated via dramatic diurnal changes of actin dynamics, leading to the rhythmic translocation of the SRF coactivator Myocardin-related transcription factor-B (MRTF-B) into the nucleus.

The Munc13 proteins differentially regulate readily releasable pool dynamics and calcium-dependent recovery at a central synapse.

The Munc13 gene family encodes molecules located at the synaptic active zone that regulate the reliability of synapses to encode information over a wide range of frequencies in response to action potentials. In the CNS, proteins of the Munc13 family are critical in regulating neurotransmitter release and synaptic plasticity. Although Munc13-1 is essential for synaptic transmission, it is paradoxical that Munc13-2 and Munc13-3 are functionally dispensable at some synapses, although their loss in other synapses leads to increases in frequency-dependent facilitation. We addressed this issue at the calyx of Held synapse, a giant glutamatergic synapse that we found to express all these Munc13 isoforms. We studied their roles in the regulation of synaptic transmission and their impact on the reliability of information transfer. Through detailed electrophysiological analyses of Munc13-2, Munc13-3, and Munc13-2-3 knock-out and wild-type mice, we report that the combined loss of Munc13-2 and Munc13-3 led to an increase in the rate of calcium-dependent recovery and a change in kinetics of release of the readily releasable pool. Furthermore, viral-mediated overexpression of a dominant-negative form of Munc13-1 at the calyx demonstrated that these effects are Munc13-1 dependent. Quantitative immunohistochemistry using Munc13-fluorescent protein knock-in mice revealed that Munc13-1 is the most highly expressed Munc13 isoform at the calyx and the only one highly colocalized with Bassoon at the active zone. Based on these data, we conclude that Munc13-2 and Munc13-3 isoforms limit the ability of Munc13-1 to regulate calcium-dependent replenishment of readily releasable pool and slow pool to fast pool conversion in central synapses.