Abundance of the benthic dinoflagellate Prorocentrum and the diversity, distribution, and diarrhetic shellfish toxin production of Prorocentrum lima complex and P. caipirignum in Japan

•Prorocentrum cell density in subtropical was higher than that in temperate Japan.•Japanese P. lima complex was separated into two clades.•The two clades and P. caipirignum were separated into several subclades.•All the subclades produced OA, whereas some produced DTX1.•Toxic subclades of P. lima complex were widespread in subtropical and temperate Japan. In the present study, the abundance of Prorocentrum and the molecular phylogeny, distribution, and DST production of P. lima complex and P. caipirignum in Japan were investigated. First, the cell densities of Prorocentrum were assessed from the temperate to subtropical zones in Japan between 2014 and 2018. The cell density in the subtropical zone [19.0 ± 40.2 cells/g wet weight (ww) algae] was significantly higher than that in the temperate zone (1.4 ± 3.4 cells/g ww algae). A total of 244 clonal strains were established from the temperate and subtropical zones. Phylogenetic analyses based on the large-subunit ribosomal DNA D1/D2 revealed that the strains were separated into four species/species complex/phylotypes (P. lima complex, P. caipirignum, and new phylotypes Prorocentrum spp. types 1 and 2). The strains of P. lima complex could be separated into two clades (1 and 3). Furthermore, the strains of clades 1 and 3 could be separated into nine subclades (1a, 1c, 1d, 1e, 1f, 1g, 1h, 1i, and 1j) and three subclades (3a, 3b, and 3c), respectively. The strains of P. caipirignum were separated into two subclades (b and e). Each phylotype/subclade showed a unique distribution pattern in Japan: P. lima complex subclades 1a, 1c, and 3a and P. caipirignum subclades b and e were widespread from the temperate to subtropical zones. On the other hand, P. lima complex subclades 1e and 1i were restricted to the temperate zone, and P. lima complex subclades 1d, 1f, 1g, 1h, 1j, 3b, and 3c and Prorocentrum spp. types 1 and 2 were restricted to the subtropical zone. Furthermore, the DST production of the 243 clonal strains was assessed by LC/MS/MS analysis. The results revealed that all strains produced okadaic acid (OA) and that the OA contents of P. lima complex subclades 1d and 1f, P. caipirignum subclades b and e, and Prorocentrum sp. type 2 tended to be higher than those of the other subclades. While P. lima complex subclades 1a, 1e, 1f, and 1i produced DTX1, the other phylotype/subclades produced either no or low quantities of DTX1. A strain of P. lima complex subclade 1e showed the highest OA and DTX1 contents (55.