Preliminary Study on Species Diversity and Community Characteristics of Gamasid Mites on Small Mammals in Three Parallel Rivers Area of China

Simple Summary In addition to causing dermatitis, some gamasid mites are also the vector or potential vector of some zoonotic diseases. Located in the northwest Yunnan Province of southwest China, the Three Parallel Rivers Area of China is one of the hotspots of biodiversity research in the world. Based on the previous field survey in the Three Parallel Rivers Area from 2001 to 2015, this paper reports the species diversity and basic community characteristics of gamasid mites on small mammals in this unique geographical area for the first time. From the body surface of 3830 small mammal hosts, 26,048 gamasid mites were collected and identified as 10 families, 21 genera and 82 species (excluding 847 unidentified specimens) with high species diversity. The species diversity of the gamasid mite community fluctuated greatly in different elevation gradients. The highest peaks of species richness and β diversity appeared at altitudes of 3000–3500 m (S = 42) and 1500–2000 m (β = 17.5), respectively. The species abundance distribution of the mites was successfully fitted by Preston’s lognormal model, and the total number of gamasid mite species in the Three Parallel Rivers Area was estimated to be 153 species. Abstract (1) Background: Gamasid mites are a large group of arthropods, and some of them are of medical importance. Besides directly biting humans and causing dermatitis, some gamasid mites are the vector of rickettsialpox and potential vector of hemorrhagic fever with renal syndrome (HFRS). The Three Parallel Rivers Area of China is one of the hotspots of biodiversity research in the world, with complicated topographic landforms, different types of vegetation, special elevation gradients and high biodiversity. (2) Methods: Species richness (S): the Shannon–Wiener diversity index (H), Simpson dominance index (D) and Pielou evenness index (E) were used to analyze the basic community structure. The β diversity (Cody index) was used to reflect the diversity difference between any two adjacent elevation gradients. The method based on Preston’s lognormal model for species abundance distribution was used to estimate the total number of gamasid mite species. (3) Results: A total of 3830 small mammal hosts captured from the nine survey sites were identified as 44 species, 27 genera and nine families in five orders. Apodemus chevrieri, Eothenomys miletus and A. draco were the dominant host species with a total constituent ratio Cr = 52.037%. From the body surface of the hosts, 26,048 gamasid mites were collected and identified as 10 families, 21 genera and 82 species (excluding 847 unidentified specimens) with high species richness (S = 82) and diversity (H = 2.33). The three dominant mite species were Dipolaelaps anourosorecis, Laelaps nuttalli and L. echidninus, with a total Cr = 64.46% (16,791/26,048). There are significant differences in the species composition, species diversity and dominant species of gamasid mites on different hosts. The species diversity of the mite community fluctuated greatly in different elevation gradients. The highest peaks of species richness and β diversity appeared at altitudes of 3000–3500 m (S = 42) and 1500–2000 m (β = 17.5), respectively. The species abundance distribution of the mites was successfully fitted by Preston’s lognormal model with S^(R)=19e−[0.22(R−0)]2 (α = 0.22, R2 = 0.9879). Based on fitting the theoretical curve by Preston’s model, the total number of gamasid mite species was estimated to be 153 species. (4) Conclusions: Gamasid mites on small mammals are abundant with complex community structures and high species diversity in the Three Parallel Rivers Area of China. There is an apparent community heterogeneity of the mites on different hosts and in different environments.


Introduction
Gamasid mites are a large group of arthropods with complex and diverse ecological behaviors. It is estimated that about 8280 species of gamasid mites have been recorded worldwide, with more than 600 species documented in China [1,2]. Gamasid mites belong to the suborder Mesostigmata and order Parasitiformes in the subclass Acari of the class Arachnoidea [3,4]. There are five primary stages in the life cycle of gamasid mites, the egg, larva, protonymph (first nymph), deutonymph (second nymph) and adult (the male and female). The majority of gamasid mites are soil-living and free-living with complex habitats and ecological behaviors. For ectoparasitic gamasid mites, nearly all stages (except the egg) can be found on the body surface of small mammals and other animal hosts which are used as food sources of the mites [5][6][7]. Some species of ectoparasitic gamasid mites can transmit some zoonoses, and some endoparasitic gamasid mites can directly parasitize the lungs of animals or people, causing pulmonary acariasis [8][9][10]. As a whole group of arthropods, ectoparasitic gamasid mites (including facultatively and occasionally ectoparasitic species) have extensive hosts, including small mammals, aves (birds), reptiles, amphibians and other arthropods. Besides directly stinging the human skin and causing dermatitis, some species of gamasid mites are suspected to be associated with the transmission of more than 20 zoonotic diseases, with some species being the transmitting vector of rickettsialpox and potential vectors of hemorrhagic fever with renal syndrome (HFRS) [11][12][13][14][15][16]. Some gamasid mites can be the intermediate host of cotton rat filarial (Litomosoides carinii Travassos, 1919), the pathogen of cotton rat filarial disease, causing harm to experimental rats and mice [17,18].
Located in the Longitudinal Range-Gorge Region of the Hengduan Mountains in the eastern Himalayas, the Three Parallel Rivers Area is a unique geographical region in the northwest of Yunnan Province, southwest China. It is a transitional zone between the Tibetan Plateau (Qinghai-Tibet Plateau) and Yunnan-Guizhou Plateau [19][20][21][22]. As the intersection of the three geographical regions of East Asia, South Asia and the Tibetan Plateau, the Three Parallel Rivers Area borders Myanmar to the west, Tibet Autonomous Region to the northwest and Sichuan Province to the northeast and east [23,24]. In the northwest of Yunnan Province, there are three rivers (Jinshajiang, Lancangjiang and Nujiang rivers) originating from the Tibetan Plateau, flowing from the north towards the south [19,24]. The three rivers parallel each other, passing through Gaoligong Mountain, Nushan Mountain, Yunling Mountain and Shaluli Mountain without intersecting [23,25]. This area has a series of deep valleys, and four central high mountains are divided by three rivers, forming a very complex longitudinal mountain-valley landscape [24,25]. The topography, landscape, ecological environment and climate type in the area are incredibly complex with considerable changes in altitudes ranging from 720 m to 6740 m, which may contribute to the complex composition and high species diversity of animals and plants. It is one of the areas with the world's most diverse species and has become one of the hotspots in biodiversity research [26,27]. Ma et al. conducted a particular study on plant species diversity in the Three Parallel Rivers Area and found 98 species of plants [22], Renee Mullen et al. conducted a regional-level assessment of the alpine flora and reported 369 species of plants [28].
Rodents and some other small mammals are not only the important infectious source of some zoonoses (zoonotic diseases), but are also the most important hosts of ectoparasitic gamasid mites [2,4,29], and therefore it is of significance to study the mites on their body surface. Some previous studies have shown that the Three Parallel Rivers Area and its adjacent areas are rich in species of aves and mammals (including small mammals), but unfortunately few studies have dealt with ectoparasitic mites (including ectoparasitic gamasid mites) in the areas [30][31][32][33]. According to the field investigation data of nine survey sites in the Three Parallel Rivers Area from 2001 to 2015, this paper conducted a preliminary study on the species diversity and basic community characteristics of gamasid mites on small mammals in this unique area for the first time.

Field Survey Sites
The original data came from a long-term field investigation of nine survey sites in the Three Parallel Rivers Area of northwest

Gamasid Mites Collection and Identification
In the involved nine survey sites, mousetraps were placed randomly in different habitats (residential areas, cultivated farmlands, grasslands, shrubs, forests, etc.) in the evening and checked the following morning. The trapped rodents and other small mammals (hosts) were separately placed in marked white cloth bags and then transported to the

Gamasid Mites Collection and Identification
In the involved nine survey sites, mousetraps were placed randomly in different habitats (residential areas, cultivated farmlands, grasslands, shrubs, forests, etc.) in the evening and checked the following morning. The trapped rodents and other small mammals (hosts) were separately placed in marked white cloth bags and then transported to the field laboratory where gamasid mites were collected. The captured hosts were anesthetized with ether and then gamasid mites on their body surface were collected one by one. According to the appearance (body sizes, shape and coat color), various comprehensive measurements (body length, body weight, tail length, ear height and hind foot length, etc.) and other morphological features, each host was identified into species [34][35][36]. The collected gamasid mites from each host were preserved in a vial containing 70% ethanol and then mounted onto a glass slide in Hoyer's mounting medium. The mounted mites were finally identified into species under a microscope based on related taxonomic literature and identification books which contain a series of identification keys [4,[37][38][39]. The use of animals (including animal euthanasia) for our research was officially approved by the Animals' Ethics Committee of Dali University, under permission number DLDXLL2020-1104. Representative specimens of animal hosts and gamasid mites were deposited in the specimen repository of the Institute of Pathogens and Vectors, Dali University, Dali, Yunnan, China.

Species Diversity and Basic Community Structure Statistics
Together with species richness (S), which was used to reflect species diversity, the Shannon-Wiener diversity index (H), Simpson dominance index (D) and Pielou evenness index (E) were used to analyze the basic community structure of the small mammal hosts and gamasid mites [40,41]. The β diversity (Cody index) was used to reflect the diversity difference between any two adjacent elevation gradients [42,43].
In the above formulas, S represents species richness (the number of species), P i represents the proportion of the individuals of gamasid mite species i (or hosts species i) to the total number of all the species, g(h) is the number of species increased along the elevation gradient h (the number of species that were absent in the previous gradient and newly emerged in the following gradient), l(h) is the number of species decreased along the elevation gradient h (the number of species present in the previous gradient but absent in the following gradient).

Species Abundance Distribution and Total Species Estimation
All gamasid mites from nine survey sites in the Three Parallel Rivers Area were taken as a gamasid mite community unit. To illustrate the species abundance distribution of the gamasid mite community on small mammals, a semi-logarithmic rectangular system was established. The X-axis was labeled with log intervals based on log 3 M, which indicates the individuals of gamasid mites, and the Y-axis was marked with arithmetic scales, representing the number of gamasid mite species. Based on the following formulae, Preston's lognormal model was used to fit the theoretical curve of species abundance distribution, and the fitting goodness (R 2 ) was calculated. According to the final fitting goodness (R 2 ), the species abundance theoretical curve equation was established and the theoretical curve was drawn [44,45]. The total number of gamasid mite species (S T ) was then estimated. According to the total number of species (S T ), the likely missed number of species (S M ) in the field survey was approximately calculated [46,47].
In the above formulas, S(R) represents the theoretical number of gamasid species at the R-th log interval, S 0 represents the number of gamasid mite species at the R 0 log interval, m is the number of log intervals and R 0 represents the mode log interval. The value of α was determined according to the best fitting goodness, R 2 . S'(R) represents the actual number of gamasid mite species at R-th log interval and S(R) represents the average number of gamasid mite species for each log interval. S T is the total theoretical number of gamasid mite species in the community (the total expected number of gamasid mite species), S M is the number of gamasid mite species likely missed in the field survey and S A is the number of gamasid mite species actually collected in the field survey.  (Table 1). Of the captured small mammal hosts, rodents (order Rodentia) accounted for 75.00% of species (33/44) and 91.83% of individuals (3517/3830), which are the most dominant category of small mammal hosts in the Three Parallel Rivers Area, and the order Eulipotyphla came next (7.57%, 290/3830). The order Carnivora was the least, accounting for only 0.03% (1/3830) ( Table 1) Of the 26,895 collected gamasid mites, 26,048 identified as comprising 10 families (Laelapidae, Dermanyssidae, Macronyssidae, Aceocejidae, Ameroseiidae, Parasitidae, Parholaspidae, Macrochelidae, Pachylaelaptidae and Blattisocidae), 21 genera and 82 species (not including 847 unidentified individuals of gamasid mites due to structural damage, dirt cover (blurry structures) and some suspected new species). Of the 10 families, the family Laelapidae had the most abundant genera, species and individuals: 11 genera, 64 species and 24,853 individuals, which accounted for 95.41% of total gamasid mite individuals. The family Parholaspidae had the fewest genera, species and individuals (one species in one genus with only one individual), which accounted for less than 0.01% of total gamasid mite individuals ( Table 2). The 82 identified gamasid mite species and their corresponding small mammal hosts in nine survey sites were given in "Appendix B: Taxonomic checklist of gamasid mites identified from nine survey sites in Three Parallel Rivers Area of northwest Yunnan, China (2001-2015)".

Structural Characteristics and Species Diversity Changes of Gamasid Mite Communities
In this study, all the captured small mammal hosts and identified gamasid mites from the nine survey sites in the Three Parallel Rivers Area were taken as a host community unit and a gamasid mite community unit, respectively, then the species richness (S), Shannon-Wiener diversity index (H), Pielou evenness index (E) and Simpson dominance index (D) of the community structure were calculated. The identified 44 species and 3830 small mammal hosts constituted the "small mammal host community in the Three Parallel Rivers Area". The species richness (S), Shannon-Wiener diversity (H), Pielou evenness (E) and Simpson dominance (D) of the host community were S = 44, H = 2.52, E = 0.67 and D = 0.12, respectively. As shown in Table 3, three species of hosts, Apodemus chevrieri (Milne-Edwards, 1868), Eothenomys miletus (Thomas, 1914) and A. draco (Barrett-Hamiliton, 1900), were the most dominant in the Three Parallel Rivers Area, and their constituent ratios (Cr) were Cr = 24.28% (930/3830), Cr = 14.41% (552/3830) and Cr = 13.34% (511/3830), respectively. The individuals of A. chevrieri, E. miletus and A. draco exceeded 500, and the total constituent ratios of the three dominant hosts were 52.04% (1993/3830) of the whole. From A. chevrieri, 30 species and 1775 individuals of gamasid mites were collected, with 34 species and 877 individuals of gamasid mites from E. miletus, and 24 species and 558 individuals of gamasid mites from A. draco. The species richness was highest in E. miletus (S = 34), and Shannon-Wiener diversity was highest in A. chevrieri (H = 2.50) ( Table 3). The species richness of the gamasid mite community on different host species is different, and the species composition of the dominant mites on different hosts is also quite different (Tables 3 and 4).  The altitudes of the nine survey sites in the Three Parallel Rivers Area are between 862 m and 3880 m. According to the recorded data, seven altitude gradients were used to analyze the variations of species diversity of gamasid mites and small mammal hosts along different altitude gradients. The species richness of small mammal hosts showed a tendency of increasing gradually in the beginning and then decreasing abruptly from the low altitude (<1000 m) to the high altitude (>3500 m) along the vertical gradients ( Figure 2). And the species richness of small mammal hosts reached the highest value (S = 32) in the high-altitude regions between 3000 and 3500 m (Figure 2, Table 5). In addition, the species richness of gamasid mites formed the first peak (S = 41) in the low-altitude regions between 1500 and 2000 m and the second peak (S = 42) in the high-altitude regions between 3000 and 3500 m (Figure 2, Table 5). This showed that the changing trend of species richness of gamasid mites is not completely consistent with the small mammal hosts. Cody diversity index was used to analyze the beta-diversity (β diversity) of the gamasid mite community in the Three Parallel Rivers Area, and the β diversity of gamasid mites fluctuated from the low latitude to the high latitude of nine survey sites, which reached the highest value (β = 17.5) in the low-altitude regions between 1500 and 2000 m (Figure 3).

Species Abundance Distribution and Total Species Estimation of Gamasid Mite Community
Of the 82 species and 26,048 individuals of identified gamasid mites, the number of gamasid mite individuals in the seventh to eighth logarithmic scale was the highest, but the number of species was minimal. When the logarithmic scale is zero (R 0 = 0), the number of gamasid mites was only one individual, but the number of species was large (S 0 = 19). Based on Preston's lognormal model, the species abundance distribution of the (S 0 = 19, R 0 = 0). The species abundance distribution curve showed that the rare species of gamasid mites with few individuals accounted for the majority of species in the community, while the dominant species with a large number of individuals were only a few. With the increase of the number of gamasid mite individuals, the number of gamasid mite species gradually decreased (Table 6, Figure 4). Of the 82 gamasid mite species, 19 species had only one individual collected, 12 species had two individuals, and one species had three individuals. On the basis of the curve fitting, the total gamasid mite species in the Three Parallel Rivers Area was estimated to be 153 species (S T = 153), 71 species more than the actual collected 82 species. ber of gamasid mites was only one individual, but the number of species was larg 19). Based on Preston's lognormal model, the species abundance distribution gamasid mite community in the Three Parallel Rivers Area was successfully fitted lognormal distribution with α = 0.22 and R 2 = 0.9879. The theoretical curve eq was ̂( ) = 19 −[0.22( −0)] 2 (S0 = 19, R0 = 0). The species abundance distribution showed that the rare species of gamasid mites with few individuals accounted f majority of species in the community, while the dominant species with a large num individuals were only a few. With the increase of the number of gamasid mite indivi the number of gamasid mite species gradually decreased (Table 6, Figure 4). Of gamasid mite species, 19 species had only one individual collected, 12 species ha individuals, and one species had three individuals. On the basis of the curve fittin total gamasid mite species in the Three Parallel Rivers Area was estimated to be 15 cies (ST = 153), 71 species more than the actual collected 82 species.

Discussion
The present study showed a great species diversity of gamasid mites and their mammal hosts in the Three Parallel Rivers Area of China. From abundant species of mammal hosts (44 species with 3830 individuals), a total of 26,895 individuals of ga mites were collected, and 26,048 of them (26,048/26,895) are identified as compris

Discussion
The present study showed a great species diversity of gamasid mites and their small mammal hosts in the Three Parallel Rivers Area of China. From abundant species of small mammal hosts (44 species with 3830 individuals), a total of 26,895 individuals of gamasid mites were collected, and 26,048 of them (26,048/26,895) are identified as comprising 82 species and 21 genera in ten families. The Three Parallel Rivers Area is only a partial area in the northwestern Yunnan Province, southwest China, in which the identified species of gamasid mites (82 species) in the present study were not only much more than those from other provinces of China, but also vastly exceeded those recorded from some other regions in the world. The identified 82 species of gamasid mites on small mammals in the present study were obviously greater than the total recorded species of the mites (excluding soil-living and free-living species) in other provinces of China. For example, there have been 53 gamasid mites species recorded from Zhejiang Province, 46 species from Hebei Province, 78 species from Hubei Province, 21 species from Chongqing municipality (a provincial administrative region), 28 species from Shandong Province and 69 species from Fujian province, respectively, according to some previous literature [48][49][50][51][52][53]. In the world, 56 species of gamasid mites were recorded from two geomorphologic complexes of South-East Slovakia: Košická kotlina basin and Východoslovenská rovina plain, 71 gamasid mite species in Asiatic Russia and 37 species in the Kurzeme Coast of Latvia (Eastern Europe) [54][55][56]. The present study involved only the ectoparasitic gamasid mites on small mammals (including facultatively and occasionally ectoparasitic species). If soil-living and free-living gamasid mites are investigated in the future, the number of species would be greatly increased. The abundance of gamasid mite species in the Three Parallel Rivers Area is related to the unique geographical environment of this area. Located in the longitudinal valley of the Hengduan Mountains, the complex and diverse variations in topographies, landscapes, ecological habitats and climate types, with huge elevation differences in the Three Parallel Rivers Area support the rich biological resources and high biodiversity together [20,57,58]. The lowest altitude of the Three Parallel Rivers Area is 720 m in Nujiang river valley, and the highest altitude is 6740 m on the top of the Meli Snow Mountains [22,59]. There are many high mountains with an altitude of more than 4000 m in the Three Parallel Rivers Area, such as Laojun Mountain, Gaoligong Mountain, Yulong Snow Mountain, Haba Snow Mountain and Meili Snow Mountain [25]. These mountain landforms usually have different vertical forest zones and climatic types, with high landscape heterogeneity and biodiversity, and the species composition of plants and animals are quite different from the bottom to the top of the mountains [20,57,58]. The species diversity of ectoparasitic gamasid mites is closely related to that of their small mammal hosts. The above complex geographical, topographical and ecological environments led to the abundant species of small mammal hosts (44 species) and gamasid mites (82 species) in the present study. The unique geographical, topographical and ecological characteristics should be the initial factors influencing the high species diversity of gamasid mite species in the Three Parallel Rivers Area.
Due to the limitations of transportation, human resources and financial funds, the nine survey sites (elevations between 862 m and 3880 m) in this study have not covered all the mountain landforms in the Three Parallel Rivers Area. The alpine zone above 4000 m has not been involved, and the number of survey sites is also limited (only nine sites). This study is only a preliminary report on gamasid species diversity and community characteristics in the Three Parallel Rivers Area of China. With the continuous expansion and deepening of field investigation in the future, more small mammal hosts and gamasid mite species will be collected, and the composition of dominant hosts and gamasid mite species may also change to some extent. The results of this study suggest that there are abundant species of small mammal hosts and gamasid mites with a complex community structure and high species diversity in the Three Parallel Rivers Area of China. Moreover, the results show that further research is of great academic value.
The variation of species diversity along elevation gradients has been one crucial topic in ecology and biodiversity research in recent years [60][61][62]. The variation patterns of species diversity along elevation gradients are often diverse in different geographical regions because of habitat heterogeneity, vegetation difference, climate instability and some other environmental factors [62][63][64]. Some previous studies showed that the species diversity of lizards decreases with increasing elevation and it peaks at middle elevations [65]. The highest species diversity of chigger mites also occurred in the intermediate altitude regions (2000-2500 m and 2000-3000 m) [41,66,67]. In this study, species richness and β diversity are used to analyze the variation of gamasid mite species diversity along altitudinal gradients. The results of this study showed that the species diversity of the gamasid mite community in the Three Parallel Rivers Area fluctuated significantly in different elevation gradients. The peak of species richness (S = 42) occurred at the altitude between 3000 and 3500 m, and the peak of β diversity (β = 17.5) occurred at the altitude between 1500 and 2000 m (Figures 2 and 3). The result in this paper is not completely consistent with the previous studies [41] and we need further studies.
In this paper, the species composition, species diversity and dominant species of gamasid mites on different host species were also significantly different (Tables 3 and 4), which reflects the heterogeneity of gamasid mite communities on different hosts. The species diversity variation and community heterogeneity of gamasid mites on different host species may be associated with different biological characteristics of different host species and the low host specificity of gamasid mites. Different species of animal hosts with different biological characteristics usually have different susceptibility to parasite infections (including ectoparasite infestations) with different parasite burdens, species composition and species diversity [40,[67][68][69]. Many ectoparasitic gamasid mites have low host specificity, and cross infestation of the mites among different host species is very common, which may also lead to the heterogeneity of a gamasid mite community on small mammals [67,70,71]. The species diversity variation and community heterogeneity of gamasid mites on different host species further suggests the extraordinary complexity of gamasid mite communities on small mammals in the Three Parallel Rivers Area of China.
Of the 82 identified gamasid mite species in this study, several species are closely related to medicine. Ornithonyssus bacoti (Hirst, 1913) often stings humans to cause dermatitis [5,7,72], and it is an effective vector of rickettsialpox and an intermediate host of cotton rat filaria, L. carinii [17,18,72]. Moreover, O. bacoti has been proved to be an important potential vector of hantavirus, the causative agent of hemorrhagic fever with renal syndrome, HFRS [2,[11][12][13]. Aside from O. bacoti, some other species of gamasid mites found in the present study, e.g., Haemolaelaps glasgowi (Ewing, 1925), Tricholaelaps myonysognathus (Grochovskaya et Nguen-Xuan-Xoe, 1961), Eulaelaps stabularis Koch, 1836, Eulaelaps shanghaiensis Wen, 1976 and L. echidninus can also invade humans to cause dermatitis, and they can be the potential vectors of hantavirus as well [3,14,73]. The occurrence of the above mite species in the Three Parallel Rivers Area may increase the potential danger to the people's health in this area.
Species abundance distribution, together with the fitting of its theoretical curve, is an important issue in community ecology, and it describes the relationship between the species and individuals in a specific community, that is, the abundance distribution of different species in the community. The species abundance distribution can directly reflect the proportion structure of the dominant, common and rare species in a community [44,74]. Based on the theoretical curve fitting of species abundance distribution, the expected total number of species (the theoretical total number of species) can be roughly estimated [45,75,76]. The species abundance distribution of gamasid mites in the present study was successfully fitted by Preston's lognormal distribution model, suggesting that it conforms to the lognormal distribution. The result showed that the expected number of gamasid mite species showed a gradual descending tendency with the increase of gamasid mite individuals. It reveals that the majority of the gamasid mite species are sporadic with few individuals, and few gamasid mite species are dominant with abundant individuals. The tendency of species abundance distribution in this paper is almost consistent with that in other ectoparasite communities, such as chigger mites in the previous papers [71,77].
How to scientifically estimate the expected total number of species (the theoretical total number of species) in a specific community is also an important issue in ecology and several methods have been previously used to estimate the total number of species [40,77,78]. As one estimation method, the method used in the present study is based on the curve fitting of species abundance, and it has been applied to study the community of chigger mites and some other ectoparasites, such as fleas and sucking lice [44,47]. The result showed that the theoretical total number of gamasid mite species in the Three Parallel Rivers Area is 153 (S T = 153), which greatly exceeded the 82 actually identified species (S A = 82). It reveals the fact that there are still many gamasid mite species (S M = 71) that have not been collected and which were missed in the actual field investigation. The results further suggest that gamasid mites on small mammals in the Three Parallel Rivers Area of China are abundant with complex community structures and high species diversity, which is of great academic value for further in-depth study.

Conclusions
Gamasid mites on small mammals are abundant with complex community structures and high species diversity in the Three Parallel Rivers Area of China. There is an apparent community heterogeneity of the mites on different hosts and in different environments. Institutional Review Board Statement: The use of animals (including animal euthanasia) for our research was officially approved by the Animals' Ethics Committee of Dali University, approval code: DLDXLL2020-1104, approval date: 4 November 2020.

Data Availability Statement:
The experimental data used to support the findings of this study are available from the corresponding author request.