Open access

Sex-biased gene expression and evolution of candidate reproductive transcripts in adult stages of salmon lice (Lepeophtheirus salmonis)

Publication: FACETS
27 June 2019

Abstract

The salmon louse Lepeophtheirus salmonis (Krøyer 1837) displays numerous sexually dimorphic characteristics. Insights into their underlying molecular components have only recently been explored, which serve to better understand both the basic biology of the louse, associated impacts on drug sensitivity, and evolution of resistance. Expression of 16 L. salmonis genes putatively involved in sexual dimorphism and reproduction were used to determine differences between sexes and better understand responses to mating using RT-qPCR of pre-adult and adult lice. Analysis of these genes revealed the dynamic nature of sex-biased expression across stages. However, female reception of a spermatophore did not appear to impact the expression of these particular genes. Furthermore six of these transcripts and 84 others identified previously in a large-scale louse transcriptomics experiment were used to estimate differences in evolutionary rates and codon-usage bias of sex-related genes using phylogenetic analysis by maximum likelihood (PAML) and codonW. Results suggest male-biased genes are evolving at significantly greater rates than female-biased and unbiased genes as evidenced by higher rates of non-synonymous substitutions and lower codon-usage bias in these genes. These analyses expand our understanding of interactions of sex-biased expression across the pre-adult and adult life stages and provide foundations for better understanding evolutionary differences in sex-biased genes of L. salmonis.

Introduction

Lepeophtheirus salmonis (Krøyer 1837), commonly referred to as salmon lice, are ectoparasitic copepods of the family Caligidae, which parasitize both wild and farmed salmonids causing major economic losses to the salmonid aquaculture industries. The Atlantic salmon (Salmo salar) sector, valued at ~$12 billion USD (FAO 2015), is most impacted by L. salmonis infections in the northern hemisphere. Costello (2009) reported a global industry loss of approximately $480 million USD from lice infections alone, whereas more recent models estimating the cost associated with treating salmon lice in Norway demonstrate that these costs continue to escalate (Liu and Bjelland 2014; Abolofia et al. 2017). Reduced treatment efficacy and increased management costs are often associated with the acquisition of resistance in lice populations to commonly used therapeutics such as organophosphates, pyrethroids, macrocyclic lactones, and hydrogen peroxide (Fallang et al. 2004; Igboeli et al. 2013; Whyte et al. 2013; Aaen et al. 2015; Sutherland et al. 2015). Development of novel alternative management strategies will be essential to combat the rise in resistant populations and maintain the sustainability of the industry.
Salmon lice are known to naturally occur on wild salmonids, attaching to the surface and feeding on the mucus, skin, and blood of the fish. In high infestation scenarios the risks for severe damage and host death greatly increase (Wagner et al. 2008). The life cycle of L. salmonis consists of eight developmental stages: nauplii I/II, copepodid, chalimus I/II, pre-adult I/II, and adult (Hamre et al. 2013). The first two stages are free-living/planktonic, the third (copepodid) are infectious, chalimus I/II are attached to the host and non-motile, and pre-adult/adults are attached to host and motile. Sexually dimorphic traits are first observed at the chalimus II stage (Eichner et al. 2014), whereas other sex differences in physiology and behavioral characteristics are observed at the pre-adult and adult stages. The rise of sexually dimorphic characteristics plays an important role in the differential development, behavior, host–parasite interaction, and evolution of male and female lice (Johnson and Albright 1991; Ritchie et al. 1996; Hull et al. 1998; Ingvarsdóttir et al. 2002; Stephenson 2012; Wotton et al. 2014). The size of the adult male cephalothorax and genital segment are approximately half the size of females or less, and the shape of genital segment varies between the sexes with that of males having an ovoid appearance and that of females having an ovoid center with prominent posterolateral lobes (Johnson and Albright 1991). Behavioral differences are also abundant in salmon lice, where males exhibit increased frequency of mobile mate-locating, host-switching, and mate-guarding behaviors in contrast with the female’s tendency to remain in preferred areas on the host surface to serially produce egg strands (Ritchie et al. 1996; Hull et al. 1998; Stephenson 2012). Host–parasite interactions are also influenced by sexual dimorphism where females tend to consume blood meals more frequently than males (Brandal et al. 1976), insinuating that they have potential to be the more damaging sex on the host. Sexual dimorphism in L. salmonis has noticeable impacts on the aquaculture industry and is best described by the reduced sensitivity of male lice compared with females (regardless of population resistance status) to emamectin benzoate (EMB; Igboeli et al. 2013; Whyte et al. 2013; Sutherland et al. 2015), a compound registered for used by nearly all Atlantic salmon-producing countries. However, with our current knowledge of sexual dimorphism in salmon lice it cannot be fully determined whether differences in resistance between sexes is driven solely by underlying molecular mechanisms, inherent physiological differences between sexes, or a combination of both.
Sexually dimorphic characters largely manifest through differences in gene expression resulting from natural selection, sexual selection, and (or) sexual antagonism whereby sex-biased or sex-specific genes are generated and further selected upon differently influencing the fitness of each sex (Ellegren and Parsch 2007; Parsch and Ellegren 2013). For the majority of cases, male-biased genes appear to have consistently faster rates of evolution than female-biased or unbiased genes, which can be measured by the ratio of non-synonymous to synonymous substitution rates (dN/dS) for protein-coding sequences between two related species (Yang and Nielsen 2000). These elevated rates of evolution are primarily seen in male reproductive tissues such as the accessory glands, seminal fluids, and testes, which are heavily influenced by sexual selection (i.e., male–male- and sperm-competition) (Swanson and Vacquier 2002; Ellegren and Parsch 2007; Haerty et al. 2007; Ram and Wolfner 2007; Avila et al. 2011; Parsch and Ellegren 2013; Ramm et al. 2014; Cassone et al. 2017). These strong sexual selection pressures on male-biased genes may lead to the fixation of alleles that are harmful to females, driving selection to minimize the associated negative effects (sexual antagonism). These types of sexually selective and antagonistic pressures can lead to a rise in the number of sex-biased genes, novel functioning genes, and increased evolutionary rates of male-biased and female-biased genes when compared with unbiased genes (Ellegren and Parsch 2007). The quantification of non-synonymous to synonymous substitutions (dN/dS), however, cannot readily discern which selective pressures may be acting on a given group of genes, or whether these genes are simply under less selective constraint and therefore passively accumulating non-synonymous substitutions. A recent paper by Poley et al. (2016) provided a large-scale identification of sex-biased transcripts involved in sexual dimorphism from a 38K microarray analysis for the pre-adult II stage of L. salmonis. The authors found a significant portion of the transcriptome (34.7%–42.7%; 1729–3303 transcripts) contained sex-biased expression, including genes with putative links to sexually dimorphic characteristics and reproductive functions. Additionally, a high-throughput sequencing of sex-specific Caligus rogercresseyi transcriptomes has been provided by Farlora et al. (2014) identifying comparable counts of sex-biased or sex-specific transcripts (1733–2206). Poley et al. (2016) also investigated differences in evolutionary rates between male-biased, female-biased, and unbiased genes by quantifying the proportions of orphans (a gene without a detectable homologue) in each category. Male-biased transcripts had the greatest proportion of orphan genes (~50%) compared with female-biased, and unbiased (~20% and ~32%, respectively) suggesting the increase in orphan male-biased genes could be the product of increased positive selection and development of novel functions.
The current study aimed to further investigate sex-biased gene expression and the differing influence of evolution, as well as investigate the enrichment of sex-biased reproductive genes by introducing a novel analysis for mating-responsive gene expression in salmon lice. The objectives of this study were to (1) confirm putative sex-biased genes previously identified in louse transcriptomic studies (Eichner et al. 2008; Poley et al. 2016) across the pre-adult II and adult life stages, (2) determine the impact of copulation/reproductive status on sex-biased gene expression of female lice, and (3) characterize evolutionary rates of sex-biased and unbiased genes of L. salmonis.

Methods

Collection/study design

Adult female L. salmonis were collected from Bay Management Area 2A (BMA2a), Bay of Fundy, New Brunswick, Canada, and returned to the Huntsman Marine Science Center in St. Andrews, New Brunswick. Egg strands were removed and hatched using a seawater hatching system at 10 °C as described by Igboeli et al. (2013). The resulting copepodid larvae were used to infect 27 Atlantic salmon (Salmo salar) smolts obtained from Cooke Aquaculture hatchery in Elmsville, New Brunswick, housed in 10 L re-circulating aquaria at 12 °C ± 0.5 °C and a salinity of 33 ± 2 ppt containing one fish per randomly selected tank. All fish were anesthetized individually using MS-222 (Sigma-Aldrich, Oakville, Ontario, Canada) at a concentration of 0.1 g/L before being placed in 10 L of fresh seawater containing L. salmonis copepodids (30 per fish) for 5 min. Resulting attached lice were allowed to develop to the pre-adult II stage at which point sex could be determined macroscopically. At 22 days post infection (dpi), 10 fish were anesthetized as before, and all pre-adult II male lice (n = 13) were removed to create a virgin female group (n = 12). Six pre-adult II females (n = 6) were also collected from the host during this sampling, flash-frozen, and stored alongside the males at −80 °C for downstream analyses. Once the lice were collected, the anesthetized fish recovered in a salt water bath and then were returned to their original tanks. Another 10 fish (i.e., those infected with both male and female lice) were allocated to the mated-female group (females, n = 15; males n = 11). The remaining seven fish were used to monitor louse development and male-to-female interactions (mate-guarding) to ensure optimal sampling points throughout the study, which was terminated at 36 dpi when all lice had developed to the adult stage. All individuals from the mated-females group had received a spermatophore; similarly, all individuals from the virgin group had an absence of a spermatophore. The presence of egg strands at various developmental stages was observed for the majority of females. Fish were euthanized using 0.25 g/L MS-222 and lice were collected and flash frozen prior to being stored at −80 °C as before.

RNA isolation and purification

Total RNA was extracted from each louse individually using TRIzol®, according to the manufacturers’ instructions (Thermo Fisher Scientific™, Burlington, Ontario, Canada; Chomczynski 1993). An in-solution DNase treatment, TURBO DNA-free DNA removal kit (Ambion™, Foster City, California, USA) was used following the manufacturer’s instructions to remove residual traces of DNA. All samples were analyzed for quantity and purity using a Nanodrop 2000 Spectrophotometer (Thermo Fisher Scientific™, Burlington, Ontario, Canada), whereas a subset of 24 samples was randomly selected for assessing RNA integrity using ExperionStdSens Analysis Kit™ (Bio-Rad Laboratories, Mississauga, Ontario, Canada) as per manufacturer’s instructions. One sample was removed from the study because of potential RNA degradation. In total, 57 RNA extractions were executed, creating a study design with stage, sex, and adult female fertilization status as explanatory variables.

RT-qPCR

DNase-treated RNA was synthesized to 2 μg of cDNA in an 40 μL solution using iScript™ Reverse Transcription Supermix for RT-qPCR (Bio-Rad) as per manufacturer’s instructions. An identical protocol was used on pooled RNA without the reverse-transcriptase enzyme (no-RT controls) to confirm an absence of genomic contamination. A pool of cDNA containing 3 μL from each sample was prepared for testing primer specificity and efficiency. Primers were designed to specifically target sex-biased sequences from Eichner et al. (2008), Poley et al. (2016), or GenBank sequences related to sexual dimorphism or reproduction (Table 1). A 1%–2% agarose gel electrophoresis procedure using either a TAE or TBE buffer was completed using PCR product to ensure primers amplified a single product from synthesized cDNA before RT-qPCR analysis. Transcript-specific primers, which did not produce a fluorescent band or had multiple bands, were removed from the study. Serial dilutions for 16 genes of interest and two reference genes, elongation factor 1α (EF1α) and ribosomal protein subunit 20 (RPS20), were generated using a 6-step 5-fold dilution series. Efficiency curves were performed in 96 well plates using an Eppendorf Realplex with a 95 °C polymerase activation step for 30 s followed by 40 cycles of 95 °C for denaturation and 60 °C for combined annealing and extension, each lasting 15 s. Three-hundred and eighty-four well plates were loaded using the Aurora Biomed VERSA 10 for RT-qPCR amplification using the Bio-Rad CFX384™ Real-Time System. SsoAdvanced SYBR Green Supermix (Bio-Rad) was used to perform 11 μL qPCR reactions with 1 μL template and a 0.1 μmol/L concentration of each primer pair following the same thermal regime described above. Melt curve analysis was performed by increasing the temperature in 0.5 °C increments every 5 s from 65 °C to 95 °C to corroborate single product formation from gel electrophoresis. Individual Cq values above 37 cycles were considered to be unreliable positives (i.e., inconsistent amplification across replicates; Ct range of –RTs, etc.) and were determined to be non-expressing for this analysis. Reference gene stability was measured using geNorm and mean normalization of RT-qPCR data was completed using qbase-PLUS (Biogazelle, Gent, Belgium) on a log2 scale relative to the geometric mean of EF1α and RPS20. Individual primer pair efficiencies were used when generating the calibrated normalized relative quantity (CNRQ) values. Gene expression data were visualized using histograms and quantile–quantile plots to confirm normal distribution within each group. In addition, a Levene’s test was used to determine homogeneity of the variances. A t test was performed to determine statistically significant differences between mated and virgin adult female lice groups, and a two-way analysis of variance (ANOVA) was performed with stage and sex as explanatory variables using a significance cut-off of p < 0.05. In cases where a significant result was obtained, Tukey’s HSD test p < 0.05 was used to determine the groups with different expression profiles.
Table 1.
Table 1. Primers for RT-qPCR using 16 genes of interest normalized to two reference genes (indicated by*).
GeneAccession No. SequenceEfficiencyProduct sizeSourcePutative function
elongation factor-1 α* FTTAAGGAAAAGGTCGACAGAC0.9777Frost and Nilsen (2003)RNA binding; transcription regulation
  RGCCGGCATCACCAGACTT    
ribosomal protein subunit 20* FGTCACCTCAACCTCCACTCC0.99274Frost and Nilsen (2003)RNA binding; cytoplasmic translation
  RTGACTTGCCTCAAAGTGAGC    
C-type kectinCDD: cd00037FAACGCTCGAAAGGTGTGTGA0.95134Poley et al. (2016)Reproduction; putative SFP; sperm storage
  RTTGGCTGAAAGATGACACCA    
peroxidaseQ01603FGGGCACTGGGACCTCTTTTG1.03113Poley et al. (2016)Putative SFP; chorion hardening
  RAACGGTCGTGAACTCAGCAA    
trypsin-like 4EF490878FTGGAGGTTGGGGTGTACTCA1.02150Eichner et al. (2008)Protease; secreted; feeding associated
  RTCTTTGTTACCCGTCCCAGC    
kunitz-like 1CDD: smart00131FAGGGGCTTCTTTCCATGCTC0.96180Poley et al. (2016)Protease inhibitor; putative SFP
  RTGGTCATCAGGCTCAATCTGG    
kunitz-like 2CDD: smart00131FTGAGCCCGTTACTGGAGAGA0.94187Eichner et al. (2008)Protease inhibitor; putative SFP
  RCACAGGAGCTTCCTCAGTGG    
carboxypeptidaseEF490891FCTTATCAAGGATCGTTTATGGAGGA0.98121Eichner et al. (2008)Protease; secreted; feeding associated
  RTATGATCTACTCCAGATTGAGCCC    
vitellogenin 1EF490954FTCAGGCTCTTCTACCCCTCC0.98200Eichner et al. (2008)Reproduction; egg yolk associated
  RCTGGAGAGGAATGCTGGTGG    
vitellogenin 2EF490955FCACCGTTGCCAAGGAAATCA1.03170Eichner et al. (2008)Reproduction; egg yolk associated
  RCTGGGGTTACAGAGCCAACA    
vitellogenin-like proteinEF490956FAGAGATGCTGTCCGTACCCT1.04163Eichner et al. (2008)Reproduction; egg yolk associated
  RTTTCGCAGGTAACTCGTGCT    
putative vitellogenin receptorP98163FAACTTGCTTCGCCCACATCA0.99166Poley et al. (2016)Reproduction; egg yolk associated
  RGGTGTAACAGAAGCCAAGCG    
nucleoplasmin-like proteinEF490859FGTGAAGGTGAGGCCGAAGAA1.06171Eichner et al. (2008)Reproduction; histone binding; sperm de-condensation during fertilization
  RTCACTGTTGTCGCAGTCCTC    
embryo cathepsin LEF490929FTTTCCTTGGGTTTGTGAGAGGG1.0098Eichner et al. (2008)Embryogenesis; cell adhesion molecules involved in egg production
  RGAGTAGGGCAATGCCAAAGTG    
odourant-binding proteinCDD: cl20701FACGGAACAGGAACCAATACAC0.98168Poley et al. (2016)Sensory; putative semiochemical associated protein
  RCAATAAGGGAAGCGTGTGAAGC    
chorion peroxidaseQ9VEG6FAGAGCAAAGGGCTGAATCCT0.99114Poley et al. (2016)Reproduction; putative SFP; chorion hardening of eggs
  RTGGACAATCAACCGGCTCTT    
akirin-2B1AXD8FCGCCATCGGACCAACAAAGC0.95148GenBankEmbryonic development; innate immune response
  RTTGCACATCCGCTCACAAATG    
protein sarahQ9XZL8FCTACGAGAGCCGATGAAGACG0.98102GenBankReproduction; meiotic progression in oocytes
  RTTGCTGTTTCCCACCCAACC    

Note: F, forward primer; RNA, ribonucleic acid; R, reverse primer; SFP, seminal fluid protein.

Selection of homologous genes and sequence alignment

To assess the evolutionary rate of L. salmonis sex-biased genes, putative full-length coding DNA sequences (FLcDNA) were obtained from a BLASTN search against L. salmonis sequences in either the NCBI Non-Redundant (NR) or Transcript Shotgun Assembly (TSA) databases using 90 contig sequences obtained from a previous microarray study on sex-biased gene expression (Poley et al. 2016). Balanced groups containing 30 FLcDNA sequences each were obtained for male-biased, female-biased, and unbiased genes selected based on largest sequence and an open reading frame (ORF) with ≥50% coverage where the transcript had to show concordant sex-bias expression across three populations of lice (Poley et al. 2016). Transcripts with ≥90% identity were selected for further analysis. An ORF finder (ncbi.nlm.nih.gov/orffinder/) was used to predict ORFs of selected transcripts. A TBLASTN using the predicted L. salmonis protein sequences was performed to obtain orthologs in C. rogercresseyi, the closest related organism to L. salmonis with sufficient genomic characterization to run this analysis. Accession numbers for L. salmonis and C. rogercresseyi orthologs can be found in Supplementary Material 2. In order for a pair of L. salmonis and C. rogercresseyi genes to be considered orthologous, results from a TBLASTN required the following criteria be met: E-value ≤10−10, similarity score ≥75%, coverage score ≥50%, and transcript length ≥300 base pairs (Leong et al. 2010). Protein sequences of orthologous pairs were aligned using EMBOSS WATER with default input settings and MARKx10 for FASTA output. The result was used as a template for manual trimming of the nucleotide sequences to match their respective protein alignment. RevTrans (1.4 server) was used to create codon-aware DNA alignments using default parameters, and the FASTA format output was selected (Wernersson and Pedersen 2003). Additionally RevTrans served as a quality assurance check for the manual nucleotide trimming as nucleotide sequences were required to be an exact match to the unaltered protein alignment supplied to RevTrans that acted as a scaffold for the nucleotide alignment.

Non-synonymous/synonymous estimation and codon usage bias

Alignments from the RevTrans output were used to analyze non-synonymous and synonymous mutation ratios (dN/dS) using the YN00 program of the phylogenetic analysis by maximum likelihood (PAML) software (Yang and Nielsen 2000; Yang 2007). Estimated evolutionary rates (dN/dS) were calculated for the 30 male-biased, female-biased, and unbiased genes from pairwise comparisons of L. salmonis and C. rogercresseyi. The YN00 program was run using the default settings with following changes; verbose: one detailed output, noisy: ≥9, and weighting—yes. Codon usage bias data for L. salmonis sequences was calculated using the program codonW (version 1.4.2 available from: codonw.sourceforge.net/) with default settings for frequency of optimal codons (Fop). Both dN/dS and codon usage bias data were visualized using histograms and quantile–quantile plots and a Levene’s test for homogeneity of variance. One-way ANOVAs were completed for both the evolutionary rate (frequency of dN/dS, dN, and dS) and codon usage bias (frequency of optimal codons) with male-biased, female-biased, and unbiased genes as explanatory variables using a significance cut-off of p < 0.05. In cases where a significant result was obtained, Tukey’s HSD test (p < 0.05) was used.

Ethical statement

All animals were cared for in accordance with the Guide to the Care and Use of Experimental Animals available at: ccac.ca/Documents/Standards/Guidelines/Experimental_Animals_Vol1.pdf. Animal care and use protocols were submitted for review and approval to the UPEI Animal Care Committee (UPEI Animal Care Protocol No. 10-014).

Results and discussion

Sex-biased gene expression

The L. salmonis transcriptome at the pre-adult II stage has been revealed to contain extensive sex-biased gene expression (ca. 35%–43%) (Poley et al. 2016). However, further characterization is required to understand the stability of sex-biased expression throughout different stages. Sex-biased genes were characterized for pre-adult II and adult L. salmonis from BMA2a in the Bay of Fundy, New Brunswick using RT-qPCR for differential gene expression of sex and stage. A total of 16 candidate transcripts (Table 1) were selected from previous salmon louse transcriptomic studies based on their high degree of sex-biased expression, putative involvement in sexually dimorphic properties and were hypothesized to have an influence on mating responsive gene expression. Putative functions include: embryogenesis, embryonic development, vitellogenesis, sperm storage, chorion hardening, feeding response, and semiochemical/pheromone detection (Eichner et al. 2008; Poley et al. 2016). Many of these genes had similar putative functions to known mating responsive genes from other arthropods with homologous sequences and, therefore, we attempted to determine whether this trend was general or taxon-specific (Mack et al. 2006; Rogers et al. 2008).
A principal component analysis (PCA) was performed to visualize the relationships between sex and stage using log2 CNRQ data for the 16 candidate genes (Fig. 1). The first and second principal components (PCs) explained 36.4% and 24.2% of the total variation in gene expression profiles, respectively. Both sex and stage were separated along PC1 (x-axis), where adult males and adult females were the most distant groups. The majority of genes assayed here were male-biased or female-biased and, therefore, contributed to the separation of samples along PC1. Genes influenced by stage only, unaffected by either sex or stage, and three male-biased genes odourant-binding protein, protein sarah, and vitellogenin-like (Table 1), controlled differences associated with PC2. Some pre-adult II female samples overlapped with male samples suggesting pre-adult II female expression of some genes is more similar to that of males, especially at the corresponding pre-adult II stage. Male samples also had less overall separation between stages compared with females along PC1, indicating stage differences in male expression was smaller than in females.
Fig. 1.
Fig. 1. Principal component analysis of pre-adult and adult stages of both female and male L. salmonis. Individual samples of pre-adult males (PM) are represented with red crosses, adult males (AM) with yellow triangles, pre-adult females (PF) with grey squares and lastly adult females (AF) with blue circles. Males and females are separated along PC1 (x-axis) and gene vectors represented with black arrows indicate association between sex and expression.
Four transcripts influenced only by sex included nucleoplasmin-like protein, vitellogenin-like protein, putative odourant-binding protein, and protein sarah (Fig. 2; Table 1). Nucleoplasmin-like protein displayed female-biased expression, whereas the remaining three genes were expressed highest in males. Sex-biased expression, however, was not always observed to be stable throughout the pre-adult II to adult life stages as displayed by the expression profiles of embryo cathepsin L, trypsin-like 4, vitellogenin 1 (LsVit1) and vitellogenin 2 (LsVit2) (Table 1; Fig. 3). All four transcripts were significantly upregulated at the adult female stage when compared with other sex and stage combinations, with embryo cathepsin L, vitellogenin 1, and vitellogenin 2 being expressed significantly lower in pre-adult II females when compared with pre-adult II males (p < 0.001; LsVit1 had no detectable expression (n.d.)) and adult males (p < 0.001) (Fig. 3). Conversely, Trypsin-4 was expressed highest in adult females followed by pre-adult II females and then males (p = 0.005) demonstrating a conservation of sex-biased expression that is further amplified across stage. Vitellogenins are precursors of egg yolk glycoproteins produced within fat bodies for many arthropods and have vital roles in progeny fitness (Schneider 1996; Eichner et al. 2008; Dalvin et al. 2009; Dalvin et al. 2011). In L. salmonis, vitellogenins (LsVit1 and LsVit2) are expressed throughout the adult female subcuticular tissue and are involved in embryogenesis. Additionally, these proteins provide an essential energy source throughout the planktonic larval life stages (Eichner et al. 2008; Dalvin et al. 2011; Edvardsen et al. 2014). Expression knock-down of a separate novel female specific L. salmonis yolk-associated protein (LsYAP) has been shown to disrupt proper embryonic development, highlighting their importance to progeny development and fitness (Dalvin et al. 2009). These results are consistent with a study by Eichner et al. (2008) who observed a gradual increase in expression of LsVit1 and LsVit2 from the PAII to adult females as well as six time points across adult stage development.
Fig. 2.
Fig. 2. L. salmonis genes with a main effect of sex determined by a one-way analysis of variance (ANOVA) with a p < 0.05 (*). A log2 calibrated normalized relative quantities (CNRQ) scale was used to represent expression differences between males and females. The circles represent data points that are outside 1.5 × the IQR of their respective boxplot.
Fig. 3.
Fig. 3. L. salmonis genes with an interaction effect of sex and stage determined by a two-way analysis of variance (ANOVA) and follow-up Tukey’s test with a p < 0.05. A log2 calibrated normalized relative quantities (CNRQ) scale was used to represent expression differences between males and females across pre-adult II and adult stages of L. salmonis. For each gene, box plots that do not share a lowercase letter are significantly different from each other. n.d. refers to transcripts which were not detected for a given sex and stage combination.
A total of six transcripts had no detectable expression in ≥2/3 of individuals from one sex across both stages. These transcripts (C-type lectin, chorion peroxidase, kunitz 1, and peroxidase) were, therefore, considered to be male-specific, whereas carboxypeptidase and vitellogenin receptor (Table 1; Table 2) were similarly classified as female-specific (Table 2). Although some individuals from the putatively non-expressing sex formed PCR products, all Cq values were high (between 34 and 39) and (or) melt curves were unreliable potentially representing artifacts rather than true expression from the opposite sex. Carboxypeptidase and vitellogenin receptor were classified as female-specific in expression and may be involved in two essential processes for female lice. Vitellogenin receptors are another egg-yolk-associated protein that facilitates the endocytosis of vitellogenins into developing oocytes (Schneider 1996; Tufail and Takeda 2009). As mentioned, successful incorporation of vitellogenins are essential to proper development of embryos (Dalvin et al. 2009) and disruption of this vital pathway may elucidate novel reproduction-associated louse management strategies. Carboxypeptidases are putatively involved in feeding and managing blood meals within L. salmonis similar to other insects and parasitic arthropods (Edwards et al. 1997; Warr et al. 2007; Isoe et al. 2009). Carboxypeptidase B proteins have been investigated in Anopheles gambiae as targets for development of a transmission-blocking vaccine for Plasmodium falciparum (Lavazec et al. 2007). An interesting secondary effect reported by Lavazec et al. (2007) was mosquitos having fed on anti-carboxypeptidase mice exhibited reduced reproductive capacity. Furthermore, a recombinant vaccine against a gut-associated glycoprotein BM86 in Boophilus microplus has been shown to be effective in controlling tick infestations (Willadsen et al. 1995). Female-specific carboxypeptidases and other gut/digestion associated proteins may provide interesting targets for further investigation to determine their viability as vaccination options in L. salmonis.
Table 2.
Table 2. Genes determined to have sex-specific expression.
GeneProportion of expressing malesProportion of expressing femalesSex-specificity
C-type lectin91.7% (22/24)23.8% (5/21)Male
peroxidase95.8% (23/24)33.3% (7/21)Male
chorion peroxidase95.8% (23/24)19.0% (4/21)Male
kunitz-like 191.7% (22/24)14.3% (3/21)Male
carboxypeptidase16.7% (4/24)71.4% (15/21)Female
putative vitellogenin receptor4.2% (1/24a)71.4% (15/21b)Female
a
No adult males had detectable expression.
b
Only one pre-adult II female had detectable expression, potential adult female-specific expression.
C-type lectin, chorion peroxidase, kunitz-1, and peroxidase are putative seminal fluid proteins (SFPs) in salmon lice providing further evidence for male-specific expression profiles. SFPs are the products of male accessory glands and ejaculatory ducts which are transferred along with sperm in the male ejaculate to a female recipient (Avila et al. 2011). SFPs have been shown to impact multiple processes such as egg production, female longevity, female behaviour, sperm storage, mating plug formation, and receptivity to re-mating across a wide range of arthropod species (Chapman 2001; Avila et al. 2011). Recently, a relevant new gene family has been identified in L. salmonis essential for proper development of the wall surrounding the spermatophore (Borchel and Nilsen 2018). Understanding the evolution of SFPs and how they elicit responses in females may provide practical implications for pest management strategies. Increased selection, or genetic engineering, of males producing SFPs that have greater association with negative impacts on female fitness such as reduced receptivity to re-mating or male sterilization may provide attractive avenues for alternative management strategies (Sirot et al. 2015). These results provided confirmation of the stability of sex-biased gene expression from the pre-adult II to adult stage of L. salmonis, but also determined transcripts may have differing expression profiles dependent on stage, similar to observations by Eichner et al. (2018) in earlier life stages (chalimus II and pre-adult I). The results provide further characterization of sexually dimorphic characteristics of L. salmonis and provide novel insights into sex-biased expression during the final life stages.

Mating-responsive gene expression

Another goal, while exploring the stability of sex-biased expression across stage, was to investigate the impact copulation/reception of a spermatophore had on female transcripts. Mated and virgin females were analyzed for differential expression for all 16 candidate transcripts, but no significant expression differences were observed between copulatory statuses (t test; p > 0.05 for all transcripts). The transfer of male SFPs are known to induce gene expression differences in Drosophila melanogaster and A. gambiae females (Soller et al. 1997; Peng et al. 2005; Mack et al. 2006; Ram and Wolfner 2007; Rogers et al. 2008; Avila et al. 2011). Reception of a spermatophore in D. melanogaster facilitates differential gene expression either transiently (1–3/4–6 h) or for longer periods (≥12 h) (McGraw et al. 2008). Candidate L. salmonis genes with similar functions to known short- and long-term mating responsive genes such as chorion protein 36/38, odourant binding protein 99a, yolk protein 1/2/3 and vitelline membrane 26Aa/32E/34Ca appear not to be impacted post-reception of a spermatophore in L. salmonis (McGraw et al. 2008). These results suggest, similar to a study by Eichner et al. (2008) that copulatory status does not appear to be a trigger for egg production despite having similar functions to genes that are sensitive to copulatory status in other arthropods.
Potential explanations for the lack of differential expression in mated female L. salmonis include, but are not limited to: (1) select candidate genes with similar functions to known mating-responsive genes within Arthropoda are potentially not mating-responsive in L. salmonis (Soller et al. 1997; McGraw et al. 2008; Rogers et al. 2008). (2) Mating-responsive gene expression is known to be transient within Arthropoda (McGraw et al. 2008; Rogers et al. 2008). Because of the limitation of lice numbers post infection, in addition to the difficulty monitoring copulatory status for each mating-pair, inclusion of multiple sampling time points post-mating wasn’t feasible. It is therefore possible the single sampling time point was either too early or late to capture differential expression. (3) Despite having mated and virgin lice isolated in separate tanks, it is possible the recirculation system was unable to remove semiochemicals produced by males from within the mated group tanks. If semiochemical cues are important to louse receptivity to mating then the virgin lice may have falsely perceived the presence of a male and, therefore, potentially influenced their reproductive mechanisms. (4) Lastly, the use of whole homogenates over specific tissues, more specifically the genital segment, may have resulted in a dilution effect of smaller scale expression differences as the greatest number of differentially expressed genes post-mating understandably appear in the reproductive tissues (Mack et al. 2006; Edvardsen et al. 2014).
Continued investigation into mating-responsive gene expression may provide insights into alternative management strategies that are becoming increasingly necessary to be effective and sustainable when managing salmon lice infestations. Reproduction targeted or disrupting strategies are not novel ideas and have been studied for decades (Cardé and Minks 1995). Synthetic semiochemicals provide an avenue to prevent or reduce mating success as a means of control instead of drug treatments susceptible to the acquisition of resistance (Smart et al. 2014). Furthermore, hormone disruption has been proposed as a viable alternative strategy to limit Anopheles spp. reproductive success and, thus, the spread of Plasmodium faliciparum (Childs et al. 2016). A recent study used the CRISPR-Cas9 system to distort the sex-ratio of A. gambiae by utilizing a specific ribosomal sequence targeted Cas9 approach during spermatogenesis to produce an extreme male-bias in the progeny (Galizi et al. 2016). Further efforts into the identification of mating-responsive gene expression, especially through means of whole transcriptome analyses of mated vs. virgin females will be beneficial to markers of insemination success, progeny development, and provide similar targets as described above for alternative management strategies.

Evolutionary rates of sex-biased genes and codon usage bias

The ratios of non-synonymous and synonymous substitutions (dN/dS) were calculated for 30 male-biased, female-biased, and unbiased transcripts identified in a previous transcriptomic study (Poley et al. 2016). Male-biased L. salmonis transcripts had significantly higher dN/dS values when compared with unbiased transcripts, whereas values for female-biased genes were between those of the two groups (Table 3; Fig. 4). This suggests that L. salmonis sex-biased genes, especially those that are male-biased, are under reduced selective constraint or greater positive selection (Yang and Nielsen 2000; Zhang et al. 2004; Ellegren and Parsch 2007). Analyses using dN/dS require a closely related species for sequence comparisons and C. rogercresseyi was selected as it is the closest related species to L. salmonis with sufficient genomic and transcriptomic information. It is important to note, however, that these species have been diverging for ~100–150 million years (Eyun 2017), which has allowed for greater accumulation of synonymous mutations over time, thereby reducing the overall sensitivity of the dN/dS analysis. When examining the dN and dS variables separately, it was observed that dS values for male-biased, female-biased, and unbiased genes were not significantly different and had similar saturation of synonymous mutations (dS). This saturation caused resulting dN/dS values to be underestimated because of the stable and consistently high dS. Similar to a study by Zhang et al. (2004) who also describe a saturation of dS even within Drosophila spp., we can consider only dN as it is the driving force on the overall dN/dS differences. Male-biased transcripts had significantly greater dN values when compared with both female-biased and unbiased transcripts (Table 3; Fig. 4). This result supports the previous analysis by Poley et al. (2016), who showed a greater proportion of orphan genes for those with male-biased expression and is consistent with results found in Drosophila spp. (Zhang et al. 2004).
Table 3.
Table 3. Rates of nucleotide substitution and frequency of optimal codon usage in sex-biased and unbiased L. salmonis genes.
 Male-biasedFemale-biasedUnbiasedPMUbPMFbPFUb
dNa0.25280.17720.13060.0003*0.0332*0.3018
dSa3.66013.38873.80550.61770.69770.9910
dN/dSa0.07060.05820.03310.0018*0.17900.1907
Fopa0.35400.40000.37300.06040.0002*0.1332

Note: dN, non-synonymous substitution rate; dS, synonymous substitution rate; Fop, frequency of optimal codon usage.

a
Each category is represented by the median value.
b
P value of a Tukey’s honest significant difference from an analysis of variance (ANOVA) for male-biased (M), female-biased (F), and unbiased (U) genes.*P < 0.05.
Fig. 4.
Fig. 4. Distributions of dN/dS, dN, and dS for male-biased (MB), female-biased (FB), and unbiased (UB) genes (n = 30/category) determined by nucleotide substitutions in protein coding sequences of L. salmonis and C. rogercresseyi.
Codon usage bias was also investigated for these gene categories to validate the dN/dS results. Codon usage bias can be described as the weak selection pressure that drives the accumulation of a specific codon over other synonymous codons (Akashi 1994; Hambuch and Parsch 2005; Behura and Severson 2013). Often these optimal codons reflect an organism’s tRNA pool and provide the benefits of increased translation efficiency and accuracy (Duret 2000). Codon usage bias is often measured by the frequency of optimal codons (Fop) where highly expressed and (or) heavily conserved genes have the greatest accumulation of optimal codons (Powell and Moriyama 1997). It was therefore expected that this weak selection pressure would have a lesser effect on L. salmonis genes with higher dN/dS values. Male-biased L. salmonis transcripts indeed had the lowest frequency of optimal codons in this analysis (Table 3), which is again consistent with results from D. melanogaster (Hambuch and Parsch 2005) and dN/dS values presented above.
Overall, these analyses indicate that male-biased L. salmonis transcripts likely follow the commonly observed trend of rapid evolution when compared with unbiased genes (Ellegren and Parsch 2007). Whether this is caused by a reduction in selective constraint or increased sexual, natural, and (or) antagonistic selection pressures is not yet fully understood (Zhang et al. 2004; Ellegren and Parsch 2007; Parsch and Ellegren 2013). The annotations of male-biased transcripts in L. salmonis are enriched for putative reproductive functions (Poley et al. 2016; Table 1) that are likely influenced by strong selective pressures, which is consistent with observed patterns of rapid evolution for reproduction associated genes (Swanson and Vacquier 2002; Haerty et al. 2007; Meisel 2011; Grath and Parsch 2012). Recently an interesting sex-associated single nucleotide polymorphism (SNP) marker was identified in salmon lice suggesting they operate using a ZW sex determining system (Carmichael et al. 2013). Global dosage compensation is a major factor in sex-biased gene expression as the homogametic sex of a given species has twice the number of available gene copies as the opposite sex for those genes present on sex chromosomes, for example D. melanogaster compensate by increasing the transcription of the single male X chromosome two-fold (Baker et al. 1994). Compensation mechanisms appear to be less effective (or selective for only a proportion of genes on the sex chromosome) in ZW systems than XY (Itoh et al. 2007; Mank 2009; Wolf and Bryk 2011). Determining the nature of L. salmonis dosage compensation will help elucidate the necessity and (or) consequences of sex-linked gene expression.
Our analyses have provided a useful starting point to better understand evolutionary aspects of male-biased, female-biased, and unbiased genes; however, they do lack some degree of specificity as transcripts were selected based on sex-biased expression and transcript length, not function. Future analyses targeting sex-biased expression with specific roles and expression breadths (i.e., reproductive genes in testes and ovaries) or sex-linked genes will provide greater insight into these mechanisms. Additionally, optimization and automation of evolutionary rate analyses for sex-specific whole transcriptome data would be useful to globally confirm elevated sex-biased gene evolution in L. salmonis. To our knowledge this is the first time these techniques have been applied to salmon lice and will serve as valuable tools for future research. These techniques can be applied to any gene groups for comparative purposes; for example, those involved in the host-parasite interaction, which are likely under similar elevated selection pressures when compared with somatic genes.

Conclusion

These analyses provide clarification and help our understanding of the existence and importance of sex-biased expression across life stage and differential evolutionary pressures on male-biased, female-biased, and unbiased genes in L. salmonis. As with any commercially important pest, understanding key biological differences involving life stage and sex (i.e., the impacts that elevated evolutionary rates have on virulence, reproductive capacity, or drug resistance) may lead to novel targeted approaches for intervention strategies and management.

Acknowledgements

This work was funded by NSERC Discovery (402288-12); ACOA-AIF TREAT2 (199308); Canada Excellence Research Chair in Aquatic Epidemiology program at UPEI (214602); Innovation PEI—Development and Commercialization grant DCFG (210205-70). The authors would like to thank the staff at Huntsman Marine Science Centre (Chris Bridger), and Cooke Aquaculture (Leighanne Hawkins, Keng Pee Ang) for aiding in site access, lice collection and culture. In addition the authors thank Jong Leong for aiding in the refinement of the methodology for sequence comparison and evolutionary analysis. Also, thanks to the aquatics staff at the Atlantic Veterinary College for help with fish husbandry and lice infections.

References

Aaen SM, Helgesen KO, Bakke MJ, Kaur K, and Horsberg TE. 2015. Drug resistance in sea lice: a threat to salmonid aquaculture. Trends in Parasitology, 31(2): 72–81.
Abolofia J, Asche F, and Wilen JE. 2017. The cost of lice: quantifying the impacts of parasitic sea lice on farmed salmon. Marine Resource Economics, 32(3): 329–349.
Akashi H. 1994. Synonymous codon usage in Drosophila melanogaster: natural selection and translational accuracy. Genetics, 136(3): 927–935.
Avila FW, Sirot LK, LaFlamme BA, Rubinstein CD, and Wolfner MF. 2011. Insect seminal fluid proteins: identification and function. Annual Review of Entomology, 56: 21–40.
Baker BS, Gorman M, and Marin I. 1994. Dosage compensation in Drosophila. Annual Review of Genetics, 28: 491–521.
Behura SK, and Severson DW. 2013. Codon usage bias: causative factors, quantification methods and genome-wide patterns: with emphasis on insect genomes. Biological Reviews, 88(1): 49–61.
Borchel A, and Nilsen F. 2018. A novel gene-family involved in spermatophore generation in the economically important salmon louse Lepeophtheirus salmonis. Molecular Reproduction and Development, 85: 478–489.
Brandal P, Egidius E, and Romslo I. 1976. Host blood: a major food component for the parasitic copepod, Lepeophtheirus salmonis Kroyer, 1838 (Crustacea caligidae). Norwegian Journal of Zoology, 24: 341–343.
Cardé RT, and Minks AK. 1995. Control of moth pests by mating disruption: successes and constraints. Annual Review of Entomology, 40(1): 559–585.
Carmichael SN, Bekaert M, Taggart JB, Christie HRL, Bassett DI., Bron JE, et al. 2013. Identification of a sex-linked SNP marker in the salmon louse (Lepeophtheirus salmonis) using RAD sequencing. PLoS ONE, 8(10): 1–8.
Cassone BJ, Kay RGG, Daugherty MP, and White BJ. 2017. Comparative transcriptomics of malaria mosquito testes: function, evolution, and linkage. G3 (Bethesda, Md.), 7(4): 1127–1136.
Chapman T. 2001. Seminal fluid-mediated fitness traits in Drosophila. Heredity, 87: 511–521.
Childs LM, Cai FY, Kakani EG, Mitchell SN, Paton D, Gabrieli P, et al. 2016. Disrupting mosquito reproduction and parasite development for malaria control. PLoS Pathogens, 12(12): 1–20.
Chomczynski P. 1993. A reagent for the single-step simultaneous isolation of RNA, DNA and proteins from cell and tissue samples. BioTechniques, 15: 532–537.
Costello MJ. 2009. The global economic cost of sea lice to the salmonid farming industry. Journal of Fish Diseases, 32: 115–118.
Dalvin S, Frost P, Biering E, Hamre LA, Eichner C, Krossøy B, et al. 2009. Functional characterisation of the maternal yolk-associated protein (LsYAP) utilising systemic RNA interference in the salmon louse (Lepeophtheirus salmonis) (Crustacea: Copepoda). International Journal for Parasitology, 39(13): 1407–1415.
Dalvin S, Frost P, Loeffen P, Skern-Mauritzen R, Baban J, Rønnestad I., et al. 2011. Characterisation of two vitellogenins in the salmon louse Lepeophtheirus salmonis: molecular, functional and evolutional analysis. Diseases of Aquatic Organisms, 94(3): 211–224.
Duret L. 2000. tRNA gene number and codon usage in the C. elegans genome are co-adapted for optimal translation of highly expressed genes. Trends in Genetics, 16(7): 287–289.
Edvardsen RB, Dalvin S, Furmanek T, Malde K, Mæhle S, Kvamme BO, et al. 2014. Gene expression in five salmon louse (Lepeophtheirus salmonis, Krøyer 1837) tissues. Marine Genomics, 18: 39–44.
Edwards MJ, Lemos FJA, Donnelly-Doman M, and Jacobs-Lorena M. 1997. rapid induction by a blood meal of a carboxypeptidase gene in the gut of the mosquito anopheles gambiae. Insect Biochemistry and Molecular Biology, 27(12): 1063–1072.
Eichner C, Frost P, Dysvik B, Jonassen I, Kristiansen B, and Nilsen F. 2008. Salmon louse (Lepeophtheirus salmonis) transcriptomes during post molting maturation and egg production, revealed using EST-sequencing and microarray analysis. BMC Genomics, 9(1): 126.
Eichner C, Hamre LA, and Nilsen F. 2014. Instar growth and molt increments in Lepeophtheirus salmonis (Copepoda: Caligidae) chalimus larvae. Parasitology International, 64(1): 86–96.
Eichner C, Dondrup M, and Nilsen F. 2018. RNA sequencing reveals distinct gene expression patterns during the development of parasitic larval stages of the salmon louse (Lepeophtheirus salmonis). Journal of Fish Diseases, 41: 1005–1029.
Ellegren H, and Parsch J. 2007. The evolution of sex-biased genes and sex-biased gene expression. Nature Reviews Genetics, 8(9): 689–698.
Eyun S. 2017. Phylogenomic analysis of Copepoda (Arthropoda, Crustacea) reveals unexpected similarities with earlier proposed morphological phylogenies. BMC Evolutionary Biology, 17(1): 1–12.
Fallang A, Ramsay JM, Sevatdal S, Burka JF, Jewess P, Hammell KL, et al. 2004. Evidence for occurrence of an organophosphate-resistant type of acetylcholinesterase in strains of sea lice (Lepeophtheirus salmonis Kroyer). Pest Management Science, 60(12): 1163–1170.
FAO. 2015. FAO yearbook. Fishery and aquaculture statistics. Food and Agriculture Organization of The United Nations, Rome, Italy [online]: Available from fao.org/3/a-i7989t.pdf.
Farlora R, Araya-Garay J, and Gallardo-Escárate C. 2014. Discovery of sex-related genes through high-throughput transcriptome sequencing from the salmon louse Caligus rogercresseyi. Marine Genomics, 15: 85–93.
Frost P, and Nilsen F. 2003. Validation of reference genes for transcription profiling in the salmon louse, Lepeophtheirus salmonis, by quantitative real-time PCR. Veterinary Parasitology, 118: 169–174.
Galizi R, Hammond A, Kyrou K, Taxiarchi C, Bernardini F, O’Loughlin SM, et al. 2016. A CRISPR-Cas9 sex-ratio distortion system for genetic control. Scientific Reports, 6: 2–6.
Grath S, and Parsch J. 2012. Rate of amino acid substitution is influenced by the degree and conservation of male-biased transcription over 50 myr of Drosophila evolution. Genome Biology and Evolution, 4: 346–359.
Haerty W, Jagadeeshan S, Kulathinal RJ, Wong A, Ram KR, Sirot LK, et al. 2007. Evolution in the fast lane: rapidly evolving sex-related genes in Drosophila. Genetics, 177: 1321–1335.
Hambuch TM, and Parsch J. 2005. Patterns of synonymous codon usage in Drosophila melanogaster genes with sex-biased expression. Genetics, 170(4): 1691–1700.
Hamre LA, Eichner C, Caipang CMA, Dalvin ST, Bron JE, Nilsen F, et al. 2013. The Salmon Louse Lepeophtheirus salmonis (Copepoda: Caligidae) life cycle has only two chalimus stages. PLoS One, 8(9): e73539.
Hull MQ, Pike AW, Mordue AJ, and Rae GH. 1998. Patterns of pair formation and mating in an ectoparasitic caligid copepod Lepeophtheirus salmonis (Kroyer 1837): implications for its sensory and mating biology. Philosophical Transactions of the Royal Society B, Biological Sciences, 353(1369): 753–764.
Igboeli OO, Burka JF, and Fast MD. 2013. Sea lice population and sex differences in P-glycoprotein expression and emamectin benzoate resistance on salmon farms in the Bay of Fundy, New Brunswick, Canada. Pest Management Science, 70: 905–914.
Ingvarsdóttir A, Birkett MA, Duce I, Mordue W, Pickett JA, Wadhams LJ, et al. 2002. Role of semiochemicals in mate location by parasitic sea louse, Lepeophtheirus salmonis. Journal of Chemical Ecology, 28: 2107–2117.
Isoe J, Zamora J, and Miesfeld RL. 2009. Molecular analysis of the Aedes aegypti carboxypeptidase gene family. Insect Biochemistry and Molecular Biology, 39(1): 68–73.
Itoh Y, Melamed E, Yang X, Kampf K, Wang S, Yehya N, et al. 2007. Dosage compensation is less effective in birds than in mammals. Journal of Biology, 6: 2.
Johnson SC, and Albright LJ. 1991. The developmental stages of Lepeophtheirus salmonis (Kroyer, 1837) (Copepoda: Caligidae). Canadian Journal of Zoology, 69: 929–950.
Lavazec C, Boudin C, Lacroix R, Bonnet S, Diop A, Thiberge S, et al. 2007. Carboxypeptidases B of Anopheles gambiae as targets for a Plasmodium falciparum transmission-blocking vaccine. Infection and Immunity, 75(4): 1635–1642.
Leong JS, Jantzen SG, Von Schalburg KR, Cooper GA, Messmer AM, Liao NY, et al. 2010. Salmo salar and Esox lucius full-length cDNA sequences reveal changes in evolutionary pressures on a post-tetraploidization genome. BMC Genomics, 11(1): 279.
Liu Y, and Bjelland HV. 2014. Estimating costs of sea lice control strategy in Norway. Preventive Veterinary Medicine, 117(3-4): 469–477.
Mack PD, Kapelnikov A, Heifetz Y, and Bender M. 2006. Mating-responsive genes in reproductive tissues of female Drosophila melanogaster. Proceedings of the National Academy of Sciences of the United States of America, 103(27): 10358–10363.
Mank JE. 2009. The W, X, Y and Z of sex-chromosome dosage compensation. Trends in Genetics, 25(5): 226–233.
McGraw LA, Clark AG, and Wolfner MF. 2008. Post-mating gene expression profiles of female Drosophila melanogaster in response to time and to four male accessory gland proteins. Genetics, 179(3): 1395–1408.
Meisel RP. 2011. Towards a more nuanced understanding of the relationship between sex-biased gene expression and rates of protein-coding sequence evolution. Molecular Biology and Evolution, 28: 1893–1900.
Parsch J, and Ellegren H. 2013. The evolutionary causes and consequences of sex-biased gene expression. Nature Reviews Genetics, 14(2): 83–87.
Peng J, Zipperlen P, and Kubli E. 2005. Drosophila sex-peptide stimulates female innate immune system after mating via the toll and Imd pathways. Current Biology, 15(18): 1690–1694.
Poley J, Sutherland B, Jones S, Koop B, and Fast M. 2016. Sex-biased gene expression and sequence conservation in Atlantic and Pacific salmon lice (Lepeophtheirus salmonis). BMC Genomics, 17(1): 483.
Powell JR, and Moriyama EN. 1997. Evolution of codon usage bias in Drosophila. Proceedings of the National Academy of Sciences of the United States of America, 94(15): 7784–7790.
Ram KR, and Wolfner MF. 2007. Seminal influences: Drosophila Acps and the molecular interplay between males and females during reproduction. Integrative & Comparative Biology, 47(3): 427–445.
Ramm SA, Schärer L, Ehmcke J, and Wistuba J. 2014. Sperm competition and the evolution of spermatogenesis. Molecular Human Reproduction, 20(12): 1169–1179.
Ritchie G, Mordue(Luntz) AJ, Pike AW, and Raec GH. 1996. Observations on mating and reproductive behaviour of Lepeophtheirus salmonis, Kroyer (Copepoda: Caligidae). Journal of Experimental Marine Biology and Ecology, 201(1-2): 285–298.
Rogers DW, Whitten MMA, Thailayil J, Soichot J, Levashina EA, and Catteruccia F. 2008. Molecular and cellular components of the mating machinery in Anopheles gambiae females. Proceedings of the National Academy of Sciences of the United States of America, 105(49): 19390–19395.
Schneider WJ. 1996. Vitellogenin receptors: oocyte-specific members of the low-density lipoprotein receptor supergene family. International Review of Cytology, 166: 103–137.
Sirot LK, Wong A, Chapman T, and Wolfner MF. 2015. Sexual conflict and seminal fluid proteins: a dynamic landscape of sexual interactions. Cold Spring Harbor Perspectives in Biology, 7(2): a017533.
Smart LE, Aradottir GI, and Bruce TJA. 2014. Chapter 6—Role of semiochemicals in integrated pest management. In Integrated pest management. Edited by DP Abrol. Academic Press, San Diego, California. pp. 93–109.
Soller M, Bownes M, and Kubli E. 1997. Mating and sex peptide stimulate the accumulation of yolk in oocytes of Drosophila melanogaster. European Journal of Biochemistry, 243(3): 732–738.
Stephenson JF. 2012. The chemical cues of male sea lice Lepeophtheirus salmonis encourage others to move between host Atlantic salmon Salmo salar. Journal of Fish Biology, 81: 1118–1123.
Sutherland BJG, Poley JD, Igboeli OO, Jantzen JR, Fast MD, Koop BF, et al. 2015. Transcriptomic responses to emamectin benzoate in Pacific and Atlantic Canada salmon lice Lepeophtheirus salmonis with differing levels of drug resistance. Evolutionary Applications, 8(2): 133–148.
Swanson WJ, and Vacquier VD. 2002. The rapid evolution of reproductive proteins. Nature Reviews Genetics, 3: 137–144.
Tufail M, and Takeda M. 2009. Insect vitellogenin/lipophorin receptors: molecular structures, role in oogenesis, and regulatory mechanisms. Journal of Insect Physiology, 55(2): 88–104.
Wagner GN, Fast MD, and Johnson SC. 2008. Physiology and immunology of Lepeophtheirus salmonis infections of salmonids. Trends in Parasitology, 24(4): 176–183.
Warr E, Aguilar R, Dong Y, Mahairaki V, and Dimopoulos G. 2007. Spatial and sex-specific dissection of the Anopheles gambiae midgut transcriptome. BMC Genomics, 8(37): 1–11.
Wernersson R, and Pedersen AG. 2003. RevTrans: multiple alignment of coding DNA from aligned amino acid sequences. Nucleic Acids Research, 31(13): 3537–3539.
Whyte SK, Westcott JD, Elmoslemany A, Hammell KL, and Revie CW. 2013. A fixed-dose approach to conducting emamectin benzoate tolerance assessments on field- collected sea lice, Lepeophtheirus salmonis. Journal of Fish Diseases, 36: 283–292.
Willadsen P, Bird P, Cobon GS, and Hungerford J. 1995. Commercialisation of a recombinant vaccine against Boophilus microplus. Parasitology, 110: S43–50.
Wolf JBW, and Bryk J. 2011. General lack of global dosage compensation in ZZ/ZW systems? Broadening the perspective with RNA-seq. BMC Genomics, 12: 91.
Wotton JH, Purcell SL, Covello JM, Koop BF, and Fast MD. 2014. Comparison of host selection and gene expression of adult Lepeophtheirus salmonis and Salmo salar during a cohabitation of initially infected and uninfected fish. Journal of Aquaculture Research and Development, 3: 226–233.
Yang Z. 2007. PAML 4: phylogenetic analysis by maximum likelihood. Molecular Biology and Evolution, 24(8): 1586–1591.
Yang Z, and Nielsen R. 2000. Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Molecular Biology and Evolution, 17(1): 32–43.
Zhang Z, Hambuch TM, and Parsch J. 2004. Molecular evolution of sex-biased genes in Drosophila. Molecular Biology and Evolution, 21(11): 2130–2139.

Supplementary Materials

Supplementary Material 1 (XLSX / 33.3 KB)
Supplementary Material 2 (XLSX / 44.6 KB)

Information & Authors

Information

Published In

cover image FACETS
FACETS
Volume 4Number 1June 2019
Pages: 254 - 274
Editor: Kristi M. Miller

History

Received: 1 June 2018
Accepted: 28 January 2019
Version of record online: 27 June 2019

Data Availability Statement

All relevant data are within the paper and in the Supplementary Material.

Key Words

  1. Copepoda
  2. Lepeophtheirus salmonis
  3. salmon lice
  4. sex-biased
  5. transcriptomics
  6. evolution

Sections

Subjects

Authors

Affiliations

Dylan R. Michaud [email protected]
Hoplite Research Group, Department of Pathology and Microbiology, Atlantic Veterinary College, UPEI, Charlottetown, PE C1A 4P3, Canada
Jordan D. Poley
Hoplite Research Group, Department of Pathology and Microbiology, Atlantic Veterinary College, UPEI, Charlottetown, PE C1A 4P3, Canada
Center for Aquaculture Technologies Canada, 20 Hope Street, Souris PE C0A 2B0, Canada
Mark D. Fast
Hoplite Research Group, Department of Pathology and Microbiology, Atlantic Veterinary College, UPEI, Charlottetown, PE C1A 4P3, Canada

Author Contributions

JDP and MDF conceived and designed the study.
DRM and JDP performed the experiments/collected the data.
DRM and JDP analyzed and interpreted the data.
MDF contributed resources.
All drafted or revised the manuscript.

Competing Interests

The authors have declared that no competing interests exist.

Metrics & Citations

Metrics

Other Metrics

Citations

Cite As

Export Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

1. Sex differences in the early life stages of the salmon louse Lepeophtheirus salmonis (Copepoda: Caligidae)
2. Sex-Biased Gene Expression and Evolution in the Cerebrum and Syrinx of Chinese Hwamei (Garrulax canorus)

View Options

View options

PDF

View PDF

Get Access

Media

Media

Other

Tables

Share Options

Share

Share the article link

Share on social media