A Strategy to Provide a Present and Future Scenario of Mexican Biodiversity of Tardigrada

: Although the number of known tardigrade taxa in Mexico has increased signiﬁcantly in the last ten years, the knowledge of their diversity faces challenges, as more than half of the Mexican territory has no records of this phylum. Thus, we developed a strategy to provide a present and future scenario for understanding the Mexican biodiversity of Tardigrada, described the distribution patterns of the current recorded species, calculated the estimated richness, and the estimated taxonomic effort needed to complete the national inventory. We obtained 474 records of 105 taxa, belonging to 42 genera and 75 species, distributed in 12 of the 14 biogeographical provinces of Mexico. We found that 54.72% of the species are present in more than three world regions and 3.79% of species that have been recorded only in Mexican provinces. Distribution patterns could be recognized for 11 species, two of which have a Nearctic distribution, seven are Neotropical and two are distributed in both regions. The Mexican biogeographical provinces with the greatest diversity of tardigrades, both at speciﬁc and generic level, were the Transmexican Volcanic Belt (TVBP) and the Sierras Madre Oriental (SMOrP) and Sierra Madre Occidental (SMOcP), which have been previously identiﬁed as particularly species-rich regions. Diversity estimation methods predict that more than 290 species of tardigrades could be found in Mexico.


Introduction
Tardigrades are micrometazoan between 50 and 1200 µm in length, ubiquitous in marine, freshwater, and terrestrial interstitial communities, and have been found in diverse habitats [1][2][3]. They are among the most desiccation-and radiation-tolerant animals surviving even to extreme levels of ionizing radiation [4][5][6][7]. Given the great variety of environments in which they can live, tardigrades have been required to allocate significant amounts of energy toward specific adaptive strategies to survive, such as resting stages, generally defined as dormant, stages involving a temporary suspension of active life, reduced or suspended metabolism, and arrested development [8]. The great resilience of these organisms has allowed them to disperse through air and water, and recently, possible evidence shows that birds could play a role as long-distance dispersers of tardigrades have been found [9].

Distribution Patterns of the Tardigrades of Mexico
The biogeographic patterns exhibited by the taxa present in Mexico show a large complexity. In this regard, various hypotheses have been raised to explain the assembly of different taxa with diverse origins [58]. Thus, to determine the biogeographic provinces in which the tardigrades are distributed, all records at the genus and species level were projected on the layer of Mexican biogeographic provinces by Morrone et al. [59], and the marine ecoregions [60]. To evaluate in a global context the distribution patterns of the taxa found, worldwide information was integrated using Morrone's biogeographic regionalization [61], which is organized into kingdoms, regions, subregions, and domains. For Mexico, the regionalization into Mexican provinces was used [59]. With this information, it was determined whether the species were found in small regions or are broadly distributed (i.e., present in more than three biogeographic regions). We also described the number of records and species per biogeographic province and for each taxon, the provinces where their records were found and identified.

Estimating Richness
There are plenty of methods to estimate species richness, making it difficult to conduct direct comparisons among them. Therefore, we decided to follow the work of Stork [62], who provided a framework and a reference of existing methods. In his work, seven estimation methods are mentioned; however, given the nature of our available dataset (a diverse range of data from literature; with no field collections was made) only two of them were applied: (1) ratio of known to unknown species, and (2) accumulation curves.

Ratio of Known to Unknown Species
This method is used to estimate the number of species of a poorly known group ("problem taxon") in a determinate geographic region, based on the proportion of species of other well-documented taxa in the same region, concerning their diversity worldwide. To do this, we used the information on the number of species of widely studied groups. A proportion or percentage of known species for those taxa in Mexico was obtained concerning the rest of the world; with the percentages obtained for each taxon, an extrapolation was made of the total number of tardigrade species recognized worldwide [62]. In this way, well studied, representative groups that have a number of species present in Mexico and a number of species worldwide were selected; in this case, we started from a Mexican and global known diversity of 15 eukaryote groups, which in turn were obtained from articles and databases: vascular plants [63,64], mosses [65,66], birds [67,68], reptiles [69,70], mammals [71,72], fishes [73,74], arthropods [75,76], proturans [77,78], diplurans [78,79], copepods [80,81], onychophorans [82,83], nematodes [84], rotifers [85] and cladocerans [86]. The number of currently described species of Tardigrada was obtained from the Actual Checklist of Tardigrade 40th ed [21]. Likewise, a global estimate of tardigrades was made by Bartels et al. [50]; these data were also taken as a second reference of estimated tardigrade diversity. In this way, two references were obtained, the first for current diversity with a value of 1380 species and the second one for estimated diversity with a value of 2654 species [50]. Hence, the use of these values were applied as follows: if Rotifera has 3000 known species worldwide, and 184 rotifer species are reported in Mexico, then 6.1% of the global rotifer diversity is present in Mexico. If this ratio is extrapolated to phylum Tardigrada, based on the 1380 described species, 6.1% of this diversity represents 84 species in Mexico, but if we use the richness estimation by Bartels et al. [50] (n = 2654), we could expect as much as 162 species. The calculation for the rest of the phyla was performed in this way. Since most of the 1380 species of tardigrades known to date have been described from morphological data, and current integrative approaches uphold that these attributes are not sufficient to support and recognize them, it is likely that there are synonymous or cryptic species in this number. For this reason, calculations obtained from methods that consider the current number of species should be taken as a conservative estimate of the possible Mexican tardigrade diversity.

Accumulation Curves
This method is based on the estimation of the predicted asymptote for the cumulative number of species over time. It uses the information of the presence and frequency of new species over time, so the closer the curve is to the asymptote, the more difficult it is to find new taxa. From the literature previously consulted in the list of Mexican tardigrades (in the first section of the methods), collection data were extracted from the sites where tardigrades were recorded, since the accumulation curve method uses standardized sampling units and in the reviewed works this unit is not defined, and so each geographical locality corresponding to a collection event was taken as a sampling unit, recognizing a total of 136 units from 25 bibliographic sources. Given that in some studies the taxonomic assignation was carried out only at a genus level, two accumulation curves were elaborated (genus and species categories). For the curve at the genus level, the data from the collections identified to species and genus were combined. However, it should also be noted that this curve was obtained using data from sources other than indexed journals, such as theses, scientific outreach, conference proceedings, and posters, which in turn tend to be perceived as taxonomically "less rigorous" [35,37,39,41,87,88]. Taxonomic modifications of the taxa, as synonymies or genus changes, were also updated from the species checklist [21].
The curves were estimated following the protocol of Jiménez-Valverde and Hortal [89], which represents the increase of the taxa in the inventory according to the performed sampling. This generates an "ideal curve", which is based on randomization of the sampling effort and the mean number of species [90]. Subsequently, the quality of the sampling n = a/(1 + b·n) 2 and the sampling effort required to record 95% of the fauna n 95 were calculated with the a and b constants of Clench's equation S n = a·n/(1 + b·n). This model assumes that the probability of adding new species to the inventory increases as more time is spent in the field, and is recommended for both large study areas and protocols [91].

Taxonomic Effort
Taxonomic effort is a function of search rate (collection effort) and species description (taxonomic description effort) [92]. The first one was calculated from accumulation curves (see above). Taxonomic description effort implies the amount of effort, after collection, invested in describing new species; this task is carried out by specialized taxonomists [93]. Therefore, in this section, the taxonomic identification effort was calculated as the number of specialists that would be needed to describe potential species. To estimate the taxonomic effort needed to complete the inventory of the Mexican tardigrades, the taxonomic effort made in the rest of the world was evaluated as a reference that could be applied to the Mexican scenario. To this end, we evaluated the relationship between the number of species described and the number of authors participating in this taxonomic activity worldwide since 1834. All the species, their year of description, and corresponding authorities were listed from the 40th edition of the actual checklist of tardigrade species [21]; from this list, we obtained the average time spent in the taxonomic field (how many years an author remains describing species) and the average number of species described by each author.
On the other hand, based on the results of the species estimation methods, the number of species missing to complete the inventory was calculated and the rate of species description in the world and Mexico was estimated. We do not attempt to assign a deterministic value to the number of authors and the time required to inventory the Mexican tardigrade fauna, since these variables depend on numerous factors that were not considered in the present study (e.g., the level of experience of taxonomists, the number of research groups and their infrastructure-national or international-, the economic resources invested, the rate of habitat loss), those are factors that could considerably affect the rate of species description. Rather, it should be considered that our estimates of potential diversity are based on biased samples from a few regions within the country and that most of the Mexican provinces have few or no records, particularly provinces characterized by their high species richness, which in turn could contribute significantly to the discovery of new taxa of Tardigrades in the country.

List of the Mexican Tardigrades
We quantified 474 tardigrade records ( Table 1). One of them was identified only at Order (Eutardigrada), one at subfamily (Florarctinae), one at the family level (Halechiniscidae), 238 at the genus level, and 233 at the species level. These records correspond to 105 taxa from 18 families, 7 subfamilies, 42 genera, and 75 species. Seven genera were recorded more than  According to the bibliographic sources, we obtained 311 of the records from theses, 151 from indexed journals, nine from conference proceedings, and three from scientific outreach journals [35,37,39,41,87,88].

Ratio of Known to Unknown Species
The Mexican and global diversity data for the 15 selected taxa for the diversity estimation of Tardigrada are shown in Table 2. Based on the percentage of known species of these 15 groups, the current [21] and estimated diversity data of tardigrades [50], it the number of species in Mexico was calculated. We excluded the values of the taxa that quantified fewer than 75 species, which is the number of current tardigrades records in Mexico and obtained in this work. When the current diversity value of tardigrades was used (1380 spp.) we could find between 77 and 262 Tardigrada species, and between 126 and 505 considering the value of estimated diversity (2654 spp.) Considering the percentage of diversity present in Mexico obtained by Martínez-Meyers et al. [24], it is estimated that between 117 and 2265 species could have occurred in the country ( Table 2). Table 2. Ratio of known to unknown species of 13 phyla in the world, species recorded, percentage of occurrence, and the number of expected tardigrades species in Mexico. Values were obtained by extrapolating the percentage of occurrence to the number of tardigrades species in the world, of each group. The numbers in bold indicate estimated species values >75, that is, the number of current tardigrades records in Mexico obtained in this work.

Accumulation Curves
The estimation based on the accumulation curves of both genus and species were adjusted to the Clench model (R 2 = 0.99; Figure 3, Tables 3 and 4). According to this, 78% of the genera and 39% of the species that could be found in Mexico have been effectively recorded, and around 572 samples would be needed to record 95% of the tardigrade genus and 2714 samples to complete 95% of the species in Mexico; according to the applied model, 292 species and 44 genera would be expected (Figure 3). The slope at the end of the curve at the genus level is lower than 0.1, which indicates that the inventory can be considered sufficiently reliable (m genera = 0.089). Although it is still incomplete, the function supports that 22% of the genera would remain undiscovered. The slope at the end of the curve at a species level is greater than 0.1, which indicates that the inventory is incomplete and unreliable (m species = 0.78) as the curve is far from reaching the asymptote, and 61% of the species have yet to be discovered. . The X-axis shows the sampling effort (each geographical locality corresponding to a collection event was taken as a sampling unit). The Y-axis represents the number of species retrieved for each given sampling level. The solid lines represent the Clench function fitted to the curve. The dotted line indicates the asymptote predicted by the model.

Taxonomic Effort
When relating the number of described species and the number of authors describing the species, we observed that the species known so far are the result of the effort of 320 authors. The histogram (Figure 4) shows an increase in the number of descriptions over time; however, there are ups and downs in the number of descriptions, one of the first increases occurred in the decade of 1904, when 57 species were described by two authors, this increase strongly was diminished by the World War I, two decades later, in which two authors described only two species.  By calculating the number of years when that the names of taxonomic authorities were recorded in the literature, we determined that one author actively describes species on average every six years (6.3 ± 0.58); however, more than half of the authors (60.31%) appeared to describe species for only one year. One author describes seven species (7.5 ± 0.99), but 44.37% of them only describe one, 14.37% describe two, and 7.5% described three species.

Discussion
Based on our strategy for the study of tardigrade diversity in Mexico, we collected 474 records for 105 taxa, belonging to 42 genera and 75 species, distributed in 12 of the 14 biogeographical provinces of Mexico [59]. With these data, it is possible to recognize that the Mexican tardigrade fauna is made up of a greater proportion of species occurring in more than three biogeographical regions of the world, followed by species only reported in the Americas and, to a lesser extent, species that have only been recorded in Mexico. Diversity estimation methods predict that more than 290 species could be found in Mexico, and 2714 samples would be required to find them.
List of the Mexican tardigrades. The study of the tardigrade fauna in Mexico began in 1911 with the work of Heinis [26] who reported four species, one new species and two genera from samples of moss and lichen from Oaxaca. Since then, to date, only two lists of the group have been published [25,33]. The first corresponds to that by Kaczmarek et al. [33] in 2011, who provided the first complete list of species documented over a century , which includes four new records and quantifies a total of 41 species. The second one is of Moreno-Talamantes et al. [25] in 2019, where 55 species were recognized over eight years, including 11 new records and three previously described species [39,40,45]. Although numerous records were generated in the period 2011-2019 in non-specialized literature, both of genera [34][35][36][37]42,44,87] and species [41,88], these were not included in the latter listing, nor in the calculations of national diversity presented in later works [43,47]. From 2019 onwards the diversity of tardigrades increased notably in Mexico, due to the description of new species [46][47][48][49] and the recognition of new records at the specific and generic level [38,43,47,49]. For example, in 2020 alone, 61 species were recognized [47] and in 2021 84 species [43]. However, these works did not provide a list of taxa in support of their calculations on the number of species.
According to our results of a detailed literature review on Mexican tardigrade fauna, 105 taxa are recognized for the country, of which 75 were "identified at species level", 27 are assigned at genus, two to a subfamily, and one to order levels. Of the 75 taxa identified at a species level, only for 62 of them were the authors able to support the specific assignment, with qualitative and quantitative characteristics, while in 13 (Dipodarctus cf. subterraneus, Echiniscus cf. tamus, Halechiniscus cf. perfectus, Pseudechiniscus cf. juanitae, Calcarobiotus cf. polygonatus, Calohypsibius cf. ornatus, Hypsibius cf. convergens, Hypsibius cf. microps, Hypsibius cf. pallidus, Macrobiotus cf. acadianus, Minibiotus cf. intermedius, Paramacrobiotus cf. klymenki and Ramazzottius cf. oberhaeuseri), due to lack of congruence in the quantitative measures, it was only possible to assign with certain probability their belonging to the identified species, and therefore they constitute putative species. Some records at the genus level, some are specimens which lack either body morphological attributes or eggs, hindering assignation to a specific level. Several of them are potentially new species represented by unique individuals that have contrasting characters with the species recognized for those genera [37]. The record at the subfamily level (Florarctinae) [37] was supported by specimens that were subsequently identified at the species level [46], and the record at the family level (Halechiniscidae) comes from a single specimen that presents unique characteristics, suggesting the occurrence of a new genus.
Based on this list, we recognized that 28.3% of the taxa have not been determined at a specific level, and in others, it has not been possible to delimit new taxa based solely on morphological characteristics, which is evidence of the absence of an integrative strategy for the study of tardigrades diversity. The use of integrative methods that include different types of observation of morphological characters, and molecular markers focused on the delimitation of taxa, has been documented in other works for the phylum [18,19,[103][104][105]. The use of a multidisciplinary approach that takes advantage of the complementarity between different techniques of analysis of morphological [106] and genetic variation analyses [107] will allow clarifying the status of the taxa in Mexico and will provide more robust identifications and hypotheses testing of species.
An example of the urgent need to use integrative methods of species description is the number of DNA sequences deposited in the GeneBank database for described tardigrade species so far (more than 1300). It is observed that molecular information can be found for only 294 species, of which 123 species include fragments of the mtDNA Cytochrome Oxidase gene, four species have fragments of the 16S rRNA, 163 species have fragments of the 18S rRNA gene, 135 species have fragments of the 28S rRNA gene, and 156 species have ITS sequences [22]. However, some of these sequences have been obtained for purposes other than taxonomy, so fewer than 294 species have been described with molecular information sources.
Distribution patterns of the Mexican tardigrades. Our records support the presence of tardigrades in 17 of the 32 states of the country. Nuevo León, Mexico State and Chihuahua, were the states with the most records (23, 20 and 14 respectively). The Mexican biogeographical provinces with the greatest diversity of tardigrades, both at species and genus level, were the Transmexican Volcanic Belt (TVBP), the Sierra Madre Oriental (SMOrP), and the Sierra Madre Occidental (SMOcP), which have been previously identified as particularly species-rich regions [108].
Looking at the distribution pattern of tardigrade taxa records in the biogeographical provinces, it is evident that exploration of the phylum has not been uniform and there are many unstudied areas (Figure 1). In the TVBP, records are concentrated in the east of the province (Popocatépetl and Iztaccíhuatl volcanoes); in the SMOrP and SMOcP, records are clustered in the northern region (Nuevo León and Tamaulipas, and Chihuahua respectively); in the Yucatán Peninsula (YPP), all records are located in the northeastern region of the province (Quintana Roo); in the Veracruzan province (VP) the records are concentrated in the southwestern region (Chiapas) and in the Tamaulipan province (TP) the records are concentrated in the northern region (Nuevo León, Tamaulipas, and Coahuila). There are also provinces with few records, such as Sierra Madre del Sur (SMSP), Pacific Lowlands (PLP), and Balsas Basin (BBP), which had less than 5 records per province. These areas have the potential to considerably increase the inventory of the Mexican tardigrade fauna, since in other groups such as vertebrates, vascular plants, and arthropods, a great richness of species is recorded in the states of Oaxaca, Veracruz, and Chiapas [109].
Of the 105 recognized taxa in Mexico, 54.72% have a wide distribution (occurring in more than three biogeographical regions of the world), 14.14% have records in two biogeographical regions and 31.12% have only been documented in American provinces.
The latter corresponds to 30 taxa, of which only 11 species distribution patterns could be recognized based on their current records. One of these patterns is the Nearctic, found in Macrobiotus anemone and Pilatobius nodulosus, which have records from Canada, USA, and Mexico. Another pattern was Neotropical, found in the species Doryphoribius quadrituberculatus, Diphascon mitrense, Echiniscus kofordi, Macrobiotus alvaroi, Macrobiotus kazmierskii, Milnesium barbadosense, and Minibiotus sidereus, with records from the Mexican Transition Zone (MTZ) [110,111], Argentina, Chile, Costa Rica, Ecuador, and Barbados. Two species also extend their distribution from the MTZ to the Nearctic and Neotropical regions: Doryphoribius dawkinsi, found in the USA, Mexico, and Costa Rica, and Diphascon chilenense, found in Canada, USA, Mexico, and Chile.
Some tardigrades with Neotropical records have distribution similarities with other groups of animals and plants. Macrobiotus alvaroi was detected in the Mesoamerican (Chiapas) and Pacific (Costa Rica) domains of the Brazilian subregion; a similar distribution has been recognized in the coleopteran Deltochilum acropyge [110], whose distribution was attributed to the Neotropical cenocron with minimal penetration; the taxa supporting this pattern have a South American origin and have entered into the large rainforest patches of the southeast, in Chiapas and Campeche [112].
Other examples are Echiniscus kofordi and Milnesium barbadosense, species that are distributed in different subregions, from the Neotropic to the Nearctic, and show similar distribution patterns with some moss species. Echiniscus kofordi is found in the Brazilian (South American transition zone domain in Ecuador and Mesoamerican domain in Chiapas) and Allhegani (Louisiana, Alabama, and Florida) subregions, a distribution that is consistent with the moss species Atractylocarpus flagellaceus, representative of the Mesoamerican element of the moss flora, which is constituted by 174 species that extend mostly from Mexico to South America, but in some cases, they present a wider distribution towards the southeastern USA [113,114].
Milnesium barbadosense is found in the Antillean (Barbados) and Brazilian subregions (Mesoamerican Domain in Chiapas and the Pacific Domain in Colombia), a distribution that is in agreement with the moss species Breutelia jamaicensis, Calyptothecium duplicatum, and Helicophyllum torquatum, representatives of the Caribbe element, which is characterized by species shared by North and South America and whose presence in the Antilles is interpreted as migration from the continental masses [114]. To confirm any of the patterns described above, more records are needed to test these hypotheses.
Based on our exhaustive literature search, 14 new records were found, mainly from sources such as scientific outreach journals, theses, conference proceedings, and posters, which could be considered "less rigorous", since these records are not subjected to the same scrutiny by specialized taxonomists, as those records are published in a scientific journals; therefore, it is necessary to perform rigorous and appropriated studies to confirm these 14 additional records.
Furthermore, the knowledge of Mexican marine tardigrades is very recent. The first record was obtained in 2019, in fact, of the total of analyzed works, 12% corresponds to the study of marine tardigrades. This is a global phenomenon, as marine tardigrades are less well known than limnoterrestrial ones [50], which comprise approximately 15% of the tardigrade fauna [21]. In Mexico, marine tardigrades have been collected at depths ranging from 0.2 m in the littoral zone to 2847 m depth in the bathypelagic zone [38], while in the rest of the world they have been found up to 4170 m [115]; therefore, in Mexico, there remain coast and deep ocean areas missing for exploration, such as the abyssal zones present in the Gulf of Mexico, Caribbean and Atlantic Ocean [116].
One of the most interesting characteristics of tardigrades is their adaptation ability to hostile environments; these traits have supported the idea that tardigrade species have been considered largely ubiquitous and cosmopolitan. Early alpha taxonomy studies using morphological characteristics as the main form of identification have led to the formation of species groups or complexes, e.g., Echiniscus arctomys, Macrobiotus hufelandi, Milnesium tardigradum, Paramacrobiotus richtersi, etc. However, results from recent integra-tive studies indicate that tardigrade species may not be as widely distributed as previously thought [117][118][119][120][121][122][123]. This idea is supported by distribution patterns reported in other metazoans [124][125][126][127][128][129].
Estimating richness. In 2014, a series of 29 publications were made to analyze the Mexican biodiversity, in which 56 biological groups quantifying the number of species occurring in the country, as well-developed predictions on how many more species remained to be found [24]. Our calculation of potential diversity with the ratio of known to unknown species indicate that 6.1% of the global diversity of the 15 eukaryotic groups analyzed are found in Mexico, a very similar value to that found by Martínez-Meyer et al. [24] of 8.5% (although the biological groups analyzed in both studies are not the same). Martinez-Meyer et al. [24] analyzed microinvertebrate groups, such as springtails [130], proturans [77], and diplurans [78], leaving aside some groups of invertebrates registered in the country (tardigrades among them they). Thus, until now, there was no formal calculation of the current and estimated richness of the Tardigrada phylum in the country.
As in all richness estimation methods, the ratio of known to unknown species method, used in the present study, has its drawbacks. In the case of tardigrades, a very wide range of richness value was obtained because groups with very little presence in Mexico were used (proturans, copepods, onychophorans, and others highly abundant, such as cladocerans). These groups included, as tardigrades are more similar in habitats and lifestyles, to Protura, Diplura, Rotifera, free-living nematodes, Collembola, and mites; therefore, more real richness ranges were expected. Another important point is that extrapolating diversity from one taxon to another may be considered guesswork.
The accumulation curve method determined that for both genera and species the inventories are incomplete because 22% of the genera and 61% of the species have yet to be discovered; however, estimates of the number of species derived in this way often lack associated margins of error, making it impossible to objectively assess their accuracy [131,132]. Moreover, for genera level, the inventory is reliable and for species it is unreliable. In other studies of invertebrates [133][134][135], the reliability of the inventories depends on the study group, their vagility, and habitat characteristics. An example of this are the families Araneidae, Thomisidae, and Salticidae (Arachnida); in the first two families, it is possible to obtain more complete inventories in plots of 1 km 2 [136], while in the family Salticidae the inventories obtained in an area of equal size are far from complete, due to the abundance of rare species [89]. In the case of the tardigrades, it is known that the species of the genus Macrobiotus are found in a wide range of habitats and hosts, being one of the most abundant genera in bryophytes [137,138], while other taxa seem to be restricted to a range of habitats and hosts. For instance, in the Iztaccíhuatl volcano in Mexico, the genera Adropion and Pilatobius, are apparently restricted to the Abies and Pinus forest, and the Hypsibiidae family are only present in the Mixed forest, Abies religiosa and Pinus hartwegii forest [47]. Therefore, to get complete the inventory, the sampling effort must be directed to a wide spectrum of environmental conditions, in addition to being directed to fulfill the objectives of the study.
Taxonomic effort. In the period from 2006 through to 2021, ten tardigrade species were described in Mexico, which correspond to 0.66 species per year. The descriptions of these species are found in seven publications and six of them were made by Mexican authors [25,[45][46][47][48][49], meaning that Mexico contributes to 0.44% of the authors and 0.43% of the tardigrades known so far. This situation can be observed in other countries, for example, in Colombia, a mega-diverse country, from 2014 to 2020 nine tardigrade species have been described, which corresponds to 0.67% of the global tardigrade diversity, contributing to 1.25% of the species authors [139][140][141][142]. Unfortunately, in other countries within the region, such as Trinidad and Tobago, Belize, Guatemala, El Salvador, and Honduras, among many others, tardigrade species remain undiscovered.
The number of tardigrade taxonomists in Mexico and Colombia is an example that in the last 10 years, more taxonomists have joined the study of tardigrades in their countries of origin. Although fortunate, the study of tardigrades from this geographic region is in the alpha phase, which implies intensive collecting, detailed observation, description, and species naming [143]. Therefore, continuity in the training of specialized taxonomists involved in the study of tardigrades is imperative, as well as the design of taxonomic work within a comprehensive research project.
Worldwide, the diversity of tardigrade species described by 10 prolific authors ranges between 50 and 204 described species per author (Pilato 204 species, Binda 111, Lisi 92, Kaczmarek 84, Michalczyk 74, Maucci 58, Murray 58, Mihelčič 57, Iharos 56 and Dastych 54 species), and in total, these authors have described 848 species (63.37%) of those known to date. Such a scenario in which prolific authors describe the bulk of species can be observed in numerous groups. A clear example is in bees, on which almost 19,508 known species in the world, and the 10 most prolific authors have described from 417 to 3394 species, and in total, these authors have described 9435 species (48.36%); the author who described the most species (3394 species) contributed more than 17% of the species so far known [144]. Comparing the number of authors who describe bees and tardigrades, it is evident that there are fewer descriptions in Tardigrada; however, the percentage described by prolific authors is considerably higher.
Finally, it is important to mention that another factor that has maintained a low rate of species description is the lack of inclusion of methods of analysis for their delimitation, particularly in taxonomically complex groups, which is reflected in the fact that 25.7% of the species recorded in Mexico are identified as "confer" (cf.); the inclusion of different methods of analysis of morphological variation, character observation techniques, and molecular data analysis in an integrative context will allow the recognition, delimitation, and generation of more robust hypotheses of the documented species.

Conclusions
Records of the tardigrade fauna in Mexico have been based on morphological characters such as size and shape of the claws, organization of the buccopharyngeal apparatus, cuticular patterns, and egg morphology. However, these structures present intra-specific variation that has not been quantified with more precise methods; it is, therefore, necessary to use integrative methods using different sources of information, such as ecological data, molecular data from nuclear and mitochondrial DNA, and morphological characters. These would be invaluable to future studies. From the information obtained by a literature review, 105 taxa were recorded, belonging to 42 genera and 75 species. According to the bibliographic source, 31.85% of the records were obtained from indexed journals, 65.6% from theses, 1.89% from conference proceedings, and 0.63% from scientific outreach journals, so many of the records reported need confirmation. The best-sampled provinces in the country were TVBP, SMOrP, and SMOcP and the least explored were SMSP, PLP, and BBP, although in all provinces, records are concentrated in one or two states, leaving much of the territory to be explored. Moreover, two biogeographic provinces and six marine ecoregions remain unexplored, which in turn should stimulate future research and Tardigrada studies in the country. Another factor to evaluate in terms of tardigrade diversity is their ocean depth distribution since in Mexico there are still coast and deeper areas to be explored, such as the abyssal zones that are present in the Gulf of Mexico, the Caribbean, and the Atlantic Ocean.
After this revision, the known distribution of six taxa was expanded as sampling efforts intensified in Mexico and other regions of the world. This pattern of spread distribution will become increasingly common. Most of the species found in the country have a wide distribution, except for species that have been recorded only from the American continent. Likewise, many species (14) have only been recorded in Mexico, so it can be said that they are temporarily limited to the type locality until new records are found that confirm their presence in other places.
For the actual detection of missing species in future studies, it is necessary to carry out systematic sampling in different types of vegetation and habitats in poorly explored areas, as well as to apply integrative methods that facilitate the delimitation of species. We also want to attract the attention of scholars of this group and encourage them to publish their results. In Mexico, this phylum requires a lot of studies, and any contribution adds to the knowledge of the group; in addition, the exercise that involves the presentation of a manuscript for its publication adds to the taxonomic and writing experience, which is essential in the formation of specialists in the group.