Yu-Mei TIAN, Jie CHEN, Yang TAO, Xia-Yun JIANG, Shu-Ming ZOU*



Molecular cloning and function analysis of insulin-like growth factor-binding protein 1a in blunt snout bream ()

Yu-Mei TIAN, Jie CHEN, Yang TAO, Xia-Yun JIANG, Shu-Ming ZOU

Key Laboratory of Freshwater Aquatic Genetic Resources, Shanghai Ocean University, Shanghai 201306, China

Insulin-like growth factor-binding protein 1 (IGFBP-1), a hypoxia-induced protein, is a member of the IGFBP family that regulates vertebrate growth and development. In this study, full-length IGFBP-1a cDNA was cloned from a hypoxia-sensitive Cyprinidae fish species, the blunt snout bream (). IGFBP-1a was expressed in various organs of adult blunt snout bream, including strongly in the liver and weakly in the gonads. Under hypoxia, IGFBP-1a mRNA levels increased sharply in the skin, liver, kidney, spleen, intestine and heart tissues of juvenile blunt snout bream, but recovered to normal levels after 24-hour exposure to normal dissolved oxygen. In blunt snout bream embryos, IGFBP-1a mRNA was expressed at very low levels at both four and eight hours post-fertilization, and strongly at later stages. Embryonic growth and development rates decreased significantly in embryos injected with IGFBP-1a mRNA. The average body length of IGFBP-1aoverexpressed embryos was 82.4% of that of the control group, and somite numbers decreased to 85.2%. These findings suggest that hypoxia-induced IGFBP-1a may inhibit growth in this species under hypoxic conditions.

; IGFBP-1a; Hypoxia; Overexpression

The growth of vertebrates is primarily regulated by the insulin-like growth factors (IGF) signaling pathway (Baker et al, 1993; Duan & Xu, 2005; Reindl & Sheridan, 2012). Insulin-like growth factor-binding proteins (IGFBPs), with molecular weights of 24000–50000, bind IGFs and regulate their binding with corresponding receptors; thus, they regulate the biological functions of the IGF signaling pathway (Hwa & Rosenfeld, 1999; Firth & Baxter, 2002; Shimizu et al, 2005; Shen et al, 2012).All six IGFBP genes are known to exist in teleosts and have been characterized individually or as a group, including their phylogenetic relationships and the impact of genome duplication (Macqueen et al, 2013).

Currently, duplicate, i.e.and, have been found in zebrafish, Nile tilapia, fugu and medaka (Kamei et al, 2008). IGFBP-1 is the only one of the six IGFBP proteins that is regulated by hypoxia. Under hypoxia, IGFBP-1 expression increases, which suppresses the IGF signaling pathway (Kajimura et al, 2006; Kajimura & Duan, 2007). Hypoxia raises the expression levels of IGFBP-1 mRNA in hypoxia-tolerant grass carp embryos, which suppresses the growth and development rates of grass carp embryos (Sun et al, 2011). Here, we describe the molecular characterization and function analysis of IGFBP-1a in blunt shout bream (), a hypoxia-sensitive and aquaculturally important species (Shen et al, 2010).

MATERIALS AND METHODS

Experimental fish

Blunt snout bream () adults and fertilized eggs were sampled from the Gen­etics and Breeding Center for Blunt Snout Bream of Shanghai Ocean University. Three ~1 000 g adult fish were sacrificed by immersion in MS-222 (tricaine methanesulfonate; Sigma, St. Louis, MO, USA) and their spleens, muscles, kidneys, eyes, hearts, gills, livers, intest­ines, brains, skins and ovaries or testes were remo­ved, cryopreserved in liquid nitrogen, and then reserved at-80 °C. Fertilized eggs were obtained from artificial fertilization and hatched in embryo rearing solution (ERS) replaced every four hours. Meanwhile, 20 viable embr­yos were chosen, observed and photographed under a stereoscopic microscope (SMZ1500; Nikon, Japan) and then reserved in RNAstore preservation solution (TIAN­GEN BIOTECH, Beijing, China) at 4 °C.

Hypoxia treatments

Totally, 30 healthy juvenile blunt snout bream fishwith weights of approximately 150–180 g were selected. After cultivation for one week in an automatic water cycle system at 20 °C, they were divided randomly into two experimental groups and a control group. The experimental facility was designed following Zhang et al (2003). With nitrogen bubbling continuously, the oxygen dissolved in the water was kept at 1.0±0.5 mg/L, a level of acute hypoxia. The control group was placed in an identical environment, that only differed in the dissolved oxygen level of 7.0±0.5 mg/L.

After experiencing hypoxia for four hours, five blunt snout bream from each group were chosen for dissection. The dissolved oxygen value for the anoxic groups was sequentially recovered over the course 60 min to a similar condition of the control group (7.0±0.5 mg/L). After 24 h later, five blunt snout bream from each oxygen-recovered group and the control group were chosen for dissection. Skin, muscles, eyes, gills, brains, intestines, spleens, livers, kidneys and hearts were extracted and then quick-frozen in liquid nitrogen and then preserved in a refrigerator at-80 °C.

Cloning of blunt snout bream IGFBP-1a cDNA

Total RNA was extracted from embryos 28 hours post-fertilization (hpf) using the TRIzol (Invitrogen) method, and randomly reverse transcribed into comp­le­mentary DNA (cDNA) using M-MLV reverse transcri­ptase (TaKaRa, Japan) following the manufac­turer’s procedures. Based on the conserved region of IGFBP-1a mRNAof zebrafish (GenBank Accession No. NM_173­283.3) and grass carp (Sun et al, 2011), a pair of specific primers (IGFBP-1a-F, 5′-GCTGCCTCGCCTGTGCGTTGAAG-3′ and IGFBP-1a-R, 5′-TCCAGGATGACACACACCAACAC-3′) were designed for blunt snout bream IGFBP-1a cDNA cloning, and reverse transcri­ptase polymerase chain reactions (RT-PCR) were perfor­med. RACE and nested PCR procedures were performed using a SMART RACE cDNA Amplification Kit (Clont­ech, TaKaRa, Japan), while the primers (3′-IGFBP-1a-GSP3r, 5′-GCATGAAATCCAAAGTCAACGCAATACG-3′, 3′-IGFBP-1a-GSP­3n,5′-AAA GTCAAACAGTGTGAATCGTCTC-3′, 5′-IGFBP-1a-GSP5r, 5′-GATTCACACTGTTTGACTTTGTAT-3′ and 5′- IGFBP-1a-GSP5n, 5′-GTGATTTAGTGATCTTGTCAAGGGC-3′; all synthes­ized by Sangon Biotech, Shanghai, China) were designed accor­ding to the kit’s instructions. Detected by 1.5% agarose gel electrophoresis and purified using TIAN gel Midi Purification Kit (Tiangen Biotech, Beijing, China), the nested PCR product was ligated into the pMD19-T (TaKaRa, Japan) vector system.

Sequence analysis

Alignments of the 5-RACE, 3-RACE and interm­ediate region were performed using ClustalW within BioEdit 7.0. Primers at both ends were cut off and full-length IGFBP-1a cDNA of the blunt snout brea­m was assembled. Similarity searches using IGFBP-1a cDNA sequences and predicted amino acids using the BLAST n/p program were performed against the GenB­ank database, while multiple alignments of the homo­logous IGFBP-1 nucleotide and amino acid sequences of the vertebrates from the search were performed using the CLUSTAL W program. A neighbor-joining phylogenetic tree was constructed using MEGA 5.0 (Tamura et al, 2011) using bootstrap values based on 1000 replications.

RT-PCR analysis

One microgram of DNAzyme digested total RNA from the tissues and different embryo developmental stages was reverse transcribed into cDNA using M-MLV reverse transcriptase (TaKaRa). The IGFBP-1aprimers (IGFBP-1a-F and IGFBP-1a-R) and β-actin gene primers (β-actin RT-F, 5′-CCGCTGCCTCTTCTTCCTC-3′ and β-actin RT-R, 5′-CTACCTCCCTTTGCCA GTTTCCGC-3′) that crossed the introns were designed based on the β-actin full-length cDNA sequence of blunt snout breamand the genomic sequence of zebrafish. RT-PCR was performed using different tissues from adult and embr­yonic blunt snout breamusing the above primer sets. The PCR conditions were as follows: an initial denaturation of 5 min at 94 °C, 30 cycles of 30 sec denaturation at 94 °C and 30 sec annealing at 55 °C followed by 72 °C for 30 sec and a final extension at 72°C for 10 min. The PCR products were detected using 1.2% agarose gel electro­phoresis and then photographed.

qRT-PCR analysis

Total RNA was extracted from tissues of the hypo­xia groups and the control group using TRIzol (Invitr­ogen) and reversed transcribed into cDNA using Oligo(dT)18 primer and M-MLV. Quantitative real-time PCR were performed in an iCycler iQ Multicolor using the cDNA template and the IGFBP-1aprimers (IGFBP-1a-F and IGFBP-1a-R) and β-actin gene primers (β-actin RT-F and β-actin RT-R). Quantitative analysis was performed using the 2method (Shen et al, 2010).

synthesis and microinjection of IGFBP-1a capping mRNA

The IGFBP-1a open reading frame (ORF) was sub-cloned into pCS2+ vector to construct the pCS2-IGFBP-1a. Capped IGFBP-1a mRNA of the blunt snout bream as well as GFP mRNA were synthesizedusing a commercial kit (mMESSAGE mMACHINE Kit; Ambion, TX, USA) with linearized plasmid DNA as the template.

TheIGFBP-1a mRNA (at a concentration of 500 pg/embryo) and EGFP mRNA (a concentration of 200 pg/embryo) were together injected into fertilized eggs at the 1–2-cells stage of blunt snout bream. The injection volume was 1–2 nL as previously described (Jiang et al, 2012). After injection, fertilized eggs were incubated in ERS at 25°C. The embryo injected with EGFP mRNA was used as a control. Growth status was observed under a Nikon SMZ1500 fluorescence microsc­ope.Phenotypes of injected eggs was observed every 4 hours.

Statistical analysis

Data from qRT-PCR, somite number and body length are expressed as mean±). Differences among groups were analyzed with a one-way ANOVA followed by Fisher’s post hoc tests or unpaired t-tests. Significance was accepted at the level of<0.01 or 0.001.

RESULTS

Cloning of full-length cDNA of blunt snout bream IGFBP-1a

The full length of the cDNA of blunt snout bream IGFBP-1a is 1 178 bp, containing 116 bp of 5′-UTR and 273 bp of 3′-UTR, which contains a microsatellite sequence of 49 CA repeats at the 913–1 030 bp region (Figure 1). The Cys residues of the protein are the same as the IGFBPs of other vertebrates in both number and position. The N-terminal of the protein contains a repres­entative IGFBP sequence motif (GCGCCXXC) and the C-terminal contains a thyroglobulin-1 sequence motif (CWCV), whereas the intermediate region is less conser­ved. Clustering analysis confirmed that duplicateds(and) exist in teleosts (Figure 2), and thegene we cloned fell into theclade.

IGFBP-1a expression features in tissues and expres­sion regulation under hypoxia

IGFBP-1a mRNAwas expressed in most tissues of adult blunt snout bream, and the transcription was highest in the liver and very weak in the testis and ovary (Figure 3A). qRT-PCR analysis of IGFBP-1a mRNA in each tissue type was performed after hypoxia and oxygen recovery (Figure 3B). The IGFBP-1a transcription level in the skin, livers, kidneys, spleens, intestines and hearts of juvenile blunt snout bream increased sharply (18.5, 19.0, 8.8, 22.1, 16.9 and 29.5 times those of the control group, respec­tively). After 24-hour exposure to normal dissolved oxy­gen, transcription in the tissues returned to normal levels. The IGFBP-1a transcriptional levels in the eyes, muscles, brains and gills were not influenced by hypoxia.

IGFBP-1a overexpression significantly reduces blunt snout bream embryonic development and growth rate

Blunt snout bream IGFBP-1a mRNA was expressed throughout the embryonic development process. The expression level was very low at the 4 and 8 hpf embr­yonic development stages, increased significantly at 12 and 16 hpf and remained stable after 20 hpf (Figure 4A).

To further determine the function of blunt snout bream IGFBP-1a, capped IGFBP-1a mRNA and EGFP mRNA were injected together into fertilized eggs at the blunt snout bream 1–2-cell stage. Capped EGFP mRNA was injected as a control. Injection success was determined by detecting the expression ofEGFP. The 24 hpf embryos injected with EGFP mRNA did not significantly differ from those of wild type embryos in terms of morphology, somite number and body length (Figure 4B, 4D). Compared with the wild type or EGFP mRNA injected group (Figure 4B, 4D, 4F), embryos injected with IGFBP-1a mRNA showed a decrease (<0.01) in body length and somite number, 30% (18/60) of which were malformed to various degrees (Figure 4C, 4E, 4G). Blunt snout bream embryos in the IGFBP-1a mRNA injected group grew more slowly (<0.01) than the wild type or EGFP mRNA injected group (Figure 4H). At 24 hpf, the average body length of the IGFBP-1a injected embryos was 2.8±0.3 mm, 82.4% of that of the EGFP mRNA injected group (3.4±0.1 mm). Additionally, 24.9±2.6 somites was found in embryos injected with IGFBP-1a mRNA at 24 hpf, which is only 85.2% of the EGFP mRNA injected control (29.2±1.6) (Figure 4I), equivalent to embryos at 17–20 hpf of the control group. These results indicate that increasing IGFBP-1a mRNA level in blunt snout bream embryos significantly decreased their growth and development rates.

Figure 1 cDNA and encoding amino acid sequences in blunt snout bream

The microsatellite sequence containing CA repeats is shown in the block. The PPD domain is shown in the dash block.

Figure 2 Clustering analysis of vertebrateamino acid sequences constructed using the neighbor-joining method

Accession numbers of sequences retrieved from GenBank are shown. Similarity searches and phylogenetic analysis revealed that the IGFBP-1a of blunt snout bream is most closely related to that of grass carp, common carp and zebrafish, with identities of 97%, 93% and 92%, respectively. The IGFBP-1a of blunt snout bream shows only about 60% identities with that of channel catfish (), Atlantic salmon (), rainbow trout () and Japanese yellowtail (), and about 50% with that of mouse () and human ().

Figure 3 Expression of blunt snout breamIGFBP-1a mRNA in adult tissues and the influence of the hypoxia treatment on juvenile blunt snout bream

A: RT-PCR analysis of blunt snout breamIGFBP-1a in adult tissues; CTR, negative control without genomic DNA; B: qRT-PCR analysis of IGFBP-1a mRNA at different tissues of juvenile blunt snout bream during normoxic (N: 7.0±0.5 mg O/L), hypoxic (H: 1.0±0.5 mg O/L) and recovery (R) conditions. IGFBP-1a mRNA copy number was normalized as a ratio to the β-actin (a house keeping gene) mRNA copy number. qRT-PCR data are presented as mean±SE; ***<0.001.

DISCUSSION

We obtained the full-length IGFBP-1acDNA in hypoxia-sensitive blunt snout bream, an aquaculturally important cyprinid fish. The 3′-UTR of the IGFBP-1a cDNA in blunt snout bream contains a microsatellite sequence of 49 CA repeats, which could be only found on IGFBP-1a in grass carp (39 CA) and IGFBP-1b in Japanese yellowtail (79 CA) and zebrafish (28 CA). Neither the cause or function of the microsatellite sequence of the CA repeats have been elucidated. Cluste­ring analysis revealed that duplicateds(and) exist in teleost fish and should result from fish-specific 3whole genome duplication. In addition, duplicatedsexists in the genome of pseudo-tetraploid Atlantic salmon and rainbow trout, which are considered pseudo-tetraploid fish species (Macqueen et al, 2013). The duplicatedgenes were supposed to be the result of the 3rd or 4th genome duplication specific to fishes (Taylor et al, 2003; Jaillon et al, 2004; Crow et al, 2006; Ocampo Daza et al, 2011).

The expression of IGFBP-1a mRNA was detected in a wide range of tissues in adult blunt snout bream. Under hypoxia, the expression levels of IGFBP-1a mRNA increased greatly in many tissues in juveniles, and recovered to normal levels after 24-hour exposure to normal dissolved oxygen. Similar to our result, average expression levels in the liver of adult zebrafish under hypoxia is approximately 280 times that of the control group (Maures & Duan, 2002, Kajimura et al, 2005). Studies of grass carp and Atlantic croaker also show increases in IGFBP-1a mRNA expression levels under hypoxia (Sun et al, 2011; Rahman & Thomas, 2011).

In blunt snout bream embryos, IGFBP-1a is highly expressed at 12 hpf and only weakly at 4 and 8 hpf. This is different to hypoxia-tolerant grass cap embryos, in which IGFBP-1a mRNA expression was not detected at 4 or 8 hpf, and only trace amounts at 12 hpf (Tao & Zou, 2011). In zebrafish,transcription of the duplication IGFBP-1b gene could not be detected until the later embryo developm­ental stage, but IGFBP-1a is present at all stages(Kamei et al, 2008), similar to the expression cha­ra­c­te­ristic of blunt snout bream. Fish eggs carrymaternal mRNA and embryonic transcription of its own mRNA dates back to the 9thcleavage (Kane & Kimmel, 1993). The IGFBP-1a mRNA detected at early stages is likely from the maternal deposit and the significant increase in mRNA during the middle and later stages results from embryonic transcription.

Figure 4 Overexpression of IGFBP-1a reduced growth and developmental rate in blunt snout bream embryos

A: RT-PCR analysis of IGFBP-1a mRNA in blunt snout bream during embryogenesis. Developmental stages are indicated above the figure. hpf, hour post-fertilization. CTR, negative control without genomic DNA; B–G: Blunt snout bream embryos are shown at the 24 hpf stage by visible light (B,C,F,G) and fluorescence (D,E). EGFP mRNA (200 pg) (B, D and F) or a mixture of EGFP mRNA (200 pg) with IGFBP-1a mRNA (500 pg) (C, E and G) was injected into each embryo at the 1–2 cell stage. (H, I) The body length (H) and somite number (I) of IGFBP-1a and/or EGFP mRNA overexpressed blunt snout bream embryos at24 hpf.The results are presented as mean ± SE; **:<0.01. Scale bar=600 μm.

Microinjection of capped IGFBP-1a mRNA caused significant developmental retardation and growth reduc­tions in blunt snout bream embryos. It will be interesting to test whether hypoxia induces IGFBP-1a expression in blunt snout bream embryos in the future. Similar to the present results, hypoxia-induced growth reduction was also caused by high IGFBP-1a expression in grass cap embryos (Sun et al, 2011). In zebrafish, hypoxia stimulates increased IGFBP-1a transcription during early embryonic develo­pment and causes embryo developm­ental retardation and growth reductions, whereas only IGFBP-1b is affected by hypoxia during later stages (Kamei et al, 2008). Acc­ordingly, IGFBP-1a likely functions as a bridge between hypoxia and the inhibition of embryonic growth in fish.

Baker J, Liu JP, Robertson EJ, Efstratiadis A. 1993. Role of insulin-like growth factors in embryonic and postnatal growth.75(1): 73-82.

Crow KD, Stadler PF, Lynch V, Amemiya C, Wagner GP. 2006. The ‘‘fishspecific’’cluster duplication is coincident with the origin of teleosts.23(1): 121-136.

Duan C, Xu Q. 2005. Roles of insulin-like growth factor (IGF) binding proteins in regulating IGF actions.142(1-2): 44-52.

Firth SM, Baxter RC. 2002. Cellular actions of the insulin-like growth factor binding proteins.23(6): 824-854.

Hwa V, Oh Y, Rosenfeld RG. 1999. The Insulin-like growth factor-binding proteins (IGFBPs) super family.20(6): 761-787.

Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, et al. 2004. Genome duplication in the teleost fishreveals the early vertebrate proto-karyotype.431(7011): 946-957.

Jiang XY, Du XD, Tian YM, Shen RJ, Sun CF, Zou SM. 2012. Goldfish transposasepresumably from recent horizontal transfer is active.26(7): 2743-2752.

Kajimura S, Aida K, Duan C. 2005. Insulin-like growth factor-binding protein-1 (IGFBP-1) mediates hypoxia-induced embryonic growth and developmental retardation.102(4): 1240-1245.

Kajimura S, Aida K, Duan C. 2006. Understanding hypoxia-induced gene expression in early development: in vitro and in vivo analysis of hypoxia inducible factor 1-regulated zebrafish insulin-like growth factor binding protein 1 gene expression.26(3): 1142-1155.

Kajimura S, Duan C. 2007. Insulin-like growth factor (IGF) binding protein-1: An evolutionarily conserved fine tuner of IGF actions under catabolic conditions.71: 309-325.

Kamei H, Lu L, Jiao S, Li Y, Gyrup C, Laursen LS, Duan C. 2008. Duplication and diversification of the hypoxia-inducible IGFBP-1 gene in zebrafish.3: e3091.

Kane DA, Kimmel CB. 1993. The zebrafish midblastula transition.119(2): 447-456.

Macqueen DJ, Garcia de la Serrana D, Johnston IA. 2013. Evolution of ancient functions in the vertebrate insulin-like growth factor system uncovered by study of duplicated salmonid fish genomes.30(5): 1060-1076.

Maures TJ, Duan C. 2002. Structure, developmental expression, and physiological regulation of zebrafish IGF binding protein-1.143(7): 2722-2731.

Ocampo Daza DO, Sundstrom G, Bergqvist CA, Duan C, Larhammar D. 2011. Evolution of the insulin-like growth factor binding protein (IGFBP) family.152(6): 2278-2289.

Rahman MS, Thomas P. 2011. Characterization of three IGFBP mRNAs in Atlantic croaker and their regulation during hypoxic stress: potential mechanisms of their upregulation by hypoxia.301(4): E637-E648.

Reindl KM, Sheridan MA. 2012. Peripheral regulation of the growth hormone-insulin-like growth factor system in fish and other vertebrates.163(3-4): 231-245.

Shen RJ, Jiang XY, Pu JW, Zou SM. 2010. HIF-1α and -2α genes in a hypoxia-sensitive teleost species: cDNA cloning, expression and different responses to hypoxia.157(3): 273-280.

Shen WY, Ren G, Zhu YR. 2012. Effects of compensatory growth on the levels of IGF-1, IGFBP-1 and expressions of IGF-1 mRNA, IGF-1R mRNA in.33(3): 298-303. (in Chinese)

Shimizu M, Dickey JT, Fukada H, Dickhoff WW. 2005. Salmon serum 22 kDa insulin-like growth factor-binding protein(IGFBP) is IGFBP-1.184(1): 267-276.

Sun CF, Tao Y, Jiang XY, Zou SM. 2011. IGF binding protein 1 is correlated with hypoxia-induced growth reduce and developmental defects in grass carp () embryos.172(3): 409-415.

Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.28(10): 2731-2739.

Tao Y, Zou SM. 2011. cDNA cloning and expression of insulin-like growth factor binding protein gene IGFBP-1 in.20(1): 11-19. (in Chinese)

Taylor JS, Braasch I, Frickey T, Meyer A, Van de Peer Y. 2003. Genome duplication, a trait shared by 22000 species of ray-finned fish.13(3): 382-390.

Zhang Z, Wu RSS, Mok HOL, Wang Y, Poon WWL, Cheng SH, Kong RYC. 2003. Isolation, characterization and expression analysis of a hypoxia-responsive glucose transporter gene from the grass carp,.270(14): 3010-3017.

29 September 2013; Accepted: 26 December 2013

This work was supported by grants from the Key Technologies Research and Development Program of China (2012B­AD26B02), National Science Foundation of China (31272633; 31201­760), Special Funding from the Ministry of Agriculture of China (2012­0­3­081, 201203086), and Shanghai University Knowledge Service Platform (ZF1206).

, E-mail:smzou@shou.edu.cn