An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony

Humble E, Paijmans A, Forcada J, Hoffman J. An 85K SNP Array Uncovers Inbreeding and Cryptic Relatedness in an Antarctic Fur Seal Breeding Colony. G3: Genes, Genomes, Genetics. 2020;10(8):2787-2799.High density single nucleotide polymorphism (SNP) arrays allow large numbers of individuals to be rapidly and cost-effectively genotyped at large numbers of genetic markers. However, despite being widely used in studies of humans and domesticated plants and animals, SNP arrays are lacking for most wild organisms. We developed a custom 85K Affymetrix Axiom array for an intensively studied pinniped, the Antarctic fur seal (Arctocephalus gazella). SNPs were discovered from a combination of genomic and transcriptomic resources and filtered according to strict criteria. Out of a total of 85,359 SNPs tiled on the array, 75,601 (88.6%) successfully converted and were polymorphic in 270 animals from a breeding colony at Bird Island in South Georgia. Evidence was found for inbreeding, with three genomic inbreeding coefficients being strongly intercorrelated and the proportion of the genome in runs of homozygosity being non-zero in all individuals. Furthermore, analysis of genomic relatedness coefficients identified previously unknown first-degree relatives and multiple second-degree relatives among a sample of ostensibly unrelated individuals. Such “cryptic relatedness” within fur seal breeding colonies may increase the likelihood of consanguineous matings and could therefore have implications for understanding fitness variation and mate choice. Finally, we demonstrate the cross-amplification potential of the array in three related pinniped species. Overall, our SNP array will facilitate future studies of Antarctic fur seals and has the potential to serve as a more general resource for the wider pinniped research community

Transcriptomic SNP discovery for custom genotyping arrays: impacts of sequence data, SNP calling method and genotyping technology on the probability of validation success

Background Single nucleotide polymorphism (SNP) discovery is an important goal of many studies. However, the number of ‘putative’ SNPs discovered from a sequence resource may not provide a reliable indication of the number that will successfully validate with a given genotyping technology. For this it may be necessary to account for factors such as the method used for SNP discovery and the type of sequence data from which it originates, suitability of the SNP flanking sequences for probe design, and genomic context. To explore the relative importance of these and other factors, we used Illumina sequencing to augment an existing Roche 454 transcriptome assembly for the Antarctic fur seal (Arctocephalus gazella). We then mapped the raw Illumina reads to the new hybrid transcriptome using BWA and BOWTIE2 before calling SNPs with GATK. The resulting markers were pooled with two existing sets of SNPs called from the original 454 assembly using NEWBLER and SWAP454. Finally, we explored the extent to which SNPs discovered using these four methods overlapped and predicted the corresponding validation outcomes for both Illumina Infinium iSelect HD and Affymetrix Axiom arrays. Results Collating markers across all discovery methods resulted in a global list of 34,718 SNPs. However, concordance between the methods was surprisingly poor, with only 51.0 % of SNPs being discovered by more than one method and 13.5 % being called from both the 454 and Illumina datasets. Using a predictive modeling approach, we could also show that SNPs called from the Illumina data were on average more likely to successfully validate, as were SNPs called by more than one method. Above and beyond this pattern, predicted validation outcomes were also consistently better for Affymetrix Axiom arrays. Conclusions Our results suggest that focusing on SNPs called by more than one method could potentially improve validation outcomes. They also highlight possible differences between alternative genotyping technologies that could be explored in future studies of non-model organisms

Chemical patterns of colony membership and mother- offspring similarity in Antarctic fur seals are reproducible

Tebbe J, Humble E, Stoffel MA, et al. Chemical patterns of colony membership and mother-offspring similarity in Antarctic fur seals are reproducible. PeerJ. 2020;8: e10131.Replication studies are essential for evaluating the validity of previous research findings. However, it has proven challenging to reproduce the results of ecological and evolutionary studies, partly because of the complexity and lability of many of the phenomena being investigated, but also due to small sample sizes, low statistical power and publication bias. Additionally, replication is often considered too difficult in field settings where many factors are beyond the investigator’s control and where spatial and temporal dependencies may be strong. We investigated the feasibility of reproducing original research findings in the field of chemical ecology by performing an exact replication of a previous study of Antarctic fur seals (Arctocephalus gazella). In the original study, skin swabs from 41 mother-offspring pairs from two adjacent breeding colonies on Bird Island, South Georgia, were analyzed using gas chromatography-mass spectrometry. Seals from the two colonies differed significantly in their chemical fingerprints, suggesting that colony membership may be chemically encoded, and mothers were also chemically similar to their pups, hinting at the possible involvement of phenotype matching in mother-offspring recognition. In the current study, we generated and analyzed chemical data from a non-overlapping sample of 50 mother-offspring pairs from the same two colonies 5 years later. The original results were corroborated in both hypothesis testing and estimation contexts, with p-values remaining highly significant and effect sizes, standardized between studies by bootstrapping the chemical data over individuals, being of comparable magnitude. However, exact replication studies are only capable of showing whether a given effect can be replicated in a specific setting. We therefore investigated whether chemical signatures are colony-specific in general by expanding the geographic coverage of our study to include pups from a total of six colonies around Bird Island. We detected significant chemical differences in all but a handful of pairwise comparisons between colonies. This finding adds weight to our original conclusion that colony membership is chemically encoded, and suggests that chemical patterns of colony membership not only persist over time but can also be generalized over space. Our study systematically confirms and extends our previous findings, while also implying more broadly that spatial and temporal heterogeneity need not necessarily negate the reproduction and generalization of ecological research findings

Genomic analysis of the carboxylesterase family in the salmon louse (Lepeophtheirus salmonis)

The pyrethroid deltamethrin and the macrocyclic lactone emamectin benzoate (EMB) are used to treat infestations of farmed salmon by parasitic salmon lice, Lepeophtheirus salmonis. While the efficacy of both compounds against Atlantic populations of the parasite has decreased as a result of the evolution of resistance, the molecular mechanisms of drug resistance in L. salmonis are currently not fully understood. The functionally diverse carboxylesterases (CaE) family includes members involved in pesticide resistance phenotypes of terrestrial arthropods. The present study had the objective to characterize the CaE family in L. salmonis and assess its role in drug resistance. L. salmonis CaE homologues were identified by homology searches in the parasite's transcriptome and genome. The transcript expression of CaEs predicted to be catalytically competent was studied using quantitative reverse-transcription PCR in drug susceptible and multi-resistant L. salmonis. The above strategy led to the identification of 21 CaEs genes/pseudogenes. Phylogenetic analyses assigned 13 CaEs to clades involved in neurodevelopmental signaling and cell adhesion, while three sequences were predicted to encode secreted enzymes. Ten CaEs were identified as being potentially catalytically competent. Transcript expression of acetylcholinesterase (ace1b) was significantly increased in multi-resistant lice compared to drug-susceptible L. salmonis, with transcript abundance further increased in preadult-II females following EMB exposure. In summary, results from the present study demonstrate that L. salmonis possesses fewer CaE gene family members than most arthropods characterized so far. Drug resistance in L. salmonis was associated with overexpression of ace1b

Genome-wide survey of cytochrome P450 genes in the salmon louse Lepeophtheirus salmonis (Krøyer, 1837)

Background The salmon louse (Lepeophtheirus salmonis) infests farmed and wild salmonid fishes, causing considerable economic damage to the salmon farming industry. Infestations of farmed salmon are controlled using a combination of non-medicinal approaches and veterinary drug treatments. While L. salmonis has developed resistance to most available salmon delousing agents, relatively little is known about the molecular mechanisms involved. Members of the cytochrome P450 (CYP) superfamily are typically monooxygenases, some of which are involved in the biosynthesis and metabolism of endogenous compounds, while others have central roles in the detoxification of xenobiotics. In terrestrial arthropods, insecticide resistance can be based on the enhanced expression of CYPs. The reported research aimed to characterise the CYP superfamily in L. salmonis and assess its potential roles in drug resistance. Methods Lepeophtheirus salmonis CYPs were identified by homology searches of the genome and transcriptome of the parasite. CYP transcript abundance in drug susceptible and multi-resistant L. salmonis was assessed by quantitative reverse transcription PCR, taking into account both constitutive expression and expression in parasites exposed to sublethal levels of salmon delousing agents, ecdysteroids and environmental chemicals. Results The above strategy led to the identification of 25 CYP genes/pseudogenes in L. salmonis, making its CYP superfamily the most compact characterised for any arthropod to date. Lepeophtheirus salmonis possesses homologues of a number of arthropod CYP genes with roles in ecdysteroid metabolism, such as the fruit fly genes disembodied, shadow, shade, spook and Cyp18a1. CYP transcript expression did not differ between one drug susceptible and one multi-resistant strain of L. salmonis. Exposure of L. salmonis to emamectin benzoate or deltamethrin caused the transcriptional upregulation of certain CYPs. In contrast, neither ecdysteroid nor benzo[a]pyrene exposure affected CYP transcription significantly. Conclusions The parasite L. salmonis is demonstrated to possess the most compact CYP superfamily characterised for any arthropod to date. The complement of CYP genes in L. salmonis includes conserved CYP genes involved in ecdysteroid biosynthesis and metabolism, as well as drug-inducible CYP genes. The present study does not provide evidence for a role of CYP genes in the decreased susceptibility of the multiresistant parasite strain studied

Time-to-response toxicity analysis as a method for drug susceptibility assessment in salmon lice

The salmon louse Lepeophtheirus salmonis (Krøyer, 1837) is an ectoparasite causing infections ofwild and farmed Atlantic salmon (Salmo salar L.) in the Northern hemisphere.While L. salmonis control at commercial mariculture sites increasingly employs non-medicinal approaches, such as cage designs reducing infection rates and biological control through cleaner fish, anti-parasitic drugs are still a requirement for effective fish health care. With only a limited range of salmon delousing agents available, all of which have been in use for more than a decade, drug resistance formation has been reported for different products. Successful resistance management requires reliable susceptibility assessment, which is usually achieved through L. salmonis bioassays. These tests involve the exposure of parasites to different drug concentrations and require significant numbers of suitable L. salmonis stages. The present study reports an alternative bioassay that is based on time-to-response toxicity analyses and can be carried outwith limited parasite numbers. The assay determines the median effective time (ET50), i.e., the time required until impaired swimming and/or attachment behaviour becomes apparent in 50% of parasites, by conducting repeated examinations of test animals starting at the timepointwhere exposure to a set drug concentration commences. This experimental approach further allows the estimation of the apparent drug susceptibility of individual L. salmonis by determining their time to response, which may prove useful in experiments designed to elucidate associations between genetic factors and the drug susceptibility phenotype of parasites. Three laboratory strains of L. salmonis differing in susceptibility to emamectin benzoate were characterised using standard 24 h bioassays and time-to-response toxicity assays. While both the median effective concentration (EC50) and the ET50 showed variability between experimental repeats, both types of bioassay consistently discriminated susceptible and drug-resistant L. salmonis laboratory strains. Statement of relevance: Infections by sea lice cause significant costs to the global salmon farming industry, which have been estimated to exceed €300 million per year worldwide. Control of sea lice still relies to a significant extent on chemical delousing; however, chemical control is threatened by resistance formation. Resistance can be combated by rotation between different drugs and strategic implementation of non-medicinal strategies. However, resistance management requires reliable and feasible methods of susceptibility assessment. The present study is a technical note introducing a novel approach to susceptibility assessments in sea lice. The method can be applied in susceptibility assessments on farms,where it offers the advantage of a reduced requirement of parasites for testing. In addition, the novel method allows deriving the times of parasite require to showa response after drug treatment has started, thus providing a variable characterizing the drug susceptibility phenotype of individual parasites. Accordingly, the bioassay approach presented here will be useful for studies aiming at unravelling the genetic determinants of drug resistance

Mutations in voltage-gated sodium channels from pyrethroid resistant salmon lice (Lepeophtheirus salmonis)

Background Parasitic salmon lice (Lepeophtheirus salmonis) cause high economic losses in Atlantic salmon farming. Pyrethroids, which block arthropod voltage‐gated sodium channels (Nav1), are used for salmon delousing. However, pyrethroid resistance is common in L. salmonis. The present study characterised Nav1 homologues in L. salmonis in order to identify channel mutations associated to resistance, called kdr (knockdown) mutations. Results Genome scans identified three L. salmonis Nav1 homologues, LsNav1.1, LsNav1.2 and LsNav1.3. Arthropod kdr mutations map to specific Nav1 regions within domains DI‐III, namely segments S5 and S6 and the linker helix connecting S4 and S5. The above channel regions were amplified by RT‐PCR and sequenced in deltamethrin‐susceptible and deltamethrin‐resistant L. salmonis. While LsNav1.1 and LsNav1.2 lacked nucleotide polymorphisms showing association to resistance, LsNav1.3 showed a non‐synonymous mutation in S5 of DII occurring in deltamethrin‐resistant parasites. The mutation is homologous to a previously described kdr mutation (I936V, numbering according to Musca domestica Vssc1) and was present in two pyrethroid‐resistant L. salmonis strains (allele frequencies of 0.800 and 0.357), but absent in two pyrethroid‐susceptible strains. Conclusions The present study indicates that a kdr‐mutation in LsNaV 1.3 may contribute to deltamethrin resistance in L. salmonis

A warmer environment can reduce sociability in an ectotherm

The costs and benefits of being social vary with environmental conditions, so individuals must weigh the balance between these trade-offs in response to changes in the environment. Temperature is a salient environmental factor that may play a key role in altering the costs and benefits of sociality through its effects on food availability, predator abundance, and other ecological parameters. In ectotherms, changes in temperature also have direct effects on physiological traits linked to social behaviour, such as metabolic rate and locomotor performance. In light of climate change, it is therefore important to understand the potential effects of temperature on sociality. Here, we took the advantage of a ‘natural experiment’ of threespine sticklebacks from contrasting thermal environments in Iceland: geothermally warmed water bodies (warm habitats) and adjacent ambient-temperature water bodies (cold habitats) that were either linked (sympatric) or physically distinct (allopatric). We first measured the sociability of wild-caught adult fish from warm and cold habitats after acclimation to a low and a high temperature. At both acclimation temperatures, fish from the allopatric warm habitat were less social than those from the allopatric cold habitat, whereas fish from sympatric warm and cold habitats showed no differences in sociability. To determine whether differences in sociability between thermal habitats in the allopatric population were heritable, we used a common garden breeding design where individuals from the warm and the cold habitat were reared at a low or high temperature for two generations. We found that sociability was indeed heritable but also influenced by rearing temperature, suggesting that thermal conditions during early life can play an important role in influencing social behaviour in adulthood. By providing the first evidence for a causal effect of rearing temperature on social behaviour, our study provides novel insights into how a warming world may influence sociality in animal populations.publishedVersio

Maternal inheritance of deltamethrin resistance in the salmon louse Lepeophtheirus salmonis (Krøyer) is associated with unique mtDNA haplotypes

Parasitic infections by the salmon louse, Lepeophtheirus salmonis (Krøyer), cause huge economic damage in salmon farming in the northern hemisphere, with combined treatment costs and production losses in 2014 having been estimated at US$ 350 million for Norway (annual production 1.25 million tonnes). The control of L. salmonis relies significantly on medicinal treatments, supplemented by non-pharmacological approaches. However, efficacy losses have been reported for several delousing agents, including the pyrethroid deltamethrin. The aim of the present study was to analyse the genetic basis of deltamethrin resistance in L. salmonis. Deltamethrin median effective concentrations (EC50) were 0.28 μg L-1 in the drug susceptible L. salmonis strain IoA-00 and 40.1 μg L-1 in the pyrethroid resistant strain IoA-02. IoA-00 and IoA-02 were crossed to produce families spanning one parental and three filial generations (P0, F1-F3). In three families derived from P0 crosses between an IoA-00 sire and an IoA-02 dam, 98.8% of F2 parasites (n = 173) were resistant, i.e. remained unaffected after exposure to 2.0 μg L-1 deltamethrin. F3 parasites from these crosses showed a deltamethrin EC50 of 9.66 μg L-1. In two families of the inverse orientation at P0 (IoA-02 sire x IoA-00 dam), 16.7% of F2 parasites were resistant (n = 84), while the deltamethrin EC50 in F3 animals was 0.26 μg L-1. The results revealed a predominantly maternal inheritance of deltamethrin resistance. The 15,947-nt mitochondrial genome was sequenced and compared among six unrelated L. salmonis strains and parasites sampled from wild salmon in 2010. IoA-02 and three further deltamethrin resistant strains, established from isolates originating from different regions of Scotland, showed almost identical mitochondrial haplotypes. In contrast, the mitochondrial genome was variable among susceptible strains and L. salmonis from wild hosts. Deltamethrin caused toxicity and depletion of whole body ATP levels in IoA-00 but not IoA-02 parasites. The maternal inheritance of deltamethrin resistance and its association with mitochondrial haplotypes suggests that pyrethroid toxicity in L. salmonis may involve molecular targets encoded by mitochondrial genes