Hatchery-Imposed Selection Does Not Impact the Genetic Diversity of Australian Farmed Blue Mussels ( Mytilus spp.)

: Australian blue mussels ( Mytilus spp.) are an increasingly important sustainable product of the Australian aquaculture industry. Although important for commercial fisheries, aquaculture may have adverse environmental and ecological impacts. This study assessed the impact of standard hatchery-imposed selection practices on the genetic diversity of farmed blue mussels. Using microsatellite markers, relatedness and genetic structure analyses showed that hatchery-reared larvae have high levels of genetic diversity without a significant decline as they move through the hatchery rearing process. Selection and/or genetic drift does appear to be operating during the hatchery rearing process, however, evidenced by an increase in relatedness among larvae over time. Significant shifts in allele frequency as well as genetic clusters provides further evidence that selection is acting on larvae due to the selection practice applied at the hatchery. Comparison of the level of genetic diversity and genetic differentiation of adults from wild and farmed populations provided no evidence that farmed mussels have lower diversity, or that they are genetically swamping local natural populations. The data suggest that careful design and implementation of mussel breeding programs can maintain high genetic diversity among larvae that does not lead to genetic swamping of natural mussel populations in the surrounding area.


Introduction
Aquaculture is one of the fastest growing food production industries in the world [1].It is an increasingly important sustainable food source, with marine bivalves production generating USD 29.8 billion worldwide [2].Bivalves (commonly separated as clam, oyster, mussel and scallop) account for approximately 14% of global marine seafood production with 89% coming from aquaculture [3].Asia is the largest contributor to bivalve production (85%) followed by the Americas (9%) and Europe (5.5%)[3].Africa and Oceania produce less than 1% of global bivalve populations [3].Although the fishery of marine bivalves is more limited in Oceania, mussels are one of the most important commercially produced bivalve species in this region [4].
Aquaculture plays a crucial role in providing food globally and this is likely to become more important as the world's population continues to grow.Aquaculture helps meet the increasing demand for seafood and helps reduce the impacts of harvesting wild populations, which is becoming progressively unsustainable [5].Unfortunately, some aspects of aquaculture can also have potential negative genetic impacts on natural populations if not carefully considered during the design and implementation of breeding programs.The choice of broodstock needs to be carefully designed in aquaculture as it can have important implications for the level of genetic diversity of animals produced.Many highly related individuals (full and half siblings) may be present under industrial scale husbandry, and subsequent selection in aquaculture systems may further erode diversity and increase relatedness [6].Thus, when animals are released back into the natural environment (either intentionally for restocking or growing out, or unintentionally through escapees), they have the potential to cause higher levels of inbreeding and genetic swamping in the wild (frequency of alleles from a farmed population to increase in the surrounding natural population) [7][8][9].For example, although selective breeding can have positive effects such as faster growth, better utilization of farming resources and disease resilience [10][11][12], there have been many cases where the release of cultured marine species resulted in genetic swamping when they hybridized with natural populations (e.g., the impact of cultured eastern oysters Crassostrea virginica [13] and European flat oysters Ostrea edulis [14] on their wild counterparts).The aquaculture process could also increase the frequency of maladaptive traits in natural populations.For example, multiple maladaptive traits have been linked to the release of salmon into natural populations including reduced developmental rate [15] and smaller egg size [16] resulting in higher mortality rates in the wild.
The taxonomic status of blue mussels around Australia remains unclear and is an active area of research [7,9,[17][18][19].Thus, in this study we use the term Mytilus spp.when talking about the study species.The Australian blue mussel (Mytilus sp.) forms the basis of an important aquaculture industry on the south coast of Australia.In 2017, the state of Victoria produced 45% of Australia's total mussel stock, worth AUD 5.2 m [20].In Victoria, the historical farming approach relied on the collection of wild spat from natural spawning grounds, which were then grown out on suspended ropes [21].However, due to high variability in spat supply from year to year, a commercial hatchery was established in 2008 to produce a more consistent and reliable spat supply for the mussel industry [21].In the commercial system, the adults are spawned, and the larvae are reared within a closed hatchery system.One potential consequence of the hatchery spawning process is that large numbers of highly related individuals are produced.This combined with selection within the hatchery (e.g., selection for larger, faster growing individuals) may favour certain family lineages and cause an overall decrease in the genetic diversity of farmed mussels [22].Once the larvae reach the settlement stage [21], they are settled onto ropes which are then hung from lines in the local environment for growing out to market size.It is during this phase that hatchery-produced and reared mussels have the potential to interbreed with surrounding natural populations.
To assess the potential impact of the hatchery process on levels of genetic diversity and relatedness of hatchery-produced larvae, microsatellite markers were used to conduct a genetic analysis of larvae throughout the hatchery rearing process.It was predicted that, due to the selection for larger larvae, genetic diversity will decrease, and the level of relatedness will increase among larvae as they move through the hatchery rearing process.Genetic diversity of mussels collected from farmed adult mussel lines in the natural environment were also established and compared to those from surrounding natural populations to identify if there was any evidence of hybridization of natural populations with farmed mussels that could lead to genetic swamping and loss of fitness in natural populations.

Sample Collection
Adult blue mussel broodstock for spawning within the hatchery were collected (N = 200) from a wild population in 2014 from the vicinity of Werribee town (Victoria, Australia 38.010 S, 144.680E) and transported back to the Victorian Shellfish Hatchery (VSH) facilities in Queenscliff.The adult mussels were spawned using standard hatchery protocols [21,23].All animals were cleaned of epiphytes and used for spawning within 24 h of collection.In total, 200 mussels were placed on a spawning table in approximately 6 cm of seawater and spawning was induced using thermal cycling.This involved progressively increasing seawater temperature from 16 to 24 • C using a water heating system over 20 min, then draining the water and repeating the procedure until the release of gametes was triggered.Males and females can readily be distinguished at the time of gamete release with females releasing orange gametes, while males release a white cloud of sperm.Individual mussels were rinsed with filtered seawater and placed into male or female communal 5 L spawning trays.Individuals were allowed to continue releasing gametes for up to 20 min.A total of 80 females were used for gamete collection, and 80 males used for sperm collection; the remaining 40 mussels failed to spawn.Egg concentrations were determined from three replicate counts using a 250 mL (5 mL eggs/245 mL sea water) sample with a coulter counter (560 micro aperture), while sperm concentrations were determined from three replicate haemocytometer counts.Fertilisations were carried out in a 5 L fertilisation chamber by adding enough eggs and sperm to achieve a final egg-sperm ratio of approximately 1:15 [21,24].Fertilisation was allowed to proceed for approximately 2 h, after which eggs were rinsed with filtered seawater onto a 30 µm mesh to remove excess sperm and prevent polyspermy.Fertilization was confirmed via microscopic examination before fertilised eggs were transferred to a 5000 L tank for larval development.
Larvae collection was done within the Queenscliff marine facility during a normal rearing run, thus larvae rearing followed VSH strict hatchery procedures as described above [21,23].In a normal rearing run, larvae are screened for size selection daily for 21 days and three sampling periods (beginning, middle, end) were chosen for analysis.Larvae were stocked into 12,400 L tanks at 16-18 million larvae per tank density and sampled on screening day 2, 12 and 20 following manual size grading.On each sampling day, larvae were washed with filtered seawater (1 µm, heated and UV-sterilised) through double filter screens at different grading mesh sizes (a 60 µm filter mesh screen on Day 2, 100 µm on Day 12 and 150 µm on Day 20 respectively) and preserved in 70% ethanol.The front filter screen removed unwanted debris.Larvae smaller than the filter mesh were considered un-fit or failed to survive the hatchery environment and discarded in a disposable container to reduce the risk of gene pollution.
A random selection of 96 mussels were collected from a wild population in Werribee (Victoria, Australia 38.010 S, 144.680E. The same location as broodstock collection) during March 2014 and 2020.A further 96 samples were purchased from mussel farms in the surrounding area that are supplied by VSH in 2020.Samples were transported live back to the facility where tissue biopsies were taken and preserved in 100% ethanol.

DNA Extraction and Microsatellite Genotyping
Individual mussel larvae were selected from the preserved larval samples under a dissection microscope using a micropipette.Each individual larva was rinsed by gently pipetting the larva into a fresh solution of 70% ethanol before being transferred into a 96-well plate (one larva per well).The PCR plate was sealed and pulse centrifuged to fix larvae to the bottom of the wells.Plates were unsealed and placed into a 56 • C incubator for approximately 30 min to remove any residual ethanol from the wells.DNA was extracted from larvae by adding 30 µL of Tween20 buffer to each well and then incubating at 56 • C for two hours followed by a 95 • C for 10 min.DNA extracts were stored at −20 • C until further use.
DNA extraction from adult tissue samples was conducted using a DNeasy Blood and Tissue Kit (Qiagen, Hilden, German) according to the manufacturer's instructions.Six microsatellite loci were selected to assess levels of genetic diversity, relatedness, and population structure.Microsatellite loci were amplified in two PCR multiplex reactions (see Supplementary Table S1).Multiplex 1 consisted of the loci MgU2-HEX, MT203-HEX and MGE005-FAM, and Multiplex 2 consisted of the loci Med744-HEX, My048-FAM and MgU3-HEX [25][26][27][28].For each multiplex, a 10 µL total reaction assay consisted of 0.2 mM dNTPs, 1 × PCR buffer, 0.25 mM MgCl 2 buffer, 0.5 U of QIAGEN Taq, 0.3 µM of each microsatellite primers and 30 ng DNA.Thermal cycling conditions for the two multiplexes varied slightly: for multiplex 1 conditions were 94 • C for 15 min, 30 cycles of 94 • C for 30 s, 55 • C for 3 s, 72 • C for 1 min.A final extension step was held at 72 • C for 15 min.Thermal cycling conditions for multiplex 2 were 94 • C for 15 min, 30 cycles of 94 • C for 30 s, 60 • C for 3 s, 72 • C for 1 min.A final extension step was held at 72 • C for 15 min.Amplified products were analysed on an ABI3730 (Applied Biosystems, Waltham, MA, USA) at the Australian Genome Facility (AGRF, Melbourne, Australia).Fragment sizes were determined with the LIZ500-250 standard using GeneMapper software (Applied Biosystems, Waltham, MA, USA).

Statistical Analysis
Null allele frequencies were estimated with the Van Oosterhout algorithm in MI-CROCHECKER v2.2.3 [29].One locus (My048) was excluded from further analysis due to low percentage of amplification.Null alleles were detected at 3 loci (MT203, MgU3, Med744).Analysis was performed with and without those loci with similar results obtained, therefore all results are presented for the 5 loci.Microsatellite data were tested for agreement with the Hardy Weinberg model in GENEPOP [30] using exact tests.All 5 loci showed significant deficits of heterozygotes for the larvae and adult samples; however, this is not surprising as a selection experiment was undertaken.
Levels of genetic diversity among larvae collected at each time point in the hatchery rearing processes were assessed as mean number of alleles and levels of expected (H E ) and observed (H O ) heterozygosity and were calculated using Arlequin v3.5.2.2 [31].Significant differences in genetic diversity between sampled time points were tested using a Friedman Chi-squared test for non-independent data, followed by pairwise paired Wilcoxon tests with false discovery rate correction for multiple testing [32,33] in R [34].To test for divergence in allele frequencies between the sampling time points, pairwise F ST values were calculated in Arlequin v3.5.2.2 [35].To further explore divergence in allele frequencies through the hatchery rearing process, a Bayesian model-based clustering analysis was used to determine the number of genetic clusters within our data set using the program STRUCTURE 2.3.4 [36].We assessed the number of genetic clusters (K) in our data for K ranging from 1 through to 10, with a burn-in length of 500,000, and with 500,000 MCMC iterations at each level and 10 replicates for each K. Genetic structuring was tested using the admixture model with allelic frequencies correlated among time points and ignoring prior population information [37].From the STRUCTURE output, the true number of clusters (K) was determined by looking at the optimum ∆K [38] using STRUCTURE HARVESTER [39].Relatedness among larvae sampled at each time point was calculated with COANCESTRY v1.0 [40], using the triadic likelihood method described by [41].This estimator was chosen as it is least biased with unrelated individuals, which is expected in wild populations [33].Significant changes in relatedness between time periods was assessed with 10,000 permutations in COANCESTRY [40].Inbreeding coefficients (F is ) were determined in GENEPOP (Rousset, 2008).
The same analyses were conducted on samples from farmed and adjacent wild adult mussels.Levels of genetic diversity between farmed and adjacent wild adult mussels was assessed as mean number of alleles and levels of expected heterozygosity (H E ), while population structure was calculated as F ST using the program Arlequin v3.5.2.2 [31].

Genetic Diversity of Hatchery-Raised Larvae
No significant change was observed in genetic diversity as larvae moved through the hatchery rearing process.Mean number of alleles was not significantly different among sampled time points (Friedman Chi-squared = 1.78; df = 2; p = 0.411).Expected heterozygosity showed a similar pattern across the sampled time points with no significant difference detected (Friedman Chi-squared = 0.4; df = 2; p = 0.819).Similarly, no significant difference was seen in observed heterozygosity (Friedman Chi-squared = 3.6; df = 2; p = 0.165).
The global F ST value indicated significant differentiation between larvae sampled at different time points in the hatchery rearing process (F ST = 0.0434, p < 0.0001).Pairwise F ST analysis revealed that this was driven by differences occurring between the Day 2 sampling period and the Day 12 and Day 20 sampling periods (Table 1).No significant genetic differentiation was observed between the Day 12 and Day 20 sampling periods (Table 1).A Bayesian STRUCTURE analyses identified three distinct genetic clusters corresponding to larvae sampled from each time point (Figure 1). 3. Results

Genetic Diversity of Hatchery-Raised Larvae
No significant change was observed in genetic diversity as larvae moved through the hatchery rearing process.Mean number of alleles was not significantly different among sampled time points (Friedman Chi-squared = 1.78; df = 2; p = 0.411).Expected heterozygosity showed a similar pa ern across the sampled time points with no significant difference detected (Friedman Chi-squared = 0.4; df = 2; p = 0.819).Similarly, no significant difference was seen in observed heterozygosity (Friedman Chi-squared = 3.6; df = 2; p = 0.165).
The global FST value indicated significant differentiation between larvae sampled at different time points in the hatchery rearing process (FST = 0.0434, p < 0.0001).Pairwise FST analysis revealed that this was driven by differences occurring between the Day 2 sampling period and the Day 12 and Day 20 sampling periods (Table 1).No significant genetic differentiation was observed between the Day 12 and Day 20 sampling periods (Table 1).A Bayesian STRUCTURE analyses identified three distinct genetic clusters corresponding to larvae sampled from each time point (Figure 1).
Estimates of relatedness among larvae within each time point also showed a significant trend of increasing relatedness through the hatchery rearing process (Figure 2, Table 2).Inbreeding coefficient values (FIS) were positive for time points Day 2 and Day 20, and negative for Day 12 (Table 2).Estimates of relatedness among larvae within each time point also showed a significant trend of increasing relatedness through the hatchery rearing process (Figure 2, Table 2).Inbreeding coefficient values (F IS ) were positive for time points Day 2 and Day 20, and negative for Day 12 (Table 2).

Genetic Diversity of Wild and Farmed Adult Mussel Populations
Estimates of genetic diversity including mean number of alleles (Friedman Chisquared = 0.44; df = 2; p = 0.8) and expected heterozygosity (Friedman Chi-squared = 3.6; df = 2; p = 0.17) did not differ between sampled wild and farmed populations.However, a significant drop in observed heterozygosity was identified in the wild population (Friedman Chi-squared = 8.4; df = 2; p = 0.015).
Significant levels of genetic differentiation were detected (FST = 0.0290, p < 0.0001) between wild and farmed populations with all pairwise FST values being significantly different (Table 1).Bayesian STRUCTURE analysis of the adult populations identified three

Genetic Diversity of Wild and Farmed Adult Mussel Populations
Estimates of genetic diversity including mean number of alleles (Friedman Chisquared = 0.44; df = 2; p = 0.8) and expected heterozygosity (Friedman Chi-squared = 3.6; df = 2; p = 0.17) did not differ between sampled wild and farmed populations.However, a significant drop in observed heterozygosity was identified in the wild population (Friedman Chi-squared = 8.4; df = 2; p = 0.015).
Significant levels of genetic differentiation were detected (F ST = 0.0290, p < 0.0001) between wild and farmed populations with all pairwise F ST values being significantly different (Table 1).Bayesian STRUCTURE analysis of the adult populations identified three potential genetic clusters within the dataset (K = 3), however the structure plot did not indicate that these were strongly aligned with source populations (Figure 3).
Hydrobiology 2024, 2, FOR PEER REVIEW 7 potential genetic clusters within the dataset (K = 3), however the structure plot did not indicate that these were strongly aligned with source populations (Figure 3).Relatedness estimates among individuals within each population were low (R ranged from 0.083 to 0.092) and not significantly different between populations (Figure 4, Table 2).Inbreeding coefficient values (FIS) were positive for all adult populations (Table 2).Relatedness estimates among individuals within each population were low (R ranged from 0.083 to 0.092) and not significantly different between populations (Figure 4, Table 2).Inbreeding coefficient values (F IS ) were positive for all adult populations (Table 2).

Genetic Diversity of Hatchery-Raised Larvae Compared to Adult Farmed and Wild Mussel Populations
Comparison of the wild and farmed populations with the hatchery-reared larvae showed that the hatchery larvae had significantly lower levels of expected heterozygosity (Friedman Chi-squared = 17.93; df = 5; p = 0.003) and mean number of alleles (Friedman Chi-squared = 17.93, p = 0.003) compared with the wild and farmed populations.The farmed and wild populations appear to have higher expected heterozygosity than the lar vae (Friedman Chi-squared = 19.29;df = 5; p = 0.0017).Levels of observed heterozygosity showed a similar pa ern to that seen for expected heterozygosity with the farmed and wild populations having a higher HO (Friedman Chi-squared = 12.31; df = 5; p = 0.031) Levels of genetic differentiation between larvae and adult populations were significantly different from zero (FST = 0.0782, p < 0.0001).Pairwise comparisons showed FST values be ing significantly different from zero between all populations (Table 1).Relatedness value of the farmed and wild populations were significantly lower than observed in larvae (Ta ble 2).

Discussion
The results from this study show that larvae do not appear to lose significant level of genetic diversity as they move through the hatchery rearing process, although selection does appear to be operating, resulting in increased relatedness among larvae.Larvae from the hatchery do appear to have lower levels of genetic diversity compared to wild and farmed adult populations of blue mussels; however, this does not appear to impact overal levels of genetic diversity and relatedness in wild and farmed adult populations.This i promising, as lower genetic diversity in adult populations can have potential fitness im pacts, including increased susceptibility to pathogens such as viruses, parasites and can cer [42,43] as well as reduced reproductive success [44].While these pa erns are positive further studies using more molecular markers would be necessary to confirm these re sults.The latest development of technologies (e.g., single cell sequencing) may be able to overcome the limitations of our study in which, due to a low amount of DNA, we could only use five microsatellite markers [45].
While no change in genetic diversity between time points was observed (indicating a

Genetic Diversity of Hatchery-Raised Larvae Compared to Adult Farmed and Wild Mussel Populations
Comparison of the wild and farmed populations with the hatchery-reared larvae showed that the hatchery larvae had significantly lower levels of expected heterozygosity (Friedman Chi-squared = 17.93; df = 5; p = 0.003) and mean number of alleles (Friedman Chi-squared = 17.93, p = 0.003) compared with the wild and farmed populations.The farmed and wild populations appear to have higher expected heterozygosity than the larvae (Friedman Chi-squared = 19.29;df = 5; p = 0.0017).Levels of observed heterozygosity showed a similar pattern to that seen for expected heterozygosity with the farmed and wild populations having a higher H O (Friedman Chi-squared = 12.31; df = 5; p = 0.031).Levels of genetic differentiation between larvae and adult populations were significantly different from zero (F ST = 0.0782, p < 0.0001).Pairwise comparisons showed F ST values being significantly different from zero between all populations (Table 1).Relatedness values of the farmed and wild populations were significantly lower than observed in larvae (Table 2).

Discussion
The results from this study show that larvae do not appear to lose significant levels of genetic diversity as they move through the hatchery rearing process, although selection does appear to be operating, resulting in increased relatedness among larvae.Larvae from the hatchery do appear to have lower levels of genetic diversity compared to wild and farmed adult populations of blue mussels; however, this does not appear to impact overall levels of genetic diversity and relatedness in wild and farmed adult populations.This is promising, as lower genetic diversity in adult populations can have potential fitness impacts, including increased susceptibility to pathogens such as viruses, parasites and cancer [42,43] as well as reduced reproductive success [44].While these patterns are positive, further studies using more molecular markers would be necessary to confirm these results.The latest development of technologies (e.g., single cell sequencing) may be able to overcome the limitations of our study in which, due to a low amount of DNA, we could only use five microsatellite markers [45].
While no change in genetic diversity between time points was observed (indicating a high standing genetic variation at the start of larvae raising), the level of relatedness increased while moving through the rearing process.Interestingly, the inbreeding coefficient (F IS ) drops from 0.350 for day 2 larvae to effectively zero for day 12 and day 20 larvae.This suggests a decrease in the heterozygous deficits detected at day 2, which may result from a decline in more homozygous larvae (family lineage).As it is a common practice for industries, including the Queenscliff Marine facility, to control their stock by selecting for desirable traits (e.g., size, growth rate), the observed increased relatedness may be the result of this selection process.This suggests that genetic drift and/or selection for familial linages may occur due to the artificial selection in the hatchery [46].For example, while the industry is selecting for larger phenotypes, these may not have underlying genetic determinants, as such the genotypes are randomly selected (i.e., genetic drift).It could also be hypothesized that the larger phenotypes are due to particular genotypes that are favoured by the selection process (i.e., selection for familial linages).
Although selecting for phenotypic traits may be beneficial for the industry, in the long term, the practice could be detrimental if it reduces overall genetic diversity.Currently, the process appears beneficial as it achieves faster growing larvae (less production time) while retaining good levels of genetic diversity.If this practice is not well maintained, however, it could have negative impacts, including higher levels of inbreeding.The observed retention of genetic diversity demonstrates the importance of rotating and introducing a diversity of broodstock between different spawning runs, the practice applied at this hatchery.Despite rotating parents in some practices, however, inbreeding can still remain (or even increase) in mussel populations, as shown by [47].The results showed similar patterns; however, this is not surprising given the many highly related individuals retained within the initial spawned larvae.While relatedness continues to increase during the rearing process, the inbreeding coefficient drops.This suggests that there is a selection for increased heterozygosity among larvae, which is known to have fitness benefits [48] and could be beneficial for both the hatchery and wild populations.
Comparison of wild and farmed adults provided no evidence that farming practices would result in a decline in genetic diversity of farmed mussels, and minimal evidence they are genetically swamping local natural populations.This was surprising as in many other bivalve species, including silver lipped pearl oysters (Pinctada maxima), eastern oyster (C.virginica), pacific oyster (C.gigas), Chinese freshwater pearl mussel (Hyriopsis cumingii) and the blue mussel (Mytilus edulis) at other locations, the artificial breeding process resulted in a significant decrease of genetic diversity (Ho, He and mean number of alleles) in the farmed populations when compared to natural ones [28,[49][50][51][52].Although the STRUCTURE analysis does not suggest any major genetic clusters relating to the populations, there is still weak (and significant) genetic differentiation based on the FST analysis but not enough for structure to distinguish different genetic clusters.This suggests that there is sufficient gene flow to prevent large amounts of genetic structuring.Levels of genetic diversity of farmed and wild populations appear high and this is likely resulting from the hatchery using a diverse set of broodstock that is changed between subsequent spawning.The use of different broodstock appears to maintain levels of diversity of the farmed mussels and there does not appear to be any evidence of genetic swamping in the wild populations.Interestingly, estimates of the inbreeding coefficient (F IS values) were higher in the natural compared to the farmed populations.In addition, when comparing the levels of heterozygosity of wild populations between 2014 and 2020, our results indicated a decline in heterozygotes in 2020 compared with 2014, potentially indicating loss of diversity in the wild populations, which may indicate some inbreeding in wild populations.
It is important to note that the farmed and wild populations reside in the same environment and hence could potentially exchange genetic materials.The observed genetic clustering in the STRUCTURE analyses and the significant level of genetic differentiation, however, suggests limited interbreeding between farmed and wild populations.Hatchery stocking can be beneficial and/or harmful for native populations.Breeding programs have been known to support the restoration of endangered species [53] and the release of some species into the wild has been seen to positively increase the number of spawning sites [54].Unfortunately, the release of hatchery stock into the wild have also had negative effects, including but not limited to lower reproductive success, behavioural changes and higher vulnerability to predation [54].Release of hatchery stock can also alter disease and parasite outcomes in wild populations, as well as cause hybridisation events.Farming involves a high density of individuals which can increase disease and parasite prevalence and risk, as seen in multiple salmonid species and their ectoparasites, sea lice [55].The use of farmed populations placed in a wild environment has also been known to decrease disease prevalence if handled correctly [56].For example, intensive oyster aquaculture (Crassostrea virginica) has been shown to reduce disease impacts by the protozoan parasite Perkinsus marinus on sympatric wild oysters [56].
The introduction of farmed animals into wild populations has been known to result in hybridisation [57], leading to genetic and phenotypic changes including epigenetic changes (DNA methylation) and the formation of maladaptive traits as seen in many salmonids [58][59][60] and mussels [57].As mentioned above, limited or no gene flow was detected between the wild and farmed populations in our study.The presence of these independent clusters may be due to the geographic distance between farmed and natural populations within the bay, or that the time to harvest does not allow for the hatchery populations to complete a full spawning cycle, as they are usually harvested at this time [61,62].
The observed low relatedness levels in the adult populations are also confirmed by the high genetic diversity of both the farmed and the natural populations, as all three populations (wild and farmed) contained mainly unrelated individuals, which decreases the risk of adverse effects of inbreeding.Interestingly, the farmed populations had a lower inbreeding coefficient than the wild ones.This may be the result of the hatchery practice where new specimens are used in each induced spawning run, while spawning in wild populations occurs en masse.Mass spawning is unpredictable and can result in the domination of certain family lineages [13,52].Induced spawning has been seen to improve reproductive success and decrease the risk of inbreeding and genetic bottleneck [13,52,63,64].

Management Implications and Conclusions
This study demonstrated the need for well-designed and carefully implemented mussel breeding programs to ensure that high levels of genetic diversity are maintained, and that genetic swamping of natural populations does not occur because of aquaculture practices.It is important to consider both the genetic implications of the larvae reared and how the farmed populations may interact and influence natural populations.The current evidence suggests that the larvae raising process at the hatchery where the study was conducted does not appear to have any associated negative effects on the natural populations or on stocks at the farms associated with their breeding program, as high levels of genetic diversity are maintained among their farmed mussels.The development of standardised breeding and rearing protocols is encouraged to ensure best practices that take into consideration the need to maintain high levels of genetic diversity of brood stocks, reared larvae and farmed adult populations.This not only provides increased resilience to future impacts such as climate change or the emergence of new pathogens, but also reduces the potential for adversely effecting and swamping natural populations.Further research should focus on investigating the effects of manual size selection and its effects on larvae phenotypic traits (growth rates, shell, and meat quality) in relation to their genotypes and genetic diversity.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/hydrobiology3010004/s1,Supplementary Table S1: A summary of all microsatellite markers used within this study.Markers were obtained from [25].

Figure 1 .Figure 1 .
Figure 1.The population genetic STRUCTURE analysis of the three larval sampling time periods (Day 2, Day 12 and Day 20) of Australian blue mussels (Mytilus spp.).Each vertical line represents a single individual and colours indicate deferred ancestry.Bar plot of K = 3, which is the most probable value based on deltaK methods.Image from CLUMPAK Server.Day 2 Day 12 Day 20 Figure 1.The population genetic STRUCTURE analysis of the three larval sampling time periods (Day 2, Day 12 and Day 20) of Australian blue mussels (Mytilus spp.).Each vertical line represents a single individual and colours indicate deferred ancestry.Bar plot of K = 3, which is the most probable value based on deltaK methods.Image from CLUMPAK Server.

Figure 2 .
Figure 2. Levels of relatedness among blue mussel larvae (Mytilus spp.) reared from three time points in the hatchery rearing process.Error bars represent 95% confidence interval.N = 96 at each time point.

Figure 2 .
Figure 2. Levels of relatedness among blue mussel larvae (Mytilus spp.) reared from three time points in the hatchery rearing process.Error bars represent 95% confidence interval.N = 96 at each time point.

Figure 3 .
Figure 3.The population genetic STRUCTURE analysis of three adult Australian blue mussels (Mytilus spp.) (Wild 2014, Farmed 2020, Wild 2020).Each vertical line represents a single individual and colours indicate deferred ancestry.Bar plot of K = 3, which is the most probable value based on deltaK methods.Image from CLUMPAK Server.

Figure 3 .
Figure 3.The population genetic STRUCTURE analysis of three adult Australian blue mussels (Mytilus spp.) (Wild 2014, Farmed 2020, Wild 2020).Each vertical line represents a single individual and colours indicate deferred ancestry.Bar plot of K = 3, which is the most probable value based on deltaK methods.Image from CLUMPAK Server.

Table 1 .
Pairwise F ST between mussel larvae sampled at different times in the hatchery rearing process and adult populations.Significance is indicated as ** p < 0.01, *** p < 0.0001.

Table 2 .
Summary of genetic diversity measures for larvae and adult mussels collected within this study.N = sample size, A x − = mean number of alleles, H O = observed heterozygosity, H E = expected heterozygosity, F IS = inbreeding coefficient.* indicates significant difference at 95% confidence.

Table 1 .
Pairwise FST between mussel larvae sampled at different times in the hatchery rearing process and adult populations.Significance is indicated as ** p < 0.01, *** p < 0.0001.

Table 2 .
Summary of genetic diversity measures for larvae and adult mussels collected within this study.N = sample size, Ax̄ = mean number of alleles, HO = observed heterozygosity, HE = expected heterozygosity, FIS = inbreeding coefficient.* indicates significant difference at 95% confidence.