Genetic fingerprinting of salmon louse (Lepeophtheirus salmonis) populations in the North-East Atlantic using a random forest classification approach

Estimating blue mussel (Mytilus edulis) connectivity and settlement capacity in mid-latitude fjord regions

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.

Aquaculture-driven evolution of the salmon louse mtDNA genome

Resistance toward the antiparasitic pyrethroid, deltamethrin, is reported in the Atlantic salmon louse (Lepeophtheirus salmonis salmonis), a persistent ectoparasite of farmed and wild salmonids. The resistance mechanism is linked to mitochondrial DNA (mtDNA), where genetic markers for resistance have been identified. Here, we investigated how widespread pyrethroid use in aquaculture may have influenced mtDNA variation in lice, and the dispersion of resistant haplotypes across the North Atlantic, using historical (2000-2002 "pre-resistance") and contemporary (2014-2017 "post-resistance") samples. To study this, we sequenced ATPase 6 and cytochrome b, genotyped two genetic markers for deltamethrin resistance, and genotyped microsatellites as "neutral" controls of potential population bottlenecks. Overall, we observed a modest reduction in mtDNA diversity in the period 2000-2017, but no reduction in microsatellite variation was observed. The reduction in mtDNA variation was especially distinct in two of the contemporary samples, fixed for one and two haplotypes, respectively. By contrast, all historical samples consisted of close to one mtDNA haplotype per individual. No population genetic structure was detected among the historical samples for mtDNA nor microsatellites. By contrast, significant population genetic differentiation was observed for mtDNA among some of the contemporary samples. However, the observed population genetic structure was tightly linked with the pattern of deltamethrin resistance, and we therefore conclude that it primarily reflects the transient mosaic of pyrethroid usage in time and space. Two historically undetected mtDNA haplotypes dominated in the contemporary samples, both of which were linked to deltamethrin resistance, demonstrating primarily two origins of deltamethrin resistance in the North Atlantic. Collectively, these data demonstrate that the widespread use of pyrethroids in commercial aquaculture has substantially altered the patterns of mtDNA diversity in lice across the North Atlantic, and that long-distance dispersion of resistance is rapid due to high level of genetic connectivity that is observed in this species.

Drug and pesticide usage for sea lice treatment in salmon aquaculture sites in a Canadian province from 2016 to 2019

We used 4 years of publicly available data (2016–2019) on chemical usage at salmon sites with information on production, stocking, locations and environmental conditions to explore patterns of anti-sea lice treatments in a Canadian province. Results show that sequential chemical treatments are prevalent, emamectin benzoate (EMB) with azamethiphos being the most used combination with a decrease in ivermectin usage. Relatively high rates of usage of EMB per fish biomass may point to potential lice resistance patterns with information needed on mechanisms and local populations. Lower or no chemical usage at some sites indicate less sea lice infestations likely influenced by localized site conditions (coves), and a lessened need for medication due to the usage of cleaner fish and possibly other non-chemical methods (not documented in this report). The year/climate influenced chemical input only in sites with higher treatment levels likely due to effects on sea lice growth and reproduction. Observed differences between years are warmer surface temperature in the fall, a higher freshwater input in spring, and stronger wind conditions for 2017 and 2018 with more medication usage for these two years. The lack of significant effect of site distances calculated in zones of influence based on 24 h potential connectivity patterns highlight the need to refine the resolution of hydrodynamic processes.

Community Parameters and Genome-Wide RAD-Seq Loci of Ceratothoa oestroides Imply Its Transfer between Farmed European Sea Bass and Wild Farm-Aggregating Fish

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.

A machine learning approach for the identification of population-informative markers from high-throughput genotyping data: application to several pig breeds

Single nucleotide polymorphisms (SNPs) able to describe population differences can be used for important applications in livestock, including breed assignment of individual animals, authentication of mono-breed products and parentage verification among several other applications. To identify the most discriminating SNPs among thousands of markers in the available commercial SNP chip tools, several methods have been used. Random forest (RF) is a machine learning technique that has been proposed for this purpose. In this study, we used RF to analyse PorcineSNP60 BeadChip array genotyping data obtained from a total of 2737 pigs of 7 Italian pig breeds (3 cosmopolitan-derived breeds: Italian Large White, Italian Duroc and Italian Landrace, and 4 autochthonous breeds: Apulo-Calabrese, Casertana, Cinta Senese and Nero Siciliano) to identify breed informative and reduced SNP panels using the mean decrease in the Gini Index and the Mean Decrease in Accuracy parameters with stability evaluation. Other reduced informative SNP panels were obtained using Delta, Fixation index and principal component analysis statistics, and their performances were compared with those obtained using the RF-defined panels using the RF classification method and its derived Out Of Bag rates and correct prediction proportions. Therefore, the performances of a total of six reduced panels were evaluated. The correct assignment of the animals to its breed was close to 100% for all tested approaches. Porcine chromosome 8 harboured the largest number of selected SNPs across all panels. Many SNPs were included in genomic regions in which previous studies identified signatures of selection or genes (e.g. ESR1, KITL and LCORL) that could contribute to explain, at least in part, phenotypically or economically relevant traits that might differentiate cosmopolitan and autochthonous pig breeds. Random forest used as preselection statistics highlighted informative SNPs that were not the same as those identified by other methods. This might be due to specific features of this machine learning methodology. It will be interesting to explore if the adaptation of RF methods for the identification of selection signature regions could be able to describe population-specific features that are not captured by other approaches.

Indications of strong adaptive population genetic structure in albacore tuna (Thunnus alalunga) in the southwest and central Pacific Ocean

Albacore tuna (Thunnus alalunga) has a distinctly complex life history in which juveniles and adults separate geographically but at times inhabit the same spaces sequentially. The species also migrates long distances and presumably experiences varied regimes of physical stress over a lifetime. There are, therefore, many opportunities for population structure to arise based on stochastic differences or environmental factors that promote local adaptation. However, with the extent of mobility consistently demonstrated by tagged individuals, there is also a strong argument for panmixia within an ocean basin. It is important to confirm such assumptions from a population genetics standpoint for this species in particular because albacore is one of the principal market tuna species that sustains massive global fisheries and yet is also a slow-growing temperate tuna. Consequently, we used 1,837 neutral SNP loci and 89 loci under potential selection to analyze population genetic structure among five sample groups collected from the western and central South Pacific. We found no evidence to challenge panmixia at neutral loci, but strong indications of structuring at adaptive loci. One population sample, from French Polynesia in 2004, was particularly differentiated. Unfortunately, the current study cannot infer whether the divergence is geographic or temporal, or possibly caused by sample distribution. We encourage future studies to include potentially adaptive loci and to continue fine-scale observations within an ocean basin, and not to assume genome-wide panmixia.

Genome-wide comparisons reveal a clinal species pattern within a holobenthic octopod-the Australian Southern blue-ringed octopus, Hapalochlaena maculosa (Cephalopoda: Octopodidae)

The southern blue-ringed octopus, Hapalochlaena maculosa (Hoyle, 1883) lacks a planktonic dispersal phase, yet ranges across Australia's southern coastline. This species' brief and holobenthic life history suggests gene flow might be limited, leaving distant populations prone to strong genetic divergence. This study used 17,523 genome-wide SNP loci to investigate genetic structuring and local adaptation patterns of H. maculosa among eight sampling sites along its reported range. Within sites, interrelatedness was very high, consistent with the limited dispersal of this taxon. However, inbreeding coefficients were proportionally lower among sites where substructuring was not detected, suggesting H. maculosa might possess a mechanism for inbreeding avoidance. Genetic divergence was extremely high among all sites, with the greatest divergence observed between both ends of the distribution, Fremantle, WA, and Stanley, TAS. Genetic distances closely followed an isolation by geographic distance pattern. Outlier analyses revealed distinct selection signatures at all sites, with the strongest divergence reported between Fremantle and the other Western Australian sites. Phylogenetic reconstructions using the described sister taxon H. fasciata (Hoyle, 1886) further supported that the genetic divergence between distal H. maculosa sites in this study was equivalent to that of between established heterospecifics within this genus. However, it is advocated that taxonomic delineations within this species should be made with caution. These data indicate that H. maculosa forms a clinal species pattern across its geographic range, with gene flow present through allele sharing between adjacent populations. Morphological investigations are recommended for a robust resolution of the taxonomic identity and ecotype boundaries of this species.