Integrative taxonomy provides evidence for the species status of the Ibero‐Maghrebian grass snake Natrix astreptophora

The grass snake (Natrix natrix) is Europe's most widely distributed and, in many regions, most common snake species, with many morphologically defined subspecies. Yet, the taxonomy of grass snakes is relatively little studied and recent work has shown major conflicts between morphologically defined subspecies and phylogeographical differentiation. Using external morphology, osteological characters, and information from 13 microsatellite loci and two mitochondrial markers, we examine differentiation of the subspecies N. n. astreptophora from the North African Maghreb region, the Iberian Peninsula and neighbouring France. According to previous studies, N. n. astreptophora corresponds to a deeply divergent mitochondrial clade and constitutes the sister taxon of all remaining grass snakes. In the French Pyrenees region, there is a contact zone of N. n. astreptophora with another subspecies, N. n. helvetica. Our analyses of microsatellites and mitochondrial DNA reveal that the distribution ranges of the two taxa abut there, but both hybridize only exceptionally. Even though many morphological characters are highly variable and homoplastic in grass snakes, N. n. astreptophora differs consistently from all other grass snakes by its reddish iris coloration and in having significantly fewer ventral scales and another skull morphology. Considering further the virtual absence of gene flow between N. n. astreptophora and N. n. helvetica, and acknowledging the morphological distinctiveness of N. n. astreptophora and its sister group relationship to all remaining subspecies of grass snakes, we conclude that Natrix astreptophora (Seoane, 1884) should be recognized as a distinct species. Further research is needed to explore whether N. astreptophora is polytypic because a single sample of N. astreptophora from Tunisia turned out to be genetically highly distinct from its European conspecifics.     

Mitochondrial phylogeography,contact zones and taxonomy of grass snakes (Natrix natrix,N. megalocephala)

Grass snakes (Natrix natrix) represent one of the most widely distributed snake species of the Palaearctic region, ranging from the North African Maghreb region and the Iberian Peninsula through most of Europe and western Asia eastward to the region of Lake Baikal in Central Asia. Within N. natrix, up to 14 distinct subspecies are regarded as valid. In addition, some authors recognize big‐headed grass snakes from western Transcaucasia as a distinct species, N. megalocephala. Based on phylogenetic analyses of a 1984‐bp‐long alignment of mtDNA sequences (ND4+tRNAs, cyt b) of 410 grass snakes, a nearly range‐wide phylogeography is presented for both species. Within N. natrix, 16 terminal mitochondrial clades were identified, most of which conflict with morphologically defined subspecies. These 16 clades correspond to three more inclusive clades from (i) the Iberian Peninsula plus North Africa, (ii) East Europe and Asia and (iii) West Europe including Corso‐Sardinia, the Apennine Peninsula and Sicily. Hypotheses regarding glacial refugia and postglacial range expansions are presented. Refugia were most likely located in each of the southern European peninsulas, Corso‐Sardinia, North Africa, Anatolia and the neighbouring Near and Middle East, where the greatest extant genetic diversity occurs. Multiple distinct microrefugia are inferred for continental Italy plus Sicily, the Balkan Peninsula, Anatolia and the Near and Middle East. Holocene range expansions led to the colonization of more northerly regions and the formation of secondary contact zones. Western Europe was invaded from a refuge within southern France, while Central Europe was reached by two distinct range expansions from the Balkan Peninsula. In Central Europe, there are two contact zones of three distinct mitochondrial clades, and one of these contact zones was theretofore completely unknown. Another contact zone is hypothesized for Eastern Europe, which was colonized, like north‐western Asia, from the Caucasus region. Further contact zones were identified for southern Italy, the Balkans and Transcaucasia. In agreement with previous studies using morphological characters and allozymes, there is no evidence for the distinctiveness of N. megalocephala. Therefore, N. megalocephala is synonymized with N. natrix.     

Mitochondrial ghost lineages blur phylogeography and taxonomy of Natrix helvetica and N. natrix in Italy and Corsica

Grass snakes are widely distributed across the Western Palearctic. Recent phylogeographic studies provided evidence that three distinct parapatric species exist. Two of these occur in Italy, Natrix helvetica and N. natrix, and a contact zone between both taxa has been suggested for north-eastern Italy. Moreover, previous investigations revealed for the Italian Peninsula a complex phylogeographic structure. Using mtDNA sequences and microsatellite loci, we examined the situation for mainland Italy, Sicily, Sardinia, and Corsica. Our study confirmed the occurrence of N. natrix in north-eastern Italy. Cline analyses revealed limited gene flow between N. helvetica and N. natrix across a narrow hybrid zone. Within N. helvetica, conflicting patterns of mitochondrial and nuclear genomic differentiation were revealed. Three nuclear genomic clusters were found; one of them corresponded to no fewer than five distinct and, in part, deeply divergent and ancient mitochondrial lineages from mainland Italy and Sicily. This cluster was paraphyletic with respect to the two remaining mitochondrial lineages, each of which matched with another nuclear genomic cluster (one from Corsica plus Sardinia and another one from western Europe north of the Alps). This unexpected pattern most likely results from mainly male-mediated gene flow and female philopatry combined with population-density-dependent processes such as ‘high-density blocking’. With respect to taxonomy, we propose to synonymize N. h. lanzai Kramer, 1970 with N. h. sicula (Cuvier, 1829), acknowledging their lacking nuclear genomic differentiation. The studied hybrid zone of N. h. helvetica and N. h. sicula in Italy is wide, with a smooth cline for nuclear markers, supporting their subspecies status. We found no evidence for the distinctiveness of the two subspecies from Corsica (N. h. corsa) and Sardinia (N. h. cetti), suggesting their synonymy, but refrain from taxonomic conclusions because of small sample sizes and the endangered status of the Sardinian taxon.     

Development of pancreatic acini in embryos of the grass snake Natrix natrix (Lepidosauria,Serpentes)

This study report about the differentiation of pancreatic acinar tissue in grass snake, Natrix natrix, embryos using light microscopy, transmission electron microscopy, and immuno-gold labeling. Differentiation of acinar cells in the embryonic pancreas of the grass snake is similar to that of other amniotes. Pancreatic acini occurred for the first time at Stage VIII, which is the midpoint of embryonic development. Two pattern of acinar cell differentiation were observed. The first involved formation of zymogen granules followed by cell migration from ducts. In the second, one zymogen granule was formed at the end of acinar cell differentiation. During embryonic development in the pancreatic acini of N. natrix, five types of zymogen granules were established, which correlated with the degree of their maturation and condensation. Within differentiating acini of the studied species, three types of cells were present: acinar, centroacinar, and endocrine cells. The origin of acinar cells as well as centroacinar cells in the pancreas of the studied species was the pancreatic ducts, which is similar as in other vertebrates. In the differentiating pancreatic acini of N. natrix, intermediate cells were not present. It may be related to the lack of transdifferentiation activity of acinar cells in the studied species. Amylase activity of exocrine pancreas was detected only at the end of embryonic development, which may be related to animal feeding after hatching from external sources that are rich in carbohydrates and presence of digestive enzymes in the egg yolk. Mitotic division of acinar cells was the main mechanism of expansion of acinar tissue during pancreas differentiation in the grass snake embryos.     

Shaped by uneven Pleistocene climate: mitochondrial phylogeographic pattern and population history of white wagtail Motacilla alba (Aves: Passeriformes)

We studied the phylogeography and population history of the white wagtail Motacilla alba, which has a vast breeding range, covering areas with different Pleistocene climatic histories. The mitochondrial NADH dehydrogenase subunit II gene (ND2) and Control Region (CR) were analyzed for 273 individuals from 45 localities. Our data comprised all nine subspecies of white wagtail. Four primary clades were inferred (M, N, SW and SE), with indications of M. grandis being nested within M. alba. The oldest split was between two haplotypes from the endemic Moroccan M. a. subpersonata (clade M) and the others, at 0.63–0.96 Mya; other divergences were at 0.31–0.38 Mya. The entire differentiation falls within the part of the Pleistocene characterized by Milankovitch cycles of large amplitudes and durations. Clade N was distributed across the northern Palearctic; clade SW in southwestern Asia plus the British Isles and was predicted by Ecological niche models (ENMs) to occur also in central and south Europe; and clade SE was distributed in central and east Asia. The deep divergence within M. a. subpersonata may reflect retention of ancestral haplotypes. Regional differences in historical climates have had different impacts on different populations: clade N expanded after the last glacial maximum (LGM), whereas milder Pleistocene climate of east Asia allowed clade SE a longer expansion time (since MIS 5); clade SW expanded over a similarly long time as clade SE, which is untypical for European species. ENMs supported these conclusions in that the northern part of the Eurasian continent was unsuitable during the LGM, whereas southern parts remained suitable. The recent divergences and poor structure in the mitochondrial tree contrasts strongly with the pronounced, well defined phenotypical differentiation, indicating extremely fast plumage divergence.     

Molecular phylogeography of common garter snakes (Thamnophis sirtalis) in western North America: implications for regional historical forces

Complete ND2 and partial ND4 and cytochrome b mitochondrial DNA (mtDNA) sequences were analysed to evaluate the phylogeographic patterns of common garter snakes (Thamnophis sirtalis) in western North America. This species is widely distributed throughout North America, and exhibits extensive phenotypic variation in the westernmost part of its range. The overall phylogeographic pattern based on mtDNA sequences is concordant with results from studies of other species in this region, implicating historical vicariant processes during the Pleistocene and indicating bottleneck effects of recent dispersal into postglacial habitat. Indeed, the topology is statistically consistent with the hypothesis of both southern (Great Basin and California) and northern (Haida Gwaii) refugia. Specifically, we identified genetic breaks among three major clades: Northwest Coastal populations, Intermountain populations, and all California populations. The California clade contained the only other well-supported branching patterns detected; relationships among populations within the two northern clades were indistinguishable. These molecular splits contrast sharply with all prior geographical analyses of phenotypic variation in T. sirtalis in this region. Our results suggest that the extensive phenotypic variation in western T. sirtalis has been shaped more by local evolutionary forces than by shared common ancestry. Consequently, we consider all morphologically based subspecies designations of T. sirtalis in this region invalid because they do not reflect reciprocal monophyly of the mtDNA sequences.     

Mitochondrial phylogeography and subspecies of the wide‐ranging sub‐Saharan leopard tortoise Stigmochelys pardalis (Testudines: Testudinidae) – a case study for the pitfalls of pseudogenes and GenBank sequences

The leopard tortoise (Stigmochelys pardalis) is the most widely distributed sub‐Saharan tortoise species, with a range extending from the Horn of Africa all over eastern Africa to the Republic of South Africa, Namibia and southernmost Angola. Using 1938 bp of mitochondrial DNA (cyt b gene, partial ND4 gene plus adjacent tRNA genes) from a nearly range‐wide sampling, we examined its phylogeographic structure and compared our findings with previously published GenBank sequences. We identified seven major clades that are largely parapatrically distributed. A few records of distinct haplotypes at the same locality or in close proximity could be the result of translocation of tortoises by man. The greatest diversity occurs in the south of the species’ range, with five out of the seven clades. Testing for isolation‐by‐distance suggests that the observed phylogeographic structure is the result of restricted geographical gene flow and not of historical vicariance. This is in sharp contrast to wide‐ranging thermophilic reptiles from the western Palaearctic, whose phylogeographic structure was significantly shaped by Pleistocene range interruptions, but also by earlier dispersal and vicariant events. Most cyt b sequences of S. pardalis from GenBank turned out to be nuclear pseudogenes, or to be of chimerical origin from such pseudogenes and authentic mitochondrial sequences, which argues for caution regarding uncritical usage of GenBank sequences. The recent revalidation of the two subspecies of S. pardalis was based on such a chimerical sequence that was erroneously identified with the subspecies S. p. babcocki. Furthermore, according to our data, the distribution of mitochondrial clades does match neither the traditional subspecies ranges nor the pronounced geographical size variation of leopard tortoises. We conclude that there is no rationale for recognizing subspecies within S. pardalis.     

Molecular systematics of the western rattlesnake, Crotalus viridis (Viperidae), with comments on the utility of the D-loop in phylogenetic studies of snakes.

Crotalus viridis, the western rattlesnake, ranges throughout western North America and has been divided into at least eight subspecies. However, the validity of and relationships among these subspecies and the monophyly of C. viridis as a whole are questionable. We used mitochondrial DNA sequence data from the D-loop region and ND2 gene to examine the relationships among 26 populations of C. viridis and to test the monophyly of this species. These data were analyzed separately and combined using maximum-likelihood and maximum-parsimony. The C. viridis group was monophyletic in all combined analyses, consisting of two strongly divergent clades. We recommend that these clades be recognized as two distinct evolutionary species: C. viridis and C. oreganus. Crotalus viridis should be restricted to the subspecies viridis and nuntius and the remaining subspecies be assigned to the species C. oreganus. Our data do not allow strong evaluation of the validity of the subspecies. We found that the ND2 gene had greater sequence divergences among closely related individuals than the D-loop region, but this relationship reversed at higher levels of divergence. This pattern is apparently due to: (1) ND2 third positions evolving faster than the D-loop but becoming saturated at higher levels of divergence, and (2) the D-loop evolving faster than ND2 second (and possibly first) positions. Our results suggest that the ND2 gene is preferable for examining intraspecific relationships and the D-loop may better resolve relationships between species of snakes. The latter result is contrary to the common perception of the phylogenetic utility of the D-loop. Another unusual result is that the 145 bp spacer region, adjacent to the 5' end of the light strand of the D-loop, provides greater phylogenetic resolution than the 1030 bp D-loop.     

Phylogeography of the Ibero-Maghrebian red-eyed grass snake (<Emphasis Type="Italic">Natrix astreptophora</Emphasis>)

We examined phylogeographic differentiation of the red-eyed grass snake (Natrix astreptophora) using 1984 bp of mtDNA and 13 microsatellite loci from specimens collected across its distribution range in southwestern Europe and northwestern Africa. Based on phylogenetic analyses of mtDNA, European N. astreptophora constituted the sister clade to a weakly supported North African clade comprised of two deeply divergent and well-supported clades, one corresponding to Moroccan snakes and the other to snakes from Algeria and Tunisia. This tripartite differentiation was confirmed by analyses of microsatellite loci. According to a fossil-calibrated molecular clock, European and North African N. astreptophora diverged 5.44 million years ago (mya), and the two Maghrebian clades split 4.64 mya. These dates suggest that the radiation of the three clades was initiated by the environmental changes related to the Messinian Salinity Crisis and the reflooding of the Mediterranean Basin. The differentiation of N. astreptophora, with distinct clades in the Iberian Peninsula and in the western and eastern Maghreb, corresponds to a general phylogeographic paradigm and resembles the differentiation found in another co-distributed Natrix species, the viperine snake (N. maura). Despite both species being good swimmers, the Strait of Gibraltar constitutes a significant biogeographic barrier for them. The discovery that North Africa harbours two endemic lineages of N. astreptophora necessitates more conservation efforts for these imperilled snakes.     

Recent rapid speciation and ecomorph divergence in Indo‐Australian sea snakes

The viviparous sea snakes (Hydrophiinae) are a young radiation of at least 62 species that display spectacular morphological diversity and high levels of local sympatry. To shed light on the mechanisms underlying sea snake diversification, we investigated recent speciation and eco‐morphological differentiation in a clade of four nominal species with overlapping ranges in Southeast Asia and Australia. Analyses of morphology and stomach contents identified the presence of two distinct ecomorphs: a ‘macrocephalic’ ecomorph that reaches >2 m in length, has a large head and feeds on crevice‐dwelling eels and gobies; and a ‘microcephalic’ ecomorph that rarely exceeds 1 m in length, has a small head and narrow fore‐body and hunts snake eels in burrows. Mitochondrial sequences show a lack of reciprocal monophyly between ecomorphs and among putative species. However, individual assignment based on newly developed microsatellites separated co‐distributed specimens into four significantly differentiated clusters corresponding to morphological species designations, indicating limited recent gene flow and progress towards speciation. A coalescent species tree (based on mitochondrial and nuclear sequences) and isolation‐migration model (mitochondrial and microsatellite markers) suggest between one and three transitions between ecomorphs within the last approximately 1.2 million to approximately 840 000 years. In particular, the macrocephalic ‘eastern’ population of Hydrophis cyanocinctus and microcephalic H. melanocephalus appear to have diverged very recently and rapidly, resulting in major phenotypic differences and restriction of gene flow in sympatry. These results highlight the viviparous sea snakes as a promising system for speciation studies in the marine environment.     

Mitochondrial DNA sequences of five squamates: phylogenetic affiliation of snakes.

Complete or nearly complete mitochondrial DNA sequences were determined from four lizards (Western fence lizard, Warren's spinytail lizard, Terrestrial arboreal alligator lizard, and Chinese crocodile lizard) and a snake (Texas blind snake). These genomes had a typical gene organization found in those of most mammals and fishes, except for a translocation of the glutamine tRNA gene in the blind snake and a tandem duplication of the threonine and proline tRNA genes in the spinytail lizard. Although previous work showed the existence of duplicate control regions in mitochondrial DNAs of several snakes, the blind snake did not have this characteristic. Phylogenetic analyses based on different tree-building methods consistently supported that the blind snake and a colubrid snake (akamata) make a sister clade relative to all the lizard taxa from six different families. An alternative hypothesis that snakes evolved from a lineage of varanoids was not favored and nearly statistically rejected by the Kishino-Hasegawa test. It is therefore likely that the apparent similarity of the tongue structure between snakes and varanoids independently evolved and that the duplication of the control region occurred on a snake lineage after divergence of the blind snake.     

Genetic variability and differentiation among the populations of Dice snake,Natrix tessellata (Sepentes,Colubridae) in the Iranian Plateau

The Dice snake, Natrix tessellate, is one of the widely distributed but poorly studied reptiles of the Iranian plateau. Here, we sequenced a partial fragment of Cyt b for 14 samples of Natrix tessellata from Iran. The data were added to the previously published data from GenBank. The phylogenetic analyses indicated that the Iranian populations of N. tessellata are well divided into two major clades. The north, eastern and NE populations of Dice snake in Iran in an association with the Central Asian populations constitute the first clade which is in turn a sister clade to the European populations. The other clade that is placed in the base of the whole Dice snake clade, is formed by the southern and western Iranian populations. These are separated from the first group with high bootstrap and posterior probability values. Position of the western and southern Iranian populations in the phylogenetic tree and its close relation to the N. natrix samples, included in the analyses, imply that the group is most ancient among all populations of Dice snake in the Palearctic region with a long history in the area.     

Length–mass allometry in snakes

Body size and body shape are tightly related to an animal's physiology, ecology and life history, and, as such, play a major role in understanding ecological and evolutionary phenomena. Because organisms have different shapes, only a uniform proxy of size, such as mass, may be suitable for comparisons between taxa. Unfortunately, snake masses are rarely reported in the literature. On the basis of 423 species of snakes in 10 families, we developed clade‐specific equations for the estimation of snake masses from snout–vent lengths and total lengths. We found that snout–vent lengths predict masses better than total lengths. By examining the effects of phylogeny, as well as ecological and life history traits on the relationship between mass and length, we found that viviparous species are heavier than oviparous species, and diurnal species are heavier than nocturnal species. Furthermore, microhabitat preferences profoundly influence body shape: arboreal snakes are lighter than terrestrial snakes, whereas aquatic snakes are heavier than terrestrial snakes of a similar length. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ●● , ●●–●●.     

Ultrastructure of endocrine pancreatic granules during pancreatic differentiation in the grass snake,Natrix natrix L. (Lepidosauria,Serpentes)

We used transmission electron microscopy to study the pancreatic main endocrine cell types in the embryos of the grass snake Natrix natrix L. with focus on the morphology of their secretory granules. The embryonic endocrine part of the pancreas in the grass snake contains four main types of cells (A, B, D, and PP), which is similar to other vertebrates. The B granules contained a moderately electron‐dense crystalline‐like core that was polygonal in shape and an electron‐dense outer zone. The A granules had a spherical electron‐dense eccentrically located core and a moderately electron‐dense outer zone. The D granules were filled with a moderately electron‐dense non‐homogeneous content. The PP granules had a spherical electron‐dense core with an electron translucent outer zone. Within the main types of granules (A, B, D, PP), different morphological subtypes were recognized that indicated their maturity, which may be related to the different content of these granules during the process of maturation. The sequence of pancreatic endocrine cell differentiation in grass snake embryos differs from that in many vertebrates. In the grass snake embryos, the B and D cells differentiated earlier than A and PP cells. The different sequence of endocrine cell differentiation in snakes and other vertebrates has been related to phylogenetic position and nutrition during early developmental stages.     

Squamate phylogeny and the relationships of snakes and mosasauroids

Cladistic analysis of extant and fossil squamates (95 characters, 26 taxa) finds the fossil squamate, Coniasaurus Owen, 1850, to be the sister-group of the Mosasauroidea (mosasaurs and aigialosaurs). This clade is supported in all 18 shortest cladograms (464 steps; CI 0.677; HI 0.772) by nine characters of the dermatocranium, maxilla, and mandible. A Strict Consensus Tree of the 18 shortest trees collapses to a basal polytomy for most major squamate clades including the clade (Coniasaurus, Mosasauroidea). A Majority Rule Consensus Tree shows that, in 12 of 18 shortest cladograms, the clade Coniasaurus- Mosasauroidea is the sister-group to snakes (Scolecophidia (Alethinophidia, Dinilysia); this entire clade, referred to as the Pythonomorpha ([[Scolecophidia [Alethinophidia, Dinilysia]], [Coniasaurus, Mosasauroidea]]) is the sister-group to all other scleroglossans. Pythonomorpha is supported in these 12 cladograms by nine characters related to the lower jaw and cranial kinesis. In 6 of 18 shortest cladograms, snakes are the sister-group to the clade (Amphisbaenia (Dibamidae (Gekkonoidea, Eublepharidae))). None of the cladograms support the hypothesis that coniasaurs and mosasauroids are derived varanoid anguimorphs. Two additional analyses were conducted: (1) manipulation and movement of problematic squamate clades while constraining ‘accepted’ relationships; (2) additional cladistic analyses beginning with extant taxa, and sequentially adding fossil taxa. From Test I, at 467 steps, Pythonomorpha can be the sister-group to the Anguimorpha, Scincomorpha, ‘scinco-gekkonomorpha’ [scincomorphs, gekkotans, and amphibaenids-dibamids]. At 471 steps Pythonomorpha can be placed within Varanoidea. Treating only mosasauroids and coniasaurs as a monophyletic group: 469 steps, mosasauroids and coniasaurs as sister-group to Anguimorpha; 479 steps, mosasauroids and coniasaurs nested within Varanoidea. Test II finds snakes to nest within Anguimorpha in a data set of only Mosasauroidea + Extant Squamates; the sistergroup to snakes + anugimorphs is (Amphisbaenia (Dibarnidae (Gekkonoidea, Eublepharidae))). No one particular taxon is identified as a keystone taxon in this analysis, though it appears truc that fossil taxa significantly alter the structure of squamate phylogenetic trees.     

Hidden support from unpromising data sets strongly unites snakes with anguimorph ‘lizards’

A combined analysis of nuclear, mitochondrial and morphological data robustly resolves snakes as the sister taxon to anguimorph ‘lizards’. Analysed in isolation, nuclear DNA (nDNA) produces a trichotomy between snakes, iguanians and anguimorphs, mitochondrial DNA (mtDNA) is largely uninformative at deeper levels, and morphology tends to nest snakes deep within anguimorphs or with various legless squamate groups. When analysed simultaneously, the nuclear signal is sufficiently strong that mtDNA and morphology are constrained to choose between alternative resolutions of the iguanian–anguimorph–snake trichotomy (generated by the nDNA) – and both support the snake–anguimorph solution. Combined analyses of fast‐evolving or idiosyncratically evolving markers (mtNDA, morphology) with conservative traits (e.g. nuclear genes) might be the best way to resolve ancient, closely spaced divergences. Fast or idiosyncratic markers potentially provide the most information about short, ancient internodes, but can converge on spurious trees if analysed in isolation. However, if constrained to only choosing between plausible trees, such data can contribute unique and valuable phylogenetic signal that resolves such problematic divergences.     

Molecular systematics,hybridization, and phylogeography of the Bufo americanus complex in Eastern North America

We reconstruct phylogenetic relationships among a well-studied group of toads and find relationships that differ greatly from the current taxonomic understanding. We use mitochondrial sequences encoding ND1, tRNA(Leu(UUR)), and part of 16S to infer relationships among members of the Bufo americanus complex. Focusing on the four taxa that historically have been most problematic due to morphological similarity and hybridization in sympatry, we sample 150 individuals from multiple populations across each species' geographic range. Our evidence conflicts with previous taxonomic hypotheses that were based on ability to hybridize, geographic distribution, and call variation. First, sequences from B. fowleri do not comprise the sister clade to sequences of B. woodhousii; therefore the previous classifications of B. fowleri as sister species to, or eastern subspecies of, B. woodhousii are both called into question. Second, sequences from B. americanus are more closely related to those of B. woodhousii than to those of B. terrestris, indicating that similar advertisement call characteristics evolved independently. Third, sequences of B. fowleri are paraphyletic, with sequences of B. terrestris embedded within. Lastly, sequences from B. fowleri cluster into three distinct mitochondrial clades, with some divergences corresponding to greater than 2mya. These clades are somewhat geographically structured, suggesting divergence in allopatry during the Pleistocene. These mitochondrial divergences are not accompanied by known phenotypic differences, however, suggesting either evolutionary stasis in morphology and behavior, cryptic phenotypic evolution, or that hybridization in secondary contact has homogenized phenotypic differences that may have arisen in allopatry.     

Molecular phylogeny and divergence dates for Australasian elapids and sea snakes (hydrophiinae): evidence from seven genes for rapid evolutionary radiations

One of the most prolific radiations of venomous snakes, the Australo-Melanesian Hydrophiinae includes approximately 100 species of Australasian terrestrial elapids plus all approximately 60 species of viviparous sea snakes. Here, we estimate hydrophiine relationships based on a large data set comprising 5800 bp drawn from seven genes (mitochondrial: ND4, cytb, 12S, 16S; nuclear: rag1, cmos, myh). These data were analysed using parsimony, likelihood and Bayesian methods to better resolve hydrophiine phylogeny and provide a timescale for the terrestrial and marine radiations. Among oviparous forms, Cacophis, Furina and Demansia are basal to other Australian elapids (core oxyuranines). The Melanesian Toxicocalamus and Aspidomorphus group with Demansia, indicating multiple dispersal events between New Guinea and Australia. Oxyuranus and Pseudonaja form a robust clade. The small burrowing taxa form two separate clades, one consisting of Vermicella and Neelaps calanotus, and the other including Simoselaps, Brachyurophis and Neelaps bimaculatus. The viviparous terrestrial elapids form three separate groups: Acanthophis, the Rhinoplocephalus group and the Notechis-Hemiaspis group. True sea snakes (Hydrophiini) are robustly united with the Notechis-Hemiaspis group. Many of the retrieved groupings are consistent with previous molecular and morphological analyses, but the polyphyly of the viviparous and burrowing groups, and of Neelaps, are novel results. Bayesian relaxed clock analyses indicate very recent divergences: the approximately 160 species of the core Australian radiation (including sea snakes) arose within the last 10 Myr, with most inter-generic splits dating to between 10 and 6 Ma. The Hydrophis sea snake lineage is an exceptionally rapid radiation, with > 40 species evolving within the last 5 Myr.     

A workflow of massive identification and application of intron markers using snakes as a model

Relative to the commonly used mitochondrial and nuclear protein‐coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single‐copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next‐generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p‐distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.     

Out of Asia: natricine snakes support the Cenozoic Beringian Dispersal Hypothesis

Based on a combination of six mitochondrial gene fragments (12S RNA, cyt b, ND1, ND2, ND4 and CO1) and one nuclear gene (c-mos) from 22 genera we infer phylogenetic relationships among natricine snakes and examine the date and area of origin of these snakes. Our phylogenetic results indicate: (1) the subfamily Natricinae is strongly supported as monophyletic including a majority of extant genera, and a poorly known and previously unassigned species Trachischium monticola, (2) two main clades are inferred within Natricinae, one containing solely taxa from the Old World (OW) and the other comprising taxa from a monophyletic New World (NW) group with a small number of OW relatives. Within the first clade, the genera Xenochrophis and Amphiesma are apparently not monophyletic. Divergence dating and ancestral area estimation indicate that the natricines originated in tropical Asia during the later Eocene or the Oligocene. We recover two major dispersals events out of Asia, the first to Africa in the Oligocene (28 Ma) and the second to the Western Palearctic and the New World at 27 Ma. This date is consistent with the dispersal of numerous other OW groups into the NW.