Congolius, a new genus of African reed frog endemic to the central Congo: A potential case of convergent evolution.

The reed frog genus Hyperolius (Afrobatrachia, Hyperoliidae) is a speciose genus containing over 140 species of mostly small to medium-sized frogs distributed in sub-Saharan Africa. Its high level of colour polymorphism, together with in anurans relatively rare sexual dichromatism, make systematic studies more difficult. As a result, the knowledge of the diversity and taxonomy of this genus is still limited. Hyperolius robustus known only from a handful of localities in rain forests of the central Congo Basin is one of the least known species. Here, we have used molecular methods for the first time to study the phylogenetic position of this taxon, accompanied by an analysis of phenotype based on external (morphometric) and internal (osteological) morphological characters. Our phylogenetic results undoubtedly placed H. robustus out of Hyperolius into a common clade with sympatric Cryptothylax and West African Morerella. To prevent the uncovered paraphyly, we place H. robustus into a new genus, Congolius. The review of all available data suggests that the new genus is endemic to the central Congolian lowland rain forests. The analysis of phenotype underlined morphological similarity of the new genus to some Hyperolius species. This uniformity of body shape (including cranial shape) indicates that the two genera have either retained ancestral morphology or evolved through convergent evolution under similar ecological pressures in the African rain forests.

The identity of emIndosylvirana nicobariensis/em (Stoliczka 1870) (Amphibia, Ranidae) from northeast India.

Stoliczka (1870) described Hylorana nicobariensis from the Nicobar Islands of India. The generic allocation of this enigmatic species is yet to be fully resolved (Chan et al. 2020b) and it has been placed in different genera (Boulenger 1885; Frost et al. 2006; Che et al. 2007; Oliver et al. 2015; Chan et al. 2020a; Chandramouli et al. 2020). As of now, the species is assigned to the genus Indosylvirana (Oliver et al. 2015; Chan et al. 2020b; Frost 2020). Apart from Nicobar archipelago, this species is distributed in the insular regions of southeast Asia (Oliver et al. 2015; Chandramouli et al. 2020; Frost 2020) and the states of Arunachal Pradesh, Assam, Mizoram, Tripura and West Bengal in northeast India (Sarkar et al. 1992, 2002; Sarkar Ray 2006; Ahmed et al. 2009; Mathew Sen 2010; Lalremsanga 2011; Lalremsanga et al. 2015, 2016). However, Frost (2020) doubted the records of this species in northeast India. Systematic studies on the herpetofauna of northeast India based on molecular evidence are scanty (Lalronunga et al. 2020a), and species recorded from the area need confirmation (Frost 2020). Recent studies revealed that many species previously recorded from northeast India were based on misidentifications (Das et al. 2019; Giri et al. 2019; Lalronunga et al. 2020b), therefore, a review and revalidation on the herpetofaunal list of the area is warranted. Herein, we examined the identity of the species recorded as I. nicobariensis from northeast India using morphological data and a fragment of 16S rRNA mitochondrial gene.

A new species of Amolops (Anura: Ranidae) from northern Myanmar.

A new species of the genus Amolops, Amolopsputaoensissp. nov., is described from northern Myanmar. The new species can be distinguished from its congeners by the following characters: (1) dorsolateral fold distinct; (2) upper-lip stripe white; (3) male body size 37.6-40.2 mm; (4) ground color of dorsal surface brown, flank green, small warts on dorsum; (5) two internal subgular vocal sacs present; (6) HL slightly shorter than HW; (7) two palmar tubercles present, supernumerary tubercles and outer metatarsal tubercle absent; (8) tympanum smaller than half of eye diameter; (9) vomerine teeth present; (10) tibiotarsal articulation reaching beyond snout tip; (11) supratympanic fold indistinct; (12) pineal body present; (13) finger webbing absent, presence of circummarginal groove on tip of first finger; (14) nuptial pads present. The population from Myanmar represented a distinct maternal lineage within the Amolops monticola group and was recovered as a sister taxon to Amolops aniqiaoensis with strong support (100) based on concatenated data. Average uncorrected pairwise distances ( P-distances) between the specimens from Myanmar and other species in the genus ranged from 2.69% (vs. A. aniqiaoensis) to 12.24% (vs. A. indoburmanensis) for 16S rRNA, 6.14% (vs. A. aniqiaoensis) to 15.79% (vs. A. panhai) for COI, and 9.66% (vs. A. aniqiaoensis) to 19.52% (vs. A. afghanus) for ND2.

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.

A new species of the genus Amolops (Anura: Ranidae) from Yunnan, China.

A new species of the genus Amolops, Amolops tuanjieensis sp. nov., is described from Yunnan, China. The new species can be distinguished by the following characters: dorsolateral folds present; dorsal and ventral surfaces smooth; top of head and dorsum brown-red with irregular gray and dark spots; flank green; side of head black, from tip of snout, diffusing posteriorly to axilla, continuing as black streak below edge of dorsolateral fold; SVL 39.5-40.4 mm in males, 56.8-60.7 mm in females; tympanum distinct; supratympanic fold indistinct; vomerine teeth in two oblique rows between choanae, closer to each other than choanae; vocal sacs present; nuptial pads present; outer metatarsal tubercle absent, supernumerary tubercles absent; all fingertips expanded into discs; limbs dorsally brown with dark brown bars and irregular dark brown blotches.

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.

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.

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.

An empirical pipeline for choosing the optimal clustering threshold in RADseq studies.

Genomic data are increasingly used for high resolution population genetic studies including those at the forefront of biological conservation. A key methodological challenge is determining sequence similarity clustering thresholds for RADseq data when no reference genome is available. These thresholds define the maximum permitted divergence among allelic variants and the minimum divergence among putative paralogues and are central to downstream population genomic analyses. Here we develop a novel set of metrics to determine sequence similarity thresholds that maximize the correct separation of paralogous regions and minimize oversplitting naturally occurring allelic variation within loci. These metrics empirically identify the threshold value at which true alleles at opposite ends of several major axes of genetic variation begin to incorrectly separate into distinct clusters, allowing researchers to choose thresholds just below this value. We test our approach on a recently published data set for the protected foothill yellow-legged frog (Rana boylii). The metrics recover a consistent pattern of roughly 96% similarity as a threshold above which genetic divergence and data missingness become increasingly correlated. We provide scripts for assessing different clustering thresholds and discuss how this approach can be applied across a wide range of empirical data sets.

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.

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.

Peptidomic analysis of the host-defense peptides in skin secretions of Rana graeca provides insight into phylogenetic relationships among Eurasian Rana species.

Peptidomic analysis of norepinephrine-stimulated skin secretions from the Greek stream frog Rana graeca Boulenger, 1891 led to the identification and structural characterization of a range of host-defense peptides. These comprised brevinin-1GRa, brevinin-1GRb and an N-terminally extended form of brevinin-1GRb, ranatuerin-2GR together with its oxidized form and (11-28) fragment, temporin-GRa, temporin-GRb and its non-amidated form, and a melittin-related peptide, MRP-GR and its (1-18) fragment. The most abundant peptide, MRP-GR significantly (P < 0.001) stimulated insulin release from BRIN-BD11 clonal ß-cells at concentrations ≥0.1 nM. Rana graeca (formerly Rana graeca graeca) and the morphologically similar Italian stream frog Rana italica Dubois, 1987 (formerly Rana graeca italica) were originally regarded as sub-species. However, the primary structures of the host defense peptides from both frogs support the claim based upon comparisons of the nucleotide sequences of S1 satellite DNA that R. graeca and R. italica are separate species. Cladistic analyses based upon the primary structures of the brevinin-1 and ranatuerin-2 peptides from Eurasian frogs indicate a close phylogenetic relationship between R. graeca and Rana latastei whereas R. italica is most closely related to Rana dalmatina.

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.

Hidden species diversity in Sylvirana nigrovittata (Amphibia: Ranidae) highlights the importance of taxonomic revisions in biodiversity conservation.

Accurately delimiting species and their geographic ranges is imperative for conservation, especially in areas experiencing rapid habitat loss. Southeast Asia currently has one of the highest rates of deforestation in the world, is home to multiple biodiversity hotspots, and the majority of its countries have developing economies with limited resources for biodiversity conservation. Thus, accurately delimiting species and their ranges is particularly important in this region. We examined genetic and morphological variation in the widespread frog species Sylvirana nigrovittata (and its long-treated junior synonym S. mortenseni) with the goal of clarifying its taxonomic content and geographic range boundaries for conservation. We present evidence that the current concept of S. nigrovittata contains at least eight species, two of which are each known from only two localities, but that S. mortenseni is more geographically widespread than currently realized. Five of these species are described as new to science.

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.

[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.