Single Nucleotide Polymorphisms Are Suitable For Assessing The Success of Restocking Efforts of The European Lobster (Homarus Gammarus, L.)

The European lobster (Homarus gammarus) forms the base of an important shery along the coasts of Europe. However, stocks have been in decline for many years, prompting new regulations in the shery and also restocking efforts. An important feature of any restocking effort is the assessment of success in the number of released juveniles that stay and become adult over time. Here, we tested the power of a SNP DNA marker panel developed for population assignment to correctly infer parentage on the maternal side of lobster larvae, in the absence of known fathers, using lobsters included in a current restocking effort on the Swedish west coast. We also examined the power to reconstruct the unknown paternal genotypes, and examined the number of fathers for each larval clutch. We found that the 96-SNP panel, despite only containing 78 informative markers, allowed us to assign all larvae to the correct mother. Furthermore, with 10 genotyped larvae or more, condent paternal genotypes could be reconstructed. We also found that 15 out of 17 clutches were full siblings, whereas two clutches had two fathers. To our knowledge, this is the rst time a SNP panel of this size has been used to assess parentage in a crustacean restocking effort. Our conclusion is that the panel works well, and that it could be an important tool for the assessment of restocking success of H. gammarus in the future.


Introduction
The world's oceans are continuously experiencing stronger and stronger shing pressure, and presently there are few sheries that are exploited less than to the maximum sustainable yield. In addition, about 30 % of the world's marine sh stocks are currently shed at biologically unsustainable levels (FAO 2020), a fraction which is currently increasing. As a consequence, a number of actions have been undertaken in order to minimize unsustainable practices, including implementation of no-take zones, shing quotas, size limits, along with stricter enforcement of regulations (FAO 2020). Many sheries have also introduced restocking efforts, where sh or shell sh are raised in aquaculture environments from eggs to juvenile life stages, and then released into the wild (reviewed in e.g. Halverson 2008; Taylor et al. 2017). In order to distinguish hatchery-reared individuals from wild ones in later life stages, a variety of different methods have been used, including different types of staining (Camp et al. 2013) and electronic tags (Nzau Matondo et al. 2019). However, the large-scale effects of restocking efforts have to date been di cult to assess, and in many small-scale efforts where evaluations have been possible the costs seem to outweigh the bene ts (Kitada 2018).
The European Lobster (Homarus gammarus) is the base for an important local shery all along the coast of Europe. Despite this, there are currently no EU-common regulations for the shery (Bryhn et al. 2020).
Along the Swedish west coast, the lobster shery is to a large extent recreational, and although there is no record of the number of landings from the recreational shery, it was estimated to account for almost 90 % of total Swedish used traps, and 75 % of total landings in 2007 (Bryhn et al. 2020 As a consequence, regulations for the recreational shery have been tightened in recent years in Sweden, with a lowered number of traps per person allowed, and higher minimum landing size limits. In addition to the stricter regulations, a restocking effort has recently been introduced on the Swedish west coast, where ovigerous females are supplied by local lobster shers, and the larvae are hatched in aquaculture facilities. After that, the lobster larvae are reared until ready to settle on the benthos, after which they are released at selected sites in the same area where their mothers were caught. Restocking has been taking place since 2018, with the expectation that the rst outplanted lobsters will have grown above the minimum landing size in the fall of 2022. Estimating the success of the restocking effort will be of crucial importance to the shery at that point. Here, we applied the SNP panel developed by Jenkins et al. (2019a) to a dataset consisting of mothers, progeny and unrelated control individuals, in order to test the power to: (1) Assign juveniles to the correct mother; (2) Not incorrectly assign unrelated individuals to a mother. As pregnant females collected for the restocking effort in the wild were used, there was no information on potential father(s), so therefore an additional goal was to assess the power of the panel to: (3) Con dently reconstruct paternal genotypes.
In addition, the panel was used to (4) Investigate whether clutches were originating from one or several fathers.

Materials And Methods
Lobster juveniles were reared at the Kristineberg Marine Research Station from wild-caught gravid females (n=17), until ca. 2 months old, during 2018-2020 (Table 1). The females were caught in the Sotefjord (Sotefjorden, SE582700-110451) and transported to Kristineberg by local lobster shermen.
They were then kept over winter, in 40 L perforated boxes in a seawater ow-through system. The water was pumped up from 32 m depth in the Gullmar fjord and varied in salinity and temperature with the season, between 30.5 -34.3 PSU and 5.4 -19.4 °C respectively. The females were fed fresh frozen shrimp (Pandalus borealis), and the embryo development in the egg clutches was checked routinely based on the eye index method developed by Perkins (1972). The females were brought over in pairs to a separate hatching system and placed in separate 40 L boxes with the water out ow leading to a larval collection tank. The temperature was incrementally raised to 18.5 ± 0.5 °C (1 °C per day) to speed up the embryonal development. Upon hatching, the larvae were counted and transferred to the larvae rearing set up, at which point larvae and tissue samples were collected.
From each family unit, a pleiopod from the mother was removed and placed in 95 % ethanol, and newly hatched larvae (n=10-30 per family) were placed whole in 95 % ethanol. Mothers were then re-released into the area where caught. Genomic DNA was then extracted from tissue dissected from the base of the pleiopods (mothers) or from whole larvae, using a Qiagen DNeasy Blood & Tissue kit following the standard protocol. Larval exoskeletons were fragmented in a tissue homogenizer (Retsch GmbH) before lysis. DNA quantity, purity and integrity were assessed using a Nanodrop ND-1000 spectrophotometer (Thermo Fisher), and through gel electrophoresis.
For maternal assignment and reconstruction of paternal genotypes, a previously published SNP panel of 96 loci was used (Jenkins et al. 2019a). Genotyping was performed at the SNP&SEQ platform at the Swedish National Genomics infrastructure (Uppsala University) through the iPLEX/MassArray system. All mothers were genotyped as two technical replicates, in order to assess the reliability of the genotype data. A total of 309 larvae from the 17 mothers were genotyped. Additionally, DNA from 10 non-related adult individuals, sampled 2007 at Kåvra on the Swedish west coast for a previous study (Huserbråten et al. 2013), was extracted and genotyped in order to test for putative false positive assignments. Four out of the 96 loci did not produce a viable assay (H_gam_15531, H_gam_4173, H_gam_31462, H_gam_35584), leaving 92 loci for genotyping. Seven loci had genotyping rates of 0 % (H_gam_11183,   H_gam_15109, H_gam_54240, H_gam_55142, H_gam_63798, H_gam_65064, H_gam_65376), and out of the well-genotyped loci, 7 were monomorphic (H_gam_25580, H_gam_58053, H_gam_32362, H_gam_53889, H_gam_32358, H_gam_8953, H_gam_22365), leaving 78 informative SNP markers.
The genotype data was analyzed using the command-line version of COLONY, Version 2.0.6.6, Build 20200830 (Wang & Santure 2009; Jones & Wang 2010), after pruning individual larvae with a 0 % genotyping rate (n=7), leaving 302 well-genotyped larvae plus the ten additional false assignment controls. The settings used were: Not updating allele frequency; dioecious species; no inbreeding; diploid species; polygamy for males, monogamy for females; no clone inference; full sibship size scaling; strong prior on sibship size, mean paternal and maternal sibship size 1.0; no known population allele frequencies; 1 run; run length 2; monitor by iteration, interval 100000; use full likelihood; high precision for full likelihood; probability of dad included 0; probability of mum included 0.95; paternity exclusion threshold 1; maternity exclusion threshold 5; no known dyads, sibships or exclusions. The last setting was used despite knowledge of known sibships, as the purpose was to test the power to detect sibships.

Results
Genotyping rates per individual and locus are given in Supplementary Tables S1-S2. Overall, the mean genotyping rate per individual was 90.2% (not including the 7 individuals with 0 % genotyping rate). The mean genotyping rate per locus was 95.7% (not including the 7 loci with 0 % genotyping rate). 27 out of the 78 polymorphic SNPs were signi cantly deviating from Hardy-Weinberg equilibrium (Online Resource 2), which was expected as the dataset does not consist of randomly mating individuals. All 17 technical replicates had identical genotypes as their counterparts at all loci with called genotypes (Online Resource 3). The full genotype dataset is given in Online Resource 4.

Maternal assignment
All larvae with genotyped mothers (100 %, n=302) could be correctly assigned to their mothers (Online Resource 5). For 14 out of the 17 mothers, larvae were assigned with >95% probability, while for mother S9 the probability was 93.5% and for mothers S15 and S17 the probabilities were 77.1% and 71.0%, respectively ( Table 1). The 10 genotyped control individuals without genotyped mothers were not assigned to any of the genotyped mothers. Rather, they were assigned to 6 different imputed maternal genotypes (ID #1-#6; Online Resource 5), indicating the possibility of some half-sibship among these individuals. However, the assignment probabilities of these individuals were generally low (Online Resource 5).

Imputation of paternal genotypes
The power to con dently impute paternal genotypes (P > 0.95) depended strongly on the number of genotyped larvae (p < 0.001), with a dramatic decrease below 9 larvae. With 20 larvae or more, essentially all paternal genotypes could be imputed ( Figure 1; Table 2). The paternal genotype imputation inferred 29 fathers for the whole dataset, of which 10 were assigned to the control lobsters with no mother genotypes present ( Table 2). For the remaining larvae, most were inferred to be full siblings. Only two mothers presented larvae from more than one father (M-S13 and M-S17; Table 1). All imputed paternal genotypes are given in Online Resource 6. In addition, the imputed maternal genotypes for the control lobsters are given in Online Resource 7.

Discussion
Despite the reduced number of SNPs that produced informative data (n=78 out of 96), we were here able to correctly assign all of our juvenile individuals to the correct mother, while not incorrectly assign any of the unrelated control individuals to any of the mothers in the dataset. Two of the mothers (S15 and S17) had lower assignment probabilities, despite not deviating in terms of genotyping rate, heterozygosity or genetic distance from others, leaving the reason unclear. Further, the success of maternal assignment was independent of the number of larvae sampled from each mother. SNP panels of similar size have successfully been used to infer parentage in restocking efforts of sh species such as bass (Zhao et al. 2018) and rainbow trout (Liu et al. 2016), and have also been shown to work well to infer parentage in selected shrimp lines in aquaculture (Sellars et al. 2014). However, this is to our knowledge the rst time a SNP panel has been used to infer motherhood in absence of known fathers in a restocking effort in any crustacean species to date. The results gained here suggest that using this SNP panel to assess parentage of lobsters caught in a restocking area would work well.
The reconstruction of the paternal genotypes, however, proved to be dependent on the number of genotyped larvae. This is not surprising, as the power to impute the paternal genotype depends on the probability of sampling an alternate allele in the larval pool. A potential confounding factor here could also be that the larvae coming from small sibling groups were also the ones with no mother present in the dataset, meaning that both maternal and paternal genotypes had to be imputed, further complicating the imputation process. Nevertheless, we here found that using ten larvae per family generate a con dent imputation of paternal genotypes. We also found that in 15 out of the 17 larval clutches examined (88 %), there was only one father, while in the remaining two clutches larvae originated from two fathers. Ellis et al. (2015a) found no multiple paternity in 34 analyzed broods, while Sørdalen et al. (2018) found 2 out of 54 broods (2 %) to have two fathers, both concluding that multiple paternal fertilization is highly unusual in H. gammarus. Here, we also nd that it is unusual, but that it does seem to occur occasionally.
Our results show that this SNP panel is a good choice for determining the success of the ongoing restocking effort on the west coast of Sweden, and possibly also in other areas in Europe. As all mothers used in the ongoing restocking effort have been sampled for DNA, it is likely that wild-caught lobsters from the fall of 2022 and onward will be able to be assigned to these mothers if they are indeed released ones, and that they will not assign to these mothers if they are wild-raised ones. Additionally, as the mothers were released back into the area, the genetic panel will be an important tool to investigate whether these individuals are re-caught in future efforts. We also conclude that as the restocking effort continues, pleiopods from all mothers used should be sampled and stored in ethanol, for future reference. Further, if there is a desire to also impute paternal genotypes from the ongoing efforts, for example to learn if the fathers of the lab-raised lobsters are being caught, then a minimum of ten larvae from each released clutch should also be preserved in ethanol.

Declarations
Funding: Funding for this research was provided by the project "Kläckning och uppfödning av utsättningsbara Hummeryngel" nanced by the Swedish board of agriculture and European Maritime and Fisheries Fund, the European Regional Development fund project "MarGen II -Kapacitetsopbygning for øget baeredygtig vaekst ved akvatiske miljøer i Kattegat-Skagerrak", and from the municipality of Sotenäs.
Competing interests: The authors declare no competing interests.
Data availability: All data generated in this project is provided in the supplementary materials.
Code availability: Not Applicable Authors' contributions: PDW and CA designed the study and performed genetic analyses. LS, IC, SPE and KS were responsible for the rearing of lobster juveniles. PDW wrote the manuscript. All co-authors contributed to the nal version of the manuscript.  Table 2. Statistics of imputed paternal genotypes ("#" denotes imputed mothers and "*" denotes imputed fathers).