Nucleotide divergence of the rp49 gene region between Drosophila melanogaster and two species of the Obscura group of Drosophila

Summary A 2.1-kb SStI fragment including the rp49 gene and the 3 end of the -serendipity gene has been cloned and sequenced in Drosophila pseudoobscura. rp49 maps at region 62 on the tip of chromosome II of this species. Both the coding and flanking regions have been aligned and compared with those of D. subobscura. There is no evidence for heterogeneity in the rate of silent substitution between the rp49 coding region and the rate of substitutions in flanking regions, the overall silent divergence per site being 0.19. Noncoding regions also differ between both species by different insertions/deletions, some of which are related to repeated sequences. The rp49 region of D. pseudoobscura shows a strong codon bias similar to those of D. subobscura and D. melanogaster. Comparison of the rates of silent (K _S ) and nonsilent (K _a ) substitutions of the rp49 gene and other genes completely sequenced in D. pseudoobscura and D. melanogaster confirms previous results indicating that rp49 is evolving slowly both at silent and nonsilent sites. According to the data for the rp49 region, D. pseudoobscura and D. subobscura lineages would have diverged some 9 Myr ago, if one assumes a divergence time of 30 Myr for the melanogaster and obscura groups.Offprint requests to: C. Segarra     

A Genomic Comparison of Faster-Sex, Faster-X, and Faster-Male Evolution Between Drosophila melanogaster and Drosophila pseudoobscura

A genomic comparison of Drosophila melanogaster and Drosophila pseudoobscura provides a unique opportunity to investigate factors involved in sequence divergence. The chromosomal arrangements of these species include an autosomal segment in D. melanogaster which is homologous to part of the X chromosome in D. pseudoobscura. Using orthologues to calculate rates of nonsynonymous (d_N) substitutions, we found genes on the X chromosome to be significantly more diverged than those on the autosomes, but it is not true for segment 3L-XR which is autosomal in D. melanogaster (3L) and X-linked in D. pseudoobscura (XR). We also found that the median d_N values for genes having reproductive functions in either the male, the female, or both sexes are higher than those for sequences without reproductive function and even higher for sequences involved in male-specific function. These estimates of divergence for male sex-related sequences are most likely underestimates, as the very rapidly evolving reproductive genes would tend to lose homology sooner and thus not be included in the comparison of orthologues. We also noticed a high proportion of male reproductive genes among the othologous genes with the highest rates of d_N. Reproductive genes with and without an orthologue in D. pseudoobscura were compared among D. melanogaster, D. simulans, and D. yakuba and it was found that there were in fact higher rates of divergence in the group without a D. pseudoobscura orthologue. These results, from widely separated taxa, bolster the thesis that sexual system genes experience accelerated rates of change in comparison to nonsexual genes in evolution and speciation. [Reviewing Editor: Dr. Willie J. Swanson]     

An evolutionary model for the duplication and divergence of esterase genes in Drosophila

Summary The esterase 5 (Est-5 = gene, EST 5 = protein) enzyme in Drosophila pseudoobscura is encoded by one of three paralogous genes, Est-5A, Est5B, and Est-5C, that are tightly clustered on the right arm of the X chromosome. The homologous Est-6 locus in Drosophila melanogaster has only one paralogous neighbor, Est-P. Comparisons of coding and flanking DNA sequences among the three D. pseudoobscura and two D. melanogaster genes suggest that two paralogous genes were present before the divergence of D. pseudoobscura from D. melanogaster and that, later, a second duplication occurred in D. pseudoobscura. Nucleotide sequences of the coding regions of the three D. pseudoobscura genes showed 78–85% similarity in pairwise comparisons, whereas the relatedness between Est-6 and Est-P was only 67%. The higher degree of conservation in D. pseudoobscura likely results from the comparatively recent divergence of Est-5B and Est-5C and from possible gene conversion events between Est-5A and Est-5B. Analyses of silent and replacement site differences in the two exons of the paralogous and orthologous genes in each species indicate that common selective forces are acting on all five loci. Further evidence for common purifying selective constraints comes from the conservation of hydropathy profiles and proposed catalytic residues. However, different levels of amino acid substitution between the paralogous genes in D. melanogaster relative to those in D. pseudoobscura suggest that interspecific differences in selection also exist.Offprint requests to: R.C. Richmond     

Historical effective size and the level of genetic diversity in Drosophila melanogaster and Drosophila pseudoobscura

We report the results of a sequential gel electrophoretic study of protein variation in Drosophila melanogaster and its comparison with D. pseudoobscura. The number of alleles and mean heterozygosity were lower in D. melanogaster than in D. pseudoobscura. On the other hand, geographical populations of Drosophila melanogaster have been shown to be much more differentiated than those of D. pseudoobscura. The results suggest that in D. melanogaster low-frequency alleles have been lost during the colonization process and that major alleles have become differentiated among populations. Population bottlenecks, due to various causes, appear to have played a significant role in the shaping of genetic variation in natural populations of many species. It is proposed that a comparison of genetic variation at homologous gene loci between related species can bring out effects of historical bottlenecks and provide an alternative approach for analyzing causes of genetic variation in natural populations.We thank the Natural Science and Engineering Research Council of Canada for financial support (Grant A0235 to R.S.S.).     

Correlated Evolution of Synonymous and Nonsynonymous Sites in <Emphasis Type="Italic">Drosophila</Emphasis>

Recent work has shown that Drosophila melanogaster genes with fast-evolving nonsynonymous sites have lower codon usage bias. This pattern has been attributed to interference between positive selection at nonsynonymous sites and weak selection on codon usage. Here we have looked for this correlation in a much larger and less biased dataset, comprising 630 gene pairs from D. melanogaster and D. yakuba. We confirmed that there is a negative correlation between the rate of nonsynonymous substitutions (d_N) and codon bias in D. melanogaster. We then tested the interference hypothesis and other alternative explanations, including one involving gene expression. We found that d_N indeed correlates with the level of gene expression. Given that gene expression is a strong determinant of codon bias, the relationship between d_N and codon bias might be a by-product of gene expression. However, our tests show that none of the hypotheses we consider seem to explain the data fully.This article contains online supplementary material.Reviewing Editor: Dr. John Huelsenbeck     

Diversification and Independent Evolution of Troponin C Genes in Insects

Troponin C (TpnC), the calcium-binding subunit of the troponin regulatory complex in the muscle thin filament, is encoded by multiple genes in insects. To understand how TpnC genes have evolved, we characterized the gene number and structure in a number of insect species. The TpnC gene complement is five genes in Drosophilidae as previously reported for D. melanogaster. Gene structures are almost identical in D. pseudoobscura, D. suboboscura, and D. virilis. Developmental patterns of expression are also conserved in Drosophila subobscura and D. virilis. Similar, but not completely equivalent, TpnC gene repertoires have been identified in the Anopheles gambiae and Apis mellifera genomes. Insect TpnC sequences can be divided into three groups, allowing a systematic classification of newly identified genes. The pattern of expression of the Apis mellifera genes essentially agrees with the pattern in Drosophilidae, providing further functional support to the classification. A model for the evolution of the TpnC genes is proposed including the most likely pathway of insect TpnC diversification. Our results suggest that the rapid increase in number and sequence specialization of the adult Type III isoforms can be correlated with the evolution of the holometabolous mode of development and the acquisition of asynchronous indirect flight muscle function in insects. This evolutionarily specialization has probably been achieved independently in different insect orders.Reviewing Editor: Dr. Rüdiger Cerff     

Sequence evolution of theGpdh gene in theDrosophila virilis species group

The nucleotide sequence of theGpdh gene from six taxa,D. virilis, D. lummei, D. novamexicana, D. a. americana, D. a. texana andD. ezoana, belonging to thevirilis species group was determined to examine details of evolutionary change in the structure of theGpdh gene. TheGpdh gene is comprised of one 5 non-translated region, eight exons, seven introns and three 3 non-translated regions. Exon/intron organization was identical in all the species examined, but different from that of mammals. Interspecific nucleotide divergence in the entireGpdh gene followed the common pattern: it was low in the exon, high in the intron and intermediate in the non-translated regions. The degree of nucleotide divergence differed within these regions, suggesting that selection exerts constraints differentially on nucleotide change of theGpdh gene. A phylogenetic tree of thevirilis phylad constructed from nucleotide variation of total sequence was consistent with those obtained from other data.Nucleotide sequences for theGpdh gene ofD. lummei, D. novamexicana, D. a. americana, D. a. texana andD. ezoana have been submitted to GenBank with accession numbers D50087, D50088, D50089, D50090 and D50091.     

Evolutionary differences between alternative and constitutive protein-coding regions of alternatively spliced genes of Drosophila

A total of 790 Drosophila melanogaster genes that are alternatively spliced in a coding region and have orthologs in Drosophila pseudoobscura were studied. It proved that nucleotide substitutions are accumulated in alternative coding regions more rapidly than in constitutive coding regions. Moreover, the evolutionary patterns of alternative regions differing in insertion-deletion mechanisms (use of alternative promoters, splicing sites, or polyadenylation sites) differ significantly. The synonymous substitution rate in coding regions of genes varies more strongly than the nonsynonymous substitution rate. The patterns of substitutions in different classes of alternative regions of Drosophila melanogaster and mammals differ considerably.     

Association Between Levels of Coding Sequence Divergence and Gene Misregulation in <Emphasis Type="Italic">Drosophila</Emphasis> Male Hybrids

Previous studies have shown widespread conservation of gene expression levels between species of the Drosophila melanogaster subgroup as well as a positive correlation between coding sequence divergence and expression level divergence between species. Meanwhile, large-scale misregulation of gene expression level has been described in interspecific sterile hybrids between D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. Using data from gene expression analysis involving D. simulans, D. melanogaster, and their hybrids, we observed a significant positive correlation between protein sequence divergence and gene expression differences between hybrids and their parental species. Furthermore, we demonstrate that underexpressed misregulated genes in hybrids are evolving more rapidly at the protein sequence level than nonmisregulated genes or overexpressed misregulated genes, highlighting the possible effects of sexual and natural selection as male-biased genes and nonessential genes are the principal gene categories affected by interspecific hybrid misregulation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Carlo G. Artieri and Wilfried Haerty contributed equally to this publication.     

Length polymorphism in the threonine-glycine-encoding repeat region of theperiod gene inDrosophila

Summary Single-fly polymerase chain reaction amplification and direct DNA sequencing revealed high levels of length polymorphism in the threonine-glycine encoding repeat region of theperiod (per) gene in natural populations ofDrosophila melanogaster. DNA comparison of two alleles of identical lengths gave a high number of synonymous substitutions suggesting an ancient time of separation. However detailed examination of the sequences of different Thr-Gly length variants indicated that this divergence could be understood in terms of four deletion/insertion events. InDrosophila pseudoobscura a length polymorphism is observed in a five-amino acid degenerate repeat, which corresponds tomelanogaster's Thr-Gly domain. In spite of the differences betweenD. melanogaster andD. pseudoobscura in the amino acid sequence of the repeats, the predicted secondary structures suggest evolutionary and mechanistic constraints on theper protein of these two species.     

Positive Selection on Nucleotide Substitutions and Indels in Accessory Gland Proteins of the Drosophila pseudoobscura Subgroup

Genes encoding reproductive proteins often diverge rapidly due to positive selection on nucleotide substitutions. While this general pattern is well established, the extent to which specific reproductive genes experience similar selection in different clades has been little explored, nor have possible targets of positive selection other than nucleotide substitutions, such as indels, received much attention. Here, we inspect for the signature of positive selection in the genes encoding five accessory gland proteins (Acps) (Acp26Aa, Acp32CD, Acp53Ea, Acp62F, and Acp70A) originally described from Drosophila melanogaster but with recognizable orthologues in the D. pseudoobscura subgroup. We compare patterns of selection within the D. psuedoobscura subgroup to those in the D. melanogaster subgroup. Similar patterns of positive selection were found in Acp26Aa and Acp62F in the two subgroups, while Acp53Ea and Acp70A experienced purifying selection in both subgroups. These proteins have thus remained targets for similar types of selection over long (>21-MY) periods of time. We also found several indel substitutions and polymorphisms in Acp26Aa and Acp32CD. These indels occur in the same regions as positively selected nucleotide substitutions for Acp26Aa in the D. pseudoobscura subgroup but not in the D. melanogaster subgroup. Rates of indel substitution within Acp26Aa in the D. pseudoobscura subgroup were up to several times those in noncoding regions of the Drosophila genome. This suggests that indel substitutions may be under positive selection and may play a key role in the divergence of some Acps. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Willis Swanson]     

Identification of One Intron Loss and Phylogenetic Evolution of Dfak Gene in the Drosophila melanogaster Species Group

Intron loss and its evolutionary significance have been noted in Drosophila. The current study provides another example of intron loss within a single-copy Dfak gene in Drosophila. By using polymerase chain reaction (PCR), we amplified about 1.3 kb fragment spanning intron 5–10, located in the position of Tyr kinase (TyK) domain of Dfak gene from Drosophila melanogaster species group, and observed size difference among the amplified DNA fragments from different species. Further sequencing analysis revealed that D. melanogaster and D. simulans deleted an about 60 bp of DNA fragment relative to other 7 Drosophila species, such as D. elegans, D. ficusphila, D. biarmipes, D. takahashii, D. jambulina, D. prostipennis and D. pseudoobscura, and the deleted fragment located precisely in the position of one intron. The data suggested that intron loss might have occurred in the Dfak gene evolutionary process of D. melanogaster and D. simulans of Drosophila melanogaster species group. In addition, the constructed phylogenetic tree based on the Dfak TyK domains clearly revealed the evolutionary relationships between subgroups of Drosophila melanogaster species group, and the intron loss identified from D. melanogaster and D. simulans provides a unique diagnostic tool for taxonomic classification of the melanogaster subgroup from other group of genus Drosophila.     

Structural and genetic studies of the proliferation disrupter genes of Drosophila simulans and D. melanogaster

To examine whether structural and functional differences exist in the proliferation disrupter (prod) genes between Drosophila simulans and D. melanogaster, we analyzed and compared both genes. The exon–intron structure of the genes was found to be the same – three exons were interrupted by two introns, although a previous report suggested that only one intron existed in D. melanogaster. The prod genes of D. simulans and D. melanogaster both turn out to encode 346 amino acids, not 301 as previously reported for D. melanogaster. The numbers of nucleotide substitutions in the prod genes was 0.0747 ±  per synonymous site and 0.0116 ± 0.0039 per replacement site, both comparable to those previously known for homologous genes between D. simulans and D. melanogaster. Genetic analysis demonstrated that D. simulans PROD can compensate for a deficiency of D. melanogaster PROD in hybrids. The PRODs of D. simulans and D. melanogaster presumably share the same function and a conserved working mechanism. The prod gene showed no significant interaction with the lethality of the male hybrid between these species. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dosage compensation of X-linked heat-shock puffs in Drophila pseudoobscura

Four major puffs are inducible by heat shock in the larval salivary gland chromosomes of D. pseudoobscura. Two of these puffs are present at 23 and 39–40 on the right arm of the X chromosome and two are present at 53 and 58 on chromosome 2. By means of in situ hybridization, residual homologies were demonstrated between the puffs at 23 in D. pseudoobscura and at 63C in D. melanogaster, and between the two chromosome 2 puffs of D. pseudoobscura and 87A and 87C of D. melanogaster. RNA synthesis was monitored as a function of ³H-uridine incorporation in the major heat-induced puffs of D. pseudoobscura and was found to be equivalent in males and females indicating dosage compensation of the two X-linked loci. The evolution of the regulatory controls of these genes is discussed.     

Analysis of a swallow homologue from Drosophila pseudoobscura

We analyzed a functional homologue of the swallow gene from Drosophila pseudoobscura. The swallow gene of D. melanogaster plays an essential role in localizing bicoid mRNA in oocytes, and swallow mutant embryos show anterior pattern defects that result from the lack of localization of the bicoid morphogen. The pseudoobscura homologue rescues the function of swallow mutants when introduced into the genome of D. melanogaster, and its expression is similar to that of the melanogaster gene. The predicted pseudoobscura and melanogaster proteins are 49% identical and 69% conserved. The coiled-coil domain previously identified in the melanogaster swallow protein is strongly conserved in the pseudoobscura homologue, but the weak similarity of the melanogaster swallow protein to the RNP class of RNA-binding proteins is not conserved in the pseudoobscura homologue. These and other observations suggest a structural role for swallow in localizing bicoid mRNA, perhaps as part of the egg cytoskeleton. Received: 3 August 1999 / Accepted: 29 September 1999     

Molecular structure of rare but geographically widespread sn-glycerol-3-phosphate dehydrogenase ‘ultra-fast’ electrophoretic alleles in Drosophila melanogaster

Six naturally occurring but rare alleles of sn-glycerol-3-phosphate dehydrogenase (Gpdh) in Drosophila melanogaster have been investigated in this study. They all belong to a class of Gpdh ^UF (ultra-fast) alleles, because their electrophoretic mobilities are faster than that of the Gpdh ^F (fast) allele. The Gpdh ^UF variants are widespread, and have been reported from five continents. DNA sequence analysis has shown that the change in electrophoretic mobility was in each allele caused by a single amino acid residue substitution in the encoded protein. In the Xiamen ^UF allele it is a substitution of lysine (AAA) to asparagine (AAT) in exon 1 (residue 3). An asparagine (AAT) to aspartate (GAT) change was found in exon 6 (residue 336) in the Iowa ^UF and Netherlands ^UF alleles. The mobility of the Raleigh ^UF allele was altered by a valine (GTG) to glutamate (GAG) substitution in exon 3 (residue 76). Two mutations were detected in the Brazzaville ^UF allele: a lysine (AAG) to methionine (ATG) substitution in exon 2 (residue 68) is responsible for the ultra-fast phenotype of this variant, while a tyrosine (TAT) to phenylalanine (TTT) substitution in exon 4 (residue 244) is not expected to alter the electrophoretic mobility of the encoded protein. These results indicate that the Gpdh ^UF alleles originate from different mutational events, and only two of them — Iowa ^UF and Netherlands ^UF — might share a common ancestry. The GPDH activity of the Iowa ^UF allele is intermediate between those of the Gpdh ^S and Gpdh ^F control stocks. The other Gpdh ^UF variants have lower activities than the controls: Xiamen ^UF -83%, Raleigh ^UF -80% and Brazzaville ^UF -73% of the Gpdh ^F control.     

Distribution and evolution of introns in drosophila amylase genes

While the two amylase genes of Drosophila melanogaster are intronless, the three genes of D. pseudoobscura harbor a short intron. This raises the question of the common structure of the Amy gene in Drosophila species. We have investigated the presence or absence of an intron in the amylase genes of 150 species of Drosophilids. Using polymerase chain reaction (PCR), we have amplified a region that surrounds the intron site reported in D. pseudoobscura and a few other species. The results revealed that most species contain an intron, with a variable size ranging from 50 to 750 bp, although the very majoritary size was around 60–80 bp. Several species belonging to different lineages were found to lack an intron. This loss of intervening sequence was likely due to evolutionarily independent and rather frequent events. Some other species had both types of genes: In the obscura group, and to a lesser extent in the ananassae subgroup, intronless copies had much diverged from intron-containing genes. Base composition of short introns was found to be variable and correlated with that of the surrounding exons, whereas long introns were all A-T rich. We have extended our study to non-Drosophilid insects. In species from other orders of Holometaboles, Lepidoptera and Hymenoptera, an intron was found at an identical position in the Amy gene, suggesting that the intron was ancestral. Received: 23 October 1995 / Accepted: 5 March 1996     

Molecular evolution of the Rh3 gene in Drosophila

Previous investigations into the evolution of the Drosophila opsin gene family are extended by inter- and intraspecific DNA sequence comparisons of the Rh3 locus in the melanogaster subgroup and D. pseudoobscura. Two separate statistical tests of the neutral-mutation hypothesis suggest that random genetic drift is responsible for virtually all of the observed amino acid replacement substitutions within the melanogaster subgroup. Analyses incorporating the D. pseudoobscura sequences are enigmatic due to the accumulation of multiple substitutions, because the McDonald-Kreitman test is not applicable to species comparisons that approach mutational saturation. However, the data from D. pseudoobscura are not inconsistent with selective neutrality. The ratio of amino acid polymorphisms within species to fixed differences between species imply that these are approximately 31 possible neutral single-step amino-acid-replacement substitutions at this locus. Synonymous substitutions are unevenly distributed among the structural domains of the Rh3 gene. Patterns of synonymous polymorphism are analyzed with respect to GC content and codon bias, and are compared to other loci from the same species.     

Molecular characterization of the zerknüllt region of the Antennapedia complex of D. subobscura

We have characterized at the molecular level the zerknüllt (zen) region of the Drosophila subobscura Antennapedia complex. The sequence comparison between D. subobscura and D. melanogaster shows an irregular distribution of the conserved and diverged regions, with the homeobox and a putative activating domain completely conserved. Comparisons of the promoter sequence and pattern of expression of the gene during development suggest that the regulation of zen has been conserved during evolution. The conservation of zen expression in a subpopulation of the polar cells indicates the existence of an important role in such cells. We describe a transitory segmented pattern of expression of zen in both species, suggesting the existence of interactions with a pair rule gene. Some indirect clues indicate that the z2 gene might be absent from the D. subobscura genome. A chromosome walk initiated to reach the proboscipedia gene of D. subobscura reveals that the distance between pb and zen is at least four times the one described for D. melanogaster and for D. pseudoobscura. Finally, we present cytological evidence showing that the ANT-C is inverted in D. subobscura as compared to D. melanogaster.