Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Porcine SOX9 Gene Expression Is Influenced by an 18bp Indel in the 5’-Untranslated Region

  • Bertram Brenig ,

    Contributed equally to this work with: Bertram Brenig, Yanyu Duan

    bbrenig@gwdg.de

    Affiliation Institute of Veterinary Medicine, Georg-August-University, Burckhardtweg 2, D-37077, Göttingen, Germany

    ⨯
  • Yanyu Duan ,

    Contributed equally to this work with: Bertram Brenig, Yanyu Duan

    Affiliation Key Laboratory for Animal Biotechnology of Jiangxi Province and the Ministry of Agriculture of China, Jiangxi Agricultural University, 330045, Nanchang, China

    ⨯
  • Yuyun Xing,

    Affiliation Key Laboratory for Animal Biotechnology of Jiangxi Province and the Ministry of Agriculture of China, Jiangxi Agricultural University, 330045, Nanchang, China

    ⨯
  • Nengshui Ding,

    Affiliation Key Laboratory for Animal Biotechnology of Jiangxi Province and the Ministry of Agriculture of China, Jiangxi Agricultural University, 330045, Nanchang, China

    ⨯
  • Lusheng Huang,

    Affiliation Key Laboratory for Animal Biotechnology of Jiangxi Province and the Ministry of Agriculture of China, Jiangxi Agricultural University, 330045, Nanchang, China

    ⨯
  • Ekkehard Schütz

    Affiliation Institute of Veterinary Medicine, Georg-August-University, Burckhardtweg 2, D-37077, Göttingen, Germany

    ⨯

Abstract

Sex determining region Y-box 9 (SOX9) is an important regulator of sex and skeletal development and is expressed in a variety of embryonal and adult tissues. Loss or gain of function resulting from mutations within the coding region or chromosomal aberrations of the SOX9 locus lead to a plethora of detrimental phenotypes in humans and animals. One of these phenotypes is the so-called male-to-female or female-to-male sex-reversal which has been observed in several mammals including pig, dog, cat, goat, horse, and deer. In 38,XX sex-reversal French Large White pigs, a genome-wide association study suggested SOX9 as the causal gene, although no functional mutations were identified in affected animals. However, besides others an 18bp indel had been detected in the 5′-untranslated region of the SOX9 gene by comparing affected animals and controls. We have identified the same indel (Δ18) between position +247bp and +266bp downstream the transcription start site of the porcine SOX9 gene in four other pig breeds; i.e., German Large White, Laiwu Black, Bamei, and Erhualian. These animals have been genotyped in an attempt to identify candidate genes for porcine inguinal and/or scrotal hernia. Because the 18bp segment in the wild type 5′-UTR harbours a highly conserved cAMP-response element (CRE) half-site, we analysed its role in SOX9 expression in vitro. Competition and immunodepletion electromobility shift assays demonstrate that the CRE half-site is specifically recognized by CREB. Both binding of CREB to the wild type as well as the absence of the CRE half-site in Δ18 reduced expression efficiency in HEK293T, PK–15, and ATDC5 cells significantly. Transfection experiments of wild type and Δ18 SOX9 promoter luciferase constructs show a significant reduction of RNA and protein levels depending on the presence or absence of the 18bp segment. Hence, the data presented here demonstrate that the 18bp indel in the porcine SOX9 5′-UTR is of functional importance and may therefore indeed be a causative variation in SOX9 associated traits.

Citation: Brenig B, Duan Y, Xing Y, Ding N, Huang L, Schütz E (2015) Porcine SOX9 Gene Expression Is Influenced by an 18bp Indel in the 5’-Untranslated Region. PLoS ONE 10(10): e0139583. https://doi.org/10.1371/journal.pone.0139583

Editor: Christoph Englert, Leibniz Institute for Age Research - Fritz Lipmann Institute (FLI), GERMANY

Received: July 10, 2015; Accepted: September 15, 2015; Published: October 2, 2015

Copyright: © 2015 Brenig et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper.

Funding: The authors have no support or funding to report.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Sry (sex determining region Y)-box 9 (SOX9) belongs to the SoxE subgroup of Sox family proteins and is expressed during embryonal development and adult life in meso-, ecto- and endoderm derived tissues [1]. It is involved in numerous cellular processes, e.g. chrondrogenesis [2], sex determination [3], pigmentation [4], organ maintenance [1], limb development [5], and cancer [6]. In the mesoderm SOX9 is involved in chondrogenesis and skeletal development [7], male gonad development, as well as development of cardiac valves and septa as well as epithelial differentiation in the pyloric sphincter [8]. However, SOX9 does not only play an important role in embryonal development. Recent data demonstrate that it is also important in the maintenance of adult organs and expression of SOX9 can be detected in adult stem and progenitor cells [9]. In mice it was shown that SOX5, SOX9, and SOX13 are expressed in adult Leydig cells and may therefore contribute to steroido- and spermatogenesis in postnatal testes [10]. In rats SOX9 expression was detected in the adult testicular cords and seminiferous tubuli suggesting a role in further germ cell differentiation [11]. SOX9 expression was also detected in goat testis in postnatal development, however, expression levels decreased to less than 50% of the concentration measured at two months of age [12].

Despite high levels of expression in both chondrogenic tissue and gonads, SOX9 is also transcribed to varying degrees in other tissues, including human adolescent heart, brain, kidney, muscle, colon, and cranial neural crest [13]. This suggests that SOX9 has other crucial functions not only in chondrogenesis and sex determination. For example, SOX9 supports tumor growth and invasion, regulates CEACAM1 expression in colon epithelium and plays a role in cranial neural crest development [14, 15].

Because of its wide range of interactions and functions, it is not surprising that mutations of the SOX9 gene locus are causative for a variety of defects in humans including campomelic dysplasia with or without sex reversal [16, 17], Pierre Robin sequence [18], Cooks syndrome [19], 46,XY gonadal dysgenesis [20], 46,XX male sex reversal, and congenital generalized hypertrichosis with or without hyperplasia [21]. Male-to-female or female-to-male sex reversal has also been described in animals [22–26]. In a Sry-negative XX European roe deer three SOX9 copies were detected leading to an incomplete male-determination. The sex reversal was presumably due to a dosage effect. In earlier studies of canine XX sex reversal SOX9 was initially excluded as candidate gene, however, recent reports show a SOX9 duplication resulting in an overexpression [25, 27]. In pigs female-to-male sex reversal of 38,XX animals has been described in several studies and it was shown that SOX9 expression is elevated in XX sex-reversed or intersex gonads [23, 24]. Recently, a genome-wide association study performed in the French Large White population demonstrated that the only significantly associated SNPs clustered around the SOX9 locus [28]. Comparative sequencing of the candidate region in affected animals and controls revealed 14 different polymorphisms. Unfortunately, all of these were located outside of the exons or splice-sites and therefore were questioned as functional candidate mutations [28]. However, at least three haplotypes were deduced that were more frequently present in the affected animals. These haplotypes included polymorphisms located in important regulatory regions as well as the 5′- and 3′-UTR of SOX9. One of the regulatory regions is the so-called TES (testis specific enhancer) approx. 13 kb upstream the transcription start site. It was shown in mice that within this region a highly conserved core element TESCO is bound by different transcription factors resulting in either up or down regulation of SOX9 [3, 29]. But SOX9 expression regulation is complex and under the control of further distant located elements. For example, eight (E1–E8) evolutionary conserved elements have been identified by comparative analysis. Five of these elements are dispersed in a region 290 kbp upstream and three up to 452 kbp downstream of SOX9. Two further elements were localized in the 3′-UTR of SOX9 [30]. In recent experiments using chromosome conformation capture-on-chip analysis even more distant regulatory regions influencing SOX9 expression have been identified [31]. These regions are located 2.46 Mb upstream as well as 1.22 Mb and 1.6 Mb downstream of SOX9. SOX9 also regulates its own expression in a positive feedback mechanism by the so-called SOM, an enhancer identified in mice 70 kb upstream of Sox9 [32]. In this respect it is noteworthy that a potential SOX9 consensus binding site [33] can be found approx. 17.8 kb upstream the porcine SOX9 transcription start site. Besides the long-range regulation, SOX9 expression is, of course, also regulated by promoter elements. In the murine Sox9 promoter an interval between -193 and -73 was identified to be essential for maximal promoter activity and tissue-specific expression in mouse cell lines [34]. Subsequently, several studies reported that functional elements in this region, i.e. two functional CCAAT boxes, a cAMP-response element (CRE) half-site and two specificity protein 1 (SP1) binding sites within 150 bp upstream the transcription start site (TSS), regulated SOX9 transcription in mouse and human chondrogenic cell lines and primary chondrocytes [35, 36]. In humans it was shown that SOX9 expression is regulated by Sp1, CREB, and CBF [37, 38]. These binding sites are located in the proximal promoter upstream the transcription start site. In vitro mutagenesis of these sites results in a reduction of SOX9 promoter activity. To define the regulatory network in response to different hormones and cellular signals, it is important to characterize the cis-regulatory elements even in the minimal promoter region. Functional elements located between -256 and +67 are important for controlling human SOX9 transcription efficiency [35]. However, transcription elements in the interval between -73 and +251 determine the minimal transcriptional activity of the mouse Sox9 promoter [34]. Whether additional regulatory elements are influencing SOX9 transcription downstream the transcription start site has not been shown so far.

In this context, it is noteworthy that in the French Large White population mentioned above, an 18bp indel was detected in the 5′-UTR in three of the sequenced control animals [28]. In this contribution we report that the same 18bp indel is also present in other pig breeds, i.e. German Large White, Laiwu Black, Bamei, and Erhualian. The occurrence of the 18bp indel in such diverse pig breeds in conjunction with the results reported about the porcine 38,XX sex reversal phenotype prompted us to address the question whether this region has an influence on SOX9 expression and might therefore at least partly explain the sex reversal. This hypothesis is further supported by the fact that the 18bp indel region harbours a CRE half-site which is usually recognized by the cAMP-response element binding protein CREB [39].

Materials and Methods

Animals

Seven Chinese domestic pig breeds (Erhualian, Laiwu Black, Bamei, Wuzhishan, Hang, Jianhexiang White, Tibetan) and different Western commercial pig breeds (Pietrian, German Landrace, German Large White, Red Duroc) were investigated. Blood samples were drawn by veterinarians as part of routine diagnostic procedures (parentage control, epidemiological testing) with informed owner consent, therefore the study was exempt of ethical approval according to the German regulations. A total of 938 animals were used for the identification of polymorphisms in the SOX9 gene. The complete SOX9 gene including 2 kb up- and downstream sequences was analysed.

PCR fragment length polymorphism (PFLP) analysis

PFLP analysis was used to screen all samples using the following PCR primers: forward 5'-GCCAGTTTTACCCCCAGGA–3' and reverse 5'-AGCGGCTCCCGGGAA GCCT–3'. The amplicon has a length of 189bp and 207bp for the Δ18 and wild type alleles, respectively. PCR was performed in a final 20 μL reaction mix contained 40 ng DNA, 1 × PCR buffer, 2.0 mM MgCl2, 0.2 mM of each dNTP, 10% DMSO, 10 μmol of each primer and 1 unit of ExTaq DNA polymerase (Takara). The PCR profile was as follows: 3 min at 94°C; 35 cycles of 30 s at 94°C and 1 min at 68°C; and a final extension of 7 min at 72°C. PCR products were evaluated by electrophoresis on 2% agarose gels.

Tissue culture

The cell lines used for the experiments were selected based on their reported transfection efficiencies and level of SOX9 expression. Nuclear extract from ESK–4 cells (female embryonal porcine kidney) were used for the immunodepletion EMSA in order to keep within the autologous porcine system. In addition it was shown previously that SOX9 is expressed in the embryonal kidney and therefore the presence of essential regulatory factors was assumed [40]. No data on SOX9 expression in PK–15 are available. Expression of SOX9 has been shown in ATDC5 [41]. In HEK293T SOX9 was not detectable using antibody NBP1-85551 (Novus Biologicals). However, due to a lower transfection efficiency of porcine epithelial PK–15 and mouse pre-chondrogenic ATDC5 cell lines, we also used human kidney HEK293T, which is known to have a higher transfection efficiency [41–44].

HEK293T, PK–15, and ATDC5 cells were cultured in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum and 2 mM glutamine at 37°C in a humidified atmosphere containing 5% CO2. Porcine kidney ESK–4 cells were purchased from the European Collection of Cell Cultures (ECACC) and cultured according to the manufacturer’s recommended protocols.

Affinity purification of CREB

DNA-binding affinity was analysed as described [45] using three oligonucleotides containing either the CRE half-site (WT: 5′-CGCTGCAGCCGAGTGACGCGCCAG GCTTCC–3′), a mutated CRE half-sitemut (MT: 5′-CGCTGCAGCCGAGCTCTACGC CAGGCTTCC–3′) or the consensus CRE binding site (CRE: 5′-AGAGGATTGCCT GACGTCAGAGAGCTAG–3′). HEK293T cells were transfected using HA-tag CREB-expression plasmid (pcDNA4-HA-CREB, see below) and nuclear and cytoplasmic extracts were prepared after 32 h. Extracts were incubated with the biotinylated probes (WT, MUT, CRE) and Streptavidin MagneSphere Paramagnetic particles (Promega) for 2 h at room temperature. The reactions were separated using a MagneSphere Magnetic Separation Stand (Promega) and washed according to the manufacturer’s instructions. Eluted extracts were separated on SDS/PAGE and analysed by Western blotting using a phospho-CREB antibody (Millipore).

Electromobility shift assay (EMSA)

Nuclear and cytoplasmic extracts were prepared from sub-confluent ESK–4 cells by NE-PER Nuclear and Cytoplasmic Extraction Kit (PIERCE) according to the manufacturer’s instructions. The DNA-protein interaction was studied by EMSA following the instructions of Chemiluminescent Nucleic Acid Detection Module (PIERCE) with minor modifications as previously described [46]. Briefly, binding reactions included 2 μL 10 x binding buffer (100 mM Tris, 500 mM KCl, 10 mM DTT; pH 7.5), 5% BSA, 0.5 mM EDTA, 5% glycerol, 1 mM MgCl2, 100 ng poly-dAdT, 0.1% NP–40 and 4–4.5 μg nuclear and cytoplasmic proteins. Binding reactions were pre-incubated for 20 min on ice followed by 20 min at room temperature after adding 20 fmol biotin-end labelled double strand oligonucleotide probe (wild type: 5′-gcgccttcctaaaagctcgccgcTGCAGCCGAGTGACGCGCcaggcttcccgggagcc–3′; Δ18: gcgccttcctaaaagctcgccgccaggcttcccgggagcc–3′). Finally, the mixtures were loaded onto 7% acrylamide–0.5x Tris-borate-EDTA gel (acrylamide/bisacrylamide = 37.5:1). For competition analysis, 2–5 pmol unlabelled competitors, i.e. the 18-bp fragment, were added to the binding reactions. For interrogation of the binding transcription factor, 3 μL phospho-CREB antibody (Millipore) or E2F4 antibody (Millipore) was pre-incubated with nuclear and cytoplasmic proteins on ice for 2 h to get the maximal signal. The mix was added to the above binding reactions for 20 min on ice prior to the addition of the labelled DNA probe.

Plasmid constructs

To analyse the effect on expression of CREB binding to the wild type 5′-UTR pcDNA4-HA-CREB, pGL4.20-WT (wild type) and pGL4.20-Δ18 (mutant) were constructed. Full-length human CREB cDNA (Sino Biological Inc.) was subcloned into the pcDNA4-HA expression vector. pGL4.20-WT and pGL4.20-Δ18 was generated by subcloning the porcine SOX9 promoter between position -100 and +399 (TSS = +1) into the KpnI/HindIII sites of the pGL4.20 vector. The promoter was obtained by artificial gene synthesis (Sangon).

To determine the effect of Δ18 on translation, two plasmids (psiCHECK-2-WT-TTG and psiCHECK-2-Δ18-TTG) were constructed by inserting the 5’-UTR plus the start codon of SOX9 into an NheI site upstream the Renilla luciferase gene. Before inserting the DNA fragments into the psiCHECK–2 vector, the start codon of the Renilla luciferase was mutated to TTG by PCR-driven overlap extension. By this way, the Renilla luciferase translation would be driven by the start codon of SOX9. The 5'-UTR of SOX9 was PCR amplified using the primers: forward 5'-TTTTGCTAG CAAGAGCCCCTGGGCTGGGAGTTGG–3' and reverse 5'-TTTTGCTAGCATCAC GCAGGCCCGGGGCAGGGGAC–3'. Amplified products were cloned and ligated into the NheI site directly preceding the Renilla luciferase gene in the plasmid psiCHECK-2-TTG. The psiCHECK2-TTG vectors are ideal for examining the effect of 5’-UTRs on gene expression as described [47, 48]. By fusing the 5’-UTR of interest to the TTG mutated Renilla luciferase gene, luciferase activity can be used as a marker for 5’-UTR regulation.

Luciferase reporter assays

Cells were seeded at 2 x 105 cells per well in 24-well dishes. When the cells were 80% confluent, the luciferase constructs were transfected using LipofectamineTM 2000 (Invitrogen) according to the manufacturer’s instructions. Thirty-six hours after transfection, cells were harvested and assayed for luciferase activity using a dual-luciferase reporter assay kit (Promega) following the manufacturer’s manual. Results are given as relative light units (RLU) that is Firefly luciferase value/Renilla luciferase value. The average RLU of cells co-transfected with pGL4.20 and pRL-TK (Renilla luciferase control vector, Promega) was considered to be 1. In each experiment 0.2 μg of pGL4.20, pGL4.20-WT, or pGL4.20-Δ18 was added together with 50 ng pRL-TK without or with 0.8 μg pcDNA4-HA-CREB. The optimal amount and ratio of reporter and control plasmids have been determined in preliminary experiments. The amount of pRL-TK was adjusted according to previous reports [49–51].

Student's t-test was used to calculate significance (* p<0.05, ** p<0.01, *** p<0.001). All experiments were performed in triplicate.

Real-time quantitative PCR assays

Total RNA was extracted with the TRIZOL Reagent (Invitrogen) according to the manufacturer’s instructions. 6 μg of RNA was transcribed into first-strand cDNA using Superscript First-Strand Synthesis System for RT-qPCR (Takara). The relative expression levels were analysed using Brilliant II FAST SYBR QPCR Master Mix (Stratagene) in an MX3000P system (ABgene) according to the manufacturer’s instructions. The following primers were used: Renilla 1, 5′-CTGATCTGATCGGAA TGGGT–3′ and 5′-GACGATGGCCTTGATCTTGT–3′; ß-actin 1, 5′-AAGGAGAAG CTGTGCTACGTCG–3′ and 5′-TGAAGGTAGTTTCGTGGATGCC–3′.

Statistical analysis

Putative transcription factor binding sites within Δ18 were analysed with MatInspector (http://www.genomatrix.de/matinspector) using default settings. The DNA sequence between the transcription start site and the start codon was used as input. The promoter alignment analysis of evolutionary relationship among pig and other vertebrates was performed using ClustalW2 (http://www.ebi.ac.uk/). Expression data were analysed using GraphPad PrismTM 5.0 for Macs (GraphPad Software, San Diego, CA).

All transfection experiments were performed in at least three independent experiments using three biological replicates and standard errors (SEM) were calculated using Microsoft Excel for Mac 2011 (version 14.5.4).

Results

CREB binds to the CRE half-site within the 5'-UTR of the porcine SOX9 gene and influences expression

In an analysis of 938 animals of different commercial and Chinese pig breeds 12 polymorphisms (SNP) were identified within in the 5′-UTR (2/12), introns (5/12), and 3′-UTR (5/12) of SOX9 (Table 1). In addition an 18bp indel (Δ18) abolishing an interval between +247bp and +266bp downstream the TSS in four of the analysed pig breeds, i.e. German Large White, Laiwu Black, Bamei, and Erhualian, was identified (Table 2). As mentioned above we previously studied Δ18 in 545 animals of the commercial breeds in an attempt to identify an association with inguinal and scrotal hernia. However, no significant association could be detected using transmission disequilibrium test (p = 0.075) [52, 53]. All SNPs including Δ18 have been already deposited either with dbSNP or described previously as mentioned in the introduction [28].

thumbnail
Download:
Table 1. Polymorphisms detected in the porcine SOX9 gene.

https://doi.org/10.1371/journal.pone.0139583.t001

thumbnail
Download:
Table 2. Genotype distribution of the 18-bp variant in different commercial breeds, seven Chinese domestic breeds and three Western commercial breeds.

https://doi.org/10.1371/journal.pone.0139583.t002

DNA sequence comparison between pig, cattle, buffalo, roe deer, and horse, demonstrates that the Δ18 region is highly conserved (Fig 1) and harbours three potential transcription factor binding sites with varying degrees of similarity to the respective consensus motives (Table 3). Due to the highly conserved sites and potential transcription factor binding sites, we hypothesized that this region might be important for SOX9 expression. To test CREB binding to the 18bp region in vivo, HEK293T cells were transfected with the CREB-expression plasmid (pcDNA4A-HA-CREB) and nuclear and cytoplasmic extracts were incubated with oligonucleotides harbouring the 18bp region (WT) and a CRE consensus sequence (CRE). As control an oligonucleotide homologous to the 18bp region with a mutated CRE half-site (MUT) was used. As shown in Fig 2A CREB binds to the WT and CRE oligonucleotides, whereas no binding was detected to the mutated control oligonucleotide. Specificity of CREB binding was analysed with competition and immunodepletion EMSAs. Nuclear and cytoplasmic extracts were prepared from ESK–4 cells and incubated with the labelled wild type DNA fragment (Fig 2B). As shown in Fig 2B, a specific band (BP) was detectable and diminished with increasing concentrations of cold unlabelled double strand homologous competitor. To exclude other factors which could bind to this region, we used antibodies directed against another potential transcription factor listed in Table 3, i.e. E2F, pre-incubated the cell extracts with the antibodies and then performed an immunodepletion EMSA as described elsewhere [54]. With the interference analysis we expected that due to binding of the antibody less transcription factor would be available to bind the labelled oligonucleotide and therefore would result in a reduced band intensity. Two antibodies were used, directed against CREB and E2F4. Fig 2C shows that when a polyclonal CREB antibody (anti-phospho-CREB) was added to the binding reaction the intensity of the protein-DNA complex was significantly reduced (lane 2) whereas an anti-E2F4 antibody did not influence the band intensity (lane 3). No binding was detected with the Δ18 fragment (lane 4). The band intensities were densitometrically quantified and differed significantly between lane 2 and 3 (Fig 2C lower panel).

thumbnail
Download:
Fig 1. Organization and location of the proximal promoter of the porcine SOX9 gene.

The 18bp sequence and the CRE half-site are indicated. The sequence alignment in the lower panel shows conserved regions between pig, cattle, buffalo, roe deer, and horse. (B) Sequence alignment of the wild type and 18bp indel (Δ18) region.

https://doi.org/10.1371/journal.pone.0139583.g001

thumbnail
Download:
Fig 2. CREB binds the CRE half-site in the 18bp region.

(A) Biotinylated oligonucleotides WT, MUT, and CRE were incubated with cell extracts of HEK293T cells transfected with a CREB expression plasmid. Bound protein was purified using Streptavidin MagneSphere Paramagnetic particles and a MagneSphere Magnetic Separation Stand (Promega), separated on an SDS/PAGE and Western blotted. Western blots were incubated with an HA-CREB antibody. Intensities of bands were densitometrically measured using ImageJ software (http://imagej.nih.gov/ij/index.html). (B) A constant amount of nuclear and cytoplasmic proteins of ESK–4 cells and labelled probe (WT) were incubated with increasing amounts of homologous unlabelled competitor. DsFP, double stranded free probe; ssFP: single stranded free probe; BP, bound protein-DNA complexes. (C) Immunodepletion assays were performed using anti-phospho-CREB and anti-E2F4 antibodies. Density of the resulting bands was determined using ImageJ software.

https://doi.org/10.1371/journal.pone.0139583.g002

thumbnail
Download:
Table 3. Putative transcription factor binding sites within the 5'-UTR between positions 66,414 to 66,446 of the porcine SOX9 gene.

https://doi.org/10.1371/journal.pone.0139583.t003

To analyse whether CREB binding and the Δ18 variant in the context of the 5'-UTR have regulatory effects on SOX9 expression the wild type and the Δ18 5′-UTR were cloned into the pGL4.20 vector and used for transfection of HEK293T cells with or without a CREB expression plasmid. As shown in Fig 3 (left panel, -CREB) the wild type 5′-UTR resulted in a significantly higher activity than the Δ18 fragment in the absence of CREB. When CREB was co-expressed (Fig 3 right panel, +CREB) the activity of the wild type 5´-UTR was reduced to the level of the Δ18 fragment and both fragments resulted in a significant lower expression than the CREB induced expression of the control vector.

thumbnail
Download:
Fig 3. Binding of CREB to the wild type 5′-UTR and Δ18 results in reduced expression.

HEK293T cells were transiently co-transfected with pGL4.20-WT, pGL4.20-Δ18 with or without pcDNA4A-HA-CREB over-expression vector. Relative luciferase activity was measured after 32 h using dual luciferase reporter assays. The results are given as relative light units (RLU). The average RLU of cells co-transfected with pGL4.20, pRL Renilla luciferase control reporter vector, and pcDNA4A-HA was considered to be 1 (error bars indicate standard deviations).

https://doi.org/10.1371/journal.pone.0139583.g003

Effects of the variable SOX9 5'-UTR on transcription and translation in different cell types

Three expression vectors were constructed, i.e. psiCHECK-2-TTG, -WT-TTG and Δ18-TTG and transfected into HEK293T, PK–15, and ATDC5 cells. As shown in Fig 4A the Renilla luciferase start codon was mutated to TTG (psiCHECK-2-TTG) to abolish translation from this site. As the 5'-UTR of porcine SOX9 with a length of 394 bp does not contain multiple open reading frames or other upstream AUGs except the natural start codon, we cloned the two different variants of the SOX9 5'-UTR into the NheI cloning site of the psiCHECK-2-TTG vector to generate psiCHECK-2-WT-TTG and psiCHECK-2-Δ18-TTG. In these constructs translation will be initiated from the original SOX9 start codon -12 bp upstream with the addition of 4 amino acids at the N-terminus of the Renilla luciferase. Although the N-terminus of the luciferase is important for the stability of the enzyme [55], it has been shown that the addition of amino acids, e.g. Myc epitope or His-tag, does not influence the activity [56, 57].

thumbnail
Download:
Fig 4. Δ18 reduces translation efficiency.

(A) Schematic diagram showing the luciferase plasmids used for translation efficiency analysis. The start codon of Renilla was mutated to TTG. Two plasmids (psiCHECK-2-WT-TTG and –Δ18-TTG) were constructed by inserting the 5’-UTR plus the start code of SOX9 into upstream of the mutated Renilla luciferase gene. By this way, the Renilla luciferase expression would be driven by the primary initiation codon of SOX9. (B) The deletion decreases the translation efficiency of porcine SOX9. The four different plasmids were transfected into HEK293T, PK–15 and ATDC5 cells. After 24 h, both Firefly and Renilla luciferase values were measured. The ratio of Renilla and Firefly luciferase value were used as RLU. The average RLU for psiCHECK-2-TTG was considered to be 1.

https://doi.org/10.1371/journal.pone.0139583.g004

In all three transfected cell types the luciferase activity was almost reduced to background levels compared to the intact expression vector demonstrating the effect of the mutated Renilla start codon (Fig 4B). Renilla luciferase activity of psiCHECK-2-TTG was then used to assess the effect of the 5’-UTR on translation. Both constructs harbouring the wild type (WT) or Δ18 were used for transfection of HEK293T, PK–15, and ATDC5 cells. In all three cells types, the Δ18 showed a significant lower luciferase activity compared to the wild type 5′-UTR construct. In addition, Renilla transcripts were analysed by qRT-PCR. As shown in Fig 5 the wild type 5′-UTR resulted in a 3-fold higher RNA concentration compared to Δ18.

thumbnail
Download:
Fig 5. Determination of wild type and Δ18 5′-UTR transcription regulation.

RT-qPCR was performed for quantifying Renilla mRNA contents. The content of Renilla was normalized to ß-actin as internal control and the relative luciferase mRNA for psiCHECK-2-TTG was fixed as 1.

https://doi.org/10.1371/journal.pone.0139583.g005

Discussion

In the present study we characterized the effects of an 18bp deletion within the 5’-UTR of porcine SOX9 on transcription and translation. A sequence comparison of the 18bp region between pig and other species showed several conserved nucleotide positions including a CRE half-site. To determine whether the CRE half-site was recognized by CREB, we performed affinity purification, competition and immunodepletion mobility shift assays and were able to show binding. CREB did functionally bind to the CRE half-site (CGTCA/TGACG), albeit with reduced binding affinity compared to the full palindrome (TGACGTCA) [58]. CREB seems to be an important transcription factor in the regulation of bone development by activating Sox9 in mouse. Overexpression of Creb enhanced the action of Sox9 in the PKA pathway in BMP-2-induced osteochondrogenic differentiation [59]. Potent dominant negative Creb transgenic mice showed short-limbed dwarfism along with a defect in chondrocyte proliferation and a delay in chondrocyte hypertrophy in endochondral bone formation [60]. Other sites could also interfere with the expression, for example, an additional conserved CRE half-site at position -147 in the human SOX9 promoter [37]. However, CREB binding in the porcine SOX9 gene 5′-UTR seems to down-regulate expression as shown in Fig 3. Interestingly, CREB co-expression clearly up-regulated the promoterless control vector pGL4.20. Although it has been reported that pRL-TK used as normalizer is up-regulated by different transcription factors, e.g. GATA transcription factors [49], 12S E1A [61], SP1 [51], Nurr1 [50], and dexamethasone [62], there are no data available that CREB influences pRL-TK expression per se. However, as it is known that the TK promoter is influenced by CREB binding [63], we also expected that pRL-TK will be up-regulated and therefore took this into account when normalizing the expression levels. Hence, the increased pGL4.20 expression in fact seems to be due to a CREB dependent induction. Although this was an unexpected and so far unknown effect of CREB on pGL4.20 expression, promotorless luciferase reporter vectors pGL3 and pGL4 have been reported to be at least responsive to steroid hormones [64].

Inhibition of expression by CREB has been shown for other genes, e.g. MuSK, AP–2α, and PPAR-γ [65–67]. Truncations or mutations in UTRs often defeat fine regulation, cause impaired protein synthesis and finally associate with various diseases or disease susceptibility in humans [68]. Two studies showed that microRNA–145 directly repressed SOX9 expression by binding a unique site in the 3’-UTR resulting in a profound change in the human chondrocyte phenotype [69, 70]. It is also known that a stable secondary structure, multiple ORFs, and uAUGs in the 5’-UTR largely determine translation efficiency [71]. However, the 5’-UTR of the porcine SOX9 has only one ORF and hence the downregulation of SOX9 expression seems to be due only to either CREB binding or the presence of Δ18. CREB binding resulted in a reduction of SOX9 transcription as shown in Fig 3, whereas Δ18 could have an effect on either RNA stability or RNA degradation. A comparative in silico analysis of the secondary structure formation of the wild type and Δ18 5′-UTR showed that wild type RNA had a lower ΔG (-170.77) than Δ18 (-161.78) (data not shown).

Although regulation of SOX9 expression is complex and the role of a single site must be interpreted with caution, our results are in agreement with the findings that the wild type 5′-UTR was identified more often in 38,XX sex-reversal pigs and SOX9 is over-expressed in XX intersex gonads [23, 28]. As shown in Figs 4 and 5 the presence of the 18bp region resulted in a significant increase of promoter activity compared to Δ18. Only, when CREB was overexpressed and bound to its binding site a reduction of expression was observed (Fig 3). From these findings different scenarios for the development of a 38,XX sex-reversal phenotype can be hypothesized. In the presence of the wild type genotype a sex reversal phenotype could result from a low CREB expression and a more stable SOX9 RNA secondary structure leading to high SOX9 concentrations and consequently to a male phenotype. On the other hand, if CREB expression is high enough to down-regulate SOX9 expression, no sex reversal will be induced. In the presence of the Δ18 genotype, where SOX9 expression was low independent of CREB and the RNA secondary structure was less stable, a normal sex development would be expected for 38,XX individuals.

Hence, in further experiments it will be interesting to analyse the expression of CREB in controls and 38,XX sex-reversal gonads.

Conclusion

We have shown that an 18bp deletion located +248 bp downstream the TSS in the 5’-UTR of the porcine SOX9 gene harbours a functional CRE half-site. The deletion segregates in several pig breeds. The 18bp region is essential for transcription and translation efficiency of porcine SOX9 and interacts with CREB. CREB binding to the wild type or presence of the 18bp deletion results in a significant reduction of SOX9 transcription and translation.

Acknowledgments

The authors are thankful to Ch. Knorr for support and encouragement at the beginning of the project. S. Pach is thanked for expert technical assistance.

Author Contributions

Conceived and designed the experiments: BB YD YX ND. Performed the experiments: YD YX ND. Analyzed the data: BB YD YX ND ES. Contributed reagents/materials/analysis tools: BB LH. Wrote the paper: BB YD ES.

References

  1. 1. Jo A, Denduluri S, Zhang B, Wang Z, Yin L, Yan Z, et al. The versatile functions of Sox9 in development, stem cells, and human diseases. Genes Dis. 2014;1(2):149–61. pmid:25685828; PubMed Central PMCID: PMC4326072.
  2. 2. Wuelling M, Vortkamp A. Chondrocyte proliferation and differentiation. Endocrine development. 2011;21:1–11. Epub 2011/08/26. pmid:21865749.
  3. 3. Jakob S, Lovell-Badge R. Sex determination and the control of Sox9 expression in mammals. FEBS J. 2011;278(7):1002–9. pmid:21281448.
  4. 4. Passeron T, Valencia JC, Bertolotto C, Hoashi T, Le Pape E, Takahashi K, et al. SOX9 is a key player in ultraviolet B-induced melanocyte differentiation and pigmentation. Proceedings of the National Academy of Sciences of the United States of America. 2007;104(35):13984–9. pmid:17702866; PubMed Central PMCID: PMC1955778.
  5. 5. Raspopovic J, Marcon L, Russo L, Sharpe J. Modeling digits. Digit patterning is controlled by a Bmp-Sox9-Wnt Turing network modulated by morphogen gradients. Science. 2014;345(6196):566–70. pmid:25082703.
  6. 6. Harris ML, Baxter LL, Loftus SK, Pavan WJ. Sox proteins in melanocyte development and melanoma. Pigment Cell Melanoma Res. 2010;23(4):496–513. pmid:20444197; PubMed Central PMCID: PMC2906668.
  7. 7. Leung VY, Gao B, Leung KK, Melhado IG, Wynn SL, Au TY, et al. SOX9 governs differentiation stage-specific gene expression in growth plate chondrocytes via direct concomitant transactivation and repression. PLoS genetics. 2011;7(11):e1002356. pmid:22072985; PubMed Central PMCID: PMC3207907.
  8. 8. Moniot B, Biau S, Faure S, Nielsen CM, Berta P, Roberts DJ, et al. SOX9 specifies the pyloric sphincter epithelium through mesenchymal-epithelial signals. Development (Cambridge, England). 2004;131(15):3795–804. pmid:15240557; PubMed Central PMCID: PMC2442161.
  9. 9. Purba TS, Haslam IS, Shahmalak A, Bhogal RK, Paus R. Mapping the expression of epithelial hair follicle stem cell-related transcription factors LHX2 and SOX9 in the human hair follicle. Exp Dermatol. 2015;24(6):462–7. pmid:25808706.
  10. 10. Daigle M, Roumaud P, Martin LJ. Expressions of Sox9, Sox5, and Sox13 transcription factors in mice testis during postnatal development. Mol Cell Biochem. 2015. pmid:26045173.
  11. 11. Frojdman K, Harley VR, Pelliniemi LJ. Sox9 protein in rat sertoli cells is age and stage dependent. Histochem Cell Biol. 2000;113(1):31–6. pmid:10664067.
  12. 12. Faucette AN, Maher VA, Gutierrez MA, Jucker JM, Yates DC, Welsh TH Jr., et al. Temporal changes in histomorphology and gene expression in goat testes during postnatal development. Journal of animal science. 2014;92(10):4440–8. pmid:25085396.
  13. 13. Akiyama H. Control of chondrogenesis by the transcription factor Sox9. Mod Rheumatol. 2008;18(3):213–9. Epub 2008/03/21. pmid:18351289.
  14. 14. Stuhlmiller TJ, Garcia-Castro MI. Current perspectives of the signaling pathways directing neural crest induction. Cell Mol Life Sci. 2012. Epub 2012/05/02. pmid:22547091.
  15. 15. Zalzali H, Naudin C, Bastide P, Quittau-Prevostel C, Yaghi C, Poulat F, et al. CEACAM1, a SOX9 direct transcriptional target identified in the colon epithelium. Oncogene. 2008;27(56):7131–8. Epub 2008/09/17. pmid:18794798.
  16. 16. Pritchett J, Athwal V, Roberts N, Hanley NA, Hanley KP. Understanding the role of SOX9 in acquired diseases: lessons from development. Trends Mol Med. 2011;17(3):166–74. pmid:21237710.
  17. 17. Lee YH, Saint-Jeannet JP. Sox9 function in craniofacial development and disease. Genesis. 2011;49(4):200–8. pmid:21309066; PubMed Central PMCID: PMC3079054.
  18. 18. Gordon CT, Tan TY, Benko S, Fitzpatrick D, Lyonnet S, Farlie PG. Long-range regulation at the SOX9 locus in development and disease. J Med Genet. 2009;46(10):649–56. pmid:19473998.
  19. 19. Brennan CB, Buehler T, Lesher JL Jr. Cooks syndrome: a case report and brief review. Pediatr Dermatol. 2013;30(4):e52–3. pmid:22329539.
  20. 20. Ostrer H. Disorders of sex development (DSDs): an update. J Clin Endocrinol Metab. 2014;99(5):1503–9. pmid:24758178.
  21. 21. Sun M, Li N, Dong W, Chen Z, Liu Q, Xu Y, et al. Copy-number mutations on chromosome 17q24.2-q24.3 in congenital generalized hypertrichosis terminalis with or without gingival hyperplasia. Am J Hum Genet. 2009;84(6):807–13. pmid:19463983; PubMed Central PMCID: PMC2694973.
  22. 22. Kropatsch R, Dekomien G, Akkad DA, Gerding WM, Petrasch-Parwez E, Young ND, et al. SOX9 duplication linked to intersex in deer. PloS one. 2013;8(9):e73734. pmid:24040047; PubMed Central PMCID: PMC3765313.
  23. 23. Pailhoux E, Parma P, Sundstrom J, Vigier B, Servel N, Kuopio T, et al. Time course of female-to-male sex reversal in 38,XX fetal and postnatal pigs. Developmental dynamics: an official publication of the American Association of Anatomists. 2001;222(3):328–40. pmid:11747069.
  24. 24. Pailhoux E, Popescu PC, Parma P, Boscher J, Legault C, Molteni L, et al. Genetic analysis of 38XX males with genital ambiguities and true hermaphrodites in pigs. Animal genetics. 1994;25(5):299–305. pmid:7818163.
  25. 25. Rossi E, Radi O, De Lorenzi L, Vetro A, Groppetti D, Bigliardi E, et al. Sox9 duplications are a relevant cause of Sry-negative XX sex reversal dogs. PloS one. 2014;9(7):e101244. pmid:25010117; PubMed Central PMCID: PMC4091935.
  26. 26. Ciotola F, Albarella S, Pasolini MP, Auletta L, Esposito L, Iannuzzi L, et al. Molecular and cytogenetic studies in a case of XX SRY-negative sex reversal in an Arabian horse. Sex Dev. 2012;6(1–3):104–7. pmid:22025175.
  27. 27. Kothapalli K, Kirkness E, Pujar S, Van Wormer R, Meyers-Wallen VN. Exclusion of candidate genes for canine SRY-negative XX sex reversal. J Hered. 2005;96(7):759–63. pmid:16267164.
  28. 28. Rousseau S, Iannuccelli N, Mercat MJ, Naylies C, Thouly JC, Servin B, et al. A genome-wide association study points out the causal implication of SOX9 in the sex-reversal phenotype in XX pigs. PloS one. 2013;8(11):e79882. Epub 2013/11/14. pmid:24223201; PubMed Central PMCID: PMC3819277.
  29. 29. Kimura R, Murata C, Kuroki Y, Kuroiwa A. Mutations in the testis-specific enhancer of SOX9 in the SRY independent sex-determining mechanism in the genus Tokudaia. PloS one. 2014;9(9):e108779. pmid:25265165; PubMed Central PMCID: PMCPMC4181316.
  30. 30. Knower KC, Kelly S, Harley VR. Turning on the male–-SRY, SOX9 and sex determination in mammals. Cytogenet Genome Res. 2003;101(3–4):185–98. 74336. pmid:14684982.
  31. 31. Smyk M, Szafranski P, Startek M, Gambin A, Stankiewicz P. Chromosome conformation capture-on-chip analysis of long-range cis-interactions of the SOX9 promoter. Chromosome Res. 2013;21(8):781–8. pmid:24254229.
  32. 32. Mead TJ, Wang Q, Bhattaram P, Dy P, Afelik S, Jensen J, et al. A far-upstream (-70 kb) enhancer mediates Sox9 auto-regulation in somatic tissues during development and adult regeneration. Nucleic acids research. 2013;41(8):4459–69. pmid:23449223; PubMed Central PMCID: PMCPMC3632127.
  33. 33. Mertin S, McDowall SG, Harley VR. The DNA-binding specificity of SOX9 and other SOX proteins. Nucleic acids research. 1999;27(5):1359–64. pmid:9973626; PubMed Central PMCID: PMC148324.
  34. 34. Kanai Y, Koopman P. Structural and functional characterization of the mouse Sox9 promoter: implications for campomelic dysplasia. Hum Mol Genet. 1999;8(4):691–6. Epub 1999/03/11. pmid:10072439.
  35. 35. Colter DC, Piera-Velazquez S, Hawkins DF, Whitecavage MK, Jimenez SA, Stokes DG. Regulation of the human Sox9 promoter by the CCAAT-binding factor. Matrix Biol. 2005;24(3):185–97. Epub 2005/05/24. pmid:15908194.
  36. 36. Piera-Velazquez S, Hawkins DF, Whitecavage MK, Colter DC, Stokes DG, Jimenez SA. Regulation of the human SOX9 promoter by Sp1 and CREB. Exp Cell Res. 2007;313(6):1069–79. Epub 2007/02/10. pmid:17289023; PubMed Central PMCID: PMC2118054.
  37. 37. Piera-Velazquez S, Hawkins DF, Whitecavage MK, Colter DC, Stokes DG, Jimenez SA. Regulation of the human SOX9 promoter by Sp1 and CREB. Experimental cell research. 2007;313(6):1069–79. Epub 2007/02/10. PubMed Central PMCID: PMCPmc2118054.
  38. 38. Colter DC, Piera-Velazquez S, Hawkins DF, Whitecavage MK, Jimenez SA, Stokes DG. Regulation of the human Sox9 promoter by the CCAAT-binding factor. Matrix biology: journal of the International Society for Matrix Biology. 2005;24(3):185–97. Epub 2005/05/24.
  39. 39. Montminy MR, Bilezikjian LM. Binding of a nuclear protein to the cyclic-AMP response element of the somatostatin gene. Nature. 1987;328(6126):175–8. pmid:2885756.
  40. 40. Kent J, Wheatley SC, Andrews JE, Sinclair AH, Koopman P. A male-specific role for SOX9 in vertebrate sex determination. Development (Cambridge, England). 1996;122(9):2813–22. pmid:8787755.
  41. 41. Murakami S, Kan M, McKeehan WL, de Crombrugghe B. Up-regulation of the chondrogenic Sox9 gene by fibroblast growth factors is mediated by the mitogen-activated protein kinase pathway. Proceedings of the National Academy of Sciences of the United States of America. 2000;97(3):1113–8. pmid:10655493; PubMed Central PMCID: PMCPMC15539.
  42. 42. Caron MM, Emans PJ, Cremers A, Surtel DA, Coolsen MM, van Rhijn LW, et al. Hypertrophic differentiation during chondrogenic differentiation of progenitor cells is stimulated by BMP–2 but suppressed by BMP–7. Osteoarthritis Cartilage. 2013;21(4):604–13. pmid:23353668.
  43. 43. Chen X, Shi X, Gan F, Huang D, Huang K. Glutamine starvation enhances PCV2 replication via the phosphorylation of p38 MAPK, as promoted by reducing glutathione levels. Vet Res. 2015;46:32. pmid:25879878; PubMed Central PMCID: PMCPMC4363047.
  44. 44. Guo H, Hao R, Wei Y, Sun D, Sun S, Zhang Z. Optimization of electrotransfection conditions of mammalian cells with different biological features. J Membr Biol. 2012;245(12):789–95. pmid:22836669.
  45. 45. Roman AC, Benitez DA, Carvajal-Gonzalez JM, Fernandez-Salguero PM. Genome-wide B1 retrotransposon binds the transcription factors dioxin receptor and Slug and regulates gene expression in vivo. Proceedings of the National Academy of Sciences of the United States of America. 2008;105(5):1632–7. pmid:18223155; PubMed Central PMCID: PMC2234196.
  46. 46. Hellman LM, Fried MG. Electrophoretic mobility shift assay (EMSA) for detecting protein-nucleic acid interactions. Nat Protoc. 2007;2(8):1849–61. Epub 2007/08/19. pmid:17703195; PubMed Central PMCID: PMC2757439.
  47. 47. Arora A, Dutkiewicz M, Scaria V, Hariharan M, Maiti S, Kurreck J. Inhibition of translation in living eukaryotic cells by an RNA G-quadruplex motif. RNA (New York, NY). 2008;14(7):1290–6. Epub 2008/06/03. PubMed Central PMCID: PMCPmc2441988.
  48. 48. Kubokawa I, Takeshima Y, Ota M, Enomoto M, Okizuka Y, Mori T, et al. Molecular characterization of the 5'-UTR of retinal dystrophin reveals a cryptic intron that regulates translational activity. Molecular vision. 2010;16:2590–7. Epub 2010/12/15. PubMed Central PMCID: PMCPmc3000234.
  49. 49. Ho CK, Strauss JF 3rd. Activation of the control reporter plasmids pRL-TK and pRL-SV40 by multiple GATA transcription factors can lead to aberrant normalization of transfection efficiency. BMC Biotechnol. 2004;4:10. pmid:15119957; PubMed Central PMCID: PMCPMC416485.
  50. 50. Matuszyk J, Ziolo E, Cebrat M, Kochel I, Strzadala L. Nurr1 affects pRL-TK but not phRG-B internal control plasmid in genetic reporter system. Biochem Biophys Res Commun. 2002;294(5):1036–9. pmid:12074581.
  51. 51. Osborne SA, Tonissen KF. pRL-TK induction can cause misinterpretation of gene promoter activity. Biotechniques. 2002;33(6):1240–2. pmid:12503307.
  52. 52. Spielman RS, McGinnis RE, Ewens WJ. Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet. 1993;52(3):506–16. pmid:8447318; PubMed Central PMCID: PMCPMC1682161.
  53. 53. Xing Y. Dissection of genetic variants affecting boar sperm quality and porcine inguinal/scrotal hernia. Göttingen: Cuvillier Verlag Göttingen; 2010. 110 p.
    • 54. Dyer RB, Herzog NK. Immunodepletion EMSA: a novel method to identify proteins in a protein-DNA complex. Nucleic acids research. 1995;23(16):3345–6. pmid:7667113; PubMed Central PMCID: PMC307198.
    • 55. Sung D, Kang H. The N-terminal amino acid sequences of the firefly luciferase are important for the stability of the enzyme. Photochem Photobiol. 1998;68(5):749–53. Epub 1998/11/24. pmid:9825705.
    • 56. Loening AM, Fenn TD, Wu AM, Gambhir SS. Consensus guided mutagenesis of Renilla luciferase yields enhanced stability and light output. Protein Eng Des Sel. 2006;19(9):391–400. Epub 2006/07/22. pmid:16857694.
    • 57. Woo J, Howell MH, von Arnim AG. Structure-function studies on the active site of the coelenterazine-dependent luciferase from Renilla. Protein Sci. 2008;17(4):725–35. Epub 2008/03/25. pmid:18359861; PubMed Central PMCID: PMC2271170.
    • 58. Conkright MD, Guzman E, Flechner L, Su AI, Hogenesch JB, Montminy M. Genome-wide analysis of CREB target genes reveals a core promoter requirement for cAMP responsiveness. Mol Cell. 2003;11(4):1101–8. Epub 2003/04/30. pmid:12718894.
    • 59. Zhao L, Li G, Zhou GQ. SOX9 directly binds CREB as a novel synergism with the PKA pathway in BMP-2-induced osteochondrogenic differentiation. Journal of bone and mineral research: the official journal of the American Society for Bone and Mineral Research. 2009;24(5):826–36. Epub 2008/12/31.
    • 60. Long F, Schipani E, Asahara H, Kronenberg H, Montminy M. The CREB family of activators is required for endochondral bone development. Development (Cambridge, England). 2001;128(4):541–50. Epub 2001/02/15.
    • 61. Thavathiru E, Das GM. Activation of pRL-TK by 12S E1A oncoprotein: drawbacks of using an internal reference reporter in transcription assays. Biotechniques. 2001;31(3):528–30, 32. pmid:11570496.
    • 62. Ibrahim NM, Marinovic AC, Price SR, Young LG, Frohlich O. Pitfall of an internal control plasmid: response of Renilla luciferase (pRL-TK) plasmid to dihydrotestosterone and dexamethasone. Biotechniques. 2000;29(4):782–4. pmid:11056808.
    • 63. Quinn PG. Distinct activation domains within cAMP response element-binding protein (CREB) mediate basal and cAMP-stimulated transcription. The Journal of biological chemistry. 1993;268(23):16999–7009. Epub 1993/08/15. pmid:8394325
    • 64. Dougherty DC, Sanders MM. Comparison of the responsiveness of the pGL3 and pGL4 luciferase reporter vectors to steroid hormones. Biotechniques. 2005;39(2):203–7. pmid:16116793.
    • 65. Herzig S, Hedrick S, Morantte I, Koo SH, Galimi F, Montminy M. CREB controls hepatic lipid metabolism through nuclear hormone receptor PPAR-gamma. Nature. 2003;426(6963):190–3. pmid:14614508.
    • 66. Kim CH, Xiong WC, Mei L. Inhibition of MuSK expression by CREB interacting with a CRE-like element and MyoD. Molecular and cellular biology. 2005;25(13):5329–38. pmid:15964791; PubMed Central PMCID: PMC1156998.
    • 67. Melnikova VO, Dobroff AS, Zigler M, Villares GJ, Braeuer RR, Wang H, et al. CREB inhibits AP-2alpha expression to regulate the malignant phenotype of melanoma. PloS one. 2010;5(8):e12452. pmid:20805990; PubMed Central PMCID: PMC2929203.
    • 68. Chatterjee S, Pal JK. Role of 5'- and 3'-untranslated regions of mRNAs in human diseases. Biology of the cell / under the auspices of the European Cell Biology Organization. 2009;101(5):251–62. Epub 2009/03/12.
    • 69. Yang B, Guo H, Zhang Y, Chen L, Ying D, Dong S. MicroRNA–145 regulates chondrogenic differentiation of mesenchymal stem cells by targeting Sox9. PloS one. 2011;6(7):e21679. Epub 2011/07/30. PubMed Central PMCID: PMCPmc3140487.
    • 70. Martinez-Sanchez A, Dudek KA, Murphy CL. Regulation of human chondrocyte function through direct inhibition of cartilage master regulator SOX9 by microRNA–145 (miRNA–145). The Journal of biological chemistry. 2012;287(2):916–24. Epub 2011/11/22. PubMed Central PMCID: PMCPmc3256897.
    • 71. Kozak M. Structural features in eukaryotic mRNAs that modulate the initiation of translation. The Journal of biological chemistry. 1991;266(30):19867–70. Epub 1991/10/25. pmid:1939050