Genetic differentiation and adaptive evolution at reproductive loci in incipient Drosophila species

Accessory gland proteins (Acps) are part of the seminal fluid of male Drosophila flies. Some Acps have exceptionally high evolutionary rates and evolve under positive selection. Proper interactions between Acps and female reproductive molecules are essential for fertilization. These observations lead to suggestions that fast evolving Acps could be involved in speciation by promoting reproductive incompatibilities between emerging species. To test this hypothesis, we used population genetics data for three sibling species: D. mayaguana, D. parisiena and D. straubae. The latter two species are morphologically very similar and show only incipient reproductive isolation. This system allowed us to examine Acp evolution at different time frames with respect to speciation and reproductive isolation. Comparing data of 14 Acp loci with data obtained for other genomic regions, we found that some Acps show extraordinarily high levels of divergence between D. mayaguana and its two sister species D. parisiena and D. straubae. This divergence was likely driven by adaptive evolution at several loci. No fixed nucleotide differences were found between D. parisiena and D. straubae, however. Nevertheless, some Acp loci did show significant differentiation between these species associated with signs of positive selection; these loci may be involved in this early phase of the speciation process.     

SHORT INDELS ARE SUBJECT TO INSERTION‐BIASED GENE CONVERSION

Recombination between homologous loci is accompanied by formation of heteroduplexes. Repairing mismatches in heteroduplexes often leads to single nucleotide substitutions in a process known as gene conversion. Gene conversion was shown to be GC‐biased in different organisms; that is, a W(A or T)→S(G or C) substitution is more likely in this process than a S→W substitution. Here, we show that the insertion/deletion ratio for short noncoding indels that reach fixation between species is positively correlated with the recombination rate in Drosophila melanogaster, Homo sapiens, and Saccharomyces cerevisiae. This correlation is both due to an increase of the fixation rate of insertions and decrease of the fixation rate of deletions in regions of high recombination. Whole‐genome data on indel polymorphism and divergence in D. melanogaster rule out mutation biases and selection as the cause of this trend, pointing to insertion‐biased gene conversion as the most likely explanation. The bias toward insertions is the strongest for single‐nucleotide indels, and decreases with indel length. In regions of high recombination rate this bias leads to an up to ～5‐fold excess of fixed short insertions over deletions, and substantially affects the evolution of DNA segments.     

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.     

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]     

Mitochondrial DNA evolution in theobscura species subgroup ofDrosophila

Summary Mitochondrial DNA (mtDNA) restriction site maps for nine species of theDrosophila obscura subgroup and forDrosophila melanogaster were established. Taking into account all restriction enzymes (12) and strains (45) analyzed, a total of 105 different sites were detected, which corresponds to a sample of 3.49% of the mtDNA genome. Based on nucleotide divergences, two phylogenetic trees were constructed assuming either constant or variable rates of evolution. Both methods led to the same relationships. Five differentiated clusters were found for theobscura subgroup species, one Nearctic, represented byDrosophila pseudoobscura, and four Palearctic, two grouping the related triads of speciesDrosophila subobscura, Drosophila madeirensis, Drosophila guanche, andDrosophila ambigua, Drosophila obscura, Drosophila subsilvestris, and two more represented by one species each,Drosophila bifasciata, andDrosophila tristis. The different Palearctic clusters are as distant between themselves as with the Nearctic one. For the related speciesD. subobscura, D. madeirensis, andD. guanche, the pairD. subobscura-D. madeirensis is the closest one. The relationships found by nucleotide divergence were confirmed by differences in mitochondrial genome size, with related species sharing similar genome lengths and differing from the distant ones. The total mtDNA size range for theobscura subgroup species was from 15.5 kb forD. pseudoobscura to 17.1 forD. tristis.     

<Emphasis Type="Italic">Saturnispora serradocipensis</Emphasis> sp. nov. and <Emphasis Type="Italic">Saturnispora gosingensis</Emphasis> sp. nov., two ascomycetous yeasts from ephemeral habitats

Two novel ascomycetous yeast species, Saturnispora serradocipensis and Saturnispora gosingensis, were isolated from leaf detritus in a tropical stream of Southeastern Brazil and a mushroom collected in Taiwan, respectively. Analysis of the D1/D2 domains of the large-subunit of the rRNA gene of these strains showed that these species are related to Saturnispora hagleri, their closest relative. Saturnispora serradocipensis and S. gosingensis differed from S. hagleri, respectively, by seven nucleotide substitutions and two indels and three nucleotide substitutions and three indels in D1/D2 rRNA sequences. The two new species differ from each another by four nucleotide substitutions and one indel in D1/D2 rRNA sequences. However, the ITS sequences of S. serradocipensis, S. gosingensis and S. hagleri were quite divergent, showing that they are genetically separate species. The type strain of S. serradocipensis is UFMG-DC-198^T (=CBS 11756^T = NRRL Y-48717^T), and of S. gosingensis GA4M05^T is (CBS 11755^T = NRRL Y-48718^T).     

Sequence and evolution of the Drosophila pseudoobscura glycerol-3-phosphate dehydrogenase locus

The Gpdh genomic region has been cloned and sequenced in Drosophila pseudoobscura. A total of 6.8 kb of sequence was obtained, encompassing all eight exons of the gene. The exons have been aligned with the sequence from D. melanogaster, and the rates of synonymous and nonsynonymous substitution have been compared to those of other genes sequenced in these two species. Gpdh has the lowest rate of nonsynonymous substitution yet seen in genes sequenced in both D. pseudoobscura and D. melanogaster. No insertion/deletion events were observed, and the overall architecture of the gene (i.e., intron sites, etc.) is conserved. An interesting amino acid reversal was noted between the D. melanogaster Fast allele and the D. pseudoobscura gene.     

The pattern of insertion/deletion polymorphism in Arabidopsis thaliana

Little is known about variation of nucleotide insertion/deletions (indels) within species. In Arabidopsis thaliana, we investigated indel polymorphism patterns between two genome sequences and among 96 accessions at 1215 loci. Our study identified patterns in the variation of indel density, size, GC content and distribution, and a correlation between indels and substitutions. We found that the GC content in indel sequences was lower than that in non-indel sequences and that indels typically occur in regions with lower GC content. Patterns of indel frequency distribution among populations were more consistent with neutral expectation than substitution patterns. We also found that the local level of substitutions is positively correlated with indel density and negatively correlated with their distance to the closed indel, suggesting that indels play an important role in nucleotide variation.     

Positive selection driving the evolution of a gene of male reproduction, Acp26Aa, of Drosophila: II. Divergence versus polymorphism

The evolution of the gene for a male ejaculatory protein, Acp26Aa, has been shown to be driven by positive selection when nonsibling species in the Drosophila melanogaster subgroup are compared. To know if selection has been operating in the recent past and to understand the details of its dynamics, we obtained DNA sequences of Acp26Aa and the nearby Acp26Ab gene from 39 D. melanogaster chromosomes. Together with the 10 published sequences, we analyzed 49 sequences from five populations in four continents. The southern African population is somewhat differentiated from all other populations, but its nucleotide diversity is lower at these two loci. We find the following results for Acp26Aa: (1) The R: S (replacement : silent changes) ratio is significantly higher in the between-species comparisons than in the within-species data by the McDonald and Kreitman test. Positive selection is probably responsible for the excess of amino acid replacements between species. (2) However, within-species nucleotide diversity is high. Neither the Tajima test nor the Fu and Li test indicates a reduction in nucleotide diversity due to positive selection in the recent past. (3) The newly derived nucleotides in D. melanogaster are at high frequency significantly more often than predicted by the neutral equilibrium. Since the nearby Acp26Ab gene does not show these patterns, these observations cannot be attributed to the characteristics of this chromosomal region. We suggest that positive selection is active, but may be weak, for each amino acid change in the Acp26Aa gene.      

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     

Connecting recombination,nucleotide diversity,and species divergence in Drosophila

The association between recombination rate and nucleotide diversity provides compelling evidence for the action of natural selection across much of the Drosophila melanogaster genome. This conclusion is further supported by the lack of association between recombination rate and nucleotide divergence between species. However, studies of other species, including other Drosophila, have not always yielded the same results. Our recent study measured these parameters within the D. pseudoobscura species group using next-generation sequencing and high-throughput genotyping technologies. We documented fine-scale variation in crossover rate within D. pseudoobscura, and we observed that crossover variation was strongly associated with nucleotide diversity only when measured at a fine-scale. We also observed associations between crossover rate and sequence differences between D. pseudoobscura and its close relatives. These latter associations could have been driven in part by mutagenic effects associated with double-strand break repair, but we cannot exclude the possibility that it results primarily from shared ancestral polymorphisms. Overall, this work strongly underscores the importance of scale in testing for associations of recombination rate with other parameters, and it brings us one small step closer to understanding the role of natural selection and other evolutionary forces in shaping divergence among genomes.     

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.     

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.     

Sources and predictors of resolvable indel polymorphism assessed using rice as a model

The principal sources of genetic variation that can be assayed with restriction enzymes are base substitutions and insertions/deletions (indels). The likelihood of detecting indels as restriction fragment length polymorphisms (RFLPs) is determined by the size and frequency of the indels, and the ability to resolve small indels as RFLPs is limited by the distribution of restriction fragment sizes. In this study, we use aligned sequences from the indica and japonica subspecies of rice ( Oryza sativa L.) to quantify and compare the ability of restriction enzymes to detect indels. We look specifically at two abundant transposable element-derived indel sources: miniature inverted repeat transposable elements (MITEs) and long terminal repeat (LTR) retroelements. From this analysis we conclude that indels rather than base substitutions are the prevailing source of the polymorphism detected in rice. We show that, although MITE derived indels are more abundant than LTR-retroelement derived indels, LTR-retroelements have a greater capacity to generate visible restriction fragment length polymorphism because of their larger size. We find that the variation in the detectability of indels among restriction enzymes can be explained by differences in the frequency and dispersion of their restriction sites in the genome. The parameters that describe the fragment size distributions obtained with the restriction enzymes are highly correlated across the sequenced genomes of rice, Arabidopsis and human, with the exception of some extreme deviations in frequency for particular recognition sequences corresponding to variations in the levels and modes of DNA methylation in the three disparate organisms. Thus, we can predict the relative ability of a restriction enzyme to detect indels derived from a specific source based on the distribution of restriction fragment sizes, even when this is estimated for a distantly related genome.Electronic Supplementary Material Supplementary Material is available in the online version of this article at Communicated by M.-A. Grandbastien     

Molecular evolution of X-linked accessory gland proteins in Drosophila pseudoobscura

In Drosophila melanogaster and Drosophila simulans, positive Darwinian selection drives high rates of evolution of male reproductive genes, and accessory gland proteins (Acps) in particular. Here, we tested whether 13 X-linked male-specific genes, 4 Acps and 9 non-Acps, are under selective forces in the Drosophila pseudoobscura species group, much as those in the D. melanogaster group. We observed a statistically significant correlation in relative rates of nonsynonymous evolution between the two species groups tested. One Acp examined had a higher rate of nonsynonymous substitution than predicted by a neutral model in both species groups, suggesting its divergence was driven by positive Darwinian selection. To further test for the signature of selection, we examined polymorphism of three Acps within D. pseudoobscura. From this test, no Acp individually bore the signature of positive selection, but the 3 Acps together possessed an excess of nonsynonymous differences between species, relative to polymorphism within species. We conclude that faster evolution of Acps in the D. pseudoobscura group appears to be driven by positive selection, as previously suggested in the D. melanogaster group.     

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     

Orthology, Function and Evolution of Accessory Gland Proteins in the Drosophila repleta Group

The accessory gland proteins (Acps) of Drosophila have become a model for the study of reproductive protein evolution. A major step in the study of Acps is to identify biological causes and consequences of the observed patterns of molecular evolution by comparing species groups with different biology. Here we characterize the Acp complement of Drosophila mayaguana, a repleta group representative. Species of this group show important differences in ecology and reproduction as compared to other Drosophila. Our results show that the extremely high rates of Acp evolution previously found are likely to be ubiquitous among species of the repleta group. These evolutionary rates are considerably higher than the ones observed in other Drosophila groups' Acps. This disparity, however, is not accompanied by major differences in the estimated number of Acps or in the functional categories represented as previously suggested. Among the genes expressed in accessory glands of D. mayaguana almost half are likely products of recent duplications. This allowed us to test predictions of the neofunctionalization model for gene duplication and paralog evolution in a more or less constrained timescale. We found that positive selection is a strong force in the early divergence of these gene pairs.     

Patterns of DNA-Sequence Divergence Between <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. pseudoobscura</Emphasis>

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d _N) and codon usage bias (F _op), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.     

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.).