Phenotypic consequences of incubation environment in the African elapid genus Aspidelaps

Eggs and embryos of the African elapid genus Aspidelaps were tested for phenotypic plasticity in response to incubation temperature and hygration. Fifty-two fertile eggs, representing both A. lubricus and A. scutatus, were utilized over two successive breeding seasons. Interspecific difference in hatchling snout-vent length (SVL) was significant, but the two species were equivalent in initial egg mass and hatchling mass. Hatchling mass was unaffected by temperature and incubation substrate moisture, but was positively correlated to initial egg mass. In A. lubricus, SVL was influenced by an incubation temperature-substrate moisture interaction. In A. scutatus, SVL was unaffected by the experimental variables. The secondary sex ratio among hatchlings was significantly female skewed under mesic incubation conditions, and male skewed under xeric conditions. One explanation for the biased sex ratios is differential mortality at wet and dry conditions. The results indicate the two species of Aspidelaps are uniform with regard to female investment in egg and hatchling biomass, yet they partition resources differently during embryogenesis to yield morphologically distinct neonates. © 1996 Wiley-Liss, Inc.     

Snake-venom resistance as a mammalian trophic adaptation: lessons from didelphid marsupials

Mammals that prey on venomous snakes include several opossums (Didelphidae), at least two hedgehogs (Erinaceidae), several mongooses (Herpestidae), several mustelids, and some skunks (Mephitidae). As a group, these taxa do not share any distinctive morphological traits. Instead, mammalian adaptations for ophiophagy seem to consist only in the ability to resist the toxic effects of snake venom. Molecular mechanisms of venom resistance (as indicated by biochemical research on opossums, mongooses, and hedgehogs) include toxin-neutralizing serum factors and adaptive changes in venom-targeted molecules. Of these, toxin-neutralizing serum factors have received the most research attention to date. All of the toxin-neutralizing serum proteins discovered so far in both opossums and mongooses are human α1B-glycoprotein homologs that inhibit either snake-venom metalloproteinases or phospholipase A(2) myotoxins. By contrast, adaptive changes in venom-targeted molecules have received far less attention. The best-documented examples include amino-acid substitutions in mongoose nicotinic acetylcholine receptor that inhibit binding by α-neurotoxins, and amino-acid substitutions in opossum von Willebrand factor (vWF) that are hypothesized to weaken the bond between vWF and coagulopathic C-type lectins. Although multiple mechanisms of venom resistance are known from some species, the proteomic complexity of most snake venoms suggests that the evolved biochemical defences of ophiophagous mammals are likely to be far more numerous than currently recognized. Whereas most previous research in this field has been motivated by the potential for medical applications, venom resistance in ophiophagous mammals is a complex adaptation that merits attention from comparative biologists. Unfortunately, evolutionary inference is currently limited by ignorance about many relevant facts that can only be provided by future research.     

Thousands of microsatellite loci from the venomous coralsnake Micrurus fulvius and variability of select loci across populations and related species

Studies of population genetics increasingly use next‐generation DNA sequencing to identify microsatellite loci in nonmodel organisms. There are, however, relatively few studies that validate the feasibility of transitioning from marker development to experimental application across populations and species. North American coralsnakes of the Micrurus fulvius species complex occur in the United States and Mexico, and little is known about their population structure and phylogenetic relationships. This absence of information and population genetics markers is particularly concerning because they are highly venomous and have important implications on human health. To alleviate this problem in coralsnakes, we investigated the feasibility of using 454 shotgun sequences for microsatellite marker development. First, a genomic shotgun library from a single individual was sequenced (approximately 7.74 megabases; 26 831 reads) to identify potentially amplifiable microsatellite loci (PALs). We then hierarchically sampled 76 individuals from throughout the geographic distribution of the species complex and examined whether PALs were amplifiable and polymorphic. Approximately half of the loci tested were readily amplifiable from all individuals, and 80% of the loci tested for variation were variable and thus informative as population genetic markers. To evaluate the repetitive landscape characteristics across multiple snakes, we also compared microsatellite content between the coralsnake and two other previously sampled snakes, the venomous copperhead (Agkistrodon contortrix) and Burmese python (Python molurus).     

Mid-Tertiary elapid snakes (Squamata, Colubroidea) from Riversleigh, northern Australia: early steps in a continent-wide adaptive radiation

Vertebral and cranial remains of elapid snakes have been collected from fossil assemblages at Riversleigh, north-west Queensland, Australia; most are Miocene but one may be late Oligocene and another as young as Pliocene. The oldest specimen (probably the oldest elapid yet known anywhere) is a vertebra that can be referred provisionally to the extant taxon Laticauda (Hydrophiinae, sensu Slowinski and Keogh, 2000), implying that the basal divergences among Australasian hydrophiine lineages had occurred by the early Miocene, in contrast to most previous estimates for the age of this geographically isolated adaptive radiation. Associated vertebrae and jaw elements from a Late Miocene deposit are described as Incongruelaps iteratus nov. gen. et sp., which has a unique combination of unusual derived characters otherwise found separately in several extant hydrophiine taxa that are only distantly related. Associated vertebrae from other sites, and two parietals from a possibly Pliocene deposit, suggest the presence of several other taxa distinct from extant forms, but the amount of material (and knowledge of variation in extant taxa) is currently insufficient to diagnose these forms. The Tertiary elapids of Riversleigh thus appear to be relatively diverse taxonomically, but low in abundance and, with one exception, not referable to extant taxa below the level of Hydrophiinae. This implies that the present diversity of hydrophiine elapids (31 recognized terrestrial genera, and approximately 16 marine) represents the result of substantial extinction as well as the “cone of increasing diversity” that could be inferred from phylogenetic studies on extant forms.     

Assembling an arsenal: origin and evolution of the snake venom proteome inferred from phylogenetic analysis of toxin sequences

We analyzed the origin and evolution of snake venom toxin families represented in both viperid and elapid snakes by means of phylogenetic analysis of the amino acid sequences of the toxins and related nonvenom proteins. Out of eight toxin families analyzed, five provided clear evidence of recruitment into the snake venom proteome before the diversification of the advanced snakes (Kunitz-type protease inhibitors, CRISP toxins, galactose-binding lectins, M12B peptidases, nerve growth factor toxins), and one was equivocal (cystatin toxins). In two others (phospholipase A(2) and natriuretic toxins), the nonmonophyly of venom toxins demonstrates that presence of these proteins in elapids and viperids results from independent recruitment events. The ANP/BNP natriuretic toxins are likely to be basal, whereas the CNP/BPP toxins are Viperidae only. Similarly, the lectins were recruited twice. In contrast to the basal recruitment of the galactose-binding lectins, the C-type lectins were shown to be Viperidae only, with the alpha-chains and beta-chains resulting from an early duplication event. These results provide strong additional evidence that venom evolved once, at the base of the advanced snake radiation, rather than multiple times in different lineages, with these toxins also present in the venoms of the "colubrid" snake families. Moreover, they provide a first insight into the composition of the earliest ophidian venoms and point the way toward a research program that could elucidate the functional context of the evolution of the snake venom proteome.     

Population systematics of the snake genus Naja (Reptilia: Serpentes: Elapidae) in Indochina: Multivariate morphometrics and comparative mitochondrial DNA sequencing (cytochrome oxidase I)

We analyze the population systematics of Asiatic cobras in Indochina, China and the Andaman Islands by means of comparative sequencing of the cytochrome oxidase subunit I gene of the mitochondrial DNA molecule and multivariate analysis of morphological characters. Canonical variate analysis and mtDNA sequence information reveal that the cobras of this region comprise four distinct species: Naja atra from China and northern Vietnam, Naja kaouthia from Burma, central Thailand, Cambodia and southern Vietnam, Naja siamensis from Thailand, Cambodia and southern Vietnam, and Naja sagittifera from the Andaman Islands. The subspecies N. kaouthia suphanensis Nutaphand 1986 shows no mtDNA sequence difference from typical N. kaouthia from central Thailand, and multivariate analysis does not reveal differences in general phenotypic profile; the subspecies is therefore synonymised with Naja kaouthia. The cytochrome oxidase subunit I gene, little used in molecular taxonomy, is shown to be well suited for studies at the species level, as it shows taxonomically useful levels of interspecific divergence but low levels of intraspecific variation; this is particularly relevant for studies of rare species, where sample size is a problem. The combination of multivariate morphometrics and molecular systematics can be particularly powerful in resolving systematic problems in such cases.     

Intraspecific variation in the direction and degree of sex-biased dispersal among sea-snake populations

Higher rates of dispersal in one sex than the other are widespread, and often attributed to the genetic advantages of reduced inbreeding. The direction of sex-biased dispersal shows strong phylogenetic conservatism (e.g. males disperse more than females in most mammals, but the reverse is true in most birds). By contrast, our genetic data reveal strong inter-population variation in the relative dispersal rates of two species of sea snakes (Laticauda saintgironsi and L. laticaudata) in the Noumea Lagoon of New Caledonia. Assignment methods using microsatellite data identified parallel variation in sex-specific dispersal in both species: dispersal was female-biased in the north-west of the sampling area (in islands far from the main island), but male-biased in the south-east (in islands closer to the main island). This flexibility may reflect sex differences in diets, with spatial variation in sex-specific resources generating spatial variation in sex-specific dispersal distances.     

Integrating phylogeny,environment and space to explore variation in macroecological traits of Viperidae and Elapidae (Squamata: Serpentes)

We used eigenvector mapping in space and phylogeny to investigate the relationships among space, phylogeny and environment on body size and range size variation across two groups of venomous snakes – Viperidae and Elapidae – from the New World. Data on species geographic range sizes, maximum body sizes and phylogenetic relationships were compiled from the available literature. The distributional data were also used to calculate the latitudinal and longitudinal midpoint and the environmental centroids for each species. The eigenvectors extracted from the pair wise spatial and phylogenetic distance matrices were integrated with environmental variables into a method of variation partitioning where the variation in each trait was quantitatively attributed to ‘pure’ and/or shared effects of phylogeny, environment and space. Our results showed that variation in body size was predominantly determined by phylogeny in both groups of snakes. For Viperidae, we found that pure ‘effects’ of phylogeny were the strongest, indicating that most of the body size evolution that was phylogenetically determined in this group occurred independently of environment and geographical proximity. Regarding range sizes, pure phylogenetic influences were very low in both groups, whereas the largest single fraction of explained variation corresponded to overlapped influences of the three sets of predictors, especially for Elapidae. Along with this, we found evidence that niche conservatism is an important processes underlying variation in body size and range size in both groups of snakes.     

Genetic surfing,not allopatric divergence,explains spatial sorting of mitochondrial haplotypes in venomous coralsnakes

Strong spatial sorting of genetic variation in contiguous populations is often explained by local adaptation or secondary contact following allopatric divergence. A third explanation, spatial sorting by stochastic effects of range expansion, has been considered less often though theoretical models suggest it should be widespread, if ephemeral. In a study designed to delimit species within a clade of venomous coralsnakes, we identified an unusual pattern within the Texas coral snake (Micrurus tener): strong spatial sorting of divergent mitochondrial (mtDNA) lineages over a portion of its range, but weak sorting of these lineages elsewhere. We tested three alternative hypotheses to explain this pattern—local adaptation, secondary contact following allopatric divergence, and range expansion. Collectively, near panmixia of nuclear DNA, the signal of range expansion associated sampling drift, expansion origins in the Gulf Coast of Mexico, and species distribution modeling suggest that the spatial sorting of divergent mtDNA lineages within M. tener has resulted from genetic surfing of standing mtDNA variation—not local adaptation or allopatric divergence. Our findings highlight the potential for the stochastic effects of recent range expansion to mislead estimations of population divergence made from mtDNA, which may be exacerbated in systems with low vagility, ancestral mtDNA polymorphism, and male‐biased dispersal.