Chris R. Feldman

BERGMANN'S RULE IN NONAVIAN REPTILES: TURTLES FOLLOW IT, LIZARDS AND SNAKES REVERSE IT

Abstract Bergmanns rule is currently defined as a within‐species tendency for increasing body size with increasing latitude or decreasing environmental temperature. This well‐known ecogeographic pattern has been considered a general trend for all animals, yet support for Bergmanns rule has only been demonstrated for mammals and birds. Here we evaluate Bergmanns rule in two groups of reptiles: chelonians (turtles) and squamates (lizards and snakes). We perform both nonphylogenetic and phylogenetic analyses and show that chelonians follow Bergmanns rule (19 of 23 species increase in size with latitude; 14 of 15 species decrease in size with temperature), whereas squamates follow the converse to Bergmanns rule (61 of 83 species decrease in size with latitude; 40 of 56 species increase in size with temperature). Size patterns of chelonians are significant using both nonphylogenetic and phylogenetic methods, whereas only the nonphylogenetic analyses are significant for squamates. These trends are consistent among major groups of chelonians and squamates for which data are available. This is the first study to document the converse to Bergmanns rule in any major animal group as well as the first to show Bergmanns rule in a major group of ectotherms. The traditional explanation for Bergmanns rule is that larger endothermic individuals conserve heat better in cooler areas. However, our finding that at least one ectothermic group also follows Bergmanns rule suggests that additional factors may be important. Several alternative processes, such as selection for rapid heat gain in cooler areas, may be responsible for the converse to Bergmanns rule in squamates.

New Chinese turtles: endangered or invalid? A reassessment of two species using mitochondrial DNA, allozyme electrophoresis and known-locality specimens

Over the past 16 years, 13 new species of geoemydid turtles have been described from China. Ten of these new species are based on specimens purchased through the Hong Kong animal trade. Unfortunately, attempts by scientists to discover wild populations of some these newly described species have failed, raising questions about the legitimacy of the type localities and concerns over the validity of the species. Here the phylogenetic and taxonomic validity of two of these species is tested. Mitochondrial DNA haplotypes and allozyme genotypes of specimens matching the descriptions of Mauremys iversoni and Cuora serrata are compared to specimens of established species collected from known localities. The available evidence is consistent with the hypothesis that the specimens represent polyphyletic, intergeneric hybrids. The systematic status of all the new forms of turtles described from pet trade specimens are critical data for conservation efforts, particularly captive breeding.

Comparative phylogeography of woodland reptiles in California: repeated patterns of cladogenesis and population expansion

The ultimate goal of comparative phylogeographical analyses is to infer processes of diversification from contemporary geographical patterns of genetic diversity. When such studies are employed across diverse groups in an array of communities, it may be difficult to discover common evolutionary and ecological processes associated with diversification. In order to identify taxa that have responded in a similar fashion to historical events, we conducted comparative phylogeographical analyses on a phylogenetically and ecologically limited set of taxa. Here, we focus on a group of squamate reptiles (snakes and lizards) that share similar ecological requirements and generally occupy the same communities in the western USA. At a gross level, deep genetic division in Contia tenuis, Diadophis punctatus, Elgaria multicarinata, the Charina bottae complex, and Lampropeltis zonata are often concordant in the Transverse Ranges, the Monterey Bay and Sacramento‐San Joaquin Delta region, and the southern Sierra Nevada in California. Molecular clock estimates suggest that major phyletic breaks within many of these taxa roughly coincide temporally, and may correspond to important geological events. Furthermore, significant congruence between the phylogeographies of E. multicarinata and L. zonata suggests that the succession of vicariance and dispersal events in these species progressed in concert. Such congruence suggests that E. multicarinata and L. zonata have occupied the same communities through time. However, across our entire multi‐taxon data set, the sequence of branching events rarely match between sympatric taxa, indicating the importance of subtle differences in life history features as well as random processes in creating unique genetic patterns. Lastly, coalescent and noncoalescent estimates of population expansion suggest that populations in the more southerly distributed clades of C. tenuis, D. punctatus, E. multicarinata, and L. zonata have been stable, while populations in more northerly clades appear to have recently expanded. This concerted demographic response is consistent with palaeontological data and previous phylogeographical work that suggests that woodland habitat has become more restricted in southern California, but more widespread in the North during Holocene warming. Future phylogeographical work focusing on allied and ecologically associated taxa may add insight into the ecological and evolutionary processes that yield current patterns of genetic diversity.

The evolutionary origins of beneficial alleles during the repeated adaptation of garter snakes to deadly prey

Where do the genetic variants underlying adaptive change come from? Are currently adaptive alleles recruited by selection from standing genetic variation within populations, moved through introgression from other populations, or do they arise as novel mutations? Here, we examine the molecular basis of repeated adaptation to the toxin of deadly prey in 3 species of garter snakes (Thamnophis) to determine whether adaptation has evolved through novel mutations, sieving of existing variation, or transmission of beneficial alleles across species. Functional amino acid substitutions in the skeletal muscle sodium channel (Nav1.4) are largely responsible for the physiological resistance of garter snakes to tetrodotoxin found in their newt (Taricha) prey. Phylogenetic analyses reject the hypotheses that the unique resistance alleles observed in multiple Thamnophis species were present before the split of these lineages, or that alleles were shared among species through occasional hybridization events. Our results demonstrate that adaptive evolution has occurred independently multiple times in garter snakes via the de novo acquisition of beneficial mutations.

Constraint shapes convergence in tetrodotoxin-resistant sodium channels of snakes

Natural selection often produces convergent changes in unrelated lineages, but the degree to which such adaptations occur via predictable genetic paths is unknown. If only a limited subset of possible mutations is fixed in independent lineages, then it is clear that constraint in the production or function of molecular variants is an important determinant of adaptation. We demonstrate remarkably constrained convergence during the evolution of resistance to the lethal poison, tetrodotoxin, in six snake species representing three distinct lineages from around the globe. Resistance-conferring amino acid substitutions in a voltage-gated sodium channel, Nav1.4, are clustered in only two regions of the protein, and a majority of the replacements are confined to the same three positions. The observed changes represent only a small fraction of the experimentally validated mutations known to increase Nav1.4 resistance to tetrodotoxin. These results suggest that constraints resulting from functional tradeoffs between ion channel function and toxin resistance led to predictable patterns of evolutionary convergence at the molecular level. Our data are consistent with theoretical predictions and recent microcosm work that suggest a predictable path is followed during an adaptive walk along a mutational landscape, and that natural selection may be frequently constrained to produce similar genetic outcomes even when operating on independent lineages.

Parallel Arms Races between Garter Snakes and Newts Involving Tetrodotoxin as the Phenotypic Interface of Coevolution

Parallel “arms races” involving the same or similar phenotypic interfaces allow inference about selective forces driving coevolution, as well as the importance of phylogenetic and phenotypic constraints in coevolution. Here, we report the existence of apparent parallel arms races between species pairs of garter snakes and their toxic newt prey that indicate independent evolutionary origins of a key phenotype in the interface. In at least one area of sympatry, the aquatic garter snake, Thamnophis couchii, has evolved elevated resistance to the neurotoxin tetrodotoxin (TTX), present in the newt Taricha torosa. Previous studies have shown that a distantly related garter snake, Thamnophis sirtalis, has coevolved with another newt species that possesses TTX, Taricha granulosa. Patterns of within population variation and phenotypic tradeoffs between TTX resistance and sprint speed suggest that the mechanism of resistance is similar in both species of snake, yet phylogenetic evidence indicates the independent origins of elevated resistance to TTX.

Genetic architecture of a feeding adaptation: garter snake (Thamnophis) resistance to tetrodotoxin bearing prey

Detailing the genetic basis of adaptive variation in natural populations is a first step towards understanding the process of adaptive evolution, yet few ecologically relevant traits have been characterized at the genetic level in wild populations. Traits that mediate coevolutionary interactions between species are ideal for studying adaptation because of the intensity of selection and the well-characterized ecological context. We have previously described the ecological context, evolutionary history and partial genetic basis of tetrodotoxin (TTX) resistance in garter snakes (Thamnophis). Derived mutations in a voltage-gated sodium channel gene (Nav1.4) in three garter snake species are associated with resistance to TTX, the lethal neurotoxin found in their newt prey (Taricha). Here we evaluate the contribution of Nav1.4 alleles to TTX resistance in two of those species from central coastal California. We measured the phenotypes (TTX resistance) and genotypes (Nav1.4 and microsatellites) in a local sample of Thamnophis atratus and Thamnophis sirtalis. Allelic variation in Nav1.4 explains 23 per cent of the variation in TTX resistance in T. atratus while variation in a haphazard sample of the genome (neutral microsatellite markers) shows no association with the phenotype. Similarly, allelic variation in Nav1.4 correlates almost perfectly with TTX resistance in T. sirtalis, but neutral variation does not. These strong correlations suggest that Nav1.4 is a major effect locus. The simple genetic architecture of TTX resistance in garter snakes may significantly impact the dynamics of phenotypic coevolution. Fixation of a few alleles of major effect in some garter snake populations may have led to the evolution of extreme phenotypes and an ‘escape’ from the arms race with newts.

Phylogenetics of the common raven complex (Corvus: Corvidae) and the utility of ND4, COI and intron 7 of the β‐fibrinogen gene in avian molecular systematics

The common raven (Corvus corax) is one of the most widely distributed and recognizable avian species in the world. Recent molecular work, however, described two mitochondrial lineages of the common raven, termed the Holarctic clade and the California clade, and questioned the monophyly of this taxon by placing the Chihuahuan raven (C. cryptoleucus) sister to the California clade. We evaluated this phylogenetic hypothesis with additional sequence data and increased taxon sampling. We used ∼3.7 kb of DNA sequence data from sections of the mitochondrial coding genes COI, cyt b and ND4, a fragment of the non‐coding mitochondrial DNA control region, and the entire intron 7 of the nuclear β‐fibrinogen gene (β‐fibint 7). We combined these DNA sequence data to erect hypotheses of relationships for lineages of the common raven and related taxa. Maximum parsimony, maximum likelihood, and Bayesian methods yield a paraphyletic common raven. These analyses nest the Chihuahuan raven within the common raven, with strong support for a sister relationship between the Chihuahuan raven and the California clade. In addition, the pied crow (C. albus) is also nested within the common raven, and is sister to the Holarctic clade. Our analyses reveal the challenge of determining phylogenetic relationships and species boundaries in this morphologically conservative genus, and suggest that future molecular work with increased taxon sampling will uncover cryptic species and novel evolutionary relationships. Lastly, this survey is one of a growing number of avian phylogenetic studies to employ either β‐fibint 7 or COI, and the first to use ND4. We developed a simple procedure for comparing rates of evolution in molecular markers, and show that in Corvus the nuclear intron β‐fibint 7 is evolving at a considerably slower pace than the mitochondrial markers, while COI is evolving at a slower rate than cyt b, and ND4 approximately the same rate as cyt b. Hence, β‐fibint 7 and other individual nuclear introns may have limited utility in resolving relationships among recently evolved taxa, whereas both COI and ND4 should be useful in a wide range of avian molecular genetic investigations.

Mitochondrial Variation in Sharp-Tailed Snakes (Contia tenuis): Evidence of a Cryptic Species

Abstract We examined genetic variation and structure in mitochondrial DNA sequences of sharp-tailed snakes (Contia tenuis) from California and southern Oregon. Maximum parsimony and maximum likelihood analyses distinguish two mitochondrial lineages: a north coast clade restricted to cool evergreen forest along the Pacific Coast; and an interior/south clade widespread throughout California. The southern limit of the north coast clade is congruent with that of several other vertebrate taxa, a historical pattern consistent with a long-term marine embayment. We interpret additional phylogeographic pattern as resulting from either gene flow or incomplete lineage sorting. Genetic, distributional, ecological, and morphological data suggest that north coast and interior/south mitochondrial lineages of C. tenuis are distinct at the species level.

Parallel Evolution of Tetrodotoxin Resistance in Three Voltage-Gated Sodium Channel Genes in the Garter Snake Thamnophis sirtalis

Members of a gene family expressed in a single species often experience common selection pressures. Consequently, the molecular basis of complex adaptations may be expected to involve parallel evolutionary changes in multiple paralogs. Here, we use bacterial artificial chromosome library scans to investigate the evolution of the voltage-gated sodium channel (Nav) family in the garter snake Thamnophis sirtalis, a predator of highly toxic Taricha newts. Newts possess tetrodotoxin (TTX), which blocks Nav’s, arresting action potentials in nerves and muscle. Some Thamnophis populations have evolved resistance to extremely high levels of TTX. Previous work has identified amino acid sites in the skeletal muscle sodium channel Nav1.4 that confer resistance to TTX and vary across populations. We identify parallel evolution of TTX resistance in two additional Nav paralogs, Nav1.6 and 1.7, which are known to be expressed in the peripheral nervous system and should thus be exposed to ingested TTX. Each paralog contains at least one TTX-resistant substitution identical to a substitution previously identified in Nav1.4. These sites are fixed across populations, suggesting that the resistant peripheral nerves antedate resistant muscle. In contrast, three sodium channels expressed solely in the central nervous system (Nav1.1–1.3) showed no evidence of TTX resistance, consistent with protection from toxins by the blood–brain barrier. We also report the exon–intron structure of six Nav paralogs, the first such analysis for snake genes. Our results demonstrate that the molecular basis of adaptation may be both repeatable across members of a gene family and predictable based on functional considerations.