Host species differences in the thermal mismatch of host-parasitoid interactions.

Extreme high temperatures associated with climate change can affect species directly, and indirectly through temperature-mediated species interactions. In most host-parasitoid systems, parasitization inevitably kills the host, but differences in heat tolerance between host and parasitoid, and between different hosts, may alter their interactions. Here, we explored the effects of extreme high temperatures on the ecological outcomes - including, in some rare cases, escape from the developmental disruption of parasitism - of the parasitoid wasp, Cotesia congregata, and two co-occurring congeneric larval hosts, Manduca sexta and M. quinquemaculata. Both host species had higher thermal tolerance than C. congregata, resulting in a thermal mismatch characterized by parasitoid (but not host) mortality under extreme high temperatures. Despite parasitoid death at high temperatures, hosts typically remain developmentally disrupted from parasitism. However, high temperatures resulted in a partial developmental recovery from parasitism (reaching the wandering stage at the end of host larval development) in some host individuals, with a significantly higher frequency of this partial developmental recovery in M. quinquemaculata than in M. sexta. Hosts species also differed in their growth and development in the absence of parasitoids, with M. quinquemaculata developing faster and larger at high temperatures relative to M. sexta. Our results demonstrate that co-occurring congeneric species, despite shared environments and phylogenetic histories, can vary in their responses to temperature, parasitism and their interaction, resulting in altered ecological outcomes.

Hypoxia tolerance, but not low pH tolerance, is associated with a latitudinal cline across populations of Tigriopus californicus.

Intertidal organisms must tolerate daily fluctuations in environmental parameters, and repeated exposure to co-occurring conditions may result in tolerance to multiple stressors correlating. The intertidal copepod Tigriopus californicus experiences diurnal variation in dissolved oxygen levels and pH as the opposing processes of photosynthesis and cellular respiration lead to coordinated highs during the day and lows at night. While environmental parameters with overlapping spatial gradients frequently result in correlated traits, less attention has been given to exploring temporally correlated stressors. We investigated whether hypoxia tolerance correlates with low pH tolerance by separately testing the hypoxia and low pH stress tolerance separately of 6 genetically differentiated populations of T. californicus. We independently checked for similarities in tolerance for each of the two stressors by latitude, sex, size, and time since collection as predictors. We found that although hypoxia tolerance correlated with latitude, low pH tolerance did not, and no predictor was significant for both stressors. We concluded that temporally coordinated exposure to low pH and low oxygen did not result in populations developing equivalent tolerance for both. Although climate change alters several environmental variables simultaneously, organisms' abilities to tolerate these changes may not be similarly coupled.

Frequent Paternal Mitochondrial Inheritance and Rapid Haplotype Frequency Shifts in Copepod Hybrids.

Mitochondria are assumed to be maternally inherited in most animal species, and this foundational concept has fostered advances in phylogenetics, conservation, and population genetics. Like other animals, mitochondria were thought to be solely maternally inherited in the marine copepod Tigriopus californicus, which has served as a useful model for studying mitonuclear interactions, hybrid breakdown, and environmental tolerance. However, we present PCR, Sanger sequencing, and Illumina Nextera sequencing evidence that extensive paternal mitochondrial DNA (mtDNA) transmission is occurring in inter-population hybrids of T. californicus. PCR on four types of crosses between three populations (total sample size of 376 F1 individuals) with 20% genome-wide mitochondrial divergence showed 2% to 59% of F1 hybrids with both paternal and maternal mtDNA, where low and high paternal leakage values were found in different cross directions of the same population pairs. Sequencing methods further verified nucleotide similarities between F1 mtDNA and paternal mtDNA sequences. Interestingly, the paternal mtDNA in F1s from some crosses inherited haplotypes that were uncommon in the paternal population. Compared to some previous research on paternal leakage, we employed more rigorous methods to rule out contamination and false detection of paternal mtDNA due to non-functional nuclear mitochondrial DNA fragments. Our results raise the potential that other animal systems thought to only inherit maternal mitochondria may also have paternal leakage, which would then affect the interpretation of past and future population genetics or phylogenetic studies that rely on mitochondria as uniparental markers.

Tolerance Patterns and Transcriptomic Response to Extreme and Fluctuating Salinities across Populations of the Intertidal Copepod Tigriopus californicus.

AbstractPopulations that tolerate extreme environmental conditions with frequent fluctuations can give valuable insights into physiological limits and adaptation. In some estuarine and marine ecosystems, organisms must adapt to extreme and fluctuating salinities, but not much is known about how varying salinities impact local adaptation across a wide geographic range. We used eight geographically and genetically divergent populations of the intertidal copepod Tigriopus californicus to test whether northern populations have greater tolerance to low salinity stresses, as they experience greater precipitation and less evaporation. We used a common-garden experiment approach and exposed all populations to acute low (1 and 3 ppt) and high (110 and 130 ppt) salinities for 24 h and to a fluctuation between baseline salinity and moderate low (7 ppt) and high (80 ppt) salinities for over 49 h. We also performed RNA sequencing at several time points during the fluctuation between baseline and salinity of 7 ppt to understand the molecular basis of divergence between two populations with differing physiological responses. We present these novel findings: (1) acute low salinity conditions caused more deaths than high salinity; (2) molecular processes that elevate proline levels increased in salinity of 7 ppt, which contrasts with other physiological studies in T. californicus that mainly associated accumulation of proline with hyperosmotic stress; and (3) tolerance to a salinity fluctuation did not follow a latitudinal trend but was instead governed by a complex interplay of factors, including population and duration of salinity stress. This highlights the importance of including a wider variety of environmental conditions in empirical studies to understand local adaptation.

The ghost of temperature past: interactive effects of previous and current thermal conditions on gene expression in Manduca sexta.

High temperatures can negatively impact the performance and survival of organisms, particularly ectotherms. While an organism's response to high temperature stress clearly depends on current thermal conditions, its response may also be affected by the temporal pattern and duration of past temperature exposures. We used RNA sequencing of Manduca sexta larvae fat body tissue to evaluate how diurnal temperature fluctuations during development affected gene expression both independently and in conjunction with subsequent heat stress. Additionally, we compared gene expression between two M. sexta populations, a lab colony and a genetically related field population that have been separated for >300 generations and differ in their thermal sensitivities. Lab-adapted larvae were predicted to show increased expression responses to both single and repeated thermal stress, whereas recurrent exposure could decrease later stress responses for field individuals. We found large differences in overall gene expression patterns between the two populations across all treatments, as well as population-specific transcriptomic responses to temperature; more differentially expressed genes were upregulated in the field compared with lab larvae. Developmental temperature fluctuations alone had minimal effects on long-term gene expression patterns, with the exception of a somewhat elevated stress response in the lab population. Fluctuating rearing conditions did alter gene expression during exposure to later heat stress, but this effect depended on both the population and the particular temperature conditions. This study contributes to increased knowledge of molecular mechanisms underlying physiological responses of organisms to temperature fluctuations, which is needed for the development of more accurate thermal performance models.

Genomic scans reveal multiple mito-nuclear incompatibilities in population crosses of the copepod Tigriopus californicus.

The evolution of intrinsic postzygotic isolation can be explained by the accumulation of Dobzhansky-Muller incompatibilities (DMI). Asymmetries in the levels of hybrid inviability and hybrid sterility are commonly observed between reciprocal crosses, a pattern that can result from the involvement of uniparentally inherited factors. The mitochondrial genome is one such factor that appears to participate in DMI in some crosses but the frequency of its involvement versus biparentally inherited factors is unclear. Here we assess the relative importance of incompatibilities between nuclear factors (nuclear-nuclear) versus those between mitochondrial and nuclear factors (mito-nuclear) in a species that lacks sex chromosomes. We used a Pool-seq approach to survey three crosses among genetically divergent populations of the copepod, Tigriopus californicus, for regions of the genome that are affected by hybrid inviability. Results from reciprocal crosses suggest that mito-nuclear incompatibilities are more common than nuclear-nuclear incompatibilities overall. These results suggest that in the presence of very high levels of nucleotide divergence between mtDNA haplotypes, mito-nuclear incompatibilities can be important for the evolution of intrinsic postzygotic isolation. This is particularly interesting considering this species lacks sex chromosomes, which have been shown to harbor a particularly high number of nuclear-nuclear DMI in several other species.

Genomic signatures of mitonuclear coevolution across populations of Tigriopus californicus.

The copepod Tigriopus californicus shows extensive population divergence and is becoming a model for understanding allopatric differentiation and the early stages of speciation. Here, we report a high-quality reference genome for one population (~190 megabases across 12 scaffolds, and ~15,500 protein-coding genes). Comparison with other arthropods reveals 2,526 genes presumed to be specific to T. californicus, with an apparent proliferation of genes involved in ion transport and receptor activity. Beyond the reference population, we report re-sequenced genomes of seven additional populations, spanning the continuum of reproductive isolation. Populations show extreme mitochondrial DNA divergence, with higher levels of amino acid differentiation than observed in other taxa. Across the nuclear genome, we find elevated protein evolutionary rates and positive selection in genes predicted to interact with mitochondrial DNA and the proteins and RNA it encodes in multiple pathways. Together, these results support the hypothesis that rapid mitochondrial evolution drives compensatory nuclear evolution within isolated populations, thereby providing a potentially important mechanism for causing intrinsic reproductive isolation.

Evolution of Melanoma Antigen-A11 (MAGEA11) During Primate Phylogeny.

Melanoma antigen-A11 (MAGE-A11) is an X-linked and primate-specific steroid hormone receptor transcriptional coregulator and proto-oncogenic protein whose increased expression promotes the growth of prostate cancer. The MAGEA11 gene is expressed at low levels in normal human testis, ovary, and endometrium, and at highest levels in castration-resistant prostate cancer. Annotated genome predictions throughout the surviving primate lineage show that MAGEA11 acquired three 5' coding exons unique within the MAGEA subfamily during the evolution of New World monkeys (NWM), Old World monkeys (OWM), and apes. MAGE-A11 in all primates has a conserved FXXIF coactivator-binding motif that suggests interaction with p160 coactivators contributed to its early evolution as a transcriptional coregulator. An ancestral form of MAGE-A11 in the more distantly related lemur has significant amino acid sequence identity with human MAGE-A11, but lacks coregulator activity based on the absence of the three 5' coding exons that include a nuclear localization signal (NLS). NWM MAGE-A11 has greater amino acid sequence identity than lemur to human MAGE-A11, but inframe premature stop codons suggest that MAGEA11 is a pseudogene in NWM. MAGE-A11 in OWM and apes has nearly identical 5' coding exon amino acid sequence and conserved interaction sites for p300 acetyltransferase and cyclin A. We conclude that the evolution of MAGEA11 within the lineage leading to OWM and apes resulted in steroid hormone receptor transcriptional coregulator activity through the acquisition of three 5' coding exons that include a NLS sequence and nonsynonymous substitutions required to interact with cell cycle regulatory proteins and transcription factors.

Using Pool-seq to Search for Genomic Regions Affected by Hybrid Inviability in the copepod T. californicus.

The formation of reproductive barriers between allopatric populations involves the accumulation of incompatibilities that lead to intrinsic postzygotic isolation. The evolution of these incompatibilities is usually explained by the Dobzhansky-Muller model, where epistatic interactions that arise within the diverging populations, lead to deleterious interactions when they come together in a hybrid genome. These incompatibilities can lead to hybrid inviability, killing individuals with certain genotypic combinations, and causing the population's allele frequency to deviate from Mendelian expectations. Traditionally, hybrid inviability loci have been detected by genotyping individuals at different loci across the genome. However, this method becomes time consuming and expensive as the number of markers or individuals increases. Here, we test if a Pool-seq method can be used to scan the genome of F2 hybrids to detect genomic regions that are affected by hybrid inviability. We survey the genome of hybrids between 2 populations of the copepod Tigriopus californicus, and show that this method has enough power to detect even small changes in allele frequency caused by hybrid inviability. We show that allele frequency estimates in Pool-seq can be affected by the sampling of alleles from the pool of DNA during the library preparation and sequencing steps and that special considerations must be taken when aligning hybrid reads to a reference when the populations/species are divergent.

Locally adapted populations of a copepod can evolve different gene expression patterns under the same environmental pressures.

As populations diverge in allopatry, but under similar thermal conditions, do similar thermal performance phenotypes evolve by maintaining similar gene expression patterns, or does genetic divergence lead to divergent patterns of gene expression between these populations? We used genetically divergent populations of the copepod Tigriopus californicus, whose performance at different thermal conditions is well characterized, to investigate transcriptome-wide expression responses under two different thermal regimes: (1) a nonvariable temperature regime and (2) a regime with variable temperature. Our results show the expression profiles of the response to these regimes differed substantially among populations, even for populations that are geographically close. This pattern was accentuated when populations were raised in the variable temperature environment. Less heat-tolerant populations mounted strong but divergent responses to the different thermal regimes, with a large heat-shock response observed in one population, and an apparent reduction in the expression of genes involved in basic cellular processes in the other. Our results suggest that as populations diverge in allopatry, they may evolve starkly different responses to changes in temperature, at the gene expression level, while maintaining similar thermal performance phenotypes.

Can we predict ectotherm responses to climate change using thermal performance curves and body temperatures?

Thermal performance curves (TPCs), which quantify how an ectotherm's body temperature (Tb ) affects its performance or fitness, are often used in an attempt to predict organismal responses to climate change. Here, we examine the key - but often biologically unreasonable - assumptions underlying this approach; for example, that physiology and thermal regimes are invariant over ontogeny, space and time, and also that TPCs are independent of previously experienced Tb. We show how a critical consideration of these assumptions can lead to biologically useful hypotheses and experimental designs. For example, rather than assuming that TPCs are fixed during ontogeny, one can measure TPCs for each major life stage and incorporate these into stage-specific ecological models to reveal the life stage most likely to be vulnerable to climate change. Our overall goal is to explicitly examine the assumptions underlying the integration of TPCs with Tb , to develop a framework within which empiricists can place their work within these limitations, and to facilitate the application of thermal physiology to understanding the biological implications of climate change.

Chromosome-Wide Impacts on the Expression of Incompatibilities in Hybrids of Tigriopus californicus.

Chromosome rearrangements such as inversions have been recognized previously as contributing to reproductive isolation by maintaining alleles together that jointly contribute to deleterious genetic interactions and postzygotic reproductive isolation. In this study, an impact of potential incompatibilities merely residing on the same chromosome was found in crosses of populations of the copepod Tigriopus californicus When genetically divergent populations of this copepod are crossed, hybrids show reduced fitness, and deviations from expected genotypic ratios can be used to determine regions of the genome involved in deleterious interactions. In this study, a set of markers was genotyped for a cross of two populations of T. californicus, and these markers show widespread deviations from Mendelian expectations, with entire chromosomes showing marked skew. Despite the importance of mtDNA/nuclear interactions in incompatibilities in this system in previous studies, in these crosses the expected patterns stemming from these interactions are not widely apparent. Females lack recombination in this species, and a striking difference is observed between male and female backcrosses. This suggests that the maintenance of multiple loci on individual chromosomes can enable some incompatibilities, perhaps playing a similar role in the initial rounds of hybridization to chromosomal rearrangements in preserving sets of alleles together that contribute to incompatibilities. Finally, it was observed that candidate pairs of incompatibility regions are not consistently interacting across replicates or subsets of these crosses, despite the repeatability of the deviations at many of the single loci themselves, suggesting that more complicated models of Dobzhansky-Muller incompatibilities may need to be considered.

In memoriam: Richard G. Harrison - his life and legacy.

Richard G. Harrison passed away unexpectedly on April 12th, 2016. In this memoriam we pay tribute to the life and legacy of an extraordinary scientist, mentor, friend, husband, and father.

Gene conversion yields novel gene combinations in paralogs of GOT1 in the copepod Tigriopus californicus.

BACKGROUND: Gene conversion of duplicated genes can slow the divergence of paralogous copies over time but can also result in other interesting evolutionary patterns. Islands of genetic divergence that persist in the face of gene conversion can point to gene regions undergoing selection for new functions. Novel combinations of genetic variation that differ greatly from the original sequence can result from the transfer of genetic variation between paralogous genes by rare gene conversion events. Genetically divergent populations of the copepod Tigriopus californicus provide an excellent model to look at the patterns of divergence among paralogs across multiple independent evolutionary lineages. RESULTS: In this study the evolution of a set of paralogous genes encoding putative aspartate transaminase proteins (called GOT1 here) are examined in populations of the copepod T. californicus. One pair of duplicated genes, GOT1p1 and GOT1p2, has regions of high divergence between the copies in the face of apparent on-going gene conversion. The GOT1p2 gene also has unique haplotypes in two populations that appear to have resulted from a transfer of genetic variation via inter-paralog gene conversion. A second pair of duplicated genes GOT1Sr and GOT1Sd also shows evidence of gene conversion, but this gene conversion does not appear to have maintained each as a functional copy in all populations. CONCLUSIONS: The patterns of conservation and sequence divergence across this set of paralogous genes among populations of T. californicus suggest that some interesting evolutionary patterns are occurring at these loci. The results for the GOT1p1/GOT1p2 paralogs illustrate how gene conversion can factor in the creation of a mosaic pattern of regions of high divergence and low divergence. When coupled with rare gene conversion events of divergent regions, this pattern can result in the formation of novel proteins differing substantially from either original protein. The evolutionary patterns across these paralogs show how gene conversion can both constrain and facilitate diversification of genetic sequences.

Quantifying the elevation of mitochondrial DNA evolutionary substitution rates over nuclear rates in the intertidal copepod Tigriopus californicus.

Mitochondrial DNA (mtDNA) genomes generally evolve rapidly in animals, but considerable variation in the rates of evolution of mtDNA occurs among taxa. Higher levels of mutation will tend to increase the amount of polymorphism, which should also scale with population size, but there are mixed signals from previous studies on the evolutionary outcomes of the interactions of these processes. The copepod Tigriopus californicus provides an interesting model in which to study the evolution of mtDNA because it has high levels of divergence among populations and there is the suggestion that this divergence could be involved in reproductive isolation. This species also appears to have an elevated mtDNA substitution rate, but previous studies did not provide an accurate measurement. This article examines the rate of mtDNA substitution versus nuclear substitution in T. californicus and finds that the mtDNA rate for synonymous sites averages 55-fold higher, a level that exceeds the rates found in most other invertebrates. Levels of polymorphism are also examined in both mtDNA and nuclear genes, and it is shown that the effective population size of mtDNA genes is much lower than that of nuclear genes. In addition, no correlation between polymorphism in mtDNA and nuclear genes is found across populations, which suggest factors other than demography may shape polymorphism in this species. The results from this study suggest that mtDNA is evolving at a very rapid rate in this copepod species, and this could increase the likelihood that mtDNA evolution is involved in the generation of reproductive isolation.

Hybrid breakdown weakens under thermal stress in population crosses of the copepod Tigriopus californicus.

The outcome of hybridization can be impacted by environmental conditions, which themselves can contribute to reproductive isolation between taxa. In crosses of genetically divergent populations, hybridization can have both negative and positive impacts on fitness, the balance between which might be tipped by changes in the environment. Genetically divergent populations of the intertidal copepod Tigriopus californicus have been shown to differ in thermal tolerance at high temperatures along a latitudinal gradient. In this study, a series of crosses were made between pairs of genetically divergent populations of T. californicus, and the thermal tolerance of these hybrids was tested. In most cases, the first-generation hybrids had relatively high thermal tolerance and the second-generation hybrids were not generally reduced below the less-tolerant parental population for high temperature tolerance. This pattern contrasts with previous studies from crosses of genetically divergent populations of this copepod, which often shows hybrid breakdown in these second-generation hybrids for other measures of fitness. These results suggest that high temperature stress could either increase the positive impacts of hybridization or decrease the negative impacts of hybridization resulting in lowered hybrid breakdown in these population crosses.

Complex deleterious interactions associated with malic enzyme may contribute to reproductive isolation in the copepod Tigriopus californicus.

Dobzhansky-Muller incompatibilities can result from the interactions of more than a single pair of interacting genes and there are several different models of how such complex interactions can be structured. Previous empirical work has identified complex conspecific epistasis as a form of complex interaction that has contributed to postzygotic reproductive isolation between taxa, but other forms of complexity are also possible. Here, I probe the genetic basis of reproductive isolation in crosses of the intertidal copepod Tigriopus californicus by looking at the impact of markers in genes encoding metabolic enzymes in F(2) hybrids. The region of the genome associated with the locus ME2 is shown to have strong, repeatable impacts on the fitness of hybrids in crosses and epistatic interactions with another chromosomal region marked by the GOT2 locus in one set of crosses. In a cross between one of these populations and a third population, these two regions do not appear to interact despite the continuation of a large effect of the ME2 region itself in both crosses. The combined results suggest that the ME2 chromosomal region is involved in incompatibilities with several unique partners. If these deleterious interactions all stem from the same factor in this region, that would suggest a different form of complexity from complex conspecific epistasis, namely, multiple independent deleterious interactions stemming from the same factor. Confirmation of this idea will require more fine-scale mapping of the interactions of the ME2 region of the genome.

The nature of interactions that contribute to postzygotic reproductive isolation in hybrid copepods.

Deleterious interactions within the genome of hybrids can lower fitness and result in postzygotic reproductive isolation. Understanding the genetic basis of these deleterious interactions, known as Dobzhansky-Muller incompatibilities, is the subject of intense current study that seeks to elucidate the nature of these deleterious interactions. Hybrids from crosses of individuals from genetically divergent populations of the intertidal copepod Tigriopus californicus provide a useful model in which to study Dobzhansky-Muller incompatibilities. Studies of the basis of postzygotic reproductive isolation in this species have revealed a number of patterns. First, there is evidence for a breakdown in genomic coadaptation between mtDNA-encoded and nuclear-encoded proteins that can result in a reduction in hybrid fitness in some crosses. It appears from studies of the individual genes involved in these interactions that although this coadaptation could lead to asymmetries between crosses, patterns of genotypic viabilities are not often consistent with simple models of genomic coadaptation. Second, there is a large impact of environmental factors on these deleterious interactions suggesting that they are not strictly intrinsic in nature. Temperature in particular appears to play an important role in determining the nature of these interactions. Finally, deleterious interactions in these hybrid copepods appear to be complex in terms of the number of genetic factors that interact to lead to reductions in hybrid fitness. This complexity may stem from three or more factors that all interact to cause a single incompatibility or the same factor interacting with multiple other factors independently leading to multiple incompatibilities.

Potential fitness trade-offs for thermal tolerance in the intertidal copepod Tigriopus californicus.

Thermal adaptation to spatially varying environmental conditions occurs in a wide range of species, but what is less clear is the nature of fitness trade-offs associated with this temperature adaptation. Here, populations of the intertidal copepod Tigriopus californicus are examined at both local and latitudinal scales to determine whether these populations have evolved differences in their survival under high temperature stress. A clear pattern of increasing high temperature stress tolerance is seen with decreasing latitude, consistent with temperature adaptation. Additionally, there is also evidence for significant variation in thermal tolerance on a smaller scale. The competitive fitness of pairs of northern and southern copepod populations were also examined under a series of lower, more moderate temperatures. These fitness assays show that the southern populations that have the best survival under extreme high temperatures have lowered competitive fitness at the lower temperatures tested, whereas the fitness of the southern populations exceeded that of the northern populations at the highest temperatures tested. Combined, these results suggest that there may be evolutionary trade-offs between performance at high and stressful temperatures and fitness at moderate temperatures in this species.

Investigations of fine-scale phylogeography in Tigriopus californicus reveal historical patterns of population divergence.

BACKGROUND: The intertidal copepod Tigriopus californicus is a model for studying the process of genetic divergence in allopatry and for probing the nature of genetic changes that lead to reproductive isolation. Although previous studies have revealed a pattern of remarkably high levels of genetic divergence between the populations of this species at several spatial scales, it is not clear what types of historical processes are responsible. Particularly lacking are data that can yield insights into population history from the finest scales of geographic resolution. RESULTS: Sequence variation in both cytochrome b (CYTB, mtDNA) and the rieske iron-sulfur protein (RISP, nuclear) are examined at a fine scale within four different regions for populations of T. californicus. High levels of genetic divergence are seen for both genes at the broader scale, and genetic subdivision is apparent at nearly all scales in these populations for these two genes. Patterns of polymorphism and divergence in both CYTB and RISP suggest that selection may be leading to non-neutral evolution of these genes in several cases but a pervasive pattern of neither selection nor coadaptation is seen for these markers. CONCLUSION: The use of sequence data at a fine-scale of resolution in this species has provided novel insights into the processes that have resulted in the accumulation of genetic divergence among populations. This divergence is likely to result from an interplay between a limited dispersal ability for this copepod and the temporal instability of copepod habitat. Both shorter-term processes such as the extinction/recolonization dynamics of copepod pools and longer-term processes such as geological uplift of coastline and sea level changes appear to have impacted the patterns of differentiation. Some patterns of sequence variation are consistent with selection acting upon the loci used in this study; however, it appears that most phylogeographic patterns are the result of history and not selection on these genes in this species.