Colonization history and genetic diversity: adaptive potential in early stage invasions.

The introduction of Anolis cristatellus from the multiple species anole community of Puerto Rico in the Greater Antilles to the island of Dominica in the Lesser Antilles, with its solitary endemic anole, provides an example of a very recent, timed, single colonization. We investigate the geographic origin and adaptive potential of the Dominican population using a range of methods including mtDNA phylogeography, nuclear microsatellite variation and multiple paternity studies, as well as heritability estimates, common garden experiments and comparative geographic studies of quantitative scalation traits. Phylogeographic analysis of NADH2 and microsatellite studies suggests that the Dominican population arose from a set of individuals from the central west area of Puerto Rico within their endemic range. The multiple-individual inoculation, together with sperm storage and evidence of multiple paternity indicate genetic variability and suggest the potential for adaptation by natural selection. Estimates of heritability, common garden experiments and broad sense Q(ST)/F(ST) ratios, linked to replicated comparisons along elevational transects go some way to suggesting that the invasive populations may be adapting by natural selection, in parallel with the endemic anole, in the brief period since their introduction.

Weak founder effect signal in a recent introduction of Caribbean Anolis.

Species introductions provide a rare opportunity to study rapid evolutionary and genetic processes in natural systems, often under novel environmental pressures. Few empirical studies have been able to characterize genetic founder effects associated with demographic bottlenecks at the earliest stages of species introductions. This study utilizes prior mitochondrial DNA information which identifies the putative source population for a recently established (c. 7 years between import and sampling) species introduction. We investigated the evidence for a founder effect in a highly successful introduction of a Puerto Rican Anolis species that has established itself on Dominica to the localized exclusion of the native, endemic anole. Five highly polymorphic microsatellite loci were used to explore the partitioning of genetic diversity within and between native source, native nonsource, and introduced populations of Anolis cristatellus. Group comparisons reveal significantly lower allelic richness and expected heterozygosity in introduced populations compared to native populations; however, tests for heterozygosity excess relative to allelic richness failed to provide consistent evidence for a founder effect within introduced populations. Significant levels of within-population genetic variation were present in both native and introduced populations. We suggest that aspects of the reproductive ecology of Anolis (high fecundity, sperm storage and multiple paternity) offer an important mechanism by which genetic variation may be maintained following demographic bottlenecks and founder events in some squamate taxa.

Adaptive radiation in Lesser Antillean lizards: molecular phylogenetics and species recognition in the Lesser Antillean dwarf gecko complex, Sphaerodactylus fantasticus.

The time associated with speciation varies dramatically among lower vertebrates. The nature and timing of divergence is investigated in the fantastic dwarf gecko Sphaerodactylus fantasticus complex, a nominal species that occurs on the central Lesser Antillean island of Guadeloupe and adjacent islands and islets. This is compared to the divergence in the sympatric anole clade from the Anolis bimaculatus group. A molecular phylogenetic analysis of numerous gecko populations from across these islands, based on three mitochondrial DNA genes, reveals several monophyletic groups occupying distinct geographical areas, these being Les Saintes, western Basse Terre plus Dominica, eastern Basse Terre, Grand Terre, and the northern and eastern islands (Montserrat, Marie Galante, Petite Terre, Desirade). Although part of the same nominal species, the molecular divergence within this species complex is extraordinarily high (27% patristic distance between the most divergent lineages) and is compatible with this group occupying the region long before the origin of the younger island arc. Tests show that several quantitative morphological traits are correlated with the phylogeny, but in general the lineages are not uniquely defined by these traits. The dwarf geckos show notably less nominal species-level adaptive radiation than that found in the sympatric southern clade of Anolis bimculatus, although both appear to have occupied the region for a broadly similar period of time. Nevertheless, the dwarf gecko populations on Les Saintes islets are the most morphologically distinct and are recognized as a full species (Sphaerodactylus phyzacinus), as are anoles on Les Saintes (Anolis terraealtae).

Evidence for a Müllerian mimetic radiation in Asian pitvipers.

Müllerian mimicry, in which toxic species gain mutual protection from shared warning signals, is poorly understood in vertebrates, reflecting a paucity of examples. Indirect evidence for mimicry is found if monophyletic species or clades show parallel geographic variation in warning patterns. Here, we evaluate a hypothesis of Müllerian mimicry for the pitvipers in Southeast Asia using a phylogeny derived from DNA sequences from four combined mitochondrial regions. Mantel matrix correlation tests show that conspicuous red colour pattern elements are significantly associated with sympatric and parapatric populations in four genera. To our knowledge, this represents the first evidence of a Müllerian mimetic radiation in vipers. The putative mimetic patterns are rarely found in females. This appears paradoxical in light of the Müllerian prediction of monomorphism, but may be explained by divergent selection pressures on the sexes, which have different behaviours. We suggest that biased predation on active males causes selection for protective warning coloration, whereas crypsis is favoured in relatively sedentary females.

The dynamics of natural selection and vicariance in the Dominican anole: patterns of within-island molecular and morphological divergence.

The larger islands of the Lesser Antilles are ecologically and geologically complex and are inhabited by single, but morphologically variable, Anolis species. Although earlier work has indicated that a large part of the morphological variation in Anolis oculatus from Dominica can be attributed to selection, a history of recurrent volcanic activity over the last few million years suggests that vicariance may have also played a significant role. We report a study of variation in the cytochrome b gene of mitochondrial DNA across the island to address this issue. We uncovered a very high degree of polymorphism, with an overall gene diversity of 0.97 and a nucleotide diversity of 0.04. Sequences, on average, differ by 3.82% and the maximum pairwise divergence (corrected for multiple hits) is 9.29%. Most haplotypes are restricted to single localities (a pattern not changed by increasing the sample size). Phylogenetic analysis revealed the presence of two distinct lineages on the island with strong phylogeographic structure. One of these is geographically restricted to a relatively small part of the central Caribbean coast. Sublineages were also discernible within the other more widely distributed lineage, but resolution within and support for these sublineages was poor. The phylogeographic pattern is not congruent with generalized body shape and scalation, but is significantly correlated with color pattern. Even when correcting for this lineage effect with partial Matrix correspondence tests, the relationship between color pattern and vegetation is reaffirmed, suggesting that although both vicariance and selection have played a role in the morphological differentiation of this species, selection for current environmental conditions has been more important. We discuss the causes of the phylogeographic structure in light of the volcanic history of the island and highlight the exceptional instance of congruence between all morphological character systems and lineage boundaries, which occurs at the transition between the northern and southern Caribbean ecotypes.

Electrophoretic profiles and biological activities: intraspecific variation in the venom of the Malayan pit viper (Calloselasma rhodostoma)

The Malayan pit viper (Calloselasma rhodostoma) is of major clinical significance both as a leading cause of snakebite and as the source of ancrod (Arvin). Although its venom has been extensively studied, the degree to which venom composition varies between individuals is poorly known. We individually analysed the venoms of over 100 C. rhodostoma using isoelectric focusing. In all populations, females produced an intense band that was absent from all males, and significant ontogenetic variation was detected. Principal components analysis of the banding profiles also revealed strong geographic variation, which was significantly congruent with variation in the biological activities of the venom (phosphodiesterase, alkalinephosphoesterase, L-amino acid oxidase, arginine ester hydrolase, 5'-nucleotidase, thrombin-like enzyme, haemorrhagic activity). Studies of captive-bred snakes indicate that the intraspecific variation in venom is genetically inherited rather than environmentally induced. The intraspecific variation in venom composition and biological activity could be of applied importance to snakebite therapy, both in correct diagnosis of the source of envenomation and in the development of a more effective antivenom. Greater attention should be given to the source of C. rhodostoma venom used in research to ensure reproducibility of results.

Molecular clocks and geological dates: cytochrome b of Anolis extremus substantially contradicts dating of Barbados emergence.

Even though molecular clocks vary in rate to some extent, they are widely used and very important in a range of evolutionary studies, not least in interpreting cause and colonization in phylogeography. Evolutionists may use island age and emergence to give the earliest possible date for colonization by a species and hence give the lower limit in a molecular clock calibration. The geology of the Lesser Antilles is well studied and Barbados, although composed of some ancient rocks, is thought to have emerged only about 1 million years ago (Ma). The cytochrome b mitochondrial gene is the most widely used gene in vertebrate phylogeography, and generally evolves at a rate of 1-2% per million years (Myr) for poikilothermic vertebrates. Divergence measured across almost all of this gene in the endemic anole (Anolis extremus) reveals a mean patristic distance of approximately 8.3% between this clade and its sister, together with distinct divergence and phylogeographical structure within Barbados. The divergence time, estimated by a range of procedures using four calibration points, is not in the least compatible with the proposed geological time of emergence of Barbados. Hence, either the molecular clock rate does not apply to the Barbadian anole population, or the geological dating of the emergence of Barbados is erroneous. The compatibility of geological times and molecular divergence of this complex on Martinique, together with relative rates tests comparing the rates on Barbados and Martinique, do not suggest atypical clock rates. The question of whether Barbados emerged much earlier than is currently thought, or whether the molecular clock assumptions are inappropriate, remains open.

GEOGRAPHIC VARIATION WITHIN AN ISLAND: UNIVARIATE AND MULTIVARIATE CONTOURING OF SCALATION, SIZE, AND SHAPE OF THE LIZARD GALLOTIA GALLOTI.

Microgeographic variation of the vegetarian lizard Gallotia galloti within the island of Tenerife is described using univariate analysis, correlation, multiple group principal component analysis, canonical analysis, transects, and contours. The size varies locally in a mosaic pattern while head shape shows a WNW-ESE cline in the south. The scalation (scale and femoral pore counts) has two facets to its geographic variation, both of which are incongruent with the primary variation in the size and shape. The scalation shows categorical variation (stepped cline) between northern and southern populations and also a strong clinal relationship with altitude. The possible causes of this variation are considered and adaptation to current ecological conditions appears to be implicated for at least the altitudinal variation, although the pertinent factors are not obvious. It is apparent that vertebrate populations distributed across small islands do not necessarily offer the advantage of a discrete homogeneous unit for evolutionary studies but can offer the opportunity for studying microgeographic variation.

Naja siamensis, a cryptic species of venomous snake revealed by mtDNA sequencing.

Because of possible variation in venom composition, an understanding of venomous snake systematics is of great importance for the optimization of antivenom treatment of snakebite patients. Intraspecific variation in the morphology of many venomous snakes complicates the definition and identification of some species when allopatric populations are involved. Selectively neutral or near-neutral mtDNA sequences can reveal evolutionary relationships obscured by ecogenetically-caused morphological variation. We use comparative sequencing of the cytochrome oxidase subunit 1 gene to reveal the existence of a widespread, cryptic species of spitting cobra from southeast Asia. This species, Naja siamensis, is widely sympatric with other Asiatic cobra species. This may be of considerable medical significance, and calls for further research into venom composition in Asiatic cobras.

Phylogeny, paraphyly and ecological adaptation of the colour and pattern in the Anolis roquet complex on Martinique.

Martinique is an environmentally heterogeneous island with a complex geological history. It is occupied by a solitary anole, Anolis roquet, showing marked geographical variation in colour and other features. Phylogenetic analysis of a segment (1 kb) of the mitochondrial cytochrome b gene across the Anolis roquet series in the southern Lesser Antilles and at 63 localities of Anolis roquet in Martinique indicate that A. roquet is paraphyletic as A. extremus (Barbados) is nested within the Martinique populations. Moreover, divergent phylogenetic lineages exist within Martinique (max. 10.6% uncorrected pairwise), and these lineages are closely associated with the geological history of this complex island. However, objective quantification of the spectroradiometric analysis of hue by delta analysis, together with analysis of the colour pattern, indicate that they are primarily determined by adaptation to environmental conditions, irrespective of these phylogenetic lineages. There is remarkable convergence in hue and pattern in both extreme xeric (dark chevrons on a dull, generally grey/brown, background), and montane conditions (black reticulation and non-UV white spots on a bright, saturated green background). Moreover, parallel trends occur between Martinique and other Lesser Antillean anoles, which further argues for adaptation (increase in green saturation in montane areas and higher levels of UV on the dewlap of some Atlantic forms). As an exception, there are two specific situations where anoles from different lineages look different. These are (i). in the low-altitude regions of the northwest where the northwestern and central lineages make contact, and (ii). in the far south of the island where the southern and central lineages meet.

Asiatic cobras: systematics and snakebite.

The population affinities of the Asiatic cobras of the genus Naja are investigated, using multivariate analysis of a range of morphological characters. This complex, which was formerly thought to be monospecific, consists of at least eight full species. In some cases, species whose bites require different antivenoms occur sympatrically. The new understanding of the systematics of the Asiatic cobra complex calls for a reappraisal of cobra antivenom use in Asia, and for more research into venom composition.

The usefulness of amplified fragment length polymorphism markers for taxon discrimination across graduated fine evolutionary levels in Caribbean Anolis lizards.

Fine-level taxon discrimination is important in biodiversity assessment and ecogeographical research. Genomic markers are often required for studies on closely related taxa, however, most existing mitochondrial and nuclear markers require prior knowledge of the genome and are impractical for use in small conservation projects. This study describes the application of amplified fragment length polymorphism (AFLP) to discriminate at four progressively finer evolutionary levels of Caribbean Anolis lizards from the central Lesser Antilles. AFLP is shown to be a rapid and effective method for discriminating between species. Separation increases with primer pair number and choice of primer combination appears to be noncritical. Initial population-level results show markedly less discriminatory power. A screening technique for the identification of population informative markers combining principal component and principal coordinate analyses is presented and assessed. Subsequent results show selected conspecific AFLP data to be remarkably congruent with those of mitochondrial DNA, microsatellite and morphological markers. The use of AFLP as a low-cost nuclear marker in species-level taxon discrimination is supported, whereas population level application demands further consideration.

Evidence that ultraviolet markings are associated with patterns of molecular gene flow.

Recent studies have shown UV vision and markings to be important in vertebrates, particularly birds, where behavioral experiments have demonstrated its potential importance in sexual selection. However, there has been no genetic evidence that UV markings determine patterns of evolution among natural populations. Here we report molecular evidence that UV markings are associated with the pattern of gene flow in the Tenerife lizard (Gallotia galloti). This species has vicariance-induced, approximate east--west lineages in Tenerife closely congruent with the primary lineages of the sympatric gecko species. Against expectations, these molecular phylogeographic lineages (representing geological history) and isolation-by-distance do not appear to influence gene flow. Sexually mature males from populations either side of a latitudinal ecotone have different UV markings and gene flow appears to be linked to this difference in UV markings. It may be that these groups with different UV sexual markings mate assortatively, restricting the gene flow between them. This has implications for debate on the relative importance of vicariance and biotopes in influencing biodiversity, with this evidence supporting the latter.

A phylogeny of the trimeresurus group of pit vipers: new evidence from a mitochondrial gene tree.

The Trimeresurus group is an important radiation of over 40 Asian pit viper species. Once considered congeneric, four genera are generally currently recognized (Trimeresurus sensu stricto, Ovophis, Protobothrops, and Tropidolaemus) but relationships within and between these are still unclear. This study, based on mitochondrial cytochrome b sequences, is the first to include a large number of species (21) and demonstrates that the current taxonomy does not adequately represent either the relationships or the genetic diversity present in the complex. Although many deeper nodes are not strongly supported, the following novel conclusions are all well supported: (1) the paraphyly of Trimeresurus sensu stricto, (2) the presence of several divergent clades within Trimeresurus sensu stricto, (3) the paraphyly of some widespread, medically significant, species, (4) the nonmonophyly of Ovophis, and (5) the monophyly of Protobothrops. Mapping of morphological characters onto the mitochondrial tree further supports the four groups proposed for Trimeresurus sensu stricto.

Can microsatellites be used to infer phylogenies? Evidence from population affinities of the Western Canary Island lizard (Gallotia galloti).

Population phylogeographic studies are generally based solely on mtDNA without corroboration, from an independent segregating unit (i.e., nuclear genes), that the mtDNA gene tree represents the organismal phylogeny. This paper attempts to evaluate the utility of microsatellites for this process by use of the Western Canary Island lacertid (Gallotia galloti) as a model. The geological times of island eruptions are known, and well-supported mtDNA phylogenies exist (corroborated as the organismal phylogeny rather than just a gene tree by nuclear random amplified polymorphic DNAs (RAPDs)). The allelic variation in 12 populations from four islands (representing five haplotype lineages) was investigated in five unlinked microsatellite loci. Analysis of molecular variance showed this data to be highly structured. A series of genetic distances among populations was computed based on both the variance in allele frequency (i.e., F(st) related) and the variance in repeat numbers (i.e., R(st) related). The genetic distances based on the former were more highly correlated with the mtDNA genetic distances than those based on the latter. All trees based on both models supported the primary division shown by mtDNA and RAPDs, which is dated at ca. 2.8 to 5.6 mybp (depending on calibration of the mtDNA clock) and which could, under the evolutionary species concept, be regarded separate species. This was achieved despite theoretical problems posed by the use of few loci, suspected bottlenecks, and large population sizes. The finer details were less consistently represented. Nevertheless, this study demonstrates that even a small number of microsatellites can be useful in corroborating the deeper divisions of a population phylogeny.

Historical biogeography of the Western Rattlesnake (Serpentes: viperidae: Crotalus viridis), inferred from mitochondrial DNA sequence information.

We infer the phylogeography of the Western Rattlesnake (Crotalus viridis) using phylogenetic analysis of mitochondrial DNA sequences from 1345 bp of the genes for cytochrome b and NADH dehydrogenase subunit 4. Two main clades are revealed: one includes populations from east and south of the Rocky Mountains (conventionally referred to as Crotalus viridis viridis and C. v. nuntius), and the other consists of populations west of the Rocky Mountains. Within the western clade, a population from southern Arizona (C. v. cerberus) represents the sister taxon to the remaining western populations. The conventional subspecies recognized in this species do not fully correspond to the phylogenetic pattern, and a review of the systematic status of several populations is needed. Our data allow the inferences that small body size evolved twice and that the ability of one population (C. v. concolor) to secrete highly lethal toxins related to Mojave toxin arose within the complex. Our phylogeny should represent the basis for further studies on the causes of geographical variation in this complex.

Ecological diversification in a group of Indomalayan pitvipers (Trimeresurus): convergence in taxonomically important traits has implications for species identification.

We analyse molecular and phenotypic evolution in a group of taxonomically problematic Indomalayan pitvipers, the Trimeresurus sumatranus group. Mitochondrial DNA sequencing provides a well-resolved phylogeny, with each species representing a distinct lineage. Multivariate morphological analysis reveals a high level of phenotypic differentiation, which is congruent between the sexes but does not reflect phylogenetic history. An adaptive explanation for the observed pattern of differentiation is supported by independent contrasts analysis, which shows significant correlations between current ecology and the characters that most account for the variation between taxa, including those that are presently used to identify the species. Reduced precipitation and altitude, and increased temperature, are correlated with higher numbers of scales on the head, body and tail. It is hypothesized that scale number plays an important role in heat and water exchange by influencing the area of exposed of interstitial skin, and that colour pattern variation reflects selection pressures involving camouflage and thermoregulation. Ecological convergence in traits used for classification is found to have important implications for species identification where taxa are distributed over varying environments.

Nuclear and mtDNA phylogenies of the Trimeresurus complex: implications for the gene versus species tree debate.

Phylogenies based on mitochondrial DNA (mtDNA) may represent gene trees that may not be congruent with the equivalent species tree. One solution to this problem is to include additional, independent loci from the nuclear genome. Sequence data from the seventh intron of the beta-fibrinogen gene were generated for 25 specimens of vipers, including 8 nominal species of the Trimeresurus complex of Asian pit vipers. Phylogenetic trees were generated using maximum-parsimony and maximum-likelihood methods. The taxonomic level at which the intron provided significant phylogenetic information was examined and the trees were compared to those produced from previously obtained mtDNA cytochrome b sequences. A variety of different approaches (separate analyses, conditional data combination, and consensus) were used in an attempt to provide a sound organismal phylogeny based on both nuclear and mtDNA data sets. We discuss the implications for the gene tree-species tree debate and its particular relevance to medically important organisms.

Phylogeography and natural selection in the Tenerife gecko Tarentola delalandii: testing historical and adaptive hypotheses.

Combining phylogeographic and matrix correspondence approaches in the analysis of geographical variation provides a fruitful approach to inferring the causes of molecular and morphological evolution within species. Here we present a study on the gecko Tarentola delalandii on the island of Tenerife, Canary Islands, which provides an outstanding model of an exceptionally high degree of phylogeographic differentiation in magnitude and pattern on a small spatial scale. We reconstruct the population history of T. delalandii using phylogeographic information, matrix correspondence tests and estimates of divergence times in conjunction with geological data. It appears that populations differentiated on three precursor islands and secondary contact followed the junction of these islands. The cytochrome b sequence appears to be evolving at least at approximately 1% per million years in this species. Matrix correspondence tests indicate that morphological character systems may reflect ecological selection regimes (colour pattern), history (body dimensions) or both (scalation). The results imply that natural selection can override a historical legacy, but also underline the potential relevance of molecular phylogenetic data for the interpretation of geographical variation in morphology.

Diet and snake venom evolution.

Venom composition within snake species can show considerable geographical variation, an important consideration because bites by conspecific populations may differ in symptomatology and require different treatments. The underlying causes of this phenomenon have never been explained. Here we present evidence that the variation in the venom of the pitviper Calloselasma rhodostoma (Serpentes: Viperidae) is closely associated with its diet. We also evaluated other possible causes of geographic variation in venom using partial Mantel tests and independent contrasts, but rejected both contemporary gene flow (estimated from geographical proximity) and the phylogenetic relationships (assessed by analysis of mitochondrial DNA) among populations as important influences upon venom evolution. As the primary function of viperid venom is to immobilize and digest prey and prey animals vary in their susceptibility to venom, we suggest that geographical variation in venom composition reflects natural selection for feeding on local prey.