| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Departments of Molecular Cell Biology (M.J.M.S., I.H.C.v.L., D.C.W.A.v.W., A.H.M., H.P.S.) and Integrative Zoology (M.K.R.), Institute of Biology, and Division of Medical Pharmacology (D.C., O.C.M.), Leiden/Amsterdam Center for Drug Research/Leiden University Medical Center, Leiden University, 2300 RA Leiden, The Netherlands
Address all correspondence and requests for reprints to: Marcel J. M. Schaaf, Department of Molecular Cell Biology, Institute of Biology, P.O. Box 9505, 2300 RA Leiden, The Netherlands. E-mail: m.j.m.schaaf{at}biology.leidenuniv.nl.
| Abstract |
|---|
|
|
|---|
and GRβ, which are identical between amino acids 1–727 and then diverge. Whereas GR
(the canonical GR) acts as a ligand-activated transcription factor, GRβ does not bind traditional glucocorticoid agonists, lacks GR
s transactivational activity, and acts as a dominant-negative inhibitor of GR
. It has been suggested that this receptor isoform is involved in the induction of glucocorticoid resistance in asthma patients. Unfortunately, a GR β-isoform has been detected in only humans, and therefore, an animal model for studies on this isoform is lacking. In the present study, we demonstrate that in zebrafish a GR isoform exists that diverges from the canonical zebrafish GR at the same position as human GRβ from human GR
. The zebrafish GR β-isoform acts as a dominant-negative inhibitor in reporter assays, and the extent of inhibition and the effective GR
/GRβ ratio is similar to studies performed with the human GR isoforms. In addition, the subcellular localization of zebrafish GRβ is similar to its human equivalent. Finally, expression levels of GR
and GRβ were determined in adult zebrafish tissues and at several developmental stages. Both receptor isoforms were detected throughout the body, and GRβ mRNA levels were relatively low compared with GR
mRNA levels, as in humans. Thus, for the first time, a GR β-isoform has been identified in a nonhuman animal species, shedding new light on the relevance of this GR splice variant and providing a versatile animal model for studies on the GR system. | Introduction |
|---|
|
|
|---|
B and activation protein-1, thereby inhibiting their activity (2).
Two splice variants of the human GR (hGR) occur: hGR
and hGRβ (3, 4). The two receptor isoforms are identical between amino acid 1–727 and then diverge. hGR
contains 50 specific C-terminal amino acid residues that form two helical structures that play an important role in ligand binding, whereas the C terminus of hGRβ contains 15 specific, nonhomologous amino acids, which are largely disordered (5). hGR
is the canonical GR, whereas hGRβ lacks GR
s transactivational activity on GRE-containing promoters and does not bind traditional glucocorticoid agonists. However, binding to the glucocorticoid antagonist RU486 has recently been reported (6). In several studies, hGRβ has been shown to act as a dominant-negative inhibitor of hGR
s transactivational properties in transiently transfected cell cultures (7, 8, 9). The mechanism of this inhibition is still unclear, but competition between hGR
and hGRβ for transcriptional coactivator proteins (10) and the formation of inactive hGR
-hGRβ heterodimers (5, 9, 11, 12) have been suggested. Specific RNA silencing of hGRβ in human monocytes and macrophages results in enhanced transactivation and transrepression induced by hGR
(13, 14), whereas overexpression of human hGRβ in mouse hybridoma cells results in decreased responsiveness to glucocorticoids in these cells (11).
In several studies, a correlation has been found between increased expression of GRβ and the occurrence of immune-related diseases. In ulcerative colitis (15), leukemia (16), and severe asthma patients (17), increased expression levels of this receptor isoform were found in various immune cells. Moreover, increased hGRβ expression correlates with resistance to glucocorticoid treatment in patients suffering from these diseases (13, 15, 16, 18, 19), suggesting that hGRβ plays an important role in glucocorticoid resistance. Interestingly, the glucocorticoid responsiveness of bronchoalveolar lavage cells from glucocorticoid-resistant asthma patients is enhanced by decreasing the hGRβ expression level using specific RNA silencing (13).
Until now, a GR β-isoform has been found only in humans. Rodents do not express a GRβ due to the absence of the alternative splice site (20). Therefore, an animal model for studies on GRβ is lacking, and studies on GRβ have been limited to human cells and tissues. In the present study, we have discovered the occurrence of a GR β-isoform in zebrafish, which is strikingly similar in structure, function, and expression level to its human equivalent. This provides a versatile animal model for studies on this pharmacologically important protein and sheds new light on the physiological relevance of hGRβ.
| Materials and Methods |
|---|
|
|
|---|
In silico prediction of the zebrafish GR gene
To identify the zebrafish GR gene(s), BLAST searches of the zebrafish genomic sequence were carried out using the TBLASTN tool at the Sanger Institute Ensembl server (www.ensembl.org). The zebrafish genome release Zv7 was used for these searches using the full-length protein sequences of Burtons mouthbrooder (21) and rainbow trout (22) GR1 and GR2 as queries. A single zebrafish GR genomic sequence was identified. Analysis of the syntenic regions of the fish GR genes in tetraodon, fugu, medaka, stickleback, and zebrafish was done using the available annotated regions of the respective genomes on the Ensembl website. The National Center for Biotechnology Information (NCBI) website (www.ncbi.nlm.nih.gov) was used to search for expressed sequence tag and cDNA sequences derived from zebrafish GR transcripts. Text searches and BLAST searches (using Burtons mouthbrooder and rainbow trout GR1 and GR2 sequences as queries) were performed.
RT-PCR analysis of zGR
and zGRβ cDNA
Total RNA was isolated from whole adult zebrafish, livers dissected from adult fish, larvae at 72 h after fertilization, or cultured zebrafish fibroblast (ZF4) cells using Trizol reagent (Invitrogen, Carlsbad, CA) according to the manufacturers instructions. Subsequently, RT-PCR analysis was performed using the Superscript III One-Step RT-PCR system with Platinum Taq (Invitrogen). Forward primers containing sequences from exon 1 (AATAACTCTGTACTGCCAAT), exon 2 (GCAAAATGGATCAAGGAGGA), and exon 7 were used (AGAGAGGATGAAGTTGCCCT), and reverse primers were based on sequences from exon 2 (CGCCTTTAATCATGGGAGAA), exon 8β (CCAAGCGGAATCACTATGACG), and exon 9 (CTGCTGTTGGGAGGAGATTC). PCR fragments were cloned into pCRII-TOPO vectors (Invitrogen) and sequenced using the sequencing service of BaseClear (Leiden, The Netherlands).
Northern blotting
Probes hybridizing in exon 8β and 9 were used (location of sequences is indicated in Fig. 1A
by solid bars). RT-PCR was performed on RNA isolated from adult zebrafish, to obtain a 519-bp amplicon from exon 8β (primers used were CCATGTCTCTCTCTGCTCCAAAGCC and TCCAATATTGTGCAGCCCTAATCCA) and a 567-bp amplicon from exon 9 (primers used were TTGGTGGGTGGACTCCTGAACTTC and GGACTGTAAGTGCTGCATCTGCTCC). Total RNA derived from whole zebrafish was reverse-transcribed using oligo-dT (Life Technologies, Inc., Rockville, MD) and Super RT (HT Biotechnology, Cambridge, UK), and PCR amplification was performed using Phusion DNA polymerase (Finnzymes, Espoo, Finland) according to the manufacturers instructions. The amplified products were cloned into pCRII-TOPO vectors and excised by restriction enzyme digestion. The DNA fragments were end-labeled with [32P]dCTP and the Klenow fragment of DNA polymerase I.
|
Luciferase reporter assays
Expression vectors for zGR
and zGRβ were constructed. Coding sequences were amplified by RT-PCR (as described under Northern blotting) using a forward primer with a linker containing a BamHI restriction site and reverse primers containing XhoI sites (forward primer, CCGGATCCGCAAAATGGATCAAGGAGGA; reverse primer for zGR
, CCCTCGAGCTGCTGTTGGGAGGAGATTC; reverse primer for zGRβ, CCCTCGAGCCAAGCGGAATCACTATGACG). Fragments were cloned into the expression vector pCS2+ (see http://sitemaker.umich.edu/dlturner.vectors/home) using BamHI and XhoI restriction sites. The resulting vectors were named pCS2+zGR
and pCS2+zGRβ. COS-1 cells (African green monkey kidney cells lacking an endogenous GR, and used in most published human GRβ studies) (5, 7, 9) were used for the luciferase assays. The day before transfection, the cells were seeded into 24-well plates at the appropriate density. Plasmids pCS2+zGR
and/or PCS2+zGRβ were transfected together with the reporter constructs pCMV-renilla (Promega, Madison, WI) and pMMTV-luciferase using Superfect reagent (QIAGEN, Valencia, CA). Different amounts of the empty vector pCS2+ were cotransfected to use equal amounts of total DNA in the individual transfections. Cells were incubated at 37 C for 24 h, a glucocorticoid ligand (dexamethasone or cortisol) was added, and cells were incubated at 37 C for 24 h. Cells were lysed and the firefly and renilla luciferase activity was measured using the Promega Dual Label Reporter assay. Bioluminescence was detected in a Berthold luminometer (LUMAT LB 9507; Berthold, Bad Wildbad, Germany). Averages were taken from at least three replicates, and the average firefly luciferase activity was divided by the average renilla luciferase activity. Data shown are means ± SE of at least three individual experiments. Data from zGR
and zGRβ coexpression studies (Fig. 3A
) were analyzed by ANOVA using the software SPSS 7.5.
|
and YFP-zGRβ were constructed to obtain YFP-GR fusion proteins with an identical linker sequence as previously made YFP-hGR constructs (23). Amplification of zGR
and zGRβ cDNA sequence was performed using primers containing XhoI and BamHI restriction sites and pCS2+zGR
and pCS2+zGRβ as templates (primer sequences: CC GCT CGA GCA GGA GCA GGC GCC GGG GCG GGT GCT ATGGATCAAGGAGGACTGGAG as forward primer, and as reverse primer, CCGGGATCCCTGCTGTTGGGAGGAGATTC for zGR
and CCGGGATCCCCAAGCGGAATCACTATGACG for zGRβ. Using XhoI and BamHI restriction sites, the amplicons were ligated into the vector pEYFP-C1 (Clontech, Mountain View, CA). The resulting plasmids were transfected into COS-1 cells (plated in four-well chambered coverglasses from Nunc, Roskilde, Denmark) using TransIt Cos (Mirus), and cells were imaged 24 h later. Imaging was performed using a Leica (Leica Microsystems, Wetzlar, Germany) SP1 confocal laser scanning microscope. Excitation was done with an argon laser at 514 nm. Images were captured using a x63 water immersion objective, a 1024 x 1024 pixel resolution was used, and images are averages of six scans. Quantitation of the fluorescence was done using LCS Lite software (Leica).
Western blotting
Lysates of COS-1 cells transfected with yellow fluorescent protein (YFP)-zGR
or YFP-zGRβ were analyzed by Western blotting. One day after transfection, cells were scraped off in ice-cold PBS, spun down, and resuspended in a 0.25 M Tris/5 mM EDTA buffer (pH 8), lysed by tip sonication, and incubated at 95 C for 10 min. Protein samples were resolved by electrophoresis through polyacrylamide gels and electrophoretically transferred to nitrocellulose membranes. Membranes were blocked at room temperature for 30 min in PBS with Tween 20 [PBST; 10 mM Tris-HCl (pH 7.4), 150 mM NaCl, and 0.1% Tween 20] containing 5% nonfat dry milk. Subsequently, the membrane was incubated overnight at 4 C with a monoclonal anti-green fluorescent protein antibody (JL-8; Clontech) at a dilution of 1:1000 in PBST containing 5% nonfat dry milk. After three 15-min washes in PBST containing 5% nonfat dry milk, the membrane was incubated with a horseradish peroxidase (HRP)-labeled antimouse secondary antibody (Cell Signaling Technology, Danvers, MA) (1:10,000 dilution in PBST containing 5% nonfat dry milk) for 1.5 h at room temperature. Finally, after three 15-min washes in PBST, immunoreactivity was visualized by enhanced chemiluminescence (GE Healthcare Life Sciences, Piscataway, NJ) according to the manufacturers instructions.
Quantitative RT-PCR analysis
Total RNA was extracted using Trizol reagent. A DNase treatment was performed, using the DNA-free kit (Ambion, Austin, TX), and cDNA was synthesized as described under Northern blotting. PCR amplification of the cDNA was performed using DNA Engine Opticon2 (Bio-Rad, Hercules, CA) and qPCR Core kit for SYBR Green I (Bio-Rad). Primers were designed to amplify a region in exon 8β (GATGAACTACGAATGTCTTA and GCAACAGACAGCCAGACAGCTCACT) or exon 9 (AACTGGCACGGTTCTATCAGCTCA and TTCTGGTGAAAGCAGCGG). For normalization, in separate reactions, primers were used to amplify β-actin mRNA (CGAGCAGGAGATGGAACC and CAACGGAAACGCTCATTGC).
To study whether genomic sequences were amplified, a control sample was used to which no reverse transcriptase was added (non-RT control). The thermocycling protocol consisted of 2 min at 50 C, 10 min at 95 C, 40 cycles of 1 sec at 95 C, 30 sec at 60 C, and 30 sec at 72 C and finished with a melting curve ranging from 65–95 C to allow discrimination of specific products. Standard curves were made using serial dilutions of known amounts of PCS2+zGR
or PCS2+zGRβ and used for quantitation of specific cDNAs in the samples. A standard curve for β-actin was made using serial dilutions of cDNA sample, so the relative amount of β-actin in each sample could be determined. Finally, the zGR
and zGRβ data were normalized for β-actin. Data shown are means ± SE of at least three individual experiments.
In situ hybridization
For both receptor isoforms, both a sense (control) and an antisense probe was constructed. Plasmids previously described under Northern blotting containing exon 9 and 8β sequences were used and linearized by enzymatic digestion, purified, and stored at –20 C until use. The RNA probe was synthesized from 1 mg of the DNA plus 2 µl 0.1 M dithiothreitol (Invitrogen), 0.5 µl RNase out (Invitrogen), 2 µl dNTPs digoxigenin labeling mix (Roche, Mannheim, Germany), 2 µl 10x transcription buffer (Ambion), and 1.5 µl of either SP6 or T7 RNA polymerase (Ambion). After a 2-h incubation period at 37 C, DNase I (Roche) was added and incubated for 15 min at 37 C. The reaction was stopped with 1 µl 0.5 M EDTA. The RNA was precipitated, air dried, and taken up in water.
Whole 24-h embryos (approximately 10 embryos were used per condition) were fixed in 4% paraformaldehyde (PFA) overnight before being washed with PBST (containing 1% Tween 20), transferred to methanol in graded steps, and then stored at –20 C until use. Before use, embryos were transferred back into PBST in graded steps. The embryos were then washed again with PBST, and all embryos were refixed with 4% PFA. After washing off the PFA with PBST, the embryos were prehybridized in hybridization buffer for 5 h before the RNA probe was added to a concentration of 1 ng/µl. The embryos were left to hybridize overnight at 70 C. Embryos were transferred in series to 2x standard saline citrate with Tween 20, then to 0.2x standard saline citrate with Tween 20, and finally to PBST. Antibody buffer was then added, and after 2 h, anti-digoxigenin alkaline phosphatase-labeled antibody (Roche) was added to the embryos in a 1:3000 dilution. The embryos were left in antibody for 2 h. The embryos were washed with PBST overnight. After additional washes with PBST, the embryos were transferred to AP buffer and then stained with purple Boehringer Mannheim AP substrate staining solution (Roche). Once staining appeared optimal, stop solution (1 mM EDTA in PBS, pH 5.5) was added before the embryos were refixed in 4% PFA for 20 min and washed with PBST. Each pair of conditions (antisense and sense for each receptor isoform) was stained for the same amount of time to validate the control in the experiment.
Images were captured using a Leica microscope (MZ FLIII) with connected camera (DC500; Leica Microsystems) and computer software (FireCam IM50).
| Results |
|---|
|
|
|---|
peptide sequence reveals that 47.1% of the amino acids are identical, with high percentages in the DNA-binding domain (96.7%, amino acids 394–454) and the ligand-binding domain (73.3%, amino acids 496–746; supplemental Fig. 1A, published as supplemental data on The Endocrine Societys Journals Online web site at http://endo.endojournals.org). In many teleostean fish species, two GR genes have been found, and the resulting receptor isoforms have been named GR1 and GR2. Our analysis of the phylogenetic tree of fish GRs (supplemental Fig. 2) confirms previously published data demonstrating that this zebrafish GR clusters within the GR2 clade of fish GRs (24, 25). We strongly suggest that the zebrafish does not contain a second GR gene (representing the GR1 ortholog), based on three lines of evidence. First, BLAST searches using GR1 peptide sequences from rainbow trout and Burtons mouthbrooder as queries returned sequences of genes encoding all zebrafish steroid receptors and many nuclear receptor but not a second GR gene. For example, when the rainbow trout GR1 sequence was used as a query, the genes for the three zebrafish estrogen receptors that have been identified before (26, 27) and the five previously described estrogen receptor-related receptor genes (28) were found, as well as genes encoding one zebrafish GR, one mineralocorticoid receptor, one progestin receptor, and one androgen receptor (supplemental Table 1). Second, BLAST searching in GenBank for cDNA and expressed sequence tag sequences representing possible transcripts from a zebrafish GR gene revealed fourteen sequences (supplemental Table 2). Sequence alignment showed that all 14 sequences turned out to be transcripts from the identified GR gene. Third, analysis of the syntenic region of GR1 and GR2 genes in fugu, tetraodon, medaka, and stickleback demonstrates that the GR2 gene and its syntenic region is intact in the zebrafish genome but that the GR1 syntenic region has undergone major rearrangement, apparently resulting in the loss of the GR1 gene (supplemental Fig. 3).
The zebrafish GR
- and β-isoform
A partial mRNA prediction had been deposited in GenBank (acc. no. XM_696233), which could be an alternative transcription product of this gene as a result of alternative splicing between exons 8 and 9. Northern blot analysis on RNA from embryonic zebrafish was performed using probes hybridizing against sequences downstream of this putative alternative splice site. Both probes stained a single band, but the size and intensity of the bands were different, indicating that alternative splicing indeed occurs (Fig. 1B
). RT-PCR analysis confirmed the occurrence and identity of the alternative splice product in adult zebrafish, in larvae, and in cultured fibroblast cells (Fig. 1C
). The full-length coding sequence of this mRNA was determined, confirming that this mRNA was the result of alternative splicing between exons 8 and 9 (Fig. 1A
). This sequence was submitted to GenBank under accession no. EF436285. The encoded protein diverges from zGR at amino acid 697 and contains a C-terminal end of 40 specific amino acids (Fig. 1D
). Protein alignments show that divergence from the zGR sequence occurs at the exact same position as divergence between human GR
and GRβ (Fig. 1D
). We therefore named this protein zGRβ.
Both in the human and in the zebrafish GR protein sequence, the point of divergence between the GR
- and β-isoform is located in the transition region between the 10th and 11th helix of the ligand-binding domain (5, 29). Both GR β-isoforms therefore lack helix 12, which contains a transactivation domain (AF-2). Analysis of the predicted secondary protein structure reveals that both hGRβ and zGRβ contain a shortened version of helix 11 and a C-terminal tail that is largely disordered (supplemental Fig. 4). However, the 40 zGRβ-specific amino acids display a very low level of homology with the 15 GRβ-specific amino acids of hGRβ (supplemental Fig. 1B), although the two-amino-acid motif (lysine-proline) that has been demonstrated to be responsible for the dominant-negative effect of hGRβ (5) is present in zGRβ and is predicted to be located at the same position in hGRβ and zGRβ, immediately behind the shortened helix 11 (supplemental Fig. 4).
Comparison between the human and zebrafish GR gene reveals that the gene organization and splicing events leading to GR
and GRβ mRNA are different (Fig. 1A
). The sequence encoding the β-specific amino acids is located in exon 9 in the human gene, whereas in the zebrafish gene, this sequence is found in exon 8. Therefore, hGR
and hGRβ mRNA are produced through alternative usage of a splice acceptor site in exon 9, whereas in zebrafish, alternative use of a splice donor site in exon 8 is the underlying mechanism. In zebrafish, using the most 5' splice donor site results in a shorter version of exon 8 and an open reading frame that includes exon 9, which encodes zGR
-specific amino acids. Using the most 3' site results in an extended version of exon 8 in zGRβ mRNA, introducing a stop codon in exon 8.
The dominant-negative activity of zGRβ
Transcriptional properties of zGR
and zGRβ were studied using luciferase assays in COS-1 cells (Fig. 2
), which are commonly used for this type of assay because they lack an endogenous GR. First, expression vectors for zGR
and zGRβ were transfected into these cells together with a vector containing a luciferase gene driven by the MMTV promoter, which contains several GREs. Luciferase activity induced by zGR
was increased 15-fold at 24 h after addition of dexamethasone (100 nM), a synthetic GR ligand. Addition of cortisol (100 nM), the main endogenous glucocorticoid hormone in fish, resulted in a 10-fold increase in luciferase activity (Fig. 2A
). In a separate experiment, dose-response curves were made for these ligands, and EC50 values were determined (0.37 and 10.1 nM for dexamethasone and cortisol respectively, Fig. 2B
) that are in line with previously found EC50 values using the human GR
in comparable assays [0.3 nM for dexamethasone (30) and 8 nM for corticosterone (31)]. No induction in transcriptional activity of zGRβ was observed using these two ligands, confirming our expectation that zGRβ has no transactivational activity.
|
and zGRβ expression vectors were cotransfected at different ratios, and luciferase activity was measured 24 h after addition of dexamethasone. Coexpression of zGRβ decreased the zGR
-induced transcription, depending on the amount of zGRβ present (Fig. 3A
and zGRβ was determined after transfection of different expression vector amounts (induced by 100 nM dexamethasone, Fig. 3B
between 30 and 300 ng DNA, whereas for zGRβ, luciferase activity still increases with higher DNA amounts in this range. These data indicate that the dominant-negative effect of zGRβ observed in this assay (Fig. 3A
Subcellular localization of zGR
and zGRβ
To study the subcellular localization of zGR
and zGRβ, expression plasmids encoding zGR
and zGRβ tagged with a YFP at the N terminus were constructed and transfected into COS-1 cells. Images of transfected cells were made using confocal microscopy, and representative images are shown in Fig. 4A
. Quantitation of the nuclear localization is shown in supplemental Table 3. YFP-zGR
was localized in both the nucleus and the cytoplasm at comparable levels. After administration of dexamethasone, complete translocation of the receptor to the nucleus was observed. YFP-zGRβ showed a slightly more nuclear localization than YFP-zGR
in the absence of ligand. This localization was not changed after addition of dexamethasone.
|
and hGRβ constructs in our study (supplemental Table 3). The linker sequence between YFP and GR in these constructs was identical to the respective sequence in the YFP-zGR constructs. Subcellular localization of the YFP-hGR constructs were in line with previously published data (6). YFP-hGRβ showed a more nuclear localization than YFP-zGR
, and this localization was not changed by dexamethasone administration. Thus, both the human and zebrafish YFP-tagged GRs show that the β-isoform displays a slightly more nuclear localization than the
-isoform (in the absence of ligand). Interestingly, both the human
-isoform and the β-isoform showed a more cytoplasmic localization in the absence of ligand than their zebrafish equivalents.
To control whether the fusion protein was correctly expressed, a Western blot was performed on cell lysates from COS-1 cells transfected with YFP-zGR
and YFP-zGRβ (Fig. 4B
). Using an antibody against YFP, a single band of approximately 120 bp was detected, which is the expected size of the fusion proteins. The expression level of YFP-zGR
was slightly higher than the level of YFP-zGRβ, probably reflecting a difference in stability between the two isoforms as observed between the human isoforms (32).
Expression of zGR
and zGRβ mRNA in vivo
To study the relative expression levels of zGR
and zGRβ mRNA, quantitative RT-PCR analysis was performed on RNA from several tissues of adult zebrafish (Fig. 5A
). In all tissues studied, zGR
and zGRβ mRNA appear to be coexpressed, with zGR
mRNA levels being significantly higher than zGRβ mRNA levels. Ratios of GR
/β mRNA levels varied between 13.3 (in muscle) and 60.6 (in spleen).
|
and zGRβ mRNA levels were measured at different stages during the first 72 h of zebrafish development (Fig. 5B
mRNA levels were significantly higher than zGRβ mRNA levels. Ratios of GR
/β mRNA ranged from 20.7 (at 72 h after fertilization) to 146 (at tail bud stage, i.e. 10 h after fertilization).
To study the localization of zGR
and zGRβ mRNA expression, whole mount in situ hybridizations were performed, using zGR
- and zGRβ-specific probes on 24-h embryos. Expression of zGR
mRNA was observed throughout the body (Fig. 6A
), confirming an in situ hybridization study previously published on the ZFIN website (33) using a probe for gene zgc:113038, which turns out to be the zGR gene. The zGRβ mRNA expression pattern looks similar, indicating that the two isoforms are colocalized (Fig. 6B
). Hybridization with sense probes did not result in any staining of the embryo (data not shown).
|
Finally, we studied whether the GR gene organization that is present in the human GR gene, i.e. the presence of an exon (9β) encoding GRβ-specific amino acids downstream of exon 9, occurs in other species as well. The exon 9 sequence of the GR gene of a wide variety of animals was studied. First we aligned exon 9 of the human GR gene with its equivalent in several nonmammalian vertebrates and nonplacental mammals. A low level of similarity was found in the region around the start of exon 9β, and the equivalent of exon 9β could not be found in these species (data not shown). Second, the human GR exon 9 sequence was aligned with its counterparts in other placental mammalian species (supplemental Fig. 6). In the rodent species studied (mouse, rat, and guinea pig), the donor splice site enabling GRβ expression is absent, as demonstrated before (20). In addition, this splice site appears to be absent in hedgehog, dog, and cat. Presence of a splice site was demonstrated in four primate species and in rabbit, horse, and cow. In both elephant and shrew, a splice site with a very low prediction score was observed. Translation products from these putative exon 9β sequences were aligned, demonstrating that this sequence is 100% identical among human, chimpanzee, and macaque and 76% identical between human and rabbit. The bushbaby, horse, cow, elephant, and shrew sequences show significantly lower similarity to the human sequence. These data show that the presence of a GR β-isoform as it occurs in humans is restricted to placental mammals and that within this group of animals, the genomic characteristics enabling GRβ expression (a specific donor splice site in exon 9 and/or GRβ-specific amino acid encoding sequences) are poorly conserved.
| Discussion |
|---|
|
|
|---|
. This splice variant was therefore suggested to be the zebrafish equivalent of hGRβ. We demonstrated that this zGRβ has a function similar to its human equivalent. In reporter assays in COS cells, it was demonstrated that zGRβ acts as a dominant-negative inhibitor of zGR
s transactivational activity. Both the extent of the inhibition (i.e. approximately 50%) and the effective zGRβ/zGR
ratio (i.e. 10) are remarkably similar to results of comparable experiments in which the transcriptional activity of hGR
and hGRβ was studied in related cell lines (5, 7, 9, 10).
Using YFP-tagged receptors, the subcellular localization of zGR
and zGRβ was studied in COS cells. The localization patterns appeared to be similar to those observed for YFP-tagged hGR
and hGRβ in our study and in previously published studies on the subcellular localization of (YFP-tagged) hGRβ (5, 6, 13, 14). Both the zebrafish and the human β-isoform show a more nuclear localization than their respective unliganded
-isoform, but less nuclear than the liganded GR
. Both human GR
and GRβ were shown to be more nuclearly localized than their zebrafish equivalents. This may be due to the fact that the zebrafish constructs were not expressed in an appropriate (zebrafish) cell line.
The relative expression levels of zGR
and zGRβ mRNA were measured, and relatively low zGRβ expression levels were observed both in tissues from adult fish and in embryos at different developmental stages. The observed zGR
/zGRβ ratios were similar to those observed in human cells and tissues (8, 15, 34, 35), in which significantly lower levels of hGRβ mRNA are observed compared with hGR
. These relatively low hGRβ expression levels raise questions about the physiological relevance of this receptor isoform. However, the low level of hGRβ mRNA expression is not reflected at the protein level (12, 36), probably due to the greater stability of the hGRβ protein (32). In addition, the expression of hGRβ appears to be inducible by proinflammatory cytokines, including IL-2 and IL-4 (18) and TNF
(32), and expression within tissues differs between cell types. Several specific cell types, mainly epithelial cells, show higher expression levels of hGRβ protein (37).
Although the human and zebrafish GR β-isoforms show remarkable similarities in structure, function, and expression level, it should be noted that there are differences in the splicing events resulting in human and zebrafish GRβ mRNA and that homology between the human and zebrafish GRβ-specific C-terminal peptide sequences is lacking. We therefore suggest that this functional outcome was achieved independently in humans and zebrafish by strong positive selection. This convergent evolution sheds new light on the physiological relevance of GRβ in the glucocorticoid system.
The presence of a GR β-isoform in other species
Our studies of the GR genes of a wide variety of animal species show that the presence of a GR β-isoform is restricted to a relatively small group of animals. The gene organization enabling GRβ expression as shown in zebrafish could not be detected in tetrapods or in the other teleost fish species tetraodon, fugu, medaka, and stickleback. Because these four fishes belong to the superorder of Acantopterygii, and the zebrafish is a member of the Ostariophysi superorder, the acquisition of GRβ appears to have occurred after the Ostariophysii superorder branched off the lineage that led to Acantopterygii (110–160 million years ago) and is therefore probably restricted to species in this superorder of teleost fish. In addition, the gene organization enabling GRβ expression as shown in humans appeared to be restricted to placental mammals, and even within this group, it was shown to be poorly conserved. Because only a limited group of animals has acquired and conserved a GR gene organization that results in the presence of a GR β-isoform, this adaptation of the GR system appears not to be essential for glucocorticoid signaling. However, the convergent evolution of GRβ in certain mammalian and fish species suggests a significant role of this receptor isoform in the regulation of glucocorticoid responsiveness in these species.
The presence of a single GR gene in zebrafish
Surprisingly, only one GR gene was identified in the zebrafish genome. This is remarkable, because in many teleostean fish species, such as rainbow trout (21), Burtons mouthbrooder (22), tetraodon, fugu and common carp (25), and sea bass (38, 39), two GR genes have been found (24, 25). These duplicate GRs appear to be a result of the genome duplication that occurred in fish between 300 and 450 million years ago, shortly after the divergence of the tetrapods from the fish lineages (40). In a few other fish species, such as Japanese flounder and brown trout (25), evidence for only one GR gene has been found thus far, but it is yet unclear whether these fishes contain a second GR gene because most of these species are poorly studied. From our analysis, it appears that the zebrafish has lost one GR gene copy during its evolution, probably as a result of rearrangements in chromosome 21. The loss of this gene appears to be limited to a small subset of fish species of the family of Cyprinidae, possibly to the Danioninae subfamily to which the zebrafish (genus Danio) belongs, because the presence of two isoforms has been reported in common carp (25), which is a member of the Cyprinidae family as well (Cyprininae subfamily, genus Cyprinus). It can be hypothesized that the zebrafish has acquired an alternative mechanism of regulation of glucocorticoid responsiveness by means of a dominant-negative GR splice variant to compensate for the loss of a second GR gene.
As a result of the loss of this GR gene, a GR splice variant, the equivalent of rainbow trout GR1b (41), which has orthologs in Burtons mouthbrooder (21), fugu, and tetraodon (25), is lacking in zebrafish. This splice variant, which contains an additional nine amino acids between the two zinc fingers of the DNA-binding domain, is a product of the GR genes of the GR1 clade, which is the GR gene that is absent in zebrafish.
Conclusions
In conclusion, during evolution, zebrafish have acquired a β-isoform of the GR, which displays similarities in structure, function, subcellular localization, and expression level to its human equivalent. However, differences in the splicing events resulting in human and zebrafish GRβ mRNA exist, and the GRβ-specific C-terminal peptide sequence displays differences. Until now, a GR β-isoform had not been observed in any species apart from humans, and its presence in humans and zebrafish is most likely a result of convergent evolution. This and the fact that only a single GR gene appears to be present in zebrafish makes it a unique animal model system for studies on the GR system, especially given the versatility and suitability of the zebrafish in genetic, endocrine, and pharmacological studies.
| Acknowledgments |
|---|
| Footnotes |
|---|
Disclosure Statement: The authors have nothing to disclose.
First Published Online December 20, 2007
Abbreviations: GR, Glucocorticoid receptor; GRE, glucocorticoid response element; hGR, human GR; PBST, PBS with Tween 20; PFA, paraformaldehyde; YFP, yellow fluorescent protein; zGR, zebrafish GR.
Received October 3, 2007.
Accepted for publication December 10, 2007.
| References |
|---|
|
|
|---|
by interfering with formation of active coactivator complexes. Mol Endocrinol 19:52–64This article has been cited by other articles:
![]() |
T. M. Rodela, M. D. McDonald, P. J. Walsh, and K. M. Gilmour The regulatory role of glucocorticoid and mineralocorticoid receptors in pulsatile urea excretion of the gulf toadfish, Opsanus beta J. Exp. Biol., June 15, 2009; 212(12): 1849 - 1858. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. H Stolte, A. F de Mazon, K. M Leon-Koosterziel, M. Jesiak, N. R Bury, A. Sturm, H. F J Savelkoul, B M L. V. van Kemenade, and G. Flik Corticosteroid receptors involved in stress regulation in common carp, Cyprinus carpio J. Endocrinol., August 1, 2008; 198(2): 403 - 417. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |