|
|
|
|||
| Home Help Feedback Subscriptions Archive Search Table of Contents | ||||
First published online September 19, 2008
Journal of Experimental Biology 211, 3077-3084 (2008)
Published by The Company of Biologists 2008
doi: 10.1242/jeb.019950
Effects of hypothermia on gene expression in zebrafish gills: upregulation in differentiation and function of ionocytes as compensatory responses
1 Institute of Cellular and Organismic Biology, Academia Sinica, Taipei, 115
Taiwan
2 Institute of Fishery Science, National Taiwan University, Taipei, 10617
Taiwan
3 Department of Bioscience Technology, Chung Yuan Christian University, Chung
Li, 32023 Taiwan
* Author for correspondence (e-mail: pphwang{at}gate.sinica.edu.tw)
Accepted 24 July 2008
| Summary |
|---|
|
|
|---|
Key words: acclimation, cold, differentiation, ionocyte, transcriptome, transporter
| INTRODUCTION |
|---|
|
|
|---|
In contrast to endothermic mammals, ectothermic vertebrates such as fish
can survive in a wide range of thermal environments, and thus fish cells may
utilize different strategies to cope with thermal fluctuations. Many
physiological responses, including changes in lipid composition
(Hazel, 1979
), increases in
pump activity and specific Na+/K+-ATPase activity
(Schwarzbaum et al., 1992a
;
Schwarzbaum et al., 1992b
) and
oxygen consumption (Raynard and Cossins,
1991
), have been extensively reported when fish are exposed to
cold stress, however there still remain several interesting issues for further
studies.
Ambient temperatures have a direct and critical impact on the physiological
functioning of fish; hence, fish are excellent experimental organisms in which
to investigate responses to environmental stresses
(Cossins and Crawford, 2005
).
To cope with fluctuating temperatures, fish can increase the quantity of
enzymes required for different physiological tasks, or can express protein
isoforms that are more suitable for changing thermal conditions
(Cossins et al., 2002
;
Driedzic and Ewart, 2004
;
Johnston et al., 1990
;
Watabe, 2002
). Cold stress
triggers a complex program of gene expression and biochemical responses in
different tissues (Gracey et al.,
2004
; Ju et al.,
2002
; Malek et al.,
2004
; Vornanen et al.,
2005
). The gills are the first target organ for environmental
stress because they are directly exposed to the ambient environment. Fish
gills have multiple functions including gas exchange, ion regulation, nitrogen
balance and acid–base adjustment. Gill ionocytes are the major response
site for osmo- and pH regulation (Evans et
al., 2005
), and the cell cycle of gill cells is about 4–5
days (Tsai and Hwang, 1998
),
indicative of a higher cell turnover rate of tissues with complex and critical
physiological activities. Moreover, gills showed great morphological
plasticity when acclimated to different ambient salinity, ionic composition,
oxygen and temperature conditions (Chang
and Hwang, 2004
; Chang et al.,
2001
; Sollid and Nilsson,
2006
). Therefore, gills may provide a suitable target for
physiological responses to environmental changes.
Zebrafish (Danio rerio), a tropical teleost species, can survive
in low-temperature (18°C) environments for 1 year and show normal motor
activity compared to control fish at 28°C
(Malek et al., 2004
). Tang and
colleagues also reported that zebrafish can survive in water less than
12°C for at least 48 h (Tang et al.,
1999
). In contrast to other species, zebrafish are small, and the
genetic resources of this model species are well known and abundant
(Malek et al., 2004
).
Zebrafish are therefore suitable for use as a model animal to study the
mechanisms of acclimation to low temperatures in fish. The purpose of this
study was to use zebrafish as a model to investigate the effects of low
temperature on the functions of gill cells. We used microarray technology to
compare the gill transcriptome between acute and chronic low-temperature
treatments. A group of upregulated genes related to ion balance and
acid–base regulation mechanisms in branchial ionocytes were identified.
We demonstrated that elevation of Ca2+ influx was accompanied by an
increase of trpv6 (epithelial Ca2+ channel) expression.
Moreover, we found cell proliferation and apoptosis in zebrafish gill cells
after cold stress. These data provide functional genomic and physiological
evidences for cold acclimation in fish gills.
| MATERIALS AND METHODS |
|---|
|
|
|---|
A total of 192 surviving fish were sacrificed for the subsequent microarray and quantitative reverse-transcription polymerase chain reaction (qRT–PCR) analysis. In order to obtain sufficient quantity of RNA, the whole gills from both sides of six individuals (three male and three female fish) were pooled as a sample. For microarray experiments, 18 fish were incubated in three different experimental tanks (i.e. N=3) at 12°C for 1 day and another 18 fish were incubated for 30 d. Control groups with 18 fish for 1 d and another 18 fish for 30 d were both incubated in 28°C. For qRT–PCR, 30 fish were incubated in 5 different experimental tanks (i.e. N=5) at 12°C for 1 day and another 30 were incubated for 30 days. Another 60 fish for the controls were treated as described above except the acclimation temperature.
For whole-mount in situ hybridization, immunohistochemistry and TUNEL assay, four fish (N=4) were used for each test group. In western blotting experiment, 48 fish, 24 for cold treatment and 24 for control (gills from six fish were pooled as a sample, i.e. N=4), were sacrificed for protein isolation. Fish were anesthetized with 100–200 mg l–1 of buffered MS-222 (3-aminobenzoic acid ethyl ester; Sigma, St Louis, MO, USA) before sampling following, guidelines of the Academia Sinica Institutional Animal Care and Utilization Committee (approval no. RFiZOOHP2007086).
RNA extraction
Dissected and pooled gill tissues were homogenized in 5 ml Trizol reagent
(Invitrogen, Carlsbad, CA, USA). After chloroform extraction, RNA
precipitation and ethanol washing, the RNA samples were purified and treated
with DNase1 to remove the genomic DNA by using RNeasy Mini Kit (Qiagen,
Huntsville, AL, USA). The quantity and quality of total RNA were assessed by
spectrophotometry and agarose gel electrophoresis, respectively.
Microarray hybridization and data analysis
The commercial zebrafish 14K oligonucleotide set (MWG Biotech AG,
Ebersbach, Germany) were obtained and were printed on an UltraGAPS-coated
slide (Corning, New York, NY, USA) with the use of the OmniGrid 100
microarrayer (Genomic Solutions, Ann Arbor, MI, USA) according to the
manufacturer's instructions. The 14 067 oligonucleotides represent 9666 genes
(7009 singlet genes and 2657 redundant genes), and the redundancy of this chip
is 31%. The detailed description of the oligonucleotide information can be
obtained on the Ocimun Biosolution website
(http://www.ocimumbio.com/web/default.asp).
cDNA probes were synthesized by reverse transcription of 20 µg total RNA using a SuperScript indirect cDNA labeling system (Invitrogen) and were labeled with Cy5 (cold treatment groups) and Cy3 (control groups) (Amersham Bioscience, Buckinghamshire, UK), respectively. The zebrafish 14K OciChip array chip was pretreated with 1% bovine serum albumin (BSA; fraction V), 4x SSC buffer and 1% sodium dodecylsulfate (SDS) at 42°C for 45 min, and then hybridized overnight in a cocktail containing 5x Denhardt's solution, 6x SSC, 0.5% SDS, 50% formamide, 50 mmol l–1 sodium phosphate and 2 µgl–1 yeast tRNA. Slides were washed with 2x SSC and 0.1% SDS (5 min), 1x SSC and 0.1% SDS (5 min), 0.5x SSC (5 min), and twice with 0.1x SSC (2 min each). Scanning was performed with a Genepix scanner (Molecular Devices, Sunnyvale, CA, USA). The acquired images were analyzed using Genepix and Genespring software (Aglient Technologies, Foster City, CA, USA). The measurements of spots were filtered by flags, and the Lowess normalization was performed after subtraction of the median background. Each experiment contained three biological replicates (including one dye swap) with different samples. In total six chips (three chips for 1 day and three chips for 30 days) were used for microarray hybridization experiments, and two (one chip for 1 day and one chip for 30 days) of the six chips were used for dye swap. Thus, 12 biological replicates were used for hybridization including dye swaps. The differentially expressed genes were selected from those with at least two of three significant signals (ratio >2 or <0.5), and then the Significant Analysis of Microarray method was used to determine statistical significances. Data were submitted to NCBI Gene Expression Omnibus (series accession no. GSE7853).
Validation of differentially expressed genes by quantitative reverse-transcription polymerase chain reaction (qRT–PCR)
We used cold-inducible RNA binding protein (cirbp,
BC057481) and high-mobility group box 1 (hmgb1, NM_199555)
to test the general responses to cold. To confirm the differentially expressed
genes collected from the microarray analysis, we selected nine transcripts
including Ca2+ transporting, cardiac muscle, fast twitch 1
(atp2a1, NM_001007029), Ca2+ ATPase, cardiac muscle,
fast twitch 1 like (atp2a1l, NM_001077533), tubulin alpha 8
like 2 (tuba8l2, NM_200691), GTP binding protein 4
(gtpbp4, NM_199851), Na+/K+-ATPase, alpha
1a.4 polypeptide (atp1a1a.4, NM_131689), glycogen synthase
kinase binding protein (gbp, NM_131442), V-ATPase subunit
A (atp6v1a, NM_201135), keratin 18 (krt18,
NM_178437) and annexin A2a (anxa2a, NM_181761) to validate
their relative expression levels by quantitative RT–PCR (N=5).
We also measured the mRNA expressional levels of ionoctye-related genes, such
as forkhead box I3a (foxi3a, NM_198917), N-myc
downstream-regulated gene 1 (ndrg1, NM_213348), carbonic
anhydrase II (ca2, NM_199215), delta C (dlc,
NM_130944), Na+/K+-ATPase, beta 1b polypeptide
(atp1b1b, NM_131671), Na+/K+-ATPase, alpha
1a.2 subunit (atp1a1a.2, NM_131687),
Na+/Cl– cotransporter (slc12a3,
EF591989), carbonic anhydrase 15a isoform (ca15a EF591981),
glial cells missing homolog 2 (gcm2, NM_001005603), and
epithelial Ca2+ channel (trpv6, NM_001001849), to
systematically analyze the expressional profiles of ionoregulatory genes. As
an internal control, primers for β-actin (NM_131031) were designed and
amplified in parallel with the genes of interest. Quantitative
reverse-transcription PCR (qRT–PCR) was carried out using a SYBR Green
dye (Qiagen, Hilden, Germany)-based assay with an ABI Prism 7000 Sequence
Detection System (Perkin-Elmer, Applied Biosystems, Wellesley, MA, USA)
according to the manufacturer's instructions. Primer targeting was designed
using Primer Express 2.0 software (Applied Biosystems). The primer sequences
are listed in supplementary material Table S1.
RNA probe synthesis
Two zebrafish nucleotide fragments from, trpv6 and ca2,
were obtained by PCR and inserted into the pGEM-T Easy vector (Promega,
Madison, WI, USA). Purified plasmids were then linearized by restriction
enzyme digestion, and in vitro transcription was carried out with T7
or SP6 RNA polymerase (Roche, Penzberg, Germany) in the presence of
digoxigenin (Dig)-UTP. Dig-labeled RNA probes were examined with RNA gels and
a dot-blot assay to confirm the quality and concentration. For the dot-blot
assay, the synthesized probes and standard RNA probes were spotted onto
nitrocellulose membranes according to the manufacturer's instructions (Dig RNA
labeling kit; Roche Diagnostics, Mannheim, Germany). After cross-linking and
blocking, the membranes were incubated with an alkaline phosphatase-conjugated
anti-dig antibody and stained with nitro blue tetrazolium (NBT; Roche) and
5-bromo-4-chloro-3-indolyl phosphate (BCIP; Roche).
Whole-mount in situ hybridization
Zebrafish gills were fixed with 4% paraformaldehyde overnight at 4°C,
and then washed several times with phosphate-buffered saline (PBS). Fixed
samples were rinsed with PBST (PBS with 0.2% Tween 20, 1.4 mmol
l–1 NaCl, 0.2 mmol l–1 KCl, 0.1 mmol
l–1 Na2HPO4, and 0.002 mmol
l–1 KH2PO4; pH 7.4). After a brief
washing with PBST, gill filaments were incubated with hybridization buffer
(HyB) containing 60% formamide, 5x SSC, and 0.1% Tween 20 for 5 min at
65°C. Prehybridization was performed in HyB+ (the hybridization
buffer supplemented with 500 µgml–1 yeast tRNA and 50
µgml–1 heparin) for 2 h at 65°C. After
prehybridization, samples were hybridized in 100 ng of the RNA probe in 200
µl of HyB+ at 65°C overnight. Gills were then washed at
65°C for 10 min in 75% HyB and 25% 2x SSC, for 10 min in 50% HyB and
50% 2x SSC, for 10 min in 25% HyB and 75% 2x SSC, for 10 min in
2x SSC, and twice for 30 min each in 0.2x SSC at 70°C. Further
washes were performed at room temperature for 5 min in 75% 0.2x SSC and
25% PBST, for 5 min in 50% 0.2x SSC and 50% PBST, for 5 min in 25%
0.2x SSC and 75% PBST, and for 5 min in PBST. After serial washings,
gill filaments were incubated in blocking solution containing 5% sheep serum
and 2 mg ml–1 BSA in PBST for 2 h and then incubated in the
1:10 000-diluted alkaline phosphatase-conjugated anti-dig antibody for another
16 h at 4°C. After the reaction, samples were washed with PBST plus
blocking reagent and then stained with NBT and BCIP.
Measurement of whole body Ca2+ influx
Whole body Ca2+ influx was measured following the method of
Chang and Hwang (Chang and Hwang,
2004
) with some modifications. Control and cold-acclimated fish
were transferred to 20 ml [45Ca2+]-containing medium for
5.5 h incubation. Following the incubation, 200 µl water medium were
sampled with addition of 2 ml counting solution (Ultima Gold, Packard,
Waltham, MA, USA), and then the radioactivities of the solutions were
determined with a LS6500 beta counter (Beckman, CA, USA). The
45Ca2+ influx was calculated using the following
formula:
![]() |
Immunohistochemistry
Zebrafish gills were fixed in 4% paraformaldehyde for 12 h at 4°C.
After being washed in PBS, fixed gills were treated with 100% methanol for 10
min at –20°C and subsequently subjected to blocking with 3% BSA at
room temperature for 60 min. Gill filaments were then incubated with a
polyclonal antibody against a short amino acid sequence containing
phosphorylated Ser 10 of histone H3 of human origin (phosphohistone H3; Santa
Cruz Biotechnology, Santa Cruz, CA, USA) diluted to 1:200, at 4°C for 16
h. Samples were washed twice in PBS for 10 min each and then incubated with
1:200 PBS-diluted Rhodamine-conjugated goat anti-rabbit IgG (Jackson
ImmunoResearch Laboratories, West Grove, PA, USA) for 1 h at room
temperature.
Western blotting
The gills were homogenized in homogenization solution (100 mmol
l–1 imidazole, 5 mmol l–1 EDTA, 200 mmol
l–1 sucrose, 0.1% sodium deoxycholate; pH 7.6) and subjected
to polyacrylamide gel electrophoresis (PAGE) in 8x10 cm sodium
dodecylsulfate (SDS)-polyacrylamide (8%) gels at 100 V for 2 h. Protein was
loaded at 30 µg per well. Separated proteins were transferred onto
polyvinylidene difluoride membranes (Millipore, Billerica, MA, USA) with a
transfer electrophoresis unit (SE 600, Hoffer) at 100 V for 2 h. After
blocking for 1.5 h in 5% nonfat dried milk, the blots were incubated with
phosphohistone H3 antibodies (overnight, diluted 1:1000). After incubation,
the membranes were washed in PBST and reacted for 90 min with an
alkaline-phosphatase-conjugated secondary antibody (Jackson ImmunoResearch
Laboratories; diluted 1:3000). Blots were extensively washed with PBST and
then developed with 0.015% NBT and 0.007% bromochloroindolyl phosphate in a
reaction buffer containing 100 mmol l–1 Tris, 100 mmol
l–1 NaCl and 5 mmol l–1 MgCl2 (pH
9.5).
TUNEL assay
The in situ cell death detection kit (Roche Diagnostics,
Indianapolis, IN, USA) was used for TUNEL detection of DNA fragmentation.
Zebrafish gills were fixed in 4% paraformaldehyde and stored in methanol at
–20°C. Zebrafish gills were incubated with 3%
H2O2 to block endogenous peroxidases, washed, and
treated with 10 g ml–1 proteinase K at, 37°C for 15 min.
Then, gill samples were incubated with fluorescein-conjugated nucleotides and
terminal deoxynucleotidyl transferase at 37°C for 1 h. After application
of the anti-fluorescein antibody conjugated with peroxidase at 37°C for 30
min, the 3,3'-diaminobenzidine (DAB) staining method was performed at
25°C for 2–5 min.
Equipment (image acquisition)
Bright-field and fluorescence images were acquired with a Zeiss Axioplan 2
Imaging MOT microscope (Carl Zeiss, Göttingen, Germany). The
z-axis series images were captured with a Zeiss Axiocam HRm cool CCD
(Carl Zeiss) with AxioVision LE Rel 4.3 software (Carl Zeiss) and then merged
by Helicon Focus 4.01.1 (Helicon Soft Ltd, Kharkov, Ukraine).
Statistical analysis
Significant analysis of microarray (SAM) method was used for microarray
analysis, and one-way ANOVA was used for the other analysis. For qRT-PCR
analysis, controls did not show significant differences between 1 day and 30
days, therefore only the 1 day data was used for the subsequent
comparisons.
| RESULTS |
|---|
|
|
|---|
|
|
There were more cold-suppressed genes than cold-induced genes. The largest GO category of both cold-suppressed and cold-induced transcripts was cellular metabolism; however, the components of cellular metabolism differed between the suppressed and induced groups.
Effects of low temperature treatment on expression of ionocyte-related genes and whole body Ca2+ influx
According to the data set of the cold-induced transcripts, atp6v1a
is involved in branchial ion transport and acid–base regulation
mechanisms in gill ionocytes. Because cold has been reported to affect
transport mechanisms and disrupt ion and acid–base homeostasis
(Hochachka, 1986
;
Hochachka, 1988
), it was
hypothesized that zebrafish may globally activate ion and acid–base
balance mechanisms to compensate for the passive lose of ions and the
imbalance in acid–base regulation. To test this hypothesis, the mRNA
expression levels of ionocyte-related genes including ca2, atp6v1a,
atp1b1b, atp1a1a.2, slc12a3, ca15a and trpv6 were examined. As
expected, all the ionocyte marker genes were significantly upregulated during
cold acclimation (summarized in Table
1). A further question asked was whether cold stress affects the
spatial distribution of ionocyte marker gene expression. In whole-mount in
situ hybridization experiments, expression of trpv6 and
ca2 mRNAs (markers of gill ionocytes) were evidently increased, and
notably, trpv6- and ca2-expressing cells were found in the
lamellae in addition to the gill filaments
(Fig. 3), whereas in the
control group, these cells only appeared in gill filaments. In further
experiments, the function of one of these genes, trpv6, was analyzed.
Whole body Ca2+ influx in zebrafish was retarded about 60% by
12°C treatment for 9 h compared with the control group at 28°C, but
was able to recover to the normal level after subsequent acclimation to
12°C for 30 days (Fig.
4).
|
|
|
The differentiation of gill ionocytes was also investigated to determine
whether they are affected during acclimation to low temperature. As shown in
Table 1, expressions of
dlc, ndrg1, gcm2 and foxi3a mRNAs, which may be involved in
ionocyte differentiation (Hsiao et al.,
2007
; Jones et al.,
1995
; Yu et al.,
1999
), were all upregulated after acclimation to low temperature
(Table 1).
Effects of temperature reduction on cell proliferation and apoptosis in zebrafish gills
For evaluation of cell proliferation rate, an M-phase cell cycle marker,
phosphohistone H3, was used for immunohistochemistry and western blot in
gills. Cell proliferation rate appeared to slow down after low temperature
treatment for 30 days (Fig.
5A). Western blot result indicated that the intensity of
phosphohistone H3 immuno-reacted band in low temperature-treated gills was
much weaker than that of control gills, supporting a lower proliferation rate
in immunohistochemistry data (Fig.
5B). We also compared TUNEL assay in control and low temperature
treated gills. The apoptotic cells in cold treatment group were much less than
in control group (Fig. 5C).
|
| DISCUSSION |
|---|
|
|
|---|
Upregulated genes in cold-acclimated gills
The largest GO category of cold-induced transcripts was cellular
metabolism, which could be further subdivided into several subgroups including
glycolysis, regulation of transcription, nucleosome assembly, proteolysis,
protein polymerization, chromatin assembly or disassembly, protein amino acid
dephosphorylation, and protein polymerization (supplementary material Table
S2). The number of upregulated genes in zebrafish acclimated to cold for 30
days were more and the functions of these genes were diverse compared to fish
acclimated for only 1 day. Most previous studies examined the time course
changes of transcriptome in fish after thermal or hypoxic stress, but
overlooked the time-dependent modulation of physiological responses
(Gracey et al., 2004
;
Gracey et al., 2001
;
Ton et al., 2003
). In annual
killifish (Austrofundulus limnaeus) liver, small and large molecular
chaperons were differentially stimulated in response to seasonal temperature
changes and daily fluctuating temperatures, respectively
(Podrabsky and Somero, 2004
).
These studies proposed the gradual enhancement and achievement of compensatory
mechanisms of gill functions following cold acclimation.
In general, the genes induced by cold in zebrafish gills are involved in
various functions including cellular metabolism, ion transport, carbohydrate
metabolism, antigen processing and presentation, immune response and signal
transduction. Similar gene expression profiles were also reported in the heart
of rainbow trout (Vornanen et al.,
2005
) and various organs of carp
(Gracey et al., 2004
) after
cold acclimation. In rainbow trout heart, genes related to protein synthesis
and intermediary metabolism were most strongly upregulated by cold, whereas
the transcription regulation-related genes were upregulated in carp gill,
kidney, brain, heart, muscle, liver and intestine. These studies suggest that
some cold responses occur in a variety of organs and are conserved among
organisms.
Downregulated genes in cold-acclimated gills
In zebrafish gills there were more genes that were suppressed under cold
stress than were upregulated (supplementary material Table S2). Based on the
GO analysis, subgroups of cold-suppressed transcripts were those involved in
DNA repair, protein folding, protein amino acid phosphorylation and response
to oxidative stress. The downregulated genes differed significantly between
the 1 day and 30 days acclimated groups. Cellular metabolism was still the
largest GO category in the downregulated gene lists of both groups. In the 1
day group, the six genes (ercc3, polb, cct3, mespb, LOC571699 and
zgc:110755) in the subgroup of DNA-dependent metabolism were the most abundant
in the cellular metabolism category. Calcium ion transport-related genes
(guca1b, actn4, calb2l and tnnc1) were also downregulated
(supplementary material Table S2).
In the 30 day group, the most intriguing finding was that six intracellular
transport genes (ndel1b, zbtb16, sec23b, ucp2, slc25a12 and
zgc:110821) were downregulated. These transporters are responsible for
trafficking of ions and proteins and important for cellular physiology. The
expression patterns of these genes was contrary to those of several cell
membrane bound transporters, which were induced by cold (discussed below).
This suggests that zebrafish activated some gill-specific responses during
cold acclimation. In addition, the immune response gene, mhclze,
which belongs to the major histocompatibility complex (MHC) family were
downregulated. MHC proteins play important roles in immune responses to
bacterial and fungal pathogens (Ojcius et
al., 1994
), and consequently, downregulation of these genes at low
temperature may increase the chance of infection by low-temperature-related
pathogens. In rainbow trout, the mRNA and protein levels of MHC II alpha and
beta were downregulated after 2°C of cold treatment
(Nath et al., 2006
). However,
rainbow trout and Atlantic salmon expressed high levels of beta-2-microgobulin
in a 2°C environment (Kales et al.,
2006
), allowing them to maintain their viral recognition machinery
at low temperatures.
In cold-exposed carp and trout only a few genes were suppressed
(Gracey et al., 2004
;
Vornanen et al., 2005
) in
contrast to zebrafish in which there were over 100 downregulated genes. The
most reasonable explanation for these differences may be to the hereditary
histories of these species. Both carp and rainbow trout are temperate species,
which have evolved sufficient mechanisms for adaptation to habitats with
temperatures near 0°C (Jain and
Farrell, 2003
; Sollid and
Nilsson, 2006
). Zebrafish, however, are considered to be a
tropical species, for which a low-temperature environment is much more
challenging. Taking all of these into account, the changes in gene expression
profiles during acclimation to low temperatures appear to reflect differential
evolutionary and environmental adaptations among species.
Gill-specific responses during acclimation to low temperature
With the aid of functional genomics, Gracey and colleagues
(Gracey et al., 2004
) extended
our knowledge and understanding of how cold exposure elicits different
responses in various tissues of fish, and other work has supported their
findings (Ju et al., 2002
;
Malek et al., 2004
;
Tang et al., 1999
;
Vornanen et al., 2005
).
However, there was very little integration of the transcriptome analysis into
physiological functions. The present study further extended the transcriptome
data to the functional analysis of fish gills, which are a multi-function
organ and the main extra-renal site responsible for ion balance and
acid–base regulation in fish (Evans
et al., 2005
). Gene expression profiling in gills indicated that
ion balance and acid–base regulation, the principal biological functions
of gill, were affected by low temperature. A group of ion and acid–base
balance-related genes, including transporters (trpv6, atp1b1b, atp1a1a.2,
atp6v1a and slc12a3) and cytosolic enzymes (ca15a and
ca2) were significantly upregulated by cold. Epithelial
Ca2+ channel (ECaC; trpv6) is the key ion channel of
Ca2+ absorption in zebrafish and rainbow trout
(Pan et al., 2005
;
Shahsavarani et al., 2006
;
Shahsavarani and Perry, 2006
),
and V-type H-ATPase (atp6v1a) plays a major role in acid secretion in
zebrafish embryos (Horng et al.,
2007
; Lin et al.,
2006
). Carbonic anhydrase (ca2 and ca15a) also
participates in acid–base balance
(Claiborne et al., 2002
;
Georgalis et al., 2006
;
Hwang and Lee, 2007
;
Lin et al., 2008
). NCC
(slc12a3) is responsible for chloride uptake in zebrafish
(Hwang and Lee, 2007
).
Na+/K+-ATPase (atp1a1a.2 and atp1b1b)
is a major driving force for other transporters in mitochondrial-rich cells
(Chang and Hwang, 2004
;
Evans et al., 2005
;
Lee et al., 1998
).
Furthermore, the Ca2+ influx, as a function of trpv6, had
recovered after 30 days cold acclimation coincidently with the elevation of
trpv6 mRNA expression, providing the evidence to integrate the
molecular findings of this study to physiological function. Taking all these
findings into account, it is apparent that ion transporters in gill ionocytes
were activated to increase ion uptake and to adjust the pH of body fluids in
zebrafish, and these may probably recover the ion and/or acid–base
balance that has been disturbed by an acute cold stress.
We examined the expression of foxi3a, dlc, gcm2 and ndrg1
mRNAs to investigate the effects of cold on ionocytes differentiation-related
genes. NDRG1 has been reported to be modulated under various conditions such
as cellular differentiation (Yu et al.,
1999
), cell cycle arrest
(Piquemal et al., 1999
) and
hypoxia (Chen et al., 2006
).
Glial cells missing (gcm), primarily identified from
Drosophila, is a genetic switch, controlling glial versus
neuronal fate (Jones et al.,
1995
). In both mouse and chicken, gcm2 is expressed in
the pharyngeal pouches and the forming parathyroid gland, whereas in teleosts,
gcm2 is expressed within the pharyngeal pouches and directs
development of the internal gill bud
(Hogan et al., 2004
;
Okabe and Graham, 2004
).
Recent studies demonstrated that a member of the forkhead transcription factor
family, foxi3a, functions as a master regulator and dlc is a
lateral inhibitor for ionocyte differentiation in zebrafish embryos
(Hsiao et al., 2007
). All
these cell differentiation-related genes, foxi3a, dlc, gcm2 and
ndrg1, were evidently stimulated after acclimation to cold. Cold
stimulates the differentiation of ionocytes in zebrafish gills, supporting our
hypothesis described above.
Increases in ionocyte numbers and expression of functional genes were observed in zebrafish gills after cold acclimation. Experiments on phosphohistone H3 and TUNEL assays indicated that both cell proliferation and apoptosis in gills were decreased after cold acclimation. These imply that the increased ionocytes may mainly originate from acceleration of the terminally differentiating pre-existing un-differentiated or immature cells located in the gill rather than from cell proliferation.
In cultured mammalian cells, cold shock retarded cell proliferation that
leads to apoptosis (Al-Fageeh et al.,
2006
; Rieder and Cole,
2002
; Sonna et al.,
2002
). By contrast, in the intestines of hibernating mammals, the
rates of both cell proliferation and apoptosis are suppressed
(Fleck and Carey, 2005
).
Similar to hibernating mammals, the overall cell proliferation and apoptosis
rates in zebrafish gills declined after cold treatment. The collective
suppression of cell proliferation (arrest of the cell cycle) and the delay of
apoptosis may result in extension of the cell lifespan. In cold environments,
zebrafish need to survive and maintain their normal physiological performance
for growth and reproduction. Zebrafish gill ionocytes displayed extended
lifespan (due to delayed apoptosis) and sustained cell functions (due to
stimulation of pre-existing undifferentiated cells into ionocytes) after cold
acclimation. These findings provide new insights into the cellular
physiological mechanisms of survival and growth of ectothermic vertebrates in
low-temperature environments.
In summary, expression of genes related to ion and acid–base regulation in the gill were stimulated during cold acclimation, suggesting an essential compensatory action for the cold-induced ion imbalance to stabilize normal physiological processes in zebrafish.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
Al-Fageeh, M. B., Marchant, R. J., Carden, M. J. and Smales, C. M. (2006). The cold-shock response in cultured mammalian cells: harnessing the response for the improvement of recombinant protein production. Biotechnol. Bioeng. 93,829 -835.[CrossRef][Medline]
Bianchi, M. E. and Manfredi, A. (2004). Chromatin and cell death. Biochim. Biophys. Acta 1677,181 -186.[Medline]
Buffenstein, R. (2005). The naked mole-rat: a new long-living model for human aging research. J. Gerontol. 60,1369 -1377.
Chang, I. C. and Hwang, P. P. (2004). Cl- uptake mechanism in freshwater-adapted tilapia (Oreochromis mossambicus). Physiol. Biochem. Zool. 77,406 -414.[CrossRef][Medline]
Chang, I. C., Lee, T. H., Yang, C. H., Wei, Y. Y., Chou, F. I. and Hwang, P. P. (2001). Morphology and function of gill mitochondria-rich cells in fish acclimated to different environments. Physiol. Biochem. Zool. 74,111 -119.[CrossRef][Medline]
Chen, B., Nelson, D. M. and Sadovsky, Y.
(2006). N-myc down-regulated gene 1 modulates the response of
term human trophoblasts to hypoxic injury. J. Biol.
Chem. 281,2764
-2772.
Claiborne, J. B., Edwards, S. L. and Morrison-Shetlar, A. I. (2002). Acid-base regulation in fishes: cellular and molecular mechanisms. J. Exp. Zool. 293,302 -319.[CrossRef][Medline]
Cossins, A. R. and Crawford, D. L. (2005). Fish as models for environmental genomics. Nat. Rev. Genet. 6, 324-333.[CrossRef][Medline]
Cossins, A. R., Murray, P. A., Gracey, A. Y., Logue, J., Polley, S., Caddick, M., Brooks, S., Postle, T. and Maclean, N. (2002). The role of desaturases in cold-induced lipid restructuring. Biochem. Soc. Trans. 30,1082 -1086.[CrossRef][Medline]
Delrow, J. J., Heath, P. J., Fujimoto, B. S. and Schurr, J. M. (1998). Effect of temperature on DNA secondary structure in the absence and presence of 0.5 M tetramethylammonium chloride. Biopolymers 45,503 -515.[CrossRef][Medline]
Driedzic, W. R. and Ewart, K. V. (2004). Control of glycerol production by rainbow smelt (Osmerus mordax) to provide freeze resistance and allow foraging at low winter temperatures. Comp. Biochem. Physiol. 139,347 -357.[CrossRef][Medline]
Evans, D. H., Piermarini, P. M. and Choe, K. P.
(2005). The multifunctional fish gill: dominant site of gas
exchange, osmoregulation, acid-base regulation, and excretion of nitrogenous
waste. Physiol. Rev. 85,97
-177.
Fleck, C. C. and Carey, H. V. (2005). Modulation of apoptotic pathways in intestinal mucosa during hibernation. Am. J. Physiol. 289,R586 -R595.
Fujita, J. (1999). Cold shock response in mammalian cells. J. Mol. Microbiol. Biotechnol. 1, 243-255.[Medline]
Georgalis, T., Perry, S. F. and Gilmour, K. M.
(2006). The role of branchial carbonic anhydrase in acid-base
regulation in rainbow trout (Oncorhynchus mykiss). J. Exp.
Biol. 209,518
-530.
Gracey, A. Y., Troll, J. V. and Somero, G. N.
(2001). Hypoxia-induced gene expression profiling in the euryoxic
fish Gillichthys mirabilis. Proc. Natl. Acad. Sci. USA
98,1993
-1998.
Gracey, A. Y., Fraser, E. J., Li, W., Fang, Y., Taylor, R. R.,
Rogers, J., Brass, A. and Cossins, A. R. (2004). Coping with
cold: An integrative, multitissue analysis of the transcriptome of a
poikilothermic vertebrate. Proc. Natl. Acad. Sci. USA
101,16970
-16975.
Guderley, H. (1990). Functional significance of metabolic responses to thermal acclimation in fish muscle. Am. J. Physiol. 259,R245 -R252.[Medline]
Hazel, J. R. (1979). Influence of thermal acclimation on membrane lipid composition of rainbow trout liver. Am. J. Physiol. 236,R91 -R101.[Medline]
Hochachka, P. W. (1986). Defense strategies
against hypoxia and hypothermia. Science
231,234
-241.
Hochachka, P. W. (1988). Metabolic-coupled, channel-coupled, and pump-coupled functions – constraints and compromises of coadaptation. Can. J. Zool. 66,1015 -1027.[CrossRef]
Hogan, B. M., Hunter, M. P., Oates, A. C., Crowhurst, M. O., Hall, N. E., Heath, J. K., Prince, V. E. and Lieschke, G. J. (2004). Zebrafish gcm2 is required for gill filament budding from pharyngeal ectoderm. Dev. Biol. 276,508 -522.[CrossRef][Medline]
Horng, J. L., Lin, L. Y., Huang, C. J., Katoh, F., Kaneko, T. and Hwang, P. P. (2007). Knockdown of V-ATPase subunit A (atp6v1a) impairs acid secretion and ion balance in zebrafish (Danio rerio). Am. J. Physiol. 292,R2068 -R2076.
Hsiao, C. D., You, M. S., Guh, Y. J., Ma, M., Jiang, Y. J. and Hwang, P. P. (2007). A positive regulatory loop between foxi3a and foxi3b is essential for specification and differentiation of zebrafish epidermal ionocytes. PLoS ONE 2, e302.[CrossRef]
Hsieh, S. L. and Kuo, C. M. (2005). Stearoyl-CoA desaturase expression and fatty acid composition in milkfish (Chanos chanos) and grass carp (Ctenopharyngodon idella) during cold acclimation. Comp. Biochem. Physiol. 141,95 -101.[CrossRef][Medline]
Hwang, P. P. and Lee, T. H. (2007). New insights into fish ion regulation and mitochondrion-rich cells. Comp. Biochem. Physiol. 148,479 -497.
Jain, K. E. and Farrell, A. P. (2003).
Influence of seasonal temperature on the repeat swimming performance of
rainbow trout Oncorhynchus mykiss. J. Exp. Biol.
206,3569
-3579.
Johnston, I. A., Fleming, J. D. and Crockford, T. (1990). Thermal acclimation and muscle contractile properties in cyprinid fish. Am. J. Physiol. 259,R231 -R236.[Medline]
Jones, B. W., Fetter, R. D., Tear, G. and Goodman, C. S. (1995). glial cells missing: a genetic switch that controls glial versus neuronal fate. Cell 82,1013 -1023.[CrossRef][Medline]
Ju, Z., Dunham, R. A. and Liu, Z. (2002). Differential gene expression in the brain of channel catfish (Ictalurus punctatus) in response to cold acclimation. Mol. Genet. Genomics 268,87 -95.[CrossRef][Medline]
Kales, S., Parks-Dely, J., Schulte, P. and Dixon, B. (2006). Beta-2-microglobulin gene expression is maintained in rainbow trout and Atlantic salmon kept at low temperatures. Fish Shellfish Immunol. 21,176 -186.[CrossRef][Medline]
Lee, T. H., Tsai, J. C., Fang, M. J., Yu, M. J. and Hwang, P. P. (1998). Isoform expression of Na+/K+-ATPase alpha-subunit in gills of the teleost Oreochromis mossambicus. Am. J. Physiol. 275,R926 -R932.[Medline]
Lin, L. Y., Horng, J. L., Kunkel, J. G. and Hwang, P. P. (2006). Proton pump-rich cell secretes acid in skin of zebrafish larvae. Am. J. Physiol. 290,C371 -C378.[CrossRef]
Lin, T. Y., Liao, B. K., Horng, J. L., Yan, J. J., Hsiao, C. D. and Hwang, P. P. (2008). Carbonic anhydrase 2-like a and 15a are involved in acid-base regulation and Na+ uptake in zebrafish H+-ATPase-rich cells. Am. J. Physiol. 294,C1250 -C1260.[CrossRef]
Malek, R. L., Sajadi, H., Abraham, J., Grundy, M. A. and Gerhard, G. S. (2004). The effects of temperature reduction on gene expression and oxidative stress in skeletal muscle from adult zebrafish. Comp. Biochem. Physiol. 138,363 -373.
Nath, S., Kales, S., Fujiki, K. and Dixon, B. (2006). Major histocompatibility class II genes in rainbow trout (Oncorhynchus mykiss) exhibit temperature dependent downregulation. Immunogenetics 58,443 -453.[CrossRef][Medline]
Ojcius, D. M., Delarbre, C., Kourilsky, P. and Gachelin, G. (1994). Major histocompatibility complex class I molecules and resistance against intracellular pathogens. Crit. Rev. Immunol. 14,193 -220.[Medline]
Okabe, M. and Graham, A. (2004). The origin of
the parathyroid gland. Proc. Natl. Acad. Sci. USA
101,17716
-17719.
Pan, T. C., Liao, B. K., Huang, C. J., Lin, L. Y. and Hwang, P. P. (2005). Epithelial Ca2+ channel expression and Ca2+ uptake in developing zebrafish. Am. J. Physiol. 289,R1202 -R1211.
Piquemal, D., Joulia, D., Balaguer, P., Basset, A., Marti, J. and Commes, T. (1999). Differential expression of the RTP/Drg1/Ndr1 gene product in proliferating and growth arrested cells. Biochim. Biophys. Acta 1450,364 -373.[Medline]
Podrabsky, J. E. and Somero, G. N. (2004).
Changes in gene expression associated with acclimation to constant
temperatures and fluctuating daily temperatures in an annual killifish
Austrofundulus limnaeus. J. Exp. Biol.
207,2237
-2254.
Raynard, R. S. and Cossins, A. R. (1991). Homeoviscous adaptation and thermal compensation of sodium pump of trout erythrocytes. Am. J. Physiol. 260,R916 -R924.[Medline]
Rieder, C. L. and Cole, R. W. (2002). Cold-shock and the Mammalian cell cycle. Cell Cycle 1, 169-175.[Medline]
Roth, G. S., Lane, M. A., Ingram, D. K., Mattison, J. A., Elahi,
D., Tobin, J. D., Muller, D. and Metter, E. J. (2002).
Biomarkers of caloric restriction may predict longevity in humans.
Science 297,811
.
Schwarzbaum, P. J., Niederstatter, H. and Wieser, W. (1992a). Effects of temperature on the Na+/K+-ATPase and oxygen-consumption in hepatocytes of 2 species of fresh-water fish, roach (Rutilus rutilus) and brook trout (Salvelinus fontinalis). Physiol. Zool. 65,699 -711.
Schwarzbaum, P. J., Wieser, W. and Cossins, A. R. (1992b). Species-specific responses of membranes and the Na+/K+ pump to temperature-change in the kidney of 2 species of fresh-water fish, roach (Rutilus rutilus) and arctic char (Salvelinus alpinus). Physiol. Zool. 65, 17-34.
Shahsavarani, A. and Perry, S. F. (2006). Hormonal and environmental regulation of epithelial calcium channel in gill of rainbow trout (Oncorhynchus mykiss). Am. J. Physiol. 291,R1490 -R1498.
Shahsavarani, A., McNeill, B., Galvez, F., Wood, C. M., Goss, G.
G., Hwang, P. P. and Perry, S. F. (2006). Characterization of
a branchial epithelial calcium channel (ECaC) in freshwater rainbow trout
(Oncorhynchus mykiss). J. Exp. Biol.
209,1928
-1943.
Sohal, R. S. and Allen, R. G. (1986). Relationship between oxygen-metabolism, aging and development. Adv. Free Radical Biol. 2,117 -160.[CrossRef]
Sollid, J. and Nilsson, G. E. (2006). Plasticity of respiratory structures-adaptive remodeling of fish gills induced by ambient oxygen and temperature. Respir. Physiol. Neurobiol. 154,241 -251.[CrossRef][Medline]
Sonna, L. A., Fujita, J., Gaffin, S. L. and Lilly, C. M.
(2002). Invited review: Effects of heat and cold stress on
mammalian gene expression. J. Appl. Physiol.
92,1725
-1742.
Sonna, L. A., Kuhlmeier, M. M., Carter, H. C., Hasday, J. D.,
Lilly, C. M. and Fairchild, K. D. (2006). Effect of moderate
hypothermia on gene expression by THP-1 cells: a DNA microarray study.
Physiol. Genomics 26,91
-98.
Tang, S. J., Sun, K. H., Sun, G. H., Lin, G., Lin, W. W. and Chuang, M. J. (1999). Cold-induced ependymin expression in zebrafish and carp brain: implications for cold acclimation. FEBS Lett. 459,95 -99.[CrossRef][Medline]
Ton, C., Stamatiou, D. and Liew, C. C. (2003).
Gene expression profile of zebrafish exposed to hypoxia during development.
Physiol. Genomics 13,97
-106.
Tsai, J. C. and Hwang, P. P. (1998). The wheat germ agglutinin binding sites and development of the mitochondria-rich cells in gills of tilapia (Oreochromis mossambicus). Fish Physiol. Biochem. 19,95 -102.[CrossRef]
Vornanen, M., Hassinen, M., Koskinen, H. and Krasnov, A. (2005). Steady-state effects of temperature acclimation on the transcriptome of the rainbow trout heart. Am. J. Physiol. 289,R1177 -R1184.
Watabe, S. (2002). Temperature plasticity of
contractile proteins in fish muscle. J. Exp. Biol.
205,2231
-2236.
Wellmann, S., Buhrer, C., Moderegger, E., Zelmer, A., Kirschner,
R., Koehne, P., Fujita, J. and Seeger, K. (2004).
Oxygen-regulated expression of the RNA-binding proteins RBM3 and CIRP by a
HIF-1-independent mechanism. J. Cell Sci.
117,1785
-1794.
Yang, H., Wang, H., Czura, C. J. and Tracey, K. J.
(2005). The cytokine activity of HMGB1. J. Leukoc.
Biol. 78,1
-8.
Yu, J., Zhang, L., Hwang, P. M., Rago, C., Kinzler, K. W. and
Vogelstein, B. (1999). Identification and classification of
p53-regulated genes. Proc. Natl. Acad. Sci. USA
96,14517
-14522.
![]()
CiteULike
Complore
Connotea
Del.icio.us
Digg
Reddit
Technorati
Twitter What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||