ABSTRACT
BACKGROUND: Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution. RESULTS: At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths. CONCLUSIONS: Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.
Subject(s)
Microsporidia , Microsporidia/genetics , Phylogeny , Evolution, Molecular , Genome , Introns , Eukaryota/geneticsABSTRACT
Environmental DNA (eDNA) metabarcoding can rapidly characterize the composition and diversity of benthic communities, thus it has high potential utility for routine assessments of benthic impacts of marine finfish farming. In this study, 126 sediment grab samples from 42 stations were collected at six salmon farms in British Columbia, Canada. Benthic community changes were assessed by both eDNA metabarcoding of metazoans and macrofaunal polychaete surveys. The latter was done by analysing 11,466 individuals using a combination of morphology-based taxonomy and DNA barcoding. Study objectives were to: (i) compare biotic signals associated with benthic impacts of salmon farming in the two data sources, and (ii) identify potential eDNA indicators to facilitate monitoring in Canada. Alpha diversity parameters were consistently reduced near fish cage edge and negatively correlated with pore-water sulphide concentration, with coefficients ranging from -0.62 to -0.48. Although Polychaeta are a common indicator group, the negative correlation with pore-water sulphide concentration was much stronger for Nematoda OTU richness (correlation coefficient: -0.86) than for Polychaeta (correlation coefficient: -0.38). Presence/absence of Capitella generally agreed well between the two methods despite that they differed in the volume of sediments sampled and the molecular marker used. Multiple approaches were used to identify OTUs related to organic enrichment statuses. We demonstrate that eDNA metabarcoding generates biotic signals that could be leveraged for environmental assessment of benthic impacts of fish farms in multiple ways: both alpha diversity and Nematoda OTU richness could be used to assess the spatial extent of impact, and OTUs related to organic enrichment could be used to develop local biotic indices.
Subject(s)
DNA Barcoding, Taxonomic , Salmon , Animals , Aquaculture , Biodiversity , British Columbia , Environmental Monitoring , Geologic Sediments , Humans , Salmon/geneticsABSTRACT
In a recent paper, "Environmental DNA: What's behind the term? Clarifying the terminology and recommendations for its future use in biomonitoring," Pawlowski et al. argue that the term eDNA should be used to refer to the pool of DNA isolated from environmental samples, as opposed to only extra-organismal DNA from macro-organisms. We agree with this view. However, we are concerned that their proposed two-level terminology specifying sampling environment and targeted taxa is overly simplistic and might hinder rather than improve clear communication about environmental DNA and its use in biomonitoring. This terminology is based on categories that are often difficult to assign and uninformative, and it overlooks a fundamental distinction within eDNA: the type of DNA (organismal or extra-organismal) from which ecological interpretations are derived.
Subject(s)
DNA, Environmental , Biodiversity , DNA/genetics , DNA Barcoding, TaxonomicABSTRACT
Latent class analysis (LCA) is a common method to evaluate the diagnostic sensitivity (DSe) and specificity (DSp) for pathogen detection assays in the absence of a perfect reference standard. Here we used LCA to evaluate the diagnostic accuracy of 3 tests for the detection of Mikrocytos mackini in Pacific oysters Crassostrea gigas: conventional polymerase chain reaction (PCR), real-time quantitative PCR (qPCR), and histopathology. A total of 802 Pacific oysters collected over 12 sampling events from 9 locations were assessed. Preliminary investigations indicated that standard LCA assumptions of test independence and constant detection accuracy across locations were likely unrealistic. This was mitigated by restructuring the LCA in a Bayesian framework to include test-derived knowledge about pathogen prevalence and load for categorizing populations into 2 classes of infection severity (low or high) and assessing separate DSe and DSp estimates for each class. Median DSp estimates were high (>96%) for all 3 tests in both population classes. DSe estimates varied between tests and population classes but were consistently highest for qPCR (87-99%) and lowest for histopathology (21-51%). Acknowledging that detection of M. mackini may be fitted to multiple diagnostic and management purposes, qPCR had the highest DSe while maintaining similar DSp to both conventional PCR and histopathology and thus is generally well-suited to most applications.
Subject(s)
Crassostrea , Animals , Bayes Theorem , Latent Class Analysis , Real-Time Polymerase Chain Reaction/veterinary , Sensitivity and SpecificityABSTRACT
Non-indigenous species (NIS) reach every corner of the world, at times wreaking havoc on ecosystems and costing the global economy billions of dollars. A rapid and accurate biosurveillance tool tailored to a particular biogeographic region is needed to detect NIS when they are first introduced into an area as traditional detection methods are expensive and require specialized expertise. Metabarcoding of environmental and community DNA meets those biosurveillance requirements; a novel tool tailored to the Northwest Pacific Ocean is presented here using an approach that could revolutionize early detection of NIS. Eight newly designed genetic markers for multiple gene regions were implemented to meet the stringent taxonomic requirements for the detection of NIS across four major marine phyla. The tool was considered highly successful because it identified 12 known NIS in the study area and a further seven species representing potential new records. Overall community composition detected here was statistically different between substrate types; zooplankton sampling accounted for significantly higher species richness than filtered sea water in most cases, but this was dominated by mollusk and arthropod species. Both substrate types sampled were required to identify the wide taxonomic breadth of known NIS in the study area. Intensive sampling is known to be paramount for the detection of rare species, including new incursions of NIS, thus it is recommended to include diverse DNA sampling protocols based on species' life-history characteristics for broad detection capacity. Application of a metabarcoding-based molecular biosurveillance tool optimized for biogeographic regions enables rapid and accurate early detection across a wide taxonomic range to allow quick implementation of eradication or control efforts and potentially mitigate some of the devastating effects of NIS worldwide.
Subject(s)
Biosurveillance , Introduced Species , Animals , Biodiversity , DNA , DNA Barcoding, Taxonomic , Ecosystem , Pacific OceanABSTRACT
Current monitoring methods to assess benthic impacts of marine finfish aquaculture are based on complex biological indices and/or geochemistry data. The former requires benthic macrofauna morpho-taxonomic characterization that is time- and cost-intensive, while the latter provides rapid assessment of the organic enrichment status of sediments but does not directly measure biotic impacts. In this study, sediment samples were collected from seven stations at six salmon farms in British Columbia, Canada, and analyzed for geochemical parameters and by eDNA metabarcoding to investigate linkages between geochemistry and foraminifera. Sediment texture across farm sites ranged from sand to silty loam, while the maximum sediment pore-water sulphide concentration at each site ranged from 1,000 to 13,000 µM. Foraminifera alpha diversity generally increased with distance from cage edge. Adonis analyses revealed that farm site explained the most variation in foraminifera community, followed by sediment type, enrichment status, and distance from cage edge. Farm-specific responses were observed in diversity analyses, taxonomic difference analyses, and correlation analyses. Results demonstrated that species diversity and composition of foraminifera characterized by eDNA metabarcoding generated signals consistent with benthic biodiversity being impacted by finfish farming activities. This substantiates the validity of eDNA metabarcoding for augmenting current approaches to benthic impact assessments by providing more cost-effective and practicable biotic measures than traditional morpho-taxonomy. To capitalize on this potential, further work is needed to design a new nomogram that combines eDNA metabarcoding data and geochemistry data to enable accurate monitoring of benthic impacts of fish farming in a time- and cost-efficient way.
Subject(s)
DNA Barcoding, Taxonomic , DNA/genetics , Foraminifera/genetics , Salmon/genetics , Animals , Aquaculture , Biodiversity , British Columbia , Environmental Monitoring , Fisheries , Geologic Sediments/chemistry , Salmon/growth & developmentABSTRACT
DNA barcoding has been used successfully for identifying specimens belonging to marine planktonic groups. However, the ability to delineate species within taxonomically diverse and widely distributed marine groups, such as the Copepoda and Thecostraca, remains largely untested. We investigate whether a cytochrome c oxidase subunit I (COI-5P) global pairwise sequence divergence threshold exists between intraspecific and interspecific divergences in the copepods plus the thecostracans (barnacles and allies). Using publicly accessible sequence data, we applied a graphical method to determine an optimal threshold value. With these thresholds, and using a newly generated planktonic marine data set, we quantify the degree of concordance using a bidirectional analysis and discuss different analytical methods for sequence-based species delimitation (e.g., BIN, ABGD, jMOTU, UPARSE, Mothur, PTP, and GMYC). Our results support a COI-5P threshold between 2.1% and 2.6% p-distance across methods for these crustacean taxa, yielding molecular groupings largely concordant with traditional, morphologically defined species. The adoption of internal methods for clustering verification enables rapid biodiversity studies and the exploration of unknown faunas using DNA barcoding. The approaches taken here for concordance assessment also provide a more quantitative comparison of clustering results (as contrasted with "success/failure" of barcoding), and we recommend their further consideration for barcoding studies.
Subject(s)
Copepoda/classification , Copepoda/genetics , DNA Barcoding, Taxonomic , Animals , Biodiversity , Canada , Cluster Analysis , Copepoda/anatomy & histology , Genetic Variation , Geography , Phenotype , PhylogenyABSTRACT
Mikrocytos mackini is an intracellular parasite of oysters and causative agent of Denman Island disease in Pacific oysters Crassostrea gigas. Although M. mackini has been investigated for decades, its natural mode of transmission, mechanism for host entry, and environmental stability are largely unknown. We explored these biological characteristics of M. mackini using a recently described quantitative PCR (qPCR) assay. We detected M. mackini in the flow-through tank water of experimentally infected oysters and during disease remission in host tissues following 6 wk of elevated water temperature. Waterborne exposure of oysters to M. mackini further confirmed the potential for extracellular seawater transmission of this parasite and also identified host gill to have the highest early and continued prevalence for M. mackini DNA compared to stomach, mantle, labial palps, or adductor muscle samples. However, infections following waterborne challenge were slow to develop despite a substantial exposure (>106 M. mackini l-1 for 24 h), and further investigation demonstrated that M. mackini occurrence and infectivity severely declined following extracellular seawater incubation of more than 24 h. This study demonstrates a potential for using qPCR to monitor M. mackini in wild or farmed oyster populations during periods of disease remission or from environmental seawater samples. This work also suggests that gill tissues may provide a primary site for waterborne entry and possibly shedding of M. mackini in oysters. Further, although extracellular seawater transmission of M. mackini was possible, poor environmental stability and infection efficiency likely restricts the geographic transmission of M. mackini between oysters in natural environs and may help to explain localized areas of infection.
Subject(s)
Crassostrea/parasitology , Eukaryota/isolation & purification , Polymerase Chain Reaction/methods , Seawater/parasitology , Animals , Filtration , Host-Parasite InteractionsABSTRACT
Potato wart, caused by the fungal pathogen Synchytrium endobioticum, is a serious disease with the potential to cause significant economic damage. The small subunit (SSU) and internal transcribed spacer (ITS) ribosomal DNA (rDNA) were sequenced for several Synchytrium spp., showing a high rate of variability for both of these markers among the different species and monophyly of the genus within phylum Chytridiomycota. The intergenic nontranscribed spacer (IGS) of rDNA was sequenced for different pathotypes and showed no intraspecific variation within S. endobioticum, similar to the other rDNA markers from this study. To facilitate screening for the pathogen in soil, three TaqMan polymerase chain reaction (PCR) assays were developed from SSU, ITS, and IGS rDNA sequences to detect S. endobioticum sporangia in the chloroform-flotation fraction of sieved soil extracts. In the screening portion of the method, a first TaqMan assay targeting the SSU rDNA was developed with positive results that were further confirmed with amplicon melt analysis. A synthetic reaction control cloned into a plasmid was incorporated into the procedure, facilitating the validation of negative results. The presence of the reaction control did not adversely affect the efficiency of the SSU target amplification. A second TaqMan assay targeting the ITS-1 region was developed as a confirmatory test. There was 100% accordance between the SSU and ITS-1 TaqMan assays. Utilizing these two assays in tandem achieved good specificity for S. endobioticum, generating negative results with the cloned SSU and ITS-1 regions from all 14 other Synchytrium spp. considered. Spike recovery experiments indicated that these assays, targeting the SSU and ITS-1 rDNA regions, developed from a phylogeny dataset of the genus, could reliably detect a single sporangium in the chloroform flotation fraction of a soil extract. Good correlation between microscopic detection of sporangia and PCR results in both positive and negative soil samples was dually demonstrated for both the SSU and ITS-1 assays.
Subject(s)
Chytridiomycota/isolation & purification , Genetic Variation , Plant Diseases/microbiology , Real-Time Polymerase Chain Reaction/methods , Soil Microbiology , Solanum tuberosum/microbiology , Base Sequence , Chytridiomycota/classification , Chytridiomycota/genetics , DNA Primers/genetics , DNA, Fungal/chemistry , DNA, Fungal/genetics , DNA, Ribosomal/chemistry , DNA, Ribosomal/genetics , DNA, Ribosomal Spacer/chemistry , DNA, Ribosomal Spacer/genetics , Molecular Sequence Data , Phylogeny , Sensitivity and Specificity , Sequence Analysis, DNAABSTRACT
The genus Mikrocytos is traditionally known for Mikrocytos mackini, the microcell parasite that typically infects Pacific oysters along the west coast of North America. Multiple factors have conspired to create difficulty for scientific research on Mikrocytos parasites. These include their tiny cell size, infections that are often of light intensity, lack of suitable cell lines and techniques for in vitro culture, and the seasonal nature of infections. The extreme rate of molecular evolution in Mikrocytos stymied new species discovery and confounded attempts to resolve its phylogenetic position for many years. Fortunately, 2 recent landmark studies have paved the way forward for future research by drastically changing our understanding of the evolution and diversity of these parasites. No longer an orphan eukaryotic lineage, the phylogenetic placement of Mikrocytos has been confidently resolved within Rhizaria and as sister taxon to Haplosporidia. The genus has also found a taxonomic home within the newly-discovered order, Mikrocytida - a globally distributed lineage of parasites infecting a wide range of invertebrate hosts. Here we review available scientific information on Mikrocytos parasites including their evolution and diversity, host and geographic ranges, epizootiology, and detection of the regulated pathogen, M. mackini. We also make recommendations towards a consistent taxonomic framework for this genus by minimally suggesting the use of 18S rDNA sequence, host species information, and histopathological presentation in new species descriptions. This is timely given that we are likely embarking on a new era of scientific advancements, including species discovery, in this genus and its relatives.
Subject(s)
Eukaryota/physiology , Ostreidae/parasitology , Animals , Biological Evolution , Eukaryota/genetics , Host-Parasite Interactions , PhylogenyABSTRACT
Mikrocytos mackini is a microcell parasite that usually infects Crassostrea gigas distributed along the Pacific Northwest coast of North America. For many years, M. mackini was the only known species in the genus, but there have been multiple recent findings of genetically divergent forms of Mikrocytos in different hosts and in distantly located geographic locations. This note describes M. boweri sp. nov. found in Olympia oysters Ostrea lurida collected from and native to British Columbia, Canada, primarily using a molecular taxonomic approach.
Subject(s)
Eukaryota/genetics , Eukaryota/physiology , Ostrea/parasitology , Animals , British Columbia , Host-Parasite Interactions , Phylogeny , Species SpecificityABSTRACT
The economic and methodological efficiencies of environmental DNA (eDNA) based survey approaches provide an unprecedented opportunity to assess and monitor aquatic environments. However, instances of inadequate communication from the scientific community about confidence levels, knowledge gaps, reliability, and appropriate parameters of eDNA-based methods have hindered their uptake in environmental monitoring programs and, in some cases, has created misperceptions or doubts in the management community. To help remedy this situation, scientists convened a session at the Second National Marine eDNA Workshop to discuss strategies for improving communications with managers. These include articulating the readiness of different eDNA applications, highlighting the strengths and limitations of eDNA tools for various applications or use cases, communicating uncertainties associated with specified uses transparently, and avoiding the exaggeration of exploratory and preliminary findings. Several key messages regarding implementation, limitations, and relationship to existing methods were prioritized. To be inclusive of the diverse managers, practitioners, and researchers, we and the other workshop participants propose the development of communication workflow plans, using RACI (Responsible, Accountable, Consulted, Informed) charts to clarify the roles of all pertinent individuals and parties and to minimize the chance for miscommunications. We also propose developing decision support tools such as Structured Decision-Making (SDM) to help balance the benefits of eDNA sampling with the inherent uncertainty, and developing an eDNA readiness scale to articulate the technological readiness of eDNA approaches for specific applications. These strategies will increase clarity and consistency regarding our understanding of the utility of eDNA-based methods, improve transparency, foster a common vision for confidently applying eDNA approaches, and enhance their benefit to the monitoring and assessment community.
ABSTRACT
Perkinsus qugwadi, a pathogenic protozoan parasite of Yesso scallops Patinopecten yessoensis, is found only in cultured populations in British Columbia, Canada. This pathogen was first identified in 1988 and caused significant mortalities at some locations during the early 1990s. Prevalence of infection decreased dramatically following 1995, and the disease was last reported in 1997, leading to speculation that the Yesso scallop stocks in Canada had developed resistance to the disease, or that P. qugwadi had disappeared. However, the present study revealed that infection with P. qugwadi and associated mortality is still occurring in scallops from at least one location in British Columbia. One of the PCR tests developed for P. qugwadi detected the parasite in a 105-fold dilution of DNA extracted from a heavily infected sample and detected 52% more positive scallops than histology; however, the assay also cross-reacted with P. honshuensis and P. olseni. The other PCR test was less sensitive and detected 34% more positives, but did not react to any of the other Perkinsus species tested, suggesting that these PCR tests are powerful tools for screening for the presence of P. qugwadi. Phylogenetic analysis of 1796 bp of SSU rRNA gene sequence clearly indicated that P. qugwadi is positioned basally to other Perkinsus species.
Subject(s)
Alveolata/isolation & purification , Pectinidae , Polymerase Chain Reaction/methods , Alveolata/classification , Alveolata/genetics , Animals , Canada , Phylogeny , RNA, Ribosomal/geneticsABSTRACT
Mikrocytos mackini, causative agent of Denman Island disease in Pacific oysters Crassostrea gigas and other oyster species, was found in 2011 in a previously unreported host, the Kumamoto oyster C. sikamea, in Humboldt Bay, California, USA. The detection was also the first reported finding of M. mackini in California. Prevalence was estimated as high as approximately 27% from pooled samples analyzed by PCR. Higher prevalence appeared related to longer residence time in the bay and somewhat colder than typical winter seawater temperatures. No M. mackini was detected in Humboldt Bay juvenile Kumamoto oysters or Pacific oyster brood or seed stock in 2011 or 2012.
Subject(s)
Crassostrea/parasitology , Parasites/classification , Parasites/isolation & purification , Animals , Host-Parasite Interactions , Species Specificity , Time FactorsABSTRACT
Advances in high-throughput sequencing (HTS) are revolutionizing monitoring in marine environments by enabling rapid, accurate and holistic detection of species within complex biological samples. Research institutions worldwide increasingly employ HTS methods for biodiversity assessments. However, variance in laboratory procedures, analytical workflows and bioinformatic pipelines impede the transferability and comparability of results across research groups. An international experiment was conducted to assess the consistency of metabarcoding results derived from identical samples and primer sets using varying laboratory procedures. Homogenized biofouling samples collected from four coastal locations (Australia, Canada, New Zealand and the USA) were distributed to 12 independent laboratories. Participants were asked to follow one of two HTS library preparation workflows. While DNA extraction, primers and bioinformatic analyses were purposefully standardized to allow comparison, many other technical variables were allowed to vary among laboratories (amplification protocols, type of instrument used, etc.). Despite substantial variation observed in raw results, the primary signal in the data was consistent, with the samples grouping strongly by geographical origin for all data sets. Simple post hoc data clean-up by removing low-quality samples gave the best improvement in sample classification for nuclear 18S rRNA gene data, with an overall 92.81% correct group attribution. For mitochondrial COI gene data, the best classification result (95.58%) was achieved after correction for contamination errors. The identified critical methodological factors that introduced the greatest variability (preservation buffer, sample defrosting, template concentration, DNA polymerase, PCR enhancer) should be of great assistance in standardizing future biodiversity studies using metabarcoding.
Subject(s)
DNA Barcoding, Taxonomic , Laboratories , Biodiversity , High-Throughput Nucleotide Sequencing , Humans , RNA, Ribosomal, 18SABSTRACT
Mikrocytos mackini is a microcell parasite of Pacific oysters only known to occur on the Pacific coast of North America. It is the only described species in the genus, although a genetically divergent Mikrocytos sp. organism has been reported once in both the Atlantic Ocean and China. We developed methods for sequencing the internal transcribed spacer (ITS) of rDNA for the purpose of characterizing extant diversity within M. mackini throughout its known geographic range, and surveying for other evidence of Mikrocytos sp. organisms. Our specific aims were to examine relatedness of M. mackini among sites to make inferences about its recent evolutionary history, and to provide baseline data for future development of a species-specific molecular detection method. We found a total lack of genetic variation within M. mackini across the complete ITS1-5.8S-ITS2 array in over 70 samples collected throughout its range. We hypothesize that this could be a result of a founder effect if the parasite had been introduced into its known range alongside its host, which was imported from Asia beginning around 1914 to about 1961. We detected a single divergent sequence at a short stretch of 18S that was identical to the Mikrocytos sp. detected elsewhere, which adds to the recent and growing body of evidence that Mikrocytos is much more broadly distributed than the limited range of M. mackini suggests. A 1903 bp section of rDNA from Mikrocytos sp. was generated that contained regions of high divergence from M. mackini (in ITS1 and ITS2) that could be exploited for molecular diagnostics.
Subject(s)
DNA, Ribosomal Spacer/genetics , Ostreidae/parasitology , Parasites/genetics , Animals , Base Sequence , Genetic Variation , Host-Parasite Interactions , Polymerase Chain ReactionABSTRACT
Incidental detection of species of concern (e.g., invasive species, pathogens, threatened and endangered species) during biodiversity assessments based on high-throughput DNA sequencing holds significant risks in the absence of rigorous, fit-for-purpose data quality and reporting standards. Molecular biodiversity data are predominantly collected for ecological studies and thus are generated to common quality assurance standards. However, the detection of certain species of concern in these data would likely elicit interest from end users working in biosecurity or other surveillance contexts (e.g., pathogen detection in health-related fields), for which more stringent quality control standards are essential to ensure that data are suitable for informing decision-making and can withstand legal or political challenges. We suggest here that data quality and reporting criteria are urgently needed to enable clear identification of those studies that may be appropriately applied to surveillance contexts. In the interim, more pointed disclaimers on uncertainties associated with the detection and identification of species of concern may be warranted in published studies. This is not only to ensure the utility of molecular biodiversity data for consumers, but also to protect data generators from uncritical and potentially ill-advised application of their science in decision-making.
ABSTRACT
Metabarcoding combines DNA barcoding with high-throughput sequencing, often using one genetic marker to understand complex and taxonomically diverse samples. However, species-level identification depends heavily on the choice of marker and the selected primer pair, often with a trade-off between successful species amplification and taxonomic resolution. We present a versatile metabarcoding protocol for biomonitoring that involves the use of two barcode markers (COI and 18S) and four primer pairs in a single high-throughput sequencing run, via sample multiplexing. We validate the protocol using a series of 24 mock zooplanktonic communities incorporating various levels of genetic variation. With the use of a single marker and single primer pair, the highest species recovery was 77%. With all three COI fragments, we detected 62%-83% of species across the mock communities, while the use of the 18S fragment alone resulted in the detection of 73%-75% of species. The species detection level was significantly improved to 89%-93% when both markers were used. Furthermore, multiplexing did not have a negative impact on the proportion of reads assigned to each species and the total number of species detected was similar to when markers were sequenced alone. Overall, our metabarcoding approach utilizing two barcode markers and multiple primer pairs per barcode improved species detection rates over a single marker/primer pair by 14% to 35%, making it an attractive and relatively cost-effective method for biomonitoring natural zooplankton communities. We strongly recommend combining evolutionary independent markers and, when necessary, multiple primer pairs per marker to increase species detection (i.e., reduce false negatives) in metabarcoding studies.
ABSTRACT
Molecular genetic approaches are playing an increasing role in conservation science by identifying biodiversity that may not be evident by morphology-based taxonomy and systematics. So-called cryptic species are particularly prevalent in freshwater environments, where isolation of dispersal-limited species, such as crayfishes, within dendritic river networks often gives rise to high intra- and inter-specific genetic divergence. We apply here a multi-gene molecular approach to investigate relationships among extant species of the crayfish genus Pacifastacus, representing the first comprehensive phylogenetic study of this taxonomic group. Importantly, Pacifastacus includes both the widely invasive signal crayfish Pacifastacus leniusculus, as well as several species of conservation concern like the Shasta crayfish Pacifastacus fortis. Our analysis used 83 individuals sampled across the four extant Pacifastacus species (omitting the extinct Pacifastacus nigrescens), representing the known taxonomic diversity and geographic distributions within this genus as comprehensively as possible. We reconstructed phylogenetic trees from mitochondrial (16S, COI) and nuclear genes (GAPDH), both separately and using a combined or concatenated dataset, and performed several species delimitation analyses (PTP, ABGD, GMYC) on the COI phylogeny to propose Primary Species Hypotheses (PSHs) within the genus. All phylogenies recovered the genus Pacifastacus as monophyletic, within which we identified a range of six to 21 PSHs; more abundant PSHs delimitations from GMYC and ABGD were always nested within PSHs delimited by the more conservative PTP method. Pacifastacus leniusculus included the majority of PSHs and was not monophyletic relative to the other Pacifastacus species considered. Several of these highly distinct P. leniusculus PSHs likely require urgent conservation attention. Our results identify research needs and conservation priorities for Pacifastacus crayfishes in western North America, and may inform better understanding and management of P. leniusculus in regions where it is invasive, such as Europe and Japan.
ABSTRACT
Successful invasion by nonindigenous species is often attributed to high propagule pressure, yet some foreign species become widespread despite showing reduced genetic variation due to founder effects. The signal crayfish (Pacifastacus leniusculus) is one such example, where rapid spread across Japan in recent decades is believed to be the result of only three founding populations. To infer the history and explore the success of this remarkable crayfish invasion, we combined detailed phylogeographical and morphological analyses conducted in both the introduced and native ranges. We sequenced 16S mitochondrial DNA of signal crayfish from across the introduced range in Japan (537 samples, 20 sites) and the native range in western North America (700 samples, 50 sites). Because chela size is often related to aggressive behavior in crayfish, and hence, their invasion success, we also measured chela size of a subset of specimens in both introduced and native ranges. Genetic diversity of introduced signal crayfish populations was as high as that of the dominant phylogeographic group in the native range, suggesting high propagule pressure during invasion. More recently established crayfish populations in Japan that originated through secondary spread from one of the founding populations exhibit reduced genetic diversity relative to older populations, probably as a result of founder effects. However, these newer populations also show larger chela size, consistent with expectations of rapid adaptations or phenotypic responses during the invasion process. Introduced signal crayfish populations in Japan originate from multiple source populations from a wide geographic range in the native range of western North America. A combination of high genetic diversity, especially for older populations in the invasive range, and rapid adaptation to colonization, manifested as larger chela in recent invasions, likely contribute to invasion success of signal crayfish in Japan.