Salmon louse (Lepeophtheirus salmonis) transcriptomes during post molting maturation and egg production, revealed using EST-sequencing and microarray analysis
- Christiane Eichner†1,
- Petter Frost†1, 5Email author,
- Bjarte Dysvik2,
- Inge Jonassen2, 4,
- Bjørn Kristiansen3 and
- Frank Nilsen1, 6
© Eichner et al; licensee BioMed Central Ltd. 2008
Received: 30 August 2007
Accepted: 10 March 2008
Published: 10 March 2008
Lepeophtheirus salmonis is an ectoparasitic copepod feeding on skin, mucus and blood from salmonid hosts. Initial analysis of EST sequences from pre adult and adult stages of L. salmonis revealed a large proportion of novel transcripts. In order to link unknown transcripts to biological functions we have combined EST sequencing and microarray analysis to characterize female salmon louse transcriptomes during post molting maturation and egg production.
EST sequence analysis shows that 43% of the ESTs have no significant hits in GenBank. Sequenced ESTs assembled into 556 contigs and 1614 singletons and whenever homologous genes were identified no clear correlation with homologous genes from any specific animal group was evident. Sequence comparison of 27 L. salmonis proteins with homologous proteins in humans, zebrafish, insects and crustaceans revealed an almost identical sequence identity with all species.
Microarray analysis of maturing female adult salmon lice revealed two major transcription patterns; up-regulation during the final molting followed by down regulation and female specific up regulation during post molting growth and egg production. For a third minor group of ESTs transcription decreased during molting from pre-adult II to immature adults. Genes regulated during molting typically gave hits with cuticula proteins whilst transcripts up regulated during post molting growth were female specific, including two vitellogenins.
The copepod L.salmonis contains high a level of novel genes. Among analyzed L.salmonis proteins, sequence identities with homologous proteins in crustaceans are no higher than to homologous proteins in humans. Three distinct processes, molting, post molting growth and egg production correlate with transcriptional regulation of three groups of transcripts; two including genes related to growth, one including genes related to egg production. The function of the regulated transcripts is discussed in relation to post molting morphological changes in adult female salmon louse. There is clear evidence that transcription of the major yolk proteins is not induced before some of the post molting growth of abdomen and the genital segment has occurred. A hallmark for the observed growth is transcription of many putative cuticula proteins prior to the size increase.
Copepods are arthropods in the aquatic environment and the extremely abundant free-living species are an essential part of the first levels of the marine food chain. Although copepods comprise the largest animal biomass on the earth, relatively limited biological information is available at the molecular level and no model species exists.
The salmon louse (L. salmonis) is an ectoparasitic copepod feeding on skin, mucous and blood from salmonid hosts. Recently it was shown that L. salmonis infections in farmed fish induce epizootics in wild fish [1, 2]. The life cycle of L.salmonis consists of 10 developmental stages separated by ecdysis [3, 4] and after the final molt, females develop into mature adults that continuously produce eggs for life. The first free-living larvae (naupli I) hatch directly from egg-strings attached to adult females and all three larval stages (naupli I, naupli II and the infectious copepidid stage) can be transported by the ocean currents over large distances depending on hydrographical conditions . After host settlement the infectious copepodids stage molt into chalimus. The four chalimus stages, all separated by molting, are anchored to the host by a frontal filament , which restricts the feeding area. However, in the succeeding pre-ad I and -II and adult stages the salmon louse can move unrestricted on the host surface resulting in increased virulence .
Sexual maturation and vitellogenesis are major physiological and behavioral changes in most animal life cycles. Germ cells are typically established early in development but arrested in development until the onset of sexual maturation. The generation of gametes is most conserved in males whereas variation is seen between different species for the development of female gametes. The ovum (i.e. the mature unfertilized egg) is a highly complex cell that is energetically expensive to produce. In order to produce high quality ova the females must undergo physiological adaptations that initiate further gamete development and maturation. Since the reproduction strategy is highly variable between different species and different life strategies (e.g. free living or parasitic) the processes of sexual maturation and the production of eggs also varies. However, there are some common hallmarks that are expected in most animals. After fertilization, the egg must contain sufficient energy to ensure development until external energy sources can be utilized. The ovarium is the site of initial development for female gametes during sexual maturation and reproduction. At ovulation, oogonia are released from the ovarium into the oviduct where growth and maturation take place. This process is typically divided into pre-vitellogenic and vitellogenic development. During vitellogenesis, yolk proteins are incorporated into the oocytes. A molecular hallmark for this process is the transcription of genes encoding egg yolk proteins like vitellogenins (Vgs). Depending on animal group, the transcription of Vgs takes place in different cell types like liver (in vertebrates), fat body (in insects) and hepatopancrease (e.g. in decapods). The Vgs are subsequently transported by the blood or hemolymph to the maturing oocytes, where they are taken up by receptor-mediated endocytosis [8, 9]. Production of vitellogenin is controlled by steroid hormones, which induce transcription of the target gene through binding of a steroid-receptor complex to the gene promoter. In arthropods, including crustaceans, ecdyson (i.e. E20) has been shown to bind to the heterodimeric ultra-spinacle (UsP) and ecdysteroid receptor (EcR) to an ecdysteroid response element (ERE) in the vitellogenin promoter . Ecdysteroids are also a key regulatory component in arthropods molting and development [11, 12].
The salmon louse reproductive systems have been described at the anatomical and histological levels  but there is no information regarding the timing of the different events during sexual maturation. It has been proposed that males depositing spermatophores triggers the egg-production in other parasitic copepods (Lepeophtheirus pectoralis) , but according to our observations using unfertilized laboratory animals, adult female L. salmonis produce eggs and external egg string also when males are not present (pers. obs., present study). Unlike crustaceans like shrimp, that produce eggs in seasons and grow/molt their entire life, salmon louse have a final molting, stop growing when egg production has started and then continuously produce eggs for life. It has been shown that L. salmonis can produce up to 11 sets of egg-strings from a single fertilization . However, immediately following the last molting, the adult female salmon louse is not fully developed. Prior to egg production the animal matures in a process that includes a large increase of the genital segment and the abdomen, whereas the frontal cephalothorax appears unchanged (present study).
In order to link transcripts to the morphological and anatomical changes that takes place during the transition from pre-adult II to egg producing females we have combined EST sequencing and microarray analysis. The microarray analysis revealed three distinct groups of transcripts that correspond to molting, post molting growth and egg production. The possible function of the regulated transcripts is discussed in relation to the anatomical and physiological changes taking place. Initial analysis of EST sequences of L. salmonis revealed a large proportion of transcripts with no significant hits in public databases. In order to obtain some initial information regarding copepod proteomes we compared at set of L. salmonis proteins to some selected crustaceans, insects and vertebrates and shows that the salmon louse proteins are equally similar to all species.
Overview of isolated clones and sequenced ESTs from the various cDNA libraries
Percent of total
Percent of total
Pre-adult II female
Female male subtracted
Summary of data obtained after EST assembly using the ContigExpress module in Vector NTI.
Number of ESTs in contigs:
Number of contigs:
Number of singletons:
Number of contigs containing:
> 20 ESTs
About 1,900 clones were not included in the contig assembly, primarily due to lack of high quality sequence data but also due to insert sizes less than 100 bp (< 1%) or empty clones/E. coli sequences (1%). Based on these considerations, the proportion of singletons (35%) in the clustering results and the proportion of mitochondrial sequences (6%), the number of additional salmon louse transcripts among the 1,900 clones without sequence data were estimated to be approximately 500. This indicates that the 7,021 clones, from which cDNA probes were printed on the microarray, represent up to 2,600 different L. salmonis transcripts.
Sequence identity analyses of L.salmonis proteins with homologous proteins in crustaceans, insects, fish and humans
40S ribosomal protein S13
60S ribosomal protein L5
60S ribosomal protein L7A
translation initiation factor 4A2 isoform 2
elongation factor 1-alpha
receptor for activated protein kinase C-like
14-3-3-like protein (Leonardo protein)
S5e ribosomal protein
phospholipid-hydroperoxide glutathione peroxidase
Adenosine kinase (AK)
ER protein disulfide isomerase
serine-type enodpeptidase (trypsin-lik)
cysteine protease (cathepsin)
heat shock protein 70
Post-molting growth of adult female L.salmonis
Microarray analysis was performed on a total of 34 female L. salmonis, representing the pre-ad II stage and 6 different time points after the last molting (Figure 2). ESTs were grouped according to expression profiles using self-organizing map (SOM) [21, 22] (25 cluster (5 × 5) neuron map). The SOM revealed three clusters with a substantial increase in gene transcription during adult development (clusters 1, 2 and 6) and two (clusters 24 and 25) with a significant decrease (see Additional file 1). All other SOM clusters showed very little variance within clusters. Although the groups sampled are not clear-bordered developmental stages but rather transitions in time, correspondence analysis (CA)  of the complete dataset (see Additional file 2) revealed consistency among the 5 biological parallels selected according to the genital segment/cephalothorax ratio (see Figure 2). With the exception of one obvious pre-ad II outlier, samples were located in adjacency with other samples from the same group or in the vicinity of a neighboring time point. The developmental group membership, and the single outlier, of individual samples were also evident in a 10 × 10 neuron heat SOM (see Additional file 3).
The majority of the annotated group 2 and 3 "early genes" (T1 or pre-ad II transcribed) were similar to cuticula proteins (22 and 40% respectively), especially Barnacle cypris larva-specific gene (BCS-1) from Balanus amphitrite , a marine crustacean, but the E-values were high (see Additional file 4). In addition, 8% of the group 2 genes were similar to genes involved in other growth related processes. Most group 2 contigs were also down regulated in adult male lice compared to pre-adII males (Figure 4), pointing towards sex independent functions. However, in group 3 (down regulated in females from pre-ad II to adult T1) all but one gene was unregulated in males.
All the salmon louse contigs homologous to cuticula-like proteins were compared to each other and based on the size of the encoded protein they separated into two sub-groups. These sub-groups could be divided into several other groups based on sequence similarity. Some short proteins (HA608, HA755, Contig137, Contig511, PU731 and PU186) are similar to BCS-1 while other cuticula proteins (Contig488, Contig34new, HA878) are similar to cuticula protein 6 in insects. Aligning these sequences to a selection of sequences from insects revealed some degree of conservation (a few stretches of conserved amino acids) indicating similar functions. The Contig166 and Contig527, similar to horseshoe crab (Tachypleus tridentatus) cuticular protein have a conserved domain, a chitin binding Peritrophin-A domain (pfam 01607.11) found in chitin binding proteins in the perotrophic matrix proteins of insects [25–27]. Interestingly, a perotrophic membrane is not described in L. salmonis.
The highest scoring sequences among the "late genes", up-regulated during adult female maturation, gave hits with genes involved in reproduction. Two large contigs (Contig552 and Contig276) showed high level of identity with two different arthropod vitellogenins [28, 29] and Contig256 had similarities to a nucleolar protein in starfish (Asterina pectinifera) only transcribed in growing oocytes . Furthermore, except for one gene (LsVit-like) all group 1 genes were up regulated in females only (Figure 4).
None of the 86 "early genes" (group 2 and 3) were in contigs consisting of more than 10 ESTs compared to 10 of the 34 "late genes" (Figure 4), including 3 genes encoding novel proteins (Contig471, Contig259 and Contig6). Five "late genes" (Contig12, Contig9, Contig241, Contig273 and LsVit-like), within the group of genes with miscellaneous functions, were also in large contigs (20–106 ESTs) indicating a high transcription level. The Contig9, that contains more ESTs (106) than both vitellogenins, gave strong hits with a Cathepsin L-associated protein in Artemia  mainly because it contains three FAS1 domains typically seen in cell adhesion molecules [32, 33].
The microarray transcription profiles were confirmed by Northern blot analysis using 2 genes in Group 1 (Contig387, Contig12) and 5 genes in Group 2 (Contig529, Contig488a, Contig34new, Contig533, Contig508) (See Additional file 5).
Salmon lice included in the present microarray study are female pre-adult II (before the final molt) and adults at different time points (T1–T6) during the post molting growth towards egg-production (see Figure 2). The microarray analysis of approximately 2,600 different transcripts revealed a total of 120 regulated genes divided into 3 patterns, down regulated from pre-ad II to adult T1 (group 3), up-regulated from pre-ad II to T1 and then down regulated (group 2), and up-regulated from T2-T6 (group 1). A typical feature of group 2 and 3 is transcripts resembling cuticula proteins, while the group 1 typically includes female specific transcripts, including vitellogenins, with many ESTs (up to 106) compared to the contigs in group 2 and 3. The diversity of transcripts is largest in group 2.
Cuticula encoding genes are transcribed prior to the final molt and prior to the post molting growth of adult female salmon louse
Animals with an exoskeleton, like the salmon louse, grow through molting and each instar typically shows limited increase in size. Developmental stages in salmon louse (3 free-living and 7 parasitic) are also separated by molting . However, the present study clearly demonstrates a substantial growth in females after the final molt. Post molt growth has been described in another parasitic copepod, Lernaeocera branchialis, where substantial growth and metamorphosis in adult females after the final molt have been reported, resulting in a 20-fold size increase of the abdomen . They observed extensive cuticular folds (4–6 um deep folds, with a density of 1–1.2 folds/um), associated with expansion of the cuticle following the final molt, which could account for some of the size increase (about 6-fold). Based on ultra structural observation they concluded that additional mechanisms (i.e. other than stretching of the cuticular folds), like large-scale cuticule secretion must account for the large size increase. In L. salmonis small folds (0.9 um) have been observed in the cuticula of chalimus larva , but currently there is no detailed description of the development of salmon louse after the final molt. Our preliminary observations on T1 and T6 lice indicate clear differences in cuticula structure, particularly on the genital segment and abdomen. Immature adult L. salmonis (T1) has cuticular folds but to what extent these folds can explain the large increases in size during the post-molt growth observed in the present study is not clear. Prior to the final molting (pre-ad II lice) and at the start of the post molting growth (T1) an increase in transcription of genes believed to be involved in growth related process in arthropods (e.g. cuticula proteins) was observed. During the post molting growth these transcripts were down regulated while female specific transcripts were up regulated. Egg strings, from which the offspring (planktonic nauphli larvae) is released, were not produced until this post molting growth was completed. This indicates that the post molting growth is connected to sexual maturation for egg production (increased genital segment) and an increased capacity for food (mainly blood) uptake (increased abdomen). The high proportion of cuticula and growth related genes identified in group-2 and -3 (25 and 40%, respectively) indicates that the group-2 and 3 genes, including the novel proteins (43% and 30%), primarily consists of genes involved in copepod growth. Many of the pre-ad II and T1 transcribed cuticula annotated genes (group 3 and 2, respectively) resemble BCS-1, known to be transcribed prior to Barnacle larval attachment and metamorphosis . For pre-ad II lice this could be related to the fact that the louse is attached to the salmon host through a frontal filament during molting. However, this is somewhat contradicted by the fact that the pre-ad II transcribed BCS-1 like genes are not regulated in males. Since both group 3 and 2 genes are down regulated after T1, when the post-molting growth of adult female salmon louse takes place (see Figure 2), growth from T1 to T6 is phenotypically delayed in time relative to transcription and probably do not include cell division growth. This is consistent with the fact that the new cuticula is produced beneath the old pre-ad II cuticula, prior to the molting when the group 3 genes are transcribed. Therefore in order to increase in size some kind of stretching/swelling must occur after the old cuticula has been shed, starting at T1 when the group 2 genes are transcribed. Based on the present study it is therefore likely that the group 2 and 3 cuticula proteins are an important part of a two-step process of cuticula formation at the final molt of L. salmonis.
Abdominal growth and gene regulation
After the final molt, the L. salmonis abdomen grows from approximately 0.8 mm to 2.8 mm (3.5 fold length increase). It is likely that this increases the capacity for food digestion, which could be visualized as increased transcriptions of digestive peptidases. The L. salmonis microarray contains more than 600 peptidase ESTs but only seven of the regulated genes (three in group 1 and four in group 2) were identified as proteases (one trypsin-like, one cysteine peptidase, two carboxypeptidases and three metallopeptidases). None of the putative digestive serine peptidases previously identified in the salmon louse [18, 19] were among the regulated transcripts. However, trypsin-like peptidases likely to be involved in digestion may also be regulated at the post-transcriptional level , and hence would not be detected in a transcriptome analysis.
Peptidases from the different families may be involved in a wide range of cellular and biological processes making it difficult to infer the function based on database searches. Host blood is a major component of the food for the salmon louse  particular for adult females. A wide range of organisms (both unicellular and metazoan species) utilize blood as nutrition. To keep blood in an easy accessible form metazoan species typically produce components with anticoagulation properties. In addition, it has been shown for several hematophagus organisms that digestion of blood and particular haemoglobin demands the action of a set of peptidases [38, 39]. It appears that aspartic- cysteine- and metallo-peptidases are key players in blood and haemoglobin digestion from different parasitic species . Of the regulated peptidases in the present study three are metallopetidases and one is cysteine peptidase. Contig127 (in group 2) is a cysteine peptidase resembling the trophozoite cysteine proteinase (TCP) from malaria (Plasmodium falciparum) . Interestingly, the TCP is believed to be involved in degradation of erythrocyte haemoglobin, which also could be the case in the salmon louse as they ingest blood from their salmon host. In addition, Contig5new, which resembles a Kuniz-like serine peptidase inhibitor, is up-regulated from T1 to T6. Serine peptidase inhibitors are involved in a wide range of biological processes including blood coagulation/anticoagulation, hence, it is possible that Contig5new could encode an anticoagulant, which is up regulated due to the abdominal growth and the corresponding increase in blood feeding during egg production.
Egg production and gene regulation
The female specific nature of all but one gene up-regulated during the T1–T6 growth of adult females (group 1) and the two identified vitellogenins, points towards genes involved in reproduction. No gene libraries were normalized and the high number of ESTs in many of the group 1 contigs, compared to group 2 and 3 contigs (Figure 4), therefore points towards functions where a high transcription level is typically seen. Both vitellogenin contigs (LsVit1 and LsVit2) contains approximately 100EST which is in compliance with the high transcription level typically seen for vitellogenins. The equally transcribed Contig9, encoding a protein containing 3 FAS1 domains, gave GenBank hits with CathepsinL-associated protein in Artemia where it is most abundantly expressed in encysted eggs and embryos . It is therefore likely that Contig9 encodes a highly expressed protein with cell adhesion properties involved in egg production. In insects, FAS1 containing proteins are primarily described as neural cell adhesion molecules but they have also been reported to be present on the surface of eggs and also to be expressed in several non-neural tissues in the embryo [32, 41]. Further information on the biological function of this gene and other genes in group 1–3 will be revealed through ongoing knock down studies.
Induction of the genes up-regulated during the post molting maturation of adult females appears to take place at time point T2 (see Figure 3). This is at the time when the first immature eggs are observed in the genital segment and 32% of the observed size increase has taken place. This, and the fact that similar growth appears in the abdomen, indicates that the final size of the genital segment in mature adults cannot be a result of stretching due to the presence of immature eggs only. From T3, the genital segment is gradually filled up with maturing eggs that are extruded as two separate strings (Figure 3). Based on these observations and the induction of a range of growth related genes (including cuticula proteins) at T1, it is likely that the observed size increase includes addition of more components into the exoskeleton rather than stretching.
Hormonal regulation and transcription profiles in the salmon louse
In the present study salmon lice go through three distinct processes, molting from pre-ad II to adult, post molting growth and egg production, three processes that correlates with regulation of the three groups of transcripts (group 3, 2 and 1, respectively). In arthropods, the ecdysteroid 20-hydroxyecdysone (20E) and juvenile hormone (JH) are key participants in regulating growth (molting), sexual maturation and egg production [12, 42]. In hematophagous mosquito a blood meal trigger the vitellogenesis, which is characterized by a specific gene expression pattern tightly controlled by the 20E . In Drosophila short day photoperiods at low temperatures induce ovarian diapause and suppressed ecdysteroid production . Temperature upshift results in an increase of steroid hormone levels, diapause termination and onset of vitellogenesis. Currently there is no available data about ecdysteroid levels in L. salmonis. In crustaceans (e.g. crabs, lobsters) information is also limited but it has been demonstrated that ecdysteroids are involved in molting, reproduction and embryogenesis . In addition, ecdysteroid levels in free living copepods are known to fluctuate during the life cycle . A recent microarray study in Drosophila assesed the genomic response of 20-hydroxyecdysone (20E) during metamorphosis . By using RNAi to knock down EcR they were able to identify EcR-dependent genes. A total of 4,188 genes were regulated in EcRi animals at one or more of the examined times points, indicating that approximately 30% of the Drosophila transcripts are affected by EcR directly or indirectly. Dana et al.  examined ~2,000 transcripts related to blood feeding in Anopheles gambiae and found that approximately 20% of the examined genes were regulated as a response to the blood meal. These regulated transcripts could be separated into early, middle and late responding genes, similar to the present study, and were suggested to be controlled by the fluctuating titers of JH and 20E either as inducers and/or repressors of transcription.
The observed post-molt growth in L. salmonis is likely to be regulated at the transcription level since distinct peaks in transcription of cuticula and growth related genes is evident in pre-ad II, where new cuticula is produced, and in newly molted adult females prior to the post molting growth. Recently, it has also been shown that transcription of cuticular proteins can be induced by 20E in the beetle Tenebrio monitor . Hence, it is likely that a 20E homolog is a key participant in inducing the transcription profiles observed in the present study, probably acting both as an inducer and repressor.
It has been proposed that fertilization is the trigger for the egg production (ecdysteriod inducer) in parasitic copepods . However, at least in L. salmonis this is not the case since in the present study normal post molt growth, including production of external egg strings (with unfertilized eggs), was observed in female only populations. Furthermore, since egg producing L. salmonis females can be found on salmon throughout the year at all temperatures and light conditions, environmental factors do not seem to be a likely candidate. Although, egg producing females feed extensively on blood, and this high energy food is probably important for egg production, blood can be part of the diet at all parasitic stages indicating that vitellogenesis is not triggered by a blood meal either.
Recently it was shown that in Drosophila growth and sexual maturation is coupled, coordinated by the DHR4 nuclear receptor . DHR4 is part of the ecdysone-triggered cascade both as a repressor (on the early 20E induced regulator genes) and an inducer of the β FTZ-F1 competence factor. By using DHR4 mutants it was suggested that DHR4 is involved in assessing critical weight and duration of larval development, linking 20E signaling to this process. It is possible that a similar molecular signaling is present in copepods and that the size of immature L. salmonis is the trigger for the onset of vitellogenesis. However, further studies are necessary to verify this hypothesis.
Many novel genes in the L. salmonisgenome
The initial analysis of individual L. salmonis EST sequences revealed that approximately 43% had no significant hits in GenBank using the nr-database. Even the re-annotation of the contigs identified by microarray analysis revealed a total of 41% novel proteins. Similar high proportion of novel genes was one of the interesting discoveries when the first non-mammalian genomes were sequenced. When the first metazoan animal (C. elegans) was completely sequenced the initial analysis revealed that more than 50% of the 19,000 genes were classified as novel . Similar figures were seen in other "early" sequenced genomes. Obviously the most important reason for the high number of novel genes was that few organisms were completely sequenced and, hence, few "complete" dataset were available for comparison. This is illustrated in the worm book , where at present only 12.3% of the C. elegans protein coding genes are classified as "function unknown". Currently, no copepod genome has been sequenced and the closest sequenced relatives to the salmon louse are insects like Drosophila and Anopheles. This is obviously the main reason for the high number of novel proteins in the salmon louse but considering the increasing number of sequenced genomes, 40–50% novelty indicates that the copepods have a distinct proteome (in terms of the number of novel proteins) compared to terrestrial arthropods. This is further indicated by the fact that searching the BaNG database  crustacean section with 35 novel proteins (see additional file 4) only resulted in six marginal (e-value 7E-6 to 3E-12) hits (data not shown). The large proportion of novel transcripts in L. salmonis was further demonstrated by the fact that of 4,586 validated nuclear L.salmonis ESTs 51% showed hits with insignificant values when compared to the only and recently available crustacean genome (D. pulex) .
Considering that reproduction processes have been intensively studied in many organisms, the number novel genes (38%) differentially transcribed during L. salmonis egg production (group 1 in the microarray study), was surprisingly high. The fact that conserved proteins in L.salmonis genes have the same low identity with homologous proteins in humans and the much closer related crustaceans and insects, further demonstrates the potential for discovery of unknown genes and "new" biological processes by sequencing genomes representing unstudied groups like copepods, instead of new species in already well studied groups.
We have analyzed 7021 ESTs from the parasitic copepod Lepeophtheirus salmonis. Annotation demonstrated high level of novel genes (40%) and sequence identity analysis demonstrated that the identity of Lepeophtheirus salmonis proteins with homologous proteins in crustaceans and terrestrial arthropods is not higher than with homologous proteins in fish and humans. We have used 7 k microarrays to identify 120 genes, of which 40% are novel, involved in sexual maturation and egg production.
Laboratory-reared strains of salmon louse were used. Infectious copepodids are produced in small buckets with flow-through water, starting with egg strings. The copepodids are allowed to infect naive Atlantic salmon kept in flow through water tanks. At relevant time points after infection, the salmon is anesthetized and salmon louse of interest collected and conserved for later use. In the present study adult fertilized individuals of L. salmonis were sampled during the period of maturation from pre-ad II to mature egg producing adults (app 14 days), photographed and stored on RNA later® according to the manufactures recommendations. In the microarray experiment adult animals were divided into six different morphological groups (T1#150;T6) based on the cephalothorax/genital complex length (mm) ratio (see Figure 2) and the degree of genital filling/eggstring development. Five (in one case four) lice from each group were used as biological parallels in the microarray analyses. In addition 5 pre-ad II females were defined as T0. Pre-ad II males and adult males were sampled for comparison between the sexes.
cDNA library construction and EST sequencing
Total RNA was extracted by Trizol (Sigma) and by RNeasy Mini Kit (Qiagen) followed by an enrichment of polyA RNA by Poly(A) Pure ™ or Poly(A) Purist ™ (Ambion) respectively. Two cDNA libraries (mRNA from female intestine and mRNA from the entire mature female) were constructed in lambda ZAP as described elsewhere . In addition, two cDNA libraries (pre-ad II female and adult male) were constructed directly in the pBluscript SK+ vector as described by the manufactures (Stratagene). A subtracted cDNA library was constructed by using the PCR-select cDNA subtraction kit (Clontech). Mature adult female lice were used as the source for tester cDNA whereas adult males were used as driver cDNA. The lambda ZAP libraries were mass excised according to the manufactures instructions (Stratagene). All clones were blue-white screened and white clones were picked randomly from all the different libraries for plasmid purification. Bacteria were grown overnight in 96 well dishes (Millipore) and plasmids were purified according to the recommendations from the manufacturer (Millipore). Clones were sequenced using vector primers T3, m13f or m13f and BigDye chemistry (Applied Biosystems). The Vector NTI software package was used for sequence handling and analysis unless otherwise stated. After retrieval, the EST sequences were trimmed for vector sequence and quality assessment using options in the ContigExpress. Each EST was trimmed until there were less than 3 ambiguous bases per 25 bases. All approved ESTs were further analyzed by NCBI BlastX and BlastN searches in GenBank using the non-redundant database. In the primary annotation process a significant hit was defined as one with an Expectation value (E-value) less than 5.5 × 10-5. Based on this analysis the ESTs were annotated as similar to the sequence showing the highest score. In addition, we searched the Crustacea section in the BaNG database  with the 27 protein sequences used for similarity calculations and we compared 4,586 validated nuclear ESTs with the recently available D. pulex genome  using Blast. All ESTs were then clustered by the Contig Express module in the Vector NTI package using a minimum overlap of 50 bp, an identity of 0.9 and a cut off score of 20. Prior to the clustering, identified mitochondrial genes were removed. For LsVit2, LsVit2 and LsVit-like, full-length sequences were obtained by supplementing EST contigs with 5'RACE clones using SMART™RACE (Clonetech) and primer walking sequencing. Sequence data from this study have been submitted to GenBank and accession numbers are available in additional file 4.
Microarray design and production
Probes were amplified from the individual cDNA clones by PCR using pBluescript-specific primers (TAATACGACTCACTATAGGGATAGGGCGAATTGGGTACCG and TAATACGACTCACTATAGGGAAAGGGAACAAAAGCTGGAGC). PCR reactions (100 μl) contained 10 μl 10 × reaction buffer (Promega), 160 μM MgCl2, 100 μM dNTPs, 0.15 μM of each primer and 2.5 U Taq U polymerase (Promega). An initial 2 min denaturation was followed by 30 PCR cycles (94°C for 30 sec, 60°C for 15 sec and 2 min elongation at 72°C) and a final 10 min elongation at 72°C. The PCR products were purified using Millipore Montage™ PCRμ96 Plate according to manufactures instructions. All purified probes were checked for size and purity by gel-electrophoresis (Invitrogen E-gel 96well 2% (GP). Furthermore, probe concentrations (70 ± 24 ng ng/μl) were estimated using Pico Green® dsDNA quantitation reagent (Molecular Probes) according to the manufactures protocol and measured on a Fluorstar Optima (BMG Labtechnoologies). Probes (70–120 ng/μl) in 50% DMSO were printed on Aminosilane coated slides (Corning® UltraGaps™) at 20–22°C and 45–55% relative humidity using a BioRobotics, Micro Grid II arrayer (Genomic Solutions®) with Mikrospot 10 K split pins. Slides were dried in a desiccator cabinet for 24–48 hours and DNA crosslinked at 350 mJ/cm2 using a UV-Stratalinker 2400, (Stratagene Inc.). The 7776 spots were printed in 24 subarrays, 7008 L. salmonis probes, 72 (triplicate on each subarray) L. salmonis reference gene EF1α , 360 negative controls mouse ESTs (15 on each subarray) and the rest were "no probe spots".
RNA isolation, cDNA labeling for microarray hybridization
RNA was isolated for individual animals using the RNAeasy Mini kit (Qiagen) according to the manufactures recommendations. The RNA was DNAse treated by TURBO DNA-free™ (Ambion) according to the supplied protocol and Superase (Ambion) was added. The RNA samples were frozen at -80°C until analysis. One aliquot was used for RNA integrity and quantity measures using the Agilent 2100 Bioanalyzer and NanoDrop Spectrophotometer (OD 260/280 and 260/230 ratios), respectively. Another aliquot was used for cDNA synthesis and labeling using Fair Play® Microarray Labeling Kit (Stratagene) according to the manufactures instructions. Typically, 10 ug total RNA was used for cDNA generation. In a few cases less RNA (down to 5 μg) were used due to lower RNA yields, typically from the smaller single pre-adult lice. Samples (Cy5 labeled) where hybridized against a common standard (Cy3 labeled) reverse transcribed and labeled in parallel. The reference RNA was a pooled and aliquoted mRNAs mix from adult females, pre-adult females, adult males and pre-adult males (1:1:1:1), tested and used as described for the samples. Labeling efficiency and quantity of labeled cDNA was determined using the NanoDrop Spectrophotometer and identically quantity of sample and standard was used in all hybridizations. Slides were pre-hybridized in 20 × SSC, 10% SDS and 1%BSA for about 45 min at 65°C followed by washing twice in water and once in isopropanol. Slides were dried by centrifugation in a mini-centrifuge. Sample and reference was unified, diluted in Tris buffer pH 8.0 and centrifuged using microcon YM-30 columns (Millipore). After sample denaturation (100°C, 2 min) hybridization was performed at 60°C over night with rotation using Agilent 2 × hybridization buffer (250 μl) in Agilent hybridization chambers. The slides were put in 2 × SSC/0.1%SDS at 65°C to remove gasket slide and then washed for 5 min in 1 × SSC at 65°C, for 5 min in 0.2 × SSC at RT, for 45 sec in 0.05 × SSC at RT, and spinned dry (mini-centrifuge). Dye swap quality control experiments were performed with two samples and technical replicates performed with five samples.
Slides were scanned directly after the washing procedure using an Agilent scanner at a resolution of 10 um with default settings. The scanned microarray images were analyzed using the GenePix Pro 6.0 software package and exported as image quantization files (gpr- and jpg-files). Each gpr- and jpg-file was then further processed in the J-Express (version 2.7) software package  and organized into a gene expression matrix where each row represents a clone, each column represents a sample (salmon louse) and each cell contains a log2 (sample vs. reference) ratio. The processing of image quantization data was performed by removing (1) all "empty" and "control" spots (696 spots) and all spots flagged by the GenePix software package or manually as "not found" or "bad" and (2) all spots with a signal to noise ratio (reported by GenePix) below 3.0 in both channels before a global lowess normalization  was applied to all remaining probes.
To prepare the gene expression matrix for further downstream analysis, we applied a KNN (k = 10) imputation to make sure an expression value existed for all genes in all samples . Prior to imputation, we removed all genes having more than 80% missing values, giving an expression matrix with 6653 clones. To remove ESTs not in the scope of this project, we used a self-organizing map to group ESTs according to expression profiles. All time courses of red/green ratios for all ESTs were sorted via self-organizing maps with Euclidian distance [21, 22] (25 cluster (5 × 5) neuron map). The self-organizing map shows clusters of ESTs sorted by cluster variance. To identify candidate gene clusters, which are most regulated over time the clusters with highest within cluster variance were selected. Using a trimmed group mean of samples from the same development stage a combined dataset was created. To further explore gene expression similarities corresponding to development stage, we applied a correspondence analysis to the complete dataset and to a group-combined expression matrix (samples within a development group was combined using a trimmed mean) . Developmental group membership was also explored by self-organizing map. Finally, probe quality of all spots from the chosen clusters was checked. All spots were the probe agarose gel electrophoresis test showed double band, smear or size below 400 bp was removed. Secondly, spots were filtered using cut off values for average signal to noise ratio (SNR 635, 532), spot intensity (Mean F635, F532) for both red and green channel, and for the percentage of feature pixels with intensities more than two standard deviations above the background intensity (at wavelength 635 or 532 (% > B635+2SD)). The cut-off values were set based on evaluation of all probes in contigs with more than 20 ESTs, setting the cut of at the level where individual probes gave inconsistent results compared to other probes representing the same contig. These cut off were: SNR< 40, F< 100 and > 10000, B+2SD< 60.
The consensus sequences of differentially expressed genes, edited contigs and re-sequenced singletons were checked for ORFs and re-annotated using BlastP (if ORF) and BlastX.
We provide MIAME-compliant description of the microarray study, available in the arrayexpress database , accession number E-BASE-5.
Total RNA from pre-adII, T1, T3, T5 (female) and pre-adII and adult male lice (3.2 μg) was mixed with Northern Max Formaldehyde Loading Dye (Ambion), denaturated (10 min, 80°C), ethidium bromide added and samples run on a 1% denaturizing agarose gel (MOPS). Quality and quantity of RNA was evaluated under UV light before RNA was blotted onto Hybond-N nylon membrane (Amersham) using standard upward blotting technique in 10 × SSC blotting buffer and crosslinked at 120 mJ/cm2 using a UV stratalinker (Stratagene). PCR probes were produced as previously described for microarray probes and analyzed on agarose gel. The individual probes were cut out of the gel and purified using DNA Gelextraction (Millipore). PCR product (25 ng) was labeled with 32P 6000 Ci/mmol (GE Healthcare) using StripEZ®DNA (Ambion) and cleaned with Quiaquick Nucleotide Removal Kit (Qiagen) according to manufactures instructions. After denaturation (90°C, 10 min) individual probe were hybridized to individual membranes at 68°C over night (Perfect Hyb 4M Plus hybridization buffer (Sigma), 5 ml per filter). The membranes were washed with 2 × SSC/0.1% SDS (2 × 5 min RT), 1 × SSC/0.1% SDS (1 × 15 min RT), 0.5 × SSC/0.1% SDS (2 × 10 min 68°C) followed by exposure on Kodak BioMax MS for 1 to 3 days. Membranes were stripped once (for re-use) with StripEZ®DNA (Ambion) according to manufactures instructions.
This project was financed by the Research Council of Norway (grants 136261/140, 153171/130) and Intervet International B.V. The project was supported by NFR's Salmon Genome Project and the functional genomics program FUGE, technology platform for microarrays. The Norwegian Microarray Platform provided the microarray printing service and also provided mouse clones for control probe production. We are grateful to Heidi Kongshaug at the Institute of Marine Research for excellent technical help in the laboratory and to Kjell Pettersen at the Computer Biology Unit for help with depositing our microarray data inn Array Express.
- Heuch PA, Bjorn PA, Finstad B, Holst JC, Asplin L, Nilsen F: A review of the Norwegian 'National Action Plan Against Salmon Lice on Salmonids': The effect on wild salmonids. Aquaculture. 2005, 246: 79-92. 10.1016/j.aquaculture.2004.12.027.View ArticleGoogle Scholar
- Krkosek M, Lewis MA, Morton A, Frazer LN, Volpe JP: Epizootics of wild fish induced by farm fish. Proc Natl Acad Sci USA. 2006, 103: 15506-15510. 10.1073/pnas.0603525103.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnson SC, Albright LJ: The developmental stages of Lepeophtheirus-salmonis (Kroyer, 1837) (Copepoda: Caligidae). Can J Zoo. 1991, 169: 929-950.View ArticleGoogle Scholar
- Schram TA: Supplementary description of the developmental stages of Lepeophtheirus salmonis (Kroyer, 1837) (Copepoda: Caligidae) on sea trout (Salmo trutta), of southern Norway. Pathogens of wild and farmed fish: sea lice. Edited by: Boxshall GA, Defaye D. 1993, Chichester, Ellis Horwood, 30-50.Google Scholar
- Brooks KM: The effects of water temperature, salinity, and currents on the survival and distribution of the infective copepodid stage of sea lice (Lepeophtheirus salmonis) originating on Atlantic salmon farms in the Broughton Archipelago of British Columbia, Canada. Rev Fish Sci. 2005, 13: 177-204. 10.1080/10641260500207109.View ArticleGoogle Scholar
- Gonzalez-Alanis P, Wright GM, Johnson SC, Burka JF: Frontal Filament Morphogenesis in the Salmon Louse Lepeophtheirus salmonis. J Parasitol. 2001, 87: 561-574.PubMedView ArticleGoogle Scholar
- Grimnes A, Jakobsen PJ: The physiological effects of salmon lice infection on post-smolt of Atlantic salmon. J Fish Biol. 1996, 48: 1179-1194. 10.1111/j.1095-8649.1996.tb01813.x.View ArticleGoogle Scholar
- Sconbaum CP, Lee S, Mahowald AP: The Drosophila yolkless gene encodes a vitellogenin receptor belonging to the low density lipoprotein receptor superfamily. Proc Natl Acad Sci USA. 1995, 92: 1485-1489. 10.1073/pnas.92.5.1485.View ArticleGoogle Scholar
- Ciudad L, Piulachs MD, Belles X: Systemic RNAi of the cockroach vitellogenin receptor results in a phenotype similar to that of the Drosophila yolkless mutant. FEBS. 2006, 273: 325-335. 10.1111/j.1742-4658.2005.05066.x.View ArticleGoogle Scholar
- Martin D, Wang SF, Raikhel AS: The vitellogenin gene of the mosquito Aedes aegypti is a direct target of ecdysteroid receptor. Mol Cell Endocrinol. 2001, 173: 75-86. 10.1016/S0303-7207(00)00413-5.PubMedView ArticleGoogle Scholar
- Chang ES, O'Connor JD: In vitro secretion and hydroxylation of alpha-ecdysone as a function of crustacean molt cycle. Gen Comp Endocrinol. 1978, 36: 151-160. 10.1016/0016-6480(78)90060-6.PubMedView ArticleGoogle Scholar
- Dubrovsky EB: Hormonal cross talk in insect development. Trends Endocrinol Metab. 2005, 16: 6-11. 10.1016/j.tem.2004.11.003.PubMedView ArticleGoogle Scholar
- Ritchie G, Mordue AJ, Pike AW, Rae GH: Morphology and ultrastructure of the reproductive system of Lepeophtheirus salmonis (Kroyer, 1837) (Copepoda: Caligidae). J Crust Biol. 1996, 16: 330-346. 10.2307/1548891.View ArticleGoogle Scholar
- Anstensrud M: Effects Of Mating On Female Behavior And Allometric Growth In The 2 Parasitic Copepods Lernaeocera-Branchialis (L 1767) (Pennellidae) And Lepeophtheirus-Pectoralis (Muller, 1776) (Caligidae). Crustaceana. 1990, 59: 245-258. 10.1163/156854090X00435.View ArticleGoogle Scholar
- Heuch PA, Nordhagen JR, Schram TA: Egg production in the salmon louse [Lepeophtheirus salmonis (Kroyer)] in relation to origin and water temperature. Aquaculture Research. 2000, 31: 805-814. 10.1046/j.1365-2109.2000.00512.x.View ArticleGoogle Scholar
- Tjensvoll K, Hodneland K, Nilsen F, Nylund A: Genetic characterization of the mitochondrial DNA from Lepeophtheirus salmonis (Crustacea; Copepoda). A new gene organization revealed. Gene. 2005, 353: 218-230. 10.1016/j.gene.2005.04.033.PubMedView ArticleGoogle Scholar
- wFleaBase, Daphnia Water Flea Genome Database. [http://wfleabase.org/genome-summaries/]
- Kvamme BO, Frost P, Nilsen F: The cloning and characterisation of full-length trypsins from the salmon louse Lepeophtheirus salmonis. Mol Biochem Parasitol. 2004, 136: 303-307. 10.1016/j.molbiopara.2003.11.017.PubMedView ArticleGoogle Scholar
- Kvamme BO, Skern R, Frost P, Nilsen F: Molecular characterisation of five trypsin-like peptidase transcripts from the salmon louse (Lepeophtheirus salmonis) intestine. Int J Parasitol. 2004, 34: 823-832. 10.1016/j.ijpara.2004.02.004.PubMedView ArticleGoogle Scholar
- Kvamme BO, Kongshaug H, Nilsen F: Organisation of trypsin genes in the salmon louse (Lepeophtheirus salmonis, Crustacea, copepoda) genome. Gene. 2005, 352: 63-74. 10.1016/j.gene.2005.03.011.PubMedView ArticleGoogle Scholar
- Toronen P, Kolehmainen M, Wong G, Castren E: Analysis of gene expression data using self-organizing maps. FEBS Lett. 1999, 451: 142-146. 10.1016/S0014-5793(99)00524-4.PubMedView ArticleGoogle Scholar
- Kohonen T: Self-organizing maps. Edited by: Kohonen T, Schroeder MR, Huang TS. 2001, New York, Springer-verlag, 3View ArticleGoogle Scholar
- Fellenberg K, Hauser NC, Brors B, Neutzner A, Hoheisel JD, Vingron M: Correspondence analysis applied to microarray data. Proc Natl Acad Sci USA. 2001, 98: 10781-10786. 10.1073/pnas.181597298.PubMedPubMed CentralView ArticleGoogle Scholar
- Okazaki Y, Shizuri Y: Effect of inducers and inhibitors on the expression of bcs genes involved in cypris larval attachment and metamorphosis of the barnacles Balanus amphitrite. Int J Dev Bio. 2000, 44: 451-456.Google Scholar
- Elvin CM, Vuocolo T, Pearson RD, East IJ, Riding GA, Eisemann CH, Tellam RL: Characterization of a major peritrophic membrane protein, peritrophin-44, from the larvae of Lucilia cuprina cDNA and deduced amino acid sequences. J Biol Chem. 1996, 271: 8925-8935. 10.1074/jbc.271.15.8925.PubMedView ArticleGoogle Scholar
- Casu R, Eisemann C, Pearson R, Riding G, East I, Donaldson A, Cadogan L, Tellam R: Antibody-mediated inhibition of the growth of larvae from an insect causing cutaneous myiasis in a mammalian host. Proc Natl Acad Sci USA. 1997, 94: 8939-8944. 10.1073/pnas.94.17.8939.PubMedPubMed CentralView ArticleGoogle Scholar
- Shen Z, Jacobs-Lorena M: A type I peritrophic matrix protein from the malaria vector Anopheles gambiae binds to chitin. Cloning, expression, and characterization. J Biol Chem. 1998, 273: 17665-17670. 10.1074/jbc.273.28.17665.PubMedView ArticleGoogle Scholar
- Trewitt PM, Heilmann LJ, Degrugillier SS, Kumaran AK: The boll weevil vitellogenin gene: nucleotide sequence, structure, and evolutionary relationship to nematode and vertebrate vitellogenin genes. J Mol Evol. 1992, 34: 478-492. 10.1007/BF00160462.PubMedView ArticleGoogle Scholar
- Nose Y, Lee JM, Ueno T, Hatakeyama M, Oishi K: Cloning of cDNA for vitellogenin of the parasitoid wasp, Pimpla nipponica (Hymenoptera: Apocrita: Ichneumonidae): vitellogenin primary structure and evolutionary considerations. Insect Biochem Mol Biol. 1997, 27: 1047-1056. 10.1016/S0965-1748(97)00091-X.PubMedView ArticleGoogle Scholar
- Nakajima H, Matoba K, Matsumoto Y, Hongo T, Kiritaka K, Sugino H, Nagamatsu Y, Hamaguchi Y, Ikegami S: Molecular characterization of a novel nucleolar protein in starfish oocytes which is phosphorylated before and during oocyte maturation. Eur J Biochem. 2000, 267: 295-304. 10.1046/j.1432-1327.2000.00931.x.PubMedView ArticleGoogle Scholar
- Warner AH, Pullumbi E, Amons R, Liu L: Characterization of a cathepsin L-associated protein in Artemia and its relationship to the FAS-I family of cell adhesion proteins. Eur J Biochem. 2004, 271: 4014-4025. 10.1111/j.1432-1033.2004.04338.x.PubMedView ArticleGoogle Scholar
- Bastiani MJ, Harrelson AL, Snow PM, Goodman CS: Expression of fasciclin I and II glycoproteins on subsets of axon pathways during neuronal development in the grasshopper. Cell. 1987, 48: 745-55. 10.1016/0092-8674(87)90072-9.PubMedView ArticleGoogle Scholar
- Clout NJ, Hohenester E: A model of FAS1 domain 4 of the corneal protein beta(ig)-h3 gives a clearer view on corneal dystrophies. Mol Vis. 2003, 9: 440-448.PubMedGoogle Scholar
- Pike A, Wadsworth SL: Sealice in salmonids: Their biology and control. Adv Parasitol. 1999, 44: 234-318.Google Scholar
- Smith JA, Whitfield PJ: Ultrastructural studies on the early cuticular metamorphosis of adult female Lernaeocera branchialis (L.) (Copepoda, Pennellidae). Hydrobiologia. 1988, 167-168. 607–616Google Scholar
- Bron JE, Shinn AP, Sommerville C: Ultrastructure of the cuticle of chalimus larva of the salmon louse Lepeophtheirus salmonis (Krøyer, 1837) (Copepoda, Caligidae). Contributions to Zoology. 2000, 69: 39-49.Google Scholar
- Brandal PO, Egidius E, Romslo I: Host blood: a major food component for the parasitic copepod Lepeophtheirus salmonis Kroyeri, 1838 (Crustacea: Caligidae). Norw J Zool. 1976, 24: 341-343.Google Scholar
- Williamson AL, Brindley PJ, Knox DP, Hotez PJ, Loukas A: Digestive proteases of blood-feeding nematodes. Trends Parasitol. 2003, 19: 417-423. 10.1016/S1471-4922(03)00189-2.PubMedView ArticleGoogle Scholar
- Williamson AL, Lecchi P, Turk BE, Choe Y, Hotez PJ, McKerrow JH, Cantley LC, Sajid M, Craik CS, Loukas A: A multi-enzyme cascade of haemoglobin proteolysis in the intestine of blood-feeding hookworms. J Biol Chem. 2004, 279: 35950-35957. 10.1074/jbc.M405842200.PubMedView ArticleGoogle Scholar
- Rosenthal PJ, Nelson RG: Isolation and characterization of a cysteine proteinase gene of Plasmodium falciparum. Mol Biochem Parasitol. 1992, 51: 143-152. 10.1016/0166-6851(92)90209-3.PubMedView ArticleGoogle Scholar
- McAllister L, Goodman CS, Zinn K: Dynamic expression of the cell adhesion molecule fasciclin I during embryonic development in Drosophila. Development. 1992, 115: 267-76.PubMedGoogle Scholar
- Kidokoro K, Iwata K, Fujiwara Y, Takeda M: Effects of juvenile hormone analogs and 20-hydroxyecdysone on diapause termination in eggs of Locusta migratoria and Oxya yezoensis. J Insect Physiol. 2006, 52: 473-479. 10.1016/j.jinsphys.2006.01.001.PubMedView ArticleGoogle Scholar
- Raikhel AS, Kokoza VA, Zhu J, Martin D, Wang SF, Li C, Sun G, Ahmed A, Dittmer N, Attardo G: Molecular biology of mosquito vitellogenesis: from basic studies to genetic engineering of antipathogen immunity. Insect Biochem Mol Biol. 2002, 32: 1275-1286. 10.1016/S0965-1748(02)00090-5.PubMedView ArticleGoogle Scholar
- Saunders DS, Henrich VC, Gilbert LI: Induction of diapause in Drosophila melanogaster: photoperiodic regulation and the impact of arrhythmic clock mutations on time measurement. Proc Natl Acad Sci USA. 1989, 86: 3748-3752. 10.1073/pnas.86.10.3748.PubMedPubMed CentralView ArticleGoogle Scholar
- Subramoniam T: Crustacean ecdysteriods in reproduction and embryogenesis. Comp Biochem Physiol C Toxicol Pharmacol. 2000, 125: 135-156.PubMedGoogle Scholar
- Block DS, Bejarano AC, Chandler GT: Ecdysteroid concentrations through various life-stages of the meiobenthic harpacticoid copepod, Amphiascus tenuiremis and the benthic estuarine amphipod, Leptocheirus plumulosus. Gen Comp Endocrinol. 2003, 132: 151-160. 10.1016/S0016-6480(03)00062-5.PubMedView ArticleGoogle Scholar
- Beckstead RB, Lam G, Thummel CS: The genomic response to 20-hydroxyecdysone at the onset of Drosophila metamorphosis. Genome Biol. 2005, 6: R99-10.1186/gb-2005-6-12-r99.PubMedPubMed CentralView ArticleGoogle Scholar
- Dana AN, Hong YS, Kern MK, Hillenmeyer ME, Harker BW, Lobo NF, Hogan JR, Romans P, Collins FH: Gene expression patterns associated with blood-feeding in the malaria mosquito Anopheles gambiae. BMC Genomics. 2005, 6: 5-10.1186/1471-2164-6-5.PubMedPubMed CentralView ArticleGoogle Scholar
- Lemoine A, Mathelin J, Braquart-Varnier C, Everaerts C, Delachambre J: A functional analysis of ACP-20, an adult-specific cuticular protein gene from the beetle Tenebrio: role of an intronic sequence in transcriptional activation during the late metamorphic period. Insect Mol Biol. 2004, 13: 481-493. 10.1111/j.0962-1075.2004.00508.x.PubMedView ArticleGoogle Scholar
- King-Jones K, Charles JP, Lam G, Thummel CS: The ecdysone-induced DHR4 orphan nuclear receptor coordinates growth and maturation in Drosophila. Cell. 2005, 121: 773-784. 10.1016/j.cell.2005.03.030.PubMedView ArticleGoogle Scholar
- C. elegans Sequencing Consortium: Genome sequence of the nematode C. elegans: a platform for investigating biology. Science. 1998, 282: 2012-2018. 10.1126/science.282.5396.2012.View ArticleGoogle Scholar
- Schwarz EM: Genomic classification of protein-coding gene families. WormBook. 2005, The C. elegans research Community, WorkBook, doi/10.1895/wormbook1.7.1, [http://www.wormbook.org]Google Scholar
- BaNG, Nematode and Neglected Genomics. [http://www.nematodes.org/index.shtml]
- Frost P, Nilsen F: Validation of reference genes for transcription profiling in the salmon louse, Lepeophtheirus salmonis, by quantitative real-time PCR. Vet Parasitol. 2003, 118: 169-174. 10.1016/j.vetpar.2003.09.020.PubMedView ArticleGoogle Scholar
- Dysvik B, Jonassen I: J-Express: exploring gene expression data using Java. Bioinformatics. 2001, 17: 369-370. 10.1093/bioinformatics/17.4.369.PubMedView ArticleGoogle Scholar
- Yang YH, Dudoit S, Luu P, Speed TP: Normalization for cDNA microarray data. Proceedings from SPIE, Microarrays optical technologies and informatics. Edited by: Bittner ML, Chen Y, Dorsel AN, Dougherty ER. 2001, , 4266: 141-152.View ArticleGoogle Scholar
- Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D, Altman RB: Missing value estimation methods for DNA microarrays. Bioinformatics. 2001, 17: 520-525. 10.1093/bioinformatics/17.6.520.PubMedView ArticleGoogle Scholar
- ArrayExpress. [http://www.ebi.ac.uk/microarray-as/aer/?#ae-main]