Coupling of GnRH concentration and the GnRH receptor-activated gene program

Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases

The pituitary gland directly regulates the reproduction of domestic animals. Research has increasingly focused on the potential regulatory mechanism of non-coding RNA in pituitary development. Little is known about the differential expression pattern of lncRNAs in Hu sheep, a famous sheep breed with high fecundity, and its role in the pituitary gland between the follicular phase and luteal phase. Herein, to identify the transcriptomic differences of the sheep pituitary gland during the estrus cycle, RNA sequencing (RNA-Seq) was performed. The results showed that 3529 lncRNAs and 16,651 mRNAs were identified in the pituitary gland. Among of them, 144 differentially expressed (DE) lncRNA transcripts and 557 DE mRNA transcripts were screened in the follicular and luteal phases. Moreover, GO and KEGG analyses demonstrated that 39 downregulated and 22 upregulated genes interacted with pituitary functions and reproduction. Lastly, the interaction of the candidate lncRNA XR_001039544.4 and its targeted gene LHB were validated in sheep pituitary cells in vitro. LncRNA XR_001039544.4 and LHB showed high expression levels in the luteal phase in Hu sheep. LncRNA XR_001039544.4 is mainly located in the cytoplasm, as determined by FISH analysis, indicating that XR_001039544.4 might act as competing endogenous RNAs for miRNAs to regulate LHB. LncRNA XR_001039544.4 knockdown significantly inhibited LH secretion and cell proliferation. LncRNA XR_001039544.4 may regulate the secretion of LH in the luteal-phase pituitary gland via affecting cell proliferation. Taken together, these findings provided genome-wide lncRNA- and mRNA-expression profiles for the sheep pituitary gland between the follicular and luteal phases, thereby contributing to the elucidation of the molecular mechanisms of pituitary function.

Single nucleus multi-omics regulatory landscape of the murine pituitary

To provide a multi-omics resource and investigate transcriptional regulatory mechanisms, we profile the transcriptome, chromatin accessibility, and methylation status of over 70,000 single nuclei (sn) from adult mouse pituitaries. Paired snRNAseq and snATACseq datasets from individual animals highlight a continuum between developmental epigenetically-encoded cell types and transcriptionally-determined transient cell states. Co-accessibility analysis-based identification of a putative Fshb cis-regulatory domain that overlaps the fertility-linked rs11031006 human polymorphism, followed by experimental validation illustrate the use of this resource for hypothesis generation. We also identify transcriptional and chromatin accessibility programs distinguishing each major cell type. Regulons, which are co-regulated gene sets sharing binding sites for a common transcription factor driver, recapitulate cell type clustering. We identify both cell type-specific and sex-specific regulons that are highly correlated with promoter accessibility, but not with methylation state, supporting the centrality of chromatin accessibility in shaping cell-defining transcriptional programs. The sn multi-omics atlas is accessible at snpituitaryatlas.princeton.edu.

Transcriptomic and Histopathological Effects of Bifenthrin to the Brain of Juvenile Rainbow Trout (Oncorhynchus mykiss)

The increased global use of pyrethroids raises concern for non-target aquatic species. Bifenthrin, among the most predominantly detected pyrethroids in the environment, is frequently measured in water samples above concentrations reported to induce neuroendocrine and neurotoxic effects to several threatened and endangered fish species, such as the Chinook salmon and steelhead trout. To better characterize the neurotoxic effect of bifenthrin to salmonids, rainbow trout were treated with environmentally relevant concentrations of bifenthrin (15 and 30 ng/L) for two weeks and assessed for changes in transcriptomic profiles and histopathological alterations. The top bioinformatic pathways predicted to be impaired in bifenthrin-exposed trout were involved in gonadotropin releasing hormone signaling, the dysregulation of iron homeostasis, reduced extracellular matrix stability and adhesion, and cell death. Subsequent histopathological analysis showed a significant increase in TUNEL positive cells in the cerebellum and optic tectum of bifenthrin-treated trout, relative to controls (p < 0.05). These findings suggest that low, ng/L concentrations of bifenthrin are capable of dysregulating proper neuroendocrine function, impair the structural integrity of the extracellular matrix and cell signaling pathways in the brain, and induce apoptosis in neurons of juvenile salmonids following bifenthrin treatment, which is consistent with metabolomic profiles demonstrating a common target and mechanism.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

Regulatory Architecture of the LβT2 Gonadotrope Cell Underlying the Response to Gonadotropin-Releasing Hormone

The LβT2 mouse pituitary cell line has many characteristics of a mature gonadotrope and is a widely used model system for studying the developmental processes and the response to gonadotropin-releasing hormone (GnRH). The global epigenetic landscape, which contributes to cell-specific gene regulatory mechanisms, and the single-cell transcriptome response variation of LβT2 cells have not been previously investigated. Here, we integrate the transcriptome and genome-wide chromatin accessibility state of LβT2 cells during GnRH stimulation. In addition, we examine cell-to-cell variability in the transcriptional response to GnRH using Gel bead-in-Emulsion Drop-seq technology. Analysis of a bulk RNA-seq data set obtained 45 min after exposure to either GnRH or vehicle identified 112 transcripts that were regulated >4-fold by GnRH (FDR < 0.05). The top regulated transcripts constitute, as determined by Bayesian massive public data integration analysis, a human pituitary-relevant coordinated gene program. Chromatin accessibility [assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq)] data sets generated from GnRH-treated LβT2 cells identified more than 58,000 open chromatin regions, some containing notches consistent with bound transcription factor footprints. The study of the most prominent open regions showed that 75% were in transcriptionally active promoters or introns, supporting their involvement in active transcription. Lhb, Cga, and Egr1 showed significantly open chromatin over their promoters. While Fshb was closed over its promoter, several discrete significantly open regions were found at -40 to -90 kb, which may represent novel upstream enhancers. Chromatin accessibility determined by ATAC-seq was associated with high levels of gene expression determined by RNA-seq. We obtained high-quality single-cell Gel bead-in-Emulsion Drop-seq transcriptome data, with an average of >4,000 expressed genes/cell, from 1,992 vehicle- and 1,889 GnRH-treated cells. While the individual cell expression patterns showed high cell-to-cell variation, representing both biological and measurement variation, the average expression patterns correlated well with bulk RNA-seq data. Computational assignment of each cell to its precise cell cycle phase showed that the response to GnRH was unaffected by cell cycle. To our knowledge, this study represents the first genome-wide epigenetic and single-cell transcriptomic characterization of this important gonadotrope model. The data have been deposited publicly and should provide a resource for hypothesis generation and further study.

Exploring Dynamics and Noise in Gonadotropin-Releasing Hormone (GnRH) Signaling

Gonadotropin-releasing hormone (GnRH) acts via G-protein coupled receptors on pituitary gonadotropes. These are G_q-coupled receptors that mediate acute effects of GnRH on the exocytotic secretion of luteinizing hormone (LH) and follicle-stimulating hormone (FSH), as well as the chronic regulation of their synthesis. FSH and LH control steroidogenesis and gametogenesis in the gonads so GnRH mediates control of reproduction by the central nervous system. GnRH is secreted in short pulses and the effects of GnRH on its target cells are dependent on the dynamics of these pulses. Here we provide a brief overview of the signaling network activated by GnRH with emphasis on the use of high content imaging for their examination. We also describe computational approaches that we have used to simulate GnRH signaling in order to explore dynamics, noise, and information transfer in this system.

Mathematical modeling of gonadotropin-releasing hormone signaling

Gonadotropin-releasing hormone (GnRH) acts via G-protein coupled receptors on pituitary gonadotropes to control reproduction. These are G_q-coupled receptors that mediate acute effects of GnRH on the exocytotic secretion of luteinizing hormone (LH) and follicle-stimulating hormone (FSH), as well as the chronic regulation of their synthesis. GnRH is secreted in short pulses and GnRH effects on its target cells are dependent upon the dynamics of these pulses. Here we overview GnRH receptors and their signaling network, placing emphasis on pulsatile signaling, and how mechanistic mathematical models and an information theoretic approach have helped further this field.

Modeling and high-throughput experimental data uncover the mechanisms underlying Fshb gene sensitivity to gonadotropin-releasing hormone pulse frequency

Neuroendocrine control of reproduction by brain-secreted pulses of gonadotropin-releasing hormone (GnRH) represents a longstanding puzzle about extracellular signal decoding mechanisms. GnRH regulates the pituitary gonadotropin's follicle-stimulating hormone (FSH) and luteinizing hormone (LH), both of which are heterodimers specified by unique β subunits (FSHβ/LHβ). Contrary to Lhb, Fshb gene induction has a preference for low-frequency GnRH pulses. To clarify the underlying regulatory mechanisms, we developed three biologically anchored mathematical models: 1) parallel activation of Fshb inhibitory factors (e.g. inhibin α and VGF nerve growth factor-inducible), 2) activation of a signaling component with a refractory period (e.g. G protein), and 3) inactivation of a factor needed for Fshb induction (e.g. growth differentiation factor 9). Simulations with all three models recapitulated the Fshb expression levels obtained in pituitary gonadotrope cells perifused with varying GnRH pulse frequencies. Notably, simulations altering average concentration, pulse duration, and pulse frequency revealed that the apparent frequency-dependent pattern of Fshb expression in model 1 actually resulted from variations in average GnRH concentration. In contrast, models 2 and 3 showed "true" pulse frequency sensing. To resolve which components of this GnRH signal induce Fshb, we developed a high-throughput parallel experimental system. We analyzed over 4,000 samples in experiments with varying near-physiological GnRH concentrations and pulse patterns. Whereas Egr1 and Fos genes responded only to variations in average GnRH concentration, Fshb levels were sensitive to both average concentration and true pulse frequency. These results provide a foundation for understanding the role of multiple regulatory factors in modulating Fshb gene activity.

GnRH Neuron Activity and Pituitary Response in Estradiol-Induced vs Proestrous Luteinizing Hormone Surges in Female Mice

During the female reproductive cycle, estradiol exerts negative and positive feedback at both the central level to alter gonadotropin-releasing hormone (GnRH) release and at the pituitary to affect response to GnRH. Many studies of the neurobiologic mechanisms underlying estradiol feedback have been done on ovariectomized, estradiol-replaced (OVX+E) mice. In this model, GnRH neuron activity depends on estradiol and time of day, increasing in estradiol-treated mice in the late afternoon, coincident with a daily luteinizing hormone (LH) surge. Amplitude of this surge appears lower than in proestrous mice, perhaps because other ovarian factors are not replaced. We hypothesized GnRH neuron activity is greater during the proestrous-preovulatory surge than the estradiol-induced surge. GnRH neuron activity was monitored by extracellular recordings from fluorescently tagged GnRH neurons in brain slices in the late afternoon from diestrous, proestrous, and OVX+E mice. Mean GnRH neuron firing rate was low on diestrus; firing rate was similarly increased in proestrous and OVX+E mice. Bursts of action potentials have been associated with hormone release in neuroendocrine systems. Examination of the patterning of action potentials revealed a shift toward longer burst duration in proestrous mice, whereas intervals between spikes were shorter in OVX+E mice. LH response to an early afternoon injection of GnRH was greater in proestrous than diestrous or OVX+E mice. These observations suggest the lower LH surge amplitude observed in the OVX+E model is likely not attributable to altered mean GnRH neuron activity, but because of reduced pituitary sensitivity, subtle shifts in action potential pattern, and/or excitation-secretion coupling in GnRH neurons.

Reactive Oxygen Species Link Gonadotropin-Releasing Hormone Receptor Signaling Cascades in the Gonadotrope

Biological rhythms lie at the center of regulatory schemes that control many aspects of living systems. At the cellular level, meaningful responses to external stimuli depend on propagation and quenching of a signal to maintain vigilance for subsequent stimulation or changes that serve to shape and modulate the response. The hypothalamus-pituitary-gonad endocrine axis that controls reproductive development and function relies on control through rhythmic stimulation. Central to this axis is the pulsatile stimulation of the gonadotropes by hypothalamic neurons through episodic release of the neuropeptide gonadotropin-releasing hormone. Alterations in pulsatile stimulation of the gonadotropes result in differential synthesis and secretion of the gonadotropins LH and FSH and changes in the expression of their respective hormone subunit genes. The requirement to amplify signals arising from activation of the gonadotropin-releasing hormone (GnRH) receptor and to rapidly quench the resultant signal to preserve an adaptive response suggests the need for rapid activation and feedback control operating at the level of intracellular signaling. Emerging data suggest that reactive oxygen species (ROS) can fulfill this role in the GnRH receptor signaling through activation of MAP kinase signaling cascades, control of negative feedback, and participation in the secretory process. Results obtained in gonadotrope cell lines or other cell models indicate that ROS can participate in each of these regulatory cascades. We discuss the potential advantage of reactive oxygen signaling for modulating the gonadotrope response to GnRH stimulation and the potential mechanisms for this action. These observations suggest further targets of study for regulation in the gonadotrope.

In vivo oestrogenic modulation of Egr1 and Pitx1 gene expression in female rat pituitary gland

EGR1 and PITX1 are transcription factors required for gonadotroph cell Lhb promoter activation. To determine changes in Egr1 and Pitx1 mRNA levels in central and peripheral pituitary stimulations, an in vivo model based on i.c.v. pulsatile (1 pulse/0.5 h over 2 h) GnRH agonist (1.5 nM buserelin) or antagonist (2 nM antide) microinjections was used. The microinjections were given to ovariectomised and 17β-oestradiol (E2) (3×20 μg), ERA (ESR1) agonist propyl pyrazole triol (PPT) (3×0.5 mg), ERB (ESR2) agonist diarylpropionitrile (DPN) (3×0.5 mg) s.c. pre-treated rats 30 min after last pulse anterior pituitaries were excised. Relative mRNA expression was determined by quantitative RT-PCR (qRT-PCR). Results revealed a gene-specific response for GnRH and/or oestrogenic stimulations in vivo. Buserelin pulses enhanced Egr1 expression by 66% in ovariectomised rats, whereas the oestradiol-supplemented+i.c.v. NaCl-microinjected group showed a 50% increase in Egr1 mRNA expression. The oestrogenic signal was transmitted via ERA (ESR1) and ERB (ESR2) activation as administration of PPT and DPN resulted in 97 and 62%, respectively, elevation in Egr1 mRNA expression. A synergistic action of GnRH agonist and 17β-oestradiol (E2) stimulation of the Egr1 gene transcription in vivo were found. GnRHR activity did not affect Pitx1 mRNA expression; regardless of NaCl, buserelin or antide i.c.v. pulses, s.c. oestrogenic supplementation (with E2, PPT or DPN) consistently decreased (by -46, -48 and -41% respectively) the Pitx1 mRNA in the anterior pituitary gland. Orchestrated Egr1 and Pitx1 activities depending on specific central and peripheral regulatory inputs could be responsible for physiologically variable Lhb gene promoter activation in vivo.

β-catenin regulates GnRH-induced FSHβ gene expression

The regulation of gonadotropin synthesis by GnRH plays an essential role in the neuroendocrine control of reproduction. The known signaling mechanisms involved in gonadotropin synthesis have been expanding. For example, involvement of β-catenin in LHβ induction by GnRH has been discovered. We examined the role of β-catenin in FSHβ gene expression in LβT2 gonadotrope cells. GnRH caused a sustained increase in nuclear β-catenin levels, which was significantly reduced by c-Jun N-terminal kinase (JNK) inhibition. Small interfering RNA-mediated knockdown of β-catenin mRNA demonstrated that induction of FSHβ mRNA by GnRH depended on β-catenin and that regulation of FSHβ by β-catenin occurred independently of the JNK-c-jun pathway. β-Catenin depletion had no impact on FSHβ mRNA stability. In LβT2 cells transfected with FSHβ promoter luciferase fusion constructs, GnRH responsiveness was conferred by the proximal promoter (-944/-1) and was markedly decreased by β-catenin knockdown. However, none of the T-cell factor/lymphoid enhancer factor binding sites in that region were required for promoter activation by GnRH. Chromatin immunoprecipitation further corroborated the absence of direct interaction between β-catenin and the 1.8-kb FSHβ promoter. To elucidate the mechanism for the β-catenin effect, we analyzed approximately 1 billion reads of next-generation RNA sequencing β-catenin knockdown assays and selected the nuclear cofactor breast cancer metastasis-suppressor 1-like (Brms1L) as one candidate for further study. Subsequent experiments confirmed that Brms1L mRNA expression was decreased by β-catenin knockdown as well as by JNK inhibition. Furthermore, knockdown of Brms1L significantly attenuated GnRH-induced FSHβ expression. Thus, our findings indicate that the expression of Brms1L depends on β-catenin activity and contributes to FSHβ induction by GnRH.

Optimized amplification and single-cell analysis identify GnRH-mediated activation of Rap1b in primary rat gonadotropes

Identifying the early gene program induced by GnRH would help understand how GnRH-activated signaling pathways modulate gonadotrope secretory response. We previously analyzed GnRH-induced early genes in LβT2 cells, however these lack GnRH self-potentiation, a physiological attribute of gonadotropes. To minimize cellular heterogeneity, rat primary pituitary cultures were enriched for gonadotropes by 40-60% using a sedimentation gradient. Given the limited number of gonadotropes, RNA was amplified prior to microarray analysis. Thirty-three genes were up-regulated 40 min after GnRH stimulation. Real-time PCR confirmed regulation of several transcripts including fosB, c-fos, egr-2 and rap1b, a small GTPase and member of the Ras family. GnRH stimulated rap1b gene expression in gonadotropes, measured by a sensitive single cell assay. Immunocytochemistry revealed increased Rap1 protein in GnRH-stimulated gonadotropes. These data establish rap1b as a novel gene rapidly induced by GnRH and a candidate to modulate gonadotropin secretion in rat gonadotropes.

Modulation of gonadotropin-releasing hormone-induced extracellular signal-regulated kinase activation by dual-specificity protein phosphatase 1 in LbetaT2 gonadotropes

As the regulator of pituitary reproductive hormone synthesis, the hypothalamic neuropeptide GnRH is the central regulator of reproduction. A hallmark of GnRH action is the differential control of gene expression in pituitary gonadotropes through varied pulsatile stimulation. Among other signaling events, GnRH activation of the ERK family of MAPKs plays a significant role in the transcriptional regulation of the luteinizing hormone β-subunit gene and regulation of cap-dependent translation. We evaluated the ERK response to different GnRH pulse amplitudes in the gonadotrope cell line LβT2. We found that low-amplitude stimulation with GnRH invokes a rapid and transient ERK activation, whereas high-amplitude stimulation invokes a prolonged activation specifically in the cytoplasm fraction of LβT2 cells. Nuclear and cytoplasmic targets of ERK, Ets-like gene 1, and eukaryotic initiation factor 4E, respectively, are similarly activated. Feedback control of ERK activation occurs mainly through the dual-specificity protein phosphatases (DUSPs). DUSP1 is localized to the nucleus in LβT2 cells but DUSP4, another member implicated in GnRH feedback, exists in both the nucleus and cytoplasm. Manipulation of nuclear DUSP activity through overexpression or knockdown of Dusp1 modulates the ERK response to low and high GnRH pulse amplitudes and activation of the Lhb promoter. Dusp1 overexpression abolishes sustained ERK activation and inhibits Lhb promoter activity induced by high amplitude pulses. Conversely, Dusp1 knockdown enhances ERK activation by low-amplitude stimulation and increases stimulation of Lhb promoter activity. We conclude that DUSP1 feedback activity modulates ERK activation and the transcriptional response to GnRH.

Plato's cave algorithm: inferring functional signaling networks from early gene expression shadows

Improving the ability to reverse engineer biochemical networks is a major goal of systems biology. Lesions in signaling networks lead to alterations in gene expression, which in principle should allow network reconstruction. However, the information about the activity levels of signaling proteins conveyed in overall gene expression is limited by the complexity of gene expression dynamics and of regulatory network topology. Two observations provide the basis for overcoming this limitation: a. genes induced without de-novo protein synthesis (early genes) show a linear accumulation of product in the first hour after the change in the cell's state; b. The signaling components in the network largely function in the linear range of their stimulus-response curves. Therefore, unlike most genes or most time points, expression profiles of early genes at an early time point provide direct biochemical assays that represent the activity levels of upstream signaling components. Such expression data provide the basis for an efficient algorithm (Plato's Cave algorithm; PLACA) to reverse engineer functional signaling networks. Unlike conventional reverse engineering algorithms that use steady state values, PLACA uses stimulated early gene expression measurements associated with systematic perturbations of signaling components, without measuring the signaling components themselves. Besides the reverse engineered network, PLACA also identifies the genes detecting the functional interaction, thereby facilitating validation of the predicted functional network. Using simulated datasets, the algorithm is shown to be robust to experimental noise. Using experimental data obtained from gonadotropes, PLACA reverse engineered the interaction network of six perturbed signaling components. The network recapitulated many known interactions and identified novel functional interactions that were validated by further experiment. PLACA uses the results of experiments that are feasible for any signaling network to predict the functional topology of the network and to identify novel relationships.

The hypothalamic GnRH pulse generator: multiple regulatory mechanisms

Pulsatile secretion of gonadotropin-releasing hormone (GnRH) release is an intrinsic property of hypothalamic GnRH neurons. Pulse generation has been attributed to multiple specific mechanisms, including spontaneous electrical activity of GnRH neurons, calcium and cAMP signaling, a GnRH receptor autocrine regulatory component, a GnRH concentration-dependent switch in GnRH receptor (GnRH-R) coupling to specific G proteins, the expression of G protein-coupled receptors (GPCRs) and steroid receptors, and homologous and heterologous interactions between cell membrane receptors expressed in GnRH neurons. The coexistence of multiple regulatory mechanisms for pulsatile GnRH secretion provides a high degree of redundancy in maintaining this crucial component of the mammalian reproductive process. These studies provide insights into the basic cellular and molecular mechanisms involved in GnRH neuronal function.

Microtranscriptome regulation by gonadotropin-releasing hormone

Gonadotropin-releasing hormone (GnRH) regulates biosynthesis in the pituitary gonadotrope via a complex signaling and gene network. Small non-coding microRNAs (miRNA) can play important roles in gene expression. We investigated the microtranscriptome in the mouse L?T2 gonadotrope cell line using microarray, single molecule coincidence detection assays, hairpin real time PCR and LNA (locked nucleic acid) primer-extension PCR. Expression of nearly 200 miRNAs were detected by array and a panel of 101 hairpin real time PCR assays. Within this broad family of expressed miRNAs, GnRH induced upregulation of two miRNA products of the same primary transcript, miR-132 and miR-212, a result confirmed by single molecule, hairpin and LNA assays. Induction peaked 6h after GnRH exposure and showed no significant frequency sensitivity. Bioinformatics analysis was used to predict potential targets of each of these GnRH-regulated miRNAs. These findings suggest the importance of the microtranscriptome in gene control in the gonadotrope and implicate miR-132 and miR-212 in the regulation of GnRH-stimulated biosynthetic response.

GnRH-regulated expression of Jun and JUN target genes in gonadotropes requires a functional interaction between TCF/LEF family members and beta-catenin

GnRH regulates gonadotrope function through a complex transcriptional network that includes three members of the immediate early gene family: Egr1, Jun, and Atf3. These DNA-binding proteins act alone or in pairs to confer hormonal responsiveness to Cga, Lhb, Fshb, and Gnrhr. Herein we suggest that the transcriptional response of Jun requires a functional interaction between the T-cell factor (TCF)/lymphoid enhancer factor (LEF) family of DNA-binding proteins and beta-catenin (officially CTNNB1), a coactivator of TCF/LEF. Supporting data include demonstration that GnRH increases activity of TOPflash, a TCF/LEF-dependent luciferase reporter, in LbetaT2 cells, a gonadotrope-derived cell line. Additional cotransfection experiments indicate that a dominant-negative form of TCF7L2 (TCFDN) that binds DNA, but not beta-catenin, blocks GnRH induction of TOPflash. Overexpression of AXIN, an inhibitor of beta-catenin, also reduces GnRH stimulation of TOPflash. Transduction of LbetaT2 cells with TCFDN adenoviruses diminishes GnRH stimulation of Jun mRNA without altering expression of Egr1 and Atf3, two other immediate early genes that confer GnRH responsiveness. Reduction of beta-catenin in LbetaT2 cells, through stable expression of short hairpin RNA, also selectively compromises GnRH regulation of Jun expression and levels of JUN protein. Finally, overexpression of TCFDN attenuates GnRH regulation of Cga promoter activity, a known downstream target of JUN. Together, these results indicate that GnRH regulation of Jun transcription requires a functional interaction between TCF/LEF and beta-catenin and that alteration of either impacts expression of JUN downstream targets such as Cga.

GNRH induces the unfolded protein response in the LbetaT2 pituitary gonadotrope cell line

The neuropeptide GNRH 1 stimulates the secretion of the reproductive hormone LH in pituitary gonadotropes. Other secretory cell types depend on the unfolded protein response (UPR) pathway to regulate protein synthesis and protect against endoplasmic reticulum (ER) stress in response to differentiation or secretory stimuli. This study investigated the role of the UPR in GNRH action within the LbetaT2 gonadotrope model. Cells were treated with GNRH, and the activation of UPR signaling components and general translational status was examined. The ER-resident stress sensors, Atf6, Eif2ak3, and Ern1, are all present, and GNRH stimulation results in the phosphorylation of eukaryotic translation initiation factor 2A kinase 3 and its downstream effector, eukaryotic translation initiation factor 2A. Additionally, activation of the UPR was confirmed both in LbetaT2 as well as mouse primary pituitary cells through identifying GNRH-induced splicing of Xbp1 mRNA, a transcription factor activated by splicing by the ER stress sensor, ER to nucleus signaling 1. Ribosome profiling revealed that GNRH stimulation caused a transient attenuation in translation, a hallmark of the UPR, remodeling ribosomes from actively translating polysomes to translationally inefficient ribonucleoprotein complexes and monosomes. The transient attenuation of specific mRNAs was also observed. Overall, the results show that GNRH activates components of the UPR pathway, and this pathway may play an important physiological role in adapting the ER of gonadotropes to the burden of their secretory demand.

Welcoming beta-catenin to the gonadotropin-releasing hormone transcriptional network in gonadotropes

GnRH binds its G-coupled protein receptor, GnRHR, on pituitary gonadotropes and stimulates transcription of Cga, Lhb, and Fshb. These three genes encode two heterodimeric glycoprotein hormones, LH and FSH, that act as gonadotropins by regulating gametogenesis and steroidogenesis in both the testes and ovary. GnRH also regulates transcription of Gnrhr. Thus, regulated expression of Cga, Lhb, Fshb, and Gnrhr provides a genomic signature unique to functional gonadotropes. Steadily increasing evidence now indicates that GnRH regulates transcription of its four signature genes indirectly through a hierarchical transcriptional network that includes distinct subclasses of DNA-binding proteins that comprise the immediate early gene (IEG) family. These IEGs, in turn, confer hormonal responsiveness to the four signature genes. Although the IEGs confer responsiveness to GnRH, they cannot act alone. Instead, additional DNA-binding proteins, including the orphan nuclear receptor steroidogenic factor 1, act permissively to allow the four signature genes to respond to GnRH-induced changes in IEG levels. Emerging new findings now indicate that beta-catenin, a transcriptional coactivator and member of the canonical WNT signaling pathway, also plays an essential role in transducing the GnRH signal by interacting with multiple DNA-binding proteins in gonadotropes. Herein we propose that these interactions with beta-catenin define a multicomponent transcriptional network required for regulated expression of the four signature genes of the gonadotrope, Cga, Lhb, Fshb, and Gnrhr.

Estrogen receptor alpha (ERalpha) and insulin-like growth factor I receptor (IGF-IR) cross-talk in the gonadotropic alphaT3-1 cell line

In reproductive tissues such as the breast and the uterus, cell proliferation and differentiation is strongly regulated by complex interactions between estrogen receptor alpha (ERalpha) and growth factor receptors. In the present study, we investigated the potential occurrence of such cross-talk in the murine, gonadotropic alphaT3-1 cell line, which expresses ERalpha and the IGF-I receptor (IGF-IR). Under estrogen-free conditions, basal cell proliferation and ER-mediated gene transcription was strongly inhibited by the selective estrogen receptor modulator (SERM) 4-hydroxy-tamoxifen (4-OH-Tam) and by the pure anti-estrogen ICI 182,780 (ICI). These effects can be reversed by either 17-beta-estradiol (E(2)) or insulin-like growth factor I (IGF-I), both exerting modest mitogenic effects in the alphaT3-1 cell line. Furthermore, IGF-I enhanced both basal and E(2)-induced ER-driven gene transcription. This may be explained, at least in part, by enhanced phosphorylation of ERalpha at serine 118, a prerequisite for the transactivation capacity of the receptor. Finally, the IGF-I-induced response on cell growth and ER-mediated transactivation can be inhibited with either ICI or 4-OH-Tam. In conclusion, our data indicate IGF-IR and ER interactions in the alphaT3-1 cell line, an in vitro model for the pituitary gonadotrophs, hereby suggesting a role of IGF-I in the regulation of gonadotropin synthesis and secretion.

Pulse sensitivity of the luteinizing hormone beta promoter is determined by a negative feedback loop Involving early growth response-1 and Ngfi-A binding protein 1 and 2

The hypothalamic-pituitary-gonadal endocrine axis regulates reproduction through estrous phase-dependent release of the heterodimeric gonadotropic glycoprotein hormones, LH and FSH, from the gonadotropes of the anterior pituitary. Gonadotropin synthesis and release is dependent upon pulsatile stimulation by the hypothalamic neuropeptide GnRH. Alterations in pulse frequency and amplitude alter the relative levels of gonadotropin synthesis and release. The mechanism of interpretation of GnRH pulse frequency and amplitude by gonadotropes is not understood. We have examined gene expression in LbetaT2 gonadotropes under various pulse regimes in a cell perifusion system by microarray and identified 1127 genes activated by tonic or pulsatile GnRH. Distinct patterns of expression are associated with each pulse frequency, but the greatest changes occur at a 60-min or less interpulse interval. The immediate early gene mRNAs encoding early growth response (Egr)1 and Egr2, which activate the gonadotropin LH beta-subunit gene promoter, are stably induced at high pulse frequency. In contrast, mRNAs for the Egr corepressor genes Ngfi-A binding protein Nab1 and Nab2 are stably induced at low pulse frequency. We show that Ngfi-A binding protein members inhibit Egr-mediated frequency-dependent induction of the LH beta-subunit promoter. This pattern of expression suggests a model of pulse frequency detection that acts by suppressing activation by Egr family members at low frequency and allowing activation at sustained high-frequency pulses.

Mixed Analog/Digital Gonadotrope Biosynthetic Response to Gonadotropin-releasing Hormone

Mammalian reproduction requires gonadotropin-releasing hormone (GnRH)-mediated signaling from brain neurons to pituitary gonadotropes. Because the pulses of released GnRH vary greatly in amplitude, we studied the biosynthetic response of the gonadotrope to varying GnRH concentrations, focusing on extracellular-regulated kinase (ERK) phosphorylation and egr1 mRNA and protein production. The overall average level of ERK activation in populations of cells increased non-cooperatively with increasing GnRH and did not show evidence of either ultrasensitivity or bistability. However, automated image analysis of single-cell responses showed that whereas individual gonadotropes exhibited two response states, inactive and active, both the probability of activation and the average response in activated cells increased with increasing GnRH concentration. These data indicate a hybrid single-cell response having both digital (switch-like) and analog (graded) features. Mathematical modeling suggests that the hybrid response can be explained by indirect thresholding of ERK activation resulting from the distributed structure of the GnRH-modulated network. The hybrid response mechanism improves the reliability of noisy reproductive signal transmission from the brain to the pituitary.

Gonadotropin-releasing hormone activation of c-jun, but not early growth response factor-1, stimulates transcription of a luteinizing hormone beta-subunit gene

Transcription of mammalian LH beta-subunit genes (LHbeta) is regulated by GnRH through activation of early growth response factor-1 (Egr-1), which interacts synergistically with steroidogenic factor-1 (Sf-1) and pituitary homeobox-1 (Pitx1) at the promoter; Egr-1 is thought to comprise the major mediator of this effect. However, the proximal promoters of LHbeta genes in lower vertebrates lack an Egr-1 response element yet are responsive to GnRH; we demonstrate here that the promoter of the Chinook salmon LHbeta (csLHbeta) gene is also unresponsive to Egr-1. The homologous LHbeta promoters in other fish contain a conserved estrogen response element-like sequence, which we recently demonstrated is not required for estrogen receptor (ER) alpha association with the csLHbeta gene. Here we show that the estrogen response element-like element is required for the GnRH effect and for a response to c-jun overexpression. Using plasmid immunoprecipitation, we show that after GnRH exposure, c-jun associates with the intact csLHbeta gene promoter through this element. We further show that the effect of c-jun requires its DNA-binding domain and that c-jun interacts with Sf-1 and ERalpha and exerts synergistic effects on promoter activity with Sf-1, ERalpha, and Pitx1. Finally, we demonstrate the role of c-jun in mediating the GnRH effect on this gene through knockdown of c-jun expression or use of a dominant negative. We conclude that c-jun mediation of the GnRH effect on the LHbeta gene may be common in lower vertebrates and may have preceded an evolutionary divergence in the cis-regulatory elements that led to its function being replaced in mammals by Egr-1.

Mechanisms for pulsatile regulation of the gonadotropin subunit genes by GNRH1

The frequency of gonadotropin-releasing hormone (GNRH1, or GnRH) pulses secreted from the hypothalamus determine the ratios of the gonadotropin subunit genes luteinizing hormone beta (Lhb), follicle-stimulating hormone beta (Fshb) and the common alpha-glycoprotein subunit gene (Cga) transcribed in the anterior pituitaries of mammals. Fshb is preferentially transcribed at slower GNRH1 pulse frequencies, whereas Lhb and Cga are preferentially transcribed at more rapid pulse frequencies. Producing the gonadotropins in the correct proportions is critical for normal fertility. Currently, there is no definitive explanation for how GNRH1 pulses differentially activate gonadotropin subunit gene transcription. Several pathways may contribute to this regulation. For example, GNRH1-regulated GNRH1-receptor concentrations may lead to variable signaling pathway activation. Several signaling pathways are activated by GnRH, including mitogen-activated protein kinase, protein kinase C, calcium influx, and calcium-calmodulin kinase, and these may be preferentially regulated under certain conditions. In addition, some signaling proteins feed back to downregulate their own levels. Other arms of gonadotroph signaling appear to be regulated by synthesis, modification, and degradation of either transcription factors or regulatory proteins. Finally, the dynamic binding of proteins to the chromatin, and how that might be regulated by chromatin-modifying proteins, is addressed. Oscillations in expression, modification, and chromatin binding of the proteins involved in gonadotropin gene expression are likely a link between GNRH1 pulsatility and differential gonadotropin transcription.

Using in vitro models for expression profiling studies on ethanol and drugs of abuse

The use of expression profiling with microarrays offers great potential for studying the mechanisms of action of drugs of abuse. Studies with the intact nervous system seem likely to be most relevant to understanding the mechanisms of drug abuse-related behaviours. However, the use of expression profiling with in vitro culture models offers significant advantages for identifying details of cellular signalling actions and toxicity for drugs of abuse. This study discusses general issues of the use of microarrays and cell culture models for studies on drugs of abuse. Specific results from existing studies are also discussed, providing clear examples of relevance for in vitro studies on ethanol, nicotine, opiates, cannabinoids and hallucinogens such as LSD. In addition to providing details on signalling mechanisms relevant to the neurobiology of drugs of abuse, microarray studies on a variety of cell culture systems have also provided important information on mechanisms of cellular/organ toxicity with drugs of abuse. Efforts to integrate genomic studies on drugs of abuse with both in vivo and in vitro models offer the potential for novel mechanistic rigor and physiological relevance.

Multiple and overlapping combinatorial codes orchestrate hormonal responsiveness and dictate cell-specific expression of the genes encoding luteinizing hormone

Normal reproductive function in mammals requires precise control of LH synthesis and secretion by gonadotropes of the anterior pituitary. Synthesis of LH requires expression of two genes [alpha-glycoprotein subunit (alphaGSU) and LHbeta] located on different chromosomes. Hormones from the hypothalamus and gonads modulate transcription of both genes as well as secretion of the biologically active LH heterodimer. In males and females, the transcriptional tone of the genes encoding alphaGSU and LHbeta reflects dynamic integration of a positive signal provided by GnRH from hypothalamic neurons and negative signals emanating from gonadal steroids. Although alphaGSU and LHbeta genes respond transcriptionally in the same manner to changes in hormonal input, different combinations of regulatory elements orchestrate their response. These hormone-responsive regulatory elements are also integral members of much larger combinatorial codes responsible for targeting expression of alphaGSU and LHbeta genes to gonadotropes. In this review, we will profile the genomic landscape of the promoter-regulatory region of both genes, depicting elements and factors that contribute to gonadotrope-specific expression and hormonal regulation. Within this context, we will highlight the different combinatorial codes that control transcriptional responses, particularly those that mediate the opposing effects of GnRH and one of the sex steroids, androgens. We will use this framework to suggest that GnRH and androgens attain the same transcriptional endpoint through combinatorial codes unique to alphaGSU and LHbeta. This parallelism permits the dynamic and coordinate regulation of two genes that encode a single hormone.

Gonadotropin-Releasing Hormone (GnRH) Receptor Expression and Membrane Signaling in Early Embryonic GnRH Neurons: Role in Pulsatile Neurosecretion

The characteristic pulsatile secretion of GnRH from hypothalamic neurons is dependent on an autocrine interaction between GnRH and its receptors expressed in GnRH-producing neurons. The ontogeny and function of this autoregulatory process were investigated in studies on the properties of GnRH neurons derived from the olfactory placode of the fetal rat. An analysis of immunocytochemically identified, laser-captured fetal rat hypothalamic GnRH neurons, and olfactory placode-derived GnRH neurons identified by differential interference contrast microscopy, demonstrated coexpression of mRNAs encoding GnRH and its type I receptor. Both placode-derived and immortalized GnRH neurons (GT1-7 cells) exhibited spontaneous electrical activity that was stimulated by GnRH agonist treatment. This evoked response, as well as basal neuronal firing, was abolished by treatment with a GnRH antagonist. GnRH stimulation elicited biphasic intracellular calcium ([Ca2+]i) responses, and both basal and GnRH-stimulated [Ca2+]i levels were reduced by antagonist treatment. Perifused cultures released GnRH in a pulsatile manner that was highly dependent on extracellular Ca2+. The amplitude of GnRH pulses was increased by GnRH agonist stimulation and was diminished during GnRH antagonist treatment. These findings demonstrate that expression of GnRH receptor, GnRH-dependent activation of Ca2+ signaling, and autocrine regulation of GnRH release are characteristics of early fetal GnRH neurons and could provide a mechanism for gene expression and regulated GnRH secretion during embryonic migration.

Gonadotropin-releasing hormone regulates expression of the DNA damage repair gene, Fanconi anemia A, in pituitary gonadotroph cells

Gonadal function is critically dependant on regulated secretion of the gonadotropin hormones from anterior pituitary gonadotroph cells. Gonadotropin biosynthesis and release is triggered by the binding of hypothalamic GnRH to GnRH receptor expressed on the gonadotroph cell surface. The repertoire of regulatory molecules involved in this process are still being defined. We used the mouse L beta T2 gonadotroph cell line, which expresses both gonadotropin hormones, as a model to investigate GnRH regulation of gene expression and differential display reverse transcription-polymerase chain reaction (RT-PCR) to identify and isolate hormonally induced changes. This approach identified Fanconi anemia a (Fanca), a gene implicated in DNA damage repair, as a differentially expressed transcript. Mutations in Fanca account for the majority of cases of Fanconi anemia (FA), a recessively inherited disease identified by congenital defects, bone marrow failure, infertility, and cancer susceptibility. We confirmed expression and hormonal regulation of Fanca mRNA by quantitative RT-PCR, which showed that GnRH induced a rapid, transient increase in Fanca mRNA. Fanca protein was also acutely upregulated after GnRH treatment of L beta T2 cells. In addition, Fanca gene expression was confined to mature pituitary gonadotrophs and adult mouse pituitary and was not expressed in the immature alpha T3-1 gonadotroph cell line. Thus, this study extends the expression profile of Fanca into a highly specialized endocrine cell and demonstrates hormonal regulation of expression of the Fanca locus. We suggest that this regulatory mechanism may have a crucial role in the GnRH-response mechanism of mature gonadotrophs and perhaps the etiology of FA.

Acute regulation of translation initiation by gonadotropin-releasing hormone in the gonadotrope cell line LbetaT2

The hypothalamic neuropeptide hormone GnRH is the central regulator of reproductive function. GnRH stimulates the synthesis and release of the gonadotropins LH and FSH by the gonadotropes of the anterior pituitary through activation of the G-protein-coupled GnRH receptor. In this study, we investigated the role of translational control of hormone synthesis by the GnRH receptor in the novel gonadotrope cell line LbetaT2. Using immunohistochemical and RIA studies with this model, we show that acute GnRH-induced synthesis and secretion of LH are dependent upon new protein synthesis but not new mRNA synthesis. We examined the response to GnRH and found that activation of cap-dependent translation occurs within 4 h. LHbeta promoter activity was also examined, and we found no increases in LHbeta promoter activity after 6 h of GnRH stimulation. Additionally, we show that increased phosphorylation of translation initiation proteins, 4E-binding protein 1, eukaryotic initiation factor 4E, and eukaryotic initiation factor 4G, occur in a dose- and time-dependent manner in response to GnRH stimulation. Quantitative luminescent image analysis of Western blots shows that 10 nm GnRH is sufficient to cause a maximal increase in factor phosphorylation, and maximal responses occur within 30 min of stimulation. Further, we demonstrate that the MAPK kinase inhibitor, PD 98059, abolishes the GnRH-mediated stimulation of a cap-dependent translation reporter. More specifically, we demonstrate that PD 98059 abolishes the GnRH-mediated stimulation of a downstream target of the ERK pathway, MAPK-interacting kinase. Based on these findings, we conclude that acute GnRH stimulation of LbetaT2 cells increases translation initiation through ERK signaling. This may contribute to the acute increases in LHbeta subunit production.

Transcriptome fingerprints distinguish hallucinogenic and nonhallucinogenic 5-hydroxytryptamine 2A receptor agonist effects in mouse somatosensory cortex

Most neuropharmacological agents and many drugs of abuse modulate the activity of heptahelical G-protein-coupled receptors. Although the effects of these ligands result from changes in cellular signaling, their neurobehavioral activity may not correlate with results of in vitro signal transduction assays. 5-Hydroxytryptamine 2A receptor (5-HT2AR) partial agonists that have similar pharmacological profiles differ in the behavioral responses they elicit. In vitro studies suggest that different agonists acting at the same receptor may establish distinct patterns of signal transduction. Testing this hypothesis in the brain requires a global signal transduction assay that is applicable in vivo. To distinguish the cellular effects of the different 5-HT2AR agonists, we developed an assay for global signal transduction on the basis of high throughput quantification of rapidly modulated transcripts. Study of the responses to agonists in human embryonic kidney 293 cells stably expressing 5-HT2ARs demonstrated that each agonist elicits a distinct transcriptome fingerprint. We therefore studied behavioral and cortical signal transduction responses in wild-type and 5-HT2AR null-mutant mice. The hallucinogenic chemicals (+/-)-2,5-dimethoxy-4-iodoamphetamine (DOI) and lysergic acid diethylamide (LSD) stimulated a head-twitch behavioral response that was not observed with the nonhallucinogenic lisuride hydrogen maleate (LHM) and was absent in receptor null-mutant mice. We also found that DOI, LSD, and LHM each induced distinct transcriptome fingerprints in somatosensory cortex that were absent in 5-HT2AR null-mutants. Moreover, DOI and LSD showed similarities in the transcriptome fingerprints obtained that were not observed with the behaviorally inactive drug LHM. Our results demonstrate that chemicals acting at the 5-HT2AR induce specific cellular response patterns in vivo that are reflected in unique changes in the somatosensory cortex transcriptome.

Validated genomic approach to study differentially expressed genes in complex tissues

Microarray-based genomic techniques allow the simultaneous determination of relative levels of expression of a large number of genes. Studies of the transcriptome in complex neurobiological systems are uniquely demanding due to the heterogeneous nature of these cells. Most brain regions contain a large variety of cell populations that are closely intermingled. The expression of any specific gene may be restricted to a subpopulation of cells, and changes in gene expression may occur in only a small fraction of the cells expressing that transcript. Due to this dilution effect, many genes of interest are expected to have relatively low levels of expression in tissue homogenates. Furthermore, biologically significant differences in expression may result in only small fold-changes. Therefore genomic approaches using brain dissections must be optimized to identify potentially regulated transcripts and differential expression should be confirmed using quantitative assays. We evaluated the effects of increasing tissue complexity on detection of regulated transcripts in focused microarray studies using a mouse cell line, mouse hypothalamus and mouse cortex. Regulated transcripts were confirmed by quantitative real-time PCR. As tissue complexity increased, distinguishing significantly regulated genes from background variation became increasingly more difficult. However, we found that cDNA microarray studies using regional brain dissections and appropriate numbers of replicates could identify genes showing less than 2-fold regulation and that most regulated genes identified fell within this range.