Cryptic vicariance in Gulf of California fishes parallels vicariant patterns found in Baja California mammals and reptiles

Comparisons across multiple taxa can often clarify the histories of biogeographic regions. In particular, historic barriers to movement should have affected multiple species and, thus, result in a pattern of concordant intraspecific genetic divisions among species. A striking example of such comparative phylogeography is the recent observation that populations of many small mammals and reptiles living on the Baja California peninsula have a large genetic break between northern and southern peninsular populations. In the present study, I demonstrate that five species of near-shore fishes living on the Baja coastline of the Gulf of California share this genetic pattern. The simplest explanation for this concordant genetic division within both terrestrial and marine vertebrates is that the Baja Peninsula was fragmented by a Plio-Pleistocene marine seaway and that this seaway posed a substantial barrier to movement for near-shore fishes. For some fish species, the signal of this vicariance in mtDNA has been eroded by gene flow and is not evident with classic, equilibrium measures of population structure. Yet, significant divisions are apparent in coalescent analyses that jointly estimate divergence with gene flow. The genetic divisions within Gulf of California fishes also coincide with recognized biogeographic regions based on fish community composition and several environmental factors. It is likely that adaptation to regional environments and present-day oceanographic circulation limit gene exchange between biogeographic regions and help maintain evidence of past vicariance.     

Population Graphs: the graph theoretic shape of genetic structure

Patterns of intraspecific genetic variation result from interactions among both historical and contemporary evolutionary processes. Traditionally, population geneticists have used methods such as F-statistics, pairwise isolation by distance models, spatial autocorrelation and coalescent models to analyse this variation and to gain insight about causal evolutionary processes. Here we introduce a novel approach (Population Graphs) that focuses on the analysis of marker-based population genetic data within a graph theoretic framework. This method can be used to estimate traditional population genetic summary statistics, but its primary focus is on characterizing the complex topology resulting from historical and contemporary genetic interactions among populations. We introduce the application of Population Graphs by examining the range-wide population genetic structure of a Sonoran Desert cactus (Lophocereus schottii). With this data set, we evaluate hypotheses regarding historical vicariance, isolation by distance, population-level assignment and the importance of specific populations to species-wide genetic connectivity. We close by discussing the applicability of Population Graphs for addressing a wide range of population genetic and phylogeographical problems.     

Contrasting population genetic patterns and evolutionary histories among sympatric Sonoran Desert cactophilic Drosophila

We studied population genetic differentiation in the sympatric Sonoran Desert cactophilic flies Drosophila pachea, D. mettleri and D. nigrospiracula across their continental and peninsular ranges. These flies show marked differences in ecology and behaviour including dispersal distances and host cactus specialization. Examination of a fragment of the mitochondrial cytochrome oxidase subunit I gene (mtCOI) reveals that the Sea of Cortez has constituted an effective dispersal barrier for D. pachea, leading to significant genetic differentiation between the continental and peninsular ranges of this species. No genetic differentiation was detected, however, within its continental and peninsular ranges. In contrast, our mtCOI-based results for D. mettleri and D. nigrospiracula are consistent with a previous allozyme-based study that showed no significant genetic differentiation between continental and peninsular ranges of these two species. For D. mettleri, we also found that the insular population from Santa Catalina Island, California, is genetically differentiated with respect to continental and peninsular localities. We discuss how differences in the genetic structure patterns of D. pachea, D. mettleri and D. nigrospiracula may correspond to differences in their dispersal abilities, host preferences and behaviour.     

Two waves of diversification in mammals and reptiles of Baja California revealed by hierarchical Bayesian analysis

Many species inhabiting the Peninsular Desert of Baja California demonstrate a phylogeographic break at the mid-peninsula, and previous researchers have attributed this shared pattern to a single vicariant event, a mid-peninsular seaway. However, previous studies have not explicitly considered the inherent stochasticity associated with the gene-tree coalescence for species preceding the time of the putative mid-peninsular divergence. We use a Bayesian analysis of a hierarchical model to test for simultaneous vicariance across co-distributed sister lineages sharing a genealogical break at the mid-peninsula. This Bayesian method is advantageous over traditional phylogenetic interpretations of biogeography because it considers the genetic variance associated with the coalescent and mutational processes, as well as the among-lineage demographic differences that affect gene-tree coalescent patterns. Mitochondrial DNA data from six small mammals and six squamate reptiles do not support the perception of a shared vicariant history among lineages exhibiting a north-south divergence at the mid-peninsula, and instead support two events differentially structuring genetic diversity in this region.     

The panmixia paradigm of eastern Pacific olive ridley turtles revised: consequences for their conservation and evolutionary biology

Previous studies of the olive ridley Lepidochelys olivacea population structure in the tropical eastern Pacific have indicated the existence of a single panmictic population ranging from Costa Rica to Mexico. This information has been used to design specific management measures to conserve primary nesting beaches in Mexico. However, little is known about olive ridleys in the Baja California Peninsula, their northernmost reproductive limit, where recent observations have shown differences in nesting female behaviour and size of hatchlings relative to other continental rookeries. We used mtDNA control region sequences from 137 turtles from five continental and four peninsular nesting sites to determine whether such differences correspond to a genetic distinction of Baja California olive ridleys or to phenotypic plasticity associated with the extreme environmental nesting conditions of this region. We found that genetic diversity in peninsular turtles was significantly lower than in continental nesting colonies. Analysis of molecular variance revealed a significant population structure (Phi ST = 0.048, P = 0.006) with the inclusion of peninsular samples. Our results: (i) suggest that the observed phenotypic variation may be associated with genetic differentiation and reproductive isolation; (ii) support the recent colonization of the eastern Pacific by Lepidochelys; (iii) reveal genetic signatures of historical expansion and colonization events; and (iv) significantly challenge the notion of a single genetic and conservation unit of olive ridleys in the eastern Pacific. We conclude that conservation measures for olive ridleys in Mexico should be revised to grant peninsular beaches special attention.     

A phytogeographical characterization of the vine flora of the Sonoran and Chihuahuan deserts

Aim This study presents a phytogeographical characterization of the vine flora of two lower North American desert regions as a biogeographical framework for further ecological inquiry into desert vines. Location The phytogeography of the vine flora of the Sonoran and Chihuahuan Deserts was c haracterized based on 263 known species. Methods Checklists of the vines of each desert were developed. Represented genera were then grouped into 10 phytogeographical elements based on worldwide distribution patterns. To compare the floristic composition of the desert floras, an index of species similarity was calculated. Results About a third more species of vines occur in the Sonoran Desert than in the Chihuahuan Desert. Based on the analysis, cosmopolitan genera are the only group more numerous in absolute terms in the Chihuahuan Desert than in the Sonoran Desert. Tropical elements are represented in about the same proportion in each desert as the number of species, however, nearly twice as many pantropical and neotropical genera are represented in the Sonoran Desert as in the Chihuahuan Desert. Proportionately, more genera of temperate elements occur in the Chihuahuan Desert than in the Sonoran desert, although the absolute number of genera is slightly higher in the latter. Main conclusions As these deserts are relatively recent ecological formations and as vines evolved in forest ecosystems, the composition of the desert vine floras is the result of the interaction between historical vegetation types, their constituent taxa and climatic and geological history. The main differences in the vining floras of the present‐day Sonoran and Chihuahuan Deserts appear to be the result of greater historical influence in the Sonoran Desert of (1) tropical vegetation types and (2) the emergence of the Gulf of California. The Chihuahuan Desert vine flora seems to be the result of (1) a more pronounced historical temperate vegetation, (2) the lack of an important isolating event, such as the creation of the Baja California peninsula, and (3) a cooler climate with shorter growing seasons.     

Genetic differentiation and demographic history in Drosophila pachea from the Sonoran Desert

Genetic variation at six microsatellite DNA loci and a segment of the mitochondrial cytochrome oxidase subunit I (COI) locus was used to estimate gene flow, population structure, and demographic history in the cactophilic Drosophila pachea from the Sonoran Desert of North America, a species that shows a strict association with its senita host cactus (genus Lophocereus). For microsatellite analyses, thirteen populations of D. pachea were sampled, five in mainland Mexico and the southwestern USA, and eight on the Baja California (Baja) peninsula, covering essentially the entire range of the species. Analysis of molecular variance (AMOVA) of microsatellite data revealed that populations from both the mainland and the Baja peninsula generally showed little structure, although there were a few exceptions, suggesting some local differentiation and restriction of gene flow within both regions. Pairwise comparisons of F(ST) among each of the mainland and Baja populations showed evidence of both panmixia and population subdivision. AMOVA performed on grouped populations from both the mainland and Baja, however, revealed significant partitioning of genetic variation among the two regions, but no partitioning among localities within each region. Bayesian skyline analyses of the COI data set, consisting of four mainland and seven peninsular populations, revealed population expansions dating to the Pleistocene or late Pliocene in D. pachea from both regions, although regional differences were seen in the estimated timing of the expansions and in changes in effective population size over time.     

The genetic structure of a columnar cactus with a disjunct distribution: Stenocereus gummosus in the Sonoran desert

Stenocereus gummosus is a columnar cactus endemic to the Sonoran desert that exhibits a disjunct distribution: it is widely distributed in Baja California and restricted to a small coastal area in mainland Sonora. In this paper, we examine the genetic structure and the mating system of this species in order to explore the origin of the disjunction and describe aspects of the pollination biology. Flowers are nocturnal, pollinated mainly by sphingids and self-incompatible. Polymorphism for allozymes (11 loci) was relatively high (P=75%) but moderate levels of heterozygosity were detected (H(o)=0.103 and H(e)=0.261). Sonoran populations exhibited higher levels of genetic variation than peninsular populations. H(e) declined with latitude when just peninsular and two island populations are included. Substantial levels of inbreeding within populations (f=0.60), moderate differentiation among populations (theta=0.10), and no evidence of isolation by distance were detected. The neighbor-joining phenogram showed Sonoran and island populations nested within peninsular populations. Mainland populations showed greater genetic similarity to island populations, supporting a dispersal hypothesis for the origin of the disjunction. Future studies using DNA markers are suggested in order to better understand the forces that have shaped the genetic structure of this species.     

Collation of data on the ploidy levels and mitochondrial DNA phylogenetic lineages in the silver crucian carp <Emphasis Type="Italic">Carassius auratus gibelio</Emphasis> from Far Eastern and Central Asian populations

The distribution of the diploid and triploid forms and the correspondence between ploidy and mitochondrial DNA (mtDNA) phylogenetic lineages of the silver crucian carp have been studied in Far Eastern water bodies and the Syr Darya River. Both diploid and triploid forms have been found in large river systems (the Amur, Suifun, Tumangan, and Syr Darya river basins). Only the diploid form has been detected in lakes of Bol’shoi Pelis Island (Peter the Great Bay of the Sea of Japan), Sakhalin Island, and the Kamchatka River basin (the Kamchatka Peninsula). It has been confirmed that there are two mtDNA phylogroups in the silver crucian carp in the area studied. Both mtDNA phylogenetic lineages are present in the Suifun and Tumangan river basins. Only one mtDNA phylogroup (characteristic of the gynogenetic form) has been detected in two samples from the Amur River and in the Syr Darya population. The other mtDNA phylogroup is predominant in insular populations and in Kamchatka. The gynogenetic form carries only mtDNA phylogroup I, whereas both phylogroups have been found in diploid bisexual fish. The existence of only two mtDNA phylogroups substantially differing from each other indicates that the gynogenetic form has emerged from the diploid form only once and evolved independently for a long time after that. The absence of haplotypes transitional between the two mtDNA phylogroups suggests that the secondary contact between the gynogenetic and bisexual forms in continental populations occurred within recent historical time. The obtained data confirm that genetic (though asymmetric) exchange between the two forms is possible, which explains the high morphological and, probably, genetic similarity between them.     

Collation of data on the ploidy levels and mitochondrial DNA phylogenetic lineages in the silver crucian carp Carassius auratus gibelio from Far Eastern and Central Asian populations

The distribution of the diploid and triploid forms and the correspondence between ploidy and mitochondrial DNA (mtDNA) phylogenetic lineages of the silver crucian carp have been studied in Far Eastern water bodies and the Syr Darya River. Both diploid and triploid forms have been found in large river systems (the Amur, Suifun, Tumangan, and Syr Darya river basins). Only the diploid form has been detected in lakes of Bol'shoi Pelis Island (Peter the Great Bay of the Sea of Japan), Sakhalin Island, and the Kamchatka River basin (the Kamchatka Peninsula). It has been confirmed that there are two mtDNA phylogroups in the silver crucian carp in the area studied. Both mtDNA phylogenetic lineages are present in the Suifun and Tumangan river basins. Only one mtDNA phylogroup (characteristic of the gynogenetic form) has been detected in two samples from the Amur River and in the Syr Darya population. The other mtDNA phylogroup is predominant in insular populations and in Kamchatka. The gynogenetic form carries only mtDNA phylogroup I, whereas both phylogroups have been found in diploid bisexual fish. The existence of only two mtDNA phylogroups substantially differing from each other indicates that the gynogenetic form has emerged from the diploid form only once and evolved independently for a long time after that. The absence of haplotypes transitional between the two mtDNA phylogroups suggests that the secondary contact between the gynogenetic and bisexual forms in continental populations occurred within recent historical time. The obtained data confirm that genetic (though asymmetric) exchange between the two forms is possible, which explains the high morphological and, probably, genetic similarity between them.     

Deep biogeographical history and cytonuclear discordance in the black-tailed brush lizard (Urosaurus nigricaudus) of Baja California

Molecular tools help us deduce historical events such as vicariance, dispersal, gene flow and speciation. However, our inferences are inevitably linked to the nature of the characters that we use to infer history. For example, the difference in inheritance patterns of mitochondrial DNA (mtDNA) and nuclear DNA (non-recombining maternal vs. recombining biparental inheritance) may lead us to propose different intraspecific histories. The peninsula of Baja California of north-western Mexico, a region affected by a complex geological history involving temporary seaways, permits evaluation of differences between these character types. We sequenced 1966 bp of mtDNA to reconstruct the genealogical history of Urosaurus nigricaudus (black-tailed brush lizard) from samples spanning the entire peninsula. The genealogy revealed several deep divergences, congruent with temporary vicariance events in the mid-peninsular, Loreto and Cape regions, as well as a major split across the Isthmus of La Paz, possibly resulting from a late Miocene seaway. The results support an emerging picture of the historical biogeography of Baja California, which suggests that key vicariance events are older than commonly perceived. The maternal genealogy of U. nigricaudus sharply contrasts with variation in allozymes that suggests very little differentiation across mitochondrial breaks, consistent with a pattern of ongoing gene flow. We interpret this cytonuclear discordance in relation to the historical biogeography of the region.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 89–104.     

Deep genealogies and the mid-peninsular seaway of Baja California

Geological forces and long-term climate changes can have profound effects on species. Such effects may be manifested in the pattern and magnitude of genealogical diversity, as revealed by mitochondrial DNA (mtDNA) lineages. The relative importance of the different forces on a regional biota must be evaluated along with a good understanding of geological and climatological history. The peninsula of Baja California of north-western Mexico is one area where both geology and climate have affected the historical biogeography of the regional biota. Molecular studies based on the genealogical relationships among mtDNA lineages have contributed greatly towards elucidating the historical biogeography of Baja California. Perhaps most noticeably, numerous concordant breaks in mtDNA genealogies half-way along the peninsula suggest a vicariant history in which the mid-peninsula was temporarily submerged. This vicariant explanation has recently been criticized, as no conclusive geological evidence exists for a continuous submergence of the mid-peninsula. As an alternative, a scenario based on climatological factors has been suggested. Here we discuss the validity of the hypothesized mid-peninsular vicariance event and the climate-based alternative in explaining the concordant genealogical breaks. We argue that, despite the significant changes in climate brought about by the glacial cycles throughout the Quaternary, a vicariant history involving a mid-peninsular seaway remains the most parsimonious explanation of the observed patterns in mtDNA genealogies.     

Model-based comparisons of phylogeographic scenarios resolve the intraspecific divergence of cactophilic Drosophila mojavensis

The cactophilic fly Drosophila mojavensis exhibits considerable intraspecific genetic structure across allopatric geographic regions and shows associations with different host cactus species across its range. The divergence between these populations has been studied for more than 60years, yet their exact historical relationships have not been resolved. We analysed sequence data from 15 intronic X-linked loci across populations from Baja California, mainland Sonora-Arizona and Mojave Desert regions under an isolation-with-migration model to assess multiple scenarios of divergence. We also compared the results with a pre-existing sequence data set of eight autosomal loci. We derived a population tree with Baja California placed at its base and link their isolation to Pleistocene climatic oscillations. Our estimates suggest the Baja California population diverged from an ancestral Mojave Desert/mainland Sonora-Arizona group around 230,000-270,000years ago, while the split between the Mojave Desert and mainland Sonora-Arizona populations occurred one glacial cycle later, 117,000-135,000years ago. Although we found these three populations to be effectively allopatric, model ranking could not rule out the possibility of a low level of gene flow between two of them. Finally, the Mojave Desert population showed a small effective population size, consistent with a historical population bottleneck. We show that model-based inference from multiple loci can provide accurate information on the historical relationships of closely related groups allowing us to set into historical context a classic system of incipient ecological speciation.     

Evidence of two Lyssavirus phylogroups with distinct pathogenicity and immunogenicity

The genetic diversity of representative members of the Lyssavirus genus (rabies and rabies-related viruses) was evaluated using the gene encoding the transmembrane glycoprotein involved in the virus-host interaction, immunogenicity, and pathogenicity. Phylogenetic analysis distinguished seven genotypes, which could be divided into two major phylogroups having the highest bootstrap values. Phylogroup I comprises the worldwide genotype 1 (classic Rabies virus), the European bat lyssavirus (EBL) genotypes 5 (EBL1) and 6 (EBL2), the African genotype 4 (Duvenhage virus), and the Australian bat lyssavirus genotype 7. Phylogroup II comprises the divergent African genotypes 2 (Lagos bat virus) and 3 (Mokola virus). We studied immunogenic and pathogenic properties to investigate the biological significance of this phylogenetic grouping. Viruses from phylogroup I (Rabies virus and EBL1) were found to be pathogenic for mice when injected by the intracerebral or the intramuscular route, whereas viruses from phylogroup II (Mokola and Lagos bat viruses) were only pathogenic by the intracerebral route. We showed that the glycoprotein R333 residue essential for virulence was naturally replaced by a D333 in the phylogroup II viruses, likely resulting in their attenuated pathogenicity. Moreover, cross-neutralization distinguished the same phylogroups. Within each phylogroup, the amino acid sequence of the glycoprotein ectodomain was at least 74% identical, and antiglycoprotein virus-neutralizing antibodies displayed cross-neutralization. Between phylogroups, the identity was less than 64.5% and the cross-neutralization was absent, explaining why the classical rabies vaccines (phylogroup I) cannot protect against lyssaviruses from phylogroup II. Our tree-axial analysis divided lyssaviruses into two phylogroups that more closely reflect their biological characteristics than previous serotypes and genotypes.     

Intraspecific phylogeography of Percina evides (Percidae: Etheostomatinae): an additional test of the Central Highlands pre-Pleistocene vicariance hypothesis

North America exhibits the most diverse freshwater fish fauna among temperate regions of the world. Species diversity is concentrated in the Central Highlands, drained by the Mississippi, Gulf Slope and Atlantic Slope river systems. Previous investigations of Central Highlands biogeography have led to conflicting hypotheses involving dispersal and vicariance to explain the diversity and distribution of the freshwater fish fauna. In this investigation predictions of the Central Highlands pre-Pleistocene vicariance hypothesis are tested with a phylogeographic analysis of the percid species Percina evides, which is widely distributed in several disjunct areas of the Central Highlands. Phylogenetic analysis of complete gene sequences of mitochondrially encoded cytochrome b recover three phylogroups, with very low levels of sequence polymorphism within groups. The two western phylogroups are monophyletic with respect to the eastern phylogroup. The recovery of two monophyletic lineages with an eastern and western distribution in the disjunct highland areas is a pattern expected from vicariance, but is not predicted by the Central Highlands pre-Pleistocene vicariance hypothesis. The recovery of very limited mitochondrial DNA polymorphism and lack of phylogeographic structuring across the entire range of the eastern clade, very shallow polymorphism between the disjunct Missouri River and upper Mississippi River populations, and lack of sequence polymorphism in the upper Mississippi River populations, support a hypothesis of dispersal during or following the Pleistocene. The present distribution of P. evides is best explained by both vicariant and dispersal events.     

Characterization of virulence factors and phylogenetic group determination of <Emphasis Type="Italic">Escherichia coli</Emphasis> isolated from diarrheic and non-diarrheic calves from Brazil

This study aimed to detect virulence factors, pathovars, and phylogenetic groups of Escherichia coli strains obtained from feces of calves with and without diarrhea up to 70 days old and to determine the association between occurrence of diarrhea, phylogenetic groups, and pathovars. Phylo-typing analysis of the 336 E. coli strains isolated from calves with Clermont method showed that 21 (6.25 %) belong to phylogroup A, 228 (67.85 %) to phylogroup B1, 2 (0.6 %) to phylogroup B2, 5 (1.49 %) to phylogroup C, 57 (16.96 %) to phylogroup E, and 3 (0.9 %) to phylogroup F. Phylogroup D was not identified and 20 strains (5.95 %) were assigned as “unknown.” The distribution of phylogenetic groups among pathovars showed that NTEC belong to phylogroups B1 (17) and C (4); EPEC to phylogroups B1 (6) and E (8); STEC to phylogroups A (5), B1 (56), B2 (2), C (1), and E (15); EHEC to phylogroups B1 (95) and E (5); and ETEC to phylogroups A (3), B1 (7), and E (10). The EAST-1 strains were phylogroups A (13), B1 (47), E (19), and F (3); E. coli strains of “unknown” phylogroups belonged to pathovars EPEC (1), EHEC (2), STEC (7), and EAST-1 strains (6). ETEC was associated with diarrhea (P = 0.002). Our study did not find association between the phylogenetic background and occurrence of diarrhea (P = 0.164) but did find some relationship in phylogenetic group and pathovar. The study showed that EHEC and STEC are classified as phylogroup B1, EAST-1 phylogroup A, ETEC, and EPEC phylogroup E.     

Use of Museum Specimens to Refine Historical Pronghorn Subspecies Boundaries

Endangered Sonoran (Antilocapra americana sonoriensis) and Peninsular (A. a. peninsularis) pronghorn persist largely because of captive breeding and reintroduction efforts. Recovery team managers want to re-establish pronghorn in their native range, but there is currently uncertainty regarding the subspecies status of extinct pronghorn populations that historically inhabited southern California, USA, and northern Baja California, Mexico. To address this uncertainty, we genotyped museum specimens and conducted phylogenetic and population genetic analyses of historical data in the context of 3 contemporary pronghorn populations. The historical northern Baja California pronghorn share the most ancestry with contemporary Peninsular pronghorn, whereas pronghorn in southern California share more ancestry with contemporary American (A. a. americana) pronghorn. For reintroductions into northern Baja California, the Peninsular subspecies is more appropriate based on museum genetic data. For reintroductions into Southern California, ecological and genetic factors are both important, as the subspecies most genetically related to historical populations (American) may not be well-adapted to the hot, low-elevation deserts of the reintroduction area. © 2019 The Wildlife Society.     

Phylogeographic structure and historical demography of the western diamondback rattlesnake (Crotalus atrox): A perspective on North American desert biogeography

The western diamondback rattlesnake (Crotalus atrox) is a prominent member of North American desert and semi-arid ecosystems, and its importance extends from its impact on the region's ecology and imagery, to its medical relevance as a large deadly venomous snake. We used mtDNA sequences to identify population genetic structure and historical demographic patterns across the range of this species, and relate these to broader patterns of historical biogeography of desert and semi-arid regions of the southwestern USA and adjacent Mexico. We inferred a Late Pliocene divergence between peninsular and continental lineages of Crotalus, followed by an Early Mid Pleistocene divergence across the continental divide within C. atrox. Within desert regions (Sonoran and Chihuahuan Deserts, Southern Plains, and Tamaulipan Plain) we observed population structure indicating isolation of populations in multiple Pleistocene refugia on either side of the continental divide, which we attempt to identify. Evidence of post-glacial population growth and range expansion was inferred, particularly in populations east of the continental divide. We observed clear evidence of (probably recent) gene flow across the continental divide and secondary contact of haplotype lineages. This recent gene flow appears to be particularly strong in the West-to-East direction. Our results also suggest that Crotalus tortugensis (Tortuga Island rattlesnake) and a population of 'C. atrox' inhabiting Santa Cruz Island (in the Gulf of California) previously suggested to be an unnamed species, are in fact deeply phylogenetically nested within continental lineages of C. atrox. Accordingly, we suggest C. tortugensis and 'C. atrox' from Santa Cruz Island be placed in the synonymy of C. atrox.     

Phylogeny, divergence times and species limits of spiny lizards (Sceloporus magister species group) in western North American deserts and Baja California

The broad distribution of the Sceloporus magister species group (squamata: phrynosomatidae) throughout western North America provides an appropriate model for testing biogeographical hypotheses explaining the timing and origins of diversity across mainland deserts and the Baja California Peninsula. We inferred concordant phylogenetic trees describing the higher-level relationships within the magister group using 1.6 kb of mitochondrial DNA (mtDNA) and 1.7 kb of nuclear DNA data. These data provide strong support for the parallel divergence of lineages endemic to the Baja California Peninsula (S. zosteromus and the orcutti complex) in the form of two sequential divergence events at the base of the magister group phylogeny. A relaxed phylogenetic analysis of the mtDNA data using one fossil and one biogeographical constraint provides a chronology of these divergence events and evidence that further diversification within the Baja California clades occurred simultaneously, although patterns of geographical variation and speciation between clades differ. We resolved four major phylogeographical clades within S. magister that (i) provide a novel phylogenetic placement of the Chihuahuan Desert populations sister to the Mojave Desert; (ii) illustrate a mixed history for the Colorado Plateau that includes Mojave and Sonoran Desert components; and (iii) identify an area of overlap between the Mojave and Sonoran Desert clades near Yuma, Arizona. Estimates of bidirectional migration rates among populations of S. magister using four nuclear loci support strong asymmetries in gene flow among the major mtDNA clades. Based on the nonexclusivity of mtDNA haplotypes, nuclear gene flow among populations and wide zones of phenotypic intergradation, S. magister appears to represent a single geographically variable and widespread species.     

A peninsula as an island: multiple forms of evidence for overwater colonization of Baja California by the gartersnake Thamnophis validus

The recent shift toward dispersal rather than vicariant explanations of disjunct distributions has been driven by the use of molecular data to estimate divergence dates between lineages. However, other kinds of evidence can also be critical in evaluating such biogeographic hypotheses. In the present study, we used a multifaceted approach employing diverse analyses of mitochondrial DNA sequences to assess explanations for the disjunct distribution of the gartersnake Thamnophis validus. The occurrence of this species in the Cape Region of the Baja California peninsula, isolated from mainland populations that occur along the west coast of Mexico, might be explained by: (1) separation of the peninsula from mainland Mexico through rifting 4–8 Mya (tectonic vicariance); (2) fragmentation of the range of this semi‐aquatic species because of post‐Pleistocene aridification (vicariance by aridification); (3) natural overwater dispersal across the Gulf of California; or (4) human introduction. Divergence dating indicates that peninsular and mainland T. validus separated from each other within the last 0.5 Myr, thus rejecting tectonic vicariance. In addition, the estimated closest mainland relatives of peninsular snakes are found farther north than expected under this hypothesis. Three findings argue against vicariance by aridification: (1) peninsular snakes and their closest mainland relatives are more genetically similar than predicted; (2) the location of closest mainland relatives is farther south than predicted; and (3) the species is absent from areas where one might expect to find relict populations. Taken together, refutation of the vicariance hypotheses and the fact that the estimated closest mainland relatives are found almost directly across the Gulf from the Cape Region supports some form of overwater colonization. Various additional arguments suggest that natural dispersal is more likely than human introduction. The present study emphasizes the need for multiple kinds of evidence, beyond divergence dates, to discriminate among hypotheses and to provide independent sources of corroboration or refutation in historical biogeography. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 409–424.