Evaluating the ribosomal internal transcribed spacer (ITS) as a candidate dinoflagellate barcode marker

Internal transcribed spacer as effective molecular marker for the detection of natural hybridization between the bivalves Pinna nobilis and Pinna rudis

The Pinna nobilis, a Mediterranean mollusc, has suffered population declines due to a massive mortality event associated with various factors including the parasite Haplosporidium pinnae. Some populations show resilience, possibly due to local environmental conditions. In this study, a molecular multiplex PCR method was developed using species-specific primers targeting Internal Transcribed Spacer (ITS) regions of P. nobilis and P. rudis, allowing accurate species identification and hybrid detection. Samples from Mediterranean areas were analysed, including putative hybrids and individuals from five other bivalve species. DNA was isolated, ITS regions were amplified and sequenced, and phylogenetic analyses confirmed species differentiation and primer specificity. The multiplex-PCR successfully identified P. nobilis, P. rudis, and their hybrids based on distinct amplicon patterns. This study highlights the value of molecular tools in species conservation, especially for monitoring and managing hybridization, supporting effective biodiversity conservation strategies.

A dinoflagellate bloom caused by multiple species of Kareniaceae in the coastal waters of Fujian in June 2022 and its adverse impacts on Brachionus plicatilis and Artemia salina

Recently, dinoflagellate blooms have frequently occurred in the coastal waters of Fujian, East China Sea. In June 2022, a fish-killing bloom of Kareniaceae species occurred in this region. In this study, four species of Kareniaceae, namely, Karenia longicanalis, K. papilionacea, Karlodinium veneficum, and Karl. digitatum were identified from this bloom event based on the results of single-cell PCR and clone libraries, and intraspecies genetic diversity was found in the Karl. veneficum population. The results of acute toxicity assays of the bloom water to two zooplankton species (Brachionus plicatilis and Artemia salina) demonstrated this bloom event strongly inhibited their swimming capacities and survival. The results of this study suggested that the bloom events caused by multiple species of Kareniaceae in the Fujian coastal waters had adverse impacts on the local fishery resources and zooplankton community.

Development of a rapid detection method for Karenia mikimotoi by using CRISPR-Cas12a

Harmful algal blooms (HABs), mainly formed by dinoflagellates, have detrimental effects on marine ecosystems and public health. Therefore, detecting HABs is crucial for early warning and prevention of HABs as well as the mitigation of their adverse effects. Although various methods, such as light microscopy, electron microscopy, real-time PCR, and microarrays, have already been established for the detection of HABs, they are still cumbersome to be exploited in the field. Therefore, rapid nucleic detection methods such as recombinase polymerase amplification (RPA) and loop-mediated isothermal amplification (LAMP)-lateral flow dipstick (LFD) have been developed for monitoring bloom-forming algae. However, the CRISPR/Cas-based detection of HABs has yet to be applied to this field. In this study, we developed a method for detecting Karenia mikimotoi (K. mikimotoi), a typical ichthyotoxic dinoflagellate responsible for global blooms. Our method utilized Cas12a from Lachnospiraceae bacterium ND2006 (LbCas12a) to target and cleave the internal transcribed spacer (ITS) of K. mikimotoi, guided by RNA. We leveraged the target-activated non-specific single-stranded deoxyribonuclease cleavage activity of LbCas12a to generate signals that can be detected using fluorescence-read machines or LFDs. By combining RPA and LbCas12a with reporters, we significantly enhanced the sensitivity, enabling the detection of ITS-harboring plasmids at concentrations as low as 9.8 aM and genomic DNA of K. mikimotoi at levels as low as 3.6 × 10-5 ng/μl. Moreover, we simplified the genomic DNA extraction method using cellulose filter paper (CFP) by directly eluting the DNA into RPA reactions, reducing the extraction time to < 30 s. The entire process, from genomic DNA extraction to result reporting, takes less than an hour, enabling the identification of nearly a single cell. In conclusion, our method provided an easy, specific, and sensitive approach for detecting K. mikimotoi, offering the potential for efficient monitoring and management of K. mikimotoi blooms.

Genetic Markers for Metabarcoding of Freshwater Microalgae: Review

The metabarcoding methods for studying the diversity of freshwater microalgae and routine biomonitoring are actively used in modern research. A lot of experience has been accumulated already, and many methodological questions have been solved (such as the influence of the methods and time of sample conservation, DNA extraction and bioinformatical processing). The reproducibility of the method has been tested and confirmed. However, one of the main problems-choosing a genetic marker for the study-still lacks a clear answer. We analyzed 70 publications and found out that studies on eukaryotic freshwater microalgae use 12 markers (different nuclear regions 18S and ITS and plastids rbcL, 23S and 16S). Each marker has its peculiarities; they amplify differently and have various levels of efficiency (variability) in different groups of algae. The V4 and V9 18S and rbcL regions are used most often. We concentrated especially on the studies that compare the results of using different markers and microscopy. We summarize the data on the primers for each region and on how the choice of a marker affects the taxonomic composition of a community.

Oleaginous Heterotrophic Dinoflagellates—Crypthecodiniaceae

The heterotrophic Crypthecodinium cohnii is a major model for dinoflagellate cell biology, and a major industrial producer of docosahexaenoic acid, a key nutraceutical and added pharmaceutical compound. Despite these factors, the family Crypthecodiniaceae is not fully described, which is partly attributable to their degenerative thecal plates, as well as the lack of ribotype-referred morphological description in many taxons. We report here significant genetic distances and phylogenetic cladding that support inter-specific variations within the Crypthecodiniaceae. We describe Crypthecodinium croucheri sp. nov. Kwok, Law and Wong, that have different genome sizes, ribotypes, and amplification fragment length polymorphism profiles when compared to the C. cohnii. The interspecific ribotypes were supported by distinctive truncation-insertion at the ITS regions that were conserved at intraspecific level. The long genetic distances between Crypthecodiniaceae and other dinoflagellate orders support the separation of the group, which includes related taxons with high oil content and degenerative thecal plates, to be ratified to the order level. The current study provides the basis for future specific demarcation-differentiation, which is an important facet in food safety, biosecurity, sustainable agriculture feeds, and biotechnology licensing of new oleaginous models.

DNA barcoding, an effective tool for species identification: a review

DNA barcoding is a powerful taxonomic tool to identify and discover species. DNA barcoding utilizes one or more standardized short DNA regions for taxon identification. With the emergence of new sequencing techniques, such as Next-generation sequencing (NGS), ONT MinION nanopore sequencing, and Pac Bio sequencing, DNA barcoding has become more accurate, fast, and reliable. Rapid species identification by DNA barcodes has been used in a variety of fields, including forensic science, control of the food supply chain, and disease understanding. The Consortium for Barcode of Life (CBOL) presents various working groups to identify the universal barcode gene, such as COI in metazoans; rbcL, matK, and ITS in plants; ITS in fungi; 16S rRNA gene in bacteria and archaea, and creating a reference DNA barcode library. In this article, an attempt has been made to analyze the various proposed DNA barcode for different organisms, strengths & limitations, recent advancements in DNA barcoding, and methods to speed up the DNA barcode reference library construction. This study concludes that constructing a reference library with high species coverage would be a major step toward identifying species by DNA barcodes. This can be achieved in a short period of time by using advanced sequencing and data analysis methods.

Isolation, identification and toxicity of three strains of Heterocapsa (Dinophyceae) in a harmful event in Fujian, China

Heterocapsa is a genus of dinoflagellates including species that can form harmful algae blooms (HABs) and cause deleterious ecological effects. To date two species ( H circularisquama and H bohaiensis) are known to be toxic. Here we isolated three strains of Heterocapsa from a multi-species dinoflagellate bloom in the aquaculture area of Fujian, China in June, 2019 that caused mass mortality of farmed abalone. Morphological analysis using light, transmission electron and scanning electron microscopy along with phylogenetic analyses with small (SSU) and large (LSU) subunit rRNA and internal transcribed spacer (ITS) gene sequences showed that these strains were H. cf. niei H. horiguchii, and H. cf. pygmaea, respectively. Furthermore, rabbit erythrocyte assay revealed hemolytic activity in all three strains in cell density dependent fashion, and only in the presence of light. In addition, the strains caused significant mortality of Artemia salina, and the toxicity was also cell density dependent. The Heterocapsa cultures and toxicity information obtained in this study expanded our knowledge of toxic species of Heterocapsa, and will facilitate further investigating the mechanism of their toxicity and developing monitoring tools for their blooms in the future.

Six marine thecate Heterocapsa (Dinophyceae) from Malaysia, including the description of three novel species and their cytotoxicity potential

Thirty-four strains of Heterocapsa were established from Malaysian waters and their morphologies were examined by light, scanning, and transmission electron microscopy. Three species, H. bohaiensis, H. huensis, and H. rotundata, and three new species, H. borneoensis sp. nov., H. limii sp. nov., and H. iwatakii sp. nov. were described in this study. The three species were differentiated morphologically by unique characteristics of cell size, shape, displacement of the cingulum, shape and position of nucleus, the number and position of pyrenoids, and body scale ultrastructure. The species delimitations were robustly supported by the molecular data. A light-microscopy-based key to species of Heterocapsa is established, with two major groups, i.e., species with a single pyrenoid, and species with multiple pyrenoids. Bioassays were conducted by exposing Artemia nauplii to Heterocapsa densities of 1-5 × 10⁵ cells mL^-1, and treatments exposed to H. borneoensis showed naupliar mortality, while no naupliar death was observed in the treatments exposed to cells of H. bohaiensis, H. huensis, H. limii, and H. iwatakii. Naupliar death was observed during the initial 24 h for both tested H. borneoensis strains, and mortality rates increased up to 50% after 72-h exposure. This study documented for the first time the diversity and cytotoxic potency of Heterocapsa species from Malaysian waters.

Bloom of Prorocentrum cordatum in Paracas Bay, Peru

During the austral winter of 2017, a bloom of Prorocentrum spp. occurred, reaching a cell density of 2.73 × 106 cells L−1, in Paracas Bay, Peru. In order to identify which, type of species generated this event and determine its toxicity, the values of the environmental parameters (temperature, winds and salinity) that induced the rapid growth of the dinoflagellate in this bloom were identified. A clonal culture was established for taxonomic (SEM), phylogenetic (ITS) and toxicological analysis via LC-MS/MS to determine the presence of tetrodotoxin (TTX) and whether the species represents a food safety hazard. This event coincided with the coastal upwelling process, which generated high concentrations of phytoplankton biomass (>10 mg m−3 chlorophyll-a) and allowed the rapid growth of P. cordatum (IMP-BG 450) in Paracas Bay. However, toxicological analyses of the IMP-BG 450 strain culture did not show the presence of TTX quantifiable through the technique used. Due to the antecedents of the presence of TTX in mollusks from other latitudes during blooms of this species, it is recommended that analyses of this toxin be carried out both in filter-feeding mollusks and in this species during a new bloom.

Using RDNA sequences to define dinoflagellate species

Dinoflagellate species are traditionally defined using morphological characters, but molecular evidence accumulated over the past several decades indicates many morphologically-based descriptions are inaccurate. This recognition led to an increasing reliance on DNA sequence data, particularly rDNA gene segments, in defining species. The validity of this approach assumes the divergence in rDNA or other selected genes parallels speciation events. Another concern is whether single gene rDNA phylogenies by themselves are adequate for delineating species or if multigene phylogenies are required instead. Currently, few studies have directly assessed the relative utility of multigene versus rDNA-based phylogenies for distinguishing species. To address this, the current study examined D1-D3 and ITS/5.8S rDNA gene regions, a multi-gene phylogeny, and morphological characters in Gambierdiscus and other related dinoflagellate genera to determine if they produce congruent phylogenies and identify the same species. Data for the analyses were obtained from previous sequencing efforts and publicly available dinoflagellate transcriptomic libraries as well from the additional nine well-characterized Gambierdiscus species transcriptomic libraries generated in this study. The D1-D3 and ITS/5.8S phylogenies successfully identified the described Gambierdiscus and Alexandrium species. Additionally, the data showed that the D1-D3 and multigene phylogenies were equally capable of identifying the same species. The multigene phylogenies, however, showed different relationships among species and are likely to prove more accurate at determining phylogenetic relationships above the species level. These data indicated that D1-D3 and ITS/5.8S rDNA region phylogenies are generally successful for identifying species of Gambierdiscus, and likely those of other dinoflagellates. To assess how broadly general this finding is likely to be, rDNA molecular phylogenies from over 473 manuscripts representing 232 genera and 863 described species of dinoflagellates were reviewed. Results showed the D1-D3 rDNA and ITS phylogenies in combination are capable of identifying 97% of dinoflagellate species including all the species belonging to the genera Alexandrium, Ostreopsis and Gambierdiscus, although it should be noted that multi-gene phylogenies are preferred for inferring relationships among these species. A protocol is presented for determining when D1-D3, confirmed by ITS/5.8S rDNA sequence data, would take precedence over morphological features when describing new dinoflagellate species. This protocol addresses situations such as: a) when a new species is both morphologically and molecularly distinct from other known species; b) when a new species and closely related species are morphologically indistinguishable, but genetically distinct; and c) how to handle potentially cryptic species and cases where morphotypes are clearly distinct but have the same rDNA sequence. The protocol also addresses other molecular, morphological, and genetic approaches required to resolve species boundaries in the small minority of species where the D1-D3/ITS region phylogenies fail.

q-PCR-based assay for the toxic dinoflagellate Karenia selliformis monitoring along the Tunisian coasts

Karenia selliformis is a marine dinoflagellate responsible for fish-kill events. Its presence has been reported along the Tunisian coasts (south-eastern Mediterranean Sea) since the 1990s. In the present study, a quantitative-PCR assay, based on the internal transcribed spacer (ITS) molecular marker, was developed to detect and quantify K. selliformis in environmental bivalve mollusk samples and in seawater samples. The assay was optimized, and its specificity was confirmed using cross-reactivity experiments against microalgal species commonly found on the Tunisian coasts and/or closely related to K. selliformis. Calibration curves were performed by tenfold dilutions of plasmid DNA harboring target sequence and genomic DNA, attaining a limit of detection of around 5 copies of target DNA per reaction, far below one K. selliformis cell per reaction. The field application of the developed assay showed a powerful detection capability. Thus, the designed assay could contribute to the deployment of in-field diagnostic tools for K. selliformis blooms monitoring.

Resolving cryptic species complexes in marine protists: phylogenetic haplotype networks meet global DNA metabarcoding datasets

Marine protists have traditionally been assumed to be lowly diverse and cosmopolitan. Yet, several recent studies have shown that many protist species actually consist of cryptic complexes of species whose members are often restricted to particular biogeographic regions. Nonetheless, detection of cryptic species is usually hampered by sampling coverage and application of methods (e.g. phylogenetic trees) that are not well suited to identify relatively recent divergence and ongoing gene flow. In this paper, we show how these issues can be overcome by inferring phylogenetic haplotype networks from global metabarcoding datasets. We use the Chaetoceros curvisetus (Bacillariophyta) species complex as study case. Using two complementary metabarcoding datasets (Ocean Sampling Day and Tara Oceans), we equally resolve the cryptic complex in terms of number of inferred species. We detect new hypothetical species in both datasets. Gene flow between most of species is absent, but no barcoding gap exists. Some species have restricted distribution patterns whereas others are widely distributed. Closely related taxa occupy contrasting biogeographic regions, suggesting that geographic and ecological differentiation drive speciation. In conclusion, we show the potential of the analysis of metabarcoding data with evolutionary approaches for systematic and phylogeographic studies of marine protists.

The parasitic dinoflagellate Hematodinium infects multiple crustaceans in the polyculture systems of Shandong Province, China

The parasitic dinoflagellates of the Hematodinium genus have impacted wild and cultured stocks of commercial crustaceans worldwide. In the past decade, outbreaks of Hematodinium epizootics resulted in substantial mortalities in cultured Chinese swimming crabs Portunus trituberculatus in the polyculture ponds located in Shandong Peninsula, whereas the source and transmission of the parasite in the polyculture pond system remains to be determined. During April to December of 2018, 2034 crabs and 108 shrimps were collected from the polyculture pond systems in the highly endemic area of Hematodinium diseases in Qingdao, Shandong Province. Among those, 188 individuals of the 6 crab species were infected by the parasite, including 4 novel host species (Uca arcuate, Hemigrapsus penicillatus, Helice wuana and Macrophthalmus japonicas). No infection was identified in Penaeus monodon. Further phylogenetic analyses indicated that the Hematodinium isolate infecting the six crab hosts, together with other isolates reported from China, composed the genotype II of Hematodinium perezi. The parasite was more infectious to cultured Portunus trituberculatus and the dominant wild crab Helice tientsinensis dwelling in the waterways connecting to the polyculture ponds, even though it was found to be a host generalist pathogen. The prevalence of Hematodinium perezi infection in Helice tientsinensis was higher than that of other wild crabs and showed significant positive correlation with that of the cultured Portunus trituberculatus. The results indicated that the wild crabs, particularly Helice tientsinensis, were the important alternate hosts closely involved in transmission and spreading of the Hematodinium disease in the polyculture pond systems.

The Application of Nanopore Sequencing Technology to the Study of Dinoflagellates: A Proof of Concept Study for Rapid Sequence-Based Discrimination of Potentially Harmful Algae

Harmful algal blooms (HABs) are a naturally occurring global phenomena that have the potential to impact fisheries, leisure and ecosystems, as well as posing a significant hazard to animal and human health. There is significant interest in the development and application of methodologies to study all aspects of the causative organisms and toxins associated with these events. This paper reports the first application of nanopore sequencing technology for the detection of eukaryotic harmful algal bloom organisms. The MinION sequencing platform from Oxford Nanopore technologies provides long read sequencing capabilities in a compact, low cost, and portable format. In this study we used the MinION to sequence long-range PCR amplicons from multiple dinoflagellate species with a focus on the genus Alexandrium. Primers applicable to a wide range of dinoflagellates were selected, meaning that although the study was primarily focused on Alexandrium the applicability to three additional genera of toxic algae, namely; Gonyaulax, Prorocentrum, and Lingulodinium was also demonstrated. The amplicon generated here spanned approximately 3 kb of the rDNA cassette, including most of the 18S, the complete ITS1, 5.8S, ITS2 and regions D1 and D2 of the 28S. The inclusion of barcode genes as well as highly conserved regions resulted in identification of organisms to the species level. The analysis of reference cultures resulted in over 99% of all sequences being attributed to the correct species with an average identity above 95% from a reference list of over 200 species (see Supplementary Material 1). The use of mock community analysis within environmental samples highlighted that complex matrices did not prevent the ability to distinguish between phylogenetically similar species. Successful identification of causative organisms in environmental samples during natural toxic events further highlighted the potential of the assay. This study proves the suitability of nanopore sequencing technology for taxonomic identification of harmful algal bloom organisms and acquisition of data relevant to the World Health Organisations "one health" approach to marine monitoring.

Establishment and application of hyperbranched rolling circle amplification coupled with lateral flow dipstick for the sensitive detection of Karenia mikimotoi

The dinoflagellate Karenia mikimotoi is a noxious and harmful algal bloom (HAB)-forming microalga. Establishing a rapid, accurate, and sensitive method of detecting this harmful alga is necessary to provide warnings of imminent HABs through field monitoring. Here, an isothermal amplification technique combined with a rapid analytical method for nucleic acid-based amplified products, i.e., hyperbranched rolling circle amplification (HRCA) coupled with lateral flow dipstick (LFD), hereafter denoted as HRCA-LFD, was established to detect K. mikimotoi. The HRCA-LFD assay relied on a padlock probe (PLP) targeting DNA template and an LFD probe targeting PLP. The sequenced internal transcribed spacer of K. mikimotoi through molecular cloning was used as the target of PLP. The optimized HRCA conditions was determined to be as follows: PLP concentration, 20 pM; ligation temperature, 65 °C; ligation time, 10 min; amplification temperature, 61 °C; and amplification time, 30 min. The developed HRCA-LFD assay was specific for K. mikimotoi, displaying no cross-reactivity with other common microalgae. Sensitivity-comparison tests indicated that HRCA-LFD assay was 100-fold more sensitive than PCR, with a detection limit of 0.1 cell mL^-1 when used to analyze spiked field samples. The analysis with field samples also indicated that HRCA-LFD assay was suitable for samples with a target cell density range of 1-1000 cells mL^-1. All of these results suggested that HRCA-LFD assay is an alternative method for the sensitive and reliable detection of K. mikimotoi from marine water samples.

Distribution of Karlodinium veneficum in the coastal region of Xiangshan Bay in the East China Sea, as detected by a real-time quantitative PCR assay of ribosomal ITS sequence

Athecate dinoflagellate Karlodinium veneficum is a universal toxic species possessing karlotoxins recognized especially as ichthyotoxic as well as cytotoxic and hemolytic. Blooms of K. veneficum, both single-species or accompanied with other species, occurred more frequently worldwide in recent years, including the coastal region of China. Normally, K. veneficum present in relatively low abundance in phytoplankton communities in estuary regions. Being small and difficult to identify with light microscopy, it has been ignored for a long time till its blooming and toxins being confirmed. How it presents in background level and what is its relationship with critical geological and hydrological environment factors are basically not clear. In this study, the paper reports the application of a real-time quantitative PCR (qPCR) method to investigate the abundance and distribution of K. veneficum in the coastal waters of Xiangshan Bay in the East China Sea (ECS), a typical bay area of harmful algae blooms and heavily affected by anthropogenic activities. The real-time qPCR assay came out being an efficient method at detecting even low cell densities of K. veneficum of different genotypes. A total of 38 field samples of surface (0.5 m) and bottom water (9-100 m in depth) were analyzed and 12 samples were found positive for K. veneficum. At least 3 genotypes of K. veneficum present in this region. Temperatures in sites of K. veneficum positive ranged from 21.7 to 23.4 °C, and salinity levels were between 21.1 and 26.3. The K. veneficum distributed quite extensively in the waters of Xiangshan Bay, cell abundance varied from a low of 4 cells/L to a maximum of 170 cells/L. Most of the samples containing K. veneficum were collected from bottom water in different sites. At three of the 19 sampling sites, K. veneficum was detected in both surface and bottom water samples. Especially at sampling site near Beilun port, where the water is typically muddy with low transparency, relative high cell numbers of K. veneficum were found in both surface and bottom waters. Mixotrophy and vertical migration of K. veneficum could be important eco-physiological factors to consider in terms of understanding these distribution characteristics. The ideal conditions for K. veneficum growth and aggregation in this area still needs further study.

Systematic Revision of Symbiodiniaceae Highlights the Antiquity and Diversity of Coral Endosymbionts

The advent of molecular data has transformed the science of organizing and studying life on Earth. Genetics-based evidence provides fundamental insights into the diversity, ecology, and origins of many biological systems, including the mutualisms between metazoan hosts and their micro-algal partners. A well-known example is the dinoflagellate endosymbionts ("zooxanthellae") that power the growth of stony corals and coral reef ecosystems. Once assumed to encompass a single panmictic species, genetic evidence has revealed a divergent and rich diversity within the zooxanthella genus Symbiodinium. Despite decades of reporting on the significance of this diversity, the formal systematics of these eukaryotic microbes have not kept pace, and a major revision is long overdue. With the consideration of molecular, morphological, physiological, and ecological data, we propose that evolutionarily divergent Symbiodinium "clades" are equivalent to genera in the family Symbiodiniaceae, and we provide formal descriptions for seven of them. Additionally, we recalibrate the molecular clock for the group and amend the date for the earliest diversification of this family to the middle of the Mesozoic Era (∼160 mya). This timing corresponds with the adaptive radiation of analogs to modern shallow-water stony corals during the Jurassic Period and connects the rise of these symbiotic dinoflagellates with the emergence and evolutionary success of reef-building corals. This improved framework acknowledges the Symbiodiniaceae's long evolutionary history while filling a pronounced taxonomic gap. Its adoption will facilitate scientific dialog and future research on the physiology, ecology, and evolution of these important micro-algae.

Delineating closely related dinoflagellate lineages using phylotranscriptomics

Recently radiated dinoflagellates Apocalathium aciculiferum (collected in Lake Erken, Sweden), Apocalathium malmogiense (Baltic Sea) and Apocalathium aff. malmogiense (Highway Lake, Antarctica) represent a lineage with an unresolved phylogeny. We determined their phylogenetic relationships using phylotranscriptomics based on 792 amino acid sequences. Our results showed that A. aciculiferum diverged from the other two closely related lineages, consistent with their different morphologies in cell size, relative cell length and presence of spines. We hypothesized that A. aff. malmogiense and A. malmogiense, which inhabit different hemispheres, are evolutionarily more closely related because they diverged from a marine common ancestor, adapting to a wide salinity range, while A. aciculiferum colonized a freshwater habitat, by acquiring adaptations to this environment, in particular, salinity intolerance. We show that phylotranscriptomics can resolve the phylogeny of recently diverged protists. This has broad relevance, given that many phytoplankton species are morphologically very similar, and single genes sometimes lack the information to determine species' relationships.

Characterising planktonic dinoflagellate diversity in Singapore using DNA metabarcoding

Dinoflagellates are traditionally identified morphologically using microscopy, which is a time-consuming and labour-intensive process. Hence, we explored DNA metabarcoding using high-throughput sequencing as a more efficient way to study planktonic dinoflagellate diversity in Singapore’s waters. From 29 minimally pre-sorted water samples collected at four locations in western Singapore, DNA was extracted, amplified and sequenced for a 313-bp fragment of the V4–V5 region in the 18S ribosomal RNA gene. Two sequencing runs generated 2,847,170 assembled paired-end reads, corresponding to 573,176 unique sequences. Sequences were clustered at 97% similarity and analysed with stringent thresholds (≥150 bp, ≥20 reads, ≥95% match to dinoflagellates), recovering 28 dinoflagellate taxa. Dinoflagellate diversity captured includes parasitic and symbiotic groups which are difficult to identify morphologically. Richness is similar between the inner and outer West Johor Strait, but variations in community structure are apparent, likely driven by environmental differences. None of the taxa detected in a recent phytoplankton bloom along the West Johor Strait have been recovered in our samples, suggesting that background communities are distinct from bloom communities. The voluminous data obtained in this study contribute baseline information for Singapore’s phytoplankton communities and prompt future research and monitoring to adopt the approach established here.

Identification, growth and toxicity assessment of Coolia Meunier (Dinophyceae) from Nova Scotia, Canada

Benthic dinoflagellates of the toxigenic genus Coolia Meunier (Dinophyceae) are known to have a global distribution in both tropical and temperate waters. The type species, C. monotis, has been reported from the Mediterranean Sea, the NE Atlantic and from Rhode Island, USA in the NW Atlantic, whereas other species in the genus have been reported from tropical locations. Coolia cells were observed in algal drift samples collected at seven sites in Nova Scotia, Canada. Clonal isolates were established from four of these locations and identified with light and scanning electron microscopy, then confirmed with genetic sequencing to be C. monotis. This is the first record of this species in Nova Scotia. The isolates were established and incubated at 18 °C under a 14:10 L:D photoperiod with an approximate photon flux density of 50-60 μmol m^-2 s^-1. Growth experiments using an isolate from Johnston Harbour (CMJH) were carried out at temperatures ranging from 5 to 30 °C under the same photoperiod with an approximate photon flux density of 45-50 μmol m^-2 s^-1. Cells tolerated temperatures from 5 to 25 °C with optimum growth and mucilage aggregate production between 15 and 20 °C. Methanol extracts of this isolate examined by Liquid Chromatography-Mass Spectrometry (LC-MS) did not show the presence of the previously reported cooliatoxin. Toxic effects were assayed using two zebrafish bioassays, the Fish Embryo Toxicity (FET) assay and the General Behaviour and Toxicity (GBT) assay. The results of this study demonstrate a lack of toxicity in C. monotis from Nova Scotia, as has been reported for other genetically-confirmed isolates of this species. Conditions in which cell growth that could potentially degrade water quality and provide substrate and dispersal mechanisms for other harmful microorganisms via mucilage production are indicated.

Comparative study of biofilm formation on biocidal antifouling and fouling-release coatings using next-generation DNA sequencing

The bacterial and eukaryotic communities forming biofilms on six different antifouling coatings, three biocidal and three fouling-release, on boards statically submerged in a marine environment were studied using next-generation sequencing. Sequenced amplicons of bacterial 16S ribosomal DNA and eukaryotic ribosomal DNA internal transcribed spacer were assigned taxonomy by comparison to reference databases and relative abundances were calculated. Differences in species composition, bacterial and eukaryotic, and relative abundance were observed between the biofilms on the various coatings; the main difference was between coating type, biocidal compared to fouling-release. Species composition and relative abundance also changed through time. Thus, it was possible to group replicate samples by coating and time point, indicating that there are fundamental and reproducible differences in biofilms assemblages. The routine use of next-generation sequencing to assess biofilm formation will allow evaluation of the efficacy of various commercial coatings and the identification of targets for novel formulations.

The diversity of ice algal communities on the Greenland Ice Sheet as revealed by oligotyping

The Arctic is being disproportionally affected by climate change compared with other geographic locations, and is currently experiencing unprecedented melt rates. The Greenland Ice Sheet (GrIS) can be regarded as the largest supraglacial ecosystem on Earth, and ice algae are the dominant primary producers on bare ice surfaces throughout the course of a melt season. Ice-algal-derived pigments cause a darkening of the ice surface, which in turn decreases albedo and increases melt rates. The important role of ice algae in changing melt rates has only recently been recognized, and we currently know little about their community compositions and functions. Here, we present the first analysis of ice algal communities across a 100 km transect on the GrIS by high-throughput sequencing and subsequent oligotyping of the most abundant taxa. Our data reveal an extremely low algal diversity with Ancylonema nordenskiöldii and a Mesotaenium species being by far the dominant taxa at all sites. We employed an oligotyping approach and revealed a hidden diversity not detectable by conventional clustering of operational taxonomic units and taxonomic classification. Oligotypes of the dominant taxa exhibit a site-specific distribution, which may be linked to differences in temperatures and subsequently the extent of the melting. Our results help to better understand the distribution patterns of ice algal communities that play a crucial role in the GrIS ecosystem.

An evaluation of the genus Amphidinium (Dinophyceae) combining evidence from morphology, phylogenetics, and toxin production, with the introduction of six novel species

The genus Amphidinium is an important group of athecated dinoflagellates because of its high abundance in marine habitats, its member's ability to live in a variety of environmental conditions and ability to produce toxins. Furthermore, the genus is of particular interest in the biotechnology field for its potential in the pharmaceutical arena. Taxonomically the there is a history of complication and confusion over the proper identities and placements of Amphidinium species due to high genetic variability coupled with high morphological conservation. Thirteen years has passed since the most recent review of the group, and while many issues were resolved, some remain. The present study used microscopy, phylogenetics of the 28S region of rDNA, secondary structure of the ITS2 region of rDNA, compensatory base change data, and cytotoxicity data from Amphidinium strains collected world-wide to elucidate remaining confusion. This holistic approach using multiple lines of evidence resulted in a more comprehensive understanding of the morphological, ecological, and genetic characteristics that are attributed to organisms belonging to Amphidinium, including six novel species: A. fijiensis, A. magnum, A. paucianulatum, A. pseudomassartii, A. theodori, and A. tomasii.

Host specificity of Symbiodinium variants revealed by an ITS2 metahaplotype approach

Analysis of the widely used ITS region is confounded by the presence of intragenomic variants (IGVs). In Symbiodinium, the algal symbionts of reef building corals, deep-sequencing analyses are used to characterise communities within corals, yet these analyses largely overlook IGVs. Here we consider that distinct ITS2 sequences could represent IGVs rather than distinct symbiont types and argue that symbionts can be distinguished by their proportional composition of IGVs, described as their ITS2 metahaplotype. Using our metahaplotype approach on Minimum Entropy Decomposition (MED) analysis of ITS2 sequences from the corals Acropora downingi, Cyphastrea microphthalma and Playgyra daedalea, we show the dominance of a single species-specific Symbiodinium C3 variant within each coral species. We confirm the presence of these species-specific symbionts using the psbA non-coding region. Our findings highlight the importance of accounting for IGVs in ITS2 analyses and demonstrate their capacity to resolve biological patterns that would otherwise be overlooked.

DNA barcoding in autotrophic euglenids: evaluation of COI and 18s rDNA

Autotrophic euglenids (Euglenophyceae) are a common and abundant group of microbial eukaryotes in freshwater habitats. They have a limited number of features, which can be observed using light microscopy, thus species identification is often problematic. Establishing a barcode for this group is therefore an important step toward the molecular identification of autotrophic euglenids. Based on the literature, we selected verified species and used a plethora of available methods to validate two molecular markers: COI and 18S rDNA (the whole sequence and three fragments separately) as potential DNA barcodes. Analyses of the COI gene were performed based on the data set of 43 sequences (42 obtained in this study) representing 24 species and the COI gene was discarded as a DNA barcode mainly due to a lack of universal primer sites. For 18S rDNA analyses we used a data set containing 263 sequences belonging to 86 taxonomically verified species. We demonstrated that the whole 18S rDNA is too long to be a useful marker, but from the three shorter analyzed variable regions we recommend variable regions V2V3 and V4 of 18S rDNA as autotrophic euglenid barcodes due to their high efficiency (above 95% and 90%, respectively).

Daphnia parasite dynamics across multiple Caullerya epidemics indicate selection against common parasite genotypes

Studies of parasite population dynamics in natural systems are crucial for our understanding of host-parasite coevolutionary processes. Some field studies have reported that host genotype frequencies in natural populations change over time according to parasite-driven negative frequency-dependent selection. However, the temporal patterns of parasite genotypes have rarely been investigated. Moreover, parasite-driven negative frequency-dependent selection is contingent on the existence of genetic specificity between hosts and parasites. In the present study, the population dynamics and host-genotype specificity of the ichthyosporean Caullerya mesnili, a common endoparasite of Daphnia water fleas, were analysed based on the observed sequence variation in the first internal transcribed spacer (ITS1) of the ribosomal DNA. The Daphnia population of lake Greifensee (Switzerland) was sampled and subjected to parasite screening and host genotyping during C. mesnili epidemics of four consecutive years. The ITS1 of wild-caught C. mesnili-infected Daphnia was sequenced using the 454 pyrosequencing platform. The relative frequencies of C. mesnili ITS1 sequences differed significantly among years: the most abundant C. mesnili ITS1 sequence decreased and rare sequences increased over the course of the study, a pattern consistent with negative frequency-dependent selection. However, only a weak signal of host-genotype specificity between C. mesnili and Daphnia genotypes was detected. Use of cutting edge genomic techniques will allow further investigation of the underlying micro-evolutionary relationships within the Daphnia-C. mesnili system.

Characterization of a Strain of Fukuyoa paulensis (Dinophyceae) from the Western Mediterranean Sea

A single cell of the dinoflagellate genus Fukuyoa was isolated from the island of Formentera (Balearic Islands, west Mediterranean Sea), cultured, and characterized by morphological and molecular methods and toxin analyses. This is the first report of the Gambierdiscus lineage (genera Fukuyoa and Gambierdiscus) from the western Mediterranean Sea, which is cooler than its eastern basin. Molecular analyses revealed that the Mediterranean strain belongs to F. paulensis and that it bears LSU rDNA sequences identical to New Zealand, Australian, and Brazilian strains. It also shared an identical sequence of the more variable ITS-rDNA with the Brazilian strain. Toxin analyses showed the presence of maitotoxin, 54-deoxyCTX1B, and gambieric acid A. This is the first observation of the two latter compounds in a Fukuyoa strain. Therefore, both Gambierdiscus and Fukuyoa should be considered when as contributing to ciguatera fish poisoning. Different strains of Fukuyoa form a complex of morphologically cryptic lineages where F. paulensis stands as the most distantly related nominal species. The comparison of the ITS2 secondary structures revealed the absence of CBCs among strains. The study of the morphological and molecular traits depicted an unresolved taxonomic scenario impacted by the low strains sampling.

Comparative study of the validity of three regions of the 18S-rRNA gene for massively parallel sequencing-based monitoring of the planktonic eukaryote community

The nuclear 18S-rRNA gene has been used as a metabarcoding marker in massively parallel sequencing (MPS)-based environmental surveys for plankton biodiversity research. However, different hypervariable regions have been used in different studies, and their utility has been debated among researchers. In this study, detailed investigations into 18S-rRNA were carried out; we investigated the effective number of sequences deposited in international nucleotide sequence databases (INSDs), the amplification bias, and the amplicon sequence variability among the three variable regions, V1-3, V4-5 and V7-9, using in silico polymerase chain reaction (PCR) amplification based on INSDs. We also examined the primer universality and the taxonomic identification power, using MPS-based environmental surveys in the Sea of Okhotsk, to determine which region is more useful for MPS-based monitoring. The primer universality was not significantly different among the three regions, but the number of sequences deposited in INSDs was markedly larger for the V4-5 region than for the other two regions. The sequence variability was significantly different, with the highest variability in the V1-3 region, followed by the V7-9 region, and the lowest variability in the V4-5 region. The results of the MPS-based environmental surveys showed significantly higher identification power in the V1-3 and V7-9 regions than in the V4-5 region, but no significant difference was detected between the V1-3 and V7-9 regions. We therefore conclude that the V1-3 region will be the most suitable for future MPS-based monitoring of natural eukaryote communities, as the number of sequences deposited in INSDs increases.

Gymnodinium catenatum Graham isolated from the Portuguese coast: Toxin content and genetic characterization

The bloom forming marine dinoflagellate Gymnodinium catenatum Graham has been linked to paralytic shellfish poisoning (PSP) outbreaks in humans. Along the Portuguese coast (NE Atlantic), G. catenatum shows a complex bloom pattern, raising questions about the origin and affinities of each bloom population. In this work, the variability within six cultured strains of G. catenatum isolated from Portuguese coastal waters (S coast, W coast and NW coast), between 1999 and 2011, was investigated. The strains were analyzed for toxin profiling and intra-specific genetic diversity. Regarding the toxin profile, differences recorded between strains could not be assigned to the time of isolation or geographical origin. The parameter that most influenced the toxin profile was the life-cycle stage that originated the culture: vegetative cell versus hypnozygote (resting cyst). At the genetic level, all strains showed similar sequences for the D1-D2 region of the large subunit (LSU) of the nuclear ribosomal DNA (rDNA) and shared complete identity with strains from Spain, Algeria, China and Australia. Conversely, we did not find a total identity match for the ITS-5.8S nuclear rDNA fragment. After sequence analysis, two guanine/adenine (R) single nucleotide polymorphisms (SNP 1 and 2) were detected for all strains, in the ITS1 region. This species has been reported to present very conservative LSU and ITS-5.8S rDNA regions, though with few SNP, including SNP1 of this study, already attributed to strains from certain locations. The SNP here described characterize G. catenatum populations from Portuguese waters and may represent valuable genetic markers for studies on the phylogeography of this species.

Cryptic infection of a broad taxonomic and geographic diversity of tadpoles by Perkinsea protists

The decline of amphibian populations, particularly frogs, is often cited as an example in support of the claim that Earth is undergoing its sixth mass extinction event. Amphibians seem to be particularly sensitive to emerging diseases (e.g., fungal and viral pathogens), yet the diversity and geographic distribution of infectious agents are only starting to be investigated. Recent work has linked a previously undescribed protist with mass-mortality events in the United States, in which infected frog tadpoles have an abnormally enlarged yellowish liver filled with protist cells of a presumed parasite. Phylogenetic analyses revealed that this infectious agent was affiliated with the Perkinsea: a parasitic group within the alveolates exemplified by Perkinsus sp., a "marine" protist responsible for mass-mortality events in commercial shellfish populations. Using small subunit (SSU) ribosomal DNA (rDNA) sequencing, we developed a targeted PCR protocol for preferentially sampling a clade of the Perkinsea. We tested this protocol on freshwater environmental DNA, revealing a wide diversity of Perkinsea lineages in these environments. Then, we used the same protocol to test for Perkinsea-like lineages in livers of 182 tadpoles from multiple families of frogs. We identified a distinct Perkinsea clade, encompassing a low level of SSU rDNA variation different from the lineage previously associated with tadpole mass-mortality events. Members of this clade were present in 38 tadpoles sampled from 14 distinct genera/phylogroups, from five countries across three continents. These data provide, to our knowledge, the first evidence that Perkinsea-like protists infect tadpoles across a wide taxonomic range of frogs in tropical and temperate environments, including oceanic islands.

An Improved DNA Extraction Method for Efficient and Quantitative Recovery of Phytoplankton Diversity in Natural Assemblages

Marine phytoplankton are highly diverse with different species possessing different cell coverings, posing challenges for thoroughly breaking the cells in DNA extraction yet preserving DNA integrity. While quantitative molecular techniques have been increasingly used in phytoplankton research, an effective and simple method broadly applicable to different lineages and natural assemblages is still lacking. In this study, we developed a bead-beating protocol based on our previous experience and tested it against 9 species of phytoplankton representing different lineages and different cell covering rigidities. We found the bead-beating method enhanced the final yield of DNA (highest as 2 folds) in comparison with the non-bead-beating method, while also preserving the DNA integrity. When our method was applied to a field sample collected at a subtropical bay located in Xiamen, China, the resultant ITS clone library revealed a highly diverse assemblage of phytoplankton and other micro-eukaryotes, including Archaea, Amoebozoa, Chlorophyta, Ciliphora, Bacillariophyta, Dinophyta, Fungi, Metazoa, etc. The appearance of thecate dinoflagellates, thin-walled phytoplankton and “naked” unicellular organisms indicates that our method could obtain the intact DNA of organisms with different cell coverings. All the results demonstrate that our method is useful for DNA extraction of phytoplankton and environmental surveys of their diversity and abundance.

Culture-Independent Study of the Late-Stage of a Bloom of the Toxic Dinoflagellate Ostreopsis cf. ovata: Preliminary Findings Suggest Genetic Differences at the Sub-Species Level and Allow ITS2 Structure Characterization

Available genomic data for the toxic, bloom-forming, benthic Ostreopsis spp. are traditionally obtained from isolates rather than from individuals originally present in environmental samples. Samples from the final phase of the first reported Ostreopsis bloom in European North Atlantic waters (Algarve, south coast of Portugal) were studied and characterized, using a culture-independent approach. In the first instance, a microscopy-based analysis revealed the intricate complexity of the samples. Then, we evaluated the adequacy of commonly used molecular tools (i.e., primers and nuclear ribosomal markers) for the study of Ostreopsis diversity in natural samples. A PCR-based methodology previously developed to identify/detect common Ostreopsis species was tested, including one new combination of existing PCR primers. Two sets of environmental rRNA sequences were obtained, one of them (1052 bp) with the newly tested primer set. These latter sequences encompass both the ITS1-5.8S-ITS2 region and the D1/D2 domain of the LSU rRNA gene, leading us to an accurate identification of ITS2. In turn, this allowed us to predict and show for the first time the ITS2 secondary structure of Ostreopsis. With 92 bp in length and a two-helix structure, the ITS2 of this genus revealed to be unique among the dinoflagellates. Both the PCR approach as the phylogenetic analyses allowed to place the Ostreopsis cells observed in the samples within the O. cf. ovata phylospecies’ complex, discarding the presence of O. cf. siamensis. The (phylo)genetic results point out a certain level of nucleotide sequence divergence, but were inconclusive in relation to a possible geographic origin of the O. cf. ovata population from the Algarve’s bloom.

Taxonomic resolutions based on 18S rRNA genes: a case study of subclass copepoda

Biodiversity studies are commonly conducted using 18S rRNA genes. In this study, we compared the inter-species divergence of variable regions (V1-9) within the copepod 18S rRNA gene, and tested their taxonomic resolutions at different taxonomic levels. Our results indicate that the 18S rRNA gene is a good molecular marker for the study of copepod biodiversity, and our conclusions are as follows: 1) 18S rRNA genes are highly conserved intra-species (intra-species similarities are close to 100%); and could aid in species-level analyses, but with some limitations; 2) nearly-whole-length sequences and some partial regions (around V2, V4, and V9) of the 18S rRNA gene can be used to discriminate between samples at both the family and order levels (with a success rate of about 80%); 3) compared with other regions, V9 has a higher resolution at the genus level (with an identification success rate of about 80%); and 4) V7 is most divergent in length, and would be a good candidate marker for the phylogenetic study of Acartia species. This study also evaluated the correlation between similarity thresholds and the accuracy of using nuclear 18S rRNA genes for the classification of organisms in the subclass Copepoda. We suggest that sample identification accuracy should be considered when a molecular sequence divergence threshold is used for taxonomic identification, and that the lowest similarity threshold should be determined based on a pre-designated level of acceptable accuracy.

Intra-genomic variation in symbiotic dinoflagellates: recent divergence or recombination between lineages?

Background

The symbiosis between corals and the dinoflagellate alga Symbiodinium is essential for the development and survival of coral reefs. Yet this fragile association is highly vulnerable to environmental disturbance. A coral’s ability to tolerate temperature stress depends on the fitness of its resident symbionts, whose thermal optima vary extensively between lineages. However, the in hospite population genetic structure of Symbiodinium is poorly understood and mostly based on analysis of bulk DNA extracted from thousands to millions of cells. Using quantitative single-cell PCR, we enumerated DNA polymorphisms in the symbionts of the reef-building coral Pocillopora damicornis, and applied a model selection approach to explore the potential for recombination between coexisting Symbiodinium populations.

Results

Two distinct Symbiodinium ITS2 sequences (denoted C100 and C109) were retrieved from all P. damicornis colonies analysed. However, the symbiont assemblage consisted of three distinct Symbiodinium populations: cells featuring pure arrays of ITS2 type C109, near-homogeneous cells of type C100 (with trace ITS2 copies of type C109), and those with co-dominant C100 and C109 ITS2 repeats. The symbiont consortia of some colonies consisted almost entirely of these putative C100 × C109 recombinants.

Conclusions

Our results are consistent with the occurrence of sexual recombination between Symbiodinium types C100 and C109. While the multiple-copy nature of the ITS2 dictates that the observed pattern of intra-genomic co-dominance may be a result of incomplete concerted evolution of intra-genomic polymorphisms, this is a less likely explanation given the occurrence of homogeneous cells of the C109 type. Conclusive evidence for inter-lineage recombination and introgression in this genus will require either direct observational evidence or a single-cell genotyping approach targeting multiple, single-copy loci.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0325-1) contains supplementary material, which is available to authorized users.

Testing ecological theories with sequence similarity networks: marine ciliates exhibit similar geographic dispersal patterns as multicellular organisms

BACKGROUND

High-throughput sequencing technologies are lifting major limitations to molecular-based ecological studies of eukaryotic microbial diversity, but analyses of the resulting millions of short sequences remain a major bottleneck for these approaches. Here, we introduce the analytical and statistical framework of sequence similarity networks, increasingly used in evolutionary studies and graph theory, into the field of ecology to analyze novel pyrosequenced V4 small subunit rDNA (SSU-rDNA) sequence data sets in the context of previous studies, including SSU-rDNA Sanger sequence data from cultured ciliates and from previous environmental diversity inventories.

RESULTS

Our broadly applicable protocol quantified the progress in the description of genetic diversity of ciliates by environmental SSU-rDNA surveys, detected a fundamental historical bias in the tendency to recover already known groups in these surveys, and revealed substantial amounts of hidden microbial diversity. Moreover, network measures demonstrated that ciliates are not globally dispersed, but are structured by habitat and geographical location at intermediate geographical scale, as observed for bacteria, plants, and animals.

CONCLUSIONS

Currently available 'universal' primers used for local in-depth sequencing surveys provide little hope to exhaust the significantly higher ciliate (and most likely microbial) diversity than previously thought. Network analyses such as presented in this study offer a promising way to guide the design of novel primers and to further explore this vast and structured microbial diversity.

Comparison of potential diatom 'barcode' genes (the 18S rRNA gene and ITS, COI, rbcL) and their effectiveness in discriminating and determining species taxonomy in the Bacillariophyta

Diatoms form an enormous group of photoautotrophic micro-eukaryotes and play a crucial role in marine ecology. In this study, we evaluated typical genes to determine whether they were effective at different levels of diatom clustering analysis to assess the potential of these regions for barcoding taxa. Our test genes included nuclear rRNA genes (the nuclear small-subunit rRNA gene and the 5.8S rRNA gene+ITS-2), a mitochondrial gene (cytochrome c-oxidase subunit 1, COI), a chloroplast gene [ribulose-1,5-biphosphate carboxylase/oxygenase large subunit (rbcL)] and the universal plastid amplicon (UPA). Calculated genetic divergence was highest for the internal transcribed spacer (ITS; 5.8S+ITS-2) (p-distance of 1.569, 85.84% parsimony-informative sites) and COI (6.084, 82.14%), followed by the 18S rRNA gene (0.139, 57.69%), rbcL (0.120, 42.01%) and UPA (0.050, 14.97%), which indicated that ITS and COI were highly divergent compared with the other tested genes, and that their nucleotide compositions were variable within the whole group of diatoms. Bayesian inference (BI) analysis showed that the phylogenetic trees generated from each gene clustered diatoms at different phylogenetic levels. The 18S rRNA gene was better than the other genes in clustering higher diatom taxa, and both the 18S rRNA gene and rbcL performed well in clustering some lower taxa. The COI region was able to barcode species of some genera within the Bacillariophyceae. ITS was a potential marker for DNA based-taxonomy and DNA barcoding of Thalassiosirales, while species of Cyclotella, Skeletonema and Stephanodiscus gathered in separate clades, and were paraphyletic with those of Thalassiosira. Finally, UPA was too conserved to serve as a diatom barcode.

ITS1: a DNA barcode better than ITS2 in eukaryotes?

A DNA barcode is a short piece of DNA sequence used for species determination and discovery. The internal transcribed spacer (ITS/ITS2) region has been proposed as the standard DNA barcode for fungi and seed plants and has been widely used in DNA barcoding analyses for other biological groups, for example algae, protists and animals. The ITS region consists of both ITS1 and ITS2 regions. Here, a large-scale meta-analysis was carried out to compare ITS1 and ITS2 from three aspects: PCR amplification, DNA sequencing and species discrimination, in terms of the presence of DNA barcoding gaps, species discrimination efficiency, sequence length distribution, GC content distribution and primer universality. In total, 85 345 sequence pairs in 10 major groups of eukaryotes, including ascomycetes, basidiomycetes, liverworts, mosses, ferns, gymnosperms, monocotyledons, eudicotyledons, insects and fishes, covering 611 families, 3694 genera, and 19 060 species, were analysed. Using similarity-based methods, we calculated species discrimination efficiencies for ITS1 and ITS2 in all major groups, families and genera. Using Fisher's exact test, we found that ITS1 has significantly higher efficiencies than ITS2 in 17 of the 47 families and 20 of the 49 genera, which are sample-rich. By in silico PCR amplification evaluation, primer universality of the extensively applied ITS1 primers was found superior to that of ITS2 primers. Additionally, shorter length of amplification product and lower GC content was discovered to be two other advantages of ITS1 for sequencing. In summary, ITS1 represents a better DNA barcode than ITS2 for eukaryotic species.

Capturing the population structure of microparasites: using ITS-sequence data and a pooled DNA approach

The internal transcribed spacer (ITS) region of nuclear ribosomal DNA is a common marker not only for the molecular identification of different taxa and strains, but also for the analysis of population structure of wild microparasite communities. Importantly, the multicopy nature of this region allows the amplification of low-quantity samples of the target DNA, a common problem in studies on unicellular, unculturable microparasites. We analysed ITS sequences from the protozoan parasite Caullerya mesnili (class Ichthyosporea) infecting waterflea (Daphnia) hosts, across several host population samples. We showed that analysing representative ITS-types [as identified by statistical parsimony network analysis (SPN)] is a suitable method to address relevant polymorphism. The spatial patterns were consistent regardless of whether parasite DNA was extracted from individual hosts or pooled host samples. Remarkably, the efficiency in detecting different sequence types was even higher after sample pooling. As shown by simulations, an easily manageable number of sequences from pooled DNA samples are sufficient to resolve the spatial population structure in this system. In summary, the ITS region analysed from pooled DNA samples can provide valuable insights into the spatial and temporal dynamics of microparasites. Moreover, the application of SPN analysis is a good alternative to the well-established neighbour-joining method (NJ) for the identification of representative ITS-types. SPN can even outperform NJ by joining most of the singleton sequences to representative sequence clusters.

DNA barcoding species in Alexandrium tamarense complex using ITS and proposing designation of five species

Alexandrium species can be very difficult to identify, with A. catenella, A. tamarense, and A. fundyense that compose "Alexandrium tamarense species complex" (Atama complex) as a distinct example. DNA barcoding is promising to offer a solution but remains to be established. In this study, we examined the utility of ITS in resolving the Atama species complex, by analyzing previously studied strains plus unstudied Chinese strains within the LSU- and SSU-rDNA based group/clade frameworks recently established. We further investigated the presence of intragenomic polymorphism and its implications in species delimitation. Similar to the previous SSU and LSU results, our ITS-based phylogenies divided the complex to five clusters, but with longer and evener branch lengths between the clusters. Based on the ITS region, the inter-cluster genetic distances (p=0.134-0.216) were consistently and substantially greater than intra-cluster genetic distances (p=0.000-0.066), with an average inter-cluster (species) distance (p=0.167) 7.6-fold of the average intraspecific difference (p=0.022), qualifying the approximately 510-520bp ITS as a DNA barcode for Atama complex. We detected varying levels of intragenomic polymorphism in ITS but found that this did not impact the taxon-resolving power of this gene. With this DNA barcode, the new East and South China Sea strains and one Antarctic strain were placed in Clade IIC/Group IV, even though there were 7-10 polymorphic sites in their ITS, in contrast to none in SSU. Furthermore, our results suggest that the five clusters are recognizable as distinct species according to the phylogenetic species concept. Based on the phylogenetic placements of the type-locality strains of the existing three morphospecies and the dominant localities of other strains, we propose that Group I/Clade I be designated as A. fundyense, Group III/Clade IIB as A. tamarense, Group IV/Clade IIC as A. catenella, Group II/Clade IIA as A. mediterranis, and Group V/Clade IID as A. australis.

An Updated List of Generic Names in the Thoracosphaeraceae

Calcareous dinophytes produce exoskeletal calcified structures during their life history (a unique character among the alveolates) and are subsumed under the Thoracosphaeraceae as part of the Peridiniales. We provide a brief synopsis about the taxonomic history of the group, from the first descriptions of fossils in the 19th century through to the results of molecular phylogenetics studies undertaken during the past two decades. Delimitation and circumscription of the Thoracosphaeraceae are challenging, as they comprise both phototrophic (presumably including endosymbiotic) as well as heterotrophic (and even parasitic) dinophytes from marine and freshwater environments, respectively. However, calcareous structures are not known from all members of the Thoracosphaeraceae, and the corresponding species and groups are considered to have lost the capacity to calcify. Five years ago, a taxonomic list of 99 generic names assigned to the Thoracosphaeraceae was published, and we update this compendium with 19 additional names based on recent studies.

Single-Cell DNA barcoding using sequences from the small subunit rRNA and internal transcribed spacer region identifies new species of Trichonympha and Trichomitopsis from the hindgut of the termite Zootermopsis angusticollis

To aid in their digestion of wood, lower termites are known to harbour a diverse community of prokaryotes as well as parabasalid and oxymonad protist symbionts. One of the best-studied lower termite gut communities is that of Zootermopsis angusticollis which has been known for almost 100 years to possess 3 species of Trichonympha (T. campanula, T. collaris, and T. sphaerica), 1 species of Trichomitopsis (T. termopsidis), as well as smaller flagellates. We have re-assessed this community by sequencing the small subunit (SSU) rRNA gene and the internal transcribed spacer (ITS) region from a large number of single Trichonympha and Trichomitopsis cells for which morphology was also documented. Based on phylogenetic clustering and sequence divergence, we identify 3 new species: Trichonympha postcylindrica, Trichomitopsis minor, and Trichomitopsis parvus spp. nov. Once identified by sequencing, the morphology of the isolated cells for all 3 new species was re-examined and found to be distinct from the previously described species: Trichonympha postcylindrica can be morphologically distinguished from the other Trichonympha species by an extension on its posterior end, whereas Trichomitopsis minor and T. parvus are smaller than T. termopsidis but similar in size to each other and cannot be distinguished based on morphology using light microscopy. Given that Z. angusticollis has one of the best characterized hindgut communities, the near doubling of the number of the largest and most easily identifiable symbiont species suggests that the diversity of hindgut symbionts is substantially underestimated in other termites as well. Accurate descriptions of the diversity of these microbial communities are essential for understanding hindgut ecology and disentangling the interactions among the symbionts, and molecular barcoding should be a priority for these systems.