POTENTIAL FITNESS TRADE-OFFS FOR THERMAL TOLERANCE IN THE INTERTIDAL COPEPOD TIGRIOPUS CALIFORNICUS

Evolution of Ecological Niche Breadth

How ecological niche breadth evolves is central to adaptation and speciation and has been a topic of perennial interest. Niche breadth evolution research has occurred within environmental, ecological, evolutionary, and biogeographical contexts, and although some generalities have emerged, critical knowledge gaps exist. Performance breadth trade-offs, although long invoked, may not be common determinants of niche breadth evolution or limits. Niche breadth can expand or contract from specialist or generalist lineages, and so specialization need not be an evolutionary dead end. Whether niche breadth determines diversification and distribution breadth and how niche breadth is partitioned among individuals and populations within a species are important but particularly understudied topics. Molecular genetic and phylogenetic techniques have greatly expanded understanding of niche breadth evolution, but field studies of how niche breadth evolves are essential for providing mechanistic details and allowing the development of comprehensive theory and improved prediction of biological responses under global change.

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.

Does increased heat resistance result in higher susceptibility to predation? A test using Drosophila melanogaster selection and hardening

Heat resistance of ectotherms can be increased both by plasticity and evolution, but these effects may have trade-offs resulting from biotic interactions. Here, we test for predation costs in Drosophila melanogaster populations with altered heat resistance produced by adult hardening and directional selection for increased heat resistance. In addition, we also tested for genetic trade-offs by testing heat resistance in lines that have evolved under increased predation risk. We show that while 35/37 °C hardening increases heat resistance as expected, it does not increase predation risk from jumping spiders or mantids; in fact, there was an indication that survival may have increased under predation following a triple 37 °C compared to a single 35 °C hardening treatment. Flies that survived a 39 °C selection cycle showed lower survival under predation, suggesting a predation cost of exposure to a more severe heat stress. There was, however, no correlated response to selection because survival did not differ between control and selected lines after selection was relaxed for one or two generations. In addition, lines selected for increased predation risk did not differ in heat resistance. Our findings suggest independent evolutionary responses to predation and heat as measured in laboratory assays, and no costs of heat hardening on susceptibility to predation.

Evolutionary and environmental determinants of freshwater fish thermal tolerance and plasticity

Understanding the extent to which phylogenetic constraints and adaptive evolutionary forces help define the physiological sensitivity of species is critical for anticipating climate-related impacts in aquatic environments. Yet, whether upper thermal tolerance and plasticity are shaped by common evolutionary and environmental mechanisms remains to be tested. Based on a systematic literature review, we investigated this question in 82 freshwater fish species (27 families) representing 829 experiments for which data existed on upper thermal limits and it was possible to estimate plasticity using upper thermal tolerance reaction norms. Our findings indicated that there are strong phylogenetic signals in both thermal tolerances and acclimation capacity, although it is weaker in the latter. We found that upper thermal tolerances are correlated with the temperatures experienced by species across their range, likely because of spatially autocorrelated processes in which closely related species share similar selection pressures and limited dispersal from ancestral environments. No association with species thermal habitat was found for acclimation capacity. Instead, species with the lowest physiological plasticity also displayed the highest thermal tolerances, reflecting to some extent an evolutionary trade-off between these two traits. Although our study demonstrates that macroecological climatic niche features measured from species distributions are likely to provide a good approximation of freshwater fish sensitivity to climate change, disentangling the mechanisms underlying both acute and chronic heat tolerances may help to refine predictions regarding climate change-related range shifts and extinctions.

Adaptation to climate change: trade-offs among responses to multiple stressors in an intertidal crustacean

Trade-offs may influence both physiological and evolutionary responses to co-occurring stressors, but their effects on both plastic and adaptive responses to climate change are poorly understood. To test for genetic and physiological trade-offs incurred in tolerating multiple stressors, we hybridized two populations of the intertidal copepod Tigriopus californicus that were divergent for both heat and salinity tolerance. Starting in the F₂ generation, we selected for increased tolerance of heat, low salinity, and high salinity in replicate lines. After five generations of selection, heat-selected lines had greater heat tolerance but lower fecundity, indicating an energetic cost to tolerance. Lines selected for increased salinity tolerance did not show evidence of adaptation to their respective environments; however, hypo-osmotic selection lines showed substantial loss of tolerance to hyperosmotic stress. Neither of the salinity selection regimes resulted in diminished heat tolerance at ambient salinity; however, simultaneous exposure to heat and hypo-osmotic stress led to decreased heat tolerance, implying a physiological trade-off in tolerance to the two stressors. When we quantified the transcriptomic response to heat and salinity stress via RNA sequencing, we observed little overlap in the stress responses, suggesting the observed synergistic effects of heat and salinity stress were driven by competing energetic demands, rather than shared stress response pathways.

Variation in tolerance to common marine pollutants among different populations in two species of the marine copepod Tigriopus

Geographical variation in chemical tolerance within a species can significantly influence results of whole animal bioassays, yet a literature survey showed that the majority of articles using bioassays did not provide detail on the original field collection site of their test specimens confounding the ability for accurate replication and comparison of results. Biological variation as a result of population-specific tolerance, if not addressed, can be misinterpreted as experimental error. Our studies of two marine copepod species, Tigriopus japonicus and Tigriopus californicus, found significant intra- and inter-specific variation in tolerance to copper and tributyltin. Because both species tolerate copper concentrations orders of magnitude higher than those found in coastal waters, difference in copper tolerance may be a by-product of adaptation to other stressors such as high temperature. Controlling for inter-population tolerance variation will greatly strengthen the application of bioassays in chemical toxicity tests.

Within-site variation of growth rates and terminal sizes in Mytilus californianus along wave exposure and tidal gradients

Mytilus californianus is a foundation species of rocky shores of western North America. Its dominance depends on rapid growth to large sizes, which confers an advantage in size-dependent species interactions. Initial rates of growth and final (terminal) sizes of the mussels depend on environmental factors. Prior comparisons of growth made over large spatial scales (tens of meters to hundreds of kilometers) indicate that temperature, submergence time, and wave exposure affect growth. However, there are few studies quantifying variation in temperature, wave force, and mussel growth parameters at small scales within local populations-that is, meter-level increments. Such measures are necessary to better understand the consequences of the complex spatial mosaic of physical factors in the intertidal zone. We measured variation in temperature, wave force, size-specific shell growth, and terminal size at 3-4-m intervals along horizontal contours within two mussel beds separated by 15 ds of latitude. Both mussel beds showed the same general trends: growth rates attenuated along gradual clines from low and wave-exposed to high shore and sheltered. For example, young adults from low and wave-exposed microhabitats grew 9- and 6-fold higher than those from high-shore-wave-sheltered points. While higher flow may promote growth by enhancing feeding, it also appears to exert a positive effect by moderating energetically costly temperature stress. Consistent with the growth rate findings, cumulative degree-hours explained 83% and 69% of the variation of terminal sizes in regressions for the two locations.

Interpopulational variation in the cold tolerance of a broadly distributed marine copepod

Latitudinal trends in cold tolerance have been observed in many terrestrial ectotherms, but few studies have investigated interpopulational variation in the cold physiology of marine invertebrates. Here, the intertidal copepod Tigriopus californicus was used as a model system to study how local adaptation influences the cold tolerance of a broadly distributed marine crustacean. Among five populations spanning 18° in latitude, the following three metrics were used to compare cold tolerance: the temperature of chill-coma onset, the chill-coma recovery time and post-freezing recovery. In comparison to copepods from warmer southern latitudes, animals from northern populations exhibited lower chill-coma onset temperatures, shorter chill-coma recovery times and faster post-freezing recovery rates. Importantly, all three metrics showed a consistent latitudinal trend, suggesting that any single metric could be used equivalently in future studies investigating latitudinal variation in cold tolerance. Our results agree with previous studies showing that populations within a single species can display strong local adaptation to spatially varying climatic conditions. Thus, accounting for local adaptation in bioclimate models will be useful for understanding how broadly distributed species like T. californicus will respond to anthropogenic climate change.

Acclimation and adaptation to common marine pollutants in the copepod Tigriopus californicus

Establishing water quality criteria using bioassays is complicated by variation in chemical tolerance between populations. Two major contributors to this variation are acclimation and adaptation, which are both linked to exposure history, but differ in how long their effects are maintained. Our study examines how tolerance changes over multiple generations of exposure to two common marine pollutants, copper (Cu) and tributyltin oxide (TBTO), in a sexually reproducing marine copepod, Tigriopus californicus. Lines of T. californicus were chronically exposed to sub-lethal levels of Cu and TBTO for 12 generations followed by a recovery period of 3 generations in seawater control conditions. At each generation, the average number of offspring produced and survived to 28 d was determined and used as the metric of tolerance. Lines exposed to Cu and TBTO showed an overall increase in tolerance over time. Increased Cu tolerance arose by generation 3 in the chronically exposed lines and was lost after 3 generations in seawater control conditions. Increased TBTO tolerance was detected at generation 7 and was maintained even after 3 generations in seawater control conditions. It was concluded from this study that tolerance to Cu is consistent with acclimation, a quick gain and loss of tolerance. In contrast, TBTO tolerance is consistent with adaptation, in which onset of tolerance was delayed relative to an acclimation response and maintained in the absence of exposure. These findings illustrate that consideration of exposure history is necessary when using bioassays to measure chemical tolerance.

Local climatic adaptation in a widespread microorganism

Exploring the ability of organisms to locally adapt is critical for determining the outcome of rapid climate changes, yet few studies have addressed this question in microorganisms. We investigated the role of a heterogeneous climate on adaptation of North American populations of the wild yeast Saccharomyces paradoxus. We found abundant among-strain variation for fitness components across a range of temperatures, but this variation was only partially explained by climatic variation in the distribution area. Most of fitness variation was explained by the divergence of genetically distinct groups, distributed along a north-south cline, suggesting that these groups have adapted to distinct climatic conditions. Within-group fitness components were correlated with climatic conditions, illustrating that even ubiquitous microorganisms locally adapt and harbour standing genetic variation for climate-related traits. Our results suggest that global climatic changes could lead to adaptation to new conditions within groups, or changes in their geographical distributions.

Climate change in the oceans: evolutionary versus phenotypically plastic responses of marine animals and plants

I summarize marine studies on plastic versus adaptive responses to global change. Due to the lack of time series, this review focuses largely on the potential for adaptive evolution in marine animals and plants. The approaches were mainly synchronic comparisons of phenotypically divergent populations, substituting spatial contrasts in temperature or CO2 environments for temporal changes, or in assessments of adaptive genetic diversity within populations for traits important under global change. The available literature is biased towards gastropods, crustaceans, cnidarians and macroalgae. Focal traits were mostly environmental tolerances, which correspond to phenotypic buffering, a plasticity type that maintains a functional phenotype despite external disturbance. Almost all studies address coastal species that are already today exposed to fluctuations in temperature, pH and oxygen levels. Recommendations for future research include (i) initiation and analyses of observational and experimental temporal studies encompassing diverse phenotypic traits (including diapausing cues, dispersal traits, reproductive timing, morphology) (ii) quantification of nongenetic trans-generational effects along with components of additive genetic variance (iii) adaptive changes in microbe-host associations under the holobiont model in response to global change (iv) evolution of plasticity patterns under increasingly fluctuating environments and extreme conditions and (v) joint consideration of demography and evolutionary adaptation in evolutionary rescue approaches.

The relative importance of predation risk and water temperature in maintaining Bergmann's rule in a marine ectotherm

Bergmann's rule-an increase in body size with latitude-correlates with latitudinal declines in ambient temperature and predation risk, but relatively few studies simultaneously explore the relative importance of these factors. Along temperate Atlantic shorelines, the isopod Idotea balthica from high latitudes are 53% longer on average than isopods from low latitudes. When reared at 6°-24°C, juveniles increased growth and development rates with temperature. Because the increase in growth rate with temperature outstripped increases in development rate, female size at maturity increased with temperature. This thermal sensitivity of growth cannot account for the latitudinal pattern in body size. Within temperature treatments, females from low latitudes reached sexual maturity at younger ages and at a smaller size than did females from higher latitudes. This shift in life-history strategy is predicted by latitudinal declines in predation pressure, which we tested using field-tethering experiments. Overall, isopods at low latitudes had a 44% greater mortality risk from daytime predators relative to isopods at higher latitudes. We conclude that a latitudinal gradient in predation risk, not temperature, is principally responsible for Bergmann's rule in I. balthica. Increases in body size during future warming of oceans may be constrained by local patterns of predation risk.

Postzygotic isolation involves strong mitochondrial and sex-specific effects in Tigriopus californicus, a species lacking heteromorphic sex chromosomes

Detailed studies of the genetics of speciation have focused on a few model systems, particularly Drosophila. The copepod Tigriopus californicus offers an alternative that differs from standard animal models in that it lacks heteromorphic chromosomes (instead, sex determination is polygenic) and has reduced opportunities for sexual conflict, because females mate only once. Quantitative trait loci (QTL) mapping was conducted on reciprocal F2 hybrids between two strongly differentiated populations, using a saturated linkage map spanning all 12 autosomes and the mitochondrion. By comparing sexes, a possible sex ratio distorter was found but no sex chromosomes. Although studies of standard models often find an excess of hybrid male sterility factors, we found no QTL for sterility and multiple QTL for hybrid viability (indicated by non-Mendelian adult ratios) and other characters. Viability problems were found to be stronger in males, but the usual explanations for weaker hybrid males (sex chromosomes, sensitivity of spermatogenesis, sexual selection) cannot fully account for these male viability problems. Instead, higher metabolic rates may amplify deleterious effects in males. Although many studies of standard speciation models find the strongest genetic incompatibilities to be nuclear-nuclear (specifically X chromosome-autosome), we found the strongest deleterious interaction in this system was mito-nuclear. Consistent with the snowball theory of incompatibility accumulation, we found that trigenic interactions in this highly divergent cross were substantially more frequent (>6×) than digenic interactions. This alternative system thus allows important comparisons to studies of the genetics of reproductive isolation in more standard model systems.

Trade-offs, geography, and limits to thermal adaptation in a tide pool copepod

Antagonistic correlations among traits may slow the rate of adaptation to a changing environment. The tide pool copepod Tigriopus californicus is locally adapted to temperature, but within populations, the response to selection for increased heat tolerance plateaus rapidly, suggesting either limited variation within populations or costs of increased tolerance. To measure possible costs of thermal tolerance, we selected for increased upper lethal limits for 10 generations in 22 lines of T. californicus from six populations. Then, for each line, we measured six fitness-related traits. Selected lines showed an overall increase in male and female body sizes, fecundity, and starvation resistance, suggesting a small benefit from (rather than costs of) increased tolerance. The effect of selection on correlated traits also varied significantly by population for five traits, indicating that the genetic basis for the selection response differed among populations. Our results suggest that adaptation was limited by the presence of variation within isolated populations rather than by costs of increased tolerance.

Multivariate analysis of adaptive capacity for upper thermal limits in Drosophila simulans

Thermal tolerance is an important factor influencing the distribution of ectotherms, but our understanding of the ability of species to evolve different thermal limits is limited. Based on univariate measures of adaptive capacity, it has recently been suggested that species may have limited evolutionary potential to extend their upper thermal limits under ramping temperature conditions that better reflect heat stress in nature. To test these findings more broadly, we used a paternal half-sibling breeding design to estimate the multivariate evolutionary potential for upper thermal limits in Drosophila simulans. We assessed heat tolerance using static (basal and hardened) and ramping assays. Our analyses revealed significant evolutionary potential for all three measures of heat tolerance. Additive genetic variances were significantly different from zero for all three traits. Our G matrix analysis revealed that all three traits would contribute to a response to selection for increased heat tolerance. Significant additive genetic covariances and additive genetic correlations between static basal and hardened heat-knockdown time, marginally nonsignificant between static basal and ramping heat-knockdown time, indicate that direct and correlated responses to selection for increased upper thermal limits are possible. Thus, combinations of all three traits will contribute to the evolution of upper thermal limits in response to selection imposed by a warming climate. Reliance on univariate estimates of evolutionary potential may not provide accurate insight into the ability of organisms to evolve upper thermal limits in nature.

Evolutionary adaptation of marine zooplankton to global change

Predicting the response of the biota to global change remains a formidable endeavor. Zooplankton face challenges related to global warming, ocean acidification, the proliferation of toxic algal blooms, and increasing pollution, eutrophication, and hypoxia. They can respond to these changes by phenotypic plasticity or genetic adaptation. Using the concept of the evolution of reaction norms, I address how adaptive responses can be unequivocally discerned from phenotypic plasticity. To date, relatively few zooplankton studies have been designed for such a purpose. As case studies, I review the evidence for zooplankton adaptation to toxic algal blooms, hypoxia, and climate change. Predicting the response of zooplankton to global change requires new information to determine (a) the trade-offs and costs of adaptation, (b) the rates of evolution versus environmental change, (c) the consequences of adaptation to stochastic or cyclic (toxic algal blooms, coastal hypoxia) versus directional (temperature, acidification, open ocean hypoxia) environmental change, and (d) the interaction of selective pressures, and evolutionary and ecological processes, in promoting or hindering adaptation.

Investigating the molecular basis of local adaptation to thermal stress: population differences in gene expression across the transcriptome of the copepod Tigriopus californicus

BACKGROUND

Geographic variation in the thermal environment impacts a broad range of biochemical and physiological processes and can be a major selective force leading to local population adaptation. In the intertidal copepod Tigriopus californicus, populations along the coast of California show differences in thermal tolerance that are consistent with adaptation, i.e., southern populations withstand thermal stresses that are lethal to northern populations. To understand the genetic basis of these physiological differences, we use an RNA-seq approach to compare genome-wide patterns of gene expression in two populations known to differ in thermal tolerance.

RESULTS

Observed differences in gene expression between the southern (San Diego) and the northern (Santa Cruz) populations included both the number of affected loci as well as the identity of these loci. However, the most pronounced differences concerned the amplitude of up-regulation of genes producing heat shock proteins (Hsps) and genes involved in ubiquitination and proteolysis. Among the hsp genes, orthologous pairs show markedly different thermal responses as the amplitude of hsp response was greatly elevated in the San Diego population, most notably in members of the hsp70 gene family. There was no evidence of accelerated evolution at the sequence level for hsp genes. Among other sets of genes, cuticle genes were up-regulated in SD but down-regulated in SC, and mitochondrial genes were down-regulated in both populations.

CONCLUSIONS

Marked changes in gene expression were observed in response to acute sub-lethal thermal stress in the copepod T. californicus. Although some qualitative differences were observed between populations, the most pronounced differences involved the magnitude of induction of numerous hsp and ubiquitin genes. These differences in gene expression suggest that evolutionary divergence in the regulatory pathway(s) involved in acute temperature stress may offer at least a partial explanation of population differences in thermal tolerance observed in Tigriopus.

Adaptive capacity of the habitat modifying sea urchin Centrostephanus rodgersii to ocean warming and ocean acidification: performance of early embryos

BACKGROUND

Predicting effects of rapid climate change on populations depends on measuring the effects of climate stressors on performance, and potential for adaptation. Adaptation to stressful climatic conditions requires heritable genetic variance for stress tolerance present in populations.

METHODOLOGY/PRINCIPAL FINDINGS

We quantified genetic variation in tolerance of early development of the ecologically important sea urchin Centrostephanus rodgersii to near-future (2100) ocean conditions projected for the southeast Australian global change hot spot. Multiple dam-sire crosses were used to quantify the interactive effects of warming (+2-4 °C) and acidification (-0.3-0.5 pH units) across twenty-seven family lines. Acidification, but not temperature, decreased the percentage of cleavage stage embryos. In contrast, temperature, but not acidification decreased the percentage of gastrulation. Cleavage success in response to both stressors was strongly affected by sire identity. Sire and dam identity significantly affected gastrulation and both interacted with temperature to determine developmental success. Positive genetic correlations for gastrulation indicated that genotypes that did well at lower pH also did well in higher temperatures.

CONCLUSIONS/SIGNIFICANCE

Significant genotype (sire) by environment interactions for both stressors at gastrulation indicated the presence of heritable variation in thermal tolerance and the ability of embryos to respond to changing environments. The significant influence of dam may be due to maternal provisioning (maternal genotype or environment) and/or offspring genotype. It appears that early development in this ecologically important sea urchin is not constrained in adapting to the multiple stressors of ocean warming and acidification. The presence of tolerant genotypes indicates the potential to adapt to concurrent warming and acidification, contributing to the resilience of C. rodgersii in a changing ocean.

Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus)

Local adaptation is often invoked to explain hybrid zone structure, but empirical evidence of this is generally rare. Hybrid zones between two poeciliid fishes, Xiphophorus birchmanni and X. malinche, occur in multiple tributaries with independent replication of upstream-to-downstream gradients in morphology and allele frequencies. Ecological niche modelling revealed that temperature is a central predictive factor in the spatial distribution of pure parental species and their hybrids and explains spatial and temporal variation in the frequency of neutral genetic markers in hybrid populations. Among populations of parentals and hybrids, both thermal tolerance and heat-shock protein expression vary strongly, indicating that spatial and temporal structure is likely driven by adaptation to local thermal environments. Therefore, hybrid zone structure is strongly influenced by interspecific differences in physiological mechanisms for coping with the thermal environment.

Biophysics, environmental stochasticity, and the evolution of thermal safety margins in intertidal limpets

As the air temperature of the Earth rises, ecological relationships within a community might shift, in part due to differences in the thermal physiology of species. Prediction of these shifts – an urgent task for ecologists – will be complicated if thermal tolerance itself can rapidly evolve. Here, we employ a mechanistic approach to predict the potential for rapid evolution of thermal tolerance in the intertidal limpet Lottia gigantea. Using biophysical principles to predict body temperature as a function of the state of the environment, and an environmental bootstrap procedure to predict how the environment fluctuates through time, we create hypothetical time-series of limpet body temperatures, which are in turn used as a test platform for a mechanistic evolutionary model of thermal tolerance. Our simulations suggest that environmentally driven stochastic variation of L. gigantea body temperature results in rapid evolution of a substantial ‘safety margin’: the average lethal limit is 5–7°C above the average annual maximum temperature. This predicted safety margin approximately matches that found in nature, and once established is sufficient, in our simulations, to allow some limpet populations to survive a drastic, century-long increase in air temperature. By contrast, in the absence of environmental stochasticity, the safety margin is dramatically reduced. We suggest that the risk of exceeding the safety margin, rather than the absolute value of the safety margin, plays an underappreciated role in the evolution of thermal tolerance. Our predictions are based on a simple, hypothetical, allelic model that connects genetics to thermal physiology. To move beyond this simple model – and thereby potentially to predict differential evolution among populations and among species – will require significant advances in our ability to translate the details of thermal histories into physiological and population-genetic consequences.

Considerations for assessing maximum critical temperatures in small ectothermic animals: insights from leaf-cutting ants

The thermal limits of individual animals were originally proposed as a link between animal physiology and thermal ecology. Although this link is valid in theory, the evaluation of physiological tolerances involves some problems that are the focus of this study. One rationale was that heating rates shall influence upper critical limits, so that ecological thermal limits need to consider experimental heating rates. In addition, if thermal limits are not surpassed in experiments, subsequent tests of the same individual should yield similar results or produce evidence of hardening. Finally, several non-controlled variables such as time under experimental conditions and procedures may affect results. To analyze these issues we conducted an integrative study of upper critical temperatures in a single species, the ant Atta sexdens rubropiosa, an animal model providing large numbers of individuals of diverse sizes but similar genetic makeup. Our specific aims were to test the 1) influence of heating rates in the experimental evaluation of upper critical temperature, 2) assumptions of absence of physical damage and reproducibility, and 3) sources of variance often overlooked in the thermal-limits literature; and 4) to introduce some experimental approaches that may help researchers to separate physiological and methodological issues. The upper thermal limits were influenced by both heating rates and body mass. In the latter case, the effect was physiological rather than methodological. The critical temperature decreased during subsequent tests performed on the same individual ants, even one week after the initial test. Accordingly, upper thermal limits may have been overestimated by our (and typical) protocols. Heating rates, body mass, procedures independent of temperature and other variables may affect the estimation of upper critical temperatures. Therefore, based on our data, we offer suggestions to enhance the quality of measurements, and offer recommendations to authors aiming to compile and analyze databases from the literature.

A gene-based SNP resource and linkage map for the copepod Tigriopus californicus

BACKGROUND

As yet, few genomic resources have been developed in crustaceans. This lack is particularly evident in Copepoda, given the extraordinary numerical abundance, and taxonomic and ecological diversity of this group. Tigriopus californicus is ideally suited to serve as a genetic model copepod and has been the subject of extensive work in environmental stress and reproductive isolation. Accordingly, we set out to develop a broadly-useful panel of genetic markers and to construct a linkage map dense enough for quantitative trait locus detection in an interval mapping framework for T. californicus--a first for copepods.

RESULTS

One hundred and ninety Single Nucleotide Polymorphisms (SNPs) were used to genotype our mapping population of 250 F2 larvae. We were able to construct a linkage map with an average intermarker distance of 1.8 cM, and a maximum intermarker distance of 10.3 cM. All markers were assembled into linkage groups, and the 12 linkage groups corresponded to the 12 known chromosomes of T. californicus. We estimate a total genome size of 401.0 cM, and a total coverage of 73.7%. Seventy five percent of the mapped markers were detected in 9 additional populations of T. californicus. Of available model arthropod genomes, we were able to show more colocalized pairs of homologues between T. californicus and the honeybee Apis mellifera, than expected by chance, suggesting preserved macrosynteny between Hymenoptera and Copepoda.

CONCLUSIONS

Our study provides an abundance of linked markers spanning all chromosomes. Many of these markers are also found in multiple populations of T. californicus, and in two other species in the genus. The genomic resource we have developed will enable mapping throughout the geographical range of this species and in closely related species. This linkage map will facilitate genome sequencing, mapping and assembly in an ecologically and taxonomically interesting group for which genomic resources are currently under development.

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.

Incorporating population-level variation in thermal performance into predictions of geographic range shifts

Determining how species' geographic ranges are governed by current climates and how they will respond to rapid climatic change poses a major biological challenge. Geographic ranges are often spatially fragmented and composed of genetically differentiated populations that are locally adapted to different thermal regimes. Tradeoffs between different aspects of thermal performance, such as between tolerance to high temperature and tolerance to low temperature or between maximal performance and breadth of performance, suggest that the performance of a given population will be a subset of that of the species. Therefore, species-level projections of distribution might overestimate the species' ability to persist at any given location. However, current approaches to modeling distributions often do not consider variation among populations. Here, we estimated genetically-based differences in thermal performance curves for growth among 12 populations of the scarlet monkeyflower, Mimulus cardinalis, a perennial herb of western North America. We inferred the maximum relative growth rate (RGR(max)), temperature optimum (T(opt)), and temperature breadth (T(breadth)) for each population. We used these data to test for tradeoffs in thermal performance, generate mechanistic population-level projections of distribution under current and future climates, and examine how variation in aspects of thermal performance influences forecasts of range shifts. Populations differed significantly in RGR(max) and had variable, but overlapping, estimates of T(opt) and T(breadth). T(opt) declined with latitude and increased with temperature of origin, consistent with tradeoffs between performances at low temperatures versus those at high temperatures. Further, T(breadth) was negatively related to RGR(max), as expected for a specialist-generalist tradeoff. Parameters of the thermal performance curve influenced properties of projected distributions. For both current and future climates, T(opt) was negatively related to latitudinal position, while T(breadth) was positively related to projected range size. The magnitude and direction of range shifts also varied with T(opt) and T(breadth), but sometimes in unexpected ways. For example, the fraction of habitat remaining suitable increased with T(opt) but decreased with T(breadth). Northern limits of all populations were projected to shift north, but the magnitude of shift decreased with T(opt) and increased with T(breadth). Median latitude was projected to shift north for populations with high T(breadth) and low T(opt), but south for populations with low T(breadth) and high T(opt). Distributions inferred by integrating population-level projections did not differ from a species-level projection that ignored variation among populations. However, the species-level approach masked the potential array of divergent responses by populations that might lead to genotypic sorting within the species' range. Thermal performance tradeoffs among populations within the species' range had important, but sometimes counterintuitive, effects on projected responses to climatic change.

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₂ 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.

Interpopulation patterns of divergence and selection across the transcriptome of the copepod Tigriopus californicus

The accumulation of genetic incompatibilities between isolated populations is thought to lead to the evolution of intrinsic postzygotic isolation. The molecular basis for these mechanisms, however, remains poorly understood. The intertidal copepod Tigriopus californicus provides unique opportunities for addressing mechanistic questions regarding the early stages of speciation; hybrids between highly divergent populations are fertile and viable, but exhibit reduced fitness at the F(2) or later generations. Given the current scarcity of genomic information in taxa at incipient stages of reproductive isolation, we utilize high-throughout 454 pyrosequencing to characterize a substantial fraction of protein-coding regions (the transcriptome) of T. californicus. Our sequencing effort was divided equally between two divergent populations in order to estimate levels of divergence and to reveal patterns of selection across the transcriptome. Assembly of sequences generated over 40,000 putatively unique transcripts (unigenes) for each population, 19,622 of which were orthologous between populations. BLAST searches of public databases determined protein identity and functional features for 15,402 and 12,670 unigenes, respectively. Based on rates of nonsynonymous and synonymous substitutions in 5897 interpopulation orthologs (those >150 bp and with at least 2X coverage), we identified 229 potential targets of positive selection. Many of these genes are predicted to be involved in several metabolic processes, and to function in hydrolase, peptidase and binding activities. The library of T. californicus coding regions, annotated with their predicted functions and level of divergence, will serve as an invaluable resource for elucidating molecular mechanisms underlying the early stages of speciation.

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.