Molecular Evaluation of Phylogenetic Significances in the Highly Divergent Karyotypes of the Genus Gonocephalus (Reptilia: Agamidae) from Tropical Asia (original) (raw)

INTRODUCTION

The agamid genus Gonocephalus Kaup, 1825, originally considered as consisting of 29 species distributed in Southeast Asia, Australia, New Guinea and the Solomons (Wermuth, 1967), was highlighted as a zoogeographical exception because of its occurrence on both sides of the Wallacia, a border of the Oriental and Australian faunal realms (e.g., Darlington, 1957). In his unpublished dissertation, Moody (1980), on the basis of cladistic analyses using morphological characters, argued that species of the genus Gonocephalus from Australia, New Guinea and the Solomons were not actually closely related to the Oriental congeners. He thus re-assigned the former to the resurrected genus Hypsilurus Peters, 1867. Such an arrangement, published by Welch et al. (1990), was favored by recent immunological (Baverstock and Donnellan, 1990; King, 1990), karyological (Ota et al., 1992), electron-microscopic (Ananjeva and Matveyeva-Dujsebayava, 1996), and molecular studies (Honda et al., 2000; Macey et al., 2000).

Recent karyological studies further revealed prominent chromosomal variation within Gonocephalus sensu stricto: four of the five species examined shared a 2n = 42 karyo-type including 22 biarmed macrochromosomes, whereas the remainder, G. robinsonii, had a 2n = 32 karyotype including no more than 12 macrochromosomes (Ota et al., 1992; Diong et al., 2000). Based on the fact that the former arrangement is exceptional as of the agamid karyotype, Diong et al. (2000) interpreted it as a synapomorph of Gonocephalus exclusive of G. robinsonii, and went so far as to argue for the possible paraphyletic nature of the genus. They, on the basis of similarities in the number of microchromosomes between G. robinsonii and a few Australian agamids, also noted that G. robinsonii may be an Asiatic representative of the Australian agamid radiation, like Physignathus cocincinus (Honda et al., 2000; Macey et al., 2000). However, these hypotheses seriously suffer due to the lack of appropriate phylogenetic analysis of chromo-some characters and of comparable karyological data for many other related genera.

In the present study, we examined the phylogenetic relationships of the two karyomorphs of Gonocephalus and other agamid genera by analyzing mitochondrial DNA sequence data. Our purposes are to examine consistency between patterns in chromosomal and phyletic diversification within the Gonocephalus, and to test the closer affinity of G. robinsonii to the Australian agamids as hypothesized by Diong et al. (2000).

MATERIALS AND METHODS

Three species of the Gonocephalus sensu stricto, G. chamaeleontinus from West Malaysia and Indonesia, G. miotympanum from Borneo and G. robinsonii from West Malaysia (Welch et al., 1990; Manthey and Grossmann, 1997), were newly examined in the present study (see Appendix for locality data of specimens examined). We incorporated published data representing the six major agamid groups of Moody (1980) (Groups I–VI), the Chamaeleonidae (Bradypodion fischeri), the Polychridae (Anolis carolinensis), and the Iguanidae (Iguana iguana) (see Appendix for accession numbers). Of these, the last three species were regarded as out-groups, because these families are considered to be basal to the Agamidae (Frost and Etheridge, 1989; Macey et al., 1997; Honda et al., 2000)

Extraction, amplification and sequencing of DNA are described in detail elsewhere (Honda et al., 1999a, b). A part of mitochondrial 12S and 16S rRNA genes consisting of approximately 860 base pairs (bp) were amplified using the polymerase chain reaction (PCR) with primers L1091, H1478, L2606, and H3056 (Kocher et al., 1989; Hedges et al., 1993). Alignments for DNA sequences were unambiguously determined based on maximum nucleotide similarity following Honda et al. (2000).

The neighbor-joining (NJ) method (Saitou and Nei, 1987) was applied to infer relationships among taxa on the basis of a pairwise matrix of the distance from Kimura's (1980) two-parameter model, using CLUSTAL X 1.8 (Thompson et al., 1994). The maximum-likelihood (ML, empirical base frequencies and equal rate substitution model) and maximum-parsimony (MP, no bias between transition and transversion) analyses were also conducted using heuristic search option of PAUP* 4.0b (Swofford, 1998). In these three analyses, gap sites were excluded, and confidences of branched were assessed by bootstrap resamplings (Felsenstein, 1985).

RESULTS

The amplified fragment of the 12S and 16S rRNA genes consisted of 862 total sites. Of these, however, 201 sites involved in the insertion or deletion were excluded from phylogenetic analyses. The NJ dendrogram derived from aligned sequences is shown in Fig. 1A. The monophyly of the Agamidae was supported in all bootstrap iterations (node 1: 100%). The ingroup portion of this dendrogram was divided into two major lineages. One of the major clusters (node 2: 83%) further split into two subclusters (nodes 4, 5), of which node 4 (80%) consisted of West Asian and North African primitive agamids (Leiolepis and Uromastyx: Group I sensu Moody, 1980), whereas node 5 (84%) accommodated Australian genera including Hypsilurus that had been previously assigned to Gonocephalus (Groups II–IV complex). The other major cluster (node 3: 71%) split into two subclusters (nodes 6, 7), of which node 6 (100%) consisted of West Asian and African derived agamids (Agama and Phrynocephalus: Group VI). Node 7 (100%), on the other hand, contained Southeast Asian derived agamids (Acanthosaura, Aphaniotis, Calotes, Draco, Gonocephalus, Japalura, Phoxophrys and Ptyctolaemus: Group V). Within this node, G. chamaeleontinus, G. grandis and G. miotympanum (node 8: 100%), and G. robinsonii and J. polygonata (node 9: 70%) constituted exclusive clusters, respectively.

Fig. 1

(A) Neighbor-joining (NJ) dendrogram derived from distance matrix from 12S and 16S rRNA sequence data. Numbers beneath branches are BPs at least 50% of the 1,000 bootstrap replications. Nodes with bold numbers indicate relationships referred to in the text. Bar equals 0.1 unit of Kimura's two-parameter distance. Published data for the diploid chromosome numbers (2n) are given in parentheses (DeSmet, 1981; Witten, 1983; Solleder and Schmid, 1988; Ota et al., 1992; Ota and Hikida, 1989; Diong et al., 2000; Ota, unpubl data). Although Groups II–IV complex is mainly distributed in Australasia, Physignathus cocincinus is occurs in the continental part of Southeast Asia. (B) Maximum-likelihood (ML) dendrogram (ln likelihood=−11727.6). Branches without BP values were not supported in≥50% of the 100 bootstrap replicates. Bar equals 0.1 unit. Bold numbers above branches are identical with those in NJ and MP dendrograms. (C) Maximum parsimony (MP) dendrogram using heuristic option (one parsimonious tree, 2,418 steps, 448 bp informative under the condition of parsimony, consistency index=0.45, homoplasy index=0.55, retention index=0.48). Branches without BP values were not supported in≥50% of the 1,000 bootstrap replicates.

Relationships depicted as a result of ML (Fig. 1B) and MP analyses (Fig. 1C) showed no inconsistency with those expressed in the NJ dendrogram in terms of topology of nodes 1–9, except for the absence of node 3 in MP. However, when the Kishino and Hasegawa (1989) test and Templeton (1983) test were applied to ML and MP dendrograms, respectively, topologies in Fig. 1B and Fig. 1C were not significantly different from those in alternative hypotheses in which Gonocephalus was constrained monophyletic (P=0.130 and 0.167, respectively).

DISCUSSION

Several authors pointed out that there are two typical karyomorphs in the family Agamidae, of which one consists of 2n=34 or 36 chromosomes including 12 biarmed macro-chromosomes and 22 or 24 microchromosomes, whereas the other of 2n=46 or 48, all telocentric chromosomes without a distinct size break (Bickham, 1984; King, 1981; Moody and Hutterer, 1978; Olmo, 1986; Ota and Hikida, 1989; Solleder and Schmid, 1988; Witten, 1983). One of these karyomorphs is considered to be derived from the other through a series of Robertsonian rearrangements of macro-chromosomes, sometimes accompanied by addition or deletion of a pair of microchromosomes (Bickham, 1984; King, 1981). Based on these empirical assumptions, Ota et al. (1992) and Diong et al. (2000) hypothesized that the 2n=42 karyomorph with 22 biarmed macrochromosomes and 20 microchromosome of Gonocephalus bellii, G. chamaeleontinus, G. grandis, G. liogaster and G. miotympanum represents a highly specialized state. The latter authors also considered such an arrangement to be a strong evidence for the monophyly of those species against G. robinsonii, a species having 2n=32 karyotype including only 12 biarmed macro-chromosomes, and further assumed the latter species to be distant from other congeners phylogenetically.

Present results, consistently indicating the monophyly of the “2n=42 species” of Gonocephalus against G. robinsonii and other related genera with high bootstrap values, offer a substantial support to these hypotheses, although the number of species examined in the present study is too small to draw any definite conclusion. The states of morphological characters shared between G. robinsonii and other congeners (i.e., diagnostic characters of Gonocephalus), such as the presence of tubercular scales among otherwise homologous granular scales (Manthey and Grossmann, 1997), thus appear to be derived from convergence or represent symplesiomorphy.

As was pointed out by Diong et al. (2000), chromosomal arrangement of G. robinsonii most resembles that of a few Australian agamids in that it has only 20 microchromosomes besides 12 biarmed macrochromosomes. Present results, however, negate their close affinity, and indicate the origin of G. robinsonii in an Asiatic agamid radiation. Simple deletion of a pair of microchromosomes seems to be responsible for the emergence of karyotype of G. robinsonii.

Acknowledgments

We thank M. Matsui, K. Araya, T. Ueda, L. David, A. A. Hamid, the staff of the Forest Research Center, Sepilok and the staff of National Park and Wildlife and Forest Research Sections, Forest Department of Sarawak for providing us with various helps and encouragements during our fieldwork, and C.-H. Diong for pulling our attention to chromosomal variation within the genus Gonocephalus. Special thanks are due N. Satoh and members of his laboratory for continuous support for our laboratory experiments. Experiments were also carried out using the facilities of the Kyoto University Museum through the courtesy of T. Nakabo and M. Motokawa.

H. Ota and T. Hikida are especially grateful to T. Hidaka and M. Matsui for providing opportunities to visit Malaysia for sampling. Our research was partially supported by Grants-in-Aid from the Japan Ministry of Education, Science, Sports and Culture (Overseas Researches Nos. 01041051 to Hidaka and 0404068 to Matsui).

REFERENCES

N. B. Ananjeva and T. N. Matveyeva-Dujsebayava . 1996. Some evidence of Gonocephalus species complex divergence based on skin sense organ morphology. Russ J Herpetol 3:82–88. Google Scholar

J. W. Bickham 1984. Patterns and models of chromosomal evolution in reptiles. In “Chromosome in Evolution of Eukaryotic Groups Vol 2”. Eds by A. K. Sharma and A. Sharma . CRC Press. Florida. pp. 13–40. Google Scholar

P. B. Baverstock and S. C. Donnellan . 1990. Molecular evolution in Australian dragons and skinks: a progress report. Mem Queensland Mus 29:323–331. Google Scholar

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

W. H. O. DeSmet 1981. Description of the orcein stained karyotypes of 27 lizards species (Lacertilia, Reptilia) belonging to the family Iguanidae, Agamidae, Chamaeleontidae and Gekkonidae (Ascalabota). Acta Zool Pathol Antverp 76:35–72. Google Scholar

C-H. Diong, M-H. Low, E-C. Tan, H-S. Yong, T. Hikida, and H. Ota . 2000. On the monophyly of the agamid genus Gonocephalus Kaup, 1925 (Reptilia: Agamidae): a chromosomal perspective. Cur Herpetol 19:71–79. Google Scholar

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

D. R. Frost and R. Etheridge . 1989. A phylogenetic analysis and taxonomy of iguanian lizards (Reptilia: Squamata). Univ Kansas Mus Nat Hist Misc Pub 81:1–65. Google Scholar

S. B. Hedges, R. A. Nussbaum, and L. R. Maxson . 1993. Caecilian phylogeny and biogeography inferred from mitochondrial DNA sequences of the 12S rRNA and 16S rRNA genes (Amphibia: Gymnophiona). Herpetol Monogr 7:64–76. Google Scholar

M. Honda, H. Ota, M. Kobayashi, J. Nabhitabhata, H-S. Yong, and T. Hikida . 1999a. Phylogenetic relationships of the flying lizards, genus Draco (Reptilia: Agamidae). Zool Sci 16:535–549. Google Scholar

M. Honda, H. Ota, M. Kobayashi, J. Nabhitabhata, H-S. Yong, and T. Hikida . 1999b. Evolution of Asian and African lygosomine skinks of the Mabuya group (Reptilia: Scincidae): a molecular perspective. Zool Sci 16:979–984. Google Scholar

M. Honda, H. Ota, M. Kobayashi, J. Nabhitabhata, H-S. Yong, S. Sengoku, and T. Hikida . 2000. Phylogenetic relationships of the family Agamidae (Reptilia: Iguania) inferred from mitochondrial DNA sequences. Zool Sci 17:527–537. Google Scholar

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

M. King 1981. Chromosomal change and speciation in lizards. In “Evolution and Speciation, Essay in Honor of MJD White”. Eds by W. R. Atchley and D. Woodruff . Cambridge Press. Cambridge. pp. 262–285. Google Scholar

M. King 1990. Chromosomal and immunogenetic data: a new respective on the origin of Australia's reptile. In “Cytogenetics of Amphibians and Reptiles”. Ed by E. Olmo Birkhauser Verlag. Basel. pp. 153–180. Google Scholar

H. Kishino and M. Hasegaw . 1989. Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in Hominoidea. J Mol Evol 29:170–179. Google Scholar

T. D. Kocher, W. K. Thomas, A. Meyer, S. V. Edwards, S. Pääbo, F. X. Villablanca, and A. C. Wilson . 1989. Dynamics of mitochondrial DNA evolution in animals: Amplifications and sequencing with conserved primers. Proc Nat Acad Sci USA 86:189–191. Google Scholar

J. R. Macey, A. Larson, N. B. Ananjeva, and T. J. Papenfuss . 1997. Evolutionary shifts in three major structural features of the mitochondrial genome among iguanian lizards. J Mol Evol 44:660–674. Google Scholar

J. R. Macey, J. A. Schulte, A. Larson, N. B. Ananjeva, Y. Wang, R. Pethiyagoda, Rasteger-Pouyani, and T. J. Papenfuss . 2000. Evaluating trans-Tethys migration: an example using Acrodont lizard phylogenetics. Syst Biol 49:233–256. Google Scholar

U. Manthey and W. Grossmann . 1997. Amphibien & Reptilien Südostasiens. Natur und Tier Verlag, Münster. Google Scholar

S. M. Moody 1980. Phylogenetic and historical biogeographical relationships of the genera in family Agamidae (Reptilia: Lacertilia). Unpubl PhD Thesis. Univ Michigan. Michigan. Google Scholar

S. M. Moody and H. R. Hutterer . 1978. Karyotypes of the agamid lizard Lyriocephalus scutatus (L 1758), with a brief review of the chromosomes of the lizard family Agamidae. Bönn Zool Beitr 29:165–170. Google Scholar

Olmo 1986. Reptilia. In “Animal Cytogenetics Vol 4 Chordata 3”. Ed by B. John Gebruder Borntraeger. Berlin and Stuttgart. pp. 1–100. Google Scholar

H. Ota and T. Hikida . 1989. Karyotypes of three species of the genus Draco (Agamidae: Lacertilia) from Sabah, Malaysia. Jpn J Herpetol 13:1–6. Google Scholar

H. Ota, M. Matsui, H. Hikida, and A. Mori . 1992. Extreme karyotypic divergence between species of the genus Gonocephalus (Reptilia: Squamata: Agamidae) from Borneo and Australia. Herpetologica 48:120–124. Google Scholar

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

E. Solleder and M. Schmid . 1988. Cytogenetic studies on Sauria (Reptilia) I. Mitotic chromosomes of the Agamidae. Amphibia-Retilia 9:301–310. Google Scholar

D. L. Swofford 1998. A Phylogenetic Analysis Using Parsimony (*and Other Method) ver 4.0. Massachusetts. Google Scholar

A. Templeton 1983. Phylogenetic inference from restrction endonuclease clavage site maps with particular reference to the evolution of humans and apes. Evolution 37:221–244. Google Scholar

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

K. R. G. Welch, P. S. Cooke, and A. S. Wright . 1990. Lizards of the Orient: A Checklist. Krieger Publ Co. Malabar, Florida. Google Scholar

G. J. Wermuth 1967. Liste der rezenten Amphibien und Reptilien. Agamidae. Das Tierreich. Walter de Gruyter. Berlin. Google Scholar

G. J. Witten 1983. Some karyotypes of Australian agamids (Reptilia: Lacertilia). Aust J Zool 31:533–540. Google Scholar

Appendices

APPENDIX.

DDBJ accession numbers of mitochondrial 12S and 16S rRNA gene sequences used in this study (given in parentheses). Localities and catalogue numbers are also given for specimens newly sequenced this time. These, as well as previously used specimens (Honda et al., 2000), were deposited in the Herpetological Collection of Department of Zoology, Kyoto University (KUZ).

Acanthosaura crucigera ( , ). Agama stellio ( , ). Aphaniotis fusca ( , ). Calotes versicolor ( , ). Chlamydosaurus kingii ( , ). Draco volans ( , ). Gonocephalus chamaeleontinus: imported by a pet dealer (detailed locality unknown), KUZ R50574 ( , ). G. grandis: ( , ). G. miotympanum: Gunong Gading N. P., Sarawak, Borneo, KUZ R27058 ( , ). G. robinsonii: Cameron highland, Peninsular Malaysia, KUZ R21343 ( , ). Hypsilurus godeffroyi ( , ). Japalura polygonata polygonata ( , ). Leiolepis belliana ( , ). Lophognathus temporalis ( , ). Phoxophrys nigrilabris ( , ). Phrynocephalus axillaris ( , ). Physignathus cocincinus ( , ). Physignathus lesueurii ( , ). Pogona vitticeps ( , ). Ptyctolaemus phuwuanensis ( , ). Uromastyx aegyptia ( , ). Bradypodion fischeri ( , ). Anolis carolinensis ( , ). Iguana iguana ( , ).