Phylogeny and phylogeography of the sword-tailed newt, Cynops ensicauda (Amphibia: Caudata), as revealed by nucleotide sequences of mitochondrial DNA

The emergence of a cryptic lineage and cytonuclear discordance through past hybridization in the Japanese fire-bellied newt, Cynops pyrrhogaster (Amphibia: Urodela)

Discrepancies in geographic variation patterns between nuclear DNA and mitochondrial DNA (mtDNA) are the result of the complicated differentiation processes in organisms and are the key to understanding their true evolutionary processes. The genetic differentiation of the northern and Southern-Izu lineages of the Japanese newt, Cynops pyrrhogaster, was investigated through their single nucleotide polymorphism variations obtained via multiplexed ISSR genotyping by sequencing (MIG-seq). We found three genetic groups (Tohoku, N-Kanto and S-Kanto), that were not detected by mtDNA variations, in the northern lineage. N-Kanto has intermediate genetic characteristics between Tohoku and S-Kanto. The genetic groups are now moderately isolated from each other and have unique genetic characteristics. An estimation of the evolutionary history using the approximate Bayesian computation (ABC) approach suggested that Tohoku diverged from the common ancestor of S-Kanto and S-Izu. Then, S-Kanto and S-Izu split, and the recent hybridization between Tohoku and S-Kanto gave rise to N-Kanto. The origin of N-Kanto through the hybridization is relatively young and seems to be related to changes in the distributions of Tohoku and S-Kanto as a result of climatic oscillation in the Pleistocene. We conclude that the mitochondrial genome of S-Kanto was captured in Tohoku and that the original mitochondrial genome of Tohoku was entirely removed through hybridization.

Distribution and Genetic Diversity of the Amphibian Chytrid in Japan

While research on frog chytrid fungus Batrachochytrium dendrobatidis (Bd), an infectious disease that threatens amphibian diversity, continues to advance worldwide, little progress has been made in Japan since around 2010. The reason for this is, which we pointed out in 2009, that the origin of frog chytrid fungus may be in the East Asian region, including Japan based on the Bd ITS-DNA variation, and as few cases of mass mortality caused by this fungus have been observed in wild amphibian populations in Japan, the interest of the Japanese government and the general public in Bd has waned. However, we believe that organizing the data obtained so far in Japan and distributing the status of frog chytrid fungus in Japan to the world will provide useful insight for future risk management of this pathogen. We collected more than 5500 swab samples from wild amphibians throughout Japan from 2009 to 2010. Then, we investigated the infection status using the Nested-PCR method. We sequenced the obtained DNA samples and constructed a maximum-parsimony (MP) tree to clarify the phylogenetic diversity of Bd. We detected Bd infection in 11 (nine native and two alien) amphibian species in Japan and obtained 44 haplotypes of Bd ITS-DNA. The MP tree showed a high diversity of Bd strains in Japan, suggesting that some strains belong to Bd-GPL and Bd-Brazil. Except for local populations of the Japanese giant salamanders Andrias japonicus in Honshu Island and the sword tail newts Cynops ensicauda in Okinawa Island, the Bd infection prevalence in native amphibian species was very low. The alien bullfrog Aquarana catesbeiana had high Bd infection rates in all areas where they were sampled. No Bd infection was detected in other native amphibians in the areas where giant salamanders, sword tail newts, and bullfrogs were collected, suggesting that many native amphibians are resistant to Bd infection. The sword tail newt of Okinawa Island had both the highest infectious incidence and greatest number of haplotypes. The giant salamanders also showed relatively high infection prevalence, but the infected strains were limited to those specific to this species. These two Caudata species are endemic to a limited area of Japan, and it was thought that they may have been refugia for Bd, which had been distributed in Japan Islands for a long time.

Unique Evolution of Hyla hallowellii Among Amphibians of the Central Ryukyus, Japan (Anura: Hylidae)

Fauna of the Central Ryukyus includes a high percentage of endemic species, and Hallowell's tree frog Hyla hallowellii Thompson, 1912 is one of such elements, occurring in a total of eight islands in the Amami and Okinawa Island groups. Using samples representing all of these eight island populations, we studied variations in morphology, karyotype, allozyme, and mtDNA, to clarify the pattern of geographic differentiation of H. hallowellii and consider factors for its formation. We could not clearly discriminate one population from another in morphology, nor could we find any interpopulation difference in karyotype. From genetic analyses, using allozymes and cyt b, we found low overall differentiations among populations. However, the southern populations from Okinawajima and Yoronjima were genetically nearly identical with the northern Amamioshima population. From that group the geographically intermediate Tokunoshima and Kakeromajima populations showed prominent differentiations. These patterns of geographical differentiation greatly differ from those known in other amphibian species of the Central Ryukyus, and suggest that H. hallowellii has evolutionary history unique to the species.

Environmental DNA metabarcoding to detect pathogenic Leptospira and associated organisms in leptospirosis-endemic areas of Japan

Leptospires, which cause the zoonotic disease leptospirosis, persist in soil and aqueous environments. Several factors, including rainfall, the presence of reservoir animals, and various abiotic and biotic components interact to influence leptospiral survival, persistence, and pathogenicity in the environment. However, how these factors modulate the risk of infection is poorly understood. Here we developed an approach using environmental DNA (eDNA) metabarcoding for detecting the microbiome, vertebrates, and pathogenic Leptospira in aquatic samples. Specifically, we combined 4 sets of primers to generate PCR products for high-throughput sequencing of multiple amplicons through next-generation sequencing. Using our method to analyze the eDNA of leptospirosis-endemic areas in northern Okinawa, Japan, we found that the microbiota in each river shifted over time. Operating taxonomic units corresponding to pathogenic L. alstonii, L. kmetyi, and L. interrogans were detected in association with 12 nonpathogenic bacterial species. In addition, the frequencies of 11 of these species correlated with the amount of rainfall. Furthermore, 10 vertebrate species, including Sus scrofa, Pteropus dasymallus, and Cynops ensicauda, showed high correlation with leptospiral eDNA detection. Our eDNA metabarcoding method is a powerful tool for understanding the environmental phase of Leptospira and predicting human infection risk.

The taxonomic position and the unexpected divergence of the Habu viper, Protobothrops among Japanese subtropical islands

There are four Habu species currently recognized in Japan: Protobothrops flavoviridis from the Amami Islands and the Okinawa Islands, P. tokarensis from the Tokara Islands, P. elegans from the Yaeyama Islands and Ovophis okinabvensis from the Amami Islands and the Okinawa Islands. To clarify their taxonomic positions, we determined the complete mitochondria genome sequence (approx. 17kb) from two specimens from two different islands each for P. flavoviridis, P. tokarensis and P. elegans as well as one specimen of O. okinavensis and reconstructed the molecular phylogeny of Protobothrops using the published sequences of related species. The maximum likelihood tree showed four major species groups within Protbothrops: Group I consisting of P. cornutus, P. dabieshanensis, P. jerdonii and P. xiangchengensis; Group II consisting of P. flavoviridis and P. tokarensis; Group III consisting of P. maolensis, P. mucrosquamatus and P. elegans; Group IV consisting of P. himalayanus and P. kaubacki. Since we observed an unexpected divergence and the paraphyly of the two samples of P. flavoviridis collected from different islands, Amami-Oshima and Okinawajima within the Group II, we expanded the analysis by increasing the number of P. flavoviridis and P. tokarensis collected from 10 islands: Amami-Oshima (5 specimens), Kakeromajima (4) and Tokunoshima (4) from the Amami Islands, Okinawajima (4), Iheyajima (4), Iejima (4), Tokashikijima (4) and Kumejima (4) from the Okinawa Islands, Kodakarajima (P. tokarensis) (4) and Takarajima (P. tokarensis) (4) from the Tokara Islands. The maximum likelihood tree of the 44 samples replicated the significant divergence of P. flavoviridis between the Amami Clade including Amami-Oshima, Kakeromajima and Tokunoshima and the Okinawa Clade including Okinawajima, Iheyajima, Iejima, Tokashikijima and Kumejima. The Amami Clade also include all specimens from the Tokara Islands currently known as an independent species, P. tokarensis, suggesting the paraphyly of the taxon, P. flavoviridis. In contrast, we observed a distinct lineage of the two specimens from the Yaeyama Islands, supporting the validity of the taxon, P. elegans as an independent species. By MCMC method, we estimated the divergence time between the Amami Clade and the Okinawa Clade to be 6.51MYA, suggesting that the vicariance of the two clades preceded the geological separation of the Amami Islands and the Okinawa Islands (∼1.5MYA). As expected from the limited mobility of terrestrial reptiles including snakes, we observed high genetic divergence in Habu mtDNA among Japanese subtropical island populations.

Contrasting genetic structure of closely related giant water bugs: phylogeography of Appasus japonicus and Appasus major (Insecta: Heteroptera, Belostomatidae)

Appasus japonicus and A. major, two belostomatid species of the giant water bug found in parts of East Asia, have very similar morphological characteristics and ecological niches, and also overlapping habitats. However, the results of our previous published study utilizing molecular phylogenetic analyses of mitochondrial and nuclear DNA revealed extensive genetic differences, which indicated that the possibility of inter-specific hybridization was extremely unlikely. We collected A. japonicus and A. major from the Japanese Archipelago, Korean Peninsula, and Russian Far East, and conducted molecular analyses of mitochondrial DNA COI and 16S rRNA to compare phylogenetic relationships between these species. Three major clades were recognized within A. japonicus. Specimens from the Korean population constituted a monophyletic clade, and were a sister group of the western region of the Japanese Archipelago. The Eastern Japanese clade was clearly differentiated. Four major clades were recognized within A. major. Specimens from the Japanese and Korean populations revealed two distinct monophyletic clades. Significant differentiation was clearly observed between their genetic structures. Furthermore, the results of mismatch distribution and Bayesian skyline plot analyses suggested the possibility of a bottleneck effect or founder effect in two of the A. major clades. Collectively, these results demonstrated both similarities and differences in these two species even though their distribution widely overlaps in East Asia, with their morphological characteristics and ecological niches being very similar. These differences in genetic structures are considered to be due to their evolutionary history.

Phylogeny and historical demography of Cynops pyrrhogaster (Amphibia: Urodela): taxonomic relationships and distributional changes associated with climatic oscillations

We investigated the phylogenetic relationships and estimated the historical demography of the Japanese fire-bellied newt, Cynops pyrrhogaster, from Japanese mainlands using 1407-bp sequences of the mitochondrial DNA (NADH6, tRNAglu, cyt b) and 1208-bp sequences of nuclear DNA (Rag-1) genes. Phylogenetic trees based on mitochondrial DNA revealed four major haplotype clades (NORTHERN, CENTRAL, WESTERN, and SOUTHERN clades) within this species. Degree of genetic differentiation among major haplotype clades was very large for intraspecific variation, suggesting this species to be composed of four species lineages that replace each other geographically. Nuclear genetic variation presented no obvious patterns of geographic structure except for the distinctness of populations diagnosed by NORTHERN clade of mitochondrial haplotypes, suggesting results of incomplete lineage sorting. Current distribution and estimated divergence times for the genus Cynops suggest that the common ancestor of two Japanese species (C. pyrrhogaster and C. ensicauda from the Ryukyu Islands) had diverged at the edge of the continent corresponding to the present East China Sea and Central Ryukyus. Subsequent range expansion to Japanese mainland seems to have occurred in the middle Miocene. Population-genetic analyses indicated that all species lineages, except for the SOUTHERN one, experienced geographic population reductions and expansions associated with glacial and postglacial climatic oscillations.

Mitochondrial genomes and divergence times of crocodile newts: inter-islands distribution of Echinotriton andersoni and the origin of a unique repetitive sequence found in Tylototriton mt genomes

Crocodile newts, which constitute the genera Echinotriton and Tylototriton, are known as living fossils, and these genera comprise many endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered taxa, we determined the complete nucleotide sequences of the mt genomes of the Japanese crocodile newt Echinotriton andersoni and Himalayan crocodile newt Tylototriton verrucosus. Although the control region (CR) is known as the most variable mtDNA region in many animal taxa, the CRs of crocodile newts are highly conservative. Rather, the genes of NADH dehydrogenase subunits and ATPase subunit 6 were found to have high sequence divergences and to be usable for population genetics studies. To estimate the inter-population divergence ages of E. andersoni endemic to the Ryukyu Islands, we performed molecular dating analysis using whole and partial mt genomic data. The estimated divergence ages of the inter-island individuals are older than the paleogeographic segmentation ages of the islands, suggesting that the lineage splits of E. andersoni populations were not caused by vicariant events. Our phylogenetic analysis with partial mt sequence data also suggests the existence of at least two more undescribed species in the genus Tylototriton. We also found unusual repeat sequences containing the 3' region of cytochrome apoenzyme b gene, whole tRNA-Thr gene, and a noncoding region (the T-P noncoding region characteristic in caudate mtDNAs) from T. verrucosus mtDNA. Similar repeat sequences were found in two other Tylototriton species. The Tylototriton taxa with the repeats become a monophyletic group, indicating a single origin of the repeat sequences. The intra-and inter-specific comparisons of the repeat sequences suggest the occurrences of homologous recombination-based concerted evolution among the repeat sequences.

Phylogeny and biogeography of the Anderson's crocodile newt, Echinotriton andersoni (Amphibia: Caudata), as revealed by mitochondrial DNA sequences

The Anderson's crocodile newt, Echinotriton andersoni, is considered a relic and endangered species distributed in the Central Ryukyus. To elucidate phylogenetic relationships and detailed genetic structures among populations, we analyzed variation in the mitochondrial cytochrome b gene. Results strongly support a primary dichotomy between populations from the Amami and Okinawa Island Groups with substantial genetic divergence, favoring a primary divergence between the two island groups. Within the latter, populations from the southern part of Okinawajima Island are shown to be more closely related to those from Tokashikijima Island than to those from the northern and central parts of Okinawajima. The prominent genetic divergence between the two island groups of the Central Ryukyus seems to have initiated in the Miocene, i.e., prior to formation of the strait that has consistently separated these island groups since the Pleistocene. The ancestor of the southern Okinawajima-Tokashikijima is estimated to have migrated from the northern and central parts of Okinawajima into southern Okinawajima at the Pleistocene, and dispersed into Tokashikijima subsequently.