Introduction

Tardigrades, also known as water bears or moss piglets, represent an invertebrate phylum closely related to onychophorans and arthropods within Panarthropoda1, which consists of ca. 1500 described species to date2. These are grouped into two classes: Heterotardigrada and Eutardigrada. Heterotardigrades are characterised by the presence of cirrus A3,4,5, whereas eutardigrades are generally distinguishable by their simplified, vermiform external morphology6. Tardigrades inhabit limno-terrestrial and aquatic, both freshwater and marine, habitats, the former typically in association with bryophytes, algae, lichens, or leaf litter5,7,8. Although the number of documented tardigrade species has increased significantly in the past three decades, most regions globally remain sparsely sampled for these diminutive animals9,10. The distributions of a few terrestrial tardigrade species are well documented for some regions of North America, e.g. see11, but the general paucity of North American tardigrade specialists has left much of the continent unexplored in this regard (most recently reviewed in12,13).

The Echiniscidae, the most speciose family within the class Heterotardigrada14, contain the morphologically odd genus Viridiscus15,16,17. The former Echiniscus viridis group was established by Ramazzotti18, who considered E. viridis Murray, 191019, E. perviridis Ramazzotti, 195918, E. viridissimus Péterfi, 195620, and E. rufoviridis du Bois-Reymond Marcus, 194421 to be closely related based on green cuticular pigments in plates forming the armour, alongside orange pigments typically present in the body cavity of echiniscids22. Two further species were added to this group many decades later: E. viridianus Pilato et al., 200723 and E. clavispinosus Fontoura et al., 201124. Viridiscus was erected based on morphological and molecular data to accommodate the distinctiveness of this group: dark body pigments, well-developed, sabre-like claws, a lack of trunk appendages except cirrus A, and typical dorsal plate sculpture composed of an intracuticular sponge layer and small, flat, densely arranged epicuticular granules15. Recently, new evidence was presented, pointing out possibly large intraspecific variability within Viridiscus, as two morphs were distinguished within V. viridissimus16,17. Moreover, V. rufoviridis was transferred to Barbaria25, rectifying the earlier assignment of this species within the viridis group18,23,26. Consequently, Viridiscus currently comprises five species: V. clavispinosus, V. perviridis, V. viridianus, V. viridis, and V. viridissimus.

Here, we describe a new species of Viridiscus found in previously described Viridiscus assemblages in Tennessee16,27 and identified as V. aff. viridianus in17. We also redescribe V. viridianus based on multiple populations collected in Alabama and Florida. Newly obtained data confirm and strengthen previous findings that Viridiscus species are variable in the dorsal plate sculpturing17 and morphology of the primary clavae. We update the genus phylogeny, and question the validity of V. clavispinosus, which we hypothesise to be a synonym of V. viridianus. We stress that there is no convincing evidence for the presence of V. viridis outside of the Hawaiian Archipelago26, and the Nearctic records28 of this species are unreliable. The case of Viridiscus variability yet again illustrates that forming taxonomic conclusions based solely on morphology is unadvised and should be renounced17, especially since DNA sequence-based tools are readily available to incorporate into the systematics of these organisms.

Results

General remarks: morphology and phylogeny

There are two general morphotypes distinguishable within Viridiscus with respect to the dorsal plate sculpturing: (I) pore-dominated, in which epicuticular granules are restricted to anterior portions of paired segmental plates (a usual morph of V. viridissimus, see Fig. 1A; however, we found several specimens with developmental aberrations embracing fusion of different elements of dorsal armour, see Fig. 1B); and (II) granule-dominated, in which pores are absent (the remaining four species, represented herein by V. viridianus: Fig. 1C). This division is also clearly visible in SEM: pores of V. viridissimus (Fig. 2A,B) are approaching the size of epicuticular granules of V. perviridis (Fig. 2C,D). However, the SEM analysis revealed also another element of sculpturing in the second Viridiscus morphotype: micropores,  1 μm in diameter and distributed irregularly between the granules (Fig. 2C,D).

Figure 1
figure 1

Two main morphotypes present in Viridiscus (PCM): (A) pores present, epicuticular granules typically reduced: the usual morph of V. viridissimus (Tennessee, dorsal view); (B) an aberrant specimen of V. viridissimus (Tennessee, dorsolateral view, the first median plate merged with the first paired segmental plate); (C) epicuticular granules dominant, pores present only in larval stage: V. viridianus (Alabama, dorsolateral view). Scale bars = 50 μm.

Figure 2
figure 2

Two main morphotypes of dorsal plate sculpturing present in Viridiscus (SEM): pores dominant, only V. viridissimus (the population from Vietnam): (A) a fragment of the scapular plate; (B) close up of the posterior portion of the second paired segmental plate; epicuticular granules dominant, micropores visible only in SEM, all remaining Viridiscus species (the population of V. perviridis from Vietnam shown): (C) a fragment of the scapular plate; (D) close up of the central portion of the scapular plate. Scale bars in μm.

Viridiscus viridissimus exhibits considerable intra-specific variation in the morphology of dorsal armour17. In the typical morph (Fig. 3A), reliably recorded from the Holarctic and Oriental regions17, epicuticular granules are limited to the anterior portion of the second median and paired segmental plates, and to the third median plate. In contrast, the much less common morph of V. viridissimus (initially described as V. miraviridis Nelson et al., 202016), so far identified only among moss samples from Tennessee, exhibits epicuticular granules in all plate portions (Fig. 3B). Having the possibility to analyse abundant populations of V. perviridis and V. viridianus, we discovered atypical morphs in these species, too. In V. perviridis, the usually well-discernible epicuticular granules (Fig. 3C) can be poorly developed (and therefore blurred with endocuticle in PCM), especially in the scapular and caudal (terminal) plates (Fig. 3D). In V. viridianus, the reduction of epicuticular granules (Fig. 3E) can be even more pronounced to the extent that granules are absent, and only the intracuticular sponge layer is identifiable in PCM (Fig. 3F). These atypical morphs were not associated with a particular life stage or sex but appeared in large monospecific populations of the analysed Viridiscus spp. No atypical morphs were observed in the case of V. celatus sp. nov. (Fig. 3G,H), but the available sample size was significantly smaller in the case of the new species compared to the other analysed Viridiscus spp. As already stated in the diagnosis of the genus15, all known larvae of Viridiscus spp., irrespectively of the adult morphotype (type I or II), possess large pores beside granules (Fig. 4). We confirmed the presence of such pores in larvae of V. perviridis, V. viridianus, and V. viridissimus, and they have been also detected in V. clavispinosus24. Larvae of V. viridis s.s. (inhabiting the Hawaiian Archipelago) have never been found19,26.

Figure 3
figure 3

Intrageneric and intraspecific variability in Viridiscus (PCM): (A) V. viridissimus (Tennessee), the typical morph; (B) V. viridissimus (Tennessee), the atypical morph with well-developed epicuticular granules16; (C) V. perviridis (Madeira), the typical morph with well-developed epicuticular granules; (D) V. perviridis (Alabama), the atypical morph with poorly delineated epicuticular granules; (E) V. viridianus (Alabama), the typical morph with well-developed epicuticular granules; (F) V. viridianus (Alabama), the atypical morph (male) with epicuticular granules absent (arrowhead points out conoid primary clava); (G,H) V. celatus sp. nov. (Tennessee), the only morph observed. Scale bars = 20 μm.

Figure 4
figure 4

Larva of V. viridianus (PCM, dorsolateral view). Scale bar = 20 μm.

The phylogeny fully conformed with the morphological analyses, indicating the presence of four species (Fig. 5), relationships among which were as follows: (((V. perviridis (V. celatus sp. nov. (V. viridianus + V. viridissimus))). Viridiscus celatus sp. nov. was classified in17 as V. aff. viridianus. Importantly, species delimitations based on COI did not distinguish between V. perviridis and V. celatus sp. nov., which are clearly separated based on both ITS markers (Fig. 5) and morphology, pinpointing the questionable utility of the COI marker in tardigrade species delimitation. Given the abundant material and the proximity of the type locality in Auburn (Alabama) to the localities sampled in this study, we redescribe V. viridianus to provide a detailed insight into its intraspecific variability.

Figure 5
figure 5

Phylogenetic relationships of the genus Viridiscus: Bayesian tree based on the concatenated ITS-1 + ITS-2 dataset (1058 bp); vertical bars denote different delineation methods used in the formulation of the primary molecular species hypotheses: ASAP, ABGD, and bPTP; (COI) refers to COI delimitation in all three methods). Asterisks indicate the maximal (1.00/100) posterior probability/bootstrap value. Echiniscus succineus was used as an outgroup. Scale bar represents substitutions per position.

Integrative redescription of Viridiscus viridianus (Pilato et al., 2007)

(Tables 1, 2, 3, 4, Figs. 6, 7, 8, 9, 10, 11, raw morphometry in Supplementary Material 1).

Table 1 Measurements (in µm) of selected morphological structures of adult females of Viridiscus viridianus mounted in Hoyer’s medium. sp the proportion between the length of a given structure and the length of the scapular plate, ? unknown.
Table 2 Measurements (in µm) of selected morphological structures of adult males of Viridiscus viridianus mounted in Hoyer’s medium. sp the proportion between the length of a given structure and the length of the scapular plate, ? unknown.
Table 3 Measurements (in µm) of selected morphological structures of juveniles of Viridiscus viridianus mounted in Hoyer’s medium. sp the proportion between the length of a given structure and the length of the scapular plate, ? unknown.
Table 4 Measurements (in µm) of selected morphological structures of larvae of Viridiscus viridianus mounted in Hoyer’s medium. sp the proportion between the length of a given structure and the length of the scapular plate, ? unknown.
Figure 6
figure 6

Habitus of V. viridianus (SEM): (A) dorsal view; (B) frontal view; (C) rear view. Scale bars = 50 μm.

Figure 7
figure 7

Anterior body portion of V. viridianus (SEM): (A) dorsal view; (B) head. Scale bars = 20 μm.

Figure 8
figure 8

Cuticle of V. viridianus (SEM): (A) epicuticular granules and micropores; (B) lateral portion of disrupted epicuticle; (C,D) close up of intracuticular pillars. Scale bars in μm.

Figure 9
figure 9

Leg morphology of V. viridianus (PCM). Arrows indicate pedal platelets in central leg portions, white asterisks indicate distinctly demarcated, central oval areas in pedal platelets, and black asterisks indicate pulvini in proximal leg portions. Scale bars = 20 μm.

Figure 10
figure 10

Morphological details of V. viridianus (PCM): (AC) conoid primary clava; (D) dactyloid primary clava; (E) typical for most echiniscids, tubby Echiniscus-like primary clava; (F) subcephalic plates; (G) claws I; (H) claws II; (I) claws III. Scale bars = 10 μm.

Figure 11
figure 11

Leg structures of V. viridianus (SEM): (A) claws I; (B) claws II; (C) claws III; (D) claws IV; (E) leg morphology. Arrows indicate pedal platelets in central leg portions, white asterisks indicate distinctly demarcated, central oval areas in pedal platelets, and black asterisks indicate pulvini in proximal leg portions. Scale bars in μm.

New material examined. Populations from Alabama and Florida, 346 specimens in total were processed for PCM, SEM and DNA analyses (Table 5); additional few hundred specimens were frozen for future analyses.

Table 5 List of newly found populations used in analyses. Types of analyses: (PCM) imaging and morphometry in PCM, (SEM) imaging in SEM, (DNA) DNA sequencing. Number in each analysis indicates how many specimens were utilised in a given method (a adults, v exuvia, j juveniles, l larvae).

Type locality. North America, USA, Alabama, Auburn.

Additional localities. North America, USA, New Mexico23; the Azores, Ribeira Fria, Lages do Pico23,29,30; North America, USA, New Jersey28; Central America, the Lesser Antilles, Antigua28. Given the reported variability in the pattern of the dorsal armour, these additional localities should be verified.

Etymology. From Latin viridianus = greenish. An adjective in nominative singular.

Animals. Females (i.e., from the third instar onwards; measurements and statistics in Table 1). Body cavity with yellowish pigments (typical for most echiniscids), whereas dorsal and pedal cuticular plates light to dark green (Figs. 1C, 3E,F). Red eyes and yellow pigments present in live specimens, but dissolve after mounting in Hoyer’s medium, thus only green pigmentation persists. Body bulky (Fig. 6), with a poorly delimited cephalic region (Fig. 7). The cervical (neck) plate is well-developed, but sculptureless (Fig. 7). Weakly developed lateralmost, rectangular portions of the scapular plate with a weak sculpturing (Fig. 1C). Dorsal plate sculpturing ordinarily comprising polygonal epicuticular granules with scarce micropores, barely identifiable, even with SEM (Fig. 8A). Lateral and ventral endocuticle with intracuticular pillars, visible in PCM as minute dark dots, but identifiable in SEM only when the thin epicuticle is ruptured (Fig. 8B–D). Pillars larger and more sclerotised in proximal and central limb portions, forming longitudinal, narrow pulvini, and pedal platelets, respectively (Fig. 9). Some specimens exhibit a differently formed central pedal portion, more convex than the remainder of each platelet (Figs. 9, 11E). Areas of more sclerotised pillars always form a pair of merged subcephalic plates (Fig. 10F). Cirrus A short (< 20% of body length) and thin. A remarkable diversity of primary clava shapes: some specimens within a population have both clavae pointed and clearly conoid (Figs. 7A, 10A–C), some have both clavae dactyloid, i.e., elongated, but without a pointed tip (Fig. 10D), and in some individuals both clavae are tubby, i.e., of a typical, Echiniscus-type shape (Fig. 10E). Importantly, numerous specimens showed a mixture of these shapes, that is the clava of one specimen differed in morphology from the other one on the same specimen. Claws massive and isonych (Figs. 10G–I, 11A–D).

Males (i.e., most probably from the third instar onwards; measurements and statistics in Table 2). No detectable sexual dimorphism besides the circular gonopore.

Juveniles (i.e., the second instar; measurements and statistics in Table 3). Smaller than adults, but qualitatively like them. Gonopore absent.

Larvae (i.e., the first instar; measurements and statistics in Table 4). Body size overlaps with juveniles. Anterior portions of paired segmental plates, and median plate 2 sculptureless. Large cuticular pores in the dorsal armour. No gonopore or anus.

Eggs. Up to five orange eggs per shed exuvia, but typically fewer (see31).

Remarks. Males were present in all examined populations of the species.

Description of Viridiscus celatus sp. nov. Momeni, Gąsiorek, Nelson & Michalczyk

(Tables 6, 7, Figs. 12, 13, 14, raw morphometry in Supplementary Material 2). ZooBank registration number: urn:lsid:zoobank.org:act:7E5416A6-E49F-46C2-8BA6-3D734B2A10A4.

Table 6 Measurements (in µm) of selected morphological structures of adult females of Viridiscus celatus sp. nov. mounted in Hoyer’s medium. sp the proportion between the length of a given structure and the length of the scapular plate, ? unknown.
Table 7 Measurements (in µm) of selected morphological structures of adult males and juveniles of Viridiscus celatus sp. nov. mounted in Hoyer’s medium. sp the proportion between the length of a given structure and the length of the scapular plate, ? unknown.
Figure 12
figure 12

Holotype of V. celatus sp. nov. (PCM, female, dorsal view). Scale bar = 20 μm.

Figure 13
figure 13

Type specimens of V. celatus sp. nov. (PCM): (A) allotype (male, dorsal view); (B) paratype (female, dorsolateral view). Arrowheads indicate rudimentary papillae on legs II–III. Scale bars = 20 μm.

Figure 14
figure 14

Morphological details of V. celatus sp. nov. (PCM): (A) subcephalic plates (holotype); (B) claws III (paratype, juvenile); (C) claws IV (holotype). Scale bars = 10 μm.

Material examined. Populations from Tennessee, 19 specimens in total processed for PCM and DNA analyses (Table 5).

Type locality. 36°18′N, 82°22′W, ca. 520 m asl: USA, Tennessee, Washington County, Johnson City. Grimmia sp. mosses from a concrete cap on a brick fence post.

Etymology. From Latin celatus = hidden, concealed. The name refers to the fact that the new species was not identified as new taxon for a long time, although the locality has been extensively sampled for tardigrades (16,27 referred to the species as V. perviridis based on the identification by Maucci32, and morphological characters, including cirrus A length). An adjective in nominative singular.

Type depositories: Type series: holotypic female (slide US.081.03), allotypic male (slide US.081.02), and nine paratypes (slides US.078.03 and US.081.03-4), are deposited at the Faculty of Biology, Jagiellonian University (Kraków, Poland).

Animals. Females (i.e., from the third instar onwards; measurements and statistics in Table 6). Body medium-sized and bulky. Body cavity with yellowish pigments (typical for most echiniscids), whereas dorsal and pedal cuticular plates olive green (Figs. 1C, 3E,F). Red eyes and yellow pigments present in live specimens, but dissolve after mounting in Hoyer’s medium, thus only green pigmentation persists (Figs. 12, 13B). Except for cirrus A, with a tubby clava near the cirrophore (Fig. 12), other dorsal and lateral trunk appendages are absent. Cephalic appendages include internal and external peribuccal cirri with tubby cephalic papillae between them (Fig. 12). Dorsal plate sculpturing comprises large epicuticular granules (Fig. 12), which may be poorly developed in central plate portions (Fig. 13B). Sponge layer identifiable beneath granules. Granules appear more convex in anterior portions of paired segmental plates than in the remainder of the armour in PCM. Micropores not visible in PCM and their presence or absence remains to be confirmed in SEM.

All plates strongly sclerotised and with clear edges. The cephalic plate with a well-marked anterior chalice-shaped incision, the cervical plate and lateral sections of the body lack dense granulation and are covered with fine regular punctuation. The scapular plate contains three portions. Only the central part is visible in the dorsal view, and two small, weakly delineated, trapezoidal sections are present on the lateral portions of the body, with intracuticular pillars visible (Fig. 13B). The first median plate is triangular and unipartite, the second median plate is subdivided into two portions, and the anterior portion lacks the sponge layer. The third median plate is absent, but the area between the paired segmental plate II and the caudal plate is covered with large granules. Paired segmental plates I and II have two clearly delineated parts. Intersegmental plate is inserted between the posterolateral edge of the paired segmental plate I and anterior margin of paired segmental plate II. The caudal incisions are unsclerotised and weakly marked (Figs. 12, 13B).

Venter densely granulated in PCM (endocuticular pillars); a pair of subcephalic plates present (Fig. 14A). Gonopore hexapartite. Pulvini (= narrow proximal bands of intracuticular pillars) and pedal platelets (= broad central bands of pillars) are visible on all legs. Dentate collar with numerous irregular teeth (Fig. 14C). Sensory organs present on all legs: a tiny spine on leg I embedded at the edge of pedal platelet; hemispherical rudimentary papillae on legs II–III, embedded in the centre of pedal platelets (identifiable only when specimens are dorsolaterally oriented); and papilla IV on hind legs (Fig. 13B). Claws anisonych; primary spurs I–III tiny and thin, positioned slightly lower on branches compared to more massive spurs IV (Fig. 14B–C).

Males (i.e., most probably from the third instar onwards; measurements and statistics in Table 7). No sexual dimorphism observable in body size or qualitative traits (Fig. 13A). Gonopore circular.

Juveniles (i.e., the second instar; measurements and statistics in Table 7). Smaller than both females and males. Qualitatively like adults; gonopore absent.

Larvae. Not found.

Eggs. Not found.

Remarks. Found only in association with large populations of V. viridissimus17.

Differential diagnosis. The new species from Tennessee is differentiated from all Viridiscus spp. based on the presence of plesiomorphic papillae on legs II–III. These structures are, however, barely identifiable in specimens oriented dorsoventrally, hence we enumerate other criteria making V. celatus sp. nov. distinct from:

  • Viridiscus clavispinosus, by the relative length of cirrus A (34–49% vs < 15% of the body length), and primary spurs IV less divergent from claw branches;

  • Viridiscus perviridis, veritably reported from the Holarctic and Oriental regions17, by the length of cirrus A (34–49% vs typically   50% of the body length, see the subsection below that addresses this character), the weakly developed caudal incisions (strongly sclerotised and well-marked in all syntypes of V. perviridis), and the body colour (light to olive green vs usually dark green to almost black in V. perviridis, also in mounted specimens);

  • Viridiscus viridianus, reliably reported only from the USA, by the relative length of cirrus A (34–49% vs < 20% of the body length), the lack of pedal platelets with a distinctly formed central portion, and the more pronounced sculpturing of the anterior portions of paired segmental plates;

  • Viridiscus viridis, reliably reported only from the Hawaiian Archipelago19,26, by the relative length of cirrus A (34–49% vs < 10% of the body length), and a different pattern of dorsal sculpturing (in general V. viridis has noticeably fewer epicuticular granules on all plates, see fig. 1 in26);

  • Viridiscus viridissimus, with a likely wide distribution in the Holarctic, Oriental, and Neotropical regions17,33, by the absence of pores in dorsal armour, and a better developed sponge layer of cuticle.

Discussion

Viridiscus viridianus vs Viridiscus clavispinosus

Viridiscus viridianus was described by Pilato et al.23 based on populations from Alabama (the origin of the holotype, thus locus typicus), New Mexico (the Nearctic region), and the Azores (Macaronesia, the westernmost part of the Palaearctic region). A few years later, Fontoura et al.24 described V. clavispinosus (from the Archipelago of Cape Verde, which is also a part of Macaronesia, but located southwards of the Azores), using a seemingly sound autapomorphy, namely conoid primary clavae (with pointed apices). However, the data gathered in the present study undermine the validity of V. clavispinosus, as the shape of primary clava is evidently a variable character in some species of Viridiscus (Fig. 10A–E), rendering it unsuitable for species delineation. The surface between the paired segmental plate II and the caudal plate can be weakly sculptured (Figs. 1C, 3E, 6A,C) or unsculptured (Fig. 3F) in a population of one species. Contrary to what was stated in24, neither V. viridianus nor V. clavispinosus exhibit well-marked epicuticular granules in the cephalic and cervical plates (compare figs. 1b, e in23, fig. 1a in24, and Fig. 7 herein). In echiniscid species, the dorsal plate sculpturing is less pronounced in these two anteriormost elements of armour. The sculpturing of the less sclerotised anterior portion of median plate 2 can also be more or less developed, and sometimes absent within a single population. The remaining characters referring to the dorsal sculpturing mentioned in24 also are invalid in light of various atypical morphs presented for Viridiscus species in this paper, including V. viridianus. We refrained however, from synonymising V. clavispinosus with V. viridianus for two reasons: (1) as rightly stated by Fontoura et al.24, primary spurs IV of V. clavispinosus are more developed than those of V. viridianus (compare fig. 2c in24 and Fig. 11D herein); (2) although V. viridianus potentially has a wide circum-Atlantic distribution (whilst molecular evidence is lacking), V. clavispinosus is known from a different part of Macaronesia. Therefore, V. clavispinosus is designated as nomen inquirendum, given the substantial doubts regarding its separateness from V. viridianus described above.

Inter- vs intraspecific variability in Viridiscus and its consequences for echiniscid taxonomy

One of the prime tasks for taxonomists is recognising the borders of inter- and intraspecific variability in a given group of organisms (for microscopic animals, see e.g.34,35,36,37). In the preliminary study that integratively addressed variability of Viridiscus17, we accentuated the role of combining molecular and morphological approaches in tardigrade taxonomy, as utilising only one line of evidence either leads to taxonomic inflation (e.g., the synonymy of V. miraviridis16,17) or deflation (e.g., the Milnesium case38). We have already underscored high variability of characters previously deemed as universally stable at the species level in the Echiniscidae (cirrus A length, details of cuticular ornamentation) in Echiniscus and Nebularmis spp.39,40. This study adds further observations on atypical morphs within Viridiscus, which show disparities with typical morphs of a given species that are larger than interspecific variability in multiple new characters, such as the general phenotype of the dorsal sculpturing, the shape of primary clava, and morphology of pedal platelets. It is gradually becoming obvious that even within such a morphologically coherent family as Echiniscidae41, phenotypic plasticity of diagnostic characters varies between lineages at various taxonomic levels (genus, species, etc.), as it happens with the shape of primary clavae in Viridiscus, which is labile in contrast to other known echiniscid genera. Consequently, conservative characters in one genus should not simply be assumed to be conservative in other taxa. In parallel, analogous variability in characters such as the shape and presence of spine I and dorsal sculpturing, has been also reported recently in Claxtonia spp.42.

On the other hand, genetic data, when taken in isolation, can also be deceiving, and this pertains specifically to the COI marker, often uncritically utilised and termed as a universal tardigrade barcode as (I) it is often difficult to amplify, and (II) may fail in delimiting species (Fig. 5;43). Multiple markers are advisable44, and both ITS markers are good predictors of intra- and interspecific differences, as they are congruent with a spectrum of morphological variability analysed on many animals. Altogether, this reinforces the necessity of the integrative approach in tardigrade studies17.

In the light of our findings, some generalisations can be made regarding the diversity and classification of Viridiscus:

  1. 1.

    The green body colour is a result of dorsal cuticular plate pigmentation, but orange and yellow (probably carotenoid-derived pigments), are present in all Viridiscus spp., as in most other echiniscids. Thus, two characteristics of the dorsal armour, the presence of the endocuticular sponge layer and green pigmentation, are advanced characters.

  2. 2.

    Although more variable than previously reported, the relative length of cirrus A with respect to the body length (bo) serves as a good criterion in morphological species delineation. The extreme elongation of cirrus A characterises V. perviridis (most often   50% of bo), followed by cirri A longer than the echiniscid average (30–50% of bo) in V. celatus sp. nov. and V. viridissimus, and cirri A that can be classified as short (< 20% of bo) in V. clavispinosus nom. inq., V. viridianus, and V. viridis23.

  3. 3.

    Atypical morphs are probably present in populations of all Viridiscus spp. We recorded them in V. perviridis, V. viridianus (this study), and V. viridissimus17. This means that morphology alone does not offer a fully credible dataset for establishing new species.

  4. 4.

    The description of V. viridis by Pilato et al.23 was based on only two individuals and likely does not reflect the spectrum of morphological variability of this species, but plainly indicates the specific nature of its dorsal sculpturing. The sculpturing embraces particularly widely spaced epicuticular granules, making it dissimilar to all other Viridiscus spp. Unless integratively redescribed from the Hawaiian Islands and later verified from other regions of the globe, it should be temporarily considered an endemic of this archipelago.

Conclusions

The extent of intraspecific morphological variability within Viridiscus is considerable and exceeds the usual variation encountered within echiniscid populations. Atypical morphs of Viridiscus species are not linked with sexual dimorphism, phenology, or ontogeny, and seem to represent idiopathic deviations from the most common morphs. Viridiscus celatus sp. nov., found only in association with large populations of V. viridissimus, is described integratively based on genetics and morphology. In addition, the validity of V. clavispinosus is questioned, and its synonymy with V. viridianus is implied. Males were reported in populations of Viridiscus celatus sp. nov., V. viridianus, and V. viridissimus so far.

Materials and methods

Sample collection and processing

Previous sampling experience31,45 showed that Viridiscus species are often found in mosses and lichens growing on rapidly draining vertical surfaces (see e.g.46) such as the sides of rocks (Fig. 15), gravestones, and tree trunks in human-altered environment. Therefore, moss and lichen samples were collected from several localities in Alabama, Florida, and Tennessee from February 2020 to May 2022 (Table 5). All samples were stored in acid-free paper envelopes and allowed to dry overnight in an air-conditioned room at 25 °C and 50% relative humidity. The air-dried samples were suspended over deionised water using the custom-built Baermann pan described by Davison47 for at least 12 h to separate the motile meiofauna from the moss and lichen particles. The deionised water was then filtered through 25 µm mesh and the sieve was rinsed into 30 mm Petri dishes with locally collected rainwater to search for tardigrades. The air-dried lichens were soaked in deionised water for at least 12 h, active tardigrades were extracted with an Irwin loop and transferred into 30 mm Petri dishes containing rainwater.

Figure 15
figure 15

Habitats harbouring Viridiscus populations in Andreaea mosses growing on boulders by Lake Harris, Tuscaloosa, Alabama: (A) an overall view of the rock substrate; (B) close up of the moss growing on the vertical surfaces of the boulder (arrow indicates the moss matt); (C) close up of the pleurocarpous moss and leaf structure when dry; (D) moss leaves seen under a stereomicroscope (arrows indicate animals on moss leaves in the tun state).

Microscopy and imaging

Individual tardigrades from the samples were visually inspected for diagnostic traits using an Olympus BX53 phase contrast microscope (PCM) associated with an Olympus DP74 digital camera. All specimens were divided into morphological and molecular analyses (details in Table 5). The excess of individuals was frozen to store for future analyses. For morphology and morphometry, individuals were permanently mounted on microscope slides in Hoyer’s medium. Dried slides were sealed with nail polish and examined under the Olympus BX53 PCM. Specimens for imaging in the SEM were CO2 critical point-dried, coated with gold and examined in the Versa 3D DualBeam SEM at the ATOMIN facility of the Jagiellonian University. All figures were assembled in Corel Photo-Paint X8. For deep structures that could not be fully focused on a single PCM photograph, a series of images was taken every ca. 0.1 mm of vertical focusing and then assembled manually in Corel Photo-Paint into a single deep-focus image.

Morphometry

All individuals of the new species from Tennessee and selected individuals of V. viridianus from Alabama were chosen for morphometry. The measurements are in micrometres (µm), and only undamaged structures with the appropriate orientation were used. The body length measurement is from the anterior to the posterior ends of the body, excluding the hind limbs. The bo% was calculated for cirrus A, which is the ratio of cirrus A length to the body, and sp% is the ratio of the length of the given structure to the length of the scapular plate48. Morphometric data were handled using the Echiniscoidea ver. 1.4 template available from the Tardigrada Register, http://www.tardigrada.net/register49.

Comparative material

The morphology of Viridiscus spp. was compared with paratypes of V. viridianus and syntypes of V. perviridis (deposited in the Civic Museum of Natural History in Verona and University of Modena and Reggio Emilia, Italy), and original descriptions or redescriptions of species17,23,24,25,26. Slides of V. viridissimus from Tennessee, including the “miraviridis morphotype” (16, synonymised with V. viridissimus by Gąsiorek et al.17), were also examined.

Genotyping

A Chelex® 100 resin (Bio-Rad) extraction method was used for DNA extraction50,51. Hologenophores were recovered after the extraction and mounted on permanent slides in Hoyer’s medium in most cases; in other cases, paragenophores were preserved52; all are deposited in the Faculty of Biology, Jagiellonian University in Kraków. Five DNA fragments (nuclear markers: 18S rRNA, 28S rRNA, ITS-1, ITS-2; mitochondrial marker: COI) were amplified and sequenced according to the protocols described in51; primers and original references for specific PCR programmes are listed in Supplementary Material 3. GenBank accession numbers for the species used in calculating phylogenies and the genetic framework (mainly based on the dataset from Gąsiorek et al.17) are provided in Table 8. ITS and COI sequences were separately aligned with sequences from Echiniscus succineus Gąsiorek & Vončina, 201953 as an outgroup using the ClustalW Multiple Alignment tool54 implemented in BioEdit ver. 7.2.555. The remaining gaps were left intact in ITS alignments. The 18S and 28S rRNA gene fragments were not used for the species delimitation purposes but are provided for future broader-scale phylogenies.

Table 8 GenBank accession numbers for the sequences analysed in this work. New sequences in bold.

Phylogeny

The sequences of the ITS fragments were concatenated to generate a matrix of 1058 bp in SequenceMatrix56. Using PartitionFinder version 2.1.157 with applied Bayesian Information Criterion (BIC) and the greedy algorithm58, the best substitution model (the lowest BIC) and partitioning scheme were chosen for posterior phylogenetic analysis. As the best-fit partitioning scheme, PartitionFinder indicated one partition with the best-fitted model GTR + G. Bayesian inference (BI) marginal posterior probabilities were calculated using MrBayes v.3.259. Random starting trees were used, and the analysis was run for ten million generations, sampling the Markov chain every 1000 generations. An average standard deviation of split frequencies of < 0.01 was used as a guide to ensure that the two independent analyses had converged. Tracer v1.360 was then used to ensure Markov chains had reached stationarity and to determine the correct ‘burn-in’ for the analysis, in this case the first 10% of generations. The Effective Sample Size values were greater than 200, and the consensus tree was obtained after summarising the resulting topologies and discarding the ‘burn-in’.

ModelFinder61 was used to choose the best-fit models for two partitions in Maximum Likelihood (ML): TIM3e + G4 (ITS-1) and K2P + G4 (ITS-2), chosen according to the BIC. W-IQ-TREE was used for ML reconstruction62,63. One thousand ultrafast bootstrap (UFBoot) replicates were applied to provide support values for branches64.

Both ITS and COI alignments were uploaded separately to the Assemble Species by Automatic Partitioning (ASAP) webpage65, and automatic barcode gap discovery (ABGD) web66 to obtain six independent marker-based primary species hypotheses using uncorrected pairwise distances. The partitions with the lowest ASAP and ABGD scores and p values < 0.05 were chosen as the best-fit hypotheses. Finally, Bayesian Poisson tree processes (bPTP67) were applied to the ML phylogenetic trees run on three markers separately in W-IQ-TREE (COI best-fitted model: HKY + F + G4). The calculations were conducted with 100,000 MCMC generations, thinning the set to 100, with 10% burn-in, and with searches for maximum likelihood and Bayesian solutions. All final consensus trees were visualised by FigTree v.1.4.3 available from https://tree.bio.ed.ac.uk/software/figtree.