Translator Disclaimer
1 November 2000 Phylogeny, Taxonomy, and Biogeography of the Oriental Pitvipers of the Genus Trimeresurus (Reptilia: Viperidae: Crotalinae): A Molecular Perspective
Author Affiliations +
Abstract

Based on sequence variation in 806 bp of the mitochondrial 12S rRNA gene, phylogenetic relationships were inferred for 14 species of Trimeresurus (sensu lato) including all East Asian members. Samples analyzed also included representatives of all assemblages of species that are frequently treated as separate genera except for T. mangshanensis, a type species of the recently described monotypic genus Ermia. Results support some previous accounts chiefly from morphological studies, such as distinct divergence of T. wagleri from the remainder, and monophyly of T. mucrosquamatus, T. flavoviridis, T. jerdonii, T. elegans and T. tokarensis. On the other hand, our results negated a putative close affinity of T. monticola and T. okinavensis, and indicated the sister relationship of the latter with T. gracilis. Phylogenetic relationships revealed in this study suggested that the genus Trimeresurus dispersed into the Ryukyu region at least three times, and that T. flavoviridis and T. tokarensis from the central Ryukyus constitute a relict clade.

INTRODUCTION

Trimeresurus Lacepede, as defined by Brattstrom (1964) (henceforth referred to as Trimeresurus [sensu lato]), is the largest group of venomous snakes in tropical-subtropical Asia, consisting of some 40 species (David and Tong, 1997; McDiarmid et al., 1999). A few species of the genus are especially dangerous and thus are of great concern to public-health authorities and medical workers in some areas (Gopalakrishnakone and Chou, 1990). Distribution of Trimeresurus (sensu lato), covering not only a broad part of the Asian continent but also a number of adjacent archipelagos with complex geohistories, makes the genus a particularly suitable model to investigate biological consequences of past geological events as expressed in a phylogeographical pattern of extant species.

A number of authors have addressed phylogenetic relationships within Trimeresurus (sensu lato) based on morphological characters, some hypothesizing scenarios for the divergence process of the genus (Brattstrom, 1964; Zhang, 1995), others proposing its taxonomic subdivisions (Burger, 1971; Hoge and Romano-Hoge, 1981, 1983; Zhang, 1993). Nevertheless, very little consensus has yet been attained to date for almost every relevant issue (e.g., see comments by Groombridge [1986], Toriba [1989] and McDiarmid et al. [1999] on the taxonomic changes proposed by Burger [1971], Hoge and Romano-Hoge [1983] and Zhang [1993], respectively).

Recently a few preliminary molecular studies, while confirming great divergences within Trimeresurus (sensu lato) predicted by previous morphological studies (see above) (Kraus et al., 1996; Knight et al., 1992), elucidated prominent convergences in some of the traditionally used taxonomic characters of this and other crotalinae genera (Cullings et al., 1997; Malhotra and Thorpe, 1997; Vidal et al., 1997; Parkinson, 1999). These findings suggest the presence of severe constraints in the use of morphological characters for the reconstruction of phylogenetic relationships in this and related genera.

For the purpose of testing some previous hypotheses regarding the phylogeny and biogeography of Trimeresurus (sensu lato), we sequenced a fragment of the mitochondrial 12S ribosomal RNA gene (12S rRNA) for representatives of the genus. Sequence variation in this domain is usually used to resolve higher level phylogenetic relationships in reptiles (Heise et al., 1995; Honda et al., 1999; Ota et al., 1999), but has also been shown to be effective to infer the infrageneric phylogeny of viperid snakes (e.g., Knight et al., 1993). Analyses of sequence data suggest relationships which substantially differ from those hypothesized for Trimeresurus (sensu lato) based on previous morphological studies. Also, our results yield new insights into the historical biogeographical analysis of terrestrial organisms in the subtropical East Asian islands.

MATERIALS AND METHODS

Eighteen specimens of 14 species and one subspecies of Trimeresurus (sensu lato) were used in this study (Table 1: Appendix 1). Of these taxa, four species and one subspecies (T. monticola monticola, T. monticola makazayazaya, T. jerdonii, T. gracilis, and T. puniceus) were examined by the molecular technique for the first time. The taxa examined include all eight species of the genus occurring in Taiwan and the Ryukyu Archipelago, as well as representatives of all species assemblages of Trimeresurus (sensu lato) but one (Ermia Zhang) that are frequently treated as separate genera (i.e., Ovophis Burger, Protobothrops Hoge et Romano-Hoge, Tropidolaemus Wagler, and Trimeresurus sensu stricto)(Hoge and Romano-Hoge, 1983; Toriba, 1993; Zhang, 1993; Wüster et al., 1997; McDiarmid et al., 1999). Specimens of Deinagkistrodon acutus and Vipera russelii that, respectively, represent another crotaline genus and a different viperid subfamily, Viperinae, were also included in the analyses as outgroups. Each specimen was designated as constituting an independent operational taxonomic unit (OTU) by itself for the phylogenetic analyses.

Table 1

Species and subspecies of Trimeresurus (sensu lato), and two outgroup species examined in this study. Information regarding their distributions was taken from Toriba (1993) and McDiarmid et al. (1999). See Appendix 1 for locality data for materials used in our analyses.

i0289-0003-17-8-1147-t01.gif

Crude DNA was extracted from fresh or 95% ethanol-preserved muscle samples following Kocher et al. (1989), but with incubation extended for approximately ten hours. After proteinase K digestion, DNA was purified by two times of phenol extraction, once or twice of phenol/chloroform/isoamyl alcohol (25:24:1) extraction, and once of chloroform/isoamyl alcohol (24:1) extraction, each followed by precipitation in ethanol with 1/10 volume NaOAC.

A portion of the 12S rRNA gene of the mitochondrial genome was amplified by the polymerase chain reaction (PCR) using two primers—SN1: 5′-AGTCTGCTCAAAAAGATTAATGTTAA-3′; and SN2: 5′-TCTTGGTCTGAAACCTCAGTTACCTA-3′ (Wang and Tu, 1997). PCR reactions were performed in 50 ul volumes consisting of 10 mM Tris-HCl (pH 9.0), 50mM KCl, 15mM MgCl2, 0.1% [w/v] gelatin, 1% Triton X-100, 0.4 pM primer, 0.2 mM dNTP, 50–100 ng of crude DNA, and 1 U Taq polymerase (InViTAQ, Germany). The temperature regimen of 35 cycles, subsequent to the two minutes of initial denaturation at 94°C, was 1 min at 94°C, 1 min at 50°C, and 1.5 min at 72°C.

PCR products were purified with the Gene Clean III elution kit (BIO 101, CA). Both DNA strands were sequenced using dye terminator cycle-sequencing reactions that were subsequently loaded on an Applied Biosystems 377A automatic sequencer. The numbering system followed Anderson et al. (1981).

DNA sequences were aligned by using the default parameters of CLUSTAL W (Thompson et al., 1994). Adjustments were made visually on the basis of maximum nucleotide similarity. Gaps and insertions (indels) were excluded from the subsequent analyses.

Based on the aligned sequences, phylogenetic relationships among OTUs were inferred using neighbor-joining (NJ: Saitou and Nei, 1987) and maximum parsimony (MP) methods. For the former, three distance models (i.e., Jukes and Cantor's [1969] one parameter model, Kimura's [1980] two parameter model, and Tamura and Nei's [1993] model) were used to correct pairwise distances for multiple hits. Resultant distance matrices were subjected to NJ analyses using PHYLIP 3.5c (Felsenstein, 1993). Degrees of supports for internal branches of the resultant tree were assessed by 1,000 boot-strap replications (Felsenstein, 1985) for each weighting scheme used. By using PAUP version 4.0 (Swofford, 1998), we further examined the content of phylogenetic information in our data set by checking the skewness of the tree distribution and the g1 values (Hillis and Huelsenbeck, 1992) for 10,000 random trees in each weighting scheme.

MP analysis was performed using PAUP*, in which each nucleotide base was regarded as an unordered character, and the four kinds of salts as different character states. A total of 1,000 bootstrap pseudoreplications (Felsenstein, 1985) were conducted using the heuristic algorithm of PAUP*. The result of this analysis was expressed as a 50% majority-rule consensus bootstrap tree. Skewness of the tree distribution and the g1 values for 10,000 random trees were also calculated.

RESULTS

Aligned sequence of 12S rRNA gene consisted of 806 bp including indels (Appendix 2). The first 230 and 200 sites of T. wagleri and D. acutus, respectively, could not be sequenced despite our several attempts. These portions were treated as gaps, and were excluded from the analyses.

Comparisons of aligned sequence revealed 362 variable sites, of which 200 were phylogenetically informative. The base composition was slightly biased, with the average nucleotide frequencies of A, T, C and G being 37.6%, 19.8%, 25.1% and 17.5%, respectively. The estimated g1 values for trees resulting from our analyses were lower than critical values (p = 0.01: Hillis and Huelsenbeck, 1992), indicating that the region sequenced contained phylogenetic information.

Distance values resulting from the three different models were very similar to each other. Likewise, resultant NJ trees showed no distinct differences in branching topology or relative branch length among the three models, either. Therefore, we provide results of NJ analysis using Kimura's distance matrix only (Fig. 1) (distance matrices and results of NJ using the remaining models are available from H.-Y. Wang upon request).

Fig. 1

Neighbor-joining (NJ) tree derived from Kimura's (1980) two parameter distance matrix for 12S rRNA sequence data. Numerals above nodes indicate proportion of support in 1000 bootstrap replications (BP). Branches without BP values were not supported in ≥ 50% of the replicates. Numerals below nodes are reference numbers used in the text.

i0289-0003-17-8-1147-f01.gif

In the NJ analysis, monophyly of the genus Trimeresurus sensu lato was supported only by a low bootstrap proportion (BP:<50%), whereas the monophyly of conspecific and consubspecific OTUs invariably received a complete bootstrap support. Five major nodes (nodes 1–5) were recognized in NJ tree, of which two consisted of single species, T. monticola (node 1) and T. wagleri (node 5), whereas another one with a complete bootstrap support (node 3) consisted of T. gracilis and T. okinavensis. Of the remaining two major nodes, one with a rather poor BP support (node 4) consisted of exclusively arboreal species, with T. popeiorum possibly initially diverging from the remainder, followed by T. trigonocephalus and T. puniceus in order. The remaining major node with a complete bootstrap support (node 2) was further divided into two nodes, one consisting of T. mucrosquamatus and T. elegans (node 6), and the other of T. jerdonii, T. flavoviridis and T. tokarensis (node 7). The last two species received a 100% bootstrap support for their sister-group relationship, whereas the sister-group relationship of these species and T. jerdonii was supported with a lower, but still substantial bootstrap proportion (67%).

The 50% majority-rule consensus tree resulting from the MP analysis did not contradict with NJ tree in branching topology (Fig. 2), but differs from the latter in the degree of supports for a few nodes. For example, monophyly of the genus Trimeresurus sensu lato was supported with a much higher BP value (85%), whereas the value in support for the node 7 (54%) was lower than the corresponding value in NJ tree (see above).

Fig. 2

Majority-rule consensus cladogram derived from maximum parsimony analysis (MP) of 12S rRNA sequence data. Numerals above nodes indicate proportion of support in 1000 bootstrap pseudoreplications. Numerals below nodes are reference numbers used in the text. Consistency intex excluding uninformative sites (Kluge and Farris, 1969) = 0.457; retention index (Farris, 1989) = 0.530.

i0289-0003-17-8-1147-f02.gif

DISCUSSION

Phylogeny

Our results corroborate the monophyly of T. mucrosquamatus, T. flavoviridis, T. jerdonii, T. elegans, and T. tokarensis, as was explicitly or implicitly predicted by a few previous authors (Hoge and Romano-Hoge, 1981; Kraus et al., 1996; David and Ineich, 1999; Parkinson, 1999). Moreover, initial divergence of T. wagleri from the remaining members of Trimeresurus (sensu lato) in our analyses (Figs. 1 and 2) agrees well with the relationships inferred by Brattstrom (1964) on the basis of morphological analysis. However, such relationships of T. wagleri and other species contradict with the hypothesis resulting from molecular analyses by Parkinson (1999), in which T. wagleri is considered to be closer to T. albolabris and T. stejnegeri than are some other species common to our study (i.e., T. elegans, T. flavoviridis, T. tokarensis, and T. okinavensis). Moreover, T. wagleri is shown to be closest to D. acutus in Parkinson's (1999) trees. Because either his trees or ours failed to receive sufficently high bootstrap supports (> 90%: Shaffer et al., 1997) for these relationships, it is obvious that unambiguous solutions of these inconsistencies on the molecular ground require additional sequence data. Nevertheless, considering their congruence with results of morphological analysis (see above), the relationships depicted in our trees seem to be more likely at present.

The remaining portions of our NJ and MP trees show considerable discrepancies with previously hypothesized relationships. For example, our results do not support at all the close affinity of Trimeresurus monticola and T. okinavensis, that were combined under a separate generic name, Ovophis (type species: monticola), by Burger (1971) and Hoge and Romano-Hoge (1981) together with three other species not studied here (i.e., T. chaseni, T. convictus, and T. tonkinensis). Such a distant relationship of monticola and okinavensis was also implied by results of comparative studies of head musculature by Groombridge (1986). Both monticola and okinavensis are stout-bodied ground dwellers, whereas most other species examined here, including the possibly most primitive T. wagleri (see above), are arboreal or semi-arboreal and have thinner bodies (Koba, 1962; Zhao et al., 1998; Gopalakrishnakone and Chou, 1990). It is thus possible that external morphological similarities between monticola and okinavensis actually represent convergence resulting from independent adaptations to similar, non-arboreal life-styles.

Brattstrom (1964), on the basis of descriptions by Maslin (1942), surmised that T. gracilis is phylogenetically closest to T. puniceus (and its sibling species, T. borneensis: see McDiarmid [1999]). Results of our analysis, however, negate such a view, and strongly suggest a much closer affinity of T. gracilis to T. okinavensis. It is likely that similarities in some scale characters between T. gracilis and T. puniceus emphasized in a key by Maslin (1942), which obviously let Brattstrom (1964) assume their close affinity, actually represent symplesiomorphy or convergence.

Based on the allozyme analyses, Toda et al. (1999) surmised a relatively close relationship between T. monticola and T. okinavensis. However, materials subjected to their analyses were limited to T. elegans, T. flavoviridis, T. mucrosquamatus, and T. tokarensis, besides those two species. Moreover, subsequent allozyme reanalyses by incorporating data for a single T. gracilis yielded results that predict the closest affinity between this species and okinavensis (Toda and Ota, unpublished data). Thus, it seems unlikely that the relationships illustrated by more comprehensive allozyme analyses substantially contradict with the relationships indicated by the present analyses.

Taxonomic implications

The relationships depicted by our analyses (Figs. 1 and 2) support the validity of Tropidolaemus, which had usually been regarded as synonymous with Trimeresurus (e.g., see a synonym list of wagleri in David and Vogel [1996]), and was resurrected to accommodate wagleri by Hoge and Romano-Hoge (1981) on the basis of Burger's (1971) morphological data. For further confirmation for the current taxonomic arrangement of Tropidolaemus, additional analyses are desired by incorporating sequence data for huttoni from India, a putative second species of the genus (David and Vogel, 1998).

Our results also support the validity of Protobothrops (node 2), which was originally described for T. flavoviridis (type species), T. mucrosquamatus and T. jerdonii, and subsequently also thought to include T. elegans and T. tokarensis, as well as T. kaulbacki, T. strigatus and T. xiangchengensis (Hoge and Romano-Hoge, 1983; Kraus et al., 1996; David and Ineich, 1999; Parkinson, 1999). In contrast, relationships revealed here negate the validity of Ovophis as arranged by Hoge and Romano-Hoge (1981). Use of the generic name, Ovophis, thus should be avoided until more plausible delimitation and morphological redefinitions are made on the basis of comprehensive phylogenetic analyses for Ovophis sensu Hoge and Romano-Hoge (1981). Considering the possible close relationships of monticola and the node 2 species as suggested in our analyses (Figs. 1 and 2), use of the generic name, Protobothrops, though most likely being monophyletic by itself (see above), should be also avoided until the relationships of other species assigned to Ovophis by Hoge and Romano-Hoge (1981) are clarified.

It is obvious that the remaining portion of Trimeresurus (sensu lato) is yet highly heterogeneous, because it includes at least two distinct lineages, one represented by two strictly terrestrial species (T. okinavensis and T. gracilis) and the other by several exclusively arboreal species (T. popeiorum, T. albolabris, T. stejnegeri, T. puniceus and T. trigonocephalus). More comprehensive analyses incorporating data for the remaining species of Trimeresurus (sensu lato), including T. gramineus (type species of the genus: McDiarmid et al., 1999) and T. mangshanensis (type species of Ermia: Zhang, 1993), are strongly desired to elucidate detailed process of divergence of those terrestrial and arboreal species, and to revise their classification at the generic level.

There is a remarkable sequence divergence between T. monticola monticola from southern continantal China and the two samples of T. m. makazayazaya, one from eastern continental China and the other from Taiwan (Fig. 1). Because values of Kimura's (1980) distance between these supposedly conspecific subspecies (0.055–0.058) are as great as or greater than those between some different species (e.g., 0.043 between T. mucrosquamatus and T. elegans, 0.035 between T. tokarensis and T. flavoviridis, and 0.055 between T. gracilis and T. okinavensis), it is probable that monticola and makazayazaya actually represent two different species. This suggests the necessity for detailed analyses of geographic variation in the highly polytypic T. monticola (McDiarmid et al., 1999; Zhao et al., 1998).

Biogeography

It is assumed that in some vertebrate examples, taxonomic diversity in Taiwan has increased through a series of in situ divergences facilitated by the large size and diverse geomorphology of this island (e.g., Yu, 1995; Ota, 1997). With respect to Trimeresurus (sensu lato), however, the diversity within Taiwan seems to be attributable to the multiple colonizations rather than to the in situ divergences, because the four Taiwanese representatives, while being distant from each other genetically, have close relatives outside Taiwan. Toda et al. (1999) analyzed allozyme variation in Trimeresurus (sensu lato), exclusive of T. gracilis and T. stejnegeri from the East Asian islands. Based on the resultant phylogeny, (monticola, okinavensis)((mucrosquamatus, elegans)(flavoviridis, tokarensis)), they assumed a landbridge dispersal of the genus into the Ryukyu Archipelago and attributed the divergence between the mucrosquamatus-elegans clade and the flavoviridis-tokarensis clade to a vicariance event involved by the initial insularization of the central Ryukyus in the late Pliocene. They also referred to the relatively large genetic distance between those clades (Nei's [1978] D = 0.198–0.289) as evidence for long isolation of the central Ryukyus from the southern Ryukyus and Taiwan.

However, present analyses with the inclusion of T. gracilis and T. jerdonii besides others that were not included in Toda et al.'s (1999) analyses yielded an essentially different picture for the historical biogeography of the East Asian Trimeresurus (sensu lato). First of all, present results suggest that T. okinavensis of the central Ryukyus was derived from a lineage distinct from T. monticola. The closest affinity of T. okinavensis with the geographically disjunct T. gracilis, a species endemic to the high altitude area of Taiwan (Ota, 1991), suggests an extremely relictual nature of the okinavensis-gracilis clade (Darlington, 1957).

More surprising is the possible sister relationship of the flavoviridis-tokarensis clade with T. jerdonii from the inland area of continental China, not with the mucrosquamatus-elegans clade from the immediate south (Fig. 3) as was postulated by Toda et al. (1999). Such a relationship suggests that T. flavoviridis and T. tokarensis were derived from an ancestral form that dispersed to the Ryukyus independently from the dispersal of the T. mucrosquamatus-elegans clade to this archipelago. This further suggests that the sister lineage of the T. flavoviridis-tokarensis clade, once occurring in the southern Ryukyus, Taiwan and eastern continent, has subsequently disappeared from these broad areas.

Fig. 3

Map of East Asia, showing geographic distributions of T. elegans (a: the Yaeyama Group of the southern Ryukyus), T. mucrosquamatus (remaining portion of stippled area), T. flavoviridis (b: the Okinawa and Amami Groups of the central Ryukyus), T. tokarensis (c: the southern Tokara islands of the central Ryukyus), and T. jerdonii (remaining portion of hatched area). Data for ranges were taken from Zhao and Adler (1993) and Toda et al. (1999).

i0289-0003-17-8-1147-f03.gif

There remain four species of Trimeresurus in continental China that have not yet been examined phylogenetically (T. mangshanensis, T. medoensis, T. tibetanus, and T. xiangchengensis: Zhao and Adler, 1993; Zhao et al., 1998). Thus, one may argue that one or more of them may be even closer to the flavoviridis-tokarensis clade than T. jerdonii is. However, even if this is actually the case, it will not lead to any substantial changes in the above-mentioned view on the historical biogeography of T. flavoviridis and T. tokarensis, because geographic ranges of those four species were also confined to far inland of the continent like that of T. jerdonii (Zhao and Adler, 1993; Zhao et al., 1998). It is thus likely that the flavoviridis-tokarensis clade of the central Ryukyus is in a relict state like several other reptiles occurring in this region (Hikida and Ota, 1997; Ota, 1998).

As such, it is likely that the current Trimeresurus assemblage of the Ryukyus was derived from three independent dispersals from Taiwan and the continent–one by the common ancestor of T. gracilis and T. okinavensis, another by the common ancestor of T. jerdonii and the T. flavoviridis-tokarensis clade, and the other by the common ancestor of T. mucrosquamatus and T. elegans. Present results also emphasize the importance of examining relevant inland representatives of the continent, if any, for appropriate biogeographical assessment of possibly relict organisms in the central Ryukyus.

Acknowledgments

We are grateful to Y. H. Chen, S. L. Chen, H. C. Lin, K. Y. Lue, J. J. Mao, G. Shang, G. Vogel and the staff of Field Museum of Natural History for supplying tissues for our study, and to M. Honda and H. Sato for technical advice. We also thank S.-H. Li, H. M. Smith and two anonymous reviewers for the critical reading of early versions of the manuscript, and P. David and I. Ineich for provision of relevant literature. Financial support was provided by National Science Council, R. O. C. (NSC 872311B003003B17), and Japan Ministry of Education, Science, Sports and Culture (Grant-in-aid C-11833013).

REFERENCES

1.

S. Anderson, A. T. Bankier, B. G. Barrell, M. H. L. de Bruijin, A. R. Coulson, J. Droun, I. C. Eperon, D. P. Nierlich, B. A. Roe, F. Sanger, P. H. Schreier, A. J. H. Smith, R. Staden, and I. G. Young . 1981. Sequence and organization of the human mitochondrial genome. Nature 290:457–465. Google Scholar

2.

B. H. Brattstrom 1964. Evolution of the pitvipers. Trans San Diego Soc Nat Hist 13:185–268. Google Scholar

3.

W. L. Burger 1971. Genera of Pitvipers (Serpentes: Crotalidae). Ph. D. Dissertation. University of Kansas. Lawrence. Google Scholar

4.

K. W. Cullings, D. J. Morafka, J. Hernandez, and J. W. Roberts . 1997. Reassessment of phylogenetic relationships among pitviper genera based on mitochondrial cytochrome b gene sequences. Copeia 1997:429–432. Google Scholar

5.

P. J. Darlington Jr 1957. Zoogeography: The Geographic Distribution of Animals. John Wiley. New York. Google Scholar

6.

P. David and H. Tong . 1997. Translations of recent descriptions of Chinese pitvipers of the Trimeresurus-complex (Serpentes, Viperidae), with a key to the complex in China and adjacent areas. Smithsonian Herpetol Inf Serv 112:1–31. Google Scholar

7.

P. David and I. Ineich . 1999. Les serpents venimeux du monde: systematique et repartition. Dumerilia 3:3–499. Google Scholar

8.

P. David and G. Vogel . 1996. The Snakes of Sumatra. Edition Chimera. Frankfurt am Main. Google Scholar

9.

P. David and G. Vogel . 1998. Redescription of Trimeresurus huttoni Smith, 1949 (Serpentes, Crotalinae), with a discussion of its relationships. Hamadryad 22:73–87. Google Scholar

10.

J. S. Farris 1989. The retention index and the rescaled consistency index. Cladistics 5:417–419. Google Scholar

11.

J. Felsenstein 1985. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 39:783–791. Google Scholar

12.

J. Felsenstein 1993. PHYLIP. 3.75c edn.Department of Genetics SK-50, University of Washington. Seattle. Google Scholar

13.

P. Gopalakrishnakone and L. M. Chou . eds. 1990. Snakes of Medical Importance (Asia-Pacific region). Venom and Toxin Research Group. Singapore. Google Scholar

14.

B. Groombridge 1986. Comments on the M. pterygoideus glandulae of crotaline snakes (Reptilia: Viperidae). Herpetologica 42:449–457. Google Scholar

15.

P. J. Heise, L. R. Maxon, H. G. Dowling, and S. B. Hedges . 1995. Higher-level snake phylogeny inferred from mitochondrial DNA sequences of 12S rRNA and 16S rRNA genes. Mol Biol Evol 12:259–265. Google Scholar

16.

D. M. Hillis and J. P. Huelsenbeck . 1992. Signal, noise, and reliability in molecular phylogenetic analysis. J Hered 83:189–195. Google Scholar

17.

T. Hikida and H. Ota . 1997. Biogeography of reptiles in the subtropical East Asian islands. In “Proceedings of the Symposium on the Phylogeny, Biogeography and Conservation of Fauna and Flora of East Asian Region”Eds by K. Y. Lue and T. H. Chen , editors. National Taiwan Normal University. Taipei. pp. 11–28. Google Scholar

18.

A. R. Hoge and S. A. L. Romano-Hoge . 1981. Poisonous snakes of the world. I. Checklist of the pitvipers: Viperoidea, Viperidae, Crotalinae. Mem Inst Butantan 42/43:179–309. Google Scholar

19.

A. R. Hoge and S. A. L. Romano-Hoge . 1983. Notes on micro and ultrastructure of “Oberhautschen” in Viperoidea. Mem Inst Butantan 44/45:81–118. Google Scholar

20.

M. Honda, H. Ota, M. Kobayashi, and T. Hikida . 1999. Phylogenetic relationships of Australian skinks of the Mabuya group (Reptilia: Scincidae) inferred from mitochondrial DNA sequences. Genes Genet Syst 74:135–139. Google Scholar

21.

T. H. Jukes and C. R. Cantor . 1969. Evolution of protein molecules. In “Mammalian Protein Metabolism”. Ed by H. N. Munro , editor. Academic Press. New York. pp. 21–132. Google Scholar

22.

M. Kimura 1980. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 16:116–120. Google Scholar

23.

A. G. Kluge and J. S. Farris . 1969. Quantitative phyletics and the evolution of anurans. Syst Zool 18:1–32. Google Scholar

24.

A. Knight, L. D. Densmore III, and E. D. Rael . 1992. Molecular systematics of the Agkistrodon complex. In “Biology of Pitvipers” Eds by J. A. Campbell, E. D. Brodie Jr, and Tyler Selva , editors. pp. 49–69. Google Scholar

25.

A. Knight, D. Styer, S. Pelikan, J. A. Campbell, L. D. Densmore III, and D. P. Mindell . 1993. Choosing among hypotheses of rattlesnake phylogeny: A best-fit rate test for DNA sequence data. Syst Biol 42:356–367. Google Scholar

26.

K. Koba 1962. Studies on the Snakes of the Genus Trimeresurus of the Amami and Tokara Islands, Japan. Japan Society for the Promotion of Science. Tokyo. in Japanese with English summary. Google Scholar

27.

T. D. Kocher, W. K. Thomas, A. Meyer, S. V. Edwards, S. Paabo, F. X. Villablanca, and A. C. Wilson . 1989. Dynamics of mitochondrial DNA evolution in animals: Amplification and sequencing with conserved primers. Proc Natl Acad Sci USA 86:6196–6200. Google Scholar

28.

F. Kraus, D. G. Mink, and W. M. Brown . 1996. Crotaline intergeneric relationships based on mitochondrial DNA sequence data. Copeia 1996:763–773. Google Scholar

29.

A. Malhotra and R. S. Thorpe . 1997. New perspectives on the evolution of south-east Asian pitvipers (genus Trimeresurus) from molecular studies. In “Venomous Snakes: Ecology, Evolution and Snakebite”. Eds by R. S. Thorpe, W. Wüster, and A. Malhotra , editors. Zoological Society of London. London. pp. 115–128. Google Scholar

30.

P. Maslin 1942. Evidence for the separation of the crotalid genera Trimeresurus and Bothrops, with a key to the genus Trimeresurus. Copeia 1942:18–24. Google Scholar

31.

R. W. McDiarmid, J. A. Campbell, and T. Toure . 1999. Snake Species of the World. A Taxonomic and Geographic References. Vol. 1. Herpetologists' League. Washington, DC. Google Scholar

32.

M. Nei 1978. Estimation of average heterozygosity and genetic distance from a small sample number of individuals. Genetics 89:583–590. Google Scholar

33.

H. Ota 1991. Systematics and biogeography of terrestrial reptiles of Taiwan. In “Proceedings of the First International Symposium on Wildlife Conservation, ROC”Eds by Y. S. Lin and K. H. Chang , editors. Council of Agriculture. Taipei. pp. 47–112. Google Scholar

34.

H. Ota 1997. Historical biogeographical implications in the variation and diversity of amphibians and reptiles in Taiwan. In “Proceedings of the Symposium on the Phylogeny, Biogeography and Conservation of Fauna and Flora of East Asian Region”Eds by K. Y. Lue and T. H. Chen , editors. National Taiwan Normal University. Taipei. pp. 75–86. Google Scholar

35.

H. Ota 1998. Geographic pattern of endemism and speciation in amphibians and reptiles of the Ryukyu Archipelago, Japan, with special reference to their paleogeographical implications. Res Popul Ecol 40:189–204. Google Scholar

36.

H. Ota, M. Honda, M. Kobayashi, S. Sengoku, and T. Hikida . 1999. Phylogenetic relationships of eublepharid geckos (Reptilia: Squamata): A molecular approach. Zool Sci 659–666. Google Scholar

37.

C. L. Parkinson 1999. Molecular systematics and biogeographical history of pitvipers as determined by mitochondrial ribosomal DNA sequences. Copeia 1999:576–586. Google Scholar

38.

N. Saitou and M. Nei . 1987. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406–425. Google Scholar

39.

H. B. Shaffer, P. Meylan, and M. L. McKnight . 1997. Tests of turtle phylogeny: molecular, morphological, and paleontological approaches. Syst Biol 46:235–268. Google Scholar

40.

D. L. Swofford 1998. PAUP: A Phylogenetic Analysis Using Parsimony, version 4.0 (preliminary release). Sinauer. Sunderland, Massachusetts. Google Scholar

41.

K. Tamura and M. Nei . 1993. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10:512–526. Google Scholar

42.

J. D. Thompson, D. G. Higgins, and T. J. Gibson . 1994. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acid Res 22:4673–4680. Google Scholar

43.

M. Toda, M. Nishida, M. C. Tu, T. Hikida, and H. Ota . 1999. Genetic variation, phylogeny and biogeography of the pitvipers of the genus Trimeresurus sensu lato (Reptilia: Viperidae) in the subtropical East Asian islands. In “Tropical Island Herpetofauna: Origin, Current Diversity, and Conservation”. Ed by H. Ota , editor. Elsevier Science. Amsterdam. pp. 249–270. Google Scholar

44.

M. Toriba 1989. Review: Snakes of the Orient, a Checklist. By R. G. Kenneth and Welch . 1988. Krieger, Melbourne, Florida. Jpn J Herpetol 13:21–25. in Japanese. Google Scholar

45.

M. Toriba 1993. The genera Ovophis and Trimeresurus. In “Endoglyphs and Other Venomous Snakes of the World: A Checklist” Eds by P. Golay, H. M. Smith, D. G. Broadley, J. R. Dixon, C. McCarthy, J. C. Rage, B. Schatti, M. Toriba, S. A. Azemiops, and Geneva , editors. pp. 81–83.94–108. Google Scholar

46.

N. Vidal, G. Lecointre, J. C. Vie, and J. P. Gasc . 1997. Molecular systematics of pitvipers: paraphyly of the Bothrops complex. CR Acad Sci III 320:95–101. Google Scholar

47.

H. Y. Wang and M. C. Tu . 1997. Molecular evidence that the endemic snake, Trimeresurus gracilis, in Taiwan is not evolved in situ. In “Proceedings of the Symposium on the Phylogeny, Biogeography and Conservation of Fauna and Flora of East Asian Region”Eds by K. Y. Lue and T. H. Chen , editors. National Taiwan Normal University. Taipei. pp. 101–108. Google Scholar

48.

W. Wüster, P. Golay, and D. A. Warrell . 1997. Synopsis of recent developments in venomous snake systematics. Toxicon 35:319–340. Google Scholar

49.

H. T. Yu 1995. Patterns of diversification and genetic population structure of small mammals in Taiwan. Biol J Linnean Soc 55:69–89. Google Scholar

50.

F. J. Zhang 1993. Division of the genus Trimeresurus (sensu lato) based on the morphology of their skulls. In “Proceedings of the First Asian Herpetological Meeting”Eds by E. M. Zhao, B. H. Chen, and T. J. Papenfuss , editors. China Forestry Press. Beijing. pp. 48–57. 373. in Chinese with English abstract. Google Scholar

51.

F. J. Zhang 1995. Study on the taxonomy and phylogeny between genera of Trimeresurus group. Acta Herpetol Sinica 4/5:105–111. in Chinese with English abstract. Google Scholar

52.

E. M. Zhao and K. Adler . 1993. Herpetology of China. Society for the Study of Amphibians and Reptiles. Athens, Oxford. Google Scholar

53.

E. M. Zhao, M. Huang, and Y. Zong . Eds. 1998. Fauna Sinica. Reptilia. Vol. 3. Squamata, Serpentes. Science Press. Beijing. in Chinese. Google Scholar

Appendices

Appendix 1

Specimens examined in this study. Catalogue numbers of voucher specimens are given in parentheses.

  • Trimeresurus mucrosquamatus: Hongya, Sichuan, continental China (uncatalogued specimen); Chishang, Kaohsiung, Taiwan (National Taiwan Normal University [NTNU] B201574).

  • T. elegans: Yaeyama Group, Okinawa, Japan (NTNU B201580).

  • T. tokarensis: Kodakarajima Island, Tokara Group, Kagoshima, Japan (Kyoto University Zoological Collection [KUZ] 21104).

  • T. flavoviridis: Kumejima Island, Okinawa Group, Okinawa, Japan (KUZ 45840).

  • T. jerdonii: Hongya, Sichuan, continental China (uncatalogued specimen).

  • T. monticola makazayazaya: Yangmingshang, Taipei, Taiwan (NTNU B200800; Yizhang, Hunan, continental China (uncatalogued specimen).

  • T. m. monticola: Kunming, Yunnan, continental China (NTNU B201401).

  • T. stejnegeri: Fushan, Ilan, Taiwan (NTNU B201588); Kunming, Yunnan, continental China (uncatalogued specimen).

  • T. albolabris: Foochow, Fujian, continental China (NTNU B201408).

  • T. puniceus: Sumatra, Indonesia (private collection of G. Vogel).

  • T. gracilis: Alishan, Chiayi, Taiwan (uncatalogued specimen).

  • T. okinavensis: Tokunoshima Island, Amami Group, Kagoshima, Japan (KUZ 45871).

  • T. popeiorum: Sakaerat, Thailand (NTNU B200511).

  • T. trigonocephalus: Sri Lanka (private collection of G. Vogel).

  • Tropidolaemus wagleri: Kuching, Sarawak, Malaysia (NTNU B200512).

  • Deinagkistrodon acutus: Hoping, Taichung, Taiwan (uncatalogued specimen)

  • Vipera russelii: Fanshan, Pintung, Taiwan (NTNU B201587).

Appendix 2

Appendix 2.

Aligned sequences of a part of the mitochondrial DNA 12S rRNA gene investigated. Dot indicates an identity with the first sequence; dash denotes an indel. Abbreviations are: MA, Trimeresurus monticola makazayazaya from Taiwan; MO, T. m. monticola; MU, T. mucrosquamatus from continental China; EL, T. elegans; TO, T. tokarensis; FL, T. flavoviridis; JE, T. jerdonii; GR, T. gracilis; OK, T. okinavensis; PO, T. popeiorum; AL, T. albolabris; ST, T. stejnegeri from Taiwan; PU, T. puniceus; TR, T. trigonocephalus, WA, T. wagleri; DA, Deinagkistrodon acutus; VR, Vipera russelii.

i0289-0003-17-8-1147-f1001.gif

Appendix 2.

Continued.

i0289-0003-17-8-1147-f1002.gif
Ming-Chung Tu, Hurng-Yi Wang, Mung-Pei Tsai, Mamoru Toda, Wen-Jen Lee, Fu-Ji Zhang, and Hidetoshi Ota "Phylogeny, Taxonomy, and Biogeography of the Oriental Pitvipers of the Genus Trimeresurus (Reptilia: Viperidae: Crotalinae): A Molecular Perspective," Zoological Science 17(8), 1147-1157, (1 November 2000). https://doi.org/10.2108/zsj.17.1147
Received: 1 March 2000; Accepted: 1 May 2000; Published: 1 November 2000
JOURNAL ARTICLE
11 PAGES


SHARE
ARTICLE IMPACT
Back to Top