2urn:lsid:arphahub.com:pub:45048D35-BB1D-5CE8-9668-537E44BD4C7Eurn:lsid:zoobank.org:pub:91BD42D4-90F1-4B45-9350-EEF175B1727AZooKeysZK1313-29891313-2970Pensoft Publishers10.3897/zookeys.365.59053065Short CommunicationBivalviaMolluscaMolecular systematicsEuropeNorth AmericaUsing DNA barcoding to differentiate invasive Dreissena species (Mollusca, Bivalvia)MarescauxJonathan1jonathan.marescaux@unamur.beVan DoninckKarine1Laboratory of Evolutionary Genetics and Ecology, Research Unit in Environmental and Evolutionary Biology, Department of Biology, University of Namur, 61 rue de Bruxelles, 5000 Namur, BelgiumUniversity of NamurNamurBelgium
Corresponding author: Jonathan Marescaux (jonathan.marescaux@unamur.be)
Academic editor: M. de Meyer
201330122013365235244FFC59237-FE0F-FF9A-FF94-FF9B31462B7457751927201316102013Jonathan Marescaux, Karine Van DoninckThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
The zebra mussel (Dreissena polymorpha) and the quagga mussel (Dreissena rostriformis bugensis) are considered as the most competitive invaders in freshwaters of Europe and North America. Although shell characteristics exist to differentiate both species, phenotypic plasticity in the genus Dreissena does not always allow a clear identification. Therefore, the need to find an accurate identification method is essential. DNA barcoding has been proven to be an adequate procedure to discriminate species. The cytochrome c oxidase subunit I mitochondrial gene (COI) is considered as the standard barcode for animals. We tested the use of this gene as an efficient DNA barcode and found that it allow rapid and accurate identification of adult Dreissena individuals.
Biological invasions are a topical issue in today’s world since they are the biggest threat to biodiversity after habitat destruction. The first, and probably the biggest, problem for scientists is to deal with widely divergent perceptions of the criteria defining “invasive” species (Colautti and MacIsaac 2004). In the management and policy field, such species are defined as “alien species whose introduction does, or is likely to, cause economic or environmental harm or harm to human health” (Invasive Species Advisory Committee 2006). By cons, from a strict scientific point of view, an invasive species is “an exotic species that present a tendency to spread with high densities from its point of introduction” (Vermeij 1996, Beisel and Lévêque 2010). A second problem for both scientists and managers is to rapidly characterize a new invasion.
The zebra mussel (Dreissena polymorpha (Pallas, 1771) and the quagga mussel (Dreissena rostriformis bugensis Andrusov, 1897) are invasive freshwater bivalves in Europe and North America (Mills et al. 1996, Son 2007). Both species are native to the Ponto-Caspian area (Son 2007) and have major negative ecological and economic impacts such as biofouling and food web alteration (Sousa et al. 2013). Several studies have shown that the newly introduced quagga mussel can often dominate well-established zebra mussel populations within only a few years and even outcompete it in some cases (Wilson et al. 2006, Heiler et al. 2012). Wilke et al. (2010) showed that, in addition to the well-known zebra and quagga mussels, two others Dreissena species native to the Balkans (Dreissena presbensis (Kobelt, 1915) and Dreissena blanci Westerlund, 1890) begin to expand in the area and may be potentially invasive in Europe.
Although Dreissena specialists may discriminate adults of the different species based on internal and external shell features (Pathy and Mackie 1993, Mills et al. 1996, Sablon et al. 2010), this task remains difficult for managers. It becomes even more problematic when identifying larvae, which is the most invasive form of Dreissena (Marescaux et al. 2012a, b). For example, the invasion of the Meuse River in Belgium by the quagga mussel remained undetected because Belgian national agencies never made the distinction with the zebra mussel. Therefore, tools for rapid identification of both adult specimens and larvae are needed in order to detect newly invaded habitats. DNA barcoding has been proven to be an effective method both for species detection and to assign new specimens to already identified species (Hebert et al. 2003a, Birky et al. 2010). Here we amplified part of the cytochrome c oxidase subunit I (COI) mitochondrial gene, the most-widely utilized gene for animal DNA barcoding (Consortium for the Barcode of Life 2013) and we tested four delimitation metrics to differentiate Dreissena species. We also demonstrate that restriction fragment length polymorphism (RFLP) could be used as an inexpensive method to distinguish between zebra and quagga mussel.
MethodsSamples collection
Dreissena samples were collected in the Meuse River (see Marescaux et al. 2012a, b for sampling protocol and locations). The mussels were collected in the littoral zone of the river bank from stones which were picked up manually from a depth of 30–40 cm.
COI sequencing
Total genomic DNA was extracted from 241 Dreissena individuals using the «DNeasy Blood and Tissue» kit (Qiagen) according to manufacturer guidelines. To minimize cost, DNA extraction with the CTAB (hexadecyltrimethylamoniumbromide) protocol proposed by Winnepenninckx et al. (1993) could also be used. A fragment of 654 base pairs (bp) of the COI mitochondrial gene was amplified using universal primers (Folmer et al. 1994). Amplifications were performed in 25 μl total volume including 0.5 or 1 μl of gDNA, 1 × GoTaq Green reaction buffer (Promega), 200 μM of dNTPs (Promega), 0.5 μM of both primers and 0.1 U of GoTaq DNA polymerase (Promega). PCR cycling conditions were as follows: an initial step of 94 °C for 4 min, followed by 30 cycles of 94 °C for 45 s, 45 °C for 45 s and 72 °C for 45 s, and then a final extension of 72 °C for 10 min. DNA sequencing was performed by the Genoscreen Company (France). Sequences were visualized and aligned using BioEdit v7.0.5.3 (Hall 1998).
Phylogenetic analysis
Sequences were collapsed into unique haplotypes using DnaSP (Librado and Rozas 2009). In order to determine the number of Dreissena species in the Meuse River we tested three barcoding methods: (i) the “Operational Taxonomic Units” (OTU) (Hebert et al. 2003a), (ii) the “Automatic Barcode Gap Discovery” (ABGD) (Puillandre et al. 2012), and (iii) the “K/θ method” (4 × rule) (Birky et al. 2010). The K/θ method specifies that if the genetic distance between clusters is higher than 4 times the genetic distance within the cluster then species are distinct (Birky et al. 2010, Tang et al. 2012). Neighbour-Joining (NJ) trees and matrix of pairwise distances were calculated using the Kimura 2-parameter (K2P) model and were generated using MEGA4 in order to define OTU’s (Tamura et al. 2007). Sequences found in GenBank (Table 1) were used to construct a haplotype network using Network v4.6 (Bandelt et al. 1999).
GenBank accession numbers and localities of Dreissena spp. sequences included in the network analysis.
GenBank
Taxon
Location
DQ840122
Dreissena polymorpha polymorpha
Black and Caspian Seas
DQ840125
Dreissena polymorpha polymorpha
Liman, Caspian Sea
DQ840123
Dreissena polymorpha polymorpha
Caspian Sea
DQ840121
Dreissena polymorpha polymorpha
Black and Caspian Seas
EF414493
Dreissena polymorpha
Turkey
U47653
Dreissena polymorpha
Lake Ontario
AF474404
Dreissena polymorpha
Poland
EU484441
Dreissena polymorpha
Lake Superior
EU484437
Dreissena polymorpha
Lake Superior
EU484448
Dreissena polymorpha
Lake Superior
EU484444
Dreissena polymorpha
Lake Superior
AM748997
Dreissena polymorpha
Italy
AM748986
Dreissena polymorpha
Germany
AM748977
Dreissena polymorpha
Italy
U47651
Dreissena bugensis
Lake Ontario
U47650
Dreissena bugensis var. profunda
Lake Ontario
DQ840132
Dreissena bugensis
Black Sea
EF080861
Dreissena rostriformis bugensis
Hollandsch Diep
AF495877
Dreissena bugensis
Ukraine
AF479637
Dreissena bugensis
Ukraine
AM748999
Dreissena polymorpha
Germany
Restriction fragment length polymorphism analysis (RFLP)
Using the restriction map application (http://www.bioinformatics.org/sms2/rest_map.html), we selected two endonucleases to differentially cut the COI gene of Dreissena species: Hinf I and Nla III. We also tested two other enzymes used in previous studies: Nla IV (Baldwin et al. 1996) and Scr FI (Claxton et al. 1998).
Restriction analysis of the amplified 654 bp COI fragment was carried out on each dreissenid haplotype (using individuals from the Meuse River). For each haplotype, the RFLP was performed in 31 μl total volume including 10 µl of PCR reaction mixture, 18 µl of distilled water, 2 µl of buffer (supplied by the manufacturer with the enzyme), and 1 µl of enzyme. Digests were incubated at 37 °C for 3 hours and then loaded on 2% agarose gels.
Results
Sequencing of the 654 bp COI fragment revealed seven haplotypes among the 241 Dreissena individuals. The OTU method revealed, by a NJ tree, two clusters separated by a genetic distance of 18.5% (Figure 1a), which is higher than the 3% threshold typically used for species delimitation with COI (Hebert et al. 2003b). This first analysis, therefore, suggests the occurrence of two species. We obtained the same results with the ABGD method. Indeed, the K2P-distances show two distinct clusters (Figure 1b). One cluster formed by haplotype 1 and 2, and a second cluster containing the five other haplotypes, all corresponding to those separated in the tree. Moreover, the genetic distances within our two clusters (0.6% and 0.2%, respectively) are four times lower than the genetic distance between them (18.5%) (Figure 1a) confirming the presence of two Dreissena species.
Barcoding analysis based on a fragment of 654 base pairs of the COI gene. a) NJ analysis of K2P-pairwise distances b) “barcoding gap” method based on the K2P-pairwise distance.
https://binary.pensoft.net/fig/20686
Our network (Figure 2) revealed that haplotypes 1 and 2 (Q1 and Q2) cluster with Dreissena rostriformis bugensis and the five other haplotypes (Z1 to Z5) cluster with Dreissena polymorpha. This, together with the three barcoding methods which each identified two clusters, shows that both Dreissena polymorpha and Dreissena rostriformis bugensis species occur in the Meuse River.
Haplotype networks based on a fragment of 654 base pairs of the COI gene. Our seven haplotypes are labelled: Q1 and Q2 for haplotypes 1 and 2 (belonging to Dreissena rostriformis bugensis) / Z1 to Z5 for the 5 other haplotypes (belonging to Dreissena polymorpha).
https://binary.pensoft.net/fig/20687
Digestion profiles for each haplotype are illustrated in Figure 3. Each of the four endonucleases tested, yielded distinct restriction patterns between both Dreissena species. Digestion with Nla IV produced four fragments in quagga mussels (Q haplotype) of approximately 70, 79, 211, and 294 bp and three distinct patterns for the zebra mussel (Z haplotype): haplotype Z1 and Z2 (91, 120, 150, and 293 bp), haplotype Z3 and Z4 (91, 150, and 413 bp), and haplotype Z5 (91, 150, 200, and 413 bp). We suggest here that the 200 bp fragment of the haplotype Z5 is an artefact, as confirmed by the restriction map, since the summed fragment lengths do not add up to the expected 654 bp. We infer that haplotype Z5 has the same pattern as haplotype Z3 and Z4. Digestion with Hinf I produced two fragments in quagga mussels of approximately 73 and 581 bp and five fragments in zebra mussels of approximately 31, 101, 114, 195, and 213 bp. The small fragments can not be distinguished on the gel but the difference between quagga and zebra is clear. Digestion with Nla III produced two fragments in quagga mussels of approximately 193 and 461 bp and three fragments in zebra mussels of approximately 193, 319, and 335 bp. Digestion with Scr FI produced five fragments in quagga mussels of approximately 42, 53, 120, 171, and 268 bp and three fragments in zebra mussels of approximately 95, 152, and 407 bp. The digestion pattern for the quagga mussel using the endonuclease Scr FI is not clearly defined (smear) since the five fragments are very short.
RFLP analysis of the COI gene to distinguish Dreissena rostriformis bugensis (Q haplotype) and Dreissena polymorpha (Z haplotype) using the endonucleases (A) Nla IV (B) Hinf I (C) Nla III and (D) Scr FI. Lane 1, 1-kb ladder; lane 2, non-digested fragment of quagga mussel; lane 3, Q1 haplotype; lane 4, Q2 haplotype; lane 5, Z1 haplotype; lane 6, Z2 haplotype; lane 7, Z3 haplotype; lane 8, Z4 haplotype; lane 9, Z5 haplotype; lane 10, 100-bp ladder.
https://binary.pensoft.net/fig/20688Discussion
On September 9 2013, the European Commission has published a proposal for a Regulation on the prevention and management of the introduction and spread of invasive alien species. This proposal highlights three types of interventions: prevention, early warning and rapid response, and then management of invasive species (European Commission 2013). In this context, rapid identification methods are needed to detect invasive species in periodic surveys, e.g. inspection of ballast water. We showed in previous work (Marescaux et al. 2012a, b) that visual identification and morphometric analyses are not always sufficient to differentiate both zebra and quagga mussel probably due to phenotypic plasticity. This is particularly true for larval identification. In addition, two other Dreissena species may become invasive and should be detected promptly.
In order to help managers and national agencies, we propose here the use of the COI mitochondrial gene as a barcode to discriminate Dreissena polymorpha and Dreissena rostriformis bugensis. Moreover, it is possible to conduct a RFLP analysis on this gene to obtain results without sequencing cost. This method could also easily be applied to Dreissena presbensis and Dreissena blanci since the COIgene have already been sequenced by Albrecht et al. (2007) and Wilke et al. (2010) and sequences are available on GenBank (accession numbers EF414478–EF414492, EF414496, HM209829–HM210081). We showed that the endonuclease Nla IV, previously used by Baldwin et al. (1996), presents different restriction patterns for the zebra mussel haplotype and not a clear distinction between some zebra mussel haplotypes (Z1 and Z2) and the quagga mussel haplotypes. Therefore, we do not recommend the use of this enzyme to discriminate between quagga and zebra mussel. The three other endonucleases tested during this study present a clear distinction between both species despite the fact that a smear appears using endonucleases Hinf I and Scr FI. Moreover, Nla III and Scr FI will produce a unique RFLP banding pattern for Dreissena blanci and Dreissena presbensis different from those observed in the zebra and quagga mussel.
This study is the first step of an extensive phylogeographical analysis on the invasion of Western Europe by the dreissenids. Further experiments will be needed to assess potential risks of both zebra and quagga mussels on native biodiversity in Western European rivers, e.g. predation on phytoplankton, infestation on native bivalves and alteration of macro-invertebrate communities.
Acknowledgements
Special thanks to Emilie Etoundi and Doctor Xiang Li for the help with the delimitation metrics. We also thank two reviewers and the editor for their helpful comments and critical reading of this manuscript. This study received financial support from the University of Namur. Jonathan Marescaux is funded by a PhD grant from the Belgian National Fund for Scientific Research (FRS-FNRS).
ReferencesAlbrechtCSchultheißRKevrekidisTStreitBWilkeT (2007) Invaders or endemics? Molecular phylogenetics, biogeography and systematics of Dreissena in the Balkans.52: 1525-1536. doi: 10.1111/j.1365-2427.2007.01784.xBaldwinBSBlackMSanjurOGustafsonRLutzRAVrijenhoekRC (1996) A diagnostic molecular marker for zebra mussels (Dreissena polymorpha) and potentially co-occurring bivalves: mitochondrial COI.5: 9-14.BandeltHJForsterPRöhlA (1999) Median-joining networks for inferring intraspecific phylogenies.16: 37-48. doi: 10.1093/oxfordjournals.molbev.a026036BeiselJNLévêqueC (2010) Introductions d’espèces dans les milieux aquatiques: Faut-il avoir peur des invasions biologiques?, 232 pp.BirkyCWAdamsJGemmelMPerryJ (2010) Using population genetic theory and DNA sequences for species detection and identification in asexual organisms.5: e10609. doi: 10.1371/journal.pone.0010609ClaxtonWTBouldingEG (1998) A new molecular technique for identifying field collections of zebra mussel (Dreissena polymorpha) and quagga mussel (Dreissena bugensis) veliger larvae applied to eastern Lake Erie, Lake Ontario, and Lake Simcoe.76: 194-198.ColauttiRIMacIsaacHJ (2004) A neutral terminology to define ‘invasive species’.10: 135-141. doi: 10.1111/j.1366-9516.2004.00061.xConsortium for the Barcode of Life (2013) . http://www.barcodeoflife.org/European Commission (2013) http://ec.europa.eu/environment/nature/invasivealien/index_en.htmFolmerOBlackMHoehWLutzRVrijenhoekR (1994) DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates.3: 294-299.HallTA (1998) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.41: 95-98.HebertPDNCywinskaABallSLdeWaardJR (2003a) Biological identifications through DNA barcodes.270: 313-321. doi: 10.1098/rspb.2002.2218HebertPDNRatnasinghamSdeWaardJR (2003b) Barcoding animal life: cytochrome c oxidase subunit I divergences among closely related species.(Supplement)270: S96-S99. doi: 10.1098/rsbl.2003.0025HeilerKCMBrandtSAlbrectCHauffeTWilkeT (2012) A new approach for dating introduction events of the quagga mussel (Dreissena rostriformis bugensis).14: 1311-1316. doi: 10.1007/s10530-011-0161-1Invasive Species Advisory Committee (2006) Invasive species definition clarification and guidance white paper.Washington, D.C.http://www.invasivespecies.gov/global/ISAC/ISAC_documents/ISAC%20Definititions%20White%20Paper%20%20-%20FINAL%20VERSION.pdfLibradoPRozasJ (2009) DnaSP v5: A software for comprehensive analysis of DNA polymorphism data.25: 1451-1452. doi: 10.1093/bioinformatics/btp187MarescauxJBij de VaateAVan DoninckK (2012a) First records of Dreissena rostriformis bugensis (Andrusov, 1897) in the Meuse River.1: 119-124. doi: 10.3391/bir.2012.1.2.05MarescauxJMolloyDPGiamberiniLAlbrechtCVan DoninckK (2012b) First records of the quagga mussel, Dreissena rostriformis bugensis (Andrusov, 1897), in the Meuse River within France.1: 273-276. doi: 10.3391/bir.2012.1.4.05MillsELRosenbergGSpidleAPLudyanskiyMPliginYMayB (1996) A review of the biology and ecology of the quagga mussel (Dreissena bugensis), a second species of freshwater dreissenid introduced to North America.36: 271-286. doi: 10.1093/icb/36.3.271PathyDAMackieGL (1993) Comparative shell morphology of Dreissena polymorpha, Mytilopsis leucophaeata and the “quagga” mussel (Bivalvia: Dreissenidae) in North America.71: 1012-1023. doi: 10.1139/z93-135PuillandreNLambertABrouilletSAchazG (2012) ABGD, Automatic Barcode Gap Discovery for primary species delimitation.21: 1864-1877. doi: 10.1111/j.1365-294X.2011.05239.xSablonRVercauterenTJacobsP (2010) De quaggamossel (Dreissena rostriformis bugensis (Andrusov, 1897)), een recent gevonden invasieve zoetwatermossel in Vlaanderen.4: 32-36.SonMO (2007) Native range of the zebra mussel and quagga mussel and new data on their invasions within the Ponto-Caspian Region.2: 174-184. doi: 10.3391/ai.2007.2.3.4SousaRNovaisACostaRStrayerDL (2013) Invasive bivalves in fresh waters: impacts from individuals to ecosystems and possible control strategies. doi: 10.1007/s10750-012-1409-1TamuraKDudleyJNeiMKumarS (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.24: 1596-1599. doi: 10.1093/molbev/msm092TangCQLeasiFOberteggerUKienekeABarracloughTGFontanetoD (2012) The widely used small subunit 18S rDNA molecule greatly underestimates true diversity in biodiversity surveys of the meiofauna.109: 16208-16212. doi: 10.1073/pnas.1209160109VermeijGJ (1996) An agenda for invasion biology.78: 3-9. doi: 10.1016/0006-3207(96)00013-4WilkeTSchultheißRAlbrechtCBornmannNTrajanovskiSKevrekidisT (2010) Native Dreissena freshwater mussels in the Balkans: in and out of ancient lakes.7: 3051–3065. doi: 10.5194/bg-7-3051-2010WilsonKAHowellETJacksonDA (2006) Replacement of zebra mussels by quagga mussels in the Canadian nearshore of Lake Ontario: the importance of substrate, round goby abundance, and upwelling frequency.32: 11-28. doi: 10.3394/0380-1330(2006)32[11:ROZMBQ]2.0.CO;2WinnepenninckxBBackeljauTDe WachterR (1993) Extraction of high molecular weight DNA from molluscs.9: 407.