Genotype-by-environment and epistatic interactions in Drosophila melanogaster: the effects of Gpdh allozymes, genetic background and rearing temperature on larval developmental time and viability

Local thermal adaptation and local temperature regimes drive the performance of a parasitic helminth under climate change: The case of Marshallagia marshalli from wild ungulates

Across a species' range, populations are exposed to their local thermal environments, which on an evolutionary scale, may cause adaptative differences among populations. Helminths often have broad geographic ranges and temperature-sensitive life stages but little is known about whether and how local thermal adaptation can influence their response to climate change. We studied the thermal responses of the free-living stages of Marshallagia marshalli, a parasitic nematode of wild ungulates, along a latitudinal gradient. We first determine its distribution in wild sheep species in North America. Then we cultured M. marshalli eggs from different locations at temperatures from 5 to 38°C. We fit performance curves based on the metabolic theory of ecology to determine whether development and mortality showed evidence of local thermal adaptation. We used parameter estimates in life-cycle-based host-parasite models to understand how local thermal responses may influence parasite performance under general and location-specific climate-change projections. We found that M. marshalli has a wide latitudinal and host range, infecting wild sheep species from New Mexico to Yukon. Increases in mortality and development time at higher temperatures were most evident for isolates from northern locations. Accounting for location-specific parasite parameters primarily influenced the magnitude of climate change parasite performance, while accounting for location-specific climates primarily influenced the phenology of parasite performance. Despite differences in development and mortality among M. marshalli populations, when using site-specific climate change projections, there was a similar magnitude of impact on the relative performance of M. marshalli among populations. Climate change is predicted to decrease the expected lifetime reproductive output of M. marshalli in all populations while delaying its seasonal peak by approximately 1 month. Our research suggests that accurate projections of the impacts of climate change on broadly distributed species need to consider local adaptations of organisms together with local temperature profiles and climate projections.

The quantitative genetic basis of clinal divergence in phenotypic plasticity

Phenotypic plasticity is thought to be an important mechanism for adapting to environmental heterogeneity. Nonetheless, the genetic basis of plasticity is still not well understood. In Drosophila melanogaster and D. simulans, body size and thermal stress resistance show clinal patterns along the east coast of Australia, and exhibit plastic responses to different developmental temperatures. The genetic basis of thermal plasticity, and whether the genetic effects underlying clinal variation in traits and their plasticity are similar, remains unknown. Here, we use line-cross analyses between a tropical and temperate population of Drosophila melanogaster and D. simulans developed at three constant temperatures (18°C, 25°C, and 29°C) to investigate the quantitative genetic basis of clinal divergence in mean thermal response (elevation) and plasticity (slope and curvature) for thermal stress and body size traits. Generally, the genetic effects underlying divergence in mean response and plasticity differed, suggesting that different genetic models may be required to understand the evolution of trait means and plasticity. Furthermore, our results suggest that nonadditive genetic effects, in particular epistasis, may commonly underlie plastic responses, indicating that current models that ignore epistasis may be insufficient to understand and predict evolutionary responses to environmental change.

Epistasis × environment interactions among Arabidopsis thaliana glucosinolate genes impact complex traits and fitness in the field

Despite the growing number of studies showing that genotype × environment and epistatic interactions control fitness, the influences of epistasis × environment interactions on adaptive trait evolution remain largely uncharacterized. Across three field trials, we quantified aliphatic glucosinolate (GSL) defense chemistry, leaf damage, and relative fitness using mutant lines of Arabidopsis thaliana varying at pairs of causal aliphatic GSL defense genes to test the impact of epistatic and epistasis × environment interactions on adaptive trait variation. We found that aliphatic GSL accumulation was primarily influenced by additive and epistatic genetic variation, leaf damage was primarily influenced by environmental variation and relative fitness was primarily influenced by epistasis and epistasis × environment interactions. Epistasis × environment interactions accounted for up to 48% of the relative fitness variation in the field. At a single field site, the impact of epistasis on relative fitness varied significantly over 2 yr, showing that epistasis × environment interactions within a location can be temporally dynamic. These results suggest that the environmental dependency of epistasis can profoundly influence the response to selection, shaping the adaptive trajectories of natural populations in complex ways, and deserves further consideration in future evolutionary studies.

ASSOCIATION BETWEEN LATITUDINAL VARIATION FOR EMBRYONIC DEVELOPMENT TIME AND CHROMOSOME STRUCTURE IN THE GRASSHOPPER CALEDIA CAPTIVA (ORTHOPTERA: ACRIDIDAE)

From southeastern Queensland to southern Victoria, over a transect of 11° latitude, the Moreton taxon of the Australian grasshopper Caledia captiva exhibits a cline in chromosome structure that involves change from a metacentric to an acrocentric genome. In this study, we show that embryonic development time covaries with chromosome structure along the transect. Both development time and chromosome short arm length exhibit an overall negative correlation with latitude, but with maxima just south of the northern limit of the taxon's distribution. Selection for such a pattern appears to arise from changes in voltinism along the cline in season length that exists along the transect. Populations with the highest temperature thresholds for avoidance of embryonic diapause also have the slowest development time and probably represent the northern extreme of a primarily univoltine life cycle. North of this region bivoltinism increases in frequency and, as expected from a split of the season length, development time decreases. Maximum chromosome short arm length occurs in the vicinity of the northern univoltine populations, rather than at the limit of distribution where bivoltinism prevails. We conclude that variation in chromosome structure could be contributing to the heritable variation for development time that forms the basis for adaptive change in this trait. These results provide justification for investigating causal relationships between chromosome structure and development time, with an ultimate aim of understanding the adaptive significance of chromosomal variation in C. captiva.

Systems genomics of metabolic phenotypes in wild-type Drosophila melanogaster

Systems biology is an approach to dissection of complex traits that explicitly recognizes the impact of genetic, physiological, and environmental interactions in the generation of phenotypic variation. We describe comprehensive transcriptional and metabolic profiling in Drosophila melanogaster across four diets, finding little overlap in modular architecture. Genotype and genotype-by-diet interactions are a major component of transcriptional variation (24 and 5.3% of the total variation, respectively) while there were no main effects of diet (<1%). Genotype was also a major contributor to metabolomic variation (16%), but in contrast to the transcriptome, diet had a large effect (9%) and the interaction effect was minor (2%) for the metabolome. Yet specific principal components of these molecular phenotypes measured in larvae are strongly correlated with particular metabolic syndrome-like phenotypes such as pupal weight, larval sugar content and triglyceride content, development time, and cardiac arrhythmia in adults. The second principal component of the metabolomic profile is especially informative across these traits with glycine identified as a key loading variable. To further relate this physiological variability to genotypic polymorphism, we performed evolve-and-resequence experiments, finding rapid and replicated changes in gene frequency across hundreds of loci that are specific to each diet. Adaptation to diet is thus highly polygenic. However, loci differentially transcribed across diet or previously identified by RNAi knockdown or expression QTL analysis were not the loci responding to dietary selection. Therefore, loci that respond to the selective pressures of diet cannot be readily predicted a priori from functional analyses.

Interactions between Controlled Atmospheres and Low Temperature Tolerance: A Review of Biochemical Mechanisms

Controlled atmosphere treatments using carbon dioxide, oxygen, and/or nitrogen, together with controlled temperature and humidity, form an important method for post-harvest sterilization against insect-infested fruit. However, in insects, the cross tolerance and biochemical interactions between the various stresses of modified gas conditions and low temperature may either elicit or block standard stress responses which can potentiate (or limit) lethal low temperature exposure. Thus, the success of such treatments is sometimes erratic and does not always result in the desired pest mortality. This review focuses on the biochemical modes of action whereby controlled atmospheres affect insects low temperature tolerance, making them more (or occasionally, less) susceptible to cold sterilization. Insights into the integrated biochemical modes of action may be used together with the pests' low temperature tolerance physiology to determine which treatments may be of value in post-harvest sterilization.

Candidate genes and thermal phenotypes: identifying ecologically important genetic variation for thermotolerance in the Australian Drosophila melanogaster cline

Clinal variation in traits often reflects climatic adaptation; in Drosophila melanogaster clinal variation provides an opportunity to link variation in chromosomal inversions, microsatellite loci and various candidate genes to adaptive variation in traits. We undertook association studies with crosses from a single population of D. melanogaster from eastern Australia to investigate the association between genetic markers and traits showing clinal variation. By genotyping parents and phenotyping offspring, we minimized genotyping costs but had the power to detect association between markers and quantitative traits. Consistent with prior studies, we found strong associations between the clinal chromosomal inversion In(3R)Payne and markers within it, as well as among these markers. We also found an association between In(3L)Payne and one marker located within this inversion. Of the five predicted associations between markers and traits, four were detected (increased heat, decreased cold resistance and body size with the heat shock gene hsr-omega S, increased cold resistance with the inversion In(3L)Payne), while one was not detected (heat resistance and the heat shock gene hsp68). In a set of eight exploratory tests, we detected one positive association (between hsp23a and heat resistance) but no associations of heat resistance with alleles at the hsp26, hsp83, Desat 2, alpha-Gpdh, hsp70 loci, while cold resistance was not associated with Frost and Dca loci. These results confirm interactions between hsr-omega and thermal resistance, as well as between In(3L)Payne and cold resistance, but do not provide evidence for associations between thermal responses and alleles at other clinically varying marker genes.

Testing the pluralist approach to sex: the influence of environment on synergistic interactions between mutation load and parasitism in Daphnia magna

Both deleterious mutations and parasites have been acknowledged as potential selective forces responsible for the evolutionary maintenance of sexual reproduction. The pluralist approach to sex proposes that these two factors may have to interact synergistically in order to stabilize sex, and one of the simplest ways this could occur is if parasites are capable of causing synergistic epistasis between mutations in their hosts. However, the effects of both deleterious mutations and parasitism are known to be influenced by a range of environmental factors, so the nature of the interaction may depend upon the organisms' environment. Using chemically mutated Daphnia magna lines, we examined the effects of mutation and parasitism under a range of temperature and food regimes. We found that although parasites were capable of causing synergistic epistasis between mutations in their hosts, these effects were dependent upon an interaction between parasite genotype and temperature.

Triglyceride pools, flight and activity variation at the Gpdh locus in Drosophila melanogaster

We have created a set of P-element excision-derived Gpdh alleles that generate a range of GPDH activity phenotypes ranging from zero to full activity. By placing these synthetic alleles in isogenic backgrounds, we characterize the effects of minor and major activity variation on two different aspects of Gpdh function: the standing triglyceride pool and glycerol-3-phosphate shuttle-assisted flight. We observe small but statistically significant reductions in triglyceride content for adult Gpdh genotypes possessing 33-80% reductions from normal activity. These small differences scale to a notable proportion of the observed genetic variation in triglyceride content in natural populations. Using a tethered fly assay to assess flight metabolism, we observed that genotypes with 100 and 66% activity exhibited no significant difference in wing-beat frequency (WBF), while activity reductions from 60 to 10% showed statistically significant reductions of approximately 7% in WBF. These studies show that the molecular polymorphism associated with GPDH activity could be maintained in natural populations by selection in the triglyceride pool.

Alternative splicing, muscle contraction and intraspecific variation: associations between troponin T transcripts, Ca2+ sensitivity and the force and power output of dragonfly flight muscles during oscillatory contraction

The flight muscles of Libellula pulchella dragonflies contain a mixture of six alternatively spliced transcripts of a single troponin T (TnT) gene. Here, we examine how intraspecific variation in the relative abundance of different TnT transcripts affects the Ca2+ sensitivity of skinned muscle fibers and the performance of intact muscles during work-loop contraction regimes that approximate in vivo conditions during flight. The relative abundance of one TnT transcript, or the pooled relative abundance of two TnT transcripts, showed a positive correlation with a 10-fold range of variation in Ca2+ sensitivity of skinned fibers (r2=0.77, P&lt;0.0001) and a threefold range in peak specific force (r2=0.74, P&lt;0.0001), specific work per cycle (r2=0.54; P&lt;0.0001) and maximum specific power output (r2=0.48, P=0.0005) of intact muscle. Using these results to reanalyze previously published data for wing kinematics during free flight, we show that the relative abundances of these particular transcripts are also positively correlated with wingbeat frequency and amplitude. TnT variation alone may be responsible for these effects, or TnT variation may be a marker for changes in a suite of co-regulated molecules. Dragonflies from two ponds separated by 16 km differed significantly in both TnT transcript composition and muscle contractile performance, and within each population there are two distinct morphs that showed different maturational trajectories of TnT transcript composition and muscle contractility. Thus, there is broad intraspecific variability and a high degree of population structure for contractile performance phenotypes, TnT ribotypes and ontogenetic patterns involving these traits that affect locomotor performance.

Variability in the size, composition, and function of insect flight muscles

In order to fly, insects require flight muscles that constitute at least 12 to 16% of their total mass, and flight performance increases as this percentage increases. However, flight muscles are energetically and materially expensive to build and maintain, and investment in flight muscles constrains other aspects of function, particularly female fecundity. This review examines ways in which insects vary the size of their flight muscles, and how variation in the relative size and composition of flight muscles affects flight performance. Sources of variability in flight muscle size and composition include genetic differences within and between species, individual phenotypic responses to environmental stimuli, and maturational changes that occur before and during the adult stage. Insects have evolved a wide variety of ways to adjust flight muscle size and contractile performance in order to meet demands imposed by variation in life history and ecology.

The evolution of recombination in a heterogeneous environment

Most models describing the evolution of recombination have focused on the case of a single population, implicitly assuming that all individuals are equally likely to mate and that spatial heterogeneity in selection is absent. In these models, the evolution of recombination is driven by linkage disequilibria generated either by epistatic selection or drift. Models based on epistatic selection show that recombination can be favored if epistasis is negative and weak compared to directional selection and if the recombination modifier locus is tightly linked to the selected loci. In this article, we examine the joint effects of spatial heterogeneity in selection and epistasis on the evolution of recombination. In a model with two patches, each subject to different selection regimes, we consider the cases of mutation-selection and migration-selection balance as well as the spread of beneficial alleles. We find that including spatial heterogeneity extends the range of epistasis over which recombination can be favored. Indeed, recombination can be favored without epistasis, with negative and even with positive epistasis depending on environmental circumstances. The selection pressure acting on recombination-modifier loci is often much stronger with spatial heterogeneity, and even loosely linked modifiers and free linkage may evolve. In each case, predicting whether recombination is favored requires knowledge of both the type of environmental heterogeneity and epistasis, as none of these factors alone is sufficient to predict the outcome.

Polymorphism at the Adh and alpha Gpdh loci in Drosophila melanogaster: effects of rearing temperature on developmental rate, body weight, and some biochemical parameters

The role of developmental time in the world-wide cline of Adh and alpha Gpdh allele frequencies of Drosophila melanogaster, and the relationship with weight and some biochemical characters, were investigated. Experimental strains were constructed with different combinations of Adh and alpha Gpdh alleles but with similar genetic background. Developmental time, adult weight, protein-and triglyceride-content, and ADH and alpha GPDH enzyme activity were measured at a rearing temperature of 20, 25 and 29 degrees C. Genotype effects were found in all studied characters. In general the developmental times of genotypes were: AdhFF less than AdhFS less than AdhSS and alpha GpdhFF greater than alpha GpdhFS = alpha GpdhSS. Developmental time and adult weight were strongly affected by rearing temperature. Triglyceride content and ADH and alpha GPDH enzyme activity were slightly affected by temperature. Interactions between genotype and temperature effects were found for developmental rate, adult weight and protein content. No trade off was observed between developmental time on the one hand and adult weight, protein- and triglyceride-content, and ADH and alpha GPDH enzyme activity on the other hand. It is argued that developmental rate differences might be one of the underlying mechanisms of the world-wide cline of the Adh and alpha GPdh allele frequencies.