Possible transmission of Sarcoptes scabiei between herbivorous Japanese serows and omnivorous Caniformia in Japan: a cryptic transmission and persistence?

Background Two transmission patterns of Sarcoptes scabiei in host mammal communities have been reported based on microsatellite-level genetic studies in the last two decades. While one involves restrictions among different host taxa, the other is associated with predator–prey interactions between different host taxa. In contrast to these observations, the present study reports a possible irregular case of transmission of S. scabiei between herbivorous Japanese serow and omnivorous Caniformia mammals in Japan, though under very weak predator–prey relationships. Methods DNA from 93 Sarcoptes mites isolated from omnivorous Caniformia (such as the domestic dog, raccoon dog, raccoon and Japanese marten), omnivorous Cetartiodactyla (wild boar) and herbivorous Cetartiodactyla (Japanese serow) in Japan were analyzed by amplifying nine microsatellite markers. Principal components analyses (PCA), Bayesian clustering analyses using STRUCTURE software, and phylogenetic analyses by constructing a NeighborNet network were applied to determine the genetic relationships among mites associated with host populations. Results In all the analyses, the genetic differentiation of Sarcoptes mites from wild boars and Japanese serows was observed. Conversely, considerably close genetic relationships were detected between Caniformia-derived and Japanese serow-derived mites. Because the predator–prey interactions between the omnivorous Caniformia and herbivorous Japanese serow are quite limited and epidemiological history shows at least a 10-year lag between the emergence of sarcoptic mange in Japanese serow and that in Caniformia, the transmission of S. scabiei from Caniformia to Japanese serow is highly suspected. Conclusions The close genetic relationships among mites beyond Host–taxon relationships and without obvious predator–prey interactions in Caniformia and Japanese serow deviate from previously reported S. scabiei transmission patterns. This type of cryptic relationship of S. scabiei populations may exist in local mammalian communities worldwide and become a risk factor for the conservation of the remnant and fragmented populations of wild mammals. Electronic supplementary material The online version of this article (10.1186/s13071-019-3630-5) contains supplementary material, which is available to authorized users.


Background
Sarcoptes scabiei is a notorious pathogenic mite species (Acari: Astigmata: Sarcoptidae) that causes a highly contagious dermatitis known as sarcoptic mange in mammalian hosts [1]. Among the mammalian species, and even within the species [2], the pathogenicity of the disease spans largely from devastating impacts on naive host individuals [1,3] to relatively low pathogenic infections in less-vulnerable host individuals [4]. The impact of sarcoptic mange on mammalian host populations changes among populations and in the trade-offs among biological factors such as population density, host mortality rate and host reproductive rate [5][6][7]. However, this has often been viewed as a possible threat for the conservation of remnant and isolated mammal populations [1]. Although emerging (or re-emerging) outbreaks of sarcoptic mange in mammalian populations have been reported worldwide [1,8], neither the origin of the epidemic lineages of S. scabiei nor the patterns of inter-species transmission of mites from the reservoir host species to the novel (or unobserved) host species have been sufficiently understood [9]. This is because the existence of a classical problem in the identification of mite strains, at least until the late 1990s, due to the nature of S. scabiei populations having different patterns of Host-specificities with indistinguishable morphology [8,10]. To resolve this problem, a method for genetic analysis using microsatellite (simple sequence repeat; SSR) markers was developed by Walton et al. [11][12][13] and, since then, this technique has been utilized as an efficient way to elucidate the hidden transmission web of Sarcoptes mites in the multi-host system [9,14].
Recent genetic studies using the SSR markers revealed two patterns of inter-species transmission of Sarcoptes mites: one is the transmission restricted by the "host taxon" [15][16][17] and the other is the transmission through "prey-predator" interaction [18][19][20]. For Host-taxon transmission restriction, Rasero et al. [15] found that Sarcoptes mite populations from ten European wild mammal species were genetically separated by three patterns of host taxon, named as carnivores, omnivores and herbivores (in the definition of Rasero et al. [15]). For predator-prey transmission, close genetic associations between mite populations derived from predators and their preferred prey were detected, e.g. between lion (Panthera leo) and wildebeest (Connochaetes taurinus), and between cheetah (Asinonyx jubatus) and Thomson's gazelle (Gazelle thompsonii) [18]. The difference between these two patterns of transmission of Sarcoptes mites suggests the existence of diversity in their transmission web, associated with the interaction within the local mammalian communities co-existing in each region [9]. In Japan, the realization of the epidemic of sarcoptic mange in wild mammal hosts was started in early 1980s with the report of an outbreak in raccoon dogs (Nyctereutes procyonoides) [21,22]. Since then, during the past 40 years, the disease has been drastically expanded in the omnivorous Carnivora species [e.g. red fox (Vulpes vulpes), masked palm civet (Paguma larvata)] and wild boars (Sus scrofa leucomystax) [21][22][23]. In addition, cases of infestation of Sarcoptes mites in Japanese serow (Capricornis crispus), an indigenous ruminant species designated as a Special National Monument in Japan, have been reported since the mid-1990s [21,24,25]. Because of the nearly simultaneous (within 10-20 years) and multiple beginnings of epidemics in wild mammals, it has been suspected that possible inter-species transmission occurred in these mammalian host communities, with the addition of transmission from domestic/companion animals (e.g. dog, swine and cattle) [21,22]. However, the genetic relationships of Sarcoptes mites with mammalian hosts have not been elucidated at the level of SSR resolution, apart from limited reports on the relationship between raccoon dog populations and domestic/companion dogs (Canis familiaris) [26].
In the present study, we aim to illuminate the neglected relationships of Sarcoptes mite populations derived from five wild mammalian species and domestic/companion dogs in Japan by SSR analysis. This enables us to explicate whether the transmission within these multi-host species follows the host taxon and predator-prey patterns. The target hosts were medium-sized omnivorous Carnivora, more specifically Caniformia [domestic dog, raccoon dog, raccoon (Procyon lotor) and Japanese marten (Martes melampus melampus)], an omnivorous Cetartiodactyla (wild boar) and a herbivorous Cetartiodactyla (Japanese serow).
To our knowledge, the predator-prey relationships between these omnivorous Caniformia and Cetartiodactyla have not been previously reported, except for a few limited cases of predation of Japanese serow by feral dogs [27] and possible occasional scavenging on carcasses of omnivorous species [28][29][30][31]. Against these inadequate predator-prey relationships and the expected Host-taxon barriers for the inter-species transmission, we report here the observation of "cryptic, " close genetic relationships among the Sarcoptes mites derived from Japanese serows and those from omnivorous species in Caniformia.

Collection of S. scabiei from skin crusts of hosts
From 1996 to 2016, we collected skin crusts from 93 mangy animals from 10 different prefectures in Japan: Tokyo, Kanagawa, Saitama, Gifu, Wakayama, Hyogo, Shimane and Yamaguchi on Honshu Island, and Oita and Saga on Kyushu Island (Fig. 1).
The animals belonged to 17 populations from 6 species that consisted of raccoon dogs (n = 60), domestic/companion dogs (n = 11), raccoons (n = 1), Japanese martens (n = 1), Japanese serows (n = 7) and wild boars (n = 13) ( Table 1). All skin crusts were stored at -30 °C or in 70% ethanol until the mite collection. We collected one Sarcoptes mite per individual animal by the postponed (postfrozen) isolation method for frozen skins proposed by Alasaad et al. [32]. All mites were confirmed to be S. scabiei from the morphological criteria defined by Fain [10].

Molecular analysis
Observed (H O ) and expected (H E ) heterozygosity, linkage disequilibrium (LD) and Hardy-Weinberg equilibrium (HWE) tests were calculated using GENEPOP v.4.0. [34]. In this software, the exact P-values of LD tests were given Fig. 1 The location of the 10 sampling areas and 93 sampled animals in the present study by Chi-square tests and those of HWE were calculated by a Markov chain algorithm [34].
While some of the analyzed mite populations significantly deviated from HWE, LD was also detected between some SSR loci. Therefore, we used two methods that do not require assumptions, like HWE and linkage equilibrium (LE) of loci, to analyze the underlying population genetic model. First, to describe differentiation between mite populations, principal components analysis (PCA) was carried out using the package adegenet [35] of R v.3.2.2 [36]. Multivariate ordination in the PCA does not require data regarding the HWE and LD [37]. In this analysis, genetic relationships among mites associated with 17 host populations (Table 1) were analyzed. Secondly, to analyze the genetic relationships between mites, we calculated genetic distances (proportion of shared alleles; D SA [38]) among mites using Populations v.1.2.32 [39]. These genetic distances were then visualized using the NeighborNet network provided by SplitsTree4 [40].
However, HWE in populations and LE between SSR loci are often violated in natural populations and are viewed as unrealistic [37]. In such circumstances, Bayesian clustering analysis by STRU CTU RE software v.2.3.4 [41] has often also been employed for populations deviated by HWE and/or for loci in LD (e.g. Rasero et al. [15]). Thus, we applied the STRU CTU RE analysis to elucidate the genetic structure of 93 Sarcoptes mites and to support the former 2 analytical methods. We performed 100,000 MCMC (Marcov Chain Monte Carlo) simulations followed by 50,000 steps of burn-in for each K value (K = 1-17), and these calculations were repeated independently 10 times. Subsequently, the most likely clustering number of K was determined by the peak of ΔK calculated by the method of Evanno et al. [42]. The website STRU CTU RE HARVESTER (taylor0.biology.ucla. edu/structureHarvester) [43] was used for this calculation. The results of 10 replicate runs for each value of K were averaged using the Greedy algorithm of CLUMPP v.1.1.2 [44], and clustering results for each value of K were displayed graphically using distruct v.1.1 [45].
The NeighborNet network composed by D SA showed the 5 genetic clusters (Clusters A-E) in 93 mites (Fig. 3). Basically, mites from the same taxonomic or close host species within the same area were included in each cluster; Cluster A comprised Caniformia-derived mites from Gifu (GCf and GN), Cluster B comprised Caniformiaderived mites from Tokyo (TCf and TN) and Kanagawa (KN), Cluster C comprised Caniformia-derived mites from Wakayama (WCf and WN), Cluster D comprised boar-derived mites (WS, HS, SmS and YS) and Cluster E comprised serow-derived mites from Oita (OCc). However, some host species of other areas and/or different host taxa were also classified in each cluster. Raccoon-and raccoon dog-derived mites (WP1, WN7 and WN8) from Wakayama, and a Japanese serow-derived mite from Saitama (StCc) were grouped in Cluster A. A domestic/companion dog-derived mite from Gifu (GCf9) was grouped in Cluster C. Raccoon dog-derived mites from Saga (SgN) and domestic/companion dog-derived mites from Gifu (GCf4 and GCf7) were grouped in Cluster E. Raccoon dog-derived mites from Gifu (GN5, GN7), a Japanese serow-derived mite and a Japanese martenderived mite from Wakayama (WCc, WM) showed just intermediate genotypes between Clusters B and C (we included them in Cluster C for the convenience in the discussion). All Japanese serow-derived mites were confirmed to be close to the Caniformia-derived mites.
The STRU CTU RE analysis revealed the highest value of ΔK at K = 9 (Additional file 2: Figure S1). At K = 9, a similar result to that obtained using the NeighborNet network was seen, wherein the wild boar-derived mite population and Caniformia-with Japanese serow-derived mite population were assigned to different clusters. Subsequently, a second run of STRU CTU RE analysis was performed for the 80 mites from the Caniformia-with Japanese serow-derived population, and the ΔK value showed the uppermost hierarchical cluster to be K = 4. The genetic structure in these 80 mites in K = 4 was also consistent with that observed using the NeighborNet network ( Fig. 3 and Additional file 2: Figure S1).

Discussion
As previously described, Rasero et al. [15] comprehensively presented the genetic differentiation of Sarcoptes mites in mammals in western Europe. Nonetheless, we have to re-attempt the interpretation of the results of their study here, before going into the profound discussion of our results. In our opinion, the results presented by Rasero et al. [15] lack accuracy because of the ambiguous definitions of terms related to "host taxon" (i.e. "carnivore", "omnivore" and "herbivore"). Although Rasero et al. [15] did not define "host taxon, " the words "carnivore, " "omnivore, " and "herbivore" are usually used in the context of the foraging ecology of animals. In their study,  Table 1 the red fox, stone marten (M. foina) and pine marten (M. martes) were defined as "carnivores. " These animals are widely recognized as omnivores based on their food habits [46,47]. Therefore, if these animals were classified by their food habits, their sorting would be incorrect. Thus, more accurately, the differentiation of Sarcoptes mite populations observed in their study did not conform to the foraging ecology of mammals, but to the taxonomy at the suborder level: (i) Caniformia in Carnivora as a substitute for "carnivore"; (ii) Suina in Cetartiodactyla as a substitute for "omnivore"; and (iii) Ruminantia in Cetartiodactyla as a substitute for "herbivore".
As far as the present study is concerned, the taxonomy of our target species was as follows: domestic dog, raccoon dog, raccoon and Japanese marten belong to the suborder Caniformia in Carnivora; wild boar belongs to the suborder Suina in Cetartiodactyla; and Japanese serow belongs to the suborder Ruminantia in Cetartiodactyla. Deficiency of observed heterozygosity and deviation from HWE in some of the SSR loci in a few mite populations were also confirmed in the previous reports of other Sarcoptes mite populations [15,16,18,20]. This phenomenon is due to a lack of random mating of mite populations between host animals, because Sarcoptes mites do not have free-living stages and individual hosts function as their ephemeral habitats in patchy environments [33]. Bearing in mind the presence of deviations in HWE, though a PCA method which will not be influenced by those genetic biases existed, we obtained the following consistent results in the multiple analyses: (i) the close genetic relationship between the mite populations derived from hosts in the same taxon, represented in omnivorous Caniformia and in wild boar populations; (ii) the separation of genetic clusters among mites derived from omnivorous Caniformia and from wild boar populations; and (iii) the close genetic relationship between mites derived from Caniformia species and from Japanese serow against the difference of host taxon. The results (i) and (ii) correspond to the previously reported genetic relationships of Sarcoptes mite populations in the European mammalian fauna and are consistent with the "Host-taxon phenomenon" formulated by Rasero et al. [15]. Meanwhile, (iii) shows the deviation from the Hosttaxon regulation in the transmission web. Those serowderived mites collected in the middle of 1990s (StCc1 in 1996 and OCc1-4 between 1996 and 1998) and in the 2010s (WCc1 in 2011 and Occ5 in 2016) from different locations highly suggested the long-time maintenance of  Table 1 gene flow between Caniformia-and serow-derived mite communities at several distant locations.
In the previous studies, the prey-predator pattern was described as "Curse of prey" [18] and tended to focus on the one-directional infections of Sarcoptes mites from the herbivorous preys to their carnivorous predators. However, this one-directional concept is not supported in case of "very weak" prey-predator relationships, the epidemiological history and the ecology of Japanese serow. Rather, infection from Caniformia mammals to serows is highly suspected. Upon reviewing the recent epidemiological history of sarcoptic mange in Japan, it was seen that its outbreaks in Japanese serows started roughly ten years later than the epidemic in raccoon dogs [21]. From a questionnaire survey for the epidemiology of sarcoptic mange in wild mammals, observations of mangy cases in hunting as well as in wildlife rescue practices began in the 1970s and 1980s in raccoon dogs and other Caniformia species [21,22]. In the mid-1990s, recordings in hunting and rescue cases in raccoon dogs had drastically increased [21]. During this time, cases of sarcoptic mange in the serows increased [21,25]. It might be possible that the epidemic in serows was associated with the increase of mangy cases in raccoon dogs (and possibly other Caniformia species).
From the view point of Japanese serow ecology, which is a territorial and monogamous species, usually living solitarily [48], the transmission of sarcoptic mange between serows has been considered to be less frequent than that in mammals which exhibit strong sociality (i.e. making herds or family groups) such as the raccoon dog, domestic dog and wild boar [2,49]. However, multi-locational outbreaks in serows were almost simultaneously recognized since the 1990s in several geographically distinct populations on the Honshu and Kyushu Islands of Japan [21]. This concurrent emergence of mange among distinct serow populations was possibly due to transmission from other Caniformia harboring S. scabiei. Based on the data obtained from a questionnaire survey conducted by the Wildlife Management Office [21], most of the epidemic areas of mange infection in Japanese serows overlapped with those of mange infection in raccoon dogs.
Given that the inter-species transmission of Sarcoptes mites arises through both direct or indirect contact [8], there are two likely explanations for the transmission from Caniformia populations to the serow populations: (i) the direct transmission upon touching the carcasses of mangy animals; and (ii) indirect transmission through the sharing of habitat intensively used by both mammal populations. For the direct contact, the probable reason is the exploring or foraging behavior for dead animals sometimes observed in ruminant species such as cattle [50] and deer [51]. There is no observation, but likelihood that serows behave similar to dead or debilitated mangy Caniformia animals. For the indirect contact, the common use of the same location (or microhabitats) by these species has often been recorded by sensor cameras [52]. Although the survival of Sarcoptes mites off the host is only a few days in the natural environment [8], infection through contaminated environments (e.g. animal trails, vegetation, soil, etc.) is not impossible [8,53]. Further study on the Sarcoptes mite's ecology and contact between host species, besides the transmission experiments, may furnish a clear answer for the observed gene flow between herbivore ruminants and Caniformia species.
Lastly, our data reveal that some Sarcoptes mites from Caniformia species and Japanese serows do not belong to dominant genotypes in the Caniformiaserow derived communities in their locations ( Fig. 3 and Additional file 2: Figure S1). If we assume that the results of the present study represent the gene pool of Sarcoptes mites in Caniformia mammals and Japanese serows in Japan, it may imply the migration of Sarcoptes mites, which is perhaps associated with (i) natural events of host migration and dispersion, or (ii) artificial translocation of Sarcoptes mites. Note that GCf4 and GCf7 from the domestic/companion dogs in Gifu belonged to the mite group from raccoon dogs and serows in Kyushu Island (OCc1-5 from Oita and SgN from Saga). Kyushu and Gifu are geographically distinct (> 550 km) and separated by sea. Thus, it is unlikely that this genetic relationship was formed by a natural event. Rather, this link could probably be owed to the artificial translocation of mangy animals, such as domestic/companion dogs. However, there are, of course, possible alternative hypotheses for the observed close genetic relationships between distant Sarcoptes mite populations, such as (i) an effect of a small sample size in mite populations between the two distant mite populations (i.e. limited observation and un-detection of similar genotypes), and (ii) a chance effect of mutation and convergence of genetic characteristics. Future research is needed for understanding of this type of genetic relationships and also the possible artificial impact for the expanding distribution of sarcoptic mange.

Conclusions
The results of our study show that in Japan, Sarcoptes mite populations follow a "Host-taxon" regulation between wild boars and Caniformia with Japanese serow communities but deviate between the Caniformia and Japanese serow. To our knowledge, this result provides the first genetic evidence for a possible hidden transmission of sarcoptic mange between host species belonging to different taxa without unambiguous predator-prey systems. In the diverted transmission patterns of S. scabiei among