ToxSci Advance Access originally published online on January 27, 2006
Toxicological Sciences 2006 92(1):254-269; doi:10.1093/toxsci/kfj118
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
HIGHLIGHTED ARTICLE |
Regenerative Growth Is Impacted by TCDD: Gene Expression Analysis Reveals Extracellular Matrix Modulation
Department of Environmental and Molecular Toxicology, Environmental Health Sciences Center, Marine and Freshwater Biomedical Sciences Center, Oregon State University, Corvallis, Oregon 97331
1 To whom correspondence should be addressed at Oregon State University, Department of Environmental and Molecular Toxicology, 1007 ALS, Corvallis, Oregon 97331-7301. Fax: (541) 737-7966. E-mail: robert.tanguay{at}oregonstate.edu.
Received November 28, 2005; accepted January 20, 2006
| ABSTRACT |
|---|
|
|
|---|
Adult zebrafish can completely regenerate their caudal fin following amputation. This complex process is initiated by the formation of an epithelial wound cap over the amputation site by 12 h post amputation (hpa). Once the cap is formed, mesenchymal cells proliferate and migrate from sites distal to the wound plane and accumulate under the epithelial cap forming the blastemal structure within 48 hpa. Blastemal cells proliferate and differentiate, replacing the amputated tissues, which are populated with angiogenic vessels and innervating nerves during the regenerative outgrowth phase which is completed around 14 days post amputation (dpa). Regenerative outgrowth does not occur in TCDD-exposed zebrafish. To identify the molecular pathways that are perturbed by TCDD exposure, male zebrafish were ip injected with 50 ng/g TCDD or vehicle and caudal fins were amputated. Regenerating fin tissue was collected at 1, 3, and 5 dpa for mRNA abundance analysis. Microarray analysis and quantitative real time PCR revealed that wound healing and regeneration alone altered the expression of nearly 900 genes by at least two-fold between 1 and 5 dpa. TCDD altered the abundance of 370 genes at least two-fold. Among these, several known aryl hydrocarbon responsive genes were identified in addition to several genes involved in extracellular matrix composition and metabolism. The profile of misexpressed genes is suggestive of impaired cellular differentiation and extracellular matrix composition potentially regulated by Sox9b.
Key Words: AHR; TCDD; regeneration; zebrafish; extracellular matrix.
| INTRODUCTION |
|---|
|
|
|---|
Adult zebrafish have the capacity to regenerate their caudal fin upon amputation, a process that is inhibited by exposure to the ubiquitous environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) (Zodrow and Tanguay, 2003
The AHR is a ligand activated, basic helix-loop-helix transcription factor and is a member of the PAS domain family of genes (reviewed in Gu et al., 2000
). The AHR signal transduction pathway in fish is similar to that in mammals, with the exception that fish have two or more AHR genes while mammals have only one (reviewed in Hahn et al., 1997
; Tanguay et al., 2003
). Similar to mammals, two dimerization partners for AHR, AHR nuclear translocator proteins 1 and 2 (ARNT1, ARNT2), have been cloned and characterized in zebrafish. AHR2 and ARNT1 are the members of the AHR pathway in zebrafish that mediate TCDD toxicity including TCDD-impaired larval fin regeneration (Andreasen et al., 2002
; Mathew et al., 2005
; Prasch et al., 2003
, 2004
, in press
). Therefore, despite the added complexity of the fish AHR phylogeny, the functionality of the fish and mammalian AHR pathways are similar.
Amputation and regeneration of the adult zebrafish caudal fin has been used as a system to study epimorphic tissue regeneration (reviewed in Poss et al., 2003
). Epimorphic regeneration involves the reprogramming and migration of cells that differentiate and restore a tissue to its original form. The caudal fin consists of a series of boney rays attached to the skeleton by ligaments. Each ray, referred to as a lepidotrichia is comprised of two hemirays. The hemirays consist of repeating segments joined by ligaments. The area between the hemirays of a lepidotrichia contains blood vessels, nerves, and mesenchymal cells. The fin grows by successive addition of hemiray segments to the most distal (end of fin) segment. After partial amputation of the caudal fin, an apical epithelial wound cap (AEC) forms over the clot at the amputation site within 12 h post amputation (hpa). The AEC cells are derived from non-proliferating local epithelial cells that migrate laterally over the wound. Once the cap is formed, mesenchymal cells proliferate and migrate from sites distal to the wound plane and accumulate under the epithelial cap forming a structure called the blastema within 48 hpa. The cells forming the blastema are believed to be de-differentiated cells that are later reprogrammed to form the regenerating tissue. Scleroblast (bone forming cells giving rise to the hemiray) and blastemal cells proliferate and differentiate, replacing the amputated tissues, which are infused with nerves and blood vessels in a process called regenerative outgrowth. Outgrowth begins at 2 days post amputation (dpa) until the regenerative event is complete at about 14 dpa.
Zebrafish caudal fin regeneration is impaired by TCDD when exposure occurs at any time during the first four days of regeneration (Zodrow and Tanguay, 2003
). An epithelial wound cap and blastema forms in the TCDD exposed regenerate, but outgrowth does not occur. The AHR pathway is active in the regenerating fin since both zfAHR2 and TCDD induced CYP1A mRNAs (a marker of an active AHR pathway) are expressed in the regenerating caudal fin (Andreasen et al., 2002
; Zodrow and Tanguay, 2003
). Several studies in zebrafish have proposed potential mechanisms that initiate AEC, blastema formation, and outgrowth (reviewed in Poss et al., 2003
). The effect of TCDD on the expression of these and other genes and pathways involved in fin regeneration was the focus of this work. Microarray analysis revealed that 370 known genes were differentially regulated by TCDD during the first five days of regeneration. The major categories of responsive genes were detoxification enzymes, extracellular matrix, development, retinoic acid signaling, and metabolism. These studies have laid the foundation for the dissection of the pathways required for fin regeneration and those that are modified by TCDD.
| MATERIALS AND METHODS |
|---|
|
|
|---|
TCDD dosing and fin amputation.
Adult male zebrafish (Danio rerio) (AB strain, Eugene, OR) were anesthetized with tricaine methanesulfonate (MS-222) and their caudal fins were amputated at the bifurcation of the fin rays. Fish were then injected with phosphatidylcholine liposomes (vehicle control) alone or liposomes incorporated with TCDD (>99% pure, Chemsyn, Lenexa, KS) prepared as previously described (Andreasen et al., 2002
Total RNA isolation.
Fins of vehicle control or TCDD exposed fish were isolated as described above. RNAlater was removed from the samples and total RNA was purified with TRI reagent (Molecular Research Laboratories, Cincinnati, OH) according to the manufacturer's instructions. Total RNA was DNase treated with RQ1 DNase (Promega, Madison, WI) according to the manufacturer's protocol and RNA quantity and quality was determined by UV absorbance. Ribosomal RNA abundance and degree of degradation was determined in electropherogram patterns using the 2100 Bioanalyzer and RNA 6000 Nano chips (Agilent Technologies, Palo Alto, CA).
Affymetrix microarray processing.
Probe synthesis, hybridization, and scanning were conducted by the Center for Gene Research at Oregon State University, Corvallis OR. 2.5 µg of total RNA was used to generate biotinylated complementary RNA (cRNA) for each treatment group using the One-Cycle Target Labeling kit (Affymetrix, Santa Clara, CA). Briefly, RNA was reverse transcribed using a T7-(dT)24 primer and Superscript II reverse transcriptase (Invitrogen, Carlsbad, CA) and double stranded cDNA was synthesized and purified. Biotinylated cRNA was synthesized from the double stranded cDNA using T7 RNA polymerase and a biotin-conjugated pseudouridine containing nucleotide mixture provided in the IVT Labeling Kit (Affymetrix, Santa Clara, CA). Prior to hybridization, the cRNA was purified, fragmented and 10 µg from each experimental sample was hybridized to zebrafish genome arrays (Zebrafish430_2) according to the Affymetrix GeneChip Expression Analysis Technical Manual (701021 Rev. 5). Affymetrix fluidics station 400 was used to wash the arrays. Arrays were scanned with an Affymetrix scanner 3000. The Affymetrix eukaryotic hybridization control kit and Poly-A RNA control kit were used to ensure efficiency of hybridization and cDNA amplification. All cRNA was synthesized at the same time. Hybridizations were conducted with one replicate of all times and treatments concurrently. Each array image was visualized to discount artifactual signals, scratches or debris. Experiments were MIAME certified. The raw data is listed at the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/projects/geo/); series record #GSE3667).
The Affymetrix zebrafish genome array was designed to detect over 14,900 transcripts accounting for approximately 33% of the predicted genome. Probe sets comprising 16 different oligonucleotides were designed to detect an individual transcript. Since the zebrafish genome has not been completely sequenced, the microarray may have different probe sets that detect the same transcript. Additionally, annotation of the microarray is incomplete. Sequence used to design the probe sets were found in dbEST (July 2003), GenBank (release 136.0, June 2003), RefSeq (July 2003), and UniGene (Build 54, June 2003). Sequence similarity to known mammalian proteins was determined by identifying full length mRNA sequence for each zebrafish probe set by conducting a BLAST search of each Affymetrix probe set against Genbank (http://www.ncbi.nlm.nih.gov/BLAST/), TIGR (http://tigrblast.tigr.org/tgi/), and Sanger (http://www.sanger.ac.uk/Projects/D_rerio/) databases. The predicted amino acid sequences from these genes were then compared to mammalian orthologs using pblast at NCBI. The deduced full-length protein predicted from each transcript's probe set was confirmed to be at least 45% identical to its mammalian ortholog.
qRT-PCR validation of mRNA abundance.
Total RNA was isolated in triplicate from the regenerating fins of vehicle control or TCDD exposed fish (n = 6/group). cDNA was prepared from 2 mg of total RNA per group using Superscript II (Life Technologies, Gaithersburg, MD) and oligo(dT) primers in a 20 ml volume. Quantitative PCR using gene specific primers was conducted using the Opticon 2 real-time PCR detection system (MJ Research, Waltham, MA). Specifically, 1 ml of each cDNA pool was used for each PCR reaction in the presence of SYBR Green, using DyNAmo SYBR green qPCR kit according to the manufacturer's instructions (Finnzymes, Espoo Finland). All experimental samples were analyzed concurrently with standard curves generated from known quantities of DNA template covering seven orders of magnitude for each gene in duplicate. Negative controls for each experiment consisted of RNA without reverse transcriptase and another without template. Agarose gel electrophoresis and thermal denaturation (melt curve analysis) were conducted to ensure formation of specific products.
Oligonucleotides.
All oligonucleotide primers were synthesized by MWG-Biotech (High Point, NC). Forward primers are prefixed with an F corresponding to sense strands and antisense reverse primers are designated with an R. Primers were designed to amplify sequence within the Affymetrix probe set sequence. Sequence for each primer can be found in Table 1.
|
Statistical analysis.
Affymetrix cell files were imported into GeneSpring 7.1 software (Agilent Technologies, Palo Alto, CA) and GC-RMA processed to correct for background signal and to normalize across arrays. Subsequently, each transcript was normalized to the median signal for each transcript across all arrays to allow for between array comparisons on a relative scale for each gene. All statistical analysis of microarray data was conducted using GeneSpring 7.1 software.
The effects of time on transcript level of vehicle control exposed animals was determined by one-way ANOVA assuming equal variance employing Benjamini and Hochberg multiple testing corrections (p < 0.05). Only genes at least two-fold differentially expressed from 1 dpa vehicle control transcript levels were considered for this analysis. Hierarchical clustering using a smooth correlation and default parameters was employed to characterize expression patterns during regeneration in vehicle exposed animals.
In order to determine the effect of TCDD on gene expression, the data set was filtered for significance by first generating a list of genes that were at least two-fold differentially expressed in comparison to 1 dpa control fins in at least one time point independent of TCDD exposure. Effects of time and treatment were first analyzed by two-way ANOVA analysis assuming equal variance and employing Benjamini and Hochberg multiple testing corrections (p < 0.05). The effects of TCDD dosing was determined by filtering control and TCDD data sets for an individual time point for genes that were at least two fold differentially expressed between control and TCDD dosed fish. Statistical significance of TCDD exposure was determined by one-way ANOVA comparing time matched controls to TCDD dosed transcript levels assuming equal variance (p < 0.05).
Significant differences of mRNA abundance determined with qRT-PCR were assessed by two-way ANOVA on log10 transformed data. Differences between groups were analyzed by the Tukey method (p < 0.05). qRT-PCR analyses were performed using SigmaStat software (Chicago, IL).
| RESULTS |
|---|
|
|
|---|
Differential Gene Expression during Fin Regeneration
Changes in gene expression were monitored at 1, 3, and 5 dpa in regenerating tissue from adult male zebrafish exposed at the time of amputation to vehicle control or 50 ng/g TCDD using Affymetrix zebrafish genome arrays. These time points were chosen because TCDD impairs regenerative fin outgrowth, a process that begins at 2 dpa. The oligonucleotide array was designed to detect approximately 14,900 transcripts accounting for approximately 33% of the estimated genome. Although each oligonucleotide (probe) on the array was designed to detect a unique transcript, the degree of redundancy is not completely known. The dose of TCDD chosen is approximately the EC90 for impairing fin regeneration. Importantly, this dose does not lead to overt toxicity such as fin necrosis, hemorrhaging wasting, or alterations in swimming behavior and feeding. Tissues isolated from 10 fish at each time and treatment, were pooled to generate an experimental group. Three independent biological replicates of each time and treatment were analyzed. A total of 18 microarrays were analyzed.
Approximately 11,700 transcripts were detected after filtering out low abundant genes. Among these, 1383 genes were differentially affected at least two fold by time and or TCDD exposure in at least one time point (Fig. 1). Wound healing and regeneration alone altered the expression of 932 genes by at least two-fold between 1 and 5 dpa (Fig. 1). Expression of 99 genes was affected by TCDD alone (Fig. 1). Time and TCDD exposure affected the abundance of 352 genes (Fig. 1). Further analysis of the change of expression of genes over time using one way ANOVA analysis of vehicle exposed animals identified 898 genes that were at least two-fold differentially regulated from 1 to 5 dpa (Fig. 2). Hierarchical clustering analysis was employed to characterize the expression patterns of these vehicle exposed animals. Five major clusters were identified and labeled AE. The ratio of genes with enhanced to repressed expression was 3 to 2. Since the Affymetrix array is not fully annotated, and the effect of TCDD was the focus of these studies, the identity of these non-annotated genes remains to be determined.
|
|
Altered Gene Expression in Response to TCDD Exposure
Gene expression altered by TCDD exposure was determined by first filtering for genes that were at least 2 fold differentially regulated from their time matched control. This gene list was then analyzed for statistical significance employing one-way ANOVA analysis comparing time matched vehicle and TCDD exposed gene profiles. The expression of 370 genes was differentially regulated by TCDD exposure (Fig. 3). The temporal response increased over time with 50, 196, and 265 being differentially expressed at 1, 3, and 5 dpa. This increasing diversity of gene expression parallels the timeline of the pathological signs of toxicity. TCDD exposure resulted in enhancing the abundance of 127 genes of which 67 were identified with mammalian orthologs (Table 2). Genes were categorized into functional groups based upon sequence homology and known function in the literature. Expression of eight genes involved in xenobiotic metabolism including genes known to be inducible by AHR activation such as cytochrome P4501A were affected by TCDD confirming that the tissues were exposed to TCDD and had an active AHR pathway as early as 1 dpa. Other major functional groups with enhanced gene expression include extracellular matrix (ECM) and metabolism, development, FGF signaling, and metabolism. TCDD exposure resulted in reduced expression of 243 genes of which 132 were identified and categorized in Table 3. Two major functional groups were identified. Thirty-six genes with reduced expression are components in ECM or metabolism of this matrix. The expression of 21 genes involved in development was negatively affected by TCDD. Other major groups of genes that were negatively affected by TCDD include retinoic acid metabolism and signaling, wound response, metabolism, and amino acid metabolism. Lack of sufficient sequence information inhibited the identification of 60 induced and 113 repressed genes.
|
|
|
Validation of Gene Expression Profiles
Separate experiments using quantitative real time PCR were conducted to measure the abundance of several of the transcripts listed in Tables 2 and 3, representing high and low responders in distinct functional categories. Separate biological samples from those used in the microarray experiments were collected to ensure that significant responses to TCDD were not experimental artifacts. Animals were dosed with vehicle or TCDD similarly to the microarray experiment and three replicates (six fins per replicate) were each isolated at 1, 3, and 5 dpa. Abundance of seven transcripts that were enhanced by TCDD exposure and four that were repressed by TCDD exposure were analyzed.
Transcripts of the AHR inducible xenobiotic detoxifying enzymes CYP1A, CYP1B1, UDPGT in vehicle exposed fish were initially decreased at 1 dpa in comparison to non-regenerating fins (Figs. 4A, 4B, 4C). However their abundance became elevated by 3 dpa and they were all significantly enhanced by TCDD dosing from 1 to 5 dpa. The low abundance of CYP1A in vehicle control animals lead to an artificially high fold induction (low denominator) and therefore was not included in Figure 4, but is listed in Table 4. Representative genes grouped into the ECM metabolism, immune and development categories also responded to TCDD. Frizzled 7a and Cadherin 19 had elevated expression in response to TCDD as early as 1 dpa and remained elevated to 5 dpa (Figs. 4D and 4E). BAMBI and CC-Chemokine 1 had transient response to TCDD with BAMBI displaying TCDD inducible expression at 3 and 5 dpa and CC-Chemokine 1 not responding to TCDD exposure until 5 dpa (Figs. 4F and 4G). Frizzled 7a, BAMBI and CC-Chemokine 1 message level in vehicle exposed fish was not repressed by regeneration alone unlike the xenobiotic metabolizing enzymes.
|
|
Four genes identified as reduced by TCDD exposure on the microarrays were quantified to determine the congruence between gene array and qRT-PCR. All four of these genes had enhanced expression during regeneration in comparison to the non-regenerating fin (Fig. 5). Sox9b was chosen because it was highly responsive to TCDD by microarray analysis. Sox9b expression was reduced by TCDD at all time points between 3 and 22 fold (Fig. 5A). Three genes that were marginally, but significantly repressed by TCDD were quantified by qRT-PCR to asses how well microarrays discern small changes in gene expression. Two genes involved in collagen maturation, lysyl oxidase 2 and prolyl 4-hydroxylase
1, each increased expression during the first five days of regeneration (Figs. 5B and 5C). Additionally both were repressed by TCDD with lysyl oxidase 2 being repressed at all time points and prolyl 4-hydroxylase
1 at 3 and 5 dpa. Procollagen c-endopeptidase enhancer 2 (PCOLCE2), another collagen modifying gene was significantly reduced by TCDD only at 1 dpa and the basal expression like the two other collagen metabolizing enzymes increased over time (Fig. 5D).
|
Comparison of mRNA abundance determined by microarray and qRT-PCR revealed similar results especially when the response to TCDD was robust (Table 4). In order to make comparisons between the two techniques, data is displayed as fold change relative to vehicle control tissues at the same time in regeneration. The expression of all genes responded to TCDD in a similar manner whether analyzed with qRT-PCR or microarray. However there were some discrepancies in response time with some genes that expressed less than a three-fold change in response to TCDD. For instance, both methods revealed significant enhancement of CC-Chemokine and BAMBI gene expression in response to TCDD exposure using microarray analysis, however they were not responsive when using qRT-PCR at 1 dpa (Table 4). Similarly, procollagen endopeptidase enhancer 2 (PCOLCE2) expression was repressed by TCDD exposure at all times investigated using microarray analysis. This gene product was only repressed 2.22.8 fold at 3 and 5 dpa, respectively, when analyzed using microarray analysis and was not found to be repressed by qRT-PCR at these time points (Table 5). Discrepancies in abundance of a transcript in Table 4 may be accounted for by differences in method, but may also be partially due to biological variability between the different biological pools of RNA that were used to conduct the different techniques.
|
| DISCUSSION |
|---|
|
|
|---|
Activation of the AHR Pathway
Fin regeneration is a complex process requiring precise regulation of several processes including wound healing, ECM production, revascularization, innervation, and bone formation. Due to the complex morphological changes during regeneration, it is not surprising that enormous transcriptional changes are required to complete this process. Since TCDD inhibits regeneration, the identification of the genes altered by TCDD at the site of toxicity over time is crucial to increase our understanding of the mechanism of toxicity. We know that this regenerative toxicity in the fin is dependent upon AHR2 and ARNT1 signaling, since toxicity is averted if either of these genes are repressed in larval zebrafish during fin regeneration (Mathew et al., 2005
Extracellular Matrix Related Transcripts
The largest category of genes misregulated by TCDD are those related to ECM composition and metabolism. TCDD exposure resulted in repression of 34 of the 41 genes in this category. This suggests that TCDD impairs the maturation of the ECM. This could occur by directly altering the expression of ECM genes or by impacting genes responsible for controlling cell differentiation and matrix maturation. Structural differences in ECM can occur due to impaired synthesis or catabolism of ECM components including collagen. This is noteworthy as pharmacological inhibition of matrix metalloproteinases (MMPs), proteins involved in metabolism of ECM constituents including collagen, inhibit blastemal proliferation and fin regeneration (Bai et al., 2005
). Moreover, proper regulation of the plasminogen and MMP pathway is known to regulate ECM composition and wound healing (reviewed in Lund et al., 1999
). Additionally, numerous drugs known to affect metabolism of collagen, one of the principal constituents of ECM, have been shown to impair collagen organization and fin regeneration (Bechara et al., 2000
). Evidence here suggests that both post-translational modifications and metabolism of collagen is impaired by TCDD.
Collagens function to provide tensile strength to tissues and are also involved in cell adhesion, chemotaxis, migration, differentiation, and wound healing. In addition to catabolism by MMPs, proper collagen structure can be compromised by altering synthesis, trafficking, post-translational modification, and fibril assembly reviewed below (reviewed in Canty and Kadler, 2005
). In addition to the repression of procollagen 1a1 and 11a1 message by TCDD, the abundance of 11 transcripts involved in collagen synthesis and catabolism were impacted. Abundance of 9 transcripts involved in post translational modification of collagen were repressed by TCDD in the fin. Prolyl-4-hydroxylase
1 and
2 (P4H) were both repressed beginning at 1 and 3 dpa respectively. P4H functions to hydroxylate proline residues of collagen in the ER thereby stabilizing the triple helical structure and increasing the denaturing temperature. The observed 3 dpa reduction in lysyl hydroxylase would be expected to impair collagen processing in the ER since it serves to generate sites of carbohydrate attachment and helps to stabilize collagen crosslinks. Impaired hydroxylation of procollagen can lead to its retention in the ER. TCDD reduced the mRNA of three genes, HSP47 and the FK506 binding proteins 9 and 14, involved in proper folding of procollagen chains in the ER, beginning at 3 dpa. FK506 binding proteins 9 and 14 contain peptidylproline cis-trans isomerase activity which converts proline to the trans conformation allowing proper protein folding. HSP47, repressed by TCDD at 3 dpa, acts as a chaperone involved in folding procollagen in the ER. Fibrillar collagens (types IIII, V, and XI) contain N and C-terminal propeptides that are removed by procollagen peptidases allowing the mature collagen to form fibrils in the ECM. Procollagen C-endopeptidase enhancers 2 and 1 (PCOLCE 2 and 1), repressed by 1 and 3 dpa respectively bind to the C-terminal propeptide enhancing the cleavage of the propeptide by proteases. Collagen fibril assembly may be misregulated since lysyl oxidases (2 and 1), repressed at 1 and 3 dpa, catalyze the first step in cross linking of collagen fibrils in the ECM.
The potential for alteration in ECM turnover by TCDD also exists since MMP13 (collagenase-3) message was greatly enhanced by TCDD exposure throughout the five days of regeneration. This is not unique to zebrafish since TCDD exposure results in elevated collagenase (MMP1 or 13) in cell culture and in the mouse heart (Partridge et al., 2000
; Thackaberry et al., 2005
; Villano et al., 2006
; Vogel et al., 2004
). TIMP2 (tissue inhibitor of metalloproteinase 2, like) was elevated at 5 dpa. TIMP2 inhibits MMP2 activation. In addition to serving as structural architecture, the ECM also serves as a source and store of biologically active molecules which can be released and or activated by MMPs affecting cell growth, cell migration, cell-cell communication and apoptosis, angiogenesis (reviewed in Mott and Werb, 2004
). By this mechanism, alterations in collagen abundance and structure could also affect availability of extracellular signaling molecules such as bFGF, TGF-ß, and VEGF.
Misexpression of a regulator of ECM production or cell differentiation fate could explain the TCDD-dependent changes in ECM transcript abundance. The HMG-box DNA-binding transcription factor SRY-box 9 (Sox9) may play this role in the regenerating fin. A functional role of Sox9 was first discovered in humans with heterozygous mutations in and around the Sox9 gene resulting in skeletal malformations known as campomelic dysplasia (CD), that is often associated with XY sex reversal and other malformations of internal organs (reviewed in Schafer et al., 1996
). Sox9 is considered essential for commitment of mesenchymal cells to chondrocyte differentiation and cartilage formation (Bi et al., 1999
; Mori-Akiyama et al., 2003
). ECM constituents enhanced by Sox9 during chondrogensis include procol2a1, procol9a1, procol11a2, CD Rap, aggrecan, cartilage link protein (reviewed in de Crombrugghe et al., 2000
). In zebrafish there are two Sox9 orthologs, Sox9a and Sox9b, which have distinct and overlapping expression patterns during development which additively are thought to duplicate the mammalian Sox9 function (Chiang et al., 2001
; Yan et al., 2005
). Both Sox9a and 9b bind specific DNA response elements and have similar transactivation properties (Chiang et al., 2001
). Studies using Sox9 mutant variant zebrafish and antisense translational repression determined that both Sox9 genes play an essential role in chondrogenesis with Sox9b involved in attaining proper numbers of chondrocytes (Yan et al., 2005
). This is striking because during development, TCDD exposure impairs chondrogenesis in the head in an AHR2 dependent manner (Prasch et al., 2003
; Teraoka et al., 2002
). The TCDD dependent morphology remarkably resembles Sox9b morphants (antisense repressed) suggesting that Sox9b could be impaired by TCDD during development (Yan et al., 2005
). It may be that Sox9b is required to lay down a provisional matrix in the regenerating fin giving it form, while the tissue is regenerating. Sox9b message was indeed enhanced by amputation 48 fold during the first five days of regeneration in comparison to non-regenerating fin tissue. However Sox9b was the most strongly reduced transcript by TCDD during the first five days of regeneration, being repressed to levels similar to non-regenerating fins. Whether Sox9b is directly regulated by AHR2 remains to be determined.
Evidence of functional repression of Sox9 can be illustrated by reduced expression of Sox9 dependent transcripts and repressed chondrogenic markers. This was observed in the fin as cartilage link protein (CRTL) the well-characterized Sox9 dependent gene (Kou and Ikegawa, 2004
) was the second strongest repressed gene by TCDD in the regenerates. CRTL plays a crucial role in chondrogenesis by stabilizing proteoglycan aggregates, giving cartilage its tensile strength, elasticity, and organization of chondrocytes. Targeted disruption of CRTL1 is lethal in mice and results in severe chondrodysplasia (Czipri et al., 2003
). Of additional importance, once CRTL is proteolytically cleaved, the N-terminus can stimulate proteoglycan and type II collagen synthesis contributing to the ECM (Liu et al., 2000
). Another indicator of impaired chondrogenesis is the reduced expression of the regulator of G-protein signaling 2 (RGS-2) at 5 dpa. RGS-2 was recently found to promote late stage chondrogenesis (James et al., 2005
).
Transcriptional regulation of Sox9 is not well understood. A potential modifier of Sox9b expression in response to TCDD is the retinoic acid pathway. Sox9 expression responds to retinoic acid in a cell or tissue specific manner, either with enhanced or repressed expression (Afonja et al., 2002
; Sekiya et al., 2000
, 2001
). TCDD is known to alter retinoid storage and metabolism of retinoic acid in liver, kidney, and plasma (Fletcher et al., 2005a
). Transcripts of several retinoic acid responsive genes including cyp26a1, cyp26b1, and short-chain dehydrogenase/reductase 3, each involved in retinoid metabolism, were repressed by TCDD, although only at 3 and 5 dpa. If a repressed retinoid pathway leads to reduced Sox9b expression as early as 1 dpa, then other RA responsive genes should be similarly down regulated. This remains to be determined. However the potential for RA regulation of Sox9b exists since the Sox9b promoter has 10 RAR elements within 4000 bp of the translation start site. Additionally Sox9a has no RAR response elements and only one RXR element suggesting that Sox9b would be preferentially regulated by the RA signaling.
Potential rationale for and consequences of reduced Sox9b expression is illustrated in Figure 6. AHR could alter Sox9 signaling by impairing its expression directly through RA signaling or other pathways such as FGF and BMP signaling known to regulate fin regeneration and Sox9 expression. AHR may interfere with Sox9b expression at its enhancer-promoter directly or through other negative regulators such as AHRR or alter Sox9b message stability. Ultimately the observed alteration in ECM constituents could be directed through AHR via Sox9b or through other AHR mediated mechanisms. The functional result of repressed Sox9b would be impaired EMC production and processing, cell adhesion, and bone formation.
|
Taken together, these results demonstrate that the regenerating fin is a dynamic tissue characterized by enormous alterations in transcript abundance within five days of amputation. The response to TCDD is drastic with complex changes in the transcriptome over time. The most striking functional category of altered transcripts are those involved in ECM composition or its regulation. It will be essential to determine the role of these responsive genes using the power of the zebrafish model, and to determine if similar TCDD responses occur in other regenerative and developmental models.
| ACKNOWLEDGMENTS |
|---|
We would like to thank Anne-Marie Girard from the Center for Gene Research and Biotechnology at Oregon State University for her valuable assistance. We would also like to thank Eric Johnson and Greg Gonnerman for their technical assistance. This work was supported in part by NIEHS grants ES10820, ES00210, and ES03850.
| REFERENCES |
|---|
|
|
|---|
Afonja, O., Raaka, B. M., Huang, A., Das, S., Zhao, X., Helmer, E., Juste, D., and Samuels, H. H. (2002). RAR agonists stimulate SOX9 gene expression in breast cancer cell lines: Evidence for a role in retinoid-mediated growth inhibition. Oncogene 21, 78507860.[CrossRef][Web of Science][Medline]
Andreasen, E. A., Hahn, M. E., Heideman, W., Peterson, R. E., and Tanguay, R. L. (2002). The zebrafish (Danio rerio) aryl hydrocarbon receptor type 1 (zfAHR1) is a novel vertebrate receptor. Mol. Pharmacol. 62, 234249.
Bai, S., Thummel, R., Godwin, A. R., Nagase, H., Itoh, Y., Li, L., Evans, R., McDermott, J., Seiki, M., and Sarras, M. P., Jr. (2005). Matrix metalloproteinase expression and function during fin regeneration in zebrafish: analysis of MT1-MMP, MMP2 and TIMP2. Matrix Biol. 24, 247260.[CrossRef][Web of Science][Medline]
Bechara, I. J., Joazeiro, P. P., Mari-Beffa, M., Becerra, J., and Montes, G. S. (2000). Collagen-affecting drugs impair regeneration of teleost tail fins. J. Submicrosc. Cytol. Pathol. 32, 273280.[Medline]
Bi, W., Deng, J. M., Zhang, Z., Behringer, R. R., and de Crombrugghe, B. (1999). Sox9 is required for cartilage formation. Nat. Genet. 22, 8589.[CrossRef][Web of Science][Medline]
Birnbaum, L. S., and Tuomisto, J. (2000). Non-carcinogenic effects of TCDD in animals. Food Addit. Contam. 17, 275288.[CrossRef][Web of Science][Medline]
Boverhof, D. R., Burgoon, L. D., Tashiro, C., Chittim, B., Harkema, J. R., Jump, D. B., and Zacharewski, T. R. (2005). Temporal and dose-dependent hepatic gene expression patterns in mice provide new insights into TCDD-mediated hepatotoxicity. Toxicol. Sci. 85, 10481063.
Canty, E. G., and Kadler, K. E. (2005). Procollagen trafficking, processing and fibrillogenesis. J. Cell Sci. 118, 13411353.
Carney, S. A., Peterson, R. E., and Heideman, W. (2004). 2,3,7,8-Tetrachlorodibenzo-p-dioxin activation of the aryl hydrocarbon receptor/aryl hydrocarbon receptor nuclear translocator pathway causes developmental toxicity through a CYP1A-independent mechanism in zebrafish. Mol. Pharmacol. 66, 512521.
Chiang, E. F., Pai, C. I., Wyatt, M., Yan, Y. L., Postlethwait, J., and Chung, B. (2001). Two sox9 genes on duplicated zebrafish chromosomes: Expression of similar transcription activators in distinct sites. Dev. Biol. 231, 149163.[CrossRef][Web of Science][Medline]
Czipri, M., Otto, J. M., Cs-Szabo, G., Kamath, R. V., Vermes, C., Firneisz, G., Kolman, K. J., Watanabe, H., Li, Y., Roughley, P. J., et al. (2003). Genetic rescue of chondrodysplasia and the perinatal lethal effect of cartilage link protein deficiency. J. Biol. Chem. 278, 3921439223.
de Crombrugghe, B., Lefebvre, V., Behringer, R. R., Bi, W., Murakami, S., and Huang, W. (2000). Transcriptional mechanisms of chondrocyte differentiation. Matrix Biol. 19, 389394.[CrossRef][Web of Science][Medline]
Fernandez-Salguero, P. M., Hilbert, D. M., Rudikoff, S., Ward, J. M., and Gonzalez, F. J. (1996). Aryl-hydrocarbon receptor-deficient mice are resistant to 2,3,7,8-tetrachlorodibenzo-p-dioxin-induced toxicity. Toxicol. Appl. Pharmacol. 140, 173179.[CrossRef][Web of Science][Medline]
Fletcher, N., Giese, N., Schmidt, C., Stern, N., Lind, P. M., Viluksela, M., Tuomisto, J. T., Tuomisto, J., Nau, H., and Hakansson, H. (2005a). Altered retinoid metabolism in female Long-Evans and Han/Wistar rats following long-term 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-treatment. Toxicol. Sci. 86, 264272.
Fletcher, N., Wahlstrom, D., Lundberg, R., Nilsson, C. B., Nilsson, K. C., Stockling, K., Hellmold, H., and Hakansson, H. (2005b). 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) alters the mRNA expression of critical genes associated with cholesterol metabolism, bile acid biosynthesis, and bile transport in rat liver: a microarray study. Toxicol. Appl. Pharmacol. 207, 124.[Web of Science][Medline]
Frueh, F. W., Hayashibara, K. C., Brown, P. O., and Whitlock, J. P., Jr. (2001). Use of cDNA microarrays to analyze dioxin-induced changes in human liver gene expression. Toxicol. Lett. 122, 189203.[CrossRef][Web of Science][Medline]
Gu, Y. Z., Hogenesch, J. B., and Bradfield, C. A. (2000). The PAS superfamily: Sensors of environmental and developmental signals. Annu. Rev. Pharmacol. Toxicol. 40, 519561.[CrossRef][Web of Science][Medline]
Guo, J., Sartor, M., Karyala, S., Medvedovic, M., Kann, S., Puga, A., Ryan, P., and Tomlinson, C. R. (2004). Expression of genes in the TGF-beta signaling pathway is significantly deregulated in smooth muscle cells from aorta of aryl hydrocarbon receptor knockout mice. Toxicol. Appl. Pharmacol. 194, 7989.[CrossRef][Web of Science][Medline]
Hahn, M. E., Karchner, S. I., Shapiro, M. A., and Perera, S. A. (1997). Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family. Proc. Natl. Acad. Sci. U.S.A. 94, 1374313748.
Handley-Goldstone, H. M., Grow, M. W., and Stegeman, J. J. (2005). Cardiovascular gene expression profiles of dioxin exposure in zebrafish embryos. Toxicol. Sci. 85, 683693.
Hanlon, P. R., Zheng, W., Ko, A. Y., and Jefcoate, C. R. (2005). Identification of novel TCDD-regulated genes by microarray analysis. Toxicol. Appl. Pharmacol. 202, 215228.[CrossRef][Web of Science][Medline]
James, C. G., Appleton, C. T., Ulici, V., Underhill, T. M., and Beier, F. (2005). Microarray analyses of gene expression during chondrocyte differentiation identifies novel regulators of hypertrophy. Mol. Biol. Cell. 16, 53165333.
Kou, I., and Ikegawa, S. (2004). SOX9-dependent and -independent transcriptional regulation of human cartilage link protein. J. Biol. Chem. 279, 5094250948.
Liu, H., McKenna, L. A., and Dean, M. F. (2000). An N-terminal peptide from link protein can stimulate biosynthesis of collagen by human articular cartilage. Arch. Biochem. Biophys. 378, 116122.[CrossRef][Web of Science][Medline]
Lund, L. R., Romer, J., Bugge, T. H., Nielsen, B. S., Frandsen, T. L., Degen, J. L., Stephens, R. W., and Dano, K. (1999). Functional overlap between two classes of matrix-degrading proteases in wound healing. Embo. J. 18, 46454656.[CrossRef][Web of Science][Medline]
Mathew, L. K., Andreasen, E. A., and Tanguay, R. L. (2005). Ahr activation inhibits regenerative growth. Mol. Pharmacol. 69, 257265.[Medline]
Mimura, J., Yamashita, K., Nakamura, K., Morita, M., Takagi, T. N., Nakao, K., Ema, M., Sogawa, K., Yasuda, M., Katsuki, M., et al. (1997). Loss of teratogenic response to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in mice lacking the Ah (dioxin) receptor. Genes Cells 2, 645654.[Abstract]
Mizutani, T., Yoshino, M., Satake, T., Nakagawa, M., Ishimura, R., Tohyama, C., Kokame, K., Kangawa, K., and Miyamoto, K. (2004). Identification of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-inducible and -suppressive genes in the rat placenta: Induction of interferon-regulated genes with possible inhibitory roles for angiogenesis in the placenta. Endocr. J. 51, 569577.[CrossRef][Web of Science][Medline]
Mori-Akiyama, Y., Akiyama, H., Rowitch, D. H., and de Crombrugghe, B. (2003). Sox9 is required for determination of the chondrogenic cell lineage in the cranial neural crest. Proc. Natl. Acad. Sci. U.S.A. 100, 93609365.
Mott, J. D., and Werb, Z. (2004). Regulation of matrix biology by matrix metalloproteinases. Curr. Opin. Cell. Biol. 16, 558564.[CrossRef][Web of Science][Medline]
Partridge, N. C., Fiacco, G. J., Walling, H. W., Barmina, O. Y., Jeffrey, J. J., and Ruh, M. F. (2000). Effects of dioxin and estrogen on collagenase-3 in UMR 10601 osteosarcoma cells. Arch. Biochem. Biophys. 382, 182188.[CrossRef][Web of Science][Medline]
Poss, K. D., Keating, M. T., and Nechiporuk, A. (2003). Tales of regeneration in zebrafish. Dev. Dyn. 226, 202210.[CrossRef][Web of Science][Medline]
Prasch, A. L., Heideman, W., and Peterson, R. E. (2004). ARNT2 is not required for TCDD developmental toxicity in zebrafish. Toxicol. Sci. 82, 250258.
Prasch, A. L., Tanguay, R. L., Mehta, V., Heideman, W., and Peterson, R. E. (2006). Identification of zebrafish ARNT1 homologs: TCDD toxicity in the developing zebrafish requires ARNT1. Mol. Pharmacol. [Epub ahead of print] PMID: 16306231.
Prasch, A. L., Teraoka, H., Carney, S. A., Dong, W., Hiraga, T., Stegeman, J. J., Heideman, W., and Peterson, R. E. (2003). Aryl hydrocarbon receptor 2 mediates 2,3,7,8-tetrachlorodibenzo-p-dioxin developmental toxicity in zebrafish. Toxicol. Sci. 76, 138150.
Schafer, A. J., Foster, J. W., Kwok, C., Weller, P. A., Guioli, S., and Goodfellow, P. N. (1996). Campomelic dysplasia with XY sex reversal: Diverse phenotypes resulting from mutations in a single gene. Ann. N.Y. Acad. Sci. 785, 137149.[Web of Science][Medline]
Sekiya, I., Koopman, P., Tsuji, K., Mertin, S., Harley, V., Yamada, Y., Shinomiya, K., Niguji, A., and Noda, M. (2001). Transcriptional suppression of Sox9 expression in chondrocytes by retinoic acid. J. Cell. Biochem. 81, 7178.
Sekiya, I., Tsuji, K., Koopman, P., Watanabe, H., Yamada, Y., Shinomiya, K., Nifuji, A., and Noda, M. (2000). SOX9 enhances aggrecan gene promoter/enhancer activity and is up-regulated by retinoic acid in a cartilage-derived cell line, TC6. J. Biol. Chem. 275, 1073810744.
Tanguay, R. L., Andreasen, E. A., Walker, M. K., and Peterson, R. E. (2003). Dioxin toxicity and aryl hydrocarbon receptor signaling in fish. In Dioxins and Health (A. Schecter, Ed.), pp. 603628. Plenum Press, New York.
Teraoka, H., Dong, W., Ogawa, S., Tsukiyama, S., Okuhara, Y., Niiyama, M., Ueno, N., Peterson, R. E., and Hiraga, T. (2002). 2,3,7,8-Tetrachlorodibenzo-p-dioxin toxicity in the zebrafish embryo: Altered regional blood flow and impaired lower jaw development. Toxicol. Sci. 65, 192199.
Thackaberry, E. A., Jiang, Z., Johnson, C. D., Ramos, K. S., and Walker, M. K. (2005). Toxicogenomic profile of 2,3,7,8-tetrachlorodibenzo-p-dioxin in the murine fetal heart: Modulation of cell cycle and extracellular matrix genes. Toxicol. Sci. 88, 231241.
Villano, C. M., Murphy, K. A., Akintobi, A., and White, L. A. (2006). 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) induces matrix metalloproteinase (MMP) expression and invasion in A2058 melanoma cells. Toxicol. Appl. Pharmacol. 210, 212224.
Vogel, C. F., Sciullo, E., and Matsumura, F. (2004). Activation of inflammatory mediators and potential role of ah-receptor ligands in foam cell formation. Cardiovasc. Toxicol. 4, 363373.[CrossRef][Medline]
Volz, D. C., Bencic, D. C., Hinton, D. E., Law, J. M., and Kullman, S. W. (2005). 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) induces organ- specific differential gene expression in male Japanese medaka (Oryzias latipes). Toxicol. Sci. 85, 572584.
Yan, Y. L., Willoughby, J., Liu, D., Crump, J. G., Wilson, C., Miller, C. T., Singer, A., Kimmel, C., Westerfield, M., and Postlethwait, J. H. (2005). A pair of Sox: Distinct and overlapping functions of zebrafish sox9 co-orthologs in craniofacial and pectoral fin development. Development 132, 10691083.
Zodrow, J. M., and Tanguay, R. L. (2003). 2,3,7,8-tetrachlorodibenzo-p-dioxin inhibits zebrafish caudal fin regeneration. Toxicol. Sci. 76, 151161.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
L. K. Mathew, S. Sengupta, J. A. Franzosa, J. Perry, J. La Du, E. A. Andreasen, and R. L. Tanguay Comparative Expression Profiling Reveals an Essential Role for Raldh2 in Epimorphic Regeneration J. Biol. Chem., November 27, 2009; 284(48): 33642 - 33653. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Jenny, S. I. Karchner, D. G. Franks, B. R. Woodin, J. J. Stegeman, and M. E. Hahn Distinct Roles of Two Zebrafish AHR Repressors (AHRRa and AHRRb) in Embryonic Development and Regulating the Response to 2,3,7,8-Tetrachlorodibenzo-p-dioxin Toxicol. Sci., August 1, 2009; 110(2): 426 - 441. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Xiong, R. E. Peterson, and W. Heideman Aryl Hydrocarbon Receptor-Mediated Down-Regulation of Sox9b Causes Jaw Malformation in Zebrafish Embryos Mol. Pharmacol., December 1, 2008; 74(6): 1544 - 1553. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. K. Mathew, S. S. Sengupta, J. LaDu, E. A. Andreasen, and R. L. Tanguay Crosstalk between AHR and Wnt signaling through R-Spondin1 impairs tissue regeneration in zebrafish FASEB J, August 1, 2008; 22(8): 3087 - 3096. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. K. Mathew, S. Sengupta, A. Kawakami, E. A. Andreasen, C. V. Lohr, C. A. Loynes, S. A. Renshaw, R. T. Peterson, and R. L. Tanguay Unraveling Tissue Regeneration Pathways Using Chemical Genetics J. Biol. Chem., November 30, 2007; 282(48): 35202 - 35210. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Chang, Y. Fan, S. Karyala, S. Schwemberger, C. R. Tomlinson, M. A. Sartor, and A. Puga Ligand-Independent Regulation of Transforming Growth Factor {beta}1 Expression and Cell Cycle Progression by the Aryl Hydrocarbon Receptor Mol. Cell. Biol., September 1, 2007; 27(17): 6127 - 6139. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. A. Andreasen, L. K. Mathew, C. V. Lohr, R. Hasson, and R. L. Tanguay Aryl Hydrocarbon Receptor Activation Impairs Extracellular Matrix Remodeling during Zebra Fish fin Regeneration Toxicol. Sci., January 1, 2007; 95(1): 215 - 226. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||










