Molecular Evolution between Drosophila Melanogaster and D. Simulans: Reduced Codon Bias, Faster Rates of Amino Acid Substitution, and Larger Proteins in D. Melanogaster

Both natural selection and mutational biases contribute to variation in codon usage bias within Drosophila species. This study addresses the cause of codon bias differences between the sibling species, Drosophila melanogaster and D. simulans. Under a model of mutation-selection-drift, variation in mutational processes between species predicts greater base composition differences in neutrally evolving regions than in highly biased genes. Variation in selection intensity, however, predicts larger base composition differences in highly biased loci. Greater differences in the G+C content of 34 coding regions than 46 intron sequences between D. melanogaster and D. simulans suggest that D. melanogaster has undergone a reduction in selection intensity for codon bias. Computer simulations suggest at least a fivefold reduction in N(e)s at silent sites in this lineage. Other classes of molecular change show lineage effects between these species. Rates of amino acid substitution are higher in the D. melanogaster lineage than in D. simulans in 14 genes for which outgroup sequences are available. Surprisingly, protein sizes are larger in D. melanogaster than in D. simulans in the 34 genes compared between the two species. A substantial fraction of silent, replacement, and insertion/deletion mutations in coding regions may be weakly selected in Drosophila.     

Comparative analysis of transposable elements in the melanogaster subgroup sequenced genomes

Transposable elements (TEs) are indwelling components of genomes, and their dynamics have been a driving force in genome evolution. Although we now have more information concerning their amounts and characteristics in various organisms, we still have little data from overall comparisons of their sequences in very closely-related species. While the Drosophila melanogaster genome has been extensively studied, we have only limited knowledge regarding the precise TE sequences in the genomes of the related species Drosophila simulans, Drosophila sechellia and Drosophila yakuba. In this study we analyzed the number and structure of TE copies in the sequenced genomes of these four species. Our findings show that, unexpectedly, the number of TE insertions in D. simulans is greater than that in D. melanogaster, but that most of the copies in D. simulans are degraded and in small fragments, as in D. sechellia and D. yakuba. This suggests that all three species were invaded by numerous TEs a long time ago, but have since regulated their activity, as the present TE copies are degraded, with very few full-length elements. In contrast, in D. melanogaster, a recent activation of TEs has resulted in a large number of almost-identical TE copies. We have detected variants of some TEs in D. simulans and D. sechellia, that are almost identical to the reference TE sequences in D. melanogaster, suggesting that D. melanogaster has recently been invaded by active TE variants from the other species. Our results indicate that the three species D. simulans, D. sechellia, and D. yakuba seem to be at a different stage of their TE life cycle when compared to D. melanogaster. Moreover, we show that D. melanogaster has been invaded by active TE variants for several TE families likely to come from D. simulans or the ancestor of D. simulans and D. sechellia. The numerous horizontal transfer events implied to explain these results could indicate introgression events between these species.     

Adaptive evolution of recently duplicated accessory gland protein genes in desert Drosophila

The relationship between animal mating system variation and patterns of protein polymorphism and divergence is poorly understood. Drosophila provides an excellent system for addressing this issue, as there is abundant interspecific mating system variation. For example, compared to D. melanogaster subgroup species, repleta group species have higher remating rates, delayed sexual maturity, and several other interesting differences. We previously showed that accessory gland protein genes (Acp's) of Drosophila mojavensis and D. arizonae evolve more rapidly than Acp's in the D. melanogaster subgroup and that adaptive Acp protein evolution is likely more common in D. mojavensis/D. arizonae than in D. melanogaster/D. simulans. These findings are consistent with the idea that greater postcopulatory selection results in more adaptive evolution of seminal fluid proteins in the repleta group flies. Here we report another interesting evolutionary difference between the repleta group and the D. melanogaster subgroup Acp's. Acp gene duplications are present in D. melanogaster, but their high sequence divergence indicates that the fixation rate of duplicated Acp's has been low in this lineage. Here we report that D. mojavensis and D. arizonae genomes contain several very young duplicated Acp's and that these Acp's have experienced very rapid, adaptive protein divergence. We propose that rapid remating of female desert Drosophila generates selection for continuous diversification of the male Acp complement to improve male fertilization potential. Thus, mating system variation may be associated with adaptive protein divergence as well as with duplication of Acp's in Drosophila.     

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.     

Wake up of transposable elements following Drosophila simulans worldwide colonization.

Transposable elements (TEs) make up around 10%-15% of the Drosophila melanogaster genome, but its sibling species Drosophila simulans carries only one third as many such repeat sequences. We do not, however, have an overall view of copy numbers of the various classes of TEs (long terminal repeat [LTR] retrotransposons, non-LTR retrotransposons, and transposons) in genomes of natural populations of both species. We analyzed 34 elements in individuals from various natural populations of these species. We show that D. melanogaster has higher average chromosomal insertion site numbers per genome than D. simulans for all TEs except five. The LTR retrotransposons gypsy, ZAM, and 1731 and the transposon bari-1 present similar low copy numbers in both species. The transposon hobo has a large number of insertion sites, with significantly more sites in D. simulans. High variation between populations in number of insertion sites of some elements of D. simulans suggests that these elements can invade the genome of the entire species starting from a local population. We propose that TEs in the D. simulans genome are being awakened and amplified as they had been a long time ago in D. melanogaster.     

Molecular Evolution of Drosophila Metallothionein Genes

The metallothionein genes of Drosophila melanogaster, Mtn and Mto, may play an important role in heavy metal detoxification. Several different tandem duplications of Mtn have been shown to increase cadmium and copper tolerance, as well as Mtn expression. In order to investigate the possibility of increased selection for duplications of these genes in natural populations exposed to high levels of heavy metals, we compared the frequencies of such duplications among flies collected from metal-contaminated and non-contaminated orchards in Pennsylvania, Tennessee and Georgia. Restriction enzyme analysis was used to screen 1666 wild third chromosomes for Mtn duplications and a subset (327) of these lines for Mto duplications. The frequency of pooled Mtn duplications found ranged from 0% to 20%, and was not significantly higher at the contaminated sites. No Mto duplications were identified. Estimates of sequence diversity at the Mtn locus among a subsample (92) of the duplication survey were obtained using four-cutter analysis. This analysis revealed a low level of polymorphism, consistent with both selection at the Mtn locus, and a fairly recent origin for the duplications. To further examine this hypothesis, we sequenced an Mtn allele of Drosophila simulans and measured the amount of nucleotide sequence divergence between D. simulans and its sibling species D. melanogaster. The levels of silent nucleotide polymorphism and divergence in the Mtn region were compared with those in the Adh region, using the neutrality test of R.R. Hudson, M. Kreitman and M. Aguadé.     

Ancestral polymorphism and recent invasion of transposable elements in Drosophila species

ABSTRACT: BACKGROUND: During the evolutionary history of transposable elements, some processes, such as ancestral polymorphisms and horizontal transfer of sequences between species, can produce incongruences in phylogenies. We investigated the evolutionary history of the transposable elements Bari and 412 in the sequenced genomes of the Drosophila melanogaster group and in the sibling species D. melanogaster and D. simulans using traditional phylogenetic and network approaches. RESULTS: The maximum likelihood (ML) phylogenetic analyses revealed incongruences and unresolved relationships for both the Bari and 412 elements. The DNA transposon Bari within the D. ananassae genome is more closely related to the element of the melanogaster complex than to the sequence in D. erecta, which is inconsistent with the species phylogeny. Divergence analysis and the comparison of the rate of synonymous substitutions per synonymous site of the Bari and host gene sequences explain the incongruence as an ancestral polymorphism inherited stochastically by the derived species. Unresolved relationships were observed in the ML phylogeny of both elements involving D. melanogaster, D. simulans and D. sechellia. A network approach was used to attempt to resolve these relationships. The resulting tree suggests recent transfers of both elements between D. melanogaster and D. simulans. The divergence values of the elements between these species support this conclusion. CONCLUSIONS: We showed that an ancestral polymorphism and recent invasion of genomes due to introgression or horizontal transfer between species occurred during the evolutionary history of the Bari and 412 elements in the melanogaster group. These invasions likely occurred in Africa during the Pleistocene, before the worldwide expansion of D. melanogaster and D. simulans.     

Transposable Element Dynamics in Two Sibling Species: Drosophila Melanogaster and Drosophila Simulans

Transposable elements (TEs) in the two sibling species, Drosophila melanogaster and D. simulans, differ considerably in amount and dynamics, with D. simulans having a smaller amount of TEs than D. melanogaster. Several hypotheses have been proposed to explain these differences, based on the evolutionary history of the two species, and claim differences either in the effective size of the population or in genome characteristics. Recent data suggest, however, that the higher amount of TEs in D. melanogaster could be associated with the worldwide invasion of D. melanogaster a long time ago while D. simulans is still under the process of such geographical spread. Stresses due to new environmental conditions and crosses between migrating populations could explain the mobilization of TEs while the flies colonize. Colonization and TE mobilization may be strong evolutionary forces that have shaped and are still shaping the eukaryote genomes.     

Mitochondrial pseudogenes in the nuclear genomes of Drosophila

Mitochondrial pseudogenes in nuclear chromosomes (numts) have been detected in the genomes of a diverse range of eukaryotic species. However, the numt content of different genomes and their properties is not uniform, and study of these differences provides insight into the mechanisms and dynamics of genome evolution in different organisms. In the genus Drosophila, numts have previously only been identified on a genome-wide scale in the melanogaster subgroup. The present study extends the identification to 11 species of the Drosophila genus. We identify a total of 302 numts and show that the numt complement is highly variable in Drosophilids, ranging from just 4 in D. melanogaster to 67 in D. willistoni, broadly correlating with genome size. Many numts have undergone large-scale rearrangements in the nucleus, including interruptions, inversions, deletions and duplications of sequence of variable size. Estimating the age of the numts in the nucleus by phylogenetic tree reconstruction reveals the vast majority of numts to be recent gains, 90% having arisen on terminal branches of the species tree. By identifying paralogs and counting duplications among the extant numts we estimate that 23% of extant numts arose through post-insertion duplications. We estimate genus average rates of insertion of 0.75 per million years, and a duplication rate of 0.010 duplications per numt per million years.     

Detection of gene duplications and block duplications in eukaryotic genomes

Several eukaryotic genomes have been completely sequenced and this provides an opportunity to investigate the extent and characteristics (e.g., single gene duplication, block duplication, etc.) of gene duplication in a genome. Detecting duplicate genes in a genome, however, is not a simple problem because of several complications such as domain shuffling, the existence of isoforms derived from alternative splicing, and annotational errors in the databases. We describe a method for overcoming these difficulties and the extents of gene duplication in the genomes of Drosophila melanogaster, Caenorhabditis elegans, and yeast inferred from this method. We also describe a method for detecting block duplications in a genome. Application of this method showed that block duplication is a common phenomenon in both yeast and nematode. The patterns of block duplication in the two species are, however, markedly different. Yeast shows much more extensive block duplication than nematode, with some chromosomes having more than 40% of the duplications derived from block duplications. Moreover, in yeast the majority of block duplications occurred between chromosomes, while in nematode most block duplications occurred within chromosomes.     

Incompatibility between X chromosome factor and pericentric heterochromatic region causes lethality in hybrids between Drosophila melanogaster and its sibling species

The Dobzhansky-Muller model posits that postzygotic reproductive isolation results from the evolution of incompatible epistatic interactions between species: alleles that function in the genetic background of one species can cause sterility or lethality in the genetic background of another species. Progress in identifying and characterizing factors involved in postzygotic isolation in Drosophila has remained slow, mainly because Drosophila melanogaster, with all of its genetic tools, forms dead or sterile hybrids when crossed to its sister species, D. simulans, D. sechellia, and D. mauritiana. To circumvent this problem, we used chromosome deletions and duplications from D. melanogaster to map two hybrid incompatibility loci in F(1) hybrids with its sister species. We mapped a recessive factor to the pericentromeric heterochromatin of the X chromosome in D. simulans and D. mauritiana, which we call heterochromatin hybrid lethal (hhl), which causes lethality in F(1) hybrid females with D. melanogaster. As F(1) hybrid males hemizygous for a D. mauritiana (or D. simulans) X chromosome are viable, the lethality of deficiency hybrid females implies that a dominant incompatible partner locus exists on the D. melanogaster X. Using small segments of the D. melanogaster X chromosome duplicated onto the Y chromosome, we mapped a dominant factor that causes hybrid lethality to a small 24-gene region of the D. melanogaster X. We provide evidence suggesting that it interacts with hhl(mau). The location of hhl is consistent with the emerging theme that hybrid incompatibilities in Drosophila involve heterochromatic regions and factors that interact with the heterochromatin.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

A Snapshot of Histone Modifications within Transposable Elements in Drosophila Wild Type Strains

Transposable elements (TEs) are a major source of genetic variability in genomes, creating genetic novelty and driving genome evolution. Analysis of sequenced genomes has revealed considerable diversity in TE families, copy number, and localization between different, closely related species. For instance, although the twin species Drosophila melanogaster and D. simulans share the same TE families, they display different amounts of TEs. Furthermore, previous analyses of wild type derived strains of D. simulans have revealed high polymorphism regarding TE copy number within this species. Several factors may influence the diversity and abundance of TEs in a genome, including molecular mechanisms such as epigenetic factors, which could be a source of variation in TE success. In this paper, we present the first analysis of the epigenetic status of four TE families (roo, tirant, 412 and F) in seven wild type strains of D. melanogaster and D. simulans. Our data shows intra- and inter-specific variations in the histone marks that adorn TE copies. Our results demonstrate that the chromatin state of common TEs varies among TE families, between closely related species and also between wild type strains.     

The correlation between intron length and recombination in drosophila. Dynamic equilibrium between mutational and selective forces

Intron length is negatively correlated with recombination in both Drosophila melanogaster and humans. This correlation is not likely to be the result of mutational processes alone: evolutionary analysis of intron length polymorphism in D. melanogaster reveals equivalent ratios of deletion to insertion in regions of high and low recombination. The polymorphism data do reveal, however, an excess of deletions relative to insertions (i.e., a deletion bias), with an overall deletion-to-insertion events ratio of 1.35. We propose two types of selection favoring longer intron lengths. First, the natural mutational bias toward deletion must be opposed by strong selection in very short introns to maintain the minimum intron length needed for the intron splicing reaction. Second, selection will favor insertions in introns that increase recombination between mutations under the influence of selection in adjacent exons. Mutations that increase recombination, even slightly, will be selectively favored because they reduce interference among selected mutations. Interference selection acting on intron length mutations must be very weak, as indicated by frequency spectrum analysis of Drosophila intron length polymorphism, making the equilibrium for intron length sensitive to changes in the recombinational environment and population size. One consequence of this sensitivity is that the advantage of longer introns is expected to decrease inversely with the rate of recombination, thus leading to a negative correlation between intron length and recombination rate. Also in accord with this model, intron length differs between closely related Drosophila species, with the longest variant present more often in D. melanogaster than in D. simulans. We suggest that the study of the proposed dynamic model, taking into account interference among selected sites, might shed light on many aspects of the comparative biology of genome sizes including the C value paradox.     

Is genome size influenced by colonization of new environments in dipteran species?

Genome size differences are usually attributed to the amplification and deletion of various repeated DNA sequences, including transposable elements (TEs). Because environmental changes may promote modifications in the amount of these repeated sequences, it has been postulated that when a species colonizes new environments this could be followed by an increase in its genome size. We tested this hypothesis by estimating the genome size of geographically distinct populations of Drosophila ananassae, Drosophila malerkotliana, Drosophila melanogaster, Drosophila simulans, Drosophila subobscura, and Zaprionus indianus, all of which have known colonization capacities. There was no strong statistical differences between continents for most species. However, we found that populations of D. melanogaster from east Africa have smaller genomes than more recent populations. For species in which colonization is a recent event, the differences between genome sizes do not thus seem to be related to colonization history. These findings suggest either that genome size is seldom modified in a significant way during colonization or that it takes time for genome size of invading species to change significantly.     

Rapid sequence turnover at an intergenic locus in Drosophila

Closely related species of Drosophila tend to have similar genome sizes. The strong imbalance in favor of small deletions relative to insertions implies that the unconstrained DNA in Drosophila is unlikely to be passively inherited from even closely related ancestors, and yet most DNA in Drosophila genomes is intergenic and potentially unconstrained. In an attempt to investigate the maintenance of this intergenic DNA, we studied the evolution of an intergenic locus on the fourth chromosome of the Drosophila melanogaster genome. This 1.2-kb locus is marked by two distinct, large insertion events: a nuclear transposition of a mitochondrial sequence and a transposition of a nonautonomous DNA transposon DNAREP1_DM. Because we could trace the evolutionary histories of these sequences, we were able to reconstruct the length evolution of this region in some detail. We sequenced this locus in all four species of the D. melanogaster species complex: D. melanogaster, D. simulans, D. sechellia, and D. mauritiana. Although this locus is similar in size in these four species, less than 10% of the sequence from the most recent common ancestor remains in D. melanogaster and all of its sister species. This region appears to have increased in size through several distinct insertions in the ancestor of the D. melanogaster species complex and has been shrinking since the split of these lineages. In addition, we found no evidence suggesting that the size of this locus has been maintained over evolutionary time; these results are consistent with the model of a dynamic equilibrium between persistent DNA loss through small deletions and more sporadic DNA gain through less frequent but longer insertions. The apparent stability of genome size in Drosophila may belie very rapid sequence turnover at intergenic loci.     

Wanderings of Hobo: A Transposon in Drosophila Melanogaster and its Close Relatives

The transposon hobo is present in the genomes of Drosophila melanogaster and Drosophila simulans (and D. mauritiana and probably D. sechellia, based on Southern blots) as full-size elements and internally deleted copies. The full-size melanogaster, simulans and mauritiana hobo elements are 99.9% identical at the DNA sequence level, and internally deleted copies in these species essentially differ only in having deletions. In addition to these, hobo-related sequences are present and detectable with a hobo probe in all these species. Those in D. melanogaster are 86-94% identical to the canonical hobo, but with many indels. We have sequenced one that appears to be inserted in heterochromatin (GenBank Acc. No. AF520587). It is 87.6% identical to the canonical hobo, but quite fragmented by indels, with remnants of other transposons inserted in and near it, and clearly is defunct. Numerous similar elements are found in the sequenced D. melanogaster genome. It has recently been shown that some are fixed in the euchromatic genome, but it is probable that still more reside in heterochromatic regions not included in the D. melanogaster genome database. They are probably all relics of an earlier introduction of hobo into the ancestral species. There appear to have been a minimum of two introductions of hobo into the melanogaster subgroup, and more likely three, two ancient and one quite recent. The recent introduction of hobo was probably followed by transfers between the extant species (whether 'horizontally' or by infrequent interspecific hybridization).     

基因倍增研究进展

基因倍增是指DNA片段在基因组中复制出一个或更多的拷贝,这种DNA片段可以是一小段基因组序列、整条染色体,甚至是整个基因组。基因倍增是基因组进化最主要的驱动力之一,是产生具有新功能的基因和进化出新物种的主要原因之一。本文综述了脊椎动物、模式植物和酵母在进化过程中基因倍增研究领域的最新进展,并讨论了基因倍增研究的发展方向。