{'fraction_non_alphanumeric': 0.036765671367420365, 'fraction_numerical': 0.02271593266629901, 'mean_word_length': 5.478394011568561, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 2, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Effects of latitude, host body size, and host trophic guild on patterns of diversity of helminths associated with humans, wild and domestic mammals of Mexico

Parasites are strictly associated with their hosts and present a great diversity of life histories, often resulting in different diversity patterns than those observed in free-living species. However, ecological approaches have detected that, in some cases, mammal-associated helminths respond similarly to non-parasitic species in terms of diversity patterns. Using 2200 recorded interactions, we analysed the diversity patterns of helminths (Acanthocephala, Nematoda and Platyhelminthes) harbored by humans, wild and domestic mammals of Mexico, depending on latitude, host body mass and trophic guild (carnivore, herbivore, insectivore, omnivore), considering helminth richness and average taxonomic distinctness, and host phylogenetic diversity and phylogenetic clustering. Latitude was positively correlated with the average taxonomic distinctness encompassing the three parasite phyla and nematodes. Northern latitudes had less taxonomically related parasite assemblages. Host body mass had a significant negative relationship with the average taxonomic distinctness of acanthocephalans and the richness of helminths associated to wild hosts. The omnivore hosts had greater parasite richness, while insectivores had a less taxonomically related parasite assemblage and herbivores had a more heterogeneous parasite assemblage. Our results highlight the importance of incorporating different dimensions of diversity, such as average taxonomic distinctness and to consider the composition of parasite assemblages to better understand their diversity patterns.


Introduction
Parasites are not distributed homogeneously in space, in time, or in their hosts and despite representing a significant portion of Earth's biodiversity, key determinants influencing parasite diversity are still not well understood. Macroecological approximations have been used to attempt to explain how patterns of parasite diversity are related to ecological, evolutionary, historical, and geographical processes (Dallas et al., 2018(Dallas et al., , 2019Kamiya et al., 2014a;Huang et al., 2015), and in some cases these studies have revealed that parasite diversity follows patterns similar to those found in non-parasitic wild species (Morand, 2015;Stephens et al., 2016).
Particularly in mammal-helminth associations, the geographical range size and population density of mammals have been positively related to the richness of parasites (Torres et al., 2006). At the same time, the latitudinal diversity gradient (LDG), a well-documented pattern of decreasing diversity with increasing latitude in several taxonomic groups, seems to be reversed (Lindenfors et al., 2007;Kamiya et al., 2014b) or nonexistent (Nunn et al., 2005;Bordes et al., 2010) when investigated in mammal-associated helminths. Host body mass also seems to influence parasite richness. Larger hosts are expected to harbor more parasites because they represent a larger "patch" of habitat, with which various parasites can be associated, and tend to have a greater probability of parasite exposure due to a higher resource consumption and a more extensive habitat range (Poulin, 1995;Ezenwa et al., 2006;Huang et al., 2015). Several authors have found a positive relationship between mammal body sizes and helminth richness (Vitone et al., 2004;Ezenwa et al., 2006;Lindefors et al., 2007;Kamiya et al., 2014b;Huang et al., 2015), although others have not (Morand and Poulin, 1998;Nunn et al., 2003;Dallas et al., 2018;Preisser, 2019).
Host diet also could play a key role in shaping helminth diversity and communities (Monello and Gompper, 2011;Locke et al., 2014;Gutiérrez et al., 2017). This is expected since many helminths require trophic transmission, and are acquired through the ingestion of infective larval stages, either directly or incidentally from the environment or from infected prey serving as intermediate or paratenic hosts (Lafferty, 1999;Poulin and Morand, 2004). Furthermore, species with a broader diet breadth are expected to harbor more parasites than hosts with a more specialized diet (Gutiérrez et al., 2019;Leung and Koprivnikar, 2019). In primate hosts diet has been found to be positively related with helminth and nematode richness (Vitone et al., 2004). This trait has also been associated to greater richness of cestodes in cricetid rodents (Preisser, 2019). However, Poulin (1995) and Nunn et al. (2003) reported no significant relationship between host diet and helminth richness in mammals.
Aspects of study design could obscure or modify the influence of predictors like latitude, host body mass and diet on parasite diversity. The geographical scale of the study (regional, continental, global) can influence the effect the predictors have on parasite richness (Kamiya et al., 2014b),similarly, lack of control in sampling effort and host phylogenetic relationships, might generate inflated estimates or mask these effects (Morand and Poulin, 1998;Nunn et al., 2003;Vitone et al., 2004;Krasnov et al., 2008b;Garrido-Olvera et al., 2012). Parasites are subject to environmental restrictions, the variation and availability of definitive or intermediate hosts across sites, which could extend or restrict their spatial distribution (Bozick and Real, 2015). Futhermore, other host traits like longevity, host density and social group behavior can also be determinants in the variation of parasite richness among hosts and act as covariables; making difficult to discern the true effect on the relationships between these predictors and parasite richness, as well as, producing different diversity patterns than those of non-parasitic organisms (Vitone et al., 2004;Nunn et al., 2003;Cooper et al., 2012;Violante-González et al., 2012).
A further potential limitation in studies of parasite diversity is the fact that studies on species diversity usually use species richness (S) as a measure. This measure is limiting because it considers all species to be equally distinct, and does not allow the incorporation of data about their ecological and evolutionary identity, or attempt to understand the structure of species assemblages more broadly (Purvis and Hector, 2000;Marrugan, 2004;Poulin, 2004;Gómez-Ortiz et al., 2017). For these reasons, the use of measures such as average taxonomic distinctness, phylogenetic and functional diversity, is currently increasing, providing information about processes or patterns that may have been overlooked (Krasnov et al., 2008a(Krasnov et al., , 2019Poulin and Leung, 2011;Rasmussen and Randhawa, 2018;Fecchio et al., 2020).
Taxonomic distinctness of parasites has been used as a complementary measure to parasite richness in studies about determinants of parasite diversity and have been shown to provide insight on how parasite assemblages are formed (Luque et al., 2004;Randhawa and Poulin, 2010;Rasmussen and Randhawa, 2018). A narrow taxonomic range in a parasite assemblage could be associated with intra-host speciation or colonization by a limited taxonomic spectrum of parasites. In contrast, when parasites in an assemblage are completely unrelated, the assemblage was likely formed by repeated colonisations by a broad taxonomic range of parasites. This allows the identification of features influencing the total number of parasites harboured by a host species, and the variety of parasite taxa a host will has; thus providing different explanations of the determinants affecting parasite biodiversity (Luque et al., 2004;Quiroz-Martínez and Salgado-Maldonado, 2013;Rasmussen and Randhawa, 2018).
Dimensions of diversity like phylogenetic or functional diversity have been used to understand the ecological mechanisms behind the correlation between host and parasite species richness. The positive correlation between these variables is frequently assumed to emerge because communities with higher host diversity offer additional niches for different parasites. Halliday et al. (2019) proposed that this relationship depends on how host communities assemble, and that shifts in host phylogenetic diversity during community assemblage can alter the relationship between host and parasite richness. More distantly related hosts are less likely to share parasites and communities with higher host phylogenetic diversity may harbor more parasite species. Nevertheless, a higher host phylogenetic diversity can also impede parasite transmission reducing parasite richness and weaken the relationship between host and parasite richness (Halliday et al., 2019(Halliday et al., , 2020. In Mexico, 564 species of mammals have been recorded (Gómez-Ortiz et al., 2017), but less than 25% of these mammal species have been examined for helminths (García-Prieto et al., 2012). This percentage reflects the lack of knowledge on helminth richness associated with different groups of mammals. In addition, the close contact between humans, wildlife and domestic species generated by the change in land use, the movement of hosts for livestock activities, and the introduction of non-native species provide opportunities for parasitic helminths to associate and interact with new hosts. Wells et al., 2019). This set of factors could alter the biological patterns and processes of associations between parasite and host (Kutz et al., 2005;Kennedy, 2006;Wells et al., 2018;Fernandez, 2019). In this study, we aim to describe the richness, average taxonomic distinctness, host phylogenetic diversity and host phylogenetic clustering of helminths associated with humans, wild and domestic mammals considering latitude, host body mass and trophic guild.

Database construction
A database was constructed using the taxonomic checklist of helminths associated with wild mammals in Mexico published by García-Prieto et al. (2012). This database was updated by gathering records published up until 2017 using search tools including the ISI Web of Science, CAB Abstracts, Biological Abstracts and Google Scholar with the key words: "mammal*", "helminth" and "Mexico". Information on records of helminths in humans and domestic animals was obtained from a bibliographic review (Lamothe-Argumedo and García-Prieto, 1988;Quiroz Romero, 1987) and electronic resources generated by the Parasitology Department of the School of Veterinary Medicine at UNAM and other departments at the UNAM (Alcalá et al., 2008; https://dato sabiertosunam.mx/; http://oreon.dgbiblio.unam.mx/). Complementary to this review, we searched the ISI Web of Science (https://apps. webofknowledge.com) using the keywords "dog", "cat", "cattle", "goat", "sheep", "rabbit", "horse", and "donkey" combined with the words "helminth" and "Mexico".
Each record in the database represented the detection of a parasite in a mammal species and contained information on the hosts (species, family) and parasites (species, family, class, order, phylum), as well as their geographical location (state, locality, geographic coordinates), and their bibliographic reference. In the records where the latitude and longitude were missing, we used the geographic coordinates provided by Google maps after entering the name of the locality and the state described for that record.
Helminths were classified into three categories: associated with wildlife mammals (W), associated with domestic mammals (D), and associated with both (D-W), based on the type of host in which the helminth species has been recorded in the literature. Since humans and domestic animals live in close proximity of each other, inhabit almost all terrestrial environments and the subsets of parasites they share increase the longer the animal species has been domesticated Wells et al., 2018), we included Homo sapiens in the domestic animals category. Synanthropic rodents were also classified within this group, because despite not going through a process of direct domestication their frequent contact and long-time commensal relationships with humans, has selected traits that differentiate them from wild rodents, and have favored their adaptation to anthropogenic environments (Weissbrod et al., 2017;Geiger et al., 2018). The domestic mammals included were: cow (Bos taurus), dog (Canis familiaris), goat (Capra hircus), donkey (Equus asinus), horse (Equus caballus), cat (Felis catus), human (Homo sapiens), mouse (Mus musculus), rabbit (Oryctolagus cuniculus), sheep (Ovis aries), brown rat (Rattus norvegicus), black rat (Rattus rattus), and pig (Sus scrofa).

Correcting sampling bias
The inventory of mammal-associated helminths distributed in the Mexican territory is far from complete; the sampling effort has been spotty and heterogeneous among wildlife hosts and geographic areas. In general, more intensely studied host species like domestic animals would be expected to have more parasite species recorded in the literature simply due to sampling effort, generating potential bias in helminth records and interactions. To address this problem, following the values of average helminth richness per mammal species presented by Dobson et al. (2008), we selected mammals with a minimum of 5 species of helminths recorded. Then, to account for sampling bias, we followed the methodology of several authors (Poulin, 1995;Morand and Poulin, 1998;Vitone et al., 2004;Rasmussen and Randhawa, 2018) to adjust diversity measure values (i.e. richness, average taxonomic distinctness) by sampling effort (i.e. number of publications). We did a linear or simple quadratic regression of the diversity measure as a function of the sampling effort (Rasmussen and Randhawa, 2018) and if significant (P < 0.05), used the residuals of the regression in place of the original diversity measure to correct for the influence of sampling effort in further analysis. The regressions were calculated using the lm() function implemented in the stats package in R version 3.5.0 (R Core Team, 2018). The results of these regressions are shown in Table 1 of the Supplementary Material.

Diversity measures
Four main dependent variables related to patterns of helminth and host diversity were analysed with respect to latitude, host body mass, and trophic guild: richness (S), host phylogenetic diversity (PD), host phylogenetic clustering (PC), and average taxonomic distinctness (Δ+). For latitude, using the geographic coordinates of the records we created a new column in the database containing the degree of latitude value without considering the decimals (i.e. 14⁰, 15⁰) as a classificatory variable. From the Phylacine and MOM v4.1 databases (Faurby et al., 2018;Smith et al., 2003) we compiled data on body mass (g), and consumed percentage of plants, vertebrates and invertebrates for each mammal host.
According to Faurvy et al. (2018) the values on host body mass were obtained in several ways. Values obtained from Smith et al. (2003) where species body mass was compiled taking into account these considerations: a) For gender-specific estimates obtained by averaging male or female values reported from across different geographic localities, these were then averaged to obtain an overall species value, b) When no geographic information was provided but gender-specific values were available, these were averaged for an overall species body mass estimate, c) If both male and female masses were not available, male mass was used preferentially, averaging over geographic localities if such information was available, d) If only female body mass was available, it was used, averaging over geographic localities as available. Values obtained from Faurby and Svenning (2016) were compiled from separate data sources. Lastly, for species for which body mass was not recorded, the assumption of isometries for related similar-size species and phylogenetic imputation were used to assign the value (Faurby et al., 2018;Faurby and Svenning, 2016).
Likewise the percentage on the diet categories, were obtained from other data sources and phylogenetic imputation (Faurby et al., 2018). These percentages were used to create a column in which we classified mammal-hosts into four trophic guilds: herbivores (species consuming 100% of plants), carnivores (species consuming 100% of vertebrates), insectivores (species consuming 100% of invertebrates) and omnivores (species consuming two or more diet categories).
Parasite richness (S) was measured as the total number of parasite species. Host phylogenetic diversity (PD) was the total sum of the branch lengths between host species present within a grouping variable (i.e. parasite phylum, type of host with which the helminth is associated) (Faith, 1992), and was measured using the pd() function implemented in the picante package (Kembel et al., 2010) in R (Team, 2018). The function ses.pd() or standardized effect size of phylogenetic diversity (SES.PD) from the picante package (Kembel et al., 2010) was used to calculate the phylogenetic clustering (PC). We generated 999 null model randomizations using the null model "taxa.labels", where taxa labels of all taxa in the phylogeny are repeatedly shuffled in order to assess whether phylogenetic diversity is high or low for a group of hosts . Negative SES. PD values and low quantiles (P < 0.05) indicate that host species within the group are more related than expected by chance, while positive values and high quantiles (P > 0.95) indicate overdispersion or the co-occurrence of distantly related hosts within the group Walker et al., 2017). The calculation of host phylogenetic diversity and phylogenetic clustering were based on the phylogeny of mammals published by Faurvy et al. (2018).
The average taxonomic distinctness (Δ+), was calculated exclusively for helminths using the taxondive() function implemented in the vegan package in R (Oksanen et al., 2016). The average taxonomic distinctness (Δ+) arranges helminth species in a taxonomic hierarchy and calculates the average number of steps that two species within a group require to reach the same level in that hierarchy .
Step lengths are standardized so parasite species connected at the highest taxonomic level (Phylum) is set to be equal to 100, with six levels require to reach species, such that each steps accumulates a value equal to 16.66 . The helminth hierarchy was Table 1 Correlations between the latitude, parasite richness (S), host phylogenetic diversity (PD), host phylogenetic clustering (PC) and parasite average taxonomic distinctness (Δ+) (dependent variables were controlled for number of publications, before calculating the correlations).

Measure
Correlation pvalue  (2013) for acanthocephalans, Jones et al. (2005) for trematodes, Caira and Jensen (2017) for cestodes and Anderson et al. (1974Anderson et al. ( -1985, Gibbons (2010) and Hodda (2011) for nematodes. The taxonomic categories recorded were: Species, Genus, Family, Order, Class and Phylum. The taxonomic distinctness is associated with a variance (Λ+), which measures the heterogeneity in the sets of parasites found in a host, since, it is possible to have two groups of hosts with the same Δ+ value but which differ in their taxonomic range Quiroz-Martínez and Salgado-Maldonado, 2013).

Data analysis
The diversity patterns of parasites and hosts were explored using the previously mentioned measures, considering overall helminths, separating the groups of helminths at the taxonomic level of phylum (Acanthocephala, Nematoda and Platyhelminthes) and at the categories of the type of host with which the parasite was commonly associated (Domestic, Wildlife, Domestic-Wildlife). For latitude, we calculated the value of the dependent variables (parasite richness, host phylogenetic diversity, host phylogenetic clustering and average taxonomic distinctness) considering the subsets of species in the before mentioned groups of phyla and type of host with which the helminth associates at each latitude degree. The association between latitude and the dependent variables was measured using a Pearson or Spearman correlation (when the values of the dependent variable did not have a normal distribution) with the stats package (R Core Team, 2018).
For host body mass, we log-transformed the values to correct the skewed distribution of this variable (Garamszegi, 2014). Then, as phylogenetically close species tend to share common traits and are not independent observations (Gay et al., 2014;Poulin and Mouillot, 2004;Vitone et al., 2004). We computed a Phylogenetic Generalized Least Squares model (PGLS) with the caper package in R (Orme et al., 2018); where host body mass was the predictor for parasite richness and parasite average taxonomic distinctness. The phylogenetic information was included from the Phylacine mammal phylogenetic tree (Faurby et al., 2018). PGLS is an extension of general linear model which allows the incorporation of phylogenetic dependence in the error term. This term was conformed by a matrix of expected trait covariances calculated using the phylogeny and the maximum likelihood (ML) estimate of Pagel's lambda λ (Cooper et al., 2012;Garamszegi, 2014). Finally, to assess whether trophic guild of the host had a relationship with parasite richness and average taxonomic distinctness, we generated a PGLS using trophic guild as predictor variable and parasite richness and average taxonomic distinctness as response variables. As described in section 2.2, whenever there was a significant relationship between sampling effort and some diversity measures, we used residuals from that regression rather than raw diversity measures to account for sampling bias (Supplementary Material Table S2).

Results
The database contained 2403 records, consisting of 118 host species, 346 helminth species, ten acanthocephalans, 84 platyhelminths, and 252 nematodes. Of the 118 nominal host species, 41 complied with the minimum of five helminth species, leaving 2200 records with 250 helminths, ten acanthocephalans, 210 nematodes, and 71 platyhelminths. The mammal order Carnivora had the greatest richness of helminths, whereas Cingulata had the lowest parasite richness (Supplementary Material Table S1). The average number of parasites per host was 6.09, with the horse (Equus caballus) having the greatest parasite richness (40). Seventeen helminths were associated with domestic animals, 158 were associated with both domestic and wild mammals, and 116 were associated with wildlife only.

Discussion
The description of the factors involved in parasite diversity patterns is of paramount importance in explaining how their diversity is distributed and improves our knowledge of parasite-host interactions. Besides understanding how parasite diversity is related to different features, it is fundamental to comprehend why some host species harbor more parasite richness than others (Luque et al., 2004;Locke et al., 2014). Due to differential sampling effort directed towards them for their importance in public health and livestock activities, hosts orders containing domestic species like Carnivora and Cetartiodactyla had more parasite richness. Our analysis also showed that nematodes and platyhelminths had higher parasite richness and host phylogenetic diversity than acanthocephalans, understanding as phylogenetic diversity, the sum of branches length of each group host associated with these two phylas. Both platyhelminths and nematodes are abundant and speciose among vertebrate hosts; they are diverse in their biology and life histories and they have evolved parasitism independently multiple times (Viney, 2018;Weinstein and Kuris, 2016). Acanthocephala, on the other hand, is a small phylum that has undergone multiple radiations from aquatic to terrestrial habitats, presents conservatism of their intermediate host and is mainly harbored by fishes and birds (Kennedy, 2006;Near, 2002).
A low value of average taxonomic distinctness (Δ+) means that on average the parasites in a group (i.e. Phylum) are closely taxonomically related (Quiroz-Martínez and Salgado-Maldonado, 2013). Nematodes of the same clade frequently parasitize different hosts, and the presence of congeneric nematode species is common in mammals (i.e., flock species) (Kennedy and Bush, 1992;Poulin, 2004). Several platyhelminths families also show diversification within mammal groups (e.g Hymenolepididae, Anoplocephalidae) (Morand et al., 2005), which could result in the observed values of Δ+. Despite having the lowest parasite richness, the phylum Acanthocephala presented the largest value of average taxonomic distinctness and their variance (Λ+). High values of Λ+ are associated with different families giving origin to single species, resulting in more heterogeneous assemblages. Since all acanthocephalans require an intermediate host (potentially including paratenic hosts) (Kennedy, 2006), the presence and consumption of different intermediate hosts such as insects or amphipods could result in the heterogenous assemblages we see in this phylum (Kennedy, 2006;Spickett et al., 2017;Prati and Henriksen, 2020).
When considering host categories, helminths associated with wildlife (W) and both domestic and wild (D-W) hosts had a greater richness and host phylogenetic diversity than those that are associated only with domestic mammals (D). Helminths associated with wildlife included marsupials, bats and armadillos among their hosts. The inclusion of these mammal orders increases the branch length measured in host phylogenetic diversity, as these organisms are more distantly related that the other mammals in this group. As for the value of average taxonomic distinctness the incorporation of families of acanthocephalans like Plagiorynchidae and Oligacanthorhynchidae harbored in marsupials, primates and carnivorans, as well as the adding of trematode species from families like Anenterotrematidae and Phaneropsolidae for bats, Rhopalidae, and the cestode order Onchoproteocephalidea for  marsupials increased the average taxonomic distinctness, generating a wider taxonomic range for the helminths associated to wildlife. Within the helminths associated with both domestic and wild mammals, although most of the wildlife hosts were phylogenetically close to some domestic mammals, we also included bats, marsupials and the armadillo (Dasypus novemcinctus) which increased their host phylogenetic diversity. Parasite recorded in these hosts like Ancylostoma braziliensis, Paragonimus mexicanus and Plagiorchis verpertilionis besides being found in domestic animals have also been found in humans (Lamothe-Argumedo and García-Prieto, 1988;Guk et al., 2007;Guillot et al., 2012).
The average taxonomic distinctness in the D-W category showed taxonomic relatedness among the parasites, with one parasite order (Rhabditida) comprising almost 50% of the parasite species in that category and parasite families like Strongylidae and Trichostrongylidae contributing several species. Both families contained a variety of mammalian nematodes that live in the gastrointestinal system of artiodactyl and perissodactyl hosts (Chilton et al., 2006). Finally, helminths associated with domestic hosts (D) had the lowest host phylogenetic diversity and average taxonomic distinctness, of the seventeen helminth species classified in this group almost 50% were comprising within two orders (Cyclophyllidea and Spirurida). The cestodes comprised within Cyclophyllidea were associated to human and rodents, similarly Spirurida was represented by nematodes harbored by human, rodents and ruminants. Furthermore, parasites in both groups (Domestic, Domestic and Wildlife) are of human, veterinary and economic importance. (Smout et al., 2017;Wells et al., 2018).
Our study registered different patterns of helminth diversity across a geographical and biological predictors. The latitudinal diversity gradient (LDG) has been described for different species of mammals, insects, fishes, and bacteria (Willig et al., 2003), and a similar pattern of diversity could be expected in parasites, since they depend on their hosts for their survival; however, this is not always the case (Kennedy, 2006;Kamiya et al., 2014b;Morand, 2015;Preisser, 2019). We did find a negative correlation between latitude and host phylogenetic diversity (rs = − 0.47, p = 0.04), but we did not find a significant correlation between parasite richness and latitude. According to several authors, factors like temperature, precipitation and the capacity of some parasites to disperse across their host ranges could be involved in latitudinal patterns of helminth diversity rather than latitude per se, since these variables influence parasite stages occurring in the environment, the transmission to their hosts, the diversity and availability of intermediate hosts, resulting in different diversity patterns than those observed in their hosts (Lindenfors et al., 2007;Morand, 2015;Dallas et al., 2018;Preisser, 2019).
We did, however found a positive correlation between latitude and average taxonomic distinctness of helminths as a group and nematodes. Although assemblages at lower latitudes were richer in helminth species, they were congregated into families within the same order, while at higher latitudes the richness was lower, but parasites were more widely distributed among families and orders. The variation observed in the Δ+ of the parasite subsets may also be a result of the host species associated with them (Krasnov et al., 2008a); at lower latitudes hosts were mainly cetartiodactyls whereas, at higher latitudes species such as carnivorans and rodents were more frequent. However, is important to consider that parasites do not distribute completely or evenly across the geographic Table 2 Phylogenetic generalized least squares model (PGLS) for host body mass (predictor), parasite richness (S) and parasite average taxonomic distinctness (Δ+) (response variables were controlled for number of publications and the predictor variable was log-transformed before calculating PGLS).  range of their hosts distribution (Cumming, 1999;Shenbrot et al., 2010); factors related to environmental conditions, competition among parasites, presence and distribution of the intermediate hosts could influence parasite species composition (Cumming, 2002;Shenbrot et al., 2010;Krasnov et al., 2010;Warburton et al., 2016). Furthermore, even though we consider sampling effort bias, most of records we used came from latitudes of 18⁰, 19⁰ and 20⁰, future studies should revisit these relationships when more data of higher latitudes have been collected. Host body mass had a significant negative relationship with average taxonomic distinctness of acanthocephalans; meaning that as host body mass increases, parasite subsets for this phylum become more similar. Our corrected database contained ten acanthocephalan species distributed among eleven hosts, but only for four hosts it was possible to estimate the average taxonomic distinctness, since the remaining hosts had only one species of acanthocephalan (Clarke and Warwick, 2001); this relationship should be re-analised in future studies when more information on acanthocephalans in mammals is available, as small sample sizes could generate unreliable results (Garamszegi, 2014).
Contrary to the works of Ezenwa et al. (2006);Lindefors et al. (2007) and Kamiya et al. (2014b), we found a significant negative relationship between host body mass and the richness of helminth associated to wildlife mammals. Within this group more than a half of the hosts were rodents and bats. A study made by Dáttilo et al., in 2020 found a similar relationship between small-bodied mammals and parasite richness.The authors mentioned that differences in sampling effort, could lead to increase richness values, as small mammals like rodents and bats tended to be more well-sampled in comparison to other mammals. Host body mass has been noted as an important predictor of parasite richness. However, the strength of its effect may vary among taxonomic groups of hosts; for example, according to Kamiya et al. (2014b), the association between host body size and metazoan richness was significantly stronger in fishes than in birds and mammals. In addition, body size variation within host species can be linked to various processes that also influence parasite richness and their opportunity to find a host, such as host density, habitat use, and host life-history traits (e.g., longevity, terrestriality), such that host body size is not a universal covariate of parasite diversity (Dallas et al., 2018;Morand and Poulin, 1998).
The omnivore trophic guild presented the highest parasite richness, while parasites associated with insectivore hosts were less taxonomically related than the parasites in other trophic guilds, and herbivore hosts had the most heterogenous parasite species assemblages.   According to Luque et al. (2004), certain host traits may favor the acquisition of a wider taxonomic spectrum of parasites, whereas the presence of other traits may result in the same number of parasite species but with narrower taxonomic spectrum. Due to the variability of their diets, omnivores have a higher probability of encountering infecting stages of parasites than species whose diets only include a few elements (Miller et al., 2009;Albuquerque et al., 2016;Leung and Koprivnikar, 2019). However, the specificity of the intermediate host that some helminths present (Locke et al., 2010;Kennedy, 2006;Near, 2002) and the consumption of several of these type of hosts as occurs in insectivore hosts could result in less taxonomically related assemblages.
As for the variance of taxonomic distinctness, the high value presented in herbivores hosts suggests that one main branch in the taxonomic tree contributes proportionally more species than the others and provides information on particular taxa are over-or under-represented (Korallo et al., 2007;Milošević et al., 2012;Rasmussen and Randhawa, 2018). In the case of this trophic guild, two families (Trichostrongylidae, Strongylidae) accounted for nearly half of the parasites within this trophic guild. On average, herbivore hosts had higher parasite richness than insectivore hosts, but their taxonomic distinctness did not reflect this. This could happen because the route by which a host acquires the parasite could be limited to a type of resource (e.g. grass, leaves) and to several parasite species, but with a narrow taxonomic breadth (e.g. strongylids in horses, oxyurids in primates) (Vitone et al., 2004;Morariu et al., 2016). As for the rest of the guilds, in some cases trophically transmitted parasite families contributed several species (e.g. Phaneropsolidae for insectivores, Oligacanthorhynchidae for omnivores), while other families like (e.g. Ascaridae for carnivores or Molineidae for omnivores) also contributed various parasites but are not likely acquired through the diet. In the same way that other factors could be confounded with latitude and host body mass, future analysis should also consider other host traits like habitat use or lifestyle in addition to trophic guild, since they could have an important role in determining the composition of parasite assemblages (Morand et al., 2005;Koprivnikar, 2016, 2019). Lastly, parasite richness among our hosts was highly variable, ranging from five to 40 species, whereas average taxonomic distinctness values were much less variable, even for species with many parasites like domestic mammals (Supplementary Material Fig. 1). In addition, this measure did not seem to be affected by sampling effort (Korallo et al., 2007;Poulin and Mouillot, 2004;Rasmussen and Randhawa, 2018). These characteristics make average taxonomic distinctness an interesting complementary measure to consider in future works looking to detect diversity patterns that cannot be observed when we only consider the raw number of species (Quiroz-Martínez and Salgado-Maldonado, 2013;Rasmussen and Randhawa, 2018).

Conclusions
Although it was possible to describe some patterns already reported in the literature, the results of these approaches should be taken with caution. In many wild mammals, only one helminth species was recorded, and there are still many mammal species that have not been sampled at all. Among our records, taxonomic orders containing domestic animals had more recorded helminths than other host species. Domestic animals are more frequently studied because parasites affect human and animal health and the economics of livestock activities (Walker and Morgan, 2014). Host species richness and sampling effort are often associated with parasite richness (Morand, 2015). These information biases could be obscuring patterns produced by factors measured in this study.
It is therefore vital to continue sampling and identifying parasites associated with mammals, allowing a description of the multiple interactions that these organisms have with their hosts. As parasites, helminths are an essential and abundant component of natural communities. According to Dobson et al. (2008), there are approximately 11, 631 helminth species associated with the 4637 mammal species on the planet. As inventories of helminth species become more complete, it will be possible to more accurately identify the factors involved in the diversity patterns of these parasites.
In general, our analysis was enriched by the calculation of the average taxonomic distinctness, since this allowed us to explore the taxonomic composition of the parasite assemblages, facilitating the understanding of the diversity patterns we found in this work.
We consider that it is only through macroecological studies that address multiple dimensions of diversity that we can more accurately describe diversity patterns, the ecological traits involved in the composition of parasite assemblages in their hosts, and how current changes could alter them.

Declaration of competing interest
We warrant that none of the authors has any conflict of interest in regard to this manuscript.