Microsatellite marker development by multiplex ion torrent PGM sequencing: a case study of the endangered Odorrana narina complex of frogs.

The endangered Ryukyu tip-nosed frog Odorrana narina and its related species, Odorrana amamiensis, Odorrana supranarina, and Odorrana utsunomiyaorum, belong to the family Ranidae and are endemically distributed in Okinawa (O. narina), Amami and Tokunoshima (O. amamiensis), and Ishigaki and Iriomote (O. supranarina and O. utsunomiyaorum) Islands. Because of varying distribution patterns, this species complex is an intrinsic model for speciation and adaptation. For effective conservation and molecular ecological studies, further genetic information is needed. For rapid, cost-effective development of several microsatellite markers for these and 2 other species, we used next-generation sequencing technology of Ion Torrent PGM™. Distribution patterns of repeat motifs of microsatellite loci in these modern frog species (Neobatrachia) were similarly skewed. We isolated and characterized 20 new microsatellite loci of O. narina and validated cross-amplification in the three-related species. Seventeen, 16, and 13 loci were cross-amplified in O. amamiensis, O. supranarina, and O. utsunomiyaorum, respectively, reflecting close genetic relationships between them. Mean number of alleles and expected heterozygosity of newly isolated loci varied depending on the size of each inhabited island. Our findings suggested the suitability of Ion Torrent PGM™ for microsatellite marker development. The new markers developed for the O. narina complex will be applicable in conservation genetics and molecular ecological studies.

Development and characterization of 14 microsatellite markers for Buergeria japonica (Amphibia, Anura, Rhacophoridae).

Buergeria japonica is a common frog species distributed throughout almost all islands in Ryukyu Archipelago. Because of their exceptionally wide distribution and higher physiological tolerance comparing to the other anurans, their demographic history and formation of distribution are intrinsic topics in the herpetological fauna of Ryukyu. Microsatellite marker is ideal genetic marker for such studies at inter- and intra-population level. We therefore developed microsatellite markers of B. japonica utilizing Ion PGM™ sequencing. As a result of the screening, we developed a total of 14 polymorphic markers. To test availabilities of these markers, we genotyped four island populations. The total number of alleles and expected hetelozygosities per locus ranged from 4 to 21 and 0.00 to 0.864, respectively. The phylogenetic relationship among the four populations based on the genetic distances of these markers was congruent with general divergence pattern of amphibians and reptiles in Ryukyu area. These markers developed in this study are considered to be useful for future studies about phylogeography and demography of this species.

Mitochondrial genomes of Japanese Babina frogs (Ranidae, Anura): unique gene arrangements and the phylogenetic position of genus Babina.

Genus Babina is a member of Ranidae, a large family of frogs, currently comprising 10 species. Three of them are listed as endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered species, we determined the complete nucleotide sequences of the mt genomes of 3 endangered Japanese Babina frogs, B. holsti, B. okinavana, and B. subaspera and 1 ranid frog Lithobates catesbeianus. The genes of NADH dehydrogenase subunit 5 (nad5) and the control region (CR) were found to have high sequence divergences and to be usable for population genetics studies. At present, no consensus on the phylogenetic position of genus Babina has been reached. To resolve this problem, we performed molecular phylogenetic analyses with the largest dataset used to date (11,345 bp from 2 ribosomal RNA- and 13 protein-encoding genes) in studies dealing with Babina phylogeny. These analyses revealed monophyly of Babina and Odorrana. It is well known that mt gene rearrangements of animals can provide usable phylogenetic information. Thus, we also compared the mt gene arrangements among Babina species and other related genera. Of the surveyed species, only L. catesbeianus manifested typical neobatrachian-type mt gene organization. In the B. okinavana, an additional pseudogene of tRNA-His (trnH) was observed in the CR downstream region. Furthermore, in the B. holsti and B. subaspera, the trnH/nad5 block was translocated from its typical position to the CR downstream region, and the translocated trnH became a pseudogene. The position of the trnH pseudogene is consistent with the translocated trnH position reported in Odorrana. Consequently, the trnH rearrangement seems to be a common ancestry characteristic (synapomorphy) of Babina and Odorrana. Based on the "duplication and deletion" gene rearrangement model, a single genomic duplication event can explain the order of derived mt genes found in Babina and Odorrana.

An Attempt at Captive Breeding of the Endangered Newt Echinotriton andersoni, from the Central Ryukyus in Japan.

Anderson's crocodile newt (Echinotriton andersoni) is distributed in the Central Ryukyu Islands of southern Japan, but environmental degradation and illegal collection over the last several decades have devastated the local populations. It has therefore been listed as a class B1 endangered species in the IUCN Red List, indicating that it is at high risk of extinction in the wild. The species is also protected by law in both Okinawa and Kagoshima prefectures. An artificial insemination technique using hormonal injections could not be applied to the breeding of this species in the laboratory. In this study we naturally bred the species, and tested a laboratory farming technique using several male and female E. andersoni pairs collected from Okinawa, Amami, and Tokunoshima Islands and subsequently maintained in near-biotopic breeding cages. Among 378 eggs derived from 17 females, 319 (84.4%) became normal tailbud embryos, 274 (72.5%) hatched normally, 213 (56.3%) metamorphosed normally, and 141 (37.3%) became normal two-month-old newts; in addition, 77 one- to three-year-old Tokunoshima newts and 32 Amami larvae are currently still growing normally. Over the last five breeding seasons, eggs were laid in-cage on slopes near the waterfront. Larvae were raised in nets maintained in a temperature-controlled water bath at 20 °C and fed live Tubifex. Metamorphosed newts were transferred to plastic containers containing wet sponges kept in a temperature-controlled incubator at 22.5 °C and fed a cricket diet to promote healthy growth. This is the first published report of successfully propagating an endangered species by using breeding cages in a laboratory setting for captive breeding. Our findings on the natural breeding and raising of larvae and adults are useful in breeding this endangered species and can be applied to the preservation of other similarly wild and endangered species such as E. chinhaiensis.

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.

Artificial production and natural breeding of the endangered frog species Odorrana ishikawae, with special reference to fauna conservation in the laboratory.

Odorrana ishikawae is listed as a class IB endangered species in the IUCN Red List and is protected by law in both Okinawa and Kagoshima Prefectures, Japan. Here, in an effort to help effectively preserve the genetic diversity of this endangered species in the laboratory, we tested a farming technique involving the artificial breeding of frogs, and also promoted natural breeding in the laboratory. Field-caught male/female pairs of the Amami and Okinawa Island populations were artificially bred using an artificial insemination method in the 2004, 2006, and 2008 breeding seasons (March to April). Although fewer than 50% of the inseminated eggs achieved metamorphosis, approximately 500, 300, and 250 offspring from the three respective trials are currently being raised in the laboratory. During the 2009 and 2010 breeding seasons, second-generation offspring were produced by the natural mating activities of the first offspring derived from the two artificial matings in 2004. The findings and the methods presented here appear to be applicable to the temporary protection of genetic diversity of local populations in which the number of individuals has decreased or the environmental conditions have worsened to levels that frogs are unable to survive by themselves.

Complete mitochondrial DNA sequence of the endangered frog Odorrana ishikawae (family Ranidae) and unexpected diversity of mt gene arrangements in ranids.

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of an endangered Japanese frog, Odorrana ishikawae (family Ranidae). We also sequenced partial mt genomes of three other Odorrana and six ranid species to survey the diversity of genomic organizations and elucidate the phylogenetic problems remaining in this frog family. The O. ishikawae mt genome contained the 37 mt genes and single control region (CR) typically found in vertebrate mtDNAs, but the region of Light-strand replication origin (OL) was triplicated in this species. Four protein-encoding genes (atp6, nd2, nd3, and nd5) were found to have high sequence divergence and to be usable for population genetics studies on this endangered species. Among the surveyed ranids, only two species (Rana and Lithobates) manifested the typical neobatrachian-type mt gene arrangement. In contrast, relatively large gene rearrangements were found in Amolops, Babina, and Staurois species; and translocations of single tRNA genes (trns) were observed in Glandirana and Odorrana species. Though the inter-generic and interspecific relationships of ranid taxa remain to be elucidated based on 12S and 16S rrn sequence data, some of the derived mt gene orders were found to have synapomorphic features useful for solving problematic ranid phylogenies. The tandem duplication and random loss (TDRL) model, the traditional model for mt gene rearrangement, failed to easily explain several of the mt gene rearrangements observed here. Indeed, the recent recombination-based gene rearrangement models seemed to be more suitable for this purpose. The high frequency of gene translocations involving a specific trn block (trnH-trnS1) and several single tRNA genes suggest that there may be a retrotranslocation in ranid mt genomes.