Phylogenetic relationships of the Chinese torrent frogs (Ranidae: Amolops) revealed by phylogenomic analyses of AFLP-Capture data.

The torrent frog genus Amolops contains nearly sixty species distributed in swift mountain streams throughout southeast Asia. The taxonomy of this genus has proven complicated due to unstable morphological diagnostic characters. The relationships of Amolops species and species groups were not readily resolved with a small number of molecular markers. Here, we applied the novel AFLP-Capture approach and acquired two large datasets (242 anonymous nuclear sequences and the mitochondrial genome) from 70 Chinese Amolops samples to study their relationships. The phylogenies inferred from the nuclear data and the mitochondrial data were both robust and revealed a primary phylogenetic split between eastern and western Chinese Amolops species. The relationships of the six species groups were clarified. While the three species groups in east China (the A. ricketti, A. daiyunensis and A. hainanensis groups) were monophyletic, the three species groups in the west (the A. mantzorum, A. monticola and A. marmoratus groups) were not monophyletic, suggesting a need for further investigation and revision. The robust phylogenies also provided new insights into species relationships, especially for the A. mantzorum group, which has been difficult to resolve due to multiple speciation events occurring approximately 7-8 million years ago. The divergence times estimated with the nuclear data indicated that the ancestor of the Chinese Amolops appeared in the late Eocene or early Oligocene, and that speciation events in the Chinese Amolops were often related to geological events (e.g. the uprising of mountains and the formation of islands). By including the mitochondrial sequences from GenBank, a more comprehensive Amolops phylogeny was constructed that reflected the origin of the Chinese Amolops. Based on all these results, a dispersal scenario of the torrent frogs was hypothesized. Our research serves as the first example of using AFLP-Capture to obtain a large amount of data for shallow-scale phylogenetic and taxonomic studies, which should be useful for other nonmodel organism groups.

Phylogeography, speciation and demographic history: Contrasting evidence from mitochondrial and nuclear markers of the Odorrana graminea sensu lato (Anura, Ranidae) in China.

Understanding the process of speciation and the factors driving the geographical distribution patterns of species is of great interest in ecology and evolutionary biology. Herein, we investigated the phylogeographic patterns, speciation, demographic history and genetic structure of the widespread endemic Odorrana graminea sensu lato in Southern China and adjacent areas. A total of 439 specimens from 68 localities were sequenced and analyzed for both mitochondrial (12S and 16S rRNA) and nuclear markers (RAG-1 and ten microsatellite loci). Phylogenetic analyses of the concatenated mtDNA data revealed five major highly divergent lineages within the O. graminea sensu lato in Southern China, and these divergent lineages were highly concordant with five geographical regions. The nuclear data showed a discordant genetic structure compared to the mtDNA lineages (Clades A, B, and C) for O. graminea sensu stricto, with an admixed pattern in the RAG-1 data and two structure clusters in the microsatellite data. The species delimitation analyses, based on three methods, supported the species status of Odorrana zhaoi and Odorrana rotodora, and revealed the existence of putative cryptic species in the O. graminea sensu stricto. In addition, one statistically significant gene flow event was detected from Clade B to Clade C based on mtDNA and RAG-1 data, and the microsatellite data suggested gene flow within the O. graminea sensu stricto. Bayesian skyline plotting analyses and ecological niche modeling supported demographic and range expansions during the LGM for Clades A and C of the O. graminea sensu stricto. In addition, ecological niche models suggested the existence of ecological divergence among the three Clades (Clades A, B and C) of the O. graminea sensu stricto. The intense uplifting of the Qinghai-Tibetan Plateau, the Quaternary climate oscillations, and drainage changes may have driven the speciation, genetic structure and phylogeoraphic patterns of the O. graminea sensu lato.

A Phylogenetically Distinct Group of Glandirana rugosa Found in Kyushu, Japan.

The Japanese wrinkled frog Glandirana rugosa is separated into five genetically different groups. One group in western Japan is further divided into three subgroups, found in Kyushu, Shikoku, and western Honshu. We collected G. rugosa frogs at 39 sites in Kyushu and determined nucleotide sequences of the mitochondrial 12S and 16S rRNA genes for phylogenetic analysis. Unexpectedly, we found a group of frogs in southeastern Kyushu that did not cluster with any of the pre-existing five groups of G. rugosa on the phylogenetic trees. The frogs in the new group and G. rugosa in Kyushu were externally similar, but there were a few significant differences in morphological features between the two populations. In addition, we observed significant differences in the frogs' calls . Thus, the group of the frogs in southeastern Kyushu may represent a new candidate species in the genus Glandirana. We discuss the possibility of a new species.

Population genomics of an exceptional hybridogenetic system of Pelophylax water frogs.

BACKGROUND: Hybridogenesis can represent the first stage towards hybrid speciation where the hybrid taxon eventually weans off its parental species. In hybridogenetic water frogs, the hybrid Pelophylax kl. esculentus (genomes RL) usually eliminates one genome from its germline and relies on its parental species P. lessonae (genomes LL) or P. ridibundus (genomes RR) to perpetuate in so-called L-E and R-E systems. But not exclusively: some all-hybrid populations (E-E system) bypass the need for their parental species and fulfill their sexual cycle via triploid hybrid frogs. Genetic surveys are essential to understand the great diversity of these hybridogenetic dynamics and their evolution. Here we conducted such study using RAD-sequencing on Pelophylax from southern Switzerland (Ticino), a geographically-isolated region featuring different assemblages of parental P. lessonae and hybrid P. kl. esculentus. RESULTS: We found two types of hybridogenetic systems in Ticino: an L-E system in northern populations and a presumably all-hybrid E-E system in the closely-related southern populations, where P. lessonae was not detected. In the latter, we did not find evidence for triploid individuals from the population genomic data, but identified a few P. ridibundus (RR) as offspring from interhybrid crosses (LR × LR). CONCLUSIONS: Assuming P. lessonae is truly absent from southern Ticino, the putative maintenance of all-hybrid populations without triploid individuals would require an unusual lability of genome elimination, namely that P. kl. esculentus from both sexes are capable of producing gametes with either L or R genomes. This could be achieved by the co-existence of L- and R- eliminating lineages or by "hybrid amphigamy", i. e. males and females producing sperm and eggs among which both genomes are represented. These hypotheses imply that polyploidy is not the exclusive evolutionary pathway for hybrids to become reproductively independent, and challenge the classical view that hybridogenetic taxa are necessarily sexual parasites.

Chromosomal Evolution in the Amolops mantzorum Species Group (Ranidae; Anura) Narrated by Repetitive DNAs.

In an attempt to analyze the organization of repetitive DNAs in the amphibian genome, 7 microsatellite motifs and a 5S rDNA sequence were synthesized and mapped in the karyotypes of 5 Amolops species. The results revealed nonrandom distribution of the microsatellite repeats, usually in the heterochromatic regions, as found in other organisms. These microsatellite repeats showed rapid changes among Amolops species, documenting the recent evolutionary history within this lineage. In contrast, 5S rDNA was localized in chromosomes 5 of all species, suggesting that these chromosomes are homologous within the monophyletic clade. Furthermore, the heteromorphic X and Y sex chromosomes (chromosomes 5) of A.mantzorum, had identical patterns of 5S rDNA, indicating that the subtelocentric Y resulted from a pericentric inversion. Several microsatellite repeats were found in the heteromorphic sex chromosomes, verifying the association of repetitive DNAs with sex chromosome differentiation in A. mantzorum.

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada, California, USA.

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross-breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short-lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.

Post-glacial colonization routes coincide with a life-history breakpoint along a latitudinal gradient.

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post-glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical QST -FST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The QST estimates between the colonization routes were four times higher than FST , indicating a prominent role for natural selection. QST within the colonization routes did not generally differ from FST , but we found temperature-dependent adaptive divergence close to the contact zone. These results indicate that lineage-specific variation can account for much of the adaptive divergence along a latitudinal gradient.

Parallel identification of novel antimicrobial peptide sequences from multiple anuran species by targeted DNA sequencing.

BACKGROUND: Antimicrobial peptides (AMPs) are multifunctional effector molecules that often combine direct antimicrobial activities with signaling or immunomodulatory functions. The skin secretions of anurans contain a variety of such bioactive peptides. The identification of AMPs from frog species often requires sacrificing several specimens to obtain small quantities of crude peptides, followed by activity based fractionation to identify the active principles. RESULTS: We report an efficient alternative approach to selectively amplify AMP-coding transcripts from very small amounts of tissue samples, based on RNA extraction and cDNA synthesis, followed by PCR amplification and high-throughput sequencing of size-selected amplicons. This protocol exploits the highly conserved signal peptide region of the AMP precursors from Ranidae, Hylidae and Bombinatoridae for the design of family-specific, forward degenerate primers, coupled with a reverse primer targeting the mRNA poly-A tail. CONCLUSIONS: Analysis of the assembled sequencing output allowed to identify more than a hundred full-length mature peptides, mostly from Ranidae species, including several novel potential AMPs for functional characterization. This (i) confirms the effectiveness of the experimental approach and indicates points for protocol optimization to account for particular cases, and (ii) encourages the application of the same methodology to other multigenic AMP families, also from other genera, sharing common features as in anuran AMPs.

Molecular phylogeny of Neotropical rock frogs reveals a long history of vicariant diversification in the Atlantic forest.

The Brazilian Atlantic coastal forest is one of the most heterogeneous morphoclimatic domains on earth and is thus an excellent region in which to examine the role that habitat heterogeneity plays in shaping diversification of lineages and species. Here we present a molecular phylogeny of the rock frogs of the genus Thoropa Cope, 1865, native to the Atlantic forest and extending to adjacent campo rupestre of Brazil. The goal of this study is to reconstruct the evolutionary history of the genus using multilocus molecular phylogenetic analyses. Our topology reveals 12 highly supported lineages among the four nominal species included in the study. Species T. saxatilis and T. megatympanum are monophyletic. Thoropa taophora is also monophyletic, but nested within T. miliaris. Populations of T. miliaris cluster in five geographically distinct lineages, with low support for relationships among them. Although all 12 lineages are geographically structured, some T. miliaris lineages have syntopic distributions with others, likely reflecting a secondary contact zone between divergent lineages. We discuss a biogeographic scenario that best explains the order of divergence and the distribution of species in Atlantic forest and adjacent areas, and outline the implications of our findings for the taxonomy of Thoropa.

Species delimitation with gene flow: A methodological comparison and population genomics approach to elucidate cryptic species boundaries in Malaysian Torrent Frogs.

Accurately delimiting species boundaries is a nontrivial undertaking that can have significant effects on downstream inferences. We compared the efficacy of commonly used species delimitation methods (SDMs) and a population genomics approach based on genomewide single-nucleotide polymorphisms (SNPs) to assess lineage separation in the Malaysian Torrent Frog Complex currently recognized as a single species (Amolops larutensis). First, we used morphological, mitochondrial DNA and genomewide SNPs to identify putative species boundaries by implementing noncoalescent and coalescent-based SDMs (mPTP, iBPP, BFD*). We then tested the validity of putative boundaries by estimating spatiotemporal gene flow (fastsimcoal2, ABBA-BABA) to assess the extent of genetic isolation among putative species. Our results show that the A. larutensis complex runs the gamut of the speciation continuum from highly divergent, genetically isolated lineages (mean Fst  = 0.9) to differentiating populations involving recent gene flow (mean Fst  = 0.05; Nm  > 5). As expected, SDMs were effective at delimiting divergent lineages in the absence of gene flow but overestimated species in the presence of marked population structure and gene flow. However, using a population genomics approach and the concept of species as separately evolving metapopulation lineages as the only necessary property of a species, we were able to objectively elucidate cryptic species boundaries in the presence of past and present gene flow. This study does not discount the utility of SDMs but highlights the danger of violating model assumptions and the importance of carefully considering methods that appropriately fit the diversification history of a particular system.

Diversity and phylogeography of Northeast Asian brown frogs allied to Rana dybowskii (Anura, Ranidae).

We investigated the species diversity and phylogeography of the Northeast Asian brown frogs allied to Rana dybowskii (the R. dybowskii species complex: R. dybowskii, R. pirica, and R. uenoi) using four mitochondrial and three nuclear loci. Phylogenetic analyses confirmed the existence of three distinct species in this complex; using extensive molecular data, we confirm the validity of Rana uenoi recognized as a distinct species, and infer R. dybowskii and R. pirica to be sister species. Also, we included populations from previously unsampled regions in Northeast China, and identified them to be R. dybowskii. While many species in Northeast Asia diverged due to Pleistocene glaciation, divergence-dating analyses inferred older, Miocene speciation in the R. dybowskii species complex. Ancestral area reconstruction identified the orogenic movement of the Changbai Mountain Range and the opening of the Sea of Japan/East Sea being major events influencing allopatric speciation.

Leapfrogging into new territory: How Mascarene ridged frogs diversified across Africa and Madagascar to maintain their ecological niche.

The Mascarene ridged frog, Ptychadena mascareniensis, is a species complex that includes numerous lineages occurring mostly in humid savannas and open forests of mainland Africa, Madagascar, the Seychelles, and the Mascarene Islands. Sampling across this broad distribution presents an opportunity to examine the genetic differentiation within this complex and to investigate how the evolution of bioclimatic niches may have shaped current biogeographic patterns. Using model-based phylogenetic methods and molecular-clock dating, we constructed a time-calibrated molecular phylogenetic hypothesis for the group based on mitochondrial 16S rRNA and cytochrome b (cytb) genes and the nuclear RAG1 gene from 173 individuals. Haplotype networks were reconstructed and species boundaries were investigated using three species-delimitation approaches: Bayesian generalized mixed Yule-coalescent model (bGMYC), the Poisson Tree Process model (PTP) and a cluster algorithm (SpeciesIdentifier). Estimates of similarity in bioclimatic niche were calculated from species-distribution models (maxent) and multivariate statistics (Principal Component Analysis, Discriminant Function Analysis). Ancestral-area reconstructions were performed on the phylogeny using probabilistic approaches implemented in BioGeoBEARS. We detected high levels of genetic differentiation yielding ten distinct lineages or operational taxonomic units, and Central Africa was found to be a diversity hotspot for these frogs. Most speciation events took place throughout the Miocene, including "out-of-Africa" overseas dispersal events to Madagascar in the East and to São Tomé in the West. Bioclimatic niche was remarkably well conserved, with most species tolerating similar temperature and rainfall conditions common to the Central African region. The P. mascareniensis complex provides insights into how bioclimatic niche shaped the current biogeographic patterns with niche conservatism being exhibited by the Central African radiation and niche divergence shaping populations in West Africa and Madagascar. Central Africa, including the Albertine Rift region, has been an important center of diversification for this species complex.

Spatiotemporal Diversification of the True Frogs (Genus Rana): A Historical Framework for a Widely Studied Group of Model Organisms.

True frogs of the genus Rana are widely used as model organisms in studies of development, genetics, physiology, ecology, behavior, and evolution. Comparative studies among the more than 100 species of Rana rely on an understanding of the evolutionary history and patterns of diversification of the group. We estimate a well-resolved, time-calibrated phylogeny from sequences of six nuclear and three mitochondrial loci sampled from most species of Rana, and use that phylogeny to clarify the group's diversification and global biogeography. Our analyses consistently support an "Out of Asia" pattern with two independent dispersals of Rana from East Asia to North America via Beringian land bridges. The more species-rich lineage of New World Rana appears to have experienced a rapid radiation following its colonization of the New World, especially with its expansion into montane and tropical areas of Mexico, Central America, and South America. In contrast, Old World Rana exhibit different trajectories of diversification; diversification in the Old World began very slowly and later underwent a distinct increase in speciation rate around 29-18 Ma. Net diversification is associated with environmental changes and especially intensive tectonic movements along the Asian margin from the Oligocene to early Miocene. Our phylogeny further suggests that previous classifications were misled by morphological homoplasy and plesiomorphic color patterns, as well as a reliance primarily on mitochondrial genes. We provide a phylogenetic taxonomy based on analyses of multiple nuclear and mitochondrial gene loci. [Amphibians; biogeography; diversification rate; Holarctic; transcontinental dispersal.

Identification and functional analyses of novel antioxidant peptides and antimicrobial peptides from skin secretions of four East Asian frog species.

In the present study, we identified 50 peptides that are classified into 21 peptide families with antioxidant and/or antimicrobial activity from Amolops daiyunensis, Pelophylax hubeiensis, Hylarana maosuoensis and Nanorana pleskei, which belong to four different genera in the Ranidae and Dicroglossidae families. These four frog species were found for the first time to express antioxidant peptides (AOPs) and antimicrobial peptides (AMPs). These peptides include seven newly discovered families daiyunin-1, daiyunin-2, daiyunin-3, maosonensis-1MS1, pleskein-1, pleskein-2, and pleskein-3. Antioxidant and antimicrobial activity assays showed that some of these peptides have good biological activities. For example, at a concentration of 50 µM, nigroain-B-MS1, and nigroain-C-MS1 both exhibited relatively strong 2,2-diphenyl-1-picrylhydrazyl (DPPH) and 2,2'-azinobis-(3-ethylbenzothiazoline-6-sulfonicacid) (ABTS) free radical scavenging ability, with eradication rates of 99.7% and 68.3% (nigroain-B-MS1), and 99.8% and 58.3% (nigroain-C-MS1), respectively. These peptides are potential candidates for the development of novel antioxidant or AMP preparations.

Analysis of fatty acid composition revealed differences in the diets of tadpoles of two amphibian species.

The study of the composition of fatty acid markers of tadpoles of cohabiting amphibian species for the first time revealed differences in their diets: the moor frog Rana arvalis prefers bacteria not associated with plant detritus, whereas the diet of the common spadefoot Pelobates fuscus is based on cyanobacteria, green algae, diatoms, and possibly higher plants. Major differences in the fatty acid composition are determined by the difference in the percentage of eicosapentaenoic and myristic acids.

[Identification of Ranae Oviductus original animal based on COI sequences].

Ranae Oviductus has a high economic and social value, but its adulterants are more numerous, which causes a great confusion to the market. Using DNA bar code technology based on COI sequence for PCR amplification and sequencing of the identified Rana dybowskii, R. chensinensis, R. huanrensis and R. amurensiss, the COI gene database of four species of Rana was established, and comparing the measured sequence with the sequence of GenBank, four kinds of Rana were identified. The MEGA (molecular evolutionary genetics analysis) 7 .0 software was used to calculate the genetic distance of K2P and construct the NJ (neighbor-joining) system cluster tree. The sequence of the four species of Rana measured were clustered into one group with the sequence of the four kinds of Rana downloaded from GenBank, but separated from the two outer groups downloaded from GenBank. The COI gene of the R. dybowskii was likely to have regional differences, however this technique failed to distinguish male and female Rana. The results showed that DNA bar code technology could accurately identify the base of original animal of R. oviductus. It indicates that DNA bar code COI provides a new method for the identification of R. oviductus.

Characterisation of major histocompatibility complex class I genes in Japanese Ranidae frogs.

The major histocompatibility complex (MHC) is a key component of adaptive immunity in all jawed vertebrates, and understanding the evolutionary mechanisms that have shaped these genes in amphibians, one of the earliest terrestrial tetrapods, is important. We characterised MHC class I variation in three common Japanese Rana species (Rana japonica, Rana ornativentris and Rana tagoi tagoi) and identified a total of 60 variants from 21 individuals. We also found evolutionary signatures of gene duplication, recombination and balancing selection (including trans-species polymorphism), all of which drive increased MHC diversity. A unique feature of MHC class I from these three Ranidae species includes low synonymous differences per site (d S) within species, which we attribute to a more recent diversification of these sequences or recent gene duplication. The resulting higher d N/d S ratio relative to other anurans studied could be related to stronger selection pressure at peptide binding sites. This is one of the first studies to investigate MHC in Japanese amphibians and permits further exploration of the polygenetic factors associated with resistance to infectious diseases.

Glaciations, gradients, and geography: multiple drivers of diversification of bush frogs in the Western Ghats Escarpment.

The historical processes underlying high diversity in tropical biodiversity hotspots like the Western Ghats of Peninsular India remain poorly understood. We sampled bush frogs on 13 massifs across the Western Ghats Escarpment and examined the relative influence of Quaternary glaciations, ecological gradients and geological processes on the spatial patterns of lineage and clade diversification. The results reveal a large in situ radiation (more than 60 lineages), exhibiting geographical structure and clade-level endemism, with two deeply divergent sister clades, North and South, highlighting the biogeographic significance of an ancient valley, the Palghat Gap. A majority of the bush frog sister lineages were isolated on adjacent massifs, and signatures of range stasis provide support for the dominance of geological processes in allopatric speciation. In situ diversification events within the montane zones (more than 1800 m) of the two highest massifs suggest a role for climate-mediated forest-grassland persistence. Independent transitions along elevational gradients among sub-clades during the Miocene point to diversification along the elevational gradient. The study highlights the evolutionary significance of massifs in the Western Ghats with the high elevations acting as centres of lineage diversification and the low- and mid-elevations of the southern regions, with deeply divergent lineages, serving as museums.

High genetic diversity but low population structure in the frog Pseudopaludicola falcipes (Hensel, 1867) (Amphibia, Anura) from the Pampas of South America.

Relative to South America's ecoregions, the temperate grasslands of the Pampas have been poorly studied from a phylogeographic perspective. Based on an intermediate biogeographic setting between subtropical forest (Atlantic Forest) and arid ecosystems (Chaco and Patagonia), Pampean species are expected to show unstable demographic histories due to the Quaternary climatic oscillations. Herein, we investigate the phylogenetic relatedness and phylogeographic history of Pseudopaludicola falcipes, a small and common frog that is widely distributed across the Pampean grasslands. First, we use molecular data to assess if P. falcipes represents a single or multiple, separately evolving cryptic lineages. Because P. falcipes is a small-size species (<20mm) with extensive coloration and morphological variation, we suspected that it might represent a complex of cryptic species. In addition, we expected strong genetic and geographic structuring within Pseudopaludicola falcipes due to its large geographic distribution, potentially short dispersal distances, and multiple riverine barriers. We found that P. falcipes is a single evolutionary lineage with poor geographic structuring. Furthermore, current populations of P. falcipes have a large effective population size, maintain ancestral polymorphisms, and have a complex network of gene flow. We conclude that the demographic history of P. falcipes, combined with its ecological attributes and the landscape features of the Pampas, favored a unique combination among anurans of small body size, large population size, high genetic variability, but high cohesiveness of populations over a wide geographic distribution.

Unusually high genetic diversity in the Bornean Limnonectes kuhlii-like fanged frogs (Anura: Dicroglossidae).

A fanged frog Limnonectes kuhlii was once thought to be wide-ranging in Southeast Asia, but is now confined to its type locality Java through recent phylogenetic studies, which clarified heterospecific status of non-Javanese populations, and monophyly of Bornean populations. However, large genetic differences among Bornean populations suggest occurrence of cryptic species, which we test using dense geographic sampling. We estimated the phylogenetic relationships among samples of Bornean populations together with their putative relatives from the continental Southeast Asia, using 2517bp sequences of the 12S rRNA, tRNA(val), and 16S rRNA of mitochondrial DNA, and 2367bp sequences of the NCX1, POMC, and RAG1 of nuclear genes. In the mtDNA trees, Bornean L. kuhlii-like frogs formed a monophyletic group split into 18 species lineages including L. hikidai, with the deepest phylogenetic split separating L. cintalubang from the remaining species. Almost all of these lineages co-occur geographically, and two to three lineages were found syntopically in each locality. Co-occurrence of more than one lineage may be maintained by differential morphology and microhabitat selection. These syntopic lineages should be regarded as distinct species. Our results clearly indicate that taxonomic revision is urgent to clarify many evolutionary problems of Bornean L. kuhlii-like frogs.