A proline repeat polymorphism of the Frost gene of Drosophila melanogaster showing clinal variation but not associated with cold resistance.

Genetic polymorphisms underlying adaptive shifts in thermal responses are poorly known even though studies are providing a detailed understanding of these responses at the cellular and physiological levels. The Frost gene of Drosophila melanogaster is a prime candidate for thermal adaptation; it is up-regulated under cold stress and knockdown of this gene influences cold resistance. Here we describe an amino-acid INDEL polymorphism in proline repeat number in the structural component of this gene. The two main repeats, accounting for more than 90% of alleles in eastern Australia, show a strong clinal pattern; the 6P allele was at a high frequency in tropical locations, and the 10P allele was common in temperate populations. However, the frequency of these alleles was not associated with three different assays of cold resistance. Adult transcription level of Frost was also unrelated to cold resistance as measured through post chill coma mobility. The functional significance of the proline repeat polymorphism therefore remains unclear despite its clinal pattern. The data also demonstrate the feasibility of using Roche/454 sequencing for establishing clinal patterns.

Widespread pyrethroid resistance in Australian diamondback moth, Plutella xylostella (L.), is related to multiple mutations in the para sodium channel gene.

Populations of Plutella xylostella, extending over 3800 km in southern Australia, show no genetic structure as assessed by microsatellite markers; yet outbreaks of pyrethroid resistance occur sporadically in cropping areas. Since mutations in the para voltage-gated sodium channel gene have been implicated in pyrethroid resistance, we looked for DNA sequence variation at this target among Australian moths. We found two resistance mutations previously reported for this species (L1014F and T929I), as well as a novel substitution (F1020S). Of the eight possible haplotypes formed by combinations of these three biallelic polymorphisms, only four were found in Australian populations: the wild-type allele (w), the kdr mutation allele (kdr) with only L1014F, the super-kdr-like combination of L1014F and T929I (skdrl), and the crashdown allele with only F1020S (cdr). Comparison of genotype frequencies among survivors of permethrin assays with those from untreated controls identified three resistant genotypes: skdrl homozygotes, cdr homozygotes and the corresponding heterozygote, cdr/skrdl - the heterozygote being at least as resistant as either homozygote. Spatial heterogeneity of allele frequencies was conspicuous, both across the continent and among local collections, consistent with reported spatial heterogeneity of pyrethroid resistance. Further, high resistance samples were sometimes associated with high frequency of cdr, sometimes high frequency of skdrl, or sometimes with a high combined cdr+skdrl frequency. The skdrl and cdr alleles explain a high proportion of the Australia-wide resistance variation. These data add to evidence that nerve insensitivity by mutations in the para-sodium channel gene is a common pyrethroid resistance mechanism in P. xylostella.

A clinally varying promoter polymorphism associated with adaptive variation in wing size in Drosophila.

Body size often shows adaptive clines in many ectotherms across altitude and latitude, but little is known about the genetic basis of these adaptive clines. Here we identify a polymorphism in the Dca (Drosophila cold acclimation) gene in Drosophila melanogaster that influences wing size, affects wing:thorax allometry and also controls a substantial proportion of the clinal wing-size variation. A polymorphism in the promoter region of Dca had two common alleles showing strong reciprocal clinal variation in frequency with latitude along the east coast of Australia. The Dca-237 allele increased towards the tropics where wing size is smaller. A within-population association study highlighted that an increase in the frequency of this allele decreased wing size but did not influence thorax size. A manipulated increase in the level of expression of Dca achieved through UAS-GAL4 was associated with a decrease in wing size but had no effect on thorax size. This was consistent with higher Dca expression levels in family lines with higher frequency of the Dca-237 allele. Genetic variation in the promoter region of the Dca gene appears to influence adaptive size variation in the eastern Australian cline of Drosophila melanogaster and accounts for more than 10% of the genetic variation in size within and between populations.

Selection for cold resistance alters gene transcript levels in Drosophila melanogaster.

Microarrays have been used to examine changes in gene expression underlying responses to selection for increased stress resistance in Drosophila melanogaster, but changes in expression patterns associated with increased resistance to cold stress have not been previously reported. Here we describe such changes in basal expression levels in replicate lines following selection for increased resistance to chill coma stress. We found significant up- or down-regulation of expression in 94 genes on the Affymetrix Genome 2.0 array. Quantitative RT-PCR was used to confirm changes in expression of six genes. Some of the identified genes had previously been associated with stress resistance but no previously identified candidate genes for cold resistance showed altered patterns of expression. Seven differentially expressed genes that form a tight chromosomal cluster and an unlinked gene AnnX may be potentially important for cold adaptation in natural populations. Artificial selection for chill coma resistance therefore altered basal patterns of gene expression, but we failed to link these changes to plastic changes in expression under cold stress or to previously identified candidate genes for components of cold resistance.

Latitudinal and cold-tolerance variation associate with DNA repeat-number variation in the hsr-omega RNA gene of Drosophila melanogaster.

An 8-bp deletion in the hsr-omega heat-stress gene of Drosophila melanogaster has previously been associated with latitude, and with heat tolerance that decreases with latitude. Here we report a second polymorphic site, at the 3'-end of hsr-omega, at which multiple alleles segregate in natural populations for copy number of a approximately 280 bp tandem repeat. On each of 3 consecutive years (2000, 2001 and 2002) among populations sampled along the Australian eastern coast, repeat number was negatively associated with latitude. Neither altitudinal association was detected in 2002 when five high-altitude sites were included, nor was a robust association detected with local temperature or rainfall measures. Although in a large number of family lines, derived from a population located centrally in the latitudinal transect, no association between hsr-omega repeat number and heat tolerance occurred, a negative association of repeat number with cold tolerance was detected. As cold tolerance also exhibits latitudinal clines we examined a set of cold-tolerant populations derived by selection and found both reduced repeat number and low constitutive levels of the omega-n repeat-bearing transcript. In a sample from the central population, linkage disequilibrium was measured between repeat number and linked markers that also cline latitudinally. However, such disequilibrium could not account for the cline in repeat number or tolerance associations. Finally, during adult recovery from cold exposure a large increase occurred in tissue levels of the omega-c transcript. Together these data suggest that a latitudinal cline in hsr-omega repeat number influences cold-tolerance variation in this species.

Robust clines and robust sampling: a reply to Kyriacou et al.

Kyriacou et al. (2007) have questioned a number of issues with our recent paper on a lack of clinal variation in the period and clock timing genes in Drosophila melanogaster from eastern Australia. Here we show why their arguments are not valid and reiterate that clinal variation in genes and molecular markers need to be assessed on field flies collected over a brief period of time.

Altitudinal patterns for latitudinally varying traits and polymorphic markers in Drosophila melanogaster from eastern Australia.

Altitudinal changes in traits and genetic markers can complement the studies on latitudinal patterns and provide evidence of natural selection because of climatic factors. In Drosophila melanogaster, latitudinal variation is well known but altitudinal patterns have rarely been investigated. Here, we examine five traits and five genetic markers on chromosome 3R in D. melanogaster collected at high and low altitudes from five latitudes along the eastern coast of Australia. Significant altitudinal differentiation was observed for cold tolerance, development time, ovariole number in unmated females, and the microsatellite marker DMU25686. Differences tended to match latitudinal patterns, in that trait values at high altitudes were also found at high latitudes, suggesting that factors linked to temperature are likely selective agents. Cold tolerance was closely associated with average temperature and other climatic factors, but no significant associations were detected for the other traits. Genes around DMU25686 represent good candidates for climatic adaptation.

In search of clinal variation in the period and clock timing genes in Australian Drosophila melanogaster populations.

Clinal variation for repeat number in the Thr-Gly region of the period circadian timing gene in Drosophila melanogaster was described in Europe and has subsequently been used as evidence of thermal selection on period alleles. To test for clinal variation in this gene along the east coast of Australia, the period polymorphism was scored on flies from multiple samples collected repeatedly over a 5-year interval, along with variation at another circadian rhythm locus, clock. For period, there was no consistent evidence of clinal variation in the 17 and/or 20 repeat alleles, although when average allele length was examined a weak consistent clinal pattern was detected. For clock there was no evidence of clinal variation in the two most common alleles or in average repeat size. These data are inconsistent with the reported patterns in Europe and suggest that clinal variation in timing genes needs to be re-examined in this region.

Microsatellites reveal a lack of structure in Australian populations of the diamondback moth, Plutella xylostella (L.).

The diamondback moth, Plutella xylostella, is renowned for developing resistance to insecticides and causing significant economic damage to Brassica vegetable crops throughout the world. Yet despite its economic importance, little is known about the population structure and movement patterns of this pest both at local and regional scales. In Australia, the movement patterns and insecticide resistance status of P. xylostella infesting canola, vegetables, forage brassicas and weeds have fundamental implications for the management of this pest. Here we use six polymorphic microsatellite loci to investigate population structure and gene flow in Australian populations of P. xylostella. Samples of P. xylostella from New Zealand, Malaysia, Indonesia and Kenya were also scored at these loci. We found no evidence of population structure within Australia, with most populations having low inbreeding coefficients and in Hardy-Weinberg equilibrium. In addition, a sample from the North Island of New Zealand was indistinguishable from the Australian samples. However, large genetic differences were found between the Australia/New Zealand samples and samples from Kenya, Malaysia and Indonesia. There was no relationship between genetic distance and geographic distance among Australian and New Zealand samples. Two of the loci were found to have null alleles, the frequency of which was increased in the populations outside the Australia/New Zealand region. We discuss these results with reference to insecticide resistance management strategies for P. xylostella in Australia.

A rapid shift in a classic clinal pattern in Drosophila reflecting climate change.

Geographical clines in genetic polymorphisms are widely used as evidence of climatic selection and are expected to shift with climate change. We show that the classic latitudinal cline in the alcohol dehydrogenase polymorphism of Drosophila melanogaster has shifted over 20 years in eastern coastal Australia. Southern high-latitude populations now have the genetic constitution of more northerly populations, equivalent to a shift of 4 degrees in latitude. A similar shift was detected for a genetically independent inversion polymorphism, whereas two other linked polymorphisms exhibiting weaker clinal patterns have remained relatively stable. These genetic changes are likely to reflect increasingly warmer and drier conditions and may serve as sensitive biomarkers for climate change.

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.

The Drosophila heat shock hsr-omega gene: an allele frequency cline detected by quantitative PCR.

The hsr-omega gene of Drosophila melanogaster produces RNA products both constitutively and at elevated levels in response to heat stress. A single-nucleotide difference in this gene that has been detected using denaturing gradient gel electrophoresis (DGGE) is responsible for an hsr-omegaa/b polymorphism, and selection experiments have indicated an association between the hsr-omegaa allele and susceptibility to heat stress. Since allele frequency estimates for population surveys using PCR and DGGE for single flies would be relatively time-consuming and expensive, we here develop a quantitative competitive-PCR method using mass-grind genomic DNA preparations for this purpose. Geographical and temporal variation of allele frequency at the hsr-omega locus in Australian populations of D. melanogaster are examined. Regular samples from a southern population through a summer season suggested stability of hsr-omegaa frequency. Field populations sampled from a approximately 2,250 km north-south transect along eastern Australia revealed a strong positive association between the frequency of hsr-omegaa and latitude, and marked spatial autocorrelation. Using appropriate analyses, strong association between population differences in hsr-omegaa frequencies and differences in temperature and rainfall measures, after controlling for latitudinal differences, support the idea that the cline in hsr-omegaa frequency may be attributable to some form of climatic selection.

Both allelic variation and expression of nuclear and cytoplasmic transcripts of Hsr-omega are closely associated with thermal phenotype in Drosophila.

Inducible heat shock genes are considered a major component of the molecular mechanisms that confer cellular protection against a variety of environmental stresses, in particular high temperature extremes. We have tested the association between expression of the heat shock RNA gene hsr-omega and thermoresistance by generating thermoresistant lines of Drosophila melanogaster after application of two distinct regimes of laboratory selection. One set of lines was selected for resistance to knockdown by heat stress and the other was similarly selected but before selection a mild heat exposure known to increase resistance (heat hardening) was applied. A cross between resistant and susceptible lines confirmed our earlier observation that increased thermal tolerance cosegregates with allelic variation in the hsr-omega gene. This cosegregating variation is attributed largely to two haplotype groups. Using quantitative reverse transcription-PCR, we find evidence for divergent phenotypic responses in the two selection regimes, involving both structural and regulatory changes in hsr-omega. Lines selected after hardening showed increased levels of the cytoplasmic transcript but decreased levels of the nuclear transcript. Lines selected without hardening showed decreased levels of the cytoplasmic transcript. The allelic frequency changes at hsr-omega could not by themselves account for the altered transcription patterns. Our results support the idea that the functional RNA molecules transcribed from hsr-omega are an important and polymorphic regulatory component of an insect thermoresistance phenotype.

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 degrees heat stress. Two selective regimes were used, one with and one without prior hardening. Mean knockdown times were increased from approximately 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.

A model of the Drosophila melanogaster glycerol-3-phosphate oxidase-encoding gene.

We have sequenced the Drosophila melanogaster gene encoding the mitochondrial (mt) enzyme glycerol-3-phosphate oxidase (GPO; EC 1.1.99.5) that is involved in flight and alcohol metabolism. Available data suggests a simple model for this gene that includes four exons. Exon I contains a mt import signal, exon II, a transmembrane segment and an FAD-binding site, and exon IV, an iron-sulfur centre.

Perturbation of gene frequencies in a natural population of Drosophila melanogaster: evidence for selection at the Adh locus.

The cellular population of Drosophila melanogaster at the Chateau Tahbilk Winery (Victoria, Australia) was perturbed for alcohol dehydrogenase (Adh) gene frequencies. Phenol oxidase (Phox) frequencies were also perturbed and monitored as a control. Subsequent gene frequency changes, together with information on population structure, indicated that selection acted on the chromosome regions of both loci. Adh gene frequencies returned to preperturbation levels in a predictable manner. A model in which the relative fitness of Adh phenotypes was determined by temperature-dependent specific activities of enzymes of Adh genotypes adequately accounts for the rate of gene frequency change at this locus. Thus temperature behaves as a selective agent in modulating Adh gene frequencies in this cellar environment.

The lepidopteran mitochondrial control region: structure and evolution.

For several species of lepidoptera, most of the approximately 350-bp mitochondrial control-region sequences were determined. Six of these species are in one genus, Jalmenus; are closely related; and are believed to have undergone recent rapid speciation. Recent speciation was supported by the observation of low interspecific sequence divergence. Thus, no useful phylogeny could be constructed for the genus. Despite a surprising conservation of control-region length, there was little conservation of primary sequences either among the three lepidopteran genera or between lepidoptera and Drosophila. Analysis of secondary structure indicated only one possible feature in common--inferred stem loops with higher-than-random folding energies--although the positions of the structures in different species were unrelated to regions of primary sequence similarity. We suggest that the conserved, short length of control regions is related to the observed lack of heteroplasmy in lepidopteran mitochondrial genomes. In addition, determination of flanking sequences for one Jalmenus species indicated (i) only weak support for the available model of insect 12S rRNA structure and (ii) that tRNA translocation is a frequent event in the evolution of insect mitochondrial genomes.

Long-chain dietary fatty acids affect the capacity of Drosophila melanogaster to tolerate ethanol.

Four-day post-hatch larvae (mid-third instar) of Drosophila melanogaster were fed an intermediate diet with or without supplement of an individual fatty acid for 2 d and then transferred to a diet with a growth-limiting level of 0.94 mol/L ethanol (5.5%, v/v) or an ethanol-free diet. The ethanol stress decreased survival and larval development rate but increased the weight of surviving adult males. Dietary long-chain fatty acids altered the fatty acid composition of tissue lipids of larvae. When an unsaturated fatty acid was fed, except for 18:2(n-6), the tissue level of total unsaturated fatty acids was markedly increased. Both saturated and unsaturated 18-carbon fatty acids shortened larval development time. Linoleic acid [18:2(n-6)] and linolenic acid [18:3(n-3)] enhanced survival overall and, together with stearic acid (18:0), gave marked protection from ethanol stress in terms of survival. Correlation analysis across the different fatty acid diets indicated a strong positive association between tissue 18-carbon fatty acid levels and ethanol tolerance and between 18-carbon fatty acid levels and development rate. No major differences were observed in the effects of the fatty acids on the Canton-S and OD4 (Tahbilk) wild-type strains. Thus, the fatty acid content of D. melanogaster larvae is important for growth and survival in ethanol-rich habitats.

Micro-evolution in a wine cellar population: an historical perspective.

The population of Drosophila melanogaster inside the wine cellar of Chateau Tahbilk of Victoria, Australia was found by McKenzie and Parsons (1974) to have undergone microevolution for greater alcohol tolerance when compared to the neighboring population outside the cellar. This triggered additional studies at Tahbilk, and at other wine cellars throughout the world. The contributions and interactions of researchers and the development of ideas on the ecology and genetics of this unique experimental system are traced. Although the ADH-F/ADH-S polymorphism was found to be maintained by selection in the Tahbilk populations, the selection is not significantly associated with alcohol tolerance. The environment inside the Tahbilk wine cellar is not as rich in ethanol as was originally anticipated, and selection that affects the alcohol dehydrogenase polymorphism may be more concerned with the relative efficiency with which ethanol is used as a nutrient by D. melanogaster. The synthesis and modification of lipids, particularly in membranes, appears to be important to alcohol tolerance. The studies of the Tahbilk population are at a crossroad. New experimental approaches promise to provide the keys to the selection that maintains the alcohol dehydrogenase polymorphism, and to factors that are important to alcohol tolerance and stress adaptation. From these research foundations at Tahbilk very significant contributions to our future understanding of the genetic processes of evolution can be made.