High diversity in Keratella cochlearis (Rotifera, Monogononta): morphological and genetic evidence

Rotifers are ubiquitous freshwater animals for which many complexes of cryptic species (i.e. distinct species that are morphologically difficult to distinguish) are described. Keratella cochlearis occurs globally and shows a wide phenotypic diversity indicating the potential presence of a species complex. We sampled lakes of the Trentino-South Tyrol region (Italy) and investigated mitochondrial genetic diversity in K. cochlearis in relation to detailed lorica measurements. We sequenced the mitochondrial cytochrome c oxidase subunit I and used the generalised mixed Yule coalescent approach, Poisson tree process model and automatic barcode gap discovery to delimit mitochondrial groups, associated with putative evolutionary significant units (ESUs). Based on 248 sequences, eight putative ESUs were indicated that could only partially be delimited by lorica morphology. Specifically, several morphological characteristics (i.e. spinelets, bended median ridge, and posterior spine) were found in specimens of different putative ESUs, and thus, these characters seem to be of poor discriminatory value. Furthermore, different putative ESUs of K. cochlearis were found in the same lake. We conclude that the high mitochondrial genetic diversity may be linked to tolerance of K. cochlearis to varying environmental conditions.


Introduction
Biodiversity is currently under threat, and our perception of species loss is highly dependent on accurate estimates of species richness. However, estimates of species richness are often impaired by the occurrence of cryptic species (i.e. species that are impossible or difficult to distinguish based on their morphology) in diverse groups such as protists (Foissner, 2006), ants (Fournier et al., 2012), harvestmen (Arthofer et al., 2013), and rotifers (Gómez & Snell, 1996). Understanding how and why species occur is one of the fundamental aspects in ecology (Gaston, 2000).
Keratella cochlearis Gosse, 1851 can be found in most freshwater lakes and ponds all over the world (Green, 1987). In fact, the whole genus Keratella is considered eurytopic and cosmopolitan (Segers & De Smet, 2008), and this makes the genus a good candidate for investigating the occurrence of cryptic species. Lauterborn (1900) described several morphotypes in K. cochlearis, and his detailed descriptions and drawings were the basis for following taxonomic work (e.g. Ahlstrom, 1943;Ruttner-Kolisko, 1974;Koste, 1978). The morphotypes described by Lauterborn (1900) encompass three series (macracantha-typica-tecta, hispida, and irregularis) and the group of robusta. These morphological varieties of K. cochlearis are different with respect to lorica length (LL), spine length, presence of spinelets on the lorica, and the course of the median ridge. Here, we give an overview of the Lauterborn (1900) series and a German to English translation of Lauterborn's (1900) descriptions. In the macracantha-typica-tecta series (Lauterborn's 1900, Figs. 1-10), the posterior spine is as long as the lorica or even longer, and the basis of the spine is so wide that it is difficult to decide where the spine begins and the lorica ends. The areolation is present on half of the spine, and only the distal part is smooth and pointed. In lateral view, the spine points to left or right, and this is according to Lauterborn (1900) not an important feature. Along the series, the reduction of the posterior spine is notable until it disappears completely. Lauterborn (1900) concluded that it is impossible to draw a line between the different morphotypes of the macracantha-typica-tecta series that only differ in size and posterior spine length (PSL). The morphotypes of the hispida (Lauterborn's 1900, Figs. 11-14) and irregularis series (Lauterborn's 1900, Figs. 15-20) show different morphological elements with respect to the macracantha-typica-tecta series, and size differences are not important. In the hispida series, small spines (called ''Pusteln'' after Lauterborn, 1900 and''spinelets'' after Ahlstrom, 1943) are present and can be so dense that the areolation and the borders of the plates become invisible. The morphotypes of the hispida series can be considered the forma punctata of the tecta series. Only for Lauterborn's (1900), Fig. 11, closely related to macracantha, and for Lauterborn's (1900), Fig. 27, closely related to tecta, the name forma punctata is given. In the irregularis series, the ridge is bended to the left in dorsal view, and a displacement of the facets is visible that leads to pointed bumps (called ''Höcker'' by Lauterborn, 1900) on the facets and an additional facet (called facet X by Lauterborn, 1900). In addition, the basal margin is divided into small posterior carinal facets. Similar to the hispida series, the lorica has small pointed spinelets on the intersection of the areolation. The robusta group (Lauterborn's 1900, Figs. 21-23) is not a series because no direction of morphological variations can be distinguished. Characteristic for this group is the wide base of the posterior spine that is the elongation of the ventral part of the lorica, the hooked form of the anterior spines, and the slightly bended median ridge.
Considering the wide morphological variability of Keratella morphotypes, Lauterborn (1900) already hypothesised a subspecies status of some morphotypes. In fact, Ahlstrom (1943) and Eloranta (1982) erected the series irregularis and hispida to separate species. However, Hofmann (1983), who did not recognise transitional forms between the morphotypes cochlearis, irregularis, and tecta as described by Lauterborn (1900), questioned the validity of the Lauterborn cycles. Especially, the presence and length of the posterior spine seems to be a morphological character whose suitability for discriminating species is questionable. In eutrophic habitats, K. cochlearis tends to be smaller and has smaller posterior spines than in oligotrophic habitats (Green, 2007). Furthermore, LL and PSL are longer with decreasing water temperature (e.g. Green, 1981Green, , 2005Bielanska-Grajner, 1995) and in the presence of predators (Conde-Porcuna et al., 1993;Green, 2005). Water conditioned with predators (i.e. Asplanchna spp., cyclopoid copepods) can induce spine formation in offspring of tecta (Stemberger & Gilbert, 1984). Derry et al. (2003) found a high mitochondrial genetic difference [4.4% cytochrome c oxidase subunit I (COI) sequence divergence] between spined and spineless individuals of K. cochlearis and hypothesised the presence of cryptic diversity within these morphotypes. Furthermore, the various morphotypes of K. cochlearis show different tolerances to temperature (Berziņ š & Pejler, 1989a), oxygen content (Berziņ š & Pejler, 1989b), trophic state, and conductivity (Berziņ š & Pejler, 1989c). The wide tolerances to environmental conditions could also indicate that K. cochlearis is a cryptic species complex composed of species with narrower ecological preferences than when taken as a complex.
Here, we identified mitochondrial DNA (mtDNA) groups and compared their lorica morphology in a complementary approach as recommended by Schlick-Steiner et al. (2006), Fontaneto et al. (2015) and Mills et al. (2016). Combining genetic information with other species-bound aspects such as species morphology and ecology or biochemistry of species habitat can result in a more robust species delimitation than when using genetic information alone (Schlick-Steiner et al., 2006;Fontaneto et al., 2015;Mills et al., 2016). We hypothesised that K. cochlearis is a complex of putative evolutionary significant units (ESUs) and that it is possible to delimit ESUs based on lorica measurements. In fact, in B. plicatilis some clusters of cryptic species [B. plicatilis (sensu stricto) L., B. rotundiformis SS, B. rotundiformis SM] can be distinguished based on body length differences (Ciros-Pérez et al., 2001b). Closely related species might have similar niches according to the phylogenetic niche conservatism theory (e.g. Wiens & Graham, 2005;Wiens et al., 2010), and this may lead to competitive exclusion (Violle et al., 2011). Thus, ESUs with their close phylogenetic relationship might be especially prone to competitive exclusion; however, co-occurrence of rotifer cryptic species has been reported (Obertegger et al., 2014). Thus, we also investigated temporal co-existence of putative ESUs of K. cochlearis and hypothesised little co-occurrence.

Sampling
From March to November 2014, six lakes in the Trentino-South Tyrol (Italy) region were sampled monthly. These lakes (called further the ''core lakes'') cover a wide range of environmental parameters (Table 1). In addition, we also sampled 11 additional lakes from Trentino-South Tyrol in the years 2010, 2013, and 2015 during summer and winter to cover a larger geographical area and altitudinal range (Table 1; Fig. 1). Environmental parameters were based on published data (IASMA, 1996(IASMA, -2000 and own analyses (Table 1). At the deepest site of each lake, plankton samples were collected with a 20 lm (Apstein) or 50 lm (Wisconsin) plankton net depending on lake depth. Both mesh sizes were small enough to effectively collect specimens of K. cochlearis (length [74 lm, width [60 lm;Lauterborn, 1900;Koste, 1978).

Measurements of specimens and morphological observations
For the core lakes and Lake Caldonazzo (July sample), single specimens of K. cochlearis were isolated under a stereomicroscope and photographed (Leica DC 300F camera, Leica IM1000 software) in dorsal and lateral view under a compound microscope. The following measurements were taken: PSL, LL excluding anterior and posterior spines, total LL (TLL) including all appendages, lorica width (LW) at its widest part, LW at the mouth opening region (''head width'', HW), anterolateral dorsal spine length (ALS), anterointermediate dorsal spine length (AIS), anteromedian dorsal spine length (AMS), and posterior spine angle (PSA, Fig. 2). For the measured specimens, we also observed the main characteristics of the dorsal plate, important to discriminate morphotypes. Measured specimens were subject to DNA extraction and sequencing. However, we could not obtain sequences for all measured specimens.

Phylogenetic reconstruction
We constructed the phylogenetic tree using a maximum likelihood (ML) and Bayesian inference (BI) The superscript c indicates the core lakes ordered by altitude (alti, m above sea level): area (910,000 m 2 ), depth (m), total phosphorus (TP, lg l -1 ) at spring overturn, nitrate (NO 3 , lg l -1 ), reactive silica (Si, mg l -1 ), sulphate (SO 4 , mg l -1 ), chloride (Cl, mg l -1 ), conductivity (cond, lS cm -1 ), mean summer surface temperature (temp), and trophic state (eu eutrophic, meso mesotrophic, oligo oligotrophic)  (Keane et al., 2006). The selected model was implemented into PhyML 3.0 (Guindon & Gascuel, 2003) to perform ML reconstruction using the approximate likelihood ratio test to evaluate node support. For BI, we used BEAST v1.8.0 (Drummond et al., 2012) with the following settings: uncorrelated lognormal relaxed clock (mean molecular clock rate set as normal), HKY ? I ? G substitution model, and the birth-death model. The posterior probability distribution was estimated with Markov chain Monte Carlo (MCMC) sampling, which was run for 100 million generations, sampling every 10,000th generation. We used Tracer v1.5 (Rambaut et al., 2014) to investigate for convergence and the correctness of the MCMC model and to determine the burnin. We used TreeAnnotator v1.7.5 to summarise trees and discard the first 2,000 trees as burn-in. As outgroup sequences, we used B. urceolaris (Genbank accession number EU499787), B. rotundiformis (JX239163), and B. plicatilis (JX293050), all belonging to the same family (i.e. Brachionidae) as Keratella.

Inference of mtDNA groups
We inferred mtDNA groups within K. cochlearis with the generalised mixed Yule coalescent (GMYC) approach (Fujisawa & Barraclough, 2013), the Poisson tree process model (PTP; Zhang et al., 2013), and the automatic barcode gap discovery (ABGD; Puillandre et al., 2012) and compared the results. For all methods, the outgroup was excluded prior to the analyses. We took the results of the GMYC approach as our baseline results because previously rotifer diversity was investigated by it for different species (Obertegger et al., , 2014Leasi et al., 2013;Malekzadeh-Viayeh et al., 2014). The GMYC approach is based on branching rates along an ultrametric tree (here from BEAST) to distinguish between species-level (Yule, slower) and populationlevel (coalescent, faster) branching rates. This model identifies GMYC ESUs. For the GMYC approach, we used R 3.0.2 (R Core Team, 2012), library splits (Ezard et al., 2009). The PTP model (http://species.hits.org) uses a phylogenetic tree as input (here the ML tree produced in PhyML 3.0.) and applies coalescent theory to distinguish between population-level and species-level processes. Similarly to GMYC, PTP assumes that there are less intraspecific substitutions than interspecific substitutions because they have less time to accumulate. This method does not require an ultrametric tree and has been shown to match other methods of species delimitation in rotifers (Tang et al., 2014) and copepods (Blanco-Bercial et al., 2014). Two types of PTP were used: ML (PTP-ML) approach and Bayesian approach (PTP-BA). The ABGD (http:// www.abi.snv.jussieu.fr/public/abgd/abgdweb.html) deliminates species without any a priori assumptions. It detects the gaps in the distribution of genetic pairwise distances. This method has been successfully used to delimit species of the meiofauna (Tang et al., 2012;Leasi et al., 2013). Here, all aligned K. cochlearis sequences were used for ABGD. We based our phylogenetic reconstructions and inference of mtDNA groups on a single mitochondrial gene (COI), and this may gave a biased estimate on genetic diversity. A higher evolutionary rate of COI with respect to other nuclear markers (Tang et al., 2012), mitochondrial introgression (reported for B. calyciflorus by Papakostas et al., 2016 but not for E. senta by Schröder & Walsh, 2010), and/or unresolved ancestral polymorphism (Funk & Omland, 2003)  could bias our inference on species diversity.
Recently, it has also been shown that the methods we used give biased results in species poor datasets (Dellicour & Flot, 2015). Thus, considering this uncertainty, our statements are about putative ESUs based on the inference of mtDNA groups.
Statistical analysis of measurements in relation to putative ESUs Green (1981Green ( , 1987 reports a positive correlation between LL and PSL in K. cochlearis from various lakes of the Auvergne region in France. To assess the general validity of this correlation, we considered only those specimens that were measured and for which we obtained COI sequences. We divided specimens into putative ESUs and investigated the sign and significance of the correlation (Pearson correlation coefficient; r P ) between LL and PSL. We performed a univariate statistical analysis and a multivariate ordination method to investigate if putative ESUs could be distinguished based on morphology. As univariate statistical analysis, we used a oneway ANOVA and post hoc Tukey multiple comparisons. We performed generalised least squares modelling to allow for dependence of measurements of ESUs coming from the same lake and checked homogeneity of residuals graphically. As multivariate ordination method, we performed non-metric multidimensional scaling (NMDS). In NMDS, Bray-Curtis distance matrix was used on centred and standardised measurement data. In NMDS, the goodness of fit was investigated by the Shepard plot that shows the relationship between the inter-object distances in NMDS and Bray-Curtis dissimilarity. The residuals of this relationship were used to calculate Kruskal's stress (S); S values \0.2 are considered statistically meaningful (Quinn & Keuogh, 2002). We, furthermore, performed a linear discriminant analysis (LDA) to investigate the discriminatory power of lorica morphology to separate ESUs. We tested for homogeneity of within-ESU covariance matrices.
We also investigated the correlation between phylogenetic and morphological diversity. Phylogenetic diversity was calculated as distance matrix based on the ultrametric tree, and morphological diversity as a distance matrix based on mean morphological values of ESUs. The correlation between both distance matrices was investigated by a Mantel test.
The ABGD and the PTP-ML grouped the same haplotypes in the same ESUs as GMYC (Fig. 3). However, PTP-BA, split GMYC ESU 3 into three and ESU 6 into five units (Fig. 3).

GMYC ESUs occurrence in lakes
GMYC ESUs 3 and 7 were found in seven lakes, ESU 8 in six, ESU 5 in five, and ESU 4 and 1 were found only in two and ESU 2 only in one lake ( Fig. 3; Tables s2, s3 supplementary material). Considering temporal co-existence of GMYC ESUs in the core lakes, no clear pattern emerged (Table s3 supplementary material). Generally, GMYC ESUs co-occurred, except for ESU 2 that was found only once in Lake Radl, despite monthly sampling during summer 2013. ESUs 3 and 7 co-occurred most often in different lakes. ESU 3 was almost always present throughout the sampling period in Lakes Kaltern and Terlago (Table s3 supplementary material); similarly, ESU 5 in Lake Glittner and ESU 6 in Lake Tovel were present throughout the sampling period (Table s3 supplementary material).

Morphology
We obtained lorica measurements from 138 individuals of K. cochlearis that could also be attributed to GMYC ESUs based on their COI sequence (Table 4;  Table s3 supplementary material). For ESUs 1 and 2, no measurements were obtained, and for ESU 7, only one specimen was measured (Table 4). All specimens of ESU 4 and three specimens of ESU 6 did not have a spine, while the other measured specimens had a spine of varying length (Table 4).
The correlation between LL and PSL was different when based on all specimens (r P = 0.68; P \ 0.001) compared to splitting it into GMYC ESUs: for ESUs 3 and 6, the correlation was higher (r P = 0.76 and 0.77, respectively; P \ 0.001) than the overall one, and no correlation was found for ESU 4 (spineless specimens), ESU 5 (r P = 0.13; P = 0.41), and ESU 8 (r P = 0.76; P = 0.13; Fig. 4).
We tested for significant differences in LL, PSL, and PSA between GMYC ESUs by ANOVA and following post hoc multiple comparisons tests by mixed modelling. LL and PSA were different between four ESUs, and PSL differed between three ESUs (Table 5). Based on all three measurements, ESU 8 was different from ESUs 3 and 5 (Table 5).
In NMDS with all measurements (S = 0.13), a gradient from specimens of ESU 5 to specimens of ESU 4 and spineless specimens of ESU 6 was evident.
To get a clearer picture on the relationships between ESUs with spines, we excluded ESU 4 and the three spineless specimens from ESU 6 from the NMDS analysis. In this NMDS with measurements of spined individuals (S = 0.17), specimens of ESUs 5 and 8 formed distinct clusters while specimens from ESUs 3 and 6 were mixed (Fig. 5). In the LDA based on PSL, LL, and PSA, the percent correct assignment of ESUs varied (ESU 3: 69%, ESU 5: 83%, ESU 6: 53%, ESU 8: 50%).
We noted the presence of spinelets (Fig. 6), additional facets, and bending of the ridge (Fig. 4) in some specimens and linked these characteristics to their association to GMYC ESUs. We observed across ESUs the presence of spinelets, additional facets, and bending of the ridge (Table 6). In addition, we observed small humps in the middle of the areolation section and the symmetrically situated lateral antenna (Fig. 6). No correlation was found between phylogenetic and morphological diversity (Mantel r = 0.07; P = 0.41).

Discussion
Our study indicated that eight putative ESUs of K. cochlearis occurred in lakes of the Trentino-South Tyrol region. This diversity may be responsible for the apparent tolerance of K. cochlearis to varying environmental conditions. The putative ESUs of K.
cochlearis had an average uncorrected genetic distance in COI between 12 and 30%, which is higher than the 3% threshold commonly used to separate species for most animals Tang et al., 2012). The general good agreement of the various methods that we used to infer putative ESUs corroborated our results. We did not consider the splitting of GMYC ESUs 3 and 6 by PTP-BA because it was not supported by the branching pattern of the tree and the other methods of species delimitation.
The wide morphological variability in K. cochlearis that led to the description of morphotypes Table 4 Length measurements of main lorica characteristics based on 138 specimens of K. cochlearis, lorica length (excluding anterior and posterior spines, LL), posterior spine length (PSL), total lorica length including all appendages (TLL), lorica width at its widest part (LW), posterior spine angle (PSA), anterolateral dorsal spine length (ALS), anterointermediate dorsal spine length (AIS), anteromedian dorsal spine length (AMS), and lorica width beneath the anterior spines (''head width'', HW) The number of individuals measured is given between brackets next to the lake name by Lauterborn (1900) has been investigated by many researchers who tried to understand factors influencing morphology such as temperature (Green, 2005), -Porcuna et al., 1993), maternal effect (Stemberger & Gilbert, 1984), or presence of distinct species (Ahlstrom, 1943;Eloranta, 1982). Our study indicated that neglecting presence of ESUs of K. cochlearis might have led to biased conclusions on their morphological variability and global distribution. For example, the correlation between LL and PSL is not always positive as stated by Green (2005) but seems to differ between ESUs showing no correlation or varying positive correlation. Furthermore, Green (2005) underlined that specimens with a LL of around 80 lm show a wide variability in PSL. We observed an overlap of specimens of different ESUs in the range of 80-90 lm. Thus, neglecting ESUs of K. cochlearis may lead to underestimating their phenotypic diversity. An important characteristic for the delimitation of K. cochlearis morphotypes is the presence and length of the posterior spine. Our study indicated that spined and unspined (=tecta) specimens occurred in the same and different ESUs (i.e. ESUs 3 and 6, respectively). Hofmann (1983) and Green (2005Green ( , 2007 noted that tecta specimens could not be explained by allometric growth because specimens with spines were smaller than those without spines. Green (2005) presented three hypotheses of the origin of spineless K. cochlearis: 1, true tecta (appearing only in colder periods of the year as the ''end'' of the posterior spine reduction); 2, aspina (truly spineless, absent in the winter, LL longer than in spined form); 3, ecaudata (the same dorsal structure, occurring in summer, LL longer than in spined form). Coherent with Green's (2005) hypothesis 1 of true tecta, our study indicated based on ESU 6 that spineless forms have the same and  Table 5 Morphological parameters showing statistical significant differences in ANOVA between different GMYC ESUs (only significant comparisons are shown), lorica length (excluding anterior and posterior spines, LL), posterior spine length (PSL), and posterior spine angle (PSA), degrees of freedom (df), 138 specimens were measured, but ESU 7 was excluded from analyses because only specimen was measured, for ANOVA on PSL, specimens without spine were excluded (7 of ESU 4 and 3 of ESU 6), in mixed modelling of ANOVA for PSL and LL, measurements from the same lakes were modelled as correlated and in mixed modelling of ANOVA for PSA, residuals were allowed to have a different spread per lake  Green's (2005) hypotheses 2 and 3, spineless specimens of ESU 4 were smaller and larger than spined morphotypes across habitats and those from the same lake. Thus, neglecting the co-occurrence of different ESUs in K. cochlearis leads to the odd situation that spineless specimens seem larger than spined ones. We suggest that tecta morphotypes can actually have at least two possible origins (Green's hypotheses 1 and 2/3) but delimiting true tecta from spineless aspina or ecaudata based on morphology seems quite tricky. We, furthermore, hypothesise that detailed SEM pictures of lorica facets might reveal features (such as the X-facet or carinal facets described by Lauterborn, 1900) helpful for delimiting putative ESUs. Spinelets and the bended ridge are other morphological features that are used in morphotype delimitation (Lauterborn, 1900) but the usefulness of spinelets was already questioned by Hofmann (1980). According to Lauterborn (1900) spinelets are characteristic for the hispida and irregularis series. However, specimens from GMYC ESUs 3 and 6 did and did not have spinelets. According to Hofmann (1980), the size of spinelets increases from spring to summer and are almost invisible during winter. In fact in our samples, specimens with spinelets occurred during summer and spring (only one was collected from Lake Terlago during November), but we cannot exclude that we missed the presence of spinelets in some specimens as they were very difficult to observe. However, it seems that spinelets are only appearing (and changing in length) in some ESUs because no spinelets were ever observed in ESU 5 regardless of sampling time. Thus, we suggest that the presence of spinelets is no valid criterion for delimitating morphotypes or putative ESUs. Ahlstrom (1943) and Eloranta (1982) already pointed out that the presence of spinelets shows high variability in most K. cochlearis species, and here we corroborated their statement with genetic data. More detailed SEM pictures of various putative ESUs taken from different seasons are, in any case, needed in order to investigate the temporal appearance of spinelets. According to the bended ridge, specimens of ESUs 4 and 5 always showed it while it was present or absent in specimens of ESU 3. We conclude that the bended ridge is also not a valid character to delimitate ESUs. In addition to spinelets and the bended ridge, we observed small humps in the middle of the areolation section. To the best of our knowledge, we do not know about any reference to these structures. We refrain from hypothesising on their function, and if they grow, they seem to be an overlooked feature of lorica morphology. Furthermore, we provided detailed SEM images on the lateral antenna that was previously only shown by Garza-Mouriño et al. (2005, their plate 1c).
Taking into account all information on lorica morphology, different ESUs showed different morphological variabilities. Both univariate and multivariate analyses indicated that ESUs 3 and 6 were not unambiguously distinguishable based on lorica measurements showing a wide phenotypic plasticity. Contrarily, ESU 8 could be distinguished from ESU 5 based on morphology based on single measurements and NMDS. In LDA, only specimens of ESU 5 were correctly assigned in most cases, while specimens of ESU 8 did not perform that well. Specimens of ESU 8 were smaller with respect to measured characters than specimens of ESU 5. Therefore, it is possible to delimit only some putative ESUs having a more restricted phenotypic plasticity with respect to other ESUs based on detailed lorica measurements. We suggest that an analysis of specimens sampled separately during cold and warm seasons in specific water layers could provide insights into the effect of water temperature on spine development of ESUs that we may have missed by our sampling strategy.
In many of our study lakes, different ESUs of K. cochlearis co-occurred. Generally, it is assumed that species with similar morphology and close phylogenetic relationship might have similar niches (e.g. Wiens & Graham, 2005;Wiens et al., 2010) and this would lead to competitive exclusion (Violle et al., 2011;Gabaldón et al., 2013). Cryptic species are not only morphologically similar but also phylogenetically closely related, and thus, the co-occurrence of cryptic species should be rarely encountered. However, cryptic species of B. plicatilis occur in temporal co-existence or in overlap, and their co-existence is mediated by disturbance and food partitioning (Ciros-Pérez et al., 2001a). Not only in the genus Brachionus but also in P. dolichoptera (Obertegger et al., 2014) co-existence of cryptic species has been observed. We found that several morphologically similar putative ESUs of K. cochlearis cooccurred but, at the moment, cannot infer their niche  Fig. 6b, ** shown in Fig. 6e partitioning because of missing information regarding their depth distribution. Furthermore, our study indicated no link between phylogenetic and morphological diversity of putative ESUs. Similarly, Gabaldón et al. (2013) found no difference between cryptic species of B. plicatilis for key parameters (i.e. clearance rates, starvation tolerance and predation susceptibility) related to body size. Recently, co-existence of cryptic species was linked to a negative feedback based on sex-based mechanisms that lead to stable co-existence (Montero-Pau et al., 2011).
In conclusion, our study indicates that K. cochlearis is composed of eight putative ESUs based on mtDNA, as indicated by three different methods. The generally good agreement between these methods enhances our inference on species diversity. Several morphological characteristics such as presence/absence of the posterior spine, spinelets, and bended ridge seem to be of poor value to discriminate ESUs. However, when all lorica measurements are taken together in a multivariate statistical approach, ESU 5 could be distinguished from ESU 8. More detailed morphological research is needed for a longer period to understand the morphological variations of K. cochlearis ESUs.