Estuaries are the important interface between the land and sea, providing significant environmental, economic, cultural and social values. However, they face unprecedented pressures including eutrophication, harmful algal blooms, habitat loss, and extreme weather due to climate change. Here we present an open access, quality-controlled water quality dataset collected from twelve diverse estuaries spanning 1000 km along the southeastern Australian coastline. Water depth, temperature and salinity data were collected across two years (2018-2021) capturing drought, wildfire and flood periods, using high accuracy Seabird MicroCAT field sensors located within oyster leases. These fully autonomous instruments collected and transmitted data every 10 minutes before downstream quality checking and uploading onto a public website. Simultaneous, high-resolution, longitudinal environmental data collected across multiple estuaries throughout a range of extreme weather events are exceptionally rare in the Southern Hemisphere, yet provide an invaluable resource for the aquaculture industry, researchers and environmental regulators alike.

Climate change is acidifying and warming our oceans, at an unprecedented rate posing a challenge for marine invertebrates vital across the globe for ecological services and food security. Here we show it is possible for resilience to climate change in an ecologically and economically significant oyster without detrimental effects to the energy budget. We exposed 24 pair-mated genetically distinct families of the Sydney rock oyster, Saccostrea glomerata to ocean acidification and warming for 4w and measured their resilience. Resilience was identified as the capacity to defend their acid-base balance without a loss of energy available for Scope for Growth (SFG). Of the 24 families, 13 were better able to defend their acid-base balance while eight had no loss of energy availability with a positive SFG. This study has found oyster families with reslience against climate change without a loss of SFG, is an essential mitigation strategy, in a critical mollusc.

The Pacific Oyster (Crassostrea gigas) is an important aquaculture species, that has been globally affected by “summer mortality” outbreaks. The causes of summer mortality are often unclear, but several contributing factors have been implicated. To better understand the factors contributing to C. gigas summer mortality, we performed a 24 week experiment in tanks filled by a continuous flow of unfiltered surface seawater during a summer mortality high-risk period. Throughout this period, we measured a suite of environmental parameters, with weekly sampling to characterise the bacterial community composition, the abundance and diversity of the Vibrio genus, and the expression of oyster immune response genes. We observed three clear phases of oyster mortality: Phase 1 (day 0 to 75) was characterised by sporadic and low levels (mean < 1% per week) of oyster mortality; Phase 2 (day 82 to 119) involved constant levels of ongoing mortality (mean < 10% per week); and Phase 3 (day 125 to 161) involved a sharp increase in mortality (mean > 10% per week) reaching up to weekly mortality of 51%. The shift between phase 2 and 3 coincided with heavy rainfall and a marine heatwave (MHW) event at our study site, which lasted 13 days, when the average water temperature was elevated by 0.8 °C when compared to the weeks prior and following the MHW, and a maximum water temperature of 27.7 °C was reached. Oyster mortality was positively correlated with the occurrence of the MHW as well as bacterial and Vibrio abundance, and negatively correlated with seawater salinity that decreased from a salinity of 34.6 to 31.4 psu following rainfall events. In addition, the relative abundance of several Vibrio species, including V. campbellii, V. rotiferianus and V. owensii displayed positive correlations to oyster mortality. Shifts in the oyster-associated bacterial community and oyster immunity gene expression profiles occurred in parallel to increased mortality levels, with genes involved in the regulation of antimicrobial peptides (IkB1 and NF-kB p105) significantly over-expressed during high mortality. Overall, we observed a C. gigas mortality event that coincided with a marine heatwave, rainfall event and a shift in bacterial community composition. Our analyses provides further evidence for the putative role of Vibrio in summer mortality events.

Aquaculture mollusc production is predominantly from Asia, with more than 80% of the total biomass produced in China. Vietnam’s annual mollusc production is growing rapidly but is comparatively small given its coastal resources. A significant challenge for future mollusc production, and oysters in particular, is the supply of high-quality spat. Most mollusc spat in Vietnam comes from local wild sources or is imported from China. Particularly in the case of oysters, where wild collection is low and importation of oyster spat into Vietnam from other jurisdictions is not controlled, supply is unreliable and presents a biosecurity risk to the industry. Controlling the life cycle in hatcheries can increase the sustainability and reliability of spat supply with the advantages of genetic enhancement in the long term. Beyond well-established selection methods based on phenotypic data for pedigreed families, the rapid development of genomic technology has enabled innovation in hatchery production based on genetic programs. This technology allows greater insight into oyster genetics in intensive aquaculture production systems. Recent reviews of selection programs demonstrate that inbreeding control via mate selection is an effective strategy for oyster species. Genetic response through combinations of individual and family-based selection can enhance morphometric traits by 10% per generation and disease resistance by 15% per generation in many aquaculture species. Genomic techniques provide information for selecting candidates at an earlier stage and improve prediction accuracy. In this paper, we review the literature on popular genomic tools and breeding techniques used for molluscs, focusing on the Portuguese oyster, Crassostrea angulata. This is to better comprehend how modern quantitative and molecular genetic technologies are being applied in mollusc breeding programs. It considers opportunities for and the feasibility of using genomic-based selection as well as the challenges that are faced in breeding programs transitioning to these new methods.

Ocean acidification will affect marine molluscs, however, transgenerational plasticity (TGP) can ameliorate some effects. Marine molluscs acquire members of their microbiome via the egg, yet we know little about how the microbiome can be influenced by transgenerational exposure to ocean acidification. We exposed adult Sydney Rock oysters (Saccostrea glomerata) from four genotypes to elevated and ambient PCO2 for nine weeks. Larvae were then raised in the same ambient and elevated PCO2 conditions. The relative abundance of bacteria in eggs and larvae were characterised using 16S RNA amplicon sequencing. Parental exposure to elevated PCO2 significantly altered the bacterial community composition of both eggs and larvae, but this was dependent on genotype. Parental exposure to elevated PCO2 caused five core Rhodobacteraceae ASVs to increase in relative abundance, and three Rhodobacteraceae ASVs to decrease in relative abundance. These findings show transfer of maternal microbiomes to larvae is altered by exposure to ocean acidification and this may play a role in TGP.

Teneurin C-terminal associated peptide (TCAP) is an ancient bioactive peptide that is highly conserved in metazoans. TCAP administration reduces cellular and behavioural stress in vertebrate and urochordate models, yet despite numerous studies in higher animals, there is limited knowledge of its role in invertebrates. In particular, there are no studies on TCAP’s effects on the heart of any metazoan, which is a critical organ in the stress response. We used the Sydney rock oyster (SRO) as an invertebrate model to investigate a potential role for sroTCAP in regulating cardiac activity, including during stress. sroTCAP is localized to the neural innervation network of the SRO heart, and suggested binding with various heart proteins related to metabolism and stress, including SOD, GAPDH and metabotropic glutamate receptor. Intramuscular injection of sroTCAP (10 pmol) significantly altered the expression of heart genes that are known to regulate remodelling processes under different conditions, and modulated several gene families responsible for stress mitigation. sroTCAP (1 and 10 pmol) was shown to cause transient bradycardia (heart rate was reduced by up to 63% and for up to 40 min post-administration), indicative of an unstressed state. In summary, this study has established a role for a TCAP in the regulation of cardiac activity through modulation of physiological and molecular components associated with energy conservation, stress and adaptation. This represents a novel function for TCAP and may have implications for higher-order metazoans.

Teneurin C-terminal associated peptide (TCAP) is an ancient bioactive peptide that is highly conserved in metazoans. TCAP administration reduces cellular and behavioral stress in vertebrate and urochordate models. There is little information for invertebrates regarding the existence or function of a TCAP. This study used the Sydney rock oyster (SRO) as a molluscan model to characterize an invertebrate TCAP, from molecular gene analysis to its physiological effects associated with hemocyte phagocytosis. We report a single teneurin gene (and 4 teneurin splice variants), which encodes a precursor with TCAP that shares a vertebrate-like motif, and is similar to that of other molluscan classes (gastropod, cephalopod), arthropods and echinoderms. TCAP was identified in all SRO tissues using western blotting at 1-2 different molecular weights (~22 kDa and ~37kDa), supporting precursor cleavage variation. In SRO hemolymph, TCAP was spatially localized to the cytosol of hemocytes, and with particularly high density immunoreactivity in granules. Based on ‘pull-down’ assays, the SRO TCAP binds to GAPDH, suggesting that TCAP may protect cells from apoptosis under oxidative stress. Compared to sham injection, the intramuscular administration of TCAP (5 pmol) into oysters modulated their immune system by significantly reducing hemocyte phagocytosis under stress conditions (low salinity and high temperature). TCAP administration also significantly reduced hemocyte reactive oxygen species production at ambient conditions and after 48 h stress, compared to sham injection. Transcriptomic hemocyte analysis of stressed oysters administered with TCAP demonstrated significant changes in expression of genes associated with key metabolic, protective and immune functions. In summary, this study established a role for TCAP in oysters through modulation of physiological and molecular functions associated with energy conservation, stress and cellular defense.

The Sydney rock oyster (SRO; Saccostrea glomerata) is the most intensively farmed oyster species in Australia however, Queensland unknown (QX) disease has resulted in substantial losses and impeded productivity. QX disease is caused by infection with the parasite Marteilia sydneyi, and like other diseases, outbreaks are driven by a series of complex environmental and host factors such as seasonality, seawater salinity and oyster genetics. A potential but understudied factor in QX disease is the SRO microbiota, which we sought to examine before and during a QX disease outbreak. Using 16S rRNA (V1 - V3 region) amplicon sequencing, we examined the microbiota of SROs deployed in an estuary where QX disease occurs, with sampling conducted fortnightly over 22 weeks. Marteilia sydneyi was detected in the SROs by PCR (QX-positive), 16 weeks after the first sampling event and sporonts were observed in the digestive gland two weeks later on. There were no apparent patterns observed between the microbiota of QX-positive SROs with and without digestive gland sporonts however, the microbiota of QX-positive SROs was significantly different from those sampled prior to detection of M. sydneyi and from those negative for M. sydneyi post detection. As a result, shifts in microbiota structure occurred before sporulation in the digestive gland and either before or shortly after pathogen colonisation. The microbiota shifts associated with QX-positive oysters were principally driven by a relative abundance increase of operational taxonomic units (OTUs) assigned to unclassified species of the Borrelia and Candidatus Hepatoplasma genera and a relative abundance decrease in an OTU assigned to an unclassified species of the Mycoplasma genus. Since Mycoplasma species are common microbiota features of SROs and other oysters, we propose that there may be an important ecological link between Mycoplasma species and the health state of SROs.

In a selective breeding program for the Portuguese oyster, Crassostrea angulata, oysters are cultured in several environments including the broodstock for breeding (Cat Ba) and production (Van Don) locations. This raises the potential that a genotype by environment interaction may occur and oyster genotypes may re-rank across environments. Whole weight and uniformity in weight at the time of harvest are two of the most important economic traits for oyster production in northern Vietnam. Limited information, particularly with respect to uniformity, exists for C. angulata. Our primary aim was to quantify the genetic variance in harvest whole weight and its uniformity and the genetic correlation between them using a double hierarchical generalized linear model. We investigated the degree of genotype by environment interaction on the uniformity of harvest whole weight in the two environments. Since the weights had a skewed distribution, we also investigated the affect of log-transformation on the estimates of the genetic variance in harvest whole weight, its uniformity, the genetic correlation between them and on the genotype by environment interactions. Results showed that heritability estimates for uniformity of harvest whole weight were low for Cat Ba (0.07) and Van Don (0.06) environments using standardized data, while these estimates were lower using log transformation for both environments. There was no re-ranking of oysters with regard to harvest whole weight and uniformity across environments; the genetic correlations were effectively one, with and without log transformation. Therefore, there is no scope to improve uniformity in the Portuguese oyster population while increasing whole weight at harvest.

Genetic improvement for quality traits, especially color and meat yield, has been limited in aquaculture because the assessment of these traits requires that the animals be slaughtered first. Genotyping technologies do, however, provide an opportunity to improve the selection efficiency for these traits. The main purpose of this study is to assess the potential for using genomic information to improve meat yield (soft tissue weight and condition index), body shape (cup and fan ratios), color (shell and mantle), and whole weight traits at harvest in the Portuguese oyster, Crassostrea angulata. The study consisted of 647 oysters: 188 oysters from 57 full-sib families from the first generation and 459 oysters from 33 full-sib families from the second generation. The number per family ranged from two to eight oysters for the first and 12–15 oysters for the second generation. After quality control, a set of 13,048 markers were analyzed to estimate the genetic parameters (heritability and genetic correlation) and predictive accuracy of the genomic selection for these traits. The multi-locus mixed model analysis indicated high estimates of heritability for meat yield traits: 0.43 for soft tissue weight and 0.77 for condition index. The estimated genomic heritabilities were 0.45 for whole weight, 0.24 for cup ratio, and 0.33 for fan ratio and ranged from 0.14 to 0.54 for color traits. The genetic correlations among whole weight, meat yield, and body shape traits were favorably positive, suggesting that the selection for whole weight would have beneficial effects on meat yield and body shape traits. Of paramount importance is the fact that the genomic prediction showed moderate to high accuracy for the traits studied (0.38–0.92). Therefore, there are good prospects to improve whole weight, meat yield, body shape, and color traits using genomic information. A multi-trait selection program using the genomic information can boost the genetic gain and minimize inbreeding in the long-term for this population.

In recent decades, developing countries with traditionally strong links with aquaculture have adopted genetic improvement programs for newly domesticated and captively bred aquaculture species. Accurate genealogical and pedigree information is an essential tool for efficient and successful genetic improvement programs and traditionally this type of information was obtained by separate rearing of families and physical tagging of the animals. More recently, molecular based pedigrees (e.g. using SNPs) is becoming a potential option, which also may permit the communal rearing of families at a stage earlier than usual, prior to physical tagging. Here we consider the utility of molecular based pedigrees in a multi-generation breeding and selection program for the Portuguese oyster (Crassostrea angulata) and assess parameters for parentage assignment.Using 18,849 SNP markers obtained from 489 sequenced oysters including parents and their progeny, our results indicated that approximately 400 SNPs with minor allele frequency ≥ 0.3 were capable of assigning offspring to their 33 parents (18 sires and 15 dams).

Microbiomes can both influence and be influenced by metabolism, but this relationship remains unexplored for invertebrates. We examined the relationship between microbiome and metabolism in response to climate change using oysters as a model marine invertebrate. Oysters form economies and ecosystems across the globe, yet are vulnerable to climate change. Nine genetic lineages of the oyster Saccostrea glomerata were exposed to ambient and elevated temperature and PCO2 treatments. The metabolic rate (MR) and metabolic by-products of extracellular pH and CO2 were measured. The oyster-associated bacterial community in haemolymph was characterised using 16 s rRNA gene sequencing. We found a significant negative relationship between MR and bacterial richness. Bacterial community composition was also significantly influenced by MR, extracellular CO2 and extracellular pH. The effects of extracellular CO2 depended on genotype, and the effects of extracellular pH depended on CO2 and temperature treatments. Changes in MR aligned with a shift in the relative abundance of 152 Amplicon Sequencing Variants (ASVs), with 113 negatively correlated with MR. Some spirochaete ASVs showed positive relationships with MR. We have identified a clear relationship between host metabolism and the microbiome in oysters. Altering this relationship will likely have consequences for the 12 billion USD oyster economy.

The aquaculture industry uses microalgae as a live feed for juvenile oysters in hatcheries to meet their nutritional requirements, including their need for several essential Poly Unsaturated Fatty Acids (PUFAs) such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA). The mass culture of microalgae is not only a major bottleneck for the production of juvenile oysters, but also a significant cost, accounting for 20-50% of hatchery operating costs. Currently, low biomass concentrations, high production costs and poor cultivation systems limit the quantity and quality of microalgae feed. This study focused on Tisochrysis lutea, a microalgae species commonly used in aquaculture, and we assessed the potential of photobioreactors with an improved light source and CO2 input to increase biomass production and improve biochemical composition of algal feed. Two photobioreactor systems were compared: the current industry set up (DPI) comprising fluorescent lighting and minimal CO2 input versus an optimized system utilising LEDs and increased CO2. Cultures of T. lutea were monitored over a 12-day growth period and harvested on day 14 for biochemical analysis. Final cell density was significantly higher in the optimized system relative to the conventional culture systems (6.2 x 10(6) cells / mL versus 3.7 x 10(6) cells / mL, respectively). The biochemical profile of T. lutea was not significantly different between the two photobioreactors systems. The algal biomass produced during this comparative experiment was used in a feeding trial on oyster spat, Saccostrea glomerata. Spat fed with algae produced in optimized vs conventional photobioreactors showed no significant difference in growth, but oyster spat fed with T. lutea grown in optimized photobioreactors did show a significant increase in their EPA content. Overall, our results contribute to our understanding of how altered culture conditions affect microalgal production and biochemical composition, and subsequently of oyster spat. This study further supports the potential of LEDs to reduce operating costs for oyster hatcheries but also to significantly improve microalgae yield with no negative effect on the condition of oyster spats production.

Genomic selection has been widely used in terrestrial animals but has had limited application in aquaculture due to relatively high genotyping costs. Genomic information has an important role in improving the prediction accuracy of breeding values, especially for traits that are difficult or expensive to measure. The purposes of this study were to (i) further evaluate the use of genomic information to improve prediction accuracies of breeding values from, (ii) compare different prediction methods (BayesA, BayesCπ and GBLUP) on prediction accuracies in our field data, and (iii) investigate the effects of different SNP marker densities on prediction accuracies of traits in the Portuguese oyster (Crassostrea angulata). The traits studied are all of economic importance and included morphometric traits (shell length, shell width, shell depth, shell weight), edibility traits (tenderness, taste, moisture content), and disease traits (Polydora sp. and Marteilioides chungmuensis). A total of 18,849 single nucleotide polymorphisms were obtained from genotyping by sequencing and used to estimate genetic parameters (heritability and genetic correlation) and the prediction accuracy of genomic selection for these traits. Multi-locus mixed model analysis indicated high estimates of heritability for edibility traits; 0.44 for moisture content, 0.59 for taste, and 0.72 for tenderness. The morphometric traits, shell length, shell width, shell depth and shell weight had estimated genomic heritabilities ranging from 0.28 to 0.55. The genomic heritabilities were relatively low for the disease related traits: Polydora sp. prevalence (0.11) and M. chungmuensis (0.10). Genomic correlations between whole weight and other morphometric traits were from moderate to high and positive (0.58–0.90). However, unfavourably positive genomic correlations were observed between whole weight and the disease traits (0.35–0.37). The genomic best linear unbiased prediction method (GBLUP) showed slightly higher accuracy for the traits studied (0.240–0.794) compared with both BayesA and BayesCπ methods but these differences were not significant. In addition, there is a large potential for using low-density SNP markers for genomic selection in this population at a number of 3000 SNPs. Therefore, there is the prospect to improve morphometric, edibility and disease related traits using genomic information in this species.

Oyster microbiomes are integral to healthy function and can be altered by climate change conditions. Genetic variation among oysters is known to influence the response of oysters to climate change and may ameliorate any adverse effects on oyster microbiome; however, this remains unstudied. Nine full-sibling selected breeding lines of the Sydney rock oyster (Saccostrea glomerata) were exposed to predicted warming (ambient = 24°C, elevated = 28°C) and ocean acidification (ambient pCO2 = 400, elevated pCO2 = 1000 µatm) for 4 weeks. The haemolymph bacterial microbiome was characterized using 16S rRNA (V3-V4) gene sequencing and varied among oyster lines in the control (ambient pCO2, 24°C) treatment. Microbiomes were also altered by climate change dependent on oyster lines. Bacterial α-diversity increased in response to elevated pCO2 in two selected lines, while bacterial β-diversity was significantly altered by combinations of elevated pCO2 and temperature in four selected lines. Climate change treatments caused shifts in the abundance of multiple amplicon sequence variants driving change in the microbiome of some selected lines. We show that oyster genetic background may influence the Sydney rock oyster haemolymph microbiome under climate change and that future assisted evolution breeding programs to enhance resilience should consider the oyster microbiome.

The wellbeing of marine organisms is connected to their microbiome. Oysters are a vital food source and provide ecological services, yet little is known about how climate change such as ocean acidification and warming will affect their microbiome. We exposed the Sydney rock oyster, Saccostrea glomerata, to orthogonal combinations of temperature (24, 28 °C) and pCO2 (400 and 1000 μatm) for eight weeks and used amplicon sequencing of the 16S rRNA (V3-V4) gene to characterise the bacterial community in haemolymph. Overall, elevated pCO2 and temperature interacted to alter the microbiome of oysters, with a clear partitioning of treatments in CAP ordinations. Elevated pCO2 was the strongest driver of species diversity and richness and elevated temperature also increased species richness. Climate change, both ocean acidification and warming, will alter the microbiome of S. glomerata which may increase the susceptibility of oysters to disease.

The Portuguese oyster Crassostrea angulata is one of the economically important mollusc species but to date there has been no genetic improvement program conducted for this species. The principal aim of this study was to assess genetic gains achieved from a selection program for improved whole weight at harvest over three generations from 2015 to 2017. A total of 30,212 oysters had performance data recorded; they were the offspring of 323 sires and 461 dams. Linear mixed models were used to estimate genetic parameters of the traits studied. Selection for improved harvest whole weight achieved significant direct response by 17.4% after three generations. There were also positively correlated changes in soft tissue weight and shell shape traits (3.67% and 3.98%, respectively). The genetic gains in harvest whole weight and favourably correlated changes in meat yield and shell shape traits were consistent with the positive genetic correlations of whole weight with these traits (0.08-0.50). There is abundant genetic variation in the selection criterion (heritability for whole weight = 0.24), suggesting that the present population of Portuguese oysters will continue showing responses to future selection. The estimates of heritability for meat yield and shell shape traits were low (0.10-0.16). Hence, a multi-trait selection program can be used to enhance overall productivity and economic return for oyster producers.

Queensland unknown (QX) disease is a significant cause of economic loss for the Sydney rock oyster (SRO) aquaculture industry. Evidence is emerging that QX disease is multi-factorial in nature, with a number of environmental and host factors contributing to disease dynamics. Efforts to mitigate the impacts of QX disease are primarily focused on breeding for disease resistance however, the mechanisms that drive disease resistance are poorly understood. One potential factor influencing disease resistance is the microbiota. To determine the influence of location, season and disease resistance on the SRO microbiota, we used 16S rRNA (V1 – V3 region) amplicon sequencing. The microbiota of six SRO families with two categorised as QX-resistant and four as QX-susceptible, deployed to two different locations (Port Stephens and Wallis Lake, NSW, Australia) and over two seasons (Austral summer and winter), were characterised. As expected, the SRO microbiota was distinct to the microbial community found in seawater. Further, the SRO microbiota was significantly influenced by location and season, with operational taxonomic units (OTUs) assigned to the Candidatus Hepatoplasma and Endozoicomonas genera identified as significant drivers of microbiota dissimilarity between locations and seasons. Disease resistance also significantly influenced the SRO microbiota but only at the winter time point which is before the typical QX disease period. Overall, OTUs assigned to the Mycoplasma, Borrelia and Endozoicomonas genera were over-represented in QX-resistant SRO microbiota, whereas members of the Pseudoalteromonas, Vibrio, and Candidatus Hepatoplasma genera were over-represented in QX-sensitive microbiota. These findings confirm the influencing role of location and season on the microbiota structure as evidenced in other molluscan species, but also provide preliminary evidence that the microbiota assemblage before the QX disease period may be important for resistance to disease and may provide new avenues for managing SRO aquaculture in the future.

Different organs of a host represent distinct microenvironments resulting in the establishment of multiple discrete bacterial communities within a host. These discrete bacterial communities can also vary according to geographical location. For the Pacific oyster, Crassostrea gigas, the factors governing bacterial diversity and abundance of different oyster microenvironments are poorly understood. In this study, the factors shaping bacterial abundance, diversity and composition associated with the C. gigas mantle, gill, adductor muscle and digestive gland were characterised using 16S (V3-V4) rRNA amplicon sequencing across six discrete estuaries. Both location and tissue-type, with tissue-type being the stronger determinant, were factors driving bacterial community composition. Bacterial communities from wave-dominated estuaries had similar compositions and higher bacterial abundance despite being geographically distant from one another, possibly indicating that functional estuarine morphology characteristics are a factor shaping the oyster bacterial community. Despite the bacterial community heterogeneity, examinations of the core bacterial community identified Spirochaetaceae bacteria as conserved across all sites and samples. Whereas members of the Vulcaniibacterium, Spirochaetaceae and Margulisbacteria, and Polynucleobacter were regionally conserved members of the digestive gland, gill and mantle bacterial communities, respectively. This indicates that baseline bacterial community profiles for specific locations are necessary when investigating bacterial communities in oyster health.

Production of the Portuguese oyster, Crassostrea angulata, is increasing rapidly in Vietnam, and consumer preferences for oyster appearance and meat quality have become more economically important. However, no studies have been conducted to estimate genetic parameters for meat colour, edible traits and the occurrence of disfiguring parasite infestation in any oyster species. In this study, we report genetic parameters for parasite disease traits, objective measurements of shell and mantle colours (L*, a* and b* colour) and edible traits (water holding capacity, taste and tenderness) for a selectively bred Portuguese oyster breeding population. We applied a restricted maximum likelihood method to estimate genetic parameters for the traits studied. Heritability estimates for the colour measurements of shell and mantle were low to high (h(2) = 0.13-0.57), whereas those obtained for the parasite diseases on shells and tissues were close to zero (h(2) = 0.011). Interestingly, there was genetic variation in edible traits (h(2) = 0.04-0.17). The genetic correlations among whole weight and parasitic shell and tissue disease traits was significant but unfavourable, whereas the genetic correlations between whole weight and water holding capacity was considerable and favourable (r(g) = 0.87 +/- 0.05). These suggest that a multi-trait restricted selection index should be used to improve productivity and economic returns for the oyster aquaculture sector.