Morphological and molecular variability of Peridinium volzii Lemmerm. (Peridiniaceae, Dinophyceae) and its relevance for infraspecific taxonomy

Contemporary delimitation of species and populations in the microbial domain relies on an integrative approach combining molecular and morphological techniques. In case of the dinophyte Peridinium volzii, a considerable number of infraspecific taxonomic entities have been reported, but it is unclear at present whether the corresponding traits are stable within reproductively isolated units or refer to intraspecific variability. We established 26 monoclonal strains from Central Europe with a morphology that is consistent for P. volzii and characterised them by sequences gained from the rRNA operon. Ten of such strains, representative for the entire diversity observed, were investigated in detail morphologically using light and electron microscopy. In the molecular tree, P. volzii was monophyletic, sister group of Peridinium willei, and three ITS ribotypes could be distinguished. Some traits corresponding to previously described varieties and forms were found in individual cells across the strains under investigation, but not as stable characters correlating to certain ribotypes. We also observed new morphological variability (e.g., unusual shape of plate 4″). Cell size and displacement of the cingulum were significantly different between certain ribotypes but in turn, such diagnostic traits are impossible to assign to already described taxa due to their ambiguity. Based on the small first apical plate as diagnostic trait and putative apomorphy, P. volzii is a characteristic species but the present data given, we are reserved to accept more than a single reproductive unit. Thus, more research is necessary, including a focus on species delimitation to putative close relatives such as Peridinium maeandricum.


Introduction
Dinophytes play an equally important role in the ecology of water bodies, although they are not as visually apparent as larger aquatic organisms. Habitat preference is clearly marine, but they are also known to be widely distributed in freshwater (Moestrup & Calado, 2018). As unicellular algae and important bioindicators in aquatic systems, they have the ability to indicate rapid environmental changes within these habitats (Camargo, 1994;Hellawell, 1986). The adaptability upon exposure to environmental and anthropogenic stressors is enabled by a high degree of plasticity in the morphology and molecular constitution (Chevin et al., 2010;Wong & Candolin, 2015) also of dinophytes. However, knowledge of intraspecific variability is currently limited, hindering proper monitoring of unicellular organisms and their ecological impact .
Generally accepted criteria for species delimitation are absent for unicellular organisms such as dinophytes (Boenigk et al., 2012;Mayr, 1982a, b). Crossing experiments are the first choice to test reproductive isolation of populations, but are challenging and therefore rare in dinophytes (Blackburn et al., 2001;Figueroa et al., 2010;Pfiester & Skvarla, 1979;Soyer-Gobillard et al., 2002). The morphological species concept is still very viable for such organisms. However, the quality of a diagnostic trait and its abundance in a certain population is not easy to work out precisely as well. This is especially true for character states that have been used to describe taxa at the various infraspecific ranks such as varieties (var.) and forms (forma). DNA sequencing has become an important tool to help resolve taxonomic issues (Blaxter, 2004;Keck et al., 2018;Kretschmann et al., 2018b;Miller, 2007). Contemporary delimitation of species and populations in the microbial domain therefore relies on an integrative approach combining molecular and morphological techniques. Anyhow, evaluation of any scientific name of microbes described prior to the last millennium has some uncertainty, as it lacks DNA sequence information (Boenigk et al., 2012;De Clerck et al., 2013;Romeikat et al., 2019).
Today, two separated species are accepted constituting the P. willei species group (Bachmann, 1911;Moestrup & Calado, 2018), and P. volzii can be distinguished from P. willei based on its considerably smaller plate 1′. Additionally, P. volzii differs from P. willei by a sulcal extension on the epitheca being longer than wide and not shorter than wide, as well as the absence of hyaline wings on the epitheca (Lindemann, 1920: 146;Moestrup & Calado, 2018). It also appears that ecological differences separate these two species, with P. volzii having more narrow seasonal and pH limits (Olrik, 1992). The occurrence of P. volzii in the plankton seems to be restricted from late-spring to summer in calcium-rich lakes and ponds, rarer in other types of waters (Höll, 1928). Moreover, P. willei occurs in ponds and lakes with transparent through polyhumic brownish water, while P. volzii is mainly found in the plankton of calcium-rich ponds and lakes avoiding humic water. The distinctiveness between P. volzii and P. willei is also displayed in molecular phylogenetic trees Hayhome & Pfiester, 1983;Ki et al., 2011).
Peridinium volzii is considered a morphologically variable species, and seven infraspecific taxa are currently accepted (Moestrup & Calado, 2018 In the original description of P. volzii from Singapore, Lemmermann (1905) noted the presence of a distinct antapical spine, which is however absent from those infraspecific taxa (or supposed synonyms such as P. guestrowiense) described from Central Europe (Lefèvre, 1927;Lindemann, 1916Lindemann, , 1919Lindemann, , 1920. Peridinium volzii appears morphologically highly variable, but no DNA sequence data are linked to the infraspecific and putatively synonymous names assigned to the species. In turn, the DNA sequence of only a single strain (NIES501; Ki et al., 2011) is available, though without knowledge of morphology. In this study, we present the first integrative approach to understand intraspecific variability in P. volzii, providing both morphological and molecular data of newly established monoclonal strains. We show already assessed morphological variation, but also new traits yet not been described in the literature, and the existence of three ITS ribotypes. In the phylogenetic tree, P. volzii and P. willei are clearly delimited species within the Peridiniaceae. We aim at a better knowledge of intraspecific variability within dinophyte species, which may have importance also for other microalgal groups.
Detailed information of the localities and collectors are compiled in Table S1.
Thecate cells were observed under an inverse light microscope (LM) CKX41 (Olympus; Hamburg, Germany). Documentation was performed with a DP73 digital camera (Olympus) and subsequent cell size measurements using the software cellSense (Olympus). Measurements were taken in LM using selected strains representative for the observed diversity. Statistically significant clusters were calculated with Tukey's Honest Significant Difference (HSD) test (p-values < 0.05) and illustrated using the R software v4.0.3 (freely available under www.r-proje ct. org).
For the preparation of permanent slides, cells were fixed with 2.5% glutar(di)aldehyde (agar scientific; Stansted, Essex, UK). Double-staining was carried out using 0.5% (water-based) astra blue in 2% tartaric acid (Fluka; Buchs, Switzerland) in WC medium and 0.1% (ethanol-based) eosin (Merck; Darmstadt, Germany) during a graded ethanol (Roth) series. Ethanol-based Technovit 7100 (Heraeus; Wehrheim, Germany) was used for embedding, following the manufacturer's instructions. For the final preparation, 30 µl aliquots of the Technovit mixture including the embedded samples were transferred to three slides. The material is deposited at the Centre of Excellence for Dinophyte Taxonomy (CEDiT; Wilhelmshaven, Germany), and duplicates are held in Berlin, B and Munich, M.
Sample preparation for scanning electron microscopy (SEM) followed a standard protocol (Janofske, 2000), applied in several previous studies (Gottschling et al., 2012;Kretschmann et al., 2018aKretschmann et al., , 2020. The cells were fixed overnight using 2.5% glutar(di)aldehyde. The cells were filtered onto an Omnipore TM -membrane filter (5 µm; Merck; Darmstadt, Germany), which was placed in a Swinnex ® filter holder (Merck), and dehydrated with a graded acetone series. Critical point drying was followed by sputtercoating (BAL-TEC SCD 050 sputter coater; Schalksmühle, Deutschland) with platinum on an aluminium stub coated with Planocarbon (Plano; Wetzlar, Germany). Subsequent imaging was performed with the SEM LEO 438VP (LEO; Cambridge, UK) at 15 kV. Image adjustment and arrangement of LM and SEM pictures were done in Photoshop and InDesign (Adobe; Munich, Germany).

Molecular phylogenetics
A systematically representative set of peridiniacean accessions was compiled, considering the information provided by dinophyte reference trees such as presented in Gottschling et al. (2020). The sample was enriched by all those sequences deposited in GenBank, which showed ultimately close relationships to the sequence data gained in the present study, as inferred from BLAST searches (Altschul et al., 1990). Full voucher information is provided in Table S1, including information of the outgroup comprising dinophytes of Heterocapsaceae and Protoperidiniaceae. The ITS sequences were also inspected for possible compensatory base changes that occur when substitutions take place in pairing regions of the molecule folded into a secondary structure (Gottschling & Plötner, 2004;Kremp et al., 2014;Thornhill & Lord, 2010).
For alignment constitution, separate matrices of the rRNA operon (i.e., SSU, ITS, LSU) were constructed, aligned using 'MAFFT' v6.502a (Katoh & Standley, 2013) and concatenated afterwards. The aligned matrices are available as *.nex files upon request. Phylogenetic analyses were carried out using maximum likelihood (ML) and the Bayes theorem, as described in detail previously  using the resources available from the CIPRES Science Gateway (Miller et al., 2010). The MPI version of 'RAxML' v8.2.4 (Stamatakis, 2014, freely available at http:// www. exeli xislab. org/) was applied using the GTR + Γ substitution model under the CAT approximation. We determined the best-scoring ML tree and performed 1,000 non-parametric bootstrap replicates (rapid analysis) in a single step. The Bayesian analysis was performed using 'MrBayes' v3.2.6 (Ronquist et al., 2012, freely available at http:// mrbay es. sourc eforge. net/ downl oad. php) under the GTR + Γ substitution model and the randomaddition-sequence method with 10 replicates. We ran two independent analyses of four chains (one cold and three heated) with 20,000,000 generations, sampled every 1,000th cycle, with an appropriate burn-in (10%) as inferred from the evaluation of the trace files using Tracer v1.5 (http:// tree. bio. ed. ac. uk/ softw are/ tracer/). Statistical support values (LBS: ML bootstrap support) were drawn on the resulting, bestscoring tree.

Morphology of Peridinium. volzii
All investigated strains exhibited cells that were either flagellated ( Fig. 1a, i) or coccoid ( Fig. 1k-m), with the motile cells being predominant. Motile cells swam homogenously within the medium or occasionally, they accumulated towards the light source. Each cell, except for necrotic cells, contained numerous gold-through olive-brownish chloroplasts and an orange to red accumulation body (Fig. 1i). An eyespot was absent. Necrotic cells included plastids that appeared grey through silver.
All strains were indistinguishable from each other in gross morphology, and motile cells were continuously covered by a theca built of cellulosic plates (as inferred from astra blue staining). The thecate cells had an ovoid shape and were slightly compressed dorso-ventrally. The hypotheca was generally smaller than the epitheca. Average cell sizes of motile cells within the ten strains under detailed study ranged from 41 to 51 µm in length and from 36 to 48 µm in width (Table 1). However, cell sizes of particular strains were significantly different from each other (Fig. 2a, b). All strains presented the thecal formula 4′, 3a, 7′′, 5c, 5s, 5′′′, 2′′′′, with few exceptions (see below). The surface of the plates was strongly reticulate, and each major plate of the epi-and hypothecae contained one or two pores.
The epithecal plate pattern was largely symmetrical, although the apical plate 2′ was always smaller than the apical plate 4′ (Fig. 1a, c-d, f). The apical plate 1′ appeared small, was rhomboid and ca 12-14 μm tall. Two significantly different suture lengths between the plates 1′ and 4′ were distinguishable across the strains (Fig. 2c), being shorter ( Fig. 1a) or longer ( Fig. 1c-d). Both apical plates 2′ and 4′ appeared hexagonal in shape and plate 3′ pentagonal. The apical plate 3′ is separated from plate 1′ by plates 2′ and 4′ and was either as wide as tall ( Fig. 1f) or slightly wider than tall (Fig. 1c). The intercalary plates 1a and 3a were pentagonal, had the same size and were smaller than the elongated plate 2a (Fig. 1c, f).
Deviations from the archetypical plate pattern could be stated (Tables 1, 2 and 3). Usually tetragonal, precingular plate 4′′ was occasionally pentagonal (Fig. 3a, b) and abutted plate 3a, from which it was otherwise separated. This conformation affected the position and size of the intercalary plate 2a, which appeared smaller and less elongated. In few individual cells of all strain, the intercalary plate 2a (Fig. 3c), the apical plate 3′ (Fig. 3m) or the postcingular plate 3′′′ (Fig. 3i) were split, and also unusual fusions of plates could be observed (Table 3; Fig. 3d, f).
The cingulum was median or sub-median and descending. The displacement was significantly different between the strains (Table 1), with either ca 1 cingulum width (Fig. 1a) or ca 1.5-2 cingulum widths (Fig. 1d). Notably, this distinction correlated with the two different suture lengths between apical plates 1′ and 4′ (Table 1). The sulcus appeared narrow and extended into the epitheca for approximately a cingulum width (Fig. 1a, d). Two flagella originated from the junction between sulcus and cingulum and were visible in LM.
Cell division was by eleutheroschisis. The resulting empty thecae were collected in high quantity at the bottom of the cultivation plate. The theca opened along the edge of the cingulum that epitheca and hypotheca were separated. The predominant opening of the epitheca was exhibited at the apical side between the three intercalary plates and the plates 2′ and 4′ (Fig. 1f). The cell opening continued at the sutures between the precingular plates 5′′ and 6′′ as well as 2′′ and 3′′, resulting in the release of the dorsal part of the epitheca. Subsequently, this thecal split led to the release of the ventral epitheca, as well as the separation of epitheca and hypotheca (Fig. 1m). This specific dehiscence line was present in all strains with a frequency of 66-82% (Table 1). Deviations from the regular opening line predominantly occurred at the precingular plates in all strains investigated. Sutures between thecal plates varied from thin lines through wide bands, the latter exhibiting cross striations (Figs. 1f,3b,e,g).
Coccoid cells varied in their average size from 39 to 50 µm in length and from 32 to 45 µm in width (Table 1) and were either thecate or athecate. Their shapes were variable ranging from spherical through obovoid. Some coccoid cells exhibited a slight dorso-ventrally flattening (Fig. 1l). Within these cells, brown granules of varying sizes were observed in addition to the plate-like chloroplasts. A single coccoid cell was observed to develop intrathecately (Fig. 1m), and the ecdysing cell was released by thecal opening along the edge of the cingulum.

Peridinium volzii is morphologically distinct from Peridinium willei
Correct and consistent species identification is the necessary prerequisite for investigations of ecosystem dynamics and conservation strategies. However, delimitation and determination of unicellular organisms at the species level are challenging due to complex identification procedures. Exemplarily, there has been some uncertainty whether the names, P. volzii and P. willei represent a single species (Playfair, 1920;Popovský & Pfiester, 1990) or in fact two distinct species. The integrative approach followed in this study, combining a molecular phylogenetic tree and detailed morphological investigations, unequivocally shows the distinctiveness of the two entities and their sister group relationship. The differences in the ecological   preferences in the two species provide further evidence for their separation (Höll, 1928;Olrek, 1992). Bachmann (1911) was probably the first to accept two distinct species. In addition to P. willei, also Lindemann (1918) recognised a second species, initially his P. guestrowiense (Lindemann, 1916) and since Lindemann (1924) at the latest the older P. volzii. The most striking diagnostic trait is the size of the first apical plate, which is large in the lineage including accessions determined as P. willei and which is small in the lineage including those determined as P. volzii. Only strain PWCL1 has been determined as P. willei (Logares et al., 2007) but is placed on the branch of P. volzii. However, no morphology is known from this strain, and it is likely a misidentification. Ling et al. (1989) reported from intermediates between P. volzii and P. willei but in the strains studied here in detail, the first apical plate was always small, without exception. From a phylogenetic perspective, the character state is presumably apomorphic for P. volzii (versus bigger in other Peridiniaceae including P. willei as ancestral state).

Peridinium volzii is morphologically variable
Despite the uniqueness of P. volzii due to the small first apical plate, the morphological plasticity within the evolutionary lineage appears great. However, all this plasticity (shown within monoclonal strains and partly formalised in numerous varieties and forms) cannot be associated with other traits such as environmental conditions, geographic occurrence or genetic constitution that evolutionary adaptation appears of minor importance in P. volzii. A considerable number of infraspecific taxa have been described (Table 2) and based on the microscopic study of monoclonal strains we confirm the existence of traits such as additional sutures (associated with, e.g., P. volzii forma complexum; Fig. 3m) and fusion of plates (associated with, e.g., P. volzii var. simplex; Fig. 3l). Other traits such as the shifted suture (of Peridinium guestrowiense var. betacollineatum Er.Lindem.) having importance in species such as P. cinctum (Izquierdo López et al., 2018;Romeikat et al., 2019) have not been found during the course of the present study. It is known today that varying widths of sutures, where the intercalary space becomes distinctively striate, are ontogenetically disposed (growth bands: Netzel, 1982). This trait (Figs. 1a, 3b, e, i, k) is associated with P. guestrowiense forma late-intercalatum and P. guestrowiense forma latissime-intercalatum Er.Lindem. and is the only one that occurs in combination with other morphological deviations.
The diagnostic trait of P. volzii var. australe G.S.West is the antapical plates of unequal size (versus equal size in P. volzii var. volzii), and this character state is frequent (> 70%) in the material under investigation (Fig. 3i). Similarly, the diagnostic trait of P. volzii var. cinctiforme M.Lefèvre is the apical plates 2′ and 4′ of unequal size (versus rather equal size in P. volzii var. volzii), and this character state (distinctive in species such as P. cinctum and Peridinium raciborskii Wołosz.) is also not rare (25-50%) in the material studied here (Fig. 3k). We even confirm the existence of Peridinium guestrowiense subvar. originale Er.Lindem. having six postcingular plates (versus five such plates usually present in peridinialean dinophytes; Fig. 3i). The reduction from six to five postcingular plates in Peridiniales may result from a fusion of the third and fourth postcingular plate  present today in dinophytes such as the Amphidomataceae (Tillmann et al., 2009). Already Lindemann (1920: 122) pointed out the possibility that phylogenetically ancestral character states occasionally recur (known as atavism), which is corroborated by our observations of individual such cells in monoclonal strains. However, all such variation occurs in individual cells and cannot be associated with other traits such as a genetic constitution or geographic occurrence. Therefore, the recognition of corresponding taxa at infraspecific ranks does not appear justified. In addition to the historically described intraspecific variabilities, we detected five cases of yet undiscovered deviation in P. volzii. The unusual epithecal configuration with plate 4″ abutting plate 3a (Fig. 3a-b) has to the best of knowledge never been noted in the literature before. More-over, alternate configurations between, for example, plates 3′ and 4″ (resp. 1a and 3a; Fig. 1c, f) within monoclonal strains are reminiscent of, for example, Parvodinium mixtum Wołosz. ex Kretschmann, Zerdoner, Owsianny & Gottschling (Kretschmann et al., 2018a). An additional suture and the split of plate 2a (Fig. 3c) have also not been reported from P. volzii so far. Fusion of plates regularly though rarely appears in the thecal pattern of various dinophytes, and corresponding observations in the present study do not come as a surprise. Usually, the fusion of plates takes place within plate series such as the intercalary (Fig. 3f), precingular (Fig. 3d) and antapical series (associated with P. volzii var. simplex; Fig. 3l). However, the fusion of plates associated with different plate series is very rare, and we found fused plates 1′″ and 1″″ (Fig. 3e) once only. Anyhow, also the newly discovered variability occurs in individual cells associated with all three ribotypes present in this study and does not make the recognition of corresponding taxa at an infraspecific rank seem justified either.

Correlations between ribotypes and morphotypes
The goal and importance of Linnaean binomials in nomenclature is to avoid different scientific names for the same biological entity (Minelli, 2019). Phenotypic variants occur at several interindividual levels, such as between reproductively isolated units (i.e., species), between populations (e.g., local communities) or within populations (e.g., alleles). There is no ultimate consensus about the usefulness of infraspecific taxonomic ranks (such as subspecies, varieties and forms) in order to provide such variants with scientific names. In practice, different taxonomists (at different points in time) treat(ed) phenotypic variation differently (Knapp et al., 2004;Mayr, 1982a, b) and created considerable complexity (Hamilton & Reichard, 1992), also in dinophytes (Izquierdo López et al., 2018;Moestrup & Calado, 2018;Romeikat et al., 2019). Infraspecific classifications should only be made with a certain reserve, as categorisation of minor morphological variability conflicts with the practical goal of stability in taxonomy (Boenigk et al., 2012;Turland, 2019;Warburton, 1967). As criteria for the acceptance of taxonomic entities, the present case studies possible correlations between genetic disposition (i.e., certain ribotypes), morphology (i.e., diagnostic traits) and/or geographic occurrence. There is considerable variation (previously already observed based on enzyme electrophoresis and flow cytometric determination of nuclear DNA quantities: Hayhome et al., 1987), and the basic question is whether the evolutionary lineage of P. volzii as displayed in the molecular tree represents a single or several species. In the previous part of the discussion, we have already argued that most of the infraspecific taxa described so far are not suitable as taxonomic subsets of P. volzii.
No morphological and ITS sequence information is available for strains PWCL1 and NIES501, that they are not worthy to be considered further here. However, there are significant correlations to the state between certain ITS ribotypes and the three traits cell size, cingulum displacement and length of the suture between plates 1′ and 4′, respectively. Data overlap is very low, which underlines the distinctiveness of certain character states. Ribotype I is associated with the smallest cell size, the single-width displacement of the cingulum and a short suture. Morphometrics are similar for ribotype II, which shows, however, greater cell size. Ribotype III is distinct by the combination of the greatest cell size, a displacement of the cingulum beyond 1.5 fold the width and a long suture. Notably, the similarity between cells assigned to ribotypes I and II is greater than any of them to ribotype III, but ribotypes II and III are closely related in the DNA tree. However, no unequivocal or operable correlations between certain ribotypes, diagnostic traits, ecological preferences and/or geographic occurrence can be drawn. Cryptic speciation (Fenchel, 2005) can certainly not be excluded as it has been shown for other dinophytes (Lajeunesse et al., 2012;Luo et al., 2017;Montresor et al., 2003;Wang et al., 2019), but a single, variable species P. volzii is more likely at this moment in time, as it was worked out also for P. cinctum of the same phylogenetic group (Izquierdo López et al., 2018).

Life-history
Life-history and metagenesis of Peridiniaceae is complex and not fully understood at present. Sexual reproduction of P. volzii is isogamous and heterothallic and includes motile and coccoid cells during the diploid stage (Hayhome et al., 1987;Kita & Fukuyo, 1995;Pfiester & Skvarla, 1979). In our monoclonal strains, we have never observed any fusion of cells (with subsequent karyogamy) or four-cell aggregations (i.e., indication for meiosis). This confirms observations of Pfiester and Skvarla (1979) that sexuality does not take place within monoclonal strains of P. volzii (other Peridiniaceae are homothallic : Pfiester, 1975: Pfiester, , 1976: Pfiester, , 1977. Three size classes are stated for the motile cells of P. volzii (Pfiester & Skvarla, 1979) and other Peridiniaceae (Pfiester, 1975(Pfiester, , 1976(Pfiester, , 1977, namely small cells < 30 µm (gametes, certainly haploid), large cells ≥ 60 µm (zygotes, certainly diploid) and mid-sized cells, which are in the range described in the present study. The ploidy level of these motile cells is unknown at present, but it is not improbable that they are haploid as the gametes. Subsequently, the various coccoid cells observed in the present study are all most likely also part of the haploid stage (as sexuality does not take place within monoclonal strains). The here reported, intrathecately formed coccoid cells are morphologically not differentiated from the coccoid cells of the diploid stage. Pfiester and Skvarla (1979) do not provide exact measurements for those cells they term hypnozygotes, but they are probably bigger than the otherwise indistinguishable cells in our study. More research is necessary to entangle all these different cell types and place them properly in the metagenesis of Peridiniaceae. Little attention has been put in the past on the dissociation of thecal elements during mitotic cell division. For the closely related P. willei, Moestrup and Calado (2018) state that the epitheca segregates into pieces, one of which is said to be composed of the plates 2′, 1a-3a, 3′′′-5′′′. In P. volzii, the lid is composed of the plates 3′, 1a-3a, 3′′-5′′, and whether the break of the epitheca is truly different between the two species, or whether the information provided by Moestrup and Calado (2018) contain typos, remains to be determined (no reference provided). The opening of the cells during regular mitosis has not been observed in P. limbatum, but the dehiscence line on the epitheca during the formation of the coccoid cell (Evitt & Wall, 1968;Wall & Dale, 1968) is similar to P. volzii. However, the hypotheca remains entire in the latter species, whereas it also splits in P. limbatum during the release of the coccoid cell. Notably, the course of the archaeopyle of P. limbatum is identical to the observations made for the epitheca of P. volzii. Alternative sheds of the theca, such as the dorso-ventral split of the epitheca in †Lingulodinium D.Wall  or the hypotheca opening in Peridiniopsidaceae (Kretschmann et al., 2018a), may indicate that this trait is phylogenetically informative, but not sufficiently investigated so far.

Conclusion
The integrative approach presented here studying both morphological and molecular data is important to clarify the boundaries between species in the microbial domain. Based on monoclonal strains established from recently collected material we confirmed a great morphological and genetic variability in P. volzii. However, the decision whether the three ITS ribotypes correspond to one, two or three species is impossible to make the present data given. Morphological data of strains NIES501 from Japan and PWCL1 from Wisconsin, or published molecular data of strain SCCAP K-1155 from Italy (Hansen & Flaim, 2007), would probably render a clearer structure to delimit species and populations within dinophytes such as P. volzii.
In the light of the present new data, it is worthy to reconsider the taxonomic status of scientific names associated with P. volzii, but also of the species itself described Fig. 3 Historically described phenotypical variations and yet undiscovered deviations in the plate pattern of P. volzii. b-c, f, i, l-m Light microscopy, a, d-e, g-h, k scanning electron microscopy. a-f Newly identified deviations a-b plate 4′′ pentagonal in strains a GeoM*793; b GeoM*788; c plate 2a split (strain GeoK*024); d plates 2′′ and 3′′ fused (strain GeoM*866); e plates 1′′′ and 1′′′′ fused (strain GeoM*788); f plates 1a and 3′ fused (strain GeoM*788). g-m Historic infraspecific taxa; g P. guestrowiense forma lineatum (strain GeoM*866); h P. guestrowiense forma compressum (strain GeoM*866); i P. guestrowiense subvar. originale (strain GeoK*024); k P. volzii var. cinctiforme (strain GeoM*793); l P. volzii var. simplex (strain GeoM*789); m P. volzii forma complexum (strain GeoM*793). Abbreviations: n′: apical plate, n′′: precingular plate, n′′′: postcingular plate, n′′′′: antapical plate, na: anterior intercalary plate, nC: cingular plate, split or fused plates are indicated by asterisks. Scale bar = 10 µm. U A = 15 kV ◂ from Singapore. In the protologue, Lemmermann (1905) highlighted the trait of an antapical spine, which appears abundant in material from Singapore, Indonesia and Australia (Lindemann, 1931;Ling et al., 1989;Playfair, 1920;West, 1909). Such spine is absent from all illustrations we are aware of from Central Europe (Baumeister, 1976;Entz, 1927;Lefèvre, 1932;Lewis & Dodge, 2011;Lindemann, 1916Lindemann, , 1920Olrik, 1992;Кpaxмaльный, 2011), and it is also absent without exception from the material studied here. Cells from Australasia are also remarkably smaller than those reported from Europe. Doubts are allowed that true P. volzii occurs in Central Europe at all, and future studies must clarify the distinctiveness from German P. guestrowiense as next-younger scientific name. Similarly, P. volzii var. maeandricum has a very distinctive ornamentation, of which it is said being consistent within strains (Hansen & Flaim, 2007). This could be indicative for the species status of Peridinium maeandricum (Lauterborn) V.Brehm rather than for the variety level, and the characteristic ornamentation is consistently absent from the material studied here. More research is necessary to entangle the complex interactions between genetic disposition, morphological variation and geographic occurrence of unicellular organisms such as P. volzii and its relatives.
Acknowledgements We sincerely thank Eva Facher (Munich) for assistance in the lab.
Funding Open Access funding enabled and organized by Projekt DEAL. This work was supported by the Deutsche Forschungsgemeinschaft (grant GO1459 10-1).

Data availability
The datasets generated during and analysed during the current study will be available in the GenBank repository after acceptance of the manuscript. The corresponding Accession numbers are listed in the Voucher List (Table S1). Comprehensive morphological data are available from the corresponding author on reasonable request.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. Fig. 4 A molecular tree of 51 systematically representative Peridiniaceae, including all 28 accessions assignable to P. volzii. Maximum Likelihood tree (-ln = 22,017.62), as inferred from a rRNA nucleotide alignment (1,129 parsimony-informative sites) and with strain number information. Numbers on branches are ML bootstrap (above) and Bayesian support values (below) for the clusters (asterisks indicate maximal support values, values under 50 and 0.90, respectively, are not shown). Clades are indicated (abbreviations: HET, Heterocapsaceae; PPE, Protoperidiniaceae) ◂