Distinguishing the effects of selection from demographic history in the genetic variation of two sister passerines based on mitochondrial–nuclear comparison

Determining the mechanisms responsible for the distribution of genetic diversity in natural populations has occupied a central role in molecular evolution. Our study was motivated by the unprecedented observation that a widespread Eurasian flycatcher, Ficedula albicilla, exhibited no variation at the mitochondrial DNA (mtDNA) ND2 gene in 75 individuals sampled over a 5000-km distance. In contrast, its sister species, F. parva, had low but considerably higher levels of mtDNA variation. We assessed whether natural selection or demographic factors could explain the absence of mtDNA variation in F. albicilla. Eighteen nuclear genes were sequenced to estimate the two species'' phylogeographic histories, and for comparison to the mtDNA data. Multilocus coalescence analyses suggested that F. albicilla experienced a population expansion perhaps following a population bottleneck. Simulations based on this demographic history, however, did not replicate the extremely low level of mtDNA variation. Historical range changes based on ecological niche models also failed to explain the observed mtDNA patterns. Neutrality tests (DHEW and ML-HKA) suggested a non-neutral pattern in the mtDNA of F. albicilla. We found a transmembrane-skewed distribution of nonsynonymous substitutions between the two species, three of which caused functional change; the results implied that positive selection could have targeted mtDNA. Several lines of evidence support selection rather than demographic history as the main force influencing the patterns of mtDNA variation. Despite the influence of natural selection, many of the phylogeographic inferences derived from mtDNA were robust, including species limits and a high level of gene flow among populations within species.     

Post-glacial partitioning of mitochondrial genetic variation in the field vole

Genetic markers are often used to examine population history. There is considerable debate about the behaviour of molecular clock rates around the population-species transition. Nevertheless, appropriate calibration is critical to any inference regarding the absolute timing and scale of demographic changes. Here, we use a mitochondrial cytochrome b gene genealogy, based entirely on modern sequences and calibrated from recent geophysical events, to date the post-glacial expansion of the Eurasian field vole (Microtus agrestis), a widespread temperate mammal species. The phylogeographic structure reflects the subsequent expansion of populations that went through bottlenecks at the time of the Younger Dryas (ca 12,000 years BP) rather than the Last Glacial Maximum (LGM, ca 24,000 years BP), which is usually seen as the time when present-day patterns were determined. The nucleotide substitution rate that was estimated here, ca 4 × 10(-7) substitutions/site/year, remains extremely high throughout the relevant time frame. Calibration with similarly high population-based substitution rates, rather than long-term rates derived from species divergence times, will show that post-LGM climatic events generated current phylogeographic structure in many other organisms from temperate latitudes.     

Patterns of genetic diversity in the critically endangered Central American river turtle: human influence since the Mayan age?

We conducted a phylogeographic analysis of the strictly aquatic and critically endangered Central American river turtle, Dermatemys mawii, as part of a conservation management program for the species. We sampled 238 individuals from 15 different localities throughout the species range. Using sequence fragments from the mtDNA Cyt b and ND4 genes, we identified 16 different haplotypes. Overall, our results reveal a signal of phylogeographic structure throughout the range, which appears to have been secondarily blurred by extensive gene flow. Notably, this also applies to genetic structuring across three major hydrological basins that pose biogeographic breaks in other aquatic taxa. Divergence times of mtDNA haplotypes in D. mawii suggest that the main lineages split in the Pliocene–Pleistocene (3.73–0.227 MA) and demographic tests indicate that the species has undergone drastic demographic size fluctuations since this time period. One ancient haplotype (1D) was found to exhibit sequence divergence of up to 2% from other haplogroups. Divergence of this magnitude is indicative of species level differentiation in other turtle genera. Haplotype 1D was found in only two localities, Sarstun and Salinas, but specimens with other haplotypes were also found in those localities. It is not known whether the individuals with the 1D haplotype interbreed with non-1D individuals. Our results suggest that human activity, such as harvesting and long distance transport of animals, may have influenced the current patterns of genetic diversity. For more than 2000 years, D. mawii has been consumed by people from Middle American cultures, and the archeological record contains strong evidence that the Mayans transported animals between villages and far away from their natural distribution range. Therefore, the large-scale pattern of haplotype sharing even across hydrological barriers, the observed low haplotype diversity in some populations and the contemporary absence of a pronounced phylogeographic pattern is likely due to a combination of population expansions, gene flow, extensive human-mediated-movements and recent bottlenecks resulting from over-harvesting.     

Selective processes and adaptive evolution of the cytochrome b gene in salamanders of the genus Salamandrella

Sequence analysis of the cytochrome b gene fragment in the salamanders of the genus Salamandrella, Siberian salamander and Schrenk salamander was performed with the purpose to elucidate the effect of natural selection on the evolution of mitochondrial DNA (mtDNA) in these species. It was demonstrated that despite of notable influence of negative selection (expressed as very low dN/dS values), speciation and intraspecific divergence in salamanders was accompanied by the appearance of radical amino acid substitutions, caused by the influence of positive (directional) selection. To examine the evolutionary pattern of synonymous mtDNA sites, distribution of conservative and non-conservative substitutions was analyzed. The rates of conservative and non-conservative substitutions were nearly equal, pointing to neutrality of mutation process at synonymous mtDNA sites of salamanders. Analysis of conservative and non-conservative synonymous substitution distributions in different parts of phylogenetic trees showed that the differences between the synonymous groups compared were statistically significant only in one phylogenetic group of Siberian salamander (haplogroup C) (P = 0.02). In the group of single substitutions, located at terminal phylogenetic branches of Siberian salamanders from group C, increased rate of conservative substitutions was observed. Based on these findings, it was suggested that selective processes could have an influence on the formation of the synonymous substitution profile in the Siberian salamander mtDNA fragment examined.     

The complex biogeographic history of a widespread tropical tree species

Many tropical forest tree species have broad geographic ranges, and fossil records indicate that population disjunctions in some species were established millions of years ago. Here we relate biogeographic history to patterns of population differentiation, mutational and demographic processes in the widespread rainforest tree Symphonia globulifera using ribosomal (ITS) and chloroplast DNA sequences and nuclear microsatellite (nSSR) loci. Fossil records document sweepstakes dispersal origins of Neotropical S. globulifera populations from Africa during the Miocene. Despite historical long-distance gene flow, nSSR differentiation across 13 populations from Costa Rica, Panama, Ecuador (east and west of Andes) and French Guiana was pronounced (F(ST)= 0.14, R(ST)= 0.39, P < 0.001) and allele-size mutations contributed significantly (R(ST) > F(ST)) to the divergences between cis- and trans-Andean populations. Both DNA sequence and nSSR data reflect contrasting demographic histories in lower Mesoamerica and Amazonia. Amazon populations show weak phylogeographic structure and deviation from drift-mutation equilibrium indicating recent population expansion. In Mesoamerica, genetic drift was strong and contributed to marked differentiation among populations. The genetic structure of S. globulifera contains fingerprints of drift-dispersal processes and phylogeographic footprints of geological uplifts and sweepstakes dispersal.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Substitution bias, rapid saturation, and the use of mtDNA for nematode systematics

Only relatively recently have researchers turned to molecular methods for nematode phylogeny reconstruction. Thus, we lack the extensive literature on evolutionary patterns and phylogenetic usefulness of different DNA regions for nematodes that exists for other taxa. Here, we examine the usefulness of mtDNA for nematode phylogeny reconstruction and provide data that can be used for a priori character weighting or for parameter specification in models of sequence evolution. We estimated the substitution pattern for the mitochondrial ND4 gene from intraspecific comparisons in four species of parasitic nematodes from the family Trichostrongylidae (38-50 sequences per species). The resulting pattern suggests a strong mutational bias toward A and T, and a lower transition/transversion ratio than is typically observed in other taxa. We also present information on the relative rates of substitution at first, second, and third codon positions and on relative rates of saturation of different types of substitutions in comparisons ranging from intraspecific to interordinal. Silent sites saturate extremely quickly, presumably owing to the substitution bias and, perhaps, to an accelerated mutation rate. Results emphasize the importance of using only the most closely related sequences in order to infer patterns of substitution accurately for nematodes or for other taxa having strongly composition-biased DNA. ND4 also shows high amino acid polymorphism at both the intra- and interspecific levels, and in higher level comparisons, there is evidence of saturation at variable amino acid sites. In general, we recommend using mtDNA coding genes only for phylogenetics of relatively closely related nematode species and, even then, using only nonsynonymous substitutions and the more conserved mitochondrial genes (e.g., cytochrome oxidases). On the other hand, the high substitution rate in genes such as ND4 should make them excellent for population genetics studies, identifying cryptic species, and resolving relationships among closely related congeners when other markers show insufficient variation.      

Speciation in an avian complex endemic to the mountains of Middle America (Ergaticus, Aves: Parulidae)

The implementation of the phylogeographic approach for the study of biodiversity is critical in poorly sampled regions like the montane systems of Middle America, as complex evolutionary histories often result in the presence of independent lineages not properly considered by traditional taxonomy. Herein we sequenced 2370 bp of mtDNA (ND2, cyt b and ATPase) from 81 individuals of Ergaticus, a complex of birds endemic to the montane forests of Middle America. Although current taxonomy recognizes two species, the results reveal considerable genetic structure with the presence of four mtDNA lineages. Two of these lineages within Ergaticus ruber evidence the need of a revaluation of the species limits for this taxon. The general phylogeographic pattern can be explained as a consequence of relative isolation of the populations in different mountain ranges separated by low elevation barriers. Most population groups did not show signals of demographic expansion with the exception of the one corresponding to clade 1. The divergence time estimates point to the Pleistocene as an important time period for the diversification of this complex.     

Nonneutral Mitochondrial DNA Variation in Humans and Chimpanzees

We sequenced the NADH dehydrogenase subunit 3 (ND3) gene from a sample of 61 humans, five common chimpanzees, and one gorilla to test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution. Within humans and within chimpanzees, the ratio of replacement to silent nucleotide substitutions was higher than observed in comparisons between species, contrary to neutral expectations. To test the generality of this result, we reanalyzed published human RFLP data from the entire mitochondrial genome. Gains of restriction sites relative to a known human mtDNA sequence were used to infer unambiguous nucleotide substitutions. We also compared the complete mtDNA sequences of three humans. Both the RFLP data and the sequence data reveal a higher ratio of replacement to silent nucleotide substitutions within humans than is seen between species. This pattern is observed at most or all human mitochondrial genes and is inconsistent with a strictly neutral model. These data suggest that many mitochondrial protein polymorphisms are slightly deleterious, consistent with studies of human mitochondrial diseases.     

Mitochondrial genetic variability of Didelphis albiventris (Didelphimorphia, Didelphidae) in Brazilian localities

Didelphis albiventris is a well-known and common marsupial. Due to its high adaptability, this very widespread generalist species occurs under various environmental conditions, this even including protected regions and disturbed urban areas. We studied a 653 bp fragment of cytochrome oxidase c (COI) from 93 biological samples from seven Brazilian localities, with linear distances ranging between 58 and about 1800 km to analyze the effects of geographic distances on variability and genetic differentiation. The haplotype network presented nine haplotypes and two genetic clusters compatible with the two most distant geographic areas of the states of Minas Gerais, in the southeast, and Rio Grande do Sul, in the extreme south. As each cluster was characterized by low nucleotide and high haplotype diversities, their populations were obviously composed of closely related haplotypes. Surprisingly, moderate to high F(ST) differentiation values and a very weak phylogeographic signal characterizes interpopulation comparisons within Minas Gerais interdemes, these being correlated with the presence of privative haplotypes. On a large rgeographic scale, a comparison between demes from Minas Gerais and Rio Grande do Sul presented high F(ST) values and a robust phylogeographic pattern. This unexpected scenario implies that mtDNA gene flow was insufficient to maintain population cohesion, reflected by the observed high differentiation.     

Phylogeography of Fischer’s blue,Tongeia fischeri,in Japan: Evidence for introgressive hybridization

The widespread lycaenid butterfly Tongeia fischeri is distributed from eastern Europe to northeastern Asia and represented by three geographically isolated populations in Japan. In order to clarify the phylogeographic history of the species, we used sequences of three mitochondrial (COI, Cyt b and ND5) and two nuclear (Rpl5 and Ldh) genes of 207 individuals collected from 55 sites throughout Japan and five sites on the Asian continent. Phylogenetic trees and the median-joining network revealed six evolutionary mitochondrial haplotype clades, which corresponded to the geographic distribution of the species. Common ancestors of Japanese T. fischeri might have come to Japan during the mid-Pleistocene by multiple dispersals of continental populations, probably via a land bridge or narrow channel between western Japan and the Korean Peninsula. The geographical patterns of variation of mitochondrial and nuclear markers are discordant in northeastern Kyushu, possibly as a result of introgressive hybridization during the ancient contact between the Kyushu and Shikoku populations in the last glacial maximum. The phylogeographic pattern of T. fischeri in Japan are probably related to the geological history, Pleistocene climatic oscillations and distribution of the host plant.     

Limited phylogeographical structure across Eurasia in two red wood ant species Formica pratensis and F. lugubris (Hymenoptera, Formicidae)

The phylogeography and demographic history of two closely related species of the red wood ant (Formica pratensis and F. lugubris) were examined across Eurasia. The phylogeny based on a 1.5-kilobase mitochondrial DNA fragment, including the cytochrome b gene and part of the ND6 gene, showed one phylogeographical division in F. pratensis. This division (0.7% of nucleotide divergence) suggests postglacial colonization of western Europe and of a wide area ranging from Sweden on the west to Lake Baikal on the east from separate forest refugia. In two localities, mitochondrial DNA has been transferred from F. lugubris to F. pratensis and all the individuals of F. pratensis sampled from the Pyrenees had haplotypes clustering in the lugubris clade. No phylogeographical divisions were detected in F. lugubris. Comparison of species-wide phylogeography between the two sympatrically distributed species of ant demonstrates a difference in phylogeographical structure that implies different vicariant histories. However, over most of the species' distribution ranges, similar signs of demographic expansion predating the last glaciation and the lack of phylogeographical structure were found in both the eastern phylogroup of F. pratensis and F. lugubris. This finding is highly consistent with the results reported for all other boreal forest animal species studied to date in Eurasia. Contraction of the distribution range of each species to a single refugial area at different times during the late Pleistocene and a subsequent population expansion seem to be an explanation for the lack of phylogeographical structure across most of Eurasia in species that are ecologically associated with the boreal forest.     

Demographic spatial genetic structure of the Neotropical tree, Jacaranda copaia

We used genotypes from six microsatellite loci and demographic data from a large mapped forest plot to study changes in spatial genetic structure across demographic stages, from seed rain to seedlings, juveniles, and adult diameter classes in the Neotropical tree, Jacaranda copaia. In pairwise comparisons of genetic differentiation among demographic classes, only seedlings were significantly differentiated from the other diameter classes; F(ST) values ranged from 0.006 to 0.009. Furthermore, only seedlings showed homozygote excess suggesting biparental inbreeding in the large diameter reproductive adults. We found very low levels of relatedness in the first distance class of trees, 1-26 cm diameter (F(ij) = 0.011). However, there was a 5- to 10-fold rise in relatedness in the smallest distance class, from the smallest to the largest tree diameter classes (F(ij) = 0.110 for individuals > 56 cm diameter). A variety of non-mutually exclusive mechanisms have been invoked perviously to explain such a pattern, including natural selection, history, or nonequilibrium population dynamics. The long-term demographic data available for this species allow us to evaluate these mechanisms. Jacaranda is a fast-growing, light-demanding species with low recruitment rates and high mortality rates in the smaller diameter classes. It successfully regenerates only in large light gaps, which occur infrequently and stochastically in space and time. These factors contribute to the nonequilibrium population dynamics and observed low genetic structure in the small size classes. We conclude that the pattern of spatial genetic transitions in Jacaranda is consistent with overlapping related generations and strong but infrequent periods of high recruitment, followed by long periods of population decline.     

Genetic diversity analysis of Arius manillensis (Siluriformes: Ariidae) using the mitochondrial control region

Arius manillensis is a Philippine endemic species and is an economically important fishery resource in Laguna de Bay, the largest lake in the country. Drastic reduction in population sizes of A. manillensis has been recorded in the past, which may have resulted in genetic bottleneck. In this study, the genetic diversity and population structure of A. manillensis in Laguna de Bay were assessed using the mitochondrial DNA control region. Specimens were obtained from three localities along Laguna de Bay, namely Binangonan (n?=?27), Tanay (n?=?29), and Calamba (n?=?30). Of the 86 DNA sequences generated, 22 distinct haplotypes were observed. There were four unique haplotypes for Binangonan, six for Calamba, and five for Tanay. There were two haplotypes common to the three sites. The maximum likelihood tree and median-joining network showed little geographic separation among the haplotypes. Chi-square test showed no significant differentiation in A. manillensis from the three sites. The overall computed F(ST) was 0.0144, indicating small genetic differentiation in A. manillensis from the three localities sampled. Likewise, analysis of molecular variance showed a greater percentage of variation within population (98.62%) than variation among populations (1.38%; P?=?0.21). Total haplotype diversity and nucleotide diversity among the specimens from the three sites were 0.775 and 0.013, respectively. The high haplotype diversity coupled with low nucleotide diversity observed in this study confirms that genetic bottleneck occurred in A. manillensis which was followed by population expansion. This is also supported by the non-significant values for both Tajima's D and Fu's F. Furthermore, multimodal mismatch distribution plots were generated, which is consistent with the model of spatial range expansion followed by demographic expansion.     

Recent habitat fragmentation due to roads can lead to significant genetic differentiation in an abundant flightless ground beetle

Although habitat fragmentation is suspected to pose a major threat to biodiversity, its impact on abundant invertebrate species has been little investigated. We assessed the genetic population structure of the flightless ground beetle Abax parallelepipedus in a forest fragmented by two main roads and a highway using five microsatellite loci. We detected low levels of genetic differentiation, which was concordant with the high population densities of 632-1707 individuals/ha estimated with a mark-recapture method. A Mantel test detected a highly significant increase of pairwise F(ST)-values with the number of roads between sampling locations. As expected, the most pronounced effect of the isolation due to roads was observed in the sample from the smallest fragment (highway exit loop), which was differentiated significantly from most other locations. However, no signs of a recent bottleneck or a loss of genetic variability were detected in this population, indicating a still relatively large effective population size (N(e)). Computer simulations confirmed that the observed F(ST)-values were indeed compatible with a N(e) of a few hundred individuals in this fragment, assuming strong or absolute isolation since the construction of the roads. We discuss the implications of our findings for the conservation of abundant but poorly dispersing species in fragmented habitats.     

Interglacial refugia preserved high genetic diversity of the Chinese mole shrew in the mountains of southwest China

The mountains of southwest China (MSC) harbor extremely high species diversity; however, the mechanism behind this diversity is unknown. We investigated to what degree the topography and climate change shaped the genetic diversity and diversification in these mountains, and we also sought to identify the locations of microrefugia areas in these mountains. For these purposes, we sampled extensively to estimate the intraspecific phylogenetic pattern of the Chinese mole shrew (Anourosorex squamipes) in southwest China throughout its range of distribution. Two mitochondrial genes, namely, cytochrome b (CYT B) and NADH dehydrogenase subunit 2 (ND2), from 383 archived specimens from 43 localities were determined for phylogeographic and demographic analyses. We used the continuous-diffusion phylogeographic model, extensive Bayesian skyline plot species distribution modeling (SDM) and approximate Bayesian computation (ABC) to explore the changes in population size and distribution through time of the species. Two phylogenetic clades were identified, and significantly higher genetic diversity was preserved in the southern subregion of the mountains. The results of the SDM, continuous-diffusion phylogeographic model, extensive Bayesian skyline plot and ABC analyses were congruent and supported that the Last Interglacial Maximum (LIG) was an unfavorable period for the mole shrews because of a high degree of seasonality; A. squamipes survived in isolated interglacial refugia mainly located in the southern subregion during the LIG and rapidly expanded during the last glacial period. These results furnished the first evidence for major Pleistocene interglacial refugia and a latitudinal effect in southwest China, and the results shedding light on the higher level of species richness in the southern subregion.     

Experimental demonstration of a causal relationship between heterogeneity of selection and genetic differentiation in quantitative traits

Comparisons of estimates of genetic differentiation at molecular markers (F(ST)) and at quantitative traits (Q(ST)) are a means of inferring the level and heterogeneity of selection in natural populations. However, such comparisons are questionable because they require that the influence of drift and selection on Q(ST) be detectable over possible background influences of environmental or nonadditive genetic effects on Q(ST)-values. Here we test this using an experimental evolution approach in metapopulations of Arabidopsis thaliana experiencing different levels of drift and selection heterogeneity. We estimated the intensity and heterogeneity of selection on morphological and phenological traits via selection differentials. We demonstrate that Q(ST)-values increased with increasing selection heterogeneity when genetic drift was limited. The effect of selection on Q(ST) was thus detectable despite significant genotype-by-environment interactions that most probably biased the estimates of genetic differentiation. Although they cannot be used as a direct validation of the conclusions of prior studies, our results strongly support both the relevance of Q(ST) as an estimator of genetic differentiation and the role of local selection in shaping the genetic differentiation of natural populations.     

Chloroplast DNA phylogeography of European ashes, Fraxinus sp. (Oleaceae): roles of hybridization and life history traits

We investigated range-wide phylogeographic variation in three European ash species (Fraxinus sp., Oleaceae). Chloroplast DNA (cpDNA) microsatellites were typed in the thermophilous Fraxinus angustifolia and Fraxinus ornus and the observed haplotypes and the geographic distribution of diversity were compared to cpDNA data previously obtained in the more cold-tolerant Fraxinus excelsior. We found wide-ranging haplotype sharing between the phylogenetically close F. angustifolia and F. excelsior, suggesting hybridization (i) in common glacial refuges in the Iberian Peninsula, northern Italy, the eastern and/or Dinaric Alps and the Balkan Peninsula, and/or (ii) during postglacial recolonization. The data allowed us to propose additional glacial refuges for F. angustifolia in southern Italy and in Turkey, and populations from the latter region were particularly polymorphic. There was evidence for refuge areas in Italy, the Balkan Peninsula and Turkey for F. ornus, which did not share any single chloroplast haplotype with the other species. In both F. angustifolia and F. ornus, cpDNA diversity (h(S) = 0.027 and h(S) = 0.009, respectively) was lower and fixation levels (G(ST) = 0.964 and G(ST) = 0.983, respectively) higher than in sympatric F. excelsior (h(S) = 0.096, G(ST) = 0.870). These diversity patterns could be due to temperature tolerance or the demographic history.     

Selective sweeps reveal candidate genes for adaptation to drought and salt tolerance in common sunflower, Helianthus annuus

Here we report the results of an analysis of variation at 128 EST-based microsatellites in wild Helianthus annuus, using populations from the species' typical plains habitat in Kansas and Colorado, as well as two arid desert and two distinct brackish marsh areas in Utah. The test statistics lnRV and lnRH were used to find regions of the genome that were significantly less variable in one population relative to the others and thus are likely to contain genes under selection. A small but detectable percentage (1.5-6%) of genes showed evidence for selection from both statistics in any particular environment, and a total of 17 loci showed evidence of selection in at least one environment. Distance-based measures provided additional evidence of selection for 15 of the 17 loci. Global F(ST)-values were significantly higher for candidate loci, as expected under divergent selection. However, pairwise F(ST)-values were lower for populations that shared a selective sweep. Moreover, while spatially separated populations undergoing similar selective pressures showed evidence of divergence at some loci, they evolved in concert at other loci. Thus, this study illustrates how selective sweeps might contribute both to the integration of conspecific populations and to the differentiation of races or species.     

Phylogeographic insights into an irruptive pest outbreak

Irruptive forest insect pests cause considerable ecological and economic damage, and their outbreaks have been increasing in frequency and severity. We use a phylogeographic approach to understand the location and progression of an outbreak by the MPB (Dendroctonus ponderosae Hopkins), an irruptive bark beetle that has caused unprecedented damage to lodgepole pine forests in western North America and is poised to expand its range across the boreal forest. We sampled MPB populations across British Columbia and Alberta and used phylogeographic methods to describe lineage diversification, characterize population structure, investigate expansion dynamics, and identify source populations of the outbreak. Using 1181 bp of mitochondrial DNA sequence from 267 individuals, we found high haplotype diversity, low nucleotide diversity, and limited lineage diversification. The overall pattern was consistent with isolation by distance at a continental scale, and with reduced diversity and population structure in the northerly, outbreak regions. Post-Pleistocene expansion was detected, however more recent expansion signals were not detected, potentially due to the size and rapid rate of range expansion. Based on the limited genetic structure, there were likely multiple source populations in southern British Columbia, although the magnitude of the demographic expansion and rate of spread have obscured the signature of these source populations. Our data highlight the need for caution in interpreting phylogeographic results for species with similar demographics.