Analysis of the Rana catesbeiana tadpole tail fin proteome and phosphoproteome during T3-induced apoptosis: identification of a novel type I keratin
© Domanski and Helbing; licensee BioMed Central Ltd. 2007
Received: 12 January 2007
Accepted: 06 August 2007
Published: 06 August 2007
Thyroid hormones (THs) are vital in the maintenance of homeostasis and in the control of development. One postembryonic developmental process that is principally regulated by THs is amphibian metamorphosis. This process has been intensively studied at the genomic level yet very little information at the proteomic level exists. In addition, there is increasing evidence that changes in the phosphoproteome influence TH action.
Here we identify components of the proteome and phosphoproteome in the tail fin that changed within 48 h of exposure of premetamorphic Rana catesbeiana tadpoles to 10 nM 3,5,3'-triiodothyronine (T3). To this end, we developed a cell and protein fractionation method combined with two-dimensional gel electrophoresis and phosphoprotein-specific staining. Altered proteins were identified using mass spectrometry (MS). We identified and cloned a novel Rana larval type I keratin, RLK I, which may be a target for caspase-mediated proteolysis upon exposure to T3. In addition, the RLK I transcript is reduced during T3-induced and natural metamorphosis which is consistent with a larval keratin. Furthermore, GILT, a protein involved in the immune system, is changed in phosphorylation state which is linked to its activation. Using a complementary MS technique for the analysis of differentially-expressed proteins, isobaric tags for relative and absolute quantitation (iTRAQ) revealed 15 additional proteins whose levels were altered upon T3 treatment. The success of identifying proteins whose levels changed upon T3 treatment with iTRAQ was enhanced through de novo sequencing of MS data and homology database searching. These proteins are involved in apoptosis, extracellular matrix structure, immune system, metabolism, mechanical function, and oxygen transport.
We have demonstrated the ability to derive proteomics-based information from a model species for postembryonic development for which no genome information is currently available. The present study identifies proteins whose levels and/or phosphorylation states are altered within 48 h of the induction of tadpole tail regression prior to overt remodeling of the tail. In particular, we have identified a novel keratin that is a target for T3-mediated changes in the tail that can serve as an indicator of early response to this hormone.
Thyroid hormones (TH) are important signaling molecules in vertebrates that regulate homeostasis, growth, and development. One developmental process that is dependent upon the presence of TH is amphibian metamorphosis. During metamorphosis the larval, aquatic, herbivorous tadpole transforms into a terrestrial, carnivorous juvenile frog. This event requires drastic changes in essentially every organ and tissue of the tadpole and includes: the resorption of larval organs and tissues, remodeling of larval organs into juvenile form, and de novo development of organs and tissues . Metamorphosis is completely controlled through the control of serum TH levels. The thyroid gland mainly produces the thyroid hormone, 3,5,3',5'-tetraiodothyronine (T4 or thyroxine), which is converted into the biologically more active form, 3,5,3'-triiodothyronine (T3), in the peripheral tissues [1, 2]. Progression through natural metamorphosis is dependent upon a tightly-controlled rise and fall in TH levels. Premetamorphic tadpoles are functionally athyroid with no measurable levels of THs . TH levels gradually increase during prometamorphosis and reach maximal levels at metamorphic climax. At this stage, overt remodelling of the tadpole rapidly ensues. Premetamorphic tadpoles can be induced to undergo precocious metamorphosis by exposure to TH .
The best understood mechanism of TH action involves TH binding to nuclear thyroid hormone receptors (TRs) which regulate gene expression [2, 4]. TH binding to TRs can either activate or repress transcription of responsive genes through recruitment of coactivators and corepressors, respectively [2, 4]. Differential nuclear TR levels and intracellular T3 levels controlled by deiodinase activity and TH-binding proteins, contribute to tissue-specific responses [1, 2, 4]. However, the response to TH also depends on the existing complement of other proteins that can influence cell fate (e.g. regression of the tail versus growth and differentiation of the hindlimb). However, the molecular mechanisms are poorly understood.
Most TH-responsive genes are up-regulated and these have been most commonly studied particularly in the tail [5–10]. The genetic program required for tail regression is established between 24 and 48 h of TH exposure at a "commitment point" after which removal of TH or exposure to transcription or protein synthesis inhibitors cannot prevent regression [5, 11]. Studies based on PCR subtractive hybridization methods and cDNA gene arrays have identified a number of possible genes involved in this process [5–9]. However, the relationship of most of these genomic findings to changes in the proteome has yet to be identified and there is growing evidence for non-classical TH action through phosphorylation signaling pathways [12–16].
Proteomic scale studies on TH-dependent tail regression have been scarce. Ray et al. identified several 35S-methionine labelled proteins from Rana catesbeiana tail fin that change during natural and precocious metamorphosis using two-dimensional (2D) gel analysis . Kobayashi et al. used 2D gel electrophoresis to analyze changes in protein expression in the back and tail skin of Xenopus laevis during metamorphosis . From the 2D protein spot patterns they could classify the back skin into larval or adult type and observe the transition. Attempts were made in these studies to identify the altered protein spots. This, however, involved identifying the spots based on position, comigration or immunological detection methods. Using 2D gel separation and mass spectrometry (MS) for protein identification, we were able to identify 9 proteins whose expression was altered in the X. laevis tadpole tail during TH-induced metamorphosis . Regardless of the approach, the lack of any sample fractionation led to the identification of only abundant proteins.
In this work we identified novel changes in the proteome and phosphoproteome associated with the induction of T3-dependent tail regression of Rana catesbeiana tadpoles. Proteins that changed in abundance were detected using two-dimensional (2D) gel electrophoresis and isobaric tags for relative and absolute quantitation (iTRAQ) methods. Alterations in phosphorylation were revealed using a phosphoprotein-specific stain. Mass spectrometry (MS) analyses were then used for protein identification. Proteome coverage was enhanced through the use of cell and protein fractionation to reveal several proteins that are altered within 48 h of T3 exposure.
Results and discussion
Fractionation of the tail fin proteome and 2D gel analyses
The subcellular fractionation method was optimized to provide highly enriched fractions of cytosolic, mitochondrial, nuclear and microsomal samples from the R. catesbeiana tail fin. The goals of the method were to limit sample complexity with effective cell fractionation with minimal cross-contamination and reasonable yields. The fractions also had to be salt and buffer-compatible with subsequent proteomic analysis.
The above methods were used to analyze the proteome and phosphoproteome of the premetamorphic R. catesbeiana tail fin undergoing precocious metamorphosis at 24 and 48 h induced with 10 nM T3. A minimum of three independent replicates allowed for the verification of changes in protein and phosphoprotein expression and MS analysis was used for protein identification.
Identification of a unique R. catesbeianakeratin fragment
MS analysis of protein spot identified to be a type I keratin fragment
Observed peptide mass (Da, [M+H]+)1
Peptide sequence from MS/MS2
Identified by MALDI-TOF-TOF3
Matched database sequence5
Keratins are expressed by epithelial cells where they impart a mechanical function. In recent years, they have also been shown to be posttranslationally modified during cell stress, apoptosis, and cell signaling [27, 28]. Keratins have been extensively studied in anurans during skin differentiation and metamorphosis. The skin in R. catesbeiana transitions from a larval type into a pre-adult type, and finally into an adult type with the onset of metamorphosis [25, 26]. These changes are associated with the differentiation and apoptosis of specific epidermal cells and changes in connective tissue [25, 29]. Each of these changes have been associated with alterations in keratin type expression in specific cells [23, 30]. RLK II expression was reduced in tadpole tail and body skin with the onset of metamorphosis while RAK expression increased . These changes are precociously induced by T3 and similar changes have been observed in X. laevis with its corresponding larval and adult keratin genes [23, 24].
Keratins range in size from 40 to 67 kDa. The keratin spot corresponding to RLK I runs at ~24 kDa and all of the peptides that we detected in the MS analysis from the protein spot mapped to the N-terminal end of the complete cloned sequence (Fig. 6). Immunoblot analyses of the 48 hour microsomal fraction using a pan-cytokeratin antibody revealed the appearance of a similarly migrating fragment in the T3 sample with a concomitant reduction of protein intensity at around 50 kDa compared to the control sample (Fig. 8C). Interestingly, the identified peptides lie just upstream of what is known as a consensus caspase cleavage site (VEMDA) identified in type I human keratins  suggesting that our observed protein spot is a caspase cleavage product of RLK I (Fig. 6). No caspase cleavage of any keratins has previously been reported in anurans during metamorphosis. However, a number of effector caspases, such as caspase 3 and 7, increase in expression and activity during metamorphosis in the tail [34–36]. Caspase 3 is the most markedly up-regulated in the tail. It is expressed in larval skin epidermal cells that undergo TH-induced cell autonomous death and is known to act on type I keratins in humans [33, 37].
RLK I was also highly similar to human K19 protein. This type I keratin and the related K18 protein form heterodimers with the type II keratin K8 in simple-type epithelia and are a prevalent and extensively-studied group of keratins [27, 28, 38]. Both K18 and K19 are known caspase 3 substrates . Type I keratin caspase cleavage occurs very early in apoptosis before the detection of DNA fragmentation. Type II keratins are not caspase substrates since they lack the caspase cleavage sequence . In addition, the onset of apoptosis is associated with rapid phosphorylation of both type I and type II keratins on their head and tail domains. For K18, Ser52 phosphorylation controls its caspase cleavage and Ser33 phosphorylation regulates binding to 14-3-3 protein [27, 33]. From our results it is possible that our keratin fragment is also phosphorylated. Like K18, it contains Ser residues at positions 34 and 52. Also, the phosphoprotein stain revealed an increased phosphoprotein spot in the nuclear fraction positioned just slightly towards the more acidic end of the gel relative to the RLK I fragment. This phosphoprotein was increased at 24 and 48 h by 1.6 and 4.0 fold respectively which matches the increase of the total-protein keratin spot (Fig. 5A and 5B). Although this phosphoprotein could not be identified by MS, its more acidic position adjacent to the RLK I fragment suggests it could be a phosphorylated form of the same fragment.
The reasons for keratin cleavage and phosphorylation are not well understood and are speculated to function as phosphate sinks, for filament re-organization during apoptosis, or mechanisms that protect cells from apoptotic damage allowing a graded sequence of events. In human cancer therapy, the appearance of the caspase cleavage products of K18 and K19 in patient serum are used as indicators of cancer prognosis and apoptotic death of tumor cells undergoing chemotherapy [39, 40]. The ease with which this fragment was detected in our immunoblot using a pan-cytokeratin-specific antibody suggests that it could be used as a protein marker for the induction of a TH response in tadpoles which may be perturbed upon exposure to disruptors of TH action. The very early appearance of this fragment during precocious tail metamorphosis before the detection of any overt apoptotic events  is very intriguing. It is currently unclear whether the appearance of this keratin fragment is a product of initiation of the apoptotic event or is required for tail regression to occur.
Phosphorylation changes in γ-interferon-inducible lysosomal thiol reductase
MS analysis of protein spot identified as GILT
Observed peptide mass (Da, [M+H]+)1
Peptide sequence from MS/MS2
Matched database sequence4
LFNLV + YK
GILT is the specific lysosomal enzyme responsible for thiol reduction of proteins in the endocytic pathway for antigen presentation . Gamma interferon (IFNγ) released by activated T cells, increases antigen presentation on antigen presenting cells (APCs) such as macrophages, dendritic cells, and B cells through the induction of major histocompatability complex II (MHCII) molecules and related proteins. In addition, IFNγ increases the constitutive expression of GILT in APCs and even induces GILT expression in non-hematopoietic cells such as fibroblasts, keratinocytes and endothelial cells . GILT is transported in its proform (30–35 kDa) from the ER and Golgi complex into the early endosomes of the endocytic pathway where it is combined with MHCII-invariant chain containing vesicles and converted into its mature form (25–30 kDa) by the removal of N- and C-terminal prosequences in the late endosomes and lysosomes [41, 43]. It is likely that the GILT protein we identified is the mature form of the protein due to its localization in the microsomal and mitochondrial fractions (which contain vesicles from the ER, Golgi and lysosomes), the size on the 2D gels, and scope of coverage of identified peptides.
Further analyses of the phosphoprotein spots using 2D immunoblots and antibodies specific for the three common phospho-amino acids, phosphoserine, phosphothreonine, and phosphotyrosine did not give any indication that the GILT protein is differentially phosphorylated on any of those residues (data not shown), while the ProQ Diamond phosphoprotein stain reproducibly detected a change in phosphorylation. As with many lysosomal enzymes, GILT is a glycoprotein with three potential N-linked glycosylation sites. Furthermore, mature GILT contains mannose-6-phosphate (M6P) on one or more of the glycan chains . Phosphorylation of mannose within the N-linked glycan chain is a signal recognized in the Golgi complex by the M6P receptor that targets lysosomal enzymes to the endocytic pathway . Therefore, it is possible that the ProQ Diamond phosphoprotein stain is detecting a phosphate in M6P and that our observed increase in phosphorylation of GILT indicates an increased amount of M6P-containing GILT being directed into the endocytic pathway due to T3. In accordance with this idea of increased endocytic pathway activity, the lysosomal endoprotease, cathepsin D, increases in expression due to T3 in our iTRAQ experiment (see below). Cathepsin D, with other lysosomal cathepsin proteases (B, L and S) also converts the proform of GILT into its mature form by cleaving off the N- and C-terminal prosequences [43, 45]. A GILT knock-out study in mice, showed that GILT is required for the presentation of disulfide-containing antigens and the resulting T cell activation . These ideas, and the additional observation of increased immunoglobulin (Ig) heavy chains from our iTRAQ experiment (see below), indicate that T3 must be stimulating antigen presentation and possibly in turn activating T cell and finally Ig-producing/possessing B cell lymphocytes.
Additional changes observed in the 2D gel analysis
Within the region of the keratin spot, additional proteins and phosphoproteins were increased in the nuclear fraction (Fig. 5A). A protein spot that could not be identified increased by 2 (data not shown) and 3.5 fold at 24 and 48 h, respectively (Fig. 5B). And two unidentified phosphoproteins were increased only at 48 h by ~5 fold (Fig. 5A and 5B). A protein spot located at 30 kDa and pI ~5.5 in the microsomal fraction was increased by 2.4 fold at 48 h upon T3 exposure [see Additional file 1]. Amino acid sequence was obtained with high confidence for three peptides from this spot but no significant homology match could be made [see Additional file 2].
Differential expression analysis using iTRAQ
2D gel analysis has the advantage of identifying specific protein isoforms and posttranslational modifications of whole proteins. This whole-protein analysis approach is limited by the inability to observe hydrophobic proteins and proteins with extremes in size or pI . In contrast, MS analysis techniques circumvent this limit by analyzing peptides derived from trypsin-cleaved protein samples. However MS peptide analysis has limitations in distinguishing between protein isoforms that potentially share identical peptides . Consequently, 2D gel analysis and MS analysis techniques are complementary methods that can be combined to study differential protein expression. Therefore, to increase the number of observed altered proteins, changes in protein expression in the R. catesbeiana tadpole tail fin were additionally assessed using the novel MS technique, iTRAQ .
Traditionally, iTRAQ data analysis is performed automatically with software, which greatly reduces the analysis by the end-user. The success of obtaining valid data in this manner, however, is first dependent on the protein identification obtained using software, such as MASCOT, where an identification is reliant on a perfect match of ion peaks in a MS/MS spectrum to an existing protein in the protein sequence database . This method works well for the analysis of organisms with extensive entries in gene and protein databases but results in many unidentified spectra when less well-defined species, such as R. catesbeiana, are analyzed. Additionally, when a high quality partial MS/MS spectrum is obtained (yielding partial peptide sequence information) software limitations do not permit identification when using ion peaks database searching.
To remedy these problems, we first analyzed the raw iTRAQ data using Excel spreadsheets to find all the MS/MS spectra that showed significant changes in expression due to the T3 treatment. The resulting MS/MS spectra were de novo sequenced manually and/or using additional software (PEAKS) . The resulting peptide sequences were then used to query protein databases using BLASTp. This resulted in substantially higher identification of proteins due to homologous matches or identical matches originally missed by MASCOT.
Summary of results for iTRAQ analysis by ESI-QqTOF
Unique spectra 2
Changing peptides 3
de novo sequenced 4
Homology matched 5
MASCOT ID 6
Total proteins ID 7
Differentially expressed proteins in the Rana catesbeiana tail fin due to T3-induction as analyzed by iTRAQ
Protein name [accession #/Species]1
Observed peptide mass (Da, iTRAQ [M+H]+)3
Observed peptide sequence4
Matched database sequence6
MGC80395 protein (Sterol regulatory element-binding transcription factor 2) [AAH72922/Xenopus laevis]
Query 1 LTPATVET 8
Sbjct 184 LTPATVQT 191
Inter-α inhibitor H4 [XP_848765/Canis familiaris]
Query 1 VTFELVYEEMLK 12
Sbjct 140 VTFELVYEELLK 151
Hemoglobin α-III chain, larval [P02011/Rana catesbeiana]
Query 1 FLSFPQTK 8
Sbjct 33 FLSFPQTK 40
Hemoglobin α-III chain, larval [P02011/Rana catesbeiana]
Query 1 FLSFPQTK 8
Sbjct 33 FLSFPQTK 40
Hemoglobin α-III chain, larval [P02011/Rana catesbeiana]
Query 1 YVPHFDLTPGSADLN 15
Sbjct 42 YFPHFDLTPGSADLN 56
MGC80107 protein (Biliverdin reductase B) [AAH72790/Xenopus laevis]
Query 1 VISTPDLSHFFLR 13
VIST DLS FFLR
Sbjct 174 VISTHDLSLFFLR 186
Immunoglobulin heavy chain [AAC12909/Hydrolagus colliei]
Query 1 VVLLPPSP 8
Sbjct 133 VILLPPSP 140
Immunoglobulin heavy prechain [CAA33212/Xenopus laevis]
Query 1 SDPDQGFDGTYTVK 14
S P++ +DGT+TVK
Sbjct 393 SAPEKAYDGTFTVK 406
Immunoglobulin heavy chain constant region [AAC12914/Hydrolagus colliei]
Query 1 LNVA--DWNSGK 10
LNV+ DW SGK
Sbjct 69 LNVSTEDWKSGK 80
Immunoglobulin M heavy chain [AAO37747/Ornithorhynchus anatinus]
Query 1 FTCTVSHSDLPAP 13
Sbjct 446 FTCTVSHADLPAP 458
Immunoglobulin heavy chain variable region [AAP41191/Lepus granatensis]
Query 3 QVVEEAGGALIK 14
Q VEE+GG LIK
Sbjct 1 QVEESGGGLIK 12
Immunoglobulin heavy prechain [CAA33212/Xenopus laevis]
Query 4 FDGTYTVK 11
Sbjct 399 YDGTFTVK 406
LOC443721 protein (cathepsin D) [AAH94178/Xenopus laevis]
Query 1 AYWQIR 6
Sbjct 257 AYWQIR 262
α-2-macroglobulin [AAY98517/Xenopus laevis]
Query 1 AYVTVLGDIMGTALENLDR 19
Sbjct 958 AYVTVLGDIMGTAMQNLDR 976
Calcium-binding protein p26olf [BAA34388/Rana catesbeiana]
Query 1 GNTTSMNFK 9
Sbjct 31 GNTTSMNFK 39
α1type I collagen [BAA29028/Rana catesbeiana]
Query 1 TGPAGAPGQDGRPGPPGPPGAR 22
Sbjct 538 TGPAGAPGQDGRPGPPGPPGAR 559
α1type I collagen [BAA29028/Rana catesbeiana]
Query 1 PPGPSGEK 8
Sbjct 912 PPGPSGEK 919
Caridac α actin 2 [AAX85445/Rana catesbeiana]
Query 1 VAPEEHPTLLTEAPLNPK 18
Sbjct 98 VAPEEHPTLLTEAPLNPK 115
myosin heavy chain (skeletal muscle MHC-3) [AAD13771/Rana catesbeiana]
Query 1 FQAALEEAEASLEHEEAK 18
Sbjct 320 FQAALEEAEASLEHEEAK 337
Mylpf-prov protein (myosin light chain 2) [AAH41503/Xenopus laevis]
Query 1 NICYVITHGED 11
Sbjct 156 NICYVITHGED 166
MGC68533 protein (coatomer protein complex, subunit γ) [AAH61661/Xenopus laevis]
Query 1 NAHSLYLAGVFR 12
Sbjct 824 NAHSLYLAGVFR 835
Cortactin [NP_005222/Homo sapiens]
Query 1 YGLFPANYVELR 12
Sbjct 538 YGLFPANYVELR 549
Triose phosphate isomerase [NP_788764/Drosophila melanogaster]
Query 1 DIGADWVILGHSER 14
Sbjct 185 DIGADWVILGHSER 198
Immune system components
T3 induces the activation of the immune system as is evident with the increased maturation of GILT as revealed in the 2D analysis. In addition, this is apparent from an increase in immunoglobulin chains and cathepsin D and is possibly related to the decrease of a coatomer protein.
Six peptides that increased 1.5–4.6 fold matched a type of immunoglobulin heavy chain (Table 4). The peptides matched sequences from X. laevis as well as from other species. The peptides were derived from constant and variable regions of the heavy chain, and one matched the IgM heavy chain isotype.
The immune system plays a major role during the metamorphic process. Macrophages increase in the regressing tadpole tail in number and phagocytic activity removing apoptotic bodies of dying muscle cells . Larval lymphocyte populations rise in the growing tadpole, decrease sharply at climax of metamorphosis, and expand again with adult lymphocytes . The exchange in lymphocyte populations appears to be required for the tolerance of new adult antigens and for the removal of larval tissues that are seen as 'non-self' in the metamorphosing animal. X. laevis tail cells possess larval antigens that stimulate the immune system of a syngeneic adult or a metamorphosing animal . Izutsu and coworkers have shown that a 59 kDa larval antigen expressed specifically by tail epidermal cells only during metamorphosis causes increased proliferation of newly arising adult-type T lymphocytes, eventually leading to the destruction of tail fin tissue by cytotoxic T lymphocytes (CTL) and natural killer cells (NK) [56, 57]. Collaboration between T- and B-cell types must occur for NK cell activation and may explain our observation of increased immunoglobulin heavy chains which can only come from B cell activation. A reported peptide sequence from the 59 kDa antigen is homologous to keratin α, and points at the possibility that the increased keratin fragment observed in the 2D gel analyses may be a larval antigen produced to direct the involvement of the immune system.
A peptide that increased by 1.7 fold was identified in the MALDI-TOF-TOF analysis corresponding to a putative Xenopus cathepsin D protein (Table 4). Cathepsin D is a lysosomal aspartyl protease that has been previously implicated in tadpole tail regression [58, 59]. The increase in cathepsin D may be related to GILT maturation inside the lysosomes and a resulting increase in antigen presentation as discussed above.
A subunit γ protein of the coatomer protein complex (COPI) corresponding to the X. laevis MGC68533 protein decreased by 1.6 fold in the MALDI-TOF-TOF analysis (Table 4). Protein transport between the ER and Golgi compartment is mediated by these non-clathrin-coated vesicular coat protein complexes . COPI is composed of seven unique subunits and coats the vesicles on the cytoplasmic side as they are transported from one organelle to another. COPI is mainly associated with retrograde transport of vesicles from the Golgi back to the ER to retrieve escaped ER-resident proteins and vesicle machinery. The γ subunit binds double lysine motifs in the transmembrane proteins targeted for return to the ER. Interestingly, GILT and other lysosomal zymogens have dilysine motifs which are used as cleavage sites during enzyme maturation. There have not been any previous associations of COPI with TH or metamorphosis, but it is possible that T3 causes a specific change in vesicle trafficking to accomodate tail regression. This change may contribute to an increase in mature GILT in the lysosomes as is observed on the 2Ds of the microsomal and mitochondrial fractions.
Extracellular matrix components and modifier proteins
The plasma proteinase inhibitor, α-2-macroglobulin (A2M), was decreased 1.5 fold upon T3 treatment (Table 4). A2M is synthesized by several cell types such as hepatocytes, lung fibroblasts and monocyte-macrophages, and functions to trap and facilitate clearing of proteases such as trypsin, thrombin and collagenases [61, 62]. Tissue-specific A2M synthesis is independent of plasma A2M and serves a compartmentalized function such as the regulation of tissue proteinases such as matrix metalloproteinases. Since there is an increase in matrix metalloprotease activity in the tadpole tail in TH-induced and natural metamorphosis , the observed decrease in the A2M levels could be a result of an organized downregulation to allow for ECM remodeling or A2M and metalloprotease covalent linkage and subsequent removal and degradation of the proteins. Recently, Lin et al. observed a T3-dependent decrease in A2M transcripts in a human hepatocellular carcinoma cell line .
A peptide matching inter-α inhibitor H4 protein increased by 2.1 fold (Table 4). The H4 polypeptide is a less described component of the trimeric plasma serine protease inhibitor, inter-α inhibitor (IaI), more commonly composed of the H1, H2 and bikunin polypeptides [64, 65]. IaI proteins are linked to inflammation processes and ECM structure. Serine protease inhibitor activity of bikunin regulates plasmin which is a known activator of matrix metalloproteases . The H4 polypeptide is also a serum marker for the acute-phase response in infected animals and humans resulting from the nonspecific immune response [66–68] and may be an additional link to our observation of the immune system activation. Additionally, the H polypeptides stabilize and protect the ECM through covalent links to ECM components and possibly through collagen binding domains .
α 1 type I collagen (α 1(I) chain) was identified by two peptides that decreased by 1.6–1.9 fold (Table 4). Type I collagen is the most abundant member of the collagen family and is the main constituent of the ECM . It is a major substrate of matrix metalloproteinase 1 (MMP1) which is upregulated in R. catesbeiana tail and back skin during natural and T3-induced metamorphosis [70, 71]. Collagen (α 1(I) and α 2(I)) gene expression also decreased under these conditions in the tail and intestine [72–74]. It is possible that the α 1(I) collagen gene is negatively regulated by a TRE-like site as in rat cardiac fibroblasts . Here we show that this change also occurs at the protein level.
The last peptide associated with ECM remodeling proteins corresponds with cortactin (Table 4). A 1.5 fold decrease in this Src homology 3 (SH3) domain-containing protein suggests altered tyrosine kinase signaling of actin cytoskeleton rearrangements . Cortactin is a filamentous-actin binding protein that stabilizes actin networks and promotes actin polymerization. Since cortactin interacts with various proteins involved in actin polymerization, endocytosis (recall an overlap with GILT and cathepsin D proteins) and cell-cell interactions, modification of any or all of these processes could contribute to tail regression.
Cell metabolism proteins
Two proteins involved in cell metabolic homeostasis were altered in abundance by T3 treatment (Table 4). The first peptide was increased by 1.6 fold matched with a sterol regulatory element binding protein-2 (SREBP-2). This membrane-bound transcription factor regulates lipid homeostasis through the control of genes encoding key enzymes important in cholesterol synthesis . The SREBP transcription factors have not been implicated in anuran metamorphosis before, however, they are involved with TH in the mammalian context [78–80]. Recently Shin and Osborne  observed that T3 regulates both SREBP-2 mRNA and nuclear protein levels, as well as SREBP-2-controlled genes such as the low density lipoprotein (LDL) receptor in mice. The SREBP-2 promoter contains at least one TR-binding site and is activated directly by TR in a ligand-dependent manner . It is therefore possible that the increased level of SREBP-2 observed in our experiments may be a direct response to T3. This remains to be determined.
The second peptide was reduced by 1.6 fold and corresponds to triose phosphate isomerase (TIM; Table 4). This enzyme is involved in glycolysis  and its decrease is compatible with a shift from glycolysis to oxidative phosphorylation observed during metamorphosis .
In addition, two proteins involved in oxygen transport were increased upon T3 treatment. Three peptides corresponding to larval hemoglobin α-III chain increased by 1.5 fold (Table 4). This protein is only found in tadpole red blood cells (RBCs) and not adult RBCs [83, 84]. In R. catesbeiana there is a change in erythropoiesis sites and hemoglobin types during tadpole growth and metamorphosis that is influenced by THs [83, 85–88]. Young tadpoles produce hemoglobin components I and II in the kidneys, older tadpoles produce components III and IV in the liver, and after metamorphosis, only adult hemoglobins are produced in the bone marrow. This sequential change produces hemoglobins with lower oxygen affinities. Our observation of an increase in hemoglobin α-III is indicative of an accelerated progression of hemoglobin switching due to T3 exposure.
Consistent with a switch in hemoglobin is an 1.5 fold increase in a peptide matching biliverdin reductase B (BVR-B; Table 4). This enzyme is important in hemoglobin metabolism and is predominantly found in human fetal liver . There, it produces bilirubin IXβ from biliverdin derived from heme. Bilirubin IXβ is the major heme catabolite during human fetal development . In adult human liver biliverdin reductase A (BVR-A) is the major form and it produces bilirubin IXα. The switch in this heme degradation pathway seems to be related to the switch in hemoglobins from fetal to adult forms.
Mechanical function and calcium modulated proteins
The abundance of peptides from muscle α-actin, myosin heavy chain (MHC) and myosin light chain (MLC) was reduced by exactly the same amount (2.6 fold; Table 4). These proteins are the fundamental components of the molecular machinery used in muscle contraction. The peptide for muscle α-actin is distinguished from the β and γ forms of cytoplasmic actins by a T residue instead of V in the HPTLL sequence within this peptide. The peptide sequence matches the R. catesbeiana cardiac α-actin 2 sequence perfectly as well as other muscle α actins in other species. The MHC peptide matches the R. catesbeiana skeletal muscle MHC-3 and is distinguishable from the other isoforms, MHCs -1, 2, 4, and 5, by three residues in this species . The MLC peptide matches the X. laevis MLC-2 (Mylpf-prov protein) as well as MLC-2 proteins from other species.
An embryonic/larval α-3 skeletal-actin gene has been found in Xenopus borealis . This gene was expressed in the skeletal muscle of embryos and tadpoles, but not in adults. However, the corresponding gene in X. laevis is highly expressed in both adult muscle and at earlier stages. The human skeletal muscle α-actin gene is a direct TH response gene in cardiocytes [93, 94]. The cardiac α-actin gene of rats is also inducible with T3 in heart myocytes .
Five skeletal muscle MHC isoforms have been cloned in R. catesbeiana . No specific larval forms have been found and all five isoforms are expressed in the tail and limb muscle of both tadpoles and adults. In the tadpole tail muscle, T3 increases the mRNA for MHCs -1, 3, 4, and 5 by two to three-fold within 48 h . However, in contrast, the rate of total MHC protein synthesis was reduced by 40% at 48 h  consistent with our observations for MHC-3. In fast twitch muscle fibers, MHC-3 is associated with MLC-2 . MLCs are Ca2+-signaling dependent regulators of muscle contraction and have been investigated in X. laevis and R. pipiens but their responses to TH are not well understood . Our data suggest that MHC and MLC expression are co-regulated.
A peptide that was decreased by 2.1 fold matched the p26olf protein of R. catesbeiana (Table 4) The p26olf protein is a S100-like Ca2+-binding protein with the very unique feature of having two S100-like regions . Thus far p26olf has only been found in the anurans R. catesbeiana, X. tropicalis and X. laevis, and appears not to have any homologues in other vertebrates. Immunohistochemical studies on R. catesbeiana frogs showed that p26olf primarily localizes in the cilia of olfactory and lung respiratory epithelia. p26olf may be involved in olfactory signal transduction and a related protein, S100/A11 (calgizzarin), is a key mediator of high Ca2+-induced growth arrest and differentiation in human keratinocytes . In addition, Makino et al.  have shown that introducing an N-terminal peptide of S100/A11 into human carcinoma cell lines induced apoptosis. Thus, it may be possible that p26olf is involved in the control of apoptotic events in the tadpole tail.
Genomic studies on tadpole tail regression have revealed a number of modulated transcripts [5–9]. In comparison, our extensive proteomic analyses showed relatively few changes within the first 48 h of T3 treatment. This provides support for the hypothesis that changes in the transcriptome are the major target for the establishment of the tail regression program. However, this does not imply that the few changes in the (phospho)proteome during this time are not important, as these changes may support the transcriptome transition. For example, the tyrosine kinase inhibitor, genistein, prevents T3-dependent tail regression likely through inhibition of protein kinase C activity and TRα phosphorylation .
Our proteomic analyses may be limited in their sensitivity to reveal phosphorylation of low-abundance signaling phosphoproteins since no enrichment was performed. Nevertheless our data suggest that RLK I could be posttranslationally modified due to a T3 signal most likely through the early action of caspases thereby indicating an early apoptotic event. Our findings with GILT maturation, Ig heavy chains and cathepsin D increases show that the involvement of the immune system in the regressing tadpole tail may be greater than what has been generally considered. Furthermore, a large proportion of the transcripts identified by genomic studies encode transcription factors [5–9] whose protein products are most likely below the detection limits of our 2D PAGE and iTRAQ analyses. An improved fractionation strategy for the iTRAQ analysis may reveal these.
The present study provides important insights into the molecular mechanism of T3-dependent tadpole tail regression during postembryonic development which includes novel components at the proteome level.
The care and treatment of animals used in this study were in accordance with the guidelines of the Animal Care Committee, University of Victoria, under the Canadian Council of Animal Care.
Rana catesbeiana tadpoles were locally caught (Victoria, BC) or purchased (Ward's Natural Science Ltd., St, Catharines, ON). Animals were housed in the University of Victoria aquatics facility and maintained in 360 L all-glass flow-through aquaria containing recirculated water at 15°C with exposure to natural daylight. Tadpoles were fed daily with spirulina (Aquatic ELO-Systems, Inc., FL). R. catesbeiana tadpoles, TK stage VI-VIII , were exposed to 10 nM T3 (80 μl of 10-3 M T3 in 2.5 mM NaOH) (Sigma-Aldrich, St. Louis MO USA) or vehicle control (80 μl 2.5 mM NaOH) in the rearing water (8 L per 6 tadpoles) for 24 or 48 h. Tadpoles were acclimatized for 24 h prior to exposures to 24°C, and were not fed during this period or during exposures. For the 2D gel analysis, each treatment or control group consisted of 18 pooled tadpoles. The 48 hour 2D analyses consisted of three, completely independent, control and T3 treatment groups, while the 24 hour analyses were from a single control and T3 treatment group. The 24 hour analyses were only done to further characterize absence or presence of proteins observed to be altered at the 48 hour time point. For the iTRAQ analysis, two independent T3-treatment groups and two independent control groups consisting of six pooled tadpoles per group were used. Tadpoles were euthanized in 0.1% tricaine methanesulfonate (Syndel Laboratories, Vancouver, BC, Canada) in 25 mM sodium bicarbonate, dorsal and ventral tail fins were removed on ice, and proteins extracted as described below. For QPCR analysis, premetamorphic tadpoles (stage VI-VIII) were exposed to 100 nM T3 or vehicle control in the rearing water for 24, 48 and 72 h and euthanized as above. Tadpoles at the indicated stages were also collected to determine gene expression in tail tissue during natural metamorphosis.
Two separate procedures were used for the isolation of nuclei and for the cytosolic-mitochondrial-microsomal fractions. All extraction procedures were performed on ice or at 4°C. Tail fin tissue was rinsed with extraction buffer lacking detergents or inhibitors (see below for composition). Three ml of buffer were used per gram of tissue, which was homogenized in a teflon-glass Dounce homogenizer (Wheaton, USA). For the nuclear isolation tissue was homogenized in 20 mM N-[2-Hydroxyethyl]piperazine-N'-[2-ethanosulfonic acid] (HEPES) pH 7.5, 0.25 M sucrose (ACP Chemicals Inc., Montreal Canada), 25 mM KCl, 5 mM MgCl2 (both from BDH chemicals, Toronto Canada), 0.5% Triton X-100, 1mM dithiothreitol (DTT), 2 mM vanadyl ribonucleoside complexes (VRCs), 100 μM phenylmethylsulfonyl fluoride (PMSF), 4 μg/ml aprotinin, 1 μg/ml leupeptin, 2 μg/ml antipain, 300 μg/ml benzamidine, 10 mM β-glycerophosphate, 0.2 mM Na3VO4, and 1 mM NaF (all remaining from Sigma-Aldrich).
Homogenization efficiency and integrity of nuclei was checked with microscopy and 4'-6-Diamidino-2-phenylindole (DAPI) staining for nuclei. The homogenate was incubated on ice with shaking for 30 min then passed through 100 μm nylon mesh cell-strainer. The homogenate was centrifuged at 1000 × g for 10 min. The nuclei were resuspended in the above extraction buffer (3 ml per g of initial tissue), layered onto a 0.88 M sucrose cushion, which contained everything the extraction buffer did except for Triton X-100 and VRCs, and centrifuged at 2000 × g for 30 min. The resulting nuclear pellet was washed in extraction buffer, without Triton X-100 and RVCs, and centrifuged at 2000 × g for 10 min. The nuclear pellet was stored at -70°C.
To obtain the cytosolic, mitochondrial and microsomal fractions, tail fin tissue was homogenized in 20 mM HEPES pH 7.5, 0.25 M sucrose, 5 mM Ethylenedinitrilo-tetraacetic acid (EDTA) (EM Science, Gibbstown NJ USA), 1mM DTT, 100 μM PMSF, 4 μg/ml aprotinin, 1 μg/ml leupeptin, 2 μg/ml antipain, 300 μg/ml benzamidine, 10 mM β-glycerophosphate, 0.2 mM Na3VO4, 1 mM NaF. The homogenate was centrifuged at 800 × g for 10 min. The supernatant was then centrifuged at 12,200 × g for 15 min. The pellet was washed twice in extraction buffer and centrifuged at 12,200 × g for 15 min. This mitochondrial pellet was stored at -70°C. The supernatant was centrifuged at 100,000 × g for 1 hour to obtain the microsomal pellet and cytosolic fractions which were then stored at -70°C.
Anion-exchange HPLC fractionation of cytosolic fraction
The cytosolic fraction was first desalted using size-exclusion chromatography with Sephadex G-25-M beads (Sigma-Aldrich). Sample was kept at 4–8°C throughout the procedure. Six ml of sample was loaded onto a 27 ml column with 60 mM ammonium bicarbonate (EM Science) as buffer. The eluted protein fractions in the first 7.5 ml were kept, frozen at -70°C, and lyophilized for 48 h. The sample was reconstituted in 40 mM ammonium bicarbonate to a concentration of 10 mg/ml. Fractionation of the desalted cytosolic sample was performed using a Waters HPLC setup with a Waters 600S controller, 626 pump and 486 tunable absorbance detector (Waters, Milford MA USA). The column used consisted of 2 ml Accell Plus QMA 300Å anion-exchange media packed into a AP Minicolumn assembly (5 × 100 mm) (Waters). Five mg of total cytosolic sample was loaded onto the column in 40 mM ammonium bicarbonate (pH 8.3), the unbound proteins were captured as the 40 mM fraction. The remaining sample was fractionated using a step-gradient of 190, 260, 340 mM and 1 M ammonium bicarbonate. The flow rate was set at 1 ml/minute and each step was 13 min. The detector was set at 280 nm. The fractions were collected on ice, frozen at -70°C, and lyophilized for 48 h.
2D polyacrylamide gel electrophoresis
Proteins from the nuclear pellet were initially extracted in an isoelectric focusing (IEF) buffer lacking ampholytes consisting of 9.5 M urea (SigmaUltra), 4% 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate (CHAPS), 1% DTT, 40 mM tris(hydroxymethyl)aminomethane (Tris-base) pH 8.5, 5 mM EDTA, 100 μM PMSF, 4 μg/ml aprotinin, 1 μg/ml leupeptin, 2 μg/ml antipain, 300 μg/ml benzamidine, 10 mM β-glycerophosphate, 0.2 mM Na3VO4, 1 mM NaF (all from Sigma-Aldrich) for 15 min at room temperature. DNA was then precipitated by adding 2% ampholytes (Pharmalyte 3–10 and 5–8 in a 2/3 and 1/3 ratio, respectively; Amersham Pharmacia Biotech AB, Uppsala Sweden), incubating for 15 min at room temperature, and pelleting the DNA at 12,000 × g for 10 min. Ampholytes were then replaced in the supernatent to 2% as indicated above. The mitochondrial and microsomal pellets and the lyophilized HPLC cytosolic fractions were all dissolved in the IEF buffer with the inclusion of 2% ampholytes as described above. Protein concentrations were assayed using the Bio-Rad Protein Assay (Hercules, CA USA) with protein standards dissolved in the IEF buffer. The amount of pellet to IEF buffer used yielded a protein extract with a protein concentration of 2–6 mg/ml.
Three hundred μg of protein for each fraction (except for the 40 mM cytosolic fraction where 100 μg of protein was used) were aliquoted, made up to 100 μl with IEF buffer, solubilized for seven hours at room temperature, centrifuged at 12,000 × g for 10 min, and loaded onto IEF tube gels used for the first dimension. The IEF tube gels (2.5 mm × 13 cm) consisted of 9.5 M urea, 3.897% acrylamide and 0.104% piperazine diacrylamide (wt/vol) (Bio-Rad), 2% CHAPS, 3.875% Pharmalyte 3–10 ampholytes, 2.375% Pharmalyte 5–8 ampholytes, and 0.05% ammonium persulphate (APS) (ACP Chemicals Inc.) and 0.07% N,N,N',N'-tetramethylethylenediamine (TEMED; Sigma-Aldrich) for polymerization.
The IEF gels were prefocused and run as described in  with the following modifications. Samples were electrophoresed for 16 h at 450 V, then hyperfocused at 800 V for 2 h (8800 Vh total). The IEF gels were then equilibrated using the same equilibration buffer as in  except the buffer contained 65 mM DTT for 15 min, and then for 15 min in 5 ml of IEF-equilibration buffer with 136 mM iodoacetamide (Sigma-Aldrich).
The IEF gels were then transferred onto the second dimension SDS-polyacrylamide gels (15 cm × 14 cm × 1.5 mm) comprised of an 11% separating gel (37.5:1 acrylamide to bis-acrylamide (Bio-Rad), 375 mM Tris-base pH 8.8, 0.1% SDS, 0.05% APS, 0.05% TEMED) and a 4% stacking gel (125 mM Tris-base pH 6.8, 0.1% SDS, 0.05% APS, 0.1% TEMED). The IEF gels were overlayed with 0.5% agarose (EM Science) in 125 mM Tris-HCl pH 6.8 with 0.05% bromophenol blue/2% SDS. Precision Plus Dual Color (Bio-Rad) molecular weight markers were used. The gels were electrophoresed at 30 mA per gel at constant current in SDS running buffer (25 mM Tris-base, 192 mM glycine (EM Science), 0.1% SDS) until the bromophenol blue ran out of the bottom of the gel,.
The 2D gels were then stained with a fluorescent phosphoprotein stain, ProQ Diamond, according to the manufacturer's instructions (Invitrogen-Molecular Probes, Burlington, Ontario, Canada). Fluorescence was detected with UV light using a ChemiImager 4000 imaging system (Alpha Innotech Corporation, San Leandro CA USA). The gels were then stained for total protein by washing in 50% ethanol (EM Science) and 3% phosphoric acid with shaking overnight, rinsing three times for 30 min in water, equilibrating in 16% ammonium sulfate (EM Science), 25% methanol (EM Science), 5% phosphoric acid for 1 hour, and subsequently staining by adding Coomassie brilliant blue G250 (Sigma-Aldrich) to 0.01% and shaking for 4 days. Images were obtained by scanning the gels with a Li-Cor Odyssey scanner (Li-Cor, Lincoln NE USA) at 700 nm. Protein spot density was analyzed using Li-Cor Odyssey Ver. 2.0.5 software (Li-Cor). Gel images were analyzed using Adobe PhotoShop Ver. 5.0 (Adobe Systems Inc., San Jose CA USA).
Mass spectrometry analysis
Protein spots of interest were excised from the 2D gels and the proteins within were reduced, alkylated, and digested with trypsin according to an in-gel digestion protocol  with a few modifications. The gel pieces were destained in 50% methanol/5% acetic acid (ACP Chemicals Inc.), dehydrated with acetonitrile (EM Science), dried, reduced with 50 mM DTT in 100 mM ammonium bicarbonate at 56°C for 30 min, alkylated with 100 mM iodoacetamide in 100 mM ammonium bicarbonate at 45°C for 30 min, dehydrated with acetonitrile, hydrated with 100 mM ammonium bicarbonate, dehydrated with acetonitrile, dried, and digested with 20 ng/μl of sequencing grade modified trypsin (Promega Corp., Madison WI USA) in 50 mM ammonium bicarbonate at 37°C overnight. The resulting peptides were extracted out of the gel pieces by successive incubation in 50 mM ammonium bicarbonate for 1 hour at 37°C and twice in 0.1% formic acid (ACP Chemicals Inc.) for 30 min at 37°C. The peptides were desalted using ZipTip pipette tips containing C18 reversed-phase media (Millipore Corp., Bedford, MA USA) by washing with 0.1% formic acid and eluting with 75% acetonitrile/0.1% formic acid. Peptides were either analyzed by MALDI-TOF MS as follows or by ESI-QqTOF or MALDI-TOF-TOF MS analysis as described in the iTRAQ section. Peptide samples were applied to the target plate with an equal volume of matrix solution (1% α-cyano-4-hydroxycynnamic acid (Sigma-Aldrich) in 50% acetonitrile/0.3% formic acid) and allowed to dry. Adrenocorticotropic hormone fragment 1–17 (FW 2093.4), bradykinin fragment 2–9 (FW 904.0) and angiotensin 1 (FW 1296.5) (Sigma-Aldrich) were used as external calibrants.
Spectra were obtained using a Voyager-DE STR Biospectrometry Workstation MALDI-TOF mass spectrometer (PE Applied Biosystems, Foster City, CA USA) operating in positive reflector mode with delayed extraction. Data were manipulated using the Voyager Version 5.1 software with Data Explorer (PE Applied Biosystems). The sample spectra were further internally calibrated using autolytic trypsin peptide peaks. The masses of the tryptic peptides were used to search a local protein database using MASCOT Peptide Mass Fingerprint software (Matrix Science Inc., Boston MA USA) using the following settings: all entries, 1 missed cleavage, cysteines modified by carbamidomethylation, oxidized methionines, deamidation of asparagine and glutamine, peptide mass tolerance of ± 50 ppm. If the tryptic peptides were analyzed by ESI-QqTOF generating MS/MS spectra, manual de novo sequencing was performed viewing the spectra on AnalystQS Ver. 1.1 software (Applied Biosystems, Foster City, CA USA) or automatically using PEAKS Studio Ver.3.0 software (Bioinformatics Solutions Inc., Waterloo ON Canada) as explained in the data analyses section for iTRAQ. The group of peptide spectra that were obtained from the RLK I spot were analyzed using SPIDER software .
Protein samples from the total extract and subcellular fractions were analyzed for lamin B1+B2 and for cytochrome c. Protein samples of the 48 hour control and treatment microsomal fraction were analyzed for cytokeratin in duplicate. Thirty μg of protein per sample was boiled for 3 min in 25 μl of 60 mM Tris-base pH 6.8, 2% SDS, 10% glycerol, 100 mM DTT, 0.03% bromophenol blue and loaded onto a 12% SDS-PAGE (see '2D gel electrophoresis' above for gel and buffer compositions). Proteins were then electrophoretically transferred onto a nitrocellulose membrane (Bio-Rad). Membranes were blocked with 5% nonfat milk in 140 mM NaCl/20 mM Tris-base pH 7.6 (TBS), for 1 hour, then probed with a pan-cytokeratin rabbit polyclonal IgG antibody (H-240; Santa Cruz Biotechnology Inc., Santa Cruz CA USA) diluted at 1/500 or an anti-lamin B1+B2 mouse monoclonal IgG1 antibody (ab4825; Cedarlane, Hornby ON Canada) diluted at 1/50 or anti-cytochrome c rabbit polyclonal antibody (H-104; Santa Cruz Biotechnology Inc.) diluted at 1/200 in 5% nonfat milk/TBS and 0.1% Tween (TBST) with gentle agitation for 1 hour at room temperature. The blots were washed with TBST for 1 hour. The membranes were then incubated with secondary antibody, IRDye-800 conjugated anti-rabbit IgG or IRDye-800 conjugated anti-mouse IgG (Rockland Inc., Gilbertsville PA USA) diluted at 1/2000 in 5% nonfat milk/TBST for 1 hour at room temperature. The blots were washed with TBST for 1 hour. The prepared blots were scanned by a Li-Cor Odyssey scanner at 800 nm.
Isolation of RNA, generation of cDNA, degenerate-primer PCR and 5'-/3'-RACE
RNA was isolated and cDNA generated from premetamorphic tadpole tail fin tissue as described in Veldhoen et al . The resulting cDNA was used for PCR using degenerate primers: forward: 5'-GAA/G GCA/C GCC AAT/C ACT/C GAA/G CT-3' and reverse: 5'-AC C/TTG T/AAG C/TTC CTC T/CTC G/ATG-3' (32 fold degeneracy each) with an initial 94°C for 10 min then 40 cycles at 94°C for 30 sec, 52°C for 30 sec, and 72°C for 30 sec, followed by 10 min at 72°C. This generated a single product at 380 bp. The PCR product was cloned using the TOPO TA Cloning Kit System with pCR2.1 TOPO vector and TOP10 cells (Invitrogen) as described previously . The insert was then sequenced in both directions using the TOPO kit M13 forward and reverse primers at the University of Victoria sequencing facility. Based on this sequence two gene specific primers (GSP) were designed: GSP1 for 5'-RACE: 5'-TGTCTCCACCGCTAACACCAACGCCAAC-3' and GSP2 for 3'-RACE: 5'-GTTGGCGTTGGTGTTAGCGGTGGAGACA-3'. RACE-ready cDNA libraries were made from total RNA using the Smart RACE cDNA kit (BD Biosciences, Burlington ON Canada). The RACE reactions used the universal primers supplied with the kit paired with the GSP1 and GSP2. Touchdown PCR was performed on the RACE-ready libraries: 94°C for 2 min; five cycles at 94°C for 30 sec, and 72°C for 3 min; five cycles at 94°C for 30 sec, 70°C for 30 sec, and 72°C for 3 min; and 25 cycles at 94°C for 30 sec, 68°C for 30 sec, and 72°C for 3 min. The 5'-RACE produced a single product of ~600 bp and the 3' RACE produced three products between ~1400 and ~1500 bp. The PCR products were cloned as above and the longest clones from each RACE reaction were sequenced as above, giving two overlapping sequences that encoded the complete open reading frame. The sequence reported in this manuscript is accessible in GenBank [GenBank: EF156435].
Quantitation of gene expression
The expression of RLK I gene transcript was quantified in tails of individual animals (n = 4 for all treatments and/or time points). Tail tissue treatment, RNA isolation and cDNA generation was performed as indicated above and expression data were collected using a real-time quantitative PCR assay on a MX4000 system (Stratagene, La Jolla, CA, USA) as described previously  using the following primers for RLK I: 5'-GTTGGCGTTGGTGTTAGCGG-3' and 5'-GGCACTGCTTCTTGCAACTTG-3'. The QPCR product from the keratin transcript was analyzed on a gel and sequenced ensuring that only this specific keratin transcript was being analyzed. The invariant L8 ribosomal protein transcript was used as a normalizer .
Differential expression analysis using iTRAQ
Peptide preparation, iTRAQ labeling, two-dimensional liquid chromatography separation and MS analysis
Proteins from tadpole tail fin tissue were extracted by homogenization in 6 M urea, 0.2% SDS, 20 mM HEPES pH 7.5 (3 ml/g tissue). The homogenate was centrifuged at 12,000 × g for 15 min at 4°C. The supernatant was adjusted to 4 M urea, 0.05% SDS, 2 mM MgCl2 and 20 mM HEPES pH 8.0, and DNA was digested with 10 U/ml of benzonase for 30 min at 8°C. The iTRAQ reagents were used according to the manufacturer's protocol with minor modifications as follows (Applied Biosystems).
Two independent control and two independent treatment protein samples were quantified using the bicinchoninic acid assay (BCA) (Sigma-Aldrich, St Louis, MO USA). One hundred μg of each sample (~55 μl) was precipitated in 1 ml cold acetone, overnight at -20°C. Proteins were then reduced in 3.3 mM tris-(2-carboxyethyl) phosphine, alkylated with 6.7 mM methyl methane thiosulfonate (MMTS) and digested overnight with 10 μg of sequencing grade modified trypsin (Promega, Madison, WI USA) at 37°C. The resulting peptides from the two control samples were labeled with iTRAQ reagents 114 and 116, while peptides from the two T3 treatment samples were labeled with iTRAQ reagents 115 and 117.
The labeled peptide samples were then pooled, adjusted to pH 2.5–3.0 with concentrated phosphoric acid (ACP Chemicals Inc., Montreal, QC Canada), and separated by strong cation-exchange chromatography (SCX) using a PolySULFOETHYL A SCX column (100 × 4.6 mm, 5μm, 300Å - PolyLC Inc., Columbia, MD USA) on the VISION Workstation (Applied Biosystems, Foster City, CA USA). Mobile phases used were: Buffer A: 10 mM monobasic potassium phosphate (Sigma-Aldrich, St Louis, MO, USA), 25% acetonitrile (ACN) (EMD Chemicals, Gibbstown, NJ, USA), pH 2.7 and Buffer B: same as A with the addition of 0.5 M potassium chloride (Sigma-Aldrich St Louis, MO, USA). Fractions of 500 μL were collected over an 80 minute gradient: 0–30 min, 5% to 35% Buffer B; 30–80 min, 35% to 100% Buffer B. Sixteen SCX fractions were reduced to 150 μL using a SpeedVac and subjected to reverse-phase chromatography using an integrated system consisting of a Famos autosampler, Switchos switching pump, UltiMate micro pump and a Probot microfraction collector (LC Packings, Amsterdam, Netherlands). The SCX fractions were first desalted on a C18 PepMap guard column (300 μm i.d. x 5 mm 5 μm, 100Å, LC Packings, Amsterdam) at 50 μl/minute flow rate for 15 min. Peptides were separated over a manually packed 75 μm × 15 cm C18 column (Magic C18Aq, 5μm, 100Å, Michrom Bioresources Inc., Auburn CA, USA) using a 85 min gradient of 5%-75% ACN in 0.1% formic acid flowing at 250 nL/min.
The eluant was either sprayed directly into an Applied Biosystems QSTAR Pulsar i (Applied Biosystems/MDS SCIEX Concord, ON Canada) (ESI-QqTOF) mass spectrometer or spotted onto a MALDI target plate by the Probot for LC-MALDI analysis. The QSTAR operating software Analyst QS v1.1 employed an information dependent acquisition (IDA) method for optimized MSMS spectra acquisition over a 6 sec cycle which was repeated over the duration of the gradient. The MS-TOF survey scan lasted for 1 sec over a range of 400–1200 m/z targeting precursor ions of charge state 2–5+ that exceeded a threshold of 20 counts. Former target ions within 100 ppm were excluded for the next 180 sec. The two most intense ions that met the IDA criteria were gated and fragmented with optimized collision energy. Each product ion scan lasted for 2.5 sec over a range of 100–1500 m/z. 'Enhance-all' was turned on for the product ion scans. LC-MALDI data was acquired on an Applied Biosystems 4800 MALDI TOF/TOF Analyzer (Applied Biosystems/MDS Sciex Foster City, CA USA). Sample spots were overlaid with 3 mg/mL α-cyano-4-hydroxycinnamic acid (CHCA) (Sigma-Aldrich, St Louis, MO USA) matrix in 50% ACN and 0.1% TFA. MS data was automatically acquired over a mass range of 800–4000 Da using fixed laser intensity for 500 shots with a uniformly random spot search pattern. In each spot, the 30 strongest peaks by MS (determined by cluster area) were selected for MS/MS using a job-wide interpretation method which excludes ions with a signal-to-noise of less than 60 and filters identical peaks detected in multiple spots. A 1 KV positive MS/MS operating mode was used with the relative precursor mass window set at 150 and collision induced dissociation turned on. MALDI-TOF-TOF data analysis was performed using GPS Explorer Ver. 3.0 (Applied Biosystems) software using the MASCOT (Matrix Science) search engine for protein database searching.
iTRAQ data analysis
All of the MS data in ProQUANT (Applied Biosystems) resulting from the ESI-QqTOF analysis was transferred to Excel (Microsoft Corporation, Redmond WA) where the data was filtered to reveal only those MS spectra where a change of 1.5 fold or higher was observed between the treatment and the control for both duplicates. An additional requirement was that the ratio of the iTRAQ tags between the two controls had to be between 0.67 and 1.5. The resulting MS spectra were then manually inspected in AnalystQS 1.1 (Applied Biosystems) software to ensure good quality of the iTRAQ tags and fragmentation information.
High quality spectra where then used for de novo peptide sequencing. Peptide sequence was obtained by either manually interpreting spectra and/or through the use of PEAKS Studio Ver.3.0 (Bioinformatics Solutions Inc.) software capable of auto de novo sequencing. The parameters were set as follows: parent mass error tolerance of 0.1 Da, fragment mass error tolerance of 0.1 Da, enzyme used was trypsin, fixed posttranslational modifications (PTMs) were MMTS on cysteines, iTRAQ tags on N-terminal and lysine, and variable PTMs were iTRAQ tag on tyrosine and oxidation of methionine histidine and tryptophan. All of the results were manually inspected to ensure correctness. The resulting peptide sequences were then queried against the Entrez Protein (NCBInr) database (Oct. 15, 2006) using BLASTp  to allow for short, nearly exact matches: no compositional adjustments, a low complexity filter, expect threshold of 20000, word size of 2, using matrix PAM 30, with gap costs of 9 for existence, and 1 for extension, searching all metazoan sequences.
Statistical analyses were performed using SPSS Ver. 12.0. (Chicago IL USA). For the 48 hour 2D gel analysis data P-values were determined using a one-way ANOVA with a test for homogeneity of variance and a Shapiro-Wilk test of normality. For the QPCR analyses P-values were determined using the Mann-Whitney U non-parametric two-tailed test.
We thank Darryl Hardie, Leanne Ohlund, Derek Smith and Monica Elliott from the University of Victoria – Genome British Columbia Proteomics Centre for assistance in the mass spectrometry analysis and helpful technical advice. This work was supported by NSERC and a MSFHR scholar award to CCH, and a NSERC PGS A and PGS D award to DD.
- Shi YB: Amphibian Metamorphosis: from morphology to molecular biology. 2000, New York, Wiley-LissGoogle Scholar
- Furlow JD, Neff ES: A developmental switch induced by thyroid hormone: Xenopus laevis metamorphosis. Trends Endocrinol Metab. 2006, 17: 40-47. 10.1016/j.tem.2006.01.007.PubMedGoogle Scholar
- Regard E, Taurog A, Nakashima T: Plasma thyroxine and triiodothyronine levels in spontaneously metamorphosing Rana catesbeiana tadpoles and in adult anuran amphibia. Endocrinology. 1978, 102: 674-684.PubMedGoogle Scholar
- Tata JR: Amphibian metamorphosis as a model for the developmental actions of thyroid hormone. Mol Cell Endocrinol. 2006, 246: 10-20. 10.1016/j.mce.2005.11.024.PubMedGoogle Scholar
- Wang Z, Brown DD: Thyroid hormone-induced gene expression program for amphibian tail resorption. J Biol Chem. 1993, 268: 16270-16278.PubMedGoogle Scholar
- Helbing CC, Werry K, Crump D, Domanski D, Veldhoen N, Bailey CM: Expression profiles of novel thyroid hormone-responsive genes and proteins in the tail of Xenopus laevis tadpoles undergoing precocious metamorphosis. Mol Endocrinol. 2003, 17: 1395-1409. 10.1210/me.2002-0274.PubMedGoogle Scholar
- Veldhoen N, Crump D, Werry K, Helbing CC: Distinctive gene profiles occur at key points during natural metamorphosis in the Xenopus laevis tadpole tail. Dev Dyn. 2002, 225: 457-468. 10.1002/dvdy.10175.PubMedGoogle Scholar
- Das B, Cai L, Carter MG, Piao YL, Sharov AA, Ko MS, Brown DD: Gene expression changes at metamorphosis induced by thyroid hormone in Xenopus laevis tadpoles. Dev Biol. 2006, 291: 342-355. 10.1016/j.ydbio.2005.12.032.PubMedGoogle Scholar
- Veldhoen N, Skirrow RC, Ji L, Domanski D, Bonfield ER, Bailey CM, Helbing CC: Use of heterologous cDNA arrays and organ culture in the detection of thyroid hormone-dependent responses in a sentinel frog, Rana catesbeiana. Comparative Biochemistry and Physiology Part D: Genomics and Proteomics. 2006, 1: 187-199. 10.1016/j.cbd.2005.10.005.Google Scholar
- Denver RJ, Pavhi S, Shi YB: Thyroid hormone-dependent gene expression program for Xenopus neural development. Journal of Biological Chemistry. 1997, 272: 8179-8188. 10.1074/jbc.272.13.8179.PubMedGoogle Scholar
- Tata JR: Requirement for RNA and protein synthesis for induced regression of the tadpole tail in organ culture. Dev Biol. 1966, 13: 77-94. 10.1016/0012-1606(66)90050-9.PubMedGoogle Scholar
- Lin HY, Davis FB, Gordinier JK, Martino LJ, Davis PJ: Thyroid hormone induces activation of mitogen-activated protein kinase in cultured cells. Am J Physiol. 1999, 276: C1014-24.PubMedGoogle Scholar
- Davis PJ, Davis FB, Cody V: Membrane receptors mediating thyroid hormone action. Trends Endocrinol Metab. 2005, 16: 429-435. 10.1016/j.tem.2005.09.007.PubMedGoogle Scholar
- Bassett JH, Harvey CB, Williams GR: Mechanisms of thyroid hormone receptor-specific nuclear and extra nuclear actions. Mol Cell Endocrinol. 2003, 213: 1-11. 10.1016/j.mce.2003.10.033.PubMedGoogle Scholar
- Skirrow RC, Helbing CC: Decreased cyclin dependent kinase activity promotes thyroid hormone-dependent tail regression in Rana catesbeiana. Cell and Tissue Research. 2007, 328: 281-289. 10.1007/s00441-006-0362-6.PubMedGoogle Scholar
- Ji L, Domanski D, Skirrow RC, Helbing CC: Genistein prevents thyroid hormone-dependent tail regression of Rana catesbeiana tadpoles by targetting protein kinase C and thyroid hormone receptor alpha. Developmental Dynamics. 2007, 236: 777-790. 10.1002/dvdy.21088.PubMedGoogle Scholar
- Ray LB, Dent JN: An analysis of the influence of thyroid hormone on the synthesis of proteins in the tail fin of bullfrog tadpoles. J Exp Zool. 1986, 240: 191-201. 10.1002/jez.1402400207.PubMedGoogle Scholar
- Kobayashi H, Sato H, Yoshizato K: Regionally regulated conversion of protein expression in the skin during anuran metamorphosis. J Exp Zool. 1996, 274: 187-192. 10.1002/(SICI)1097-010X(19960215)274:3<181::AID-JEZ5>3.0.CO;2-K.PubMedGoogle Scholar
- Steinberg TH, Agnew BJ, Gee KR, Leung WY, Goodman T, Schulenberg B, Hendrickson J, Beechem JM, Haugland RP, Patton WF: Global quantitative phosphoprotein analysis using Multiplexed Proteomics technology. Proteomics. 2003, 3: 1128-1144. 10.1002/pmic.200300434.PubMedGoogle Scholar
- Miyatani S, Winkles JA, Sargent TD, Dawid IB: Stage-specific keratins in Xenopus laevis embryos and tadpoles: the XK81 gene family. J Cell Biol. 1986, 103: 1957-1965. 10.1083/jcb.103.5.1957.PubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralPubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Research. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralPubMedGoogle Scholar
- Watanabe Y, Kobayashi H, Suzuki K, Kotani K, Yoshizato K: New epidermal keratin genes from Xenopus laevis: hormonal and regional regulation of their expression during anuran skin metamorphosis. Biochim Biophys Acta. 2001, 1517: 339-350.PubMedGoogle Scholar
- Watanabe Y, Tanaka R, Kobayashi H, Utoh R, Suzuki K, Obara M, Yoshizato K: Metamorphosis-dependent transcriptional regulation of xak-c, a novel Xenopus type I keratin gene. Dev Dyn. 2002, 225: 561-570. 10.1002/dvdy.10196.PubMedGoogle Scholar
- Suzuki K, Utoh R, Kotani K, Obara M, Yoshizato K: Lineage of anuran epidermal basal cells and their differentiation potential in relation to metamorphic skin remodeling. Dev Growth Differ. 2002, 44: 225-238. 10.1046/j.1440-169X.2002.00637.x.PubMedGoogle Scholar
- Ishida Y, Suzuki K, Utoh R, Obara M, Yoshizato K: Molecular identification of the skin transformation center of anuran larval skin using genes of Rana adult keratin (RAK) and SPARC as probes. Dev Growth Differ. 2003, 45: 515-526. 10.1111/j.1440-169X.2003.00719.x.PubMedGoogle Scholar
- Marceau N, Loranger A, Gilbert S, Daigle N, Champetier S: Keratin-mediated resistance to stress and apoptosis in simple epithelial cells in relation to health and disease. Biochem Cell Biol. 2001, 79: 543-555. 10.1139/bcb-79-5-543.PubMedGoogle Scholar
- Oshima RG: Apoptosis and keratin intermediate filaments. Cell Death Differ. 2002, 9: 486-492. 10.1038/sj.cdd.4400988.PubMedGoogle Scholar
- Kawai A, Ikeya J, Kinoshita T, Yoshizato K: A three-step mechanism of action of thyroid hormone and mesenchyme in metamorphic changes in anuran larval skin. Dev Biol. 1994, 166: 477-488. 10.1006/dbio.1994.1331.PubMedGoogle Scholar
- Suzuki K, Sato K, Katsu K, Hayashita H, Kristensen DB, Yoshizato K: Novel Rana keratin genes and their expression during larval to adult epidermal conversion in bullfrog tadpoles. Differentiation. 2001, 68: 44-54. 10.1046/j.1432-0436.2001.068001044.x.PubMedGoogle Scholar
- Taylor AC, Kollros JJ: Stages in the normal development of Rana pipiens larvae. Anat Rec. 1946, 7-23. 10.1002/ar.1090940103.Google Scholar
- Gilbert LI, Tata JR, Atkinson BG: Metamorphosis: postembryonic reprogramming of gene expression in amphibian and insect cells. Metamorphosis: postembryonic reprogramming of gene expression in amphibian and insect cells. 1996, San Diego, Academic PressGoogle Scholar
- Ku NO, Omary MB: Effect of mutation and phosphorylation of type I keratins on their caspase-mediated degradation. J Biol Chem. 2001, 276: 26792-26798. 10.1074/jbc.M103315200.PubMedGoogle Scholar
- Sachs LM, Abdallah B, Hassan A, Levi G, De Luze A, Reed JC, Demeneix BA: Apoptosis in Xenopus tadpole tail muscles involves Bax-dependent pathways. Faseb J. 1997, 11: 801-808.PubMedGoogle Scholar
- Sachs LM, Le Mevel S, Demeneix BA: Implication of bax in Xenopus laevis tail regression at metamorphosis. Dev Dyn. 2004, 231: 671-682. 10.1002/dvdy.20166.PubMedGoogle Scholar
- Rowe I, Coen L, Le Blay K, Le Mevel S, Demeneix BA: Autonomous regulation of muscle fibre fate during metamorphosis in Xenopus tropicalis. Dev Dyn. 2002, 224: 381-390. 10.1002/dvdy.10117.PubMedGoogle Scholar
- Schreiber AM, Brown DD: Tadpole skin dies autonomously in response to thyroid hormone at metamorphosis. Proc Natl Acad Sci U S A. 2003, 100: 1769-1774. 10.1073/pnas.252774999.PubMed CentralPubMedGoogle Scholar
- Omary MB, Ku NO, Liao J, Price D: Keratin modifications and solubility properties in epithelial cells and in vitro. Subcell Biochem. 1998, 31: 105-140.PubMedGoogle Scholar
- Ueno T, Toi M, Linder S: Detection of epithelial cell death in the body by cytokeratin 18 measurement. Biomed Pharmacother. 2005, 59 Suppl 2: S359-62. 10.1016/S0753-3322(05)80078-2.PubMedGoogle Scholar
- Linder S, Havelka AM, Ueno T, Shoshan MC: Determining tumor apoptosis and necrosis in patient serum using cytokeratin 18 as a biomarker. Cancer Lett. 2004, 214: 1-9. 10.1016/j.canlet.2004.06.032.PubMedGoogle Scholar
- Arunachalam B, Phan UT, Geuze HJ, Cresswell P: Enzymatic reduction of disulfide bonds in lysosomes: characterization of a gamma-interferon-inducible lysosomal thiol reductase (GILT). Proc Natl Acad Sci U S A. 2000, 97: 745-750. 10.1073/pnas.97.2.745.PubMed CentralPubMedGoogle Scholar
- Luster AD, Weinshank RL, Feinman R, Ravetch JV: Molecular and biochemical characterization of a novel gamma-interferon-inducible protein. J Biol Chem. 1988, 263: 12036-12043.PubMedGoogle Scholar
- Phan UT, Arunachalam B, Cresswell P: Gamma-interferon-inducible lysosomal thiol reductase (GILT). Maturation, activity, and mechanism of action. J Biol Chem. 2000, 275: 25907-25914. 10.1074/jbc.M003459200.PubMedGoogle Scholar
- Kornfeld S, Mellman I: The biogenesis of lysosomes. Annu Rev Cell Biol. 1989, 5: 483-525. 10.1146/annurev.cb.05.110189.002411.PubMedGoogle Scholar
- Honey K, Duff M, Beers C, Brissette WH, Elliott EA, Peters C, Maric M, Cresswell P, Rudensky A: Cathepsin S regulates the expression of cathepsin L and the turnover of gamma-interferon-inducible lysosomal thiol reductase in B lymphocytes. J Biol Chem. 2001, 276: 22573-22578. 10.1074/jbc.M101851200.PubMedGoogle Scholar
- Maric M, Arunachalam B, Phan UT, Dong C, Garrett WS, Cannon KS, Alfonso C, Karlsson L, Flavell RA, Cresswell P: Defective antigen processing in GILT-free mice. Science. 2001, 294: 1361-1365. 10.1126/science.1065500.PubMedGoogle Scholar
- Patton WF, Schulenberg B, Steinberg TH: Two-dimensional gel electrophoresis; better than a poke in the ICAT?. Curr Opin Biotechnol. 2002, 13: 321-328. 10.1016/S0958-1669(02)00333-6.PubMedGoogle Scholar
- Nesvizhskii AI, Aebersold R: Interpretation of shotgun proteomic data: the protein inference problem. Mol Cell Proteomics. 2005, 4: 1419-1440. 10.1074/mcp.R500012-MCP200.PubMedGoogle Scholar
- Ross PL, Huang YN, Marchese JN, Williamson B, Parker K, Hattan S, Khainovski N, Pillai S, Dey S, Daniels S, Purkayastha S, Juhasz P, Martin S, Bartlet-Jones M, He F, Jacobson A, Pappin DJ: Multiplexed protein quantitation in Saccharomyces cerevisiae using amine-reactive isobaric tagging reagents. Mol Cell Proteomics. 2004, 3: 1154-1169. 10.1074/mcp.M400129-MCP200.PubMedGoogle Scholar
- Perkins DN, Pappin DJ, Creasy DM, Cottrell JS: Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis. 1999, 20: 3551-3567. 10.1002/(SICI)1522-2683(19991201)20:18<3551::AID-ELPS3551>3.0.CO;2-2.PubMedGoogle Scholar
- Ma B, Zhang K, Hendrie C, Liang C, Li M, Doherty-Kirby A, Lajoie G: PEAKS: powerful software for peptide de novo sequencing by tandem mass spectrometry. Rapid Commun Mass Spectrom. 2003, 17: 2337-2342. 10.1002/rcm.1196.PubMedGoogle Scholar
- Chong PK, Gan CS, Pham TK, Wright PC: Isobaric tags for relative and absolute quantitation (iTRAQ) reproducibility: Implication of multiple injections. J Proteome Res. 2006, 5: 1232-1240. 10.1021/pr060018u.PubMedGoogle Scholar
- Nishikawa A, Murata E, Akita M, Kaneko K, Moriya O, Tomita M, Hayashi H: Roles of macrophages in programmed cell death and remodeling of tail and body muscle of Xenopus laevis during metamorphosis. Histochem Cell Biol. 1998, 109: 11-17. 10.1007/s004180050197.PubMedGoogle Scholar
- Rollins-Smith LA: Metamorphosis and the amphibian immune system. Immunol Rev. 1998, 166: 221-230. 10.1111/j.1600-065X.1998.tb01265.x.PubMedGoogle Scholar
- Izutsu Y, Tochinai S, Iwabuchi K, Onoe K: Larval antigen molecules recognized by adult immune cells of inbred Xenopus laevis: two pathways for recognition by adult splenic T cells. Dev Biol. 2000, 221: 365-374. 10.1006/dbio.2000.9681.PubMedGoogle Scholar
- Izutsu Y, Tochinai S, Maeno M, Iwabuchi K, Onoe K: Larval antigen molecules recognized by adult immune cells of inbred Xenopus laevis: partial characterization and implication in metamorphosis. Dev Growth Differ. 2002, 44: 477-488. 10.1046/j.1440-169X.2002.00660.x.PubMedGoogle Scholar
- Watanabe M, Ohshima M, Morohashi M, Maeno M, Izutsu Y: Ontogenic emergence and localization of larval skin antigen molecule recognized by adult T cells of Xenopus laevis: Regulation by thyroid hormone during metamorphosis. Dev Growth Differ. 2003, 45: 77-84. 10.1046/j.1440-169X.2003.00676.x.PubMedGoogle Scholar
- Seshimo H, Ryuzaki M, Yoshizato K: Specific inhibition of Triiodothyronine-induced tadpole tail-fin regression by cathepsin D-inhibitor pepstatin. Dev Biol. 1977, 59: 96-100. 10.1016/0012-1606(77)90244-5.PubMedGoogle Scholar
- Mukai M, Obara M, Yoshizato K: Characterization of gene of anuran cathepsin D as a metamorphosis-associated enzyme. Development, Growth and Differentiation. 1995, 463-477. 10.1046/j.1440-169X.1995.t01-4-00001.x.Google Scholar
- Bethune J, Wieland F, Moelleken J: COPI-mediated Transport. J Membr Biol. 2006, 211: 65-79. 10.1007/s00232-006-0859-7.PubMedGoogle Scholar
- Borth W: Alpha 2-macroglobulin, a multifunctional binding protein with targeting characteristics. Faseb J. 1992, 6: 3345-3353.PubMedGoogle Scholar
- Dodds AW, Law SK: The phylogeny and evolution of the thioester bond-containing proteins C3, C4 and alpha 2-macroglobulin. Immunol Rev. 1998, 166: 15-26. 10.1111/j.1600-065X.1998.tb01249.x.PubMedGoogle Scholar
- Lin KH, Lee HY, Shih CH, Yen CC, Chen SL, Yang RC, Wang CS: Plasma protein regulation by thyroid hormone. J Endocrinol. 2003, 179: 367-377. 10.1677/joe.0.1790367.PubMedGoogle Scholar
- Fries E, Kaczmarczyk A: Inter-alpha-inhibitor, hyaluronan and inflammation. Acta Biochim Pol. 2003, 50: 735-742.PubMedGoogle Scholar
- Josic D, Brown MK, Huang F, Lim YP, Rucevic M, Clifton JG, Hixson DC: Proteomic characterization of inter-alpha inhibitor proteins from human plasma. Proteomics. 2006, 6: 2874-2885. 10.1002/pmic.200500563.PubMedGoogle Scholar
- Pineiro M, Andres M, Iturralde M, Carmona S, Hirvonen J, Pyorala S, Heegaard PM, Tjornehoj K, Lampreave F, Pineiro A, Alava MA: ITIH4 (inter-alpha-trypsin inhibitor heavy chain 4) is a new acute-phase protein isolated from cattle during experimental infection. Infect Immun. 2004, 72: 3777-3782. 10.1128/IAI.72.7.3777-3782.2004.PubMed CentralPubMedGoogle Scholar
- Harraghy N, Mitchell TJ: Isolation and characterization of the promoter and partial enhancer region of the porcine inter-alpha-trypsin inhibitor heavy chain 4 gene. Clin Diagn Lab Immunol. 2005, 12: 1336-1339. 10.1128/CDLI.12.11.1336-1339.2005.PubMed CentralPubMedGoogle Scholar
- Pineiro M, Alava MA, Gonzalez-Ramon N, Osada J, Lasierra P, Larrad L, Pineiro A, Lampreave F: ITIH4 serum concentration increases during acute-phase processes in human patients and is up-regulated by interleukin-6 in hepatocarcinoma HepG2 cells. Biochem Biophys Res Commun. 1999, 263: 224-229. 10.1006/bbrc.1999.1349.PubMedGoogle Scholar
- Ricard-Blum S, Ruggiero F: The collagen superfamily: from the extracellular matrix to the cell membrane. Pathol Biol (Paris). 2005, 53: 430-442.Google Scholar
- Sawada T, Oofusa K, Yoshizato K: Thyroid hormone response element-like sequence in anuran matrix metalloproteinase 1 gene is responsive to in vivo thyroid hormone administration. J Endocrinol. 2001, 169: 477-486. 10.1677/joe.0.1690477.PubMedGoogle Scholar
- Oofusa K, Yomori S, Yoshizato K: Regionally and hormonally regulated expression of genes of collagen and collagenase in the anuran larval skin. Int J Dev Biol. 1994, 38: 345-350.PubMedGoogle Scholar
- Utoh R, Asahina K, Suzuki K, Kotani K, Obara M, Yoshizato K: Developmentally and regionally regulated participation of epidermal cells in the formation of collagen lamella of anuran tadpole skin. Dev Growth Differ. 2000, 42: 571-580. 10.1046/j.1440-169x.2000.00543.x.PubMedGoogle Scholar
- Asahina K, Oofusa K, Obara M, Yoshizato K: Cloning and characterization of the full length cDNA encoding alpha2 type I collagen of bullfrog Rana catesbeiana. Gene. 1997, 194: 283-289. 10.1016/S0378-1119(97)00210-2.PubMedGoogle Scholar
- Asahina K, Utoh R, Obara M, Yoshizato K: Cell-type specific and thyroid hormone-dependent expression of genes of alpha1(I) and alpha2(I) collagen in intestine during amphibian metamorphosis. Matrix Biol. 1999, 18: 89-103. 10.1016/S0945-053X(99)00005-0.PubMedGoogle Scholar
- Chen WJ, Lin KH, Lee YS: Molecular characterization of myocardial fibrosis during hypothyroidism: evidence for negative regulation of the pro-alpha1(I) collagen gene expression by thyroid hormone receptor. Mol Cell Endocrinol. 2000, 162: 45-55. 10.1016/S0303-7207(00)00203-3.PubMedGoogle Scholar
- Daly RJ: Cortactin signalling and dynamic actin networks. Biochem J. 2004, 382: 13-25. 10.1042/BJ20040737.PubMed CentralPubMedGoogle Scholar
- Eberle D, Hegarty B, Bossard P, Ferre P, Foufelle F: SREBP transcription factors: master regulators of lipid homeostasis. Biochimie. 2004, 86: 839-848. 10.1016/j.biochi.2004.09.018.PubMedGoogle Scholar
- Shin DJ, Osborne TF: Thyroid hormone regulation and cholesterol metabolism are connected through Sterol Regulatory Element-Binding Protein-2 (SREBP-2). J Biol Chem. 2003, 278: 34114-34118. 10.1074/jbc.M305417200.PubMedGoogle Scholar
- Hashimoto K, Yamada M, Matsumoto S, Monden T, Satoh T, Mori M: Mouse sterol response element binding protein-1c gene expression is negatively regulated by thyroid hormone. Endocrinology. 2006, 147: 4292-4302. 10.1210/en.2006-0116.PubMedGoogle Scholar
- Zhang Y, Yin L, Hillgartner FB: SREBP-1 integrates the actions of thyroid hormone, insulin, cAMP, and medium-chain fatty acids on ACCalpha transcription in hepatocytes. J Lipid Res. 2003, 44: 356-368. 10.1194/jlr.M200283-JLR200.PubMedGoogle Scholar
- Sullivan DT, MacIntyre R, Fuda N, Fiori J, Barrilla J, Ramizel L: Analysis of glycolytic enzyme co-localization in Drosophila flight muscle. J Exp Biol. 2003, 206: 2031-2038. 10.1242/jeb.00367.PubMedGoogle Scholar
- Smith-Gill SJ, Carver V: Biochemical characterization of organ differentiation and maturation. Metamorphosis: A problem in developmental biology 2nd Edition. Edited by: Gilbert LI and Frieden E. 1981, New York, Plenum Press, 491-544.Google Scholar
- Maruyama T, Watt KW, Riggs A: Hemoglobins of the tadpole of the bullfrog, Rana catesbeiana. Amino acid sequence of the alpha chain of a major component. J Biol Chem. 1980, 255: 3285-3293.PubMedGoogle Scholar
- Smith DJ, Zhu H, Kolatkar PR, Tam LT, Baldwin TO, Roe BA, Broyles RH, Riggs AF: The hemoglobins of the bullfrog, Rana catesbeiana. The cDNA-derived amino acid sequences of the alpha chains of adult hemoglobins B and C: their roles in deoxygenation-induced aggregation. J Biol Chem. 1993, 268: 26961-26971.PubMedGoogle Scholar
- Hasebe T, Oshima H, Kawamura K, Kikuyama S: Rapid and selective removal of larval erythrocytes from systemic circulation during metamorphosis of the bullfrog, Rana catesbeiana. Dev Growth Differ. 1999, 41: 639-643. 10.1046/j.1440-169x.1999.00461.x.PubMedGoogle Scholar
- Just JJ, Atkinson BG: Hemoglobin transitions in the bullfrog, Rana catesbeiana, during spontaneous and induced metamorphosis. Journal of Experimental Zoology. 1972, 182: 271-280. 10.1002/jez.1401820210.PubMedGoogle Scholar
- Watt KW, Maruyama T, Riggs A: Hemoglobins of the tadpole of the bullfrog, Rana catesbeiana. Amino acid sequence of the beta chain of a major component. J Biol Chem. 1980, 255: 3294-3301.PubMedGoogle Scholar
- Forman LJ, Just JJ: Cellular quantitation of hemoglobin transition during natural and thyroid-hormone-induced metamorphosis of the bullfrog, Rana catesbeiana. Gen Comp Endocrinol. 1981, 44: 1-12. 10.1016/0016-6480(81)90349-X.PubMedGoogle Scholar
- Yamaguchi T, Komoda Y, Nakajima H: Biliverdin-IX alpha reductase and biliverdin-IX beta reductase from human liver. Purification and characterization. J Biol Chem. 1994, 269: 24343-24348.PubMedGoogle Scholar
- Yamaguchi T, Nakajima H: Changes in the composition of bilirubin-IX isomers during human prenatal development. Eur J Biochem. 1995, 233: 467-472. 10.1111/j.1432-1033.1995.467_2.x.PubMedGoogle Scholar
- Hu H, Merrifield P, Atkinson BG: Expression of the myosin heavy chain genes in the tail muscle of thyroid hormone-induced metamorphosing Rana catesbeiana tadpoles. Dev Genet. 1999, 24: 151-164. 10.1002/(SICI)1520-6408(1999)24:1/2<151::AID-DVG14>3.0.CO;2-4.PubMedGoogle Scholar
- Boardman M, Cross GS, Jones EA, Woodland HR: Regulation of expression of a Xenopus borealis embryonic/larval alpha 3 skeletal-actin gene. Eur J Biochem. 1992, 208: 241-249. 10.1111/j.1432-1033.1992.tb17179.x.PubMedGoogle Scholar
- Collie ES, Muscat GE: The human skeletal alpha-actin promoter is regulated by thyroid hormone: identification of a thyroid hormone response element. Cell Growth Differ. 1992, 3: 31-42.PubMedGoogle Scholar
- Muscat GE, Griggs R, Downes M, Emery J: Characterization of the thyroid hormone response element in the skeletal alpha-actin gene: negative regulation of T3 receptor binding by the retinoid X receptor. Cell Growth Differ. 1993, 4: 269-279.PubMedGoogle Scholar
- Gustafson TA, Bahl JJ, Markham BE, Roeske WR, Morkin E: Hormonal regulation of myosin heavy chain and alpha-actin gene expression in cultured fetal rat heart myocytes. J Biol Chem. 1987, 262: 13316-13322.PubMedGoogle Scholar
- Lutz GJ, Lieber RL: Myosin isoforms in anuran skeletal muscle: their influence on contractile properties and in vivo muscle function. Microsc Res Tech. 2000, 50: 443-457. 10.1002/1097-0029(20000915)50:6<443::AID-JEMT3>3.0.CO;2-5.PubMedGoogle Scholar
- Miwa N, Kawamura S: Frog p26olf, a molecule with two S100-like regions in a single peptide. Microsc Res Tech. 2003, 60: 593-599. 10.1002/jemt.10301.PubMedGoogle Scholar
- Sakaguchi M, Miyazaki M, Takaishi M, Sakaguchi Y, Makino E, Kataoka N, Yamada H, Namba M, Huh NH: S100C/A11 is a key mediator of Ca(2+)-induced growth inhibition of human epidermal keratinocytes. J Cell Biol. 2003, 163: 825-835. 10.1083/jcb.200304017.PubMed CentralPubMedGoogle Scholar
- Makino E, Sakaguchi M, Iwatsuki K, Huh NH: Introduction of an N-terminal peptide of S100C/A11 into human cells induces apoptotic cell death. J Mol Med. 2004, 82: 612-620. 10.1007/s00109-004-0560-1.PubMedGoogle Scholar
- Kinter M, Sherman NE: The preparation of protein digests for mass spectrometric sequencing experiments. Protein sequencing and identification using tandem mass spectrometry. 2000, New York, Wiley-Interscience Inc, 152-160.Google Scholar
- Han Y, Ma B, Zhang K: SPIDER: software for protein identification from sequence tags with de novo sequencing error. J Bioinform Comput Biol. 2005, 3: 697-716. 10.1142/S0219720005001247.PubMedGoogle Scholar
- Wagner MJ, Helbing C: Multiple variants of the ING1 and ING2 tumor suppressors are differentially expressed and thyroid hormone-responsive in Xenopus laevis. General and Comparative Endocrinology. 2005, 144: 38-50. 10.1016/j.ygcen.2005.04.008.PubMedGoogle Scholar
- Veldhoen N, Skirrow RC, Osachoff H, Wigmore H, Clapson DJ, Gunderson MP, van Aggelen G, Helbing CC: The bactericidal agent Triclosan modulates thyroid hormone-associated gene expression and disrupts postembryonic anuran development. Aquatic Toxicology. 2006, 80: 217-227. 10.1016/j.aquatox.2006.08.010.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.