Skip to main content

Phylogenomic and functional analyses of salmon lice aquaporins uncover the molecular diversity of the superfamily in Arthropoda

Abstract

Background

An emerging field in biomedical research is focusing on the roles of aquaporin water channels in parasites that cause debilitating or lethal diseases to their vertebrate hosts. The primary vectorial agents are hematophagous arthropods, including mosquitoes, flies, ticks and lice, however very little is known concerning the functional diversity of aquaporins in non-insect members of the Arthropoda. Here we conducted phylogenomic and functional analyses of aquaporins in the salmon louse, a marine ectoparasitic copepod that feeds on the skin and body fluids of salmonids, and used the primary structures of the isolated channels to uncover the genomic repertoires in Arthropoda.

Results

Genomic screening identified 7 aquaporin paralogs in the louse in contrast to 42 in its host the Atlantic salmon. Phylogenetic inference of the louse nucleotides and proteins in relation to orthologs identified in Chelicerata, Myriapoda, Crustacea and Hexapoda revealed that the arthropod aquaporin superfamily can be classified into three major grades (1) classical aquaporins including Big brain (Bib) and Prip-like (PripL) channels (2) aquaglyceroporins (Glp) and (3) unorthodox aquaporins (Aqp12-like). In Hexapoda, two additional subfamilies exist as Drip and a recently classified entomoglyceroporin (Eglp) group. Cloning and remapping the louse cDNAs to the genomic DNA revealed that they are encoded by 1–7 exons, with two of the Glps being expressed as N-terminal splice variants (Glp1_v1, −1_v2, −3_v1, −3_v2). Heterologous expression of the cRNAs in amphibian oocytes demonstrated that PripL transports water and urea, while Bib does not. Glp1_v1, −2, −3_v1 and −3_v2 each transport water, glycerol and urea, while Glp1_v2 and the Aqp12-like channels were retained intracellularly. Transcript abundance analyses revealed expression of each louse paralog at all developmental stages, except for glp1_v1, which is specific to preadult and adult males.

Conclusions

Our data suggest that the aquaporin repertoires of extant arthropods have expanded independently in the different lineages, but can be phylogenetically classified into three major grades as opposed to four present in deuterostome animals. While the aquaporin repertoire of Atlantic salmon represents a 6-fold redundancy compared to the louse, the functional assays reveal that the permeation properties of the different crustacean grades of aquaporin are largely conserved to the vertebrate counterparts.

Background

Aquaporins are protein channels, which are assembled as tetramers in biological membranes to facilitate the rapid, yet selective conductance of water and other small solutes and gasses down their concentration gradients [1]. Each protomer is comprised of a conserved tertiary structure consisting of 6 transmembrane domains with two hemihelices that typically project opposing Asn-Pro-Ala (NPA) motifs as part of a central proton-excluding selectivity filter [2, 3]. In some channels, the NPA motifs may be modified for other functions [4, 5]. An additional selectivity filter located in the outer channel vestibule typically consists of aromatic amino acids and an Arg residue (ar/R), which determines the molecular substrate specificity of the pore [6, 7]. Previous research has shown that aquaporins are encoded in the genomes of nearly all organisms studied to date, with diverse subfamilies expressed in prokaryotic and eukaryotic organisms [8, 9]. Although high gene copy numbers are found in plants and deuterostome animals, the superfamily has been categorised into four major grades consisting of classical water-selective type aquaporins (Aqp0, −1, −2, −4, −5, -5-like, −6, −14, and −15), Aqp8-type aquaamoniaporins (Aqp8 and −16), aquaglyceroporins (Glp, Aqp3, −7, 9, −10 and −13) and unorthodox aquaporins (Aqp11 and −12) [913].

Recent research in the biomedical field has begun to focus on the roles of aquaporin water channels in drug-resistant parasites as potential therapeutic targets [14, 15]. This includes aquaporins expressed in protozoan parasites causing malaria, toxoplasmosis, trypanosomiasis and leishmaniasis, and metazoan parasites causing fascioliasis and schistosomiasis in mammals [14, 1620]. Experiments have shown that double allelic disruption of the major glycerol-transporting channel (AQP9) in mouse erythrocytes inhibits the virulence of the malaria causing Plasmodium species of parasite by attenuating its growth and lipogenesis [17, 18]. Other studies have investigated aquaporins in the hexapod mosquitoes (Diptera) and arachnid ticks (Acari) that act as vectors for such parasites or other bacterial and viral pathogens [2129]. Unlike the protozoan parasites, which infect endogenous tissues, however, mosquitoes and ticks are ectoparasites, which gorge on the blood of mammals as part of their reproductive cycle [30, 31]. This hematophagous behaviour induces osmoregulatory stress due to the manifold quantities of fluid imbibed [24, 25, 32, 33], and a number of aquaporins have been identified that are considered to be involved in the fluid homeostasis of ingestion, excretion and chemoreception. For dipteran mosquitoes, chironomids and tsetse flies, this includes up to six aquaporins variously labelled as AeaAQP or AQP1 to −6 [21, 22, 2527, 29, 32, 34], which may be confused with the mammalian aquaporins (AQP1 to −6) encompassing both water-selective channels and Glps. For other insects, an alternative nomenclature has been partially applied [35], based upon their phylogeny and function, including water-selective channels similar to the Drosophila intrinsic protein (Drip) [3640], a water and urea-transporting channel similar to Pyrocoelia rufa integral protein (Prip) [41, 42], and a cell-adhesive, cation-permeating neurogenic protein termed big brain (Bib) [43, 44]. The aquaporin nomenclature for ticks is also confusingly annotated as AQP1 or −2 [23, 24, 28], which are Glps rather than the water-selective channels found in mammalian genomes. As for the mouse AQP9 trial, it is one of the tick Glps that has been targeted for vaccine development. In this instance, recombinant regions of the Glp (RmAQP1) were intramuscularly injected into cattle, and provided sufficient antigenic capacity to reduce tick infection by up to 75 % [28].

An equally problematic blood-feeding ectoparasite with broad resistance to available drugs is the saltwater salmon louse (Lepeophtheirus salmonis), which is a caligid crustacean copepod (class: Maxillipoda) causing severe infections of salmonids across the Atlantic and Pacific regions [45, 46]. It has been established that in contrast to its salmonid host, which is an osmoregulator, the ionic concentration of the louse’s hemolymph conforms to full strength seawater and remains significantly hyperosmotic in dilute seawater, however, cell volume regulation breaks down in fresh and brackish waters [47]. These features suggest that aquaporins likely play important physiological roles in the louse’s fluid homeostasis, which, given the divergent osmoregulatory capabilities to its salmonid host, could represent eminent therapeutic targets. To date, however, almost nothing is known of the genomic diversity of these channels in crustaceans, which are represented by tens of thousands of species [48]. To redress this lack of knowledge we screened the available genome of the salmon louse, and identified the full complement of water channel genes and cloned the corresponding cDNAs. We then used maximum likelihood and Bayesian inference of the louse nucleotides and deduced proteins to classify the genomic repertoires of aquaporins in the four extant lineages of Arthropoda, the Chelicerata (scorpions, spiders, ticks), Myriapoda (centipedes), Crustacea (water fleas, copepods) and Hexapoda (insects). We subsequently functionally characterised the permeability properties and stage specific-expression patterns of each channel, and used the integrated results to propose a uniform aquaporin nomenclature for arthropods and vertebrates.

Results and discussion

The aquaporin superfamily of the salmon louse

The screening of the salmon louse genome identified nine putative sequences with homologies to the aquaporin family of proteins. Each putative mRNA sequence was isolated and cloned using gene-specific primers and rapid amplification of cDNA ends (RACE) analyses. This resulted in 9 cDNAs with open reading frames ranging in length from 777–1992 bp, encoding proteins of 258–663 amino acids, which, following phylogenetic analyses (see below), we termed Bib, Prip-like (PripL), Glp1-3, Aqp12-like1 (Aqp12L1) and Aqp12-like 2 (Aqp12L2) (Fig. 1a). We found no evidence for the existence of a putative lsal-evj-541-119 Aqp (BT120567) deposited in GenBank, either experimentally using PCR with gene specific primers or in silico via searches of the salmon louse genome.

Fig. 1
figure1

Molecular structures of salmon lice aquaporins. a Multiple sequence alignment of the full-length proteins showing secondary structural features including the N-terminal regions (NT), α-helices (H), transmembrane domains (TMD), loops A-E and C-termini (CT). Unique (white text on black background) and common (boxed residues) putative phosphorylation sites are shown for the N- and C-teminal regions. The two Asn-Pro-Ala (NPA) motifs are highlighted in green and yellow on hemihelices 3 and 7, respectively. Aromatic-arginine (ar/R) residues associated with the selectivity filter in the outer channel vestibule are highlighted in magenta, organge, red and blue, respectively. The position of Ser172, which normally corresponds to the ar/R residue on H2 is highlighted in gray for Bib. Residues predicted to confer glycerol selectivity (P1–P5) are highlighted in gray [50]. b, c and d Extracelluar views of cartoon renders of Bib, PripL and Glp1_v1 with ar/R residues (spacefill) coloured according to highlighted residues in (A). Helices and extracellular loops are annotated

Alignment and in silico analyses of the secondary and tertiary structures showed that each louse protein consists of six transmembrane domains, two Asn-Pro-Ala (NPA)-like motifs at the termini of two hemihelices, and five loops A – E typically found in mammalian aquaporins [49]. The largest protein, which has a deduced length of 633 amino acids and a predicted molecular mass of 72 kD is Bib, and harbours a long N-terminus (135 amino acids) and an exceptionally long C-terminus (329 amino acids), similar to the subdomain structure of the Drosphila Bib molecule [43]. The N-terminus of the louse Bib is thus equivalent in length to the membrane-spanning domains (169 amino acids), while the C-terminus is twice the size. Three dimensional reconstruction of the louse Bib revealed that in contrast to the majority of aquaporins displaying an ar/R residue on transmembrane domain 2 (TMD2), the putative constriction residue is Phe146 on TMD1, rather than Ser172 on TMD2 (Fig. 1b). Conversely, PripL, which represents the smallest aquaporin in the louse superfamily with 258 amino acids and a predicted molecular mass of 27.8 kD, displays the more canonical ar/R constriction and functional (P1–P5) residues associated with the water-selective branch of aquaporins (Fig. 1a, c) [7, 50]. The tertiary structures of the Glps, which range in length from 266–332 amino acids with predicted molecular masses from 28.6–35.8 kD, appear less constricted in the central pore due to a generally more relaxed fold conformation and an uncharged Thr on TMD5 of Glp1_v1 and −1_v2 (Fig. 1d, Additional file 1: Figure S1A). Interestingly, although Glp2, −3_v1 and −3_v2 harbour a His on TMD5 that is typical of water-selective channels [7], the in silico models indicate that the TMD5 His does not strongly protrude into the central channel as in PripL and thus may not present a hindrance to glycerol conductance (Additional file 1: Figure S1B-D). The remaining sequences, Aqp12L1 and -12L2, were more divergent (11–18 % identity over the membrane-spanning domains compared to the other louse aquaporins) with lengths of 262–263 amino acids and predicted molecular masses of 29.1 and 28.8 kD, respectively. These latter channels have an ungapped alignment identity of 61.8 % and encode non-canonical CPY and NPV motifs associated with unorthodox members of the aquaporin superfamily [4]. Three-dimensional models based upon the structure mask of AqpZ (2ABM, Z-score = 33–34) indicate that Aqp12L1 is highly constricted due to a Phe on TMD2, while Aqp12L2 has a more open conformation due to the substitution of the Phe for an Ile in this position (Additional file 1: Figure S1E-F).

Molecular phylogeny of salmon louse aquaporins uncovers the genomic repertoires in Arthropoda

To understand the evolutionary history of the salmon louse aquaporins, each paralog was subjected to phylogenetic analyses in relation to 435 orthologs retrieved from Chelicerata, Myriapoda, Crustacea and Hexapoda using maximum likelihood and Bayesian protocols. In addition, 21 prokaryotic AqpZ, GlpF and AqpM orthologs were included to ensure that none of the isolated louse aquaporins were derived from prokaryotic endosymbionts. Two main arthropod data sets were constructed consisting of 350 non-redundant orthodox aquaporins (Fig. 2a) and 94 non-redundant unorthodox aquaporins (Fig. 2b). These analyses revealed that two of the louse paralogs, Bib and PripL are members of the classical water-selective branch of aquaporins as judged by their derivative clustering above bacterial AqpZ. Although the salmon louse Bib clustered as an arthropod Bib, it did not cluster together with the other malacostracan (crab and crayfish) or branchiopod (water flea) crustaceans. This is primarily due to the derived nature of the louse sequence, which has a membrane-spanning identity of 27–28 % compared to the 46 other arthropod Bibs analysed. In contrast the water flea (Daphnia pulex) Bib is less derived with 32–38 % identity to the membrane-spanning regions of other arthropod Bib-like proteins, and clusters as a sister branch to Hexapoda consistent with the pancrustacean model of arthropod evolution [51, 52]. The present phylogenetic data thus reveal that neurogenic Bib channels are not only found in dipteran or hemipteran insects [25, 32, 43, 5355], but are encoded in the genomes of all extant lineages of Arthropoda.

Fig. 2
figure2

Phylogeny of the genomic repertoires of arthropod aquaporins compared to the tetraploid complement in Atlantic salmon. a Maximum likelihood of arthropod orthodox aquaporins. The tree is rooted with AqpM and was inferred from 30 million MCMC generations on 484,924 nucleotide sites from a codon alignment of 350 non-redundant aquaporins. Posterior probabilities of the codon/amino acid analyses are shown at each node, with “-“ representing a polytomy (<50 %). The numbers of taxa are indicated in the collapsed branch triangles. Scale bars indicate the rate of nucleotide substitution per site. b Bayesian majority-rule phylogeny of arthropod unorthodox aquaporins. The tree is mid-point rooted and was inferred from 5 million MCMC generations on 87,613 nucleotide sites in a codon alignment of 94 non-redundant aquaporins. Annotations as in (a). c Genomic organisation of the salmon louse aquaporins. Colored and gray boxes indicate coding and non-coding exons, respectively. Long introns are annotated in kb, and gray backgrounds linking paralogs highlight conserved exons. Exon numbers are annotated with roman numerals, with the first proximal and distal exons of splice variants labelled A and B, respectively. d Gene copy numbers and orthology of salmon louse and Atlantic salmon aquaporins

The second water-selective type classical aquaporin isolated from the louse clustered together with other maxillipod and crustacean water channels that are closely related to insect Prips. Previously Prip-like channels, variously named as Prip, AQP, AQP1, AQP2, AQP4A, AQP4B have only been identified in Hexapoda [2527, 32, 34, 35, 41, 42, 5558], however the present findings reveal that Prip-like water channels exist throughout the Arthropoda. For most species, including the salmon louse, single copy genes were identified, although spiders (Araneae) have at least two, and the centipede species studied has specifically expanded this class of aquaporins with 11 divergent paralogs (27–28 % identity). Very recently we have also identified multiple prip-like genes in the oldest hexapod lineages, including Archeognatha, Diplura, Collembola and Protura [59]. In addition to the prip genes, the phylogenetic analyses show that two insect-specific classes of aquaporin exist as drip and the recently identified glycerol transporters termed entomoglyceroporins (eglps [59]). We therefore conclude that arthropods have two ubiquitous classes of classical aquaporin proteins termed Bib and PripL, while while hexapods have four, Bib, Prip, Drip and the glycerol-transporting Eglp group.

The remaining orthodox aquaporins isolated from the salmon louse clustered as canonical glps, as opposed to orthologs in bacterial symbionts (Fig. 2a). Interestingly the phylogenetic signal revealed that glp evolution in arthropods appears to be polyphyletic, with multiple gene copies existing in the different lineages. For the louse, two subclusters are observed, with the first including glp1 co-clustering with other crustacean glps, while the second, which includes glp2 and −3, clusters with chelicerate glps. These data suggest that glp1-type genes share a common ancestor in Crustacea, while glp2- and −3-type genes have a different origin. In contrast to the polyphyletic evolution of maxillipod and chelicerate glps, however, the data for Branchiopoda indicate that the glp repertoires specifically expanded within this lineage with 6 gene products on the same branch.

The present data reveal that aqp12-related genes are encoded in the genomes of each of the major lineages of arthropod (Fig. 2b). This latter finding is therefore consistent with previous studies that identified Aqp12-like orthologs in Diptera [32, 60], Tardigrada (water bears) and Cnidaria (corals and sea anemones), the latter of which represent basal members of the Panarthropoda and Metazoa, respectively [9, 61]. Unlike the other arthropod species studied here, however, we show that the salmon louse encodes two unorthodox aquaporins (Aqp12L1 and -12L2), which co-cluster with other maxillipod Aqp12-like sequences as a sister branch to hexapods. Consequently, despite the relatively low amino acid conservation between the Aqp12L1 and -12L2 paralogs (61.8 %), the phylogenetic data imply that the genes may have diverged within parasitic caligid copepods.

Genomic structure and N-terminal splice variants

To decipher the genomic organization of the salmon louse aquaporins, we mapped each mRNA sequence, including the 5′ and 3′untranslated regions (UTR), to the genomic DNA. This revealed that Although the bib transcript includes a second exon for the 3′UTR, the Bib protein is encoded within a single exon, while the coding regions of all other paralogs are located within 3–4 exons (Fig. 2c). The single coding exon structure of the Bib protein is unusual compared to other arthropod orthologs, which are encoded by five exons in Drosophila, up to nine exons in the assassin bug (Rhodnius prolixus), and based upon the tblastn searches of whole genome shotgun sequences at least three exons in the common house spider (Parasteatoda tepidariorum) and at least five exons in the coastal European centipede (Strigamia maritima). The intronless structure of the coding region of the louse bib gene thus appears to be specific to the salmon louse, and may have evolved via reverse transcription and replacement of an older split gene as suggested for human G-protein-coupled receptors [62].

The pripL gene is by far the longest spanning >43 kb with four coding exons and three non-coding exons in the 5′ UTR. Unlike the prip of the malaria mosquito (Anopheles gambiae), which differentially expresses two splice variants in the ovary and gut [27], no splice variants were detected for the louse pripL. By contrast two of the louse glps (glp1 and −3), which are encoded by 2–3 exons are transcribed as splice variants. In both cases the variation occurs at the N-terminus whereby either a proximal exon (v1) or a distal exon (v2) is spliced. We therefore named these isoforms glp1_v1, glp1_v2 and glp3_v1, glp3_v2, respectively (Fig. 2c). Interestingly, a glp1_v1-like transcript described here for the salmon louse has recently been identified through high-throughput transcriptome sequencing as a sex-specific variant termed AQP3 in a different parasitic copepod, Caligus rogercresseyi [63]. To determine whether related splice variants exist in other copepods, we analyzed the available transcriptomes of five species, including the salmon louse, and found that N-terminal splice variation of both glp1 and −3-type genes appears to be a conserved feature within the maxillopod class of Copepoda (see Additional file 1: Figure S2).

Genomic analysis of the glp2 gene revealed that it is closely related to glp3 with 77 % nucleotide identity for the overlapping regions of exons II and III. Since this gene has not previously been reported for any organism, and glp3-type genes are found in the transcriptome shotgun assemblies of other species of copepod, it seems likely that glp2 may have recently evolved as a duplicate of glp3 in the salmon louse. By contrast both of the louse aqp12-like genes are expressed in the transcriptomes of other caligid parasitic copepods, such as Caligus rogercresseyi, but apparently not in other members of the Crustacea. The genomic organisation of the louse aqp12L1 and -12L2 paralogs is conserved with 4 exons with an extended first intron in the aqp12L1 gene (Fig. 2c). Although the Drosphila aqp12-like gene is also split into 4 exons, this is not a conserved feature of unorthodox aquaporins in Arthropoda as the number of exons varies between species with up to 8 exons in butterflies (Lepidoptera).

Aquaporin orthology between parasite and host

To identify potential therapeutic targets, we compared the seven paralogs isolated from the louse to the 42 orthologs recently identified in its tetraploid host, the Atlantic salmon [9, 6467], we aligned the two superfamilies and examined their phylogenetic distributions using Bayesian inference. These results provided robust statistical evidence that the louse bib and pripL gene products are closely related to the salmon aqp4 paralogs, and thus members of the classical grade of aquaporins (Additional file 1: Figure S3). Similarly the louse glp and aqp12-like gene products are related to the salmon aquaglyceroporins and unorthodox aquaporins, respectively, revealing that the salmon louse superfamily consists of 3 major grades of aquaporins, while that of the Atlantic salmon consists of four. Consequently due to serial rounds of whole genome duplication in tetraploid Salmonidae [68, 69], there is a 6-fold redundancy of aquaporins encoded in the genome of Atlantic salmon compared to the salmon louse. A summary of these findings is given in Fig. 2d which illustrates the absence of vertebrate aqp8-type aquaammoniaporins in the louse. Based upon the broader molecular phylogenetic analysis (Fig. 2a), the present data further indicate that this is also true for each of the major lineages of arthropod.

Channel permeation properties of the salmon louse aquaporins

To corroborate the phylogenetic and structural analyses, the functional permeation properties of the salmon louse aquaporins were examined by heterologous expression in Xenopus laevis oocytes followed by measurement of the oocyte osmotic water permeability (P f) in hypotonic media or uptake of radioactive solutes (glycerol and urea). The subcellular localization of the aquaporins in the injected oocytes was monitored by immunofluorescence microscopy and immunoblotting using specific affinity-purified antibodies raised in rabbits against synthetic C-terminal peptides (see Additional file 1: Table S1). Injection of 5–10 ng of cRNA and subsequent expression of PripL, Glp3_v1 and Glp3_v2 targeted the proteins to the oocyte plasma membrane (Fig. 3) and induced a 2.5–4.5-fold increase in the oocyte water permeability compared to water injected controls (Fig. 4). Although antibody development for the Glp2 paralog was not successful, the heterologous expression experiments showed that it functions as a water transporter. By contrast, oocytes expressing high cRNA titres (20 ng) of the Glp1_v1 splice variant were not permeable to water due to the retention of the channel just below the oocyte plasma membrane (see Additional file 1: Figure S4). Addition of dibutyryl cyclic adenosine monophosphate (cAMP) in the bathing media was necessary to traffic the channel to the plasma membrane (Fig. 3) and induce the increase in water permeability (Fig. 4). This suggested that a cAMP-dependent phosphorylation mechanism may be involved in the trafficking of the Glp_v1 channel to the plasma membrane as reported for some vertebrate [10, 12, 7074] and plant aquaporins [7577]. An in silico analysis of the putative phosphorylation sites indicates that five residues in the N-terminus of Glp1_v1 (Ser2, Thr3, Tyr9, Thr14and Ser21) differ in phosphorylation potential compared to the Glp1_v2 isoform (Fig. 1a), implying that a casein type kinase, a protein kinase type C or a receptor tyrosine kinase [78] may be involved in the trafficking of this splice variant to the oocyte plasma membrane. By contrast, the 17 amino acid shorter Glp1_v2 isoform, which retains only one unique putative phosphorylation site, was not permeable to water in Xenopus oocytes, even with cAMP added to the media. The immunolocalization and Western blot experiments revealed that Glp1_v2 was retained intracellularly (Fig. 3a) and was quickly degraded in the oocytes (see ~20 kD band in Fig. 3b).

Fig. 3
figure3

Immunolocalisation of salmon louse aquaporins expressed in X. laevis oocytes. a Immunofluorescence photomicrographs of paraffin sections of water- (Ctrl) and salmon louse cRNA-injected oocytes probed with paralog-specific antisera followed by Cy3-labeled anti-rabbit IgG. Arrows point to the plasma membrane. Sections shown are representative of five different sections per treatment showing identical staining. b Representative immunoblots of total oocyte membranes (0.5-3 oocyte equivalents per lane) with molecular markers shown to the left. Arrows point to bands that match the predicted molecular mass of each aquaporin paralog, while lower bands likely repesent misfolded and degraded proteins. For Aqp12L2 a cross-reaction is seen against a ~45 kD X. laevis protein

Fig. 4
figure4

Functional characterization of salmon louse aquaporins. Osmotic water permeability in the presence and absence of HgCl2 and βME (P f, left panel), and urea and glycerol uptake (right panel) of X. laevis oocytes expressing the salmon louse aquaporins. Values (mean ± SEM; n = 10–15 oocytes per treatment) are from representative experiments. Statistical differences (* = p < 0.05, ** = p < 0.01, ***p < 0.001) are indicated with brackets

While cAMP-dependent membrane trafficking involving the phosphorylation of C-terminal aquaporin residues is well established [10, 12, 7077], the role of the N-terminus in this pathway is less common. Nevertheless, the N-terminus has been associated with trafficking of some mammalian channels, including AQP4 and AQP6 [79, 80]. In the present context our data suggest that only the N-terminus of the Glp1 isoforms retains putative phosphorylation sites which may be the targets of kinases in the cAMP-dependent trafficking mechanism. Further research will be necessary to decipher the precise mechanism involved.

Addition of 0.3 mM HgCl2 abolished the water transport function of the PripL, Glp1_v1 and −2 channels and significantly (p < 0.001) inhibited the water permeability of oocytes expressing Glp3_v1 and −3_v2 (Fig. 4). This inhibition was reversed by the addition of β-mercaptoethanol (βME) to Glp1_v1, −2 and −3_v1, but not to PripL or Glp3_v2. Studies of mammalian AQP1 channels have established that the mercury-induced inhibition is associated with the binding of Hg2+ to a Cys189 residue upstream of the second NPA motif resulting in a collapse in the orientation of several residues at the ar/R region [8183]. Consequently, the irreversible mercury sensitivity of the louse PripL channel may therefore be associated with three Cys residues (Cys77, Cys184 and Cys185) located upstream upstream of the two NPA motifs (see Fig. 1a). However, as for certain piscine [11], mammalian [84, 85] and plant [8688] aquaporins, the mercury sensitivity of the other louse paralogs was independent of any Cys residue upstream of either NPA motif. This observation reaffirms that the mechanism of mercurial inhibition is complex and likely paralog-specific due to the position of the Cys residues in the primary structures [11, 85, 88].

Heterologous expression of the Bib, Aqp12L1 and -12L2 channels in X. laevis oocytes did not induce any significant swelling compared to the water injected controls (Fig. 4). The immunolocalization and Western blot experiments revealed that despite being translated and expressed as 27–29 kD proteins, both of the unorthodox aquaporins were retained intracellularly (Fig. 3a, b). Interestingly, although a cross reaction was noted for the Aqp12L1 antibody in the Western blot, it was not observed in the control immunohistochemical analysis, indicating that it is the Aqp12L2 protein located in the yolk vesicles of the cRNA-injected oocytes. In contrast to these observations, the affinity-purified antibody developed against a C-terminal peptide of Bib revealed that it is expressed as a 72 kD protein and trafficked to the oocyte plasma membrane (Fig. 3a, b). Consequently the absence of water transport through the louse Bib paralog is potentially due to the unusual constriction residue (Phe146) on TMD1 (Fig. 1b), and may support its role as an ion channel and/or a cell adhesion molecule rather than a water transporter as reported for the Drosophila ortholog [44, 89]. Further research will be necessary to identify the molecular and cell physiological function of the louse Bib.

The uptake assays for urea and glycerol revealed that PripL transports a significant amount of urea compared to the controls, but not glycerol, while all of the Glps expressed in the plasma membrane transport both urea and glycerol (Fig. 4). These findings are consistent with the open pore configurations of the Glps (Fig. 1d; Additional file 1: Figure S1), and a recent report demonstrating that a Prip ortholog isolated from the German cockroach (Blattella germanica) transports water and urea, but not glycerol [42].

Expression analysis of salmon louse aquaporins

To gain insight into the developmental importance of the louse aquaporins, the stage-specific expression patterns were evaluated in whole individuals in relation to elongation factor 1α (ef1α) by RT-PCR. With the exception of glp1_v1, each paralog was detected in all developmental stages from egg to adult (Fig. 5). The most striking expression pattern was observed for the glp1_v1 isoform, which is expressed in pre-adult II and adult males, while the glp1_v2 isoform is expressed in the pre-adult and adult stages of both sexes. This finding thus supports the recent observation of the differential expression of an aqp3-annotated gene in the parasitic copepod C. rogercresseyi [63]. It is not clear, however, whether C. rogercresseyi also expresses a shorter splice variant of this gene. To date, very few arthropod aquaporins have been identified as splice variants, and sex-specific splicing has until now only been reported in the malaria mosquito Anopheles gambiae [27]. The present study reveals that the functional repertoire of crustacean aquaporins is increased through splice-variation, which can contribute to the sex-specific fluid homeostasis of crustaceans. Further studies will be necessary to understand how the aquaporin genes are regulated and what are the subcellular sites of expression, in order to better understand the differential roles of the parasite channels compared those in the host.

Fig. 5
figure5

Developmental expression of salmon louse aquaporins. RT-PCR data are shown against expression of elongation factor 1α for egg, nauplius (N), the infective copepodite (Cop) stage and the obligate chalimus (Ch) pre-adult (P and adult parasitic stages for males and females. No template controls (NTC) were included to exclude the possibility of non-specific amplification due to primer-dimer formation or contamination

Conclusions

In the present study we identified seven aquaporin genes in the salmon louse compared to 42 in its host, the Atlantic salmon. However, in contrast to the salmon repertoire, which is phylogenetically divided into four major grades of classical-type water channels, Aqp8-type aquaamoniaporins, aquaglyceroporins and unorthodox aquaporins, Bayesian inference revealed that the salmon louse and other members of the Arthropoda have lost the Aqp8-type orthologs. The remaining classes of louse and non-hexapod aquaporins are classified as Bib, Prip-like, Glps and Aqp12-like channels, while hexapod insects have an additional Drip channel and a specific group of Eglps, which supplanted the ancestral Glps as the major vehicles of glycerol conductance. Heterologous expression of the cloned louse aquaporins in X. laevis oocytes showed that the Prip-like channel is a water and urea transporter, while Bib is not, and all but one of the Glps transports water, glycerol and urea. Expression data showed that each aquaporin is expressed throughout the life cycle of the salmon louse, with one exception where an N-terminal splice variant of Glp1_v1 is specific to pre-adult and adult males. Taken together, these data are the first to reveal the substrate permeability properties of aquaporins in a crustacean organism and uncover the broader genomic diversity of the superfamily in the four extant lineages of Arthropoda.

Methods

Genome analysis and cDNA cloning

Salmon louse aquaporin-coding sequences were identified by a basic local alignment search tool (BLAST) in the salmon louse genomic database (sealouse.imr.no). Based on these sequences, we designed gene-specific primers to amplify the 3′- and 5′-terminus by rapid amplification of cDNA ends (RACE) following the manufacturers instructions (SMARTer TM RACE, Clontech). All experimental procedures were done in accordance with national legislation. RNA was extracted with Trizol (Life Technologies) from laboratory-reared strains of salmon lice (copepodites and adult females), stored on RNA-later® (Life Technologies) and was DNAse treated (Ambion) prior to preparation of 5′and 3′ RACE-cDNA. After sequencing the 5′ and 3′ ends of each aquaporin paralog, new primers with restriction enzyme sites were respectively designed in the 5′UTR and 3′UTR regions in order to produce full length clones that included the open reading frames (ORF). cDNA template was prepared from DNAse treated RNA as described above with the Affinityscript cDNA synthesis kit (Agilent Technologies) followed by a ten-fold dillution of the cDNA. The PCR reaction for producing the aquaporin-clones was carried out using a high fidelity polymerase (Phusion, New England Biosystems, Canada). The resulting PCR product was run on a 1 % agarose gel, purified with a gel extraction kit (Quiagen gel extraction kit), digested with EcoRV/SpeI estriction enzymes and ligated into the X. laevis oocyte expression vector pT7Ts – vector [90] using a T4 ligase (Quick ligation kit, New England Biosystems). Finally, the pT7Ts vector containing the ligated aquaporin paralogue was cloned into E. coli (Invitrogen, USA), and sequenced plasmids from overnight cultures from single bacterial clones were used in later experiments. The cloned aquaporin mRNA sequences were further used to identify the intron-exon borders of the respective genes by alignment to the genomic sequences of the salmon louse. The experiments were performed in accordance with Norwegian legislation and approved by the Norwegian Animal Research Authority (permits nr. 2010/ 245410 and 2009/186329).

Phylogenetic and structural analyses

Phylogenetic analyses were performed as described previously [9, 9193] using maximum likelihood (PAUP v4b10-x86-macosx) and Bayesian inference (Mr Bayes v3.2.2) of codon and deduced amino acid alignments of 435 non-redundant arthropod aquaporins. Orthologous sequences were obtained from open sources (Ensembl, Genbank, whole genome shotgun contigs, transcriptome shotgun assemblies, and expressed sequence tag databases) and assembled as described previously [9]. A full list of the accession numbers of sequences used in the study is provided in Additional file 1: Table S2. Multiple sequence alignments were constructed using default t-coffee v9.01 [94] or L-INS-I MAFFT v7.058b [95] algorithms, and Bayesian model parameters were nucmodel = 4by4, nst = 2, rates = gamma for codon alignments, and aamodel = mixed for amino acid alignments. Markov chain Monte Carlo (MCMC) algorithms were run with 3 heated and 1 cold chain with resulting probability distributions examined for convergence using Tracer version 1.6 (tree.bio.ed.ac.uk/software/tracer/), and majority rule consensus trees summarized with a burnin of 25 %. All trees generated were processed with Archaeopteryx [96] and rendered with Geneious (Biomatters Ltd, New Zealand). A nexus alignment and tree file is deposited in the Dryad Digital Repository. Doi:10.5061/dryad.kd742.

The three-dimensional structures of the salmon louse aquaporins were built using the model leverage option in the Modeller server (modbase.compbio.ucsf.edu). The best scoring models were selected using the slow (Seq-Prf, PSI-BLAST) assignment method and rendered with MacPymol (pymol.org).

Functional expression in X. laevis oocytes

Capped RNAs (cRNAs) were synthesized in vitro with T7 RNA Polymerase (Roche) from XbaI-linearized pT7Ts vector containing the different aquaporin sequences. The isolation and microinjection of stage V-VI oocytes was performed as described previously [90]. Oocytes were injected with 50 nl of water alone (control) or containing 1–25 ng cRNA and incubated in Barth’s culture medium (MBS; 0.33 mM Ca(NO3)2, 0.4 mM CaCl2, 88 mM NaCl, 1 mM KCl, 2.4 mM NaHCO3, 10 mM Hepes, 0.82 mM MgSO4, pH 7.5) at 18 °C over night. The following day the oocytes were de-folliculated and the third day the oocytes were tested for osmotic permeability in a swelling assay.

Swelling assays

The P f was measured in control and injected oocytes from the time course of oocyte swelling in a standard assay. Oocytes were transferred from 200 mOsm MBS to 20 mOsm MBS at room temperature. Oocyte swelling was followed by video microscopy using serial images at 2 s intervals during the first 20 s period. The P f values were calculated taking into account the time-course changes in relative oocyte volume [d(V/Vo)/ dt], the molar volume of water (Vw = 18 cm3/ml) and the oocyte surface area (S) using the formula Vo [d(V/ Vo )/dt]/[SVw (Osmin - Osmout )]. To examine the inhibitory effect of mercury on P f oocytes were pre-incubated for 15 min in MBS containing 0.3 mM HgCl2 before and during the swelling assays. To determine the reversibility of the inhibition, the oocytes were rinsed 3 times with fresh MBS and incubated for another 15 min with 5 mM β-mercaptoethanol before being subjected to swelling assays. In some experiments, the cell-permeable dibutyryl cAMP analog (N6,2′-O-Dibutyryladenosine 3′,5′-cyclic monophosphate sodium salt; Sigma-Aldrich) was added to the MBS for 30 min prior to and during the hypotonic shock.

Radioactive solute uptake assays

To determine the uptake of [3H]-glycerol (60 Ci/mmol) and [14C]-urea (52 mCi/mmol) groups of 6 oocytes, control oocytes or oocytes injected with cRNA, were incubated in 200 μl of MBS containing 20 μCi of the radiolabeled solute (cold solute was added to give 1 mM final concentration) at room temperature. After 10 min (including zero time for subtraction of the signal from externally bound solute), oocytes were washed rapidly in ice-cold MBS three times, and individual oocytes were dissolved in 5 % SDS for scintillation counting.

Statistical analyses of P f and solute uptake

The P f and solute uptake data are expressed as mean ± SEM. Data shown are from a representative experiment out of 2–3 different trials (n = 10–15 oocytes per treatment). The measured values were analysed by one-way analysis of variance (ANOVA), after log-transformation of the data when necessary, followed by Tukey‘s pairwise comparison with a 95 % confidence interval (p < 0.05).

Antiobody production

N-terminal or C-terminal peptide sequences of each aquaporin-paralog (see Additional file 1: Table S1) were injected in rabbits to raise paralog specific polyclonal antibodies (Agrisera AB, Sweden). The antisera were affinity-purified against the synthetic peptides as described previously [92], and their specificity confirmed by ELISA, as well as by immunofluorescence microscopy and immunoblotting on X. laevis oocytes (see below).

Immunohistochemistry and immunoblot analysis

For immunohistochemistry, oocytes were fixed for 6 h in 4 % formaldehyde (4 % PFA dissolved in PBS), washed in PBS before being dehydrated in steps by ascending concentrations of ethanol (50, 70, 96 100 % ethanol in H2O) and finally in xylene, before being embedded in paraffin. The paraffin embedded oocytes were cut in 7 μm thick sections and the slides de-paraffinised in xylene and rehydrated in a series of descending ethanol baths. Epitope retrieval was done by either a 10 min incubation in 0.1 % Triton X-100, or a 10 min incubation in acetone at −20 °C, followed by three washes with 95 °C citrate buffer (10 mM sodium citrate, pH 6.0). Sections were blocked with 5 % goat serum and 0.1 % BSA for 20 min before the sections were incubated with the paralog-specific antibodies (1:400 dilution). The primary antibodies were detected using FITC anti-rabbit secondary antibodies (1:800 dilution, Sigma-Aldrich). Sections were mounted with Fluoromount aqueous mounting medium (Sigma-Aldrich), and immunofluorescence was observed and documented with a Zeiss imager.z1 microscope (Carl Zeiss MicroImaging, S.L.).

For Western blot analyses, aquaporin-injected oocytes were homogenized in HbA-buffer (20 mM Tris, pH 7.4, 5 mM MgCl2, 5 mM NaH2PO4, 1 mM EDTA, 80 mM sucrose, and protease inhibitor cocktail [Mini EDTA-free; Roche], 20 μL/oocyte), followed by two centrifugations for 5 min at 500 g where the supernatant was taken. Finally, a spin for 20 min at maximum speed pelleted the total membranes and this pellet was dissolved in 1x Laemmli buffer (5 μl per oocyte) and heated 10 min at 95 °C. The total membrane extract was loaded on a 12 % SDS gel and run at 90 V for 1.5 h. The separated peptides were transferred to nitro cellulose membranes and blocked for 1 h at room temperature in Tris buffered saline with Tween 20 (TBST; 20mM Tris, 140 mM NaCl, 0.1 % Tween, pH 7.6) with 5 % non-fat milk powder and subsequently incubated with primary antibodies (1:1000–1:2000 dilutions). Bound antibodies were detected by horseradish peroxidase-coupled anti-rabbit IgG and visualized by using enhanced chemiluminescence (ECL; Picomax, Rockland) as described previously [92].

RT-PCR analyses of aquaporin expression

Total RNA was extracted from RNA later stabilized L. salmonis from different developmental stages; fertilized eggs, nauplia 1 and 2, copepods, chalimus 1 and 2, preadult 1 and 2 males and females, and adult males and females with the RNAeasy Mini Kit (Qiagen). The RNA was DNAse treated before 80 ng was used for cDNA synthesis (as described above). 2 ul of 20-fold dilluted cDNA was used for template in 25 μL RT-PCR reactions using Gotaq DNA polymerase (Promega) and primer concentrations of 0.4 μM. (Primers are listed in Additional file 1: Table S3). 5uL of each PCR reaction (35 cycles) was run on 1 % agarose gels stained with GelredTM (Biotum) in TAE-buffer, and PCR expression products were visualised with UV-light in a gel doc imaging system (Biorad). To verify specific amplification the PCR products were sequenced following the BigDye protocol (v3.1).

Availability of supporting data

All the supporting data are included as additional files.

Abbreviations

ANOVA:

Analysis of variance

Aqp:

Aquaporin

ar/R:

Aromatic-arginine constriction

Bib:

Big brain

BLAST:

Basic local alignment search tool

cAMP:

Cyclic adenosine monophosphate

Drip:

Drosophila intrinsic protein

Eglp:

Entomoglyceroporin

ELISA:

Enzyme linked immunosorbant assay

Glp:

Aquaglyceroporin

MAFFT:

Multiple alignment using fast fourier transform

MBS:

Barth’s culture medium

MCMC:

Markov chain Monte Carlo

NPA:

Asparagine-Proline-Alanine motif

ORF:

Open reading frame

PBS:

Phosphate-buffered saline

PCR:

Ploymerase chain reaction

PripL:

Pyrocoelia rufa integral protein-like

RACE:

Rapid amplification of cDNA ends

RT:

Reverse transcriptase

TBST:

Tris buffered saline with Tween 20

t-coffee:

Tree-based consistency objective function for alignment evaluation

TMD:

Transmembrane domain

UTR:

Untranslated region.

References

  1. 1.

    King LS, Kozono D, Agre P. From structure to disease: the evolving tale of aquaporin biology. Nat Rev Mol Cell Bio. 2004;5:687–98.

    CAS  Article  Google Scholar 

  2. 2.

    Murata K, Mitsuoka K, Hirai T, Walz T, Agre P, Heymann JB, et al. Structural determinants of water permeation through aquaporin-1. Nature. 2000;407:599–605.

  3. 3.

    Agre P, King LS, Yasui M, Guggino WB, Ottersen OP, Fujiyoshi Y, et al. Aquaporin water channels - from atomic structure to clinical medicine. J Physiol-London. 2002;542:3–16.

  4. 4.

    Yakata K, Hiroaki Y, Ishibashi K, Sohara E, Sasaki S, Mitsuoka K, et al. Aquaporin-11 containing a divergent NPA motif has normal water channel activity. Biochim Biophys Acta-Biomembr. 2007;1768:688–93.

  5. 5.

    Ikeda M, Andoo A, Shimono M, Takamatsu N, Taki A, Muta K, et al. The NPC motif of aquaporin-11, unlike the NPA motif of known aquaporins, is essential for full expression of molecular function. J Biol Chem. 2011;286:3342–50.

  6. 6.

    Fu DX, Libson A, Miercke LJW, Weitzman C, Nollert P, Krucinski J, et al. Structure of a glycerol-conducting channel and the basis for its selectivity. Science. 2000;290:481–6.

  7. 7.

    Sui HX, Han BG, Lee JK, Walian P, Jap BK. Structural basis of water-specific transport through the AQP1 water channel. Nature. 2001;414:872–8.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Abascal F, Irisarri I, Zardoya R. Diversity and evolution of membrane intrinsic proteins. Biochim Biophys Acta-Gen Subj. 2014;1840:1468–81.

    CAS  Article  Google Scholar 

  9. 9.

    Finn RN, Chauvigné F, Hlidberg JB, Cutler CP, Cerdà J. The lineage-specific evolution of aquaporin gene clusters facilitated tetrapod terrestrial adaptation. PLoS One. 2014;9(11):e113686.

    Article  PubMed Central  PubMed  Google Scholar 

  10. 10.

    Cerdà J, Finn RN. Piscine aquaporins: an overview of recent advances. J Exp Zool Part A. 2010;313A:623–50.

    Article  Google Scholar 

  11. 11.

    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.

    Article  PubMed Central  PubMed  Google Scholar 

  12. 12.

    Finn RN, Cerdà J. Aquaporin evolution in fishes. Front. Physiol. 2011;2:441–6

  13. 13.

    Soto G, Alleva K, Amodeo G, Muschietti J, Ayub ND. New insight into the evolution of aquaporins from flowering plants and vertebrates: orthologous identification and functional transfer is possible. Gene. 2012;503:165–76.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Song J, Mak E, Wu B, Beitz E. Parasite aquaporins: current developments in drug facilitation and resistance. Biochim Biophys Acta-GenSubj. 2014;1840:1566–73.

    CAS  Article  Google Scholar 

  15. 15.

    Verkman AS, Anderson MO, Papadopoulos MC. Aquaporins: important but elusive drug targets. Nat Rev Drug Discov. 2014;13:259–77.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  16. 16.

    Beitz E. Aquaporins from pathogenic protozoan parasites: structure, function and potential for chemotherapy. Biol Cell. 2005;97:373–83.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Beitz E. Jammed traffic impedes parasite growth. Proc Natl Acad Sci U S A. 2007;104:13855–6.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  18. 18.

    Liu Y, Promeneur D, Rojek A, Kumar N, Frokiaer J, Nielsen S, et al. Aquaporin 9 is the major pathway for glycerol uptake by mouse erythrocytes, with implications for malarial virulence. Proc Natl Acad Sci U S A. 2007;104:12560–4.

  19. 19.

    Biyani N, Mandal S, Seth C, Saint M, Natarajan K, Ghosh I, et al. Characterization of Leishmania donovani aquaporins shows presence of subcellular aquaporins similar to tonoplast Intrinsic proteins of plants. PLoS One. 2011;6(9):e24820.

  20. 20.

    Mandal G, Orta JF, Sharma M, Mukhopadhyay R. Trypanosomatid aquaporins: roles in physiology and drug response. Diseases. 2014;22:3–23.

    Google Scholar 

  21. 21.

    Pietrantonio PV, Jagge C, Keeley LL, Ross LS. Cloning of an aquaporin-like cDNA and in situ hybridization in adults of the mosquito Aedes aegypti (Diptera:Culicidae). Insect Mol Biol. 2000;9:407–18.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Duchesne L, Hubert JF, Verbavatz JM, Thomas D, Pietrantonio PV. Mosquito (Aedes aegypti) aquaporin, present in tracheolar cells, transports water, not glycerol, and forms orthogonal arrays in Xenopus oocyte membranes. Eur J Biochem. 2003;270:422–9.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Ball A, Campbell EM, Jacob J, Hoppler S, Bowman AS. Identification, functional characterization and expression patterns of a water-specific aquaporin in the brown dog tick, Rhipicephalus sanguineus. Insect Biochem Mol Biol. 2009;39:105–12.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Campbell EM, Burdin M, Hoppler S, Bowman AS. Role of an aquaporin in the sheep tick Ixodes ricinus: assessment as a potential control target. Int J Parasit. 2010;40:15–23.

    CAS  Article  Google Scholar 

  25. 25.

    Drake LL, Boudko DY, Marinotti O, Carpenter VK, Dawe AL, Hansen IA. The aquaporin gene family of the yellow fever mosquito, Aedes aegypti. PLoS One. 2010;5(12):e15578.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  26. 26.

    Liu K, Tsujimoto H, Cha S-J, Agre P, Rasgon JL. Aquaporin water channel AgAQP1 in the malaria vector mosquito Anopheles gambiae during blood feeding and humidity adaptation. Proc Natl Acad Sci U S A. 2011;108:6062–6.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  27. 27.

    Tsujimoto H, Liu K, Linser PJ, Agre P, Rasgon JL. Organ-specific splice variants of aquaporin water channel AgAQP1 in the malaria vector Anopheles gambiae. PLoS One. 2013;8(9):e75888.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  28. 28.

    Guerrero FD, Andreotti R, Bendele KG, Cunha RC, Miller RJ. Yeater K, de L, Adalberto A Perez: Rhipicephalus (Boophilus) microplus aquaporin as an effective vaccine antigen to protect against cattle tick infestations. Parasit Vectors. 2014;7:475.

    PubMed Central  PubMed  Google Scholar 

  29. 29.

    Drake LL, Rodriguez SD, Hansen IA. Functional characterization of aquaporins and aquaglyceroporins of the yellow fever mosquito, Aedes aegypti. Sci Rep. 2015;5:7795.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  30. 30.

    Fuchs S, Behrends V, Bundy JG, Crisanti A, Nolan T. Phenylalanine metabolism regulates reproduction and parasite melanization in the malaria mosquito. PLoS One. 2014;9(1):e84865.

    Article  PubMed Central  PubMed  Google Scholar 

  31. 31.

    Potter CJ. Stop the biting: targeting a mosquito's sense of smell. Cell. 2014;156:878–81.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Benoit JB, Hansen IA, Attardo GM, Michalkova V, Mireji PO, Bargul JL, et al. Aquaporins are critical for provision of water during lactation and intrauterine progeny hydration to maintain Tsetse fly reproductive success. Plos Neglect Trop Dis. 2014;8(4):e2517.

  33. 33.

    Benoit JB, Hansen IA, Szuter EM, Drake LL, Burnett DL, Attardo GM. Emerging roles of aquaporins in relation to the physiology of blood-feeding arthropods. J Comp Physiol B-Biochem Syst Environ Physiol. 2014;184:811–25.

    Article  Google Scholar 

  34. 34.

    Kikawada T, Saito A, Kanamori Y, Fujita M, Snigorska K, Watanabe M, et al. Dehydration-inducible changes in expression of two aquaporins in the sleeping chironomid, Polypedilum vanderplanki. Biochim Biophys Acta-Biomembr. 2008;1778:514–20.

  35. 35.

    Campbell EM, Ball A, Hoppler S, Bowman AS. Invertebrate aquaporins: a review. J Comp Physiol B-Biochem Syst Environ Physiol. 2008;178:935–55.

    CAS  Article  Google Scholar 

  36. 36.

    Beuron F, Lecahérec F, Guillam MT, Cavalier A, Garret A, Tassan JP, et al. Structural analysis of a MIP family protein from the digestive tract of Cicadella viridis. J Biol Chem. 1995;270:17414–22.

  37. 37.

    Dow JAT, Kelly DC, Davies SA, Maddrell SHP, Brown D. A novel member of the major intrinsic protein family in Drosophila: are aquaporins involved in insect Malpighian (renal) tubule secretion? J Physiol. 1995;489:P110–1.

    Google Scholar 

  38. 38.

    Le Cahérec F, Deschamps S, Delamarche C, Pellerin I, Bonnec G, Guillam MT, et al. Molecular cloning and characterization of an insect aquaporin-Functional comparison with aquaporin 1. Eur J Biochem. 1996;241:707–15.

  39. 39.

    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–407.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Shakesby AJ, Wallace IS, Isaacs HV, Pritchard J, Roberts DM, Douglas AE. A water-specific aquaporin involved in aphid osmoregulation. Insect Biochem Mol Biol. 2009;39:1–10.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Lee KS, Kim SR, Lee SM, Lee KR, Sohn HD, Jin BR. Molecular cloning and expression of a cDNA encoding the aquaporin homologue from the firefly, Pyrocoelia rufa. Korean J Entomol. 2001;31:269–79.

    CAS  Google Scholar 

  42. 42.

    Herraiz A, Chauvigné F, Cerdà J, Belles X, Piulachs M-D. Identification and functional characterization of an ovarian aquaporin from the cockroach Blattella germanica L. (Dictyoptera, Blattellidae). J Exp Biol. 2011;214:3630–8.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Rao Y, Jan LY, Jan YN. Similarity of the product of the Drosphila neurogenic gene big brain to transmembrane channel proteins. Nature. 1990;345:163–7.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Yanochko GM, Yool AJ. Regulated cationic channel function in Xenopus oocytes expressing Drosophila big brain. J Neurosci. 2002;22:2530–40.

    CAS  PubMed  Google Scholar 

  45. 45.

    Frazer IN. Sea-cage aquaculture, sea lice, and declines of world fish. Conserv Biol. 2008;23:599–607.

    Article  PubMed  Google Scholar 

  46. 46.

    Rosenberg AA. Aquaculture - the price of lice. Nature. 2008;451:23–4.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Hahnenkamp L, Fyhn HJ. The osmotic response of salmon louse Lepeoptheirus salmonis (Copepoda:Caligidae), during the transition from seawater to fresh water. J Comp Physiol. 1985;155B:357–65.

  48. 48.

    Schram FR. Comments on Crustacean Biodiversity and Disparity of Body Plans. In: Watling L, Thiel M, editors. The Natural History of Crustacea. Functional Morphology and Diversity. Oxford: Oxford University Press; 2013. p. 1–33.

    Google Scholar 

  49. 49.

    Verkman AS, Mitra AK. Structure and function of aquaporin water channels. Am J Physiol-Renal Physiol. 2000;278:F13–28.

    CAS  PubMed  Google Scholar 

  50. 50.

    Froger A, Tallur B, Thomas D, Delamarche C. Prediction of functional residues in water channels and related proteins. Protein Sci. 1998;7:1458–68.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  51. 51.

    Regier JC, Shultz JW, Kambic RE. Pancrustacean phylogeny: hexapods are terrestrial crustaceans and maxillopods are not monophyletic. Proc R Soc B-Biol Sci. 2005;272:395–401.

    Article  Google Scholar 

  52. 52.

    Von R, Bjoern M, Jenner RA, Wills MA, Dell'Ampio E, Pass G, et al. Pancrustacean phylogeny in the light of new phylogenomic data: support for Remipedia as the possible sister group of Hexapoda. Mol Biol Evol. 2012;29:1031–45.

    Article  Google Scholar 

  53. 53.

    Rao Y, Bodmer R, Jan LY, Jan YN. The big brain gene of Drosophila functions to control the number of neuronal precursors in the peripheral nervous system. Development. 1992;116:31–3.

    CAS  PubMed  Google Scholar 

  54. 54.

    Burris PA, Zhang Y, Rusconi JC, Corbin V. The pore-forming and cytoplasmic domains of the neurogenic gene product, BIG BRAIN, are conserved between Drosophila virilis and Drosophila melanogaster. Gene. 1998;206:69–76.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Fabrick JA, Pei J, Hull JJ, Yool AJ. Molecular and functional characterization of multiple aquaporin water channel proteins from the western tarnished plant bug, Lygus hesperus. Insect Biochem Mol Biol. 2014;45:125–40.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Goto SG, Philip BN, Teets NM, Kawarasaki Y, Lee Jr RE, Denlinger DL. Functional characterization of an aquaporin in the Antarctic midge Belgica antarctica. J Insect Physiol. 2011;57:1106–14.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Azuma M, Nagae T, Maruyama M, Kataoka N, Miyake S. Two water-specific aquaporins at the apical and basal plasma membranes of insect epithelia: molecular basis for water recycling through the cryptonephric rectal complex of lepidopteran larvae. J Insect Physiol. 2012;58:523–33.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Ibanez F, Hancock J, Tamborindeguy C. Identification and expression analysis of aquaporins in the potato psyllid, Bactericera cockerelli. PLoS One. 2014;9(10):e111745.

    Article  PubMed Central  PubMed  Google Scholar 

  59. 59.

    Finn RN, Chauvigné F, Stavang JA, Belles X, Cerdà J. Insect glycerol transporters evolved by functional co-option and gene replacement. Nat Commun. 2015;6:7814.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  60. 60.

    Zardoya R. Phylogeny and evolution of the major intrinsic protein family. Biol Cell. 2005;97:397–414.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Grohme MA, Mali B, Wełnicz W, Michel S. The aquaporin channel repertoire of the tardigrade Milnesium tardigradum. Bioinformatics Biol Insights. 2013;7:153–65.

    Article  Google Scholar 

  62. 62.

    Gentles AJ, Karlin S. Why are human G-protein-coupled receptors predominantly intronless. Trends Genet. 1999;15:47–9.

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Farlora R, Araya-Garay J, Gallardo-Escarate C. Discovery of sex-related genes through high-throughput transcriptome sequencing from the salmon louse Caligus rogercresseyi. Mar Genom. 2014;15:85–93.

    Article  Google Scholar 

  64. 64.

    Tipsmark CK, Sorensen KJ, Madsen SS. Aquaporin expression dynamics in osmoregulatory tissues of Atlantic salmon during smoltification and seawater acclimation. J Exp Biol. 2010;213:368–79.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Madsen SS, Olesen JH, Bedal K, Engelund MB, Velasco-Santamaria YM, Tipsmark CK. Functional characterization of water transport and cellular localization of three aquaporin paralogs in the salmonid intestine. Front Physiol. 2011;2:56.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  66. 66.

    Engelund MB, Chauvigné F, Christensen BM, Finn RN, Cerdà J, Madsen SS. Differential expression and novel permeability properties of three aquaporin 8 paralogs from seawater-challenged Atlantic salmon smolts. J Exp Biol. 2013;216:3873–85.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Chauvigné F, Zapater C, Stavang JA, Taranger GL, Cerdà J, Finn RN. The pH sensitivity of Aqp0 channels in tetraploid and diploid teleosts. FASEB J. 2015;29:2172–84.

    Article  PubMed Central  PubMed  Google Scholar 

  68. 68.

    Berthelot C, Brunet F, Chalopin D, Juanchich A, Bernard M, Noel B, et al. The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates. Nat Commun. 2014;5:3657.

  69. 69.

    Macqueen DJ, Johnston IA. A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification. Proc R Soc B-Biol Sci. 2014;281:20132881.

    Article  Google Scholar 

  70. 70.

    Kuwahara M, Fushimi K, Terada Y, Bai LQ, Marumo F, Sasaki S. cAMP-dependent phosphorylation stimulates water permeability of aquaporin-collecting duct water channel protein expressed in Xenopus oocytes. J Biol Chem. 1995;270:10384–7.

    CAS  Article  PubMed  Google Scholar 

  71. 71.

    van Balkom BWM, Savelkoul PJM, Markovich D, Hofman E, Nielsen S, van der Suijs P, et al. The role of putative phosphorylation sites in the targeting and shuttling of the aquaporin-2 water channel. J Biol Chem. 2002;277:41473–9.

  72. 72.

    Valenti G, Procino G, Tamma G, Carmosino M, Svelto M. Minireview: Aquaporin 2 trafficking. Endocrinology. 2005;146:5063–70.

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Suzuki M, Tanaka S. Molecular and cellular regulation of water homeostasis in anuran amphibians by aquaporins. Comp Biochem Physiol A-Mol Integr Physiol. 2009;153:231–41.

    Article  PubMed  Google Scholar 

  74. 74.

    Saitoh Y, Ogushi Y, Shibata Y, Okada R, Tanaka S, Suzuki M. Novel vasotocin-regulated aquaporins expressed in the ventral skin of semiaquatic anuran amphibians: Evolution of cutaneous water-absorbing mechanisms. Endocrinology. 2014;155:2166–77.

    Article  PubMed  Google Scholar 

  75. 75.

    Johansson I, Karlsson M, Johanson U, Larsson C, Kjellbom P. The role of aquaporins in cellular and whole plant water balance. Biochim Biophys Acta-Biomembr. 2000;1465:324–42.

    CAS  Article  Google Scholar 

  76. 76.

    Prak S, Hem S, Boudet J, Viennois G, Sommerer N, Rossignol M, et al. Multiple phosphorylations in the C-terminal tail of plant plasma membrane aquaporins. Mol Cell Proteomics. 2008;7:1019–30.

  77. 77.

    Verdoucq L, Rodrigues O, Martiniere A, Luu DT, Maurel C. Plant aquaporins on the move: reversible phosphorylation, lateral motion and cycling. Curr Opin Plant Biol. 2014;22:101–7.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Blom N, Sicheritz-Ponten T, Gupta R, Gammeltoft S, Brunak S. Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics. 2004;4:1633–49.

    CAS  Article  PubMed  Google Scholar 

  79. 79.

    Neely JD, Christensen BM, Nielsen S, Agre P. Heterotetrameric composition of aquaporin-4 water channels. Biochemistry. 1999;38:11156–63.

    CAS  Article  PubMed  Google Scholar 

  80. 80.

    Beitz E, Liu K, Ikeda M, Guggino WB, Agre P, Yasui M. Determinants of AQP6 trafficking to intracellular sites versus the plasma membrane in transfected mammalian cells. Biol Cell. 2006;98:101–9.

    CAS  Article  PubMed  Google Scholar 

  81. 81.

    Preston GM, Jung JS, Guggino WB, Agre P. The mercury sensitive residue at Cysteine 189 in the CHIP28 water channel. J Biol Chem. 1993;268:17–20.

    CAS  PubMed  Google Scholar 

  82. 82.

    Zhang RB, Vanhoek AN, Biwersi J, Verkman AS. A point mutation at cysteine 189 blocks the water permeability of rat kidney water channel CHIP28k. Biochemistry. 1993;32:2938–41.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Hirano Y, Okimoto N, Kadohira I, Suematsu M, Yasuoka K, Yasui M. Molecular mechanisms of how mercury inhibits water permeation through aquaporin-1: understanding by molecular dynamics simulation. Biophys J. 2010;98:1512–9.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  84. 84.

    Echevarria M, Windhager EE, Tate SS, Frindt G. Cloning and expression of AQP3, a water channel from the medullary collecting duct of rat kidney. Proc Natl Acad Sci U S A. 1994;91:10997–1001.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  85. 85.

    Yukutake Y, Tsuji S, Hirano Y, Adachi T, Takahashi T, Fujihara K, et al. Mercury chloride decreases the water permeability of aquaporin-4-reconstituted proteoliposomes. Biol Cell. 2008;100:355–63.

  86. 86.

    Maurel C, Reizer J, Schroeder JI, Chrispeels MJ. The vacuolar membrane protein gamma-TIP creates water specific channels in Xenopus oocytes. Embo J. 1993;12:2241–7.

    CAS  PubMed Central  PubMed  Google Scholar 

  87. 87.

    Daniels MJ, Mirkov TE, Chrispeels MJ. The plasma-membrane of Arabidopsis thaliana contains a mercury-insensitive aquaporin that is a homolog of the tonoplast water channel protein TIP. Plant Physiol. 1994;106:1325–33.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  88. 88.

    Daniels MJ, Chaumont F, Mirkov TE, Chrispeels MJ. Characterization of a new vacuolar membrane aquaporin sensitive to mercury at a unique site. Plant Cell. 1996;8:587–99.

    CAS  Article  PubMed Central  PubMed  Google Scholar 

  89. 89.

    Tatsumi K, Tsuji S, Miwa H, Morisaku T, Noriya M, Orihara M, et al. Drosophila big brain does not act as a water channel, but mediates cell adhesion. FEBS Lett. 2009;583:2077–82.

  90. 90.

    Deen PMT, Verdijk MAJ, Knoers NVAM, Wieringa B, Monnens LAH, Vanos CH, et al. Requirement of human renal water channel aquaporin-2 for vasopressin-dependent concentration of urine. Science. 1994;264:92–5.

  91. 91.

    Finn RN, Kristoffersen BA. Vertebrate vitellogenin gene duplication in relation to the "3R hypothesis": correlation to the pelagic egg and the oceanic radiation of teleosts. PLoS One. 2007;2:e169.

    Article  PubMed Central  PubMed  Google Scholar 

  92. 92.

    Zapater C, Chauvigné F, Norberg B, Finn RN, Cerdà J. Dual neofunctionalization of a rapidly evolving aquaporin-1 paralog resulted in constrained and relaxed traits controlling channel function during meiosis resumption in teleosts. Mol Biol Evol. 2011;28:3151–69.

    CAS  Article  PubMed  Google Scholar 

  93. 93.

    Applebaum SL, Finn RN, Faulk CK, Holt GJ, Nunez BS. Developmental expression, differential hormonal regulation and evolution of thyroid and glucocorticoid receptor variants in a marine acanthomorph teleost (Sciaenops ocellatus). Gen Comp Endocrinol. 2012;176:39–51.

    CAS  Article  PubMed  Google Scholar 

  94. 94.

    Notredame C, Higgins DG, Heringa J. T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000;302:205–17.

    CAS  Article  PubMed  Google Scholar 

  95. 95.

    Katoh K, Toh H. Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform. 2008;9:286–98.

    CAS  Article  PubMed  Google Scholar 

  96. 96.

    Han MV, Zmasek CM.phyloXML: XML for evolutionary biology and comparative genomics. BMC Bioinformatics. 2009;10:356.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the Research Council of Norway (RCN) project 204813/F20 to R.N.F., the SFI-Sea Lice Research Centre, Grant No. 203513 to F.N., and the Spanish Ministry of Economy and Competitiveness (MINECO), project AGL2013-41196-R to JC. J.A.S. and F.C. were supported by postdoctoral fellowships from the RCN (204813/F20 and 224816/E40, respectively).

Author information

Affiliations

Authors

Corresponding author

Correspondence to Roderick Nigel Finn.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JC, FN and RNF, designed the experiments. RNF and JAS conducted the bioinformatic and phylogenetic analyses. JAS and HK conducted the cloning and expression experiments. JAS and FC performed the X. laevis swelling experiments. JAS, FC, HK, JC, FN and RNF analysed the data. FN contributed new reagents/analytical tools. JAS, JC and RNF wrote the manuscript. All authors discussed the data and commented on the manuscript. All authors read and approved the final manuscript.

Additional file

Additional file 1: Figure S1.

Extracellular views of cartoon renders of salmon louse aquaporins illustrating the ar/R constriction residues (spacefill) that determine the molecular selectivity of the channel. Labels are annotated as for Fig. 1 in the main text. Figure S2. Bayesian majority-rule mid-point rooted tree of an amino acid alignment of expressed transcripts in Maxillipoda. The tree was inferred from 500,000 MCMC generations on 8,085 amino acid sites. Posterior probabilities are shown at each node. Scale bar represents the number of amino acid substitutions per site. Grey bars highlight species with N-terminal splice-variants. Figure S3. Bayesian majority-rule mid-point rooted tree of a codon alignment of the aquaporin superfamilies in the Atlantic salmon and salmon louse. The tree was inferred from 1 million MCMC generations on 57,075 nucleotide sites. Posterior probabilities of the codon/amino acid analyses are shown at each node. Scale bar represents the number of nucleotide substitutions per site. Figure S4. In the absence of cAMP, salmon louse Glp1_v1 is non functional in X. laevis oocytes. (A) Osmotic water permeability (P f) of Glp1_v1 without the addition of cAMP compared to water- and glp_3v2-injected oocytes as negative and positve controls, respectively. (B) Immunofluorescence micrographs of paraffin sections of water- (Ctrl) and glp1_v1-injected oocytes in the absence of cAMP probed with paralog-specific antisera followed by Cy3-labeled anti-rabbit IgG. Arrows point to the plasma membrane. Inset: Magnified region shows the retention of Glp1_v1 (green) just below the oocyte plasma membrane. Table S1. C-terminal peptides used to immunize rabbits to raise salmon louse affinity-purified antibodies. Table S2. List of aquaporin accession numbers used in the study. Table S3. Oligonucleotide primers used for RT-PCR analysis. Nucleic acid sequences for primers specific for each salmon louse aquaporin mRNA, ef1α and expected product size. (PDF 2063 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Stavang, J.A., Chauvigné, F., Kongshaug, H. et al. Phylogenomic and functional analyses of salmon lice aquaporins uncover the molecular diversity of the superfamily in Arthropoda. BMC Genomics 16, 618 (2015). https://doi.org/10.1186/s12864-015-1814-8

Download citation

Keywords

  • Aquaporin
  • Aquaglyceroporin
  • Arthropod
  • Crustacea
  • Copepod
  • Salmon louse
  • Lepeoptheirus
  • Parasite
  • Atlantic salmon
  • Evolution
  • Osmoregulation
  • Fluid homeostasis
  • Permeability
  • Selectivity