SNP discovery and genetic structure in blue mussel species using low coverage sequencing and a medium density 60 K SNP-array.

Blue mussels from the genus Mytilus are an abundant component of the benthic community, found in the high latitude habitats. These foundation species are relevant to the aquaculture industry, with over 2 million tonnes produced globally each year. Mussels withstand a wide range of environmental conditions and species from the Mytilus edulis complex readily hybridize in regions where their distributions overlap. Significant effort has been made to investigate the consequences of environmental stress on mussel physiology, reproductive isolation, and local adaptation. Yet our understanding on the genomic mechanisms underlying such processes remains limited. In this study, we developed a multi species medium-density 60 K SNP-array including four species of the Mytilus genus. SNPs included in the platform were called from 138 mussels from 23 globally distributed mussel populations, sequenced using a whole-genome low coverage approach. The array contains polymorphic SNPs which capture the genetic diversity present in mussel populations thriving across a gradient of environmental conditions (~59 K SNPs) and a set of published and validated SNPs informative for species identification and for diagnosis of transmissible cancer (610 SNPs). The array will allow the consistent genotyping of individuals, facilitating the investigation of ecological and evolutionary processes in these taxa. The applications of this array extend to shellfish aquaculture, contributing to the optimization of this industry via genomic selection of blue mussels, parentage assignment, inbreeding assessment and traceability. Further applications such as genome wide association studies (GWAS) for key production traits and those related to environmental resilience are especially relevant to safeguard aquaculture production under climate change.

Two parallel chromosome-level reference genomes to support restoration and aquaculture of European flat oyster Ostrea edulis.

This volume of Evolutionary Applications sees the publication of two genomes for the European native flat oyster Ostrea edulis, a species of significant evolutionary, ecological and commercial value. Each is a highly contiguous chromosome-level assembly from individuals of different genetic backgrounds, which have been benchmarked against one another. This situation has resulted from the serendipitous discovery that two independent research groups were both deep into the process of building, annotating and investigating separately produced assemblies. Due to constraints with funder requirements and the need to recognize early career researchers for their work, alongside the technical challenge of integrating assemblies from two very different genomes, there was limited capacity to merge the sequences into one publication at the stage of discovery. This issue is likely to become very common over the next few years until the technologies for working with multiple genomes at once, for example, graph genomes, become commonplace in nonmodel species. Consequently, both of our teams have decided to collaborate rather than compete, recognizing the benefit to copublishing two separate genome resources for the research community, each with distinct scientific investigations, and working collaboratively to benchmark the assemblies.

Chromosome-level reference genome for European flat oyster (Ostrea edulis L.).

The European flat oyster (Ostrea edulis L.) is a bivalve naturally distributed across Europe, which was an integral part of human diets for centuries, until anthropogenic activities and disease outbreaks severely reduced wild populations. Despite a growing interest in genetic applications to support population management and aquaculture, a reference genome for this species is lacking to date. Here, we report a chromosome-level assembly and annotation for the European Flat oyster genome, generated using Oxford Nanopore, Illumina, Dovetail OmniC™ proximity ligation and RNA sequencing. A contig assembly (N50: 2.38 Mb) was scaffolded into the expected karyotype of 10 pseudochromosomes. The final assembly is 935.13 Mb, with a scaffold-N50 of 95.56 Mb, with a predicted repeat landscape dominated by unclassified elements specific to O. edulis. The assembly was verified for accuracy and completeness using multiple approaches, including a novel linkage map built with ddRAD-Seq technology, comprising 4016 SNPs from four full-sib families (eight parents and 163 F1 offspring). Annotation of the genome integrating multitissue transcriptome data, comparative protein evidence and ab-initio gene prediction identified 35,699 protein-coding genes. Chromosome-level synteny was demonstrated against multiple high-quality bivalve genome assemblies, including an O. edulis genome generated independently for a French O. edulis individual. Comparative genomics was used to characterize gene family expansions during Ostrea evolution that potentially facilitated adaptation. This new reference genome for European flat oyster will enable high-resolution genomics in support of conservation and aquaculture initiatives, and improves our understanding of bivalve genome evolution.

A single genomic region involving a putative chromosome rearrangement in flat oyster (Ostrea edulis) is associated with differential host resilience to the parasite Bonamia ostreae.

European flat oyster (Ostrea edulis) is an ecologically and economically important marine bivalve, that has been severely affected by the intracellular parasite Bonamia ostreae. In this study, a flat oyster SNP array (~14,000 SNPs) was used to validate previously reported outlier loci for divergent selection associated with B. ostreae exposure in the Northeast Atlantic Area. A total of 134 wild and hatchery individuals from the North Sea, collected in naïve (NV) and long-term affected (LTA) areas, were analysed. Genetic diversity and differentiation were related to the sampling origin (wild vs. hatchery) when using neutral markers, and to bonamiosis status (NV vs. LTA) when using outlier loci for divergent selection. Two genetic clusters appeared intermingled in all sampling locations when using outlier loci, and their frequency was associated with their bonamiosis status. When both clusters were compared, outlier data sets showed high genetic divergence (F ST > 0.25) unlike neutral loci (F ST not ≠ 0). Moreover, the cluster associated with LTA samples showed much higher genetic diversity and significant heterozygote excess with outlier loci, but not with neutral data. Most outliers mapped on chromosome 8 (OE-C8) of the flat oyster genome, supporting a main genomic region underlying resilience to bonamiosis. Furthermore, differentially expressed genes previously reported between NV and LTA strains showed higher mapping density on OE-C8. A range of relevant immune functions were specifically enriched among genes annotated on OE-C8, providing hypotheses for resilience mechanisms to an intracellular parasite. The results suggest that marker-assisted selection could be applied to breed resilient strains of O. edulis to bonamiosis, if lower parasite load and/or higher viability of the LTA genetic cluster following B. ostreae infection is demonstrated.

Genome-Wide Association and Genomic Prediction of Growth Traits in the European Flat Oyster (Ostrea edulis).

The European flat oyster (Ostrea edulis) is a bivalve mollusc that was once widely distributed across Europe and represented an important food resource for humans for centuries. Populations of O. edulis experienced a severe decline across their biogeographic range mainly due to overexploitation and disease outbreaks. To restore the economic and ecological benefits of European flat oyster populations, extensive protection and restoration efforts are in place within Europe. In line with the increasing interest in supporting restoration and oyster farming through the breeding of stocks with enhanced performance, the present study aimed to evaluate the potential of genomic selection for improving growth traits in a European flat oyster population obtained from successive mass-spawning events. Four growth-related traits were evaluated: total weight (TW), shell height (SH), shell width (SW) and shell length (SL). The heritability of the growth traits was in the low-moderate range, with estimates of 0.45, 0.37, 0.22, and 0.32 for TW, SH, SW and SL, respectively. A genome-wide association analysis revealed a largely polygenic architecture for the four growth traits, with two distinct QTLs detected on chromosome 4. To investigate whether genomic selection can be implemented in flat oyster breeding at a reduced cost, the utility of low-density SNP panels was assessed. Genomic prediction accuracies using the full density panel were high (> 0.83 for all traits). The evaluation of the effect of reducing the number of markers used to predict genomic breeding values revealed that similar selection accuracies could be achieved for all traits with 2K SNPs as for a full panel containing 4,577 SNPs. Only slight reductions in accuracies were observed at the lowest SNP density tested (i.e., 100 SNPs), likely due to a high relatedness between individuals being included in the training and validation sets during cross-validation. Overall, our results suggest that the genetic improvement of growth traits in oysters is feasible. Nevertheless, and although low-density SNP panels appear as a promising strategy for applying GS at a reduced cost, additional populations with different degrees of genetic relatedness should be assessed to derive estimates of prediction accuracies to be expected in practical breeding programmes.

Evaluation of interactive effects of phosphorus-32 and copper on marine and freshwater bivalve mollusks.

PURPOSE: Contaminants seldom occur in isolation in the aquatic environment. While pollution of coastal and inland water bodies has received considerable attention to date, there is limited information on potential interactive effects between radionuclides and metals. Whether by accidental or controlled release, such contaminants co-exist in aquatic ecosystems and can pose an enhanced threat to biota. Using a range of biological responses, the study aimed to evaluate relative interactive effects on representative freshwater and marine bivalve species. METHODS: An integrated, multi-biomarker approach was adopted to investigate response to copper (Cu, 18 µg L-1), a known environmentally relevant genotoxic metal and differing concentrations of phosphorus-32 (32P; 0.1 and 1 mGy d-1), alone and in combination in marine (Mytilus galloprovincialis) and freshwater (Dreissena polymorpha) mussels. Genetic and molecular biomarkers were determined post-exposure and included DNA damage (as measured by the comet assay), micronuclei (MN) formation, γ-H2AX foci induction and the expression of key stress-related genes (i.e. hsp70/90, sod, cat, gst). RESULTS: Overall, using a tissue-specific (i.e. gill and digestive gland) approach, genotoxic response was reflective of exposures where Cu had a slight additive effect on 32P-induced damage across the species (but not all), cell types and dose rates. Multivariate analysis found significant correlations between comet and γ-H2AX assays, across both the tissues. Transcriptional expression of selected genes were generally unaltered in response to contaminant exposures, independent of species or tissues. CONCLUSIONS: Our study is the first to explore the interactive effects of ionizing radiation (IR) and Cu on two bivalve species representing two ecological habitats. The complexity of IR-metal interactions demonstrate that extrapolation of findings obtained from single stressor studies into field conditions could be misrepresentative of real-world environments. In turn, environmental protective strategies deemed suitable in protecting biota from a single, isolated stressor may not be wholly adequate.

Identification and Full Characterisation of Two Novel Crustacean Infecting Members of the Family Nudiviridae Provides Support for Two Subfamilies.

Multiple enveloped viruses with rod-shaped nucleocapsids have been described, infecting the epithelial cell nuclei within the hepatopancreas tubules of crustaceans. These bacilliform viruses share the ultrastructural characteristics of nudiviruses, a specific clade of viruses infecting arthropods. Using histology, electron microscopy and high throughput sequencing, we characterise two further bacilliform viruses from aquatic hosts, the brown shrimp (Crangon crangon) and the European shore crab (Carcinus maenas). We assembled the full double stranded, circular DNA genome sequences of these viruses (~113 and 132 kbp, respectively). Comparative genomics and phylogenetic analyses confirm that both belong within the family Nudiviridae but in separate clades representing nudiviruses found in freshwater and marine environments. We show that the three thymidine kinase (tk) genes present in all sequenced nudivirus genomes, thus far, were absent in the Crangon crangon nudivirus, suggesting there are twenty-eight core genes shared by all nudiviruses. Furthermore, the phylogenetic data no longer support the subdivision of the family Nudiviridae into four genera (Alphanudivirus to Deltanudivirus), as recently adopted by the International Committee on Taxonomy of Viruses (ICTV), but rather shows two main branches of the family that are further subdivided. Our data support a recent proposal to create two subfamilies within the family Nudiviridae, each subdivided into several genera.

Ancestral Physical Stress and Later Immune Gene Family Expansions Shaped Bivalve Mollusc Evolution.

Bivalve molluscs comprise 20,000 species occupying a wide diversity of marine habitats. As filter feeders and detritivores they act as ecosystem engineers clarifying water, creating reefs, and protecting coastlines. The global decline of natural oyster reefs has led to increased restoration efforts in recent years. Bivalves also play an important role in global food security contributing to >20% of worldwide aquaculture production. Despite this importance, relatively little is known about bivalve evolutionary adaptation strategies. Difficulties previously associated with highly heterozygous and repetitive regions of bivalve genomes have been overcome by long-read sequencing, enabling the generation of accurate bivalve assemblies. With these resources we have analyzed the genomes of 32 species representing each molluscan class, including 15 bivalve species, to identify gene families that have undergone expansion during bivalve evolution. Gene family expansions across bivalve genomes occur at the point of evolutionary pressures. We uncovered two key factors that shape bivalve evolutionary history: expansion of bivalvia into environmental niches with high stress followed by later exposure to specific pathogenic pressures. The conserved expansion of protein recycling gene families we found across bivalvia is mirrored by adaptations to a sedentary lifestyle seen in plants. These results reflect the ability of bivalves to tolerate high levels of environmental stress and constant exposure to pathogens as filter feeders. The increasing availability of accurate genome assemblies will provide greater resolution to these analyses allowing further points of evolutionary pressure to become clear in other understudied taxa and potentially different populations of a single species.

Genomic Diversity of the Ostreid Herpesvirus Type 1 Across Time and Location and Among Host Species.

The mechanisms underlying virus emergence are rarely well understood, making the appearance of outbreaks largely unpredictable. This is particularly true for pathogens with low per-site mutation rates, such as DNA viruses, that do not exhibit a large amount of evolutionary change among genetic sequences sampled at different time points. However, whole-genome sequencing can reveal the accumulation of novel genetic variation between samples, promising to render most, if not all, microbial pathogens measurably evolving and suitable for analytical techniques derived from population genetic theory. Here, we aim to assess the measurability of evolution on epidemiological time scales of the Ostreid herpesvirus 1 (OsHV-1), a double stranded DNA virus of which a new variant, OsHV-1 µVar, emerged in France in 2008, spreading across Europe and causing dramatic economic and ecological damage. We performed phylogenetic analyses of heterochronous (n = 21) OsHV-1 genomes sampled worldwide. Results show sufficient temporal signal in the viral sequences to proceed with phylogenetic molecular clock analyses and they indicate that the genetic diversity seen in these OsHV-1 isolates has arisen within the past three decades. OsHV-1 samples from France and New Zealand did not cluster together suggesting a spatial structuration of the viral populations. The genome-wide study of simple and complex polymorphisms shows that specific genomic regions are deleted in several isolates or accumulate a high number of substitutions. These contrasting and non-random patterns of polymorphism suggest that some genomic regions are affected by strong selective pressures. Interestingly, we also found variant genotypes within all infected individuals. Altogether, these results provide baseline evidence that whole genome sequencing could be used to study population dynamic processes of OsHV-1, and more broadly herpesviruses.

What do the terms resistance, tolerance, and resilience mean in the case of Ostrea edulis infected by the haplosporidian parasite Bonamia ostreae.

The decline of the European flat oyster Ostrea edulis represents a loss to European coastal economies both in terms of food security and by affecting the Good Environmental Status of the marine environment as set out by the European Council's Marine Strategy Framework Directive (2008/56/EC). Restoration of O. edulis habitat is being widely discussed across Europe, addressing key challenges such as the devastating impact of the haplosporidian parasite Bonamia ostreae. The use of resistant, tolerant, or resilient oysters as restoration broodstock has been proposed by restoration practitioners, but the definitions and implications of these superficially familiar terms have yet to be defined and agreed by all stakeholders. This opinion piece considers the challenges of differentiating Bonamia resistance, tolerance, and resilience; challenges which impede the adoption of robust definitions. We argue that, disease-resistance is reduced susceptibility to infection by the parasite, or active suppression of the parasites ability to multiply and proliferate. Disease-tolerance is the retention of fitness and an ability to neutralise the virulence of the parasite. Disease-resilience is the ability to recover from illness and, at population level, tolerance could be interpreted as resilience. We concede that further work is required to resolve practical uncertainty in applying these definitions, and argue for a collaboration of experts to achieve consensus. Failure to act now might result in the future dispersal of this disease into new locations and populations, because robust definitions are important components of regulatory mechanisms that underpin marine management.

Potential of genomic technologies to improve disease resistance in molluscan aquaculture.

Molluscan aquaculture is a major contributor to global seafood production, but is hampered by infectious disease outbreaks that can cause serious economic losses. Selective breeding has been widely used to improve disease resistance in major agricultural and aquaculture species, and has clear potential in molluscs, albeit its commercial application remains at a formative stage. Advances in genomic technologies, especially the development of cost-efficient genomic selection, have the potential to accelerate genetic improvement. However, tailored approaches are required owing to the distinctive reproductive and life cycle characteristics of molluscan species. Transgenesis and genome editing, in particular CRISPR/Cas systems, have been successfully trialled in molluscs and may further understanding and improvement of genetic resistance to disease through targeted changes to the host genome. Whole-organism genome editing is achievable on a much greater scale compared to other farmed species, making genome-wide CRISPR screening approaches plausible. This review discusses the current state and future potential of selective breeding, genomic tools and genome editing approaches to understand and improve host resistance to infectious disease in molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

A chromosome-level genome assembly for the Pacific oyster Crassostrea gigas.

BACKGROUND: The Pacific oyster (Crassostrea gigas) is a bivalve mollusc with vital roles in coastal ecosystems and aquaculture globally. While extensive genomic tools are available for C. gigas, highly contiguous reference genomes are required to support both fundamental and applied research. Herein we report the creation and annotation of a chromosome-level assembly for C. gigas. FINDINGS: High-coverage long- and short-read sequence data generated on Pacific Biosciences and Illumina platforms were used to generate an initial assembly, which was then scaffolded into 10 pseudo-chromosomes using both Hi-C sequencing and a high-density linkage map. The assembly has a scaffold N50 of 58.4 Mb and a contig N50 of 1.8 Mb, representing a step advance on the previously published C. gigas assembly. Annotation based on Pacific Biosciences Iso-Seq and Illumina RNA-Seq resulted in identification of ∼30,000 putative protein-coding genes. Annotation of putative repeat elements highlighted an enrichment of Helitron rolling-circle transposable elements, suggesting their potential role in shaping the evolution of the C. gigas genome. CONCLUSIONS: This new chromosome-level assembly will be an enabling resource for genetics and genomics studies to support fundamental insight into bivalve biology, as well as for selective breeding of C. gigas in aquaculture.

A Single-Tube HNB-Based Loop-Mediated Isothermal Amplification for the Robust Detection of the Ostreid herpesvirus 1.

The Ostreid herpesvirus 1 species affects shellfish, contributing significantly to high economic losses during production. To counteract the threat related to mortality, there is a need for the development of novel point-of-care testing (POCT) that can be implemented in aquaculture production to prevent disease outbreaks. In this study, a simple, rapid and specific colorimetric loop-mediated isothermal amplification (LAMP) assay has been developed for the detection of Ostreid herpesvirus1 (OsHV-1) and its variants infecting Crassostrea gigas (C. gigas). The LAMP assay has been optimized to use hydroxynaphthol blue (HNB) for visual colorimetric distinction of positive and negative templates. The effect of an additional Tte UvrD helicase enzyme used in the reaction was also evaluated with an improved reaction time of 10 min. Additionally, this study provides a robust workflow for optimization of primers for uncultured viruses using designed target plasmid when DNA availability is limited.

Ephemeral detection of Bonamia exitiosa (Haplosporida) in adult and larval European flat oysters Ostrea edulis in the Solent, United Kingdom.

The haplosporidian parasite Bonamia exitiosa was detected using PCR in four adult and six larval brood samples of the European flat oyster Ostrea edulis from the Solent, UK. This represents the second reported detection of this parasite along the south coast of England. Adult oysters were collected and preserved from seabed populations or restoration broodstock cages between 2015 and 2018. The larvae within brooding adults sampled during 2017 and 2018 were also preserved. Molecular analysis of all samples was performed in 2019. The DNA of B. exitiosa was confirmed to be present within the gill tissue of one oyster within the Portsmouth wild fishery seabed population (n = 48), sampled in November 2015; the congeneric parasite Bonamia ostreae was not detected in this individual. This is the earliest record of B. exitiosa in the Solent. Concurrent presence of both B. ostreae and B. exitiosa, determined by DNA presence, was confirmed in the gill and heart tissue of three mature individuals from broodstock cages sampled in October 2017 (n = 99), two from a location on the River Hamble and one from the Camber Dock in Portsmouth Harbour. B. exitiosa was not detected in the November 2018 broodstock populations. A total of six larval broods were positive for B. exitiosa, with five also positive for B. ostreae. None of the brooding adults were positive for B. exitiosa suggesting that horizontal transmission from the surrounding environment to the brooding larvae is occurring. Further sampling of broodstock populations conducted by the Fish Health Inspectorate at the Centre for Environment, Fisheries and Aquaculture Science in June 2019 did not detect infection of O. edulis by B. exitiosa. These findings together suggest that the pathogen has not currently established in the area.

Developments in marine invertebrate primary culture reveal novel cell morphologies in the model bivalve Crassostrea gigas.

Cell culture provides useful model systems used in a wide range of biological applications, but its utility in marine invertebrates is limited due to the lack of immortalised cell lines. Primary cell and tissue cultures are typically used but remain poorly characterised for oysters, which can cause issues with experimental consistency and reproducibility. Improvements to methods of repeatable isolation, culture, and characterisation of oyster cells and tissues are required to help address these issues. In the current study, systematic improvements have been developed to facilitate the culture of primary cells from adult Pacific oyster tissues and identify novel cell morphologies that have not been reported previously. Cultures analysed by light microscopy, qPCR, and live cell imaging demonstrated maintenance of live, metabolically active Pacific oyster cells for several weeks post-explant. Interestingly, whole hearts dissected from adult oysters were found to continue contracting rhythmically up to 8 weeks after being transferred to a tissue culture system. Mantle tissue explants were also actively moving in the culture system. These improvements in primary cell culture of bivalves may be beneficial for research in ecotoxicology, virology, immunology, and genetic resistance to disease.

Insights into the development of hepatocellular fibrillar inclusions in European flounder (Platichthys flesus) from UK estuaries.

Hepatocellular fibrillar inclusions (HFI) are an unusual pathology of unknown aetiology affecting European flounder (Platichthys flesus), particularly from estuaries historically impacted by pollution. This study demonstrated that the HFI prevalence range was 6-77% at several UK estuaries, with Spearman rank correlation analysis showing a correlation between HFI prevalence and sediment concentrations of ∑PBDEs and ∑HBCDs. The data showed that males exhibit higher HFI prevalence than females, with severity being more pronounced in estuaries exhibiting higher prevalence. HFI were not age associated indicating a subacute condition. Electron microscopy confirmed that HFI were modified proliferating rough endoplasmic reticulum (RER), whilst immunohistochemistry provided evidence of VTG production in HFI of male P. flesus. Despite positive labelling of aberrant VTG production, we could not provide additional evidence of xenoestrogen exposure. Gene transcripts (VTG/CHR) and plasma VTG concentrations (>1 µg ml-1), were only considered elevated in four male fish showing no correlation with HFI severity. Further analysis revealed that reproductively mature female P. flesus i.e. >3-year-old, did not exhibit HFI, whereas males of all ages were affected. This, combined with previous reports that estradiol (E2) can impair mixed function oxygenase activity, supports a hypothesis that harmful chemical metabolites (following phase 1 metabolism of their parent compounds) are potentially responsible for HFIs observed in male and ≤ 3-year-old female fish. Consequently, HFI and xenoestrogenic induced VTG production could be independent of each other resulting from different concurrent toxicopathic mechanisms, although laboratory exposures will likely be the only way to determine the true aetiology of HFI.

Harnessing genomics to fast-track genetic improvement in aquaculture.

Aquaculture is the fastest-growing farmed food sector and will soon become the primary source of fish and shellfish for human diets. In contrast to crop and livestock production, aquaculture production is derived from numerous, exceptionally diverse species that are typically in the early stages of domestication. Genetic improvement of production traits via well-designed, managed breeding programmes has great potential to help meet the rising seafood demand driven by human population growth. Supported by continuous advances in sequencing and bioinformatics, genomics is increasingly being applied across the broad range of aquaculture species and at all stages of the domestication process to optimize selective breeding. In the future, combining genomic selection with biotechnological innovations, such as genome editing and surrogate broodstock technologies, may further expedite genetic improvement in aquaculture.

Assessing relative biomarker responses in marine and freshwater bivalve molluscs following exposure to phosphorus 32 (32P): Application of genotoxicological and molecular biomarkers.

Anthropogenic radionuclides can enter water bodies through accidental or controlled discharges. In order to assess their potential impact, understanding the link between exposure, tissue specific bioaccumulation and radiation dose rate, to biological or biomarker responses in aquatic biota is required. Adopting an integrated, multi-biomarker, multi-species approach, we have investigated potential biological responses induced by short-lived radionuclide, phosphorus-32 (32P, radiophosphorus) in two ecologically important mussel species, the freshwater Dreissena polymorpha (DP) and marine Mytilus galloprovincialis (MG). Adult individuals were exposed to 32P for 10 days, to acquire nominal whole-body average dose rates of 0.10, 1 and 10 mGy d-1, which encompass a screening value of 10 µGy h-1 (0.24 mGy d-1), in accordance with the ERICA tool. Following exposure, a suite of genotoxic biomarkers (DNA damage, γ-H2AX induction and micronucleus [MN] formation) were measured in gill and digestive gland tissues, along with transcriptional expression of selected stress-related genes in both the species (i.e. hsp70/90, sod, cat and gst). Our results demonstrate the relationship between tissue specific dosimetry, where 32P induced a dose-dependent increase, and biological responses independent of species. Gene expression analysis revealed little significant variation across species or tissues. Overall, MG appeared to be more sensitive to short-term damage (i.e. high DNA damage and γ-H2AX induction), particularly in digestive gland. This study contributes to limited knowledge on the transfer and biological impact of radionuclides within differing aquatic systems on a tissue specific level, aiding the development of adequate management and protective strategies.

De novo transcriptome assembly of the Qatari pearl oyster Pinctada imbricata radiata.

The pearl oyster Pinctada imbricata radiata is an iconic species in Qatar, representing an integral part of the nation's cultural heritage and one of the main economic foundations upon which the nation developed. During the early part of the 20th century, nearly half the Qatar population was involved in the pearl oyster industry. However, the fishery has undergone steady decline since the 1930s, and the species is now under threat due to multiple confounding pressures. This manuscript presents the first de novo transcriptome of the Qatari pearl oyster assembled into 30,739 non-redundant coding sequences and with a BUSCO completeness score of 98.4%. Analysis of the transcriptome reveals the close evolutionary distance to the conspecific animal Pinctada imbricata fucata but also highlights differences in immune genes and the presence of distinctive transposon families, suggesting recent adaptive divergence. This data is made available for all to utilise in future studies on the species.

Signatures of selection for bonamiosis resistance in European flat oyster (Ostrea edulis): New genomic tools for breeding programs and management of natural resources.

The European flat oyster (Ostrea edulis) is a highly appreciated mollusk with an important aquaculture production throughout the 20th century, in addition to playing an important role on coastal ecosystems. Overexploitation of natural beds, habitat degradation, introduction of non-native species, and epidemic outbreaks have severely affected this important resource, particularly, the protozoan parasite Bonamia ostreae, which is the main concern affecting its production and conservation. In order to identify genomic regions and markers potentially associated with bonamiosis resistance, six oyster beds distributed throughout the European Atlantic coast were sampled. Three of them have been exposed to this parasite since the early 1980s and showed some degree of innate resistance (long-term affected group, LTA), while the other three were free of B. ostreae at least until sampling date (naïve group, NV). A total of 14,065 SNPs were analyzed, including 37 markers from candidate genes and 14,028 from a medium-density SNP array. Gene diversity was similar between LTA and NV groups suggesting no genetic erosion due to long-term exposure to the parasite, and three population clusters were detected using the whole dataset. Tests for divergent selection between NV and LTA groups detected the presence of a very consistent set of 22 markers, located within a putative single genomic region, which suggests the presence of a major quantitative trait locus associated with B. ostreae resistance. Moreover, 324 outlier loci associated with factors other than bonamiosis were identified allowing fully discrimination of all the oyster beds. A practical tool which included the 84 highest discriminative markers for tracing O. edulis populations was developed and tested with empirical data. Results reported herein could assist the production of stocks with improved resistance to bonamiosis and facilitate the management of oyster beds for recovery production and ecosystem services provided by this species.