- Research article
- Open Access
Molecular evolution of sex-biased genes in the Drosophila ananassaesubgroup
© Grath et al; licensee BioMed Central Ltd. 2009
Received: 3 July 2009
Accepted: 16 December 2009
Published: 16 December 2009
Genes with sex-biased expression often show rapid molecular evolution between species. Previous population genetic and comparative genomic studies of Drosophila melanogaster and D. simulans revealed that male-biased genes have especially high rates of adaptive evolution. To test if this is also the case for other lineages within the melanogaster group, we investigated gene expression in D. ananassae, a species that occurs in structured populations in tropical and subtropical regions. We used custom-made microarrays and published microarray data to characterize the sex-biased expression of 129 D. ananassae genes whose D. melanogaster orthologs had been classified previously as male-biased, female-biased, or unbiased in their expression and had been studied extensively at the population-genetic level. For 43 of these genes we surveyed DNA sequence polymorphism in a natural population of D. ananassae and determined divergence to the sister species D. atripex and D. phaeopleura.
Sex-biased expression is generally conserved between D. melanogaster and D. ananassae, with the majority of genes exhibiting the same bias in the two species. However, about one-third of the genes have either gained or lost sex-biased expression in one of the species and a small proportion of genes (~4%) have changed bias from one sex to the other. The male-biased genes of D. ananassae show evidence of positive selection acting at the protein level. However, the signal of adaptive protein evolution for male-biased genes is not as strong in D. ananassae as it is in D. melanogaster and is limited to genes with conserved male-biased expression in both species. Within D. ananassae, a significant signal of adaptive evolution is also detected for female-biased and unbiased genes.
Our findings extend previous observations of widespread adaptive protein evolution to an independent Drosophila lineage, the D. ananassae subgroup. However, the rate of adaptive evolution is not greater for male-biased genes than for female-biased or unbiased genes, which suggests that there are differences in sex-biased gene evolution between the two lineages.
Sex-biased genes, i.e. those that differ in expression level between males and females, may be subject to differing selective constraints depending on the sex in which they are expressed or they may experience conflicting selective pressures in males and females (reviewed in ). Previous studies of Drosophila melanogaster have shown that male-biased genes, especially those expressed in reproductive tissues, consistently exhibit high levels of adaptive protein evolution [2, 3]. Genome-wide comparisons of the ratio of the nonsynonymous substitution rate to the synonymous substitution rate (dN/dS) also indicate that male-biased genes are more functionally divergent between closely-related Drosophila species and are less likely to have identifiable orthologs between distantly-related species than genes with female-biased or unbiased expression [2, 4]. A limitation of these studies is that they rely on gene expression data from D. melanogaster for sex-bias classifications and, thus, are not informative with respect to differences in sex-biased gene expression or sex-biased gene evolution between lineages.
In contrast to the studies that focused primarily on D. melanogaster, a recent SAGE (serial analysis of gene expression) study found no accelerated rate of protein evolution for male-biased genes in D. pseudoobscura . This study measured the rate of evolution by the proportion of nonsynonymous substitutions (d N ) between species and confirmed a higher rate of protein evolution in genes that had male-biased expression in both D. melanogaster and D. pseudoobscura, but found no evidence for an increased rate of evolution of genes that had male-biased expression only in D. pseudoobscura. The latter genes were only about half as divergent as the former and showed evolutionary rates similar to those of female-biased and unbiased genes. These results suggest that patterns of sex-biased gene evolution may have changed since the split of the D. melanogaster and D. pseudoobscura lineages. To further investigate this possibility, we analyzed sex-biased gene expression and DNA sequence polymorphism in D. ananassae, a species within the melanogaster group that serves as an outgroup to all species within the melanogaster subgroup, but is more closely related to D. melanogaster than D. pseudoobscura [6, 7].
D. ananassae is distributed throughout the subtropical and tropical regions of the world. In contrast to D. melanogaster, D. ananassae is a species displaying significant population structure [8, 9]. The demographic history has been investigated previously, with the ancestral range of the species being defined as a region of Southeast Asia that existed as a single landmass (Sundaland) during the late Pleistocene around 18,000 years ago [8–10]. Ancestral populations are expected to have colonized Asia and the South Pacific since the last glaciation and human migration to Oceania.
In this paper, we use species-specific microarrays to investigate sex-biased gene expression in D. ananassae for a set of genes previously investigated in D. melanogaster. We find that 60% of these genes show conserved sex-bias, while 40% differ in their sex-bias classification between the two species. Using multilocus statistical tests that compare ratios of polymorphism and divergence at synonymous and nonsynonymous sites, we detect a general signal of adaptive protein evolution in D. ananassae. However, this signal is not stronger for male-biased genes than for female-biased or unbiased genes, which is consistent with there being differences in sex-biased gene evolution between the D. melanogaster and D. ananassae lineages.
Sex-biased gene expression in D. ananassae
To investigate sex-biased gene expression in D. ananassae, we designed a species-specific microarray of PCR-amplified exon sequences from 148 genes (Additional file 1) whose orthologs had been previously classified as male-biased, female-biased or unbiased in D. melanogaster (see Methods). This set of genes was of particular interest because the vast majority (136 genes) had been studied at the population-genetic level in D. melanogaster and estimates of the rate of adaptive evolution in the melanogaster subgroup were already available [2, 3]. After quality control, we obtained sufficient hybridization signal to reliably classify 60 genes (22 male-biased, 13 female-biased, and 25 unbiased). Of the remaining genes, we were able to classify 69 (21 male-biased, 10 female-biased, and 38 unbiased) using the D. ananassae whole-genome microarray data of Zhang et al. . In total, 129 of the 148 genes (87%) were classified, with 43 male-biased, 23 female-biased, and 63 unbiased (Additional file 2).
Phylogenetic relationship of the focal species
Levels of polymorphism and divergence
Evolutionary rates of genes with male-, female-, and unbiased expression
d N /d S
Dmel vs. Dsim
Dana vs. Dat
Dana vs. Dph
Comparison of the divergence between D. ananassae and D. atripex/D. phaeopleura with the divergence between D. melanogaster and D. simulans for the different sex-bias classes of genes suggests that the accelerated rate of evolution of male-biased genes is limited to the melanogaster subgroup (Table 1). When comparing D. melanogaster and D. simulans, values of d N and dN/dS are consistently higher for male-biased genes than for female-biased or unbiased genes, which is not the case for comparisons of D. ananassae to either D. atripex or D. phaeopleura (Table 1).
However, there is no significant difference in d N /d S between male-biased genes on the different lineages (Mann-Whitney test, P > 0.05). It should be noted, however, that the number and identity of the genes are not constant among the three comparisons.
Intraspecific polymorphism in D. melanogaster and D. ananassae
π s a
π n b
π n /π s
θ s c
θ n d
D. melanogaster (Zimbabwe, Africa)
D. ananassae (Bangkok, Thailand)
Synonymous polymorphism and divergence at X-linked and autosomal loci
θ s /d S
Genes with significant McDonald-Kreitman tests for positive selection
D S b
P S c
D N d
D N e
D S b
P S c
D N d
P N e
Average values of Tajima's D
A comparison of our MK test results to those from D. melanogaster/D. simulans [2, 3] revealed that, in all but one case (CG10750), the individual genes showing significant evidence for positive selection (i.e. those with MK-test P-values less than 0.05) differed between the two lineages (Table 4). Most genes significant for positive selection were male-biased in at least one of the species. However, the proportion of positively-selected male-biased genes in D. ananassae (1/17 = 6%) was not as high as in previous studies of D. melanogaster where 7/33 (21%) of autosomal male-biased , and 7/17 (41%) of X-linked male-biased genes  gave significant MK tests for positive selection. Of the three genes with significant MK tests for positive selection in D. ananassae, one (CG6980) showed conserved male-biased expression between D. ananassae and D. melanogaster, while two (CG14717 and CG10750) were male-biased in D. melanogaster, but unbiased in D. ananassae (Table 4).
Summary of polymorphic and divergent sites
P S b
P N b
α a, b
Estimation of the proportion of adaptive amino acid replacements
We also used a multi-locus version of the MK test to estimate α, the fraction of amino acid replacements between species that were fixed by positive selection . For this, we divided the genes into six different groups. First, we considered male-, female-, and unbiased genes that had conserved sex-biased expression between D. melanogaster and D. ananassae. Second, we considered the genes according only to their bias in D. ananassae. We could not consider all 12 groups of genes mentioned above, because some of the groups contained too few genes to allow maximum likelihood analysis (e.g. there were only two genes with male-biased expression private to D. ananassae). Because the segregation of slightly deleterious nonsynonymous mutations can lead to a downward bias in the estimate of α, we excluded all polymorphisms (both synonymous and nonsynonymous) segregating at frequency ≤ 15% from our analysis .
Although sex-biased expression was conserved between D. ananassae and D. melanogaster for the majority of the genes we analyzed, a large minority (40%) showed a gain, loss, or reversal of sex-biased expression between species (Figure 1). In a previous microarray study, Ranz et al.  found that 20% of genes showed a gain, loss, or reversal of sex-biased expression between D. melanogaster and D. simulans. The higher percentage observed in the present study is likely attributable to the greater evolutionary distance separating D. melanogaster and D. ananassae (Figure 2), which provides more opportunity for expression changes. Consistent with this, a previous SAGE study found that 34% of genes changed their sex-biased expression pattern between D. melanogaster and the more distantly-related D. pseudoobscura .
Several aspects of our experimental design might also contribute to the gene expression patterns described above. First, the genes on our array were not a random set, but instead were enriched for those showing strong sex-biased expression in multiple, independent experiments in D. melanogaster. Thus, we might expect to see an overrepresentation of genes that are sex-biased only in D. melanogaster, and an underrepresentation of genes that are sex-biased only in D. ananassae. Indeed, this is what we find: of the 48 genes that are sex-biased in only one species, 33 (69%) are sex-biased in D. melanogaster, while only 15 (31%) are sex-biased in D. ananassae (Figure 1B). A second factor that may influence our results is that our D. ananassae classifications come from a combination of our own experimental data and those from a published whole-genome microarray study  and the two studies differ in their design and replication. Our own microarray experiments examined only a small set of genes, which allowed for many replicate probes per gene to be present on each array. We also enforced strict quality control to exclude genes with weak hybridization signal (which would otherwise be classified as unbiased). The use of the published whole-genome arrays allowed us to classify more genes, but the classification may be less reliable due to the lower replication and the increased problem of multiple testing that weakens the statistical analysis. Indeed, the whole-genome data classify a higher fraction of unbiased genes (68% vs. 33% for our custom arrays), which would be expected since the null hypothesis of these comparisons is that there is no difference in expression level between males and females. Also, it should be noted that the whole-genome arrays of Zhang et al.  only identified ~12% of the D. ananassae transcriptome as sex-biased, while the previous studies that our D. melanogaster classifications are based on identified 20-70% of the transcriptome as sex-biased [19–22]. This difference is most likely the result of differences in experimental design and statistical power between the studies, not an underlying difference in the amount of sex-biased expression between the two species, as Zhang et al.  also reported similarly low percentages of sex-biased genes in six other Drosophila species, including D. melanogaster.
Despite the presumably weaker power to detect sex-biased expression in D. ananassae, we observed 15 genes that were sex-biased in D. ananassae, but unbiased in D. melanogaster (Figure 1B). For 12 of these genes, the ancestral expression state could be inferred using microarray data from the outgroup species D. pseudoobscura  (Additional file 4). Interestingly, 11 of these 12 genes showed a match between the D. melanogaster and the D. pseudoobscura classification, suggesting that sex-biased expression was gained on the D. ananassae lineage in the vast majority of cases. The ancestral expression state could also be inferred for four of the five genes that showed a reversal of sex-bias between D. melanogaster and D. ananassae. Three of these genes were female-biased in D. melanogaster and D. pseudoobscura, but male-biased in D. ananassae. One gene was male-biased in D. pseudoobscura and D. melanogaster, but female-biased in D. ananassae. However, this gene (CG7387) also differed in its sex-bias classification between our custom microarrays (female-biased) and the whole-genome D. ananassae microarrays (male-biased) . Thus, it is possible that the sex-biased expression of this gene is strain- or condition-dependent. This was the only such conflict between the two D. ananassae expression datasets, although there were 26 cases where a gene was classified as sex-biased in one dataset and unbiased in the other.
Our survey of DNA sequence polymorphism in D. ananassae is the largest performed to date in terms of number of loci investigated and the first to examine genes with sex-biased expression. Overall, the level of polymorphism in the D. ananassae population from Bangkok, Thailand is similar to that in an ancestral African D. melanogaster population, which is consistent with Bangkok being within the ancestral species range of D. ananassae . When considering all loci we find that D. ananassae has slightly more synonymous polymorphism than D. melanogaster (Table 2), which suggest that the former has a larger Ne. Consistent with this interpretation, nonsynonymous polymorphism and the ratio of nonsynonymous to synonymous polymorphism are lower in D. ananassae (Table 2). This is expected if most segregating nonsynonymous mutations are slightly deleterious, as is suggested by their negative values of Tajima's D (Table 4), because purifying selection is more effective at removing deleterious mutations when Ne is large.
Multi-locus analyses of polymorphism and divergence indicate that adaptive protein evolution is prevalent in D. ananassae, with estimates of α in the range of 50-60% (Table 3 and Figure 3). These α values are remarkably similar to those estimated for other Drosophila species (reviewed in ), which suggests that there is a consistently high rate of adaptive protein evolution throughout the genus. However, despite the strong overall signal of positive selection acting on D. ananassae proteins, we do not see clear differences in the prevalence of adaptive protein evolution among male-, female-, and unbiased genes. This contrasts with previous results from D. melanogaster that indicated an increased rate of adaptive evolution in male-biased genes [2, 3], but is consistent with the results of Metta et al. , who found accelerated rates of evolution (as measured by d N ) for D. melanogaster male-biased genes, but not D. pseudoobscura male-biased genes. We see the same general trend for d N (and d N /d S ) values in our data, which are highest for the male-biased genes of D. melanogaster, but lower for the male-biased genes of D. ananassae (Table 1). Taken together, these results suggest that male-biased genes have experienced an increase in the rate of adaptive protein evolution since the divergence of the ananassae and melanogaster subgroups (Figure 2).
The above results could be explained if sexual selection drives the evolution of many male-biased genes and the prevalence of sexual selection differs between the two subgroups. In ancestral D. melanogaster populations, X-linked synonymous polymorphism is significantly greater than 3/4 of autosomal synonymous polymorphism , which is expected if sexual selection acts on males (reviewed in ). We observe a similar pattern for D. ananassae (Table 3), suggesting that sexual selection does not differ greatly between the species. However, due to the relatively small number of D. ananassae loci, we cannot detect a significant difference between the observed X:autosome diversity ratio and the expected ratio of 3/4, nor can we detect a significant difference in this ratio between D. ananassae and D. melanogaster. Thus, we cannot exclude the possibility that sexual selection is weaker in the D. ananassae lineage. Mate-choice experiments, however, indicate that there is significant female mate preference in D. ananassae  and that male courtship song plays a major role in female mate discrimination . Furthermore, it has been shown that D. ananassae males are subject to stronger intrasexual selection than females and that male mating success is correlated with morphological traits, such as body size and sternopleural bristle number . These findings suggest that sexual selection plays an important role in the species' evolution.
The observed differences in male-biased gene evolution between D. ananassae and D. melanogaster could also be influenced by the particular genes that were investigated. Our initial gene set was enriched for genes that showed strong and consistent sex-biased expression across multiple D. melanogaster microarray experiments. On average, the subset of male-biased genes showed a 12-fold male-bias in D. melanogaster, but only a 5-fold male bias in D. ananassae. While it is difficult to directly compare male/female expression ratios across different experiments and microarray platforms, it is likely that most of the genes investigated had a stronger male-bias in D. melanogaster than in D. ananassae. If the degree of male bias is correlated with the rate of adaptive evolution across the genus, as it appears to be in the melanogaster subgroup , it could explain differences in male-biased gene evolution between the lineages. This possibility could be addressed in future studies focusing on genes with exceptionally strong male-biased expression in D. ananassae. However, there is not a significant correlation between d N /d S and the male/female expression ratio within our D. ananassae dataset (Pearson's R = 0.10, P = 0.52). In general, the female-biased genes included in our study showed weaker sex-biased expression than the male-biased genes in D. melanogaster (4-fold vs. 12-fold). This might explain why expression conservation between D. melanogaster and D. ananassae was greater for male-biased genes than female-biased genes (Figure 1), as highly biased genes are more likely to de detected as significant, while weakly biased genes are more likely to be non-significant and classified as unbiased.
Although sex-biased gene expression is abundant in Drosophila species, the sex-biased expression pattern of many genes differs between species. Species-specific microarray data indicate that 26% of genes with a strong sex bias in D. melanogaster show no detectable sex bias in D. ananassae, while 12% of genes with unbiased expression in D. melanogaster show significantly sex-biased expression in D. ananassae. The accelerated rate of adaptive evolution seen for male-biased genes in D. melanogaster is not observed for male-biased genes in D. ananassae, which suggests that there are differences in sex-biased gene evolution between the two lineages. This is in agreement with a previous study on rates of protein evolution in D. pseudoobscura  and suggests that the rapid adaptive evolution of male-biased genes is unique to the melanogaster subgroup and not a general pattern in Drosophilids. Despite these differences, the overall signal of adaptive protein evolution is strong in D. ananassae (α ≈ 50%) and is consistent with previous estimates throughout the genus.
To analyze sex-biased gene expression in D. ananassae, we designed a species-specific PCR-amplicon microarray. We began with a set of 148 genes that had previously been studied in D. melanogaster [2, 3] and used the available D. ananassae genome sequence (Assembly August 2005; http://genome.ucsc.edu/) to identify their orthologs and design PCR primers that amplify exonic sequences (mean length = 458 bp). A complete list of genes and PCR primers is provided in (Additional file 1). PCR products were purified using genPURE PCR 96 well filter plates (Genetix) and spotted on UltraGAPS slides (Corning) using a GeneMachines OmniGrid Accent microarrayer. D. ananassae genomic DNA (gDNA) was also spotted as a control. Eight replicates of each gene probe were spotted per array, with each replicate in a different subarray. 12-14 control gDNA probes were spotted per subarray.
For hybridization, we extracted total RNA from four- to five-day old males and females using Trizol reagent (Invitrogen) and the manufacturer's protocol. Two inbred strains of D. ananassae from Kota Kinabalu, Borneo were used for RNA extraction . Reverse transcription was performed using 25 μg of total RNA and an anchored oligo(dT) primer. cDNA was labelled with fluorescent dyes (Alexa Fluor 555 and 647) using the amino-allyl labelling system (Invitrogen). Labelled male and female cDNA was competitively hybridized to the arrays for 20 hours at 42 C. Arrays were scanned using a Genetix aQuire 2-laser microarray scanner and the Genetix Qscan software. In total, we performed 12 replicate hybridizations. Six of these were biological replicates, i.e. from different RNA extractions performed at different times, while the other six were "dye-swap" replicates where we exchanged the fluorescent dyes used to label male and female cDNA within each biological replicate.
Prior to statistical analysis, normalization of the two dye channels was performed. Since our arrays contained only a subset of the genome with a non-random distribution of sex-biased genes, we used the gDNA spots as controls for normalization. The signal for each spot was calculated by subtracting the median local background from the mean spot intensity. If this result was negative for a given channel, that channel was assigned a value of 0.5 according to the MINIMUM approach of LIMMA . For each replicate subarray, the red/green ratio was calculated for each of the gDNA control spots and the mean of these ratios was taken as the normalization factor. The raw red value of each spot was then multiplied by this factor to obtain the normalized red value used to calculate the corrected red/green ratio. Note that this is a local normalization, with a different normalization factor for each replicate subarray. As a quality control measure, only spots with a mean intensity at least 20% above the mean local background in at least one of the channels were used. For any one array, we required that at least half of the replicate spots per gene (in this case four) displayed adequate signal as defined above. The median red/green ratio (after normalization) of these replicate spots was used for statistical analysis using the BAGEL software [31, 32]. To determine the P-value cut-off used to define significant sex-biased genes, we performed randomizations of the BAGEL input file to estimate the false-discovery rate (FDR) for the dataset. For classification of sex-biased genes, we used a P-value cut-off of 0.01, which corresponds to a FDR of 10%. All new microarray data have been deposited in NCBI's Gene Expression Omnibus  and are accessible through GEO Series accession number GSE19096 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE19096.
Fly strains, polymerase chain reaction and DNA sequencing
For the polymorphism survey, we used 12 inbred strains of D. ananassae from Bangkok, Thailand. Like the Kota Kinabalu strains used for the microarray experiments, these strains were included in a previous population genetic survey of intronic loci and found to be in the ancestral range of D. ananassae . Thus, we expect little differentiation between these populations. We chose the Bangkok population for the polymorphism survey because a larger number of strains were available.
One strain each of D. atripex and D. phaeopleura (kindly provided by M. Schug) were used as outgroups. We used the D. ananassae genome (Assembly August 2005; http://genome.ucsc.edu/) to design PCR primers flanking the coding sequence of 43 target genes. These genes were chosen to represent different categories of sex-biased gene expression on the basis of the microarray data described above. Following PCR, the amplified products were purified with ExoSAP-IT (USB, Cleveland, OH) and sequenced on both strands using BigDye version 1.1 chemistry and a 3730 automated sequencer (Applied Biosystems, Foster City, CA). The PCR primers were also used as sequencing primers in the polymorphism survey. A complete list of primers is provided in Additional file 5. For some genes, we were unable to get successful PCR or DNA sequence from all 12 D. ananassae strains (Additional file 3). The average number of strains sequenced per gene was 11 and we required at least eight sequences for a gene to be included in our analysis. All new DNA sequences have been submitted to the GenBank/EMBL databases under accession numbers FN546265-FN546780.
Analysis of polymorphism and divergence
Sequences were edited using DNAstar (Madison, WI) and multiple alignments were calculated using MUSCLE http://www.ebi.ac.uk/Tools/muscle/. Polymorphism and divergence statistics were calculated using DnaSP 4.5 . For McDonald-Kreitman (MK) table data we used the number of segregating mutations instead of the number of segregating sites as some genes had sites with three segregating variants. For divergence, we considered only sites with fixed differences between all D. ananassae lines and a single D. atripex (or D. phaeopleura) sequence. Multilocus Tajima's D tests were performed using the HKA program http://lifesci.rutgers.edu/~heylab/heylabsoftware.htm. To calculate the fraction of positively selected amino acid substitutions, a, with the method of , we used the DoFE program http://www.lifesci.sussex.ac.uk/home/Adam_Eyre-Walker/Website/Software.html. For MK tests and a calculations, the D. atripex sequence was used as the outgroup whenever available. When the D. atripex sequence was not available, D. phaeopleura was used as the outgroup.
For 13 genes we were able to get sequences from both D. atripex and D. phaeopleura [Additional file 3]. In addition, we downloaded the amino acid sequences of D. melanogaster, D. simulans, D. ananassae, and D. pseudoobscura for the entire set of genes analyzed from Flybase http://flybase.org/. For one gene (CG18418), we were unable to get the orthologous sequence from D. pseudoobscura. This left us with 12 genes that could be aligned across all species (CG2577, CG3004, CG3024, CG4593, CG6459, CG7840, CG8277, CG9135, CG10853, CG11379, CG15336, and CG15717). Amino acid sequences of these 12 genes were concatenated for phylogenetic reconstruction Bayesian inference (BI), as implemented in MrBayes version 3.1.1 , was used to reconstruct the phylogeny with D. pseudoobscura set as the outgroup. For the BI analysis, two distinct 10,000 generation runs were conducted (three incrementally heated chains with model jumping between fixed-rate amino acid models were used and trees were saved to a file every 10 generations). Identical topologies were recovered from both runs. A burn-in period of 2,500 generations was determined graphically by plotting likelihood values for each sample. The results were presented in the form of a 50% majority-rule consensus tree in which trees corresponding to the burn-in period were discarded. Support for the nodes was given by posterior probability estimates of clades. Using model jumping, only the Jones (JTT) model  contributed to the final result with 100% posterior probability. This was also supported by the program ProtTest version 2.0 , which determines the best-fit substitution model for amino acid data under a likelihood framework. Within this framework, the Akaike information criterion  selected the JTT  model including the proportion of invariable sites and the gamma distribution of rate variation among sites (JTT + I + G) as the best-fitting model with an Akaike weight of 0.54. The second-best model was JTT+G with an Akaike weight of 0.46. Trees were visualized with TreeView .
All additional material is provided on a small webpage that can be downloaded as a zip-archive. All files can be accessed through index.html, which can be viewed with any standard internet browser (Additional File 6)
We thank A. Das and W. Stephan for providing D. ananassae strains and M. Schug for providing D. atripex and D. phaeoplura strains. This research was supported by Deutsche Forschungsgemeinschaft grant PA 903/4-1.
- Ellegren H, Parsch J: The evolution of sex-biased genes and sex-biased gene expression. Nat Rev Genet. 2007, 8: 689-698. 10.1038/nrg2167.View ArticlePubMedGoogle Scholar
- Baines JF, Sawyer SA, Hartl DL, Parsch J: Effects of X-linkage and sex-biased gene expression on the rate of adaptive protein evolution in Drosophila. Mol Biol Evol. 2008, 25: 1639-1650. 10.1093/molbev/msn111.PubMed CentralView ArticlePubMedGoogle Scholar
- Pröschel M, Zhang Z, Parsch J: Widespread adaptive evolution of Drosophila genes with sex-biased expression. Genetics. 2006, 174: 893-900. 10.1534/genetics.106.058008.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Z, Parsch J: Positive correlation between evolutionary rate and recombination rate in Drosophila genes with male-biased expression. Drosophila. 2005, 22: 1945-1947.Google Scholar
- Metta M, Gudavalli R, Gibert JM, Schlötterer C: No accelerated rate of protein evolution in male-biased Drosophila pseudoobscura genes. Genetics. 2006, 174: 411-420. 10.1534/genetics.106.057414.PubMed CentralView ArticlePubMedGoogle Scholar
- Drosophila 12 Genomes Consortium, et al: Evolution of genes and genomes on the Drosophila phylogeny. Nature. 2007, 450: 203-218. 10.1038/nature06341.View ArticleGoogle Scholar
- Larracuente AM, Sackton TB, Greenberg AJ, Wong A, Singh ND, Sturgill D, Zhang Y, Oliver B, Clark AG: Evolution of protein-coding genes in Drosophila. Trends in Genetics. 2008, 24-3: 114-123. 10.1016/j.tig.2007.12.001.View ArticleGoogle Scholar
- Das A, Mohanty S, Stephan W: Inferring population structure and demography of Drosophila ananassae from multilocus data. Genetics. 2004, 168: 1975-1985. 10.1534/genetics.104.031567.PubMed CentralView ArticlePubMedGoogle Scholar
- Schug MD, Smith SG, Tozier-Pearce A, McEvey SF: The genetic structure of Drosophila ananassae populations from Asia, Australia and Samoa. Genetics. 2007, 175: 1429-1440. 10.1534/genetics.106.066613.PubMed CentralView ArticlePubMedGoogle Scholar
- Das A: Population genomics and bioinformatic studies reveal evolutionary history of Drosophila ananassae. Current Science. 2005, 89: 8-Google Scholar
- Zhang Y, Sturgill D, Parisi M, Kumar S, Oliver B: Constraint and turnover in sex-biased gene expression in the genus Drosophila. Nature. 2007, 450: 233-238. 10.1038/nature06323.PubMed CentralView ArticlePubMedGoogle Scholar
- Da Lage JL, Kergoat GJ, Maczkowiak F, Silvain JF, Cariou ML, Lachaise D: A phylogeny of Drosophilidae using the Amyrel gene: questioning the Drosophila melanogaster species group boundaries. J Zool Syst Evol Res. 2007, 45: 47-63. 10.1111/j.1439-0469.2006.00389.x.View ArticleGoogle Scholar
- Matsuda M, Ng C-S, Doi M, Kopp A, Tobari YN: Evolution in the Drosophila ananassae species subgroup. Fly. 2009, 3 (2): 157-169.View ArticlePubMedGoogle Scholar
- McDonald JH, Kreitman M: Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991, 351: 652-654. 10.1038/351652a0.View ArticlePubMedGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMed CentralPubMedGoogle Scholar
- Fay JC, Wyckoff GJ, Wu CI: Positive and negative selection on the human genome. Genetics. 2001, 158: 1227-1234.PubMed CentralPubMedGoogle Scholar
- Bierne N, Eyre-Walker A: The genomic rate of adaptive amino acid substitution in Drosophila. Mol Biol Evol. 2004, 21: 1350-1360. 10.1093/molbev/msh134.View ArticlePubMedGoogle Scholar
- Charlesworth J, Eyre-Walker A: The McDonald-Kreitman test and slightly deleterious mutations. Mol Biol Evol. 2008, 25: 1007-1015. 10.1093/molbev/msn005.View ArticlePubMedGoogle Scholar
- Ranz JM, Castillo-Davis CI, Meiklejohn CD, Hartl DL: Sex-dependent gene expression and evolution of the Drosophila transcriptome. Science. 2003, 300: 1742-1745. 10.1126/science.1085881.View ArticlePubMedGoogle Scholar
- Gibson G, Riley-Berger R, Harshman L, Kopp A, Vacha S, Nuzhdin S, Wayne M: Extensive sex-specific non-additivity of gene expression in Drosophila melanogaster. Genetics. 2004, 167: 1791-1799. 10.1534/genetics.104.026583.PubMed CentralView ArticlePubMedGoogle Scholar
- Parisi M, Nuttall R, Naiman D, Bouffard G, Malley J, Andrews J, Eastman S, Oliver B: Paucity of genes on the Drosophila X chromosome showing male-biased expression. Science. 2003, 299: 697-700. 10.1126/science.1079190.PubMed CentralView ArticlePubMedGoogle Scholar
- Gnad F, Parsch J: Sebida: a database for the functional and evolutionary analysis of genes with sex-biased expression. Bioinformatics. 2006, 22: 2577-2579. 10.1093/bioinformatics/btl422.View ArticlePubMedGoogle Scholar
- Sella G, Petrov DA, Przeworski M, Andolfatto P: Pervasive Natural Selection in the Drosophila Genome?. PLoS Genetics. 2009, 5: e1000495-10.1371/journal.pgen.1000495.PubMed CentralView ArticlePubMedGoogle Scholar
- Parsch J, Zhang Z, Baines JF: The influence of demography and weak selection on the McDonald-Kreitman test: an empirical study in Drosophila. Mol Biol Evol. 2009, 26: 691-698. 10.1093/molbev/msn297.View ArticlePubMedGoogle Scholar
- Ellegren H: The different levels of genetic diversity in sex chromosomes and autosomes. Trends in Genetics. 2009, 25: 278-284. 10.1016/j.tig.2009.04.005.View ArticlePubMedGoogle Scholar
- Schug MD, Baines JF, Killon-Atwood A, Mohanty S, Das A, Grath S, Smith SG, Zargham S, McEvey SF, Stephan W: Evolution of mating isolation between populations of D. ananassae. Molecular Ecology. 2008, 17: 2706-2721. 10.1111/j.1365-294X.2008.03770.x.View ArticlePubMedGoogle Scholar
- Doi M, Matsuda M, Tomaru M, Matsubayashi H, Oguma Y: A locus for female discrimination behavior causing sexual isolation in Drosophila. Proc Natl Acad Sci. 2001, 98: 6714-6719. 10.1073/pnas.091421598.PubMed CentralView ArticlePubMedGoogle Scholar
- Singh SR, Singh BN: Behavioral genetics of Drosophila ananassae. Genetics and Molecular Research. 2003, 2: 394-409.PubMedGoogle Scholar
- FlyBase: A Database of Drosophila Genes & Genomes. [http://flybase.org/]
- Smyth GK, Speed TP: Normalization of cDNA microarray data. Methods. 2002, 31: 265-273. 10.1016/S1046-2023(03)00155-5.View ArticleGoogle Scholar
- Townsend JP, Hartl DL: Bayesian analysis of gene expression levels: statistical quantification of relative mRNA level across multiple strains or treatments. Genome Biol. 2002, 3: 12-10.1186/gb-2002-3-12-research0071.View ArticleGoogle Scholar
- Townsend JP: Resolution of large and small differences in gene expression using models for the Bayesian analysis of gene expression levels and spotted DNA microarrays. BMC Bioinformatics. 2004, 5: 54-10.1186/1471-2105-5-54.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMed CentralView ArticlePubMedGoogle Scholar
- UCSC Genome Browser. [http://genome.ucsc.edu]
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5: 113-10.1186/1471-2105-5-113.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MrBayes: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8: 275-282.PubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: Selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.View ArticlePubMedGoogle Scholar
- Akaike H: Information theory and an extension of the Maximum Likelihood principle. Second International Symposium on Information Theory. Edited by: Petrov BN, Csaki F. 1973, Budapest: Akademiai Kiado, 267-281.Google Scholar
- Page RDM: TREEVIEW: An application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.