Molecular evolution in the Drosophila melanogaster species subgroup: frequent parameter fluctuations on the timescale of molecular divergence

Although mutation, genetic drift, and natural selection are well established as determinants of genome evolution, the importance (frequency and magnitude) of parameter fluctuations in molecular evolution is less understood. DNA sequence comparisons among closely related species allow specific substitutions to be assigned to lineages on a phylogenetic tree. In this study, we compare patterns of codon usage and protein evolution in 22 genes (>11,000 codons) among Drosophila melanogaster and five relatives within the D. melanogaster subgroup. We assign changes to eight lineages using a maximum-likelihood approach to infer ancestral states. Uncertainty in ancestral reconstructions is taken into account, at least to some extent, by weighting reconstructions by their posterior probabilities. Four of the eight lineages show potentially genomewide departures from equilibrium synonymous codon usage; three are decreasing and one is increasing in major codon usage. Several of these departures are consistent with lineage-specific changes in selection intensity (selection coefficients scaled to effective population size) at silent sites. Intron base composition and rates and patterns of protein evolution are also heterogeneous among these lineages. The magnitude of forces governing silent, intron, and protein evolution appears to have varied frequently, and in a lineage-specific manner, within the D. melanogaster subgroup.     

Evolution of nested genes with special reference to cuticle proteins inDrosophila melanogaster

Summary One of the pupal cuticle protein (PCP) genes has been found within an intron of aDrosophila housekeeping gene (theGart locus) that encodes three enzymes involved in the purine pathway. This intronic gene has been described as a gene within a gene, and the gene is now called a “nested” gene. Because the intronic PCP gene has sequence similarity with the larval cuticle protein (LCP) gene, it may have been derived from one of the LCP genes or their ancestral gene. We have studied possible phylogenetic relationships among these five genes by comparing nucleotide sequences of four LCP genes with that of the PCP gene. The results obtained suggest that the PCP gene may have originated from an ancestral gene before duplication of the LCP genes occurred. Using the number of synonymous (silent) substitutions, we then estimated the divergence time between the PCP gene and the LCP genes to be about 70 million years (Myr). The divergence time estimated is much larger than that for the sibling species ofD. melanogaster (about 2.5 Myr), indicating that the “nested” gene structure can be seen not only inDrosophila melanogaster, but also in other distantly relatedDrosophila species.     

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     

INFERRING THE RATES OF BRANCHING AND EXTINCTION FROM MOLECULAR PHYLOGENIES

Molecular techniques provide ancestral phylogenies of extant taxa with estimated branching times. Here we studied the pattern of ancestral phylogeny of extant taxa produced by branching (or cladogenesis) and extinction of taxa, assuming branching processes with time-dependent rates. (1) If the branching rate b and extinction rate c are constant, the semilog plot of the number of ancestral lineages over time is not a straight line but is curvilinear, with increasing slope toward the end, implying that ancestral phylogeny shows apparent increase in the branching rate near the present. The estimate of b and c based on nonlinear fitting is examined by computer simulation. The estimate of branching rate can be usable for a large phylogeny if b is greater than c, but the estimate of extinction rate c is unreliable because of large bias and variance. (2) Gradual decrease in the slope of the semilog plot of the number of ancestral lineages over time, as was observed in a phylogeny of bird families based on DNA hybridization data, can be explained equally well by either the decreasing branching rate or the increasing extinction rate. Infinitely many pairs of branching and extinction rates as functions of time can produce the same ancestral phylogeny. (3) An explosive branching event in the past would appear as a quick increase in the number of ancestral lineages. In contrast, mass extinction occurring in a brief period, if not accompanied by an increase in branching rate, does not produce any rapid change in the number of ancestral lineages at the time. (4) The condition in which the number of ancestral lineages of extant species changes in parallel with the actual number of species in the past is derived.     

Latitudinal clines for nucleotide polymorphisms in the Esterase 6 gene of <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Previous studies have found non-neutral patterns of nucleotide polymorphism in the promoter and coding regions of Est6 in D. melanogaster. Coding region polymorphism peaks around two closely linked replacement differences associated with the EST6-F/EST6-S allozyme polymorphism. The promoter contains two common, highly diverged haplotype groups, P1 and P7, that differentially affect Est6 expression. Allozyme studies have also revealed latitudinal clines in EST6-F and EST6-S frequencies that recur across continents. Here we analyse nucleotide polymorphisms across the promoter and the region of peak coding sequence polymorphism in 10 Australian populations along a 25° latitudinal gradient in order to examine the basis for the allozyme clines. As with the earlier studies, we find an excess of intermediate to high frequency variants in both the P1/P7 region and around the two EST6-F/EST6-S replacements in some populations. The two EST6-F/EST6-S replacement polymorphisms show latitudinal clines whereas the P1 and P7 groups of promoter haplotypes do not. However the strongest clines are for three co-segregating silent site polymorphisms in a 4 bp stretch at the 3′ end of the sequenced region. Monte Carlo simulations show that the clines for those three sites can explain all others in the data but none of the others can explain those three. Thus the allozyme clines may not reflect selection on either the P1/P7 polymorphism or the two replacements previously associated with the EST6-F/EST-S difference.     

Evolutionary History and Mode of the <Emphasis Type="Italic">amylase</Emphasis> Multigene Family in <Emphasis Type="Italic">Drosophila</Emphasis>

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 ± 1.93 and 12.38 ± 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.     

Molecular Evolution of Cytochrome c Oxidase Subunit IV: Evidence for Positive Selection in Simian Primates

Cytochrome c oxidase (COX) is a multi-subunit enzyme complex that catalyzes the final step of electron transfer through the respiratory chain on the mitochondrial inner membrane. Up to 13 subunits encoded by both the mitochondrial (subunits I, II, and III) and nuclear genomes occur in eukaryotic organisms ranging from yeast to human. Previously, we observed a high number of amino acid replacements in the human COX IV subunit compared to mouse, rat, and cow orthologues. Here we examined COX IV evolution in the two groups of anthropoid primates, the catarrhines (hominoids, cercopithecoids) and platyrrhines (ceboids), as well as one prosimian primate (lorisiform), by sequencing PCR-amplified portions of functional COX4 genes from genomic DNAs. Phylogenetic analysis of the COX4 sequence data revealed that accelerated nonsynonymous substitution rates were evident in the early evolution of both catarrhines and, to a lesser extent, platyrrhines. These accelerated rates were followed later by decelerated rates, suggesting that positive selection for adaptive amino acid replacement became purifying selection, preserving replacements that had occurred. The evidence for positive selection was especially pronounced along the catarrhine lineage to hominoids in which the nonsynonymous rate was first faster than the synonymous rate, then later much slower. The rates of three types of ``neutral DNA' nucleotide substitutions (synonymous substitutions, pseudogene nucleotide substitutions, and intron nucleotide substitutions) are similar and are consistent with previous observations of a slower rate of such substitutions in the nuclear genomes of hominoids than in the nuclear genomes of other primate and mammalian lineages. Received: 22 May 1996 / Accepted: 24 November 1996     

Rates of Transition and Transversion in Coding Sequences since the Human-Rodent Divergence

Protein-coding sequences of 337 human genes were compared with those of homologous genes from rodent (mouse or rat). A composite alignment containing 477,189 nucleotide positions was constructed, and 21,570 amino acid replacements were inferred. The rates of transitional and transversional silent substitutions in fourfold degenerate sites are estimated as 1.71 × 10^-9 and 1.22 × 10^-9 site^-1 year^-1, respectively. Rates of substitutions in replacement sites, subject to selective constraints mediated by the genetic code, are lower, but also reflect a transitional bias. The amino acid exchange rejected least often during evolution is Asp/Glu, which is fixed at 30% the rate of transversions in silent sites. The most mutable amino acids in this survey are threonine and serine; serine coded by AGY is more mutable than serine coded by TCN. A scoring matrix for evaluating amino acid similarity was derived from this study.     

Comparative Genomics of Mitochondrial DNA in Members of the Drosophila melanogaster Subgroup

In this study, a comparative genomics approach is employed to investigate the forces that shape evolutionary change in the mitochondrial DNA (mtDNA) of members of the Drosophila melanogaster subgroup. This approach facilitates differentiation of the patterns of variation resulting from processes acting at a higher level from those acting on a single gene. The mitochondrial genomes of three isofemale lines of D. simulans (siI, -II, and -III), two of D. melanogaster (Oregon R and a line from Zimbabwe), and D. mauritiana (maI and -II), and one of D. sechellia were sequenced and compared with that derived from D. yakuba. Data presented here indicate that at least three broad mechanisms shape the evolutionary dynamics of mtDNA in these taxa. The first set of mechanisms is intrinsic to the molecule. Dominant processes may be interpreted as selection for an increased rate of replication of the mtDNA molecule, biases in DNA repair, and differences in the pattern of nucleotide substitution among strands. In the genes encoded on the major strand (62% of the coding DNA) changes to or from C predominate, whereas on the minor changes to or from G predominate. The second set of mechanisms affects distinct lineages. There are evolutionary rate differences among lineages, possibly owing to population demographic changes or changes in mutational biases. This is supported by the heterogeneity found in synonymous, nonsynonymous, and silent substitutions. The third set of mechanisms differentially affects distinct genes. A maximum-likelihood sliding-window analysis detected four disjunct regions that have a significantly different nucleotide substitution process from that derived from the complete sequence. These data show the potential for comparative genomics to tease apart subtle forces that shape the evolution of DNA. Received: 30 July 1999 / Accepted: 16 March 2000     

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]     

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.     

Basal relationships in the Drosophila melanogaster species group

The Drosophila melanogaster species group is a popular model for evolutionary studies due to its morphological and ecological diversity and its inclusion of the model species D. melanogaster. However, phylogenetic relationships among major lineages within this species group remain controversial. In this report, the phylogeny of 10 species representing each of the well-supported monophyletic clades in the melanogaster group was studied using the sequences of 14 loci that together comprise 9493 nucleotide positions. Combined Bayesian analysis using gene-specific substitution models produced a 100% credible set of two trees. In the strict consensus of these trees, the ananassae subgroup branches first in the melanogaster species group, followed by the montium subgroup. The remaining lineages form a monophyletic clade in which D. ficusphila and D. elegans branch first, followed by D. biarmipes, D. eugracilis, and the melanogaster subgroup. This strongly supported phylogeny resolves most basal relationships in the melanogaster species group, and provides a framework that can be extended in the future to encompass more species.     

Parallel evolution of site‐specific changes in divergent caribou lineages

The parallel evolution of phenotypes or traits within or between species provides important insight into the basic mechanisms of evolution. Genetic and genomic advances have allowed investigations into the genetic underpinnings of parallel evolution and the independent evolution of similar traits in sympatric species. Parallel evolution may best be exemplified among species where multiple genetic lineages, descended from a common ancestor, colonized analogous environmental niches, and converged on a genotypic or phenotypic trait. Modern North American caribou (Rangifer tarandus) originated from three ancestral sources separated during the Last Glacial Maximum (LGM): the Beringian–Eurasian lineage (BEL), the North American lineage (NAL), and the High Arctic lineage (HAL). Historical introgression between the NAL and the BEL has been found throughout Ontario and eastern Manitoba. In this study, we first characterized the functional differentiation in the cytochrome‐b (cytB) gene by identifying nonsynonymous changes. Second, the caribou lineages were used as a direct means to assess site‐specific parallel changes among lineages. There was greater functional diversity within the NAL despite the BEL having greater neutral diversity. The patterns of amino acid substitutions occurring within different lineages supported the parallel evolution of cytB amino acid substitutions suggesting different selective pressures among lineages. This study highlights the independent evolution of identical amino acid substitutions within a wide‐ranging mammal species that have diversified from different ancestral haplogroups and where ecological niches can invoke parallel evolution.     

Evidence for interspecific transfer of the transposable element mariner betweenDrosophila andZaprionus

Summary The transposable element mariner occurs widely in themelanogaster species group ofDrosophila. However, in drosophilids outside of themelanogaster species group, sequences showing strong DNA hybridization with mariner are found only in the genusZaprionus. the mariner sequence obtained fromZaprionus tuberculatus is 97% identical with that fromDrosophila mauritiana, a member of themelanogaster species subgroup, whereas a mariner sequence isolated fromDrosophila tsacasi is only 92% identical with that fromD. mauritiana. BecauseD. tsacasi is much more closely related toD. mauritiana than isZaprionus, the presence of mariner inZaprionus may result from horizontal transfer. In order to confirm lack of a close phylogenetic relationship between the genusZaprionus and themelanogaster species group, we compared the alcohol dehydrogenase (Adh) sequences among these species. The results show that the coding region of Adh is only 82% identical betweenZ. tuberculatus andD. mauritiana, as compared with 90% identical betweenD. tsacasi andD. mauritiana. Furthermore, the mariner gene phylogeny obtained by maximum likelihood and maximum parsimony analyses is discordant with the species phylogeny estimated by using the Adh genes. The only inconsistency in the mariner gene phylogeny is in the placement of theZaprionus mariner sequence, which clusters with mariner fromDrosophila teissieri andDrosophila yakuba in themelanogaster species subgroup. These results strongly suggest horizontal transfer.     

Parallel amino acid replacements in the rhodopsins of the rockfishes (Sebastes spp.) associated with shifts in habitat depth

Among various groups of fishes, a shift in peak wavelength sensitivity has been correlated with changes in their photic environments. The genus Sebastes is a radiation of marine fish species that inhabit a wide range of depths from intertidal to over 600 m. We examined 32 species of Sebastes for evidence of adaptive amino acid substitution at the rhodopsin gene. Fourteen amino acid positions were variable among these species. Maximum likelihood analyses identify several of these to be targets of positive selection. None of these correspond to previously identified critical amino acid sites, yet they may in fact be functionally important. The occurrence of independent parallel changes at certain amino acid positions reinforces this idea. Reconstruction of habitat depths of ancestral nodes in the phylogeny suggests that shallow habitats have been colonized independently in different lineages. The evolution of rhodopsin appears to be associated with changes in depth, with accelerated evolution in lineages that have had large changes in depth.     

Nucleotide and repeat length variation at the nonA gene of the Drosophila virilis group species and its effects on male courtship song

Genes found to affect male courtship song characters in Drosophila melanogaster are good candidates when tracing genes responsible for species-specific songs in other Drosophila species. It has previously been shown that Thr–Gly repeat length variation at the period gene affects song traits in D. melanogaster, which gives the repetitive regions a special interest. In this work, we have characterised the patterns of nucleotide variation for gene regions containing two Gly and one Gln–Ala repeat in another D. melanogaster song gene, no-on-transient A, in D. virilis group species. The levels of nucleotide variability in D. virilis nonA were similar to those found for other genes of the species, and the gene sequences showed no signs of deviation from neutrality. The Gly 2 repeat preceding the central domain of the gene exhibited length variation, which did not, however, correlate with song variation either within D. virilis or between the species of D. virilis group. The Gly 3 repeat located on the other side of the central domain showed amino acid divergence parallel to the consensus phylogeny of the D. virilis group species. The species of the virilis subgroup having Asn after the first three glycines in this repeat have simple songs with no species-specificity, while the species of the montana subgroup having two Gly or Asn–Ser in this site have unique courtship songs. Amino acid differences between the species in this repeat may, however, reflect species phylogeny rather than have an effect on song divergence per se.     

Lessons from the deep: mechanisms behind diversification of eukaryotic protein complexes

Genetic variation is the major mechanism behind adaptation and evolutionary change. As most proteins operate through interactions with other proteins, changes in protein complex composition and subunit sequence provide potentially new functions. Comparative genomics can reveal expansions, losses and sequence divergence within protein-coding genes, but in silico analysis cannot detect subunit substitutions or replacements of entire protein complexes. Insights into these fundamental evolutionary processes require broad and extensive comparative analyses, from both in silico and experimental evidence. Here, we combine data from both approaches and consider the gamut of possible protein complex compositional changes that arise during evolution, citing examples of complete conservation to partial and total replacement by functional analogues. We focus in part on complexes in trypanosomes as they represent one of the better studied non-animal/non-fungal lineages, but extend insights across the eukaryotes by extensive comparative genomic analysis. We argue that gene loss plays an important role in diversification of protein complexes and hence enhancement of eukaryotic diversity.     

Maximum likelihood inference of protein phylogeny and the origin of chloroplasts

Summary A maximum likelihood method for inferring protein phylogeny was developed. It is based on a Markov model that takes into account the unequal transition probabilities among pairs of amino acids and does not assume constancy of rate among different lineages. Therefore, this method is expected to be powerful in inferring phylogeny among distantly related proteins, either orthologous or parallogous, where the evolutionary rate may deviate from constancy. Not only amino acid substitutions but also insertion/deletion events during evolution were incorporated into the Markov model. A simple method for estimating a bootstrap probability for the maximum likelihood tree among alternatives without performing a maximum likelihood estimation for each resampled data set was developed. These methods were applied to amino acid sequence data of a photosynthetic membrane protein,psbA, from photosystem II, and the phylogeny of this protein was discussed in relation to the origin of chloroplasts.     

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.