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 log
2 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.
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.
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.
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).
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 (F
op) 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.