Search

Background: The salmon louse (Lepeophtheirus salmonis) is a parasite of wild and farmed salmonid fish, causing huge economic damage to the commercial farming of Atlantic salmon (Salmo salar) in the northern hemisphere. The avermectin emamectin benzoate (EMB) is widely used for salmon delousing. While resistance to EMB is widespread in Atlantic populations of L. salmonis, the molecular mechanisms of resistance remain to be elucidated. The aim of the present work was to obtain insights into potential EMB resistance mechanisms by identifying genetic and transcriptomic markers associated with EMB resistance., Results: Crosses were performed between EMB-susceptible and -resistant L. salmonis, sourced from two parental strains isolated in Scotland, producing fully pedigreed families. The EMB susceptibility of individual parasites was characterised using time-to-response bioassays. Parasites of two families were subjected to double digest restriction site-associated DNA sequencing (ddRAD-seq) for simultaneous discovery of single nucleotide polymorphisms (SNPs) and genotyping. Data analysis revealed that EMB resistance is associated with one quantitative trait locus (QTL) region on L. salmonis chromosome 5. Marker-trait association was confirmed by genotyping assays for 7 SNPs in two additional families. Furthermore, the transcriptome of male parasites of the EMB-susceptible and -resistant L. salmonis parental strains was assessed. Among eighteen sequences showing higher transcript expression in EMB-resistant as compared to drug-susceptible lice, the most strongly up-regulated gene is located in the above QTL region and shows high homology to β spectrin, a cytoskeleton protein that has roles in neuron architecture and function. Further genes differentially regulated in EMB-resistant lice include a glutathione S-transferase (GST), and genes coding for proteins with predicted roles in mitochondrial function, intracellular signalling or transcription., Conclusions: Major determinants of EMB resistance in L. salmonis are located on Chromosome 5. Resistance can be predicted using a limited number of genetic markers. Genes transcriptionally up-regulated in EMB resistant parasites include a β spectrin, a cytoskeletal protein with still incompletely understood roles in neuron structure and function, as well as glutathione S-transferase, an enzyme with potential roles in the biochemical defence against toxicants., Competing Interests: Declarations. Ethics approval and consent to participate: All research projects involving the Institute of Aquaculture at the University of Stirling are subjected to a thorough Ethical Review Process prior to any work being undertaken. The present research was assessed and approved by the Institute of Aquaculture Ethical Review Committee. Laboratory infections of Atlantic salmon with L. salmonis were carried out under UK Home Office license. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Polyploidy occurs naturally across eukaryotic lineages and has been harnessed in the domestication of many crops and vertebrates. In aquaculture, triploidy can be induced as a biocontainment strategy, as it creates a reproductive barrier preventing farm-to-wild introgression, which is currently a major conservation issue for the industry. However, recent work suggests that triploidisation protocols may, on occasion, produce 'failed triploids' displaying diploidy, aneuploidy and aberrant inheritance. The potentially negative consequences for conservation and animal welfare motivate the need for methods to evaluate the success of ploidy-manipulation protocols early in the production process. We developed a semi-automated version of the MAC-PR (microsatellite DNA allele counting - peak ratios) method to resolve the allelic configuration of large numbers of individuals across a panel of microsatellite markers that can be used to infer ploidy, pedigree and inheritance aberrations. We demonstrate an application of the approach using material from a series of Atlantic salmon (Salmo salar) breeding experiments where ploidy was manipulated using a hydrostatic pressure treatment. We validated the approach to infer ploidy against blood smears, finding a > 99% agreement between these methods, and demonstrate its potential utility to infer ploidy as early as the embryonic stage. Furthermore, we present tools to assign diploid and triploid progeny to families and to detect aberrant inheritance, which may be useful for breeding programmes that utilise ploidy manipulation techniques. The approach adds to the ploidy verification toolbox. The increased precision in detecting ploidy and inheritance aberrations will facilitate the ability of triploidisation programmes to prevent farm-to-wild introgression., (© 2024 The Author(s). Molecular Ecology Resources published by John Wiley & Sons Ltd.)

The mussel industry faces challenges such as low and inconsistent levels of larvae settlement and poor-quality spat, leading to variable production. However, mussel farming remains a vital sustainable and environmentally responsible method for producing protein, fostering ecological responsibility in the aquaculture sector. We investigate the population connectivity and larval dispersion of blue mussels (Mytilus edulis) in Scottish waters, as a case study, using a multidisciplinary approach that combined genetic data and particle modelling. This research allows us to develop a thorough understanding of blue mussel population dynamics in mid-latitude fjord regions, to infer gene-flow patterns, and to estimate population divergence. Our findings reveal a primary south-to-north particle transport direction and the presence of five genetic clusters. We discover a significant and continuous genetic material exchange among populations within the study area, with our biophysical model's outcomes aligning with our genetic observations. Additionally, our model reveals a robust connection between the southwest coast and the rest of the west coast. This study will guide the preservation of mussel farming regions, ensuring sustainable populations that contribute to marine ecosystem health and resilience., (© 2024. The Author(s).)

Cleaner fish species have gained great importance in the control of sea lice, among them, lumpfish (Cyclopterus lumpus) has become one of the most popular. Lumpfish life cycle has been closed, and hatchery reproduction is now possible, however, current production is reliant on wild caught broodstock to meet the increasing demand. Selective breeding practices are called to play an important role in the successful breeding of most aquaculture species, including lumpfish. In this study we analysed a lumpfish population for the identification of genomic markers linked to production traits. Sequencing of RAD libraries allowed us to identify, 7193 informative markers within the sampled individuals. Genome wide association analysis for sex, weight, condition factor and standard length was performed. One single major QTL region was identified for sex, while nine QTL regions were detected for weight, and three QTL regions for standard length. A total of 177 SNP markers of interest (from QTL regions) and 399 high F st SNP markers were combined in a low-density panel, useful to obtain relevant genetic information from lumpfish populations. Moreover, a robust combined subset of 29 SNP markers (10 associated to sex, 14 to weight and 18 to standard length) provided over 90% accuracy in predicting the animal's phenotype by machine learning. Overall, our findings provide significant insights into the genetic control of important traits in lumpfish and deliver important genomic resources that will facilitate the establishment of selective breeding programmes in lumpfish., Competing Interests: Declaration of Competing Interest T.C.R. and E.J.F. are employed by Otter Ferry SeaFish Ltd., a funding partner of the research project. The other authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

The current methods used for producing triploid Atlantic salmon are generally reliable but not infallible, and each batch of triploids must be validated to ensure consumer trust and licensing compliance. Microsatellites have recently been shown to offer a cheaper and more convenient alternative to traditional flow cytometry for triploidy validation in a commercial setting. However, incubating eggs to at least the eyed stage for microsatellite validation poses challenges, such as reduced quality and performance of triploids produced from later eggs in the stripping season. To address these issues, we propose another option: extracting DNA from recently fertilised eggs for use in conjunction with microsatellite validation. To achieve this, we have developed an optimized protocol for HotSHOT extraction that can rapidly and cheaply extract DNA from Atlantic salmon eggs, which can then be used for triploidy validation through microsatellites. Our approach offers a simpler and more cost-effective way to validate triploidy, without the need for skilled dissection or expensive kits., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Howard et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

A non-synonymous single nucleotide polymorphism (SNP) underlies a diallelic allozyme polymorphism at the mitochondrial NADP-dependent mMEP-2* locus in Atlantic salmon (Salmo salar L.). The resultant amino acid substitution, which alters the charge of the allelic products, matches the differential mobility of the two allozyme alleles, whereas allozyme and SNP assays revealed genotyping concordance in 257 of 258 individuals. A single mismatch, homozygous allozyme vs. heterozygote SNP, suggests the presence of a second, less common null allele., (© 2022 The Authors. Journal of Fish Biology published by John Wiley & Sons Ltd on behalf of Fisheries Society of the British Isles.)

Genetic resistance to infectious pancreatic necrosis virus (IPNV) in Atlantic salmon is a rare example of a trait where a single locus (QTL) explains almost all of the genetic variation. Genetic marker tests based on this QTL on salmon chromosome 26 have been widely applied in selective breeding to markedly reduce the incidence of the disease. In the current study, whole genome sequencing and functional annotation approaches were applied to characterise genes and variants in the QTL region. This was complemented by an analysis of differential expression between salmon fry of homozygous resistant and homozygous susceptible genotypes challenged with IPNV. These analyses pointed to the NEDD-8 activating enzyme 1 (nae1) gene as a putative functional candidate underlying the QTL effect. The role of nae1 in IPN resistance was further assessed via CRISPR-Cas9 knockout of the nae1 gene and chemical inhibition of the nae1 protein activity in Atlantic salmon cell lines, both of which resulted in highly significant reduction in productive IPNV replication. In contrast, CRISPR-Cas9 knockout of a candidate gene previously purported to be a cellular receptor for the virus (cdh1) did not have a major impact on productive IPNV replication. These results suggest that nae1 is the causative gene underlying the major QTL affecting resistance to IPNV in salmon, provide further evidence for the critical role of neddylation in host-pathogen interactions, and highlight the value in combining high-throughput genomics approaches with targeted genome editing to understand the genetic basis of disease resistance., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Tambaqui (Colossoma macropomum, Cuvier, 1818) is the most economically important native freshwater fish species in Brazil. It can reach a total length of over 1 m and a weight of over 40 kg. The species displays a clear sex dimorphism in growth performance, with females reaching larger sizes at harvest. In aquaculture, the production of monosex populations in selective breeding programmes has been therefore identified as a key priority., Results: In the present study, a genetic linkage map was generated by double digest restriction-site associated DNA (ddRAD) sequencing from 248 individuals sampled from two F1 families. The map was constructed using 14,805 informative SNPs and spanned 27 linkage groups. From this, the tambaqui draft genome was improved, by ordering the scaffolds into chromosomes, and sex-linked markers were identified. A total of 235 markers on linkage group 26 showed a significant association with the phenotypic sex, supporting an XX/XY sex determination system in the species. The four most informative sex-linked markers were validated on another 206 sexed individuals, demonstrating an accuracy in predicting sex ranging from 90.0 to 96.7%., Conclusions: The genetic mapping and novel sex-linked DNA markers identified and validated offer new tools for rapid progeny sexing, thus supporting the development of monosex female production in the industry while also supporting breeding programmes of the species., (© 2021. The Author(s).)

Wild fish assemblages that aggregate within commercial marine aquaculture sites for feeding and shelter have been considered as a primary source of pathogenic parasites vectored to farmed fish maintained in net pens at an elevated density. In order to evaluate whether Ceratothoa oestroides (Isopoda, Cymothoidae), a generalist and pestilent isopod that is frequently found in Adriatic and Greek stocks of farmed European sea bass ( Dicentrarchus labrax ), transfers between wild and farmed fish, a RAD-Seq (restriction-site-associated DNA sequencing)-mediated genetic screening approach was employed. The double-digest RAD-Seq of 310 C. oestroides specimens collected from farmed European sea bass (138) and different wild farm-aggregating fish (172) identified 313 robust SNPs that evidenced a close genetic relatedness between the "wild" and "farmed" genotypes. ddRAD-Seq proved to be an effective method for detecting the discrete genetic structuring of C. oestroides and genotype intermixing between two populations. The parasite prevalence in the farmed sea bass was 1.02%, with a mean intensity of 2.0 and mean abundance of 0.02, while in the wild fish, the prevalence was 8.1%; the mean intensity, 1.81; and the mean abundance, 0.15. Such differences are likely a consequence of human interventions during the farmed fish's rearing cycle that, nevertheless, did not affect the transfer of C. oestroides .

Gilthead sea bream is an important target for both recreational and commercial fishing in Europe, where it is also one of the most important cultured fish. Its distribution ranges from the Mediterranean to the African and European coasts of the North-East Atlantic. Until now, the population genetic structure of this species in the wild has largely been studied using microsatellite DNA markers, with minimal genetic differentiation being detected. In this geographically widespread study, 958 wild gilthead sea bream from 23 locations within the Mediterranean Sea and Atlantic Ocean were genotyped at 1159 genome-wide SNP markers by RAD sequencing. Outlier analyses identified 18 loci potentially under selection. Neutral marker analyses identified weak subdivision into three genetic clusters: Atlantic, West, and East Mediterranean. The latter group could be further subdivided into an Ionian/Adriatic and an Aegean group using the outlier markers alone. Seascape analysis suggested that this differentiation was mainly due to difference in salinity, this being also supported by preliminary genomic functional analysis. These results are of fundamental importance for the development of proper management of this species in the wild and are a first step toward the study of the potential genetic impact of the sea bream aquaculture industry., Competing Interests: The authors and the commercial company involved declare that they have no competing financial interests or personal relationships that could have influenced the work reported in this paper. The commercial affiliation with Fios Genomics Ltd. does not alter our adherence to PLOS ONE policies on sharing data and materials.

Background: Triploid organisms have three sets of chromosomes. In Atlantic salmon, hydrostatic pressure treatment of newly fertilized eggs has been extensively used to produce triploids which are functionally sterile due to their unpaired chromosomes. These fish often perform poorly on commercial farms, sometimes without explanation. Inheritance patterns in individuals subjected to pressure treatment have not been investigated in Atlantic salmon thus far. However, work on other species suggests that this treatment can result in aberrant inheritance. We therefore studied this in Atlantic salmon by genotyping 16 polymorphic microsatellites in eyed eggs and juveniles which had been subjected to pressure-induction of triploidy. Communally reared juveniles including fish subjected to pressure-induction of triploidy and their diploid siblings were included as a control., Results: No diploid offspring were detected in any of the eggs or juveniles which were subjected to hydrostatic pressure; therefore, the induction of triploidy was highly successful. Aberrant inheritance was nevertheless observed in 0.9% of the eggs and 0.9% of the juveniles that had been subjected to pressure treatment. In the communally reared fish, 0.3% of the fish subjected to pressure treatment displayed aberrant inheritance, while their diploid controls displayed 0% aberrant inheritance. Inheritance errors included two eyed eggs lacking maternal DNA across all microsatellites, and, examples in both eggs and juveniles of either the maternal or paternal allele lacking in one of the microsatellites. All individuals displaying chromosome aberrations were otherwise triploid., Conclusions: This is the first study to document aberrant inheritance in Atlantic salmon that have been subjected to pressure-induction of triploidy. Our experiments unequivocally demonstrate that even when induction of triploidy is highly successful, this treatment can cause chromosome aberrations in this species. Based upon our novel data, and earlier studies in other organisms, we hypothesize that in batches of Atlantic salmon where low to modest triploid induction rates have been reported, aberrant inheritance is likely to be higher than the rates observed here. Therefore, we tentatively suggest that this could contribute to the unexplained poor performance of triploid salmon that is occasionally reported in commercial aquaculture. These hypotheses require further investigation.

Background: Domestication is the process by which organisms become adapted to the human-controlled environment. Since the selection pressures that act upon cultured and natural populations differ, adaptations that favour life in the domesticated environment are unlikely to be advantageous in the wild. Elucidation of the differences between wild and domesticated Atlantic salmon may provide insights into some of the genomic changes occurring during domestication, and, help to predict the evolutionary consequences of farmed salmon escapees interbreeding with wild conspecifics. In this study the transcriptome of the offspring of wild and domesticated Atlantic salmon were compared using a common-garden experiment under standard hatchery conditions and in response to an applied crowding stressor., Results: Transcriptomic differences between wild and domesticated crosses were largely consistent between the control and stress conditions, and included down-regulation of environmental information processing, immune and nervous system pathways and up-regulation of genetic information processing, carbohydrate metabolism, lipid metabolism and digestive and endocrine system pathways in the domesticated fish relative to their wild counterparts, likely reflective of different selection pressures acting in wild and cultured populations. Many stress responsive functions were also shared between crosses and included down-regulation of cellular processes and genetic information processing and up-regulation of some metabolic pathways, lipid and energy in particular. The latter may be indicative of mobilization and reallocation of energy resources in response to stress. However, functional analysis indicated that a number of pathways behave differently between domesticated and wild salmon in response to stress. Reciprocal F1 hybrids permitted investigation of inheritance patterns that govern transcriptomic differences between these genetically divergent crosses. Additivity and maternal dominance accounted for approximately 42 and 25% of all differences under control conditions for both hybrids respectively. However, the inheritance of genes differentially expressed between crosses under stress was less consistent between reciprocal hybrids, potentially reflecting maternal environmental effects., Conclusion: We conclude that there are transcriptomic differences between the domesticated and wild salmon strains studied here, reflecting the different selection pressures operating on them. Our results indicate that stress may affect certain biological functions differently in wild, domesticated and hybrid crosses and these should be further investigated.

Background: Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The 'genetically improved farmed tilapia' (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins., Results: A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284)., Conclusions: Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

Sustainable fisheries management requires detailed knowledge of population genetic structure. The European sprat is an important commercial fish distributed from Morocco to the Arctic circle, Baltic, Mediterranean, and Black seas. Prior to 2018, annual catch advice on sprat from the International Council for the Exploration of the Sea (ICES) was based on five putative stocks: (a) North Sea, (b) Kattegat-Skagerrak and Norwegian fjords, (c) Baltic Sea, (d) West of Scotland-southern Celtic Seas, and (e) English Channel. However, there were concerns that the sprat advice on stock size estimates management plan inadequately reflected the underlying biological units. Here, we used ddRAD sequencing to develop 91 SNPs that were thereafter used to genotype approximately 2,500 fish from 40 locations. Three highly distinct and relatively homogenous genetic groups were identified: (a) Norwegian fjords; (b) Northeast Atlantic including the North Sea, Kattegat-Skagerrak, Celtic Sea, and Bay of Biscay; and (c) Baltic Sea. Evidence of genetic admixture and possibly physical mixing was detected in samples collected from the transition zone between the North and Baltic seas, but not between any of the other groups. These results have already been implemented by ICES with the decision to merge the North Sea and the Kattegat-Skagerrak sprat to be assessed as a single unit, thus demonstrating that genetic data can be rapidly absorbed to align harvest regimes and biological units., Competing Interests: None declared, (© 2020 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Tilapias (family Cichlidae) are of importance in aquaculture and fisheries. Hybridisation and introgression are common within tilapia genera but are difficult to analyse due to limited numbers of species-specific genetic markers. We tested the potential of double digested restriction-site associated DNA (ddRAD) sequencing for discovering single nucleotide polymorphism (SNP) markers to distinguish between 10 tilapia species. Analysis of ddRAD data revealed 1,371 shared SNPs in the de novo-based analysis and 1,204 SNPs in the reference-based analysis. Phylogenetic trees based on these two analyses were very similar. A total of 57 species-specific SNP markers were found among the samples analysed of the 10 tilapia species. Another set of 62 species-specific SNP markers was identified from a subset of four species which have often been involved in hybridisation in aquaculture: 13 for Oreochromis niloticus, 23 for O. aureus, 12 for O. mossambicus and 14 for O. u. hornorum. A panel of 24 SNPs was selected to distinguish among these four species and validated using 91 individuals. Larger numbers of SNP markers were found that could distinguish between the pairs of species within this subset. This technique offers potential for the investigation of hybridisation and introgression among tilapia species in aquaculture and in wild populations.

Local adaptations can determine the potential of populations to respond to environmental changes, yet adaptive genetic variation is commonly ignored in models forecasting species vulnerability and biogeographical shifts under future climate change. Here we integrate genomic and ecological modeling approaches to identify genetic adaptations associated with climate in two cryptic forest bats. We then incorporate this information directly into forecasts of range changes under future climate change and assessment of population persistence through the spread of climate-adaptive genetic variation (evolutionary rescue potential). Considering climate-adaptive potential reduced range loss projections, suggesting that failure to account for intraspecific variability can result in overestimation of future losses. On the other hand, range overlap between species was projected to increase, indicating that interspecific competition is likely to play an important role in limiting species' future ranges. We show that although evolutionary rescue is possible, it depends on a population's adaptive capacity and connectivity. Hence, we stress the importance of incorporating genomic data and landscape connectivity in climate change vulnerability assessments and conservation management., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Background: Arapaima gigas (Schinz, 1822) is the largest freshwater scaled fish in the world, and an emerging species for tropical aquaculture development. Conservation of the species, and the expansion of aquaculture requires the development of genetic tools to study polymorphism, differentiation, and stock structure. This study aimed to investigate genomic polymorphism through ddRAD sequencing, in order to identify a panel of single nucleotide polymorphisms (SNPs) and to simultaneously assess genetic diversity and structure in wild (from rivers Amazon, Solimões, Tocantins and Araguaia) and captive populations., Results: Compared to many other teleosts, the degree of polymorphism in A. gigas was low with only 2.3% of identified RAD-tags (135 bases long) containing SNPs. A panel of 393 informative SNPs was identified and screened across the five populations. Higher genetic diversity indices (number of polymorphic loci and private alleles, Shannon's Index and H O ) were found in populations from the Amazon and Solimões, intermediate levels in Tocantins and Captive, and very low levels in the Araguaia population. These results likely reflect larger population sizes from less urbanized environments in the Amazon basin compared to Araguaia. Populations were significantly differentiated with pairwise F ST values ranging from 0.086 (Amazon × Solimões) to 0.556 (Amazon × Araguaia). Mean pairwise relatedness among individuals was significant in all populations (P < 0.01), reflecting a degree of inbreeding possibly due to severe depletion of natural stocks, the species sedentary behaviour and possible sampling biases. Although Mantel test was not significant (P = 0.104; R 2  = 0.65), Bayesian analysis in STRUCTURE and discriminant analysis of principal components (DAPC) showed populations of Amazon and Solimões to be genetically differentiated from Araguaia, with Tocantins comprising individuals from both identified stocks., Conclusions: This relatively rapid genotyping by sequencing approach proved to be successful in delineating arapaima stocks. The approach and / or SNP panels identified should prove valuable for more detailed genetic studies of arapaima populations, including the elucidation of the genetic status of described discrete morphotypes and aid in delivery of conservation programs to maintain genetic diversity in reservoirs across the Amazon region.

MicroRNAs (miRNAs) are important post-transcriptional gene expression regulators. Here, 448 different miRNA genes, including 17 novel miRNAs, encoding for 589 mature Atlantic salmon miRNAs were identified after sequencing 111 samples (fry, pathogen challenged fry, various developmental and adult tissues). This increased the reference miRNAome with almost one hundred genes. Prior to isomiR characterization (mature miRNA variants), the proportion of erroneous sequence variants (ESVs) arising in the analysis pipeline was assessed. The ESVs were biased towards 5' and 3' end of reads in unexpectedly high proportions indicating that measurements of ESVs rather than Phred score should be used to avoid misinterpreting ESVs as isomiRs. Forty-three isomiRs were subsequently discovered. The biological effect of the isomiRs measured as increases in target diversity was small (<3%). Five miRNA genes showed allelic variation that had a large impact on target gene diversity if present in the seed. Twenty-one miRNAs were ubiquitously expressed while 31 miRNAs showed predominant expression in one or few tissues, indicating housekeeping or tissue specific functions, respectively. The miR-10 family, known to target Hox genes, were highly expressed in the developmental stages. The proportion of miR-430 family members, participating in maternal RNA clearance, was high at the earliest developmental stage.

In order to obtain an insight into genomic changes and associated evolution and adaptation of Infectious Pancreatic Necrosis Virus (IPNV), the complete coding genomes of 57 IPNV isolates collected from Scottish aquafarms from 1982 to 2014 were sequenced and analysed. Phylogenetic analysis of the sequenced IPNV strains showed separate clustering of genogroups I, II, III and V. IPNV isolates with genetic reassortment of segment A/B of genogroup III/II were determined. About 59 % of the IPNV isolates belonged to the persistent type and 32 % to the low-virulent type, and only one highly pathogenic strain (1.79 %) was identified. Codon adaptation index calculations indicated that the IPNV major capsid protein VP2 has adapted to its salmonid host. Under-representation of CpG dinucleotides in the IPNV genome to minimize detection by the innate immunity receptors, and observed positive selection in the virulence determination sites of VP2 embedded in the variable region of the main antigenic region, suggest an immune escape mechanism driving virulence evolution. The prevalence of mostly persistent genotypes, together with the assumption of adaptation and immune escape, indicates that IPNV is evolving with the host.

The development of Genotyping-By-Sequencing (GBS) technologies enables cost-effective analysis of large numbers of Single Nucleotide Polymorphisms (SNPs), especially in "non-model" species. Nevertheless, as such technologies enter a mature phase, biases and errors inherent to GBS are becoming evident. Here, we evaluated the performance of double digest Restriction enzyme Associated DNA (ddRAD) sequencing in SNP genotyping studies including high number of samples. Datasets of sequence data were generated from three marine teleost species (>5500 samples, >2.5 × 10 12 bases in total), using a standardized protocol. A common bioinformatics pipeline based on STACKS was established, with and without the use of a reference genome. We performed analyses throughout the production and analysis of ddRAD data in order to explore (i) the loss of information due to heterogeneous raw read number across samples; (ii) the discrepancy between expected and observed tag length and coverage; (iii) the performances of reference based vs. de novo approaches; (iv) the sources of potential genotyping errors of the library preparation/bioinformatics protocol, by comparing technical replicates. Our results showed use of a reference genome and a posteriori genotype correction improved genotyping precision. Individual read coverage was a key variable for reproducibility; variance in sequencing depth between loci in the same individual was also identified as an important factor and found to correlate to tag length. A comparison of downstream analysis carried out with ddRAD vs single SNP allele specific assay genotypes provided information about the levels of genotyping imprecision that can have a significant impact on allele frequency estimations and population assignment. The results and insights presented here will help to select and improve approaches to the analysis of large datasets based on RAD-like methodologies., (Crown Copyright © 2018. Published by Elsevier B.V. All rights reserved.)

Unraveling adaptive genetic variation represents, in addition to the estimate of population demographic parameters, a cornerstone for the management of aquatic natural living resources, which, in turn, represent the raw material for breeding programs. The turbot ( Scophthalmus maximus ) is a marine flatfish of high commercial value living on the European continental shelf. While wild populations are declining, aquaculture is flourishing in southern Europe. We evaluated the genetic structure of turbot throughout its natural distribution range (672 individuals; 20 populations) by analyzing allele frequency data from 755 single nucleotide polymorphism discovered and genotyped by double-digest RAD sequencing. The species was structured into four main regions: Baltic Sea, Atlantic Ocean, Adriatic Sea, and Black Sea, with subtle differentiation apparent at the distribution margins of the Atlantic region. Genetic diversity and effective population size estimates were highest in the Atlantic populations, the area of greatest occurrence, while turbot from other regions showed lower levels, reflecting geographical isolation and reduced abundance. Divergent selection was detected within and between the Atlantic Ocean and Baltic Sea regions, and also when comparing these two regions with the Black Sea. Evidence of parallel evolution was detected between the two low salinity regions, the Baltic and Black seas. Correlation between genetic and environmental variation indicated that temperature and salinity were probably the main environmental drivers of selection. Mining around the four genomic regions consistently inferred to be under selection identified candidate genes related to osmoregulation, growth, and resistance to diseases. The new insights are useful for the management of turbot fisheries and aquaculture by providing the baseline for evaluating the consequences of turbot releases from restocking and farming.

Genome size varies significantly across eukaryotic taxa and the largest changes are typically driven by macro-mutations such as whole genome duplications (WGDs) and proliferation of repetitive elements. These two processes may affect the evolutionary potential of lineages by increasing genetic variation and changing gene expression. Here, we elucidate the evolutionary history and mechanisms underpinning genome size variation in a species-rich group of Neotropical catfishes (Corydoradinae) with extreme variation in genome size-0.6 to 4.4 pg per haploid cell. First, genome size was quantified in 65 species and mapped onto a novel fossil-calibrated phylogeny. Two evolutionary shifts in genome size were identified across the tree-the first between 43 and 49 Ma (95% highest posterior density (HPD) 36.2-68.1 Ma) and the second at approximately 19 Ma (95% HPD 15.3-30.14 Ma). Second, restriction-site-associated DNA (RAD) sequencing was used to identify potential WGD events and quantify transposable element (TE) abundance in different lineages. Evidence of two lineage-scale WGDs was identified across the phylogeny, the first event occurring between 54 and 66 Ma (95% HPD 42.56-99.5 Ma) and the second at 20-30 Ma (95% HPD 15.3-45 Ma) based on haplotype numbers per contig and between 35 and 44 Ma (95% HPD 30.29-64.51 Ma) and 20-30 Ma (95% HPD 15.3-45 Ma) based on SNP read ratios. TE abundance increased considerably in parallel with genome size, with a single TE-family (TC1-IS630-Pogo) showing several increases across the Corydoradinae, with the most recent at 20-30 Ma (95% HPD 15.3-45 Ma) and an older event at 35-44 Ma (95% HPD 30.29-64.51 Ma). We identified signals congruent with two WGD duplication events, as well as an increase in TE abundance across different lineages, making the Corydoradinae an excellent model system to study the effects of WGD and TEs on genome and organismal evolution., (© 2018 The Authors.)

The continuing decline in forest elephant ( Loxodonta cyclotis ) numbers due to poaching and habitat reduction is driving the search for new tools to inform management and conservation. For dense rainforest species, basic ecological data on populations and threats can be challenging and expensive to collect, impeding conservation action in the field. As such, genetic monitoring is being increasingly implemented to complement or replace more burdensome field techniques. Single-nucleotide polymorphisms (SNPs) are particularly cost-effective and informative markers that can be used for a range of practical applications, including population census, assessment of human impact on social and genetic structure, and investigation of the illegal wildlife trade. SNP resources for elephants are scarce, but next-generation sequencing provides the opportunity for rapid, inexpensive generation of SNP markers in nonmodel species. Here, we sourced forest elephant DNA from 23 samples collected from 10 locations within Gabon, Central Africa, and applied double-digest restriction-site-associated DNA (ddRAD) sequencing to discover 31,851 tags containing SNPs that were reduced to a set of 1,365 high-quality candidate SNP markers. A subset of 115 candidate SNPs was then selected for assay design and validation using 56 additional samples. Genotyping resulted in a high conversion rate (93%) and a low per allele error rate (0.07%). This study provides the first panel of 107 validated SNP markers for forest elephants. This resource presents great potential for new genetic tools to produce reliable data and underpin a step-change in conservation policies for this elusive species.

Climate change is a major threat to global biodiversity that will produce a range of new selection pressures. Understanding species responses to climate change requires an interdisciplinary perspective, combining ecological, molecular and environmental approaches. We propose an applied integrated framework to identify populations under threat from climate change based on their extent of exposure, inherent sensitivity due to adaptive and neutral genetic variation and range shift potential. We consider intraspecific vulnerability and population-level responses, an important but often neglected conservation research priority. We demonstrate how this framework can be applied to vertebrates with limited dispersal abilities using empirical data for the bat Plecotus austriacus. We use ecological niche modelling and environmental dissimilarity analysis to locate areas at high risk of exposure to future changes. Combining outlier tests with genotype-environment association analysis, we identify potential climate-adaptive SNPs in our genomic data set and differences in the frequency of adaptive and neutral variation between populations. We assess landscape connectivity and show that changing environmental suitability may limit the future movement of individuals, thus affecting both the ability of populations to shift their distribution to climatically suitable areas and the probability of evolutionary rescue through the spread of adaptive genetic variation among populations. Therefore, a better understanding of movement ecology and landscape connectivity is needed for predicting population persistence under climate change. Our study highlights the importance of incorporating genomic data to determine sensitivity, adaptive potential and range shift potential, instead of relying solely on exposure to guide species vulnerability assessments and conservation planning., (© 2017 The Authors. Molecular Ecology Resources Published by John Wiley & Sons Ltd.)

Fish species exhibit substantial variation in the degree of genetic differentiation between sex chromosome pairs, and therefore offer the opportunity to study the full range of sex chromosome evolution. We used restriction-site associated DNA sequencing (RAD-seq) to study the sex chromosomes of Characidium gomesi, a species with conspicuous heteromorphic ZW/ZZ sex chromosomes. We screened 9863 single-nucleotide polymorphisms (SNPs), corresponding to ~1 marker/100 kb distributed across the genome for sex-linked variation. With this data set, we identified 26 female-specific RAD loci, putatively located on the W chromosome, as well as 148 sex-associated SNPs showing significant differentiation (average F ST =0.144) between males and females, and therefore in regions of more recent divergence between the Z and W chromosomes. In addition, we detected 25 RAD loci showing extreme heterozygote deficiency in females but which were in Hardy-Weinberg equilibrium in males, consistent with degeneration of the W chromosome and therefore female hemizygosity. We validated seven female-specific and two sex-associated markers in a larger sample of C. gomesi, of which three localised to the W chromosome, thereby providing useful markers for sexing wild samples. Validated markers were evaluated in other populations and species of the genus Characidium, this exploration suggesting a rapid turnover of W-specific repetitive elements. Together, our analyses point to a complex origin for the sex chromosome of C. gomesi and highlight the utility of RAD-seq for studying the composition and evolution of sex chromosomes systems in wild populations.

With increasing interest in the use of triploid salmon in commercial aquaculture, gaining an understanding of how economically important pathogens affect triploid stocks is important. To compare the susceptibility of diploid and triploid Atlantic salmon (Salmo salar L.) to viral pathogens, fry were experimentally infected with Salmonid alphavirus sub-type 1 (SAV1), the aetiological agent of pancreas disease (PD) affecting Atlantic salmon aquaculture in Europe. Three groups of fry were exposed to the virus via different routes of infection: intraperitoneal injection (IP), bath immersion, or cohabitation (co-hab) and untreated fry were used as a control group. Mortalities commenced in the co-hab challenged diploid and triploid fish from 11 days post infection (dpi), and the experiment was terminated at 17 dpi. Both diploid and triploid IP challenged groups had similar levels of cumulative mortality at the end of the experimental period (41.1% and 38.9% respectively), and these were significantly higher (p < 0.01) than for the other challenge routes. A TaqMan-based quantitative PCR was used to assess SAV load in the heart, a main target organ of the virus, and also liver, which does not normally display any pathological changes during clinical infections, but exhibited severe degenerative lesions in the present study. The median viral RNA copy number was higher in diploid fish compared to triploid fish in both the heart and the liver of all three challenged groups. However, a significant statistical difference (p < 0.05) was only apparent in the liver of the co-hab groups. Diploid fry also displayed significantly higher levels of pancreatic and myocardial degeneration than triploids. This study showed that both diploid and triploid fry are susceptible to experimental SAV1 infection. The lower virus load seen in the triploids compared to the diploids may possibly be related to differences in cell metabolism between the two groups, however, further investigation is necessary to confirm this and also to assess the outcome of PD outbreaks in other developmental stages of the fish when maintained in commercial production systems.

Background: Fully isogenic lines in fish can be developed using "mitotic" gynogenesis (suppression of first zygotic mitosis following inactivation of the sperm genome). However, genome-wide verification of the steps in this process has seldom been applied. We used ddRADseq to generate SNP markers in a meiotic gynogenetic family of European seabass (Dicentrarchus labrax): (i) to verify the lack of paternal contribution in a meiotic gynogenetic family; (ii) to generate a gene-centromere map from this family; (iii) to identify telomeric markers that could distinguish mitotic gynogenetics from meiotic gynogenetics, which sometimes arise spontaneously in mitotic gynogenetic families., Results: From a single meiotic gynogenetic family consisting of 79 progeny, 42 million sequencing reads (Illumina, trimmed to 148 bases) resolved 6866 unique RAD-tags. The 340 male-informative SNP markers that were identified confirmed the lack of paternal contribution. A gene-centromere map was constructed based on 804 female-informative SNPs in 24 linkage groups (2n = 48) with a total length of 1251.02 cM (initial LG assignment was based on the seabass genome assembly, dicLab v1). Chromosome arm structure could be clearly discerned from the pattern of heterozygosity in each linkage group in 18 out of 24 LGs: the other six showed anomalies that appeared to be related to issues in the genome assembly., Conclusion: Genome-wide screening enabled substantive verification of the production of the gynogenetic family used in this study. The large number of telomeric and subtelomeric markers with high heterozygosity values in the meiotic gynogenetic family indicate that such markers could be used to clearly distinguish between meiotic and mitotic gynogenetics.

Genomic selection uses genome-wide marker information to predict breeding values for traits of economic interest, and is more accurate than pedigree-based methods. The development of high density SNP arrays for Atlantic salmon has enabled genomic selection in selective breeding programs, alongside high-resolution association mapping of the genetic basis of complex traits. However, in sibling testing schemes typical of salmon breeding programs, trait records are available on many thousands of fish with close relationships to the selection candidates. Therefore, routine high density SNP genotyping may be prohibitively expensive. One means to reducing genotyping cost is the use of genotype imputation, where selected key animals ( e.g. , breeding program parents) are genotyped at high density, and the majority of individuals ( e.g. , performance tested fish and selection candidates) are genotyped at much lower density, followed by imputation to high density. The main objectives of the current study were to assess the feasibility and accuracy of genotype imputation in the context of a salmon breeding program. The specific aims were: (i) to measure the accuracy of genotype imputation using medium (25 K) and high (78 K) density mapped SNP panels, by masking varying proportions of the genotypes and assessing the correlation between the imputed genotypes and the true genotypes; and (ii) to assess the efficacy of imputed genotype data in genomic prediction of key performance traits (sea lice resistance and body weight). Imputation accuracies of up to 0.90 were observed using the simple two-generation pedigree dataset, and moderately high accuracy (0.83) was possible even with very low density SNP data (∼250 SNPs). The performance of genomic prediction using imputed genotype data was comparable to using true genotype data, and both were superior to pedigree-based prediction. These results demonstrate that the genotype imputation approach used in this study can provide a cost-effective method for generating robust genome-wide SNP data for genomic prediction in Atlantic salmon. Genotype imputation approaches are likely to form a critical component of cost-efficient genomic selection programs to improve economically important traits in aquaculture., (Copyright © 2017 Tsai et al.)

High density linkage maps are useful tools for fine-scale mapping of quantitative trait loci, and characterization of the recombination landscape of a species' genome. Genomic resources for Atlantic salmon (Salmo salar) include a well-assembled reference genome, and high density single nucleotide polymorphism (SNP) arrays. Our aim was to create a high density linkage map, and to align it with the reference genome assembly. Over 96,000 SNPs were mapped and ordered on the 29 salmon linkage groups using a pedigreed population comprising 622 fish from 60 nuclear families, all genotyped with the 'ssalar01' high density SNP array. The number of SNPs per group showed a high positive correlation with physical chromosome length (r = 0.95). While the order of markers on the genetic and physical maps was generally consistent, areas of discrepancy were identified. Approximately 6.5% of the previously unmapped reference genome sequence was assigned to chromosomes using the linkage map. Male recombination rate was lower than females across the vast majority of the genome, but with a notable peak in subtelomeric regions. Finally, using RNA-Seq data to annotate the reference genome, the mapped SNPs were categorized according to their predicted function, including annotation of ∼2500 putative nonsynonymous variants. The highest density SNP linkage map for any salmonid species has been created, annotated, and integrated with the Atlantic salmon reference genome assembly. This map highlights the marked heterochiasmy of salmon, and provides a useful resource for salmonid genetics and genomics research., (Copyright © 2016 Tsai et al.)

Aquaculture production of cod has decreased from over 20,000 tonnes in 2009 to less than 2,000 tonnes in 2014 and the industry faces many challenges, one of which is high and unpredictably variable mortality rates in the early life stages. Hence, full-cycle farming with hatchery produced juveniles is still considered unprofitable compared to fisheries and on-growing of wild cod. In the present study, potential batch differences in progeny survival of wild-caught, hatchery-spawned Faroe Bank cod (Gadus morhua L.) were investigated at two defined periods during early life history; i) the embryo stage (60 day degrees post fertilisation) and ii) the fry stage (110 days post hatch), post metamorphosis. The fry stage experiment was conducted in three replicates (N = 300 per replicate), and a panel of three polymorphic microsatellite markers was used for parental analysis. Mean survival rate at the embryo stage was 69% (± 20% SD). Survival was positively associated with egg diameter (P < 0.01), explaining 90% of the variation in egg survival rates. The data were too scarce to conclude either way concerning a possible correlation between survival rates between the two periods (P < 0.10). Offspring from three batches (from a total of eight) dominated in the fry stage, contributing over 90% of the progeny, and results were consistent over all three replicate tanks. The skewed batch representation observed may be of relevance to the effective management of selective breeding programmes for cod.

Background: Sea lice have significant negative economic and welfare impacts on marine Atlantic salmon farming. Since host resistance to sea lice has a substantial genetic component, selective breeding can contribute to control of lice. Genomic selection uses genome-wide marker information to predict breeding values, and can achieve markedly higher accuracy than pedigree-based methods. Our aim was to assess the genetic architecture of host resistance to sea lice, and test the utility of genomic prediction of breeding values. Individual lice counts were measured in challenge experiments using two large Atlantic salmon post-smolt populations from a commercial breeding programme, which had genotypes for ~33 K single nucleotide polymorphisms (SNPs). The specific objectives were to: (i) estimate the heritability of host resistance; (ii) assess its genetic architecture by performing a genome-wide association study (GWAS); (iii) assess the accuracy of predicted breeding values using varying SNP densities (0.5 to 33 K) and compare it to that of pedigree-based prediction; and (iv) evaluate the accuracy of prediction in closely and distantly related animals., Results: Heritability of host resistance was significant (0.22 to 0.33) in both populations using either pedigree or genomic relationship matrices. The GWAS suggested that lice resistance is a polygenic trait, and no genome-wide significant quantitative trait loci were identified. Based on cross-validation analysis, genomic predictions were more accurate than pedigree-based predictions for both populations. Although prediction accuracies were highest when closely-related animals were used in the training and validation sets, the benefit of having genomic-versus pedigree-based predictions within a population increased as the relationships between training and validation sets decreased. Prediction accuracy reached an asymptote with a SNP density of ~5 K within populations, although higher SNP density was advantageous for cross-population prediction., Conclusions: Host resistance to sea lice in farmed Atlantic salmon has a significant genetic component. Phenotypes relating to host resistance can be predicted with moderate to high accuracy within populations, with a major advantage of genomic over pedigree-based methods, even at relatively sparse SNP densities. Prediction accuracies across populations were low, but improved with higher marker densities. Genomic selection can contribute to lice control in salmon farming.

Background: Hāpuku (Polyprion oxygeneios) is a member of the wreckfish family (Polyprionidae) and is highly regarded as a food fish. Although adults grow relatively slowly, juveniles exhibit low feed conversion ratios and can reach market size in 1-2 years, making P. oxygeneios a strong candidate for aquaculture. However, they can take over 5 years to reach sexual maturity in captivity and are not externally sexually dimorphic, complicating many aspects of broodstock management. Understanding the sex determination system of P. oxygeneios and developing accurate assays to assign genetic sex will contribute significantly towards its full-scale commercialisation., Results: DNA from parents and sexed offspring (n = 57) from a single family of captive bred P. oxygeneios was used as a template for double digestion Restriction-site Associated DNA (ddRAD) sequencing. Two libraries were constructed using SbfI - SphI and SbfI - NcoI restriction enzyme combinations, respectively. Two runs on an Illumina MiSeq platform generated 70,266,464 raw reads, identifying 19,669 RAD loci. A combined sex linkage map (1367 cM) was constructed based on 1575 Single Nucleotide Polymorphism (SNP) markers that resolved into 35 linkage groups. Sex-specific linkage maps were of similar size (1132 and 1168 cM for male and female maps respectively). A single major sex-determining locus, found to be heterogametic in males, was mapped to linkage group 14. Several markers were found to be in strong linkage disequilibrium with the sex-determining locus. Allele-specific PCR assays were developed for two of these markers, SphI6331 and SphI8298, and demonstrated to accurately differentiate sex in progeny within the same pedigree. Comparative genomic analyses indicated that many of the linkage groups within the P. oxygeneios map share a relatively high degree of homology with those published for the European seabass (Dicentrarchus labrax)., Conclusion: P. oxygeneios has an XX/XY sex determination system. Evaluation of allele-specific PCR assays, based on the two SNP markers most closely associated with phenotypic sex, indicates that a simple molecular assay for sexing P. oxygeneios should be readily attainable. The high degree of synteny observed with D. labrax should aid further molecular genetic study and exploitation of hāpuku as a food fish.

Background: Infectious Pancreatic Necrosis (IPN) is a highly contagious birnavirus disease of farmed salmonid fish, which often causes high levels of morbidity and mortality. A large host genetic component to resistance has been previously described for Atlantic salmon (Salmo salar L.), which mediates high mortality rates in some families and zero mortality in others. However, the molecular and immunological basis for this resistance is not yet fully known. This manuscript describes a global comparison of the gene expression profiles of resistant and susceptible Atlantic salmon fry following challenge with the IPN virus., Results: Salmon fry from two IPNV-resistant and two IPNV-susceptible full sibling families were challenged with the virus and sampled at 1 day, 7 days and 20 days post-challenge. Significant viral titre was observed in both resistant and susceptible fish at all timepoints, although generally at higher levels in susceptible fish. Gene expression profiles combined with gene ontology and pathway analyses demonstrated that while a clear immune response was observed in both resistant and susceptible fish, there were striking differences between the two phenotypes. The susceptible fish showed marked up-regulation of genes related to cytokine activity and inflammatory response that evidently failed to protect against the virus. In contrast, the resistant fish demonstrated a less pronounced immune response including up-regulation of genes relating to the M2 macrophage system., Conclusions: While only the susceptible phenotype shows appreciable mortality levels, both resistant and susceptible fish can become infected with IPNV. Susceptible fish are characterized by a much larger, yet ineffective, immune response, largely related to cytokine and inflammatory systems. Resistant fish demonstrate a more moderate, putative macrophage-mediated inflammatory response, which may contribute to their survival.

Background: Due to selective breeding, domesticated and wild Atlantic salmon are genetically diverged, which raises concerns about farmed escapees having the potential to alter the genetic composition of wild populations and thereby disrupting local adaptation. Documenting transcriptional differences between wild and domesticated stocks under controlled conditions is one way to explore the consequences of domestication and selection. We compared the transcriptomes of wild and domesticated Atlantic salmon embryos, by using a custom 44k oligonucleotide microarray to identify perturbed gene pathways between the two stocks, and to document the inheritance patterns of differentially-expressed genes by examining gene expression in their reciprocal hybrids., Results: Data from 24 array interrogations were analysed: four reciprocal cross types (W♀ × W♂, D♀ × W♂; W♀ × D♂, D♀ × D♂) × six biological replicates. A common set of 31,491 features on the microarrays passed quality control, of which about 62 % were assigned a KEGG Orthology number. A total of 6037 distinct genes were identified for gene-set enrichment/pathway analysis. The most highly enriched functional groups that were perturbed between the two stocks were cellular signalling and immune system, ribosome and RNA transport, and focal adhesion and gap junction pathways, relating to cell communication and cell adhesion molecules. Most transcripts that were differentially expressed between the stocks were governed by additive gene interaction (33 to 42 %). Maternal dominance and over-dominance were also prevalent modes of inheritance, with no convincing evidence for a stock effect., Conclusions: Our data indicate that even at this relatively early developmental stage, transcriptional differences exist between the two stocks and affect pathways that are relevant to wild versus domesticated environments. Many of the identified differentially perturbed pathways are involved in organogenesis, which is expected to be an active process at the eyed egg stage. The dominant effects are more largely due to the maternal line than to the origin of the stock. This finding is particularly relevant in the context of potential introgression between farmed and wild fish, since female escapees tend to have a higher spawning success rate compared to males.

Common pandora (Pagellus erythrinus) is a benthopelagic marine fish belonging to the teleost family Sparidae, and a newly recruited species in Mediterranean aquaculture. The paucity of genetic information relating to sparids, despite their growing economic value for aquaculture, provides the impetus for exploring the genomics of this fish group. Genomic tool development, such as genetic linkage maps provision, lays the groundwork for linking genotype to phenotype, allowing fine-mapping of loci responsible for beneficial traits. In this study, we applied ddRAD methodology to identify polymorphic markers in a full-sib family of common pandora. Employing the Illumina MiSeq platform, we sampled and sequenced a size-selected genomic fraction of 99 individuals, which led to the identification of 920 polymorphic loci. Downstream mapping analysis resulted in the construction of 24 robust linkage groups, corresponding to the karyotype of the species. The common pandora linkage map showed varying degrees of conserved synteny with four other teleost genomes, namely the European seabass (Dicentrarchus labrax), Nile tilapia (Oreochromis niloticus), stickleback (Gasterosteus aculeatus), and medaka (Oryzias latipes), suggesting a conserved genomic evolution in Sparidae. Our work exploits the possibilities of genotyping by sequencing to gain novel insights into genome structure and evolution. Such information will boost the study of cultured species and will set the foundation for a deeper understanding of the complex evolutionary history of teleosts., (Copyright © 2016 Manousaki et al.)

Background: The largest of the tuna species, Atlantic bluefin tuna (Thunnus thynnus), inhabits the North Atlantic Ocean and the Mediterranean Sea and is considered to be an endangered species, largely a consequence of overfishing. T. thynnus aquaculture, referred to as fattening or farming, is a capture based activity dependent on yearly renewal from the wild. Thus, the development of aquaculture practices independent of wild resources can provide an important contribution towards ensuring security and sustainability of this species in the longer-term. The development of such practices is today greatly assisted by large scale transcriptomic studies., Results: We have used pyrosequencing technology to sequence a mixed-tissue normalised cDNA library, derived from adult T. thynnus. A total of 976,904 raw sequence reads were assembled into 33,105 unique transcripts having a mean length of 893 bases and an N50 of 870. Of these, 33.4% showed similarity to known proteins or gene transcripts and 86.6% of them were matched to the congeneric Pacific bluefin tuna (Thunnus orientalis) genome, compared to 70.3% for the more distantly related Nile tilapia (Oreochromis niloticus) genome. Transcript sequences were used to develop a novel 15 K Agilent oligonucleotide DNA microarray for T. thynnus and comparative tissue gene expression profiles were inferred for gill, heart, liver, ovaries and testes. Functional contrasts were strongest between gills and ovaries. Gills were particularly associated with immune system, signal transduction and cell communication, while ovaries displayed signatures of glycan biosynthesis, nucleotide metabolism, transcription, translation, replication and repair., Conclusions: Sequence data generated from a novel mixed-tissue T. thynnus cDNA library provide an important transcriptomic resource that can be further employed for study of various aspects of T. thynnus ecology and genomics, with strong applications in aquaculture. Tissue-specific gene expression profiles inferred through the use of novel oligo-microarray can serve in the design of new and more focused transcriptomic studies for future research of tuna physiology and assessment of the welfare in a production environment.

Pangasianodon hypophthalmus (striped catfish) is an important aquaculture species and intensification of farming has increased disease problems, particularly Edwardsiella ictaluri. The effects of feeding β-glucans on immune gene expression and resistance to E. ictaluri in P. hypophthalmus were explored. Fish were fed 0.1% fungal-derived β-glucan or 0.1% commercial yeast-derived β-glucan or a basal control diet without glucan. After 14 days of feeding, the mRNA expression of immune genes (transferrin, C-reactive protein, precerebellin-like protein, Complement C3 and factor B, 2a MHC class II and interleukin-1 beta) in liver, kidney and spleen were determined. Following this fish from each of the three diet treatment groups were infected with E. ictaluri and further gene expression measured 24 h post-infection (h.p.i.), while the remaining fish were monitored over 2 weeks for mortalities. Cumulative percentage mortality at 14 days post-infection (d.p.i.) was less in β-glucan fed fish compared to controls. There was no difference in gene expression between dietary groups after feeding for 14 days, but there was a clear difference between infected and uninfected fish at 24 h.p.i., and based on principal component analysis β-glucans stimulated the overall expression of immune genes in the liver, kidney and spleen at 24 h.p.i., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Salmon lice, Lepeophtheirus salmonis (Krøyer, 1837), are fish ectoparasites causing significant economic damage in the mariculture of Atlantic salmon, Salmo salar Linnaeus, 1758. The control of L. salmonis at fish farms relies to a large extent on treatment with anti-parasitic drugs. A problem related to chemical control is the potential for development of resistance, which in L. salmonis is documented for a number of drug classes including organophosphates, pyrethroids and avermectins. The ATP-binding cassette (ABC) gene superfamily is found in all biota and includes a range of drug efflux transporters that can confer drug resistance to cancers and pathogens. Furthermore, some ABC transporters are recognised to be involved in conferral of insecticide resistance. While a number of studies have investigated ABC transporters in L. salmonis, no systematic analysis of the ABC gene family exists for this species. This study presents a genome-wide survey of ABC genes in L. salmonis for which, ABC superfamily members were identified through homology searching of the L. salmonis genome. In addition, ABC proteins were identified in a reference transcriptome of the parasite generated by high-throughput RNA sequencing (RNA-seq) of a multi-stage RNA library. Searches of both genome and transcriptome allowed the identification of a total of 33 genes / transcripts coding for ABC proteins, of which 3 were represented only in the genome and 4 only in the transcriptome. Eighteen sequences were assigned to ABC subfamilies known to contain drug transporters, i.e. subfamilies B (4 sequences), C (11) and G (2). The results suggest that the ABC gene family of L. salmonis possesses fewer members than recorded for other arthropods. The present survey of the L. salmonis ABC gene superfamily will provide the basis for further research into potential roles of ABC transporters in the toxicity of salmon delousing agents and as potential mechanisms of drug resistance.

Background: The present study aimed to explore the impact of dietary docosahexaenoic acid (DHA) on aspects of the metabolism of Atlantic salmon (Salmo salar). The effects of diets containing increasing levels of DHA (1 g kg(-1), 3 g kg(-1), 6 g kg(-1), 10 g kg(-1) and 13 g kg(-1)) on the liver transcriptome of post-smolt salmon was examined to elucidate patterns of gene expression and responses of specific metabolic pathways. Total RNA was isolated from the liver of individual fish and analyzed using a custom gene expression 44K feature Atlantic salmon oligo-microarray., Results: The expression of up to 911 unique annotated genes was significantly affected by dietary DHA inclusion relative to a low DHA reference diet. Analysis of a total of 797 unique genes were found with a significant linear correlation between expression level and dietary DHA. Gene-Set Enrichment Analysis (GSEA) identified a range of pathways that were significantly affected by dietary DHA content., Conclusions: Pathways that showed a significant response to dietary DHA level included those for long-chain polyunsaturated fatty acid biosynthesis, fatty acid elongation, steroid biosynthesis, glycan biosynthesis, protein export and protein processing in the endoplasmic reticulum. These findings suggest that in addition to clear roles in influencing lipid metabolic pathways, DHA might also have key functional roles in other pathways distinct from lipid metabolism.

Background: European sea bass (Dicentrarchus labrax) is one of the most important farmed species in Mediterranean aquaculture. The observed sexual growth and maturity dimorphism in favour of females adds value towards deciphering the sex determination system of this species. Current knowledge indicates the existence of a polygenic sex determining determination system that interacts with temperature. This was explored by restriction-site associated DNA (RAD) marker analysis in a test panel of 175 offspring that originated from a factorial cross between two dams and four sires from a single full-sib family., Results: The first high-density single nucleotide polymorphism (SNP) based linkage map for sea bass was constructed, consisting of 6706 SNPs on 24 linkage groups. Indications for putative sex-determining QTL (quantitative trait loci) that were significant at the genome-wide threshold were detected on linkage groups 6, 11 and 18 to 21, although a genome-wide association study (GWAS) did not identify individual significant SNPs at a genome-wide threshold. A preliminary genomic prediction approach that tested the efficiency of SNP-based selection for female sea bass showed a slight advantage compared to traditional pedigree-based selection. However, when the same models were tested on the same animals for selection for greater length, a clear advantage of the SNP-based selection was observed., Conclusions: Overall, the results of this study provide additional support to the polygenic sex determination hypothesis in sea bass. In addition, identification of sex-ratio QTL may provide new opportunities for sex-ratio control in sea bass.

The aim of the present study was to generate an experimental model to characterize the nutrigenomic profile of a plant-derived nutritional stress. Atlantic salmon (Salmo salar) was used as the model species. The nutritional stress was induced by inclusion of dietary defatted soybean meal (SBM), as this ingredient had been previously demonstrated to induce enteropathy in the distal intestine and reduce growth in salmon. Triplicate groups of Atlantic salmon were fed concentrations of 0, 100, 200 and 300 g kg(-1) SBM for 12 weeks and reduced growth performance was used as the indicator of nutritional stress. The transcriptome was analyzed in two tissues, liver and distal intestine, with the hypothesis being that the liver transcriptome would be characterized by gene expression responses related to overall growth and health performance, whereas intestinal gene expression would be dominated by specific responses to SBM. A set of 133 genes was differentially expressed in liver including 44 genes in common with the intestinal response. The liver-specific response included up-regulation of genes involved in protein digestion, energy metabolism and immune functions, whereas genes in other metabolic pathways were generally anabolic and down-regulated. These responses may be more related to general nutritional stress than to SBM per se. The transcriptomic profile in the distal intestine was consistent with the enteritis response as described previously. This study provides a comprehensive report on the profiles of liver and distal intestine transcriptomes, specifically highlighting the role of the liver in fish undergoing SBM-induced nutritional stress., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Background: Fish species often exhibit significant sexual dimorphism for commercially important traits. Accordingly, the control of phenotypic sex, and in particular the production of monosex cultures, is of particular interest to the aquaculture industry. Sex determination in the widely farmed Nile tilapia (Oreochromis niloticus) is complex, involving genomic regions on at least three chromosomes (chromosomes 1, 3 and 23) and interacting in certain cases with elevated early rearing temperature as well. Thus, sex ratios may vary substantially from 50%., Results: This study focused on mapping sex-determining quantitative trait loci (QTL) in families with skewed sex ratios. These included four families that showed an excess of males (male ratio varied between 64% and 93%) when reared at standard temperature (28°C) and a fifth family in which an excess of males (96%) was observed when fry were reared at 36°C for ten days from first feeding. All the samples used in the current study were genotyped for two single-nucleotide polymorphisms (rs397507167 and rs397507165) located in the expected major sex-determining region in linkage group 1 (LG 1). The only misassigned individuals were phenotypic males with the expected female genotype, suggesting that those offspring had undergone sex-reversal with respect to the major sex-determining locus. We mapped SNPs identified from double digest Restriction-site Associated DNA (ddRAD) sequencing in these five families. Three genetic maps were constructed consisting of 641, 175 and 1,155 SNPs from the three largest families. QTL analyses provided evidence for a novel genome-wide significant QTL in LG 20. Evidence was also found for another sex-determining QTL in the fifth family, in the proximal region of LG 1., Conclusions: Overall, the results from this study suggest that these previously undetected QTLs are involved in sex determination in the Nile tilapia, causing sex reversal (masculinisation) with respect to the XX genotype at the major sex-determining locus in LG 1.

Background: Atlantic salmon have been subject to domestication for approximately ten generations, beginning in the early 1970s. This process of artificial selection will have created various genetic differences between wild and farmed stocks. Each year, hundreds of thousands of farmed fish escape into the wild. These escapees may interbreed with wild conspecifics raising concerns for both the fish-farming industry and fisheries managers. Thus, a better understanding of the interactions between domesticated and wild salmon is essential to the continued sustainability of the aquaculture industry and to the maintenance of healthy wild stocks., Results: We compared the transcriptomes of a wild Norwegian Atlantic salmon population (Figgjo) and a Norwegian farmed strain (Mowi) at two life stages: yolk sac fry and post first-feeding fry. The analysis employed 44 k oligo-microarrays to analyse gene expression of 36 farmed, wild and hybrid (farmed dam x wild sire) individuals reared under identical hatchery conditions. Although some of the transcriptional differences detected overlapped between sampling points, our results highlighted the importance of studying various life stages. Compared to the wild population, the Mowi strain displayed up-regulation in mRNA translation-related and down regulation in nervous and immune system -related pathways in the sac fry, whereas up-regulation of digestive and endocrine activities, carbohydrate, energy, amino acid and lipid metabolism and down-regulation of environmental information processing and immune system pathways were evident in the feeding fry. Differentially regulated pathways that were common among life stages generally belonged to environmental information processing and immune system functional groups. In addition, we found indications of strong maternal effects, reinforcing the importance of including reciprocal hybrids in the analysis., Conclusions: In agreement with previous studies we showed that domestication has caused changes in the transcriptome of wild Atlantic salmon and that many of the affected pathways are life-stage specific We highlighted the importance of reciprocal hybrids to the deconvolution of maternal/paternal effects and our data support the view that the genetic architecture of the strains studied highly influences the genes differentially expressed between wild and domesticated fish.

Concerns are being expressed recently over possible environmental effects of human pharmaceuticals. Although the likelihood of acute toxicity is low, the continuous discharge of pharmaceuticals into the aquatic environment means that sublethal effects on non-target organisms need to be seriously considered. One-year-old Atlantic salmon parr were exposed to 7.85±0.13μgL(-1) of the antidepressant drug Carbamazepine (CBZ) for five days to investigate changes of mRNA expression in the brain by means of a custom 17k Atlantic salmon cDNA microarray. The selected concentration is similar to upper levels that can be found in hospital and sewage treatment plant effluents. After treatment, 373 features were differently expressed with 26 showing up- or down-regulation of ≥2-fold (p≤0.05). Among the mRNAs showing the highest change were the pituitary hormones encoding features somatolactin, prolactin and somatotropin, or growth hormone. Functional enrichment and network analyses of up- and down-regulated genes showed that CBZ induced a highly different gene expression profile in comparison to untreated organisms. CBZ induced expression of essential genes of the focal adhesion and extracellular matrix - receptor interaction pathways most likely through integrin alpha-6 (itga6) activation. Negative regulation of apoptotic process, extracellular matrix organization and heme biosynthesis were the most enriched biological process related GO-terms, with the simultaneous enrichment of collagen and extracellular region related cellular component GO-terms, and extracellular matrix structural constituent, hormone activity and chromatin binding molecular function related GO-terms. These results show that relatively low doses of CBZ may affect brain physiology in exposed salmon parr, targeting similar processes as in human, indicating a high degree of conservation of targets of CBZ action. However, and since the mRNAs showing most changes in expression are critical for adaptation to different stressors and life history transitions in Atlantic salmon, more research should be undertaken to assess CBZ effects to avoid impairment of normal development and maintenance of natural populations., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Background: Dense single nucleotide polymorphism (SNP) genotyping arrays provide extensive information on polymorphic variation across the genome of species of interest. Such information can be used in studies of the genetic architecture of quantitative traits and to improve the accuracy of selection in breeding programs. In Atlantic salmon (Salmo salar), these goals are currently hampered by the lack of a high-density SNP genotyping platform. Therefore, the aim of the study was to develop and test a dense Atlantic salmon SNP array., Results: SNP discovery was performed using extensive deep sequencing of Reduced Representation (RR-Seq), Restriction site-Associated DNA (RAD-Seq) and mRNA (RNA-Seq) libraries derived from farmed and wild Atlantic salmon samples (n = 283) resulting in the discovery of > 400 K putative SNPs. An Affymetrix Axiom® myDesign Custom Array was created and tested on samples of animals of wild and farmed origin (n = 96) revealing a total of 132,033 polymorphic SNPs with high call rate, good cluster separation on the array and stable Mendelian inheritance in our sample. At least 38% of these SNPs are from transcribed genomic regions and therefore more likely to include functional variants. Linkage analysis utilising the lack of male recombination in salmonids allowed the mapping of 40,214 SNPs distributed across all 29 pairs of chromosomes, highlighting the extensive genome-wide coverage of the SNPs. An identity-by-state clustering analysis revealed that the array can clearly distinguish between fish of different origins, within and between farmed and wild populations. Finally, Y-chromosome-specific probes included on the array provide an accurate molecular genetic test for sex., Conclusions: This manuscript describes the first high-density SNP genotyping array for Atlantic salmon. This array will be publicly available and is likely to be used as a platform for high-resolution genetics research into traits of evolutionary and economic importance in salmonids and in aquaculture breeding programs via genomic selection.

The salmon louse (Lepeophtheirus salmonis (Krøyer, 1837)) is a parasitic copepod that can, if untreated, cause considerable damage to Atlantic salmon (Salmo salar Linnaeus, 1758) and incurs significant costs to the Atlantic salmon mariculture industry. Salmon lice are gonochoristic and normally show sex ratios close to 1:1. While this observation suggests that sex determination in salmon lice is genetic, with only minor environmental influences, the mechanism of sex determination in the salmon louse is unknown. This paper describes the identification of a sex-linked Single Nucleotide Polymorphism (SNP) marker, providing the first evidence for a genetic mechanism of sex determination in the salmon louse. Restriction site-associated DNA sequencing (RAD-seq) was used to isolate SNP markers in a laboratory-maintained salmon louse strain. A total of 85 million raw Illumina 100 base paired-end reads produced 281,838 unique RAD-tags across 24 unrelated individuals. RAD marker Lsa101901 showed complete association with phenotypic sex for all individuals analysed, being heterozygous in females and homozygous in males. Using an allele-specific PCR assay for genotyping, this SNP association pattern was further confirmed for three unrelated salmon louse strains, displaying complete association with phenotypic sex in a total of 96 genotyped individuals. The marker Lsa101901 was located in the coding region of the prohibitin-2 gene, which showed a sex-dependent differential expression, with mRNA levels determined by RT-qPCR about 1.8-fold higher in adult female than adult male salmon lice. This study's observations of a novel sex-linked SNP marker are consistent with sex determination in the salmon louse being genetic and following a female heterozygous system. Marker Lsa101901 provides a tool to determine the genetic sex of salmon lice, and could be useful in the development of control strategies.

Background: Atlantic halibut (Hippoglossus hippoglossus) is a high-value, niche market species for cold-water marine aquaculture. Production of monosex female stocks is desirable in commercial production since females grow faster and mature later than males. Understanding the sex determination mechanism and developing sex-associated markers will shorten the time for the development of monosex female production, thus decreasing the costs of farming., Results: Halibut juveniles were masculinised with 17 α-methyldihydrotestosterone (MDHT) and grown to maturity. Progeny groups from four treated males were reared and sexed. Two of these groups (n = 26 and 70) consisted of only females, while the other two (n = 30 and 71) contained balanced sex ratios (50% and 48% females respectively). DNA from parents and offspring from the two mixed-sex families were used as a template for Restriction-site Associated DNA (RAD) sequencing. The 648 million raw reads produced 90,105 unique RAD-tags. A linkage map was constructed based on 5703 Single Nucleotide Polymorphism (SNP) markers and 7 microsatellites consisting of 24 linkage groups, which corresponds to the number of chromosome pairs in this species. A major sex determining locus was mapped to linkage group 13 in both families. Assays for 10 SNPs with significant association with phenotypic sex were tested in both population data and in 3 additional families. Using a variety of machine-learning algorithms 97% correct classification could be obtained with the 3% of errors being phenotypic males predicted to be females., Conclusion: Altogether our findings support the hypothesis that the Atlantic halibut has an XX/XY sex determination system. Assays are described for sex-associated DNA markers developed from the RAD sequencing analysis to fast track progeny testing and implement monosex female halibut production for an immediate improvement in productivity. These should also help to speed up the inclusion of neomales derived from many families to maintain a larger effective population size and ensure long-term improvement through selective breeding.

Atlantic salmon (Salmo salar L.), a member of the family Salmonidae, is a totemic species of ecological and cultural significance that is also economically important in terms of both sports fisheries and aquaculture. These factors have promoted the continuous development of genomic resources for this species, furthering both fundamental and applied research. MicroRNAs (miRNA) are small endogenous non-coding RNA molecules that control spatial and temporal expression of targeted genes through post-transcriptional regulation. While miRNA have been characterised in detail for many other species, this is not yet the case for Atlantic salmon. To identify miRNAs from Atlantic salmon, we constructed whole fish miRNA libraries for 18 individual juveniles (fry, four months post hatch) and characterised them by Illumina high-throughput sequencing (total of 354,505,167 paired-ended reads). We report an extensive and partly novel repertoire of miRNA sequences, comprising 888 miRNA genes (547 unique mature miRNA sequences), quantify their expression levels in basal conditions, examine their homology to miRNAs from other species and identify their predicted target genes. We also identify the location and putative copy number of the miRNA genes in the draft Atlantic salmon reference genome sequence. The Atlantic salmon miRNAs experimentally identified in this study provide a robust large-scale resource for functional genome research in salmonids. There is an opportunity to explore the evolution of salmonid miRNAs following the relatively recent whole genome duplication event in salmonid species and to investigate the role of miRNAs in the regulation of gene expression in particular their contribution to variation in economically and ecologically important traits.

Sex in Oreochromis niloticus (Nile tilapia) is principally determined by an XX/XY locus but other genetic and environmental factors also influence sex ratio. Restriction Associated DNA (RAD) sequencing was used in two families derived from crossing XY males with females from an isogenic clonal line, in order to identify Single Nucleotide Polymorphisms (SNPs) and map the sex-determining region(s). We constructed a linkage map with 3,802 SNPs, which corresponded to 3,280 informative markers, and identified a major sex-determining region on linkage group 1, explaining nearly 96% of the phenotypic variance. This sex-determining region was mapped in a 2 cM interval, corresponding to approximately 1.2 Mb in the O. niloticus draft genome. In order to validate this, a diverse family (4 families; 96 individuals in total) and population (40 broodstock individuals) test panel were genotyped for five of the SNPs showing the highest association with phenotypic sex. From the expanded data set, SNPs Oni23063 and Oni28137 showed the highest association, which persisted both in the case of family and population data. Across the entire dataset all females were found to be homozygous for these two SNPs. Males were heterozygous, with the exception of five individuals in the population and two in the family dataset. These fish possessed the homozygous genotype expected of females. Progeny sex ratios (over 95% females) from two of the males with the "female" genotype indicated that they were neomales (XX males). Sex reversal induced by elevated temperature during sexual differentiation also resulted in phenotypic males with the "female" genotype. This study narrows down the region containing the main sex-determining locus, and provides genetic markers tightly linked to this locus, with an association that persisted across the population. These markers will be of use in refining the production of genetically male O. niloticus for aquaculture.

Background: Caligid copepods, also called sea lice, are fish ectoparasites, some species of which cause significant problems in the mariculture of salmon, where the annual cost of infection is in excess of €300 million globally. At present, caligid control on farms is mainly achieved using medicinal treatments. However, the continued use of a restricted number of medicine actives potentially favours the development of drug resistance. Here, we report transcriptional changes in a laboratory strain of the caligid Lepeophtheirus salmonis (Krøyer, 1837) that is moderately (~7-fold) resistant to the avermectin compound emamectin benzoate (EMB), a component of the anti-salmon louse agent SLICE® (Merck Animal Health)., Results: Suppression subtractive hybridisation (SSH) was used to enrich transcripts differentially expressed between EMB-resistant (PT) and drug-susceptible (S) laboratory strains of L. salmonis. SSH libraries were subjected to 454 sequencing. Further L. salmonis transcript sequences were available as expressed sequence tags (EST) from GenBank. Contiguous sequences were generated from both SSH and EST sequences and annotated. Transcriptional responses in PT and S salmon lice were investigated using custom 15 K oligonucleotide microarrays designed using the above sequence resources. In the absence of EMB exposure, 359 targets differed in transcript abundance between the two strains, these genes being enriched for functions such as calcium ion binding, chitin metabolism and muscle structure. γ-aminobutyric acid (GABA)-gated chloride channel (GABA-Cl) and neuronal acetylcholine receptor (nAChR) subunits showed significantly lower transcript levels in PT lice compared to S lice. Using RT-qPCR, the decrease in mRNA levels was estimated at ~1.4-fold for GABA-Cl and ~2.8-fold for nAChR. Salmon lice from the PT strain showed few transcriptional responses following acute exposure (1 or 3 h) to 200 μg L-1 of EMB, a drug concentration tolerated by PT lice, but toxic for S lice., Conclusions: Avermectins are believed to exert their toxicity to invertebrates through interaction with glutamate-gated and GABA-gated chloride channels. Further potential drug targets include other Cys-loop ion channels such as nAChR. The present study demonstrates decreased transcript abundances of GABA-Cl and nAChR subunits in EMB-resistant salmon lice, suggesting their involvement in avermectin toxicity in caligids.