- Open Access
Cloning of aquaporin-1 of the blue crab, Callinectes sapidus: its expression during the larval development in hyposalinity
Aquatic Biosystems volume 8, Article number: 21 (2012)
Ontogenetic variation in salinity adaptation has been noted for the blue crab, Callinectes sapidus, which uses the export strategy for larval development: females migrate from the estuaries to the coast to spawn, larvae develop in the ocean, and postlarvae (megalopae) colonize estuarine areas. We hypothesized that C. sapidus larvae may be stenohaline and have limited osmoregulatory capacity which compromises their ability to survive in lower salinity waters. We tested this hypothesis using hatchery-raised larvae that were traceable to specific life stages. In addition, we aimed to understand the possible involvement of AQP-1 in salinity adaptation during larval development and during exposure to hyposalinity.
A full-length cDNA sequence of aquaporin (GenBank JQ970426) was isolated from the hypodermis of the blue crab, C. sapidus, using PCR with degenerate primers and 5′ and 3′ RACE. The open reading frame of CasAQP-1 consists of 238 amino acids containing six helical structures and two NPA motifs for the water pore. The expression pattern of CasAQP-1 was ubiquitous in cDNAs from all tissues examined, although higher in the hepatopancreas, thoracic ganglia, abdominal muscle, and hypodermis and lower in the antennal gland, heart, hemocytes, ovary, eyestalk, brain, hindgut, Y-organs, and gill. Callinectes larvae differed in their capacity to molt in hyposalinity, as those at earlier stages from Zoea (Z) 1 to Z4 had lower molting rates than those from Z5 onwards, as compared to controls kept in 30 ppt water. No difference was found in the survival of larvae held at 15 and 30 ppt. CasAQP-1 expression differed with ontogeny during larval development, with significantly higher expression at Z1-2, compared to other larval stages. The exposure to 15 ppt affected larval-stage dependent CasAQP-1 expression which was significantly higher in Z2- 6 stages than the other larval stages.
We report the ontogenetic variation in CasAQP-1 expression during the larval development of C. sapidus and the induction of its expression at early larval stages in the exposure of hyposalinity. However, it remains to be determined if the increase in CasAQP-1 expression at later larval stages may have a role in adaptation to hyposalinity.
Ontogenetic variation in salinity tolerance and osmoregulatory capacity may be directly related to patterns of dispersal and recruitment of animals in various aquatic habitats. In decapod crustaceans, adults of the blue crab, Callinectes sapidus, the green shore crab, Carcinus maenas, and the American lobster, Homarus americanus[2, 3] are known to be strong hyper- and hypo-osmoregulators and inhabit a wide range of salinities. On the other hand, their embryonic and larval stages require high salinity water, possibly due to a limited osmoregulatory capacity . Consequently, larvae are typically exported to higher salinity waters for larval development either by migration of females prior to spawning or rapid transport of larvae out of estuaries during ebb tides.
In estuaries such as the Chesapeake Bay, life stage-dependent osmoregulatory capacity and salinity tolerance may be the driving force underlying population structures of C. sapidus, resulting from migration of adult females to high salinity waters for spawning and the return migration of postlarvae (megalopae). First, adult females migrate to higher salinity areas near coastal waters after the pubertal- terminal molting and mating, where they spawn and release pelagic larvae . These larvae largely spend seven-eight zoeal stages in coastal ocean waters . However, upon molting to the megalopa stage, they migrate back to the coast and invade lower salinity estuarine areas where they metamorphose to the first crab stage . Thus, the life cycle of C. sapidus presents a typical ontogenetic variation in salinity adaptation as osmoregulatory capacity and salinity tolerance are acquired during late larval development or the megalopal stage.
Salinity adaptation involves a complex process that entails dramatic changes in cell volume, ion transport, cellular metabolism, and whole-scale tissue remodeling. A large number of genes are involved in this osmoregulatory process in Carcinus maenas. The aquaporin (AQP) family of water channels, small and very hydrophobic intrinsic membrane proteins, is critical in the physiological processes of water and solute transport for salinity adaptation . AQPs are ubiquitous, being present in bacteria, plants, and animals. To date, 13 isoforms of the AQP family can be grouped into three subfamilies: aquaporins, aquaglyceroporins, and superaquaporins . Among these three subfamilies, the aquaporin subfamily including AQPs 0, 1, 2, 4, 5, 6, and 8 is selective for water transport .
The involvement of aquaporins in salinity adaptation has been most studied for AQP-1 in teleosts, although other aquaporins have been identified in this process. The expression of AQP-1 is found in most organs of fish with high expression in gill, intestine, and kidney, where the expression levels change in response to different salinities. Acclimation to hyposalinity up-regulated AQP 1a expression in the gill of Atlantic salmon and black porgy [12, 13]. On the other hand, acclimation to hypersalinity increased AQP-1 expression in the intestines and kidneys of Atlantic salmon , and in the intestines of European eels , Japanese eels , and sea bass . These studies demonstrate the potential for the involvement of aquaporins in the adaptation of C. sapidus to hyposalinity.
In view of the fact that adult females migrate to high salinity waters for spawning and that high salinity is required for larval development, we hypothesized that C. sapidus larvae may be stenohaline and have limited osmoregulatory capacity which compromises their ability to survive in lower salinity waters. We tested this hypothesis using hatchery-raised larvae that were traceable to specific life stages. In addition, we aimed to understand the possible involvement of AQP-1 in salinity adaptation during larval development and during exposure to hyposalinity. We present evidence that larval stages Z2-6 exposed to artificial seawater (ASW) at 15 ppt showed significantly higher expression of the blue crab aquaporin orthologue CasAQP-1 as compared to those exposed to 30 ppt. Molting percentage is much lower in ASW-exposed Z2-6 larvae as opposed to those reared at 15 ppt, suggesting that energy reserves are diverted to survival through osmoregulation under these conditions.
Sequence analyses of C. sapidus aquaporin 1 (CasAQP-1)
The nucleotide and deduced aa sequences of C. sapidus aquaporin-1 (CasAQP-1: GenBank JQ970426) are presented in Figure 1A. Both the 5′ and 3′ UTRs (italicized in Figure 1A) contained three terminal oligopyrimidine tracts (TOP) as a translation regulatory site: two located in the 5′ UTR and one located in the 3′ UTR (highlighted in bold and underlined in Figure 1A). The deduced amino acid sequence of CasAQP-1 does not contain the signal peptide (P = 0.068 by http://www.cbs.dtu.dk/services/SignalP). The Conserved Domain database (http://www.ncbi.nlm.nih.gov/Structure/cdd/) identified from the deduced amino a region from amino acid E15 to V214 of CasAQP-1 as a putative major intrinsic protein (MIP) superfamily member (Figure 1A, marked with arrows). The two highly conserved hydrophobic stretch regions, with two NPA boxes (boxed) that are involved in forming the water pore, are underlined. Four putative phosphorylation sites were predicted by NetPhos 2.0 Server (http://www.cbs.dtu.dk/services/NetPhos/) with a value >0.9 at three serine residues (S93, 199, and 224) and of 0.6 at one threonine residue (T119).
The 3D structure of CasAQP-1 (Figure 1B) was obtained using the 3D structure of c1ymgA (PDB) as the template showing that 91% of 217 deduced aa of CasAQP-1 were modeled with 100% confidence by the single highest scoring template. 10 α helices including six transmembrane helices, one β strand and 12 random coil secondary structures were predicted. Color rainbows indicates N- to C-termini of CasAQP-1.
A phylogram was generated with the deduced aa sequences of 15 different AQP-1 including seven vertebrates and 8 invertebrates (Figure 1C). The tree contains two separate clades: one for vertebrates and the other for invertebrates. CasAQP-1 was located close to AQP-1 of the cockroach, Blattella germanica, both of which were separate from the rest of the branch of insects and the water flea, Daphnia pulex.
Spatial expression of C. sapidus aquaporin I (CasAQP-1)
The expression pattern of CasAQP-1 in various tissue cDNAs prepared from a juvenile female C. sapidus at intermolt was examined with a pair of primers (CasAQP-1-3 F1 and -5R1, Table 1) amplifying the MIP domain (Figure 2). The ubiquitous CasAQP-1 expression was found in all tissues tested, although there were some differences in expression levels. The thoracic ganglia complex, hepatopancreas, abdominal muscle, and hypodermis had the highest levels of expression, followed by progressively lower levels in the antennal gland, heart, hemocytes, ovary, eyestalk, brain, hindgut, Y-organs, and gill.
Effects of salinity on molting and survival rate during the development of C. sapiduslarvae
The exposure of larvae to ASW at 15 ppt did not affect survival, but negatively affected the molting rate of C. sapidus larvae, as compared to controls held in ASW at 30 ppt (Figure 3A). In general, all larval stages (Z1-Z8) in ASW at 30 ppt showed molting rates varying from 12 ± 9% (Z1) to 47 ± 8% (Z4), in which Z1 had significantly lower molting rates than the rest of the larval stages. Larvae exposed to ASW at 15 ppt displayed a similar trend to those of the control but had much lower molting rates with 1 ± 1% at the Z1 stage, which increased significantly to a high rate of 13% (P <0.05) at Z2 onwards (Figure 3A). The statistical significance between the two groups is noted only at Z4 (P <0.001) and Z7 (P <0.05) larvae, due to a large variation at other larval stages.
Interestingly, however, there were no overall statistical differences in the survival of larvae among those exposed to ASW at 30 and 15 ppt (Figure 3B). At 30 ppt, the larvae at the first seven stages had 90- 97% survival, while Z8 to megalopae had survival of 78-97%. The larvae at Z1 to Z5 had 75-86% survival at 15 ppt, which was not significantly different from the controls with survival of 90-98%. The larvae at Z5-6 to Z8 also did not differ in survival at 15 ppt (83-97%) when compared to controls (78-97%).
Effect of salinity on the expression of CasAQP-1 and arginine kinase (CasAK) in C. sapiduslarvae
Experimental animals exposed to ASW at 30 ppt showed varying degrees of CasAQP-1 expression from 0.7 ± 0.4 to 6.0 ± 0.9 x 106 copies/μg total RNA (Figure 4). Larvae at stages Z1-2 and Z7-8 contained significantly greater expression of CasAQP-1: 5.1 ± 1.0 (n = 9) and 6.0 ± 0.9 x 106 copies/μg total RNA (n = 6), compared to the larvae at stage Z2-3 (0.7 ± 0.4 x 106 copies/μg total RNA, n = 7), Z3-4 (2.6 ± 0.5 x 106 copies/μg total RNA, n = 5), and Z5-6 (0.6 ± 0.2 x 106 copies/μg total RNA, n = 6).
Animals exposed to ASW at 15 ppt exhibited high levels of CasAQP-1 expression ranging from 11.6 ± 3.2 to 2.4 ± 0.9 x 106 copies/μg total RNA. Larvae at stage Z3-4 expressed the greatest amounts of CasAQP-1 with 11.6 ± 3.2 x 106 copies/μg total RNA (n = 5), which was followed by larvae at stage Z1-2 with 6.7 ± 1.7 x 106 copies/μg total RNA (n = 7). At stages Z2-3 and Z5-6, the larvae show markedly less expression of CasAQP-1 with 2.7 ± 0.7 x 106 (n = 7) and 2.4 ± 0.9 x 106 copies/μg total RNA (n = 5), respectively. The expression levels of CasAQP-1 at stage Z7-8 and Z8-megalopa were slightly higher with 3.3 ± 1.8 x 106 copies/μg total RNA (n = 6) and 3.6 ± 1.3 x 106 copies/μg total RNA (n = 6). These values did not significantly differ from those measured at stages Z2-3 and Z5-6. Larvae at stages Z2-3, Z3-4, and Z5-6 exposed to ASW at 15 ppt showed significantly higher expression of CasAQP-1 compared to those exposed to 30 ppt.
The expression of AK in the larvae exposed to ASW at 30 ppt was relatively consistent throughout the larval stages with a two-fold difference ranging from the lowest level of 0.9 ± 0.3 to the highest level of 1.9 ± 0.9 x 106 copies/μg total RNA (Figure 5). The expression of AK in the larvae exposed to ASW at 15 ppt was also similarly consistent throughout the larval stages with a three-fold difference with the lowest of 0.8 ± 0.4 to the highest amounts of 2.4 ± 1.2 x 106 copies/μg total RNA. Indeed, AK expression did not differ among the larval stages exposed to ASW at 30 ppt and 15 ppt.
In this study, we isolated the full length cDNA of AQP-1 from the hypodermis of C. sapidus and examined its mRNA expression level in various larval stages in response to hyposalinity.
CasAQP- 1 cDNA is a 5′ TOP mRNA with two TOP sites, one in the 5′ UTR and another in the 3′ UTR. These have been shown to be critical for translational control. The 5′ TOP found in ribosomal elongation factors is known to be the target of ripampicin, resulting in repression under the stressful, suboptimal growth conditions . Thus, the presence of a TOP site in the 5′UTR of CasAQP-1 implies that upstream signaling pathways and other trans-acting factors regulate its translation . The functional significance of 3′TOP is not yet understood.
The ORF of Callinectes AQP-1 encodes a deduced protein of 238 aa containing two NPA motifs forming the water pore. Four putative phosphorylation sites are predicted in the deduced CasAQP-1, whereas S228 is the only site, known to be phophorylated and involved in AQP trafficking [19, 20]. The expression pattern of CasAQP-1 is ubiquitous. In decapod crustaceans, the hepatopancreas and gills are known as the water uptake and osmoregulatory sites, respectively, and these tissues displayed differential expression of CasAQP-1 with higher expression in the former and lower expression in the latter. Due to the ubiquitous expression of CasAQP-1 in all tissues examined, and also its small size, we used whole larvae for determining its expression in response to exposure of larvae to different salinities.
The larvae exposed to a constant salinity of 30 ppt exhibited changes in CasAQP-1 expression changes during development, compared to a relatively consistent AK expression. Z1-2 and 7–8 larvae have the highest expression of CasAQP-1, while Z3-4 and 5–6 larvae showed the lowest expression. Our data imply that C. sapidus larvae may undergo endogenous changes in the expression level of CasAQP-1 during their larval development.
In addition, each of these larval stages shows differential responses when exposed to hyposalinity (15 ppt). Z2-3, Z3-4, and Z5-6 larvae respond by over-expression of CasAQP-1, whereas the rest of the larval stages remained at similar levels of expression, compared to the controls exposed to 30 ppt. Because our expression data were obtained after 96 hrs exposure to hyposalinity, we are not certain when the initial response of the CasAQP-1 upregulation occurs. In contrast to our data, AQP expression in the gills of the green crab, C. maenas, exposed to 10–15 ppt for 15 days steadily declined, and did not significantly differ from that of the control at 32 ppt . Interestingly, a similar discrepancy was noted in fish in that the expression of AQP 1 of the kidney of animals diverges in response to exposure to hypersalinity and hyposalinity. On the other hand, European eels and black porgy exposed to hypersalinity reduce AQP-1 expression in the kidney [21, 22], whereas Atlantic salmon increase its expression in the kidney .
AQPs belonging to a superfamily of MIPs channels passive permeation of water molecules across cellular membranes of bacteria, plants, and animals [23, 24]. Thirteen paralogs of AQPs found in mammals show often very specific to a particular cell type, while AQP 1, 3, 7 and 9 are found in various organs including kidney, red blood cells, eyes, ears and lungs . Water permeability is endowed by activation or deactivation of AQPs through phosphorylation or translocation in and out of the cell membrane . Most studies have focused on vertebrate AQPs, particularly euryhaline teleosts’ AQP 1 as they are capable of maintaining fluid homeostasis against fluctuating salinity conditions. In most fish, two paralogs of AQP 1 are found: AQP 1aa and 1ab in the kidney, although zebrafish kidney possesses only AQP 1 aa .
Much less is known about invertebrate AQPs, although putative AQP homologues are found in genomic databases. In insects, AQPs include the Drosophila integral protein (DRIPs) family with two NPA motifs specific for water transport , the Drosophila Big Brain gene (DmBiB) family with an extended C-terminal tail similar to human AQP 4, and the PRIP family closely related to DRIPs . Genomic data of Caenorhabditis elegans reveal the presence of eight AQPs  including three exclusive for water transport, one for glycerol, two for both water and glycerol, and two others not involved in the transport of any solutes examined . Considering the variety of AQPs in other invertebrates, we expect that C. sapidus also possesses more than one AQP.
A species-dependent ontogenetic variation in salinity tolerance can drive migratory patterns of dispersal and recruitment of animals. Juvenile and adult crustaceans with high osmoregulatory capacity often occur in estuarine conditions and appear to display two alternative strategies of dispersal and recruitment: retention and export strategies . Those species retaining their larvae in low salinity adult habitats such as Armases miersii, Sesarma curacaoense, Palaemonetes argentines, and Astacus leptodactylus show a strong osmoregluatory capacity during the larval stages. On the other hand, C. sapidus adopts an export strategy for its larvae with limited osmoreglatory capacity through development of larval stages in coastal shelf or oceanic waters . Consequently, later life stages (megalopae and juveniles) of C. sapidus migrate to low salinity areas and colonize estuaries.
The larval development of C. sapidus is complex, undergoing 7–8 zoeal stages. C. sapidus larvae in this study survived at 15 ppt water for 96 hrs but showed a stunted growth rate at early stages. This indicates that even larvae at earlier zoeal stages do have osmoregulatory capacity, but it is limited. Considering that osmoregulation is an energy-dependent process, it appears that larvae at earlier stages (Z1-4) exposed to 15 ppt compromise growth for survival. However, we are not certain that animals held under different conditions consumed the same amount of food. For this study, we did not directly measure the food intake of each animal, or whether hyposalinity might affect food consumption.
For C. sapidus larvae at later stages (Z5 to megalopae), growth increased by stage when larvae were exposed to hyposalinity at 15 ppt. This suggests that larvae at later stages that were exposed to hyposalinity utilized the energy gained from food consumption for growth whereas younger larvae may divert their energy largely for osmoregulation for survival, compromising their growth in the process. This difference in energy diversion could be critical to those larvae at earlier stages with smaller body sizes and hemolymph volumes. Our findings are consistent with reports that C. sapidus megalopae collected in the field exhibit a significant osmoregulatory capacity . On the other hand, our data contradict some prior studies which suggest that zoeae have limited osmoregulatory capacity and that metamorphosis results in the immediate appearance of adult-type osmoregulation .
In juveniles and adults of decapod crustaceans, gills and guts are the main tissues of iono-and osmoregulation that are mainly regulated by the CHH neuropeptide family . Animals at intermolt stage utilize pleiotrophic CHH and its isoform to increase Na+ influx through gill epithelial cells of Pachygraphsus marmoratus and hemolymph osmolarity and N+ influx in eyestalk-ablated crayfish, Astacus leptodactylus. On the other hand, animals during and immediately after ecdysis uptake iso-osmotic water through guts by drinking, the process of which is exerted by the release of gut CHH . The amount of water uptake seems closely associated with the level of CHH concentration in hemolymph, directly resulting in 20-50% molt-related somatic growth over a period of 1–3 hr.
Water uptake occurring during the molting process is recapitulated throughout the life cycle of crustaceans, starting from hatching that coincides with the first molting process. The onset of CHH expression occurs early in the developmental stage and its presence lasts throughout the life cycle [38, 42]. CHH mRNA in X-organ cells and neuropeptide in the sinus gland appears at an early eye anlage stage during embryogenesis . As opposed to juveniles and adults, however, the water uptake during hatching is also driven by the release of CHH but not from the endocrine cells located in the fore-and hindguts but those of the abdominal segments of embryos.
The expression pattern of CasAQP-1 is not determined during embryonic development prior to hatching. CasAQP-1 expression found in the larvae at Z1-2 implies a possible occurrence of its expression during embryonic development. The activation of vertebrate and mammalian aquaporins by phosphoryation is under the control of various hormones. Thus, it will be interesting to study whether the dramatic water uptake during hatching and molting may be driven by the cascade events of the release of CHH activating crustacean aquaporins.
We report the ontogenic variation in CasAQP-1 expression during the larval development of C. sapidus at 30 ppt and the induction of its expression at early larval stages in the exposure of hyposalinity. Hyposalinity (15 ppt) compromises only the growth rate during larval development, but not their survival rate. The expression of CasAQP-1 was common in all tissues obtained from a juvenile crab. It remains to be determined if the increase in CasAQP-1 expression at later larval stages may have a role in adaptation to hyposalinity.
Materials and methods
C. sapidus larvae were collected on the day of hatching and reared in a tank holding 1.5 m3 of artificial seawater (ASW) at 30 ppt at 22°C as described previously . During larval rearing, the density of larvae was ~100/L.
Abrupt exposure of larvae to the ASW at 30 and 15 ppt for 96 hrs
Acute static bioassays were employed to examine salinity tolerance. Two -seven larvae were placed directly into individual wells of a 24-well plate with each well containing 1 ml of filtered ASW at 30 and 15 ppt and kept for 96 hrs at 22-24°C. Larval stages Z1-Z4 were fed 10–15 rotifers/well/day. The larvae at Z5-Z8 were fed 10–15 newly-hatched fresh Artemia nauplii daily. The larvae at Z2-3 stages were fed 10–15 one-day old Artemia nauplii daily. Z7-8 larvae were fed two days old Artemia nauplii daily. The animals were monitored daily for survival and molting for 96 hrs, the time by which the most significant changes in gene expression have been noted in the gills of the green crab (Carcinus maenas) exposed to hyposalinity . During the exposure of 96 hrs, some of the larvae underwent molting. Therefore, larval stages were presented as 1–2, 2–3, 3–4, 5–6, 7–8, 8-megalopa in Figures 4 and 5. At the end of exposure, all the larvae retrieved from each well were gently blotted on a tissue paper and placed in a tube. Each larval stage was identified as described [44, 45]. Each treatment was hexaplicated and all experiments were replicated using larvae from a total of five independent spawns.
PCR with degenerate primers
Four degenerate primers were produced (Invitrogen) based on the conserved amino acid sequences of the following aquaporins found in GenBank (Dmel CG9023; AgaP 14229; hsa361_AQP4; bta_281008_AQP4; mmu_11829_AQP4; rno_25293_AQP4; mdo_100010466; gga_421088_AQP4; dre_445293; xla_443817; xla_495037; xtr_448309; dre_335821; dre_559284; hsa_359_AQP2; mcc_711719; cfa486552; bta_539870; xla_378655).
Total RNA was extracted from the tissues using TRIZOl by following the manufacturer’s protocol (Invitrogen) and quantified using a NanoDrop spectrometer (Fisher Scientific). Initially, the total RNA of the hypodermis of C. sapidus at premolt was subjected to 5′ and 3′ RACE cDNA synthesis using the GeneRacer® kit (Invitrogen) following the manufacturer’s protocol. The first amplicon (the expected size, 330 bp) was obtained using a two-step PCR method: 1) touch-down (TD)-PCR and 2) nested PCR, as described [46–48]. In brief, for the TD-PCR, hypodermis cDNA was amplified using Advantage Taq (BD Biosciences) with a combination of primers of AQPdF1 and CasAQPdR1. PCR conditions were as follows: initial denaturation at 94°C for 2.5 min; 94°C, 30 sec, 44°C, 30 sec, 72°C, 1 min, repeated for 3 cycles; 94°C, 30 sec, 44°C, 30 sec, 72°C, 1 min, repeated 3 cycles; 94°C, 30 sec, 42°C, 30 sec., 72°C, 1 min, repeated 3 cycles; 94°C, 30 sec, 40°C, 30 sec, 72°C, 1 min, repeated 3 cycles; 94°C, 30 sec, 45°C, 30 sec, 72°C, 1 min, repeated for 25 cycles and extended 7 min at 72°C. The TD-PCR products were diluted 10 times in sterilized water and then served as the template for a nested PCR with CasAQPdF2 and CasAQPdR2. Eppendorf Taq polymerase was used for the nested PCR at the following conditions: initial denaturation at 94°C for 2.5 min; 94°C, 30 sec, 50°C, 30 sec, 70°C, 1 min, repeated for 35 cycles and extended 7 min at 70°C. The procedures of DNA extraction and subcloning into a pGEM®-T Easy vector (Promega) for sequencing were described previously [46–48].
5′ and 3′ Rapid Amplification of cDNA Ends (RACE) of C. sapidusaquaporin
The TD-PCRs were performed with primers CasAQP-1-5R1 and CasAQP-1-3 F1 (Table 1) and the corresponding manufacturer’s primers respectively (Invitrogen), using conditions stated above except for the annealing temperatures: 54o, 52o, and 50o, at a final temperature of 55°C for 1 min extension. The nested PCRs for 5′ and 3′ RACE were carried out with 10-fold diluted TD-PCR product as template(s) and a combination of primers of CasAQP-1-5R2 and CasAQP-1-3 F2 (Table 1) with 5′ nested primer and 3′ nested primer of the manufacturer’s primers, respectively. The reaction was amplified at 58°C for annealing and 70°C for 1 min extension. The remaining cloning and sequencing procedures were as described above.
Spatial expression pattern of CasAQP-1in various tissue cDNAs of a juvenile female
Sample cDNAs were prepared from the dissection of a juvenile female at intermolt stage that was raised in 15 ppt salinity. Each tissue cDNA (1.5 μg total RNA) was diluted with sterilized water at the final concentration of 12.5 ng total RNA/ μl. One μl was amplified with CasAQP-1-3F1 and CasAQP-1- 5R1 (listed in Table 1) at 60°C annealing temperature for 35 cycles with Taq polymerase (Eppendorf). The arginine kinase (AK) gene was used as a reference gene, as reported elsewhere [47, 49]. After electrophoresis on a 1.5% agarose gel containing ethidium bromide, the gel was visualized and digitally-photographed using a Kodak-gel documentation system (Kodak).
The cDNA sequence was analyzed using an ORF finder program (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). The RNA regulatory motifs terminal oligopyrimidine tract (TOP) and upstream ORF (uORF) were predicted using RegRNA  (http://regrna.mbc.nctu.edu.tw). The signal peptide of the deduced aa of CasAQP was examined using SignalP 3.0 Server (http://www.cbs.dtu.dk/services/SignalP). The prediction of the phosphorylation sites was performed using NetPhos 2. 0 Server (http://www.cbs.dtu.dk/services/NetPhos/). Potential kinase specific phosphorylation was predicted using NetPhosK 1.0 Server (http://www.cbs.dtu.dk/services/NetPhosK/). The sequence homology was examined using the BLAST network server (blast.ncbi.nih.gov/Blast.cgi). Multiple protein sequences were aligned using ClustalW (http://www.genome.ad.jp).
The deduced aa sequence of CasAQP-1 was formatted into PDB sequences using the Phyre program (http://www.sbg.bio.ic.ac.uk/Phyre², c1mgA (PDB)  and viewed using the Jmol program (Version 1.2r3pre) for tertiary structure prediction.
A phylogram was constructed using the neighbor-joining method with 15 deduced amino acid sequences of aquaporin (http://www.phylogeny.fr) . The conserved domain was searched using the conserved domain database (http://www.ncbi.nlm.nih.gov/Structure/cdd) .
Quantitative RT-PCR (qRT-PCR) analysis
The expression of CasAQP-1 in the samples was determined using a qRT-PCR assay with each sample cDNA containing 25 ng total RNA. Primers for the assays are listed in Table 1. The qRT-PCR standard of CasAQP-1 and AK was prepared as described [47–49, 54–56]. The level of AK expression was examined as a reference gene with an AK standard that was generated similarly to that of CasAQP-1. The data were presented as mean ± 1SE copies/μg total RNA.
All results represent mean ± 1SE (n), in which n is the number of replicates. GraphPad InStat 3 program (GraphPad Software, Inc) was used to evaluate the statistical significance of the data. Statistical significance among the different larval stages at the exposure to 30 or 15 ppt was determined using one-way ANOVA with post-hoc Tukey-Kramer multiple comparison tests and was accepted at P <0.05. The difference of the same larval stage at two different salinities was calculated using Student’s t test.
Cieluch U, Anger K, Aujoulat F, Buchlolz F, Charmantier-Daures M, Charmantier G: Ontogeny of osmoregulatory structures and functions in the green crab Carinus maenas (Crustacea, Decapoda). J Exp Biol. 2004, 207: 325-336. 10.1242/jeb.00759.
Charmantier G, Charmantier-Daures M, Bouaricha N, Thuet P, Trilles J-P, Aiken DE: Ontogeny of osmoregulation and salinity tolerance in two decapod crustaceans: Homams americanus and Penaeus japonicus. Biol Bull. 1988, 175: 102-110. 10.2307/1541897.
Charmantier G, Haond C, Lignot J-H, Charmantier-Daures M: Ecophysiological adaptation to salinity throughout a life cycle: a review in homarid lobsters. J Exp Biol. 2001, 204: 967-977.
Anger K, Spivak E, Luppi T: Effects of reduced salinities on development and bioenergetics of early larval shore crab. Carcinus maenas J Exp Mar Biol Ecol. 1998, 220: 287-304. 10.1016/S0022-0981(97)00110-X.
Van Engel WA: The blue crab and its fishery in Chesapeake Bay: Part I. Reproduction, early development, growth, and migration. Commer Fish Rev. 1958, 20: 6-17.
Costlow JD, Bookhout CG: The larval development of Callinectes sapidus Rathbun reared in the laboratory. Biol Bull. 1959, 116: 373-396. 10.2307/1538947.
Epifanio CE, Garvine RW: Larval transport on the Atlantic Continental Shelf of North America: a review. Estuar Coast Shelf Sci. 2001, 52: 51-77. 10.1006/ecss.2000.0727.
Towle DW, Henry RP, Terwilliger NB: Microarray-detected changes in gene expression in gills of green crabs (Carcinus maenas) upon dilution of environmental salinity. Comp Biochem Physiol Part D Genomics Proteomics. 2011, 6 (2): 115-125. 10.1016/j.cbd.2010.11.001.
Borgnia M, Nielsen S, Engel A, Agre P: Cellular and molecular biology of the aquaporin. Ann Rev Biochem. 1999, 68: 425-458. 10.1146/annurev.biochem.68.1.425.
Ishibashi K, Kondo S, Hara S, Mprishita Y: The evolutionary aspects of aquaporin family. Amer J Physiol Regul Interg Comp Physiol. 2011, 300: R566-R576. 10.1152/ajpregu.90464.2008.
Verkman AS, Mitra AK: Structure and function of aquaporin water channel. Am J Physiol Renal Physiol. 2000, 278: F13-F28.
An KW, Kim NN, Shin HS, Kil G-S, Choi CY: Profiles of antioxidant gene expression and physiological changes by thermal and hypoosmotic stresses in black porgy (Acanthopagrus schlegeli). Comp Biochem Physiol A Mol Integr Physiol. 2010, 156: 262-268. 10.1016/j.cbpa.2010.02.013.
Tipsmark CK, Sorensen KJ, Madsen SS: Aquaporin expression dynamics in osmoregulatory tissues of Atlantic salmon during smoltification and sewater acclimation. J Exp Biol. 2010, 213: 368-379. 10.1242/jeb.034785.
Martinez AS, Cutler CP, Wilson GD, Philips C, Hazon N, Cram G: Regulation of expression of two aquaporin homologues in the intestine of the European eel: effects of seawater acclimation and contrisol treatment. Am J Physiol. 2005, 57: 1733-1743.
Aoki M, Kaneko T, Katoh F, Hasegawa S, Tsutsui N, Aida K: Intestinal water absorption through aquaporin 1 expressed in the apical membrane of mucosal epithelial cells in seawater-adapted Japanese eel. J Exp Biol. 2003, 206: 3495-3505. 10.1242/jeb.00579.
Giffard-Mena I, Boulo V, Aujoulat F, Fowden H, Castile R, Charmantier G: Aquaporin molecular characterization inthe sea-bass (Dicentrarchus labrax): the effect of salinity of AQP1 and AQP3 expression. Comp Bioch Physiol. 2007, 148: 430-444. 10.1016/j.cbpa.2007.06.002.
Livingstone M, Atas E, Meller A, Sonenberg N: Mechanisms governing the control of mRNA translation. Phys Biol. 2010, 7: 021001-10.1088/1478-3975/7/2/021001.
Hamilton TL, Stoneley M, Spriggs KA, Bushell M: TOPs and their regulation. Biochem Soc Trans. 2006, 34: 12-16. 10.1042/BST0340012.
Boassa D, Yool AJ: A fascinating tail: cGMP activation of aquaporin-1 ion channels. Trends Pharmacol Sci. 2002, 23: 558-562. 10.1016/S0165-6147(02)02112-0.
Han Z, Patil RV: Protein kinase A-dependent phosphorylation of aquaporin-1. Biochem Biophys Res Commun. 2000, 273: 328-332. 10.1006/bbrc.2000.2944.
An KW, Kim NN, Choi CY: Cloning and expression of aquaporin 1 and arginine vasotocin receptor mRNA from the black porgy, Acanthopagrus schlegeli: effect of freshwater acclimation. Fish Physiol Biochem. 2008, 34: 185-194. 10.1007/s10695-007-9175-0.
Martinez A-S, Cutler CP, Wilson GD, Phillips C, Hazon N, Cramb G: Regulation of expression of two aquaporin homologs in the intestine of the European eel: effects of seawater acclimation and cortisol treatment. Am J Physiol Regul Interg Comp Physiol. 2005, 288: R1733-R1743. 10.1152/ajpregu.00747.2004.
King LS, Agre P: From structure to disease: the evolving tale of aquaporin biology. Nat Rev Mol Cell Biol. 2004, 5: 687-698. 10.1038/nrm1469.
Zardoya R: Phyology and evolution of the major intrinsic protein family. Biol Cell. 2005, 97: 397-414. 10.1042/BC20040134.
Avonce N, Mendoza-Vargas A, Morett E, Iturriaga G: Insights on the evolution of trehalose biosynthesis. BMC Evol Biol. 2006, 19: 119-
Nejsum LN: The renal plumbing system: aquaporin water channels. Cell Mol Life Sci. 2005, 62 (15): 1692-1706. 10.1007/s00018-005-4549-x.
Tingaud-Sequeira A, Calusinska M, Finn RN, Chauvigné F, Lozano J, Cerdà J: The zebrafish genome encodes the largest vertebrate repertoire of functional aquaporins with dual paralogy and substrate specificities similar to mammals. BMC Evol Biol. 2010, 10: 38-10.1186/1471-2148-10-38.
Kaufmann N, Mathai JC, Hill WG, Dow JAT, Zeidel ML, Brodsky JL: Developmental expression and biophysical characterization of a Drosophila melanogaster aquaporin. Am J Physiol Cell Physiol. 2005, 289: C397-C407. 10.1152/ajpcell.00612.2004.
Campbell EM, Ball A, Hoppler S, Bowman AS: Invertebrate aquaporins: a review. J Comp Physiol B. 2008, 178: 935-955. 10.1007/s00360-008-0288-2.
Kuwahara M, Asai T, Sato K, Shinbo I, Terada Y, Marumo F, Sasaki S: Functional characterization of a water channel of the nematode Caenorhabditis elegans. Biochem Biophys Acta. 2000, 1517: 107-112. 10.1016/S0167-4781(00)00268-2.
Huang CG, Lamitina T, Arge P, Strange K: Functional analysis of the aquaporin gene family in Caenorhabditis elegans. Am J Physiol Cell Physiol. 2007, 292: C1867-C1873.
Selection for retention or export of larvae in estuaries. Edited by: Strathmann RR. 1982, New York: Academic Press
Charmantier G, Charmantier-Daures M, Anger K: Ontogeny of osmoregulation in the grapsid crab Armases miersii (Crustacea, Decapoda). Mar Ecol Prog Ser. 1998, 164: 285-292.
Anger K, Charmantier G: Ontogeny of osmoregulation and salinity tolerance in a mangrove crab, Sesarma curacaoense (Decapoda: grapsidae). J Exp Mar Biol Ecol. 2000, 251: 265-274. 10.1016/S0022-0981(00)00223-9.
Charmantier G, Anger K: Ontogeny of osmoregulation in the palaemonid shrimp Palaemonetes argentinus (Crustacea: Decapoda). Mar Ecol Prog Ser. 1999, 181: 125-129.
Susanto GN, Charmantier G: Ontogeny of osmoregulation in the crayfish Astacus leptodactylus. Physiol Biochem Zool. 2000, 73: 169-176. 10.1086/316736.
Ogburn MB, Jackson JL, Forward RB: Comparision of low salinity tolerance in Callinectes sapidus Rathbun and Callinectes similis Williams postlarvae upon entry into an estuary. J Exp Mar Biol Ecol. 2007, 352: 343-350. 10.1016/j.jembe.2007.08.008.
Chung JS, Zmora N, Tsutsui N, Katayama H: Crustacean hyperglycemic hormone (CHH) neuropeptides family: function, titer, and binding to target tissues. Gen Comp Endocrinol. 2010, 166: 447-454. 10.1016/j.ygcen.2009.12.011.
Spanings-Pierrot C, Soyez D, Van Herp F, Gompel M, Skaret G, Grousset E, Charmantier G: Involvement of crustacean hyperglycemic neurohormone in the control of gill ion transport in the crab Pachygrapsus marmoratus. General and Comparative Endocrinolgy. 2000, 119: 340-350. 10.1006/gcen.2000.7527.
Serrano L, Blanvilain G, Soyez D, Charmantier G, Grousset E, Aujoulat F, Spanings-Pierrot C: Putative involvement of crustacean hyperglycemic hormone isoforms in the neuroendocrine mediation of osmoregulation in the crayfish Astacus leptodactylus. J Exp Biol. 2003, 2003: 979-988.
Chung JS, Dircksen H, Webster SG: A remarkable, precisely timed release of hyperglycemic hormone from endocrine cells in the gut is associated with ecdysis in the crab Carcinus maenas. Proc Nat Acad Sci USA. 1999, 96: 13013-13107.
Webster SG, Keller R, Dircksen H: The CHH-superfamily of multifunctional peptide hormones controlling crustacean metabolism, osmoregulation, moulting and reproduction. Gen Comp Endocrinol. 2012, 175: 217-233. 10.1016/j.ygcen.2011.11.035.
Chung JS, Webster SG: Expression and release patterns of neuropeptides during embryonic development and hatching in the green crab, Carcinus maenas. Development. 2004, 131: 4751-4761. 10.1242/dev.01312.
Zmora O, Findiesen A, Stubblefield J, Fraenkel V, Zohar Y: Large-scale juvenile production of the blue crab Callinectes sapidus. Aquaculture. 2005, 244: 129-139. 10.1016/j.aquaculture.2004.11.012.
Kennedy VS: External anatomy of blue crab larvae. The Blue Crab. Edited by: Kennedy VS, Cronin LE. 2007, College Park: Maryland Sea Grant, 23-54.
Chung JS, Bembe S, Tamone S, Andrews E, Thomas H: Molecular cloning of the crustacean hyperglycemic hormone (CHH) precursor from the X-organ and the identification of the neuropeptide from sinus gland of the Alaskan Tanner crab, Chionoecetes bairdi. Gen Comp Endocrinol. 2009, 162: 129-133. 10.1016/j.ygcen.2009.03.012.
Chung JS, Wilcockson DC, Zmora N, Zohar Y, Dircksen H, Webster SG: Identification and developmental expression of mRNA encoding crustacean cardioactive peptide (CCAP) in decapod crustaceans. J Exp Biol. 2006, 209: 3862-3872. 10.1242/jeb.02425.
Chung JS, Zmora N: Functional studies of crustacean hyperglycemic hormone (CHHs) of the blue crab, Callinectes sapidus- the expression and release of CHH in eyestalk and pericardial organ in response to environmental stress. FEBS J. 2008, 275: 693-704. 10.1111/j.1742-4658.2007.06231.x.
Chung JS, Bachvaroff TR, Trant J, Place A: A second copper zinc superoxide dismutase (CuZnSOD) in the blue crab Callinectes sapidus: Cloning and up-regulated expression in the hemocytes after immune challenge. Fish Shellfish Immunol. 2012, 32 (1): 16-25. 10.1016/j.fsi.2011.08.023.
Huang HY, Chien CH, Jen KH, Huang HD: RegRNA: A regulatory RNA motifs and elements finder. Nucleic Acids Res. 2006, 34: W429-W434. 10.1093/nar/gkl333.
Kelley LA, Sternberg MJE: Protein structure prediction on the web: A case study using the Phyre server. Nat Protoc. 2009, 4: 363-371.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, et al: Phyologeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36: W465-W459. 10.1093/nar/gkn180.
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, et al: CDD: a Conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: 225-229. 10.1093/nar/gkq769.
Chung JS, Manor R, Sagi A: Cloning of an insulin-like androgenic gland factor (IAG) from the blue crab, Callinectes sapidus: Implications for eyestalk regulation of IAG expression. Gen Comp Endocrinol. 2011, 173: 4-10. 10.1016/j.ygcen.2011.04.017.
Chung JS, Webster SG: Moult cycle-related changes in biological activity of moult-inhibiting hormone (MIH) and crustacean hyperglycaemic hormone (CHH) in the crab, Carcinus maenas. Eur J Biochem. 2003, 270: 3280-3288. 10.1046/j.1432-1033.2003.03720.x.
Chung JS: A trehalose 6-phosphate synthase gene of the hemocytes of the blue crab, Callinectes sapidus: cloning, the expression, its enzyme activity and relationship to hemolymph trehalose levels. Saline Systems. 2008, 4: 18-10.1186/1746-1448-4-18.
Authors thank the personnel in the blue crab hatchery program for the larval rearing. This article is contribution No. 4651 of the University of Maryland Center for Environmental Science and contribution No. 12–232 of the Institute of Marine and Environmental Science. This project was supported by the NOAA Living Marine Resources Cooperative Science Center (NOAA Award No. NA11SEC4810002) and NSF-UMES CREST center grant (1036586). NIH MBRS Award R25GM063775 supported MB’s summer internship at Chung’s laboratory at the Institute of Marine and Environmental Technology (Baltimore, MD).
Authors declare no conflicting interests.
JSC for cloning AQP-1, qRT-PCR assays and data analyses; JSC, JSP, and MBO for larval experimental design; MB for salinity exposure study; LM for salinity exposure study, RNA extraction and cDNA synthesis; JSC & JSP for funding; JSC, JSP & MBO for writing the manuscript. All authors read and approved the final manuscript.