Evolutionary Relationships and Reproductive Isolating Mechanisms in the Rice Frog (Fejervarya limnocharis) Species Complex from Sri Lanka, Thailand, Taiwan and Japan, Inferred from mtDNA Gene Sequences, Allozymes, and Crossing Experiments

Postmating isolation and evolutionary relationships among Fejervarya species from Lesser Sunda, Indonesia and other Asian countries revealed by crossing experiments and mtDNA Cytb sequence analyses

To evaluate the degree of postmating isolation and the evolutionary relationships among frog species in the genus Fejervarya from Indonesia (Lesser Sunda), Bangladesh, China, and Japan, crossing experiments and molecular phylogenetic analyses were carried out. Crossing experiments revealed that reciprocal hybrids among F. iskandari, F. verruculosa, and F. sp. large type and between F. multistriata and F. kawamurai are viable through metamorphosis, while those between the F. iskandari group and F. limnocharis group were completely or partially inviable at the tadpole stage and those between Southeast Asian and South Asian Fejervarya groups were completely inviable at the embryonic stage. The mature reciprocal hybrids between F. iskandari and F. verruculosa from Lesser Sunda, Indonesia, showed some degree of abnormality in spermatogenesis. In phylogenetic analyses based on mtDNA Cytb sequences, F. iskandari formed a sister clade with F. verruculosa from Lesser Sunda, Indonesia, with 8.1% sequence divergence. F. multistriata from China formed a clade with populations of F. limnocharis in Thailand, Malaysia, and Indonesia (topotype), and these taxa showed sister relationships to F. kawamurai from Japan with 8.9% sequence divergence. Fejervarya sp. small type from Bangladesh formed a clade with the other South Asian members of the Fejervarya group and formed a sister clade with the Southeast Asian Fejervarya group, with 23.1% sequence divergence in the Cytb gene. These results showed that the degree of postmating isolation reflects molecular phylogenetic relationships and that F. iskandari and F. verruculosa from Indonesia (Lesser Sunda) are reproductively isolated by abnormalities in spermatogenesis and show genetic differentiation.

Multiple Mechanistic Action of Brevinin-1FL Peptide against Oxidative Stress Effects in an Acute Inflammatory Model of Carrageenan-Induced Damage

Amphibian skin is acknowledged to contain an antioxidant system composed of various gene-encoded antioxidant peptides, which exert significant effects on host defense. Nevertheless, recognition of such peptides is in its infancy so far. Here, we reported the antioxidant properties and underlying mechanism of a new antioxidant peptide, brevinin-1FL, identified from Fejervarya limnocharis frog skin. The cDNA sequence encoding brevinin-1FL was successfully cloned from the total cDNA of F. limnocharis and showed to contain 222 bp. The deduced mature peptide sequence of brevinin-1FL was FWERCSRWLLN. Functional analysis revealed that brevinin-1FL could concentration-dependently scavenge ABTS⁺, DPPH, NO, and hydroxyl radicals and alleviate iron oxidation. Besides, brevinin-1FL was found to show neuroprotective activity by reducing contents of MDA and ROS plus mitochondrial membrane potential, increasing endogenous antioxidant enzyme activity, and suppressing H₂O₂-induced death, apoptosis, and cycle arrest in PC12 cells which were associated with its regulation of AKT/MAPK/NF-κB signal pathways. Moreover, brevinin-1FL relieved paw edema, decreased the levels of TNF-α, IL-1β, IL-6, MPO, and malondialdehyde (MDA), and restored catalase (CAT) and superoxide dismutase (SOD) activity plus glutathione (GSH) contents in the mouse injected by carrageenan. Together, these findings indicate that brevinin-1FL as an antioxidant has potent therapeutic potential for the diseases induced by oxidative damage. Meanwhile, this study will help us further comprehend the biological functions of amphibian skin and the mechanism by which antioxidants protect cells from oxidative stress.

Cytonuclear discordance, reticulation and cryptic diversity in one of North America's most common frogs

Complicated phylogenetic histories benefit from diverse sources of inference. Pseudacris crucifer (spring peeper) spans most of eastern North America and comprises six mtDNA lineages that form multiple contact zones. The putative Miocene or early Pliocene origins of the oldest lineages within Pseudacris crucifer imply sufficient time for species-level divergence. To understand why this species appears unified while congeners have radiated, we analyze and compare male advertisement calls, mitochondrial, and nuclear markers and speak to the complex processes that have potentially influenced its contemporary patterns. We find extensive geographic and topological mitonuclear discordance, with three nuclear lineages containing 6 more-structured mtDNA lineages, and nuclear introgression at some contact zones. Male advertisement call differentiation is incongruent with the genetic structure as only one lineage appears differentiated. Occupying the Interior Highlands of the central United States, this Western lineage also has the most concordant mitochondrial and nuclear geographic patterns. Based on our findings we suggest that the antiquity of common ancestors was not as important as the maintenance of allopatry in the divergence in P. crucifer genetic lineages. We use multiple lines of evidence to generate hypotheses of isolation, reticulation, and discordance within this species and to expand our understanding of the early stages of speciation.

Filling the BINs of life: Report of an amphibian and reptile survey of the Tanintharyi (Tenasserim) Region of Myanmar, with DNA barcode data

Despite threats of species extinctions, taxonomic crises, and technological advances in genomics and natural history database informatics, we are still distant from cataloguing all of the species of life on earth. Amphibians and reptiles are no exceptions; in fact new species are described nearly every day and many species face possible extinction. The number of described species continues to climb as new areas of the world are explored and as species complexes are examined more thoroughly. The use of DNA barcoding provides a mechanism for rapidly estimating the number of species at a given site and has the potential to record all of the species of life on Earth. Though DNA barcoding has its caveats, it can be useful to estimate the number of species in a more systematic and efficient manner, to be followed in combination with more traditional, morphology-based identifications and species descriptions. Herein, we report the results of a voucher-based herpetological expedition to the Tanintharyi (Tenasserim) Region of Myanmar, enhanced with DNA barcode data. Our main surveys took place in the currently proposed Tanintharyi National Park. We combine our results with photographs and observational data from the Chaung-nauk-pyan forest reserve. Additionally, we provide the first checklist of amphibians and reptiles of the region, with species based on the literature and museum. Amphibians, anurans in particular, are one of the most poorly known groups of vertebrates in terms of taxonomy and the number of known species, particularly in Southeast Asia. Our rapid-assessment program combined with DNA barcoding and use of Barcode Index Numbers (BINs) of voucher specimens reveals the depth of taxonomic diversity in the southern Tanintharyi herpetofauna even though only a third of the potential amphibians and reptiles were seen. A total of 51 putative species (one caecilian, 25 frogs, 13 lizards, 10 snakes, and two turtles) were detected, several of which represent potentially undescribed species. Several of these species were detected by DNA barcode data alone. Furthermore, five species were recorded for the first time in Myanmar, two amphibians (Ichthyophis cf. kohtaoensis and Chalcorana eschatia) and three snakes (Ahaetulla mycterizans, Boiga dendrophila, and Boiga drapiezii).

Changes of diet and dominant intestinal microbes in farmland frogs

BACKGROUND

Agricultural activities inevitably result in anthropogenic interference with natural habitats. The diet and the gut microbiota of farmland wildlife can be altered due to the changes in food webs within agricultural ecosystems. In this work, we compared the diet and intestinal microbiota of the frog Fejervarya limnocharis in natural and farmland habitats in order to understand how custom farming affects the health of in vivo microbial ecosystems.

RESULTS

The occurrence, abundance, and the numbers of prey categories of stomach content were significantly different between the frogs inhabiting natural and farmland habitats. In addition, differences in the abundance, species richness, and alpha-diversity of intestinal microbial communities were also statistically significant. The microbial composition, and particularly the composition of dominant microbes living in intestines, indicated that the land use practices might be one of factors affecting the gut microbial community composition. Although the first three dominant microbial phyla Bacteroidetes, Firmicutes, and Proteobacteria found in the intestines of frogs were classified as generalists among habitats, the most dominant gut bacterial phylum Bacteroidetes in natural environments was replaced by the microbial phylum Firmicutes in farmland frogs. Increased intestinal microbial richness of the farmland frogs, which is mostly contributed by numerous microbial species of Proteobacteria, Actinobacteria, Acidobacteria, and Planctomycetes, not only reflects the possible shifts in microbial community composition through the alteration of external ecosystem, but also indicates the higher risk of invasion by disease-related microbes.

CONCLUSIONS

This study indicates that anthropogenic activities, such as the custom farming, have not only affected the food resources of frogs, but also influenced the health and in vivo microbial ecosystem of wildlife.

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group)

Background

Hybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group.

Results

We found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30 km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts.

Conclusions

Our findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0385-2) contains supplementary material, which is available to authorized users.

Inferring the degree of incipient speciation in secondary contact zones of closely related lineages of Palearctic green toads (Bufo viridis subgroup)

Reproductive isolation between lineages is expected to accumulate with divergence time, but the time taken to speciate may strongly vary between different groups of organisms. In anuran amphibians, laboratory crosses can still produce viable hybrid offspring >20 My after separation, but the speed of speciation in closely related anuran lineages under natural conditions is poorly studied. Palearctic green toads (Bufo viridis subgroup) offer an excellent system to address this question, comprising several lineages that arose at different times and form secondary contact zones. Using mitochondrial and nuclear markers, we previously demonstrated that in Sicily, B. siculus and B. balearicus developed advanced reproductive isolation after Plio-Pleistocene divergence (2.6 My, 3.3-1.9), with limited historic mtDNA introgression, scarce nuclear admixture, but low, if any, current gene flow. Here, we study genetic interactions between younger lineages of early Pleistocene divergence (1.9 My, 2.5-1.3) in northeastern Italy (B. balearicus, B. viridis). We find significantly more, asymmetric nuclear and wider, differential mtDNA introgression. The population structure seems to be molded by geographic distance and barriers (rivers), much more than by intrinsic genomic incompatibilities. These differences of hybridization between zones may be partly explained by differences in the duration of previous isolation. Scattered research on other anurans suggests that wide hybrid zones with strong introgression may develop when secondary contacts occur <2 My after divergence, whereas narrower zones with restricted gene flow form when divergence exceeds 3 My. Our study strengthens support for this rule of thumb by comparing lineages with different divergence times within the same radiation.

Cryptic anuran biodiversity in Bangladesh revealed by mitochondrial 16S rRNA gene sequences

To survey the diversity of anuran species in Bangladesh, we compared mitochondrial 16S rRNA gene sequences (approximately 1.4 kbp) from 107 Bangladesh frog specimens. The results of genetic divergence and phylogenetic analyses incorporating data from related species revealed the occurrence of at least eight cryptic species. Hoplobatrachus tigerinus from two districts diverged considerably, indicating the involvement of a cryptic species. Two Fejervarya sp. (large and medium types) and Hylarana cf. taipehensis formed lineages distinct from related species and are probably new species. Microhyla cf. ornata differed from M. ornata with respect to type locality area and involved two distinct species. In addition, we found that Hylarana sp. and Microhyla sp. did not match congeners examined to date in either morphology or 16S rRNA sequence. The occurrence of M. fissipes was tentatively suggested. Consequently, at least, 19 species were found from Bangladesh in this study. These findings revealed a rich anuran biodiversity in Bangladesh, which is unexpected considering the rather simple topographic features of the country.

A new species of the Fejervarya limnocharis complex from Japan (Anura, Dicroglossidae)

We describe a new species of dicroglossid frog of the Fejervarya limnocharis complex from western Honshu, Japan Mainland. The new species, Fejervarya kawamurai, is genetically closer to F. sakishimensis than to F. limnocharis. It differs from F. sakishimensis by smaller tympanum, head, forelimb, hindlimb, foot, and tibia lengths, all relative to snout-vent length, and from F. multistriata by relatively shorter forelimb, hindlimb, foot, and tibia. From F. limnocharis and F. iskandari, it is differentiated by relatively smaller forelimb, hindlimb, foot, and tibia lengths. Taxonomic problems of Fejervarya populations occurring in Central Ryukyus, continental China, and Taiwan are discussed.

Irreversibility of a bad start: early exposure to osmotic stress limits growth and adaptive developmental plasticity

Harsh environments experienced early in development have immediate effects and potentially long-lasting consequences throughout ontogeny. We examined how salinity fluctuations affected survival, growth and development of Fejervarya limnocharis tadpoles. Specifically, we tested whether initial salinity effects on growth and rates of development were reversible and whether they affected the tadpoles' ability to adaptively accelerate development in response to deteriorating conditions later in development. Tadpoles were initially assigned to either low or high salinity, and then some were switched between salinity levels upon reaching either Gosner stage 30 (early switch) or 38 (late switch). All tadpoles initially experiencing low salinity survived whereas those initially experiencing high salinity had poor survival, even if switched to low salinity. Growth and developmental rates of tadpoles initially assigned to high salinity did not increase after osmotic stress release. Initial low salinity conditions allowed tadpoles to attain a fast pace of development even if exposed to high salinity afterwards. Tadpoles experiencing high salinity only late in development metamorphosed faster and at a smaller size, indicating an adaptive acceleration of development to avoid osmotic stress. Nonetheless, early exposure to high salinity precluded adaptive acceleration of development, always causing delayed metamorphosis relative to those in initially low salinity. Our results thus show that stressful environments experienced early in development can critically impact life history traits, having long-lasting or irreversible effects, and restricting their ability to produce adaptive plastic responses.

Impacts of the herbicide butachlor on the larvae of a paddy field breeding frog (Fejervarya limnocharis) in subtropical Taiwan

Butachlor is the most commonly used herbicide on paddy fields in Taiwan and throughout Southeast Asia. Since paddy fields provide habitat for pond breeding amphibians, we examined growth, development, time to metamorphosis, and survival of alpine cricket frog tadpoles (Fejervarya limnocharis) exposed to environmentally realistic concentrations of butachlor. We documented negative impacts of butachlor on survival, development, and time to metamorphosis, but not on tadpole growth. The 96 h LC(50) for tadpoles was 0.87 mg/l, much lower than the 4.8 mg/l recommended dosage for application to paddy fields. Even given the rapid breakdown of butachlor, tadpoles would be exposed to concentrations in excess of their 96 h LC(50) for an estimated 126 h. We also documented DNA damage (genotoxicity) in tadpoles exposed to butachlor at concentrations an order of magnitude less than the 4.8 mg/l recommended application rate. We did not find that butachlor depressed cholinesterase activity of tadpoles, unlike most organophosphorus insecticides. We conclude that butachlor is likely to have widespread negative impacts on amphibians occupying paddy fields with traditional herbicide application.

Complete mitochondrial genomes and novel gene rearrangements in two dicroglossid frogs, Hoplobatrachus tigerinus and Euphlyctis hexadactylus, from Bangladesh

We determined the complete nucleotide sequences of mitochondrial (mt) genomes from two dicroglossid frogs, Hoplobatrachus tigerinus (Indian Bullfrog) and Euphlyctis hexadactylus (Indian Green frog). The genome sizes are 20462 bp in H. tigerinus and 20280 bp in E. hexadactylus. Although both genomes encode the typical 37 mt genes, the following unique features are observed: 1) the ND5 genes are duplicated in H. tigerinus that have completely identical sequences, whereas duplicated ND5 genes in E. hexadactylus possessed dissimilar substitutions; 2) duplicated control region (CR) in H. tigerinus has almost identical sequences whereas single control region (CR) was found in E. hexadactylus; 3) the tRNA-Leu (CUN) gene is translocated from the LTPF tRNA cluster to downstream of ND5-1 in H. tigerinus, and the tRNA-Pro gene is translocated from the LTPF tRNA cluster to downstream of CR in E. hexadactylus; 4) pseudo tRNA-Leu (CUN) and tRNA-Pro genes are observed in E. hexadactylus; and 5) two tRNA-Met genes are encoded in both species, as observed in the previously reported dicroglossid mt genomes. Almost all observed gene rearrangements in H. tigerinus and E. hexadactylus can be explained by the tandem duplication and random loss model, except translocation of tRNA-Pro in E. hexadactylus. The novel mt genomic features found in this study may be useful for future phylogenetic studies in the dicroglossid taxa. However, the mt genome with interesting features found in the present study reveal a high level of variation of gene order and gene content, inspiring more research to understand the mechanisms behind gene and genome evolution in the dicroglossid and as well as in the amphibian taxa in future studies.

Strong reproductive barriers in a narrow hybrid zone of West-Mediterranean green toads (Bufo viridis subgroup) with Plio-Pleistocene divergence

BACKGROUND

One key question in evolutionary biology deals with the mode and rate at which reproductive isolation accumulates during allopatric speciation. Little is known about secondary contacts of recently diverged anuran species. Here we conduct a multi-locus field study to investigate a contact zone between two lineages of green toads with an estimated divergence time of 2.7 My, and report results from preliminary experimental crosses.

RESULTS

The Sicilian endemic Bufo siculus and the Italian mainland-origin B. balearicus form a narrow hybrid zone east of Mt. Etna. Despite bidirectional mtDNA introgression over a ca. 40 km North-South cline, no F1 hybrids could be found, and nuclear genomes display almost no admixture. Populations from each side of the contact zone showed depressed genetic diversity and very strong differentiation (FST = 0.52). Preliminary experimental crosses point to a slightly reduced fitness in F1 hybrids, a strong hybrid breakdown in backcrossed offspring (F1 x parental, with very few reaching metamorphosis) and a complete and early mortality in F2 (F1 x F1).

CONCLUSION

Genetic patterns at the contact zone are molded by drift and selection. Local effective sizes are reduced by the geography and history of the contact zone, B. balearicus populations being at the front wave of a recent expansion (late Pleistocene). Selection against hybrids likely results from intrinsic genomic causes (disruption of coadapted sets of genes in backcrosses and F2-hybrids), possibly reinforced by local adaptation (the ranges of the two taxa roughly coincide with the borders of semiarid and arid climates). The absence of F1 in the field might be due to premating isolation mechanisms. Our results, show that these lineages have evolved almost complete reproductive isolation after some 2.7 My of divergence, contrasting sharply with evidence from laboratory experiments that some anuran species may still produce viable F1 offspring after > 20 My of divergence.

Genetic differentiation of the Fejervarya limnocharis complex from Bangladesh and other Asian countries elucidated by allozyme analyses

The present study was conducted to elucidate the genetic divergence and the phylogenetic relationships in the F. limnocharis complex from Bangladesh and other Asian countries such as Sri Lanka, Thailand, Malaysia, Taiwan and Japan by allozyme analyses. We used a total of 95 frogs of the F. limnocharis complex from these countries and F. cancrivora from the Philippines as an outgroup. Based on body size, the F. limnocharis complex from Bangladesh was divided into three distinct groups: large, medium and small types. Allozyme analyses were carried out with 28 loci encoding 20 enzymes and two blood proteins by horizontal starch-gel electrophoresis. When genetic distance was calculated, distinct divergence was found among the three types: mean genetic distance was 0.782 between the small and medium types, 1.458 between the large and medium types, and 1.520 between the large and small types. Phylogenetic trees based on genetic distance showed that all populations of Bangladesh small type strongly formed a cluster and were found to be most closely related to the Sri Lanka population; that all populations of Bangladesh large type formed a very strong cluster and were grouped with several populations from Thailand, Malaysia, Japan, and Taiwan; and that the medium type was segregated from all other groups. This may imply that each of the three types is a different species, and that the medium type is possibly an undescribed taxon.

Phylogeography of the rice frog, Fejervarya multistriata (Anura: Ranidae), from China based on mtDNA D-loop sequences

The rice frog, Fejervarya multistriata, is an amphibian widely distributed in China. In this study, we sampled the species across its distributional area in China and sequenced the mtDNA D-loop to investigate the genetic diversity and geographical pattern of the frog population. The results revealed 38 haplotypes in the population, with K2P values varying from 0.19% to 4.22%. Both a phylogenetic analysis and a nested clade analysis (NCA) detected two geographically isolated lineages respectively distributed around the Yangtze drainage (Yangtze lineage) and the south of China (southern lineage). NCA inferred a contiguous range expansion within the Yangtze lineage and allopatric fragmentation within the southern lineage, which might be partly due to the limited samples from this lineage. Accordingly, Fu's Fs test also indicated a population expansion after glacial movement. Therefore, we assumed that the species history responding to glacial events shaped the present population pattern of F. multistriata on the Chinese mainland.

Genetic divergence and evolutionary relationships in six species of genera Hoplobatrachus and Euphlyctis (Amphibia: Anura) from Bangladesh and other Asian countries revealed by mitochondrial gene sequences

To elucidate the species composition, genetic divergence, evolutionary relationships, and divergence time of Hoplobatrachus and Euphlyctis frogs (subfamily Dicroglossinae, family Ranidae) in Bangladesh and other Asian countries, we analyzed the mitochondrial Cyt b, 12S, and 16S rRNA genes of 252 specimens. Our phylogenetic analyses showed 13 major clades corresponding to several cryptic species as well as to nominal species in the two genera. The results suggested monophyly of Asian Hoplobatrachus species, but the position of African Hoplobatrachus occipitalis was not clarified. Nucleotide divergence and phylogenetic data suggested the presence of allopatric cryptic species allied to Euphlyctis hexadactylus in Sundarban, Bangladesh and several parapatric cryptic species in the Western Ghats, India. The presence of at least two allopatric cryptic species among diverged Euphlyctis cyanophlyctis in Bangladesh, India, and Sri Lanka was also suggested. In some cases, our estimated divergence times matched the paleogeological events of South and Southeast Asian regions that may have led to the divergence of Hoplobatrachus and Euphlyctis taxa. Especially, land formation at Bangladesh (15-10Ma) may have allowed the spread of these frog taxa to Southeast Asian areas, and the aridification of central India (5.1-1.6Ma) might have affected the gene flow of widely distributed species. The present study revealed prior underestimation of the richness of the amphibian fauna in this region, indicating the possible occurrence of many cryptic species among these groups.

Morphological divergence, reproductive isolating mechanism, and molecular phylogenetic relationships among Indonesia, Malaysia, and Japan populations of the Fejervarya limnocharis complex (Anura, Ranidae)

In order to elucidate the taxonomic status of the Fejervarya limnocharis complex relative to Malaysia and Japan populations, morphological observations and molecular phylogenetic analysis were carried out using three populations from Indonesia (type locality), Malaysia, and Japan. In addition, we conducted histological and spermatogenic observations using hybrids among these populations. Principal component and cluster analyses demonstrated that these populations could be clearly separated from one another. Abnormal testes were found in the hybrids between the Japan and Indonesia populations and between the Japan and Malaysia populations, but testes of the controls and hybrids between the Malaysia and Indonesia populations were quite normal. The mean number of univalents per cell was 5.42, 4.58, and 0.20 in hybrids between the Indonesia and Japan populations, Malaysia and Japan populations, and Indonesia and Malaysia populations, respectively. Sequence divergences in 16S rRNA and Cyt b genes were 0-0.4% (xbar=0.2%) and 0.3-1.5% (xbar=1.0%), respectively, between the Malaysia and Indonesia populations, and 2.4-2.6% (xbar=2.5%) and 11.0-12.0% (xbar=11.5%) between the Japan population and F. limnocharis complex, including the Malaysia and Indonesia populations and F. multistriata from China. This study indicated that the Malaysia population and F. multistriata from China should be designated as a subspecies of topotypic F. limnocharis, and that the Japan population should be regarded as a distinct species.