Molecular evolution of the period gene in Drosophila athabasca

We measured nucleotide variability within and between the three semispecies of the Drosophila athabasca complex, at the period (per) gene by using a polymerase chain reaction-based four-cutter restriction- enzyme analysis. The levels of polymorphism varied considerably between the three semispecies. Our results for per, combined with previous data for X-linked allozymes, suggest that the X chromosome in the western- northern semispecies is less variable than expected under an equilibrium-neutral model. Both the pattern of divergence between the semispecies and a cladistic clustering of per haplotypes support the previously hypothesized grouping of eastern A and eastern B as the two most recently diverged semispecies. A 21-bp in-frame segment in the region of per which shares sequence similarity with the neuronal development gene single minded is deleted in all eastern A and eastern B flies examined but is present in all of the western-northern flies and all other published per sequences. Despite these hints that there may be significant differences at the per gene between the semispecies, especially the western-northern group versus the two eastern groups, there is no compelling evidence that per is involved in the mating song differences between the semispecies.      

Molecular evolution of a repetitive region within the per gene of Drosophila

The clock gene period (per) controls a number of biological rhythms in Drosophila. In D. melanogaster, per has a repetitive region that encodes a number of alternating threonine-glycine residues. We sequenced and compared this region from several different Drosophila species belonging to various groups within the Drosophila and Sophophora subgenera. This part of per shows a great variability in both DNA sequence and length. Furthermore, analysis of the data suggests that changes in the length of this variable region might be associated with amino acid replacements in the more conserved flanking sequences.      

Molecular evolution of a Y chromosome to autosome gene duplication in Drosophila

In contrast to the rest of the genome, the Y chromosome is restricted to males and lacks recombination. As a result, Y chromosomes are unable to respond efficiently to selection, and newly formed Y chromosomes degenerate until few genes remain. The rapid loss of genes from newly formed Y chromosomes has been well studied, but gene loss from highly degenerate Y chromosomes has only recently received attention. Here, we identify and characterize a Y to autosome duplication of the male fertility gene kl-5 that occurred during the evolution of the testacea group species of Drosophila. The duplication was likely DNA based, as other Y-linked genes remain on the Y chromosome, the locations of introns are conserved, and expression analyses suggest that regulatory elements remain linked. Genetic mapping reveals that the autosomal copy of kl-5 resides on the dot chromosome, a tiny autosome with strongly suppressed recombination. Molecular evolutionary analyses show that autosomal copies of kl-5 have reduced polymorphism and little recombination. Importantly, the rate of protein evolution of kl-5 has increased significantly in lineages where it is on the dot versus Y linked. Further analyses suggest this pattern is a consequence of relaxed purifying selection, rather than adaptive evolution. Thus, although the initial fixation of the kl-5 duplication may have been advantageous, slightly deleterious mutations have accumulated in the dot-linked copies of kl-5 faster than in the Y-linked copies. Because the dot chromosome contains seven times more genes than the Y and is exposed to selection in both males and females, these results suggest that the dot suffers the deleterious effects of genetic linkage to more selective targets compared with the Y chromosome. Thus, a highly degenerate Y chromosome may not be the worst environment in the genome, as is generally thought, but may in fact be protected from the accumulation of deleterious mutations relative to other nonrecombining regions that contain more genes.     

Molecular evolution of sex-biased genes in Drosophila

Studies of morphology, interspecific hybridization, protein/DNA sequences, and levels of gene expression have suggested that sex-related characters (particularly those involved in male reproduction) evolve rapidly relative to non-sex-related characters. Here we report a general comparison of evolutionary rates of sex-biased genes using data from cDNA microarray experiments and comparative genomic studies of Drosophila. Comparisons of nonsynonymous/synonymous substitution rates (d(N)/d(S)) between species of the D. melanogaster subgroup revealed that genes with male-biased expression had significantly faster rates of evolution than genes with female-biased or unbiased expression. The difference was caused primarily by a higher d(N) in the male-biased genes. The same pattern was observed for comparisons among more distantly related species. In comparisons between D. melanogaster and D. pseudoobscura, genes with highly biased male expression were significantly more divergent than genes with highly biased female expression. In many cases, orthologs of D. melanogaster male-biased genes could not be identified in D. pseudoobscura through a Blast search. In contrast to the male-biased genes, there was no clear evidence for accelerated rates of evolution in female-biased genes, and most comparisons indicated a reduced rate of evolution in female-biased genes relative to unbiased genes. Male-biased genes did not show an increased ratio of nonsynonymous/synonymous polymorphism within D. melanogaster, and comparisons of polymorphism/divergence ratios suggest that the rapid evolution of male-biased genes is caused by positive selection.     

Molecular evolution of phosphoprotein phosphatases in Drosophila

Phosphoprotein phosphatases (PPP), these ancient and important regulatory enzymes are present in all eukaryotic organisms. Based on the genome sequences of 12 Drosophila species we traced the evolution of the PPP catalytic subunits and noted a substantial expansion of the gene family. We concluded that the 18-22 PPP genes of Drosophilidae were generated from a core set of 8 indispensable phosphatases that are present in most of the insects. Retropositons followed by tandem gene duplications extended the phosphatase repertoire, and sporadic gene losses contributed to the species specific variations in the PPP complement. During the course of these studies we identified 5, up till now uncharacterized phosphatase retrogenes: PpY+, PpD5+, PpD6+, Pp4+, and Pp6+ which are found only in some ancient Drosophila. We demonstrated that all of these new PPP genes exhibit a distinct male specific expression. In addition to the changes in gene numbers, the intron-exon structure and the chromosomal localization of several PPP genes was also altered during evolution. The G-C content of the coding regions decreased when a gene moved into the heterochromatic region of chromosome Y. Thus the PPP enzymes exemplify the various types of dynamic rearrangements that accompany the molecular evolution of a gene family in Drosophilidae.     

Molecular evolution of the histone 3 multigene family in the Drosophila melanogaster species subgroup

Molecular evolution of the histone multigene family was studied by cloning and sequencing regions of the histone 3 gene in the Drosophila melanogaster species subgroup. Analysis of the nucleotide substitution pattern showed that in the coding region synonymous changes occurred more frequently to A or T in contrast to the GC-rich base composition, while in the 3' region the nucleotide substitutions were most likely in equilibrium. These results suggested that the base composition at the third codon position of the H3 gene, i.e., codon usage, has been changing to A or T in the Drosophila melanogaster species subgroup.     

Molecular evolution of glutathione S-transferases in the genus Drosophila

As classical phase II detoxification enzymes, glutathione S-transferases (GSTs) have been implicated in insecticide resistance and may have evolved in response to toxins in the niche-defining feeding substrates of Drosophila species. We have annotated the GST genes of the 12 Drosophila species with recently sequenced genomes and analyzed their molecular evolution. Gene copy number variation is attributable mainly to unequal crossing-over events in the large delta and epsilon clusters. Within these gene clusters there are also GST genes with slowly diverging orthologs. This implies that they have their own unique functions or have spatial/temporal expression patterns that impose significant selective constraints. Searches for positively selected sites within the GSTs identified G171K in GSTD1, a protein that has previously been shown to be capable of metabolizing the insecticide DDT. We find that the same radical substitution (G171K) in the substrate-binding domain has occurred at least three times in the Drosophila radiation. Homology-modeling places site 171 distant from the active site but adjacent to an alternative DDT-binding site. We propose that the parallel evolution observed at this site is an adaptive response to an environmental toxin and that sequencing of historical alleles suggests that this toxin was not a synthetic insecticide.     

Molecular evolution of the sex peptide network in Drosophila

Successful reproduction depends on interactions between numerous proteins beyond those involved directly in gamete fusion. Although such reproductive proteins evolve in response to sexual selection pressures, how networks of interacting proteins arise and evolve as reproductive phenotypes change remains an open question. Here, we investigated the molecular evolution of the ‘sex peptide network’ of Drosophila melanogaster, a functionally well‐characterized reproductive protein network. In this species, the peptide hormone sex peptide (SP) and its interacting proteins cause major changes in female physiology and behaviour after mating. In contrast, females of more distantly related Drosophila species do not respond to SP. In spite of these phenotypic differences, we detected orthologs of all network proteins across 22 diverse Drosophila species and found evidence that most orthologs likely function in reproduction throughout the genus. Within SP‐responsive species, we detected the recurrent, adaptive evolution of several network proteins, consistent with sexual selection acting to continually refine network function. We also found some evidence for adaptive evolution of several proteins along two specific phylogenetic lineages that correspond with increased expression of the SP receptor in female reproductive tracts or increased sperm length, respectively. Finally, we used gene expression profiling to examine the likely degree of functional conservation of the paralogs of an SP network protein that arose via gene duplication. Our results suggest a dynamic history for the SP network in which network members arose before the onset of robust SP‐mediated responses and then were shaped by both purifying and positive selection.     

Molecular evolution of sex-biased genes in the Drosophila ananassae subgroup

Background  

Genes with sex-biased expression often show rapid molecular evolution between species. Previous population genetic and comparative genomic studies of Drosophila melanogaster and D. simulans revealed that male-biased genes have especially high rates of adaptive evolution. To test if this is also the case for other lineages within the melanogaster group, we investigated gene expression in D. ananassae, a species that occurs in structured populations in tropical and subtropical regions. We used custom-made microarrays and published microarray data to characterize the sex-biased expression of 129 D. ananassae genes whose D. melanogaster orthologs had been classified previously as male-biased, female-biased, or unbiased in their expression and had been studied extensively at the population-genetic level. For 43 of these genes we surveyed DNA sequence polymorphism in a natural population of D. ananassae and determined divergence to the sister species D. atripex and D. phaeopleura.     

Molecular evolution of Drosophila odorant receptor genes

A total of 752 odorant receptor (Or) genes, including pseudogenes, were identified in 11 Drosophila species and named after their orthologs in Drosophila melanogaster. The 813 Or genes, including 61 from D. melanogaster, were classified into 59 orthologous groups that are well supported by gene phylogeny. By reconciling with the gene family phylogeny, we estimated the number of gene duplication/loss events and intron gain/loss events in the species phylogeny. We found that these events are particularly frequent in Drosophila grimshawi, Drosophila willistoni, and obscura group. More than half of the duplicated genes stay as tandem arrays, whose size range from 2 to 8. These genes vary in sequence and some likely underwent positive selection, indicating that the gene duplication was important for flies to acquire new olfactory functions. We hypothesize that Or genes conferred the basic olfactory repertoire to ancestral flies before the speciation of the Drosophila and Sophophora subgenera about 40 Mya. This repertoire has been largely maintained in the current species, whereas lineage-specific gene duplication seems to have led to additional specialization in some species in response to specific ecological conditions.     

Concerted evolution within the Drosophila dumpy gene

We have determined by reverse Southern analysis and direct sequence comparisons that most of the dumpy gene has evolved in the dipteran and other insect orders by purifying selection acting on amino acid replacements. One region, however, is evolving rapidly due to unequal crossing over and/or gene conversion. This region, called "PIGSFEAST," or PF, encodes in D. melanogaster 30-47 repeats of 102 amino acids rich in serines, threonines, and prolines. We show that the processes of concerted evolution have been operating on all species of Drosophila examined to date, but that an adjacent region has expanded in Anopheles gambiae, Aedes aegypti, and Tribolium castaneum, while the PF repeats are reduced in size and number. In addition, processes of concerted evolution have radically altered the codon usage patterns in D. melanogaster, D. pseudoobscura, and D. virilis compared with the rest of the dumpy gene. We show also that the dumpy gene is expressed on the inner surface of the micropyle of the mature oocyte and propose that, as in the abalone system, concerted evolution may be involved in adaptive changes affecting Dumpy's possible role in sperm-egg recognition.     

Neutral evolution of the sex-determining gene transformer in Drosophila

The amino acid sequence of the transformer (tra) gene exhibits an extremely rapid rate of evolution among Drosophila species, although the gene performs a critical step in sex determination. These changes in amino acid sequence are the result of either natural selection or neutral evolution. To differentiate between selective and neutral causes of this evolutionary change, analyses of both intraspecific and interspecific patterns of molecular evolution of tra gene sequences are presented. Sequences of 31 tra alleles were obtained from Drosophila americana. Many replacement and silent nucleotide variants are present among the alleles; however, the distribution of this sequence variation is consistent with neutral evolution. Sequence evolution was also examined among six species representative of the genus Drosophila. For most lineages and most regions of the gene, both silent and replacement substitutions have accumulated in a constant, clock-like manner. In exon 3 of D. virilis and D. americana we find evidence for an elevated rate of nonsynonymous substitution, but no statistical support for a greater rate of nonsynonymous relative to synonymous substitutions. Both levels of analysis of the tra sequence suggest that, although the gene is evolving at a rapid pace, these changes are neutral in function.     

Molecular evolution of the calmodulin gene

On the basis of the intron/exon organization and the intramolecular homology of DNA sequences, I propose a novel model for genesis of the calmodulin gene. A primordial calmodulin gene consisting of 51 base pairs (17 amino acids) was subjected to three-fold duplication to create modern calmodulin with four calcium-binding subdomains. The model elucidates the seemingly enigmatic positions of splice junctions observed in calmodulin genes.     

Rh gene evolution in primates: study of intron sequences

By amplification and sequencing of RH gene intron 4 of various primates we demonstrate that an Alu-Sx-like element has been inserted in the RH gene of the common ancestor of humans, apes, Old World monkeys, and New World monkeys. The study of mouse and lemur intron 4 sequences allowed us to precisely define the insertion point of the Alu-Sx element in intron 4 of the RH gene ancestor common to Anthropoidea. Like humans, chimpanzees and gorillas possess two types of RH intron 4, characterized by the presence (human RHCE and ape RHCE-like genes) or absence (human RHD and ape RHD-like genes) of the Alu-Sx element. This led us to conclude that in the RH common ancestor of humans, chimpanzees, and gorillas, a duplication of the common ancestor gene gave rise to two genes, one differing from the other by a 654-bp deletion encompassing an Alu-Sx element. Moreover, most of chimpanzees and some gorillas posses two types of RHD-like intron 4. The introns 4 of type 1 have a length similar to that of human RHD intron 4, whereas introns 4 of type 2 display an insertion of 12 bp. The latest insertion was not found in the human genome (72 individuals tested). The study of RH intron 3 length polymorphism confirmed that, like humans, chimpanzees and gorillas possess two types of intron 3, with the RHD-type intron 3 being 289 bases shorter than the RHCE intron 3. By amplification and sequencing of regions encompassing introns 3 and 4, we demonstrated that chimpanzee and gorilla RH-like genes displayed associations of introns 3 and 4 distinct to those found in man. Altogether, the results demonstrate that, as in humans, chimpanzee and gorilla RH genes experienced intergenic exchanges.     

Molecular evolution of the Sex-Ratio inversion complex in Drosophila pseudoobscura: analysis of the Esterase-5 gene region

The Sex-Ratio chromosome in Drosophila pseudoobscura is subject to meiotic drive. It is associated with a series of three nonoverlapping paracentric inversions on the right arm of the X chromosome. The esterase-5 gene region has been localized to section 23 within the subbasal inversion of the Sex-Ratio inversion complex, making esterase- 5 a convenient locus for molecular evolutionary analyses of the Sex- Ratio inversion complex and the associated drive system. A 504-bp fragment of noncoding, intergenic DNA from the esterase-5 gene region was amplified and sequenced from 14 Sex-Ratio and 14 Standard X chromosomes of D. pseudoobscura, and from 9 X chromosomes of its two sibling species, Drosophila persimilis and Drosophila miranda. There is extensive sequence differentiation between the Sex-Ratio and Standard chromosomal types. The common Standard chromosome is highly polymorphic, while, as expected from either the neutral mutation theory or the selective sweep hypothesis, the rarer Sex-Ratio chromosome has much less within-chromosome nucleotide polymorphism. We estimate that the Standard and Sex-Ratio chromosomes in D. pseudoobscura diverged between 700,000 and 1.3 Mya, or at least 2 million generations ago. The clustering of D. pseudoobscura Sex-Ratio chromosomes in a neighbor- joining phylogeny indicates a fairly old, monophyletic origin in this species. It appears from these data that Sex-Ratio genes were present prior to the divergence of D. pseudoobscura and D. persimilis and that both the Standard and Sex-Ratio chromosomes of D. persimilis were derived from the Standard chromosome of D. pseudoobscura after the inversion events that isolated the D. pseudoobscura Sex-Ratio chromosome.      

Phylogenetic utility of the mitochondrial cytochrome oxidase gene: Molecular evolution of the Drosophila buzzatii species complex

Phylogenetic relationships among eight species of the Drosophila buzzatii species complex (D. mulleri subgroup; D. repleta species group) and D. hamatofila were determined by sequencing the mitochondrial cytochrome oxidase subunit 1, II, and III genes. The species examined included members of the martensis cluster (D. martensis, D. starmeri, D. venezolana), the buzzatii cluster (D. buzzatii, D. serido, D. borborema), and the stalkeri cluster (D. stalkeri, D. richardsoni). The molecular phylogeny was found to be congruent with the chromosomal inversion phylogeny. Analyzing the cytochrome oxidase subunits separately revealed that not all the subunits seem to have the same phylogenetic information content. Parameters are discussed that might explain these differences. Correspondence to: G. Spicer     

Molecular evolution of the period gene in sandflies

The molecular evolution of the clock gene period was studied in Phlebotomine sandflies (Diptera: Psychodidae). Comparison of the synonymous and nonsynonymous substitution rates between sandflies and Drosophila revealed a significantly higher evolutionary rate in the latter in three of the four regions analyzed. The differences in rate were higher in the sequences flanking the Thr–Gly repetitive domain, a region that has expanded in Drosophila but remained stable and short in sandflies, a result consistent with the coevolutionary scenario proposed for this region of the gene. An initial phylogenetic analysis including eight neotropical sandfly species and one from the Old World was also carried out. The results showed that only the subgenus Nyssomyia is well supported by distance (neighbor-joining) and maximum parsimony analysis. The grouping of the other species from the subgenus Lutzomyia and Migonei group shows very low bootstrap values and is not entirely consistent with classical morphological systematics of the genus Lutzomyia.     

Molecular evolution of a duplication: the sex-peptide (Acp70A) gene region of Drosophila subobscura and Drosophila madeirensis

In Drosophila melanogaster, the Acp70A gene, which is involved in the postmating reactions of the female, is a single-copy gene. However, in Drosophila subobscura, the gene is duplicated and both copies are transcribed. To study the molecular evolution of the duplication, a 2.1- kb fragment encompassing both copies of the duplication was sequenced for 10 lines of D. subobscura and one line of Drosophila madeirensis. Estimates of the divergence between the two copies of the duplicated region and between the two species studied, D. subobscura and D. madeirensis, revealed that both copies of the Acp70a gene had evolved independently since their duplication. The ratio of nonsynonymous to silent divergence between copies was generally higher than one. The McDonald and Kreitman test revealed an excess of nonsynonymous changes fixed since the duplication and before the split of the D. subobscura and D. madeirensis lineages. These results point to natural selection driving protein evolution after the duplication. Specifically, adaptive evolution appears to have caused the initial differentiation between copies of the N-terminal parts of the proteins, while purifying selection could be responsible for the high conservation of the C- terminal parts.