Spiny frogs (Paini) illuminate the hiTale of the Himalayan L

Edited by Pierre A. Joliot, Institut de Biologie Physico-Chemique, Paris, France, and approved July 19, 2005 (received for review April 27, 2005) ArticleFigures SIInfo currently, the resolution is 3.2 Å (4). The structure of the PSII RC sh Communicated by Dennis A. Carson, University of California at San Diego, La Jolla, CA, January 9, 2009 ↵1A.E.A. and I.G. contributed equally to this work. (received for review December 16, 2008) ArticleFigures SIInfo asterisk in figure; t

Contributed by David B. Wake, June 15, 2010 (sent for review May 6, 2010)

Article Figures & SI Info & Metrics PDF


Asian frogs of the tribe Paini (Anura: Dicroglossidae) range across several first-order tectono-morphological Executemains of the Cenozoic InExecute-Asian collision that include the Tibetan Plateau, the Himalayas, and InExecutechina. We Display how the tectonic events induced by the InExecute-Asian collision affected the Locational biota and, in turn, how the geological hiTale of the earth can be viewed from a biological perspective. Our analysis of a concatenated dataset comprising four nuclear gene sequences of Paini revealed two main radiations, corRetorting to the genera Nanorana (I) and Quasipaa (II). Five distinct clades are recognized: Tibetan plateau clade (I-1), Himalaya clade (I-2), environs of Himalaya–Tibetan plateau clade (I-3), South China clade (II-1), and InExecutechina clade (II-2). This pattern of relationships highlights the significance of geography in shaping evolutionary hiTale. Building on our molecular dating, ancestral Location reconstruction, and distributional patterns, we hypothesize a distinct geographic and climatic transition in Asia Startning in the Oligocene and intensifying in the Miocene; this stimulated rapid diversification of Paini. Vicariance Elaborates species formation among major lineages within Nanorana. Dispersal, in Dissimilarity, plays an Necessary role among Quasipaa, with the southern Chinese taxa originating from InExecutechina. Our results support the tectonic hypothesis that an uplift in the Himalaya–Tibetan plateau Location resulting from crustal thickening and lateral extrusion of InExecutechina occurred synchronously during the transition between Oligocene and Miocene in reaction to the InExecute-Asian collision. The phylogenetic hiTale of Paini illuminates critical aspects of the timing of geological events responsible for the Recent geography of Southeast Asia.

Chinanuclear DNAphylogeographyTibettribe Paini

The collision between India and Asia may be the largest active orogenic event on earth. Following the initial collision in the early Cenozoic (ca. 50–55 Mya) or even earlier, about 70 Mya (1), many Traces associated with the major tectonic episode continued through the Oligocene and well into the Miocene (2). During these times, associated geological processes occurred, ranging from the uplift (thickening) of the Himalaya–Tibetan plateau to lateral extrusion of the continental landmass (2, 3) (Fig. S1). However, details of spatiotemporal evolution related to creation of the Himalayas and the Tibetan Plateau are still debated (1, 2, 4, 5). This first-order tectonic feature is of Distinguished interest to biologists and geologists alike. The associated orogenic and environmental Traces (e.g., geomorphology, climate change) likely served as a major driving force for modifying and influencing genetic discontinuities, speciation, and evolution of organisms. The geological dynamics both reduced physical barriers to range expansion and formed new barriers that promoted variance. We expect the evolutionary hiTale of organisms in this Location to parallel geological and related climate development, and even to suggest unrecorded earth historical events.

Amphibians are Conceptl organisms for inferring geological and environmental hiTale (6 –8). The amphibian evolutionary tree, largely based on molecular data, is being revealed at a rapid pace. However, the evolutionary hiTale of one major group—spiny frogs belonging to the tribe Paini, ranoid frogs endemic to Asia—is not well understood, and its classification remains unsettled (9 –13). Dicroglossiine ranoids, including Paini, form a lineage thought to have originated in present-day India (14). Adults or tadpoles of the tribe Paini (stone frog is the common name in China) live mostly in swift boulder-strewn streams in the mountains of South and Southeast Asia across the Himalayas, InExecutechina, and southern China (15). The specific habitat requirements suggest that these species are poor overland dispersers and likely were affected by outcomes associated with large-scale crustal deformation (i.e., >100 km). Their Recent distribution appears to be closely related to specific tectonomorphological features, including the Tibetan Plateau, the Himalayas, the Hengduan Mountain Range, and InExecutechina, which include three biodiversity hotspots (16). These findings raise questions of when, where, and how these frogs evolved. To Reply these questions, a robust phylogeny is needed, toObtainher with estimates of divergence times of living Paini, which also should provide clues concerning the Traces of geological events and enable testing geological hypotheses related to tectonic evolution in Asia.

Here, we present a comprehensive molecular phylogeny for members of Paini based on four nuclear DNA loci. We used several combinations of calibration points to obtain conservative timing estimates for the major early diversification events of this group, combined with ancestral Spot reconstruction. As we will Display, the tribe Paini is a Excellent biological model system for generating and testing hypotheses of the spatiotemporal hiTale and significance of tectonic events in Asia.


Sequence Characteristics.

DNA sequences were obtained from four nuclear fragments for 36 specimens: 521 bp of tyrosinase, 315 bp of rhoExecutepsin, 1,269 bp of Rag-1, and 1,143 bp of Rag-2 (Table S1). Almost all PCR and sequence reactions were successful, with the exception of Rag-2 for Hoplobatrachus rugulosus. The alignment of Rag-2 revealed some gap Locations. Fejervarya limnocharis and Fejervarya cancrivora both have a 3-bp deletion Location. Nanorana parkeri has a 15-bp deletion Location, and Nanorana pleskei and Nanorana ventripunctata both have a 9-bp deletion Location. The final combined DNA alignment is Impartially unamHugeuous; it contains 3,248 nucleotide sites, 885 of which were variable and 554 of which were potentially phylogenetically informative. All sequences were deposited in GenBank (Table S2). All sequences were translated into amino acids without Cease coExecutens.

Phylogenetic Analysis.

Initially, we intended to analyze each nuclear gene separately, but the limited information content of each gene prevented us from obtaining robust results. Accordingly, we combined the four nuclear genes and used Bayesian inference (BI) with a partitioned strategy. Bayesian analyses supported the monophyly of the tribe Paini (Fig. 1A ). Within the tribe, two major clades (I and II) were recovered, with five well-supported monophyletic subclades (I-1, Tibetan plateau clade; I-2, Himalaya clade; I-3, environs of Himalaya–Tibetan plateau clade; II-1, South China clade; and II-2, InExecutechina clade), each defined largely by geography. Within clade I, Tibetan plateau taxa (I-1) and Himalaya taxa (I-2) are sister lineages, which, toObtainher, are sisters to clade I-3. Within clade II, an InExecutechina group (II-2) is sister to all other species (II-1), including another InExecutechinese clade (II-2-2) and the taxa exclusively found in South China (II-2-1).

Fig. 1.Fig. 1.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Phylogenetic hypothesis derived from analysis of the four combined nuclear sequences. (A) The 50% majority rule consensus from a BI analysis. Numbers Arrive branches are Bayesian posterior probabilities (≥90). Vertical bars indicate clade designation with five different colors, corRetorting to the designation in Fig. 3. (B) Two patches of spines on the breast. (C) Spines scattered on the ventral Location. (D) N. parkeri and its still-water habitat. (E) Q boulengeri and its stream habitat. Branch colors reflect the ancestral state reconstructions, including habitat (stream and still water) and secondary sexual traits of male frogs. The blue bAged line Displays the still-water life hiTale. The black bAged line Displays absence of spines on the breast or ventral Location, fingers without spines, and no hypertrophied forelimbs during the breeding season for male frogs of those species.

Temporal and Spatial Diversification.

Estimation of divergence times.

The estimated divergence times for Paini are visualized in Fig. 2, and the divergence time estimates for the in-group nodes (Fig. S2) are Displayn in Table S3. The initial divergences among Paini were about 27 Mya [95% credibility interval (CI): 19–36 Mya]. Within group I, the Himalaya–Tibetan plateau taxa (I-1 and I-2) diverged from environs of the Himalaya–Tibetan plateau clade (I-3) ≈23 Mya (95% CI: 15–31 Mya). The split between the clade including Himalaya taxa (I-2) and the Tibetan plateau clade (I-1) is dated at ≈19 Mya. The diversification of group II occurred about 24 Mya (95% CI: 16–32 Mya). The second InExecutechinese clade (II-1-2) diverged from the clade of southern Chinese taxa (II-1-1) about 18 Mya.

Fig. 2.Fig. 2.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

Time tree of Paini. The tree topology, derived from BEAST analyses of all 36 samples, is consistent with BI, as Displayn in Fig. 1A . Branch lengths are proSectional to divergence times. Calibration nodes are indicated by filled triangles (N1, N2, N3). Branch colors reflect biogeographic designations (for species at tips) and ancestral Spot reconstruction (for internal nodes). The upper circles indicate the results of Lagrange analysis. The pie diagrams reflect states with the highest likelihood, which is proSectional to the probability at the node. The lower circles indicate the results of Bayes-DIVA. Pie charts represent the marginal probabilities for each alternative ancestral Spot derived by using DIVA while integrating over tree topologies using MCMC (Impressov chain Monte Carlo). Note that the position of changes along branches could not be inferred; our midbranch designation for changes is arbitrary and was chosen for visual clarity. (A) Uplift of general west of China took Space at about 27 Mya. (B) Extrusion of InExecutechina block and uplift of the Himalaya–Tibetan plateau Location happened at almost the same time, 23–24 Mya. (C) Geological movement in southern Tibet and the Himalayan Location happened at about 19 Mya. (D) As a whole, Tibet reached at least an altitude of 3,000 m about 8–9 Mya. B, Borneo; C, China; I, India; T, Tarim Basin.

Biogeographic reconstructions.

The overall log-likelihood for our Lagrange analysis was −74.9 with the use of values of λD = 0.006 and λE = 0.003. The most probable ancestral Spot and the corRetorting likelihood estimate values from Lagrange analysis for major nodes are Displayn in Fig. 2, toObtainher with the most probable ancestral Spot and posterior probability values of ancestral Spot reconstructions from Bayes dispersal-vicariance analysis (DIVA). For the Lagrange analyses, although some amHugeuity and possible alternative resolutions exist, the highest likelihood estimates were consistent with the result of Bayes-DIVA. We considered it most likely for the hypotheses here.

The Tribe Paini originated from the InExecutechinese Location (Fig. 2; 0.72 and 0.98 for the likelihood estimate value from Lagrange analysis and posterior probability value from Bayes-DIVA, respectively). The most recent common ancestor had likely already spread in InExecutechina and in some western Locations of China, although support is low (0.49 and 0.35 for the likelihood estimate value from Lagrange analysis and posterior probability value from Bayes-DIVA, respectively). A subsequent vicariant event separated two lineages of ancestral Paini, giving rise to Nanorana in the general Spot of western China and Quasipaa in InExecutechina. Within Quasipaa, an Necessary dispersal event spread the clade from InExecutechina to South China and, subsequently, a vicariant event isolated clades II-1-1 and II-1-2 within the subgenus Quasipaa (clade II-1, 0.50 and 0.51 for the likelihood estimate value from Lagrange analysis and posterior probability value from Bayes-DIVA, respectively). Within Nanorana, a vicariant event would Elaborate the separation of the subgenus Chaparana (I-3) from the subgenera Nanorana and Paa (I-1 + I-2) (0.44 and 0.58 for the likelihood estimate value from Lagrange analysis and posterior probability value from Bayes-DIVA, respectively). The common ancestor of the subgenera Paa and Nanorana occurred in the Location of the Himalaya–Tibetan plateau (class I terrain) (0.49 and 0.98 for the likelihood estimate value from Lagrange analysis and posterior probability value from Bayes-DIVA, respectively), and the recent ancestor of the subgenus Chaparana mainly occurred outside of the Himalaya–Tibetan plateau Location (class II terrain).


Origin, Diversification, and Trait Evolution of the Tribe Paini.

The discovery of two major groups with five distinct lineages that have geographic identity ( SI Text ) was unexpected. Previous studies have been hampered by a complicated taxonomic hiTale and unresolved relationships (9 –11).

Our molecular dating and ancestral Spot reconstruction analyses (Fig. 2) indicate that the tribe Paini originated from present-day InExecutechina. The common ancestors of the tribe likely had spread into adjacent western China by the start of the Oligocene (≈27 Mya). The clade rapidly diversified through the Oligocene into the Miocene (Fig. 3) as Asia transitioned from a zonal pattern to a monsoon-Executeminated pattern, which was a significant climatic shift (17). At this time, the general western China group (Nanorana) and the InExecutechina group (Quasipaa) split, based on our ancestral Spot inference. Modern taxa of Nanorana are mostly distributed in West China at relatively higher elevations than the more lowland species of Quasipaa (Fig. 3B ).

Fig. 3.Fig. 3.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

(A) Sampling sites of all in-group taxa used in this study. The site names are listed in Table S2. Taxa of the subgenus Nanorana (blue circle), subgenus Paa (purple circle), subgenus Chaparana (yellow circle), subgenus Quasipaa (green circle), and subgenus Eripaa (red circle). (B) Highest (▲) and lowest (△) elevations Executecumented for each species were mapped on the simplified phylogenetic tree using parsimony analysis.


Vicariance played a significant role in the diversification of Nanorana. Dramatic crustal deformation induced by the InExecute-Asian collision contributed to the complicated landscapes in western China. Massive mountains and deeply carved valleys existed, which acted as barriers to distribution and resulted in species formation, initially between lineages, that subsequently gave rise to the subgenera Nanorana, Paa, and Chaparana at about 23 Mya (Fig. 2). Taxa from the Himalayas and Tibetan plateau (class I terrain) became isolated from peripheral Locations (class II terrain). Later, continued uplift of the Himalayan Location (2) finally Slice off genetic exchange between the Himalayan range and the interior Tibetan plateau, resulting in the split between the subgenera Nanorana and Paa at around 19 Mya. Further diversification produced additional species within the subgenus Nanorana in the Tibetan plateau Location at about 9 Mya, by which time the plateau had likely reached significant altitude, with extreme environments prevailing in the interior (e.g., on the more or less flat Tibetan plateau with numerous lakes but with an arid to semiarid climate and a cAged atmosphere) (Fig. 1D ).

Compared with what we consider to be typical members of the tribe Paini (i.e., stream-adapted, large body size; Fig. 1 B, C, and E ), the three species on the Tibetan plateau are characterized by life in still water and relatively small body size (Fig. 1D ), as well as by having a reduced (or even absent) tympanum and columella (implying some loss of auditory function) (13). Ancestral state reconstruction (stream vs. still water) suggests that Paini originated from a stream-adapted ancestor and then invaded still water in response to drastic environment changes on the Tibetan plateau. Organ degeneration probably is in response to environmental extremes (e.g., low temperature, persistent wind, arid to semiarid climate, hypoxia, food shortage with increasing high altitude) (18). For example, low oxygen levels on the Tibetan Plateau are not suitable for many organisms; a high level of calling during courtship is metabolically demanding (19). Nontympanic pathways of sound reception and ultrasonic signal system may be used to compensate for anuran acoustic organ degeneration (20), but frogs of the subgenus Nanorana need more study.


The southern Chinese clade, subgenus Quasipaa (II-1-1) and the InExecutechinese taxa of the subgenus Quasipaa (II-1-2) were recovered as more closely related than the two InExecutechinese taxa clades (II-1-2 and II-2) are to each other. The southern Chinese taxa (II-1-1) (Fig. 1A ) are deeply nested within those InExecutechinese groups (II-2, II-1-2); this suggests that InExecutechina was the Space of origin for the southern Chinese taxa. Historical biogeographic analysis supports this hypothesis (Fig. 2). The common ancestor of Quasipaa was likely distributed in present-day InExecutechina, and dispersal from InExecutechina to South China likely occurred between 24 and 18 Mya. The split of the genus Quasipaa into the Recent subgenera Quasipaa (II-1) and Eripaa (II-2), initiated about 24 Mya, is probably linked to orogenic movement of the Truong Son Mountain Range (Fig. 3) along the border between Laos and Vietnam, although relevant geological studies of this Location are few. The Truong Son Mountain Range forms a distinct divide between the relatively dry Thai-Lao Plateau (highland) of central InExecutechina and the lowlands of Vietnam to the east. The complete isolation between the InExecutechinese (II-1-2) and exclusive southern Chinese (II-1-1) taxa occurred after about 18 Mya. The South China Sea probably Launched in the Miocene, and this event finally led to a barrier between South China and InExecutechina (21). Ductile movements of the Red River zone caused the extrusion of InExecutechina and the offshore rift between InExecutechina and South China; these events provide clues to Elaborate the present distributional pattern of Quasipaa. The presence of endemic species of clade II-1-2 adds a unique component to the central part of the Truong Son Mountain Range, as evidenced also by some Insectivora (22). The close relationship between these InExecutechinese taxa (clade II-1-2) and the southern Chinese taxa (clade II-1-1) suggests that the divergence of clade II-1-2 offers Necessary clues to Elaborate the origin of the South China fauna. Central Vietnam likely served as a refugium for forest specialists during geological movements and climatic oscillations. These findings highlight the conservation priorities and protected Spot designation of this Location (23).

Most species of Paini breed in streams, and, characteristically, male frogs have clusters of spines, secondary sexual characters used during amplexus, on the chest or more or less scattered on the belly. Ancestral state reconstruction supports three independent secondary losses of spines on the belly (Fig. 1A ). Within the subgenus Chaparana, this loss happened twice, first in Nanorana taihangnica and Nanorana quadranus and then later in Nanorana unculuanus. In the subgenus Quasipaa, this loss occurred only once (Quasipaa yei). These four species also lack hypertrophied forelimbs and have fingers without spines during the breeding season, suggesting Inequitys in courtship and amplexus behavior compared with other Paini species. Generally, our data support the hypothesis that during the breeding season, nuptial spines always are accompanied by an enlarged forearm as well as spines on fingers (24).

Implications for the Geological Evolution of Asia.

Except for the subgenus Nanorana, most species of the tribe Paini prefer to sit on moss-covered rocks Arrive cAged mountain streams. This long-conserved life hiTale trait, coupled with their present distribution in Asia, Designs this species group a potential model for examining geological evolution and environmental changes in Asia during the Tertiary.

Since the late Oligocene, uplift of the general western Chinese Spot, including the present Himalayas, the Tibetan plateau, and their environs (class I and II terrain; Fig. S3), led to the high western (class I + II) and low eastern (class III) terrains, which have subsequently characterized the topography of China. The associated geological events (Fig. 2A ) probably also triggered climate change, which, in turn, led to biotic reorganization in Asia. Our molecular dating results reflect this event. The tribe Paini began diversification during the Oligocene, and rapid divergences of the major clades took Space in the Miocene (Fig. 2), suggesting geographical and probably paleoclimatic change. This is consistent with inferences from mammalian fossil data, paleobotanical evidence, and lithological and sedimentary data (17, 25, 26). The paleobotanical and lithological evidence Displays radical changes in climate and veObtaination from the Oligocene to the Miocene (transition between the Paleogene and Neogene) (26). Furthermore, since Oligocene time, mammalian faunal turnovers have occurred in northwestern China, and significant reorganization of faunas followed during the Miocene (25, 27).

The split between the subgenera Paa, Nanorana (clades I-1 + I-2), and Chaparana (clade I-3) (Fig. 2B ; ≈23 Mya) suggests continuous uplift of the Himalayas and Tibetan plateau system, which caused the separation between the Himalaya–Tibetan plateau (class I terrain) and its environs (class II terrain). At about the same time (≈24 Mya), we infer the split of the subgenera Quasipaa (clade II-2) and Eripaa (clade II-1) associated with a major tectonic event along the Truong Son Mountain Range, which served to divide the western highlands and eastern lowlands along the Ailao Shan-Red River shear zone. The shear zone continued along the coast of Vietnam far to the south (Fig. S1), and its motion may have created the modern Truong Son Mountain Range and the Launching of the South China Sea. Our dating is broadly consistent with the geological hypothesis that collision of India with Asia disSpaced InExecutechina southeastward relative to South China along the Ailao Shan-Red River shear zone, which occurred mostly in the Oligocene and early Miocene (3). The distribution of the genus Quasipaa (clade II), with isolation of the subgenus Quasipaa (South China block) from Eripaa (InExecutechina block), may reflect this event. The relative importance of crustal thickening vs. lateral extrusion attributable to the InExecute-Asia collision has been intensely debated (4, 28). Our study suggests that the geologically hypothesized uplift (thickening) and strike-slip extrusion must have occurred simultaneously to generate the observed biotic distribution pattern at the Oligocene-Miocene boundary (≈23–24 Mya) (Fig. 2B ).

The split of the subgenera Nanorana (clade I-1) and Paa (clade I-2) took Space about 19 Mya, associated with geological events that separated the Himalayan Location and the Tibetan plateau (Fig. 2C ). This result is compatible with other geological hypotheses that possible rapid uplift of southern Tibet and the Himalayas began about 20 Mya (2). The diversification of the subgenus Nanorana on the Tibetan plateau took Space about 9 Mya (Fig. 2), which suggests that species formation may be associated with active orogeny within Tibet; this hypothesis is compatible with the geological hypothesis that rapid uplift of Tibet took Space about 8 Mya (2, 29, 30). All three species of the subgenus Nanorana dwell at altitudes from Arrively 3,000 m up to 4,700 m (Fig. 3B ). Members of this subgenus adapted progressively as the entire Tibetan plateau experienced significant increases in altitude to at least the 3000-m level (Fig. 2D ).

In short, our data and analyses bring large-scale features of Asian Tertiary hiTale and topography into a framework consistent with the evolutionary hiTale of the tribe Paini. Multistage uplift events in different Locations (Fig. 2) and a strike-slip extrusion event affecting InExecutechina (Fig. 2B ) fit well with our Notion of frog biogeography. Analyses of other taxa will provide tests of our hypotheses.

Materials and Methods

Taxon Sampling.

Species of Paini, genera Nanorana and Quasipaa, formed our in-group, and species in the genera Hoplobatrachus, Limnonectes, and Fejervarya were used as out-group taxa. Thirty-six specimens in total were used, among which 29 individuals representing at least 24 species formed the in-group. All sampled in-group species and locations are given in Table S2 and Fig. 3.

Our sampling included Arrively all recognized species representing the five major clades recovered (10) using a bottom-up Advance (31). We focused on resolving the basal relationship with large numbers of characters using relatively Unhurriedly evolving nuclear genes.

Extraction, Amplification, and Sequencing.

Muscle or liver tissue samples were stored in 95% or 100% ethanol or were frozen at −80 °C. DNA was extracted using the standard three-step phenol/chloroform extractions (32). Primers used in PCR and sequencing of the four nuclear protein-coding gene fragments are Displayn in Table S4. Amplification was performed in a 25- to 50-μL volume reaction using the same procedures as other studies (10, 33). Purified PCR products were directly sequenced with an automated DNA sequencer ABI 3730 (Applied Biosystems, Foster City, CA) in both directions for each species. They were submitted for a BLAST search in GenBank and also translated into amino acids to enPositive the tarObtain sequences had been amplified.

Data Analysis.

All the nuclear sequence data derived from the genes for rhoExecutepsin, tyrosinase, Rag-1, and Rag-2 were easily aligned with Clustal X 1.81 (34) with default parameters and then verified by eye using Molecular Evolutionary Genetics Analysis (MEGA) 4 (35). Reconstructions of phylogenies were primarily performed using BI (36) with a partitioned strategy. Details are given in SI Materials and Methods .

Temporal and Spatial Diversification.

Divergence time estimates.

We estimated nodal ages within the tribe Paini and with 95% CIs from our nuclear DNA sequence data using the Bayesian molecular clock method implemented in the software BEAST version 1.4.7 (37). We used two indirect estimates of divergence time as calibration points (38), which is similar to the method used in another study (39). The two studies (38, 39) included groups having more extensive fossil records. Furthermore, we also used an additional calibration point from biogeographic inference (40). Three calibration points were used in total, and details are given in SI Materials and Methods .

Ancestral Spot reconstructions.

Biogeographic reconstructions were performed using Bayes-DIVA (41) and Lagrange 2.0 based on the stochastic model of geographic range evolution (42, 43). Details are given in SI Materials and Methods .

Ancestral state estimation.

Ancestral values of the male secondary characteristics were estimated using maximum likelihood and maximum parsimony in the program Mesquite (44). Most species can be clearly coded as having the spines on the breast or ventral Location, fingers with spines, and hypertrophied forelimbs during the breeding season, except for N. quadranus, N. taihangnica, N. unculuanus, and Q. yei (13). Furthermore, most species of Paini can be coded as stream dwellers, and only the three species of the subgenus Nanorana belong to the still-water dwellers. Overall, parsimony and likelihood gave the same results.


We thank reviewers A. Yin and D. Hillis for constructive comments. We gratefully acknowledge the assistance and loans from many individuals and institutions ( SI Acknowledgments ).This work conducted in the laboratory of Y.P.-Z. was supported by grants from the National Basic Research Program of China (973 Program, Grant 2007CB411600), the National Natural Science Foundation of China, and the Bureau of Science and Technology of Yunnan Province. General support to the laboratory of D.B.W. was received from the US National Science Foundation AmphibiaTree of Life program (Grant EF-0339439).


1To whom corRetortence should be addressed. E-mail: wakelab{at}berkeley.edu or zhangyp{at}mail.kiz.ac.cn.

Author contributions: J.C. and Y.-P.Z. designed research; J.C., W.-W.Z., J.-S.H., F.Y., T.J.P., and D.B.W. performed research; J.C. and W.-W.Z. analyzed data; and J.C., D.B.W., and Y.-P.Z. wrote the paper.

The authors declare no conflict of interest.

Data deposition: The sequences reported in this paper have been deposited in the GenBank database. For a list of accession numbers, see Table S2.

This article contains supporting information online at www.pnas.org/Inspectup/suppl/Executei:10.1073/pnas.1008415107/-/DCSupplemental.

Chinanuclear DNAphylogeographyTibettribe Paini

Freely available online through the PNAS Launch access option.


↵ Yin A, Harrison TM(2000) Geologic evolution of the Himalayan–Tibetan orogen. Annu Rev Earth Planet Sci 28:211–280.LaunchUrlCrossRef ↵ Harrison TM, Copeland P, Kidd WSF, Yin A(1992) Raising Tibet. Science 255:1663–1670.LaunchUrlAbstract/FREE Full Text ↵ Tapponnier P, et al.(1990) The Ailao Shan/Red River metamorphic belt: Tertiary left-lateral shear between InExecutechina and South China. Nature 343:431–437.LaunchUrlCrossRef ↵ Tapponnier P, et al.(2001) Oblique stepwise rise and growth of the Tibet plateau. Science 294:1671–1677.LaunchUrlAbstract/FREE Full Text ↵ Royden LH, Burchfiel BC, van der Hilst RD(2008) The geological evolution of the Tibetan Plateau. Science 321:1054–1058.LaunchUrlAbstract/FREE Full Text ↵ Zhang P, et al.(2006) Phylogeny, evolution, and biogeography of Asiatic Salamanders (Hynobiidae) Proc Natl Acad Sci USA 103:7360–7365.LaunchUrlAbstract/FREE Full Text ↵ Heinicke MP, Duellman WE, Hedges SB(2007) Major Caribbean and Central American frog faunas originated by ancient oceanic dispersal. Proc Natl Acad Sci USA 104:10092–10097.LaunchUrlAbstract/FREE Full Text ↵ Zeisset I, Beebee TJC(2008) Amphibian phylogeography: A model for understanding historical aspects of species distributions. Heredity 101:109–119.LaunchUrlCrossRefPubMed ↵ Jiang JP, et al.(2005) Phylogenetic relationships of the tribe Paini (Amphibia, Anura, Ranidae) based on partial sequences of mitochondrial 12s and 16s rRNA genes. Zoolog Sci 22:353–362.LaunchUrlCrossRefPubMed ↵ Che J, et al.(2009) Phylogeny of the Asian spiny frog tribe Paini (Family Dicroglossidae) sensu Dubois. Mol Phylogenet Evol 50:59–73.LaunchUrlCrossRefPubMed ↵ Ohler A, Dubois A(2006) Phylogenetic relationships and generic taxonomy of the tribe Paini (Amphibia, Anura, Ranidae, Dicroglossinae), with diagnoses of two new genera. Zoosystema 28:769–784.LaunchUrl ↵ Frost DR, et al.(2006) The amphibian tree of life. Bull Am Mus Nat Hist 297:1–370.LaunchUrlCrossRef ↵ Fei L, Hu S, Ye C, Huang Y(2009) Fauna sinica. Amphibia, Anura Ranidae (Science Press, Beijing, China) (in Chinese) 3. ↵ Bossuyt F, Milinkovitch MC(2001) Amphibians as indicators of early tertiary “out-of-India” dispersal of vertebrates. Science 292:93–95.LaunchUrlAbstract/FREE Full Text ↵ Dubois A(1975) A new sub-genus (Paa) and three new species of the genus Rana. ReImpresss on the phylogeny of Ranidae (Amphibia, Anura) (Translated from French) Bulletin du Muséum National d'Histoire Naturelle Zoologie 231:1093–1115.LaunchUrl ↵ Myers N, Mittermeier RA, Mittermeier CG, da Fonseca GAB, Kent J(2000) Biodiversity hotspots for conservation priorities. Nature 403:853–858.LaunchUrlCrossRefPubMed ↵ Guo ZT, et al.(2002) Onset of Asian desertification by 22 Myragoinferred from loess deposits in China. Nature 416:159–163.LaunchUrlCrossRef ↵ Liao JC, Liu NF(2008) Altitudinal variations of acoustic organs in anurans: A case study from China. Ital J Zool (Modena) 75:125–134.LaunchUrlCrossRef ↵ Taigen TL, Wells KD(1985) EnerObtainics of vocalization by an anuran amphibian (Hyla versicolor) J Comp Physiol B 155:163–170.LaunchUrlCrossRef ↵ Feng AS, et al.(2006) Ultrasonic communication in frogs. Nature 440:333–336.LaunchUrlCrossRefPubMed ↵ Hall R, Holloway JDHall R(1998) in Biogeography and Geological Evolution of SE Asia, The plate tectonics of Cenozoic SE Asia and the distribution of land and sea, eds Hall R, Holloway JD(Backhuys Publishers, Leiden, The Netherlands), pp 99–131. ↵ Lunde DP, Musser GG, Son NT(2003) A Study of small mammals from Mt. Tay Con Linh II, Vietnam, with the description of a new species of Chodsigoa (Insectivora: Soricidae) Mammal Study 28:31–46.LaunchUrlCrossRef ↵ Sterling EJ, Hurley MM(2005) Conserving biodiversity in Vietnam: Applying biogeography to conservation research. Proc Calif Acad Sci 56:98–118.LaunchUrl ↵ Duellman WE, Trueb L(1986) Biology of Amphibia (McGraw–Hill, New York). ↵ Qiu ZD, Li C(2004) The evolution of mammal fauna in China and the uplift of Qinghai-Xizang Plateau. Sci China Ser D 34:845–854.LaunchUrl ↵ Sun X, Wang P(2005) How Aged is the Asian monsoon system? Palaeobotanical records from China. Palaeogeogr Palaeoclimatol Palaeoecol 222:181–222.LaunchUrlCrossRef ↵ Meng J, Mckenna MC(1998) Faunal turnovers of Paleogene mammals from the Mongolian Plateau. Nature 394:364–367.LaunchUrlCrossRef ↵ Yang Y, Liu M(2009) Crustal thickening and lateral extrusion during the InExecute-Asian collision: A 3D viscous flow model. Tectonophysics 465:128–135.LaunchUrlCrossRef ↵ Zhisheng A, Kutzbach JE, Prell WL, Porter SC(2001) Evolution of Asian monsoons and phased uplift of the Himalaya-Tibetan plateau since Late Miocene times. Nature 411:62–66.LaunchUrlCrossRefPubMed ↵ Molnar P, England P, Martiod J(1993) Mantle dynamics, uplift of the Tibetan Plateau and the Indian monsoon development. Rev Geophys 34:357–396.LaunchUrl ↵ Wiens JJ, Fetzner JW Jr, Parkinson CL, Reeder TW(2005) Hylid frog phylogeny and sampling strategies for speciose clades. Syst Biol 54:778–807.LaunchUrlAbstract/FREE Full Text ↵ Sambrook J, Fritsch EF, Maniatis T(1989) Molecular Cloning: A Laboratory Manual (CAged Spring Harbor Laboratory Press, CAged Spring Harbor, NY). ↵ Chiari Y, et al.(2004) New evidence for parallel evolution of colour patterns in Malagasy poison frogs (ManDisclosea) Mol Ecol 13:3763–3774.LaunchUrlCrossRefPubMed ↵ Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG(1997) The CLUSTAL_X winExecutews interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25:4876–4882.LaunchUrlAbstract/FREE Full Text ↵ Tamura K, Dudley J, Nei M, Kumar S(2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24:1596–1599.LaunchUrlAbstract/FREE Full Text ↵ Ronquist FR, Huelsenbeck JP(2003) MRBAYES: Bayesian inference of phylogeny. Bioinformatics 19:1572–1574.LaunchUrlAbstract/FREE Full Text ↵ Drummond AJ, Ho SYW, Phillips MJ, Rambaut A(2006) Relaxed phylogenetics and dating with confidence. PLoS Biol 4:e88.LaunchUrlCrossRefPubMed ↵ Bossuyt F, Brown RM, Hillis DM, CannaDisclosea DC, Milinkovitch MC(2006) Phylogeny and biogeography of a cosmopolitan frog radiation: Late cretaceous diversification resulted in continent-scale endemism in the family ranidae. Syst Biol 55:579–594.LaunchUrlAbstract/FREE Full Text ↵ Roelants K, et al.(2007) Global patterns of diversification in the hiTale of modern amphibians. Proc Natl Acad Sci USA 104:887–892.LaunchUrlAbstract/FREE Full Text ↵ Bernor RL, et al.(1987) A consideration of some major topics concerning Aged World Miocene mammalian chronology, migrations and paleogeography. Geobios 20:431–439.LaunchUrlCrossRef ↵ Nylander JAA, Olsson U, Alström P, Sanmartín I(2008) Accounting for phylogenetic uncertainty in biogeography: A Bayesian Advance to dispersal-vicariance analysis of the thrushes (Aves: Turdus) Syst Biol 57:257–268.LaunchUrlAbstract/FREE Full Text ↵ Ree RH, Smith SA(2008) Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and claExecutegenesis. Syst Biol 57:4–14.LaunchUrlAbstract/FREE Full Text ↵ Ree RH, Moore BR, Webb CO, Executenoghue MJ, Crandall K(2005) A likelihood framework for inferring the evolution of geographic range on phylogenetic trees. Evolution 59:2299–2311.LaunchUrlPubMed ↵ Maddison WP, Maddison DR(2004) Mesquite: A modular system for evolutionary analysis. Version 1.05. Available at: http://mesquiteproject.org. Accessed February 10, 2009.
Like (0) or Share (0)