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.     

Microsatellite variation among diverging populations of Drosophila mojavensis

Divergence and speciation may occur by various means, depending on the particular history, selective environments, and genetic composition of populations. In Drosophila mojavensis, a good model of incipient speciation, understanding the population genetic structure within this group facilitates our ability to understand the context in which reproductive isolation among populations is developing. Here we report the genetic structure and relationships of D. mojavensis populations at nuclear loci. We surveyed 29 populations throughout the distribution of D. mojavensis for four microsatellite loci to differentiation among populations of this species. These loci reveal four distinct geographical regions of D. mojavensis populations in the south-western United States and north-western Mexico--(i) Baja California peninsula (Baja), (ii) Sonora, Mexico-southern Arizona, United States (Sonora), (iii) Mojave Desert and Grand Canyon (Mojave), and (iv) Santa Catalina Island (Catalina). While all regions show strong isolation, Mojave and Catalina are highly diverged from other regions. Within any region, populations are largely homogenous over broad geographical distances. Based on the population structure, we find clear geographical barriers to gene flow appear to have a strong effect in isolating populations across regions for this species.     

Evolutionary relationships of Drosophila mojavensis geographic host races and their sister species Drosophila arizonae

The cactophilic Drosophila mojavensis species group living in the deserts and dry tropical forests of the southwestern United States and Mexico provides a valuable system for studies in diversification and speciation. Rigorous studies of the relationships between host races of D. mojavensis and the relationships among the members of the species group (D. mojavensis, Drosophila arizona, and Drosophila navojoa) are lacking. We used mitochondrial CO1 sequence data to address the phylogenetics and population genetics of this species group. In this study we have found that the sister species D. mojavensis and D. arizonae share no mitochondrial haplotypes and thus show no evidence for recent introgression. We estimate the divergence time between D. mojavensis and D. arizonae to be between 1.91 and 2.97 million years ago. D. arizonae shows little structure in our population genetic analyses but there is phylogenetic differentiation between southeastern and northern populations of D. arizonae. Drosophila mojavensis shows significant population and phylogenetic structure across the four geographic regions of its distribution. The mitochondrial data support an origin of D. mojavensis on the mainland with early differentiation into the populations now found in the Mojave Desert and the Mainland Sonoran Desert and later colonization of the Baja Peninsula, in contrast to previous models. Also, the sister clade to D. mojavensis/D. arizonae includes D. navojoa and Drosophila huaylasi. By defining the genetic relationships among these populations, we provide a foundation for more sophisticated hypothesis testing regarding the timing of early speciation events and host switches in this species group.     

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.     

As the raven flies: using genetic data to infer the history of invasive common raven (Corvus corax) populations in the Mojave Desert

Common raven (Corvus corax) populations in Mojave Desert regions of southern California and Nevada have increased dramatically over the past five decades. This growth has been attributed to increased human development in the region, as ravens have a commensal relationship with humans and feed extensively at landfills and on road-killed wildlife. Ravens, as a partially subsidized predator, also represent a problem for native desert wildlife, in particular threatened desert tortoises (Gopherus agassizii). However, it is unclear whether the more than 15-fold population increase is due to in situ population growth or to immigration from adjacent regions where ravens have been historically common. Ravens were sampled for genetic analysis at several local sites within five major areas: the West Mojave Desert (California), East Mojave Desert (southern Nevada), southern coastal California, northern coastal California (Bay Area), and northern Nevada (Great Basin). Analyses of mtDNA control region sequences reveal an increased frequency of raven 'Holarctic clade' haplotypes from south to north inland, with 'California clade' haplotypes nearly fixed in the California populations. There was significant structuring among regions for mtDNA, with high F(ST) values among sampling regions, especially between the Nevada and California samples. Analyses of eight microsatellite loci reveal a mostly similar pattern of regional population structure, with considerably smaller, but mostly significant, values. The greater mtDNA divergences may be due to lower female dispersal relative to males, lower N(e), or effects of high mutation rates on maximal values of F(ST). Analyses indicate recent population growth in the West Mojave Desert and a bottleneck in the northern California populations. While we cannot rule out in situ population growth as a factor, patterns of movement inferred from our data suggest that the increase in raven populations in the West Mojave Desert resulted from movements from southern California and the Central Valley. Ravens in the East Mojave Desert are more similar to ones from northern Nevada, indicating movement between those regions. If this interpretation of high gene flow into the Mojave Desert is correct, then efforts to manage raven numbers by local control may not be optimally effective.     

Tempo and mode of speciation in the Baja California disjunct fish species Anisotremus davidsonii

The Baja California region provides a natural setting for studying the early mechanisms of allopatric speciation in marine systems. Disjunct fish populations from several species that occur in the northern Gulf of California and northern Pacific coast of Baja California, but are absent from its southern shores, were previously shown to be genetically isolated, making them excellent candidates for studying allopatry. In addition, one of these species, the sargo Anisotremus davidsonii, has two pairs of congeneric Panamic trans-isthmian geminate species that allow for internal molecular clock calibration. Phylogeographic and demographic approaches based on mitochondrial (cytochrome b) and nuclear (S7 ribosomal protein) sequences showed that A. davidsonii entered the gulf from the south, and later colonized the Pacific coast, approximately 0.6-0.16 million years ago. Pacific coast colonization may have used a route either around the southern cape of Baja California or across the peninsula through a natural seaway. However, while several seaways have been described from different geological times, none matches the dates of population disjunction, yet much geological work remains to be done in that area. At the present time, there is no evidence for dispersal around the southern tip of the Baja California Peninsula. Signatures of incipient allopatric speciation were observed, such as the reciprocal monophyly of disjunct populations for the mitochondrial marker. However, other characteristics were lacking, such as a strong difference in divergence and coalescence times. Taken together, these results suggest that disjunct populations of A. davidsonii may be consistent with the earliest stages of allopatric speciation.     

Colour variation is incongruent with mitochondrial lineages: cryptic speciation and subsequent diversification in a Gulf of California reef fish (Teleostei: Blennioidei)

The Gulf of California endemic reef fish, Acanthemblemaria crockeri (Blennioidei, Chaenopsidae), reportedly has two colour morphs, one with melanic lateral spots ('Gulf' morph) and one with orange spots ('Cape' morph). In this study, we recorded colour morph in both males and females and collected mitochondrial DNA sequence data for cytochrome c oxidase I (COI) and tRNA-Pro/D-loop of specimens from throughout the Gulf to explore the genetic basis of the colour morphs. Two highly divergent (HKY + I distance = 11.9% for COI), reciprocally monophyletic lineages were identified, consistent with the presence of two parapatric species. A 30-km gap between the distributions of mitochondrial lineages roughly corresponds to a hypothesized former seaway across the Baja California peninsula north of La Paz, although the estimated divergence time (1.84 million years ago) is more recent than the hypothetical seaway (3–4 million years ago). Surprisingly, the distribution of mitochondrial species is not congruent with the distribution of either male or female colour morphs. Our analysis also revealed significant population differentiation within both species and no shared haplotypes among populations. The northern Gulf species includes four populations (NB, CB, NM and CM) corresponding to northern and central Baja and northern and central mainland sites, while the Cape species includes two populations (SB and SM) corresponding to the Baja and mainland sides of the southern Gulf. The NB/CB division corresponds to a hypothesized Plio–Pleistocene mid-peninsular seaway. The level of genetic divergence documented in this lineage is extraordinary for a marine fish with a pelagic larval stage within a semi-enclosed basin.     

Genealogical concordance and the specific status of Peromyscus sejugis

Peromyscus sejugis, a peripheral isolate of Peromyscus maniculatus, is a threatened taxon endemic to 2 small islands in the Sea of Cortés. Although its insularity makes the specific recognition of P. sejugis inherently problematic, resolution of this problem has important conservation implications. To evaluate the specific validity and evolutionary history of P. sejugis, we compared sequence variation (ND3/ND4L/ND4) in mtDNA for both island populations of P. sejugis with that for 8 populations of P. maniculatus from mainland Baja California. Each island population of P. sejugis had a single haplotype (0.7% sequence divergence), whereas 11 different haplotypes (mean sequence divergence = 0.68%) were obtained for the populations of P. maniculatus. The mean sequence divergence between the populations of the 2 species was 2.0%. Nested clade analysis supports the conclusion that P. sejugis is an insular isolate of P. maniculatus from mainland Baja California. Although our analysis confirms a low level of mtDNA divergence between P. sejugis and P. maniculatus from Baja California, the genealogical concordance of morphological, chromosomal, microsatellite, and mtDNA haplotype distinctiveness supports the conclusion that the 2 island populations of P. sejugis constitute independent evolutionarily significant units and together represent a phylogenetic species distinct from the P. maniculatus from Baja California.     

Historical connectivity,contemporary isolation and local adaptation in a widespread but discontinuously distributed species endemic to Taiwan,Rhododendron oldhamii (Ericaceae)

Elucidation of the evolutionary processes that constrain or facilitate adaptive divergence is a central goal in evolutionary biology, especially in non-model organisms. We tested whether changes in dynamics of gene flow (historical vs contemporary) caused population isolation and examined local adaptation in response to environmental selective forces in fragmented Rhododendron oldhamii populations. Variation in 26 expressed sequence tag-simple sequence repeat loci from 18 populations in Taiwan was investigated by examining patterns of genetic diversity, inbreeding, geographic structure, recent bottlenecks, and historical and contemporary gene flow. Selection associated with environmental variables was also examined. Bayesian clustering analysis revealed four regional population groups of north, central, south and southeast with significant genetic differentiation. Historical bottlenecks beginning 9168–13,092 years ago and ending 1584–3504 years ago were revealed by estimates using approximate Bayesian computation for all four regional samples analyzed. Recent migration within and across geographic regions was limited. However, major dispersal sources were found within geographic regions. Altitudinal clines of allelic frequencies of environmentally associated positively selected outliers were found, indicating adaptive divergence. Our results point to a transition from historical population connectivity toward contemporary population isolation and divergence on a regional scale. Spatial and temporal dispersal differences may have resulted in regional population divergence and local adaptation associated with environmental variables, which may have played roles as selective forces at a regional scale.     

Population genetics of <Emphasis Type="Italic">Penstemon albomarginatus</Emphasis> (Plantaginaceae), a rare Mojave Desert species of conservation concern

Penstemon albomarginatus is a psammophytic endemic of the Mojave Desert, found only in deep sand and dune habitats of San Bernardino County, California, Mohave County, Arizona, and Clark and Nye Counties in Nevada. We used six microsatellite loci to assess genetic differentiation and diversity for 228 individuals across the 12 known populations of this rare species. A slight heterozygote deficiency was found in two populations, but most populations show no signs of inbreeding. Results show a geographic pattern of northern populations being more closely related to one another compared to all other geographic regions. Genetic diversity was greatest in the southern populations, with decreasing amounts of diversity observed with latitude. In general, the geographic pattern of genetic diversity among all populations suggests a post-glacial dispersal from south-to-north. Our results are discussed in the framework of anthropogenic pressures on deep sand habitats of the Mojave Desert.     

Analysis of microsatellite variation in Pinus radiata reveals effects of genetic drift but no recent bottlenecks

Most conifer species occur in large continuous populations, but radiata pine, Pinus radiata, occurs only in five disjunctive natural populations in California and Mexico. The Mexican island populations were presumably colonized from the mainland millions of years ago. According to Axelrod (1981), the mainland populations are relicts of an earlier much wider distribution, reduced some 8,000 years ago, whereas according to Millar (1997, 2000), the patchy metapopulation-like structure is typical of the long-term population demography of the species. We used 19 highly polymorphic microsatellite loci to describe population structure and to search for signs of the dynamics of population demography over space and time. Frequencies of null alleles at microsatellite loci were estimated using an approach based on the probability of identity by descent. Microsatellite genetic diversities were high in all populations [expected heterozygosity (H(e)) = 0.68-0.77], but the island populations had significantly lower estimates. Variation between loci in genetic differentiation (F(ST)) was high, but no locus deviated statistically significantly from the rest at an experiment wide level of 0.05. Thus, all loci were included in subsequent analysis. The average differentiation was measured as F(ST) = 0.14 (SD 0.012), comparable with earlier allozyme results. The island populations were more diverged from the other populations and from an inferred common ancestral gene pool than the mainland ones. All populations showed a deficiency of expected heterozygosity given the number of alleles, the mainland populations more so than the island ones. The results thus do not support a recent important contraction in the mainland range of radiata pine.     

Speciation and demographic history of the Cortez bonefish, Albula sp. A (Albuliformes: Albulidae), in the Gulf of California inferred from mitochondrial DNA

Nucleotide sequence data from a segment of the mitochondrial cytochrome b (Cyt b ) gene were used to infer demographic history and examine conditions that may have led to speciation in the Cortez bonefish ( Albula sp. A) in the Gulf of California, Mexico, a currently undescribed species of bonefish in the Albula vulpes complex. Analysis of molecular variance in 39 individuals collected from three localities along the eastern gulf coast, over c. 850 km, revealed a lack of population structure among localities (overall Φ _ST=−0·015), with 100% of the genetic variation distributed within populations. Analysis of combined sequences from these individuals using neutrality tests and the mismatch distribution provided evidence of a population expansion dating to the Pleistocene. The population expansion was supported by maximum likelihood estimates of changes in long-term female effective population size ( N _ef). A molecular clock estimate of divergence, provisionally calibrated for the Cyt b gene in Albula , indicates that Albula sp. A and its sister species in the eastern Pacific, Albula esuncula , diverged from a common ancestor c. 5·0–8·8 million years ago. This date is about the time the Baja California peninsula separated from mainland Mexico during the formation of the Gulf of California. Oceanographic and ecological changes associated with the opening of the gulf likely provided conditions favourable for adaptive radiation and reproductive isolation, ultimately resulting in the allopatric formation of two separate lineages. The co-occurrence of Albula sp. A and A. esuncula found today in the coastal Pacific waters of northern Mexico is most probably the result of secondary contact after speciation.     

Phylogeography and systematics of the Peromyscus eremicus species group and the historical biogeography of North American warm regional deserts

Phylogeographic relationships among 26 populations from throughout the geographic range of the Peromyscus eremicus species group are described based on sequence data for a 699-bp fragment of the mitochondrial DNA COIII gene. Distance, maximum-likelihood, and maximum-parsimony analyses of phylogenetic trees generated under four separate character-weighting strategies and representing five alternative biogeographic hypotheses revealed the existence of a cryptic species (Peromyscus fraterculus, previously included under P. eremicus) on the Baja California Peninsula and adjacent southwestern California and two distinct forms of P. eremicus, one from the Mojave, Sonoran, and northwestern Chihuahuan regional deserts (West) and one from the remainder of the Chihuahuan Desert (East). Distinctiveness of P. fraterculus is supported by previous morphometric and allozyme analyses, including comparisons with neighboring P. eremicus and parapatric P. eva, with which P. fraterculus shares a sister taxon relationship. Divergence of the eva + fraterculus, West + East eremicus, and P. merriami haplotype lineages likely occurred in the late Neogene (3 Ma), in response to northern extension of the Sea of Cortéz and elevation of the Sierra Madre Occidental; divergence of eva from fraterculus is concordant with the existence of a trans-Peninsular seaway during the Pleistocene (1 Ma); and divergence of West from East eremicus occurred during the Pleistocene pluvial-interpluvial cycles, but well before the Wisconsinan glacial interval. The sequence of divergence within the eremicus species group and causal association of geological events of the Neogene and Holocene provide a working hypothesis against which phylogeographic patterns among other arid-adapted species of the warm regional deserts of North America may be compared.     

Phylogeography of the ground squirrel subgenus Xerospermophilus and assembly of the Mojave Desert biota

Aim The Mohave ground squirrel (Xerospermophilus mohavensis) is one of a few endemic species of the Mojave Desert of south‐western North America. We describe phylogeographic patterns within this species and its sister taxon (Xerospermophilus tereticaudus) and test hypotheses concerning their biogeographical history using genetic signatures of stable versus expanding populations. We compare these patterns with those of other Mojave species to evaluate the role of vicariance in producing phylogeographic structure during the assembly of the Mojave Desert biota. Location The Mojave Desert and adjacent desert regions of south‐western North America. Methods Complete cytochrome b gene sequences of X. mohavensis (46 individuals representing 11 localities) and X. tereticaudus (38 individuals representing 14 localities) were analysed using Bayesian methods to infer phylogenetic relationships. Genetic signals of stable or expanding populations were examined based on the distribution of recent mutations and pairwise differences, as well as with a coalescent‐based approach. Results The two species are reciprocally monophyletic and may have diverged in response to the late Pliocene–early Pleistocene uplift of the Transverse Ranges and Mojave block. Little phylogeographic structure is evident within X. mohavensis, but there is a signature of northern expansion from a presumably full‐pluvial refugium in the Mojave River basin. Four geographic subgroups are evident within X. tereticaudus, and there is a signature of northern expansion from a presumably full‐pluvial refugium in the Sonoran coastal plains. Roughly congruent phylogeographic patterns are found within five arid‐adapted taxa, indicating a strong element of vicariance during the assembly of the generally transitional Mojave Desert biota. Main conclusions We present a preliminary model for the historical assembly of the Mojave Desert biota that indicates a strong vicariant element producing autochthonous lineages (including X. mohavensis) that diverged during the major geological and climatic events of the last 5 Myr. Phylogeographic partitioning within the Mojave Desert underscores the necessity of immediate conservation measures for this unique and fragile arid ecosystem that is locked between two large metropolitan population centres and is the target of continued adverse environmental impact.     

GENETICS OF INCIPIENT SPECIATION IN DROSOPHILA MOJAVENSIS. III. LIFE‐HISTORY DIVERGENCE IN ALLOPATRY AND REPRODUCTIVE ISOLATION

We carried out a three‐tiered genetic analysis of egg‐to‐adult development time and viability in ancestral and derived populations of cactophilic Drosophila mojavensis to test the hypothesis that evolution of these life‐history characters has shaped premating reproductive isolation in this species. First, a common garden experiment with 11 populations from Baja California and mainland Mexico and Arizona reared on two host species revealed significant host plant X region and population interactions for viability and development time, evidence for host plant adaptation. Second, replicated line crosses with flies reared on both hosts revealed autosomal, X chromosome, cytoplasmic, and autosome X cactus influences on development time. Viability differences were influenced by host plants, autosomal dominance, and X chromosomal effects. Many of the F₁, F₂, and backcross generations showed evidence of heterosis for viability. Third, a QTL analysis of male courtship song and epicuticular hydrocarbon variation based on 1688 Baja × mainland F₂ males also revealed eight QTL influencing development time differences. Mainland alleles at six of these loci were associated with longer development times, consistent with population‐level differences. Eight G × E interactions were also detected caused by longer development times of mainland alleles expressed on a mainland host with smaller differences among Baja genotypes reared on the Baja host plant. Four QTL influenced both development time and epicuticular hydrocarbon differences associated with courtship success, and there was a significant QTL‐based correlation between development time and cuticular hydrocarbon variation. Thus, the regional shifts in life histories that evolved once D. mojavensis invaded mainland Mexico from Baja California by shifting host plants were genetically correlated with variation in cuticular hydrocarbon‐based mate preferences.     

Historical Divergence and Gene Flow in the Genus Zea

Gene flow plays a fundamental role in plant evolutionary history, yet its role in population divergence—and ultimately speciation—remains poorly understood. We investigated gene flow and the modalities of divergence in the domesticate Zea mays ssp. mays and three wild Zea taxa using sequence polymorphism data from 26 nuclear loci. We described diversity across loci and assessed evidence for adaptive and purifying selection at nonsynonymous sites. For each of three divergence events in the history of these taxa, we used approximate Bayesian simulation to estimate population sizes and divergence times and explicitly compare among alternative models of divergence. Our estimates of divergence times are surprisingly consistent with previous data from other markers and suggest rapid diversification of lineages within Zea in the last ~150,000 years. We found widespread evidence of historical gene flow, including evidence for divergence in the face of gene flow. We speculate that cultivated maize may serve as a bridge for gene flow among otherwise allopatric wild taxa.     

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.     

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

ABSTRACT: BACKGROUND: Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated ('allopatric') populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations. RESULTS: As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation. CONCLUSIONS: Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of 'genomic islands of speciation'.     

Population history of the Hispaniolan hutia Plagiodontia aedium (Rodentia: Capromyidae): testing the model of ancient differentiation on a geotectonically complex Caribbean island

Hispaniola is a geotectonically complex island consisting of two palaeo-islands that docked c. 10 Ma, with a further geological boundary subdividing the southern palaeo-island into eastern and western regions. All three regions have been isolated by marine barriers during the late Cenozoic and possess biogeographically distinct terrestrial biotas. However, there is currently little evidence to indicate whether Hispaniolan mammals show distributional patterns reflecting this geotectonic history, as the island's endemic land mammal fauna is now almost entirely extinct. We obtained samples of Hispaniolan hutia (Plagiodontia aedium), one of the two surviving Hispaniolan land mammal species, through fieldwork and historical museum collections from seven localities distributed across all three of the island's biogeographic regions. Phylogenetic analysis using mitochondrial DNA (cytochrome b) reveals a pattern of historical allopatric lineage divergence in this species, with the spatial distribution of three distinct hutia lineages biogeographically consistent with the island's geotectonic history. Coalescent modelling, approximate Bayesian computation and approximate Bayes factor analyses support our phylogenetic inferences, indicating near-complete genetic isolation of these biogeographically separate populations and differing estimates of their effective population sizes. Spatial congruence of hutia lineage divergence is not however matched by temporal congruence with divergences in other Hispaniolan taxa or major events in Hispaniola's geotectonic history; divergence between northern and southern hutia lineages dates to c. 0.6 Ma, significantly later than the unification of the palaeo-islands. The three allopatric Plagiodontia populations should all be treated as distinct management units for conservation, with particular attention required for the northern population (low haplotype diversity) and the south-western population (high haplotype diversity but highly threatened).     

Characterization of 10 microsatellite loci in the spiny‐tailed iguana Ctenosaura hemilopha

Ten microsatellite loci were developed for the Baja California spiny‐tailed iguana Ctenosaura hemilopha, using an enriched genomic DNA library. In the Cerralvo Island population, seven loci were polymorphic and presented moderate levels of variation. Number of alleles ranged from two to eight (average 4.43), and observed heterozygosity from 0.150 to 0.857 (average 0.492). Polymorphism was detected at six loci on C. hemilopha individuals from the southern Baja California Peninsula. These markers will be useful to study familiar relationships and behaviour on the Cerralvo Island population.