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.     

Genetic identification of bucktooth parrotfish Sparisoma radians (Valenciennes, 1840) (Labridae,Scarinae) by chromosomal and molecular markers

Parrotfishes (Labridae, Scarinae) comprise a large marine fish group of difficult identification, particularly during juvenile phase when the typical morphology and coloration of adults are absent. Therefore, the goal of this study was to test cytogenetic markers and DNA barcoding in the identification of bucktooth parrtotfish Sparisoma radians from the northeastern coast of Brazil. Sequencing of cytochrome c oxidase subunit I (COI) confirmed all studied samples as S. radians, and all showed high similarity (99–100%) with Caribbean populations. The karyotype of this species was divergent from most marine Perciformes, being composed of 2n = 46 chromosomes. These consisted of a large number of metacentric and submetacentric pairs with small amounts of heterochromatin and GC-rich single nucleolar organizer regions (NORs) not syntenic to 5S rDNA clusters. These are the first data about DNA barcoding in parrotfish from the Brazilian province and the first refined chromosomal analysis in Scarinae, providing useful data to a reliable genetic identification of S. radians.     

rpoB-Based Microbial Community Analysis Avoids Limitations Inherent in 16S rRNA Gene Intraspecies Heterogeneity

Contemporary microbial community analysis frequently involves PCR-amplified sequences of the 16S rRNA gene (rDNA). However, this technology carries the inherent problem of heterogeneity between copies of the 16S rDNA in many species. As an alternative to 16S rDNA sequences in community analysis, we employed the gene for the RNA polymerase beta subunit (rpoB), which appears to exist in one copy only in bacteria. In the present study, the frequency of 16S rDNA heterogeneity in bacteria isolated from the marine environment was assessed using bacterial isolates from the red alga Delisea pulchra and from the surface of a marine rock. Ten strains commonly used in our laboratory were also assessed for the degree of heterogeneity between the copies of 16S rDNA and were used to illustrate the effect of this heterogeneity on microbial community pattern analysis. The rock isolates and the laboratory strains were also used to confirm nonheterogeneity of rpoB, as well as to investigate the versatility of the primers. In addition, a comparison between 16S rDNA and rpoB PCR-DGGE (denaturing gradient gel electrophoresis)-based community analyses was performed using a DNA mixture of nine isolates from D. pulchra. Eight out of 14 isolates from D. pulchra, all rock isolates, and 6 of 10 laboratory strains displayed multiple bands for 16S rDNA when analyzed by DGGE. There was no indication of heterogeneity for either the rock isolates or the laboratory strains when rpoB was used for PCR-DGGE analysis. Microbial community pattern analysis using 16S rDNA PCR-DGGE showed an overestimation of the number of laboratory strains in the sample, while some strains were not represented. Therefore, the 16S rDNA PCR-DGGE-based community analysis was proven to be severely limited by 16S rDNA heterogeneity. The mixture of isolates from D. pulchra proved to be more accurately described using rpoB, compared to the 16S rDNA-based PCR-DGGE.     

DISCRIMINATIVE POWER OF NUCLEAR rDNA SEQUENCES FOR THE DNA TAXONOMY OF THE DINOFLAGELLATE GENUS PERIDINIUM (DINOPHYCEAE)1

The genus Peridinium Ehrenb. comprises a group of highly diversified dinoflagellates. Their morphological taxonomy has been established over the last century. Here, we examined relationships within the genus Peridinium, including Peridinium bipes F. Stein sensu lato, based on a molecular phylogeny derived from nuclear rDNA sequences. Extensive rDNA analyses of nine selected Peridinium species showed that intraspecies genetic variation was considerably low, but interspecies genetic divergence was high (>1.5% dissimilarity in the nearly complete 18S sequence; >4.4% in the 28S rDNA D1/D2). The 18S and 28S rDNA Bayesian tree topologies showed that Peridinium species grouped according to their taxonomic positions and certain morphological characters (e.g., epithecal plate formula). Of these groups, the quinquecorne group (plate formula of 3′, 2a, 7″) diverged first, followed by the umbonatum group (4′, 2a, 7″) and polonicum group (4′, 1a, 7″). Peridinium species with a plate formula of 4′, 3a, 7″ diverged last. Thus, 18S and 28S rDNA D1/D2 sequences are informative about relationships among Peridinium species. Statistical analyses revealed that the 28S rDNA D1/D2 region had a significantly higher genetic divergence than the 18S rDNA region, suggesting that the former as DNA markers may be more suitable for sequence‐based delimitation of Peridinium. The rDNA sequences had sufficient discriminative power to separate P. bipes f. occultaum (Er. Lindem.) M. Lefèvre and P. bipes f. globosum Er. Lindem. into two distinct species, even though these taxa are morphologically only marginally discriminated by spines on antapical plates and the shape of red bodies during the generation of cysts. Our results suggest that 28S rDNA can be used for all Peridinium species to make species‐level taxonomic distinctions, allowing improved taxonomic classification of Peridinium.     

Development of a Real-Time PCR Assay for Rapid Detection and Quantification of Alexandrium minutum (a Dinoflagellate)

The marine dinoflagellate genus Alexandrium includes a number of species which produce neurotoxins responsible for paralytic shellfish poisoning (PSP), which in humans may cause muscular paralysis, neurological symptoms, and, in extreme cases, death. A. minutum is the most widespread toxic PSP species in the western Mediterranean basin. The monitoring of coastal waters for the presence of harmful algae also normally involves microscopic examinations of phytoplankton populations. These procedures are time consuming and require a great deal of taxonomic experience, thus limiting the number of specimens that can be analyzed. Because of the genetic diversity of different genera and species, molecular tools may also help to detect the presence of target microorganisms in marine field samples. In this study, we developed a real-time PCR-based assay for rapid detection of all toxic species of the Alexandrium genus in both fixative-preserved environmental samples and cultures. Moreover, we developed a real-time quantitative PCR assay for the quantification of A. minutum cells in seawater samples. Alexandrium genus-specific primers were designed on the 5.8S rDNA region. Primer specificity was confirmed by using BLAST and by amplification of a representative sample of the DNA of other dinoflagellates and diatoms. Using a standard curve constructed with a plasmid containing the ITS1-5.8S-ITS2 A. minutum sequence and cultured A. minutum cells, we determined the absolute number of 5.8S rDNA copies per cell. Consequently, after quantification of 5.8S rDNA copies in samples containing A. minutum cells, we were also able to estimate the number of cells. Several fixed A. minutum bloom sea samples from Arenys Harbor (Catalan Coast, Spain) were analyzed using this method, and quantification results were compared with standard microscopy counting methods. The two methods gave comparable results, confirming that real-time PCR could be a valid, fast alternative procedure for the detection and quantification of target phytoplankton species during coastal water monitoring.     

Identification of species in the genus Pichia by restriction of the internal transcribed spacers (ITS1 and ITS2) and the 5.8S ribosomal DNA gene

In this study, the variability within the ribosomal DNA region spanning the internal transcribed spacers ITS1 and ITS2 and the 5.8S gene (5.8S-ITS rDNA) was used to differentiate species in the genus Pichia. The 5.8S-ITS rDNA region was PCR-amplified and the PCR product digested with the enzymes CfoI, HinfI, and HaeIII. The variability in the size of the amplified product and in the restriction patterns enabled differentiation between species in the genus Pichia, and between Pichia species and yeast species from other genera in the Yeast-id database (). Moreover, the restriction fragment length polymorphism (RFLP) patterns of the 5.8S-ITS enabled misidentified strains to be detected and revealed genetic heterogeneity between strains within the Pichia membranifaciens and Pichia nakazawae species. Ultimately, the RFLP patterns of the 5.8S-ITS rDNA failed to differentiate between some Pichia and Candida species that could be distinguished on the basis of the sequence of the 5.8S-ITS rRNA region or the sequence of the D1/D2 domain of the 26S rDNA gene.     

Characterization of Trichodesmium spp. by Genetic Techniques

The genetic diversity of Trichodesmium spp. from natural populations (off Bermuda in the Sargasso Sea and off North Australia in the Arafura and Coral Seas) and of culture isolates from two regions (Sargasso Sea and Indian Ocean) was investigated. Three independent techniques were used, including a DNA fingerprinting method based on a highly iterated palindrome (HIP1), denaturing gradient gel electrophoresis of a hetR fragment, and sequencing of the internal transcribed spacer (ITS) of the 16S-23S rDNA region. Low genetic diversity was observed in natural populations of Trichodesmium spp. from the two hemispheres. Culture isolates of Trichodesmium thiebautii, Trichodesmium hildebrandtii, Trichodesmium tenue, and Katagnymene spiralis displayed remarkable similarity when these techniques were used, suggesting that K. spiralis is very closely related to the genus Trichodesmium. The largest genetic variation was found between Trichodesmium erythraeum and all other species of Trichodesmium, including a species of Katagnymene. Our data obtained with all three techniques suggest that there are two major clades of Trichodesmium spp. The HIP1 fingerprinting and ITS sequence analyses allowed the closely related species to be distinguished. This is the first report of the presence of HIP1 in marine cyanobacteria.     

Description and genetic characterisation of Hysterothylacium (Nematoda: Raphidascarididae) larvae parasitic in Australian marine fishes

Nematodes belonging to the genus Hysterothylacium (family Raphidascarididae) infect various species of marine fish in both the larval and adult stages. Humans can be accidentally infected upon eating infected seafood. In spite of their importance, relatively little is known of their occurrence and systematics in Australia. An examination of various species of marine teleosts in Australian waters revealed a high prevalence of Hysterothylacium larval types. In the present study, seven previously undescribed Hysterothylacium larval morphotypes (V to VII and IX to XII) were discovered. In total we found 10 different morphotypes and we genetically characterised nine morphotypes identified. A morphological dichotomous identification key has been established to differentiate these morphotypes. Since some larvae of Hysterothylacium from marine fishes cannot be differentiated morphologically from other nematode larvae, such as Paraheterotyphlum, Heterotyphlum, Iheringascaris and Lapetascaris, the first and second internal transcribed spacers (ITS-1 and ITS-2) of nuclear ribosomal DNA (rDNA) of these larvae were characterised to confirm their taxonomic status. This genetic characterisation implied that some distinct morphotypes belong to different developmental stages of the same species. In addition, it revealed that some morphotypes can comprise distinct genotypes. No match was found between ITS-1 and ITS-2 sequences obtained from larvae in the present study and those from adults available in the GenBank, highlighting the lack of knowledge on occurrence of adult nematodes infecting Australian fish.     

Real-Time Quantitative PCR for Assessment of Abundance of Pseudoalteromonas Species in Marine Samples

A real-time quantitative PCR (RTQ-PCR) method for measuring the abundance of Pseudoalteromonas species in marine samples is presented. PCR primers targeting a Pseudoalteromonas-specific region of the 16S rRNA gene were tested at three different levels using database searches (in silico), a selection of pure cultures (in vitro), and a combined denaturing gradient gel electrophoresis and cloning approach on environmental DNA (in situ). The RTQ-PCR method allowed for the detection of SYBR Green fluorescence from double-stranded DNA over a linear range spanning six orders of magnitude. The detection limit was determined as 1.4 fg of target DNA (1,000 gene copies) measured in the presence of 20 ng of nontarget DNA from salmon testes. In this study, we discuss the importance of robust post-PCR analyses to overcome pitfalls in RTQ-PCR when samples from different complex marine habitats are analyzed and compared on a nonroutine basis. Representatives of the genus Pseudoalteromonas were detected in samples from all investigated habitats, suggesting a widespread distribution of this genus across many marine habitats (e.g., seawater, rocks, macroalgae, and marine animals). Three sample types were analyzed by RTQ-PCR to determine the relative abundance of Pseudoalteromonas ribosomal DNA (rDNA) compared to the total abundance of eubacterial rDNA. The rDNA fractions of Pseudoalteromonas compared to all Eubacteria were 1.55% on the green alga Ulva lactuca, 0.10% on the tunicate Ciona intestinalis, and 0.06% on the green alga Ulvaria fusca.     

Construction and characterization of a BAC library from the Coffea arabica genotype Timor Hybrid CIFC 832/2

Most of the world’s coffee production originates from Coffea arabica, an allotetraploid species with low genetic diversity and for which few genomic resources are available. Genomic libraries with large DNA fragment inserts are useful tools for the study of plant genomes, including the production of physical maps, integration studies of physical and genetic maps, genome structure analysis and gene isolation by positional cloning. Here, we report the construction and characterization of a Bacterial Artificial Chromosome (BAC) library from C. arabica Timor Hybrid CIFC 832/2, a parental genotype for several modern coffee cultivars. The BAC library consists of 56,832 clones with an average insert size of 118 kb, which represents a dihaploid genome coverage of five to sixfold. The content of organellar DNA was estimated at 1.04 and 0.5 % for chloroplast and mitochondrial DNA, respectively. The BAC library was screened for the NADPH-dependent mannose-6-phosphate reductase gene (CaM6PR) with markers positioned on four linkage groups of a partial C. arabica genetic map. A mixed approach using PCR and membrane hybridization of BAC pools allowed for the discovery of nine BAC clones with the CaM6PR gene and 53 BAC clones that were anchored to the genetic map with simple sequence repeat markers. This library will be a useful tool for future studies on comparative genomics and the identification of genes and regulatory elements controlling major traits in this economically important crop species.     

Nontranscribed spacers in Drosophila ribosomal DNA

Ribosomal DNA nontranscribed spacers in Drosophila virilis DNA have been examined in some detail by restriction site analysis of cloned segments of rDNA, nucleic acid hybridizations involving unfractionated rDNA, and base composition estimates. The overall G+C content of the spacer is 27–28%; this compares with 39% for rDNA as a whole, 40% for main band DNA, and 26% for the D. virilis satellites. Much of the spacer is comprised of 0.25 kb repeats revealed by digestion with Msp I, Fnu DII or Rsd I, which terminate very near the beginning of the template for the ribosomal RNA precursor. The spacers are heterogeneous in length among rDNA repeats, and this is largely accounted for by variation among rDNA units in the number of 0.25 kb elements per spacer. Despite its high A+T content and the repetitive nature of much of the spacer, and the proximity of rDNA and heterochromatin in Drosophila, pyrimidine tract analysis gave no indication of relatedness between the spacer and satellite DNA sequences. Species of Drosophila closely related to D. virilis have rDNA spacers that are homologous with those in D. virilis to the extent that hybridization of a cloned spacer segment of D. virilis rDNA to various DNA is comparable with hybridization to homologous DNA, and distributions of restriction enzyme cleavage sites are very similar (but not identical) among spacers of the various species. There is spacer length heterogeneity in the rDNA of all species, and each species has a unique major rDNA spacer length. Judging from Southern blot hybridization, D. hydei rDNA spacers have 20–30% sequence homology with D. virilis rDNA spacers, and a repetitive component is similarly sensitive to Msp I and Fnu DII digestion, D. melanogaster rDNA spacers have little or no homology with counterparts in D. virilis rDNA, despite a similar content of 0.25 kb repetitive elements. In contrast, sequences in rDNA that encode 18S and 28S ribosomal RNA have been highly conserved during the divergence of Drosophila species; this is inferred from interspecific hybridizations involving ribosomal RNA and a comparison of distributions of restriction enzyme cleavage sites in rDNA.Dedicated to Professor Wolfgang Beermann on the occasion of his sixtieth birthday     

Phylogenetic analysis of the pathogenic bacteria Spiroplasma penaei based on multilocus sequence analysis

A pathogenic Spiroplasma penaei strain was isolated from the hemolymph of moribund Pacific white shrimp, Penaeus vannamei. The shrimp sample originated from a shrimp farm near Cartagena, Colombia, that was suffering from high mortalities in ponds with very low salinity and high temperatures. This new emerging disease in a marine crustacean in the Americas is described as a systemic infection. The multilocus phylogenetic analysis suggests that S. penaei strain has a terrestrial origin. Evolutionary relationship trees, based on five partial DNA sequences of 16S rDNA, 23S rDNA, 5S rDNA, gyrB, rpoB genes and two complete DNA sequences of 16S-23S rDNA and 23S-5S rDNA intergenic spacer region, were reconstructed using the distance-based Neighboring-Joining (NJ) method with Kimura-2-parameter substitution model. The NJ trees based on all DNA sequences investigated in this study positioned S. penaei in the Citri-Poulsonii clade and corroborate the observations by other investigators using the 16S rDNA gene. Pairwise genetic distance calculation between sequences of spiroplasmas showed S. penaei to be closely related to Spiroplasma insolitum and distantly related to Spiroplasma sp. SHRIMP from China.     

Differentiation of Species and Populations of Aphelenchoides and of Ditylenchus angustus Using a Fragment of Ribosomal DNA

The polymerase chain reaction (PCR) was used to amplify a fragment of the ribosomal DNA (rDNA) from species and undescribed populations of Aphelenchoides and Ditylenchus angustus. The PCR primers used were based on conserved sequences in the 18S and 26S ribosomal RNA genes of Caenorhabditis elegans. In C. elegans, these primers amplify a 1,292 base pair (bp) fragment, which consists of the two internal transcribed spacers and the entire 5.8S gene. Amplification products from crude DNA preparations of 12 species and populations of Aphelenchoides and from D. angustus ranged in size from approximately 860-1,100bp. Southern blots probed with a cloned ribosomal repeat from C. elegans confirmed the identity of these amplified bands as ribosomal fragments. In addition to the differing sizes of the amplified rDNA fragments, the relative intensity of hybridization with the C. elegans probe indicated varying degrees of sequence divergence between species and populations. In some cases, amplified rDNA from the fungal host was evident. Storage of A. composticola at - 45 C for 2 years did not affect the ability to obtain appropriate amplified products from crude DNA preparations. Amplified rDNA fragments were cut with six restriction enzymes, and the restriction fragments produced revealed useful diagnostic differences between species and some undescribed populations. These results were consistent with previous studies based on morphology and isoenzymes. Three undescribed populations of Aphelenchoides were found to be different from all the species examined and from each other.     

Detection and quantification of cultured marine Alexandrium species by real-time PCR

The occurrence of harmful algal blooms (HABs) throughout the world has increased and poses a large threat to human health, fishery resources and tourism industries. The genus Alexandrium includes a number of toxic species associated with HABs. Therefore, it is very important to rapidly detect and monitor the harmful algae, such as Alexandrium genus. In this study, a standard curve of plasmid containing 18S rDNA-28S rDNA region from Alexandrium catenella was constructed and 5.8S rDNA sequence served as the primer of the real-time PCR. Cultured A. catenella, Alexandrium affine, Alexandrium lusitanicum and Alexandrium minutum samples were analyzed by real-time PCR using the same set of primers simultaneously. Using microscopy cells counts, 5.8S rDNA copies per cell and total DNA per cell were estimated. This assay method is promising for rapid detection of large number of Alexandrium samples.     

Oramoeba fumarolia gen. nov., sp. nov., a new marine heterolobosean amoeboflagellate growing at 54 °C

An amoeba strain was isolated from marine sediment taken from the beach near a fumarole in Italy. The trophozoites of this new marine species transforms into flagellates with variable numbers of flagella, from 2 to 10. The strain forms round to oval cysts. This thermophilic amoeboflagellate grows at temperatures up to 54 °C. Molecular phylogenetic analysis of the small subunit ribosomal DNA (SSU rDNA) places the amoeboflagellate in the Heterolobosea. The closest relatives are Stachyamoeba sp. ATCC50324, a strain isolated from an ocean sample, and Vrihiamoeba italica, a recent isolate from a rice field. Like some other heterolobosean species, this new isolate has a group I intron in the SSU rDNA. Because of the unique place in the molecular phylogenetic tree, and because there is no species found in the literature with similar morphological and physiological characteristics, this isolate is considered to be a new genus and a new species, Oramoeba fumarolia gen. nov., sp. nov.     

Characterization and Sequence Variation in the rDNA Region of Six Nematode Species of the Genus Longidorus (Nematoda)

Total DNA was isolated from individual nematodes of the species Longidorus helveticus, L. macrosoma, L. arthensis, L. profundorum, L. elongatus, and L. raskii collected in Switzerland. The ITS region and D1-D2 expansion segments of the 26S rDNA were amplified and cloned. The sequences obtained were aligned in order to investigate sequence diversity and to infer the phylogenetic relationships among the six Longidorus species. D1-D2 sequences were more conserved than the ITS sequences that varied widely in primary structure and length, and no consensus was observed. Phylogenetic analyses using the neighbor-joining, maximum parsimony and maximum likelihood methods were performed with three different sequence data sets: ITS1-ITS2, 5.8S-D1-D2, and combining ITS1-ITS2+5.8S-D1-D2 sequences. All multiple alignments yielded similar basic trees supporting the existence of the six species established using morphological characters. These sequence data also provided evidence that the different regions of the rDNA are characterized by different evolution rates and by different factors associated with the generation of extreme size variation.     

Molecular cytogenetic characterisation of Salix viminalis L. using repetitive DNA sequences

Salix viminalis L. (2n?=?38) is a diploid dicot species belonging to the Salix genus of the Salicaceae family. This short-rotation woody crop is one of the most important renewable bioenergy resources worldwide. In breeding for high biomass productivity, limited knowledge is available on the molecular cytogenetics of willow, which could be combined with genetic linkage mapping. The present paper describes the adaptation of a fluorescence in situ hybridisation (FISH) protocol as a new approach to analyse the genomic constitution of Salix viminalis using the heterologous DNA clones pSc119.2, pTa71, pTa794, pAs1, Afa-family, pAl1, HT100.3, ZCF1 and the GAA microsatellite marker. Three of the nine probes showed unambiguous signals on the metaphase chromosomes. FISH analysis with the pTa71 probe detected one major 18S-5.8S-26S rDNA locus on the short arm of one chromosome pair; however, the pTa794 rDNA site was not visible. One chromosome pair showed a distinct signal around the centromeric region after FISH with the telomere-specific DNA clone HT100.3. Two chromosome pairs were found to have pAs1 FISH signals, which represent a D-genome-specific insert from Aegilops tauschii. Based on the FISH study, a set of chromosomes with characteristic patterns is presented, which could be used to establish the karyotype of willow species.     

Henneguya tunisiensis n. sp. (Myxosporea: Bivalvulida), a new gill parasite of Symphodus tinca (L.) (Teleostei: Labridae) off Tunisia

Henneguya tunisiensis n. sp., a new myxosporean, is described from the gill-arches of the East Atlantic peacock wrasse Symphodus tinca (L.) collected from off the Kerkennah Islands, Tunisia. It is characterised by the presence of elongate white plasmodia of 1–1.5 × 1.5–2 mm in size. The mature spores are rounded in frontal view and have two identical polar capsules and two caudal appendages which taper considerably at the end. Both light and electron microscopical data show that this species differs in several morphological features from all previously described Henneguya spp. A molecular analysis, based on 18S rDNA sequence data, indicates that H. tunisiensis n. sp. is readily distinguishable from other myxozoan DNA sequences in GenBank. Phylogenetically, the new species is placed in the marine Henneguya clade, which is a sister group of marine Myxobolus spp. from perciform fishes in Tunisian waters.     

Genetic Diversity and Distribution of the Ciguatera-Causing Dinoflagellate Gambierdiscus spp. (Dinophyceae) in Coastal Areas of Japan

Background

The marine epiphytic dinoflagellate genus Gambierdiscus produce toxins that cause ciguatera fish poisoning (CFP): one of the most significant seafood-borne illnesses associated with fish consumption worldwide. So far, occurrences of CFP incidents in Japan have been mainly reported in subtropical areas. A previous phylogeographic study of Japanese Gambierdiscus revealed the existence of two distinct phylotypes: Gambierdiscus sp. type 1 from subtropical and Gambierdiscus sp. type 2 from temperate areas. However, details of the genetic diversity and distribution for Japanese Gambierdiscus are still unclear, because a comprehensive investigation has not been conducted yet.

Methods/Principal Finding

A total of 248 strains were examined from samples mainly collected from western and southern coastal areas of Japan during 2006–2011. The SSU rDNA, the LSU rDNA D8–D10 and the ITS region were selected as genetic markers and phylogenetic analyses were conducted. The genetic diversity of Japanese Gambierdiscus was high since five species/phylotypes were detected: including two reported phylotypes (Gambierdiscus sp. type 1 and Gambierdiscus sp. type 2), two species of Gambierdiscus (G. australes and G. cf. yasumotoi) and a hitherto unreported phylotype Gambierdiscus sp. type 3. The distributions of type 3 and G. cf. yasumotoi were restricted to the temperate and the subtropical area, respectively. On the other hand, type 1, type 2 and G. australes occurred from the subtropical to the temperate area, with a tendency that type 1 and G. australes were dominant in the subtropical area, whereas type 2 was dominant in the temperate area. By using mouse bioassay, type 1, type 3 and G. australes exhibited mouse toxicities.

Conclusions/Significance

This study revealed a surprising diversity of Japanese Gambierdiscus and the distribution of five species/phylotypes displayed clear geographical patterns in Japanese coastal areas. The SSU rDNA and the LSU rDNA D8–D10 as genetic markers are recommended for further use.