Proof of the environmental circulation of veterinary drug albendazole in real farm conditions.

Anthelmintics, drugs against parasitic worms, are frequently used in livestock and might act as danger environmental microcontaminants. The present study was designed to monitor the possible circulation of common anthelmintic drug albendazole (ABZ) and its metabolites in the real agriculture conditions. The sheep were treated with the recommended dose of ABZ. Collected faeces were used for the fertilization of a field with fodder plants (alfalfa and clover) which served as feed for sheep from a different farm. The selective ultrasensitive mass spectrometry revealed surprisingly high concentrations of active ABZ metabolite (ABZ-sulphoxide) in all samples (dung, plants, ovine plasma, rumen content and faeces). Our results prove for the first time an undesirable permeation of ABZ metabolites from sheep excrement into plants (used as fodder) and subsequently to other sheep in real agricultural conditions. This circulation causes the permanent exposition of the ecosystems and food-chain to the drug and can promote the development of drug resistance in helminths.

Sub-lethal doses of albendazole induce drug metabolizing enzymes and increase albendazole deactivation in Haemonchus contortus adults.

The efficacy of anthelmintic therapy of farm animals rapidly decreases due to drug resistance development in helminths. In resistant isolates, the increased expression and activity of drug-metabolizing enzymes (DMEs), e.g. cytochromes P450 (CYPs), UDP-glycosyltransferases (UGTs) and P-glycoprotein transporters (P-gps), in comparison to sensitive isolates have been described. However, the mechanisms and circumstances of DMEs induction are not well known. Therefore, the present study was designed to find the changes in expression of CYPs, UGTs and P-gps in adult parasitic nematodes Haemonchus contortus exposed to sub-lethal doses of the benzimidazole anthelmintic drug albendazole (ABZ) and its active metabolite ABZ-sulfoxide (ABZSO). In addition, the effect of ABZ at sub-lethal doses on the ability to deactivate ABZ during consequent treatment was studied. The results showed that contact of H. contortus adults with sub-lethal doses of ABZ and ABZSO led to a significant induction of several DMEs, particularly cyp-2, cyp-3, cyp-6, cyp-7, cyp-8, UGT10B1, UGT24C1, UGT26A2, UGT365A1, UGT366C1, UGT368B2, UGT367A1, UGT371A1, UGT372A1 and pgp-3, pgp-9.1, pgp-9.2, pgp-10. This induction led to increased formation of ABZ metabolites (especially glycosides) and their increased export from the helminths' body into the medium. The present study demonstrates for the first time that contact of H. contortus with sub-lethal doses of ABZ (e.g. during underdose treatment) improves the ability of H. contortus adults to deactivate ABZ in consequent therapy.

Ivermectin-induced changes in the expression of cytochromes P450 and efflux transporters in Haemonchus contortus female and male adults.

Haemonchus contortus, one of the most pathogenic of all small ruminant parasites, have developed resistance to all used anthelmintics. Detoxification enzymes, e.g. cytochromes P450 (CYPs) and efflux transporters P-glycoproteins (P-gps), which represent the main defense system against harmful xenobiotics, have been suggested to contribute to drug resistance development. The present study was designed to compare the constitutive expression of individual CYPs and P-gps in females and males of H. contortus adults and to follow up on the changes in expression of these genes in nematodes exposed to sub-lethal concentrations of ivermectin (IVM), which might occur during inaccurate treatment. The adults of inbred susceptible-Edinburgh strain (ISE, MHco3) of H. contortus were used for this purpose. The nematodes were incubated ex vivo with or without IVM (1, 10 and 100 nM) in culture medium for 4, 12 and 24 h. After incubation, total RNA was isolated and expression levels of individual CYPs and P-gps were analyzed using qPCR. Our results showed a great variability in the constitutive expression of individual CYPs and P-gps in H. contortus adults. The constitutive expression as well as the inducibility of CYPs and P-gps significantly differed in males and females. Contact of adult nematodes with sub-lethal IVM concentrations led to only minor changes in expression of CYPs, while expression of several P-gps, particularly pgp-9.2 in males and pgp-10, pgp-11 in females was increased significantly in IVM-exposed nematodes. In conclusion, inaccurate treatment of sheep with IVM might contribute to drug resistance development via increased expression of efflux transporters in H. contortus adults.

UDP-glycosyltransferase family in Haemonchus contortus: Phylogenetic analysis, constitutive expression, sex-differences and resistance-related differences.

UDP-glycosyltransferases (UGT), catalysing conjugation of UDP-activated sugar donors to small lipophilic chemicals, are widespread in living organisms from bacteria to fungi, plant, or animals. The progress of genome sequencing has enabled an assessment of the UGT multigene family in Haemonchus contortus (family Trichostrongylidae, Nematoda), a hematophagous gastrointestinal parasite of small ruminants. Here we report 32 putative UGT genes divided into 15 UGT families. Phylogenetic analysis in comparison with UGTs from Caenorhabditis elegans, a free-living model nematode, revealed several single member homologues, a lack of the dramatic gene expansion seen in C. elegans, but also several families (UGT365, UGT366, UGT368) expanded in H. contortus only. The assessment of constitutive UGT mRNA expression in H. contortus adults identified significant differences between females and males. In addition, we compared the expression of selected UGTs in the drug-sensitive ISE strain to two benzimidazole-resistant strains, IRE and WR, with different genetic backgrounds. Constitutive expression of UGT368B2 was significantly higher in both resistant strains than in the sensitive strain. As resistant strains were able to deactivate benzimidazole anthelmintics via glycosylation more effectively then the sensitive strain, UGT368B2 enhanced constitutive expression might contribute to drug resistance in H. contortus.

Metabolism of albendazole, ricobendazole and flubendazole in Haemonchus contortus adults: Sex differences, resistance-related differences and the identification of new metabolites.

Haemonchus contortus (family Trichostrongylidae, Nematoda), a hematophagous gastrointestinal parasite found in small ruminants, has a great ability to develop resistance to anthelmintic drugs. We studied the biotransformation of the three benzimidazole anthelmintics: albendazole (ABZ), ricobendazole (albendazole S-oxide; RCB) and flubendazole (FLU) in females and males of H. contortus in both a susceptible ISE strain and resistant IRE strain. The ex vivo cultivation of living nematodes in culture medium with or without the anthelmintics was used. Ultrasensitive UHPLC/MS/MS analysis revealed 9, 7 and 12 metabolites of ABZ, RCB and FLU, respectively, with most of these metabolites now described in the present study for the first time in H. contortus. The structure of certain metabolites shows the presence of biotransformation reactions not previously reported in nematodes. There were significant qualitative and semi-quantitative differences in the metabolites formed by male and female worms. In most cases, females metabolized drugs more extensively than males. Adults of the IRE strain were able to form many more metabolites of all the drugs than adults of the ISE strain. Some metabolites were even found only in adults of the IRE strain. These findings suggest that increased drug metabolism may play a role in resistance to benzimidazole drugs in H. contortus.

Seroprevalence of Neospora caninum and Toxoplasma gondii in exotic ruminants and camelids in the Czech Republic.

Neospora caninum and Toxoplasma gondii are the protozoan parasites with definitive hosts from order Carnivora. Due to vertical transmission, both parasites can cause abortions and neonatal mortality that lead to significant productive and economic losses in the domestic ruminants. The aim of this study was to describe N. caninum and T. gondii seroprevalence in the group of frequently farmed captive exotic ruminants (n = 184) including Bovidae (barbary sheep, bezoar goat, common eland, American bison, water buffalo, and yak) and Camelidae (bactrian camel, guanaco, llama, and alpaca). Antibodies were tested by indirect fluorescent antibody test (IFAT) and enzyme-linked immunosorbent assay (ELISA). Higher prevalence of T. gondii antibodies (31% in IFAT and 52% in ELISA) was detected compared to N. caninum (24% in IFAT and 17% in cELISA). Mixed infection was found in 18 (10%) and 22 (12%) animals by IFAT and ELISA, respectively. Higher seroprevalence of both N. caninum and T. gondii was found in Camelidae compared to Bovidae. To author knowledge, this is the first detection of T. gondii and N. caninum antibodies in common elands and bezoar goats.

Serological Prevalence of Enteropathogenic Yersinia spp. in Pigs and Wild Boars from Different Production Systems in the Moravian Region, Czech Republic.

Human yersiniosis caused by pathogenic Yersinia spp. is one of the most common reported zoonoses in the European Union and pigs are considered as the major reservoir of these bacteria. Serological testing represents a suitable method to obtain information about the prevalence of enteropathogenic Yersinia spp. in food animals. The prevalence of antibodies against enteropathogenic Yersinia spp. was studied in 319 slaughtered pigs and 135 wild boars from different production systems in the Moravian region (Czech Republic) using a commercially available ELISA test (an apparent prevalence). The seroprevalence was significantly associated with the type of breeding system, with the lowest seroprevalence being observed in household-raised pigs (13/29, 44.8%). No significant difference between the prevalence of anti-Yersinia antibodies in conventional (146/180, 81.1%) and organic pigs (92/110, 83.6%) was found. Antibodies were found in 65.9% (89/135) of wild boars without a significant difference between adult (23/41, 56.1%) and young (66/94, 70.2%) animals. Seropositivity was significantly higher in domestic (251/319, 78.7% in total) compared to feral pigs. A Bayesian approach taking into account the sensitivity and specificity of the ELISA test was used to estimate the true prevalence of anti-Yersinia antibodies in pigs and wild boars. According to our results, domestic pigs and wild boars proved to be an important reservoir of enteropathogenic Yersinia in the Czech Republic. Attention should be paid to good hygienic practice during slaughtering and handling of meat to prevent meat contamination and subsequently human infection.

Albendazole in environment: faecal concentrations in lambs and impact on lower development stages of helminths and seed germination.

Albendazole (ABZ), widely used benzimidazole anthelmintic, administered to animals enters via excrements into environment and may impact non-target organisms. Moreover, exposure of lower development stages of helminths to anthelmintics may also encourage the development of drug-resistant strains of helminths. In present project, the kinetics of ABZ (10 mg kg(-1) p.o.) and its metabolite (ABZ.SO, ABZSO2) elimination in faeces from treated Texel lambs were studied using UHPLC/MS/MS with the aim to find out their concentrations achievable in the environment. Consequently, the effect of these compounds on lower development stages of Barber's pole worm (Haemonchus contortus) and on germination of white mustard (Sinapis alba) seeds was evaluated. The results showed that ABZ concentrations in faeces excreted in 4-60 h after treatment were above the concentrations lethal for H. contortus eggs. Moreover, pre-incubation with sub-lethal doses of ABZ and ABZ.SO did not increase the resistance of H. contortus eggs and larvae to anthelmintics. On the other hand, concentrations of ABZ and ABZ.SO in faeces are so high that might have negative influence on non-target soil invertebrates. As neither ABZ nor its metabolites affect the germination of mustard seeds, phytoremediation could be considered as potential tool for detoxification of ABZ in the environment.

The Role of Xenobiotic-Metabolizing Enzymes in Anthelmintic Deactivation and Resistance in Helminths.

Xenobiotic-metabolizing enzymes (XMEs) modulate the biological activity and behavior of many drugs, including anthelmintics. The effects of anthelmintics can often be abolished by XMEs when the drugs are metabolized to an inefficient compound. XMEs therefore play a significant role in anthelmintic efficacy. Moreover, differences in XMEs between helminths are reflected by differences in anthelmintic metabolism between target species. Taking advantage of the newly sequenced genomes of many helminth species, progress in this field has been remarkable. The present review collects up to date information regarding the most important XMEs (phase I and phase II biotransformation enzymes; efflux transporters) in helminths. The participation of these XMEs in anthelmintic metabolism and their possible roles in drug resistance are evaluated.

Comparison of biotransformation and efficacy of aminoacetonitrile anthelmintics in vitro.

The present in vitro study was designed to test and compare anthelmintic activity, hepatotoxicity, and biotransformation of four selected aminoacetonitrile derivatives (AADs): monepantel (MOP, anthelmintic approved for the treatment), AAD-970, AAD-1154, and AAD-1336. Micro-agar larval development test, MTT test of cytotoxicity, and biotransformation study coupled with Ultra high performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS) technique were used for this purpose. Larvae of two Haemonchus contortus strains (drug susceptible and multi-drug resistant) and primary cultures of rat and ovine hepatocytes served as model systems. All AADs (including MOP) exhibited significant larvicidal effect in H. contortus susceptible as well as multi-resistant strains, much higher than those of reference anthelmintics thiabendazole and flubendazole. AAD-1154 provides the best results for most tested parameters among all AADs in this study. The cytotoxicity test showed that all AADs can be considered as nontoxic for hepatocytes. In the biotransformation study, Phase I and Phase II metabolites of AADs were identified and schemes of possible metabolic pathways in ovine hepatocytes were proposed. Biotransformation of MOP was much more extensive than biotransformation of other AADs. Based on obtained results, AAD-1154 and AAD-1336 can be considered as promising candidates for further in vivo testing.

Metabolism of drugs and other xenobiotics in giant liver fluke (Fascioloides magna).

1. Giant liver fluke Fascioloides magna is a dangerous parasite, which infects herbivores. It was imported to Europe from North America and started to spread. Benzimidazoles like albendazole, mebendazole, triclabendazole and salicylanilides closantel and rafoxanide are the most used anthelmintics to control fascioloidosis. However their effect might be altered via drug-metabolizing enzymes of this parasite. 2. The aim of our study was to determine the activities of drug-metabolizing enzymes in F. magna and the metabolism of above mentioned anthelmintics. 3. Activities of several oxidative, reductive and conjugative enzymes towards various model xenobiotic substrates were found in F. magna subcellular fractions. 4. Subcellular fractions from F. magna oxidized albendazole to its sulphoxide metabolite and reduced mebendazole to hydroxyl-mebendazole. Under ex vivo conditions, only very-low concentrations of these compounds were detected using high-performance liquid chromatography/mass spectrometry. 5. The results indicate that the giant liver fluke possesses the active xenobiotic-metabolizing system. The overexpression of this system may play an important role in parasite resistance against these anthelmintics.

Seasonal Dynamics, Parity Rate, and Composition of Culicoides (Diptera: Ceratopogonidae) Occurring in the Vicinity of Wild and Domestic Ruminants in the Czech Republic.

In the light of the emergence of bluetongue and Schmallenberg viruses in northern and central Europe, an extensive entomological survey within the framework of a bluetongue control program was undertaken from 2008 to 2013 in the Czech Republic to investigate Culicoides biting midges (Diptera: Ceratopogonidae) collected in close proximity of domestic livestock and semiwild ruminants. Insects were sampled using CDC black-light suction traps placed overnight near ruminants in farms or in forest game preserves to provide data on Culicoides fauna collected near these two groups of hosts inhabiting different environments. From almost a half million biting midge specimens collected at 41 sampling sites, 34 species were identified including three species newly recorded for the Czech Republic: Culicoides (Oecacta) clastrieri Callot, Kremer & Deduit, Culicoides (Oecacta) odiatus Austen, and Culicoides (Pontoculicoides) saevus Kieffer. The Culicoides obsoletus species group, incriminated as a bluetongue virus vector, was predominant in both domestic livestock (91%) and semiwild game (52%). A relatively high proportion (around 30%) of C. obsoletus Meigen females with pigmented abdomen (= more likely parous) was observed from spring till autumn. In contrast, adult biting midges were found to be largely absent during at least three winter months, approximately December till March, which could be considered as the biting midge vector-free period.

Reliable reference gene selection for quantitative real time PCR in Haemonchus contortus.

The aim of this work was to identify reliable reference genes for expression studies in adult Haemonchus contortus. Eleven candidate genes were identified and the stability of their expression was assessed in adult males and females of two genetically divergent H. contortus isolates: drug-susceptible (ISE) and multi-drug-resistant (WR). Five genes with the most stable expression patterns were further assessed for suitability as reference genes in anthelmintic-treated H. contortus adults versus non-treated controls. We identified important differences in the expression of a number of candidate genes in anthelmintic-treated samples, confirming the need for careful validation of control genes for such experiments. We propose the use of multiple reference genes for expression studies in this species and found gpd, ama and far most suitable for adult H. contortus.

Prevalence of Hepatitis E Virus in Populations of Wild Animals in Comparison with Animals Bred in Game Enclosures.

Hepatitis E virus (HEV) is now accepted as a zoonotic virus, and domestic pigs, wild boars and deer are recognised as natural reservoirs of the pathogen. In this study, 762 animals (wild boars, fallow deer, red deer, sika deer, roe deer and mouflons) originating from the wild and from game enclosures were tested for the presence of HEV RNA by qRT-PCR. HEV RNA was detected in wild boars (96/450), red deer (2/169), roe deer (1/30) and mouflons (5/39). The sequence relationship between HEV isolates from wild boars and domestic pigs or humans indicate a circulation of HEV in the Czech Republic.

Monepantel induces hepatic cytochromes p450 in sheep in vitro and in vivo.

Monepantel (MOP), a new amino-acetonitrile anthelmintic for the treatment and control of gastrointestinal nematode infections and associated diseases in sheep, is approved and marketed as oral solution under the trade name Zolvix® (Novartis Animal Health Inc., Switzerland). The effect of MOP on hepatic cytochromes P450 (CYP) has been investigated in sheep. In an in vivo experiment, castrated rams (9-months old) were treated with the recommended therapeutic dose of MOP. Non-treated animals represented the controls. After 24 h, the animals were stunned and exsanguinated. Microsomal fractions and total RNA were prepared from liver homogenates. The activities towards alkyloxyresorufins, 7-methoxy-4-trifluoromethylcoumarin and midazolam were assayed and mRNAs of individual CYP isoforms were quantified. In an in vitro procedure, primary cultures of ovine hepatocytes were incubated with or without MOP (10 µM) for 24 h and then expression levels of individual CYP isoforms were analyzed. Results showed that MOP significantly increased all CYP-related activities and CYP3A24 mRNA in sheep. The induction effect of MOP on CYP3A was similar or even higher than those of dexamethasone and rifampicin, well-known CYP3A inducers. As CYP3A enzymes belongs to the most important biotransformation enzymes, their induction may have serious pharmacological and/or toxicological consequences. These facts should be taken into account when other drugs together with or after MOP (Zolvix®) are administered to sheep.

Biotransformation of anthelmintics and the activity of drug-metabolizing enzymes in the tapeworm Moniezia expansa.

The sheep tapeworm Moniezia expansa is very common parasite, which affects ruminants such as sheep, goats as well as other species. The benzimidazole anthelmintics albendazole (ABZ), flubendazole (FLU) and mebendazole (MBZ) are often used to treat the infection. The drug-metabolizing enzymes of helminths may alter the potency of anthelmintic treatment. The aim of our study was to assess the activity of the main drug-metabolizing enzymes and evaluate the metabolism of selected anthelmintics (ABZ, MBZ and FLU) in M. expansa. Activities of biotransformation enzymes were determined in subcellular fractions. Metabolites of the anthelmintics were detected and identified using high performance liquid chromatography/ultra-violet/VIS/fluorescence or ultra-high performance liquid chromatography/mass spectrometry. Reduction of MBZ, FLU and oxidation of ABZ were proved as well as activities of various metabolizing enzymes. Despite the fact that the conjugation enzymes glutathione S-transferase, UDP-glucuronosyl transferase and UDP-glucosyl transferase were active in vitro, no conjugated metabolites of anthelmintics were identified either ex vivo or in vitro. The obtained results indicate that sheep tapeworm is able to deactivate the administered anthelmintics, and thus protects itself against their action.

Influence of Stress Connected with Moving to a New Farm on Potentially MAP-Infected Mouflons.

There is no European legislation concerning paratuberculosis that requires that imported animals be kept in quarantine and commonly they are directly released into areas with other animals. In this study, detection of latent infection of paratuberculosis in healthy mouflons previously diagnosed as paratuberculosis-free, but originating from a real time quantitative PCR- (qPCR-) positive herd, occurred after their transport to a new farm. During a twelve-day quarantine period, all mouflons irregularly shed Mycobacterium avium subsp. paratuberculosis (MAP) in faeces, and in a small number of cases also in milk. After the animals were released from quarantine, MAP was detected for a further two days, after which, testing was negative, except in one case. Therefore, the stress connected with transport, novel environment, dietary change, or limited area with high density of animals might have contributed to the induction of paratuberculosis and the shedding of MAP from the animals, previously diagnosed as MAP-negative. According to these results, the keeping of imported animals in quarantine and their examination for MAP presence not only before the transport but also afterwards should be recommended. The designation of a particular area of a farm as a quarantine enclosure could help to mitigate the impact of stress caused by a confined space with a high density of animals.

Metabolic pathways of anthelmintic drug monepantel in sheep and in its parasite (Haemonchus contortus).

Monepantel (MOP) is a new anthelmintic drug intended for the treatment and control of gastrointestinal roundworms (nematodes) infection and associated disease in sheep. The aim of our study was to find out metabolic pathways of MOP in sheep in vivo and in its parasite Haemonchus contortus ex vivo. MOP biotransformation in two H. contortus strains with different sensitivity to anthelmintics was also compared. Ultra high-performance liquid chromatography with tandem mass spectrometry (UHPLC-MS/MS) technique is used for the identification of MOP metabolites in ovine urine, faeces, and nematodes. MOP biotransformation study in sheep in vivo led to the identification of 13 MOP metabolites; 7 of them have not been described previously in in vitro study. The study of MOP biotransformation in H. contorus ex vivo reveals four MOP metabolites. The nitrile hydrolysis as a new biotransformation pathway in helminths ex vivo is reported here for the first time. Unlike sheep, H. contorus nematodes are not able to metabolize MOP via phase II biotransformation. Nematodes of resistant White river (WR) strain form more types of MOP metabolites than nematodes of sensitive inbred susceptible Edinburgh (ISE) strain. Based on obtained results, schemes of metabolic pathways of MOP in sheep and nematodes are proposed.

Biotransformation of albendazole and activities of selected detoxification enzymes in Haemonchus contortus strains susceptible and resistant to anthelmintics.

The increased activity of drug-metabolizing enzymes can protect helminths against the toxic effect of anthelmintics. The aim of this study was to compare the metabolism of the anthelmintic drug albendazole (ABZ) and the activities of selected biotransformation and antioxidant enzymes in three different strains of Haemonchus contortus: the ISE strain (susceptible to common anthelmintics), the BR strain (resistant to benzimidazole anthelmintics) and the WR strain (multi-resistant). H. contortus adults were collected from the abomasum of experimentally infected lambs. In vitro (subcellular fractions of H. contortus homogenate) as well as ex vivo (living nematodes cultivated in flasks with medium) experiments were performed. HPLC with spectrofluorimetric and mass-spectrometric detection was used in the analysis of ABZ metabolites. The in vitro activities of oxidation/antioxidation and conjugation enzymes toward model substrates were also assayed. The in vitro data showed significant differences between the susceptible (ISE) and resistant (BR, WR) strains regarding the activities of peroxidases, catalase and UDP-glucosyltransferases. S-oxidation of ABZ was significantly lower in BR than in the ISE strain. Ex vivo, four ABZ metabolites were identified: ABZ sulphoxide and three ABZ glucosides. In the resistant strains BR and WR, the ex vivo formation of all ABZ glucosides was significantly higher than in the susceptible ISE strain. The altered activities of certain detoxifying enzymes might partly protect the parasites against the toxic effect of the drugs as well as contribute to drug-resistance in these parasites.

Dybowski's sika deer (Cervus nippon hortulorum): genetic divergence between natural primorian and introduced Czech populations.

Dybowski's sika deer (Cervus nippon hortulorum) originally inhabited the majority of the Primorsky Krai in Far Eastern Russia, north-eastern China, and Korean Peninsula. At present, only the Russian population seems to be stable, even though this taxon is still classified as endangered by the Russian Federation. Almost 100 years ago, this subspecies, among others, was imported to several European countries including the Czech Republic. We used both mitochondrial (mtDNA; the cytochrome b gene and the control region) and nuclear DNA markers to examine the actual taxonomic status of modern Czech Dybowski's sika population and to compare the genetic diversity between the introduced and the native populations. Altogether, 124 Czech samples and 109 Primorian samples were used in the analyses. Within the samples obtained from individuals that were all morphologically classified as Dybowski's sika, we detected mtDNA haplotypes of Dybowski's sika (84 samples), as well as those belonging to other sika subspecies: northern Japanese sika (25 samples), southern Japanese sika (6 samples), and south-eastern Chinese sika (8 samples). Microsatellite analysis revealed a certain level of heterozygote deficiency and a high level of inbreeding in both populations. The high number of private alleles, factorial correspondence analysis, and Bayesian clustering analysis indicate a high level of divergence between both populations. The large degree of differentiation and the high number of population-specific alleles could be a result of a founder effect, could be a result of a previously suggested bottleneck within the Primorian population, and could also be affected by the crossbreeding of captive individuals with other sika subspecies.