Genetic Evidence for the Introduction of Rhagoletis pomonella (Diptera: Tephritidae) into the Northwestern United States.

The apple maggot fly, Rhagoletis pomonella Walsh (Diptera: Tephritidae), is a serious quarantine pest in the apple-growing regions of central Washington and Oregon. The fly is believed to have been introduced into the Pacific Northwest via the transport of larval-infested apples near Portland, Oregon, within the last 40 yr. However, R. pomonella also attacks native black hawthorn, Crataegus douglasii Lindley (Rosales: Rosaceae), and introduced ornamental hawthorn, Crataegus monogyna Jacquin, in the region. It is, therefore, possible that R. pomonella was not introduced but has always been present on black hawthorn. If true, then the fly may have independently shifted from hawthorn onto apple in the Pacific Northwest within the last 40 yr after apples were introduced. Here, we test the introduction hypothesis through a microsatellite genetic survey of 10 R. pomonella sites in Washington and 5 in the eastern United States, as well as a comparison to patterns of genetic variation between populations of Rhagoletis cingulata Loew and Rhagoletis indifferens Curran, two sister species of cherry-infesting flies known to be native to the eastern and western United States, respectively. We report results based on genetic distance networks, patterns of allelic variation, and estimated times of population divergence that are consistent with the introduction hypothesis for R. pomonella. The results have important implications for R. pomonella management, suggesting that black hawthorn-infesting flies near commercial apple-growing regions of central Washington may harbor sufficient variation to utilize apple as an alternate host, urging careful monitoring, and possible removal of hawthorn trees near orchards.

Gene expression in closely related species mirrors local adaptation: consequences for responses to a warming world.

Local adaptation of populations could preclude or slow range expansions in response to changing climate, particularly when dispersal is limited. To investigate the differential responses of populations to changing climatic conditions, we exposed poleward peripheral and central populations of two Lepidoptera to reciprocal, common-garden climatic conditions and compared their whole-transcriptome expression. We found evidence of simple population differentiation in both species, and in the species with previously identified population structure and phenotypic local adaptation, we found several hundred genes that responded in a synchronized and localized fashion. These genes were primarily involved in energy metabolism and oxidative stress, and expression levels were most divergent between populations in the same environment in which we previously detected divergence for metabolism. We found no localized genes in the species with less population structure and for which no local adaptation was previously detected. These results challenge the assumption that species are functionally similar across their ranges and poleward peripheral populations are preadapted to warmer conditions. Rather, some taxa deserve population-level consideration when predicting the effects of climate change because they respond in genetically based, distinctive ways to changing conditions.

Thermal variability increases the impact of autumnal warming and drives metabolic depression in an overwintering butterfly.

Increases in thermal variability elevate metabolic rate due to Jensen's inequality, and increased metabolic rate decreases the fitness of dormant ectotherms by increasing consumption of stored energy reserves. Theory predicts that ectotherms should respond to increased thermal variability by lowering the thermal sensitivity of metabolism, which will reduce the impact of the warm portion of thermal variability. We examined the thermal sensitivity of metabolic rate of overwintering Erynnis propertius (Lepidoptera: Hesperiidae) larvae from a stable or variable environment reared in the laboratory in a reciprocal common garden design, and used these data to model energy use during the winters of 1973-2010 using meteorological data to predict the energetic outcomes of metabolic compensation and phenological shifts. Larvae that experienced variable temperatures had decreased thermal sensitivity of metabolic rate, and were larger than those reared at stable temperatures, which could partially compensate for the increased energetic demands. Even with depressed thermal sensitivity, the variable environment was more energy-demanding than the stable, with the majority of this demand occurring in autumn. Autumn phenology changes thus had disproportionate influence on energy consumption in variable environments, and variable-reared larvae were most susceptible to overwinter energy drain. Therefore the energetic impacts of the timing of entry into winter dormancy will strongly influence ectotherm fitness in northern temperate environments. We conclude that thermal variability drives the expression of metabolic suppression in this species; that phenological shifts will have a greater impact on ectotherms in variable thermal environments; and that E. propertius will be more sensitive to shifts in phenology in autumn than in spring. This suggests that increases in overwinter thermal variability and/or extended, warm autumns, will negatively impact all non-feeding dormant ectotherms which lack the ability to suppress their overwinter metabolic thermal sensitivity.

Population-level transcriptome sequencing of nonmodel organisms Erynnis propertius and Papilio zelicaon.

BACKGROUND: Several recent studies have demonstrated the use of Roche 454 sequencing technology for de novo transcriptome analysis. Low error rates and high coverage also allow for effective SNP discovery and genetic diversity estimates. However, genetically diverse datasets, such as those sourced from natural populations, pose challenges for assembly programs and subsequent analysis. Further, estimating the effectiveness of transcript discovery using Roche 454 transcriptome data is still a difficult task. RESULTS: Using the Roche 454 FLX Titanium platform, we sequenced and assembled larval transcriptomes for two butterfly species: the Propertius duskywing, Erynnis propertius (Lepidoptera: Hesperiidae) and the Anise swallowtail, Papilio zelicaon (Lepidoptera: Papilionidae). The Expressed Sequence Tags (ESTs) generated represent a diverse sample drawn from multiple populations, developmental stages, and stress treatments. Despite this diversity, > 95% of the ESTs assembled into long (> 714 bp on average) and highly covered (> 9.6x on average) contigs. To estimate the effectiveness of transcript discovery, we compared the number of bases in the hit region of unigenes (contigs and singletons) to the length of the best match silkworm (Bombyx mori) protein--this "ortholog hit ratio" gives a close estimate on the amount of the transcript discovered relative to a model lepidopteran genome. For each species, we tested two assembly programs and two parameter sets; although CAP3 is commonly used for such data, the assemblies produced by Celera Assembler with modified parameters were chosen over those produced by CAP3 based on contig and singleton counts as well as ortholog hit ratio analysis. In the final assemblies, 1,413 E. propertius and 1,940 P. zelicaon unigenes had a ratio > 0.8; 2,866 E. propertius and 4,015 P. zelicaon unigenes had a ratio > 0.5. CONCLUSIONS: Ultimately, these assemblies and SNP data will be used to generate microarrays for ecoinformatics examining climate change tolerance of different natural populations. These studies will benefit from high quality assemblies with few singletons (less than 26% of bases for each assembled transcriptome are present in unassembled singleton ESTs) and effective transcript discovery (over 6,500 of our putative orthologs cover at least 50% of the corresponding model silkworm gene).

Translocation experiments with butterflies reveal limits to enhancement of poleward populations under climate change.

There is a pressing need to predict how species will change their geographic ranges under climate change. Projections typically assume that temperature is a primary fitness determinant and that populations near the poleward (and upward) range boundary are preadapted to warming. Thus, poleward, peripheral populations will increase with warming, and these increases facilitate poleward range expansions. We tested the assumption that poleward, peripheral populations are enhanced by warming using 2 butterflies (Erynnis propertius and Papilio zelicaon) that co-occur and have contrasting degrees of host specialization and interpopulation genetic differentiation. We performed a reciprocal translocation experiment between central and poleward, peripheral populations in the field and simulated a translocation experiment that included alternate host plants. We found that the performance of both central and peripheral populations of E. propertius were enhanced during the summer months by temperatures characteristic of the range center but that local adaptation of peripheral populations to winter conditions near the range edge could counteract that enhancement. Further, poleward range expansion in this species is prevented by a lack of host plants. In P. zelicaon, the fitness of central and peripheral populations decreased under extreme summer temperatures that occurred in the field at the range center. Performance in this species also was affected by an interaction of temperature and host plant such that host species strongly mediated the fitness of peripheral individuals under differing simulated temperatures. Altogether we have evidence that facilitation of poleward range shifts through enhancement of peripheral populations is unlikely in either study species.

The response of two butterfly species to climatic variation at the edge of their range and the implications for poleward range shifts.

To predict changes in species' distributions due to climate change we must understand populations at the poleward edge of species' ranges. Ecologists generally expect range shifts under climate change caused by the expansion of edge populations as peripheral conditions increasingly resemble the range core. We tested whether peripheral populations of two contrasting butterflies, a small-bodied specialist (Erynnis propertius) and a large-bodied generalist (Papilio zelicaon), respond favorably to warmer conditions. Performance of populations related to climate was evaluated in seven peripheral populations spanning 1.2 degrees latitude (160 km) using: (1) population density surveys, an indirect measure of site suitability; and (2) organismal fitness in translocation experiments. There was evidence that population density increased with temperature for P. zelicaon whose population density declined with latitude in 1 of 3 sample years. On the other hand, E. propertius showed a positive relationship of population density with latitude, apparently unrelated to climate or measured habitat variables. Translocation experiments showed increased larval production at increased temperatures for both species, and in P. zelicaon, larval production also increased under drier conditions. These findings suggest that both species may increase at their range edge with warming but the preference for core-like conditions may be stronger in P. zelicaon. Further, populations of E. propertius at the range boundary may be large enough to act as sources of colonists for range expansions, but range expansion in this species may be prevented by a lack of available host plants further north. In total, the species appear to respond differently to climate and other factors that vary latitudinally, factors that will likely affect poleward expansion.