Skip to main content

The development of multiplex PCR assays for the rapid identification of multiple Saccostrea species, and their practical applications in restoration and aquaculture

Abstract

Background

The ecology and biology of oysters (Ostreidae) across the tropics is poorly understood. Morphological plasticity and shared characteristics among oysters have resulted in the misidentification of species, creating challenges for understanding basic species-specific biological information that is required for restoration and aquaculture. Genetic barcoding has proven essential for accurate species identification and understanding species geographic ranges. To reduce the costs of molecular species identification we developed multiplex assays using the cytochrome c oxidase subunit I (COI or cox1) barcoding gene for the rapid identification of five species of oysters within the genus Saccostrea that are commonly found in Queensland, Australia: Saccostrea glomerata, Saccostrea lineage B, Saccostrea lineage F, Saccostrea lineage G, and Saccostrea spathulata (lineage J).

Results

Multiplex assays were successful in species-specific amplification of targeted species. The practical application of these primers was tested on wild spat collected from a pilot restoration project in Moreton Bay, Queensland, with identified species (S. glomerata, lineage B and lineage G) validated by Sanger sequencing. DNA sampling by extraction of oyster pallial fluid was also tested on adult oysters collected from the Noosa estuary in Queensland to assess whether oysters were able to be identified non-destructively. DNA concentrations as low as 1 ng/ μL still amplified in most cases, allowing for identification, and mortality at 6 weeks post pallial fluid collection was low (3 out of 104 sampled oysters).

Conclusion

These multiplex assays will be essential tools for species identification in future studies, and we successfully demonstrate their practical application in both restoration and aquaculture contexts in Queensland. The multiplex assays developed in this study outline easily replicable methods for the development of additional species-specific primer sets for the rapid identification of other species of Saccostrea found across the Indo-Pacific, which will be instrumental in unravelling the taxonomic ambiguities within this genus in tropical regions.

Peer Review reports

Background

The global decline of oyster reefs is unprecedented compared to other intertidal estuarine marine habitats, and they are now considered functionally extinct in many regions [1]. There is a widely recognised need to restore oyster reefs and the ecosystem services they once provided, including water filtration, nutrient cycling, shoreline protection, and habitat for fish and marine invertebrates [2,3,4,5].

Globally, tropical regions contain over four times more reef-building species of oysters than temperate regions [6]. Molecular techniques are required for accurate identification, however, as shell morphology is highly variable which has resulted in species being visually misidentified [7, 8]. Misidentification is especially prevalent within the genus Saccostrea, where multiple species occur across Australia and the Indo-Pacific. Multiple species within this genus have been collectively lumped into two “superspecies”; Saccostrea scyphophilla (previously Saccostrea mordax) which is comprised of at least three genetically distinct lineages, designated A – C, and Saccostrea cuccullata which contains at least ten genetically distinct lineages (designated A – J), with most lineages likely representing a separate species [7, 9,10,11,12]. For example, Saccostrea scyphophilla lineage C has now been described as a new species Saccostrea mordoides [13], Saccostrea cuccullata lineage J has recently been assigned to the species Saccostrea spathulata [12], Saccostrea cuccullata lineage F may correspond to the species Saccostrea malabonensis [10], and accidental sampling of Saccostrea cuccullata lineage B and G during a Saccostrea glomerata population study showed that these three lineages are likely separate species through nuclear DNA markers [14].

Tropical Australia has been identified as having high potential for the establishment and expansion of oyster aquaculture and oyster reef restoration, however the high diversity of species and taxonomic uncertainties pose a challenge [6, 12, 15]. For example, Saccostrea echinata is a small spiky oyster endemic to the Indo-Pacific and likely Australia [7, 11, 16], however, the name continues to be incorrectly applied to the tropical blacklip oyster (now known to be S. spathulata [12]), a promising species for aquaculture in the Australian tropics [17,18,19,20]. The use of “superspecies” names in scientific literature, especially S. cuccullata, is also problematic as it is unclear which specific lineage research is reporting on [21,22,23,24]. The past use of what we now know to be incorrect or “superspecies” names has trickled through to government and industry, causing difficulty for summarising species-specific biological information necessary for advancing restoration and aquaculture [6]. There is a need for all stakeholders to adopt the lineage naming system proposed by Lam and Morton, 2006 to delineate species and species-specific information until species boundaries have been ascertained and scientific names can be assigned [9, 10, 12, 15].

Molecular research has been essential in the accurate identification of morphologically similar oyster species and resolving taxonomic uncertainties. The cost of molecular work can present a barrier to research, however, the development of tools such as species-specific primers can drastically reduce costs by allowing identification by gel electrophoresis rather than nucleotide sequencing. For example, in China species-specific primer sets were developed using the mitochondrial cytochrome oxidase I (COI) gene for use within a multiplex PCR for the identification of five co-occuring oysters within the genus Magallana (previously Crassostrea): Magallana angulata, Magallana ariakensis, Magallana gigas, Magallana hongkongensis and Magallana sikamea [25]. The development of species-specific primers can also allow for the analysis of larger sample sizes for ecological studies. Primer sets using the mitochondrial ND5 gene were developed for the identification of M. sikamea and M. ariakensis which are known to create mixed species oyster reefs in China [26]. The development of these primer sets facilitated research that explored settlement preferences between the two species and highlighted potential reasons for limited restoration success [27].

While most lineages within the S. cucullata complex are yet to be assigned to a new or existing species, molecular work has enabled a better understanding of habitat preferences and geographic ranges which will aid in species identification and understanding aquaculture and restoration potential in localised regions [7, 9, 15, 28]. Five species belonging to the S. cuccullata complex have recently been documented in the Queensland tropics: lineage B, lineage F, lineage G, lineage I and Saccostrea spathulata (lineage J) [12]. Most of these species are sympatric with one another as well as with S. glomerata (in subtropical locations), and in most cases are unable to be distinguished from one another morphologically. To enable the rapid identification of these species we developed species-specific primers for use in multiplex PCRs and gel electrophoresis and demonstrate their practical applications in both restoration and aquaculture.

Methods

Primer design

The species chosen for primer development were S. glomerata, lineage B, lineage F, lineage G and S. spathulata. Lineage I was not selected as this species has several distinctive characteristics including broad but delicate purple lobe-like squamae (lamellae) on the lower/bottom (left) valve, internal shell iridescence, and preferred settlement on the roots of mangroves, enabling straightforward identification [12, 16]. S. scyphophilla and S. mordoides were also excluded as these oysters are not part of the S. cuccullata complex and can be distinguished morphologically by their distinctive shape, deep longitudinal ridges, and purple colouration [12].

The cytochrome c oxidase subunit I (COI or cox1) barcoding gene was chosen as the target for primer design based on the availability of sequences for all species of interest, and the presence of sufficient nucleotide diversity to enable the design of species-specific primers. All publicly submitted COI gene sequences from Australia and the Indo-Pacific were compiled for each target species from NCBI (https://www.ncbi.nlm.nih.gov), and one sequence representing the consensus was used for primer design. Several primer sets were designed for each species using Primer3 (https://primer3.ut.ee) specifying the following parameters: optimal primer size of 20 base pairs (minimum 18 and maximum 23); optimal annealing temperature of 59 °C (minimum 57 °C and maximum 62 °C); and 50 % optimal GC content (minimum 30 % and maximum 70 %). Final primers were then selected based on the uniqueness of nucleotide positions compared to the same gene regions of other species, to ensure that polymorphic sites within species were avoided, and to ensure that there was a minimum difference of 55 base pairs between the product size of each species for readability in gel electrophoresis. Details of species-specific primer sets, product sizes and sequence of origin for primer design can be found in Table 1. The number of unique nucleotide positions between each species/lineage at each region of all species-specific forward and reverse primers can be found in Table 2. Nucleotide sequences among all species and lineages at primer-annealing regions can be found in Fig. 1.

Table 1 Summary of primers designed for morphologically similar Saccostrea species in Queensland
Table 2 The number of unique nucleotide positions between each species/lineage at each region of species-specific forward and reverse primers
Fig. 1
figure 1

Nucleotide sequences among all species and lineages at primer-annealing regions. Forward primers are highlighted in green and reverse primers are highlighted in purple

Primers were tested on oyster DNA from specimens collected across Queensland that had previously been unambiguously identified by 16S and/or COI barcoding [12]. Species-specific primer sets were first tested individually on confirmed DNA from 20 individuals of the target species to confirm functionality, except for lineage F where only 10 samples were used due to low sample availability. In each reaction the following quantities were added to the PCR mixture; 2 μL of New England Biolabs ThermoPol® 10X buffer (Catalog # B9004); 1 μL of each primer (forward and reverse) at a concentration of 20 μM; 0.4 μL of New England Biolabs® Deoxynucleotide (dNTP) Solution Mix (Catalog #N0447) at a concentration of 10 mM; and 0.2 μL of New England Biolabs® Taq DNA polymerase (Catalog # M0267). All DNA used for testing was added at a concentration of 25 – 50 ng/ μL, and PCR mixtures were then brought to 20 μL with molecular grade water. A Bio-Rad T-100™ thermal cycler was used for PCR reactions and the thermoprofile consisted of 5 min at 94 °C as initial denaturation, followed by 25 cycles of 30 sec at 94 °C for denaturation, 30 sec at 62 °C for annealing, 30 sec at 68 °C for extension, and final extension at 68°C for 5 min. 10 μL of each PCR reaction was then visualised by gel electrophoresis at 110 V for 30 minutes in a 2.5 % TAE agarose gel with a 100 bp DNA ladder.

We then tested each primer set on confirmed DNA from the seven remaining Saccostrea lineages and species found in Queensland (including lineage I, S. scyphophilla and S. mordoides) to check for amplification of non-targeted species using the same thermoprofile. Where required, gradient PCRs (with annealing temperatures ranging from 62 – 70 °C) were conducted for further optimisation. Once primers were optimized to amplify targeted species, we tested all five primer pairs together in a single multiplex PCR. Unfortunately, successful multiplex PCR of all five primers was unable to be achieved as this decreased the positive detection rate of all species. Instead, four-species multiplex assays using all combinations of primer sets using the same initial thermoprofile were successful.

Practical applications – restoration

Gabion baskets filled with sterilised oyster shell are becoming an increasingly popular method of oyster reef restoration in the US [29, 30]. This method has now been adopted in Queensland as a part of OzFish Unlimited’s “Shellfish Revolution” restoration initiative (see https://ozfish.org.au/programs/shellfish-restoration/). To test our primers on wild oysters recruited in a restoration pilot project, we deployed 15 gabion baskets built by OzFish (Fig. 2A) at a commercial oyster lease in Moreton Bay, Queensland (-27° 27’ 31.9644” N and 153° 23’ 45.4236” S) in November 2022 to coincide with known S. glomerata spawning season [31]. To account for other species of Saccostrea that may have had different spawning triggers/seasons [32,33,34], baskets were left to establish for one year. In November 2023, five spat per basket (n = 75) were collected for identification, with care taken to sample a broad range of morphologies and sizes (Fig. 2B).

Fig. 2
figure 2

A Gabion baskets designed and used by OzFish for oyster reef restoration in southeast Queensland; and B an example of different spat morphologies recruiting to oyster shell in gabion baskets. Photo credit: Robert Porter

To extract DNA to enable species identifications to be made, a piece of mantle tissue was sampled from spat over 2 cm in size and all tissue was collected for smaller spat. Tissue samples were placed in sterilised microcentrifuge tubes upon collection and stored in 70% ethanol for later analysis. DNA was then extracted using the DNeasy Blood and Tissue Kit (Qiagen), and yields were measured in a Qubit Fluorometer and diluted to concentrations of 20 – 50 ng/ μL.

Although S. glomerata is the recognised reef-building species in southeast Queensland, lineage B and lineage G are also known to occur in this region [12, 14]. Lineage B, lineage G and S. glomerata primer sets were therefore selected for use in a multiplex PCR for species identification. DNA was amplified in a multiplex PCR using the same PCR mixture and thermoprofile described above, and 10 μL of PCR product was then visualised via agarose gel electrophoresis. Although lineage F and lineage J do not have known distributions in Moreton Bay [12], all unamplified DNA was run in a multiplex PCR using these primer sets as a precaution, however no DNA amplified. A subset of samples from each identified species as well as all unidentified samples were then amplified in a PCR using 16S ribosomal RNA (16S) primer sets described in McDougall and Walker, 2021 [35] and sent to Macrogen for Sanger sequencing to validate results.

Practical applications – aquaculture

Tissue sampling from oysters for the purpose of DNA extractions often results in mortality. This is problematic for many applications, including aquaculture which relies on the identification of live specimens for future hatchery production. Lineage G is a strong candidate for aquaculture given its morphological similarity to S. glomerata and widespread distribution across Australia and the Indo-Pacific [12]. To positively identify live lineage G specimens, we collected 20 oysters from the Noosa estuary in southeast Queensland (-26° 23’ 37.3776”N and 153° 2’ 30.6888” E) where populations of lineage G and S. glomerata are found [12]. To non-destructively sample DNA, we first trialled magnesium chloride (MgCl2) baths following the methods of Nowland et al., 2021. We had limited success obtaining DNA samples with this method which we attribute to the small size of our collected oysters compared to the tropical blacklip.

To trial a second method, 104 oysters were collected from the same estuary and transported to the Bribie Island Research Centre. A mechanical handheld 1 mm drill was used to make a small hole in the ventral half of the right (upper) valve, with care to avoid contact with the animal inside. A 23G 0.65x32mm needle on a 1ml syringe was inserted and 250 – 500 μL of pallial fluid was extracted, placed in a microcentrifuge tube and stored at 4 °C until processing. The hole was then plugged with beeswax and covered with quick dry waterproof superglue, and the oyster retained in aquaria at the facility to monitor recovery.

Prior to DNA extractions, microcentrifuge tubes containing pallial fluid were spun down for 5 minutes at 13,000 rpm to create a pellet. Due to expected low DNA yields from pallial fluid, the DNeasy Blood and Tissue Kit (Qiagen) protocol was modified as DNA binding in silica-solution versus silica spin columns can increase DNA yields [36]. In summary, 180 μL of Buffer ATL (Qiagen) and 10 μL of ProK buffer was added to the pallial fluid and incubated at 56 °C for 3 – 6 hours. 400 μL of Buffer AL was added to the solution and vortexed for 15 seconds. 400 μL of 100% ethanol was added to the solution and vortexed for 15 seconds, followed by incubation at room temperature for 5 minutes. 20 μL of silica solution was added, mixed, and incubated at 50 °C for 15 minutes. The solution was then centrifuged for 30 seconds at 13,000 rpm, and supernatant was discarded. 500 μL of Buffer AW1 (Qiagen) was added and the silica pellet was resuspended. The solution was centrifuged at 13,000 rpm for 30 seconds, supernatant was discarded, and the previous wash step was repeated with Buffer AW2 (Qiagen). Once supernatant from the AW2 wash step was discarded, the silica pellet was left to air dry for 5 minutes. The silica pellet was then resuspended in 10 μL of H2O and centrifuged at 500 rpm for 2 minutes. DNA yields were measured in a Qubit Fluorometer and ranged from approximately 15 ng/ μL to 50 ng/ μL, however some samples were less than 1 ng/ μL. Because S. glomerata and lineage G were the only species of interest for this experiment, these were the only primers used in the multiplex PCR.

Results

All species-specific primer pairs successfully and exclusively amplified targeted species with the exception of the lineage B primer set, which amplified a product of the same expected size from all lineage I DNA. Six alternative lineage B primer sets were designed and tested but none showed species specificity. A gradient PCR was then performed from 62 °C – 68 °C using the original primer set to determine the optimal annealing temperature for species specificity; this was achieved with an annealing temperature of 68°C. Occasional non-specific amplification of lineage G samples was also observed with the lineage B primers at the expected size for lineage B products (four out of 55 tested lineage G DNA samples). We repeated the gradient PCR, however, two bands were still faintly present at 68 °C (Fig. 3a). This was resolved by raising the annealing temperature to 70 °C, however this decreased the positive detection rate of the Lineage B primers (Fig. 3b). Species specificity of the Lineage B primer sets is shown in Fig. 3c.

Fig. 3
figure 3

Gel electrophoresis of PCR products using the lineage B primer set (380 bp): a lineage B and lineage G DNA amplified at an annealing temperature of 68 °C (* denotes non-targeted amplification of lineage G DNA); b lineage B and lineage G DNA amplified at an annealing temperature of 70 °C; and c PCR products of all Queensland Saccostrea species DNA amplified at an annealing temperature of 70 °C (lineage G DNA amplified in this PCR had previously produced non-targeted bands). NC = negative control. See Supplementary Figures S1 – S3 for original photographs

Four-species multiplex assays using all combinations of primer sets were successful in species-specific amplification of targeted species (Fig. 4). Although a single five-species multiplex assay was unable to be achieved, S. glomerata has a subtropical distribution and does not overlap with lineage F and has very little overlap with S. spathulata. It is therefore highly unlikely that a five-species multiplex assay would be required for species identification in Queensland. The S. glomerata primer showed some non-specific amplification, producing both the expected product at 185 bp, as well as a larger band of approximately 600 bp (Fig. 4a). This non-specific amplification will not impact species identification as the largest product size for primers designed in this study is 380 bp.

Fig. 4
figure 4

Gel electrophoresis of PCR products demonstrating species-specific amplification of targeted species: a Saccostrea multiplex 1 PCR products of targeted Saccostrea species DNA amplified at an annealing temperature of 62 °C; b Saccostrea multiplex 1 PCR products of non-targeted Saccostrea species DNA (same PCR as 3a but run in a different gel); and c Saccostrea multiplex 2 PCR products of targeted Saccostrea species DNA amplified at an annealing temperature of 62 °C. NC = negative control. See Supplementary Figures S4 – S6 for original photographs

Spat identification

Three species of Saccostrea were identified from gabion baskets deployed in Moreton Bay.

Almost all spat were identified as S. glomerata (n = 63) and only two lineage B and one lineage G were identified. Out of the 75 samples we collected, only nine did not amplify. These nine samples, as well as five of the samples identified as S. glomerata and the lineage B and G samples were then barcoded via 16S Sanger sequencing. This confirmed that identifications made via multiplex PCRs were correct. Of the nine samples for which multiplex PCRs failed, five were also unable to be identified via sequencing. This was likely due to either inadequate amounts of template DNA or the presence of DNA inhibitors which is a common issue when working with mollusc DNA [37,38,39], as some of the bands produced for sequencing were weak. Three were successfully identified as S. glomerata, while the fourth sample was identified as a species belonging to the genus Nanostrea. Our sequence aligned with Nanostrea fluctigera sequences from Japan [40], however the identification of these species was based on morphology, which we know to be unreliable [10]. Given that this species was originally described from the Red Sea [41], it is unlikely that the Nanostrea species identified in this study and by Hamaguchi et al., 2017 are Nanostrea fluctigera. We have therefore given this species the temporary name Nanostrea QLD sp. 1. The 16S sequence generated for this specimen has been submitted to GenBank (accession number PP116107). To our knowledge this species has not previously been documented in Moreton Bay. All other accession numbers for generated sequences can be found in “Availability of data and materials”.

Non-destructive DNA sampling

Out of the 104 oysters collected from Noosa, only three mortalities were recorded post pallial fluid collection after six weeks. Using the S. glomerata and lineage G primer sets, we were able to positively identify 33 live lineage G oysters which have been retained at Bribie Island Research Centre for future spawning trials. 64 oysters were identified as S. glomerata and the remaining 7 samples did not amplify and were therefore unable to be identified. Despite some samples having yields less than 5 ng/ μL, including one sample which was less than 1 ng/ μL, DNA successfully amplified and subsequent identifications were able to be made. Unamplified samples were therefore likely due to the presence of DNA inhibitors, or could have been species other than S. glomerata or lineage G.

Discussion

In this study we successfully developed a series of multiplex assays to facilitate rapid and cost-effective identification of Saccostrea species. Species specificity was achieved for all primer sets at an annealing temperature of 62 °C, except for the lineage B primer set which required a higher annealing temperature of 70 °C. This high temperature is not suitable for multiplex PCRs and will increase false negative lineage B detections. Using Lineage B primers in multiplex PCRs at 62 °C can be problematic, however, as they will amplify lineage I DNA and occasionally lineage G DNA. To deal with this, we recommend using the lineage B and lineage G primer sets together (along with up to two other primer sets) in a multiplex PCR, with any double bands produced at 244 bp and 380 bp identified as lineage G (see Supplementary Figure S7). Although the lineage B primer set amplified lineage I DNA at annealing temperatures below 68 °C, lineage I oysters are found almost exclusively on the aerial roots of Rhizophora mangroves and have several morphologically distinct features, making them one of the only species of Saccostrea that can be visually identified [12, 16] (Fig. 5). We therefore suggest that Lineage B primers can be used in multiplex assays (at annealing temperatures of 62 °C), so long as lineage I oysters are visually identified prior to DNA extractions, and lineage B and G primer sets are used conjunctively in PCR reactions. Given that lineage B and lineage G have closely overlapping distributions in Queensland [12], these primers should always be used together in both subtropical regions along with S. glomerata primers, and in tropical regions along with lineage F and S. spathulata primers. In tropical and subtropical Western Australia, however, all primers need to be tested for non-targeted amplification of lineage A before multiplex assays can be used for reliable identification in this region [15].

Fig. 5
figure 5

Lineage I oysters photographed on the roots and trunks of mangrove trees in Queensland, displaying characteristic broad but delicate purple lobe-like squamae (lamellae) on the lower/bottom (left) valve. A Gladstone; B and C Keppel Islands; and D Turkey Beach. Photo credit: A – C – Marina Richardson; D – Robert Porter

In this study we demonstrated the utility of our developed species-specific primer sets in multiplex assays in both restoration and aquaculture contexts. The results from our spat identifications in Moreton Bay showed that S. glomerata was the most common Saccostrea species at this site. Although lineage B and lineage G were also identified, they were not observed in high enough numbers to be classified as reef-building. S. cuccullata is often cited as a reef-building species in the tropics [22, 42, 43], however reef-building capabilities of specific lineages within this “superspecies” remains unknown. S. glomerata does not occur in the tropics [12], so at least one other ‘S. cuccullata’ lineage must form reefs. Both lineage B and lineage G have widespread distributions across tropical Australia and the Indo-Pacific [7, 11], and it is possible that either (or both) of these lineages could be reef-building in these areas. The development of these multiplex assays will enable the identification of large numbers of oysters in future studies to determine this rapidly and cost-effectively.

The distributions of multiple oyster species should also be monitored to assess whether there will be range shifts in response to climate change, and what potential problems or opportunities may arise. Southern shifts into colder waters have already been documented for S. glomerata [44], and as the effects of climate change are intensified in subtropical regions it is possible that S. glomerata reefs could eventually be replaced by tropical reef-building species [45,46,47]. The development of these multiplex assays will be an essential tool for future monitoring as the visual misidentification of oysters could lead to the mistaken detection of range shifts. For example, there are existing records for S. glomerata in far north Queensland and the Indo-Pacific [22, 23], however these specimens have almost certainly been misidentified, given that surveys based on genetic identification limit S. glomerata to the subtropics [7, 12]. Future assessments of range shifts should consider the robustness of species records that are based on morphology alone.

In this study we also demonstrate how the development of these multiplex assays coupled with successful non-destructive DNA sampling methods have made it possible to identify live oysters. Sanger sequencing for unambiguous identification generally requires higher DNA input for reliable reads, which can be challenging for non-destructive sampling methods where yields can be too low. In most cases we were still able to identify oysters by amplifying DNA with yields as low as 1 – 5 ng/ μL using our primer sets. The positive identification of lineage G has been a breakthrough for aquaculture, and lineage G spawning runs and grow out trials are currently underway at the Bribie Island Research Centre (Fig. 6). With no pre-existing information on the spawning or seasonality of lineage G, future research investigating the biology of the species including spawning triggers, larval development, settlement cues and substrates, growth, susceptibility to disease, potential hybridisation and preferred temperature and salinity will be essential in assessing their performance and potential for aquaculture and/or shellfish reef restoration.

Fig. 6
figure 6

A Lineage G samples positively identified non-destructively; B a close-up image of a positively identified lineage G oyster; and C lineage G spat following a successful spawning trial. Photo credit: Max Wingfield

The COI marker gene has become one of the most popular genes for species identifications and metabarcoding [48, 49]. Although sequencing gives more reliable results, sequencing of this region of the gene is notoriously difficult in invertebrates and especially bivalves, and often PCR bands are not obtained, or sequencing results from obtained bands are poor [37,38,39, 50, 51]. The multiplex primer assays designed in this study enables species identification in this eventuality, and identifications can be undertaken in a standard laboratory without the additional cost of capillary sequencing. The methods developed in this study will be most beneficial when species identifications are required for large numbers of individual specimens, however, we recommend that a subset is sequenced to confirm accuracy. The successful development of multiplex species-specific primer sets depends on the availability of diagnostic single nucleotide polymorphism [25]. Although the COI/cox1 gene is highly variable, available sequences for Saccostrea species and lineages are generally less than 600 bp long, making the development of multiple species-specific primers of unique product lengths challenging. The diversity of Saccostrea species and lineages in the Indo-Pacific is even higher than the diversity in Australia [6, 11]. Generating longer sequences for individual Saccostrea species and lineages occurring throughout region will therefore be instrumental to the successful development of future species-specific primers.

Conclusion

The development of multiplex primer assays will enable the rapid, cost-effective, and accurate identification of Saccostrea species and lineages across tropical Australia. All tropical species and lineages of Saccostrea occurring in Queensland also have broad distributions across the Indo-Pacific [7, 10, 11, 40]. The assays developed in this study present easily replicable methods for the development of additional species-specific primer sets across the tropics, and future studies can therefore build on this work by testing these primers on additional Saccostrea species lineages found in the Indo-Pacific. Visual misidentifications and the use of incorrect species names have perpetuated a state of confusion for understanding the true diversity, ecology, and distribution of tropical oysters. Incorporating these primers into future studies will therefore ensure that species are correctly identified, aiding in the accurate collation of species-specific information to advance both restoration and aquaculture in tropical and subtropical regions.

Availability of data and materials

The sequences generated during this study are available at GenBank. The accession numbers are: Nanostrea sp1. QLD = PP116107; Saccostrea glomerata = PP261151, PP261152, PP261153, PP261154, PP261155, PP261156, PP261157, PP261158; Saccostrea sp. non-mordax lineage B = PP261159, PP261160; and Saccostrea sp. non-mordax lineage G = PP261161.

References

  1. Beck MW, Brumbaugh RD, Airoldi L, Carranza A, Coen LD, Crawford C, et al. Oyster reefs at risk and recommendations for conservation, restoration, and management. BioScience. 2011;61:107–16.

    Article  Google Scholar 

  2. Coen LD, Brumbaugh RD, Bushek D, Grizzle R, Luckenbach MW, Posey MH, et al. Ecosystem services related to oyster restoration. Mar Ecol Prog Ser. 2007;341:303–7.

    Article  Google Scholar 

  3. Grabowski JH, Brumbaugh RD, Conrad RF, Keeler AG, Opaluch JJ, Peterson CH, et al. Economic Valuation of ecosystem services provided by oyster reefs. BioScience. 2012;62:900–9.

    Article  Google Scholar 

  4. Humphries AT, La Peyre MK. Oyster reef restoration supports increased nekton biomass and potential commercial fishery value. PeerJ. 2015;3:e1111.

    Article  PubMed  PubMed Central  Google Scholar 

  5. zu Ermgassen PSE, Hancock B, DeAngelis B, Greene J, Schuster E, Spalding M, et al. Setting objectives for oyster habitat restoration using ecosystem services: A manager’s guide. Arlington VA 76pp: The Nature Conservancy; 2016.

    Google Scholar 

  6. Richardson MA, Zhang Y, Connolly RM, Gillies CL, McDougall C. Some Like it Hot: The ecology, ecosystem benefits and restoration potential of oyster reefs in tropical waters. Front Mar Sci. 2022;9:873768.

    Article  Google Scholar 

  7. Lam K, Morton B. Morphological and mitochondrial-DNA analysis of the Indo-West Pacific rock oysters (Ostreidae: Saccostrea species). J Mollusc Stud. 2006;72:235–45.

    Article  Google Scholar 

  8. Lam K, Morton B. Oysters (Bivalvia: Ostreidae and Gryphaeidae) recorded from Malaysia and Singapore. Raffles Bull Zool. 2009;57:481–94.

    Google Scholar 

  9. Hamaguchi M, Shimabukuro H, Usuki H, Hori M. Occurrences of the Indo-West Pacific rock oyster Saccostrea cucullata in mainland Japan. Marine Biodivers Rec. 2014;7:e84.

    Article  Google Scholar 

  10. Sekino M, Yamashita H. Mitochondrial and nuclear DNA analyses of Saccostrea oysters in Japan highlight the confused taxonomy of the genus. J Mollusc Stud. 2016;82:492–506.

    Article  Google Scholar 

  11. Guo X, Li C, Wang H, Xu Z. Diversity and evolution of living oysters. J Shellfish Res. 2018;37:755–71.

    Article  Google Scholar 

  12. McDougall C, Nenadic N, Richardson M, Healy JM. Molecular identification of intertidal rock oyster species in North-Eastern Australia reveals new candidates for aquaculture. Aquac. 2024;587:740838.

    Article  CAS  Google Scholar 

  13. Cui Z, Hu L, Li C, Zhang Z, Guo X, Wang H. Identification of Saccostrea mordax and a New Species Saccostrea mordoides sp. nov. (Bivalvia: Ostreidae) from China. J Shellfish Res. 2021;40(1):63–75.

    Article  Google Scholar 

  14. O’Hare JA, Momigliano P, Raftos DA, Stow AJ. Genetic structure and effective population size of Sydney rock oysters in eastern Australia. Conserv Genet. 2021. https://doi.org/10.1007/s10592-021-01343-4.

    Article  Google Scholar 

  15. Snow M, Fotedar S, Wilson NG, Kirkendale LA. Clarifying the natural distribution of Saccostrea (edible rock oyster) species in Western Australia to guide development of a fledgling aquaculture industry. Aquac. 2023;566:739202.

    Article  CAS  Google Scholar 

  16. McDougall C, Nenadic N, Healy J. Guide to Queensland’s intertidal oyster. Molecular identification of intertidal rock oyster species in north-eastern Australia reveals new candidates for aquaculture. 2020.

  17. Thomson JM. The genera of oysters and the Australian species. Aust J Mar Freshw Res. 1953;5(1):132–68.

    Article  Google Scholar 

  18. Braley RD. Mariculture potential of introduced oysters Saccostrea cucullata tuberculata and Crassostrea echinata, and a histoogical study of reproduction of C. echinata. Aust J Mar Freshw Res. 1984;35:129–41.

    Article  Google Scholar 

  19. Southgate PC, Lee PS. Hatchery rearing of the tropical blacklip oyster Saccostrea echinata (Quoy and Gaimard). Aquac. 1998;169:275–81.

    Article  Google Scholar 

  20. Nowland SJ, O’Connor WA, Osborne MWJ, Southgate PC. Current status and potential of tropical rock oyster aquaculture. Rev Fish Sci Aquac. 2019;28:57–70.

    Article  Google Scholar 

  21. Davenport J, Wong TM. Effects of temperature and aerial exposure on three tropical oyster species, Crassostrea belcheri, Crassostrea iradelei and Saccostrea cucullata. J Therm Biol. 1992;17:135–9.

    Article  Google Scholar 

  22. Alleway HK, Crawford C, Creighton C, Diggles B, Ford J, Gillies CL, et al. Shellfish reef habitats: a synopsis to underpin the repair and conservation of Australia’s environmentally, socially and economically important bays and estuaries. Townsville: Centre for Tropical Water and Aquatic Ecosystem Research (TropWATER) Publication, James Cook University; 2015.

    Google Scholar 

  23. do Amaral VS, Simone LRL. Comparative anatomy of five species of Saccostrea Dollfus and Dautzenberg, 1920 (Bivalvia: Ostreidae) from the Pacific Ocean. Nautilus. 2016;130:53–71.

    Google Scholar 

  24. Chowdhury MSN, Wijsman JWM, Shahadat Hossain M, Ysebaert T, Smaal AC. A verified habitat suitability model for the intertidal rock oyster, Saccostrea cucullata. Plos One. 2019;14(6):e0217688.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang H, Guo X. Identification of Crassostrea ariakensis and related oysters by multiplex species-specific PCR. J Shellfish Res. 2008;27:481–7.

    Article  Google Scholar 

  26. Wang T, Quan W, Cheng Q, Fan R. The development of two novel species-specific primers for identifying genus Crassostrea, with focus on C. sikamea and C. ariakensis. Mitochondrial DNA B Resour. 2018;3:892–5.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Wang T, Fan R, Cheng Q, Sun Z, Fan X, Li N, et al. Intertidal zonation of the Suminoe oyster Crassostrea ariakensis and the Kumamoto oyster Crassostrea sikamea on the coast of the northern East China Sea. J Shellfish Res. 2020;39:31–41.

    Article  Google Scholar 

  28. Chowdhury MS, Reece KS, Joy SS, Coen LD, Smaal AC. Morphological and mitochondrial dna analyses of oysters in the northern Bay of Bengal. Bangladesh J Shellfish Res. 2021;40:213–29.

    Google Scholar 

  29. Safak I, Norby PL, Dix N, Grizzle RE, Southwell M, Veenstra JJ, et al. Coupling breakwalls with oyster restoration structures enhances living shoreline performance along energetic shorelines. Ecol Eng. 2020;158.

    Article  Google Scholar 

  30. Walters LJ, Roddenberry A, Crandall C, Wayles J, Donnelly M, Barry SC, et al. The use of non-plastic materials for oyster reef and shoreline restoration: understanding what is needed and where the field is headed. Sustain. 2022;14:8055.

    Article  CAS  Google Scholar 

  31. Diggles BK. Annual pattern of settlement of Sydney rock oyster (Saccostrea glomerata) spat in Pumicestone Passage, Moreton Bay. Proc R Soc Queen. 2017;122:17–33.

    Google Scholar 

  32. Angell CL. The biology and culture of tropical oysters. international center for living aquatic resources management. Manila, Philippines: ICLARM Studies and Reviews; 1986.

    Google Scholar 

  33. Nowland SJ, O’Connor WA, Penny SS, Southgate PC. Monsoonally driven reproduction in the tropical black-lip oyster Saccostrea echinata (Quoy & Gaimard, 1835) in northern Australia. J Shellfish Res. 2019;38:89–100.

    Article  Google Scholar 

  34. Nowland SJ, O’Connor WA, Elizur A, Southgate PC. Evaluating spawning induction methods for the tropical black-lip rock oyster Saccostrea echinata. Aquac Rep. 2021;20:100676.

    Article  Google Scholar 

  35. McDougall C, Walker S. Biosphere to restore a functional estuary: final report. Brisbane: Australian Rivers Institute, Griffith University; 2021.

    Google Scholar 

  36. Armbrecht L, Herrando-Pérez S, Eisenhofer R, Hallegraeff GM, Bolch CJS, Cooper A. An optimized method for the extraction of ancient eukaryote DNA from marine sediments. Mol Ecol Res. 2020;20:906–19.

    Article  CAS  Google Scholar 

  37. Pereira JC, Chaves R, Bastos E, Leitão A, Guedes-Pinto H. An Efficient method for genomic DNA Extraction from Different Molluscs Species. IJMS. 2011;12:8086–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Adema CM. Sticky problems: extraction of nucleic acids from molluscs. Phil Trans R Soc B. 2021;376:20200162.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Layton KKS, Martel AL, Hebert PDn. Patterns of DNA barcode variation in canadian marine molluscs. PLoS ONE. 2014;9:e95003.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Hamaguchi M, Manabe M, Kajihara N, Shimabukuro H, Yamada Y, Nishi E. DNA barcoding of flat oyster species reveals the presence of Ostrea stentina Payraudeau, 1826 (Bivalvia: Ostreidae) in Japan. Mar Biodivers Rec. 2017;10:4.

    Article  Google Scholar 

  41. Blatterer H. Mollusca of the Dahab region (Gulf of Aqaba, Red Sea). Denisia. 2019;43:1–480.

    Google Scholar 

  42. McAfee D, Bishop MJ, Yu T-N, Williams GA. Structural traits dictate abiotic stress amelioration by intertidal oysters. Funct Ecol. 2018;32:2666–77.

    Article  Google Scholar 

  43. Al-Kandari M, Oliver PG, Salvi D. Molecular and morphological systematics of a new, reef forming, cupped oyster from the northern Arabian Gulf: Talonostrea salpinx. ZooKeys. 2021;1043:1–20.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Crawford C. Report on the occurrence of Sydney rock oysters on the east coast of Tasmania. Hobart, Tasmania: Institute for Marine and Antarctic Studies, University of Tasmania; 2015.

    Google Scholar 

  45. Watson S, Southgate PC, Tyler PA, Peck LS. Early larval development of the Sydney rock oyster Saccostrea glomerata under near-future predictions of CO2-driven ocean acidification. J Shellfish Res. 2009;28:431–7.

    Article  Google Scholar 

  46. Parker LM, O’Connor WA, Byrne M, Dove M, Coleman RA, Portner HO, et al. Ocean acidification but not warming alters sex determination in the Sydney rock oyster. Saccostrea glomerata Proc Biol Sci. 1872;2018(285):20172869.

    Google Scholar 

  47. Scanes E, Parker LM, O’Connor WA, Dove MC, Ross PM. Heatwaves alter survival of the Sydney rock oyster. Saccostrea glomerata Mar Pollut Bull. 2020;158:111389.

    Article  CAS  PubMed  Google Scholar 

  48. Porter TM, Hajibabaei M. Over 2.5 million COI sequences in GenBank and growing. PLoS ONE. 2018;13:e0200177.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome C oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3:294–9.

    CAS  PubMed  Google Scholar 

  50. Geller J, Meyer C, Parker M, Hawk H. Redesign of primers for mitochondrial cytochrome oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys. Molecular Ecology Resources. 2013;13:851–61.

    Article  CAS  PubMed  Google Scholar 

  51. Sharma P, Kobayashi T. Are, “universal” DNA primers really universal? J Appl Genetics. 2014;55:485–96.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Thank you to OzFish Unlimited and Robert Porter for assisting in the deployment and collection of samples in Moreton Bay, Aiden Mellor and Luke Dutney at the Department of Agriculture and Fisheries Bribie Island Aquaculture facility for their assistance and contributions collecting and maintaining oysters, and to The Nature Conservancy for funding this research through their Reef Builder program.

Funding

The Nature Conservancy 2020 Reef-Builder program.

Author information

Authors and Affiliations

Authors

Contributions

MR and CM conceptualized the study. MR undertook the primer design and testing, collected and identified all oysters in Moreton Bay, prepared all figures and tables, and wrote the initial draft. NN and CM developed methods for non-destructive DNA extraction, MW collected and maintained oysters in aquaria from Noosa and NN performed species identifications. All authors contributed to the article and approved the submitted version.

Corresponding authors

Correspondence to Marina A. Richardson or Carmel McDougall.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Richardson, M.A., Nenadic, N., Wingfield, M. et al. The development of multiplex PCR assays for the rapid identification of multiple Saccostrea species, and their practical applications in restoration and aquaculture. BMC Ecol Evo 24, 67 (2024). https://doi.org/10.1186/s12862-024-02250-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-024-02250-1

Keywords