Poult. Sci.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Poult Sci 2006. 85:1746-1754
© 2006 Poultry Science Association
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zhan, K.
Right arrow Articles by Yang, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zhan, K.
Right arrow Articles by Yang, N.

MOLECULAR, CELLULAR, AND DEVELOPMENTAL BIOLOGY

Molecular Cloning and Expression of the Duplicated Thyroid Hormone Responsive Spot 14 (THRSP) Genes in Ducks1

K. Zhan, Z. C. Hou, H. F. Li, G. Y. Xu, R. Zhao and N. Yang2

Department of Animal Genetics and Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100094

2 Corresponding author: nyang{at}cau.edu.cn


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Thyroid hormone responsive Spot 14 (THRSP) is suggested as a transcription factor involved in the regulation of adipogenic enzymes by 3 thyroid response elements in the promoter region. In the chicken genome, THRSP gene was identified to duplicate into 2 paralogs, THRSP{alpha} and THRSPß. In the current study, cDNA sequences of the duplicated duck THRSP genes were cloned by real-time PCR and rapid amplification of cDNA ends. Duck THRSP{alpha} and THRSPß were predicted to encode peptides with 133 amino acids, which had 74 and 68% sequence identity at cDNA level, 78 and 74% identity at amino acid level to the chicken counterparts, respectively. A high percentage (73.1%) of G and C nucleotides were found in the 3' untranslated region of duck THRSPß cDNA. Although a low similarity of peptide composition was shared between ducks and mammals, and a moderate similarity was shared between ducks and chickens, many predicted properties of THRSP, including the pI, subcellular localization and functional domains seemed to be highly conserved. The present study demonstrated that the duck THRSP gene duplicates into the 2 paralogs as in chickens. Phylogenetic analysis indicated that the duplication for THRSP paralogs appeared to have taken place preceding the chicken–duck split, and the diverging rate between THRSP paralogs seemed faster in the chicken genome than that in the duck genome. Expression analysis by real-time quantitative PCR showed that THRSP paralogs in ducks were more actively transcribed in fat tissues (i.e., s.c. fat and abdominal fat) than in liver, and the mRNA concentrations of THRSPß were higher than that of THRSP{alpha} in liver and s.c. fat.

Key Words: thyroid hormone responsive Spot 14 • paralog • complementary deoxyribonucleic acid • expression • duck


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
The thyroid hormone responsive Spot 14 (THRSP) gene, also known as Spot 14, encodes a small acidic protein expressed predominately in the adipogenic tissues, such as the lactating mammary gland, fat, and the liver (Seelig et al., 1981; Jump et al., 1984; Jump and Oppenheimer, 1985). In humans, THRSP was shown to be associated with obesity (Chagnon et al., 1998), growth, and differentiation of breast cancer cells (Moncur et al., 1998; Sanchez-Rodriguez et al., 2005). As a transcription factor with 3 thyroid response elements in the promoter region, THRSP is postulated to govern the expression of a cascade of enzymes in the adipogenic pathway of rat liver by responding to the stimulation of triiodothyronine, as well as other factors (i.e., carbohydrate, polyunsaturated fatty acids, insulin, and glucagon) related to hepatic lipid metabolism (Jump et al., 1993; Liu and Towle, 1994).

The chicken THRSP gene was found by microarray analysis as a differentially regulated EST in livers of chickens divergently selected for fast or slow growth rates (Cogburn et al., 2000), and its expression was regulated by thyroid hormone status (Wang et al., 2002). The gene was mapped to a chromosomal region (1q41–44; Carre et al., 2001) harboring QTL for s.c. fatness (Ikeobi et al., 2002) and abdominal fatness (Lagarrigue et al., 2003). Recently, the chicken THRSP gene was cloned by in silico EST assembling and was identified to duplicate into 2 paralogs, THRSP{alpha} and THRSPß (Wang et al., 2004). Furthermore, Cogburn et al. (2003) and Wang et al. (2004) observed that the polymorphisms in the coding regions of the paralogous genes were associated with the abdominal fatness of chicken. Sequence analysis of mammal and chicken THRSP genes demonstrated that both of them shared a similar gene organization, with 2 exons and an intron, but the deduced chicken THRSP{alpha} peptide had a low similarity (less than 30% of identities and 50% of positives) to the THRSP amino acid sequences of mammals and fishes (Grillasca et al., 1997; Wang et al., 2004). In addition, the ortholog of chicken THRSPß hasn’t been identified in other vertebrates.

Ducks and chickens belong to sister avian orders, and both of them share de novo fatty acid biosynthesis, mainly in the liver. Although comparative fluorescence in situ hybridization mapping suggests poor conservation between the chicken and duck genes (Yuan et al., 2005), it is of interest to explore if the gene organization and duplication episode of the chicken THRSP are shared in the duck genome. In addition, THRSP may play an important role in avian adipogenesis, and it is necessary to get cDNA sequence and expression profiles of duck THRSP for future functional genomic investigations. To this end, we carried out the present study to identify and characterize the duplicated THRSP genes in the duck genome.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Tissue Sampling and Total RNA Isolation
Six Pekin ducks (3 males and 3 females) at 6 wk of age were obtained from Jinxing Pekin Duck Breeding Center in Beijing, China. The birds were killed by decapitation for the collection of the blood and the tissues of interest, including abdominal fat, s.c. fat, liver, uropygial gland, muscle, heart, pancreas, kidney, brain, pituitary, thymus, spleen, testes, and ovary. The China Agricultural University Animal Care and Use Committee approved the protocol for bird treatments. The tissues were collected and snap-frozen in liquid N and stored at –80°C until the extraction of RNA. Total RNA was prepared using Trizol LS reagent (Gibco BRL, Gaithersburg, MD) according to the manufacturer’s guidelines. Genomic DNA was isolated from blood samples by the phenol–chloroform method. Ribonucleic concentration and purity were quantified in the samples by densitometry at 260 and 280 nm, and the integrity was examined by electrophoresis in formaldehyde-agarose gels. Finally, each RNA was diluted to 1 µg/µl with RNase-free water and stored at –80°C.

Primer Design
As low sequence similarities in both nucleotides and amino acids are shared by chickens and mammals, primers were designed only based on the available sequences of chicken THRSP{alpha} and THRSPß cDNA, with GenBank Accession no. AY568628 [GenBank] for THRSP{alpha} and AY568630 [GenBank] for THRSPß (Wang et al., 2004). The gene-specific primers to clone the 3' ends of THRSP{alpha} and THRSPß cDNA using the rapid amplification of cDNA ends (RACE) procedure were first designed from the nucleotides encoding the conservative domains of the chicken THRSP peptides. The primers for 5'-RACE of THRSP{alpha} or amplifying open reading frame (ORF) of THRSPß were designed from the nucleotides sequenced from the 3'-RACE products and those encoding conservative N-terminus of THRSPß in chickens. The primers amplifying the full-length cDNA of THRSP paralogs and intron sequence of THRSP{alpha} were designed according to the assembled sequences of RACE products. Primers for real-time PCR were designed using Primer Express 2.0 software (Applied Biosystems, Foster City, CA). All primers were synthesized with ABI 3900 DNA/RNA Synthesizer (Applied Biosystems), and they are listed in Table 1Go.


View this table:
[in this window]
[in a new window]
 
Table 1. Primers for the cloning and expression of duck thyroid hormone responsive Spot 14 (THRSP) genes
 
Cloning of Duck THRSP cDNA
Attempts to obtain the full length of duck THRSP mRNA sequences were carried out by the technique of rapid amplification of cDNA ends; 3'- and 5'-RACE kits (Invitrogen, Carlsbad, CA) were used according to the procedures provided by the manufacturer. For 3'-RACE, reverse transcription was conducted to synthesize RACE-ready cDNA using Moloney murine leukemia virus reverse transcriptase (SuperScript II reverse transcription, Invitrogen). The 3' adapter primer (500 nmol/L) and the duck liver RNA (2 µg) were heated to 70°C for 10 min. Then, 200 U of SuperScript II reverse transcription was added to each reaction, and the reactions were incubated for 50 min at 42°C. Two microliters of RACE-ready cDNA product in a 50-µL volume were amplified with the sense primers ({alpha} 3 and ß 3) and abridge universal amplification primer. A touchdown PCR was used with a PTC-200 thermocycler (MJ Research, Watertown, MA), and temperature profile was set as follows: 1 cycle of 95°C for 5 min; 5 cycles of 95°C for 45 s, 70°C for 45 s (–2°C per cycle), and 72°C for 1 min; 20 cycles of 95°C for 45 s, 60°C for 45 s, and 72°C for 1 min; 1 cycle of 72°C for 10 min. Due to the high GC content in the 3' untranslated region (UTR) of chicken THRSPß, Takara LA Taq with GC buffer (Takara, Tokyo, Japan) was used in performing the hot start condition, with 98°C at the first denaturation step to clone the duck THRSPß 3'-UTR.

For 5'-RACE of THRSP{alpha}, 2 µg of liver RNA were primed by using the obtained cDNA sequence (primer {alpha} 5–1) and reverse transcribed as described above. After the cDNA synthesis, RNase H was used to degrade the remainder of the RNA template. The first-strand cDNA was purified by spin cartridge column. A homopolymeric C tail was then added to the 3' end of the cDNA, using deoxycytidine triphosphate (dCTP) and terminal deoxynucleotide transferase. The first-round PCR amplification was performed with 5' abridged anchor primer and the nested primer {alpha} 5–2. The obtained primary PCR product was diluted and reamplified using abridge universal amplification primer and the nested primer {alpha} 5–3. The thermal profile for the first-round PCR process consisted of an initial 94°C denaturation step for 2 min followed by 35 cycles of denaturation (94°C for 30 s), annealing (55°C for 50 s), and extension (72°C for 1 min). Finally, a 7-min 72°C step was used before holding at 4°C. The nested PCR amplification for the final 5'-RACE result was conducted as the touchdown protocol outlined above. A degenerated sense primer ß ICF For was designed to amplify the ORF and partial 3'-UTR of THRSPß with antisense primer ß ICF Rev together, and the initial annealing temperature was set at 59°C in touchdown PCR. The PCR products were separated by electrophoresis in 1.5% agarose, stained with ethidium bromide and visualized with a UV transilluminator gel documentation system (BioRad, Richmond, CA). The separated PCR products were purified with the QIAquick DNA purification kit (Qiagene, Hilden, Germany) and directly sequenced or ligated into the pMD18-T Vector (Takara). Multiple clones for each fragment were bidirectionally sequenced using Big Dye Terminator v 3.1 reagents on an Applied Biosystems 3730 DNA analyzer (Applied Biosystems) by the dideoxy-mediated chain-termination method (Sanger et al., 1977).

After the overlapping cDNA fragments were assembled according to the sequencing results, the full-length or complete cDNA sequences of THRSP{alpha} and THRSPß were con-firmed by amplifying the cDNA templates with the primer sets PS {alpha} For and PS {alpha} Rev and PS ß For and PS ß Rev, respectively. The primers of {alpha} I For and {alpha} I Rev were used to amplify the intron sequence of THRSP{alpha} in genomic level. Simultaneously, the duck THRSP{alpha} gene containing 2 exons and 1 intron was amplified using genomic DNA templates and the primers of PS {alpha} For and PS {alpha} Rev. A 50-µL volume of Takara LA Taq with GC buffer (Takara) was chosen to perform the similar touchdown PCR protocol depicted above with the different initial annealing temperatures displayed in Table 1Go. Following the consecutive steps of purification and clone sequencing for the PCR products, the integrity of cDNA of THRSP paralogs was confirmed by the alignment of the full-length nucleotides sequenced with the assembling nucleotides and exhibiting the gel picture of the PCR products.

Sequence and Phylogenetic Analysis
Nucleotide sequences of 3'- and 5'-RACE clones were assembled and identified by the NCBI BLAST search program, and the overall cDNA sequences were aligned using the DNAMAN software package (Lynnon Biosoft, Quebec, Canada). The amino acid similarity, including identities and positives, was analyzed using the BLASTp program with BLOSUM62 scoring matrix. The alignment of multiple peptides was created with the ClustalW multiple sequences alignment program and displayed using the BOX-SHADE 3.21 program (http://www.ch.embnet.org). After the ClustalW alignment using amino acid sequences of THRSP, the phylogenetic tree was constructed by MEGA software version 3.0,using maximum parsimony methods (Kumar et al., 2004). The p-distances (proportion of differences) with a high resolution of branching pattern were calculated (Nei and Kumar, 2000). The reliability of the constructed tree was tested by bootstrap analysis implemented in MEGA.

Expression Profiling by Real-Time PCR
The mRNA levels of THRSP paralogs were measured using the QuantiTect SYBR Green PCR Kit (Qiagen, Tokyo, Japan) and the ABI PRISM 7900 Sequence Detection System (Applied Biosystems). The quantities of THRSP{alpha} and THRSPß messages from different tissues were normalized with the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) mRNA to compensate for variations in input RNA amounts. Eight serial dilutions (from 1 x 102 to 1 x 109 copies per µl) of a plasmid from the GAPDH clone were selected as the templates to generate a standard curve. All reactions were carried out in a 15-µL volume containing 7.5 µL of SYBR Green Master Mix (Qiagen), 5.5 µL of water, 5 pmol of each primer, 1 µL of serial diluted GAPDH plasmids, or 1 µL of cDNA derived from reverse transcription of each tissue RNA using Moloney murine leukemia virus reverse transcriptase (Promega Corp., Madison, WI). Real-time PCR was performed to amplify each cDNA sample in duplicate, the negative control, and 8 diluted GAPDH plasmids within the same 96-well microplates. The amplification program recommended by the manufacturer was used. Specificity of the real-time PCR was evaluated by both high-resolution gel electrophoresis and melting-curve analysis. Normalization was conducted by dividing the average copies of each isoform of THRSP paralogs by the average copies of GAPDH in each tissue. The relative mRNA concentrations of THRSP paralogs were determined by transformation of common logarithm for the 1,000-fold normalized values and represented in arbitrary units.

Statistical Analysis
Relative mRNA levels of THRSP genes in the duck adipogenic tissues were subject to ANOVA. The expression of THRSP paralogs in same tissue was analyzed by Student’s t-test. Analysis of variance and Student’s t-test were performed using SAS 8.2 package (SAS Institute Inc., Cary, NC). The data were presented as means ± SEM.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Molecular Cloning of Duck THRSP Paralogs
Two fragments of 669 bp for THRSP{alpha} and 402 bp for THRSPß were generated by 3'-RACE, and 5'-RACE for THRSP{alpha} yielded a 282-bp fragment. The full-length THRSP{alpha} cDNA, with 817 nucleotides, was obtained by assembling the overlapping 3'-RACE and 5'-RACE sequences. The ORF for THRSP{alpha} cDNA comprises 402 nucleotides deducing a peptide with 133 amino acids (Figure 1Go, panel A).The THRSP{alpha} peptide was predicted to have a molecular weight of 14.7 kDa and an isoelectric point of 4.44 by using Compute pI/Mw tool (http://ca.expasy.org/tools/pi_tool.html). The ORF is flanked by a 5'-UTR, which is 21 nucleotides long, and a 3'-UTR, which is 394 nucleotides long. The integrity of the full-length THRSP{alpha} cDNA was confirmed by real-time PCR using the primers of PS {alpha} For and PS {alpha} Rev (Figure 2Go, lane 2) and sequencing. A pair of primers of {alpha} I For and {alpha} I Rev allowed amplification of a 1,008-bp fragment containing the intron sequence of duck THRSP{alpha} based on the overall organization of the chicken THRSP{alpha} gene (Figure 2Go, lane 4). As a result, a 739 bp of intron of duck THRSP{alpha} gene was obtained, corresponding to a homolog of 637 bp found in chicken genome by searching in Genome Browser Gateway (http://genome.ucsc.edu/cgi-bin/hgGateway?org=chicken). The boundaries between the exon and intron conform to the GT-AG rule (data not shown). A 1,556-bp fragment (Figure 2Go, lane 5), which was amplified using the primers of PS {alpha} For and PS {alpha} Rev in genomic DNA level, was found to be composed of a 817-bp transcript and a 739-bp intron of the duck THRSP{alpha} by sequencing and alignment.


Figure 1
View larger version (30K):
[in this window]
[in a new window]
 
Figure 1. Nucleotide and deduced amino acid sequences of the duck thyroid hormone responsive Spot 14 (THRSP) {alpha} (panel A) and ß (panel B). Nucleotides are numbered on the left side of the sequences. The bold underlined bases denote the primer sequences for rapid amplification of cDNA ends and PCR amplification. Amino acids are shown in signal letters under the cDNA sequences. Polyadenylation signals (AATAAA) and stop codons are indicated by short underlines and asterisks, respectively. The boundary between exon1 and exon2 in THRSP{alpha} is shown by the solid triangle. The alignment of peptide sequences of the duck THRSP{alpha} and THRSPß with the + sign representing positive amino acid substitutions (panel C). The cDNA sequences are available from the GenBank/EMBL/DDBJ DNA databases under Accession no.: DQ227766 and DQ227767 for THRSP{alpha} and THRSPß, respectively.

 

Figure 2
View larger version (78K):
[in this window]
[in a new window]
 
Figure 2. Gel picture of the fragments of thyroid hormone responsive Spot 14 (THRSP) paralogs amplified using complementary and genomic DNA templates: lane 1 = 100 bp of DNA ladder I; lane 2 = full-length cDNA of THRSP{alpha} (817 bases); lane 3 = complete cDNA of THRSPß (711 bases); lane 4 = fragment containing THRSP{alpha} intron (1,008 bases); lane 5 = overall THRSP{alpha} fragment consisting of 2 exons and 1 intron (1,556 bases); and lane 6 = marker DGL 2000.

 
The duck THRSPß isoform has a high GC content in the 3'-UTR (73.1%), confirmed by the 3'-RACE product, similar to the chicken ortholog (74.7%; Wang et al., 2004). A 446-bp fragment containing ORF and partial 3'-UTR of the duck THRSPß was obtained using the degenerated sense primer (ß ICF For) and the antisense primer (ß ICF Rev) targeting to the 3'-RACE product. Of the 711 nucleotides, which were sequenced from THRSPß cDNA products, 309 and 402 bases represent the 3'-UTR and ORF, which also predicts a peptide of 133 amino acids, respectively (Figure 1Go, panel B). The attempt to amplify the 5'-UTR of duck THRSPß cDNA by 5'-RACE assay failed to yield the predicted products due to high GC content and high similarity in ORF (91% identity) between THRSP paralogs. However, the 711-bp fragment of THRSPß cDNA, which was amplified using the primers of PS ß For and PS ß Rev (Figure 2Go, lane 3) and sequenced, demonstrated that the complete cDNA sequence of the duck THRSPß was as desired. Similar to the duck THRSP{alpha} isoform, the predicted THRSPß protein is acidic (pI of 4.53), with a molecular weight of 14.8 kDa. Without exception, the conserved AATAAA hexanucleotide polyadenylation signal (Proudfoot and Brownlee, 1976) was detected 18 nucleotides away from the poly(dA) tracts of both THRSP{alpha} and THRSPß cDNA.

The duck THRSP{alpha} and THRSPß were found to have 74 and 68% nucleotide sequence identities at the cDNA level as compared with the chicken orthologs, respectively. Higher similarities of nucleotide sequence (91%) and amino acid sequence (85% identities; 90% positives) in the putative coding region were detected between the duck THRSP paralogs (Figure 1Go, panel C) corresponding to those of nucleotide (80%) and amino acid (70% identities; 79% positives) between the chicken THRSP paralogs (Wang et al., 2004). Analysis with BLASTp revealed that the duck THRSP{alpha} peptide shared a moderate similarity (78% identities; 83% positives) to the chicken THRSP{alpha} protein and low similarity to that of humans (33% identities; 50% positives; Grillasca et al., 1997) and zebrafish (30% identities; 45% positives; zTC192887 found in TIGR: http://www.tigr.org/tdb/tgi/zgi). The duck THRSPß peptide is similar to the homologous amino acid sequences of the chicken (74% identities; 80% positives), the human (34% identities; 46% positives), and the zebrafish (33% identities; 45% positives). Both THRSP{alpha} and THRSPß peptides in ducks were predicted to localize in the nucleus and have a Leu zipper motif in the C-terminus by the PSORT II program (http://psort.ims.u-tokyo.ac.jp/), which is similar to the THRSP proteins of the chicken and other vertebrates.

The alignment of multiple amino acid sequences exhibited structural similarity of THRSP orthologs among the zebrafish, duck, chicken, bovine, human, mouse, and rat (Figure 3Go). The THRSP amino acid sequence is weakly conserved in vertebrates, as described above, but both isoforms of duck THRSP also share 3 conserved domains with the orthologs of the chicken and other vertebrates, which include a highly hydrophobic domain 1 near N-terminus, a second hydrophobic domain 2 in the middle, and a Leu zipper domain 3 in the carboxyl terminus. No difference of amino acids was found in the 2 hydrophobic regions of either THRSP isoform between duck and chicken. Leucine zipper motif, which confers DNA-binding ability of THRSP and conducts homodimerization of THRSP with the first hydrophobic domain (Cunningham et al., 1997), is more classical in ducks and chickens than in mammals and fish in terms of the number of Leu. As suggested by Compe et al. (2001), homodimers of THRSP interact with chicken ovalbumin upstream promoter-transcription factor 1 and upregulate the expression of L-type pyruvate kinase gene through a specificity protein 1 binding site localized in the proximal promoter.


Figure 3
View larger version (84K):
[in this window]
[in a new window]
 
Figure 3. Sequences alignment of orthologous thyroid hormone responsive Spot 14 (THRSP) proteins among different vertebrate species. Protein sequences for duck THRSP{alpha} (ABB16437) and THRSPß (ABB16438), chicken THRSP{alpha} (AAS77855) and THRSPß (AAS77857), bovine THRSP (AAT79485), human THRSP (AAH31989), rat THRSP (P04143), mouse THRSP (Q62264), and zebrafish THRSP (zTC192887) are aligned using ClustalW (http://www.ch.embnet.org) with default parameters and BLOSUM scoring matrix. Fraction of sequences is 0.5, in which identical amino acid residues are shown in black, and similar amino acid residues are shown in gray; the hyphens denote gaps. The amino acid residues in 3 frames indicate the putative domain 1, domain 2, and Leu zipper domain 3, respectively.

 
Duplication of THRSP Paralogs in the Duck Genome
A phylogenetic tree of the orthologous peptides for the vertebrate THRSP was constructed based on their homology in amino acid sequence using the zebrafish THRSP (zTC192887) as the outgroup (Figure 4Go). The duck THRSP paralogs were clustered into the same orthologous group with the chicken THRSP{alpha} and THRSPß. The duplication episode of THRSP occurred only in the genomes of both duck and chicken among different vertebrate species. Duplicated THRSP paralogs were not found in the mammalian genome when performing homologous searching by using the BLAST program in GenBank. No counterpart of THRSPß isoform was found in other vertebrate genomes by searching for the homologs in the protein database. Comparison of human genome sequence to that of chicken showed a similar number of older duplications in both species genomes, but there were significantly fewer paralogs specific to the chicken lineage than those specific to the human lineage (International Chicken Genome Sequencing Consortium, 2004). The overall stability of avian genomes is indicated by a low rate of gene duplication (Ellegren, 2005). Although many extant orders of birds date the divergence between chicken and duck back to 89.8 ± 6.97 million years (Tuinen and Hedges, 2001), Lynch and Conery (2000) estimated that gene duplication occurred and was fixed in populations at an average rate of about 1 per gene per 100 million years in eukaryotes. Thus, the genetic event generating the THRSP paralogs appeared to have taken place in the last common progenitor of both the duck and the chicken.


Figure 4
View larger version (6K):
[in this window]
[in a new window]
 
Figure 4. A dendrogram displays the phylogenetic relationship of orthologous thyroid hormone responsive Spot 14 (THRSP) proteins in vertebrates. The phylogenetic tree is constructed using the maximum parsimony method to calculate p-distances. The duck and the chicken have both THRSP{alpha} and THRSPß, whereas other species have only 1 of them. Bootstrap support values based on 1,000 replicates for the phylogenetic tree are shown the nodes above.

 
The phylogenetic analysis also suggests that THRSP{alpha} and THRSPß evolve independently in chickens and ducks. A larger genetic distance between duplicated THRSP paralogs in chickens (0.185) than that in ducks (0.111) was observed, implying a faster divergence for duplicated THRSP genes in the chicken genome than in the duck genome. Different functional constraints might have contributed to the phenomenon that 1 copy evolves faster than the other at the amino acid level after the gene duplication (Zhang et al., 2003). Gene duplication is an important mechanism that enables the morphological and functional innovation in evolution, and it may lead to distinct function for each copy by accumulating mutations that arose widely in the genome evolution (Tatusov et al., 1997).

Expression Analysis of THRSP Paralogs in Ducks
The expression profiles of the duck THRSP genes were examined by real-time PCR using the specific {alpha} For/{alpha} Rev primer set for THRSP{alpha}, and ß For/ß Rev for THRSPß (Figure 5Go). Gel electrophoresis revealed single bands of expected sizes for THRSP{alpha} (72 bp), THRSPß (137 bp), and GAPDH (108 bp) products amplified from adipogenic tissues or liver (Figure 5Go, panel D). The expression data of THRSP{alpha} and THRSPß subjected to the normalization and logarithm transformation showed the positive values only in adipogenic tissues including abdominal fat, s.c. fat, and the liver. The remaining tissues expressed THRSP genes at extremely low rates and showed negative values of relative expression (data not shown). The highest expression levels of both THRSP{alpha} and THRSPß mRNA were found identically in the fat tissues (i.e., abdominal fat or s.c. fat) of 6-wk-old Pekin ducks. The mRNA concentrations of the duck THRSP{alpha} and THRSPß were 1.5 to 2 times greater in fat tissues than that in liver (P < 0.01, Figure 5Go, panels A and B), but there was no difference between s.c. and abdominal fat tissues for THRSP paralog mRNA (P > 0.05). Expression concentrations of THRSP{alpha} and THRSPß in the duck abdominal fat were almost equal.


Figure 5
View larger version (37K):
[in this window]
[in a new window]
 
Figure 5. Expression of thyroid hormone responsive Spot 14 (THRSP) paralogs in abdominal fat (AF), s.c. fat (SF), and liver (L) from 6-wk-old Pekin ducks. Relative mRNA levels of THRSP{alpha} (panel A) and THRSPß (panel B) were measured by real-time PCR. The expression differences of THRSP paralogs in the same tissues were also analyzed (panel C). Each value represented the mean ± SEM of 6 individuals. **P < 0.01. The electrophoretic image of real-time PCR products is shown in panel D. Real-time PCR products of THRSP paralogs in abdominal fat (lane 1 = THRSP{alpha}; lane 2 = THRSPß), s.c. fat (lane 3 = THRSP{alpha}; lane 4 = THRSPß), liver (lane 5 = THRSP{alpha}, lane 6 = THRSPß), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) in fat tissues (lane 7, lane 8) were separated on a 3.5% agarose gel to document the amplification specificity of the real-time PCR. Single bands of a 72-bp fragment for THRSP{alpha}, a 137-bp fragment for THRSPß, and a 108-bp fragment for GAPDH were observed. Lane 9 = SSR DNA mark II, a DNA molecular weight standard with the fragments ranging from 41 to 533 bp. au = arbitrary units.

 
The level of THRSP mRNA in chicken liver was first found to rise by feeding a diet containing triiodothyronine (Wang et al., 2002). Although the expression of THRSP mRNA in the chicken s.c. fat was not reported, the mRNA abundances of THRSP paralogs were found consistently higher in liver than in abdominal fat tissue in 5-wk-old broiler chickens. At the same time, the mRNA concentration of THRSP{alpha} were 2 to 3 times greater than that of THRSPß in liver and abdominal fat, respectively (Wang et al., 2004). In the current study, intensive expression of THRSP mRNA was observed in adipogenic tissues of 6-wk-old ducks, as reported for the 5-wk-old chickens. However, the THRSP paralogs in 6-wk-old ducks were more actively transcribed in adipose tissues than in liver. In contrast to the expression patterns of THRSP orthologs in 5-wk-old chickens, the relative mRNA levels of THRSPß in liver and s.c. fat of 6-wk-old ducks were significantly higher than those of THRSP{alpha} (P < 0.01, Figure 5Go, panel C). In mammals, the presence of lineage-specific gene duplications could increase the human–mouse divergence of the expression patterns in the corresponding tissues of both species (Huminiecki and Wolfe, 2004). The present study on THRSP paralogs exhibited the similar picture of expression differences of duplicate genes in avian species. There is a distinct character on fat deposition and distribution in ducks when compared with the chickens, although both of them belong to the avian class. A large mass of s.c. fat can be deposited in ducks, corresponding to the excessive abdominal fat in chickens. The changes in gene expression were proposed to play a key role in the evolution course of vertebrate species (Schulte, 2004). Gu et al. (2004) pointed out that the expression patterns of duplicate genes were, in general, expected to diverge faster between species than those of single-copy gene, because changes in gene expression may lead to important phenotypic evolution. To date, the exact mechanism of THRSP involved in mammalian and avian lipid metabolism has not been well characterized. Therefore, the expression difference of THRSP mRNA in adipogenic tissues of avian species and their relationship with gene duplication awaits further investigation.

In conclusion, the current study demonstrated that the duplication of the THRSP gene is not a unique case in chickens, but it also exists in ducks and is possible in other avian species. Comparison of nucleotide and amino acid sequences in the putative coding region between the 2 paralogs showed that less difference was present in the duck THRSP compared with that of the chicken. Furthermore, expression differences of duplicated THRSP genes were presented in adipogenic tissues of the duck and the chicken. It is, therefore, the objective for further functional genomics studies to determine the physiological roles of the 2 THRSP paralogs.


    ACKNOWLEDGMENTS
 
This work is supported in part by a grant from the National Outstanding Youth Science Foundation of China (no. 30225032).


    FOOTNOTES
 
1 The cDNA sequences of duck THRSP paralogs and the intron sequence of THRSP{alpha} reported in this paper have been submitted to Gen-Bank nucleotide sequence database and assigned Accession no. DQ227766 [GenBank] for THRSP{alpha}, DQ227767 [GenBank] for THRSPß, and DQ334339 [GenBank] for the intron of THRSP{alpha}. Back

Received for publication January 9, 2006. Accepted for publication June 1, 2006.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Carre, W., C. Diot, V. Fillon, R. P. Crooijmans, S. Lagarrigue, M. Morrisson, A. Vignal, M. A. Groenen, and M. Douaire. 2001. Development of 112 unique expressed sequence tags from chicken liver using an arbitrarily primed reverse transcriptase-polymerase chain reaction and single strand conformation gel purification method. Anim. Genet. 32:289–297.[ISI][Medline]

Chagnon, Y. C., L. Perusse, and C. Bouchard. 1998. The human obesity gene map: The 1997 update. Obes. Res. 6:76–92.[ISI][Medline]

Cogburn, L. A., J. Tang, J. Cui, L. Sofer, B. Leclercq, J. Simon, and J. Burnside. 2000. DNA microarray analysis of gene expression in the liver of broiler chickens divergently selected for growth rate. Poult. Sci. 79(Suppl. 1):72.

Cogburn, L. A., X. Wang, W. Carre, L. Rejto, T. E. Porter, and S. E. Aggrey. 2003. Systems-wide chicken DNA microarrays, gene expression profiling, and discovery of functional genes. Poult. Sci. 82:6378–6383.

Compe, E., G. de Sousa, K. Francois, R. Roche, R. Rahmani, J. Torresani, M. Raymondjean, and R. Planells. 2001. Spot 14 protein interacts and co-operates with chicken ovalbumin upstream promoter transcription factor 1 in the transcription of the L-type pyruvate kinase gene through a specificity protein 1 (Sp1) binding site. Biochem. J. 358:175–183.[ISI][Medline]

Cunningham, B. A., M. Maloney, and W. B. Kinlaw. 1997. Spot 14 protein–protein interactions: Evidence for both homo- and heterodimer formation in vivo. Endocrinology 138:5184–5188.[Abstract/Free Full Text]

Ellegren, H. 2005. The avian genome uncovered. Trends Ecol. Evol. 20:180–186.[Medline]

Grillasca, J. P., M. Gastaldi, H. Khiri, A. Dace, N. Peyrol, P. Reynier, J. Torresani, and R. Planells. 1997. Cloning and initial characterization of human and mouse Spot 14 genes. FEBS Lett. 401:38–42.[ISI][Medline]

Gu, Z. L., S. A. Rifkin, K. P. White, and W. H. Li. 2004. Duplicate genes increase gene expression diversity within and between species. Nat. Genet. 36:577–579.[ISI][Medline]

Huminiecki, L., and K. H. Wolfe. 2004. Divergence of spatial gene expression profiles following species-specific gene duplications in human and mouse. Genome Res. 14:1870–1879.[Abstract/Free Full Text]

Ikeobi, C. O., J. A. Woolliams, D. R. Morrice, A. Law, D. Windsor, D. W. Burt, and P. M. Hocking. 2002. Quantitative trait loci affecting fatness in the chicken. Anim. Genet. 33:428–435.[ISI][Medline]

International Chicken Genome Sequencing Consortium. 2004. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature 432:695–716.[Medline]

Jump, D. B., S. D. Clarke, O. MacDougald, and A. Thelen. 1993. Polyunsaturated fatty acids inhibit S 14 gene transcription in rat liver and cultured hepatocytes. Proc. Natl. Acad. Sci. USA 90:8454–8458.[Abstract/Free Full Text]

Jump, D. B., P. Narayan, H. Towle, and J. H. Oppenheimer. 1984. Rapid effects of triiodothyronine on hepatic gene expression. Hybridization analysis of tissue-specific triiodothyronine regulation of mRNAS14. J. Biol. Chem. 259:2789–2797.[Abstract/Free Full Text]

Jump, D. B., and J. H. Oppenheimer. 1985. High basal expression and 3,5,3'-triiodothyronine regulation of messenger ribonucleic acid S14 in adipogenic tissues. Endocrinology 117:2259–2266.[Abstract]

Kumar, S., K. Tamura, and M. Nei. 2004. MEGA3: Integrated Software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief. Bioinform. 5:150–163.[Abstract/Free Full Text]

Lagarrigue, S., F. Pitel, W. Carre, P. Le Roy, A. Neau, Y. Amigues, J. Simon, A. Vignal, B. Leclercq, L. A. Cogburn, and M. Douaire. 2003. An initial QTL scan for abdominal fatness and breast muscle weight in broiler chickens. Plant Anim. Genome XI Conf., San Diego, CA. http://www.intl-pag.org/11/abstracts/P5l_P595_XI.htm/ Accessed Aug 2006.

Liu, H. C., and H. C. Towle. 1994. Functional synergism between multiple thyroid hormone response elements regulates hepatic expression of the rat S14 gene. Mol. Endocrinol. 8:1021–1037.[Abstract]

Lynch, M., and J. S. Conery. 2000. The evolutionary fate and consequences of duplicate genes. Science 290:1151–1155.[Abstract/Free Full Text]

Moncur, J. T., J. P. Park, V. A. Memoli, T. K. Mohandas, and W. B. Kinlaw. 1998. The "Spot 14" gene resides on the telomeric end of the 11q13 amplicon and is expressed in adipogenic breast cancers: Implications for control of tumor metabolism. Proc. Natl. Acad. Sci. USA 95:6989–6994.[Abstract/Free Full Text]

Nei, M., and S. Kumar. 2000. Molecular Evolution and Phylogenetics. Oxford Univ. Press, New York, NY

Proudfoot, N. J., and G. G. Brownlee. 1976. 3' non-coding region sequences in eukaryotic messenger RNA. Nature 263:211–214.[Medline]

Sanchez-Rodriguez, J., J. P. Kaninda-Tshilumbu, A. Santos, and A. Perez-Castillo. 2005. The Spot 14 protein inhibits growth and induces differentiation and cell death of human MCF-7 breast cancer cells. Biochem. J. 390:57–65.[ISI][Medline]

Sanger, F., S. Nicklen, and A. R. Coulson. 1977. DNA sequencing with chain terminating inhibitors. Proc. Natl. Acad. Sci. USA 74:5463–5467.[Abstract/Free Full Text]

Schulte, P. M. 2004. Changes in gene expression as biochemical adaptations to environmental change: A tribute to Peter Hochachka. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 139:519–529.[Medline]

Seelig, S., C. Liaw, H. C. Towle, and J. H. Oppenheimer. 1981. Thyroid hormone attenuates and augments hepatic gene expression at a pretranslational level. Proc. Natl. Acad. Sci. USA 78:4733–4737.[Abstract/Free Full Text]

Tatusov, R. L., E. V. Koonin, and D. J. Lipman. 1997. A genomic perspective on protein families. Science 278:631–637.[Abstract/Free Full Text]

Tuinen, M. V., and S. B. Hedges. 2001. Calibration of avian molecular clocks. Mol. Biol. Evol. 18:206–213.[Abstract/Free Full Text]

Wang, X. F., W. Carre, L. Rejto, and L. A. Cogburn. 2002. Global gene expression profiling in liver of thyroid manipulated and/or growth hormone injected broiler chickens. Poult. Sci. 81(Suppl. 1):63.

Wang, X. F., W. Carre, H. J. Zhou, S. J. Lamont, and L. A. Cogburn. 2004. Duplicated Spot 14 genes in the chicken: Characterization and identification of polymorphisms associated with abdominal fat traits. Gene 332:79–88.[ISI][Medline]

Yuan, J., B. Tang, L. Zhang, and N. Li. 2005. Comparative FISH mapping suggests poor conservation between chicken and duck genes. Anim. Genet. 36:436–437.

Zhang, P., Z. Gu, and W. H. Li. 2003. Different evolutionary patterns between young duplicate genes in the human genome. Genome Biol. 4:R56.[Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zhan, K.
Right arrow Articles by Yang, N.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zhan, K.
Right arrow Articles by Yang, N.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS