Extensive amino acid polymorphism at the pgm locus is consistent with adaptive protein evolution in Drosophila melanogaster

Balanced Polymorphism at the Pgm-1 Locus of the Pompeii Worm Alvinella pompejana and Its Variant Adaptability Is Only Governed by Two QE Mutations at Linked Sites

The polychaete Alvinella pompejana lives exclusively on the walls of deep-sea hydrothermal chimneys along the East Pacific Rise (EPR), and displays specific adaptations to withstand the high temperatures and hypoxia associated with this highly variable habitat. Previous studies have revealed the existence of a balanced polymorphism on the enzyme phosphoglucomutase associated with thermal variations, where allozymes 90 and 100 exhibit different optimal activities and thermostabilities. Exploration of the mutational landscape of phosphoglucomutase 1 revealed the maintenance of four highly divergent allelic lineages encoding the three most frequent electromorphs over the geographic range of A. pompejana. This polymorphism is only governed by two linked amino acid replacements, located in exon 3 (E155Q and E190Q). A two-niche model of selection, including ‘cold’ and ‘hot’ conditions, represents the most likely scenario for the long-term persistence of these isoforms. Using directed mutagenesis and the expression of the three recombinant variants allowed us to test the additive effect of these two mutations on the biochemical properties of this enzyme. Our results are coherent with those previously obtained from native proteins, and reveal a thermodynamic trade-off between protein thermostability and catalysis, which is likely to have maintained these functional phenotypes prior to the geographic separation of populations across the Equator about 1.2 million years ago.

Molecular population genetics and selection in the glycolytic pathway

In this review, I discuss the evidence for differential natural selection acting across enzymes in the glycolytic pathway in Drosophila. Across the genome, genes evolve at very different rates and possess markedly varying levels of molecular polymorphism, codon bias and expression variation. Discovering the underlying causes of this variation has been a challenge in evolutionary biology. It has been proposed that both the intrinsic properties of enzymes and their pathway position have direct effects on their molecular evolution, and with the genomic era the study of adaptation has been taken to the level of pathways and networks of genes and their products. Of special interest have been the energy-producing pathways. Using both population genetic and experimental approaches, our laboratory has been engaged in a study of molecular variation across the glycolytic pathway in Drosophila melanogaster and its close relatives. We have observed a pervasive pattern in which genes at the top of the pathway, especially around the intersection at glucose 6-phosphate, show evidence for both contemporary selection, in the form of latitudinal allele clines, and inter-specific selection, in the form of elevated levels of amino acid substitutions between species. To further explore this question, future work will require corroboration in other species, expansion into tangential pathways, and experimental work to better characterize metabolic control through the pathway and to examine the pleiotropic effects of these genes on other traits and fitness components.

Flux control and excess capacity in the enzymes of glycolysis and their relationship to flight metabolism in Drosophila melanogaster

An important question in evolutionary and physiological genetics is how the control of flux-base phenotypes is distributed across the enzymes in a pathway. This control is often related to enzyme-specific levels of activity that are reported to be in excess of that required for demand. In glycolysis, metabolic control is frequently considered vested in classical regulatory enzymes, each strongly displaced from equilibrium. Yet the contribution of individual steps to control is unclear. To assess enzyme-specific control in the glycolytic pathway, we used P-element excision-derived mutagenesis in Drosophila melanogaster to generate full and partial knockouts of seven metabolic genes and to measure tethered flight performance. For most enzymes, we find that reduction to half of the normal activity has no measurable impact on wing beat frequency. The enzymes catalyzing near-equilibrium reactions, phosphoglucose isomerase, phosphoglucomutase, and triosephosphate isomerase fail to show any decline in flight performance even when activity levels are reduced to 17% or less. At reduced activities, the classic regulatory enzymes, hexokinase and glycogen phosphorylase, show significant drops in flight performance and are nearer to saturation. Our results show that flight performance is canalized or robust to the activity variation found in natural populations. Furthermore, enzymes catalyzing near-equilibrium reactions show strong genetic dominance down to low levels of activity. This implies considerable excess enzyme capacity for these enzymes.

Towards genetic markers in animal populations as biomonitors for human-induced environmental change

Genetic markers provide potentially sensitive indicators of changes in environmental conditions because the genetic constitution of populations is normally altered well before populations become extinct. Genetic indicators in populations include overall genetic diversity, genetic changes in traits measured at the phenotypic level, and evolution at specific loci under selection. While overall genetic diversity has rarely been successfully related to environmental conditions, genetically based changes in traits have now been linked to the presence of toxins and both local and global temperature shifts. Candidate loci for monitoring stressors are emerging from information on how specific genes influence traits, and from screens of random loci across environmental gradients. Drosophila research suggests that chromosomal regions under recent intense selection can be identified from patterns of molecular variation and a high frequency of transposable element insertions. Allele frequency changes at candidate loci have been linked to pesticides, pollutants and climate change. Nevertheless, there are challenges in interpreting allele frequencies in populations, particularly when a large number of loci control a trait and when interactions between alleles influence trait expression. To meet these challenges, population samples should be collected for longitudinal studies, and experimental programmes should be undertaken to link variation at candidate genes to ecological processes.

Molecular identification and expression of the phosphoglucomutase (PGM) gene from the Pacific oyster Crassostrea gigas

Phosphoglucomutase is a key enzyme in glycolysis and has been widely studied in vertebrates and some invertebrates but no molecular information is available in marine invertebrates despite the importance of this marker in ecological and genetical studies. In this work, we isolated a cDNA and the corresponding genomic sequence that encode PGM-2 locus in the Pacific oyster Crassostrea gigas. We used sequences drawn from the database to construct an evolutionary framework for examining the position of mollusc PGM sequences among prokaryotic and eukaryotic homologues and showed that oyster PGM gene organization was closer to vertebrates PGM genes than other invertebrates as previously found in other Lophotrochozoa species. We also investigated PGM mRNA expression in oyster tissues in response to xenobiotics (i.e hydrocarbons and pesticides). The results obtained showed that PGM mRNA expression is mostly up-regulated in the first steps of the response to pollutant exposure and is xenobiotic-dependent.

Natural and synthetic alleles provide complementary insights into the nature of selection acting on the Men polymorphism of Drosophila melanogaster

Two malic enzyme alleles, Men(113A) and Men(113G), occur at approximately equal frequency in North American populations of Drosophila melanogaster, while only Men(113A) occurs in African populations. We investigated the population genetics, biochemical characteristics, and selective potential of these alleles. Comparable levels of nucleotide polymorphism in both alleles suggest that the Men(113G) allele is not recently derived, but we find no evidence in the DNA sequence data for selection maintaining the polymorphism. Interestingly, the alleles differ in both V(max) and K(m) for the substrate malate. Triglyceride concentration and isocitrate dehydrogenase (IDH) and glucose-6-phosphate dehydrogenase (G6PD) activities are negatively correlated with the in vivo activities of the Men alleles. We examined the causality of the observed correlations using P-element excision-derived knockout alleles of the Men gene and found significant changes in the maximum activities of both IDH and G6PD, but not in triglyceride concentration, suggesting compensatory interactions between MEN, IDH, and G6PD. Additionally, we found significantly higher than expected levels of MEN activity in knockout heterozygotes, which we attribute to transvection effects. The distinct differences in biochemistry and physiology between the naturally occurring alleles and between the engineered alleles suggest the potential for selection on the Men locus.

Multiple signatures of positive selection downstream of notch on the X chromosome in Drosophila melanogaster

To identify genomic regions affected by the rapid fixation of beneficial mutations (selective sweeps), we performed a scan of microsatellite variability across the Notch locus region of Drosophila melanogaster. Nine microsatellites spanning 60 kb of the X chromosome were surveyed for variation in one African and three non-African populations of this species. The microsatellites identified an approximately 14-kb window for which we observed relatively low levels of variability and/or a skew in the frequency spectrum toward rare alleles, patterns predicted at regions linked to a selective sweep. DNA sequence polymorphism data were subsequently collected within this 14-kb region for three of the D. melanogaster populations. The sequence data strongly support the initial microsatellite findings; in the non-African populations there is evidence of a recent selective sweep downstream of the Notch locus near or within the open reading frames CG18508 and Fcp3C. In addition, we observe a significant McDonald-Kreitman test result suggesting too many amino acid fixations species wide, presumably due to positive selection, at the unannotated open reading frame CG18508. Thus, we observe within this small genomic region evidence for both recent (skew toward rare alleles in non-African populations) and recurring (amino acid evolution at CG18508) episodes of positive selection.

The role of natural selection in genetic differentiation of worldwide populations of Drosophila ananassae

The main evolutionary forces leading to genetic differentiation between populations are generally considered to be natural selection, random genetic drift, and limited migration. However, little empirical evidence exists to help explain the extent, mechanism, and relative role of these forces. In this study, we make use of the differential migration behavior of genes located in regions of low and high recombination to infer the role and demographic distribution of natural selection in Drosophila ananassae. Sequence data were obtained from 13 populations, representing almost the entire range of cosmopolitan D. ananassae. The pattern of variation at a 5.1-kb fragment of the furrowed gene, located in a region of very low recombination, appears strikingly different from that of 10 noncoding DNA fragments (introns) in regions of normal to high recombination. Most interestingly, two main haplotypes are present at furrowed, one being fixed in northern populations and the other being fixed or in high frequency in more southern populations. A cline in the frequency of one of these haplotypes occurs in parallel latitudinal transects. Taken together, significant clinal variation and a test against alternative models of natural selection provide evidence of two independent selective sweeps restricted to specific regions of the species range.

Multiple mutations and gene duplications conferring organophosphorus insecticide resistance have been selected at the Rop-1 locus of the sheep blowfly, Lucilia cuprina

Sequences of the esterase gene alpha E7 were compared across 41 isogenic (IV) strains of the sheep blowfly, Lucilia cuprina, and one strain of the sibling species, L. sericata. The 1.2-kb region sequenced includes sites of two insecticide resistance mutations. Gly137Asp confers resistance to organophosphorus insecticides (OPs), particularly preferring diethyl OPs such as diazinon, while Trp251Leu prefers dimethyl OPs, and particularly malathion, with the additional presence of carboxylester moieties. We found that there are just eight haplotypes among the 41 chromosomes studied: two Gly137Asp containing haplotypes, two Trp251Leu containing haplotypes, and four susceptible haplotypes, including the L. sericata sequence. While phylogenetic analysis of these haplotypes suggests that the Asp137 and Leu251 mutations each arose at least twice, evidence for recombination was detected across the region, therefore single origins for these resistance mutations cannot be ruled out. Levels of linkage disequilibrium in the data are high and significant hitchhiking is indicated by Fay and Wu' s H test but not the Tajima test. A test of haplotype diversity indicates a paucity of diversity compared with neutral expectations. Both these results are consistent with a very recent selective sweep at the Lc alphaE7 locus. Interestingly, gene duplications of three different combinations of OP resistant haplotypes were identified in seven of the isogenic (IV) strains. All three types of duplication involve an Asp137 and a Trp251 haplotype. To examine whether more haplotypes existed before the hypothesised selective sweep, fragments of alpha E7 surrounding the resistance mutations were amplified from pinned material dating back to before OPs were used. Four new sequence haplotypes, not sampled in the survey of extant haplotypes, were obtained that are all associated with susceptibility. This is suggestive of a higher historical level of susceptible allelic diversity at this locus.

The structure and population genetics of the breakpoints associated with the cosmopolitan chromosomal inversion In(3R)Payne in Drosophila melanogaster

We report here the breakpoint structure and sequences of the Drosophila melanogaster cosmopolitan chromosomal inversion In(3R)P. Combining in situ hybridization to polytene chromosomes and long-range PCR, we have identified and sequenced the distal and proximal breakpoints. The breakpoints are not simple cut-and-paste structures; gene fragments and small duplications of DNA are associated with both breaks. The distal breakpoint breaks the tolkin (tok) gene and the proximal breakpoint breaks CG31279 and the tolloid (tld) gene. Functional copies of all three genes are found at the opposite breakpoints. We sequenced a representative sample of standard (St) and In(3R)P karyotypes for a 2-kb portion of the tok gene, as well as the same 2 kb from the pseudogene tok fragment found at the distal breakpoint of In(3R)P chromosomes. The tok gene in St arrangements possesses levels of polymorphism typical of D. melanogaster genes. The functional tok gene associated with In(3R)P shows little polymorphism. Numerous single-base changes, as well as deletions and duplications, are associated with the truncated copy of tok. The overall pattern of polymorphism is consistent with a recent origin of In(3R)P, on the order of Ne generations. The identification of these breakpoint sequences permits a simple PCR-based screen for In(3R)P.

Thermal selection of PGM allozymes in newly founded populations of the thermotolerant vent polychaete Alvinella pompejana

Alvinella pompejana lives on the top of chimneys at deep-sea hydrothermal vents of the East Pacific Rise. It is thought to be one of the most thermotolerant and eurythermal metazoans. Our experimental approach combines methods of population genetics and biochemistry, considering temperature as a potential selective factor. Phosphoglucomutase (Pgm-1 locus) is one of the most polymorphic loci of A. pompejana and exhibits four alleles, from which alleles 90 and 100 dominate with frequencies of approximately 0.5 in populations. Results from previous studies suggested that allele 90 might be more thermostable than allele 100. Significant genetic differentiation was found by comparing contrasted microhabitats, especially the young, still hot, versus older and colder chimneys, with allele 90 being at highest frequency on young chimneys. Moreover the frequency of allele 90 was positively correlated with mean temperature at the opening of Alvinella tubes. In parallel, thermostability and thermal optimum experiments demonstrated that allele 90 is more thermostable and more active at higher temperatures than allele 100. This dataset supports an additive model of diversifying selection in which allele 90 is favoured on young hot chimneys but counterbalanced over the whole metapopulation by the dynamics of the vent ecosystem.

Single-locus latitudinal clines and their relationship to temperate adaptation in metabolic genes and derived alleles in Drosophila melanogaster

We report a study in Drosophila melanogaster of latitudinal clines for 23 SNPs embedded in 13 genes (Pgi, Gapdh1, UGPase, Pglym78, Pglym87, Eno, Men, Gdh, Sod, Pgk, Mdh1, TreS, Treh) representing various metabolic enzymes. Our samples are from 10 populations spanning latitude from southern Florida to northern Vermont. Three new clines with latitude were detected. These are the amino acid polymorphisms in the NAD-dependent glutamate dehydrogenase (Gdh) and trehalase (Treh) genes, and a silent site polymorphism in the UDP-glucose pyrophosphorylase gene (UGPase). The result, when combined with the overall incidence and pattern of reports for six other genes (Adh, Gpdh, Pgm, G6pd, 6Pgd, Hex-C), presents a picture of latitudinal clines in metabolic genes prevalent around the branch point of competing pathways. For six of the seven amino acid polymorphisms showing significant latitudinal clines in North America, the derived allele is the one increasing with latitude, suggesting temperate adaptation. This is consistent with a model of an Afrotropical ancestral species adapting to temperate climates through selection favoring new mutations.

Patterns of Polymorphism and Divergence from Noncoding Sequences of Drosophila melanogaster and D. simulans: Evidence for Nonequilibrium Processes

Despite the fact that D. melanogaster and D. simulans have been the central model system for molecular population genetics, few data are available for noncoding regions. Here, we present an analysis of population genetic data from intergenic regions and comparisons of these data to previously collected data from introns and exons. Polymorphisms and fixations were categorized as A/T to G/C or G/C to A/T changes and were polarized by inferring the ancestral state using both parsimony and maximum likelihood. Noncoding fixations in both D. melanogaster and D. simulans were consistent with equilibrium base-composition evolution. However, polarized noncoding polymorphisms, revealed a different pattern. Although A/T to G/C and G/C to A/T polymorphisms in D. simulans were consistent with equilibrium, we observed a highly significant dearth of A/T to G/C polymorphisms in D. melanogaster introns but not in intergenic sequences. Such data could be explained by recent evolution of mutational biases associated with transcription or by lineage-specific selection on base composition. These data reveal the complexity of evolutionary processes acting even on noncoding DNA in Drosophila.

Natural selection drives altitudinal divergence at the albumin locus in deer mice, Peromyscus maniculatus

In populations that are distributed across steep environmental gradients, the potential for local adaptation is largely determined by the spatial scale of fitness variation relative to dispersal distance. Since altitudinal gradients are generally characterized by dramatic ecological transitions over relatively short linear distances, adaptive divergence across such gradients will typically require especially strong selection to counterbalance the homogenizing effect of gene flow. Here we report the results of a study that was designed to test for evidence of adaptive divergence across an altitudinal gradient in a natural population of deer mice, Peromyscus maniculatus. We conducted a multilocus survey of allozyme variation across a steep altitudinal gradient in the southern Rocky Mountains that spanned several distinct biomes, from prairie grassland to alpine tundra. As a control for the effects of altitude, we also surveyed the same loci in mice sampled along a latitudinal transect through the prairie grassland that ran perpendicular to the east-west altitudinal transect. We used a coalescent-based simulation model to identify loci that deviated from neutral expectations, and we then assessed whether locus-specific patterns of variation were nonrandom with respect to altitude. Results indicated that the albumin locus (Alb) reflects a history of diversifying selection across the altitudinal gradient. This conclusion is supported by two main lines of evidence: (1) Alb was characterized by levels of divergence across the altitudinal transect that exceeded neutral expectations in two consecutive years of sampling (in contrast to the spatial pattern of variation across the latitudinal transect), and (2) levels of divergence at the Alb locus exhibited a positive association with altitudinal distance in both years (in contrast to the pattern observed at unlinked loci). We conclude that clinal variation at the Alb locus reflects a balance between gene flow and diversifying selection that results from elevational changes in fitness rankings among alternative genotypes.

The evolution of molecular genetic pathways and networks

There is growing interest in the evolutionary dynamics of molecular genetic pathways and networks, and the extent to which the molecular evolution of a gene depends on its position within a pathway or network, as well as over-all network topology. Investigations on the relationships between network organization, topological architecture and evolutionary dynamics provide intriguing hints as to how networks evolve. Recent studies also suggest that genetic pathway and network structures may influence the action of evolutionary forces, and may play a role in maintaining phenotypic robustness in organisms.

Influence of Two Wolbachia Strains on Population Structure of East AfricanDrosophila simulans

Drosophila simulans is hypothesized to have originated in continental East Africa or Madagascar. In this study, we investigated evolutionary forces operating on mitochondrial DNA (mtDNA) in populations of D. simulans from Zimbabwe, Malawi, Tanzania, and Kenya. Variation in mtDNA may be affected by positive selection, background selection, demographic history, and/or any maternally inherited factor such as the bacterial symbiont Wolbachia. In East Africa, the wRi and wMa Wolbachia strains associate with the siII or siIII mitochondrial haplogroups, respectively. To ask how polymorphism relates to Wolbachia infection status, we sequenced 1776 bp of mitochondrial DNA and 1029 bp of the X-linked per locus from 79 lines. The two southern populations were infected with wRi and exhibited significantly reduced mtDNA variation, while Wolbachia-uninfected siII flies from Tanzania and Kenya showed high levels of mtDNA polymorphism. These are the first known populations of D. simulans that do not exhibit reduced mtDNA variation. We observed no mitochondrial variation in the siIII haplogroup regardless of Wolbachia infection status, suggesting positive or background selection. These populations offer a unique opportunity to monitor evolutionary dynamics in ancestral populations that harbor multiple strains of Wolbachia.

Latitudinal clines inDrosophila melanogaster: Body size, allozyme frequencies, inversion frequencies, and the insulin-signalling pathway

Many latitudinal clines exist in Drosophila melanogaster: in adult body size, in allele frequency at allozyme loci, and in frequencies of common cosmopolitan inversions. The question is raised whether these latitudinal clines are causally related. This review aims to connect data from two very different fields of study, evolutionary biology and cell biology, in explaining such natural genetic variation in D. melanogaster body size and development time. It is argued that adult body size clines, inversion frequency clines, and clines in allele frequency at loci involved in glycolysis and glycogen storage are part of the same adaptive strategy. Selection pressure is expected to differ at opposite ends of the clines. At high latitudes, selection on D. melanogaster would favour high larval growth rate at low temperatures, and resource storage in adults to survive winter. At low latitudes selection would favour lower larval critical size to survive crowding, and increased male activity leading to high male reproductive success. Studies of the insulin-signalling pathway in D. melanogaster point to the involvement of this pathway in metabolism and adult body size. The genes involved in the insulin-signalling pathway are associated with common cosmopolitan inversions that show latitudinal clines. Each chromosome region connected with a large common cosmopolitan inversion possesses a gene of the insulin transmembrane complex, a gene of the intermediate pathway and a gene of the TOR branch. The hypothesis is presented that temperate D. melanogaster populations have a higher frequency of a 'thrifty' genotype corresponding to high insulin level or high signal level, while tropical populations possess a more 'spendthrift' genotype corresponding to low insulin or low signal level.

Linkage Disequilibrium Patterns Across a Recombination Gradient in African Drosophila melanogaster

Previous multilocus surveys of nucleotide polymorphism have documented a genome-wide excess of intralocus linkage disequilibrium (LD) in Drosophila melanogaster and D. simulans relative to expectations based on estimated mutation and recombination rates and observed levels of diversity. These studies examined patterns of variation from predominantly non-African populations that are thought to have recently expanded their ranges from central Africa. Here, we analyze polymorphism data from a Zimbabwean population of D. melanogaster, which is likely to be closer to the standard population model assumptions of a large population with constant size. Unlike previous studies, we find that levels of LD are roughly compatible with expectations based on estimated rates of crossing over. Further, a detailed examination of genes in different recombination environments suggests that markers near the telomere of the X chromosome show considerably less linkage disequilibrium than predicted by rates of crossing over, suggesting appreciable levels of exchange due to gene conversion. Assuming that these populations are near mutation-drift equilibrium, our results are most consistent with a model that posits heterogeneity in levels of exchange due to gene conversion across the X chromosome, with gene conversion being a minor determinant of LD levels in regions of high crossing over. Alternatively, if levels of exchange due to gene conversion are not negligible in regions of high crossing over, our results suggest a marked departure from mutation-drift equilibrium (i.e., toward an excess of LD) in this Zimbabwean population. Our results also have implications for the dynamics of weakly selected mutations in regions of reduced crossing over.

Mapping Determinants of Variation in Energy Metabolism, Respiration and Flight in Drosophila

We employed quantitative trait locus (QTL) mapping to dissect the genetic architecture of a hierarchy of functionally related physiological traits, including metabolic enzyme activity, metabolite storage, metabolic rate, and free-flight performance in recombinant inbred lines of Drosophila melanogaster. We identified QTL underlying variation in glycogen synthase, hexokinase, phosphoglucomutase, and trehalase activity. In each case variation mapped away from the enzyme-encoding loci, indicating that trans-acting regions of the genome are important sources of variation within the metabolic network. Individual QTL associated with variation in metabolic rate and flight performance explained between 9 and 35% of the phenotypic variance. Bayesian QTL analysis identified epistatic effects underlying variation in flight velocity, metabolic rate, glycogen content, and several metabolic enzyme activities. A region on the third chromosome was associated with expression of the glucose-6-phosphate branchpoint enzymes and with metabolic rate and flight performance. These genomic regions are of special interest as they may coordinately regulate components of energy metabolism with effects on whole-organism physiological performance. The complex biochemical network is encoded by an equally complex network of interacting genetic elements with potentially pleiotropic effects. This has important consequences for the evolution of performance traits that depend upon these metabolic networks.

DNA sequence variation and latitudinal associations in hsp23, hsp26 and hsp27 from natural populations of Drosophila melanogaster

Heat shock genes are considered to be likely candidate genes for environmental stress resistance. Nucleotide variation in the coding sequence of the small heat shock genes (hsps) hsp26 and hsp27 from Drosophila melanogaster was studied in flies originating from the Netherlands and eastern Australia. The hsp26 gene was polymorphic for an insertion/deletion of three extra amino acids and two nonsynonymous changes in all populations. The hsp27 gene exhibited two nonsynonymous changes and three synonymous mutations. The hsp26 polymorphism showed a latitudinal cline along the east coast of Australia. This pattern was not confounded by the fact that the shsps are located in the inversion In(3 L)P which also shows a latitudinal cline in eastern Australia. A similar latitudinal cline was found for the previously described variation in hsp23, while frequencies of hsp27 alleles did not change with latitude. These findings suggest that variation at two of the shsps or closely linked loci are under selection in natural populations of D. melanogaster.

Drosophila melanogaster and D. simulans rescue strains produce fit offspring, despite divergent centromere-specific histone alleles

The interaction between rapidly evolving centromere sequences and conserved kinetochore machinery appears to be mediated by centromere-binding proteins. A recent theory proposes that the independent evolution of centromere-binding proteins in isolated populations may be a universal cause of speciation among eukaryotes. In Drosophila the centromere-specific histone, Cid (centromere identifier), shows extensive sequence divergence between D. melanogaster and the D. simulans clade, indicating that centromere machinery incompatibilities may indeed be involved in reproductive isolation and speciation. However, it is presently unclear whether the adaptive evolution of Cid was a cause of the divergence between these species, or merely a product of postspeciation adaptation in the separate lineages. Furthermore, the extent to which divergent centromere identifier proteins provide a barrier to reproduction remains unknown. Interestingly, a small number of rescue lines from both D. melanogaster and D. simulans can restore hybrid fitness. Through comparisons of cid sequence between nonrescue and rescue strains, we show that cid is not involved in restoring hybrid viability or female fertility. Further, we demonstrate that divergent cid alleles are not sufficient to cause inviability or female sterility in hybrid crosses. Our data do not dispute the rapid divergence of cid or the coevolution of centromeric components in Drosophila; however, they do suggest that cid underwent adaptive evolution after D. melanogaster and D. simulans diverged and, consequently, is not a speciation gene.

Enhanced phosphorylation and enzymatic activity of phosphoglucomutase by the Btk29A tyrosine kinase in Drosophila

The Drosophila Btk29A tyrosine kinase is suggested to be involved in diverse processes, although its target proteins are unknown. In the present study, we investigated substrates of Btk29A tyrosine kinase by expressing a catalytically activated form of Btk29A-P1 (Btk-EG) in Drosophila compound eyes. Expression in eye disks led to the development of the rough-eye phenotype and increased tyrosine phosphorylation of a 65-kDa protein. Partial amino acid sequence analysis of this protein showed that it was phosphoglucomutase. Phosphoglucomutase activity in heads from Btk-EG-expressing flies was higher than that in controls, suggesting that the levels of tyrosine phosphorylation and activity of the enzyme are associated with Btk29A tyrosine kinase activity.

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.

Applications of selective neutrality tests to molecular ecology

This paper reviews how statistical tests of neutrality have been used to address questions in molecular ecology are reviewed. The work consists of four major parts: a brief review of the current status of the neutral theory; a review of several particularly interesting examples of how statistical tests of neutrality have led to insight into ecological problems; a brief discussion of the pitfalls of assuming a strictly neutral model if it is false; and a discussion of some of the opportunities and problems that molecular ecologists face when using neutrality tests to study natural selection.

DNA sequence variation at a duplicated gene: excess of replacement polymorphism and extensive haplotype structure in the Drosophila melanogaster bicoid region

The bicoid (bcd) gene of Drosophila has played an important role in understanding the system of developmental regulatory genes that controls segmentation in the fruit fly. Several studies in Drosophila and closely related insects suggest that bcd may be the result of a gene duplication in the Dipteran lineage. In addition, the presence of a large, conserved secondary structure in the 3' untranslated region (UTR) makes the bcd gene a good candidate for studying compensatory evolution and the relationship between RNA secondary structure and patterns of standing variation in natural populations. Despite these interesting aspects, a population-level analysis has until now not been performed on bcd. In this study, DNA sequence variation was examined for a 4-kb region of the bcd gene, including a portion of the 5' UTR, the entire coding region, and the 3' UTR, for 25 Drosophila melanogaster isofemale lines from Zimbabwe and one allele from D. simulans. Statistical tests revealed a significant excess of replacement polymorphisms in the D. melanogaster lineage that are clustered in two putative linker regions of the Bicoid protein. This result is consistent with a relaxation of selective constraints in these regions. In addition, we found a distinct haplotype structure and a significantly smaller number of haplotypes than predicted by the standard neutral model. It is unlikely that the haplotype structure is maintained by epistatic selection acting on the secondary structure in the 3' UTR or by the association of the bcd gene with polymorphic inversions. Instead, our two main observations, namely the occurrence of a haplotype structure and the excess of replacement polymorphisms, may indicate that the selective history of this gene is rather complex, involving both the relaxation of purifying selection in some parts of the protein and the action of positive selection in other parts of the gene region.

Distinguishing between selection and population expansion in an experimental lineage of bacteriophage T7

Experimental evolution of short-lived organisms offers the opportunity to study the dynamics of polymorphism over time in a controlled environment. Here, we characterize DNA polymorphism data over time for four genes in bacteriophage T7. Our experiment ran for 2500 generations and populations were sampled after 500, 2000, and 2500 generations. We detect positive selection, purifying ("negative") selection, and population expansion in our experiment. We also present a statistical test that is able to distinguish demographic from selective events, processes that are hard to identify individually because both often produce an excess of rare mutations. Our "heterogeneity test" modifies common statistics measuring the frequency spectrum of polymorphism (e.g., Fu and Li's D) by looking for processes producing different patterns on nonsynonymous and synonymous mutations. Test results agree with the known conditions of the experiment, and we are therefore confident that this test offers a tool to evaluate natural populations. Our results suggest that instances of segregating deleterious mutations may be common, but as yet undetected, in nature.

Doxycycline-induced expression of sense and inverted-repeat constructs modulates phosphogluconate mutase (Pgm) gene expression in adult Drosophila melanogaster

BACKGROUND

A tetracycline-regulated (conditional) system for RNA interference (RNAi) would have many practical applications. Such a strategy was developed using RNAi of the gene for phosphogluconate mutase (Pgm). Pgm is a candidate lifespan regulator: PgmS allele frequency is increased by selection for increased lifespan, whereas PgmM and PgmF allele frequencies are decreased.

RESULTS

The Pgm alleles were cloned and sequenced and were found to differ by amino-acid substitutions consistent with the relative electrophoretic mobilities of the proteins. The 'tet-on' doxycycline-regulated promoter system was used to overexpress PgmS in a wild-type (PgmM) background. Enzyme activity increases of two- to five-fold were observed in five independent transgenic lines. Tet-on was also used to drive expression of an inverted-repeat fragment of Pgm coding region. The inverted-repeat transcript was expected to form a dsRNA hairpin, induce RNAi, and thereby reduce endogenous Pgm gene expression at the RNA level. Endogenous Pgm RNA levels in adult flies were found to be reduced or eliminated by doxycycline treatment in five independent inverted-repeat transgenic lines. Our results show that doxycycline-regulated expression of inverted-repeat constructs can cause a conditional reduction in specific gene expression. The effect of sense and inverted-repeat construct expression on lifespan was assayed in multiple transgenic lines. Under the conditions tested, altered Pgm gene expression had no detectable effect on adult Drosophila lifespan.

CONCLUSIONS

A system for conditional RNAi in Drosophila adults shows promise for assay of gene functions during aging. Our results indicate that Pgm does not have a simple strong effect on longevity.

Genome-wide variation in the human and fruitfly: a comparison

Average levels of nucleotide diversity are ten-fold lower in humans than in the fruitfly, Drosophila melanogaster. Despite this difference, apparently as a result of a lower population size, patterns of genomic diversity are strikingly similar in being correlated with local rates of recombination, and influenced by similar interactions between positive natural selection and recombination. Both species also show lower levels of variation on average in non-African compared to African populations, reflecting a similar evolutionary history and perhaps both natural selection and founder effects in new environments.

The functional impact of Pgm amino acid polymorphism on glycogen content in Drosophila melanogaster

Earlier studies of the common PGM allozymes in Drosophila melanogaster reported no in vitro activity differences. However, our study of nucleotide variation observed that PGM allozymes are a heterogeneous mixture of amino acid polymorphisms. In this study, we analyze 10 PGM protein haplotypes with respect to PGM activity, thermostability, and adult glycogen content. We find a twofold difference in activity among PGM protein haplotypes that is associated with a threefold difference in glycogen content. The latitudinal clines for several Pgm amino acid polymorphisms show that high PGM activity, and apparently higher flux to glycogen synthesis, parallel the low activity clines at G6PD for reduced pentose shunt flux in northern latitudes. This suggests that amino acid polymorphism is under selection at this branch point and may be favored for increased metabolic storage associated with stress resistance and adaptation to temperate regions.

Clinal variation for amino acid polymorphisms at the Pgm locus in Drosophila melanogaster

Clinal variation is common for enzymes in the glycolytic pathway for Drosophila melanogaster and is generally accepted as an adaptive response to different climates. Although the enzyme phosphoglucomutase (PGM) possesses several allozyme polymorphisms, it is unique in that it had been reported to show no clinal variation. Our recent DNA sequence investigation of Pgm found extensive cryptic amino acid polymorphism segregating with the allozyme alleles. In this study, we characterize the geographic variation of Pgm amino acid polymorphisms at the nucleotide level along a latitudinal cline in the eastern United States. A survey of 15 SNPs across the Pgm gene finds significant clinal differentiation for the allozyme polymorphisms as well as for many of the cryptic amino acid polymorphisms. A test of independence shows that pervasive linkage disequilibrium across this gene region can explain many of the amino acid clines. A single Pgm haplotype defined by two amino acid polymorphisms shows the strongest correlation with latitude and the steepest change in allele frequency across the cline. We propose that clinal selection at Pgm may in part explain the extensive amino acid polymorphism at this locus and is consistent with a multilocus response to selection in the glycolytic pathway.

Contrasting molecular population genetics of four hexokinases in Drosophila melanogaster, D. simulans and D. yakuba

As part of a larger study contrasting patterns of variation in regulatory and nonregulatory enzymes of the central metabolic pathways we have examined the molecular variation in four uncharacterized hexokinase genes unique to muscle, fat body, and testis in Drosophila melanogaster, D. simulans, and D. yakuba. Earlier isoenzyme studies had designated these genes as Hex-A, Hex-C, and Hex-t. There are two tightly linked testes-specific genes designated here as Hex-t1 and Hex-t2. Substantial and concordant differences across species are seen in levels of both amino acid and silent polymorphism. The flight muscle form Hex-A is the most conserved followed by the fat body hexokinase Hex-C and testis-specific hexokinases Hex-t1 and Hex-t2. While constraints acting at the amino acid level are expected, the silent polymorphisms follow this pattern as well. All genes are in regions of normal recombination, therefore hitchhiking and background selection are not likely causes of interlocus differences. In D. melanogaster latitudinal clines are seen for amino acid polymorphisms at the Hex-C and Hex-t2 loci. There is evidence for accelerated amino acid substitution in Hex-t1 that has lost residues known to be associated with glucose and glucose-6-phosphate binding. D. simulans shows substantial linkage phase structuring that suggests historical population subdivision.