Expression and water calcium dependence of calcium transporter isoforms in zebrafish gill mitochondrion-rich cells
BMC Genomics volume 8, Article number: 354 (2007)
Freshwater fish absorb Ca2+ predominantly from ambient water, and more than 97% of Ca2+ uptake is achieved by active transport through gill mitochondrion-rich (MR) cells. In the current model for Ca2+ uptake in gill MR cells, Ca2+ passively enters the cytosol via the epithelium Ca2+ channel (ECaC), and then is extruded into the plasma through the basolateral Na+/Ca2+ exchanger (NCX) and plasma membrane Ca2+-ATPase (PMCA). However, no convincing molecular or cellular evidence has been available to support the role of specific PMCA and/or NCX isoforms in this model. Zebrafish (Danio rerio) is a good model for analyzing isoforms of a gene because of the plentiful genomic databases and expression sequence tag (EST) data.
Using a strategy of BLAST from the zebrafish genome database (Sanger Institute), 6 isoforms of PMCAs (PMCA1a, PMCA1b, PMCA2, PMCA3a, PMCA3b, and PMCA4) and 7 isoforms of NCXs (NCX1a, NCX1b, NCX2a, NCX2b, NCX3, NCX4a, and NCX4b) were identified. In the reverse-transcriptase polymerase chain reaction (RT-PCR) analysis, 5 PMCAs and 2 NCXs were ubiquitously expressed in various tissues including gills. Triple fluorescence in situ hybridization and immunocytochemistry showed the colocalization of zecac, zpmca2, and zncx1b mRNAs in a portion of gill MR cells (using Na+-K+-ATPase as the marker), implying a subset of ionocytes specifically responsible for the transepithelial Ca2+ uptake in zebrafish gills. The gene expressions in gills of high- or low-Ca2+-acclimated zebrafish by quantitative real-time PCR analysis showed that zecac was the only gene regulated in response to environmental Ca2+ levels, while zpmcas and zncxs remained steady.
The present study provides molecular evidence for the specific isoforms of Ca2+ transporters, zECaC, zPMCA2, and zNCX1b, supporting the current Ca2+ uptake model, in which ECaC may play a role as the major regulatory target for this mechanism during environmental challenge.
Ca2+ is an essential element for almost all organisms and plays comprehensive regulatory roles in eukaryotic cells; therefore, Ca2+ must be maintained within a narrow concentration range in organisms (cellular free [Ca2+]: 100 nM; plasma total [Ca2+]: 2~3 mM). In terrestrial vertebrates, the whole-body Ca2+ balance is mainly achieved by intestinal absorption and kidney reabsorption. Ca2+ uptake mechanisms in the mammalian kidney have been most extensively studied, while fish gills, a specialized organ in non-mammalian vertebrates, is another model for studying Ca2+ uptake mechanisms . In fish, gills are the main site (> 97% of the whole body) of Ca2+ uptake from the aquatic environment to maintain the Ca2+ balance , and skin takes the place of gills during early developmental stages when the gills are not yet developed and functioning [3–5]. Compared with terrestrial vertebrates, Ca2+ regulation mechanisms in aquatic animals are probably more complicated and challenging because of the dramatic fluctuations in ambient Ca2+ concentrations which can occur (seawater, around 10 mM; fresh water, 0.01~3 mM).
According to the current model in mammals, active and transcellular Ca2+ transport is carried out as a 3-step process . Following entry of Ca2+ through apical epithelial Ca2+ channels (ECaC, TRPV5, and/or TRPV6), Ca2+ is bound to calbindins that facilitate diffusion to the basolateral membrane, and then it is extruded via the plasma membrane Ca2+-ATPase (PMCA) and/or Na+/Ca2+ exchanger (NCX). In this way, net transepithelial Ca2+ absorption from the luminal compartment (or environment) to the plasma is accomplished. A model of the capability of physiological regulation has also been similarly proposed in fish gills . However, only very little molecular evidence is currently available to support this model in fish gills, which is a specialized organ for Ca2+ uptake. The ecac gene has recently been cloned and sequenced in fugu, zebrafish, and trout [4, 7, 8]. In zebrafish embryos, low-Ca2+ fresh water causes upregulation of the whole-body Ca2+ influx and zECaC expression in mitochondrion-rich (MR) cells of both gills and skin, providing molecular evidence for the role of ECaC in fish Ca2+ uptake . On the other hand, some biochemical and physiological studies have investigated the possibility of PMCA and NCX's involvement in fish gill Ca2+ uptake [9, 10]. Nevertheless, no convincing evidence has been presented demonstrating the existence and involvement of PMCA and NCX in fish Ca2+ uptake mechanisms.
Three NCX genes (SLC8A1, SLC8A2, and SLC8A3) and 4 PMCA genes (ATP2B1, ATP2B2, ATP2B3, and ATP2B4) have been identified in mammals so far. The fundamental work of PMCA is the highly regulated active extrusion of Ca2+ from cells for maintaining a gradient across the plasma membrane . PMCA belongs to the P-type primary ion transport ATPase superfamily with 10 transmembrane domains and 1 calmodulin-binding site. At least 20 alternative splicing variants of PMCA have been found in mammals . In mammals, PMCA1 and PMCA4 are ubiquitously expressed, whereas PMCA2 and PMCA3 are more tissue specific; it has been suggested that PMCA1 and PMCA4 are housekeeping isoforms involved in the maintenance of cellular Ca2+ homeostasis . However, no molecular data about PMCAs are available in fish. On the other hand, at least 32 alternatively spliced isoforms of NCX1 gene products have been identified . Knock-out of NCX1 in mice is lethal due to cardiac defects before birth , while many physiological studies have indicated that vitamin D3 modulates Ca2+ balance by upregulating NCX1 expression . Recently, NCX1-3 and NCX4, another putative prototype of NCX, have been reported in fish, indicating the existence of NCXs in non-mammals as well [17, 18]. However, there is no molecular physiological evidence to support the involvement of NCX in fish Ca2+ uptake mechanisms.
Because of a theoretical process of whole-genome duplication (WGD) from ancestor vertebrates [19, 20], the numbers of gene isoforms found in teleosts are usually more than those found in mammals and tetrapods. While different isoforms coded from distinct genomic loci tend to have tissue-specific expressions, there are some cases of isoforms being differentially regulated upon facing environmental changes (reviewed by Schulte . Isoforms also serve compensatory functions; upregulation of TRPV6 was found in the kidneys of TRPV5-deficient mice . In biochemical studies on fish, the gill activities of PMCA and NCX, which are altered based upon physiological needs and environmental changes, have been proposed as being associated with transepithelial Ca2+ transport [9, 23]. However, nothing is known about which specific isoforms are responsible for Ca2+ uptake mechanisms in teleost gills. Therefore, using a functional genomic approach to investigate all isoforms of PMCA and NCX in a species may provide both physiological and genomic insights into this issue.
In the present study, a strategy of a whole-genome survey was used to uncover any unidentified PMCA and NCX isoforms in zebrafish, and then these candidates were refined with the EST database and cloned. Six isoforms of PMCAs (PMCA1a, PMCA1b, PMCA2, PMCA3a, PMCA3b, and PMCA4) and 7 isoforms of NCXs (NCX1a, NCX1b, NCX2a, NCX2b, NCX3, NCX4a, and NCX4b) were identified. Moreover, we report that some duplication of PMCA and NCX occurred, and the tissue distributions and expression patterns of these isoforms were analyzed by reverse-transcriptase polymerase chain reaction (RT-PCR) and triple fluorescence labeling, respectively. The gene expressions in gills of zebrafish acclimated to high- or low-Ca2+ environments were also examined by quantitative real-time PCR analysis. The present study provides molecular evidence to support the current Ca2+ uptake model in fish gill cells .
Molecular cloning and bioinformatics analysis of pmca and ncx
In this whole-genome survey, 6 distinct PMCAs (Figure 1) and 7 distinct NCXs (Figure 2) were found in zebrafish. According to the phylogenetic analysis, these transporters were annotated as zPMCA1a, zPMCA1b, zPMCA2, zPMCA3a, zPMCA3b, zPMCA4 (Figure 3a, and in standard nomenclature ATP2B1a, ATP2B1b, ATP2B2, ATP2B3a, ATP2B3b, and ATP2B4, respectively) and zNCX1a, zNCX1b, zNCX2a, zNCX2b, zNCX3, zNCX4a, and zNCX4b (Figure 3b, in standard nomenclature SLC8a1a, SLC8a1b, SLC8a2a, SLC8a2b, SLC8a3, SLC8a4a, and SLC8a4b, respectively). zNCX1a and zNCX1b were previously described and named zNCX1h and zNCX1n, respectively , and zNCX2a, zNCX3, zNCX4a, and zNCX4b have also been published [17, 18]. However, zNCX2b and the other zPMCAs are first described in this study.
Based on the topological structures of PMCA and NCX summarized in mammals [12, 26], zebrafish PMCAs and NCXs share similar patterns as mammals at the deduced amino acid sequences. According to the hydropathy analysis, 10 putative transmembrane domains were predicted for zPMCA (Figures 1, 4a), whereas NCX has 11 hydrophobic segments and 9 putative transmembrane domains (Figures 2, 4b) due to an intracellular helix front transmembrane domain 6 and a P loop-like structure of hydrophobic segment 8, which contained the amino acid sequence "GIG" . PMCA and sarco-endoplasmic reticulum calcium ATPase (SERCA), the closest kin to PMCA, share a similar topological structure: hydrophobic segments span the membrane 10 times, and 3 intracellular loops consist of Ca2+ binding and transporting, and calmodulin-regulating functions [12, 27]; the zPMCAs were also found to have 3 cytoplasmic domains 1 each between TM2 and TM3, TM4 and TM5, and TM10 and the C-terminal (Figure 1). As for NCX, a large intracellular loop for ion translocation [28, 29] was located between TM5 and TM6 (Figure 2). There were 20 alternative splicing isoforms of zPMCAs found when cloning was conducted, while 9 alternative splicing isoforms were found in zNCXs (data not shown). The hot spots of alternative splicing regions in the zPMCAs are at the 5' untranslated regions (UTRs), the fragments around the 300th amino acid including from 1 to 3 small exons (less than 100 bp each in those cases), and 3' amino acid tails that usually stretch to the 3' UTR and alter the location of stop codons. On the other hand, a region of alternative splicing sites was found in zNCX1b and zNCX3 around the 610th amino acid, which is a part of the intracellular loop.
zPMCAs share 63%~86% identities with each other, while zNCXs exhibit higher variety among isoforms (i.e., 36%~76% identity). According to the phylogenetic analysis, zebrafish and 2 pufferfish species so far all have 6 PMCA and 7 NCX isoforms (Figure 3a,b). Compared with the isoforms annotated in tetrapods, duplication from the ancestor ray-finned fish occurred in PMCA1, PMCA3, NCX1, NCX2, and NCX4, according to the latest version of the genomic database utilized. Moreover, in the karyotype, zPMCAs and zNCXs are located on different chromosomes (see additional file 1), suggesting that these paralogous groups may have originated from genome duplication but not from tandem duplication. Relationships among these pairs of duplicates were evident not only in the phylogenetic tree, but also in their gene structures (Figure 4). All teleost PMCAs and NCXs examined in this study showed an outgroup topology to their orthologues of mammals (Figure 3a,b), and the paralogous teleosts formed an inner-group topology. The 2 prototypes of NCXs, NCX4a and NCX4b, found in teleosts were most closely related to NCX1 and NCX3 in both phylogenetic trees (Figure 3) and gene structures (Figure 4b). In gene structures, all of them had a relatively large first exon, even though some small introns were found in the split exons in zNCX4s (Figure 4b).
mRNA expressions of pmca and ncx in various tissues
RT-PCR analysis revealed that zpmca2 and zpmca4 were ubiquitously expressed in various tissues compared with zpmca1a, zpmca1b, zpmca3a, and zpmca3b (Figure 5). All zpmcas were abundant in the brain and eye except zpmca1b, which was mainly expressed in the eye. These findings agree with the Ca2+ homeostasis function of PMCAs found in the mammalian neuron system . zncx1b and zncx4a were also ubiquitously expressed, while zncx2b and zncx3 were only detected in the brain and eyes. zpmca1a, zpmca1b, zpmca2, zpmca3a, and all the zncxs were expressed in ovary and 1-cell embryo (Figure 5), suggesting that these genes may regulate intracellular Ca2+ homeostasis at earlier stages. zpmca1a, zpmca2, zpmca4, zncx1, and zncx4a were expressed in gills (Figure 5).
Effects of environmental Ca2+ levels on ecac, pmca, and ncx mRNA expressions
According to the quantitative real time-PCR analysis, low-Ca2+ FW-acclimated zebrafish expressed significantly higher levels of zecac in gills than did high-Ca2+ FW-acclimated fish (Figure 6). However, no significant difference was found between the 2 groups in the expressions of any of the zpmcas and zncxs (Figure 6). Consistent with the results of the tissue scans shown in Figure 5, zpmca4, zpmca1a, and zpmca2 exhibited more than 20-fold higher expression levels than did the other zpmca isoforms in gills, whereas zncx1b was found to have an approximately 5-times stronger signal compared with the other zncxs (Figure 6).
Colocalization of ecac, pmca, and ncx mRNAs and Na+-K+-ATPase
Fluorescence in situ hybridization and immunocytochemistry were used to determine the isoforms of zpmca and zncx that were specifically expressed in gill ionocytes. Specific mRNA probes of all of the zpmca and zncx isoforms were used to conduct the in situ hybridizations for whole-mount embryos and adult gills. Among the isoforms, only zpmca2 and zncx1b mRNA signals were found in specific groups of cells in embryonic skin and adult gills (data not shown).
Subsequently, double-fluorescence in situ hybridizations were used for zecac and zpmca2 (Figure 7A–C), and zecac and zncx1b (Figure 7D–F), respectively. More gill cells expressed zpmca2 and zncx1b mRNAs than zecac mRNA (Figure 7A–F). Notably, zpmca2 and zncx1b were also abundantly expressed in gill lamellae, where fewer MR cells are usually found . zecac and zpmca2 mRNAs were colocalized only in a portion of gill cells, but were not colocalized all of the time. Around half of the zecac-positive cells co-expressed zpmca2 signals (51.7 +/- 3.50%, n = 3), while only 40.3 +/- 1.35% (n = 3) of zpmca2-positive cells showed zecac signals (Figures 7A–C, 8). Similar results were found in the case of double labeling of zecac and zncx1b (Figures 7D–F, 8). The proportion of zncx1b- positive cells that co-expressed zecac was 77.1 +/- 1.35% (n = 3), while zncx1b signals could be detected in all zecac-positive cells (Figures 7D–F, 8).
The subsequent triple-labeling experiments further demonstrated the colocalization of zecac mRNA/zpmca2 mRNA/Na+-K+-ATPase (Figure 8A–E), and zecac mRNA/zncx1b mRNA/Na+-K+-ATPase (Figure 8F–J), respectively, in specific groups of gill ionocytes. Taken all together, only a portion of gill ionocytes co-expressed the mRNA of zecac, zpmca2, and zncx1b, and Na+-K+-ATPase.
So far, 4 PMCA genes have been identified in mammals, and PMCA1 and PMCA4 were suggested to be housekeeping isoforms. Mammalian PMCA2 was found to be expressed in the hair bundle of hair cells of the inner ear , suggesting that it is responsible for endolymph Ca2+ homeostasis . A null mutation of PMCA2 caused balance and hearing defects in mice . On the other hand, mammalian PMCA1 was deemed a major Ca2+ extrusion transporter for Ca2+ absorption in intestine and kidney [32–34]. Our data demonstrated opposite phenomena in fish PMCAs. zpmca2 was expressed with a housekeeping pattern (Figure 5), while zpmca1a showed a hair cell expression pattern during embryonic stages (data not shown). This functional substitution of paralogues found in Ca2+ transporters is not a unique example. In zebrafish, the expression profiles suggest the hypothesis that zebrafish hoxB1a and hoxB1b are functional equivalents of mouse Hoxb1 and Hoxa1, respectively. On the contrary, the zebrafish orthologue of mouse Hoxa1, zebrafish hoxA1a, is not expressed in a similar pattern . Further examination of the promoter regions of these genes may provide insights into this phenomenon. As far as our data, zPMCA2 was proposed to participate in the transcellular Ca2+ uptake pathway in fish gills. NCX1a and NCX1b from teleosts including zebrafish and pufferfish were not clustered within an inner group on the phylogenetic tree (Figure 3). Further collection is necessary to enlarge the sequence pool, which may allow us to reconstruct the phylogenetic tree or to rename members of the gene family. In the present study, zNCX1b was demonstrated to be responsible for the Ca2+ extrusion mechanism in fish gill cells, while the paralogue, zNCX1a, was proven to be a heart-specific isoform that causes cardiac fibrillations in zNCX1a mutants . The combined functions from both zNCX1s are equivalent to those of mammalian NCX1; however, investigations of the direct function of NCX1 in Ca2+ absorption were not possible due to fatal heart failure in knockout mice. Therefore, teleosts may serve as a good model to examine the role of NCX in Ca2+ uptake without resulting heart failure.
As described above, Flik et al.  proposed a model of a fish gill Ca2+uptake mechanism, which is similar to that for Ca2+ reabsorption in mammalian kidneys . There has been no convincing or comprehensive molecular evidence to support the existence of the 3 major transporters, ECaC, PMCA, and NCX, in fish gill ionocytes until the present study. Pan et al.  for the first time demonstrated the expression of zecac in a subset of MR cells in zebrafish skin/gills and a correlation between zecac expression and in vivo Ca2+ uptake function. Subsequently in rainbow trout, Perry and colleagues found the expression of ECaC mRNA and/or protein in both pavement and MR cells in rainbow trout [8, 35], which is inconsistent with the current view that MR cells are the predominant site for branchial Ca2+ uptake , and the authors also reported that hypercapnia, implantation with cortisol, and infusion with CaCl2 affected ECaC expression in MR cells and/or pavement cells . In a recent study on isolated trout gill cells, PNA+ MR cells showed an over 3-fold higher 45Ca2+ uptake capacity compared to either PNA- MR or pavement cells . The subset of gill MR cells and/or pavement cells that expressed ECaC is not necessarily involved in transepithelial Ca2+ uptake as Hwang and Lee  claimed, because those previous studies lacked molecular and physiological data to demonstrate the transport pathway across the basolateral membrane. The present triple in situ hybridization and immunocytochemical experiments demonstrate the co-expressions of zecac, zpmca2, and zncx1b in a specific group of gill MR cells, providing comprehensive molecular evidence for the involvement of these transporters in a fish gill Ca2+ uptake mechanism.
Double in situ hybridization and immunocytochemistry experiments indicated that about 90% of gill MR cells expressed zecac mRNA, although this may have been underestimated because the antigenicity of Na+-K+-ATPase (the marker for MR cells) is attenuated during the high-temperature treatments used for in situ hybridization. All zecac-expressing cells co-expressed zncx1b, and about half of the zecac-expressing cells co-expressed zpmca2, indicating that about 50% of zecac-expressing cells co-expressed both zncx1b and zpmca2. Taking all these findings into account, only about 45% of the entire population of MR cells co-express zecac, zpmca2, and zncx1b, and this subset of ionocytes may specifically carry out the function of transepithelial Ca2+ uptake in zebrafish gills. It was noted that some gill cells express only 1 or 2 (never all) of the 3 major Ca2+ transporters. These cells may be in the process of terminal differentiation as Hsiao et al.  reported in skin ionocytes in zebrafish embryos. Alternatively, the Ca2+ transporter(s) expressed in those cells may be involved in intracellular Ca2+ homeostasis or other cellular events [4, 39, 40]. This remains to be confirmed in future studies.
ECaC is a vitamin D-sensitive transporter in mammals , and the expression of PMCA1 was also found to be stimulated by 1,25(OH)2D3 in mammals and chicken [34, 41, 42]. 1,25(OH)2-vitamin D3 was shown to have a hypercalcemic effect on both freshwater and saltwater teleosts [6, 43], whereas stanniocalcin, a predominant hypocalcemic hormone produced by Stannius bodies, was found to show opposite effect on fish . Stanniocalcin did not show a significant effect on basolateral Ca2+ transport , but quickly reduced the permeability of Ca2+ at the apical membrane [6, 46]. It has been proposed that ECaC is the rate-limiting step and the gatekeeper channel for active Ca2+ transport . The present study provides further molecular physiological evidence to support this notion. The present study for the first time demonstrates that acclimation to a low-Ca2+ environment, which would stimulate Ca2+ uptake capacity in fish [3, 48, 49], caused upregulation of mRNA expression of zECaC but not zPMCA2 or zNCX1b, implying that the steady-state expressions of zPMCA2 and zNCX1b may fulfill the requirement for the transepithelial transport machinery under all situations. Moreover, this study provides molecular evidence to support a previous biochemical study by Flik et al. , who found that the basolateral Ca2+ extrusion mechanisms by PMCA and NCX is far below their maximum capacity in fish gills based on the kinetic properties of the 2 enzymes.
The major findings of the present study were that (1) 6 PMCAs and 7 NCXs were identified from zebrafish tissues; (2) differential expressions of these genes were found in various tissues including gills; (3) acclimation to low-Ca2+ environment stimulated the mRNA expression of ecac in gills, but not those of pmcas or ncxs; and (4) only a portion of gill MR cells co-expressed zecac, zpmca2, and zncx1b mRNAs.
Animals and embryos
Zebrafish (Danio rerio) brood stocks at the Institute of Cellular and Organismic Biology, Academia Sinica were kept in fresh water (local tap water, FW) at 28.5°C under a 14-h:10-h light: dark photoperiod. Embryos were collected within 30 min after fertilization and incubated in a Petri dish until the desired developmental stages. For whole-mount in situ hybridization, PTU (1-phenyl 2-thiourea) at a final concentration of 0.003% was added to prevent melanogenesis. Fish were anesthetized with buffered MS222 before sampling following the guidelines of the Academia Sinica Institutional Animal Care and Utilization Committee (approval no.: RFiZOOHP2006086).
Following a previous process in our laboratory , high-Ca2+ (2.00 mM, control) and low-Ca2+ (0.02 mM) artificial fresh waters were prepared with double-deionized water (Milli-RO60, Millipore, Billerica, MA, USA) supplemented with adequate CaSO4·2H2O, MgSO4·7H2O, NaCl, K2HPO4, and KH2PO4. Nominal Ca2+ concentrations of the high- and low-Ca2+ media were 2.00 and 0.02 mM, respectively, but the other ion concentrations of the media were the same ([Na+], 0.5 mM; [Mg2+], 0.16 mM; and [K+], 0.3 mM) as those in the local tap water. Variations in the ion concentrations were maintained within 10% of the predicted values by examination with an atomic absorption spectrophotometer (Hitachi Z-8000, Tokyo, Japan). Zebrafish were transferred to high- and low-Ca2+ media, respectively, for 2 wk. At the end of acclimation, fish gills were sampled for quantitative real-time PCR analysis.
Molecular cloning and sequences analysis
The peptide sequences from other species (teleost had a higher priority) were used to BLAST the genome databases (of NCBI and Ensembl) for zebrafish, pufferfish (Fugu rubripes), and tetradon (Tetraodon nigroviridis). The putative full-length or partial open reading frames of PMCAs and NCXs obtained were confirmed using the EST database, and/or were used to design primers for cloning and the RT-PCR analysis. PCR products thus obtained were subcloned into a pGEM-T Easy vector (Promega, Madison, WI, USA), and the nucleotide sequences were determined with an ABI 377 sequencer (Applied Biosystems, Warrington, UK). Sequence analysis was conducted with a BLASTx program (NCBI). The specific primers of 5' and 3'-rapid amplification of cDNA ends (RACE) were designed from the partial sequences obtained from the PCR with degenerate primers. The RACE PCR program followed a commercial protocol (Clontech, Mountain View, CA, USA), and the RACE PCR products were also subcloned into pGEM-T Easy vectors and sequenced. The full-length deduced amino-acid sequences were aligned with ClustalX and then phylogenetically analyzed with Mega3.1 . For the proteomic structure analysis, the hydropathy plot and transmembranes were predicted with EMBOSS (the European Molecular Biology Open Software Suite, version 2.8.0  algorithm according to Persson and Argos .
An appropriate amount of zebrafish adult tissues was collected for the total RNA preparation with TRIZol reagent (Invitrogen, Carlsbad, CA, USA). The amount and quality of total RNA were determined by measuring the absorbances at 260 and 280 nm with a spectrophotometer (Hitachi U-2000) and RNA denatured gels. For cDNA synthesis, 5 μg of total RNA was reverse-transcribed in a final volume of 20 μL containing 0.5 mM dNTPs, 2.5 μM oligo (dT)18, 5 mM dithiothreitol, and 200 units superscript reverse transcriptase III (Invitrogen) for 1.5 h at 42°C, followed by a 15-min incubation at 70°C. For the PCR amplification, 1 μL cDNA was used as a template in a 25-μL final reaction volume containing 0.25 mM dNTP, 1.25 units Gen-Taq polymerase (Genemark, Taipei, Taiwan), and 0.2 μM of each primer. The primer sets for the PCR of tissue scans were zPMCA1a (428-bp fragment) forward 5'-AGACAGGGTGGATAGAAGGT-3', reverse 5'-CCCCAATAATGTGAAGATGA-3'; zPMCA1b (556-bp fragment), forward 5'-AGCCCCTTATTTCCCGCAC-3', reverse 5'-GCCCCCTTCTCAGCTCCATT-3'; zPMCA2 (444-bp fragment) forward 5'-GGGATCGGGATGAGATGGTA-3', reverse 5'-GCGTGTGTTTGTCTGTGGGT-3'; zPMCA3a (457-bp fragment) forward 5'-TTGCTGGTACAGATGTGGC-3', reverse 5'-GGAGGAGAGTGAAGCGGAG-3'; zPMCA3b (614-bp fragment) forward 5'-GGTCTGCTGGGTTTCCTACT-3', reverse 5'-TCCTGCTCAATCTCTCCTTT-3'; zPMCA4 (388-bp fragment) forward 5'-GTCAAGGCTGTGATGTGGGC-3', reverse 5'-TGGTTGGGAGTGGAGAAGGG-3'; zNCX1b (358-bp fragment) forward 5'-AGAGACGAGGAGAAGGAGGT-3', reverse 5'-GCACGAAAGCAAAGAGAACT-3'; zNCX2a (276-bp fragment) forward 5'-TGATGAAGAAGGAGGTGAAC-3', reverse 5'-CTTGCGAATGTGTCTGGTAT-3'; zNCX2b (477-bp fragment) forward 5'-GCAGCCGCATTTCCTTCG-3', reverse 5'-CGCCAGAGTTTCGGACCAC-3'; zNCX3 (244-bp fragment) forward 5'-GAGGAAGCAAGAAGAATAGC-3', reverse 5'-AGTCAAAACAAGATGGCAGA-3'; zNCX4a (462-bp fragment) forward 5'-GGAAGAGAGCGGAGAAGAGC-3', reverse 5'- ATGGTGAAGAGGGTGACGGA-3'; and zNCX4b (350-bp fragment) forward 5'-ATGCAGCAGCAGAAGAGC-3', reverse 5'-ACATTCTCGCCAGGTTTG-3'. A 514-bp fragment of zebrafish beta-actin was used as an internal control to evaluate the relative amounts of complementary DNAs (cDNAs), and the design of the primer pairs followed Hsiao et al. . The amplicons were all sequenced to make sure that the PCR products were the desired gene fragments.
Fluorescence double in situ hybridization
Fragments of the target genes obtained by PCR were inserted into pGEM-T Eeasy vectors. Digoxigenin- (Dig) (Roche, Penzberg, Germany) or dinitrophenol (DNP)-labeled (Perkin-Elmer, Boston, MA, USA) RNA probes were synthesized by in vitro transcription with T7 and SP6 RNA polymerase (Takara, Shiga, Japan). The qualities of the probes were examined using RNA gels, and the concentrations were determined by a dot-blot assay with standard DIG-labeled RNA (100 ng/μl) (Roche). Excised gills were fixed with 4% paraformaldehyde overnight at 4°C and then washed several times with PBS. For the cryosections, fixed samples were immersed in PBS containing 30% sucrose overnight, and embedded in Optimal Cutting Temperature (OCT) compound embedding medium (Sakura, Tokyo, Japan) at -20°C, and 10-μm frozen cross-sections were cut with a CM 1900 rapid sectioning cryostat (Leica, Heidelberg, Germany) and attached to poly-L-lysine coated slides (Electron Microscopy Sciences, Ft. Washington, PA, USA). Prepared samples from either cryosections or methanol-dehydrated whole gills were washed several times with phosphate-buffered saline with 0.1% tween-20 (PBST). After a brief washing with PBST, samples were incubated with hybridization buffer (HyB) containing 50% formamide, 5× SSC, and 0.1% Tween-20 for 5 min at 65°C. Prehybridization was performed for 2 h at 65°C with HyB+, which is the hybridization buffer supplemented with 500 ng/mL yeast tRNA and 50 μg/mL heparin. For hybridization, samples were incubated in 100 ng of the RNA probe in 200 μL HyB+ at 65°C overnight. Then, slides were washed at 65°C for 10 min in 75% HyB and 25% 2× SSC, 10 min in 50% HyB and 50% 2× SSC, 10 min in 25% HyB and 75% 2× SSC, 10 min in 2× SSC, and 30 min for 2 times in 0.2× SSC at 70°C. Further washes were performed at room temperature for 5 min in 75% 0.2× SSC and 25% PBST, 5 min in 50% 0.2× SSC and 50% PBST, 5 min in 25% 0.2× SSC and 75% PBST, and 5 min in PBST. Fluorescence staining was conducted with a commercial kit, TSA Plus Fluorescence Systems (Perkin-Elmer). The hybridization signals detected by the DIG-labeled RNA probes were amplified through fluorescein-tyramide signal amplification (TSA), while cyanine 3-TSA was used for the DNP-labeled probes. For triple labeling, samples were first subjected to in situ hybridizations of the 2 transporter genes and then to Na+/K+-ATPase immunocytochemistry (see below).
Zebrafish embryos and gills were fixed in 4% paraformaldehyde for 10 min at 4°C. After washing in PBS, fixed embryos were treated with 100% ethanol for 10 min at -20°C and subsequently subjected to blocking with 3% BSA at room temperature for 30 min. Embryos were then incubated with 1:200 PBS-diluted mouse anti-chicken Na+/K+-ATPase α subunit (cytosolic epitope for all isoforms) monoclonal immunoglobulin G (IgG) (Developmental Studies Hybridoma Bank, University of Iowa, Iowa City, IA, USA) at room temperature for 2 h. Samples were washed twice in PBS for 10 min each, and then incubated with 1:200 PBS-diluted goat anti-mouse IgG conjugated with Cy5 (Jackson Immunoresearch Laboratories, West Grove, PA, USA) for 2 h at room temperature. Images were acquired with a confocal laser scanning microscope (TCS-SP5, Leica Lasertechnik, Heidelberg, Germany).
Quantitative real-time PCR
Quantitative real-time PCR (qPCR) was carried out using a SYBR green dye (Applied Biosystems)-based assay with an ABI Prism 7000 Sequence Detection System (Applied Biosystems) according to the manufacturer's instructions. Primers targeting Ca2+ transporters (see additional file 1) and the endogenous control gene, beta-actin, were designed using the Primer Express 2.0 software (Applied Biosystems). In each assay, 25 ng cDNA was amplified in a 20-μL reaction containing 2x SYBR green master mix, 100 nM of forward and reverse primers, and nuclease-free water. Beta-actin, with a 151-bp amplicon (forward 5'-CACCTTCCAGCAGATGTGGA-3' and reverse 5'-AAAAGCCATGCCAATGTTGTC-3'), was used as an internal control to construct the standard curves.
Values are presented as the mean ± SE and were compared by the 2-sample t-test.
bovine serum albumin
epithelium Ca2+ channel
expressed sequence tag
Optimal Cutting Temperature
phosphate-buffered saline with 0.1% tween-20
plasma membrane Ca2+-ATPase
peanut lectin agglutinin
Quantitative real-time PCR
rapid amplification of cDNA ends
reverse-transcriptase polymerase chain reaction
sarco-endoplasmic reticulum calcium ATPase
SSC, saline-sodium citrate
transient receptor potential cation channel
fluorescein-tyramide signal amplification
Hoenderop JG, Nilius B, Bindels RJ: Calcium absorption across epithelia. Physioll Rev. 2005, 85: 373-422. 10.1152/physrev.00003.2004.
Flik G, Verbost PM, Wendelaar Bonga SE: Calcium transport process in fishes. Cellular and molecular approaches to fish ionic regulation. Edited by: Wood CM, Shuttleworth TJ. 1995, San Diego: Academic, 317-342.
Hwang PP, Tung YC, Chang MH: Effect of environmental calcium levels on calcium uptake in tilapia larvae (Oreochromis mossambicus). Fish Physiol Biochem. 1996, 15: 363-370. 10.1007/BF01875578.
Pan TC, Liao BK, Huang CJ, Lin LY, Hwang PP: Epithelial Ca2+ channel expression and Ca(2+) uptake in developing zebrafish. Am J Physiol Regul Integr Comp Physiol. 2005, 289: R1202-1211.
Hwang PP, Tsai YN, Tung YC: Calcium Balance in Embryos and Larvae of the Fresh-Water-Adapted Teleost, Oreochromis-Mossambicus. Fish Physio Biochem. 1994, 13: 325-333. 10.1007/BF00003437.
Flik G, Verbost PM, Wendelaar Bongar SE: Calcium transport process in fishes. Cellular and molecular approaches to fish ionic regulation. 1995, San Diego: Academic Press
Qiu A, Hogstrand C: Functional characterisation and genomic analysis of an epithelial calcium channel (ECaC) from pufferfish, Fugu rubripes. Gene. 2004, 342: 113-123. 10.1016/j.gene.2004.07.041.
Shahsavarani A, McNeill B, Galvez F, Wood CM, Goss GG, Hwang PP, Perry SF: Characterization of a branchial epithelial calcium channel (ECaC) in freshwater rainbow trout (Oncorhynchus mykiss). J Exp Biol. 2006, 209: 1928-1943. 10.1242/jeb.02190.
Verbost PM, Schoenmakers TJ, Flik G, Wendelaar Bonga SE: Kinetics of ATP- and Na+-gradient driven Ca2+ transport in basolateral membranes from gills of freshwater- and seawater-adapted tilapia. J Exp Biol. 1994, 186: 95-108.
van der Heijden AJ, Verbost PM, Bijvelds MJ, Atsma W, Wendelaar Bonga SE, Flik G: Effects of sea water and stanniectomy on branchial Ca2+ handling and drinking rate in eel (Anguilla anguilla L.). J Exp Biol. 1999, 202: 2505-2511.
Shull GE: Gene knockout studies of Ca2+-transporting ATPases. Eur J Biochem. 2000, 267: 5284-5290. 10.1046/j.1432-1327.2000.01568.x.
Strehler EE, Zacharias DA: Role of alternative splicing in generating isoform diversity among plasma membrane calcium pumps. Physiol Rev. 2001, 81: 21-50.
Guerini D: The significance of the isoforms of plasma membrane calcium ATPase. Cell Tissue Res. 1998, 292: 191-197. 10.1007/s004410051050.
Blaustein MP, Lederer WJ: Sodium/calcium exchange: its physiological implications. Physiol Rev. 1999, 79: 763-854.
Wakimoto K, Kobayashi K, Kuro OM, Yao A, Iwamoto T, Yanaka N, Kita S, Nishida A, Azuma S, Toyoda Y: Targeted disruption of Na+/Ca2+ exchanger gene leads to cardiomyocyte apoptosis and defects in heartbeat. J Biol Chem. 2000, 275: 36991-36998. 10.1074/jbc.M004035200.
Hoenderop JG, Dardenne O, Van Abel M, Van Der Kemp AW, Van Os CH, St-Arnaud R, Bindels RJ: Modulation of renal Ca2+ transport protein genes by dietary Ca2+ and 1,25-dihydroxyvitamin D3 in 25-hydroxyvitamin D3-1alpha-hydroxylase knockout mice. Faseb J. 2002, 16: 1398-1406. 10.1096/fj.02-0225com.
Marshall CR, Fox JA, Butland SL, Ouellette BF, Brinkman FS, Tibbits GF: Phylogeny of Na+/Ca2+ exchanger (NCX) genes from genomic data identifies new gene duplications and a new family member in fish species. Physiol Genomics. 2005, 21: 161-173. 10.1152/physiolgenomics.00286.2004.
Shu X, Huang J, Dong Y, Choi J, Langenbacher A, Chen JN: Na,K-ATPase alpha2 and Ncx4a regulate zebrafish left-right patterning. Development. 2007, 134: 1921-1930. 10.1242/dev.02851.
Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431: 946-957. 10.1038/nature03025.
Prince VE, Pickett FB: Splitting pairs: the diverging fates of duplicated genes. Nat Rev. 2002, 3: 827-837. 10.1038/nrg928.
Schulte PM: Changes in gene expression as biochemical adaptations to environmental change: a tribute to Peter Hochachka. Comp Biochem Physiol B Biochem Mol Biol. 2004, 139: 519-529. 10.1016/j.cbpc.2004.06.001.
Hoenderop JG, van Leeuwen JP, van der Eerden BC, Kersten FF, van der Kemp AW, Merillat AM, Waarsing JH, Rossier BC, Vallon V, Hummler E: Renal Ca2+ wasting, hyperabsorption, and reduced bone thickness in mice lacking TRPV5. J Clin Invest. 2003, 112: 1906-1914. 10.1172/JCI200319826.
Bijvelds M, Heijden A, Flik G, Verbost P, Kolar Z, Bonga S: Calcium pump activities in the kidneys of Oreochromis mossambicus. J Exp Biol. 1995, 198: 1351-1357.
Hwang PP, Lee T: New insights into fish ion regulation and mitochondria-rich cells. Comp Biochem Physiol A Mol Integr Physiol. 2007, 148: 479-497. 10.1016/j.cbpa.2007.06.416.
Langenbacher AD, Dong Y, Shu X, Choi J, Nicoll DA, Goldhaber JI, Philipson KD, Chen JN: Mutation in sodium-calcium exchanger 1 (NCX1) causes cardiac fibrillation in zebrafish. P Natl Acad Sci USA. 2005, 102: 17699-17704. 10.1073/pnas.0502679102.
Philipson KD, Nicoll DA: Sodium-calcium exchange: a molecular perspective. Annu Rev Physiol. 2000, 62: 111-133. 10.1146/annurev.physiol.62.1.111.
Verma AK, Filoteo AG, Stanford DR, Wieben ED, Penniston JT, Strehler EE, Fischer R, Heim R, Vogel G, Mathews S: Complete primary structure of a human plasma membrane Ca2+ pump. J Biol Chem. 1988, 263 (28): 14152-14159.
Iwamoto T, Nakamura TY, Pan Y, Uehara A, Imanaga I, Shigekawa M: Unique topology of the internal repeats in the cardiac Na+/Ca2+ exchanger. FEBS letters. 1999, 446: 264-268. 10.1016/S0014-5793(99)00218-5.
Nicoll DA, Ottolia M, Lu L, Lu Y, Philipson KD: A new topological model of the cardiac sarcolemmal Na+-Ca2+ exchanger. J Biol Chem. 1999, 274: 910-917. 10.1074/jbc.274.2.910.
Yamoah EN, Lumpkin EA, Dumont RA, Smith PJ, Hudspeth AJ, Gillespie PG: Plasma membrane Ca2+-ATPase extrudes Ca2+ from hair cell stereocilia. J Neurosci. 1998, 18: 610-624.
Kozel PJ, Friedman RA, Erway LC, Yamoah EN, Liu LH, Riddle T, Duffy JJ, Doetschman T, Miller ML, Cardell EL: Balance and hearing deficits in mice with a null mutation in the gene encoding plasma membrane Ca2+-ATPase isoform 2. J Biol Chem. 1998, 273: 18693-18696. 10.1074/jbc.273.30.18693.
Armbrecht HJ, Boltz MA, Kumar VB: Intestinal plasma membrane calcium pump protein and its induction by 1,25-OH2D3 decrease with age. Am J Physiol. 1999, 277: G41-47.
Cai Q, Chandler JS, Wasserman RH, Kumar R, Penniston JT: Vitamin D and adaptation to dietary calcium and phosphate deficiencies increase intestinal plasma membrane calcium pump gene expression. P Natl Acad Sci USA. 1993, 90: 1345-1349. 10.1073/pnas.90.4.1345.
Zelinski JM, Sykes DE, Weiser MM: The effect of vitamin D on rat intestinal plasma membrane Ca-pump mRNA. Biochem Bioph Res Co. 1991, 179: 749-755. 10.1016/0006-291X(91)91880-L.
Shahsavarani A, Perry SF: Hormonal and environmental regulation of epithelial calcium channel in gill of rainbow trout (Oncorhynchus mykiss). Am J Physiol Regul Integr Comp Physiol. 2006, 291: R1490-1498.
Evans DH, Piermarini PM, Choe KP: The multifunctional fish gill: dominant site of gas exchange, osmoregulation, acid-base regulation, and excretion of nitrogenous waste. Physiol Rev. 2005, 85: 97-177. 10.1152/physrev.00050.2003.
Galvez F, Wong D, Wood CM: Cadmium and calcium uptake in isolated mitochondria-rich cell populations from the gills of the freshwater rainbow trout. Am J Physiol Regul Integr Comp Physiol. 2006, 291: R170-176.
Hsiao CD, You MS, Guh YJ, Ma M, Jiang YJ, Hwang PP: A Positive Regulatory Loop between foxi3a and foxi3b is Essential for Specification and Differentiation of Zebrafish Epidermal Ionocytes. PLoS ONE. 2007, 21: e302-10.1371/journal.pone.0000302.
Prasad V, Okunade G, Liu L, Paul RJ, Shull GE: Distinct phenotypes among plasma membrane Ca2+-ATPase knockout mice. Ann NY Acad Sci. 2007, 1099: 276-286. 10.1196/annals.1387.029.
Reppel M, Sasse P, Malan D, Nguemo F, Reuter H, Bloch W, Hescheler J, Fleischmann BK: Functional expression of the Na+/Ca2+ exchanger in the embryonic mouse heart. J Mol Cell Cardiol. 2007, 42: 121-132. 10.1016/j.yjmcc.2006.10.008.
Wasserman RH, Smith CA, Brindak ME, De Talamoni N, Fullmer CS, Penniston JT, Kumar R: Vitamin D and mineral deficiencies increase the plasma membrane calcium pump of chicken intestine. Gastroenterology. 1992, 102: 886-894.
Pannabecker TL, Chandler JS, Wasserman RH: Vitamin-D-dependent transcriptional regulation of the intestinal plasma membrane calcium pump. Biochem Bioph Res Co. 1995, 213: 499-505. 10.1006/bbrc.1995.2159.
Marshall WS, Bryson SE: Transport mechanisms of seawater teleost chloride cells: an inclusive model of a multifunctional cell. Comp Biochem Physiol A Mol Integr Physiol. 1998, 119: 97-106. 10.1016/S1095-6433(97)00402-9.
Srivastav AK, Flik G, Wendelaar Bonga SE: Plasma calcium and stanniocalcin levels of male tilapia, Oreochromis mossambicus, fed calcium-deficient food and treated with 1,25 dihydroxyvitamin D3. Gen Compa Endocr. 1998, 110: 290-294. 10.1006/gcen.1998.7074.
Verbost PM, Flik G, Fenwick JC, Greco AM, Pang PKT, Bonga SEW: Branchial Calcium-Uptake – Possible Mechanisms of Control by Stanniocalcin. Fish Physiol Biochem. 1993, 11: 205-215. 10.1007/BF00004568.
Verbost PM, Butkus A, Atsma W, Willems P, Flik G, Bonga SE: Studies on stanniocalcin: characterization of bioactive and antigenic domains of the hormone. Mol Cell Endocrinol. 1993, 93: 11-16. 10.1016/0303-7207(93)90133-5.
Vennekens R, Hoenderop JG, Prenen J, Stuiver M, Willems PH, Droogmans G, Nilius B, Bindels RJ: Permeation and gating properties of the novel epithelial Ca2+ channel. J Biol Chem. 2000, 275: 3963-3969. 10.1074/jbc.275.6.3963.
Chang IC, Lee TH, Yang CH, Wei YY, Chou FI, Hwang PP: Morphology and function of gill mitochondria-rich cells in fish acclimated to different environments. Physiol Biochem Zool. 2001, 74: 111-119. 10.1086/319304.
Chen YY, Lu FI, Hwang PP: Comparisons of calcium regulation in fish larvae. J Exp Zool. 2003, 295: 127-135. 10.1002/jez.a.10195.
Flik G, Kaneko T, Greco AM, Li J, Fenwick JC: Sodium dependent ion transporters in trout gills. Fish Physiol Biochem. 1997, 17: 385-396. 10.1023/A:1007768825043.
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Briefings in bioinformatics. 2004, 5: 150-163. 10.1093/bib/5.2.150.
Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.
Persson B, Argos P: Prediction of transmembrane segments in proteins utilising multiple sequence alignments. J Mol Biol. 1994, 237: 182-192. 10.1006/jmbi.1994.1220.
Hsiao CD, Tsai WY, Horng LS, Tsai HJ: Molecular structure and developmental expression of three muscle-type troponin T genes in zebrafish. Dev Dyn. 2003, 227: 266-279. 10.1002/dvdy.10305.
This study was financially supported by grants to P.P.H. from the National Science Council, and Academia Sinica, Taiwan, R.O.C. We extend our thanks to Ms. Y. C. Tung and Mr. J. Y. Wang for their assistance during the experiments, and to the Core Facility of the Institute of Cellular and Organismic Biology, the Institute of Botany and Microbiology, Academia Sinica, for assistance with array printing and scanning, sequencing and microscopy.
PPH, the correspondence author, organized the whole project and manuscript. BKL designed and conducted the experiments and the data analysis; AND performed molecular cloning experiments; SCC, and MYC performed other biochemical experiments. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Primers for quantitative real-time PCR and the chromosome loci of zebrafish plasma membrane Ca2+-ATPase (PMCA) and Na+/Ca2+ exchanger (NCX). The chromosome number and loci were obtained by the SSAHA program from the Ensembl database (zebrafish genome assemble version 6). (XLS 31 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Liao, BK., Deng, AN., Chen, SC. et al. Expression and water calcium dependence of calcium transporter isoforms in zebrafish gill mitochondrion-rich cells. BMC Genomics 8, 354 (2007). https://doi.org/10.1186/1471-2164-8-354