Phylogeography of Korean raccoon dogs: implications of peripheral isolation of a forest mammal in East Asia

Genetic tracing of market wildlife and viruses at the epicenter of the COVID-19 pandemic

Zoonotic spillovers of viruses have occurred through the animal trade worldwide. The start of the COVID-19 pandemic was traced epidemiologically to the Huanan Wholesale Seafood Market, the site with the most reported wildlife vendors in the city of Wuhan, China. Here, we analyze publicly available qPCR and sequencing data from environmental samples collected in the Huanan market in early 2020. We demonstrate that the SARS-CoV-2 genetic diversity linked to this market is consistent with market emergence, and find increased SARS-CoV-2 positivity near and within a particular wildlife stall. We identify wildlife DNA in all SARS-CoV-2 positive samples from this stall. This includes species such as civets, bamboo rats, porcupines, hedgehogs, and one species, raccoon dogs, known to be capable of SARS-CoV-2 transmission. We also detect other animal viruses that infect raccoon dogs, civets, and bamboo rats. Combining metagenomic and phylogenetic approaches, we recover genotypes of market animals and compare them to those from other markets. This analysis provides the genetic basis for a short list of potential intermediate hosts of SARS-CoV-2 to prioritize for retrospective serological testing and viral sampling.

Impacts of Quaternary glaciation, geological history and geography on animal species history in continental East Asia: A phylogeographic review

Continental East Asia has a mild Pleistocene climate and a complex recent geological history. Phylogeographic studies of animals over the last 30 years have produced several distinctive patterns. Glaciation refugia are numerous and are not restricted to any particular regions. Most of them are localized and species-specific, although several large refugia, for example the mountains of SW China, are shared by multiple species and have refugia-within-refugia. Furthermore, postglaciation range expansion events vary greatly in time, scale and direction. Large-scale south-to-north post-LGM expansions are few and mostly occurred in the northern regions. Additionally, several unique geographic features, including the three-step terrain of China and the northern arid belt, have significant impacts on many species histories. Overall, the impacts of Pleistocene glaciations, particularly the LGM, on species history vary drastically from nondetectable to significant. The impacts are the least for species from the southwestern region and are most dominant for species from the north. Geological events play a more significant role in shaping species history than Pleistocene climatic changes. Phylogeographic patterns among animals species are highly consistent with those of plants. Future phylogeographic endeavour in East Asia should be hypothesis-driven and seek processes that underlie common patterns. The wide use of genomic data allow accurate estimates of historical population processes and exploration of older history beyond the Pleistocene.

Genome-scale phylogeography resolves the native population structure of the Asian longhorned beetle, Anoplophora glabripennis (Motschulsky)

Human-assisted movement has allowed the Asian longhorned beetle (ALB, Anoplophora glabripennis (Motschulsky)) to spread beyond its native range and become a globally regulated invasive pest. Within its native range of China and the Korean peninsula, human-mediated dispersal has also caused cryptic translocation of insects, resulting in population structure complexity. Previous studies used genetic methods to detangle this complexity but were unable to clearly delimit native populations which is needed to develop downstream biosurveillance tools. We used genome-wide markers to define historical population structure in native ALB populations and contemporary movement between regions. We used genotyping-by-sequencing to generate 6102 single-nucleotide polymorphisms (SNPs) and amplicon sequencing to genotype 53 microsatellites. In total, we genotyped 712 individuals from ALB's native distribution. We observed six distinct population clusters among native ALB populations, with a clear delineation between northern and southern groups. Most of the individuals from South Korea were distinct from populations in China. Our results also indicate historical divergence among populations and suggest limited large-scale admixture, but we did identify a restricted number of cases of contemporary movement between regions. We identified SNPs under selection and describe a clinal allele frequency pattern in a missense variant associated with glycerol kinase, an important enzyme in the utilization of an insect cryoprotectant. We further demonstrate that small numbers of SNPs can assign individuals to geographic regions with high probability, paving the way for novel ALB biosurveillance tools.

Phylogeographic study of the Bufo gargarizans species complex, with emphasis on Northeast Asia

We conduct a phylogeographic and population genetic study of the Asiatic toad (Bufo gargarizans) to understand its evolutionary history, and the influence of geology and climate. A total of 292 individuals from 94 locations were genotyped for two mitochondrial loci (cytb, ND2) and five nuclear introns (Sox9-2, Rho-3, CCNB2-3, UCH-2, and DBI-2), and we performed a suite of phylogenetic, population genetic, and divergence dating analyses. The phylogenetic trees constructed using mitochondrial loci inferred B. gargarizans being divided into two major groups: China mainland and Northeast Asia (Northeast China, Russia, and Korean Peninsula). As with previous studies of this species, we recover population genetic structure not tied to geographic region. Additionally, we discover a new genetic clade restricted to Northeast Asia that points towards the Korean Peninsula being a glacial refugium during the Pleistocene. The weak phylogeographic pattern of B. gargarizans is likely the result of multiple biological, anthropogenic, and historical factors – robust dispersal abilities as a consequence of physiological adaptations, human translocation, geologic activity, and glacial cycles of the Pleistocene. We highlight the complex geologic and climatic history of Northeast Asia and encourage further research to understand its impact on the biodiversity in the region.

Update on Distribution and Conservation Status of Amphibians in the Democratic People's Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties

Determining the range, status, ecology and behaviour of species from areas where surveys and samplings are uncommon or difficult to conduct is a challenge, such as in the Democratic People's Republic of Korea (DPR Korea). Here, we used genetic samples, field surveys, call recordings, photographic identification and a literature review to estimate the presence, range and status of amphibians in the DPR Korea. From our combined results and based on the IUCN Red List categories and criteria, we were able to estimate the national threat levels for most species. Our results demonstrated the presence of 18 native species and the suspected presence of Karsenia koreana and two Onychodactylus species. We reported the first record for Rana uenoi in the vicinity of Pyongyang using molecular tools and similarly confirmed the presence of Dryophytes japonicus at the same location. Based on distribution and modelling, we can expect the contact zone between species within the Rana and Onychodactylus genera to be located along the Changbai Massif, a mountain range that marks a shift in ecoregions and acts as a barrier to dispersion. The species richness was higher in the lowlands and at lower latitudes, with such areas populated by up to 11 species, while more northern regions were characterised by species richness of about half of that value. The combination of ecological models and known threats resulted in the recommendation of ten species as threatened at the national level following the IUCN Red List categories and criteria. This high number of threatened species was anticipated based on the high threat level to amphibians in bordering nations and globally. While the ecology of species in the DPR Korea is still understudied, we argue that species relying on agricultural wetlands such as rice paddies are not under imminent threat due to the enduring presence of extensive agricultural landscapes with low rates of chemical use and mechanisation. The maintenance of such landscapes is a clear benefit to amphibian species, in contrast to more industrialised agricultural landscapes in neighbouring nations. In comparison, the status of species dependent on forested habitats is unclear and threat levels are likely to be higher because of deforestation, as in neighbouring nations.

Evolution of MHC class I genes in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Major histocompatibility complex (MHC) genes have been widely studied to assess the immunological fitness and evolutionary adaptation of animal populations. Among the Canidae, the raccoon dog's adventurous nature, omnivorous behavior, and high variability of intracellular pathogens make it ideal to study selection on MHC class I in a non-model canid species. Here, we examined allelic diversity and evolutionary patterns of MHC class I genes in the raccoon dog (Nyctereutes procyonoides). We identified 48 novel MHC class I alleles from 31 raccoon dogs from Japan and Russia. Some alleles were geographically restricted, whereas others were widely distributed across the species' range. The rate of non-synonymous substitutions was greater than that of synonymous substitutions for both exon 2 and exon 3 encoding α1 and α2 domains, respectively, in the α chain of the MHC class I protein. Positively selected sites at the amino acid level were evident in both the α1 and α2 domains, and a recombination breakpoint was found in exon 3. Bayesian phylogenetic trees showed no evidence of trans-species polymorphism (TSP) with alleles from carnivoran species in other families but did detect TSP between raccoon dogs and the domestic dog, Canis familiaris, indicative of long-term balancing selection in canids. Our results indicate that the extensive allelic diversity of MHC class I in Japanese and Russian raccoon dogs has been influenced and maintained by pathogen-driven positive selection, recombination, and long-term balancing selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13364-021-00561-y.

Phylogenetic relationships between different raccoon dog (Nyctereutes procyonoides) populations based on four nuclear and Y genes

BACKGROUND

The raccoon dog (Nyctereutes procyonoides), endemic to East Asia, is classified as six subspecies according to their geographical distribution including a population introduced to Europe. Studies on phylogenetic relationship or population genetics in both native and introduced areas have been carried out recently. Lately, opinions that Japanese raccoon dogs should be classified as a different species were asserted based on several studies using karyotypes, morphometric characters, mtDNA, and microsatellites analysis. However, no data pertaining to the nuclear DNA (nDNA) or Y chromosome are available.

OBJECTIVE

To estimate the relationship among the species using different genes is necessary in understanding of the history of this species.

METHOD

Therefore, we investigated nDNA and Y chromosomes in our study to define relationships: (1) between continental raccoon dog populations, (2) between original and introduced groups, and (3) between continental and Japanese groups.

RESULTS

The analysis of four nuclear (CHRNA1, VTN, TRSP, WT1) and ZFY genes indicated that there had been no genetic differentiation among the continental populations. However, significant differences were observed between continental and Japanese raccoon dogs in VTN and ZFY genes implying genetic differentiation has been going between them.

CONCLUSION

To better understand the phylogenetic relationship among raccoon dog populations, further study will be necessary.

Molecular evolution of MHC class II DRB exon 2 in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Raccoon dogs, Nyctereutes procyonoides, are native to East Asia, but have been introduced into western Russia and eastern Europe. To determine allelic diversity and elucidate the evolution of major histocompatibility complex (MHC) class II genes in the raccoon dog, we analysed a 237-bp region of DRB exon 2 from 36 individuals of native and introduced populations from Japan and Russia. We detected 23 DRB alleles (Nypr-DRBs), 22 of which were novel. Some alleles were found across the species’ range, while others were geographically restricted. For both native and introduced populations, the ratio of non-synonymous to synonymous substitution rates for codons at predicted antigen-binding sites was significantly greater than 2, indicating that Nypr-DRBs have evolved under positive selection. Mixed effect model evolution analysis and an algorithm to detect recombination showed five positively selected codons and one recombination breakpoint, respectively. Overall, our results suggest that the diversity of MHC class II DRB in N. procyonoides was influenced and maintained by recombination, pathogen-driven positive selection, geographical barriers and the founder effect. A Bayesian phylogenetic tree revealed no evidence of trans-species polymorphism (TSP), but instead showed monophyly for the Nypr-DRB alleles within a larger clade of canid sequences. The lack of TSP may have been due to long-term divergence of raccoon dogs from other canids, or to their having encountered different sets of pathogens due to occupying a different ecological niche.

Comparison of feeding habits and habitat use between invasive raccoons and native raccoon dogs in Hokkaido, Japan

Background

In Japan, invasive raccoons cause severe ecological and social problems by transmitting pathogens to humans, livestock, and native species, causing substantial crop damage, and competing with native species. Possible competition between invasive raccoons and native raccoon dogs is of concern in Japan because Japanese raccoon dogs have a limited distribution and are native only to Japan and the two species have similar characteristics. We assessed potential competition between raccoons and raccoon dogs by comparing feeding habits and habitat use.

Results

Both species were captured in Hokkaido, Japan from 2004 to 2017. More raccoons were captured close to agricultural land at the forest periphery (70.1%, 358/511); conversely, more raccoon dogs were captured in the forest core (74.9%, 253/338). Feeding habits were then examined by fecal analysis and stable isotope analyses. Fecal analysis revealed both species to be opportunistic omnivores that consumed easily found food items. However, raccoon feces contained more crops, whereas raccoon dog feces contained more insects, reflecting the different locations in which the species were trapped. Moreover, stable isotope ratios were significantly higher in raccoons than raccoon dogs (Corn has the highest carbon stable isotope (δ¹³C) value, and amphibians and reptiles are high in nitrogen stable isotope (δ¹⁵N); forest resources such as insects and wild fruits are low in δ¹³C and δ¹⁵N).

Conclusions

We conclude that both species ate similar food types, but their food preferences appeared to differ. Raccoon and raccoon dog habitat use also differed, possibly because the two species inhabited areas where they could easily obtain their preferred foods. Therefore, the current feeding habits and habitat use of raccoons do not appear to overlap sufficiently with those of raccoon dogs to impact the latter. The results of this study, particularly the stable isotope data, may provide a useful precedent for future studies of competition in medium-sized mammals, particularly canids.

What is the taxonomic status of East Asian otter species based on molecular evidence?: focus on the position of the Japanese otter holotype specimen from museum

The Japanese otter (Lutra nippon), once inhabited in most islands of Japan, is now considered as an extinct species. Although the Japanese otter is regarded as a distinct species from the Eurasian otter (L. lutra), its phylogeny and taxonomic status are based on limited information on morphological and genetic data, and thus further clarification is required. Here, we assessed the phylogenetic relationship among the genus Lutra and taxonomic status of L. nippon by using the complete sequences of cytochrome b gene of its holotype. The present phylogenic trees supported that the genus Lutra specimens largely formed monophyletic group, with L. sumatrana as a basal to other Lutra species. Within Lutra species, L. nippon was distantly related with L. lutra. The European otter population of L. l. lutra were clustered together with its subspecies, L. l. chinensis rather than the same subspecies, Korean otter population. The discrepancy between the genetic data and traditional taxonomy justifies the necessity of reexamination of the current subspecific classification system of Eurasian otters. Level of genetic divergence between the holotype of L. nippon and L. lutra was two to three-fold lower than those among the other sister species of the Lutrinae. Based on the level of divergence between the L. nippon and L. lutra, and insufficient evidence of morphological difference between them, it is suggested that designation of Japanese otter as a separate species from L. lutra will be reconsidered.

Genetic Diversity and Population Structure of East Asian Raccoon Dog (Nyctereutes procyonoides): Genetic Features in Central and Marginal Populations

The raccoon dog (Nyctereutes procyonoides) is endemic to East Asia but has been introduced in Europe. Its high adaptability enabled its rapid colonization of European countries, where population growth has been raising concerns regarding ecosystem disturbance and the spread of zoonotic diseases. The genetic diversity and structure of endemic, source, and introduced populations from seven locations across South Korea, China, Russian Far East, Finland (spread to Finland after introduction to European part of Russia from Russian Far East), Vietnam, and Japan (Honshu and Hokkaido) were examined based on 16 microsatellite loci. Two major and significantly different (F_ST = 0.236) genetic clusters were found: continental (South Korean, Chinese, Russian, Finnish, and Vietnamese) and island (Japanese) populations. The continental raccoon dog population comprises three subpopulations (Chinese_Russian_Finnish, South Korean, and Vietnamese) and the Japanese population consists of Honshu and Hokkaido subpopulations. The genetic diversity and geographic structure of raccoon dogs in East Asia has been influenced by natural barriers to gene flow and reveals a typical central-marginal trend in genetic diversity (continental vs. island, and central vs. marginal or source vs. introduced within continental populations). The detected differences between continental and island populations agree with those reported in previous studies that considered these populations as different species.

Deep Phylogeographic Breaks in Magpie Pica pica Across the Holarctic: Concordance with Bioacoustics and Phenotypes

We examined sequences of the mitochondrial control region in magpies (Pica pica) from the entire distribution range and found deep genetic splits into four major lineages: (1) group West (Europe-Siberia), (2) group East (southern Far East), (3) P. p. mauritanica (North Africa), and (4) P. p. hudsonia (North America). These lineages show a geographic pattern corresponding to known subspecies or subspecies groups. Genetic variation within the widely-distributed group West is low and neutrality tests supported a recent expansion scenario. The haplotypes from Kamchatka, representing a separated sublineage with clear affinity to the European-Siberian group, are almost identical, implying a recent bottleneck. Group East contained two subclades without clear geographic pattern, presumably due to admixing of populations that had diverged in Pleistocene refuges. The homogeneity of the Kyushu population supports historical reports of introduction of the species from Korea. In contrast, the high variation in the recently established Hokkaido population may reflect an ongoing invasion from several populations of the Far Eastern mainland. Bioacoustic data based on chatter call differentiate groups of subspecies and reflect phylogeographic patterns, i.e., mitochondrial lineages. Furthermore, we report the fast spreading of P. p. jankowskii towards the west along the upper Amur River, and a slower shifting of P. p. leucoptera in the opposite direction thus yielding a new contact zone. Overall, our data support a scenario of divergence in geographic isolation, but the ongoing expansion of distribution ranges may lead to major changes in phylogeographic patterns.

Biogeography of Korea's top predator, the yellow-throated Marten: evolutionary history and population dynamics

BACKGROUND

Peninsulas often harvest high genetic diversity through repeated southward migrations of species during glacial maxima. Studies addressing within-species evolutionary responses to climate fluctuations in northeast Asia are limited compared to other regions of the world, and more so in the Korean Peninsula. In this study, we conducted the first population-level study of the yellow-throated marten, Martes flavigula, from the Korean Peninsula, Russian, Taiwanese and Chinese localities in a biogeographic framework using mitochondrial (cyt-b, nd2, cr) and nuclear gene sequencing (ghr).

RESULTS

Bayesian analyses revealed a rather young population, with a split from the most recent common ancestor at around 125 kya. Martes flavigula likely colonized the Korean Peninsula from Mainland China through the Yellow Sea twice, ca. 60 kya and 20 kya. Korean martens diversified during the Late Pleistocene with at least two dispersal events out of Korea, towards the southwest to Taiwan (ca. 80 kya) and towards the North into Russia and eastern China; most likely after the Last Glacial Maxima (ca. 20 kya). We argue that the lack of population structure and mixed populations is possibly a consequence of the high dispersal capability of the species. The Bayesian skyline plot revealed a population decline within the last 5000 years, suggesting potential negative biotic and anthropogenic effects in the area. We find that local populations are not genetically differentiated, therefore no perceptible population structure within Korea was found.

CONCLUSIONS

The topography and geography of the Korean Peninsula has played a pivotal role in its colonization. Connections between the Korean Peninsula and the Mainland through sea-level drops of the Yellow Sea at times of glacial maxima and the high dispersal capability of M. flavigula adds to the lack of geographical structure in this species and the paraphyly of Korean lineages.

Identification of novel haplotypes and interpretation of gene flow of mitochondrial DNA control region of Eurasian otter (Lutra lutra) for the effective conservation

The number and distribution of Eurasian otters have declined during twentieth century due to human activity and water pollution. The global conservation status of Eurasian otter is presently ‘Near Threatened (NT)’ and strictly protected by being listed on the international legislation and conventions. A number of studies using the mitochondrial DNA (mtDNA) control region (CR) have been conducted in order to effectively apply conservation and reintroduction programs, especially in Europe. However, aside from Europe, there have been few studies concerning genetic diversity and phylogeny of Eurasian otters. Therefore, in this study, we sequenced partial mtDNA CR sequences (232 bp) from five South Korean Eurasian otters and analyzed 27 otters originating from parts of northeast Asia (South Korea, China, Japan and Russia (Sakhalin)), and Europe. Out of 232 bp partial mtDNA CR sequences, 13 polymorphic sites (5.6%) were identified and 4 novel mtDNA CR haplotypes (Lut16–19) were discovered from 12 Eurasian otters originating from northeast Asian region. In this study, a comprehensive analysis of genetic diversity and population structure of Eurasian otter between Europe and northeast Asia continents were conducted. Of these, different past demographic histories in Pleistocene period might have largely impacted the genetic structure of each population differently. In addition, low degree of gene flow, isolation by distance (IBD) pattern from geographically wide distanced dataset and analysis of molecular variance (AMOVA) also represented distinct genetic characteristics of Eurasian otter between Europe and northeast Asia.

Population structure of the raccoon dog (Nyctereutes procyonoides) using microsatellite loci analysis in South Korea: Implications for disease management

The prevention and control of infectious diseases transmitted by wildlife are gaining importance. To establish effective management strategies, it is essential to understand the population structure of animals. Raccoon dogs (Nyctereutes procyonoides) in South Korea play a key role in the maintenance of food web stability and possess genetic compositions that are unique compared to those in other areas. However, wild raccoon dogs play another role as the main host of various infectious diseases. To establish long-term strategies for disease management, we investigated the genetic structure and possible geographic barriers that influence the raccoon dog population in South Korea by analyzing 16 microsatellite loci. The present study showed that mountains were the major factors responsible for genetic structuring, along with distance. We proposed potential management units (MUs) for raccoon dogs based on the genetic structuring and gene-flow barrier data obtained in this study. Four MUs were suggested for the Korean raccoon dog population (Northern, Central, Southwestern, and Southeastern). The Korean raccoon dog population structure determined in this study and the proposed MUs will be helpful to establish pragmatic strategies for managing Korean raccoon dog population and for preventing the transmission of infectious diseases.

Phylogeographic and population insights of the Asian common toad (Bufo gargarizans) in Korea and China: population isolation and expansions as response to the ice ages

The effects of ice ages on speciation have been well documented for many European and North American taxa. In contrast, very few studies have addressed the consequences of such environmental and topographical changes in North East Asian species. More precisely, the Korean Peninsula offers a unique model to assess patterns and processes of speciation as it hosts the northern- and eastern-most distribution limit of some widespread Asian taxa. Despite this, studies addressing phylogeographic patterns and population genetics in the peninsula and surrounding countries are few and studies for most families are lacking. Here we inferred the phylogenetic relationships of the common toad (Bufo gargarizans) from South Korea and their North East Asian counterpart populations, based on mitochondrial data. Korean B. gargarizans GenBank BLASTs matched few individuals from nearby China, but the presence of a Korean clade suggests isolation on the Korean Peninsula, previous to the last glacial maximum, linked to sea level resurgence. Molecular clock calibrations within this group were used to date the divergence between clades and their relationship to paleo-climatic events in the area. Lack of genetic structure among South Korean populations and strong homogeneity between the Korean and some Chinese localities suggest weak isolation and recent expansion. Geographical projection of continuous coalescent maximum-clade-credibility trees shows an original Chinese expansion towards the Korean Peninsula through the Yellow Sea circa two million years ago with colonisation events dating circa 800 thousand years ago (K. y. a.). Following this colonisation, the data point to outgoing Korean Peninsula dispersal events throughout different periods, towards the North through land, and West through land bridge formations over the Yellow Sea during sea level falls. In accordance, demographic analyses revealed a population expansion in the Koran Peninsula circa 300 K. y. a., likely attributed to glacial cycle fluctuations.

Population Structure of the Raccoon Dog on the Grounds of the Imperial Palace, Tokyo, Revealed by Microsatellite Analysis of Fecal DNA

The raccoon dog (Nyctereutes procyonoides, Canidae, Carnivora) is highly adaptable to urban environments. Populations of carnivorans inhabiting urban areas sometimes differ ecologically and genetically from those in rural areas. However, there is little information on urban raccoon dogs. This study focused on raccoon dog populations in Tokyo, one of the most highly urbanized cities in the world. We examined the genotypes of 10 microsatellites for 101 fecal samples from raccoon dogs inhabiting the grounds of the Imperial Palace, a green space in central Tokyo. We successfully genotyped 58 samples originating from 31 individuals. We also analyzed muscle tissue samples from raccoon dogs from the grounds of the Imperial Palace, the Akasaka Imperial Grounds (a green space close to the Imperial Palace), and the surrounding urban area, and then investigated the genetic structure and diversity of these populations, and the genetic differentiation among them. The population on the grounds of the Imperial Palace was genetically differentiated from that in the Akasaka Imperial Grounds, suggesting that the roads and buildings act as barriers to gene flow. In addition, the population on the grounds of the Imperial Palace showed greater genetic difference from that in the surrounding area than that in the Akasaka Imperial Grounds. We speculate that the moats around the Imperial Palace restrict individual ranges within the palace grounds and limit migration and gene flow to other areas.

Intraspecific comparison of complete mitogenome sequences from two Asian raccoon dogs (Canidae: Nyctereutes procyonoides)

We determined the complete mitochondrial genome (GenBank accession number: KF709435) of the Korean raccoon dog Nyctereutes procyonoides koreensis and compared it with a previously published mitogenome (GenBank accession number: GU256221) of a Chinese raccoon dog. The total length of N. p. koreensis mitogenome is 16,802 bp, with a base composition of 32.1% A, 26.9% T, 26.8% C and 14.2% G. High similarity of 98.7% was found between the complete mitogenome sequences of Korean and Chinese raccoon dogs. Sequence similarity of the two mitogenomes was 99.3% in the other gene regions except for D-loop. The sequence similarity of 99.1% was found in the 13 protein-coding gene regions, whereas 99.6% was identical in mtDNA regions covering all the 22 tRNA genes. There was no variation between 12S rRNAs, whereas 0.5% difference was found between 16S rRNAs.