PMCCPMCCPMCC

Search tips
Search criteria 

Advanced


 
Formats
Article sections
Authors
Related links
Logo of plosgenPLoS GeneticsSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Genet. 2012 October; 8(10): e1003013.
Published online 2012 October 11. doi:  10.1371/journal.pgen.1003013
PMCID: PMC3469423
Copyright notice
Faster-X Evolution of Gene Expression in Drosophila
Richard P. Meisel,1* John H. Malone,2 and Andrew G. Clark1
1Department of Molecular Biology and Genetics, Cornell University, Ithaca, New York, United States of America
2Department of Biological Science, Florida State University, Tallahassee, Florida, United States of America
Doris Bachtrog, Editor
University of California Berkeley, United States of America
* E-mail: meisel/at/cornell.edu
The authors have declared that no competing interests exist.
Conceived and designed the experiments: RPM JHM AGC. Performed the experiments: RPM JHM. Analyzed the data: RPM JHM. Wrote the paper: RPM JHM AGC.
Received March 28, 2012; Accepted August 22, 2012.
This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Small right arrow pointing to: See "An Excess of Gene Expression Divergence on the X Chromosome in Drosophila Embryos: Implications for the Faster-X Hypothesis" , e1003200.
 
Small right arrow pointing to: This article has been cited by other articles in PMC.
Abstract
DNA sequences on X chromosomes often have a faster rate of evolution when compared to similar loci on the autosomes, and well articulated models provide reasons why the X-linked mode of inheritance may be responsible for the faster evolution of X-linked genes. We analyzed microarray and RNA–seq data collected from females and males of six Drosophila species and found that the expression levels of X-linked genes also diverge faster than autosomal gene expression, similar to the “faster-X” effect often observed in DNA sequence evolution. Faster-X evolution of gene expression was recently described in mammals, but it was limited to the evolutionary lineages shortly following the creation of the therian X chromosome. In contrast, we detect a faster-X effect along both deep lineages and those on the tips of the Drosophila phylogeny. In Drosophila males, the dosage compensation complex (DCC) binds the X chromosome, creating a unique chromatin environment that promotes the hyper-expression of X-linked genes. We find that DCC binding, chromatin environment, and breadth of expression are all predictive of the rate of gene expression evolution. In addition, estimates of the intraspecific genetic polymorphism underlying gene expression variation suggest that X-linked expression levels are not under relaxed selective constraints. We therefore hypothesize that the faster-X evolution of gene expression is the result of the adaptive fixation of beneficial mutations at X-linked loci that change expression level in cis. This adaptive faster-X evolution of gene expression is limited to genes that are narrowly expressed in a single tissue, suggesting that relaxed pleiotropic constraints permit a faster response to selection. Finally, we present a conceptional framework to explain faster-X expression evolution, and we use this framework to examine differences in the faster-X effect between Drosophila and mammals.
Author Summary
As species diverge over evolutionary time, they accumulate differences in the sequences of their genes and how those genes are expressed. We show that gene expression changes accumulate faster for genes on the X chromosome than for genes on the other chromosomes (autosomes) in Drosophila (the “faster-X” effect). The X chromosome is only found in a single copy in males, whereas the autosomes are found in two copies in both sexes. To compensate for the reduced dosage of X-linked genes in males, a molecular complex binds the Drosophila X chromosome to upregulate gene expression in males. We demonstrate that genes that escape this dosage compensation process have faster evolving expression levels. X-linked genes are inherited in a unique manner, and we hypothesize that this permits a faster rate of adaptive evolution, thereby driving the faster-X evolution of gene expression. We compare these observations with the recently described faster-X evolution of gene expression in mammals, and we explain how differences in dosage compensation, mutation rate, and population size could affect the extent of the faster-X effect.
Comparing the evolutionary rates of X-linked (or Z-linked) and autosomal genes can be informative of the nature of allelic dominance [1], the type of variation acted upon by natural selection [2], [3], the mutational process [4]–[8], and the effect of differences in population size on the efficacy of natural selection across taxa [9], [10]. Notably, DNA sequences on X (or Z) chromosomes often evolve faster than autosomal sequences (i.e., the “faster-X” effect). This may be a result of the adaptive fixation of recessive beneficial mutations in X-linked genes [1], [11]–[14], mutational biases associated with dosage compensation [15], or the smaller effective population size (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e001.jpg) of sex chromosomes [9], [10]. The faster-X effect is especially pronounced in the protein coding sequences of genes with male-biased expression (i.e., genes expressed higher in males than females) or genes specifically expressed in male reproductive tissues in male heterogametic (XY) taxa [16]–[20]. These results support the theoretical prediction that the adaptive fixation of recessive X-linked male-beneficial mutations in hemizygous males can drive faster-X evolution [1].
Comparisons of expression divergence between X-linked and autosomal genes are not as prevalent as analyses of DNA sequences. Some experiments have suggested that the disproportionate effect of X-linked loci on interspecific hybrid fitness (the “large X” effect [21]) is the result of divergence in the regulation of gene expression. For example, gene expression from the X chromosome may be misregulated in the male germline of interspecific hybrids [22]–[24], and dosage compensation of the X chromosome could also be affected in hybrids [25]–[27]. With the advent of high throughput technologies to measure expression in multiple species we can now directly test whether the rate of expression evolution differs between X-linked and autosomal genes. The first such analysis did indeed find evidence for the faster-X evolution of gene expression shortly following the creation of the therian X chromosome [28].
Gene expression is determined by an interaction of cis regulatory elements and the proteins that bind to them (e.g., transcription factors, histones) to either promote or inhibit transcription. X chromosomes often have a unique chromatin environment because of the need to dosage compensate X-linked genes in males. In mammals, this is hypothesized to be accomplished by the upregulation of X-linked gene expression in both sexes, followed by random silencing of one X chromosome in females [29]–[31] (although this model is not universally accepted [32]). Drosophila compensate for reduced X chromosome dose in males by modifying the chromatin structure of the X in a male-specific manner. The dosage compensation complex (DCC; or male-specific lethal [MSL] complex), a ribonucleoprotein structure, binds the X chromosome in males, acetylating histone H4 at lysine 16 [33]–[35]. This is thought to promote the expression of X-linked genes via some combination of relaxing compacted chromatin [36], [37], enhancing recruitment of RNA polymerase II [38], and/or increasing transcriptional elongation [39]. The DCC only assembles in males because one of the essential proteins, MSL-2, is not produced in females [40]–[42].
Recently, chromatin immunoprecipitation (ChIP) experiments followed by microarrays (ChIP-chip) or sequencing (ChIP-seq) have revealed regions of the Drosophila melanogaster X chromosome that are enriched with DCC binding and bound by the DCC in the absence of essential DCC components [43], 44. These chromatin entry or high affinity sites (HASs) contain a DNA sequence motif that is thought to direct the DCC to the Drosophila X chromosome [43], [44]. After initially binding to the 100–300 HASs, the DCC is hypothesized to spread in cis to promote the upregulation of expression by inducing transcriptionally activating chromatin marks [45]–[49].
To examine how X-linkage, chromatin environment, and breadth of expression affect the evolution of gene expression, we calculated expression differences between Drosophila species using data collected from male and female whole flies and heads using microarrays and high throughput RNA sequencing (RNA-seq) [50]–[52]. We detect a robust signal of faster-X evolution of gene expression. This faster-X effect is most pronounced in genes that are located in transcriptionally repressive chromatin in cell culture and genes that are narrowly expressed in a limited number of tissues. In addition, we analyzed measurements of intraspecific variation in gene expression, and we found that the faster-X effect cannot be explained by relaxed selective constraints. Our results suggest that the faster-X evolution of gene expression is the result of the adaptive fixation of X-linked mutations that affect gene expression in cis.
Faster-X evolution of gene expression
We analyzed expression measurements in six Drosophila species (D. melanogaster, Drosophila yakuba, Drosophila ananassae, Drosophila pseudoobscura, Drosophila mojavensis, and Drosophila virilis) collected from whole females and males separately using microarrays. Following the method of Brawand et al. [28], we calculated the similarity in expression between pairs of species, within each sex, using Spearman's rank correlation coefficient (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e002.jpg) sampling only genes present as 1[ratio]1[ratio]1[ratio]1[ratio]1[ratio]1 orthologs across all six species [53]. To determine if the expression levels of X-linked genes diverge faster than autosomal genes, we compared An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e003.jpg across the five major chromosome arms (also known as Muller elements [54]). In every pairwise comparison, the correlation of expression levels of X-linked genes (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e004.jpg), in both females and males, is significantly lower than that of autosomal genes (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e005.jpg) (Figure 1). We confirmed this pattern using a different pipeline to handle the microarray data (Figure S1) and with RNA-seq data from three species (D. melanogaster, D. pseudoobscura, and D. mojavensis; Figure S2). These results suggest that X-linked gene expression levels diverge faster than the expression levels of autosomal genes.
Figure 1
Figure 1
Faster-X evolution of gene expression.
While there is congruence between the microarray and RNA-seq data with regards to the faster-X evolution of expression, we observe two notable differences between these data sets. First, expression levels estimated using RNA-seq are more highly correlated than those estimated from microarray data (Figure 1, Figure S2), possibly because of the increased dynamic range of RNA-seq [55]. Second, the magnitude of the difference between An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e006.jpg and An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e007.jpg is greater in the microarray data than the RNA-seq data (Figure 1, Figure S2). We reanalyzed the microarray data using only the genes present in the RNA-seq dataset, and these correlations resemble the microarray analysis more than the RNA-seq (Figure S3). Therefore, the difference in magnitude of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e008.jpg in the microarray and RNA-seq analyses is not attributable to differences in the gene sets analyzed. Regardless of the cause of these differences, both methodologies provide evidence in support of the faster-X evolution of gene expression (Figure 1, Figure S2).
Faster-X expression evolution is not limited to genes with male-biased expression or new genes
The faster-X evolution of Drosophila protein coding sequences is especially pronounced in genes with male-biased expression that are expressed in male reproductive tissues [18]–[20], possibly because the hemizygosity of the X chromosome favors the adaptive fixation of recessive male-beneficial mutations in X-linked genes [1]. Additionally, genes with male-biased expression tend to have more divergent expression between species than genes with female-biased or non-sex-biased expression [56], [57]. The faster-X evolution of gene expression, however, is not limited to genes expressed in male reproductive tissues: we detect the faster-X effect when gene expression is measured in females (Figure 1) or heads (Figure 2, Figure S4), although the pattern is not as striking as when whole fly data are used.
Figure 2
Figure 2
Faster-X evolution of gene expression measured in head.
To further examine the effect of expression in male-reproductive tissues on the faster-X effect, we excluded genes with male-biased expression in D. melanogaster (765 genes at a false discovery rate [FDR] of 0.05), male-biased expression in all of the species (2027 genes at a FDR of 0.20), or biased expression in male reproductive tissues in D. melanogaster (439 genes based on expression data from FlyAtlas [58]). In all cases, we detect the faster-X evolution of gene expression even when genes with male-biased expression are removed (Figures S5, S6, S7). In addition, genes that are narrowly expressed in non-reproductive tissues exhibit a faster-X effect comparable to genes with biased expression in male reproductive tissues (Figure 3). On the other hand, we fail to detect the faster-X effect when we consider only genes with biased expression in female-limited reproductive tissues (Figure 3). We therefore conclude that the faster-X evolution of gene expression requires expression in males but not necessarily male-biased expression.
Figure 3
Figure 3
Faster-X expression evolution for genes expressed in male, but not female, reproductive tissues.
Genes that arose recently by duplication tend to have male-biased expression [59], [60]. Many new Drosophila genes are located on the X chromosome and show evidence of a faster-X effect in their protein coding sequences [61]. We find that new genes (those that arose after the split between the D. melanogaster and D. virilis lineages) do not exhibit evidence of faster-X expression evolution, while old genes shared by the entire genus do (Figure S8). Therefore, the faster-X evolution of male reproductive genes, genes with male-biased expression, or new genes are not entirely responsible for the faster-X evolution of gene expression in Drosophila.
Faster-X expression evolution along both internal and tip lineages of the Drosophila phylogeny
Previous work in mammals found evidence for faster-X evolution of gene expression that was limited to evolutionary lineages closely following the creation of the therian X chromosome [28]. To test for a similar lineage-specific faster-X effect in Drosophila, we used our calculations of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e009.jpg from whole fly expression measurements between species to estimate branch lengths along the Drosophila phylogeny. This approach assumes that there is a phylogenetic signal in these pairwise correlations. Some pairwise correlations are lower for more closely related species than more distantly related ones (Figure 1), suggesting that the correlations may not have an underlying phylogenetic signal. To further test for phylogenetic signal, we estimated the divergence in expression between species as An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e010.jpg [28]. We were indeed able to reconstruct the evolutionary relationships of the six species using this distance matrix (Figure S9), demonstrating that the expression correlations contain a phylogenetic signal.
We next tested whether the faster-X evolution of gene expression is limited to particular branches in the Drosophila phylogeny using matrices of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e011.jpg and An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e012.jpg to estimate branch lengths along the known topology (Figure 4). These branch lengths represent the contribution that each lineage makes toward An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e013.jpg, and larger values indicate a lower correlation in expression. Phylogenies constructed using X-linked gene expression from both males and females have longer internal and terminal branch lengths (Figure 4), unlike in mammals where a faster-X effect is only detected on internal branches [28].
Figure 4
Figure 4
Branch length estimates and bootstrap support from expression level divergence.
Interestingly, branch length estimates closest to the root do not show evidence for a faster-X effect in Drosophila (Figure 4). This is not necessarily evidence against the faster-X evolution of gene expression along these internal branches. We instead hypothesize that it is the result of low power to resolve deep branching orders using these correlation matrices, which leads to poor estimates of branch lengths around deep nodes. Supporting this hypothesis, when we use the correlation matrices to estimate the tree topology, some of the deep nodes have the lowest bootstrap support for the correct topology (Figure 4, Figure S9). In addition, the bootstrap support for these nodes is lower for X-linked gene expression levels than autosomal expression (Figure 4). Furthermore, when we exclude genes on the X chromosome, we observe an increase in bootstrap support for the correct branching order between D. pseudoobscura and the melanogaster group (Figure 4, bottom number in bootstrap boxes, in italics). We therefore hypothesize that the faster-X evolution of gene expression complicates the inference of the correct branching order more for X-linked genes than autosomal genes at these deep nodes, leading to a flawed measurement of the branch lengths. In summary, depending on the branch in question, either longer branch lengths or lower bootstrap values support the hypothesis that X-linked gene expression levels diverge faster than autosomal expression levels across most of the phylogeny.
D. pseudoobscura and Drosophila willistoni each have an independently derived neo-X chromosome arm (Muller element D) that is autosomal in all other species [62]. If the faster evolution of gene expression closely follows the creation of an X chromosome, we would expect to detect a faster-X effect in genes on these neo-X chromosome arms. Using available RNA-seq data collected from D. pseudoobscura, D. willistoni, D. melanogaster, and D. mojavensis heads [52], we find some evidence for the faster evolution of gene expression on the neo-X chromosome arms (Figure S4). However, we fail to detect evidence for faster-X expression evolution in genes on the D. pseudoobscura neo-X chromosome when expression is measured in whole fly (Figure 1, Figure 4). The latter result may be because of low power to detect a faster-neo-X effect; the X-autosome fusion giving rise to the D. pseudoobscura neo-X occurred recently relative to the pseudoobscura-melanogaster common ancestor [60], [63].
Faster-X expression evolution of genes unbound by the DCC and further from HASs
The Drosophila X chromosome has a unique chromatin environment because of the need to compensate dosage in hemizygous males [33], [34], [36], [64], and these histone modifications are correlated with gene expression levels [65]. We therefore considered whether DCC binding is associated with the faster-X effect. To do so, we calculated pairwise expression divergence for each 1[ratio]1[ratio]1 ortholog between D. melanogaster, D. yakuba, and D. ananassae. We selected these three closely related species because DCC binding and chromatin states have only been inferred for D. melanogaster [45], [46], and these inferences are less likely to be accurate in more distantly related species. In addition, these gene-wise estimates of expression divergence differ from the correlations in expression levels across entire chromosomes (see Methods). We chose this approach because, as we introduce more parameters into our analysis, gene-wise expression divergence is easier to interpret than correlations of chromosome-wide expression between species.
Using the gene-wise estimates of expression divergence between D. melanogaster and either D. yakuba or D. ananassae, we found that X-linked genes that are unbound by the DCC have greater expression divergence than DCC bound X-linked genes (Figure 5A). Additionally, in the comparison between D. melanogaster and D. ananassae, unbound X-linked genes have greater expression divergence than autosomal genes (Figure 5A). Furthermore, the expression levels of DCC bound X-linked genes are more evolutionarily conserved than autosomal genes (Figure 5A). DCC bound genes tend to be in close proximity to HASs [48], and HASs have the highest concentration of DCC binding [43], suggesting that proximity to an HAS may also predict expression divergence. Distance from the nearest HAS is indeed positively correlated with expression divergence (Figure 5B). We observe these patterns when expression is measured in either females or males (Figure 5A–5B).
Figure 5
Figure 5
Faster expression evolution of X-linked genes not directly regulated by the DCC.
Highly expressed genes tend to have more conserved protein coding sequences [19], [66], and there may be a positive correlation between the rate of protein coding sequence evolution and divergence in gene expression [67]–[69]. Genes bound by the DCC have higher expression levels than unbound genes [48], suggesting that the relationship between DCC binding and expression divergence (Figure 5A–5B) may be a byproduct of highly expressed genes with less expression divergence. We found a negative correlation between expression level and expression divergence for both X-linked and autosomal genes (Figure 5C), demonstrating that highly expressed genes have more conserved expression levels.
To test whether the relationship between DCC binding and expression divergence is merely an artifact of the correlation between expression level and expression divergence, we calculated partial correlations between expression divergence, distance from the nearest HAS, and expression level. If DCC binding and expression divergence are directly related, genes further from an HAS should have elevated expression divergence even when expression level is taken into account. Distance from the nearest HAS is positively correlated with expression divergence in most of our partial correlations (Figure 5D), demonstrating that genes that are not directly regulated by the DCC have faster evolving expression levels. In addition, expression level and expression divergence are negatively correlated (Figure 5D), supporting the hypothesis that highly expressed X-linked genes have more conserved expression levels independent of DCC binding. Lastly, distance from an HAS is negatively correlated with expression level (Figure 5D), providing additional evidence that highly expressed genes are more directly regulated by the DCC [48].
Faster expression evolution of X-linked genes located in transcriptionally repressive chromatin
The DCC is both attracted to and promotes chromatin modifications associated with transcriptional activity [35], [49], suggesting that genes unbound by the DCC are in transcriptionally repressive chromatin. The faster expression evolution of X-linked genes that are unbound by the DCC could therefore be a general property of genes associated with repressive chromatin. To test this hypothesis, we obtained mapped chromatin states in the D. melanogaster genome from two different cell lines [65], and we used these data to assign genes to one of two chromatin states: transcriptionally active or repressive. We found that X-linked genes that are bound by the DCC are indeed almost always (97.8–100%) associated with active chromatin, while unbound genes tend to be in repressive chromatin (Table S1). In addition, genes in transcriptionally active chromatin have higher expression levels than genes in repressive chromatin (Figure S10).
Both autosomal and X-linked genes associated with transcriptionally repressive chromatin have more divergent expression levels than genes associated with active chromatin regardless of which cell type is used to infer chromatin states (Figure 6). However, X-linked genes that are located in repressive chromatin have more divergent expression between D. melanogaster and D. ananassae than autosomal genes in repressive chromatin (Figure 6). Furthermore, X-linked genes associated with active chromatin tend to have less expression divergence than other genes in the genome (Figure 6). We observe similar patterns when we use DCC binding as a proxy for transcriptionally active chromatin in X-linked genes (Figure S11). These results provide further support for the hypothesis that the faster-X evolution of gene expression is driven by genes that are not directly regulated by the DCC.
Figure 6
Figure 6
Faster expression evolution of genes associated with transcriptionally repressive chromatin.
Expression breadth, chromatin environment, and the faster-X effect
Genes expressed narrowly (i.e., in a limited number of tissues) tend to have rapidly evolving protein coding sequences [19], [66], which raises the possibility that expression breadth may also affect expression divergence. We find that narrowly expressed genes do tend to have elevated expression divergence (Figure 7A, Figure S12). In addition, DCC bound genes tend to be broadly expressed, and genes further from an HAS are more narrowly expressed (Figure 7B, Figure S13). Furthermore, genes that are in transcriptionally active chromatin tend to be broadly expressed, while genes in repressive chromatin tend to be narrowly expressed (Figure 7C, Figure S14). This raises the possibility that the association between chromatin environment and the faster-X effect (Figure 6) could be an artifact of the correlation between expression breadth and expression divergence.
Figure 7
Figure 7
Association between expression breadth, DCC binding, chromatin state, and expression divergence.
If the faster-X evolution of gene expression is affected by expression breadth and not chromatin environment, we expect to only detect the faster-X effect in narrowly expressed genes. Consistent with this prediction, we detect the strongest evidence for faster-X evolution in narrowly expressed genes (Figure 8, Figure S15). The faster-X effect is, however, limited to narrowly expressed genes in transcriptionally repressive chromatin (Figure 8, Figure S15), suggesting that narrow expression and transcriptionally repressive chromatin environment both promote faster-X expression evolution. Narrowly expressed genes in transcriptionally repressive chromatin are more likely to have low expression levels [19], [66] (Figure S10), which could increase the error in expression level measurements. However, measurement error is unlikely to explain the association of expression breadth and chromatin environment with the faster-X effect for two reasons. First, experimental and biological variance should not produce the consistent signal of faster-X evolution. Second, we still detect the faster-X effect when genes with low expression levels are excluded (Figure S16). x
Figure 8
Figure 8
The faster-X effect is limited to narrowly expressed genes in transcriptionally repressive chromatin.
The faster-X effect could be a result of differences in gene content between the X chromosome and the autosomes if, for example, X-linked genes were more narrowly expressed or more likely to be in transcriptionally repressive chromatin than autosomal genes. The D. melanogaster X chromosome, however, harbors a deficiency of narrowly expressed genes [52], [70], and there is a paucity of X-linked genes in repressive chromatin (Figure 6; An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e034.jpg using Fisher's exact test). In addition, we fail to detect a significant difference in expression breadth between X-linked and autosomal genes in repressive chromatin (Figure S17). It is therefore unlikely that the unique gene content of the Drosophila X chromosome is responsible for the faster-X effect. Our power to detect the faster-X effect is also limited by the small sample size of narrowly expressed genes in repressive chromatin on the X chromosome, demonstrating that our results are conservative.
The faster-X evolution of protein-coding sequences is most pronounced for genes that are narrowly expressed in male reproductive tissues [19], [20]. We showed, however, that expression in male reproductive tissues is not solely responsible for the faster-X evolution of gene expression (Figure 2, Figure 3; Figures S4, S5, S6, S7). This does not exclude the possibility that X-linked genes expressed in male reproductive tissues have faster evolving expression levels than autosomal genes (e.g., Figure 3). We do find some support for the faster-X evolution of male expression levels among genes in repressive chromatin that are expressed narrowly in male reproductive tissues, but the evidence is not exceedingly strong (Figure S18). Most notably, we fail to detect the faster-X effect when we limit the analysis to genes in repressive chromatin that are narrowly expressed in female reproductive tissues (Figure S18), consistent with our earlier analysis of chromosome-wide correlations of expression (Figure 3). Therefore, genes with limited expression in females do not experience faster-X expression evolution.
Faster-X expression evolution is not driven by relaxed constraints
Accelerated evolutionary divergence can be the result of relaxed selective constraints or an elevated rate of adaptive evolution. To distinguish between these two explanations for increased divergence in gene expression one can examine intraspecific variation in expression levels [57], [71], [72]. If relaxed selective constraints were responsible for greater divergence, we would expect increased intraspecific variation in genes with rapidly evolving expression levels. Conversely, if the fast evolution of gene expression is driven by positive selection, we expect rapidly evolving genes to have equivalent (or less) expression polymorphism when compared to non-rapidly evolving genes. In making these interpretations we assume that expression variation segregating in natural populations has neutral or slightly deleterious fitness effects, an assumption common to the interpretation of DNA sequence polymorphism and divergence data [73], [74].
One way to estimate intraspecific variation is to compare expression levels between females and males of the same species. Higher correlation between sexes suggests greater constraint on gene expression. We find that X-linked expression levels are often more correlated between the sexes than autosomal expression levels (Figure 9A), suggesting that X-linked expression levels are not under relaxed constraints.
Figure 9
Figure 9
Intraspecific expression variation, X-linkage, DCC binding, expression breadth, and chromatin state.
We also used available calculations of the broad sense heritability (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e044.jpg) of gene expression measured in whole flies from 40 inbred D. melanogaster lines [75] as an estimate of the intraspecific variation in gene expression contributed by genetic factors. Higher An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e045.jpg implies greater genetic variation underlying gene expression levels, which suggests relaxed selective constraints. In the results presented below, all genes with estimates of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e046.jpg were included, but we obtain similar results if we limit ourselves to only genes included in our analysis of expression divergence. Consistent with a previous analysis of the same data [75] and independent experiments in Drosophila simulans [76], we detect significantly reduced An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e047.jpg among X-linked genes (Figure 9B). This provides further evidence that the expression regulation of X-linked genes is not under relaxed selective constraints and that the faster-X effect is not a result of relaxed constraints.
The faster-X evolution of expression is most pronounced for genes that are unbound by the DCC, in transcriptionally repressive chromatin, or narrowly expressed (Figure 5A, Figure 6, Figure 8). If the faster-X effect were the result of relaxed selective constraints, we should observe increased An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e048.jpg values in genes with the most pronounced faster-X effect. Consistent with this prediction, X-linked genes that are unbound by the DCC have higher An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e049.jpg than bound genes (Figure 9C), suggesting that unbound genes are under relaxed constraints. Comparisons between X-linked and autosomal genes with similar expression breadth or in similar chromatin environments, however, suggest that the faster-X effect is not the result of relaxed constraints. For example, while narrowly expressed genes have elevated An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e050.jpg values, there is not a significant difference in An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e051.jpg between X-linked and autosomal narrowly expressed genes (Figure 9D). Additionally, X-linked genes in repressive chromatin tend to have lower An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e052.jpg than autosomal genes in repressive chromatin (Figure 9E). Faster-X expression evolution is therefore unlikely to be a result of relaxed selective constraints on X-linked expression levels.
We showed that X-linked gene expression levels in Drosophila have more interspecific divergence than autosomal expression levels (Figure 1, Figure 2, Figure 4), demonstrating faster-X evolution of gene expression in this genus. A similar faster-X effect has been observed for expression levels in Drosophila embryos (Kayserili, Gerrard, Tomancak, and Kalinka, in review, http://arxiv.org/abs/1209.0968). The faster-X effect is most pronounced for genes that are unbound by the DCC (Figure 5), in transcriptionally repressive chromatin (Figure 6), and narrowly expressed in a limited number of tissues (Figure 8), as summarized in Table 1. The expression levels of X-linked genes are not more variable within species (Figure 9), suggesting that the faster-X evolution of expression is not the result of relaxed selective constraints. We therefore hypothesize that the faster-X effect is driven by a higher rate of adaptive substitutions that affect the expression of X-linked genes relative to those that affect autosomal gene expression levels. Below, we expand on this hypothesis, explaining how allelic dominance, dosage compensation, and population size may contribute to the faster-X evolution of gene expression in Drosophila.
Table 1
Table 1
Factors that contribute to the faster-X evolution of gene expression.
Adaptive evolution of DCC proteins and HASs
Our analysis relies on inferences of DCC binding and HASs based on experiments performed in D. melanogaster cells [43], [45]. DCC proteins have experienced adaptive evolution along the D. melanogaster lineage [27], [77], as have three HASs on the X chromosome [78]. Despite the potential for enhanced expression divergence because of this rapid evolution, we see greater conservation of expression levels associated with DCC bound genes (Figure 5A). This result implies that, despite the accelerated evolution of DCC components and their binding sites, DCC binding is likely to be conserved across species. Furthermore, if DCC binding sites are turning over, this makes our discovery of a relationship between DCC binding and expression divergence conservative.
Chromatin environment, expression breadth, and the faster-X effect
DCC binding and chromatin states were inferred in a limited number of D. melanogaster cell lines. Genes that were never identified as bound by the DCC are either never compensated (because their dose does not need to be tightly controlled) or are compensated in cell types other than those studied thus far. Similarly, genes categorized as in regions of transcriptionally repressive chromatin are likely to be transcriptionally activated in a tissue-specific manner that differs from their state in S2 or BG3 cells. In this way, chromatin state can be used as a second, independent measurement of expression breadth: genes inferred to be in repressive chromatin can be assumed to be narrowly expressed, while genes in active chromatin are likely to be broadly expressed (Figure 7C).
Narrowly expressed genes tend to have rapidly evolving protein coding sequences, possibly because they are under fewer evolutionary constraints [19], [66]. Not only can these relaxed constraints permit faster evolution by purely neutral processes, but genes that are under fewer constraints are also expected to have a higher likelihood of adaptive fixations because they have less pleiotropic restrictions on their evolution [79]. Similarly, genes that are unbound by the DCC, genes that are in transcriptionally repressive chromatin, and narrowly expressed genes have rapidly evolving expression levels (Figure 5, Figure 6, Figure 7). In addition, these genes also have more intraspecific variation in expression (Figure 9), as do genes with fewer annotated functions [80], suggesting that the regulation of their expression is under relaxed constraints. If the faster-X evolution of gene expression is driven by positive selection, the faster-X effect should be most pronounced in genes that are most likely to experience adaptive substitutions. Genes in transcriptionally repressive chromatin and narrowly expressed genes do indeed have the most robust evidence for a faster-X effect (Figure 6, Figure 8), supporting an adaptive model of faster-X expression evolution.
A framework for adaptation-driven faster-X evolution of gene expression
The canonical model of faster-X evolution driven by positive selection posits that X-linked recessive beneficial mutations will be exposed to selection in hemizygous males, this will lead to an increased probability of invasion for X-linked beneficial alleles, and there will be a higher rate of adaptive evolution in X-linked genes [1], [3]. Before we can apply this model to the faster-X evolution of gene expression, we must determine if two assumptions are met: 1) mutations that affect the expression of X-linked genes are themselves X-linked; 2) mutations that affect expression level are recessive. We consider each of these assumptions below, and then we develop a conceptual framework to explain the faster-X evolution of expression.
Gene expression is inherited in a polygenic manner [81], and both cis and trans acting factors are responsible for expression differences between Drosophila species [82], [83]. The expression divergence of X-linked genes is therefore determined by substitutions in X-linked cis regulatory sequences and the trans acting proteins that bind to them. While the cis regulatory elements are all X-linked, the trans factors can be encoded by X-linked or autosomal genes. There are trans factors that preferentially affect the expression of X-linked genes (e.g., some nuclear pore proteins [84], [85], the DCC [35], and chromatin modifications that are enriched on the X chromosome [64]), but these are unlikely to be the norm [65]. In addition, the preferential targeting of certain trans factors to X-linked loci is ultimately attributable to sequences that are enriched on the X chromosome—either X-linked motifs direct trans factors to the X chromosome [43], [44] or trans factors are attracted by other proteins that are enriched on the X chromosome because they themselves were directed there by X-linked sequences [85]. Similarly, the expression of autosomal genes is determined by X-linked and autosomal trans acting factors, but the cis regulatory sequences are all autosomal. It is well documented that cis regulatory changes play an important role in gene expression divergence [82], [83], [86]–[88]. Therefore, expression changes in X-linked genes are more likely to be the result of mutations in X-linked loci when compared to similar expression changes in autosomal genes. This supports the hypothesis that the faster-X evolution of gene expression is the result of an increased rate of X-linked substitutions affecting expression levels.
If the faster-X evolution of expression is driven by adaptive substitutions in the cis regulatory sequences of X-linked genes, we would expect to detect signatures of positive selection near genes on the X chromosome. X-linked loci in D. melanogaster do tend to have reduced genetic variation, and this can be attributed to genetic hitchhiking because of selection at loci within or near X-linked genes [13], [89], [90]. In addition, DNA sequence variation is positively correlated with intraspecific expression variation in D. simulans [76], and sequence divergence upstream of coding sequences is correlated with expression divergence between D. melanogaster and a close relative [83]. These patterns further support the hypothesis that the faster-X evolution of gene expression is driven by X-linked substitutions that affect expression level in cis.
While it is clear how X-linked mutations can affect the expression of X-linked genes, it is not obvious why those mutations would be recessive. Non-additive inheritance of gene expression levels is common [91], [92], but cis regulatory differences between species are more likely to be inherited in an additive manner [83], [93]. These results suggest that the phenotypic effects of mutations that affect expression in cis are not likely to be recessive, but what is more important is whether the fitness effects of the mutations are recessive. It is reasonable to assume that the fitness landscape near an optimum is concave, which implies that mutations that push the expression level of a gene toward the optimum will be dominant [94], [95]. Therefore, empirical results and theoretical predictions appear to challenge the assumption that beneficial mutations that affect expression in cis will be recessive.
Recently, however, it has been demonstrated that beneficial mutations with additive phenotypic effects can increase fitness in heterozygotes while at the same time being less fit when homozygous because they overshoot the adaptive peak [96], [97] (Figure 10A). Therefore, mutations with additive phenotypic effects can have over-dominant fitness effects as a consequence of diploidy [97]. While this could impede adaptation at autosomal loci, the dynamics of this process are likely to differ at X-linked genes because they are effectively haploid in males in the absence of dosage compensation (Figure 10B). Notably, we detect the faster-X effect in genes that appear not be dosage compensation (Figure 5). Beneficial mutations with additive phenotypic effects on the expression of these uncompensated X-linked genes may therefore be more likely to fix because selection in males does not run the risk of overshooting the fitness optimum as a consequence of diploidy (Figure 10). Further theoretical work is needed to determine whether this intuitive prediction is a feasible adaptive explanation for the faster-X evolution of gene expression.
Figure 10
Figure 10
Additive phenotypic effects, fitness landscapes, and X-linkage.
If the aforementioned model could explain the faster-X evolution of gene expression, we would expect the faster-X effect to be limited to genes expressed in males because they would be present in the hemizygous (i.e., haploid) state. Consistent with this hypothesis, the expression levels of X-linked genes transcribed primarily in female-limited tissues do not evolve faster than autosomal genes with equivalent expression profiles (Figure 3). We do, however, detect the faster-X effect when expression is measured in either males or females (Figure 1, Figure 2, Figure 4, Figure 6, Figure 8). This is because male and female phenotypes are correlated so that selection on expression levels in males will affect the expression levels of genes that are also expressed in females [98], [99].
Lineage-specific patterns of faster-X expression evolution in Drosophila
Within the Drosophila genus, we observe the most pronounced faster-X effect along the lineage leading to D. ananassae (Figure 4, Figure 6, Figure 8). The Painting of fourth (POF) protein localizes specifically to the diminutive D. melanogaster fourth chromosome (Muller element F, or dot chromosome) [100]. Numerous lines of evidence suggest that POF promotes the transcriptional output of genes on the fourth chromosome by an unknown RNA-binding mechanism [101]–[103]. Intriguingly, while POF is dot-chromosome-specific in most Drosophila species, it also co-localizes with the DCC on the X chromosome in D. ananassae males [104].
POF localization to the D. ananassae X chromosome in males suggests that X-linked gene expression is uniquely affected in D. ananassae. This could contribute to the increased expression divergence of the X chromosome along the D. ananassae lineage either by directly affecting the expression levels of D. ananassae genes or by creating unique selection pressures on X-linked gene expression. We detect a substantial faster-X effect in female expression along the D. ananassae lineage (Figure 4A), despite the fact that POF does not localize to the X chromosome in D. ananassae females [104]. Therefore, the accentuated faster-X effect along the lineage leading to D. ananassae is unlikely to be a direct result of POF modifying expression levels. It is instead more likely that POF localization to the X chromosome in D. ananassae creates novel selection pressures on X-linked expression levels, which leads to a more pronounced faster-X effect.
Faster-X expression evolution in Drosophila and mammals
Expression levels of X-linked genes diverge faster than those of autosomal genes along both internal and terminal branches of the Drosophila phylogeny (Figure 4). The faster-X effect in mammals, on the other hand, is limited to only some deep lineages [28]. If our conceptual framework for understanding the faster-X evolution of gene expression is correct, we should be able to use it to explain differences in the faster-X effect between Drosophila and mammals. We consider four hypotheses that could explain the extent of faster-X expression evolution in the two taxa, and we examine how each could contribute to the observed incongruities.
First, X chromosome gene content differs between Drosophila and mammals [52], [105]. These differences, however, are unlikely to be responsible for the differences in the faster-X effect between mammals and Drosophila. The mammalian X chromosome harbors an excess of narrowly expressed genes [52], [106], i.e., the same type of genes with the most pronounced faster-X effect in Drosophila (Figure 8). Therefore, we would expect an even more substantial faster-X effect in mammals if differences in X chromosome gene content were an important contributor to taxon-specific faster-X expression evolution.
Second, Drosophila and mammals deal with the haploid dose of the male X chromosome in different ways [29], [37]. Faster-X gene expression evolution in Drosophila is most pronounced for genes that are unbound by the DCC (Figure 5), and we hypothesize that the effective haploidy of X-linked alleles in uncompensated Drosophila genes promotes the faster-X effect (Figure 10). Mammalian dosage compensation, on the other hand, is thought to be a two step process: X chromosome gene expression is upregulated in both sexes, followed by random silencing of one X chromosome in females [29]–[31]. The specific mechanisms of mammalian X chromosome upregulation are not understood, and the phenomenon itself remains controversial [32]. Regardless of the details of mammalian dosage compensation, if the allelic dominance of fitness effects for mutations that change gene expression are affected by the mechanism of dosage compensation, the differences in dosage compensation between Drosophila and mammals could be responsible for the taxon-specific patterns of faster-X expression evolution [1], [3].
Third, the rate of evolution depends on mutational input and the fixation rate of those mutations. A higher autosomal mutation rate could therefore counteract a higher fixation rate on the X chromosome [107]. The mutation rate in the germline of many mammals is higher in males than females (“male mutation bias”) [4]–[8]. Because the X chromosome is transmitted through females 2/3 of the time, the population mutation rate is lower for the X chromosome than the autosomes in species with male mutation bias. This downward biased mutation rate of the X chromosome in some mammalian lineages could therefore be responsible for the lineage-specificity of the faster-X effect in mammals [107], [108].
Fourth, if the faster-X evolution of gene expression is driven by adaptive substitutions, as we propose, it is likely to be sensitive to An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e056.jpg [9], [10]. In small populations a larger fraction of mutations will be effectively neutral [109], which will decrease the number of beneficial mutations fixed by positive selection. The higher An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e057.jpg of Drosophila relative to mammals may therefore be more permissive of adaptive faster-X evolution [10].
In summary, we conclude that the difference in the extent of the faster-X evolution of gene expression between Drosophila and mammals could be a result of the unique mechanism of dosage compensation in Drosophila, the pervasiveness of male mutation bias, and/or the differences in An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e058.jpg between taxa. Determining which factor is most important will require additional theoretical and empirical work to identify the key determinants of gene expression evolution, the nature of selection on expression, and the effects of gene dosage on the dominance of fitness effects.
Expression measurements
We obtained microarray measurements of expression from whole fly or head from previously published results [50]–[52]. We calculated the expression level of each gene by first taking the median signal across all probes for each gene within each replicate, and then calculating the median for each gene across all replicates. As an alternative approach, we used expression levels estimated in the LIMMA package of Bioconductor [110], as described previously [52].
We tested for significant differences in expression between males and females (i.e., sex-biased expression) using moderated t-tests implemented in the LIMMA package with the empirical Bayes function to pool sample variances toward a common value [110], as described previously [52]. We corrected for multiple tests using a FDR [111] of 0.05 when only sex-biased expression in D. melanogaster was considered, and with a FDR of 0.20 when sex-biased expression in all species was considered. Genes with significantly higher expression in males were classified as having male-biased expression, and those with higher expression in females as female-biased.
RNA-seq data collected from whole D. melanogaster, D. pseudoobscura, and D. mojavensis males and female or D. melanogaster, D. pseudoobscura, D. willistoni, and D. mojavensis heads were obtained from previously published results [52]. Reads longer than 36 bases (bp) were trimmed to 36 bp so that all reads were the same length, and reads were then mapped to the transcriptome using BWA [112]. Any read mapping to multiple locations in the genome was discarded, and genes with fewer than 50 mapped reads were excluded from the subsequent analysis. The expression level of each gene was estimated as the number of unique reads mapping to the gene standardized by the total number of mapped reads and the transcript length.
We extracted only those genes present as 1[ratio]1 orthologs on the same chromosome arm in all species under consideration, and we quantile normalized the expression levels so that they are identically distributed across all species. We next calculated SpearmanÕs An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e059.jpg between all pairwise comparisons of expression levels from the species under consideration. This was repeated for each chromosome arm. Confidence intervals (CIs) of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e060.jpg were estimated by bootstrapping the data 1000 times in the R statistical computing environment [113]. We also calculated correlations of expression between sexes within each species.
Microarray measurements of expression were obtained for 14 adult D. melanogaster tissues from FlyAtlas [58], and the expression breadth was determined for each gene as described previously [19]. Briefly, we calculated An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e061.jpg [114], a metric that ranges between 0 (for broadly expressed genes) to 1 (for narrowly expressed genes). Genes were said to be narrowly expressed in a tissue if An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e062.jpg, and genes with An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e063.jpg were classified as broadly expressed.
Phylogenetic analysis
We used An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e064.jpg as an estimate of the pairwise divergence in expression between D. melanogaster, D. yakuba, D. ananassae, D. pseudoobscura, D. mojavensis, and D. virilis. We then reconstructed the phylogenetic relationships using the method of Fitch and Margoliash [115] implemented in the PHYLIP software package [116]. Bootstrap support for phylogenetic nodes was estimated by resampling the 1[ratio]1[ratio]1[ratio]1[ratio]1[ratio]1 orthologs 1000 times. Branch lengths were estimated using the method of Fitch and Margoliash [115] with a fixed tree topology implemented in the PHYLIP software package [116]. CIs of the branch lengths were calculated by bootstrapping the data 1000 times. Bootstrap support and branch lengths were estimated for all 1[ratio]1[ratio]1[ratio]1[ratio]1[ratio]1 orthologs, and this was repeated for genes on each chromosome arm separately. All bootstrapping was performed using the R statistical computing environment [113].
DCC binding
We obtained a list of genes bound by the DCC identified using ChIP-chip in three cell types: SL2 embryonic cell culture, larval wing imaginal disc cell culture, and late embryo [45]. A gene was said to be bound by the DCC if it is bound in at least one cell type. In addition, a second list of DCC bound genes was kindly provided by D. Bachtrog [48]. Our results are robust to the gene list used in our analysis. X-linked regions identified as HASs were obtained from previously published results [43].
Gene-wise expression divergence and distance to nearest HAS
We calculated the gene-wise expression divergence between 1[ratio]1[ratio]1 orthologs in D. melanogaster, D. yakuba, and D. ananassae as:
A mathematical equation, expression, or formula. Object name is pgen.1003013.e065.jpg
where An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e066.jpg and An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e067.jpg are the expression levels of ortholog An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e068.jpg in species An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e069.jpg and An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e070.jpg, respectively. We then calculated Spearman's An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e071.jpg between An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e072.jpg, An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e073.jpg, and distance to the nearest HAS (where An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e074.jpg = D. melanogaster). From these pairwise correlations, we calculated partial correlations to determine the direct relationship of each pair of values [117]. The CIs of the pairwise and partial correlations were estimated by bootstrapping the 1[ratio]1[ratio]1 orthologs 1000 times in the R statistical programming environment [113].
Chromatin states
Kharchenko et al. [65] analyzed ChIP-chip results for multiple histone modifications and DNA binding proteins in two D. melanogaster cell lines (S2 and BG3), and they used a hidden Markov model to assign each region of the genome to one of nine chromatin states. States 1–5 are associated with transcriptionally active chromatin marks and states 6–9 with repressive marks. We used the overlap of these regions with annotated protein coding genes to determine whether each D. melanogaster gene is associated with a region of active or repressive chromatin marks. A gene was considered to be in active chromatin if An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e075.jpg of the gene body overlaps with regions identified as containing active marks, and, conversely, a gene was considered to be in repressive chromatin in An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e076.jpg of the gene body overlaps with regions identified as harboring repressive marks. We obtain similar results when use overlap cutoffs of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e077.jpg and An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e078.jpg.
Intraspecific expression variation
We obtained estimates of broad sense heritability (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e079.jpg) for D. melanogaster genes from a published analysis of microarray expression measurements in females and males from 40 inbred lines [75]. These estimates were calculated using an analysis of variance (ANOVA), and we considered only genes in which the line term (the estimate of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e080.jpg) was significant at a FDR of 0.05.
Figure S1
Faster-X evolution of gene expression using an alternative analysis of microarray data. Pairwise correlations of gene expression are shown for genes on each chromosome arm, using expression measurements from females (top) or males (bottom). Expression levels were estimated using the LIMMA package of Bioconductor [110] (see Methods). See Figure 1 for a description of the graphs.
(PDF)
Click here for additional data file.(257K, pdf)
Figure S2
Faster-X evolution of gene expression using RNA-seq data. Pairwise correlations of gene expression are shown for genes on each chromosome arm, using expression measurements from females (top) or males (bottom) collected by RNA-seq. See Figure 1 for a description of the graphs.
(PDF)
Click here for additional data file.(134K, pdf)
Figure S3
Faster-X evolution of gene expression using microarray data with only the genes included in the RNA-seq data set. Pairwise correlations of gene expression are shown for genes on each chromosome arm, using expression measurements from females (top) or males (bottom) collected by microarray. Only those genes present in both the microarray and RNA-seq data sets are included. See Figure 1 for a description of the graphs.
(PDF)
Click here for additional data file.(121K, pdf)
Figure S4
Faster-X evolution of gene expression measured in head using RNA-seq. Pairwise correlations of gene expression are shown for genes on each chromosome arm, using expression measurements from female (top) or male (bottom) heads collected with RNA-seq. See Figure 1 for a description of the graphs.
(PDF)
Click here for additional data file.(162K, pdf)
Figure S5
Faster-X evolution of gene expression with D. melanogaster male-biased genes removed. Pairwise correlations of gene expression are shown for genes on each chromosome arm, using expression measurements from females (top) or males (bottom). Genes with male-biased expression in D. melanogaster were excluded. See Figure 1 for a description of the graphs.
(PDF)
Click here for additional data file.(265K, pdf)
Figure S6
Faster-X evolution of gene expression with male-biased genes removed. Pairwise correlations of gene expression are shown for genes on each chromosome arm, using expression measurements from females (top) or males (bottom). Genes with male-biased expression in any of the six species were excluded. See Figure 1 for a description of the graphs.
(PDF)
Click here for additional data file.(267K, pdf)
Figure S7
Faster-X evolution of gene expression with genes with biased expression in male reproductive tissues removed. Pairwise correlations of gene expression are shown for genes on each chromosome arm, using expression measurements from females (top) or males (bottom). Genes with biased expression in either testis or accessory gland were excluded. See Figure 1 for a description of the graphs.
(PDF)
Click here for additional data file.(267K, pdf)
Figure S8
New genes do not exhibit faster-X evolution of gene expression Pairwise correlations of gene expression between D. melanogaster and D. yakuba (mel-yak) and D. melanogaster and D. ananassae (mel-ana) for genes on each chromosome arm. Correlations were calculated using genes shared by all species in the Drosophila genus and genes that arose along the lineage leading to D. melanogaster after the split with the Drosophila subgenus [61].
(PDF)
Click here for additional data file.(132K, pdf)
Figure S9
Phylogenetic reconstruction using the correlation of expression. We reconstructed the evolutionary relationships of the six species using the pairwise correlation of expression levels in (A) females and (B) males. A distance matrix of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e081.jpg was used to build the phylogenies using the Fitch and Margoliash [115] method. We bootstrap sampled the genes 1000 times to estimate the support for each node, and the percent of bootstrap replicates supporting each node is given on the tree.
(PDF)
Click here for additional data file.(257K, pdf)
Figure S10
Expression levels of genes in transcriptionally active and repressive chromatin Boxplots show the distribution of expression levels for D. melanogaster genes in transcriptionally active and repressive (repress) chromatin. Horizontal dashed lines represent the genome-wide expression level. Chromatin states were measured in two cell lines (BG3 and S2), and expression levels were measured in either females or males. In addition, genes were divided into those that are autosomal and those that are X-linked. Significant differences between the expression levels of genes in active and repressive chromatin are indicated by asterisks (*** An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e082.jpg, ***** An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e083.jpg).
(PDF)
Click here for additional data file.(119K, pdf)
Figure S11
Faster expression evolution of un-dosage-compensated X-linked genes not associated with active chromatin. This graph is the same as the one in Figure 6 except that X-linked genes are now divided into those that are both bound by the DCC and in active chromatin or unbound by the DCC and in repressive chromatin.
(PDF)
Click here for additional data file.(202K, pdf)
Figure S12
Narrowly expressed genes have greater expression divergence. Spearman's rank order correlation between An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e084.jpg and expression divergence is plotted, along with the 95% CI of the correlation. Larger values of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e085.jpg indicate narrower expression breadth. Divergence was measured between D. melanogaster (mel) and both D. yakuba (yak) and D. ananassae (ana) using measurements from both females and males. The dashed line indicates the null expectation of no correlation.
(PDF)
Click here for additional data file.(244K, pdf)
Figure S13
DCC bound genes tend to be broadly expressed. X-linked genes were divided into those that are bound and unbound by the DCC. (A) Expression breadth (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e086.jpg) was compared between bound and unbound genes using a Mann-Whitney U test (***** An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e087.jpg). Larger values of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e088.jpg indicate narrower expression breadth. (B) SpearmanÕs rank order correlation between distance from the nearest HAS and An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e089.jpg is plotted, along with the 95% CI of the correlation. The dashed line indicates the null expectation of no correlation.
(PDF)
Click here for additional data file.(237K, pdf)
Figure S14
Genes in repressive chromatin are more narrowly expressed. Genes were divided into those that are in transcriptionally active and repressive chromatin using data from two cell lines: BG3 and S2. (A) Expression breadth (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e090.jpg) was compared between genes in active and repressive chromatin using a Mann-Whitney U test (***** An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e091.jpg). Larger values of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e092.jpg indicate narrower expression breadth. (B) Genes were additionally divided into those that are broadly and narrowly expressed, and the counts of genes in each expression breadth and chromatin state are plotted for data collected in BG3 cells (S2 data are available in Figure 7). Fisher's exact test was used to determine if there is a significantly non-random distribution of genes in the four classes.
(PDF)
Click here for additional data file.(257K, pdf)
Figure S15
The faster-X effect is limited to narrowly expressed genes in transcriptionally repressive chromatin, with chromatin environment measured in BG3 cells. This figure is the same as Figure 8 except chromatin state was measured in BG3 cells.
(PDF)
Click here for additional data file.(156K, pdf)
Figure S16
The faster-X effect is limited to narrowly expressed genes in transcriptionally repressive chromatin, with chromatin environment measured in BG3 cells. This figure is the same as Figure 8 and Figure S15 except that genes with the lowest 5% expression levels were excluded.
(PDF)
Click here for additional data file.(307K, pdf)
Figure S17
Expression breadth, chromatin environment, and X-linkage. The expression breadth (An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e093.jpg) of genes in transcriptionally active and repressive chromatin on the X chromosome and autosomes is plotted. Larger values of An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e094.jpg indicate narrower expression breadth. Chromatin state was inferred in BG3 and S2 cells. Significant differences between X-linked and autosomal genes in the same chromatin context are indicated by asterisks (* An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e095.jpg, *** An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e096.jpg, **** An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e097.jpg). (A) Median expression breadth across the entire genome is indicated by a dashed line. (B) Genes were additionally divided into narrowly and broadly expressed.
(PDF)
Click here for additional data file.(278K, pdf)
Figure S18
No faster-X effect for genes narrowly expressed in female reproductive tissues. Boxplots show the pairwise divergence in expression between 1[ratio]1[ratio]1 orthologs in the D. melanogaster (mel) and D. yakuba (yak) or D. ananassae (ana) genomes measured in whole females and males (see Figure 6). Only genes that are narrowly expressed in non-reproductive tissues (non-repro), female reproductive tissues (female repro), or male reproductive tissues (male repro) are included. X-linked (X, red) and autosomal (A) genes were assigned to transcriptionally active and repressive chromatin based on the results of experiments in BG3 (top) and S2 (bottom) cells. The horizontal dashed line indicates the genome-wide average pairwise divergence for narrowly expressed genes. Subsets of narrowly expressed genes whose pairwise divergence is significantly different than all other narrowly expressed genes are marked (N). The same was done for subsets of narrowly expressed genes that reside in repressive chromatin (R). Significant differences between X-linked and autosomal genes in the same chromatin state and with the same expression breadth are marked with asterisks. Mann-Whitney U tests were used to assess significant differences (one symbol = An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e098.jpg, two symbols = An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e099.jpg, three symbols = An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e100.jpg, four symbols = An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e101.jpg, five symbols = An external file that holds a picture, illustration, etc. Object name is pgen.1003013.e102.jpg).
(PDF)
Click here for additional data file.(315K, pdf)
Table S1
DCC binding and chromatin state. Genes were called as bound by the DCC if they were bound in S2 cells (SL2); either SL2, larval wing imaginal disc, or late embryonic cells (Any); or based on a separate analysis of the data (Bachtrog) [48]. Chromatin states (chr.state) were inferred in S2 cells or BG3 cells. The number of genes in each DCC binding and chromatin state class (num_genes) were tested for independence using Fisher's exact test, and the P value of this test is reported (FET.p).
(PDF)
Click here for additional data file.(43K, pdf)
Acknowledgments
We are extremely grateful to Tim Connallon for discussions about the theory underlying faster-X evolution. We thank Peter Park and Doris Bachtrog for sharing their data on DCC binding, Peter Kharchenko and Peter Park for sharing their chromatin state calls, and Julien Ayroles for sharing the estimates of expression heritability. This manuscript benefited from comments by Tim Connallon, Julian Ayroles, other members of the Clark lab, Erin Kelleher, and four anonymous reviewers.
Funding Statement
RPM was supported by National Institutes of Health fellowship F32GM087611. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
References
1. Charlesworth B, Coyne JA, Barton NH (1987) The relative rates of evolution of sex chromosomes and autosomes. Am Nat 130: 113–146.
2. Orr HA, Betancourt AJ (2001) Haldane's sieve and adaptation from the standing genetic variation. Genetics 157: 875–884. [PubMed]
3. Connallon T, Singh ND, Clark AG (2012) Impact of genetic architecture on the relative rates of X versus autosomal adaptive substitution. Mol Biol Evol 29: 1933–1942. [PubMed]
4. Haldane JBS (1935) The rate of spontaneous mutation of a human gene. J Genet 31: 317–326.
5. Miyata T, Hayashida H, Kuma K, Mitsuyasu K, Yasunaga T (1987) Male-driven molecular evolution: a model and nucleotide sequence analysis. Cold Spring Harb Symp Quant Biol 52: 863–867. [PubMed]
6. Li WH, Yi S, Makova K (2002) Male-driven evolution. Curr Opin Genet Dev 12: 650–656. [PubMed]
7. Ellegren H (2007) Characteristics, causes and evolutionary consequences of male-biased mutation. Proc R Soc Lond, B, Biol Sci 274: 1–10. [PMC free article] [PubMed]
8. Wilson Sayres MA, Makova KD (2011) Genome analyses substantiate male mutation bias in many species. BioEssays 33: 938–945. [PubMed]
9. Vicoso B, Charlesworth B (2009) Effective population size and the faster-X effect: an extended model. Evolution; international journal of organic evolution 63: 2413–26. [PubMed]
10. Mank JE, Vicoso B, Berlin S, Charlesworth B (2010) Effective population size and the Faster-X effect: empirical results and their interpretation. Evolution 64: 663–674. [PubMed]
11. Lu J, Wu CI (2005) Weak selection revealed by the whole-genome comparison of the X chromosome and autosomes of human and chimpanzee. Proc Natl Acad Sci USA 102: 4063–4067. [PubMed]
12. Nielsen R, Bustamante C, Clark AG, Glanowski S, Sackton TB, et al. (2005) A scan for positively selected genes in the genomes of humans and chimpanzees. PLoS Biol 3: e170 doi:10.1371/journal.pbio.0030170. [PMC free article] [PubMed]
13. Mackay TFC, Richards S, Stone EA, Barbadilla A, Ayroles JF, et al. (2012) The Drosophila melanogaster genetic reference panel. Nature 482: 173–178. [PubMed]
14. Hvilsom C, Qian Y, Bataillon T, Li Y, Mailund T, et al. (2012) Extensive X-linked adaptive evolution in central chimpanzees. Proc Natl Acad Sci USA 109: 2054–2059. [PubMed]
15. Begun DJ, Holloway AK, Stevens K, Hillier LW, Poh YP, et al. (2007) Population genomics: whole-genome analysis of polymorphism and divergence in Drosophila simulans. PLoS Biol 5: e310 doi:10.1371/journal.pbio.0050310. [PMC free article] [PubMed]
16. Torgerson DG, Singh RS (2003) Sex-linked mammalian sperm proteins evolve faster than autosomal ones. Mol Biol Evol 20: 1705–1709. [PubMed]
17. Torgerson DG, Singh RS (2006) Enhanced adaptive evolution of sperm-expressed genes on the mammalian X chromosome. Heredity 96: 39–44. [PubMed]
18. Baines JF, Sawyer SA, Hartl DL, Parsch J (2008) Effects of X-linkage and sex-biased gene expression on the rate of adaptive protein evolution in Drosophila. Mol Biol Evol 25: 1639–1650. [PMC free article] [PubMed]
19. Meisel RP (2011) Towards a more nuanced understanding of the relationship between sex-biased gene expression and rates of protein coding sequence evolution. Mol Biol Evol 28: 1893–1900. [PMC free article] [PubMed]
20. Grath S, Parsch J (2012) Rate of amino acid substitution is inuenced by the degree and conservation of male-biased transcription over 50 million years of Drosophila evolution. Genome Biol Evol 4: 346–359. [PMC free article] [PubMed]
21. Coyne JA, Orr HA (1989) The two rules of speciation. In: Otte D, Endler J, editors, Speciation and its Consequences, Sinauer Associates. pp. 180–207.
22. Masly JP, Presgraves DC (2007) High-Resolution Genome-Wide Dissection of the Two Rules of Speciation in Drosophila. PLoS Biol 5: e243 doi:10.1371/journal.pbio.0050243. [PMC free article] [PubMed]
23. Presgraves DC (2008) Sex chromosomes and speciation in Drosophila. Trends Genet 24: 336–343. [PMC free article] [PubMed]
24. Good JM, Giger T, Dean MD, Nachman MW (2010) Widespread over-expression of the X chromosome in sterile F1 hybrid mice. PLoS Genet 6: e1001148 doi:10.1371/journal.pgen.1001148. [PMC free article] [PubMed]
25. Pal Bhadra M, U B, Birchler JA (2006) Misregulation of sex-lethal and disruption of male-specific lethal complex localization in Drosophila species hybrids. Genetics 174: 1151–1159. [PubMed]
26. Chatterjee RN, Chatterjee P, Pal A, Pal-Bhadra M (2007) Drosophila simulans Lethal hybrid rescue mutation (Lhr) rescues inviable hybrids by restoring X chromosomal dosage compensation and causes uctuating asymmetry of development. J Genet 86: 203–215. [PubMed]
27. Rodriguez MA, Vermaak D, Bayes JJ, Malik HS (2007) Species-specific positive selection of the male-specific lethal complex that participates in dosage compensation in Drosophila. Proc Natl Acad Sci USA 104: 15412–15417. [PubMed]
28. Brawand D, Soumillon M, Necsulea A, Julien P, Csárdi G, et al. (2011) The evolution of gene expression levels in mammalian organs. Nature 478: 343–348. [PubMed]
29. Straub T, Becker PB (2007) Dosage compensation: the beginning and end of generalization. Nat Rev Genet 8: 47–57. [PubMed]
30. Payer B, Lee JT (2008) X chromosome dosage compensation: how mammals keep the balance. Annu Rev Genet 42: 733–772. [PubMed]
31. Deng X, Hiatt JB, Nguyen DK, Ercan S, Sturgill D, et al. (2011) Evidence for compensatory upregulation of expressed X-linked genes in mammals, Caenorhabditis elegans and Drosophila melanogaster. Nat Genet 43: 1179–1185. [PubMed]
32. Lin F, Xing K, Zhang J, He X (2012) Expression reduction in mammalian X chromosome evolution refutes Ohno's hypothesis of dosage compensation. Proc Natl Acad Sci USA 109: 11752–11757. [PubMed]
33. Bone JR, Lavender J, Richman R, Palmer MJ, Turner BM, et al. (1994) Acetylated histone H4 on the male X chromosome is associated with dosage compensation in Drosophila. Genes Dev 8: 96–104. [PubMed]
34. Gelbart ME, Larschan E, Peng S, Park PJ, Kuroda MI (2009) Drosophila MSL complex globally acetylates H4K16 on the male X chromosome for dosage compensation. Nat Struct Mol Biol 16: 825–832. [PMC free article] [PubMed]
35. Conrad T, Akhtar A (2012) Dosage compensation in Drosophila melanogaster: epigenetic finetuning of chromosome-wide transcription. Nat Rev Genet 13: 123–134. [PubMed]
36. Park Y, Kuroda MI (2001) Epigenetic aspects of X-chromosome dosage compensation. Science 293: 1083–1085. [PubMed]
37. Vicoso B, Bachtrog D (2009) Progress and prospects toward our understanding of the evolution of dosage compensation. Chromosome Res 17: 585–602. [PMC free article] [PubMed]
38. Conrad T, Cavalli FMG, Vaquerizas JM, Luscombe NM, Akhtar A (2012) Drosophila dosage compensation involves enhanced Pol II recruitment to male X-linked promoters. Science In press. [PubMed]
39. Larschan E, Bishop EP, Kharchenko PV, Core LJ, Lis JT, et al. (2011) X chromosome dosage compensation via enhanced transcriptional elongation in Drosophila. Nature 471: 115–118. [PMC free article] [PubMed]
40. Bashaw GJ, Baker BS (1995) The msl-2 dosage compensation gene of Drosophila encodes a putative DNA-binding protein whose expression is sex specifically regulated by Sex-lethal. Development 121: 3245–3258. [PubMed]
41. Kelley RL, Solovyeva I, Lyman LM, Richman R, Solovyev V, et al. (1995) Expression of Msl-2 causes assembly of dosage compensation regulators on the X chromosomes and female lethality in Drosophila. Cell 81: 867–877. [PubMed]
42. Zhou S, Yang Y, Scott MJ, Pannuti A, Fehr KC, et al. (1995) Male-specific lethal 2, a dosage compensation gene of Drosophila, undergoes sex-specific regulation and encodes a protein with a RING finger and a metallothionein-like cysteine cluster. EMBO J 14: 2884–2895. [PubMed]
43. Alekseyenko AA, Peng S, Larschan E, Gorchakov AA, Lee OK, et al. (2008) A sequence motif within chromatin entry sites directs MSL establishment on the Drosophila X chromosome. Cell 134: 599–609. [PMC free article] [PubMed]
44. Straub T, Grimaud C, Gilfillan GD, Mitterweger A, Becker PB (2008) The chromosomal highaffinity binding sites for the Drosophila dosage compensation complex. PLoS Genet 4: e1000302 doi:10.1371/journal.pgen.1000302. [PMC free article] [PubMed]
45. Alekseyenko AA, Larschan E, Lai WR, Park PJ, Kuroda MI (2006) High-resolution ChIP-chip analysis reveals that the Drosophila MSL complex selectively identifies active genes on the male X chromosome. Genes Dev 20: 848–857. [PubMed]
46. Gilfillan GD, Straub T, De Wit E, Greil F, Lamm R, et al. (2006) Chromosome-wide gene-specific targeting of the Drosophila dosage compensation complex. Genes Dev 20: 858–870. [PubMed]
47. Gorchakov AA, Alekseyenko AA, Kharchenko P, Park PJ, Kuroda MI (2009) Long-range spreading of dosage compensation in Drosophila captures transcribed autosomal genes inserted on X. Genes Dev 23: 2266–2271. [PubMed]
48. Bachtrog D, Toda NR, Lockton S (2010) Dosage compensation and demasculinization of X chromosomes in Drosophila. Curr Biol 20: 1476–1481. [PubMed]
49. Alekseyenko AA, Ho JWK, Peng S, Gelbart M, Tolstorukov MY, et al. (2012) Sequence-specific targeting of dosage compensation in Drosophila favors an active chromatin context. PLoS Genet 8: e1002646 doi:10.1371/journal.pgen.1002646. [PMC free article] [PubMed]
50. Parisi M, Nuttall R, Naiman D, Bouffard G, Malley J, et al. (2003) Paucity of genes on the Drosophila X chromosome showing male-biased expression. Science 299: 697–700. [PMC free article] [PubMed]
51. Zhang Y, Sturgill D, Parisi M, Kumar S, Oliver B (2007) Constraint and turnover in sex-biased gene expression in the genus Drosophila. Nature 450: 233–237. [PMC free article] [PubMed]
52. Meisel RP, Malone JH, Clark AG (2012) Disentangling the relationship between sex-biased gene expression and X-linkage. Genome Res 22: 1255–1265. [PubMed]
53. Drosophila 12 Genomes Consortium (2007) Evolution of genes and genomes on the Drosophila phylogeny. Nature 450: 203–218. [PubMed]
54. Muller HJ (1940) Bearings of the ‘Drosophila’ work on systematics. In: Huxley J, editor, The New Systematics, Oxford: Clarendon Press, Oxford. pp. 185–268.
55. Wang Z, Gerstein M, Snyder M (2009) RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10: 57–63. [PMC free article] [PubMed]
56. Ranz JM, Castillo-Davis CI, Meiklejohn CD, Hartl DL (2003) Sex-dependent gene expression and evolution of the Drosophila transcriptome. Science 300: 1742–1745. [PubMed]
57. Meiklejohn CD, Parsch J, Ranz JM, Hartl DL (2003) Rapid evolution of male-biased gene expression in Drosophila. Proc Natl Acad Sci USA 100: 9894–9899. [PubMed]
58. Chintapalli VR, Wang J, Dow JAT (2007) Using FlyAtlas to identify better Drosophila melano-gaster models of human disease. Nat Genet 39: 715–720. [PubMed]
59. Betrán E, Thornton K, Long M (2002) Retroposed new genes out of the X in Drosophila. Genome Res 12: 1854–1859. [PubMed]
60. Meisel RP, Han MV, Hahn MW (2009) A complex suite of forces drives gene traffic from Drosophila X chromosomes. Genome Biol Evol 1: 176–188. [PMC free article] [PubMed]
61. Zhang YE, Vibranovski MD, Krinsky BH, Long M (2010) Age-dependent chromosomal distribution of male-biased genes in Drosophila. Genome Res 20: 1526–1533. [PubMed]
62. Stone WS (1955) Genetic and chromosomal variability in Drosophila. Cold Spring Harb Symp Quant Biol 20: 256–270. [PubMed]
63. Tamura K, Subramanian S, Kumar S (2004) Temporal patterns of fruit y (Drosophila) evolution revealed by mutation clocks. Mol Biol Evol 21: 36–44. [PubMed]
64. Zhang Y, Oliver B (2010) An evolutionary consequence of dosage compensation on Drosophila melanogaster female X-chromatin structure? BMC Genomics 11: 6. [PMC free article] [PubMed]
65. Kharchenko PV, Alekseyenko AA, Schwartz YB, Minoda A, Riddle NC, et al. (2011) Comprehensive analysis of the chromatin landscape in Drosophila melanogaster. Nature 471: 480–485. [PMC free article] [PubMed]
66. Larracuente AM, Sackton TB, Greenberg AJ, Wong A, Singh ND, et al. (2008) Evolution of protein-coding genes in Drosophila. Trends Genet 24: 114–123. [PubMed]
67. Nuzhdin SV, Wayne ML, Harmon KL, McIntyre LM (2004) Common pattern of evolution of gene expression level and protein sequence in Drosophila. Mol Biol Evol 21: 1308–1317. [PubMed]
68. Lemos B, Bettencourt BR, Meiklejohn CD, Hartl DL (2005) Evolution of proteins and gene expression levels are coupled in Drosophila and are independently associated with mRNA abundance, protein length, and number of protein-protein interactions. Mol Biol Evol 22: 1345–1354. [PubMed]
69. Good JM, Hayden CA, Wheeler TJ (2006) Adaptive protein evolution and regulatory divergence in Drosophila. Mol Biol Evol 23: 1101–1103. [PubMed]
70. Mikhaylova LM, Nurminsky DI (2011) Lack of global meiotic sex chromosome inactivation, and paucity of tissue-specific gene expression on the Drosophila X chromosome. BMC Biol 9: 29. [PMC free article] [PubMed]
71. Rifkin SA, Kim J, White KP (2003) Evolution of gene expression in the Drosophila melanogaster subgroup. Nat Genet 33: 138–144. [PubMed]
72. Lemos B, Meiklejohn CD, Caceres M, Hartl DL (2005) Rates of divergence in gene expression profiles of primates, mice, and ies: stabilizing selection and variability among functional categories. Evolution 59: 126–137. [PubMed]
73. Hudson RR, Kreitman M, Aguade M (1987) A test of neutral molecular evolution based on nucleotide data. Genetics 116: 153–159. [PubMed]
74. McDonald JH, Kreitman M (1991) Adaptive protein evolution at the Adh locus in Drosophila. Nature 351: 652–654. [PubMed]
75. Ayroles JF, Carbone MA, Stone EA, Jordan KW, Lyman RF, et al. (2009) Systems genetics of complex traits in Drosophila melanogaster. Nat Genet 41: 299–307. [PMC free article] [PubMed]
76. Lawniczak MK, Holloway AK, Begun DJ, Jones CD (2008) Genomic analysis of the relationship between gene expression variation and DNA polymorphism in Drosophila simulans. Genome Biol 9: R125. [PMC free article] [PubMed]
77. Levine MT, Holloway AK, Arshad U, Begun DJ (2007) Pervasive and largely lineage-specific adaptive protein evolution in the dosage compensation complex of Drosophila melanogaster. Genetics 177: 1959–1962. [PubMed]
78. Bachtrog D (2008) Positive selection at the binding sites of the male-specific lethal complex involved in dosage compensation in Drosophila. Genetics 180: 1123–1129. [PubMed]
79. Fisher RA (1930) The Genetical Theory of Natural Selection. Oxford: Clarendon Press.
80. Hutter S, Saminadin-Peter SS, Stephan W, Parsch J (2008) Gene expression variation in African and European populations of Drosophila melanogaster. Genome Biol 9: R12. [PMC free article] [PubMed]
81. Gibson G, Weir B (2005) The quantitative genetics of transcription. Trends Genet 21: 616–623. [PubMed]
82. Wittkopp PJ, Haerum BK, Clark AG (2004) Evolutionary changes in cis and trans gene regulation. Nature 430: 85–88. [PubMed]
83. McManus CJ, Coolon JD, Duff MO, Eipper-Mains J, Graveley BR, et al. (2010) Regulatory divergence in Drosophila revealed by mRNA-seq. Genome Res 20: 816–825. [PubMed]
84. Mendjan S, Taipale M, Kind J, Holz H, Gebhardt P, et al. (2006) Nuclear pore components are involved in the transcriptional regulation of dosage compensation in Drosophila. Mol Cel 21: 811–823. [PubMed]
85. Vaquirizas JM, Suyama R, Kind J, Miura K, Luscombe NM, et al. (2010) Nuclear pore proteins Nup153 and Megator define transcriptionally active regions in the Drosophila genome. PLoS Genet 6: e1000846 doi:10.1371/journal.pgen.1000846. [PMC free article] [PubMed]
86. Lemos B, Araripe LO, Fontanillas P, Hartl DL (2008) Dominance and the evolutionary accumulation of cis- and trans-effects on gene expression. Proceedings of the National Academy of Sciences of the United States of America 105: 14471–14476. [PubMed]
87. Wittkopp PJ, Haerum BK, Clark AG (2008) Regulatory changes underlying expression differences within and between Drosophila species. Nat Genet 40: 346–350. [PubMed]
88. Graze RM, McIntyre LM, Main BJ, Wayne ML, Nuzhdin SV (2009) Regulatory divergence in Drosophila melanogaster and D. simulans, a genomewide analysis of allele-specific expression. Genetics 183: 547–561. [PubMed]
89. Hutter S, Li H, Beisswanger S, De Lorenzo D, Stephan W (2007) Distinctly different sex ratios in African and European populations of Drosophila melanogaster inferred from chromosomewide single nucleotide polymorphism data. Genetics 177: 469–480. [PubMed]
90. Singh N, Macpherson JM, Jensen J, Petrov D (2007) Similar levels of X-linked and autosomal nucleotide variation in African and non-African populations of Drosophila melanogaster. BMC Evol Biol 7: 202. [PMC free article] [PubMed]
91. Gibson G, Riley-Berger R, Harshman L, Kopp A, Vacha S, et al. (2004) Extensive sex-specific nonadditivity of gene expression in Drosophila melanogaster. Genetics 167: 1791–1799. [PubMed]
92. Wayne ML, Pan YJ, Nuzhdin SV, McIntyre LM (2004) Additivity and trans-acting effects on gene expression in male Drosophila simulans. Genetics 168: 1413–1420. [PubMed]
93. Lemos B, Araripe LO, Hartl DL (2008) Polymorphic Y chromosomes harbor cryptic variation with manifold functional consequences. Science 319: 91–93. [PubMed]
94. Fry JD (2010) The genomic location of sexually antagontic variation: some cautionary comments. Evolution 64: 1510–1516. [PubMed]
95. Connallon T, Clark AG (2010) Sex linkage, sex-specific selection, and the role of recombination in the evolution of sexually dimorphic gene expression. Evolution 64: 3417–3442. [PMC free article] [PubMed]
96. Manna F, Martin G, Lenormand T (2011) Fitness landscapes: an alternative theory for the dominance of mutation. Genetics 189: 923–937. [PubMed]
97. Sellis D, Callahan BJ, Petrov DA, Messer PW (2011) Heterozygote advantage as a natural consequence of adaptation in diploids. Proc Natl Acad Sci USA 108: 20666–20671. [PubMed]
98. Lande R (1980) Sexual dimorphism, sexual selection, and adaptation in polygenic characters. Evolution 34: 292–305.
99. Rice WR, Chippindale AK (2001) Intersexual ontogenetic conict. J Evol Biol 14: 685–693.
100. Larsson J, Chen JD, Rasheva V, Rasmuson-Lestander As, Pirrotta V (2001) Painting of fourth, a chromosome-specific protein in Drosophila. Proc Natl Acad Sci USA 98: 6273–6278. [PubMed]
101. Johansson AM, Stenberg P, Bernhardsson C, Larsson J (2007) Painting of fourth and chromosomewide regulation of the 4th chromosome in Drosophila melanogaster. EMBO J 26: 2307–2316. [PubMed]
102. Johansson AM, Stenberg P, Pettersson F, Larsson J (2007) POF and HP1 bind expressed exons, suggesting a balancing mechanism for gene regulation. PLoS Genet 3: e209 doi:10.1371/journal.pgen.0030209. [PubMed]
103. Johansson AM, Stenberg P, Allgardsson A, Larsson J (2012) POF regulates the expression of genes on the fourth chromosome in Drosophila melanogaster by binding to nascent RNA. Mol Cell Biol 32: 2121–2134. [PMC free article] [PubMed]
104. Larsson J, Svensson MJ, Stenberg P, Mäkitalo M (2004) Painting of the fourth in genus Drosophila suggests autosome-specific gene regulation. Proc Natl Acad Sci USA 101: 9728–9733. [PubMed]
105. Vicoso B, Charlesworth B (2006) Evolution on the X chromosome: unusual patterns and processes. Nat Rev Genet 7: 645–653. [PubMed]
106. Lercher MJ, Urrutia AO, Hurst LD (2003) Evidence that the human X chromosome is enriched for male-specific but not female-specific genes. Mol Biol Evol 20: 1113–1116. [PubMed]
107. Kirkpatrick M, Hall DW (2004) Male-biased mutation, sex linkage, and the rate of adaptive evolution. Evolution 58: 437–440. [PubMed]
108. Xu K, Oh S, Park T, Presgraves DC, Yi SV (In press) Lineage-speci_c variation in slow- and fast-X evolution in primates. Evolution doi:10.1111/j.1558-5646.2011.01556.x. [PubMed]
109. Ohta T (1973) Slightly deleterious mutant substitutions in evolution. Nature 246: 96–98. [PubMed]
110. Smyth GK (2004) Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3: 3. [PubMed]
111. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy, Statist Soc Ser B 57: 289–300.
112. Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754–1760. [PMC free article] [PubMed]
113. R Development Core Team (2011) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.
114. Yanai I, Benjamin H, Shmoish M, Chalifa-Caspi V, Shklar M, et al. (2005) Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification. Bioinformatics 21: 650–659. [PubMed]
115. Fitch WM, Margoliash E (1967) Construction of phylogenetic trees. Science 155: 279–284. [PubMed]
116. Felsenstein J (2005) PHYLIP (Phylogeny Inference Package) version 3.6. Department of Genome Sciences, University of Washington, Seattle.
117. Schäffer J, Strimmer K (2005) A shrinkage approach to large-scale covariance matrix estimation and implications for functional genomics. Statist Appl Genet Mol Biol 4: 32. [PubMed]
Articles from PLoS Genetics are provided here courtesy of
Public Library of Science