Lessons in Détente or know thy host: The immunomodulatory gene products of myxoma virus

The poxvirus, myxoma virus, encodes within its genome at least eleven different proteins that compromise, skew, or disable the innate and adaptive responses of its hosts. In the laboratory rabbit, Oryctolagus cuniculus, these effects result in myxomatosis, a fatal condition characterized by skin lesions and systemic immunosuppression. Interestingly, while myxoma infection also causes skin lesions in its natural host and in natural populations of O. cuniculus in Australia where this novel host and the virus have co-evolved, the condition of myxomatosis does not ensue and infection is not fatal. In this review I discuss the biochemical properties of the characterized immunomodulatory proteins of myxoma virus, and their pathogenic effects in laboratory rabbits. Disruption of any one myxoma immunomodulatory gene diminishes the severity of the infection without compromising infectivity. Thus, the characterized immunomodulatory genes appear not to be required for a productive infection in vivo. The differences in the severity of their effects in laboratory-bred versus wild O. cuniculus suggest that the outcome of myxoma infection is a consequence of the interplay between the viral immunomodulatory gene products and the cells and molecules of the host immune system.     

Evolution of viral sensing RIG-I-like receptor genes in Leporidae genera Oryctolagus, Sylvilagus, and Lepus

One of the most severe European rabbit (Oryctolagus cuniculus) pathogens is myxoma virus (MYXV), a rabbit-specific leporipoxvirus that causes the highly lethal disease myxomatosis. Other leporid genera, Sylvilagus and Lepus, encompass species with variable susceptibilities to MYXV, but these do not develop the lethal form of the disease. The protective role of the retinoic acid-inducible gene-I (RIG-I/DDX58) in sensing MYXV in nonpermissive human myeloid cells prompted the study of the RIG-I-like receptor (RLR) family evolution in the three leporid genera. This viral-sensor family also includes the melanoma differentiation-associated factor 5 (MDA5/IFIH1), and the laboratory of genetics and physiology 2 (LGP2/DHX58). Considering specifically the MYXV susceptible host (European rabbit) and one of the virus natural long-term hosts (Sylvilagus bachmani, brush rabbit), the amino acid differences of positively selected sites in RIG-I between the two species were located in the protein region responsible for viral RNA recognition and binding, the repressor domain. Such differences might play a determinant role in how MYXV is sensed. When looking for episodic selection on MDA5 and LGP2 of the eastern cottontail (Sylvilagus floridanus), we also uncovered evidence of selective pressures that might be exerted by a species-specific leporipoxvirus, the Shope fibroma virus. Finally, a putative alternative splicing case was identified in Oryctolagus and Lepus MDA5 isoforms, corresponding to the deletion of one specific exon. This study provided the first insights into the evolution of the leporid RLR gene family that helps illuminate the origins of the species-specific innate responses to pathogens and more specifically to MYXV.     

European rabbit research in the Iberian Peninsula: state of the art and future perspectives

The European rabbit (Oryctolagus cuniculus) is a high-profile prey, native from the Iberian Peninsula, the only region in the world where the two rabbit subspecies (O. cuniculus algirus and O. cuniculus cuniculus) currently co-exist in natural conditions. In this area, this important prey represents a keystone species and ecosystem engineer of Mediterranean landscapes, being also the most harvested and one of the most managed small-game species. Additionally, the species can create damage to crops in some parts of the Iberian Peninsula where it is regarded as an agricultural pest. The scientific interest towards the species is becoming increasingly apparent most likely as a repercussion of declining population trends over the last decades. The latter has been the result of the impact of habitat deterioration, viral diseases, unsustainable hunting, and predation. In this paper, I present a review of the scientific literature currently available on the European rabbit in the Iberian Peninsula. I discuss knowledge gaps and highlight priority research guidelines to suppress them, in an attempt to provide a general perspective to target research efforts more effectively. This analysis is particularly relevant due to the current vulnerability of rabbit populations in Iberia and to the recent news of cuts in scientific funding in most Mediterranean countries.     

Convergent evolution of IL-6 in two leporids (Oryctolagus and Pentalagus) originated an extended protein

Interleukin 6 (IL-6) is a class-I helical cytokine with a broad spectrum of biological activities and a gene structure conserved throughout vertebrates, with five coding exons. IL-6 from European rabbits belonging to the subspecies Oryctolagus cuniculus cuniculus was previously shown to differ from other mammals by extending an additional 27 amino acids. However, in other leporids (Sylvilagus spp and Lepus spp) that diverged from the European rabbit ~12 million years ago this mutation was not present. In this study, we extended the study of IL-6 for the Oryctolagus cuniculus algirus subspecies and five additional lagomorphs’ genera (Brachylagus, Bunolagus, Pentalagus, Romerolagus, and Ochotona). We confirmed the presence of the mutated stop codon in both O. c. cuniculus and O. c. algirus. We found that the typical stop codon is present in Sylvilagus bachmani and Lepus europaeus, in agreement with previous reports, but also in Bunolagus, Brachylagus, and Ochotona. Remarkably, in Pentalagus we detected a deletion of the stop codon causing an extension of IL-6 for 17 extra residues. Our results indicate that the IL-6 extension in those species occurred by two independent events: one occurred between 2 and 8 million years ago in the ancestral of the Oryctolagus subspecies, and the other occurred in a Pentalagus ancestral at a maximum of 9 million years ago. The absence of this IL-6 extension in Bunolagus, sister genus of Oryctolagus, shows that this evolutionary event happened by convergence suggesting some functional relevance.     

Biometrical analysis reveals major differences between the two subspecies of the European rabbit

The climatic oscillations that have occurred in the last few million years have strongly affected species distribution ranges. Highly divergent genetic lineages arose, some of which correspond to recognized subspecies that currently occupy small geographical areas. Understanding the implications of the genetic differences between these subspecies is crucial for proper conservation of Evolutionarily Significant Units. We use the two European rabbit subspecies, Oryctolagus cuniculus cuniculus and O. c. algirus, in the Iberian Peninsula as a model to investigate the repercussions at the biometric level of their largely recognized genetic differentiation. To accomplish this we analysed the ear and hind foot length, and the body mass of 999 adult rabbits from 27 locations across the distribution range of both subspecies in their native range, the Iberian Peninsula. Our results show biometric differences between the two subspecies, also explained by geographical location and sex, O. c. algirus being lighter and having shorter ear and hind foot lengths. We examine these findings under an evolutionary framework, and discuss their implications for current conservation efforts. Future research should focus on the ecological implications of these biometric differences, namely potential different habitat use and anti‐predatory strategies in the species' native range.     

Serp2, an Inhibitor of the Interleukin-1β-Converting Enzyme, Is Critical in the Pathobiology of Myxoma Virus

Recently, myxoma virus was shown to encode an additional member of the serpin superfamily. The viral gene, called serp2, was cloned, and the Serp2 protein was shown to specifically bind to interleukin-1β (IL-1β)-converting enzyme (ICE), thus inhibiting the cleavage of pro-IL-1β by the protease (F. Petit, S. Bertagnoli, J. Gelfi, F. Fassy, C. Boucraut-Baralon, and A. Milon, J. Virol. 70:5860–5866, 1996). Here, we address the role of Serp2 in the development of myxomatosis, a lethal infectious disease of the European rabbit. A Serp2 mutant myxoma virus was constructed by disruption of the single-copy serp2 gene and insertion of the Escherichia coli gpt gene serving as the selectable marker. A revertant virus was obtained by replacing the E. coli gpt gene by the intact serp2 open reading frame. The Serp2⁻ mutant virus replicated with wild-type kinetics both in rabbit fibroblasts and a rabbit CD4⁺ T-cell line (RL5). Moderate reduction of cell surface levels of major histocompatibility complex I was observed after infection with wild-type or Serp2⁻ mutant myxoma virus, and both produced white pocks on the chorioallantoic membrane of the chick embryo. After the infection of European rabbits, the Serp2⁻ mutant virus proved to be highly attenuated compared to wild-type myxoma virus, as demonstrated by the clinical course of myxomatosis and the survival rates of infected animals. Pathohistological examinations revealed that infection with wild-type myxoma virus resulted in a blockade of the inflammatory response at the vascular level. In contrast, rapid inflammatory reactions occurred upon infection with the Serp2⁻ mutant virus. Furthermore, lymphocytes in lymph nodes derived from animals inoculated with Serp2 mutant virus were shown to rapidly undergo apoptosis. We postulate that the virulence of myxoma virus in the European rabbit can be partially attributed to an impairment of host inflammatory processes and to the prevention of apoptosis in lymphocytes. The weakening of host defense is directly linked to serp2 gene function and is likely to involve the inhibition of IL-1β-converting-enzyme-dependent pathways.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Local feeding specialization by badgers (Meles meles) in a mediterranean environment

A case of local feeding specialization in the European badger (Meles meles), a carnivore species with morphological, physiological and behavioural traits proper to a trophic generalist, is described. For the first time, we report a mammalian species, the European rabbit (Oryctolagus cuniculus), as the preferred prey of badgers. Secondary prey are consumed according to their availability, compensating for temporal fluctuations in the abundance of rabbit kittens. We discuss how both predator (little ability to hunt) and prey (profitability and predictability) features, may favour the observed specialization, as predicted by foraging theory. Badgers show a trend to specialize on different prey in different areas throughout the species range. It is suggested that changes in prey features can reverse the badger feeding strategy at the population level. Such dynamic behavioural responses make difficult to label badgers as generalists or specialists at the species level.     

Recombination and Speciation: Loci Near Centromeres Are More Differentiated Than Loci Near Telomeres Between Subspecies of the European Rabbit (Oryctolagus cuniculus)

Recent empirical and theoretical studies suggest that regions of restricted recombination play an important role in the formation of new species. To test this idea, we studied nucleotide variation in two parapatric subspecies of the European rabbit (Oryctolagus cuniculus). We surveyed five loci near centromeres, where recombination is expected to be suppressed, and five loci near telomeres, where recombination is expected to be higher. We analyzed this multilocus data set using a divergence-with-gene flow framework and we report three main findings. First, we estimated that these subspecies diverged ~1.8 MYA and maintained large effective population sizes (O. c. algirus N_e ≈ 1,600,000 and O. c. cuniculus N_e ≈ 780,000). Second, we rejected a strict allopatric model of divergence without gene flow; instead, high rates of gene flow were inferred in both directions. Third, we found different patterns between loci near centromeres and loci near telomeres. Loci near centromeres exhibited higher levels of linkage disequilibrium than loci near telomeres. In addition, while all loci near telomeres showed little differentiation between subspecies, three of five loci near centromeres showed strong differentiation. These results support a view of speciation in which regions of low recombination can facilitate species divergence in the presence of gene flow.     

Influence of climate and community composition on the population demography of pasture species in semi-arid Australia

Substantial recruitment of Callitris glaucophylla in woodland, Sclerolaena birchii in cleared woodland, and Astrebla lappacea in grassland is related to catastrophic events of the past century in the form of interactions between climate, the impact of European land use (sheep, cattle, rabbits) and the rabbit myxoma epizootic. The direct effect of rainfall on the demography of these species and its indirect effect through competition via suites of accompanying plant species are examined. Major long-term changes in plant populations are generated by extreme sequential events rather than by random isolated events. One of the most potent climatic agents for change in eastern Australia is the El Niño/Southern Oscillation phenomenon.     

Involvement of the second extracellular loop and transmembrane residues of CCR5 in inhibitor binding and HIV-1 fusion: insights into the mechanism of allosteric inhibition

C-C chemokine receptor 5 (CCR5), a member of G-protein-coupled receptors, serves as a coreceptor for human immunodeficiency virus type 1 (HIV-1). In the present study, we examined the interactions between CCR5 and novel CCR5 inhibitors containing the spirodiketopiperazine scaffolds AK530 and AK317, both of which were lodged in the hydrophobic cavity located between the upper transmembrane domain and the second extracellular loop (ECL2) of CCR5. Although substantial differences existed between the two inhibitors—AK530 had 10-fold-greater CCR5-binding affinity (K_d = 1.4 nM) than AK317 (16.7 nM)—their antiviral potencies were virtually identical (IC₅₀ = 2.1 nM and 1.5 nM, respectively). Molecular dynamics simulations for unbound CCR5 showed hydrogen bond interactions among transmembrane residues Y108, E283, and Y251, which were crucial for HIV-1-gp120/sCD4 complex binding and HIV-1 fusion. Indeed, AK530 and AK317, when bound to CCR5, disrupted these interhelix hydrogen bond interactions, a salient molecular mechanism enabling allosteric inhibition. Mutagenesis and structural analysis showed that ECL2 consists of a part of the hydrophobic cavity for both inhibitors, although AK317 is more tightly engaged with ECL2 than AK530, explaining their similar anti-HIV-1 potencies despite the difference in K_d values. We also found that amino acid residues in the β-hairpin structural motif of ECL2 are critical for HIV-1-elicited fusion and binding of the spirodiketopiperazine-based inhibitors to CCR5. The direct ECL2-engaging property of the inhibitors likely produces an ECL2 conformation, which HIV-1 gp120 cannot bind to, but also prohibits HIV-1 from utilizing the “inhibitor-bound” CCR5 for cellular entry—a mechanism of HIV-1's resistance to CCR5 inhibitors. The data should not only help delineate the dynamics of CCR5 following inhibitor binding but also aid in designing CCR5 inhibitors that are more potent against HIV-1 and prevent or delay the emergence of resistant HIV-1 variants.     

European rabbit survival and recruitment are linked to epidemiological and environmental conditions in their exotic range

The European rabbit, Oryctolagus cuniculus, is threatened within its native range, yet it is a highly successful colonizing pest species across its worldwide introduced range, causing large economic losses and widespread environmental degradation. To date, there has been no long‐term empirical evidence documenting the relative roles of climatic, epidemiological and biological factors in limiting life‐history determinants of rabbit range and abundance. Using 12 years of capture–mark–recapture data from their exotic range in Australia, we constructed candidate Cormack–Jolly–Seber models to test the influence of environmental, competition and disease conditions on rabbit survival and recruitment. Our results show that: (i) population‐level disease infection rate has the largest overall impact on rabbit survival, explaining 80% of variance in survival rates; (ii) environmental as well as epidemiological conditions constrain rabbit survival, especially for younger animals; (iii) temporal variation in rabbit kitten recruitment patterns are best described by a combination of climate, competition and disease settings (accounting for 68% of variance), while temperature alone has a strong negative influence on kitten recruitment; and (iv) recruitment responds positively to rabbit haemorrhagic disease, but negatively to myxomatosis – the former, probably being mediated through a disease driven effect on intraspecific competition for food. A strengthened understanding of climate change impacts on rabbit range and abundance can be achieved by accounting explicitly for potential synergisms between disease dynamics and climate. In this analysis, we provide the first step towards such an attempt for this important mammal species. Integrated approaches of this kind are essential for future forecasts of rabbit range and abundance, offsetting the conservation threat faced by O. cuniculus in its native range, and achieving effective management in exotic habitats.     

Oxidative stress in wild European rabbits naturally infected with myxoma virus and rabbit haemorrhagic disease virus

The European rabbit (Oryctolagus cuniculus) is one of the most important vertebrate species in the Mediterranean Basin ecosystem. Over the last 60 years, the arrival of two viral diseases, myxomatosis and rabbit haemorrhagic disease, have led to dramatic declines in wild rabbit populations across the Iberian Peninsula. These diseases are currently endemic. Periodic outbreaks occur and have significant impacts on wild populations. Both infection types have diverse physiological effects on their hosts that are rooted in aerobic metabolic processes. To fight off these viruses, rabbits activate their immune systems. However, the production of immune defences generates reactive oxygen species that may consequently damage host tissues. Hypothesising that immune responses increase oxidative stress, we examined whether wild rabbits naturally infected with myxoma virus (MV) and rabbit haemorrhagic disease virus (RHDV) had high oxidative stress. Using blood samples, we measured anti-MV and anti-RHDV antibody concentrations and different oxidative stress markers (i.e., glutathione peroxidase, glutathione reductase, superoxide dismutase, catalase, and malondialdehyde). Our results show that rabbits that were seropositive for both MV and RHDV had high concentrations of malondialdehyde. Age and body condition were also positively related to dual seropositivity. No significant relationships were observed between serostatus and the concentrations of the other oxidative stress markers. Although we expected infection with MV and RHDV to be correlated with oxidative stress, the influence of external sources of oxidative stress (e.g., climatic conditions) likely made it more difficult to detect such relationships in wild rabbits.     

Rabbit (Oryctolagus cuniculus) abundance and protected areas in central-southern Spain: why they do not match?

European rabbits (Oryctolagus cuniculus), a keystone species in the Iberian Mediterranean ecosystem, are the staple prey of the Iberian lynx (Lynx pardinus) and the Spanish imperial eagle (Aquila adalberti). These predators require medium to high rabbit densities and a low degree of human disturbance. We compared rabbit abundances in areas of central-southern Spain under three levels of protection and management: protected areas, intensively managed (nonprotected) hunting estates, and other nonprotected areas. We used pellet abundance indices to estimate rabbit density in 118 surveys conducted during the summers of 2002 and 2003. We observed greater rabbit abundance in intensively managed hunting estates compared to protected areas and other nonprotected areas, perhaps because policy makers did not consider rabbit numbers when selecting priority areas. Alternatively, differences in game management practices (e.g., predator control or habitat management) may explain the higher rabbit densities observed in managed hunting estates. Our results suggest that the best feeding conditions for the Iberian lynx and the Spanish imperial eagle occur in intensively managed hunting areas, where such predators are frequently persecuted. The conservation of these endangered predators may require efforts to increase rabbit densities in protected areas.     

Extensive Genetic Polymorphism of Peptidases A, B, C, and D, in Wild Rabbit (Oryctolagus cuniculus) Populations from the Iberian Peninsula

Genetic polymorphism of peptidases A, B, C, andD in the wild rabbit (Oryctolagus cuniculus)was assessed by horizontal starch gel electrophoresis,in a total of 1003 individuals comprising 16 Iberian and 2 French populations and 1 domestic stock.Twenty-five different alleles were detected: 6 in PEPA,4 in PEPB, 8 in PEPC, and 7 in PEPD. The highest valuesof observed heterozygosity averaged over the four loci were obtained for the southwestern Iberianpopulations and a clinal loss of variability in anortheastern direction was detected. A clear separationbetween the two putative subspecies O. c.cuniculus and O. c. algirus was notobtained.     

A 7-bp insertion in the 3' untranslated region suggests the duplication and concerted evolution of the rabbit SRY gene

In this work we report the genetic polymorphism of a 7-bp insertion in the 3'' untranslated region of the rabbit SRY gene. The polymorphic GAATTAA motif was found exclusively in one of the two divergent rabbit Y-chromosomal lineages, suggesting that its origin is more recent than the separation of the O. c. algirus and O. c. cuniculus Y-chromosomes. In addition, the remarkable observation of haplotypes exhibiting 0, 1 and 2 7-bp inserts in essentially all algirus populations suggests that the rabbit SRY gene is duplicated and evolving under concerted evolution.     

Contrasting patterns of introgression at X-linked loci across the hybrid zone between subspecies of the European rabbit (Oryctolagus cuniculus)

Hybrid zones provide an excellent opportunity for studying the consequences of genetic changes between closely related taxa. Here we investigate patterns of genetic variability and gene flow at four X-linked loci within and between the two subspecies of European rabbit (Oryctolagus cuniculus cuniculus and O. c. algirus). Two of these genes are located near the centromere and two are located near the telomeres. We observed a deep split in the genealogy of each gene with the root located along the deepest branch in each case, consistent with the evolution of these subspecies in allopatry. The two centromeric loci showed low levels of variability, high levels of linkage disequilibrium, and little introgression between subspecies. In contrast, the two telomeric loci showed high levels of variability, low levels of linkage disequilibrium, and considerable introgression between subspecies. These data are consistent with suppression of recombination near the centromere of the rabbit X chromosome. These observations support a view of speciation where genomic incompatibilities at different loci in the genome create localized differences in levels of gene flow between nascent species.     

Oryctolagus giberti n. sp. (Lagomorpha, Mammalia) du Pléistocène inférieur de Cueva Victoria (Murcie, Espagne)

A new species of leporid (Lagomorpha, Mammalia) is described based on the material from the Lower Pleistocene of Cueva Victoria (Murcia, Spain). This species, named Oryctolagus giberti n. sp., presents several intermediate characters between the first known representative of the genus, O. laynensis, from the Middle Pliocene of Spain, and the modern European rabbit (O. cuniculus), specially at the level of the palate, the mandible, the third lower premolar, the ulna, the coxal and the femur. In this respect, Oryctolagus giberti n. sp. is a firm candidate to occupy the vacant phylogenetic place between these two species. It also questions the results of molecular studies that date the divergence of groups A and B of Oryctolagus cuniculus at about 2 Myr ago.     

First description of sarcoptic mange in wild European rabbit (Oryctolagus cuniculus)

The European wild rabbit (Oryctolagus cuniculus) is a keystone species in the Iberian Mediterranean ecosystems being the staple prey of protected and endangered predators. Wild rabbits were once widespread, but the introduction of two viral diseases, myxomatosis in the 1950s and the rabbit hemorrhagic disease in 1989, resulted into a dramatic decline of its populations. Sarcoptic mange is a highly contagious parasitic infection caused by Sarcoptes scabiei. The first cases of sarcoptic mange in a wild rabbit population are recorded from a hunting area in Majorca (Balearic Islands, Spain). Five of 32 inspected rabbits (14.7%) were affected, with similar prevalences in summer and autumn. Sarcoptes scabiei were obtained from the edge of the lesions of two of the rabbits. The most frequently observed lesions were small areas of alopecia and crusts in the limbs. Affected limbs presented also a marked hypertrichosis and an apparent abnormal growth of the nails. One rabbit also presented lesions around mouth and nostrils. Parasitized rabbits were in significant lower body condition than healthy ones. According to previous reports of mange epizooties in other naïve wild species that caused marked short-term effects in their populations, this disease may pose a risk for the conservation of wild rabbit and their predators. Whether mange is endemic in rabbits from Majorca or it has been introduced is unknown. Game managers are encouraged to be aware of introducing sarcoptic mange during rabbit translocations.

     

Worldwide migration of parasitic mites as a result of bumblebee commercialization

We investigated the status of infestation by a tracheal mite, Locustacarus buchneri, in natural populations of a Japanese native bumblebee species, Bombus hypocrita, collected on Hokkaido Island and in the Aomori prefecture between 1997 and 2001. We also investigated mite infestation in commercial colonies of the European bumblebee, Bombus terrestris, imported from the Netherlands and Belgium, and the Japanese native species, B. ignitus, imported from the Netherlands, between 1997 and 2001. We detected mites in natural populations of the two B. hypocrita subspecies and in the commercial colonies. Analysis of variations in 535 bp sequences of the mitochondrial cytochrome oxidase subunit 1 (CO1) gene showed that the mite haplotypes in the native populations and in the imported colonies did not overlap in 1997–1999, but in 2000–2001 some mites possessing European CO1 haplotypes were detected in the natural populations of Japanese native bumblebees. In addition, many mites possessing Japanese haplotypes were detected in the imported commercial colonies from Europe. Considering the fact that the Japanese native bumblebees, B. hypocrita, were once exported to Europe for commercialization, these results suggest that bumblebee commercialization has caused overseas migration and cross-infestation of parasitic mites among natural and commercial colonies. However, because the Japanese and European CO1 haplotypes were closely related, there was a possibility that the European haplotypes found in the mites in the Hokkaido Island revealed native variation. To clarify the status of mite invasion, further detailed analysis of genetic variation of the mite, using other genetic markers on additional samples, need to be performed.