Altered Expression of Genes Encoding Neurotransmitter Receptors in GnRH Neurons of Proestrous Mice

Gonadotropin-releasing hormone (GnRH) neurons play a key role in the central regulation of reproduction. In proestrous female mice, estradiol triggers the pre-ovulatory GnRH surge, however, its impact on the expression of neurotransmitter receptor genes in GnRH neurons has not been explored yet. We hypothesized that proestrus is accompanied by substantial changes in the expression profile of genes coding for neurotransmitter receptors in GnRH neurons. We compared the transcriptome of GnRH neurons obtained from intact, proestrous, and metestrous female GnRH-GFP transgenic mice, respectively. About 1500 individual GnRH neurons were sampled from both groups and their transcriptome was analyzed using microarray hybridization and real-time PCR. In this study, changes in mRNA expression of genes involved in neurotransmitter signaling were investigated. Differential gene expression was most apparent in GABA-ergic (Gabbr1, Gabra3, Gabrb3, Gabrb2, Gabrg2), glutamatergic (Gria1, Gria2, Grin1, Grin3a, Grm1, Slc17a6), cholinergic (Chrnb2, Chrm4) and dopaminergic (Drd3, Drd4), adrenergic (Adra1b, Adra2a, Adra2c), adenosinergic (Adora2a, Adora2b), glycinergic (Glra), purinergic (P2rx7), and serotonergic (Htr1b) receptors. In concert with these events, expression of genes in the signaling pathways downstream to the receptors, i.e., G-proteins (Gnai1, Gnai2, Gnas), adenylate-cyclases (Adcy3, Adcy5), protein kinase A (Prkaca, Prkacb) protein kinase C (Prkca) and certain transporters (Slc1a4, Slc17a6, Slc6a17) were also changed. The marked differences found in the expression of genes involved in neurotransmitter signaling of GnRH neurons at pro- and metestrous stages of the ovarian cycle indicate the differential contribution of these neurotransmitter systems to the induction of the pre-ovulatory GnRH surge, the known prerequisite of the subsequent hormonal cascade inducing ovulation.


INTRODUCTION
Gonadotropin-releasing hormone (GnRH) neurons play a fundamental role in the maintenance of reproduction (Knobil and Neill, 2006). GnRH axons project to the median eminence where they release (Merchenthaler et al., 1980) GnRH into the portal circulation to regulate the pituitary-gonadal axis (Carmel et al., 1976). GnRH neurons are controlled by different neuronal networks via specific membrane receptors for neurotransmitters and neuropeptides released from the presynaptic neuronal afferents (Smith and Jennes, 2001;Campbell, 2007) that modulate the synthesis rate and release pattern of GnRH and other coproduced neuromodulators (Finn et al., 1998;Christian and Moenter, 2010). The operation of the hypothalamo-pituitary gonadal (HPG) axis is cyclic including the cellular activity of GnRH neurons (Plant, 2015). Gonadal hormones have a substantial role in the modulation of GnRH neurons and their neuronal afferents (Radovick et al., 2012).
In female rodents, estradiol (E2) exerts biphasic effects on GnRH neurons and the release of the decapeptide (Sarkar and Fink, 1980;Herbison, 1998). It predominantly suppresses the GnRH system via negative feedback effects with the exception of proestrus when the rising level of E2 primes the system for the preovulatory GnRH surge (Sarkar et al., 1976). This process is driven by E2 acting on estrogen receptors (ERα, ERβ, GPR30 and STX-sensitive membrane receptors; Chu et al., 2009;Terasawa et al., 2009;Kenealy et al., 2011;Moenter and Chu, 2012). While the afferent systems of GnRH neurons are known to be regulated predominantly by ERα (Wintermantel et al., 2006;Yeo and Herbison, 2014;Cheong et al., 2015), GnRH neurons express exclusively the beta subtype of the nuclear receptor (Hrabovszky et al., 2000(Hrabovszky et al., , 2001. The positive E2 feedback, therefore, can target the widespread neuronal regulators of the GnRH system and also the GnRH neurons themselves via direct actions. The neuronal networks mediating the negative and positive feedback effects of E2 to GnRH neurons have been extensively studied by morphological and functional tools (Wintermantel et al., 2006;Yeo and Herbison, 2014). In the preovulatory GnRH surge period, GnRH neurons undergo activation exemplified by expression of the immediate early gene c-Fos (Lee et al., 1990), increased transcriptional activity (Chiu et al., 1988;Wang et al., 1995), induction of hormone synthesis (Gore and Roberts, 1997;Finn et al., 1998) and altered firing pattern (Christian et al., 2005;Farkas et al., 2013).
The classical neurotransmitter systems of the brain are potent regulators of the GnRH system as reviewed earlier (Smith and Jennes, 2001). Powerful regulatory role has been revealed for gamma-aminobutyric acid (GABA) (Herbison and Moenter, 2011), glutamate (Iremonger et al., 2010) dopamine (DA) (Liu and Herbison, 2013), norepinephrine (NE) (Hosny and Jennes, 1998), serotonin (Bhattarai et al., 2014), acetylcholine (ACh) (Turi et al., 2008), and histamine (H) (Fekete et al., 1999). In concert with the rich communication of GnRH neurons with diverse transmitter systems of the brain, the expression of genes encoding for neurotransmitter receptors (Todman et al., 2005) and ion channels (Bosch et al., 2013;Norberg et al., 2013) in GnRH neurons has also been verified. In the present study, the proestrus-associated changes in the expression of neurotransmitter receptor genes of GnRH neurons have been challenged. To achieve this goal, we carried out microarray-and PCR-based transcriptome analysis of GnRH neurons harvested from regularly cycling GnRH-GFP transgenic mice at proestrus and metestrus stages of the ovarian cycle. The comparative study revealed a differential expression of neurotransmitter receptor genes of GnRH neurons in proestrous mice providing novel data to have a better understanding of the communication and plasticity between GnRH neurons and their afferents under the positive feedback action of estradiol.

Ethics Statement
All experiments were performed with permissions from the Animal Welfare Committee of the Institute of Experimental Medicine Hungarian Academy of Sciences (Permission Number: A5769-01) and in accordance with legal requirements of the European Community (Decree86/609/EEC). All animal experimentation described was conducted in accordance with accepted standards of humane animal care and all efforts were made to minimize suffering.

Animals
Adult, gonadally intact female mice were used from local colonies bred at the Medical Gene Technology Unit of the Institute of Experimental Medicine (IEM). They were housed in light (12:12 light-dark cycle, lights on at 06:00 h)-and temperature (22 ± 2 • C) controlled environment, with free access to standard food and tap water. GnRH-green-fluorescent protein (GnRH-GFP) transgenic mice (Suter et al., 2000) bred on a C57BL/6J genetic background were used. In this animal model, a GnRH promoter segment drives selective GFP expression in the majority of GnRH neurons. The estrous cycle was monitored daily between 9 and 10 a.m. by microscopic evaluation of vaginal cytology (Nelson et al., 1982;Byers et al., 2012;Cora et al., 2015). Proestrous (n = 6) and metestrous (n = 6) female mice with at least two consecutive, regular estrous cycles were used. In order to avoid the possible circadian effect, animals were sacrificed at the same period of the day, between 16:00 and 18:00 h. Those animals were considered to be in the proestrus stage that fulfilled the following criteria: (1) vaginal smear staining with predominance of nucleated epithelial cells (Byers et al., 2012); (2) LH serum concentrations >5 mg/L (15.11 ± 3.4 mg/L); (3) uterus wet weights >0.15 g (0.19 ± 0.01 g). Accordingly, the following criteria were applied for the metestrous cycle phase: (1) vaginal smears consisting of the three cell types: leukocytes, cornified and nucleated epithelial cells (Byers et al., 2012); (2) serum LH levels <0.5 mg/L (0.35 ± 0.02 mg/L); (3) uterus wet weights <0.1 g (0.08 ± 0.01 g). Serum LH concentrations and uterus weight data are presented in Supplementary Figure 1.

Serum LH Measurements
Blood samples were collected from the heart of deeply anesthetized mice immediately before the brain fixation step.
The samples were chilled on ice, centrifuged at 1300 g for 3 min at 4 • C. Plasma was aspirated then frozen and stored at −80 • C until further use. Serum LH concentrations were measured with a rodent LH ELISA kit #ERK R7010 from Endocrine Technologies Inc. (Newark, CA, USA) according to manufacturers' instructions.

Laser Capture Microdissection
Brain fixation, preparation of sections for the subsequent laser capture microdissection (LCM) and microarray profiling were performed as reported elsewhere (Khodosevich et al., 2007;Vastagh et al., 2015). Briefly, metestrous (n = 6) and proestrous female (n = 6) mice were deeply anesthetized and perfused transcardially with 80 ml 0.5% paraformaldehyde followed by 20% sucrose. For microdissection, 7 µm thick coronal brain sections were cut. Sections were mounted on PEN-membrane slides (Zeiss, Jena, Germany), processed further for laser microdissection. Uniform and representative sampling of the entire GnRH neuronal population was performed using LCM performed on a PALM Microbeam system (Carl Zeiss Microimaging Gmbh, Jena, Germany) which was equipped with an epifluorescent setup. About 250 GFP-positive neurons were dissected per animal from 80 to 100 consecutive sections to generate GnRH cell samples from each brain.

RNA Isolation
GnRH cell samples (metestrous: n = 6, proestrous: n = 6) collected with LCM were incubated in 200 ml lysis buffer at 56 • C for 3 h. RNA was isolated from the lysate by proteinase K/acid phenol method (Khodosevich et al., 2007), then precipitated by adding isopropanol (Sigma-Aldrich) and 20 µg glycogen (Thermo Fischer Scientific, Waltham, MA, USA) at −20 • C for 30 min followed by centrifugation at 14,000 g at 4 • C. The pellet was washed in 70% ethanol, air-dried and resuspended in water. The remaining genomic DNA was eliminated by treatment the mixture with 1U of RNase-free DNase-I (Thermo Fischer Scientific). RNA was purified using RNeasy MinElute Cleanup kit (Qiagen, Hilden, Germany). Total RNA was eluted with 14 µl of ribonuclease-free water. For the analysis of RNA integrity, RNA was isolated from the tissue of the medial preoptic area, and measured using RNA Pico Chip on the 2100 Bioanalyzer (Agilent, Santa Clara, CA, US; Supplementary Figure 2A).

Whole Transcriptome Amplification (WTA)
Library preparation and amplification were performed according to the manufacturer's instructions for the WTA2 kit (Sigma-Aldrich). In the first step of the WTA2 protocol, the RNA was reverse transcribed implementing non self-complementary and quasi-random 3 ′ and universal 5 ′ primers. The newly synthetized single strands were the templates for annealing and extension steps. When the SYBR Green signal reached a plateau, the reaction was stopped. The resultant library contained cDNA fragments with size distribution between 100 and 1000 base pairs (data obtained from Agilent Bioanalyzer using DNA 1000 microfluidic chip; Supplementary Figure 2B). The amplified double-stranded cDNA was purified and quantified on a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific).

Mouse Genome 430 PM Arrays
Eight microgram of cDNA was fragmented by DNase I and biotinylated by terminal transferase obtained from the GeneChip Mapping 250K Nsp Assay Kit (Affymetrix Inc., Santa Clara, CA, USA). Hybridization, washing, staining, and scanning of Affymetrix Mouse Genome 430 PM Strip arrays were performed following the manufacturer's recommendations. The Mouse Genome 430 PM Strip array allows the analysis of 34,325 wellannotated genes using 45,123 distinct probe sets. Scanned images (DAT files) were transformed into intensities (CEL files) using the AGCC software (Affymetrix). RMA analysis was performed by means of the statistical analysis software Partek Genomics Suite (Partek Inc., St. Louis, MO, USA) to obtain probe set level expression estimates.

Bioinformatics and Pathway Analysis
All statistical and data mining works were performed in R-environment (R Core Team, 2016) with Bioconductor packages (Huber et al., 2015). Quality assessment of microarrays (n = 12) was performed using affyQCReport (Kauffmann et al., 2009). For quality control data regarding RNA integrity, see Supplementary Figure 3. Raw microarray data were preprocessed for analysis using RMA (Robust Multi-Array Average (Irizarry et al., 2003). Fold change (FC) estimation and difference analysis of gene expression were based on linear models combined with Bayesian methods using limma package (Ritchie et al., 2015). FC was calculated from normalized and log 2 transformed gene expression microarray data for each probe sets. The obtained p-values were adjusted by the FDR-based method. The following cut-off criteria were applied on the differentially expressed genes (DEG): |FC| >1.6; and adjusted p-value (p adj ) <0.05. The top 10 differentially expressed genes (up-and downregulated in proestrus, respectively) are shown in Supplementary  Table 1. For over-representation analysis, the list of DEGs was analyzed using the web-based functional annotation tool DAVID Bioinformatics Resources 6.7 (https://david.ncifcrf.gov) at default settings (Huang da et al., 2009a,b).
For visualization of biological pathways, the data matrix of log 2 transformed and RMA-normalized microarray data were analyzed using a public web server Graphite Web (Sales et al., 2013). Of the available analytical methods, Signaling Pathway Impact Analysis (SPIA) was used. Over the classical probabilistic components, the SPIA also applies statistical systems biology approaches: it takes into consideration (1) the extent of the expression changes of each gene (2) the position of the DEG within a pathway, (3) the topology of the pathway, (4) the type of interaction between the genes to identify significantly impacted pathways where the total net accumulated perturbation in the pathway (tA) can be calculated (Tarca et al., 2009) The level of significance was 0.05 using FDR correction. Interacting genes of the significant neurotransmitter signaling pathways (pGFDR < 0.05) were visualized using Cytoscape (ver. 3.2.1) open source software platform.

Quantitative Real-Time PCR Studies
For quantitative real-time PCR (qPCR) investigations of LCMderived GnRH samples (proestrous females n = 6, metestrous females n = 5) RNA isolation and WTA were performed as described in the previous section. Amplified and column-purified cDNA was used as template for qPCR. Whole transcriptomeamplified cDNA from LCM samples were diluted in 0.1x TE buffer for qPCR investigation. Inventoried TaqMan assays were used to confirm microarray results by qPCR. Each assay consisted of a FAM dye-labeled TaqMan MGB probe and two PCR primers. Thermal cycling conditions of the qPCR were as follows: 2 min at 50 • C and 10 min at 95 • C, followed by 40 cycles of 15 s at 95 • C and 1 min at 60 • C using ViiA 7 real-time PCR platform (Thermo Fisher Scientific). Differential expression of genes was calculated by the 2 − Ct method (Livak and Schmittgen, 2001) using GAPDH gene as a reference. Student's t-test were used as a statistical method in comparison of gene expression of the two groups (proestrus: n = 6; metestrus: n = 5).

RESULTS
In this study, we analyzed characteristic changes in neurotransmitter signaling of GnRH neurons-regulation of receptor subunits and downstream elements of the signaling cascades at the level of gene expression-by comparison of the neurotransmitter signaling-related elements of the whole transcriptome of GnRH neurons obtained from intact, metestrous and proestrous GnRH-GFP transgenic mice, respectively.

Differential Expression of Genes Involved in Neurotransmitter Signaling
Analysis of microarray data revealed differentially expressed genes associated with various neurotransmitter signaling mechanisms. Receptor subunits of GABAergic, glutamatergic, and cholinergic communication, as well as receptors participating in adrenergic, serotonergic, dopaminergic, glycinergic, adenosinergic, and purinergic neurotransmission exhibited significant changes (FC > 1.6; adjusted p < 0.05) in gene expression levels. In addition, genes encoding G-proteins and downstream effectors as well as neurotransmitter/amino acid transporters of the solute carrier family also showed differential expression (Figure 1 and Table 1).
Significant downregulation was found in gene expression of adenosinergic (Adora2a, Adora2b), adrenergic (Adra1b, Adra2a, Adra2c) serotonergic (Htr1b), purinergic (P2rx7), and cholinergic (Chrm4, Chrne) receptors or their subunits, respectively, by microarray data analysis. Similarly, the dopaminergic receptor Drd3 and Drd4 exhibited lower mRNA level in proestrus. The Gabbr1 coding for GABA-B receptor 1 and the GABA-A receptor subunit Gabra3, Gabrb2, Gabrb3, Gabrd, Gabrg2 were upregulated. The gene of the GABA-A receptor associated protein, Gabarap were also upregulated. Other GABA-A receptor subunit, the Gabrd showed downregulated gene expression. Receptor subunits of the glutamatergic signaling were upregulated (Grin3a, Gria1, Gria2, and Grin1) whereas Grin2d, the metabotropic glutamate receptor Grm1 and glutaminase (Gls) were downregulated. The expression level of the nicotinic cholinergic receptor beta polypeptide 2 (Chrnb2), as well as several heterotrimeric G-protein alpha subunits (Gnai1, Gnai2, and Gnas) were increased. Genes of downstream primary effector proteins were expressed at either higher (Adcy5) or lower level (Adcy3, Plcb2) in proestrus. The secondary effector Prkaca, Prkacb, and Prkca were upregulated. As far as the solute carrier families are concerned, the glutamate Slc1a4, the vesicular glutamate transporter Slc17a6 and sodium-dependent vesicular transporter Slc6a17 exhibited higher expression level in GnRH neurons in proestrus. Expression of Clcn6 and Glra1 genes, both involved in cellular chloride ion transport, were up-or down-regulated, respectively.

Validation of the Microarray Data
TaqMan real-time PCR was used for validation of microarray data (Figure 2A). Out of the 9 selected target genes, representing GABA, glutamate and acetylcholine neurotransmission, the differential expression of 8 genes was confirmed (Slc17a6, Grin1, Gria1, Chrnb2, Gabrb2, Gabbr1, Gabrb3, Gabra3) and 6 of them exhibited FC >1.5 (Slc17a6, Grin1, Gria1, Chrnb2, Gabrb2, and Gabbr1). We found significant differences in gene expression of GABAB1 and Slc17a6 genes by qPCR. Differences between other targets were not significant due to the high standard deviation of individual samples. The correlation coefficient (Pearson's r = −0.8445) indicates the quantitative interrelation between microarray and real-time PCR data ( Figure 2B).

Pathway Analyses
High throughput analysis of microarray data from 6 proestrous and 6 metestrous females by Signaling Pathway Impact Analysis (False Discovery Rate, FDR of DEGs < 0.05) identified 32 KEGG and 34 REACTOME pathways (adjusted pG using FDR, pGFdr < 0.05). Three of the KEGG pathways (GABAergic, cholinergic, and dopaminergic signaling) and two of the REACTOME pathways (Neurotransmitter receptor binding and downstream transmission in the postsynaptic cell; Activation of NMDA receptor upon glutamate binding and postsynaptic events) were strongly associated with neurotransmitter signaling and/or its postsynaptic events ( Table 2).

Cholinergic Signaling
Changes in gene expression levels of the "cholinergic synapse" (pGFdr = 3.80E-02) indicate an activation (tA = 11.92, Table 2) of this KEGG pathway that involved 29 differentially expressed genes ( Figure 2B). Nicotinic (Chrne) and muscarinic (Chrm4) cholinergic receptors were down-regulated in GnRH neurons, with the exception of the nicotinic acetylcholine receptor subunit Chrnb2 (Table 1) which was up-regulated at proestrus. Upregulation of genes in the signaling cascade (Gnai1, Gnai2, Gnb2, Gng2, Gng3, Gnao1, Prkaca, and Prkacb) was observed. The Gna11 and Gnaq subunits that form the Gq alpha complex and activate the phospholipase C beta 2 (Plcb2), were up-regulated, whereas Plcb2 was down-regulated in proestrus. The expression level of the calcium/calmodulindependent protein kinase II alpha, beta and gamma (Camk2a, Camk2g, and Camk2b, respectively) and Fos were significantly higher in proestrus ( Figure 3B, Supplementary Figure 4 and Supplementary Table 2).

Dopaminergic Signaling
Dopaminergic signaling pathway (pGFdr = 3.55E-02) in GnRH neurons was stated as inhibited (tA = −0.18) in proestrus revealed by SPIA (Table 2, Figure 3C, and Supplementary  Figures 4, 5). According to the microarray data, the expression level of the dopaminergic receptor D3 and D4 (Drd3 and Drd4) was significantly decreased as compared to metestrus. Up-regulation of G-proteins (Gnai1, Gnai2, Gnb2, Gng2, Gng3, and Gnao1) and protein kinases (Prkaca and Prkacb) have been already shown at the GABAergic and cholinergic signaling cascades. Adenylate cyclase (Adcy5), inhibited by the upstream G-protein complex, was up-regulated, whereas Plbc2 activated at protein level, was down-regulated. The glutamatergic pathway interferes with the dopaminergic signals: the upregulation of AMPA 1 and 2 (Gria1, Gria2) receptors are linked with the upregulation of the Akt/Gsk3 pathway: kinesin family members Kif5c, Kif5b, and Kif5 indirectly, whereas Akt3, from the dopaminergic side, directly inhibits Gsk3a and Gsk3b.

Gabaergic Signaling
In the previous decades, it was shown that hypothalamic GnRH neurons receive functional GABAergic primary afferent inputs (Leranth et al., 1985;Adler and Crowley, 1986;Herbison and Dyer, 1991;Kimura and Jinnai, 1994;Jung et al., 1998;Sim et al., 2000;DeFazio et al., 2002;Temple and Wray, 2005;Zhang et al., 2009). The estrous cycle is characterized by profound changes in GABA transmission to GnRH neurons (Herbison, 1997): prior the LH surge the mRNA level of the GABA synthetizing enzyme GAD67 is significantly reduced (Herbison et al., 1992) which exhibits E2-dependent temporal pattern in the territory of the AVPV (Curran- Rauhut and Petersen, 2002). The release of GABA also falls significantly prior to and during the time of estrogen-induced LH surges (Jarry et al., 1992;Tin-Tin-Win-Shwe et al., 2004). This drop of GABA level before the surge seems essential for the positive feedback (Herbison and Dyer, 1991;Kimura and Jinnai, 1994). Expression of the α1-3, α5, β1-3, and γ2 subunits of the GABA A receptor were detected in the GnRH neurons of adult female mice earlier (Sim et al., 2000). We found an increased mRNA expression level of the α1, β1-3, and γ2 subunits and in the proestrous female. Although a GnRH neuron-specific γ2 knock-out mice exhibited normal fertility (Lee et al., 2010), higher level of this subunit can increase the conductance and sensitivity of the GABA A receptors (Brickley et al., 1999). The increased mRNA level of the GABA B receptor 1 and that of several GABA A subunits might also indicate the upregulation of these receptors in response of the dropping GABA level in proestrus. However, the observed negative value of the perturbation score (tA = −4.55) indicates a repressed status of the complete GABAergic synapse pathway in GnRH neurons, probably due to downregulation of genes (effector proteins i.e., Adcy4, Adcy3, Plcb2), downstream to the GABA receptors in the signaling pathway. Recent functional data support the view that net GABA effects might be determined by the balance of the excitatory GABA A (Herbison and Moenter, 2011) and inhibitory GABA B (Zhang et al., 2009;Liu et al., 2011) receptor-mediated tones on GnRH neurons, even in vivo (Constantin et al., 2013). We showed higher expression of both type of the GABA receptors and a repressed status of the GABAergic pathway suggesting an increased sensitivity for the ligand and inhibited signaling cascade that might play a role in the facilitation of the GnRH surge.

Glutamatergic Signaling
GnRH neurons express ionotropic AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) and NMDA (N-methyl-Daspartate) receptors (Gore et al., 1996;Spergel et al., 1999;Ottem et al., 2002) and their antagonists block the LH surge (López et al., 1990Brann and Mahesh, 1991;Ping et al., 1997) indicating a prominent role of glutamate in the positive feedback regulation. At the time of the GnRH/LH surge, the number of contacts that express both glutamate and GABA transporters (Ottem et al., 2004), as well as the level of glutamate increases in the neighborhood of the GnRH cell bodies (Ping et al., 1994;Jarry et al., 1995). The expression of AMPA receptors changes at the level of the somata in steroid-induced LH surge model of young female rats, although, very few (1%) GnRH neurons cosynthetize AMPA receptors and Fos in the morning of proestrus. In contrast, the increase in the number of Gria1/Fos positive (25%) and Gria2/Fos positive (71%) GnRH neurons coincides with the increase and peak of LH (Bailey et al., 2006). In concert with these data, we found that expression of the AMPA1 and AMPA2 receptor subunits (Gria1 and Gria2, respectively) increased by the time of LH surge. Though the percentage of NMDAR1-immunopositive GnRH neurons does not change in the mouse (Adjan et al., 2008), we found that the expression level of NMDA subunits Grin1 and Grin3a increased significantly in proestrus. The higher expression level of these receptor subunits and enrichment of genes in the ionotropic glutamate receptor complex suggest an enhanced glutamate receptor signaling in GnRH neurons during the preovulatory GnRH/LH surge.

Cholinergic Signaling
In the past decades, acetylcholine (ACh) has been shown to exert a stimulatory effect on the release of the GnRH peptide in vitro (Fiorindo and Martini, 1975;Richardson et al., 1982). Cholinergic afferents to hypothalamic GnRH neurons have been detected at both light and electron microscopic levels in the rat (Turi et al., 2008). Studies using perifused hypothalamic and immortalized GnRH (GT1-7) neurons indicated that selective activation of the nicotinic ACh receptor stimulated, in contrast, muscarinic receptor-specific activation inhibited GnRH release (Krsmanovic et al., 1998). In vivo pharmacological experiments in rats suggested an estrous cycle-dependent stimulation of GnRH release by selective muscarinic antagonists (Koren et al., 1992) induced via the M4 receptor subtype encoded by the Chrm4 gene. We found a significantly downregulated Chrm4 expression in proestrus and an upregulated level of Chrnb2. The altered regulation of the two subunits and the diverse action of the nicotinic and muscarinic-type ACh receptor activation can promote GnRH release and surge in proestrus.

Dopaminergic Signaling
GnRH neurons receive tyrosine-hydroxylase (TH) positive terminals demonstrated both at light microscopic (Jennes et al., 1983) and ultrastructural levels (Leranth et al., 1988;Chen et al., 1989;Horvath et al., 1993). Dual phenotype kisspeptin-TH positive fibers of AVPV origin apposed to the somata of GnRH neurons can be one of the sources of these dopaminergic  inputs (Clarkson and Herbison, 2011). Dopamine inhibited the firing and AVPV-evoked GABA/glutamate postsynaptic currents in ∼50% of the GnRH neurons in vitro mediated by D1 and D2like receptors in male and female mice (Liu and Herbison, 2013) showing no evidence of an estrous cycle-dependent modulation when diestrous, proestrous, and estrous female mice were compared in the morning period. In the current study, D2-like dopamine receptors encoded by Drd3 and Drd4 genes were found to be down-regulated in proestrous GnRH neurons, suggesting that expression of these genes are affected by the estrous cycle. However, further functional studies are needed to clarify if downregulation in the expression of the D3 and D4 receptor genes can decrease the inhibitory effect of dopamine on the excitability of the GnRH neurons in the late afternoon period of proestrus.

Adrenergic Signaling
Numerous studies have supported the view for a long time that norepinephrine (NE) exerts a stimulatory effect on GnRH release in the presence of E2 in OVX+E replacement models (Gallo and Drouva, 1979;Ferris et al., 1984;Wise, 1984). Later, NE was found to suppress the excitability of the GnRH neurons acting through both alpha 1 and beta adrenergic receptors with reduced rate of responsive GnRH neurons during proestrus, otherwise exerting a similar effect across the estrous cycle (Han and Herbison, 2008). We found that adrenergic receptors Adra1b, Adra2a, and Adra2c are significantly downregulated in proestrus compared to metestrus. The decreased mRNA expression level of these receptors might indicate an attenuated catecholamine sensitivity of GnRH neurons specifically in the late afternoon period; prior to and during LH surge.

Serotonergic Signaling
Earlier investigations provided evidence that 5-HT neurons project directly to GnRH neurons in rodents (Jennes et al., 1982;Kiss and Halász, 1985;Wada et al., 2006;Campbell and Herbison, 2007) In the mouse, serotonin (5-HT) exerts a biphasic effect on the excitability of GnRH neurons: both inhibitory and excitatory responses have been demonstrated, mediated by activation of 5-HT1A and 5-HT2A receptors, respectively (Bhattarai et al., 2014). The serotonergic input to GnRH neurons showed an estrous cycle-dependent regulation in the adult female mouse: during proestrus, the serotonergic inhibition was reduced, accompanied by the absence of the biphasic responses. The 5-HT1B receptor also mediates inhibitory neurotransmission (Mathur et al., 2011;Huang et al., 2013) in the CNS. Although the function of the 5-HT1B receptor was not investigated in GnRH neurons to date, we found that it was down-regulated in proestrous GnRH neurons. Compared to other stages of the estrous cycle, 5-HT1 receptor expression is decreased up to 40% in the basal forebrain during proestrus (Biegon et al., 1980) when peaking level of E2 significantly down-regulates 5-HT receptors (Biegon and McEwen, 1982). The current study further supports that a suppression of the inhibitory 5-HT tone takes place in GnRH neurons in proestrus.

P2X Receptors
P2X receptors are trimeric, ligand-gated cation channels activated upon binding of extracellular ATP. They are highly permeable to Ca 2+ and expressed in diverse organs including the brain (Soto et al., 1996;Burnstock, 2013). Purinergic signaling plays a fundamental role in the regulation of hypothalamopituitary functions (Stojilkovic, 2009). Expression of P2X2 and P2X4 has been confirmed in GnRH neurons in olfactory placode cultures (Terasawa et al., 2005) where application of ATP helps to synchronize [Ca 2+ ] i oscillations. Furthermore, P2X5 and P2X6 receptors are also expressed in GnRH neurons of mice as revealed by double-labeling immunofluorescence (Fu et al., 2009). Interestingly, though the presence of P2X7 receptor has not been confirmed in GnRH neurons to date, the mRNA level of P2xr7 is significantly down-regulated in proestrous GnRH neurons. A possible explanation of this phenomenon may be the E2-dependent regulation of P2X7R expression during the estrous cycle. Data from literature indicate that treatment with E2 decreases the expression of P2X7 receptor both in vitro (Cario-Toumaniantz et al., 1998) and in vivo (Xu et al., 2016).

P1 Receptors
The adenosine receptor 2A and 2B receptors are members of the G protein-coupled receptor (GPCR) family. The A2 adenosine receptors are coupled to G s and G olf family of G proteins stimulating adenylate cyclase activity to increase the intracellular cAMP level (Schwindinger et al., 2010). Expression of Adora2b has been shown in hypothalamic GnRH neurons (Todman et al., 2005). In this study we found significant down-regulation of the Adora2a and Adora2b transcripts in proestrous GnRH neurons. The role and significance of the regulation of adenosine receptor 2 are unknown at present, they may be linked to the maintenance of intracellular cAMP levels along the estrous cycle.

Methodological Considerations
Sampling by LCM LCM is among the best approaches available for cell-type-specific microarray gene expression profiling (Lin et al., 2007;Okaty et al., 2011;Demarest et al., 2012). This methodology helps to preserve the integrity of the GnRH transcriptome and ensures the systematic sampling of the hypothalamic GnRH neurons. Our strategy aimed to minimize the extent of contamination from surrounding cells was as follows: (1) brain sections were cut at 7 µm thickness; (2) during microdissection, the cutting laser were directed strictly along the perikaryon of the GnRH-GFP neurons. While it is not possible to eliminate contamination during sample preparation completely, careful implementation of the steps above results in highly enriched mRNAs isolated from GnRH neurons.

Linearity of the RNA Amplification
In this work we implemented a relatively new approach for gene expression profiling validated and published earlier (Gonzalez-Roca et al., 2010) showing that "pico profiling" allows an accurate measure of transcript levels from populations as low as 10 cells with insignificant number of false positive or negative hits, using a commercially available TransPlex WTA2 kit from Sigma-Aldrich. Importantly, at the reverse transcription step the quasi-random 3 ′ and universal 5 ′ primers resulted in DNA fragments flanked by universal end sequences to ensure the linear amplification of the expressed genes without 3 ′ and 5 ′ bias.

Validation of the Data
High throughput gene expression analysis of GnRH neurons across the estrous cycle was performed for the first time in the present study. A previous study has addressed ion channel expression and E2 regulations in GnRH neurons during negative and positive feedback (Bosch et al., 2013). They showed that E2 treatment up-regulates voltage-gated calcium channel (VGCC) Cav1.3, Cav2.2, and Cav2.3 mRNAs in GnRH neurons. In consistence with their results regarding the expressional changes in the positive feedback period, we have also found a significant increase in the expression of the Cav 2.2, the voltage-dependent N-type calcium channel alpha 1B subunit encoded by Cacna1b gene (FC 1.7; p adj < 0.01) in proestrus whereas the increase in mRNA expression the Cav1.3 [Cacna1d] and Cav2.3 [Cacna1e] subunits was not significant. This discrepancy may be due to the different experimental paradigm: Bosch et al. used ex-vivo sorted GnRH neurons from OVX mice primed with E2 dose followed by a surge-inducing E2 dose, whereas in the current study naturally cycling intact mice were used.

Correlation of mRNA and Protein Abundances
In this study, differential expression of genes was investigated as alterations in the mRNA abundances. In general, changes of the protein levels should not be well correlated with that of the mRNA levels due to many factors that influence the expression of proteins (ribosomal density, protein half-life, etc.) and the experimental error and noise (Maier et al., 2009;Vogel and Marcotte, 2012). Genome-wide correlations between mRNA and protein expression levels are usually weak (de Sousa Abreu et al., 2009). However, the level of mRNAs expressed differentially significantly correlate better with their proteins compared to genes that are not expressed differentially (Koussounadis et al., 2015) further increasing the confidence of studies applying differential mRNA expression measurements.

CONCLUSION
We found considerable differences in the expression of genes encoding neurotransmitter receptors and their effectors in downstream signaling cascades by comparison of GnRH neurons obtained from pro-and metestrous mice, respectively. Signaling systems known to facilitate the GnRH surge (glutamatergic, GABA-A receptor-mediated inputs and cholinergic neurotransmission via nicotinic receptors) become activated, while neuronal inputs that exert inhibitory effects on GnRH release (dopaminergic, serotonergic, adrenergic systems) seem to be negatively regulated at the level of transcripts in proestrus. These complex changes in the gene expression of proestrous GnRH neurons may alter diverse intracellular mechanisms that culminate in the preovulatory GnRH surge.

AUTHOR CONTRIBUTIONS
CV designed and performed the experiments, analyzed the data and wrote the manuscript. AR executed hybridization and scanning of microarrays. NS carried out the bioinformatical analysis of the microarray data. ZL designed and supervised the project, and wrote the manuscript.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fncel. 2016.00230