Dynamic patterns of growth hormone gene transcription in individual living pituitary cells

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.

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

Single-Cell Isolation and Gene Analysis: Pitfalls and Possibilities

During the last two decades single-cell analysis (SCA) has revealed extensive phenotypic differences within homogenous cell populations. These phenotypic differences are reflected in the stochastic nature of gene regulation, which is often masked by qualitatively and quantitatively averaging in whole tissue analyses. The ability to isolate transcripts and investigate how genes are regulated at the single cell level requires highly sensitive and refined methods. This paper reviews different strategies currently used for SCA, including harvesting, reverse transcription, and amplification of the RNA, followed by methods for transcript quantification. The review provides the historical background to SCA, discusses limitations, and current and future possibilities in this exciting field of research.

From analog to digital models of gene regulation

Recently, major progress has been made to develop computational models to predict and explain the mechanisms and behaviors of gene regulation. Here, we review progress on how these mechanisms and behaviors have been interpreted with analog models, where cell properties continuously modulate transcription, and digital models, where gene modulation involves discrete activation and inactivation events. We introduce recent experimental approaches, which measure these gene regulatory behaviors at single-cell and single-molecule resolution, and we discuss the integration of these approaches with computational models to reveal biophysical insight. By analyzing simple toy models in the context of existing experimental capabilities, we discuss the interplay between different experiments and different models to measure and interpret gene regulatory behaviors. Finally, we review recent successes in the development of predictive computational models for the control of gene regulation behaviors.

The role of master regulators in the metabolic/transcriptional coupling in breast carcinomas

Metabolic transformations have been reported as involved in neoplasms survival. This suggests a role of metabolic pathways as potential cancer pharmacological targets. Modulating tumor's energy production pathways may become a substantial research area for cancer treatment. The significant role of metabolic deregulation as inducing transcriptional instabilities and consequently whole-system failure, is thus of foremost importance. By using a data integration approach that combines experimental evidence for high-throughput genome wide gene expression, a non-equilibrium thermodynamics analysis, nonlinear correlation networks as well as database mining, we were able to outline the role that transcription factors MEF2C and MNDA may have as main master regulators in primary breast cancer phenomenology, as well as the possible interrelationship between malignancy and metabolic dysfunction. The present findings are supported by the analysis of 1191 whole genome gene expression experiments, as well as probabilistic inference of gene regulatory networks, and non-equilibrium thermodynamics of such data. Other evidence sources include pathway enrichment and gene set enrichment analyses, as well as motif comparison with a comprehensive gene regulatory network (of homologue genes) in Arabidopsis thaliana. Our key finding is that the non-equilibrium free energies provide a realistic description of transcription factor activation that when supplemented with gene regulatory networks made us able to find deregulated pathways. These analyses also suggest a novel potential role of transcription factor energetics at the onset of primary tumor development. Results are important in the molecular systems biology of cancer field, since deregulation and coupling mechanisms between metabolic activity and transcriptional regulation can be better understood by taking into account the way that master regulators respond to physicochemical constraints imposed by different phenotypic conditions.

Non-equilibrium hyperbolic transport in transcriptional regulation

In this work we studied memory and irreversible transport phenomena in a non-equilibrium thermodynamical model for genomic transcriptional regulation. Transcriptional regulation possess an extremely complex phenomenology, and it is, of course, of foremost importance in organismal cell development and in the pathogenesis of complex diseases. A better understanding of the way in which these processes occur is mandatory to optimize the construction of gene regulatory networks, but also to connect these networks with multi-scale phenomena (e.g. metabolism, signalling pathways, etc.) under an integrative Systems Biology-like vision. In this paper we analyzed three simple mechanisms of genetic stimulation: an instant pulse, a periodic biochemical signal and a saturation process with sigmoidal kinetics and from these we derived the system's thermodynamical response, in the form of, for example, anomalous transcriptional bursts.

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.

Quantitative cell-based reporter gene assays using droplet-based microfluidics

We used a droplet-based microfluidic system to perform a quantitative cell-based reporter gene assay for a nuclear receptor ligand. Single Bombyx mori cells are compartmentalized in nanoliter droplets which function as microreactors with a >1000-fold smaller volume than a microtiter-plate well, together with eight or ten discrete concentrations of 20-hydroxyecdysone, generated by on-chip dilution over 3 decades and encoded by a fluorescent label. The simultaneous measurement of the expression of green fluorescent protein by the reporter gene and of the fluorescent label allows construction of the dose-response profile of the hormone at the single-cell level. Screening approximately 7500 cells per concentration provides statistically relevant data that allow precise measurement of the EC(50) (70 nM +/- 12%, alpha = 0.05), in agreement with standard methods as well as with literature data.

Dynamic organisation of prolactin gene expression in living pituitary tissue

Gene expression in living cells is highly dynamic, but temporal patterns of gene expression in intact tissues are largely unknown. The mammalian pituitary gland comprises several intermingled cell types, organised as interdigitated networks that interact functionally to generate co-ordinated hormone secretion. Live-cell imaging was used to quantify patterns of reporter gene expression in dispersed lactotrophic cells or intact pituitary tissue from bacterial artificial chromosome (BAC) transgenic rats in which a large prolactin genomic fragment directed expression of luciferase or destabilised enhanced green fluorescent protein (d2EGFP). Prolactin promoter activity in transgenic pituitaries varied with time across different regions of the gland. Although amplitude of transcriptional responses differed, all regions of the gland displayed similar overall patterns of reporter gene expression over a 50-hour period, implying overall co-ordination of cellular behaviour. By contrast, enzymatically dispersed pituitary cell cultures showed unsynchronised fluctuations of promoter activity amongst different cells, suggesting that transcriptional patterns were constrained by tissue architecture. Short-term, high resolution, single cell analyses in prolactin-d2EGFP transgenic pituitary slice preparations showed varying transcriptional patterns with little correlation between adjacent cells. Together, these data suggest that pituitary tissue comprises a series of cell ensembles, which individually display a variety of patterns of short-term stochastic behaviour, but together yield long-range and long-term coordinated behaviour.

Circadian timekeeping and multiple timescale neuroendocrine rhythms

The hormones secreted by the anterior pituitary gland regulate major functions such as reproduction, as well as body growth and metabolism. Their efficiency of action highly depends on their temporal profile of release in the blood stream. This review summarises the recent evidence suggesting that the circadian clock genes that pace our daily rhythms may also contribute to the regulation of pituitary pulsatility, even in the non 24-h range. This inter-relation between molecular circadian oscillators and endocrine rhythmicities is discussed in light of the longstanding literature that has considered the involvement of the central circadian pacemaker located within the suprachiasmatic nuclei. Other arguments that suggest a role for circadian clock genes outside the suprachiasmatic nuclei are also presented, with a special emphasis on endocrine pituitary cells and hypothalamic neuroendocrine neurones that directly pace pituitary secretion rates.

Developmental programming: impact of prenatal testosterone excess on pre- and postnatal gonadotropin regulation in sheep

The goal of this study was to explore mechanisms that mediate hypersecretion of LH and progressive loss of cyclicity in female sheep exposed during fetal life to excess testosterone. Our working hypothesis was that prenatal testosterone excess, by its androgenic action, amplifies GnRH-induced LH (but not FSH) secretion and, thus, hypersecretion of LH in adulthood, and that this results from altered developmental gene expression of GnRH and estradiol (E2) receptors, gonadotropin subunits, and paracrine factors that differentially regulate LH and FSH synthesis. We observed that, relative to controls, females exposed during fetal life to excess testosterone, as well as the nor-aromatizable androgen dihydrotestosterone, exhibited enhanced LH but not FSH responses to intermittent delivery of GnRH boluses under conditions in which endogenous LH (GnRH) pulses were suppressed. Luteinizing hormone hypersecretion was more evident in adults than in prepubertal females, and it was associated with development of acyclicity. Measurement of pituitary mRNA concentrations revealed that prenatal testosterone excess induced developmental changes in gene expression of pituitary GnRH and E2 receptors and paracrine modulators of LH and FSH synthesis in a manner consistent with subsequent amplification of LH release. Together, this series of studies suggests that prenatal testosterone excess, by its androgenic action, amplifies GnRH-induced LH response, leading to LH hypersecretion and acyclicity in adulthood, and that this programming involves developmental changes in expression of pituitary genes involved in LH and FSH release.

Stochastic mRNA synthesis in mammalian cells

Individual cells in genetically homogeneous populations have been found to express different numbers of molecules of specific proteins. We investigated the origins of these variations in mammalian cells by counting individual molecules of mRNA produced from a reporter gene that was stably integrated into the cell's genome. We found that there are massive variations in the number of mRNA molecules present in each cell. These variations occur because mRNAs are synthesized in short but intense bursts of transcription beginning when the gene transitions from an inactive to an active state and ending when they transition back to the inactive state. We show that these transitions are intrinsically random and not due to global, extrinsic factors such as the levels of transcriptional activators. Moreover, the gene activation causes burst-like expression of all genes within a wider genomic locus. We further found that bursts are also exhibited in the synthesis of natural genes. The bursts of mRNA expression can be buffered at the protein level by slow protein degradation rates. A stochastic model of gene activation and inactivation was developed to explain the statistical properties of the bursts. The model showed that increasing the level of transcription factors increases the average size of the bursts rather than their frequency. These results demonstrate that gene expression in mammalian cells is subject to large, intrinsically random fluctuations and raise questions about how cells are able to function in the face of such noise.

pTcINDEX: a stable tetracycline-regulated expression vector for Trypanosoma cruzi

BACKGROUND

Trypanosoma cruzi is a protozoan pathogen of major medical importance in Latin America. It is also an early diverging eukaryote that displays many unusual biochemical features. The completion of the T. cruzi genome project has highlighted the need to extend the range of techniques available to study gene function. To this end we report the development of a stable tetracycline-dependent expression vector applicable to this parasite and describe in detail the parameters of the system.

RESULTS

We first produced T. cruzi cell lines that constitutively expressed bacteriophage T7 RNA polymerase and the tetracycline repressor protein from a multicopy episome. An integrative vector with an inducible expression site under the control of a tetracycline-regulatable T7 promoter (pTcINDEX) was targeted to the transcriptionally silent ribosomal RNA spacer region of these parasites and transformants selected using a T7 RNA polymerase-dependent hygromycin resistance gene. To test the system we used two marker proteins, luciferase and red fluorescent protein (RFP), and an endogenous parasite protein (a mitochondrial superoxide dismutase). In each case we found that induction was both time and dose-dependent. Luciferase mRNA could be induced by at least 100-fold, and luciferase activity up to 60-fold, within 24 hours of the addition of tetracycline. When we examined RFP induction by confocal microscopy and fluorescence activated cell sorter, we observed very high levels of expression (>1000-fold increase in fluorescence intensity), although this was not synchronous throughout clonal populations. Induction of superoxide dismutase resulted in an 18-fold increase in cellular activity. The observation that a tagged version of the enzyme was correctly targeted to the mitochondrion demonstrates that our expression system may also provide a high-throughput strategy for subcellular localisation.

CONCLUSION

Our results show that pTcINDEX represents a valuable addition to the genetic tools available for T. cruzi. The vector system is sufficiently flexible that it should have widespread uses including inducible expression of tagged proteins, generation of conditional knockout cell lines and the application of dominant-negative approaches.

Mice, melatonin and the circadian system

Melatonin effects are discussed by reviewing results from mice with intact or disrupted melatonin signaling. Melatonin, the neuroendocrine hand of the clock produced in the pineal gland during night, acts upon two receptor subtypes. Melatonin receptors are found in the suprachiasmatic nuclei (SCN), hypophysial pars tuberalis (PT) and adrenal gland. In SCN, melatonin interacts with PACAP, a neuropeptide of the retinohypothalamic tract. Moreover, melatonin acts on the SCN to modulate the activity of the sympathetic nervous system. Melatonin is not required to maintain rhythmic clock gene expression in SCN. By contrast, the rhythmic clock gene expression in PT depends on a melatonin signal interacting with adenosine. Melatonin may also affect clock gene protein levels in the adrenal cortex and influence adrenal functions. In conclusion, melatonin may serve the synchronization of peripheral oscillators by interacting with other neuroactive substances. A stress-reducing potency of melatonin needs to be explored in further studies.

Tumor necrosis factor-alpha activates the human prolactin gene promoter via nuclear factor-kappaB signaling

Pituitary function has been shown to be regulated by an increasing number of intrapituitary factors, including cytokines. Here we show that the important cytokine TNF-alpha activates prolactin gene transcription in pituitary GH3 cells stably expressing luciferase under control of 5 kb of the human prolactin promoter. Similar regulation of the endogenous rat prolactin gene by TNF-alpha in GH3 cells was confirmed using real-time PCR. Luminescence microscopy revealed heterogeneous dynamic response patterns of promoter activity in individual cells. In GH3 cells treated with TNF-alpha, Western blot analysis showed rapid inhibitory protein kappaB (IkappaBalpha) degradation and phosphorylation of p65. Confocal microscopy of cells expressing fluorescence-labeled p65 and IkappaBalpha fusion proteins showed transient cytoplasmic-nuclear translocation and subsequent oscillations in p65 localization and confirmed IkappaBalpha degradation. This was associated with increased nuclear factor kappaB (NF-kappaB)-mediated transcription from an NF-kappaB-responsive luciferase reporter construct. Disruption of NF-kappaB signaling by expression of dominant-negative variants of IkappaB kinases or truncated IkappaBalpha abolished TNF-alpha activation of the prolactin promoter, suggesting that this effect was mediated by NF-kappaB. TNF-alpha signaling was found to interact with other endocrine signals to regulate prolactin gene expression and is likely to be a major paracrine modulator of lactotroph function.

Diurnal variation in CREB phosphorylation and PER1 protein levels in lactotroph cells of melatonin-proficient C3H and melatonin-deficient C57BL mice: similarities and differences

The pineal hormone melatonin plays an important role in the maintenance of rhythmic functions of the hypophyseal pars tuberalis, which controls the lactotroph cells of the pars distalis. To analyze the effects of melatonin deficiency on the activity state of these cells, we have investigated the levels of Ser133-phosphorylated (p)CREB and PER1 protein in immunocytochemically identified lactotroph cells of melatonin-proficient C3H and melatonin-deficient C57BL mice at four different time points of a 12/12 LD cycle. At night, the percentage of lactotroph cells showing a positive nuclear pCREB and PER1 immunoreaction is significantly smaller in C57BL than in C3H mice. In both mouse strains, the percentage of pCREB-immunoreactive cells is minimal in the early morning and gradually increases to reach a maximum in the late night. PER1 levels show a parallel temporal variation in C3H, but in C57BL, they are drastically reduced in the early afternoon. The observation that, during darkness, the percentage of lactotroph cells with nuclear pCREB immunoreaction is significantly higher in C3H than in C57BL mice suggests the existence of a distinct cell population that is under the control of melatonin-dependent intrapituitary signaling. Interestingly, the percentage of pCREB- and PER1-immunoreactive lactotroph cells reaches minimal and maximal values at the same time points. This suggests that the correlation between CREB phosphorylation and PER1 induction differs between these cells and other neuroendocrine centers, e.g., the pineal organ and suprachiasmatic nucleus, displaying a temporal gap between CREB phosphorylation and PER1 induction.

Effects of hormones targeting nuclear receptors on transcriptional regulation of the growth hormone gene in the MtT/S rat somatotrope cell line

We have examined the effects of nuclear receptor hormones such as glucocorticoid, gonadal steroid hormones, thyroid hormone and retinoids on the transcriptional regulation of the 5'-promoter activity of growth hormone (GH) gene using the MtT/S rat pure somatotrope cell line or MtT/SGL, a subclone of MtT/S in which the rat GH gene 5'-promoter (1.7 Kb)-luciferase fusion gene was stably incorporated. RT-PCR analyses revealed that receptors for all the hormones except androgen receptor were expressed in the cell line. Triiodothyronine (T(3), 10 nM) transiently but significantly stimulated the promoter activity of GH gene, whereas retinoic acids (9-cis and all-trans, both 1 microM) showed sustained stimulation. There were no additive effects among the T(3), all-trans, and9-cis retinoic acids. Synthetic glucocorticoid hormone dexamethasone (100 nM) showed an inhibitory effect but, interestingly, significantly enhanced T(3)-stimulated GH promoter activity during long-term incubation. Among the gonadal steroid hormones tested, estradiol and estriol had significant stimulatory effects, and deletion analysis showed that the estrogen effect was maintained with the shortest construct examined (-150 to +6, +1 denotes the transcription start site). These results suggest that thyroid hormone and retinoids stimulate the transcription of GH gene, probably through a common response element, whereas glucocorticoid has both negative and positive effects on GH expression, depending on the combination with other hormones and the time of exposure. Estrogens also have direct stimulatory effects through the proximal promoter region of GH gene.

Heterogeneous regulation of individual lactotroph cells by photoperiod in the Syrian hamster (Mesocricetus auratus)

In many mammals, changes in daylength (photoperiod) regulate multiple aspects of physiology, including the synthesis and secretion of the anterior pituitary hormone prolactin. Here, we tested the hypothesis that individual lactotroph cells exhibit a heterogeneous response to changes in photoperiod, by exploiting a recently developed assay for prolactin gene expression in single pars distalis (PD) cells. Male Syrian hamsters were exposed to either long (LD; 16 h light: 8 h dark) or short (SD; 8 h light: 16 h dark) photoperiods for 12 weeks. Response of the lactotrophic axis to photoperiod was confirmed by the significantly (P<0.01) lower plasma prolactin concentrations in SD than LD hamsters. Analysis of freshly dispersed PD cells by in situ hybridisation demonstrated that photoperiod has no effect (P>0.05) on the proportion of PD cells (approximately 25%) that expressed prolactin mRNA. Heterogeneity of prolactin mRNA expression was observed in both LD and SD. A similar proportion of cells expressed low levels of prolactin mRNA in both photoperiods, suggesting that they may be unresponsive to photoperiod change. In contrast, the remaining PD cells that expressed prolactin mRNA exhibited markedly increased gene expression in LD, consistent with the selective recruitment of a lactotroph subpopulation to a more transcriptionally active state in this photoperiod.