The evolution of antifungal peptides in Drosophila

An essential component of the immune system of animals is the production of antimicrobial peptides (AMPs). In vertebrates and termites the protein sequence of some AMPs evolves rapidly under positive selection, suggesting that they may be coevolving with pathogens. However, antibacterial peptides in Drosophila tend to be highly conserved. We have inferred the selection pressures acting on Drosophila antifungal peptides (drosomycins) from both the divergence of drosomycin genes within and between five species of Drosophila and polymorphism data from Drosophila simulans and D. melanogaster. In common with Drosophila antibacterial peptides, there is no evidence of adaptive protein evolution in any of the drosomycin genes, suggesting that they do not coevolve with pathogens. It is possible that this reflects a lack of specific fungal and bacterial parasites in Drosophila populations. The polymorphism data from both species differed from neutrality at one locus, but this was not associated with changes in the protein sequence. The synonymous site diversity was greater in D. simulans than in D. melanogaster, but the diversity both upstream of the genes and at nonsynonymous sites was similar. This can be explained if both upstream and nonsynonymous mutations are slightly deleterious and are removed more effectively from D. simulans due to its larger effective population size.     

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.     

Molecular population genetics of male accessory gland proteins in the Drosophila simulans complex

Accessory gland proteins are a major component of Drosophila seminal fluid. These proteins have a variety of functions and may be subject to sexual selection and/or antagonistic evolution between the sexes. Most population genetic data from these proteins are from D. melanogaster and D. simulans. Here, we extend the population genetic analysis of Acp genes to the other simulans complex species, D. mauritiana and D. sechellia. We sequenced population samples of seven Acp's from D. mauritiana, D. sechellia, and D. simulans. We investigated the population genetics of these genes on individual simulans complex lineages and compared Acp polymorphism and divergence to polymorphism and divergence from a set of non-Acp loci in the same species. Polymorphism and divergence data from the simulans complex revealed little evidence for adaptive protein evolution at individual loci. However, we observed a dramatically inflated index of dispersion for amino acid substitutions in the simulans complex at Acp genes, but not at non-Acp genes. This pattern of episodic bursts of protein evolution in Acp's provides the strongest evidence to date that the population genetic mechanisms driving Acp divergence are different from the mechanisms driving evolution at most Drosophila genes.     

Inferring Weak Selection from Patterns of Polymorphism and Divergence at ``silent' Sites in Drosophila DNA

Patterns of codon usage and ``silent'''' DNA divergence suggest that natural selection discriminates among synonymous codons in Drosophila. ``Preferred'''' codons are consistently found in higher frequencies within their synonymous families in Drosophila melanogaster genes. This suggests a simple model of silent DNA evolution where natural selection favors mutations from unpreferred to preferred codons (preferred changes). Changes in the opposite direction, from preferred to unpreferred synonymous codons (unpreferred changes), are selected against. Here, selection on synonymous DNA mutations is investigated by comparing the evolutionary dynamics of these two categories of silent DNA changes. Sequences from outgroups are used to determine the direction of synonymous DNA changes within and between D. melanogaster and Drosophila simulans for five genes. Population genetics theory shows that differences in the fitness effect of mutations can be inferred from the comparison of ratios of polymorphism to divergence. Unpreferred changes show a significantly higher ratio of polymorphism to divergence than preferred changes in the D. simulans lineage, confirming the action of selection at silent sites. An excess of unpreferred fixations in 28 genes suggests a relaxation of selection on synonymous mutations in D. melanogaster. Estimates of selection coefficients for synonymous mutations (3.6 <|N(e)s| < 1.3) in D. simulans are consistent with the reduced efficacy of natural selection (|N(e)s| < 1) in the three- to sixfold smaller effective population size of D. melanogaster. Synonymous DNA changes appear to be a prevalent class of weakly selected mutations in Drosophila.     

Adaptive protein evolution of X-linked and autosomal genes in Drosophila: implications for faster-X hypotheses

Patterns of sex chromosome and autosome evolution can be used to elucidate the underlying genetic basis of adaptative change. Evolutionary theory predicts that X-linked genes will adapt more rapidly than autosomes if adaptation is limited by the availability of beneficial mutations and if such mutations are recessive. In Drosophila, rates of molecular divergence between species appear to be equivalent between autosomes and the X chromosome. However, molecular divergence contrasts are difficult to interpret because they reflect a composite of adaptive and nonadaptive substitutions between species. Predictions based on faster-X theory also assume that selection is equally effective on the X and autosomes; this might not be true because the effective population sizes of X-linked and autosomal genes systematically differ. Here, population genetic and divergence data from Drosophila melanogaster, Drosophila simulans, and Drosophila yakuba are used to estimate the proportion of adaptive amino acid substitutions occurring in the D. melanogaster lineage. After gene composition and effective population size differences between chromosomes are controlled, X-linked and autosomal genes are shown to have equivalent rates of adaptive divergence with approximately 30% of amino acid substitutions driven by positive selection. The results suggest that adaptation is either unconstrained by a lack of beneficial genetic variation or that beneficial mutations are not recessive and are thus highly visible to natural selection whether on sex chromosomes or on autosomes.     

Nucleotide variation in the triosephosphate isomerase (Tpi) locus of Drosophila melanogaster and Drosophila simulans

DNA sequence variation in a 1.1-kb region including the coding portion of the Tpi locus was examined in 25 homozygous third-chromosome lines of Drosophila melanogaster, nine lines of Drosophila simulans, and one line of Drosophila yakuba. Our data show that the widespread allozyme polymorphism observed in cosmopolitan D. melanogaster is due to a glutamic acid substitution occurring in a phylogenetically conserved lysine that has been identified as part of the "hinged-lid" active site of the enzyme. This observation suggests that the replacement polymorphism may have important functional consequences. One replacement polymorphism was also observed in D. simulans, although its functional relevance is more difficult to assess, since it affects a site that is not strongly conserved. This amino acid change in D. simulans is associated with a single lineage possessing seven unique silent substitutions, which may be indicative of balancing selection or population subdivision. The absence of fixed amino acid differences between D. melanogaster and D. simulans and only a single difference with D. yakuba suggests that triose phosphate isomerase is under strong functional constraint. Silent variation is slightly higher for D. melanogaster than for D. simulans. Finally, we outline the general lack of evidence for old balanced polymorphisms at allozyme loci in D. melanogaster.      

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

In silico identification of new secretory peptide genes in Drosophila melanogaster

Bioactive peptides play critical roles in regulating most biological processes in animals. The elucidation of the amino acid sequence of these regulatory peptides is crucial for our understanding of animal physiology. Most of the (neuro)peptides currently known were identified by purification and subsequent amino acid sequencing. With the entire genome sequence of some animals now available, it has become possible to predict novel putative peptides. In this way, BLAST (Basic Local Alignment Searching Tool) analysis of the Drosophila melanogaster genome has allowed annotation of 36 secretory peptide genes so far. Peptide precursor genes are, however, poorly predicted by this algorithm, thus prompting an alternative approach described here. With the described searching program we scanned the Drosophila genome for predicted proteins with the structural hallmarks of neuropeptide precursors. As a result, 76 additional putative secretory peptide genes were predicted in addition to the 43 annotated ones. These putative (neuro)peptide genes contain conserved motifs reminiscent of known neuropeptides from other animal species. Peptides that display sequence similarities to the mammalian vasopressin, atrial natriuretic peptide, and prolactin precursors and the invertebrate peptides orcokinin, prothoracicotropic hormones, trypsin modulating oostatic factor, and Drosophila immune induced peptides (DIMs) among others were discovered. Our data hence provide further evidence that many neuropeptide genes were already present in the ancestor of Protostomia and Deuterostomia prior to their divergence. This bioinformatic study opens perspectives for the genome-wide analysis of peptide genes in other eukaryotic model organisms.     

The evolution of an alpha-esterase pseudogene inactivated in the Drosophila melanogaster lineage

Previous analyses of the alpha-esterase cluster of Drosophila melanogaster revealed 10 active genes and the DmalphaE4a-Psi pseudogene. Here, we reconstruct the evolution of the pseudogene from the sequences of 12 alleles from widely scattered D. melanogaster populations and single alleles from Drosophila simulans and Drosophila yakuba. All of the DmalphaE4a-Psi alleles contain numerous inactivating mutations, suggesting that pseudogene alleles are fixed in natural populations. Several lines of evidence also suggest that DmalphaE4a is now evolving without selective constraint in the D. melanogaster lineage. There are three polymorphic indels which result in frameshifts; a key nucleotide of the intron splice acceptor is polymorphic; the neutral mutation parameter is the same for replacement and silent sites; one of the nonsilent polymorphisms results in a stop codon; only 1 of the 13 replacement polymorphisms is biochemically conservative; residues that are conserved among active esterases have different states in DmalphaE4a-Psi; and there are about half as many transitional polymorphisms as transversional ones. In contrast, the D. simulans and D. yakuba orthologs DsalphaE4a and DyalphaE4a do not have the inactivating mutations of DmalphaE4a-Psi and appear to be evolving under the purifying selection typical of protein- encoding genes. For instance, there have been more substitutions in the introns than in the exons, and more in silent sites than in replacement sites. Furthermore, most of the amino acid substitutions that have occurred between DyalphaE4a and DsalphaE4a are located in sites that typically vary among active alpha-esterases rather than those that are usually conserved. We argue that the original alphaE4a gene had a function which it has lost since the divergence of the D. melanogaster and D. simulans lineages.     

Primary structure and in vitro antibacterial properties of the Drosophila melanogaster attacin C Pro-domain

In Drosophila melanogaster, seven distinct families of antimicrobial peptides with different structures and specificities are synthesized by the fat body and released into the hemolymph during the immune response. Using microscale high performance liquid chromatography, matrix-assisted laser desorption/ionization time-of-flight mass spectrometry, and Edman degradation, we have isolated and characterized from immune-challenged Drosophila two novel induced molecules, under the control of the Imd pathway, that correspond to post-translationally modified antimicrobial peptides or peptide fragments. The first molecule is a doubly glycosylated form of drosocin, an O-glycosylated peptide that kills Gram-negative organisms. The second molecule represents a truncated form of the pro-domain of the Drosophila attacin C carrying two post-translational modifications and has significant structural similarities to proline-rich antibacterial peptides including drosocin. We have synthesized this peptide and found that it is active against Gram-negative bacteria. Furthermore, this activity is potentiated when the peptide is used in combination with the Drosophila antimicrobial peptide cecropin A. The synergistic action observed between these two molecules suggests that the truncated post-translationally modified pro-domain of attacin C by itself may play an important role in the antimicrobial defense of Drosophila.     

Rapid evolution of the sex-determining gene,transformer: structural diversity and rate heterogeneity among sibling species of Drosophila

While developmentally regulated genes are generally conserved, transformer (tra), a key locus involved in the regulation of sexual differentiation, is highly diverged between species of Drosophila. With an aim to understand its divergence between sibling species, we investigated tra sequence variation among members of the Drosophila melanogaster species complex, D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In this species group, tra divergence is rapid yet clocklike and exhibits large differences in protein size. D. melanogaster contains a 13-amino acid tandem duplication, whereas D. sechellia possesses a 72-amino acid tandem duplication representing a 30% increase in total amino acid residues. We also found evidence of a nonrandom distribution of replacement substitutions and heterogeneity in substitution rates using clustering statistics and a codon substitution model. We show that tra's rapid divergence in this species complex is the result of generally lower selective constraints around regions that encode arginine-serine (RS) domains and a significantly higher rate of substitutions around the insertion site of D. sechellia's large duplication. The proximity of rapidly diverged regions to sites of nucleotide insertion suggests that higher local rates of mutation may provide a causal mechanism for TRA's rapid divergence in this subgroup. A comparison of tra orthologs across the genus Drosophila suggest that TRA maintains an assortment of RS domains for proper sex determining function while much of the protein evolves relatively unconstrained.     

The effect of developmental temperature on the genetic architecture underlying size and thermal clines in Drosophila melanogaster and D. simulans from the east coast of Australia

Body size and thermal tolerance clines in Drosophila melanogaster occur along the east coast of Australia. However the extent to which temperature affects the genetic architecture underlying the observed clinal divergence remains unknown. Clinal variation in these traits is associated with cosmopolitan chromosome inversions that cline in D. melanogaster. Whether this association influences the genetic architecture for these traits in D. melanogaster is unclear. Drosophila simulans shows linear clines in body size, but nonlinear clines in cold resistance. Clinally varying inversions are absent in D. simulans. Line-cross and clinal analyses were performed between tropical and temperate populations of D. melanogaster and D. simulans from the east coast of Australia to investigate whether clinal patterns and genetic effects contributing to clinal divergence in wing centroid size, thorax length, wing-to-thorax ratio, cold and heat resistance differed under different developmental temperatures (18 °C, 25 °C, and 29 °C). Developmental temperature influenced the genetic architecture in both species. Similarities between D. melanogaster and D. simulans suggest clinally varying inversion polymorphisms have little influence on the genetic architecture underlying clinal divergence in size in D. melanogaster. Differing genetic architectures across different temperatures highlight the need to consider different environments in future evolutionary and molecular studies of phenotypic divergence.     

Novel odorant-binding proteins expressed in the taste tissue of the fly

A taste tissue cDNA library of the fleshfly Boettcherisca peregrina was screened with a subtracted cDNA probe enriched with taste-receptor-tissue-specific cDNA. Seven genes were identified with sequence similarity to insect odorant-binding protein (OBP) genes. The predicted amino acid sequences of the genes contain the putative signal peptide sequence at the N-terminal and most of them conserve the six cysteines common to known insect OBPs. These genes show a high degree of sequence divergence with approximately 20% amino acid identity. The most striking feature was that all seven of these genes are expressed mainly in the taste tissues, such as the labellum and tarsus, unlike the known insect OBP genes expressed in olfactory tissue. The predicted amino acid sequences had the highest degree of sequence similarity to the Drosophila melanogaster OBPs named pheromone binding protein-related proteins (PBPRPs). These gene products are here referred to as gustatory PBP-related proteins (GPBPRPs) 1-7. Homologous GPBPRP genes were found also in D. melanogaster by database search and are shown to be expressed in Drosophila taste tissues.     

Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans

Surveys of molecular variation in Drosophila melanogaster and Drosophila simulans have suggested that diversity outside of Africa is a subset of that within Africa. It has been argued that reduced levels of diversity in non-African populations reflect a population bottleneck, adaptation to temperate climates, or both. Here, I summarize the available single-nucleotide polymorphism data for both species. A simple "out of Africa" bottleneck scenario is consistent with geographic patterns for loci on the X chromosome but not with loci on the autosomes. Interestingly, there is a trend toward lower nucleotide diversity on the X chromosome relative to autosomes in non-African populations of D. melanogaster, but the opposite trend is seen in African populations. In African populations, autosomal inversion polymorphisms in D. melanogaster may contribute to reduced autosome diversity relative to the X chromosome. To elucidate the role that selection might play in shaping patterns of variability, I present a summary of within- and between-species patterns of synonymous and replacement variation in both species. Overall, D. melanogaster autosomes harbor an excess of amino acid replacement polymorphisms relative to D. simulans. Interestingly, range expansion from Africa appears to have had little effect on synonymous-to-replacement polymorphism ratios.