Historical demography and genetic structure of sister species: deermice (Peromyscus) in the North American temperate rain forest

Phylogeography within the Peromyscus maniculatus species group: Understanding past distribution of genetic diversity and areas of refugia in western North America

The effects of anthropogenic climate change on biodiversity have been recognized on every continent, ocean, and across different taxonomic groups. Here, we study the range dynamics and demography of a cosmopolitan species: the deer mouse, Peromyscus maniculatus. We generated a multilocus SNP dataset using the ddRADseq protocol for 218 individuals across the geographic range within three western North American lineages of this species group. We evaluated population structure using several methods and explored the correlation between geographic and genetic distances. We modeled the demographic history using a site frequency spectrum approach and used a machine learning algorithm to infer current and past (Last Glacial Maximum; LGM) environmental suitability. Lastly, we explored the origin of population expansion for the identified lineages. The genome-wide SNP dataset was able to identify-three regionally distinct groups- 1) P. m. gambelii (southern California); 2) P. keeni (Pacific Northwest); 3) P. m. sonoriensis (a broad population spanning the Pacific Northwest through central California and across the Rocky Mountains into the Great Plains). Demographic analysis indicated the splits between the three populations occurred within the last 500 thousand years, with one very recent (late Holocene) split. Ecological niche models for each of these lineages predicted suitable environment present throughout their known ranges for current conditions, and a severe reduction of northern habitat in the past. The deer mouse has responded to past climate changes by expanding its range during interglacial periods and contracting its range during glacial periods leading to strong population differentiation. But lower magnitude climate change or other processes within the Holocene interglacial period led to population differentiation as well, which is likely still ongoing today given the substantial anthropogenic climate change and other landscape transformations caused by humans during the Anthropocene. By understanding the historical processes that led to the contemporary geographic distribution of biodiversity, we can determine the relative importance of different factors that shape biodiversity, now and into the future.

Short-term responses of small mammal diversity to varying stand-scale patterns of retention tree patches

Retention forestry is a common practice for biodiversity conservation in forests managed for wood production. Retention forestry often leaves unharvested patches of trees that vary in size and spatial pattern but experiments evaluating the effects of different retention patch configurations at a constant level of retention are lacking for many regions and taxonomic groups. We implemented an experimental study in clearcut conifer stands with retention across the U.S. Pacific Northwest region. The study consisted of five stand-level (11-55 ha) experimental treatments each replicated 10 times within a randomized complete block design, resulting in 50 treated stands. Retained tree density was comparable across treatments but size, number, and location (upland or riparian) of patches within stands varied among the five treatments. Within experimental treatments, we measured small mammal (<1kg) species and functional trait (i.e., body size, diet, activity stratum) richness in retention patches, surrounding harvested portions of stands, and nearby unharvested stands. We evaluated species and functional trait richness by treatment using generalized linear mixed-effects models and species-specific responses to retention placement using a community occupancy model. We obtained repeat captures of 21 species of small mammals but found limited evidence of a treatment effect on species richness, and no differences in functional trait richness. Species richness was highest where all retained trees were aggregated into one patch placed adjacent to a forested riparian buffer (mean = 6.6 species, 95% CI = 5.7-7.5), and lowest in the treatment containing one retention patch in the upland portion of a harvested stand (mean = 4.7 species, 95% CI = 3.8-5.6). Furthermore, estimates of species richness within retention patches of harvested stands (i.e., not considering species in harvested areas) did not differ among treatments, indicating that the slightly elevated species richness in riparian-associated retention results from 1-2 species in these patches that do not occur in adjacent harvested portions of each treated stand. Patch occupancy of several species was higher in riparian patches than harvested portions of the treated stands, and fewer species had higher occupancy in upland patches compared to harvested portions of treated stands. Our results indicated that at retention densities currently required in Oregon and Washington, the location of retention patches had a small influence on stand-scale measures of small mammal diversity, but local increases in species richness may be obtained by retaining trees adjacent to riparian buffers.

The genetics of morphological and behavioural island traits in deer mice

Animals on islands often exhibit dramatic differences in morphology and behaviour compared with mainland individuals, a phenomenon known as the 'island syndrome'. These differences are thought to be adaptations to island environments, but the extent to which they have a genetic basis or instead represent plastic responses to environmental extremes is often unknown. Here, we revisit a classic case of island syndrome in deer mice (Peromyscus maniculatus) from British Columbia. We first show that Saturna Island mice and those from neighbouring islands are approximately 35% (approx. 5 g) heavier than mainland mice and diverged approximately 10 000 years ago. We then establish laboratory colonies and find that Saturna Island mice are heavier both because they are longer and have disproportionately more lean mass. These trait differences are maintained in second-generation captive-born mice raised in a common environment. In addition, island-mainland hybrids reveal a maternal genetic effect on body weight. Using behavioural testing in the laboratory, we also find that wild-caught island mice are less aggressive than mainland mice; however, laboratory-raised mice born to these founders do not differ in aggression. Together, our results reveal that these mice have different responses to the environmental conditions on islands-a heritable change in a morphological trait and a plastic response in a behavioural trait.

Living on the edge: Exploring the role of coastal refugia in the Alexander Archipelago of Alaska

Although islands are of long-standing interest to biologists, only a handful of studies have investigated the role of climatic history in shaping evolutionary diversification in high-latitude archipelagos. In this study of the Alexander Archipelago (AA) of Southeast Alaska, we address the impact of glacial cycles on geographic genetic structure for three mammals co-distributed along the North Pacific Coast. We examined variation in mitochondrial and nuclear loci for long-tailed voles (Microtus longicaudus), northwestern deermice (Peromyscus keeni), and dusky shrews (Sorex monticola), and then tested hypotheses derived from Species Distribution Models, reconstructions of paleoshorelines, and island area and isolation. In all three species, we identified paleoendemic clades that likely originated in coastal refugia, a finding consistent with other paleoendemic lineages identified in the region such as ermine. Although there is spatial concordance at the regional level for endemism, finer scale spatial and temporal patterns are less clearly defined. Demographic expansion across the region for these distinctive clades is also evident and highlights the dynamic history of Late Quaternary contraction and expansion that characterizes high-latitude species.

The ultimate and proximate mechanisms driving the evolution of long tails in forest deer mice

Understanding both the role of selection in driving phenotypic change and its underlying genetic basis remain major challenges in evolutionary biology. Here, we use modern tools to revisit a classic system of local adaptation in the North American deer mouse, Peromyscus maniculatus, which occupies two main habitat types: prairie and forest. Using historical collections, we find that forest‐dwelling mice have longer tails than those from nonforested habitat, even when we account for individual and population relatedness. Using genome‐wide SNP data, we show that mice from forested habitats in the eastern and western parts of their range form separate clades, suggesting that increased tail length evolved independently. We find that forest mice in the east and west have both more and longer caudal vertebrae, but not trunk vertebrae, than nearby prairie forms. By intercrossing prairie and forest mice, we show that the number and length of caudal vertebrae are not correlated in this recombinant population, indicating that variation in these traits is controlled by separate genetic loci. Together, these results demonstrate convergent evolution of the long‐tailed forest phenotype through two distinct genetic mechanisms, affecting number and length of vertebrae, and suggest that these morphological changes—either independently or together—are adaptive.

Evidence for Host-Genotype Associations of Borrelia burgdorferi Sensu Stricto

Different genotypes of the agent of Lyme disease in North America, Borrelia burgdorferi sensu stricto, show varying degrees of pathogenicity in humans. This variation in pathogenicity correlates with phylogeny and we have hypothesized that the different phylogenetic lineages in North America reflect adaptation to different host species. In this study, evidence for host species associations of B. burgdorferi genotypes was investigated using 41 B. burgdorferi-positive samples from five mammal species and 50 samples from host-seeking ticks collected during the course of field studies in four regions of Canada: Manitoba, northwestern Ontario, Quebec, and the Maritimes. The B. burgdorferi genotypes in the samples were characterized using three established molecular markers (multi-locus sequence typing [MLST], 16S-23S rrs-rrlA intergenic spacer, and outer surface protein C sequence [ospC] major groups). Correspondence analysis and generalized linear mixed effect models revealed significant associations between B. burgdorferi genotypes and host species (in particular chipmunks, and white-footed mice and deer mice), supporting the hypotheses that host adaptation contributes to the phylogenetic structure and possibly the observed variation in pathogenicity in humans.

Phylogeographic Structure of the White-Footed Mouse and the Deer Mouse, Two Lyme Disease Reservoir Hosts in Québec

Modification of a species range is one of many consequences of climate change and is driving the emergence of Lyme disease in eastern Canada. The primary reservoir host of the bacteria responsible for Lyme disease, Borrelia burgdorferi, is the white-footed mouse (Peromyscus leucopus), whose range is rapidly shifting north into southern Québec. The deer mouse, P. maniculatus, is occurring over most Québec province and is a less competent host for B. burgdorferi. Here, we compared the phylogeographic structure of both Peromyscus species in Québec. Using a combination of multiple mitochondrial DNA markers and phylogeographic methods, we detected an ongoing and rapid expansion of P. leucopus, while P. maniculatus appears more stable. Haplotype and populations networks indicated that populations of P. maniculatus exhibit more genetic structure than P. leucopus across the study area. Furthermore, significant and consistent genetic divergences between populations of the two species on both sides of the St. Lawrence River suggest that distinct lineages of P. leucopus and P. maniculatus with different ancestral origins colonized Southern Québec following the Last Glacial Maximum. The phylogeographic structure of pathogens is expected to mirror the structure observed in their reservoir hosts. As different strains of Borrelia burgdorferi may be associated with different levels of pathogenicity and immune responses of their hosts, our results are helpful at better understanding the pattern of spread of Lyme disease in a zone of emergence, and associated risk for human populations.

The endemic insular and peninsular species Chaetodipus spinatus (Mammalia, Heteromyidae) breaks patterns for Baja California

The Baja California peninsula is the second longest, most geographically isolated peninsula on Earth. Its physiography and the presence of many surrounding islands has facilitated studies of the underlying patterns and drivers of genetic structuring for a wide spectrum of organisms. Chaetodipus spinatus is endemic to the region and occurs on 12 associated islands, including 10 in the Gulf of California and two in the Pacific Ocean. This distribution makes it a model species for evaluating natural historical barriers. We test hypotheses associated with the relationship between the range of the species, patterns in other species, and its relationship to Pleistocene-Holocene climatic changes. We analyzed sequence data from mtDNA genes encoding cytochrome b (Cytb) and cytochrome c oxidase subunits I (COI) and III (COIII) in 26 populations including all 12 islands. The matrilineal genealogy, statistical parsimony network and Bayesian skyline plot indicated an origin of C. spinatus in the southern part of the peninsula. Our analyses detected several differences from the common pattern of peninsular animals: no mid-peninsula break exists, Isla Carmen hosts the most divergent population, the population on an ancient southern Midriff island does not differ from peninsular populations, and a mtDNA peninsular discordance occurs near Loreto.

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high-amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model-averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.

Nuclear and mitochondrial DNA reveal contrasting evolutionary processes in populations of deer mice (Peromyscus maniculatus)

We investigated a major geographic break in mitochondrial DNA (mtDNA) haplotypes in deer mice, Peromyscus maniculatus, by analysing spatial variation in a 491-bp fragment of the mtDNA control region from 455 samples distributed across a north-south transect of 2000 km in Western North America. To determine whether the mtDNA break was reflected in the nuclear genome, we then compared spatial variation in 13 nuclear microsatellites of 95 individuals surrounding the mtDNA break. Using a canonical correlation analysis we found that nuclear genomic variation was not correlated with mtDNA differentiation. The contrasting patterns of variation in mtDNA and nuclear DNA are consistent with a hypothesis of historic genetic drift that occurred in isolated refugia combined with recent gene flow between the formerly isolated refugial populations. A Mantel test of genetic vs. geographic distance revealed that recent gene flow between deer mouse populations has been high. We conclude that past vicariant events associated with Pleistocene climate changes together with recent gene flow have created the observed intra-specific cytonuclear discordance in Western North America.

High altitude adaptation and phylogenetic analysis of Tibetan horse based on the mitochondrial genome

To investigate genetic mechanisms of high altitude adaptations of animals living in the Tibetan Plateau, three mitochondrial genomes (mt-genome) of Tibetan horses living in Naqu (4,500 m) of Tibetan, Zhongdian (3,300 m) and Deqin (3,100 m) of Yunnan province were sequenced. The structures and lengths of these three mt-genomes are similar to the Cheju horse, which is related to Tibetan horses, but little shorter than the Swedish horse. The pair-wise identity of these three horses on nucleotide level is more than 99.3%. When the gene encoding the mitochondrial protein of Tibetan horses was analyzed, we found that NADH6 has higher non-synonymous mutation rate in all of three Tibetan horses. This implies that NADH6 may play a role in Tibetan horses' high altitude adaptation. NADH6 is one of the subunits of the complex I in the respiratory chain. Furthermore, 7 D-loop sequences of Tibetan horse from different areas were sequenced, and the phylogeny tree was constructed to study the origin and evolutionary history of Tibetan horses. The result showed that the genetic diverse was high among Tibetan horses. All of Tibetan horses from Naqu were clustered into one clade, and Tibetan horses from Zhongdian and Deqin were clustered into others clades. The first molecular evidence of Tibetan horses indicated in this study is that Tibetan horse population might have multiple origins.

Intraspecific patterns of mitochondrial variation in natural population fragments of a localized desert dung beetle species, Pachysoma gariepinum (Coleoptera: Scarabaeidae)

Phylogenetic, population, and coalescent methods were used to examine the genetic structuring of Pachysoma gariepinum, a flightless dung beetle species endemic to the arid west coast of southern Africa that exhibits interrupted south to north morphological clinal variation along a distributional gradient. Mitochondrial cytochrome oxidase I sequence data of 67 individuals from 5 localities revealed the presence of 3 geographically distinct evolutionary lineages (with an overall nucleotide divergence of 5.7% and a per-locality divergence of 1.9-3.8%) which display significant levels of genetic structuring. The separation of the lineages was estimated to have occurred between 2.2 and 5.7 million years ago--which is the late Miocene, early Plio-Pleistocene era--possibly in response to the ebb and flow of the Orange and Holgat River systems as well as the interactions between the moving and stable sand dune systems. Moreover the species' current range appears to have been influenced by the formation of advective fog resulting in a constant source of water in an area with low precipitation thereby allowing for the beetles to radiate to areas that were previously inhospitable. Fu's F-statistics and population parameters based on recent mutations indicated that little to no recent population growth has occurred. This together with changing anthropogenic factors and the recovery of 3 geographically discrete management units, points to a need for census data in order to monitor and conserve the genetic diversity of this species.

Forest-obligate Sabethes mosquitoes suggest palaeoecological perturbations

The origin of tropical forest diversity has been hotly debated for decades. Although specific mechanisms vary, many such explanations propose some vicariance in the distribution of species during glacial cycles and several have been supported by genetic evidence in Neotropical taxa. However, no consensus exists with regard to the extent or time frame of the vicariance events. Here, we analyse the cytochrome oxidase II mitochondrial gene of 250 Sabethes albiprivus B mosquitoes sampled from western Sao Paulo in Brazil. There was very low population structuring among collection sites (Phi(ST)=0.03, P=0.04). Historic demographic analyses and the contemporary geographic distribution of genetic diversity suggest that the populations sampled are not at demographic equilibrium. Three distinct mitochondrial clades were observed in the samples, one of which differed significantly in its geographic distribution relative to the other two within a small sampling area (approximately 70 x 35 km). This fact, supported by the inability of maximum likelihood analyses to achieve adequate fits to simple models for the population demography of the species, suggests a more complex history, possibly involving disjunct forest refugia. This hypothesis is supported by a genetic signal of recent population growth, which is expected if population sizes of this forest-obligate insect increased during the forest expansions that followed glacial periods. Although a time frame cannot be reliably inferred for the vicariance event leading to the three genetic clades, molecular clock estimates place this at approximately 1 Myr before present.

Coevolutionary relationship between helminth diversity and MHC class II polymorphism in rodents

Parasite-mediated selection on major histocompatibility complex (MHC) genes has mainly been explored at the intraspecific level, although many molecular studies have revealed trans-species polymorphism. Interspecific patterns of MHC diversity might reveal factors responsible for the long-term evolution of MHC polymorphism. We hypothesize that host taxa harbouring high parasite diversity should exhibit high levels of MHC genetic diversity. We test this assumption using data on rodent species and their helminth parasites compiled from the literature. Controlling for similarity due to common descent, we present evidence indicating that high helminth species richness in rodent species is associated with increased MHC class II polymorphism. Our results are consistent with the idea that parasites sharing a long-term coevolutionary history with their hosts are the agents of selection explaining MHC polymorphism.

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.

Cytochrome-b haplotypes suggest an undescribed Peromyscus species from the Yukon

Peromyscus maniculatus (Wagner, 1845) and Peromyscus keeni Merriam, 1897 are two species of deer mouse currently recognized in the Yukon. Phylogenetic analyses (Kimura two-parameter and maximum parsimony) of cytochrome-b sequences (560 base pairs) from deer mouse specimens (n = 4) collected near Haines Junction, Yukon, resulted in a monophyletic clade genetically distant from the two currently recognized species. We suggest that the Haines Junction specimens may represent a previously undescribed peromyscine species. Peromyscus arcticus is an available name.

Comparative phylogeography of eastern chipmunks and white-footed mice in relation to the individualistic nature of species

Palaeoecological studies have demonstrated that ecological communities as a whole did not remain stable throughout the climatic fluctuations of the Quaternary. The result is that long-term associations of species cannot be inferred by contemporary associations in ecological communities. Therefore, the evolutionary significance of any contemporary ecological interactions among species and of the biotic community within which species have evolved also cannot be assumed from contemporary conditions. Comparative phylogeographic data provide a method to identify species within ecological communities that have shared biogeographic histories. We present an example of a long-term association between populations of two mammalian species, eastern chipmunks (Tamias striatus) and white-footed mice (Peromyscus leucopus), which are commonly associated with deciduous forest habitats. The distribution of mitochondrial DNA variation in T. striatus and P. leucopus from previously glaciated regions of the eastern United States support the hypothesis that, in at least part of their range, genetic lineages of the two species have expanded from similar population sources since the Last Glacial Maximum. In addition, the spatial concordance of genetic lineages of T. striatus and P. leucopus with the oak-savannah forest formations of Wisconsin and Illinois, suggest that populations associated with this community colonized the area in association with a set of arboreal species that comprise their deciduous forest habitat.

Effects of historical climate change, habitat connectivity, and vicariance on genetic structure and diversity across the range of the red tree vole (Phenacomys longicaudus) in the Pacific Northwestern United States

Phylogeographical analyses conducted in the Pacific Northwestern United States have often revealed concordant patterns of genetic diversity among taxa. These studies demonstrate distinct North/South genetic discontinuities that have been attributed to Pleistocene glaciation. We examined phylogeographical patterns of red tree voles (Phenacomys longicaudus) in western Oregon by analysing mitochondrial control region sequences for 169 individuals from 18 areas across the species' range. Cytochrome b sequences were also analysed from a subset of our samples to confirm the presence of major haplotype groups. Phylogenetic network analyses suggested the presence of two haplotype groups corresponding to northern and southern regions of P. longicaudus' range. Spatial genetic analyses (samova and Genetic Landscape Shapes) of control region sequences demonstrated a primary genetic discontinuity separating northern and southern sampling areas, while a secondary discontinuity separated northern sampling areas into eastern and western groups divided by the Willamette Valley. The North/South discontinuity likely corresponds to a region of secondary contact between lineages rather than an overt barrier. Although the Cordilleran ice sheet (maximum approximately 12,000 years ago) did not move southward to directly affect the region occupied by P. longicaudus, climate change during glaciation fragmented the forest landscape that it inhabits. Signatures of historical fragmentation were reflected by positive associations between latitude and variables such as Tajima's D and patterns associated with location-specific alleles. Genetic distances between southern sampling areas were smaller, suggesting that forest fragmentation was reduced in southern vs. northern regions.

A mitochondrial genome sequence of the Tibetan antelope (Pantholops hodgsonii)

To investigate genetic mechanisms of high altitude adaptations of native mammals on the Tibetan Plateau, we compared mitochondrial sequences of the endangered Pantholops hodgsonii with its lowland distant relatives Ovis aries and Capra hircus, as well as other mammals. The complete mitochondrial genome of P. hodgsonii (16,498 bp) revealed a similar gene order as of other mammals. Because of tandem duplications, the control region of P. hodgsonii mitochondrial genome is shorter than those of O. aries and C. hircus, but longer than those of Bos species. Phylogenetic analysis based on alignments of the entire cytochrome b genes suggested that P. hodgsonii is more closely related to O. aries and C. hircus, rather than to species of the Antilopinae subfamily. The estimated divergence time between P. hodgsonii and O. aries is about 2.25 million years ago. Further analysis on natural selection indicated that the COXI (cytochrome c oxidase subunit I) gene was under positive selection in P. hodgsonii and Bos grunniens. Considering the same climates and environments shared by these two mammalian species, we proposed that the mitochondrial COXI gene is probably relevant for these native mammals to adapt the high altitude environment unique to the Tibetan Plateau.

Phylogeny and historical biogeography of African ground squirrels: the role of climate change in the evolution of Xerus

We used phylogenetic and phylogeographical methods to infer relationships among African ground squirrels of the genus Xerus. Using Bayesian, maximum-parsimony, nested clade and coalescent analyses of cytochrome b sequences, we inferred interspecific relationships, evaluated the specific distinctness of Cape (Xerus inauris) and mountain (Xerus princeps) ground squirrels, and tested hypotheses for historical patterns of gene flow within X. inauris. The inferred phylogeny supports the hypothesized existence of an 'arid corridor' from the Horn of Africa to the Cape region. Although doubts have been raised regarding the specific distinctness of X. inauris and X. princeps, our analyses show that each represents a distinct well-supported, monophyletic lineage. Xerus inauris includes three major clades, two of which are geographically restricted. The distributions of X. inauris populations are concordant with divergences within and disjunctions between other taxa, which have been interpreted as results of Plio-Pleistocene climate cycles. Nested clade analysis, coalescent analyses, and analyses of genetic structure support allopatric fragmentation as the cause of the deep divergences within this species.

Genetic structure, phylogeography and demography of two ground-beetle species endemic to the Tenerife laurel forest (Canary Islands)

The volcanic island of Tenerife (Canary archipelago) was formerly covered at 600-1200 m above sea level on most of its northern side by a cloud forest holding much of the endemic insect fauna. In the most significant surviving patches of this laurel forest at the eastern and western tips of the island occur two forest-specialist, closely related species of Eutrichopus (Coleoptera, Carabidae); here we present data on mitochondrial DNA variation among populations of these species. In total, 116 individuals from 16 localities were sampled and a 638 bp fragment of the cytochrome oxidase subunit II gene was sequenced, obtaining evidence for two distinct evolutionary lineages, in accordance with morphological and biogeographical data. Volcanic events at approximately 0.7 Ma might be responsible for vicariance and the fragmentation of the geographical range of an ancestral species, causing the establishment of two matrilineal lineages. Using nested clade and historical demography analyses we infer past cycles of demographic bottlenecks followed by population expansion, mostly in agreement with the geological time scale of volcanic events. Recent trends, however, refer to fragmentation of the cloud forest due to human intervention.

Genetic structure of desert ground squirrels over a 20-degree-latitude transect from Oregon through the Baja California peninsula

The genetic structure of populations over a wide geographical area should reflect the demographic and evolutionary processes that have shaped a species across its range. We examined the population genetic structure of antelope ground squirrels (Ammospermophilus leucurus) across the complex of North American deserts from the Great Basin of Oregon to the cape region of the Baja California peninsula. We sampled 73 individuals from 13 major localities over this 2500-km transect, from 43 to 22 degrees north. Our molecular phylogeographical analysis of 555 bp of the mitochondrial cytochrome b gene and 510 bp of the control region revealed great genetic uniformity in a single clade that extends from Oregon to central Baja California. A second distinct clade occupies the southern half of the peninsula. The minimal geographical structure of the northern clade, its low haplotype diversity and the distribution of pairwise differences between haplotypes suggest a rapid northward expansion of the population that must have followed a northward desert habitat shift associated with the most recent Quaternary climate warming and glacial retreat. The higher haplotype diversity within the southern clade and distribution of pairwise differences between haplotypes suggest that the southern clade has a longer, more stable history associated with a southern peninsular refugium. This system, as observed, reflects both historical and contemporary ecological and evolutionary responses to physical environmental gradients within genetically homogeneous populations.

Geographic genetic architecture of pocket gopher (Thomomys bottae) populations in Baja California, Mexico

Phylogenetic analyses of complete mitochondrial cytochrome b sequences support the monophyly of pocket gopher (Thomomys bottae) populations from the 1000 km length of the Baja California peninsula of Mexico, relative to other geographical segments of the species range in western North America. The Baja California peninsula is an area that encompasses considerable ecomorphological and infraspecific diversity within this pocket gopher species. However, detailed population analyses encompassing 35 localities distributed over the southern half of the peninsula reveal only trivial phylogeographical structure. Rather, most of the 72 unique 500-base pair haplotypes examined from 142 individuals is restricted to single populations, although a few haplotypes are shared broadly across geography. Individual populations are typically comprised of haplotype sets from different branches in a network of relationships. Analysis of molecular variance (amova) indicates that approximately half of the total pool of variation is contained among individuals within local populations, and that only about 25% can be explained by the regional subdivisions of current subspecies distributions or physiographic realms. A hypothesized historical vicariant event that has been causally linked to the phylogeographical structure of other, codistributed species has had little influence on these pocket gopher populations, explaining only 13% of the total variation. The temporal depth, estimated by coalescence parameters, of the haplotype lineage in Baja California is relatively recent, approximately 300,000 generations; both the mismatch distribution of pairwise comparisons and a significantly positive exponential growth estimate support a recent history of expanding populations; but current, or recent past, migration estimates have remained small, are largely unidirectional from north to south, and weak isolation by distance is present. All data suggest that pocket gophers have relatively recently invaded the southern half of peninsular Baja California, with the genetic signature of expansion still evident but with sufficient time having lapsed to result in a weak isolation by distance pattern. The geographical assemblage of sampled populations thus appears as a meta-population, with limited gene flow contrasting with random haplotype loss due to drift in small, localized populations.

Testing for asymmetrical gene flow in a Drosophila melanogaster body-size cline

Asymmetrical gene flow is an important, but rarely examined genetic parameter. Here, we develop a new method for detecting departures from symmetrical migration between two populations using microsatellite data that are based on the difference in the proportion of private alleles. Application of this approach to data collected from wild-caught Drosophila melanogaster along a latitudinal body-size cline in eastern Australia revealed that asymmetrical gene flow could be detected, but was uncommon, nonlocalized, and occurred in both directions. We also show that, in contrast to the findings of a previous study, there is good evidence to suggest that the cline experiences significant levels of gene flow between populations.

Phylogeography of the coronulid barnacle, Chelonibia testudinaria, from loggerhead sea turtles, Caretta caretta

The barnacle, Chelonibia testudinaria, is a common inhabitant of the marine turtle epibiont community and plays a key role in the development of this community. Phylogeographic analysis of 79 cytochrome c oxidase I (COX1) sequences for barnacles collected from five populations found contrasting patterns of divergence for populations in the Atlantic vs. the Pacific Ocean. Our analysis indicates that the two Pacific populations, Senri Beach, Japan and Bahia Magdalena, Mexico, were not only highly divergent from the Atlantic populations but are highly divergent from one another. We suggest that barnacles from these populations may represent cryptic species. In contrast, sequence divergence was greatly reduced among barnacles collected from Wassaw Island, GA, USA, Keewaydin, FL, USA, and Kyparissia, Pèloponnésus Island, Greece. A reduction in sequence diversity at the latter site was attributed to a recent range expansion into the Mediterranean Sea. We examined historical patterns of migration among the three Atlantic and Mediterranean populations using the program migrate. This analysis indicates a high rate of migration from Keewaydin to Wassaw Island, contrasted with a much lower rate of migration in the opposite direction. The estimated migration rate from Kyparissia to Keewaydin was also non-negligible. We suggest that the association between C. testudinaria and loggerhead turtles and the patterns of turtle migration have played key roles in the expansion of the range of C. testudinaria into the Mediterranean Sea and the subsequent patterns of barnacle migration. We further propose that the difference between ocean basins, with respect to the impact of host migration on barnacle gene flow, probably stems from the fact that host-mediated dispersal in the Atlantic depends on advanced stage juveniles and adults while any host-mediated dispersal in the Pacific would have to involve early "pelagic" stage juvenile loggerheads.