Dynamic organisation of prolactin gene expression in living pituitary tissue

Dynamic Ins2 Gene Activity Defines β-Cell Maturity States

Transcriptional and functional cellular specialization has been described for insulin-secreting β-cells of the endocrine pancreas. However, it is not clear whether β-cell heterogeneity is stable or reflects dynamic cellular states. We investigated the temporal kinetics of endogenous insulin gene activity using live cell imaging, with complementary experiments using FACS and single-cell RNA sequencing, in β-cells from Ins2GFP knockin mice. In vivo staining and FACS analysis of islets from Ins2GFP mice confirmed that at a given moment, ∼25% of β-cells exhibited significantly higher activity at the evolutionarily conserved insulin gene, Ins2. Live cell imaging over days captured Ins2 gene activity dynamics in single β-cells. Autocorrelation analysis revealed a subset of oscillating cells, with mean oscillation periods of 17 h. Increased glucose concentrations stimulated more cells to oscillate and resulted in higher average Ins2 gene activity per cell. Single-cell RNA sequencing showed that Ins2(GFP)HIGH β-cells were enriched for markers of β-cell maturity. Ins2(GFP)HIGH β-cells were also significantly less viable at all glucose concentrations and in the context of endoplasmic reticulum stress. Collectively, our results demonstrate that the heterogeneity of insulin production, observed in mouse and human β-cells, can be accounted for by dynamic states of insulin gene activity.

Cell Networks in Endocrine/Neuroendocrine Gland Function

Reproduction, growth, stress, and metabolism are determined by endocrine/neuroendocrine systems that regulate circulating hormone concentrations. All these systems generate rhythms and changes in hormone pulsatility observed in a variety of pathophysiological states. Thus, the output of endocrine/neuroendocrine systems must be regulated within a narrow window of effective hormone concentrations but must also maintain a capacity for plasticity to respond to changing physiological demands. Remarkably most endocrinologists still have a "textbook" view of endocrine gland organization which has emanated from 20^th century histological studies on thin 2D tissue sections. However, 21^st -century technological advances, including in-depth 3D imaging of specific cell types have vastly changed our knowledge. We now know that various levels of multicellular organization can be found across different glands, that organizational motifs can vary between species and can be modified to enhance or decrease hormonal release. This article focuses on how the organization of cells regulates hormone output using three endocrine/neuroendocrine glands that present different levels of organization and complexity: the adrenal medulla, with a single neuroendocrine cell type; the anterior pituitary, with multiple intermingled cell types; and the pancreas with multiple intermingled cell types organized into distinct functional units. We give an overview of recent methodologies that allow the study of the different components within endocrine systems, particularly their temporal and spatial relationships. We believe the emerging findings about network organization, and its impact on hormone secretion, are crucial to understanding how homeostatic regulation of endocrine axes is carried out within endocrine organs themselves. © 2022 American Physiological Society. Compr Physiol 12:3371-3415, 2022.

Transcription Factor Pit-1 Affects Transcriptional Timing in the Dual-Promoter Human Prolactin Gene

Gene transcription occurs in short bursts interspersed with silent periods, and these kinetics can be altered by promoter structure. The effect of alternate promoter architecture on transcription bursting is not known. We studied the human prolactin (hPRL) gene that contains 2 promoters, a pituitary-specific promoter that requires the transcription factor Pit-1 and displays dramatic transcriptional bursting activity and an alternate upstream promoter that is active in nonpituitary tissues. We studied large hPRL genomic fragments with luciferase reporters, and used bacterial artificial chromosome recombineering to manipulate critical promoter regions. Stochastic switch mathematical modelling of single-cell time-lapse luminescence image data revealed that the Pit-1-dependent promoter showed longer, higher-amplitude transcriptional bursts. Knockdown studies confirmed that the presence of Pit-1 stabilized and prolonged periods of active transcription. Pit-1 therefore plays an active role in establishing the timing of transcription cycles, in addition to its cell-specific functions.

Calcium dynamics and chromatin remodelling underlie heterogeneity in prolactin transcription

Pituitary cells have been reported to show spontaneous calcium oscillations and dynamic transcription cycles. To study both processes in the same living cell in real time, we used rat pituitary GH3 cells stably expressing human prolactin-luciferase or prolactin-EGFP reporter gene constructs loaded with a fluorescent calcium indicator and measured activity using single-cell time-lapse microscopy. We observed heterogeneity between clonal cells in the calcium activity and prolactin transcription in unstimulated conditions. There was a significant correlation between cells displaying spontaneous calcium spikes and cells showing spontaneous bursts in prolactin expression. Notably, cells showing no basal calcium activity showed low prolactin expression but elicited a significantly greater transcriptional response to BayK8644 compared to cells showing basal calcium activity. This suggested the presence of two subsets of cells within the population at any one time. Fluorescence-activated cell sorting was used to sort cells into two populations based on the expression level of prolactin-EGFP however, the bimodal pattern of expression was restored within 26 h. Chromatin immunoprecipitation showed that these sorted populations were distinct due to the extent of histone acetylation. We suggest that maintenance of a heterogeneous bimodal population is a fundamental characteristic of this cell type and that calcium activation and histone acetylation, at least in part, drive prolactin transcriptional competence.

Patterns of prolactin secretion

Prolactin is pleotropic in nature affecting multiple tissues throughout the body. As a consequence of the broad range of functions, regulation of anterior pituitary prolactin secretion is complex and atypical as compared to other pituitary hormones. Many studies have provided insight into the complex hypothalamic-pituitary networks controlling prolactin secretion patterns in different species using a range of techniques. Here, we review prolactin secretion in both males and females; and consider the different patterns of prolactin secretion across the reproductive cycle in representative female mammals with short versus long luteal phases and in seasonal breeders. Additionally, we highlight changes in the pattern of secretion during pregnancy and lactation, and discuss the wide range of adaptive functions that prolactin may have in these important physiological states.

Discontinuous transcription

Numerous studies based on new single-cell and single-gene techniques show that individual genes can be transcribed in short bursts or pulses accompanied by changes in pulsing frequencies. Since so many examples of such discontinuous or fluctuating transcription have been found from prokaryotes to mammals, it now seems to be a common mode of gene expression. In this review we discuss the occurrence of the transcriptional fluctuations, the techniques used for their detection, their putative causes, kinetic characteristics, and probable physiological significance.

Disentangling juxtacrine from paracrine signalling in dynamic tissue

Prolactin is a major hormone product of the pituitary gland, the central endocrine regulator. Despite its physiological importance, the cell-level mechanisms of prolactin production are not well understood. Having significantly improved the resolution of real-time-single-cell-GFP-imaging, the authors recently revealed that prolactin gene transcription is highly dynamic and stochastic yet shows space-time coordination in an intact tissue slice. However, it still remains an open question as to what kind of cellular communication mediates the observed space-time organization. To determine the type of interaction between cells we developed a statistical model. The degree of similarity between two expression time series was studied in terms of two distance measures, Euclidean and geodesic, the latter being a network-theoretic distance defined to be the minimal number of edges between nodes, and this was used to discriminate between juxtacrine from paracrine signalling. The analysis presented here suggests that juxtacrine signalling dominates. To further determine whether the coupling is coordinating transcription or post-transcriptional activities we used stochastic switch modelling to infer the transcriptional profiles of cells and estimated their similarity measures to deduce that their spatial cellular coordination involves coupling of transcription via juxtacrine signalling. We developed a computational model that involves an inter-cell juxtacrine coupling, yielding simulation results that show space-time coordination in the transcription level that is in agreement with the above analysis. The developed model is expected to serve as the prototype for the further study of tissue-level organised gene expression for epigenetically regulated genes, such as prolactin.

New Understanding of β-Cell Heterogeneity and In Situ Islet Function

Insulin-secreting β-cells are heterogeneous in their regulation of hormone release. While long known, recent technological advances and new markers have allowed the identification of novel subpopulations, improving our understanding of the molecular basis for heterogeneity. This includes specific subpopulations with distinct functional characteristics, developmental programs, abilities to proliferate in response to metabolic or developmental cues, and resistance to immune-mediated damage. Importantly, these subpopulations change in disease or aging, including in human disease. Although discovering new β-cell subpopulations has substantially advanced our understanding of islet biology, a point of caution is that these characteristics have often necessarily been identified in single β-cells dissociated from the islet. β-Cells in the islet show extensive communication with each other via gap junctions and with other cell types via diffusible chemical messengers. As such, how these different subpopulations contribute to in situ islet function, including during plasticity, is not well understood. We will discuss recent findings revealing functional β-cell subpopulations in the intact islet, the underlying basis for these identified subpopulations, and how these subpopulations may influence in situ islet function. Furthermore, we will discuss the outlook for emerging technologies to gain further insight into the role of subpopulations in in situ islet function.

Asymmetry between Activation and Deactivation during a Transcriptional Pulse

Transcription in eukaryotic cells occurs in gene-specific bursts or pulses of activity. Recent studies identified a spectrum of transcriptionally active “on-states,” interspersed with periods of inactivity, but these “off-states” and the process of transcriptional deactivation are poorly understood. To examine what occurs during deactivation, we investigate the dynamics of switching between variable rates. We measured live single-cell expression of luciferase reporters from human growth hormone or human prolactin promoters in a pituitary cell line. Subsequently, we applied a statistical variable-rate model of transcription, validated by single-molecule FISH, to estimate switching between transcriptional rates. Under the assumption that transcription can switch to any rate at any time, we found that transcriptional activation occurs predominantly as a single switch, whereas deactivation occurs with graded, stepwise decreases in transcription rate. Experimentally altering cAMP signalling with forskolin or chromatin remodelling with histone deacetylase inhibitor modifies the duration of defined transcriptional states. Our findings reveal transcriptional activation and deactivation as mechanistically independent, asymmetrical processes.

•

Gene transcription switches between variable rates

•

Single-cell microscopy and mathematical modeling quantifies switch dynamics

•

We observe an asymmetry in the activation/deactivation of transcriptional bursts

Fluctuations of pol I and fibrillarin contents of the nucleoli

Nucleoli are formed on the basis of ribosomal DNA (rDNA) clusters called Nucleolus Organizer Regions (NORs). Each NOR contains multiple genes coding for RNAs of the ribosomal particles. The prominent components of the nucleolar ultrastructure, fibrillar centers (FC) and dense fibrillar components (DFC), together compose FC/DFC units. These units are centers of rDNA transcription by RNA polymerase I (pol I), as well as the early processing events, in which an essential role belongs to fibrillarin. Each FC/DFC unit probably corresponds to a single transcriptionally active gene. In this work, we transfected human-derived cells with GFP-RPA43 (subunit of pol I) and RFP-fibrillarin. Following changes of the fluorescent signals in individual FC/DFC units, we found two kinds of kinetics: 1) the rapid fluctuations with periods of 2-3 min, when the pol I and fibrillarin signals oscillated in anti-phase manner, and the intensities of pol I in the neighboring FC/DFC units did not correlate. 2) fluctuations with periods of 10 to 60 min, in which pol I and fibrillarin signals measured in the same unit did not correlate, but pol I signals in the units belonging to different nucleoli were synchronized. Our data indicate that a complex pulsing activity of transcription as well as early processing is common for ribosomal genes.

An updated view of hypothalamic-vascular-pituitary unit function and plasticity

The discoveries of novel functional adaptations of the hypothalamus and anterior pituitary gland for physiological regulation have transformed our understanding of their interaction. The activity of a small proportion of hypothalamic neurons can control complex hormonal signalling, which is disconnected from a simple stimulus and the subsequent hormone secretion relationship and is dependent on physiological status. The interrelationship of the terminals of hypothalamic neurons and pituitary cells with the vasculature has an important role in determining the pattern of neurohormone exposure. Cells in the pituitary gland form networks with distinct organizational motifs that are related to the duration and pattern of output, and modifications of these networks occur in different physiological states, can persist after cessation of demand and result in enhanced function. Consequently, the hypothalamus and pituitary can no longer be considered as having a simple stratified relationship: with the vasculature they form a tripartite system, which must function in concert for appropriate hypothalamic regulation of physiological processes, such as reproduction. An improved understanding of the mechanisms underlying these regulatory features has implications for current and future therapies that correct defects in hypothalamic-pituitary axes. In addition, recapitulating proper network organization will be an important challenge for regenerative stem cell treatment.

Cell signaling as a cognitive process

Cellular identity as defined through morphology and function emerges from intracellular signaling networks that communicate between cells. Based on recursive interactions within and among these intracellular networks, dynamical solutions in terms of biochemical behavior are generated that can differ from those in isolated cells. In this way, cellular heterogeneity in tissues can be established, implying that cell identity is not intrinsically predetermined by the genetic code but is rather dynamically maintained in a cognitive manner. We address how to experimentally measure the flow of information in intracellular biochemical networks and demonstrate that even simple causality motifs can give rise to rich, context-dependent dynamic behavior. The concept how intercellular communication can result in novel dynamical solutions is applied to provide a contextual perspective on cell differentiation and tumorigenesis.

A continuum model of transcriptional bursting

Transcription occurs in stochastic bursts. Early models based upon RNA hybridisation studies suggest bursting dynamics arise from alternating inactive and permissive states. Here we investigate bursting mechanism in live cells by quantitative imaging of actin gene transcription, combined with molecular genetics, stochastic simulation and probabilistic modelling. In contrast to early models, our data indicate a continuum of transcriptional states, with a slowly fluctuating initiation rate converting the gene between different levels of activity, interspersed with extended periods of inactivity. We place an upper limit of 40 s on the lifetime of fluctuations in elongation rate, with initiation rate variations persisting an order of magnitude longer. TATA mutations reduce the accessibility of high activity states, leaving the lifetime of on- and off-states unchanged. A continuum or spectrum of gene states potentially enables a wide dynamic range for cell responses to stimuli.

DOI:http://dx.doi.org/10.7554/eLife.13051.001

Understanding how gene activity is regulated relies on accurate measurements of the output of genes. Proteins are generated from genes via a multi-step process. In the first step, called transcription, the DNA of a gene is copied by complex cell machinery to create molecules of mRNA. Subsequently, these mRNA molecules are ‘translated’ into proteins.

Previous studies have assayed gene transcription by measuring mRNA production in millions of cells at the same time. The resulting measurements give the impression that transcription occurs as a continuous, smooth process. However, when individual gene transcription is measured in single cells, mRNA production between cells is unexpectedly variable. This challenged the view that transcription is a continuous process.

One idea that explains this variability – the "two-state" or "bursting" model – proposes that genes switch between "on" and "off" states with a certain probability. Thus, at any one time, a gene will be off in many cells and on in others. However, the methods used in these experiments measure mRNA in dead cells, and so the dynamic switching of genes between on and off states was presumed, but not accurately measured.

Corrigan et al. have now imaged the transcription of a single gene – a gene for a protein called actin – in living cells of an amoeba called Dictyostelium. Genetic techniques and computational modeling were then used to explore what affects the variability in this gene’s activity. These approaches revealed that transcription occurs across a spectrum of activity, rather than in rigid on or off states. The transcription process itself may also contribute to where a gene’s activity sits on this spectrum. Furthermore, Corrigan et al. found that a specific DNA sequence found at the start of the actin gene, that is also found in many genes in complex life-forms, is required for the gene to reach the highest levels of activity on the spectrum.

This spectrum of activity states could allow cells to finely tune their responses to the signals they receive. A future challenge will be to assess how the activity of other genes compare to the actin gene and to discover what underlies the variation in the timing of transcription’s different stages.

DOI:http://dx.doi.org/10.7554/eLife.13051.002

Spatially coordinated dynamic gene transcription in living pituitary tissue

Transcription at individual genes in single cells is often pulsatile and stochastic. A key question emerges regarding how this behaviour contributes to tissue phenotype, but it has been a challenge to quantitatively analyse this in living cells over time, as opposed to studying snap-shots of gene expression state. We have used imaging of reporter gene expression to track transcription in living pituitary tissue. We integrated live-cell imaging data with statistical modelling for quantitative real-time estimation of the timing of switching between transcriptional states across a whole tissue. Multiple levels of transcription rate were identified, indicating that gene expression is not a simple binary 'on-off' process. Immature tissue displayed shorter durations of high-expressing states than the adult. In adult pituitary tissue, direct cell contacts involving gap junctions allowed local spatial coordination of prolactin gene expression. Our findings identify how heterogeneous transcriptional dynamics of single cells may contribute to overall tissue behaviour.

Role of Estrogen Response Element in the Human Prolactin Gene: Transcriptional Response and Timing

The use of bacterial artificial chromosome (BAC) reporter constructs in molecular physiology enables the inclusion of large sections of flanking DNA, likely to contain regulatory elements and enhancers regions that contribute to the transcriptional output of a gene. Using BAC recombineering, we have manipulated a 160-kb human prolactin luciferase (hPRL-Luc) BAC construct and mutated the previously defined proximal estrogen response element (ERE) located -1189 bp relative to the transcription start site, to assess its involvement in the estrogen responsiveness of the entire hPRL locus. We found that GH3 cell lines stably expressing Luc under control of the ERE-mutated hPRL promoter (ERE-Mut) displayed a dramatically reduced transcriptional response to 17β-estradiol (E2) treatment compared with cells expressing Luc from the wild-type (WT) ERE hPRL-Luc promoter (ERE-WT). The -1189 ERE controls not only the response to E2 treatment but also the acute transcriptional response to TNFα, which was abolished in ERE-Mut cells. ERE-WT cells displayed a biphasic transcriptional response after TNFα treatment, the acute phase of which was blocked after treatment with the estrogen receptor antagonist 4-hydroxy-tamoxifen. Unexpectedly, we show the oscillatory characteristics of hPRL promoter activity in individual living cells were unaffected by disruption of this crucial response element, real-time bioluminescence imaging showed that transcription cycles were maintained, with similar cycle lengths, in ERE-WT and ERE-Mut cells. These data suggest the -1189 ERE is the dominant response element involved in the hPRL transcriptional response to both E2 and TNFα and, crucially, that cycles of hPRL promoter activity are independent of estrogen receptor binding.

A stochastic transcriptional switch model for single cell imaging data

Gene expression is made up of inherently stochastic processes within single cells and can be modeled through stochastic reaction networks (SRNs). In particular, SRNs capture the features of intrinsic variability arising from intracellular biochemical processes. We extend current models for gene expression to allow the transcriptional process within an SRN to follow a random step or switch function which may be estimated using reversible jump Markov chain Monte Carlo (MCMC). This stochastic switch model provides a generic framework to capture many different dynamic features observed in single cell gene expression. Inference for such SRNs is challenging due to the intractability of the transition densities. We derive a model-specific birth–death approximation and study its use for inference in comparison with the linear noise approximation where both approximations are considered within the unifying framework of state-space models. The methodology is applied to synthetic as well as experimental single cell imaging data measuring expression of the human prolactin gene in pituitary cells.

Anterior pituitary cell networks

Both endocrine and non-endocrine cells of the pituitary gland are organized into structural and functional networks which are formed during embryonic development but which may be modified throughout life. Structural mapping of the various endocrine cell types has highlighted the existence of distinct network motifs and relationships with the vasculature which may relate to temporal differences in their output. Functional characterization of the network activity of growth hormone and prolactin cells has revealed a role for cell organization in gene regulation, the plasticity of pituitary hormone output and remarkably the ability to memorize altered demand. As such, the description of these endocrine cell networks alters the concept of the pituitary from a gland which simply responds to external regulation to that of an oscillator which may memorize information and constantly adapt its coordinated networks' responses to the flow of hypothalamic inputs.

Interrogation of inhibitor of nuclear factor κB α/nuclear factor κB (IκBα/NF-κB) negative feedback loop dynamics: from single cells to live animals in vivo

Full understanding of the biological significance of negative feedback processes requires interrogation at multiple scales as follows: in single cells, cell populations, and live animals in vivo. The transcriptionally coupled IκBα/NF-κB negative feedback loop, a pivotal regulatory node of innate immunity and inflammation, represents a model system for multiscalar reporters. Using a κB(5)→IκBα-FLuc bioluminescent reporter, we rigorously evaluated the dynamics of ΙκBα degradation and subsequent NF-κB transcriptional activity in response to diverse modes of TNFα stimulation. Modulating TNFα concentration or pulse duration yielded complex, reproducible, and differential ΙκBα dynamics in both cell populations and live single cells. Tremendous heterogeneity in the transcriptional amplitudes of individual responding cells was observed, which was greater than the heterogeneity in the transcriptional kinetics of responsive cells. Furthermore, administration of various TNFα doses in vivo generated ΙκBα dynamic profiles in the liver resembling those observed in single cells and populations of cells stimulated with TNFα pulses. This suggested that dose modulation of circulating TNFα was perceived by hepatocytes in vivo as pulses of increasing duration. Thus, a robust bioluminescent reporter strategy enabled rigorous quantitation of NF-κB/ΙκBα dynamics in both live single cells and cell populations and furthermore, revealed reproducible behaviors that informed interpretation of in vivo studies.

The prolactin gene: a paradigm of tissue-specific gene regulation with complex temporal transcription dynamics

Transcription of numerous mammalian genes is highly pulsatile, with bursts of expression occurring with variable duration and frequency. The presence of this stochastic or 'noisy' expression pattern has been relatively unexplored in tissue systems. The prolactin gene provides a model of tissue-specific gene regulation resulting in pulsatile transcription dynamics in both cell lines and endocrine tissues. In most cell culture models, prolactin transcription appears to be highly variable between cells, with differences in transcription pulse duration and frequency. This apparently stochastic transcription is constrained by a transcriptional refractory period, which may be related to cycles of chromatin remodelling. We propose that prolactin transcription dynamics result from the summation of oscillatory cellular inputs and by regulation through chromatin remodelling cycles. Observations of transcription dynamics in cells within pituitary tissue show reduced transcriptional heterogeneity and can be grouped into a small number of distinct patterns. Thus, it appears that the tissue environment is able to reduce transcriptional noise to enable coordinated tissue responses to environmental change. We review the current knowledge on the complex tissue-specific regulation of the prolactin gene in pituitary and extra-pituitary sites, highlighting differences between humans and rodent experimental animal models. Within this context, we describe the transcription dynamics of prolactin gene expression and how this may relate to specific processes occurring within the cell.

A tridimensional view of pituitary development and function

Recent advances in tridimensional (3D) tissue imaging have considerably enriched our view of the pituitary gland and its development. Whereas traditional histology of the pituitary anterior lobe portrayed this tissue as a patchwork of cells, 3D imaging revealed that cells of each lineage form extensive and structured homotypic networks. In the adult gland these networks contribute to the robustness and coordination of the cell response to secretagogs. In addition, the network organization adapts to changes in endocrine environment, as revealed by the sexually dimorphic growth hormone (GH) cell network. Further work is required to establish better the molecular basis for homotypic and heterotypic interactions in the pituitary as well as the implications of these interactions for pituitary function and dysfunction in humans.

Micro-endoscope for in vivo widefield high spatial resolution fluorescent imaging

In this paper we report the design, testing and use of a scannerless probe specifically for minimally invasive imaging of deep tissue in vivo with an epi-fluorescence modality. The probe images a 500 μm diameter field of view through a 710 μm outer diameter probe with a maximum tissue penetration depth of 15 mm specifically configured for eGFP imaging. Example results are given from imaging the pituitary gland of rats and zebrafish hearts with lateral resolution of 2.5 μm.

Peritonitis activates transcription of the human prolactin locus in myeloid cells in a humanized transgenic rat model

Prolactin (PRL) is mainly expressed in the pituitary in rodents, whereas in humans, expression is observed in many extrapituitary sites, including lymphocytes. Due to the lack of adequate experimental models, the function of locally produced PRL in the immune system is largely unknown. Using transgenic rats that express luciferase under the control of extensive human PRL regulatory regions, we characterized immune cell responses to thioglycollate (TG)-induced peritonitis. Resident populations of myeloid cells in the peritoneal cavity of untreated rats expressed barely detectable levels of luciferase. In contrast, during TG-induced peritonitis, cell-specific expression in both neutrophils and monocytes/macrophages in peritoneal exudates increased dramatically. Elevated luciferase expression was also detectable in peripheral blood and bone marrow CD11b(+) cells. Ex vivo stimulation of primary myeloid cells showed activation of the human extrapituitary promoter by TNF-α, lipopolysaccharide, or TG. These findings were confirmed in human peripheral blood monocytes, showing that the transgenic rat provided a faithful model for the human gene. Thus, the resolution of an inflammatory response is associated with dramatic activation of the PRL gene promoter in the myeloid lineage.

Live imaging of nascent RNA dynamics reveals distinct types of transcriptional pulse regulation

Transcription of genes can be discontinuous, occurring in pulses or bursts. It is not clear how properties of transcriptional pulses vary between different genes. We compared the pulsing of five housekeeping and five developmentally induced genes by direct imaging of single gene transcriptional events in individual living Dictyostelium cells. Each gene displayed its own transcriptional signature, differing in probability of firing and pulse duration, frequency, and intensity. In contrast to the prevailing view from both prokaryotes and eukaryotes that transcription displays binary behavior, strongly expressed housekeeping genes altered the magnitude of their transcriptional pulses during development. These nonbinary "tunable" responses may be better suited than stochastic switch behavior for housekeeping functions. Analysis of RNA synthesis kinetics using fluorescence recovery after photobleaching implied modulation of housekeeping-gene pulse strength occurs at the level of transcription initiation rather than elongation. In addition, disparities between single cell and population measures of transcript production suggested differences in RNA stability between gene classes. Analysis of stability using RNAseq revealed no major global differences in stability between developmental and housekeeping transcripts, although strongly induced RNAs showed unusually rapid decay, indicating tight regulation of expression.

Quantitative measurement of single cell dynamics

Over the past 20 years luminescent and fluorescent imaging assays have been developed to report on the dynamics of transcription and protein translocation in single cells. The combination of these measurements with mathematical analysis is having an increasingly significant impact on cell biology. There is an urgent need to translate these assays to the study of cells and tissues in vivo, which requires new tools and technologies. Emergence of these new tools and techniques will further the understanding of the role of signalling and transcriptional dynamics in the generation of cellular heterogeneity and the control of cell fate.

Imaging of endocrine gene expression in a humanized transgenic rat

Reporter gene imaging has revealed cyclical patterns of gene expression in living cells. Transgenic animal studies show that these patterns are modified by tissue architecture.

Tracing cells for tracking cell lineage and clonal behavior

Reconstructing the lineage of cells is central to understanding development and is now also an important issue in stem cell research. Technological advances in genetically engineered permanent cell labeling, together with a multiplicity of fluorescent markers and sophisticated imaging, open new possibilities for prospective and retrospective clonal analysis.

Pulsatile patterns of pituitary hormone gene expression change during development

Important questions in biology have emerged recently concerning the timing of transcription in living cells. Studies on clonal cell lines have shown that transcription is often pulsatile and stochastic, with implications for cellular differentiation. Currently, information regarding transcriptional activity at cellular resolution within a physiological context remains limited. To investigate single-cell transcriptional activity in real-time in living tissue we used bioluminescence imaging of pituitary tissue from transgenic rats in which luciferase gene expression is driven by a pituitary hormone gene promoter. We studied fetal and neonatal pituitary tissue to assess whether dynamic patterns of transcription change during tissue development. We show that gene expression in single cells is highly pulsatile at the time endocrine cells first appear but becomes stabilised as the tissue develops in early neonatal life. This stabilised transcription pattern might depend upon tissue architecture or paracrine signalling, as isolated cells, generated from enzymatic dispersion of the tissue, display pulsatile luminescence. Nascent cells in embryonic tissue also showed coordinated transcription activity over short distances further indicating that cellular context is important for transcription activity. Overall, our data show that cells alter their patterns of gene expression according to their context and developmental stage, with important implications for cellular differentiation.

Molecular morphology of pituitary cells, from conventional immunohistochemistry to fluorescein imaging

In situ hybridization (ISH) at the electron microscopic (EM) level is essential for elucidating the intracellular distribution and role of mRNA in protein synthesis. EM-ISH is considered to be an important tool for clarifying the intracellular localization of mRNA and the exact site of pituitary hormone synthesis on the rough endoplasmic reticulum. A combined ISH and immunohistochemistry (IHC) under EM (EM-ISH&IHC) approach has sufficient ultrastructural resolution, and provides two-dimensional images of the subcellular localization of pituitary hormone and its mRNA in a pituitary cell. The advantages of semiconductor nanocrystals (quantum dots, Qdots) and confocal laser scanning microscopy (CLSM) enable us to obtain three-dimensional images of the subcellular localization of pituitary hormone and its mRNA. Both EM-ISH&IHC and ISH & IHC using Qdots and CLSM are useful for understanding the relationships between protein and mRNA simultaneously in two or three dimensions. CLSM observation of rab3B and SNARE proteins such as SNAP-25 and syntaxin has revealed that both rab3B and SNARE system proteins play important roles and work together as the exocytotic machinery in anterior pituitary cells. Another important issue is the intracellular transport and secretion of pituitary hormone. We have developed an experimental pituitary cell line, GH3 cell, which has growth hormone (GH) linked to enhanced yellow fluorescein protein (EYFP). This stable GH3 cell secretes GH linked to EYFP upon stimulation by Ca²⁺ influx or Ca²⁺ release from storage. This GH3 cell line is useful for the real-time visualization of the intracellular transport and secretion of GH. These three methods from conventional immunohistochemistry and fluorescein imaging allow us to consecutively visualize the process of transcription, translation, transport and secretion of anterior pituitary hormone.

Dynamic analysis of stochastic transcription cycles

In individual mammalian cells the expression of some genes such as prolactin is highly variable over time and has been suggested to occur in stochastic pulses. To investigate the origins of this behavior and to understand its functional relevance, we quantitatively analyzed this variability using new mathematical tools that allowed us to reconstruct dynamic transcription rates of different reporter genes controlled by identical promoters in the same living cell. Quantitative microscopic analysis of two reporter genes, firefly luciferase and destabilized EGFP, was used to analyze the dynamics of prolactin promoter-directed gene expression in living individual clonal and primary pituitary cells over periods of up to 25 h. We quantified the time-dependence and cyclicity of the transcription pulses and estimated the length and variation of active and inactive transcription phases. We showed an average cycle period of approximately 11 h and demonstrated that while the measured time distribution of active phases agreed with commonly accepted models of transcription, the inactive phases were differently distributed and showed strong memory, with a refractory period of transcriptional inactivation close to 3 h. Cycles in transcription occurred at two distinct prolactin-promoter controlled reporter genes in the same individual clonal or primary cells. However, the timing of the cycles was independent and out-of-phase. For the first time, we have analyzed transcription dynamics from two equivalent loci in real-time in single cells. In unstimulated conditions, cells showed independent transcription dynamics at each locus. A key result from these analyses was the evidence for a minimum refractory period in the inactive-phase of transcription. The response to acute signals and the result of manipulation of histone acetylation was consistent with the hypothesis that this refractory period corresponded to a phase of chromatin remodeling which significantly increased the cyclicity. Stochastically timed bursts of transcription in an apparently random subset of cells in a tissue may thus produce an overall coordinated but heterogeneous phenotype capable of acute responses to stimuli.

The living microarray: a high-throughput platform for measuring transcription dynamics in single cells

Background

Current methods of measuring transcription in high-throughput have led to significant improvements in our knowledge of transcriptional regulation and Systems Biology. However, endpoint measurements obtained from methods that pool populations of cells are not amenable to studying time-dependent processes that show cell heterogeneity.

Results

Here we describe a high-throughput platform for measuring transcriptional changes in real time in single mammalian cells. By using reverse transfection microarrays we are able to transfect fluorescent reporter plasmids into 600 independent clusters of cells plated on a single microscope slide and image these clusters every 20 minutes. We use a fast-maturing, destabilized and nuclear-localized reporter that is suitable for automated segmentation to accurately measure promoter activity in single cells. We tested this platform with synthetic drug-inducible promoters that showed robust induction over 24 hours. Automated segmentation and tracking of over 11 million cell images during this period revealed that cells display substantial heterogeneity in their responses to the applied treatment, including a large proportion of transfected cells that do not respond at all.

Conclusions

The results from our single-cell analysis suggest that methods that measure average cellular responses, such as DNA microarrays, RT-PCR and chromatin immunoprecipitation, characterize a response skewed by a subset of cells in the population. Our method is scalable and readily adaptable to studying complex systems, including cell proliferation, differentiation and apoptosis.

Oscillatory control of signalling molecules

The emergence of biological function from the dynamic control of cellular signalling molecules is a fundamental process in biology. Key questions include: How do cells decipher noisy environmental cues, encode these signals to control fate decisions and propagate information through tissues? Recent advances in systems biology, and molecular and cellular biology, exemplified by analyses of signalling via the transcription factor Nuclear Factor kappaB (NF-κB), reveal a critical role of oscillatory control in the regulation of these biological functions. The emerging view is that the oscillatory dynamics of signalling molecules and the epigenetically regulated specificity for target genes contribute to robust regulation of biological function at different levels of cellular organisation through frequency-dependent information encoding.

Noise-driven heterogeneity in the rate of genetic-variant generation as a basis for evolvability

Molecular biologists have long searched for molecular mechanisms responsible for tuning the rate of genetic-variant generation (RGVG) in fluctuating environments. In spite of several bacterial examples, no regulated variation in the RGVG has been identified in eukaryotic systems. Based notably on the example of industrial and pathogenic yeasts, this article proposes a nonregulated molecular evolutionary mechanism for the appearance of the transient increase of the RGVG in eukaryotic cell populations facing challenging environments. The stochastic nature of gene expression allows a model in which the RGVG in the population can be rapidly tuned as a result of a simple Darwinian process acting on noise-driven heterogeneity in the RGVG from cell to cell. The high flexibility conferred through this model could resolve paradoxical situations, especially concerning the mutator phenotype in cancer cells.

Measurement of single-cell dynamics

Populations of cells are almost always heterogeneous in function and fate. To understand the plasticity of cells, it is vital to measure quantitatively and dynamically the molecular processes that underlie cell-fate decisions in single cells. Early events in cell signalling often occur within seconds of the stimulus, whereas intracellular signalling processes and transcriptional changes can take minutes or hours. By contrast, cell-fate decisions, such as whether a cell divides, differentiates or dies, can take many hours or days. Multiparameter experimental and computational methods that integrate quantitative measurement and mathematical simulation of these noisy and complex processes are required to understand the highly dynamic mechanisms that control cell plasticity and fate.