Divergent activity of the gonadotropin-releasing hormone receptor gene promoter among genetic lines of pigs is partially conferred by nuclear factor (NF)-kB, specificity protein (SP)1-like and GATA-4 binding sites
- Emily A. McDonald†1, 2,
- Jacqueline E. Smith†1, 3,
- Rebecca A. Cederberg1 and
- Brett R. White1Email author
© The Author(s). 2016
Received: 2 June 2016
Accepted: 22 June 2016
Published: 29 June 2016
Binding of gonadotropin-releasing hormone (GnRH) to its receptor (GnRHR) on gonadotropes within the anterior pituitary gland is essential to reproduction. In pigs, the GnRHR gene is also located near a genetic marker for ovulation rate, a primary determinant of prolificacy. We hypothesized that pituitary expression of the GnRHR gene is alternatively regulated in genetic strains with elevated ovulation rates (Chinese Meishan and Nebraska Index) vs. standard white crossbred swine (Control).
Luciferase reporter vectors containing 5118 bp of GnRHR gene promoter from either the Control, Index or Meishan swine lines were generated. Transient transfection of line-specific, full length, deletion and mutation constructs into gonadotrope-derived αT3-1 cells were performed to compare promoter activity and identify regions necessary for divergent regulation of the porcine GnRHR gene. Additionally, transcription factors that bind the GnRHR promoter from each line were identified with electrophoretic mobility shift assays (EMSA).
Dramatic differences in luciferase activity among Control, Index and Meishan promoters (19-, 27- and 49-fold over promoterless control, respectively; P < 0.05) were established. A single bp substitution (−1690) within a previously identified upstream enhancer (−1779/−1667) bound GATA-4 in the Meishan promoter and the p52/p65 subunits of nuclear factor (NF)-κB in the homologous Control/Index promoters. Transient transfection of vectors containing block replacement mutations of either the GATA-4 or NF-κB binding sites within the context of their native promoters resulted in a 50 and 60 % reduction of luciferase activity, respectively (P < 0.05). Furthermore, two single-bp substitutions in the Meishan compared to Control/Index promoters resulted in binding of the p52 and p65 subunits of NF-κB and a specificity protein 1 (SP1)-like factor (−1235) as well as GATA-4 (−845). Vectors containing the full-length Meishan promoter harboring individual mutations spanning these regions reduced luciferase activity by 25 and 20 %, respectively, compared to native sequence (P < 0.05).
Elevated activity of the Meishan GnRHR gene promoter over Control/Index promoters in αT3-1 cells is partially due to three single nucleotide polymorphisms resulting in the unique binding of GATA-4 (−1690), the p52/p65 subunits of NF-kB in combination with a SP1-like factor (−1235), and GATA-4 (−845).
KeywordsGnRH receptor Transcriptional regulation NF-κB GATA-4 SP1-like factor Porcine Single nucleotide polymorphism Ovulation rate Anterior pituitary Gonadotrope
The hypothalamic decapeptide, gonadotropin-releasing hormone (GnRH), binds to its cognate receptor (GnRHR) on the surface of gonadotrope cells within the anterior pituitary gland, stimulating the synthesis and secretion of the gonadotropins, follicle stimulating hormone (FSH) and luteinizing hormone (LH) [1, 2]. The gonadotropins subsequently act on the gonads to trigger secretion of steroid hormones, which feedback at the level of the hypothalamus and anterior pituitary gland to regulate GnRH and gonadotropin levels, respectively. Binding of GnRH to its seven transmembrane, G-protein-coupled receptor activates multiple signal transduction cascades, ultimately resulting in up-regulation of the genes that encode the common α- and unique β-subunits of FSH and LH [3, 4] as well as the GnRHR itself . Therefore, the interaction between GnRH and GnRHR represents a crucial point for regulation of reproductive function in mammals. Furthermore, the porcine GnRHR gene is located on chromosome 8, in close proximity to a quantitative trait locus for ovulation rate, a primary determinant of litter size . Consistent with this, a C/G substitution in the 3’ untranslated region was shown to be significantly associated with ovulation rate . Consequently, the GnRHR gene represents both a physiological and positional candidate for genes influencing prolificacy in pigs.
The GnRHR gene promoter has been extensively studied in the mouse, rat, human and sheep [8, 9]. In the mouse, gonadotrope-specific expression is conferred by 500 bp of GnRHR gene promoter  comprised of binding sites for steroidogenic factor 1 (SF1), activator protein 1 (AP1) and a GnRHR activating sequence (GRAS) [11, 12], although pituitary homeobox (Pitx)-1 and a member of the LIM homeodomain family, Lhx3, have also been implicated [13, 14]. Additionally, protein kinase C (PKC) activation of AP1 is critical for GnRH responsiveness of the murine GnRHR gene promoter , whereas the GRAS binding site acts in conjunction with a downstream activin regulatory element to control activin responsiveness . Subsequent studies revealed that the GRAS element binds a complex of transcription factors including AP1, Smad3 and 4 and FOXL2, a member of the forkhead family of transcription factors . Furthermore, an enhancer element, sequence underlying responsiveness to GnRH (SURG)-1 , binds octamer transcription factor-1 (OCT1) and nuclear factor (NF)-Y  for basal and maximal GnRH stimulation of GnRHR gene transcription. A proximal homeodomain (Hbox) binding motif also binds OCT1, indicating the transcription factor acts at multiple TAAT sites to direct basal expression . Interestingly, the CLOCK and BMAL1 drive activity of the murine GnRHR promoter through their interaction with E-box enhancer sequences . The importance of protein kinase A (PKA) signaling was also illustrated by the role of a cAMP responsive element (CRE) in activation of the GnRHR gene .
The rat and mouse GnRHR gene promoters appear consistent, both containing SF1, AP1 and GRAS binding sites . However, the GRAS element in the rat GnRHR promoter harbors an A → G bp alteration, dramatically reducing effectiveness of this binding site . Furthermore, the rat promoter contains CRE-like and SF1 adjacent protein (SAP) binding sites involved in basal activity  and a GnRHR specific enhancer (GnSE)  that interacts with GATA factors and the LIM-related factors, Isl-1 and Lhx3, to promote maximal basal activity . Although activity of the human GnRHR promoter in cell lines of non-gonadotrope origin and in response to hormones have been characterized , investigation into gonadotrope-specific activity has only elucidated a SF1 binding site . Additional studies revealed that AP1 confers down-regulation following GnRH stimulation  and OCT1 serves as a strong constitutive repressor . Elements conferring gonadotrope-specific expression of the ovine GnRHR gene remain to be elucidated; however, Duval and coworkers  reported a SF1 binding site that mediates basal expression.
Our laboratory has previously shown that gonadotrope-specific activity of the porcine GnRHR gene is partially conferred by a SF1 binding site positioned within a 112-bp upstream enhancer (−1779/−1667), as well as two additional SF1 and one retinoid X receptor (RXR) binding sites located within 315 bp of proximal promoter . In order to compare transcriptional regulation of the GnRHR gene among pig lines with divergent ovulation rates, we constructed luciferase reporter constructs containing 5118 bp of 5’ flanking sequence from three genetic lines of swine: a Control white-crossbred line; a Nebraska Index line selected for over 14 generations based on an index of ovulation rate and embryonic survival ; and the Chinese Meishan breed, a line with increased prolificacy over white-crossbred lines, largely due to a greater ovulation rate [32, 33]. Previously, our laboratory has shown that anterior pituitary levels of GnRHR mRNA were highest in Meishan, intermediate in Index and lowest in Control . Herein, we demonstrate differential activity among these line-specific GnRHR promoters utilizing transient transfections assays in gonadotrope-derived αT3-1 cells. In addition, we identified three bp substitutions at −1690 (T → C), −1235 (C → G) and −845 (G → T) of proximal promoter that allow GATA-4, the p52 and p65 subunits of nuclear factor (NF)-κB as well as a specificity protein (SP)1-like factor, and GATA-4, respectively, to preferentially bind the Meishan compared to Index or Control GnRHR gene promoters.
Experiments involving the use of recombinant DNA have been approved by the UNL Institutional Biosafety Committee under Protocol ID # 12 entitled: Functional Analysis of GnRHR I and II in Swine. The UNL Radiation Safety Office has approved the use of isotopes in the following experiments via AU License # I-387.
Sense strand of EMSA oligosa
Primers used to generate reporter vectors
proximal A OF
proximal B OF
proximal C OF
Cell culture and transient transfections
Cultures of αT3-1 cells (Dr. Pam Mellon, Salk Institute, La Jolla, CA) were maintained at 37 °C in a humidified 5 % CO2 in air atmosphere. The αT3-1 cells were cultured in high-glucose DMEM (4.5 g/L; Mediatech, Herndon, VA) supplemented with 5 % fetal bovine serum, 5 % horse serum, 2 mM glutamine, 100 U/ml penicillin and 100 μg/ml streptomycin sulfate (Gibco, Grand Island, NY). Transient transfections were carried out using a liposome-mediated protocol (Fugene6, Roche Diagnostics Corp., Indianapolis, IN) according to manufacturer’s instructions. Briefly, 2 × 106 cells were plated in 6-well culture dishes 1 d prior to transfection. Cells were transfected with a 3:1 Fugene6 to DNA ratio. A total of 1 μg of DNA, 0.9 μg of luciferase test vector and 0.1 μg of RSV-βgal control vector were used per well. Approximately 20–24 h post-transfection, cells were washed twice with ice-cold PBS and harvested with 200 μl of lysis buffer [100 mM potassium phosphate (pH 7.8), 0.2 % Triton X-100 and 1 mM dithiothreitol (DTT)]. Lysates were cleared by centrifugation at 14,000 X g for 2 min at 4 °C. Lysates (20 μl) were immediately analyzed according to manufacturer’s instructions for both luciferase (Promega Corp.) and β-gal (Applied Biosystems, Bedford, MA) activity using a Wallac Victor2 microplate reader (PerkinElmer Life Sciences, Boston, MA). Luciferase values were divided by β-gal values to adjust for transfection efficiency. The raw data for all transfections utilized in this study have been included (Additional file 1).
Nuclear protein extracts were obtained from approximately 2.8 × 108 αT3-1 cells using the NE-PER® Nuclear and Cytoplasmic Extraction Reagents Kit (Pierce Biotechnology, Rockford, IL). The nuclear extracts were treated with protease (catalog no. P8340; Sigma Chemical Co., St. Louis, MO) and phosphatase (catalog no. 524625; Calbiochem, La Jolla, CA) inhibitor cocktail solutions to prevent enzymatic degradation of proteins. The amount of protein present in the extracts was determined using bicinchoninic acid (BCA assay, Pierce Biotechnology). Oligonucleotides were end-labeled with [γ-32P]ATP using T4 polynucleotide kinase (Fermentas Inc., Hanover, MD) and purified using sephadex G-25 spin columns (Amersham Biosciences Corp., Piscataway, NJ). EMSAs were completed through incubation of nuclear extracts (5 μg) in 20 μl reactions containing 4 μl of Dignam D buffer (20 mM HEPES, 20 % glycerol, 0.1 M potassium chloride, 0.2 mM EDTA and 0.5 mM DTT), 1 mM DTT, 2 μg of poly(dI•dC) (Amersham Biosciences) and, where indicated, a rabbit polyclonal antibody directed against the p65 (Calbiochem), p52 (Upstate, Lake Placid, NY), or p50 (Santa Cruz Biotechnology) subunits of NF-κB; GATA-1, -2 and -4 (Santa Cruz Biotechnology); SP1 (Upstate), SP1, SP2, SP3 and SP4 (Santa Cruz Biotechnology) or an equal mass of rabbit IgG (Santa Cruz Biotechnology) at 4 °C for 2 h. Following incubation, radiolabeled probe (100,000 cpm) and 50-fold molar excess of either homologous or heterologous unlabeled competitor was added. Where indicated, 50-fold molar excess of unlabeled oligonucleotides containing consensus binding sequences for AP2, NF-κB, SP1, glucocorticoid receptor (GR), nuclear factor (NF)-1 or GATA-4 were also added. The final reactions were incubated at 25 °C for 20 min before bound probe was separated from free at 30 mA for 1.5 h on a 5 % polyacrylamide gel that had been prerun at 100 V for 1 h in 1X TGE [25 mM Tris (pH 8.3), 190 mM glycine and 1 mM EDTA]. Gels were transferred to blotting paper, dried, and exposed to Biomax MS film (Eastman Kodak Co., Rochester, NY) for 20–24 h at −80 °C before being developed.
Nuclear proteins from αT3-1 cells were extracted using the Nuclear Complex Co-IP kit from Active Motif (Carlsbad, CA), quantitated with a BCA protein assay kit (Pierce) and stored at −80 °C. Protein samples (40 μg) were boiled for 5 min in a 2X reducing loading buffer (130 mM Tris pH 6.8, 4 % SDS, 0.02 % Orange G, 20 % glycerol, 100 mM DTT), cooled to room temperature (RT) and loaded onto an SDS polyacrylamide gel (PAGE) with a 5 % stacking and 10 % resolving gel. Gels were run at 40 mA for approximately 90 min and electrophoresed proteins were transferred to polyvinylidene difluoride (PVDF, Immobilon -FL, Millipore, Billerica, MA) membrane with a semi-dry electroblotter (Panther, Owl Separation Systems, Portsmouth, NH). Briefly, PVDF membrane was pre-wetted in 100 % methanol and soaked with the gel in transfer buffer (25 mM Tris pH 8.3, 192 mM glycine, 0.1 % SDS, 20 % methanol) for 15 min. The proteins were transferred at 200 mA for 1 h. Membranes were blocked with StartingBlock™ TBS buffer (Pierce) for 30 min at room temperature with agitation. Incubation of primary antibodies directed against the p50 (Santa Cruz Biotechnology), p52 (Upstate) and p65 (Calbiochem) subunits of NF-κB were performed in StartingBlock™ TBS buffer supplemented with 0.05 % Tween-20. Antibodies were used at 1:500 (p50) or 1:5000 (p52 and p65) dilutions. Blots were incubated with primary antibody overnight at 4 °C with gentle shaking. After incubation, the blots were washed four times with TBS-T (20 mM Tris pH 7.6, 137 mM sodium chloride, 0.1 % Tween-20). Each wash was performed for 5 min with gentle agitation. The secondary antibody, Alexa Fluor 680 goat anti-rabbit IgG (A21076, Invitrogen, Carlsbad, CA) was diluted 1:15,000 in StartingBlock™ TBS buffer (Pierce) supplemented with 0.01 % SDS and 0.05 % Tween-20. The incubation was performed at RT for 1 h with gentle shaking. Blots were washed four times with TBS-T for 5 min with gentle agitation. After a final rinse with TBS, blots were scanned on the 700 channel of the Odyssey Infrared Imaging System (LI-COR, Lincoln, NE) following manufacturer’s instructions.
Data were analyzed using the general linear models (GLM) procedure of the Statistical Analysis System (SAS, version 8.2, Cary, NC). To control for transfection efficiency, the arbitrary light value for each replicate was divided by the respective β-gal value. These values were then divided by the mean of the empty vector and reported as fold activity over pGL3. All transfections were performed a minimum of three times, with samples in triplicate using different plasmid preparations for each transfection. Individual values from all the replicates were used to generate the mean ± SEM. Comparisons between pGL3 and test vectors were evaluated with Dunnett’s t-test. Least squares means for luciferase activity were compared among test vectors using least significant differences.
The −5118 bp Control, Index and Meishan GnRHR promoters display divergent activity in αT3-1 cells, which is maintained upon promoter reduction to −1915 bp
Previous work in our laboratory indicated that gonadotrope-specific expression of the porcine GnRHR gene was conferred by elements located within 1915 bp of proximal promoter . To determine if this region contained elements conferring line-specific expression, we reduced the 5’ flanking region for all three pig lines from −5118 to −1915 bp. Transient transfection of plasmids containing either −5118 (described above) or −1915 GnRHR promoters for the Meishan (−M1915pGL3), Index (−I1915pGL3) and Control (−C1915pGL3) lines resulted in maintenance of significant differences in luciferase activity among lines (Fig. 1). Therefore, the elements conferring line-specific expression of the porcine GnRHR gene reside within 1915 bp of proximal promoter.
Sequence alignment of the Control, Index and Meishan GnRHR promoters identified polymorphisms at several locations within −1915 bp of promoter
Base changes between lines within the −1915 bp promotersa
C to Tb,c
A to G
T to C
G to Tb
T to C
T to G
C to T
A to Gc
C to Gc
A to Cc
C to Ac
A to Gc
T to Cc
22 bp deletionc
A to T
A to G
Differential GnRHR promoter activity among genetic lines is not attributable to a SF1 binding site, an essential element to basal activity of the promoter in αT3-1 cells
The region between −1915/−1431 of the porcine GnRHR gene promoter is not alone responsible for line divergent promoter activity
Enhanced activity of the Meishan GnRHR gene promoter is diminished upon reduction of the 5’ flanking region to −1004 bp
A T → C bp substitution at −1690 confers divergent binding between the Meishan and Control/Index GnRHR gene promoters
Block replacement mutations of both the GATA-4 and NF-κB binding sites in the Meishan and Control/Index promoters, respectively, attenuated luciferase activity
A C → G bp substitution at −1235 promotes binding of transcription factors unique to the Meishan GnRHR promoter
Block replacement mutations of the NF-κB and SP1 binding sites located between −1240/−1230 of 5’ flanking region reduce activity of the Meishan GnRHR promoter
A third bp substitution (G → T) located at −845 bp confers divergent GATA-4 binding to the Meishan GnRHR promoter
Block replacement mutation of the GATA-4 binding site at −845/−840 of 5’ flanking region reduces activity of the Meishan GnRHR promoter
Here, we have demonstrated a dramatic variance in activity among reporter vectors containing GnRHR promoters from the Control, Index and Meishan lines of pigs using transient transfections in the gonadotrope-derived αT3-1 cell line. Divergent activity of the Control, Index, and Meishan promoters illustrates a potential alteration in the mechanisms underlying transcriptional regulation of the porcine GnRHR gene among genetic strains. Ultimately, differential regulation of GnRHR gene expression may be correlated with divergent ovulation rates as well as other reproductive traits observed among the Control, Index and Meishan pig lines. This characteristic difference in activity among pig strains was lost after promoter constructs were reduced from −1431 to −1004 bp of 5’ flanking region, inferring that the element(s) responsible for the elevated luciferase activity of the Meishan promoter construct is located within the proximal 1431 bp of GnRHR promoter. Within gonadotrope-derived cell lines, the majority of elements responsible for basal and hormonally-induced expression of the GnRHR gene in other species are located within 1000 bp of proximal promoter [10–19, 22–28, 36, 37], although placental-, granulosa/luteal cell- and neuronal-specific promoters have been identified further upstream in the human [38–40]. On the other hand, the spatial arrangement of the porcine GnRHR promoter is somewhat unique because it requires approximately 1800 bp of 5’ flanking sequence for basal activity in αT3-1 cells . Despite the fact that the elements which confer line-specific activity of the porcine GnRHR promoter lie further upstream than 1000 bp, they still remain within the boundaries of the gonadotrope-specific promoter.
Transcriptional regulation of the GnRHR gene in different species is achieved through a variety of mechanisms and a number of different transcription factors [8, 9]. One of the most characterized mechanisms for transcriptional regulation of the GnRHR gene involves the orphan nuclear receptor, SF1, known to be vital for gonadotrope-specific expression of the GnRHR gene in the human , mouse , rat , sheep  and pig . SF1 binding is also known to confer expression of the gonadotropin subunit genes within gonadotrope cells [41–43]. Previously, our laboratory demonstrated that three SF1 binding sites are involved in transcriptional regulation of the porcine GnRHR gene promoter in αT3-1 cells. In reporter assays, mutation of the proximal SF1 binding site at −179/−171 bp resulted in complete ablation of promoter activity, indicating that this site is essential for gene expression . In this study, however, we have shown that the critical proximal SF1 binding site is not responsible for line-specific regulation of the GnRHR gene in swine (Fig. 2).
Identification of an NF-κB site involved in line-specific expression of the porcine GnRHR gene in gonadotrope cells is unique. While NF-κB has traditionally been associated with genes involved in immune and inflammatory responses , it has been implicated in mediating the apoptotic effects of GnRH in ovarian cancer cells . As reviewed by Hayden and Ghosh , NF-κB typically exists as a heterodimer composed of p65 (RelA), p50, p52, RelB and/or c-Rel. Although p52/p65 and p52/RelB heterodimers are prevalent, most commonly the heterodimer is comprised of the p65 and p50 subunits. Previously we reported that p52/p65 subunits are involved in the transcriptional regulation of the porcine GnRHR2 gene in the testis . Regarding gonadotrope cells of the anterior pituitary gland, however, NF-κB remains noticeably absent from the cast of transcription factors known to regulate the GnRHR or gonadotropin subunit genes. In this study, we detected p52 and p65 proteins within αT3-1 nuclear extracts and established the role of p52/p65 heterodimer binding in transcriptional regulation of the GnRHR gene. Consistent with this, NF-κB is an important regulator of Cox-2 promoter activity in the gonadotrope-derived cell line, LβT2, and treatment with GnRH stimulated phosphorylation of the p65 subunit by 22-fold . In addition, NF-κB binding sites have also been associated with regulation of genes expressed in other pituitary cell types including somatotropes  and corticotropes .
The −1240/−1230 region of the Meishan GnRHR promoter also appears to bind another transcription factor capable of interacting with the SP1 consensus binding site. Although the use of an SP1 consensus oligonucleotide revealed competition for DNA-protein binding (Fig. 4), the inability of SP1-specific antibodies acquired from two separate commercial vendors to bind the DNA-protein complex (Fig. 4) suggests that an alternative transcription factor recognizes the SP1 consensus binding sequence. To further verify the integrity of the antibodies directed against SP1, we also performed EMSAs with αT3-1 nuclear extracts and radiolabeled oligonucleotide containing consensus binding sites for SP1. In this instance, a specific complex was formed and the addition of SP1-specific antibodies resulted in supershifted DNA-protein complexes (data not shown), confirming the presence of SP1 in nuclei of αT3-1 cells and effectiveness of the SP1 antibodies. Another concern was whether the SP1 consensus oligonucleotides actually contained one or more NF-κB elements and therefore, was merely mimicking the NF-κB consensus oligonucleotide. However, sequence analysis of the oligonucleotide containing consensus SP1 binding sites did not reveal any NF-κB elements. Next, we examined other members of the SP1 family of transcription factors. Of primary interest, SP3 and SP4 recognize the same binding sequence with similar affinities . However, the addition of antibodies directed against SP2, SP3 and SP4 were unable to bind to the specific complex. Thus, we were able to eliminate SP1-4 as potential transcription factors binding to the SP1 element located within −1240/−1230 of the Meishan GnRHR promoter. Currently, nine members of the SP1 family (SP1-9) have been identified [52, 53], and represent a subgroup of a larger class of transcription factors, the SP1-like/Krüppel-like factor (KLF) family . These factors share a highly conserved DNA-binding domain containing three Cys2/His2 zinc fingers , which is the most abundant transcription factor motif in the human genome . In fact, over 25 SP1-like/KLF genes have been reported in mammals . The nine SP transcription factors can further be divided into the SP1-like family (SP1-4) and SP8-like family (SP5-9) . While SP1 and SP3 are ubiquitously expressed, SP5-9 expression patterns are more specific and temporal. SP5 and SP8 expression has been linked with Wnt activity and are important for stem cell differentiation and early embryonic development [58, 59]. SP6, also known as epiprofin, is involved in epidermal differentiation  whereas SP7, known as osterix, is expressed in osteoblasts . Additionally, SP9 is also expressed embryonically and affects limb outgrowth . Given the developmental roles that SP5-9 frequently play, it is unlikely that any of these factors are binding the GnRHR gene promoter. Thus, more studies are required to determine the identity of the factor(s) binding to the SP1 site within the −1240/−1230 promoter region of the Meishan GnRHR gene.
Mutation of the individual binding sites for NF-κB and the SP1-like factor within the Meishan GnRHR promoter demonstrated a significant loss of promoter activity. In addition, mutation of both the NF-κB and the SP1 recognition sequences (−Mμ1240pGL3) diminished luciferase activity to approximately the same level as either single block mutation. Due to the lack of further reduction of luciferase activity by the double block replacement mutation, it would appear that the transcription factors binding to the two elements work synergistically, but not additively, to stimulate increased GnRHR promoter activity in the Meishan line of swine. Indeed, NF-κB interacts with a variety of other transcription factors including AP1, estrogen receptor α, C/EBP, SF1 and SP1 [63–66]. The Nabel laboratory determined that SP1 interacts with p65 via the DNA binding regions of each factor and that this interaction is necessary for activation of the HIV-1 gene . Despite the importance of the complex of transcription factors binding at −1240/−1230, none of the mutations lowered Meishan GnRHR promoter activity to that of the native full-length Control promoter suggesting another element(s) within the Meishan GnRHR promoter also confers line-specific expression.
We also identified 2 additional single-bp alterations located at −845 (G → T) and at −1690 (T → C) within the Meishan GnRHR promoter region that allowed binding of GATA-4 to the recognition sites. While GATA-4 has not previously been implicated in transcriptional regulation of the GnRHR gene, it is involved in the expression of other gonadotropic genes. A GATA motif detected within the human α-subunit gene promoter binds GATA-2 and a GATA-4-related protein in αT3-1 cells . These investigators were unable to confirm GATA-4 binding because specific antibodies directed against GATA-4 were not commercially available at the time. Consistent with these results, our study confirmed the presence of GATA-4 in gonadotrope-derived αT3-1 cells and implicated its importance in regulation of gonadotropic gene expression. Steger and coworkers  also reported that the same GATA element binds GATA-2 and -3, but not GATA-4-related protein, in placental-derived cell lines. Binding of GATA-4 also regulates other genes essential to reproduction . In the neuronal GT1-7 cell line, GATA-4 binds recognition sites within the promoter for the GnRH gene . Transcription of the Müllerian inhibiting substance (MIS) gene in Sertoli cells is enhanced by the direct interaction of GATA-4 and SF1, although GATA-4 binding to the DNA is not required for this synergistic effect [71, 72]. Additionally, adrenal-specific transcription of the human P450c17 gene is regulated by the interaction of GATA-4 or GATA-6 with SP1 . However, despite its contribution to enhanced promoter activity of the Meishan GnRHR gene, the GATA-4 binding sites, like the NF-κB and SP1 elements, do not fully explain the increased activity of reporter constructs containing the Meishan compared to Control/Index promoters. Therefore, future studies in our laboratory will focus on identifying the remaining elements and corresponding binding factors that contribute to the enhanced activity of the Meishan GnRHR promoter in αT3-1 cells.
AP1, activator protein 1; CRE, cAMP responsive element; EMSA, electrophoretic mobility shift assay; FOXL2, forkhead box L2; FSH, follicle stimulating hormone; GnRH, gonadotropin-releasing hormone; GnRHR, gonadotropin-releasing hormone receptor; GnSE, GnRHR specific enhancer; GRAS, GnRH receptor activating sequence; GSE, gonadotrope specific element; JNK, c-Jun N-terminal kinase; LH, luteinizing hormone; LHX3, LIM homeobox 3; NF-Y, nuclear factor Y; NF-κB, nuclear factor-kB; OCT1, octamer transcription factor 1; PCR, polymerase chain reaction; Pitx-1, pituitary homeobox 1; PKA, protein kinase A; PKC, protein kinase C; RSV, Rous sarcoma virus; RXR, retinoid X receptor; SAP, SF1 adjacent protein; SF1, steroidogenic factor 1; SP, specificity protein; SURG-1, sequence underlying responsiveness to GnRH
The αT3-1 cells were a generous gift from Dr. Pamela Mellon (Salk Institute, La Jolla, CA). We would like to thank Dr. Colin Clay (Colorado State University, Ft Collins, CO) for the –m600pGL3 plasmid. Genomic DNA from the Meishan line of swine was kindly provided by Dr. Matthew B. Wheeler (University of Illinois, Urbana, IL). The authors would like to thank Rachel Friedrich, Chanho Lee and Jocelyn Wiarda for their assistance with completion of this study, as well as Dr. Janos Zempleni, Department of Nutrition and Health Sciences at UNL, for the use of his film developer. A contribution of the University of Nebraska Agricultural Research Division, Lincoln, Nebraska 68583. Journal Series No. 15039.
This project was supported by National Research Initiative Grant No. 2004-35205-14204 from the USDA Cooperative State Research, Education, and Extension Service Animal Genome program, the USDA National Institute of Food and Agriculture Hatch funds (NEB-13-159) and a UNL Layman Award to BRW.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional file.
EAM and JES contributed equally to this work, performed EMSAs, constructed reporter vectors and carried out transient transfections. EAM drafted the manuscript. RAC assisted in vector design, data interpretation and manuscript preparation. BRW conceived of the study, participated in its design and coordinated and assisted in the drafting of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Clayton RN, Catt KJ. Gonadotropin-releasing hormone receptors: Characterization, physiological regulation, and relationship to reproductive function. Endocr Rev. 1981;2:186–209.View ArticlePubMedGoogle Scholar
- Clarke IJ, Cummins JT, de Kretser DM. Pituitary gland function after disconnection from direct hypothalamic influences in the sheep. Neuroendocrinology. 1983;36:376–84.View ArticlePubMedGoogle Scholar
- Hamernik DL, Nett TM. Gonadotropin-releasing hormone increases the amount of messenger ribonucleic acid for gonadotropins in ovariectomized ewes after hypothalmic-pituitary disconnection. Endocrinology. 1988;122:959–66.View ArticlePubMedGoogle Scholar
- Gharib SD, Wierman ME, Shupnik MA, Chin WW. Molecular biology of the pituitary gonadotropins. Endocr Rev. 1990;11:177–90.View ArticlePubMedGoogle Scholar
- Sealfon SC, Millar RP. The gonadotropin-releasing hormone receptor: Structural determinants and regulatory control. Hum Reprod Update. 1995;1:216–30.View ArticlePubMedGoogle Scholar
- Wilkie PJ, Paszek AA, Beattie CW, Alexander LJ, Wheeler MB, Schook LB. A genomic scan of porcine reproductive traits reveals possible quantitative trait loci (QTLs) for number of corpora lutea. Mamm Genome. 1999;10:573–8.View ArticlePubMedGoogle Scholar
- Jiang Z, Gibson JP, Archibald AL, Haley CS. The porcine gonadotropin-releasing hormone receptor gene (GnRHR): Genomic organization, polymorphisms, and association with the number of corpora lutea. Genome. 2001;44:7–12.View ArticlePubMedGoogle Scholar
- Cheng CK, Leung PC. Molecular biology of gonadotropin-releasing hormone (GnRH)-I, GnRH-II, and their receptors in humans. Endocr Rev. 2005;26:283–306.View ArticlePubMedGoogle Scholar
- Hapgood JP, Sadie H, van Biljon W, Ronacher K. Regulation of expression of mammalian gonadotrophin-releasing hormone receptor genes. J Neuroendocrinol. 2005;17:619–38.View ArticlePubMedGoogle Scholar
- Clay CM, Nelson SE, DiGregorio GB, Campion CE, Wiedemann AL, Nett TM. Cell-specific expression of the mouse gonadotropin-releasing hormone (GnRH) receptor gene is conferred by elements residing within 500 bp of proximal 5' flanking region. Endocrine. 1995;3:615–22.View ArticlePubMedGoogle Scholar
- Duval DL, Nelson SE, Clay CM. A binding site for steroidogenic factor-1 is part of a complex enhancer that mediates expression of the murine gonadotropin-releasing hormone receptor gene. Biol Reprod. 1997;56:160–8.View ArticlePubMedGoogle Scholar
- Duval DL, Nelson SE, Clay CM. The tripartite basal enhancer of the gonadotropin-releasing hormone (GnRH) receptor gene promoter regulates cell-specific expression through a novel GnRH receptor activating sequence. Mol Endocrinol. 1997;11:1814–21.View ArticlePubMedGoogle Scholar
- McGillivray SM, Bailey JS, Ramezani R, Kirkwood BJ, Mellon PL. Mouse GnRH receptor gene expression is mediated by the LHX3 homeodomain protein. Endocrinology. 2005;146:2180–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Jeong KH, Chin WW, Kaiser UB. Essential role of the homeodomain for pituitary homeobox 1 activation of mouse gonadotropin-releasing hormone receptor gene expression through interactions with c-jun and DNA. Mol Cell Biol. 2004;24:6127–39.View ArticlePubMedPubMed CentralGoogle Scholar
- White BR, Duval DL, Mulvaney JM, Roberson MS, Clay CM. Homologous regulation of the gonadotropin-releasing hormone receptor gene is partially mediated by protien kinase c activation of an activator protein-1 element. Mol Endocrinol. 1999;13:566–77.View ArticlePubMedGoogle Scholar
- Duval DL, Ellsworth BS, Clay CM. Is gonadotrope expression of the gonadotropin releasing hormone receptor gene mediated by autocrine/paracrine stimulation of an activin response element? Endocrinology. 1999;140:1949–52.View ArticlePubMedGoogle Scholar
- Ellsworth BS, Burns AT, Escudero KW, Duval DL, Nelson SE, Clay CM. The gonadotropin releasing hormone (GnRH) receptor activating sequence (GRAS) is a composite regulatory element that interacts with multiple classes of transcription factors including smads, AP-1 and a forkhead DNA binding protein. Mol Cell Endocrinol. 2003;206:93–111.View ArticlePubMedGoogle Scholar
- Norwitz ER, Cardona GR, Jeong KH, Chin WW. Identification and characterization of the gonadotropin-releasing hormone response elements in the mouse gonadotropin-releasing hormone receptor gene. J Biol Chem. 1999;274:867–80.View ArticlePubMedGoogle Scholar
- Kam KY, Jeong KH, Norwitz ER, Jorgensen EM, Kaiser UB. Oct-1 and nuclear factor Y bind to the SURG-1 element to direct basal and gonadotropin-releasing hormone (GnRH)-stimulated mouse GnRH receptor gene transcription. Mol Endocrinol. 2005;19:148–62.View ArticlePubMedGoogle Scholar
- Lents CA, Farmerie TA, Cherrington BD, Clay CM. Multiple core homeodomain binding motifs differentially contribute to transcriptional activity of the murine gonadotropin-releasing hormone receptor gene promoter. Endocrine. 2009;35:356–64.View ArticlePubMedGoogle Scholar
- Resuehr D, Wildemann U, Sikes H, Olcese J. E-box regulation of gonadotropin-releasing hormone (GnRH) receptor expression in immortalized gonadotrope cells. Mol Cell Endocrinol. 2007;278:36–43.View ArticlePubMedGoogle Scholar
- Maya-Nunez G, Conn PM. Transcriptional regulation of the gonadotropin-releasing hormone receptor gene is mediated in part by a putative repressor element and by the cyclic adenosine 3',5'-monophosphate response element. Endocrinology. 1999;140:3452–8.PubMedGoogle Scholar
- Pincas H, Forraı̈ Z, Chauvin S, Laverrière J, Counis R. Multiple elements in the distal part of the 1.2 kb 5′-flanking region of the rat GnRH receptor gene regulate gonadotrope-specific expression conferred by proximal domain. Mol Cell Endocrinol. 1998;144:95–108.View ArticlePubMedGoogle Scholar
- Pincas H, Amoyel K, Counis R. Laverriere J-. Proximal cis-acting elements, including steroidogenic factor-1, mediate the efficiency of a distal enhancer in the promoter of the rat gonadotropin-releasing hormone receptor gene. Mol Endocrinol. 2001;15:319–37.View ArticlePubMedGoogle Scholar
- Granger A, Bleux C, Kottler ML, Rhodes SJ, Counis R, Laverriere JN. The LIM-homeodomain proteins isl-1 and Lhx3 act with steroidogenic factor 1 to enhance gonadotrope-specific activity of the gonadotropin-releasing hormone receptor gene promoter. Mol Endocrinol. 2006;20:2093–108.View ArticlePubMedGoogle Scholar
- Ngan ESW, Cheng KW, Leung PCK, Chow BKC. Steroidogenic factor-1 interacts with a gonadotrope-specific element within the first exon of the human gonadotropin-releasing hormone receptor gene to mediate gonadotrope-specific expression. Endocrinology. 1999;140:2452–62.PubMedGoogle Scholar
- Cheng KW, Ngan ESW, Kang SK, Chow BKC, Leung PCK. Transcriptional down-regulation of human gonadotropin-releasing hormone (GnRH) receptor gene by GnRH: Role of protein kinase C and activating protein 1. Endocrinology. 2000;141:3611–22.PubMedGoogle Scholar
- Cheng CK, Yeung CM, Hoo RLC, Chow BKC, Leung PCK. Oct-1 is involved in the transcriptional repression of the gonadotropin-releasing hormone receptor gene. Endocrinology. 2002;143:4693–701.View ArticlePubMedGoogle Scholar
- Duval DL, Farris AR, Quirk CC, Nett TM, Hamernik DL, Clay CM. Responsiveness of the ovine gonadotropin-releasing hormone receptor gene to estradiol and gonadotropin-releasing hormone is not detectable in vitro but is revealed in transgenic mice. Endocrinology. 2000;141:1001–10.PubMedGoogle Scholar
- Cederberg RA, Smith JE, McDonald EA, Lee C, Perkins AR, White BR. Activity of the porcine gonadotropin-releasing hormone receptor gene promoter is partially conferred by a distal gonadotrope specific element (GSE) within an upstream enhancing region, two proximal GSEs and a retinoid X receptor binding site. Reprod Biol Endocrinol. 2015;13:45.View ArticlePubMedPubMed CentralGoogle Scholar
- Johnson RK, Zimmerman DR, Kittok RJ. Selection for components of reproduction in swine. Livest Prod Sci. 1984;11:541–58.View ArticleGoogle Scholar
- Christenson RK, Vallet JL, Leymaster KA, Young LD. Uterine function in meishan pigs. J Reprod Fertil Suppl. 1993;48:279–89.PubMedGoogle Scholar
- White BR, McLaren DG, Dziuk PJ, Wheeler MB. Age at puberty, ovulation rate, uterine length, prenatal survival and litter size in chinese meishan and yorkshire females. Theriogenology. 1993;40:85–97.View ArticlePubMedGoogle Scholar
- Bass BE. Expression analysis of GnRH responsive anterior pituitary genes in lines of swine divergent for ovulation rate, M.S. Thesis. Department of Animal Science, University of Nebraska-Lincoln; 2005. Accessed from http://digitalcommons.unl.edu/animalscidiss/9/.
- Ho SN, Hunt HD, Horton RM, Pullen JK, Pease LR. Site-directed mutagenesis by overlap extension using the polymerase chain reaction. Gene. 1989;77:51–9.View ArticlePubMedGoogle Scholar
- Ellsworth BS, White BR, Burns AT, Cherrington BD, Otis AM, Clay CM. C-jun N-terminal kinase (JNK) activation of AP-1 underlies homologous regulation of the gonadotropin-releasing hormone receptor gene in αT3-1 cells. Endocrinology. 2003;144:839–49.View ArticlePubMedGoogle Scholar
- Cherrington BD, Farmerie TA, Lents CA, Cantlon JD, Roberson MS, Clay CM. Activin responsiveness of the murine gonadotropin-releasing hormone receptor is mediated by a composite enhancer containing spatially distinct regulatory elements. Mol Endocrinol. 2005;19:898–912.View ArticlePubMedGoogle Scholar
- Cheng CK, Chow BKC, Leung PCK. Functional mapping of a placenta-specific upstream promoter for human gonadotropin-releasing hormone receptor gene. Endocrinology. 2001;142:1506–16.PubMedGoogle Scholar
- Cheng CK, Yeung CM, Chow BKC, Leung PCK. Characterization of a new upstream GnRH receptor promoter in human ovarian granulosa-luteal cells. Mol Endocrinol. 2002;16:1552–64.View ArticlePubMedGoogle Scholar
- Yeung CM, An BS, Cheng CK, Chow BK, Leung PC. Expression and transcriptional regulation of the GnRH receptor gene in human neuronal cells. Mol Hum Reprod. 2005;11:837–42.View ArticlePubMedGoogle Scholar
- Barnhart KM, Mellon PL. The orphan nuclear receptor, steroidogenic factor-1, regulates the glycoprotein hormone alpha-subunit gene in pituitary gonadotropes. Mol Endocrinol. 1994;8:878–85.PubMedGoogle Scholar
- Keri RA, Nilson JH. A steroidogenic factor-1 binding site is required for activity of the luteinizing hormone ß subunit promoter in gonadotropes of transgenic mice. J Biol Chem. 1996;271:10782–5.View ArticlePubMedGoogle Scholar
- Brown P, McNeilly AS. Steroidogenic factor-1 (SF-1) and the regulation of expression of luteinising hormone and follicle stimulating hormone b-subunits in the sheep anterior pituitary in vivo. Int J Biochem Cell Biol. 1997;29:1513–24.View ArticlePubMedGoogle Scholar
- Hayden MS, Ghosh S. Signaling to NF-kappaB. Genes Dev. 2004;18:2195–224.View ArticlePubMedGoogle Scholar
- Grundker C, Schulz K, Gunthert AR, Emons G. Luteinizing hormone-releasing hormone induces nuclear factor kB-activation and inhibits apoptosis in ovarian cancer cells. J Clin Endocrinol Metab. 2000;85:3815–20.PubMedGoogle Scholar
- Hayden MS, Ghosh S. NF-kappaB, the first quarter-century: Remarkable progress and outstanding questions. Genes Dev. 2012;26:203–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Brauer VM, Bell JR, Desaulniers AT, Cederberg RA, White BR. Functional activity of the porcine Gnrhr2 gene promoter in testis-derived cells is partially conferred by nuclear factor-kappaB, specificity protein 1 and 3 (SP1/3) and overlapping early growth response 1/SP1/3 binding sites. Gene. 2016;587:137–46.View ArticlePubMedGoogle Scholar
- Naidich M, Shterntal B, Furman R, Pawson AJ, Jabbour HN, Morgan K, et al. Elucidation of mechanisms of the reciprocal cross talk between gonadotropin-releasing hormone and prostaglandin receptors. Endocrinology. 2010;151:2700–12.View ArticlePubMedGoogle Scholar
- Karalis KP, Venihaki M, van Vlerken LE, Chandras C. NF-kB participates in the corticotropin-releasing, hormone-induced regulation of the pituitary proopiomelanocortin gene. J Biol Chem. 2004;279:10837–40.View ArticlePubMedGoogle Scholar
- Parnet P, Pousset F, Laye S. NFkB activation in mouse pituitary: Comparison of response to interleukin-1ß and lipopolysaccharide. J Neuroendocrinol. 2003;15:304–14.View ArticlePubMedGoogle Scholar
- Philipsen S, Suske G. A tale of three fingers: The family of mammalian sp/XKLF transcription factors. Nucleic Acids Res. 1999;27:2991–3000.View ArticlePubMedPubMed CentralGoogle Scholar
- Suske G. The sp-family of transcription factors. Gene. 1999;238:291–300.View ArticlePubMedGoogle Scholar
- Suske G, Bruford E, Philipsen S. Mammalian SP/KLF transcription factors: Bring in the family. Genomics. 2005;85:551–6.View ArticlePubMedGoogle Scholar
- Lomberk G, Urrutia R. The family feud: Turning off Sp1 by Sp1-like KLF proteins. Biochem J. 2005;392:1–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Bieker JJ. Kruppel-like factors: Three fingers in many pies. J Biol Chem. 2001;276:34355–8.View ArticlePubMedGoogle Scholar
- van Vliet J, Crofts LA, Quinlan KG, Czolij R, Perkins AC, Crossley M. Human KLF17 is a new member of the sp/KLF family of transcription factors. Genomics. 2006;87:474–82.View ArticlePubMedGoogle Scholar
- Sahara S, Kawakami Y, Izpisua Belmonte JC, O'Leary DD. Sp8 exhibits reciprocal induction with Fgf8 but has an opposing effect on anterior-posterior cortical area patterning. Neural Dev. 2007;2:10.View ArticlePubMedPubMed CentralGoogle Scholar
- Harrison SM, Houzelstein D, Dunwoodie SL, Beddington RS. Sp5, a new member of the Sp1 family, is dynamically expressed during development and genetically interacts with brachyury. Dev Biol. 2000;227:358–72.View ArticlePubMedGoogle Scholar
- Kennedy MW, Chalamalasetty RB, Thomas S, Garriock RJ, Jailwala P, Yamaguchi TP. Sp5 and Sp8 recruit beta-catenin and Tcf1-Lef1 to select enhancers to activate wnt target gene transcription. Proc Natl Acad Sci U S A. 2016;113:3545–50.View ArticlePubMedGoogle Scholar
- Nakamura T, Yoshitomi Y, Sakai K, Patel V, Fukumoto S, Yamada Y. Epiprofin orchestrates epidermal keratinocyte proliferation and differentiation. J Cell Sci. 2014;127:5261–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakashima K, Zhou X, Kunkel G, Zhang Z, Deng JM, Behringer RR, et al. The novel zinc finger-containing transcription factor osterix is required for osteoblast differentiation and bone formation. Cell. 2002;108:17–29.View ArticlePubMedGoogle Scholar
- Kawakami Y, Esteban CR, Matsui T, Rodriguez-Leon J, Kato S, Izpisua Belmonte JC. Sp8 and Sp9, two closely related buttonhead-like transcription factors, regulate Fgf8 expression and limb outgrowth in vertebrate embryos. Development. 2004;131:4763–74.View ArticlePubMedGoogle Scholar
- Stein B, Baldwin Jr AS. Distinct mechanisms for regulation of the interleukin-8 gene involve synergism and cooperativity between C/EBP and NF-kappa B. Mol Cell Biol. 1993;13:7191–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Perkins ND, Edwards NL, Duckett CS, Agranoff AB, Schmid RM, Nabel GJ. A cooperative interaction between NF-kappa B and Sp1 is required for HIV-1 enhancer activation. EMBO J. 1993;12:3551–8.PubMedPubMed CentralGoogle Scholar
- Hong CY, Park KH, Kim J, Im SY, Lee JW, Choi H, et al. Expression of MIS in the testis is downregulated by tumor necrosis factor alpha through the negative regulation of SF-1 transactivation by NF-kB. Mol Cell Biol. 2003;23:6000–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Wissink S, van der Burg B, Katzenellenbogen BS, van der Saag PT. Synergistic activation of the serotonin-1A receptor by nuclear factor-kappa B and estrogen. Mol Endocrinol. 2001;15:543–52.PubMedGoogle Scholar
- Perkins ND, Agranoff AB, Pascal E, Nabel GJ. An interaction between the DNA-binding domains of RelA(p65) and Sp1 mediates human immunodeficiency virus gene activation. Mol Cell Biol. 1994;14:6570–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Steger DJ, Hecht JH, Mellon PL. GATA-binding proteins regulate the human gonadotropin alpha-subunit gene in the placenta and pituitary gland. Mol Cell Biol. 1994;14:5592–602.View ArticlePubMedPubMed CentralGoogle Scholar
- LaVoie HA. The role of GATA in mammalian reproduction. Exp Biol Med (Maywood). 2003;228:1282–90.Google Scholar
- Lawson MA, Whyte DB, Mellon PL. GATA factors are essential for activity of the neuron-specific enhancer of the gonadotropin-releasing hormone gene. Mol Cell Biol. 1996;16:3596–605.View ArticlePubMedPubMed CentralGoogle Scholar
- Watanabe K, Clarke TR, Lane AH, Wang X, Donahoe PK. Endogenous expression of mullerian inhibiting substance in early postnatal rat sertoli cells requires multiple steroidogenic factor-1 and GATA-4-binding sites. Proc Natl Acad Sci U S A. 2000;97:1624–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Tremblay JJ, Viger RS. Transcription factor GATA-4 enhances mullerian inhibiting substance gene transcription through a direct interaction with the nuclear receptor SF-1. Mol Endocrinol. 1999;13:1388–401.PubMedGoogle Scholar
- Fluck CE, Miller WL. GATA-4 and GATA-6 modulate tissue-specific transcription of the human gene for P450c17 by direct interaction with Sp1. Mol Endocrinol. 2004;18:1144–57.View ArticlePubMedGoogle Scholar