Response of Two Heat Shock Genes to Selection for Knockdown Heat Resistance in Drosophila Melanogaster

To identify genes involved in stress resistance and heat hardening, replicate lines of Drosophila melanogaster were selected for increased resistance to knockdown by a 39° heat stress. Two selective regimes were used, one with and one without prior hardening. Mean knockdown times were increased from ~5 min to >20 min after 18 generations. Initial realized heritabilities were as high as 10% for lines selected without hardening, and crosses between lines indicated simple additive gene effects for the selected phenotypes. To survey allelic variation and correlated selection responses in two candidate stress genes, hsr-omega and hsp68, we applied denaturing gradient gel electrophoresis to amplified DNA sequences from small regions of these genes. After eight generations of selection, allele frequencies at both loci showed correlated responses for selection following hardening, but not without hardening. The hardening process itself was associated with a hsp68 frequency change in the opposite direction to that associated with selection that followed hardening. These stress loci are closely linked on chromosome III, and the hardening selection established a disequilibrium, suggesting an epistatic effect on resistance. The data indicate that molecular variation in both hsr-omega and hsp68 contribute to natural heritable variation for hardened heat resistance.     

The latitudinal cline in the In(3R)Payne inversion polymorphism has shifted in the last 20 years in Australian Drosophila melanogaster populations

Clinal variation has been described in a number of inversions in Drosophila but these clines are often characterized by cytological techniques using small sample sizes, and associations with specific genes are rarely considered. Here we have developed a molecular assay for In(3R)Payne in Drosophila melanogaster from eastern Australia populations. It shows in repeated samples that the inversion cline is very tightly associated with latitude and is almost fixed in tropical populations while relatively rare in temperate populations. This steep cline has shifted in position in the last 20 years. The heat shock gene, hsr-omega, located centrally inside the inversion sequence, shows a different clinal pattern to In(3R)Payne. These results suggest strong ongoing selection on In(3R)Payne over the last 100 years since the colonization of Australia that is partly independent of hsr-omega.     

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.     

The Impact of Divergence Time on the Nature of Population Structure: An Example from Iceland

The Icelandic population has been sampled in many disease association studies, providing a strong motivation to understand the structure of this population and its ramifications for disease gene mapping. Previous work using 40 microsatellites showed that the Icelandic population is relatively homogeneous, but exhibits subtle population structure that can bias disease association statistics. Here, we show that regional geographic ancestries of individuals from Iceland can be distinguished using 292,289 autosomal single-nucleotide polymorphisms (SNPs). We further show that subpopulation differences are due to genetic drift since the settlement of Iceland 1100 years ago, and not to varying contributions from different ancestral populations. A consequence of the recent origin of Icelandic population structure is that allele frequency differences follow a null distribution devoid of outliers, so that the risk of false positive associations due to stratification is minimal. Our results highlight an important distinction between population differences attributable to recent drift and those arising from more ancient divergence, which has implications both for association studies and for efforts to detect natural selection using population differentiation.     

Thermal tolerance trade-offs associated with the right arm of chromosome 3 and marked by the hsr-omega gene in Drosophila melanogaster

Drosophila melanogaster occurs in diverse climatic regions and shows opposing clinal changes in resistance to heat and resistance to cold along a 3000 km latitudinal transect on the eastern coast of Australia. We report here on variation at a polymorphic 8 bp-indel site in the heat shock hsr-omega gene that maps to the right arm of chromosome 3. The frequency of the genetic element marked by the L form of the gene was strongly and positively associated with latitude along this transect, and latitudinal differences in L frequency were robustly associated with latitudinal differences in maximum temperature for the hottest month. On a genetic background mixed for genes from each end of the cline a set of 10 lines was derived, five of which were fixed for the L marker, the absence of In(3R)P and 12 kb of repeats at a second polymorphic site at the 3' end of hsr-omega, and five that were fixed for the S marker, In(3R)P and 15 kb of hsr-omega repeats. For two different measures of heat tolerance S lines outperformed L lines, and for two different measures of cold tolerance L lines outperformed S lines. These data suggest that an element on the right arm of chromosome 3, possibly In(3R)P, confers heat resistance but carries the trade-off of also conferring susceptibility to cold. This element occurs at high frequency near the equator. The alternate element on the other hand, at high frequency at temperate latitudes, confers cold resistance at the cost of heat susceptibility.     

High resolution mapping of candidate alleles for desiccation resistance in Drosophila melanogaster under selection

The ability to counter periods of low humidity is an important determinant of distribution range in Drosophila. Climate specialists with low physiological tolerance to desiccation stress are restricted to the tropics and may lack the ability to further increase resistance through evolution. Although the physiological adaptations to desiccation stress are well studied in Drosophila and other ectotherms, factors underlying evolutionary responses remain unknown because of a paucity of genetic data. We address this issue by mapping evolutionary shifts in D. melanogaster under selection for desiccation resistance. Genomic DNA from five independent replicate selected, and control lines were hybridized to high density Affymetrix Drosophila tiling arrays resulting in the detection of 691 single feature polymorphisms (SFPs) differing between the treatments. While randomly distributed throughout the genome, the SFPs formed specific clusters according to gene ontology. These included genes involved in ion transport and respiratory system development that provide candidates for evolutionary changes involving excretory and respiratory water balance. Changes to genes related to neuronal control of cell signaling, development, and gene regulation provide candidates to explore novel biological processes in stress resistance. Sequencing revealed the nucleotide shifts in a subset of the SFPs and highlighted larger regions of genomic diversity surrounding SFPs. The association between natural desiccation resistance and a 463-bp region of the 5' promoter region of the Dys gene undergoing allele frequency changes in response to selection in the experimental evolution lines was tested in an independent population from Coffs Harbour, Australia. The allele frequencies of 23 SNPs common to the two populations were inferred from the parents of the 10% most and 10% least resistant Coffs Harbour flies. The frequencies of the selected alleles were higher at all sites, with three sites significantly associated with the resistant Coffs Harbour flies. This study illustrates how rapid mapping can be used for discovering natural molecular variants associated with survival to low humidity and provides a wealth of candidate alleles to explore the genetic basis of physiological differences among resistant and susceptible Drosophila populations and species.     

Genetic Divergence under Uniform Selection. II. Different Responses to Selection for Knockdown Resistance to Ethanol among DROSOPHILA MELANOGASTER Populations and Their Replicate Lines

We have tested the hypothesis that genetic differences among conspecific populations may result in diverse responses to selection, using natural populations of Drosophila melanogaster. Selection for ethanol tolerance in a tube measuring knockdown resistance was imposed on five West Coast populations. In 24 generations the selected lines increased their mean knockdown times, on average, by a factor of 2.7. An initially weak latitudinal cline was steepened by selection. The two southernmost populations showed the same increases in the selected character, but differed consistently in their correlated responses in characters related to ethanol tolerance. This result indicates that the populations responded to selection by different genetic changes. Selection decreased female body weight and increased resistance to acetone, suggesting components of the response unrelated to ethanol metabolism. The Adhs allele was favored by selection in all populations at the onset, but increased in frequency only in the selected lines of the southernmost population. There was a correlation between latitude and Adh frequency changes, suggesting that fitnesses of the Adh alleles were dependent on the genetic background. Genetic background also had a large effect on the loss of fitness due to selection. Genetic drift between replicate lines caused more variation in selection response than initial genetic differences between populations. This result demonstrates the importance of genetic drift in divergence among natural populations undergoing uniform selection, since the effective population sizes approached those of small natural populations. Drift caused greater divergence between selected replicates than control replicates. Implications of this result for the genetic model of selection response are discussed.     

The genetic basis of adaptive pigmentation variation in Drosophila melanogaster

In a broad survey of Drosophila melanogaster population samples, levels of abdominal pigmentation were found to be highly variable and geographically differentiated. A strong positive correlation was found between dark pigmentation and high altitude, suggesting adaptation to specific environments. DNA sequence polymorphism at the candidate gene ebony revealed a clear association with the pigmentation of homozygous third chromosome lines. The darkest lines sequenced had nearly identical haplotypes spanning 14.5 kb upstream of the protein-coding exons of ebony. Thus, natural selection may have elevated the frequency of an allele that confers dark abdominal pigmentation by influencing the regulation of ebony.     

Worldwide population differentiation at disease-associated SNPs

Background

Recent genome-wide association (GWA) studies have provided compelling evidence of association between genetic variants and common complex diseases. These studies have made use of cases and controls almost exclusively from populations of European ancestry and little is known about the frequency of risk alleles in other populations. The present study addresses the transferability of disease associations across human populations by examining levels of population differentiation at disease-associated single nucleotide polymorphisms (SNPs).

Methods

We genotyped ~1000 individuals from 53 populations worldwide at 25 SNPs which show robust association with 6 complex human diseases (Crohn's disease, type 1 diabetes, type 2 diabetes, rheumatoid arthritis, coronary artery disease and obesity). Allele frequency differences between populations for these SNPs were measured using Fst. The Fst values for the disease-associated SNPs were compared to Fst values from 2750 random SNPs typed in the same set of individuals.

Results

On average, disease SNPs are not significantly more differentiated between populations than random SNPs in the genome. Risk allele frequencies, however, do show substantial variation across human populations and may contribute to differences in disease prevalence between populations. We demonstrate that, in some cases, risk allele frequency differences are unusually high compared to random SNPs and may be due to the action of local (i.e. geographically-restricted) positive natural selection. Moreover, some risk alleles were absent or fixed in a population, which implies that risk alleles identified in one population do not necessarily account for disease prevalence in all human populations.

Conclusion

Although differences in risk allele frequencies between human populations are not unusually large and are thus likely not due to positive local selection, there is substantial variation in risk allele frequencies between populations which may account for differences in disease prevalence between human populations.     

A new test for detecting ongoing selection

Much effort has been made to search for signatures of past natural selection in DNA sequences. However, currently acting selection is rarely detected in natural populations because of its rarity, low detection power of available methods, or both. Here, we develop a new test to detect viability selection over a single generation. In this test, one specific type of chromosomes is chosen as a reference, while all other chromosomes are designated as "focal". The test compares measures of variation between two groups of "focal" chromosomes: those found in reference/focal heterozygous individuals and those found in focal/focal homozygous individuals. In the absence of selection, we do not expect differences between these two groups as long as mating is random. On the other hand, currently acting selection can cause differences in some measures of variation. We applied this test to typing data for In(2L)t inversion polymorphism in a Drosophila melanogaster population, using "standard" (non-inverted) chromosomes as the focal class. Although the frequencies of In(2L)t and standard chromosomes did not deviate from the Hardy-Weinberg equilibrium, we found differences in allele frequency and the number of haplotypes between the two groups of standard chromosomes. This new test, in conjunction with the Hardy-Weinberg test, may shed light on how often strong selection is operating in extant populations.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

An Experimental Test of Adaptive Introgression in Locally Adapted Populations of Splash Pool Copepods

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT₅₀) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.     

The perturbation ofDrosophila melanogaster esterase-6 allozyme frequencies by selection for malathion resistance

The effects of selection for resistance to the organophosphate insecticide malathion on esterase-6 polymorphism was studied in laboratory populations ofD. melanogaster. A genetically well-mixed population was constructed from 40 locally-caught, iso-female lines, divided into control and malathion-selected lines and the frequency of the majorEst-6 alleles was followed for more than 100 generations. The main findings were: (1) The allele frequency in control replicates remained stable for over one hundred generations; (2) the allele frequency in the populations exposed to malathion changed dramatically and in the opposite direction between replicates during the early generations of the selection experiment; (3) all selected populations eventually returned to the controlEst-6 allele frequencies. This return was more rapid in populations exposed to lower selection intensity compared to those exposed to higher selection intensity; (4) the convergence ofEst-6 allele frequencies to control values was also observed in three populations obtained by mixing flies of appropriate genotypes from the control population to give different initial frequencies. These results have been interpreted to mean that esterase-6 does not have a direct role in malathion resistance, and that theEst-6 polymorphism in our experimental population was maintained by balancing selection.     

S-element insertions are associated with the evolution of the Hsp70 genes in Drosophila melanogaster

The "selfish DNA" theory postulates that transposable elements (TEs) are intragenomic parasites, and that natural selection against deleterious effects associated with their presence is the main force preventing their genomic spread in natural populations. In agreement with this model, TEs in Drosophila melanogaster populations are usually found at low frequencies in most genomic locations. Only a few cases of fixation of TE insertions have been reported, usually in regions of low recombination, where selection is expected to be less effective. Here, we report a population genetics study on the apparent fixation of an S-element in a highly recombining region in two natural populations of D. melanogaster. Three similar fragments of an S-element are inserted into the 5' regions of three members of a heat shock gene family, Hsp70 (Hsp70Aa and Hsp70Ab in polytene chromosome band 87A, and Hsp70Bb in 87C). A PCR-based analysis suggests that the insertions are fixed or at high frequencies in the entire species. A population survey of the levels of nucleotide sequence variation at the insertion site in 87C in two natural populations of D. melanogaster provided evidence for reduced levels of variation in the region, normal levels of recombination, and selection, reflected in a significant departure from neutrality of the variant frequency spectrum. This was particularly strong for the S-element inverted repeats (IRs) and suggests that these are of functional significance for the host.     

Physiological stress and the maintenance of adaptive genetic variation in a livebearing fish

The importance of genetic variation in evolution is well established. Yet, the mechanisms by which genetic variation—particularly variation in traits under selection—is maintained in natural populations has long been an evolutionary puzzle. Understanding individual variables driving selection and their functional mechanisms is increasingly important in the context of global change and its potential consequences for biodiversity. Here we examined intra-population performance among allelic variants of a pleiotropic locus in response to thermal stress in the variable platyfish, Xiphophorus variatus. The wild-type tailspot allele exhibited significantly lower heat tolerance than all three pattern alleles found in the population, conforming to predictions based on previously observed correlations between temperature and pattern frequencies in the wild. Furthermore, differences between tailspot pattern frequencies in adults and juveniles were broadly consistent with this trend. Thus, it appears that physiological stress and reduced performance of the wild-type allele at higher relative temperatures is a mechanism balancing its frequency in natural populations. Temperature variation and not dissolved oxygen alone, as previously reported, is likely a important abiotic variable contributing to the maintenance of adaptive polymorphism. Furthermore, our findings underscore the potential implications of rising temperatures and physiological stress for levels of genetic variation in natural populations.     

Dieldrin and diazinon resistance in populations of the Australian sheep blowfly, Lucilia cuprina, from sheep-grazing areas and rubbish tips

Populations of L. cuprina collected from adjacent sheep-grazing areas and rubbish tips in Victoria (Mansfield and Warrnambool) and New South Wales (Lismore) were tested for resistance to the insecticides diazinon and dieldrin. Populations from sheep-grazing areas had a significantly higher diazinon Rop-1 allele frequency than those from adjacent tips with the Victorian populations being more resistant than those from Lismore. Victorian sheep and tip populations had similar gene frequencies at the dieldrin resistance locus, but the Rdl allele frequency was significantly greater in the population at the tip than in the population from the sheep-grazing area at Lismore. The Rdl allele is at a higher frequency in flies from the Lismore area than in Victorian populations. The results at both loci are explained by a balance of selection and gene flow between sheep and tip populations and by selective differences between geographical areas. The exceptionally high frequency of the dieldrin Rdl allele in populations at the Lismore tip may be partially explained by the use of dichlorvos for fly control. Dosage mortality curve and genetic analyses suggest that dichlorvos (an organophosphorus compound) may select at the dieldrin resistance locus. Possible mechanisms for this are discussed. The consequences of genetic differentiation between L. cuprina populations within a region for an autocidal control program are considered.     

ALLOZYMIC DIFFERENTIATION IN RESPONSE TO LABORATORY DEMOGRAPHIC SELECTION OF DROSOPHILA MELANOGASTER

Drosophila melanogaster populations that exhibit constrasting life histories as a result of laboratory selection were compared at several potentially relevant enzyme loci. Selection regimes included postponed reproduction, accelerated development, and intermediate generation time. Each selection regime was represented by fivefold replicated populations maintained for between 50 and 500 generations. For each population, allele frequencies were calculated from frequencies of electrophoretically distinguishable allozymes of alcohol dehydrogenase, α-glycerol-3-phosphate dehydrogenase, phosphoglucomutase, and CuZn-superoxide dismutase. Based on allozyme frequency changes consistent across replicate populations, two of the studied loci responded to both selection for postponed reproduction and selection for accelerated development. The responses to contrasting selection regimes were in opposing directions, suggesting antagonistic pleiotropy.