Expression of Hoxa-11 and Hoxa-13 in the pectoral fin of a basal ray-finned fish, Polyodon spathula: implications for the origin of tetrapod limbs

Summary Paleontological and anatomical evidence suggests that the autopodium (hand or foot) is a novel feature that distinguishes limbs from fins, while the upper and lower limb (stylopod and zeugopod) are homologous to parts of the sarcopterygian paired fins. In tetrapod limb development Hoxa-11 plays a key role in differentiating the lower limb and Hoxa-13 plays a key role in differentiating the autopodium. It is thus important to determine the ancestral functions of these genes in order to understand the developmental genetic changes that led to the origin of the tetrapod autopodium. In particular it is important to understand which features of gene expression are derived in tetrapods and which are ancestral in bony fishes. To address these questions we cloned and sequenced the Hoxa-11 and Hoxa-13 genes from the North American paddlefish, Polyodon spathula, a basal ray-finned fish that has a pectoral fin morphology resembling that of primitive bony fishes ancestral to the tetrapod lineage. Sequence analysis of these genes shows that they are not orthologous to the duplicated zebrafish and fugu genes. This implies that the paddlefish has not duplicated its HoxA cluster, unlike zebrafish and fugu. The expression of Hoxa-11 and Hoxa-13 in the pectoral fins shows two main phases: an early phase in which Hoxa-11 is expressed proximally and Hoxa-13 is expressed distally, and a later phase in which Hoxa-11 and Hoxa-13 broadly overlap in the distal mesenchyme of the fin bud but are absent in the proximal fin bud. Hence the distal polarity of Hoxa-13 expression seen in tetrapods is likely to be an ancestral feature of paired appendage development. The main difference in HoxA gene expression between fin and limb development is that in tetrapods (with the exception of newts) Hoxa-11 expression is suppressed by Hoxa-13 in the distal limb bud mesenchyme. There is, however, a short period of limb bud development where Hoxa-11 and Hoxa-13 overlap similarly to the late expression seen in zebrafish and paddlefish. We conclude that the early expression pattern in tetrapods is similar to that seen in late fin development and that the local exclusion by Hoxa-13 of Hoxa-11 from the distal limb bud is a derived feature of limb developmental regulation.     

The tetrapod limb: a hypothesis on its origin

A wrist joint and structures typical of the hand, such as digits, however, are absent in [Eustenopteron] (Andrews and Westoll, '68, p 240). Great changes must have been undergone during evolution of the ankle joint; the small number of large bones in the fin must somehow have developed into a large number of small bones, and it is very difficult to draw homologies in this region, or even be certain of what is being compared (Andrews and Westoll, '68, p 268). The tetrapod limb is one of the major morphological adaptations that facilitated the transition from an aquatic to a terrestrial lifestyle in vertebrate evolution. We review the paleontological evidence for the fin-limb transition and conclude that the innovation associated with evolution of the tetrapod limb is the zeugopodial-mesopodial transition, i.e., the evolution of the developmental mechanism that differentiates the distal parts of the limb (the autopodium, i.e., hand or foot) from the proximal parts. Based on a review of tetrapod limb and fish fin development, we propose a genetic hypothesis for the origin of the autopodium. In tetrapods the genes Hoxa-11 and Hoxa-13 have locally exclusive expression domains along the proximal-distal axis of the limb bud. The junction between the distal limit of Hoxa-11 expression and of the proximal limit of Hoxa-13 expression is involved in establishing the border between the zeugopodial and autopodial anlagen. In zebrafish, the expression domains of these genes are overlapping and there is no evidence for an autopodial equivalent in the fin skeleton. We propose that the evolution of the derived expression patterns of Hoxa-11 and Hoxa-13 may be causally involved in the origin of the tetrapod limb.     

Hoxa-11 and Hoxa-13 are involved in repression of MyoD during limb muscle development

Under the influence of the limb mesenchyme, Hoxa-11 is expressed in migrating and proliferating premyoblasts in the limb field and Hoxa-13 is induced in subdomains of congregated limb muscle masses. To evaluate the roles of Hoxa-11 and Hoxa-13 in myogenesis of the limb, we performed electroporation in ovo to force expression of these Hox genes in limb muscle precursors. In the presence of ectopic Hoxa-11, expression of MyoD was blocked transiently. In C2C12 myoblasts, transfection of Hoxa-11 also repressed the expression of endogenous MyoD. Forced expression of Hoxa-13 resulted in more pronounced repression of MyoD in both limb and C2C12 myoblasts. In contrast, targeted disruption of Hoxa-13 gave rise to enhanced expression of MyoD in the flexor carpi radialis muscle, a forearm muscle that normally expressed Hoxa-13. These results suggest that Hoxa-11 and Hoxa-13 are involved in the negative regulation of MyoD expression in limb muscle precursors.     

Phylogenetic relationships and limb loss in sub-Saharan African scincine lizards (Squamata: Scincidae)

Skinks are the largest family of lizards and are found worldwide in a diversity of habitats. One of the larger and more poorly studied groups of skinks includes members of the subfamily Scincinae distributed in sub-Saharan Africa. Sub-Saharan African scincines are one of the many groups of lizards that show limb reduction and loss, and the genus Scelotes offers an excellent opportunity to look at limb loss in a phylogenetic context. Phylogenetic relationships were reconstructed for a total of 52 taxa representing all subfamilies of skinks as well as other Autarchoglossan families using sequence from six gene regions including; 12S, 16S, and cytochrome b (mitochondrial), as well as alpha-Enolase, 18S, and C-mos (nuclear). The family Scincidae is recovered as monophyletic and is the sister taxon to a (Cordylidae+Xantusiidae) clade. Within skinks the subfamily Acontinae is monophyletic and sister group to all remaining skinks. There is no support for the monophyly of the subfamilies Lygosominae and Scincinae, but sub-Saharan African scincines+Feylinia form a well supported monophyletic group. The monophyly of Scelotes is confirmed, and support is found for two geographic groups within the genus. Reconstructions of ancestral states for limb and digital characters show limited support for the reversal or gain of both digits and limbs, but conservative interpretation of the results suggest that limb loss is common, occurring multiple times throughout evolutionary history, and is most likely not reversible.     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

Temporal and Spatial Expression of Hoxa-2 During Murine Palatogenesis

1. Mice homozygous for a targeted mutation of the Hoxa-2 gene are born with a bilateral cleft of the secondary palate associated with multiple head and cranial anomalies and these animals die within 24 hr of birth (Gendron-Maguire et al., 1993; Rijli et al., 1993; Mallo and Gridley, 1996). We have determined the spatial and temporal expression of the Hoxa-2 homeobox protein in the developing mouse palate at embryonic stages E12, E13, E13.5, E14, E14.5, and E15.2. Hoxa-2 is expressed in the mesenchyme and epithelial cells of the palate at E12, but is progressively restricted to the tips of the growing palatal shelves at E13.3. By the E13.5 stage of development, Hoxa-2 protein was found to be expressed throughout the palatal shelf. These observations correlate with palatal shelf orientation and Hoxa-2 protein may play a direct or indirect role in guiding the palatal shelves vertically along side the tongue, starting with the tips of the palatal shelves at E13, followed by the entire palatal shelf at E13.5.4. As development progresses to E14, the stage at which shelf elevation occurs, Hoxa-2 protein is downregulated in the palatal mesenchyme but remains in the medial edge epithelium. Expression of Hoxa-2 continues in the medial edge epithelium until the fusion of opposing palatal shelves.5. By the E15 stage of development, Hoxa-2 is downregulated in the palate and expression is localized in the nasal and oral epithelia.6. In an animal model of phenytoin-induced cleft palate, we report that Hoxa-2 mRNA and protein expression were significantly decreased, implicating a possible functional role of the Hoxa-2 gene in the development of phenytoin-induced cleft palate.7. A recent report by Barrow and Capecchi (1999), has illustrated the importance of tongue posture during palatal shelf closure in Hoxa-2 mutant mice. This along with our new findings of the expression of the Hoxa-2 protein during palatogenesis has shed some light on the putative role of this gene in palate development.     

Limb development in the gekkonid lizard gonatodes albogularis: A reconsideration of homology in the lizard carpus and tarsus

Despite the attention squamate lizards have received in the study of digit and limb loss, little is known about limb morphogenesis in pentadactyl lizards. Recent developmental studies have provided a basis for understanding lizard autopodial element homology based on developmental and comparative anatomy. In addition, the composition and identity of some carpal and tarsal elements of lizard limbs, and reptiles in general, have been the theme of discussions about their homology compared to non‐squamate Lepidosauromorpha and basal Amniota. The study of additional embryonic material from different lizard families may improve our understanding of squamate limb evolution. Here, we analyze limb morphogenesis in the gekkonid lizard Gonatodes albogularis describing patterns of chondrogenesis and ossification from early stages of embryonic development to hatchlings. Our results are in general agreement with previous developmental studies, but we also show that limb development in squamates probably involves more chondrogenic elements for carpal and tarsal morphogenesis, as previously recognized on the grounds of comparative anatomy. We provide evidence for the transitory presence of distal carpale 1 and intermedium in the carpus and tibiale, intermedium, distal centralia, and distal tarsale 2 in the tarsus. Hence, we demonstrate that some elements that were believed to be lost in squamate evolution are conserved as transitory elements during limb development. However, these elements do not represent just phylogenetic burden but may be important for the morphogenesis of the lizard autopodium. J. Morphol., 2010. © 2010 Wiley‐Liss, Inc.     

Partial tail loss has no severe effects on energy stores and locomotor performance in a lacertid lizard, <Emphasis Type="Italic">Takydromus septentrionalis</Emphasis>

Many species of lizards use caudal autotomy as a defense strategy to avoid predation, but tail loss entails costs. These topics were studied experimentally in the northern grass lizard, Takydromus septentrionalis. We measured lipids in the three-tail segments removed from each of the 20 experimental lizards (adult females) initially having intact tails to evaluate the effect of tail loss on energy stores; we obtained data on locomotor performance (sprint speed, the maximal length traveled without stopping and the number of stops in the racetrack) for these lizards before and after the tail-removing treatments to evaluate the effect of tail loss on locomotor performance. An independent sample of 20 adult females that retained intact tails was measured for locomotor performance to serve as controls for successive measurements taken for the experimental lizards. The lipids stored in the removed tail was positively correlated with tailbase width when holding the tail length constant, indicating that thicker tails contained more lipids than did thinner tails of the same overall length. Most of the lipids stored in the tail were concentrated in the proximal portion of the tail. Locomotor performance was almost unaffected by tail loss until at least more than 71% of the tail (in length) was lost. Our data show that partial tail loss due to predatory encounters or other factors may not severely affect energy stores and locomotor performance in T. septentrionalis.     

Road edge effect on the abundance of the lizard <Emphasis Type="Italic">Gallotia galloti</Emphasis> (Sauria: Lacertidae) in two Canary Islands forests

Transportation infrastructure is a main cause of environmental change in forest landscapes worldwide. In the Canary Islands, a dense road system fragment the native Canarian pine and laurel forests causing potential changes in population densities of endemic lacertid lizards (genus Gallotia). Our aim was to assess road edge effects on relative abundance patterns of the endemic Gallotia galloti in both forests. We also explored the species–habitat relationships in this road-fragmentation context. We found that lizard relative density in relation to road edges differed between forests. Lizards were more abundant along edges and leeward interior, but virtually absent from the interior of the windward laurel forest. In the pine forest, lizards were present at three distances from edge, with a net decrease in abundance from edge to interior. These patterns may be explained partly by differences in vegetation structure regarding road proximity in each forest that potentially affect the helio- and thigmothermic character of G. galloti, and thus its habitat use. A general suggestion of this study is that road margins create corridors that may be used by native lizards for dispersal through inhospitable forest matrix. The high road density in Tenerife may have negative implications for the conservation of the genetic variability of G. galloti. At the island scale, increased communication between lizard populations through road corridors might increase homogenization of the gene pool. Ecological processes in which this lizard plays important roles may also experience changes along road edges.     

矛尾鱼Hoxa—11基因克隆及5′端序列分析

Hoxa-11基因调节鱼类鳍和四足动物肢的发育,在脊椎动物进化过程中起着重要的作用,利用人和鼠的Hoxa-11基因保守序列设计了两个兼并引物,通过PCR扩增到了矛尾鱼的Hoxa-11基因,经克隆和DNA序列分析,该片段为2065bp,包括绝大部分外显子Ⅰ,内含子和部分外显子Ⅱ,编码204个氨基酸,其氨基酸序列与人、鼠、鸡、蛙和斑马鱼的同源性分别为66.0%、67.6%、74.4%、72.8%和59.7%。外显子Ⅰ的长度从矛尾鱼到蛙、鸡、鼠和人呈现逐步上升趋势,人比矛尾鱼增长了16%,进一步分析,外显子Ⅰ可分为4个区域;两个高度保守区域,1个中度保守区域和1个可变区域,外显子Ⅰ的长度变化主要是由于可变区域内丙氨酸同聚物以及两侧富含甘氨酸和丝氨酸序列的累积。矛尾鱼只有1个由两个丙氨酸组成的同类物,蛙有1个由5个连续丙氨酸组成的同聚物,而鸡、鼠和人有3个丙氨酸同聚物,其中最大的同聚物由7个连续丙氨酸组成,而且在同聚物两侧出现了富含甘氨酸和丝氨酸序列。这表明可变区域可能与脊椎动物进化和鳍-肢转换过程中新功能的获得有关。同源异型盒所在的外显子Ⅱ区和剪接位点是高度保守的。内含子的长度变化较大,但在其内部也发现了两个高度保守的35bp和16bp的DNA片段,这两个片段在人、鼠、鸡、蛙和矛尾鱼中是完全相同的,这些序列的高度保守性提示其功能上的重要性。     

Allelic Combinations of Promoter and Exon 2 in <Emphasis Type="Italic">DQB1</Emphasis> in Dogs and Wolves

Polymorphism of PBRs of the major histocompatibility complex (MHC) genes is well recognized, but the polymorphism also extends to proximal promoter regions. Examining DQB1 variability in dogs and wolves, we identified 7 promoter variants and 13 exon 2 alleles among 89 dogs, including a previously unknown DQB1 exon 2 allele, and 8 promoter variants and 9 exon 2 alleles among 85 wolves. As expected from previous studies and from a close chromosomal location, strong linkage disequilibrium was demonstrated in both wolves and dogs by having significantly fewer promoter/exon 2 combinations than expected from simulations of randomized data sets. Interestingly, we noticed weaker haplotypic associations in dogs than in wolves. Dogs had twice as many promoter/exon 2 combinations as wolves and an almost 2-fold difference in the number of exon 2 alleles per promoter variant. This difference was not caused by an admixture of breeds in our group of dogs because the high ratio of observed to expected number of haplotypes persisted within a single dog breed, the German Shepherd. Ewens-Watterson tests indicated that both the promoter and exon 2 are under the balancing selection, and both regions appear to be more recently derived in the dog than in the wolf. Hence, although reasons for the differences are unknown, they may relate to altered selection pressure on patterns of expression. Deviations from normal MHC expression patterns have been associated with autoimmune diseases, which occur frequently in several dog breeds. Further knowledge about these deviations may help us understand the source of such diseases.     

Limbs in whales and limblessness in other vertebrates: mechanisms of evolutionary and developmental transformation and loss

We address the developmental and evolutionary mechanisms underlying fore- and hindlimb development and progressive hindlimb reduction and skeletal loss in whales and evaluate whether the genetic, developmental, and evolutionary mechanisms thought to be responsible for limb loss in snakes "explain" loss of the hindlimbs in whales. Limb loss and concurrent morphological and physiological changes associated with the transition from land to water are discussed within the context of the current whale phylogeny. Emphasis is placed on fore- and hindlimb development, how the forelimbs transformed into flippers, and how the hindlimbs regressed, leaving either no elements or vestigial skeletal elements. Hindlimbs likely began to regress only after the ancestors of whales entered the aquatic environment: Hindlimb function was co-opted by the undulatory vertical axial locomotion made possible by the newly evolved caudal flukes. Loss of the hindlimbs was associated with elongation of the body during the transition from land to water. Limblessness in most snakes is also associated with adoption of a new (burrowing) lifestyle and was driven by developmental changes associated with elongation of the body. Parallels between adaptation to burrowing or to the aquatic environment reflect structural and functional changes associated with the switch to axial locomotion. Because they are more fully studied and to determine whether hindlimb loss in lineages that are not closely related could result from similar genetically controlled developmental pathways, we discuss developmental (cellular and genetic) processes that may have driven limb loss in snakes and leg-less lizards and compare these processes to the loss of hindlimbs in whales. In neither group does ontogenetic or phylogenetic limb reduction result from failure to initiate limb development. In both groups limb loss results from arrested development at the limb bud stage, as a result of inability to maintain necessary inductive tissue interactions and enhanced cell death over that seen in limbed tetrapods. An evolutionary change in Hox gene expression--as occurs in snakes--or in Hox gene regulation--as occurs in some limbless mutants--is unlikely to have initiated loss of the hindlimbs in cetaceans. Selective pressures acting on a wide range of developmental processes and adult traits other than the limbs are likely to have driven the loss of hindlimbs in whales.     

<Emphasis Type="Italic">Borrelia lusitaniae</Emphasis> OspA Gene Heterogeneity in Mediterranean Basin Area

In this study, Borrelia lusitaniae DNA extracted from ticks and lizards was used to amplify the outer surface protein A (OspA) gene in order to increase knowledge about sequence variability in the Mediterranean basin area, to better understand how Borrelia lusitaniae has evolved and how its distribution has expanded. Phylogenetic trees including Italian and reference sequences showed a clear separation of B. lusitaniae OspA strains in two different major clades. North African isolates form a clade with Portuguese POTIB strains, whereas Italian samples are grouped with German strains and a human Portuguese strain. This subdivision was supported by very high posterior probability values in the trees, by both analysis of molecular variance and selective pressure. These results, based on phylogenetic information contained in the OspA gene sequences, show the presence of two different B. lusitaniae strains circulating in the Mediterranean basin area, suggesting two different evolution paths.     

Molecular evidence and marine snake origins

A molecular phylogeny was used to refute the marine scenario for snake origins. Nuclear gene sequences suggested that snakes are not closely related to living varanid lizards, thus also apparently contradicting proposed relationships between snakes and marine mosasaurs (usually considered to be varanoids). However, mosasaurs share derived similarities with both snakes and living varanids. A reanalysis of the morphological data suggests that, if the relationships between living taxa are constrained to the proposed molecular tree, with fossil forms allowed to insert in their optimal positions within this framework, mosasaurs cluster with snakes rather than with varanids. Combined morphological and molecular analyses also still unite marine lizards with snakes. Thus, the molecular data do not refute the phylogenetic evidence for a marine origin of snakes.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Species richness and composition of amphibians and reptiles in a fragmented forest landscape in northeastern Bolivia

We quantified patterns of species richness and species composition of frogs and reptiles (lizards and snakes) among three habitats (continuous forest, forest islands, and a seasonally flooded savannah) and between forest island size and isolation classes in a floristic transition zone in northeastern Santa Cruz Department, Bolivia. Species richness was similar across macrohabitats, as was faunal composition of forested habitats, although savannah harbored a distinct herpetofauna. On forest islands, richness and composition of forest frogs was largely related to isolation, whereas reptiles were affected by both isolation and habitat. The observation that isolation rather than area was the primary driver of distribution patterns on forest islands stands in contrast to many studies, and may be a function of (1) the greater range in forest island isolation values compared to area or (2) the long history of isolation in this landscape.     

Retinal projections in the Northern Water snake Natrix sipedon sipedon (L.)

Retinal projections were studied experimentally in the Northern water snake using modifications of the Nauta silver impregnation technique. Contralaterally, the retina projects to nucleus geniculatus lateralis pars dorsalis and pars ventralis, nucleus lentiformis mesencephali and nucleus geniculatus pretectalis. A sparse projection was also observed to nucleus ovalis. An additional afferent thalamic projection to nucleus ventrolateralis was found in two cases. The retina projects ipsilaterally to the dorsolateral portion of nucleus geniculatus lateralis pars dorsalis, and sparsely to nucleus lentiformis mesencephali and nucleus geniculatus pretectalis. Nucleus posterodorsalis receives dense bilateral retinal projections. Contralaterally, the retina also projects to the superficial layers of the tectum (layers 8–13 of Ramón) and to nucleus opticus tegmenti. Armstrong's findings that the retinal projections in Natrix are qualittatively similar to those in lizards were confirmed. However there are marked quantitative differences among the various pathways and their corresponding nuclei. These differences are particularly striking in comparing the visual projections to the dorsal thalamus, the retino-tecto-rotundal and the retino-geniculate systems. The first is reduced in volume and the second is markedly increased in volume in comparison with lizards. These data lend support to the theories of Walls that snakes evolved from fossorial lizards and of Underwood that the eyes of these lizards underwent reduction but not complete degeneration. Qualitatively the retinal projections are conservative among lizards and snakes, but a history of reduction of these pathways in ancestral snakes with a selective increase in the retino-geniculate system as a surface niche was reattained is reflected in the anatomy of this ophidian visual system.     

Ventral coloration differentiation and mitochondrial sequences of the Chinese Cobra (<Emphasis Type="Italic">Naja atra</Emphasis>) in Taiwan

Differences in coloration between eastern and western populations of the Chinese cobra (Naja atra) in Taiwan have been noted by snake collectors, snake keepers, and users of Chinese traditional medicine, but have never been verified by scientific research. In this study, we quantified the amount of black pigment on ventral scales, and found prominent differences in ventral coloration of populations across Taiwan; populations in eastern Taiwan have black ventral scales while populations in the west are predominantly white. Previous studies have shown a similar east-west population differentiation in regards to venom components. We supplement these data with mitochondrial control region sequences, which show extremely low nucleotide diversity. Black-ventral and white-ventral snakes share major haplotypes and show no genetic differentiation. Nevertheless, moderate F_st and low N_m values between populations indicate low levels of gene flow. With a morphological fixation earlier than mitochondrial sequences on a relatively short evolutionary time scale, ventral coloration is potentially a result of local adaptation. Based upon the results of this study, along with traditional observations, we strongly recommend treating each of the four populations of the Chinese cobra in Taiwan as distinct ESUs. Reintroducing confiscated snakes from the illegal trade back into the wild needs to be halted to prevent artifical gene flow.