The unique genomic properties of sex-biased genes: Insights from avian microarray data
© Mank et al; licensee BioMed Central Ltd. 2008
Received: 03 December 2007
Accepted: 31 March 2008
Published: 31 March 2008
In order to develop a framework for the analysis of sex-biased genes, we present a characterization of microarray data comparing male and female gene expression in 18 day chicken embryos for brain, gonad, and heart tissue.
From the 15982 significantly expressed coding regions that have been assigned to either the autosomes or the Z chromosome (12979 in brain, 13301 in gonad, and 12372 in heart), roughly 18% were significantly sex-biased in any one tissue, though only 4 gene targets were biased in all tissues. The gonad was the most sex-biased tissue, followed by the brain. Sex-biased autosomal genes tended to be expressed at lower levels and in fewer tissues than unbiased gene targets, and autosomal somatic sex-biased genes had more expression noise than similar unbiased genes. Sex-biased genes linked to the Z-chromosome showed reduced expression in females, but not in males, when compared to unbiased Z-linked genes, and sex-biased Z-linked genes were also expressed in fewer tissues than unbiased Z coding regions. Third position GC content, and codon usage bias showed some sex-biased effects, primarily for autosomal genes expressed in the gonad. Finally, there were several over-represented Gene Ontology terms in the sex-biased gene sets.
On the whole, this analysis suggests that sex-biased genes have unique genomic and organismal properties that delineate them from genes that are expressed equally in males and females.
Many genes are more actively transcribed in one sex than the other, and this sex-biased expression pattern is a mechanism by which heritable sexual dimorphisms can arise from a genome that is largely identical in males and females . Sex-biased gene expression is relatively common in metazoans [2–5], and has important evolutionary [1, 6–9], medical [10–12], and genomic [13, 14] implications.
Additionally, sex-biased genes are by their very nature often linked to reproduction, and this makes many sex-biased genes subject not only to natural selection, but the powerful pressures of sexual selection as well [15, 16]. The reproductive role of sex-biased gene expression is most easily observed in the global transcription profiles of the gonad, which shows the highest degree of sex-biased gene expression of all organs [17, 18]. In analysis where the entire animal was used for microarray analysis [2, 7, 19, 20], the gonad is typically the highest contributor to sex-biased gene lists. However, somatic tissues can exhibit remarkably high sex-biased expression patterns as well [5, 17], and these may produce secondary sexual characteristics and behaviors, or result from metabolic differences between males and females.
Microarray technologies have recently made it possible to study expression patterns for entire annotated transcriptomes, which has the potential to vastly increase our understanding of the underlying genomic mechanisms that give rise to recognizably dimorphic sexes. Initial evolutionary studies of transcriptome-wide sex-biased gene expression in Drosophila melanogaster suggested that sex-biased genes have peculiar evolutionary properties. Specifically, male-sex-biased genes, particularly those associated with gonad, evolve at a faster functional rate [9, 21–23], have higher rates of recombination that are tied to higher GC levels [14, 24], and are non-randomly distributed among the fly chromosomes [4, 25].
Recent results have suggested that the syndrome of sex-biased genes observed in D. melanogaster may not extend to all metazoans. Work in other Drosophila species did not recover an elevated rate of evolution for male-biased genes , and recent work in birds indicated that female-biased genes exhibit the highest rate of functional change . Birds are particularly of interest due to their female heterogametic (ZZ-ZW) system of sex chromosome inheritance, which allows for the analytical partitioning of maleness from heterogamety. Genomes with female heterogamety are subjected to different evolutionary forces than male heterogametic genomes , which can make birds a revealing contrast to male heterogametic mammalian and Drosophila lineages.
Elucidating underlying genomic and molecular differences between sex-biased and unbiased genes can give insight into the functionality and constraints for these expression classes, and ultimately help us understand the forces that produce their unique evolutionary properties. Global analysis of sex-biased gene expression in female heterogametic systems have so far been limited to birds [17, 28] and frogs , and a detailed understanding of the genomic properties of sex-biased genes in a female heterogametic system is warrented in order to understand how this sex chromosome system can effect sex-specific evolutionary pathways.
Here we study several expression (tissue specificity, expression level variance, gene function) and genomic (chromosome distribution, GC content, codon bias) properties of gene expressed in somatic (brain and heart) and gonad tissue of embryonic chickens, and investigate how these properties relate to sex-biased genes.
Distribution of microarray targets showing significant (p < 0.05) expression differences between male and female 18 day embryos in brain, heart, and gonad.
Significantly Expressed Targets
Absolute log2 Fold Change
Total Biased Targets
Clustering of genes with similar patterns of expression
Genes with similar breadth and level of expression have been shown to cluster in eukaryotic genomes [33–35], probably due to passive co-regulation of neighboring genes. We tested this possibility by calculating the average difference between both normalized level of expression between neighboring autosomal genes for which we had expression data, taking data from each tissue separately. We calculated the significance of local similarities by comparing these statistics with those derived from 10000 replicates where gene position was randomized on each chromosome. Neighboring genes showed significantly correlated levels of expression in gonad (p = 0.0002) and heart (p = 0.0114) although not in brain (p = 0.1848).
We also looked for clustering of sex-biased genes. Fold-change in expression, based only upon difference in expression between the sexes without a statistical cut-off p-value to define sex-bias, was significantly correlated in gonad (p < 0.0001), heart (p = 0.0175) and brain (p = 0.0131). All these correlations remained significant at the same level upon removal of comparisons of neighboring paralogs. This suggests that tandem duplications are not responsible for the observed local similarities, and that genes with similar expression patterns for some tissues are more likely than chance to be proximately located.
We also examined the genome for neighboring genes that both had significant evidence for biased expression in the same sex using the absolute value of log2 fold change > 1 with p adj < 0.05 to denote sex-bias. There were 157 such pairs using expression in gonad, compared with 100.9 expected by chance. When tandem duplicates were excluded, we observed 144 neighbor pairs, compared with 97.8 expected by chance alone. Both these associations are statistically significant (p < 0.0001 in both cases). However, when brain or heart expression was considered, there were no such examples of neighboring sex-biased genes.
GC3 and codon bias
Over-represented gene ontology terms for sex-biased categories (log2 fold-change. < 1,-> 1) within heart-expressed genes. Significantly (p < 0.01) over-represented terms are shown in term-for-term comparisons. Subsequent adjustments with the Bonferroni correction for multiple comparisons are shown.
Intrinsic to Golgi membrane
Amine metabolic process
Amino acid metabolic process
Amino acid and derivative process
Organic acid metabolic process
Carboxylic acid metabolic process
Nitrogen compound metabolic process
Over-represented gene ontology terms for sex-biased categories (log2 fold-change < 1,-> 1) within brain-expressed genes. Significantly (p < 0.01) over-represented terms are shown in term-for-term comparisons. Subsequent adjustments with the Bonferroni correction for multiple comparisons are shown.
Homocysteine S- methyltransferase activity
Voltage-gated calcium channel complex
Base excision repair
Fructose bisphosphate aldolase activity
External encapsulating structure
External encapsulating structure part
Cell outer membrane
Over-represented gene ontology terms for sex-biased categories (log2 fold-change < 1,-> 1) within gonad-expressed genes. Significantly (p < 0.01) over-represented terms are shown in term-for-term comparisons. Subsequent adjustments with the Bonferroni correction for multiple comparisons are shown.
Calcium ion binding
Cytoskeletal-dependent intracellular transport
Iron ion transport
Enzyme regulator activity
Voltage-gated ion channel activity
Voltage-gated potassium channel activity
Iron-sulfur cluster binding
Metal cluster binding
Protein-tyrosine kinase activity
Electron carrier activity
Integral to membrane
Intrinsic to membrane
Galactosylgalactosylxylosylprotein 3-beta glucuronosyltransferase activity
Cellular metabolic process
Primary metabolic process
Acid-amino acid ligase activity
Small conjugating protein ligase activity
Ligase activity, forming carbon-nitrogen bonds
Macromolecule metabolic process
Protein metabolic process
Cellular catabolic process
Nucleic acid binding
Biopolymer metabolic process
Cellular protein metabolic process
Intracellular membrane-bound organelle
Cellular macromolecule metabolic process
The data in this study was collected from both somatic and gonad tissue of embryonic day 18 chicks, after the circulating testosterone and estrogen levels diverge between males and females and initiate sexual differentiation . This suggests that the sex-biased genes identified in this study are involved in development of sexually dimorphic anatomy, physiology, and behavior. However, gametogenesis does not commence until well after hatching in chicken, therefore these experiments are more focused on the development of the gonad than on genes specific to gametogenesis. It is therefore difficult to compare the degree of sex-bias observed here to studies of other organisms taken at different stages of the life-cycle [5, 23, 29]. Sex-biased gene expression will increase with the onset of sexual maturity, due both to the expression of gamete-specific transcripts and genes with sex hormone-receptors. However, in general it would be expected that juvenile individuals would show lower levels of sex-biased gene expression than adults in all tissues, and this is indeed what we observe. The degree of sex-bias for autosomal genes is less in our juvenile gonad tissue when compared both to Xenopus laevis  and Drosophila melanogaster  adult tissue at similar fold-change cut-off levels. Chicken embryonic brain tissue also showed less degree of sex-bias for autosomal genes than adult mouse brains  at similar fold-change levels.
Sex-biased genes were present in both somatic and gonad tissues, though the gonad showed the highest proportion of sex-biased gene expression across the range of fold-change cutoff values. This likely stems from the fact that the gonad is the most sexually dimorphic organ, and although the testes and ovary share a common precursor, in many ways they represent two distinct tissue sets. It is therefore logical that sex-biased gene expression would be the most pronounced in this region.
Sex-biased genes share several expression properties that have implications to their evolution. First, sex-biased genes were more tissue specific in our analysis (Fig. 1), which is consistent with previous work suggesting sex-biased genes are more narrowly expressed than unbiased genes [5, 39]. Coding regions that are broadly expressed face many, sometimes contradictory, evolutionary constraints on gene expression evolution , and it may be that genes with extensive expression throughout an organism are less able to evolve sex-biased expression patterns without disrupting critical functions in other tissues.
In addition to being expressed at lower levels (Fig. 2), sex-biased genes also showed more variance in expression (Fig. 3), at least for genes mapped to autosomal locations. Variance is an indicator of expression noise, and genes with larger amounts of expression noise have been shown to be less important in organismal survival and functionality [30–32]. Somatic organs are functionally identical in males and females, therefore sex-biased autosomal genes expressed in these tissues theoretically do not have important developmental or regulatory roles. The reduced difference between autosomal biased and unbiased genes in the gonad may reflect the sex-specific nature of the ovary and testes, where the pattern of sex-biased expression observed in the gonad is required for the formation of these distinct tissues. The lack of dosage compensation of the Z chromosome means that most genes have a default male-biased expression pattern, and this confounds the analysis of male-biased Z-genes and essentiality. The reduced variance observed across all tissues for male-biased Z-linked genes when compared to male-biased autosomal genes likely stems from the default male-biased expression pattern, as male-biased genes may have more important roles in the organism, but still exhibit male-biased expression simply due to Z-linkage. Overall, these findings suggest that somatic sex-biased genes, particularly those attached to the autosomes, may have less critical roles for organismal survival than unbiased genes.
There are indications  that avian sex-biased genes show an altered form of the accelerated rates of protein evolution seen in other animals[4, 9], with female biased genes in birds having the highest rate of functional change. The above-described connections between tissue specificity, expression level, and expression variance and sex-bias may offer explanations for this rapid rate of evolution. Rates of evolution have been positively correlated with tissue specificity [41, 42], and negatively correlated with overall expression level [43, 44]. Additionally, increased variance in expression level may indicate that a gene is not critical , and therefore subject to less evolutionary constraint, which could theoretically manifest in an accelerated rate of protein change for genes with high expression variance. This would mean that genes which are tissue specific, have low expression levels, as well as high expression variances, are also likely to exhibit accelerated rates of protein evolution. Sex-biased genes possess all the characteristics, and it seems therefore likely that these characteristics play an important role in the rapid rate of evolution for sex-biased genes.
We found significant clustering of similarly expressed genes. Neighboring genes in both the heart and gonad showed significant correlations in sex-averaged expression level. Additionally, when assessed just by fold-change, sex-biased genes in all three tissues showed significant clustering, though only gonad showed significant non-random distribution after adjusted-p-values were used to denote significantly sex-biased genes. All the observed clustering patterns, for both expression level and for sex-biased expression may best be explained by local effects such as bidirectional promoters and transcriptional read-through, which cause the transcription of genes to influence expression of their neighbors .
GC3 and codon bias
Work in D. melanogaster has demonstrated different levels of codon-bias among sex-biased expression classes , and we recovered a similar pattern in chicken. Specifically, we observed codon-bias in female-biased genes expressed in the gonad. This could be explained be due to selection for specific codon usage patterns, or to differing GC3 levels that we observed across sex-biased categories (Fig. 4), as GC3 patterns can influence codon bias . When we compared the codon-bias values of each sex bias category with those predicted by GC3 from the model fitting, we found no significant differences (one-way ANOVA on residuals), suggesting that there no evidence for any differences in selection for codon-usage between sex-biased and unbiased genes. Our observations therefore suggest that in chicken, as in mammals, codon-usage bias is mainly determined by the GC content of the isochore where each gene is situated, rather than natural selection acting on some aspect of gene expression .
We found some over-representation of Gene Ontology terms in our sex-biased expression categories, particularly for those sex-biased genes expressed in the gonad. It is of course tempting to speculate about the possible reasons that over-represented GO terms are associated with sex-bias expression classes. While significantly over-represented GO terms from the male-and female-biased lists lack an immediately obvious sex-specific functionality, the unbiased gonad list is striking in the number of basal metabolic terminologies. Housekeeping and basal-functioning genes would not be expected to exhibit sex-biased expression, partly because they are important to both sexes and partly because they are generally broadly expressed and therefore less likely to exhibit sex-biased expression, as described above. This may explain the over-abundance of such terms as cellular metabolic process, mitochondrion, and ligase activity, which are theoretically common to all organs and tissues.
We have presented a characterization of several properties of sex-biased genes identified from chicken microarray experiments. Sex-biased genes in chicken have a larger variance in expression level, are expressed in fewer tissues, and at lower overall levels than unbiased genes. Additionally, sex-biased genes have somewhat different GC3, and codon-bias properties. These qualities may help explain why sex-biased genes have unique evolutionary properties, and present a useful framework for future analyses of sex-biased expression patterns.
Sample collection and preparation
White Leghorn embryos were euthanized after 18 days of incubation (ed18), and the brain, heart, and left gonad collected from four male and four female individuals. This data was originally collected for an analysis of dosage compensation, and further details about study design, tissue preparation, RNA isolation, and microarray hybridization and analysis have been previously published . We note only major experimental details and where the analysis presented here differs. All microarray data is available at the Gene Expression Omnibus (project GSE 8693).
Annotation of probe sets
Annotations and genomic locations for the microarray probe sets were extracted from Ensembl  using version 2.1 of the chicken genome (WASHUC2 May 2006). All genes that were not assigned to a specific chromosome in this build were removed from further analysis, as many of these genes may be sex-linked and show dosage effects in expression. Additionally, as the W chromosome assembly is incomplete, the few W-assigned gene targets were also removed. The remaining dataset was comprised of coding regions that have been definitively mapped to either the autosomes or the Z chromosome. This dataset was then parsed into autosomal and Z components for all subsequent analyses, as the lack of effective and complete dosage compensation on the chicken Z chromosome [17, 28] results in male-biased expression for Z-linked genes due to simple dosage effects, which confounds subsequent studies of expression bias.
Identification of sex-biased genes
For each actively transcribed target, the log2 (average male/average female) fold-change value was calculated. This treatment results in negative values for genes that are female-biased, positive values for genes that are male-biased, and values approaching 0 for genes with similar expression in both sexes. Genes were further divided into fold-change cutoff categories, based on the Wilcoxon test statistic of the log2 values, and adjusted for a 5% false discovery rate using the Benjamini Hochberg correction. Therefore, absolute value of log2 fold change > 1 corresponds to a 2× difference in expression in one sex compared to the other, absolute value of log2 fold change > 1.5 represents roughly a 3× increase, and absolute value of log2 fold change > 2 corresponds to a 4× difference. Relative expression was calculated across within-sex replicates for sex-biased (p adj < 0.05, absolute fold-change >1) and unbiased gene classes for each tissue separately, and 95% confidence intervals were determined via bootstrapping (1000 repetitions).
In order to investigate the relationship between expression variance, an indicator of noise, and fold-change, we normalized the expression level for each gene target such that the averaged expression level was equal to one within each set of sex-specific replicates. This reduces the problem where genes with low expression levels will have associated small variance estimates, which could create a false signal. We then computed average within-sex variance for each sex-bias expression class (male biased log2 fold change > 1, p adj > 0.05; female biased genes where log2 fold change < 1- and p adj < 0.05; and unbiased probe sets), and computed 95% confidence intervals via bootstrapping (1000 replicates).
Clustering of genes with similar patterns of expression
For all genes that were significantly expressed in our study, we calculated the difference in average level of expression between it and its nearest neighbor for which gene expression had been measured in the same tissue. We used sex-averaged expression to correct for sex-bias effects. We also measured the difference in fold-change between neighboring gene pairs for which we had expression data, and counted the number of times gene pairs both had evidence for significant biased expression in the same sex (absolute log2 fold change > 1, p < 0.05). We examined the significance of these values by comparing them with datasets where gene order was randomized on each chromosome using 10000 replicates.
GC3 and bias in codon usage
Third codon GC (GC3) content has been linked to codon-bias [36, 47], which has in turn been linked to sex-biased gene expression patterns in Drosophila . Therefore, in order to investigate the relationship between sex-biased gene expression and codon bias, we calculated GC3 and the Effective Number of Codons (ENC)  from the Ensembl annotation of the chicken genome. ENC can range from 20, in the case of extreme bias where only one codon is used for each amino acid, to 61, where all alternative codons are equally likely. We calculated the average GC3 and ENC for each sex-bias expression category within each tissue at the fold-change < 1- or > 1, p adj < 0.05 cutoff. 95% confidence intervals were computed with bootstrapping (1000 replicates). Within each tissue, we compared the GC3 and ENC among autosomal expression categories using one-way ANOVA (2 d.f.), and between Z-linked categories using a one-tailed t-test assuming unequal variance. As ENC is strongly governed by GC content in vertebrates, we also tested whether differences in ENC values between categories were significantly different from that predicted by the observed differences in GC3. To do this we fitted a smoothed cubic spline to the relationship between ENC and GC3, and used this to estimate predicted ENC values from GC3. Significance was tested performing a one-way ANOVA on the residuals of model fitting.
Gene Ontology analysis
We compared the functional role of biased and unbiased gene sets for each tissue using Ontologizer  and the Gallus Gene Ontology (GO)  database to identify over-represented GO categories. For each tissue, the actively expressed genes were used as the population gene set in term-for-term comparisons with each of the sex-biased expression categories as separate study sets (absolute fold change > 1, p adj < 0.05 cutoffs to denote sex-bias). Significant over-representation was determined for each term with Fisher's exact test, and the Bonferroni correction was used for multiple comparisons. As the Bonferroni method can be unnecessarily over-conservative, we report over-represented terms that are significant before Bonferroni correction as well.
We thank the Swedish Research Council and the Wenner-Gren Foundation for support for this work.
- Ellegren H, Parsch J: The evolution of sex-biased genes and sex-biased gene expression. Nature Reviews Genetics. 2007, 8: 689-698. 10.1038/nrg2167.PubMedView Article
- Eads BD, Colbourne JK, Bohuski E, Andrews J: Profiling sex-biased gene expression during parthenogenetic reproduction in Daphnia pulex. BMC Genomics. 2007,
- Marinotti O, Calvo E, Nguyen QK, Dissanayake S, Ribeiro JMC, James AA: Genome-wide analysis of gene expression in adult Anopheles gambiae. Insect Molecular Biology. 2006, 15 (1): 1-12. 10.1111/j.1365-2583.2006.00610.x.PubMedView Article
- Ranz JM, Castillo-Davis CI, Meiklejohn CD, Hartl DL: Sex-dependent gene expression and evolution of the Drosophila transcriptome. Science. 2003, 300 (5626): 1742-1745. 10.1126/science.1085881.PubMedView Article
- Yang X, Schadt EE, Wang S, Wang H, Arnold AP, Ingram-Drake L, Drake TA, Lusis AJ: Tissue-specific expression and regulation of sexually dimorphic genes in mice. Genome Research. 2006, 16 (8): 995-1004. 10.1101/gr.5217506.PubMedPubMed CentralView Article
- Mank JE, Hultin-Rosenberg L, Axelsson E, Ellegren H: Rapid evolution of female-biased, but not male-biased, genes expressed in avian brain. Molecular Biology and Evolution. 2007, 24: 2698-2706. 10.1093/molbev/msm208.PubMedView Article
- Meiklejohn CD, Parsch J, Ranz JM, Hartl DL: Rapid evolution of male-biased gene expression in Drosophila. Proceedings of the National Academy of Sciences, USA. 2003, 100 (17): 9894-9899. 10.1073/pnas.1630690100.View Article
- Pröschel M, Zhang Z, Parsch J: Widespread adaptive evolution of Drosophila genes with sex-biased expression. Genetics. 2005, 174: 893-900. 10.1534/genetics.106.058008.View Article
- Zhang Z, Hambuch TM, Parsch J: Molecular evolution of sex-biased genes in Drosophila. Molecular Biology and Evolution. 2004, 21 (11): 2130-2139. 10.1093/molbev/msh223.PubMedView Article
- Ivakine EA, Fox CJ, Pateron AD, Mortin-Toth SM, Canty A, Walton DS, Aleksa K, Ito S, FDanska JS: Sex-specific effect of insulin-dependent diabetes 4 on regulation of diabetes pathogenesis in the nonobese diabetic mouse. Journal of Immunology. 2005, 174: 7129-7140.View Article
- Naugler WE, Sakurai T, Kim S, SMaeda S, Kim K, Elsharkawy AM, LKarin M: Gender disparity in liver cancer due to sex differences in MyD88-dependent IL-6 production. Science. 2007, 317: 121-124. 10.1126/science.1140485.PubMedView Article
- Vandenbroeck K, Cunningham S, Goris A, Alloza I, Heggarty S, Graham C, Bell A, Rooney M: Polymorphisms in the interferon-gamma/interleukin-26 gene region contribute to bias in susceptibility to rheumatoid arthritis. Arthritis and Rheumatism. 2003, 48: 2773-2778. 10.1002/art.11236.PubMedView Article
- Connallon T, Knowles LL: Intergenomic conflict revealed by patterns of sex-biased gene expression. Trends in Genetics. 2005, 21: 495-499. 10.1016/j.tig.2005.07.006.PubMedView Article
- Hambuch TM, Parsch J: Patterns of synonymous codon usage in Drosophila melanogaster genes with sex-biased expression. Genetics. 2005, 170 (4): 1691-1700. 10.1534/genetics.104.038109.PubMedPubMed CentralView Article
- Civetta A, Singh RS: Broad-sense sexual selection, sex gene pool evolution, and speciation. Genome. 1999, 42 (6): 1033-1041. 10.1139/gen-42-6-1033.PubMedView Article
- Swanson WJ, Vacquier VD: The rapid evolution of reproductive proteins. Nature Reviews Genetics. 2002, 3: 137-144. 10.1038/nrg733.PubMedView Article
- Ellegren H, Hultin-Rosenberg L, Brunström B, Dencker L, Kultima K, Scholtz B: Faced with inequality: Chicken does not have general dosage compensation of sex-linked genes. BMC Biology. 2007
- Parisi M, Nuttall R, Edwards P, Minor J, Naiman D, Lu JN, Doctolero M, Vainer M, Chan C, Malley J: A survey of ovary-, testis-, and soma-biased gene expression in Drosophila melanogaster adults. Genome Biology. 2004, 5: R40-10.1186/gb-2004-5-6-r40.PubMedPubMed CentralView Article
- Cutter AD, Ward S: Sexual and temporal dynamics of molecular evolution in C. elegans development. Molecular Biology and Evolution. 2005, 22 (1): 178-188. 10.1093/molbev/msh267.PubMedView Article
- Reinke V, Gil IS, Ward S, Kazmer K: Genome-wide germline-enriched and sex-biased expression profiles in Caenorhabditis elegans. Development. 2004, 131 (2): 311-323. 10.1242/dev.00914.PubMedView Article
- Haerty W, Jagadeeshan S, LKulathinal RJ, Wong AR, K R, Sirot LK, Levesque L, Artieri CG, Wolfner MF, Civetta A, Singh RS: Evolution in the fast lane: Rapidly evolving sex-related genes in Drosophila. Genetics. 2007, 177: 1321-1335. 10.1534/genetics.107.078865.PubMedPubMed CentralView Article
- Jagadeeshan S, Singh RS: Rapidly evolving genes of Drosophila: Differing levels of selective pressure in testis, ovary, and head tissues between sibling species. Molecular Biology and Evolution. 2005, 22 (9): 1793-1801. 10.1093/molbev/msi175.PubMedView Article
- Zhang Y, Sturgill D, Parisi M, Kumar S, Oliver B: Constraint and turnover in sex-biased gene expression in the genus Drosophila. Nature. 2007, 450: 233-238. 10.1038/nature06323.PubMedPubMed CentralView Article
- Zhang Z, Parsch J: Positive correlation between evolutionary rate and recombination rate in Drosophila genes with male-biased expression. Molecular Biology and Evolution. 2005, 22 (10): 1945-1947. 10.1093/molbev/msi189.PubMedView Article
- Parisi M, Nuttall R, Naiman D, Bouffard G, Malley J, Andrews J, Eastman S, Oliver B: Paucity of genes on the Drosophila X chromosome showing male-biased expression. Science. 2003, 299 (5607): 697-700. 10.1126/science.1079190.PubMedPubMed CentralView Article
- Metta M, Gudavalli R, Gibert JM, Schlotterer C: No accelerated rate of protein evolution in male-biased Drosophila pseudoobscura genes. Genetics. 2006, 174 (1): 411-420. 10.1534/genetics.106.057414.PubMedPubMed CentralView Article
- Axelsson E, Smith NGC, Sundström H, Berlin S, Ellegren H: Male-biased mutation rate and divergence in autosomal, Z-linked and W-linked introns of chicken and turkey. Molecular Biology and Evolution. 2004, 21 (8): 1538-1547. 10.1093/molbev/msh157.PubMedView Article
- Itoh Y, Melamed E, Yang X, Ka pK, Wang S, Yehya N, van Nas A, Replogle K, Band MR, Clayton DF: Dosage compensation is less effective in birds than in mammals. Journal of Biology. 2007, 6: 2-10.1186/jbiol53.PubMedPubMed CentralView Article
- Malone JH, Hawkins DL, Michalak P: Sex-biased gene expression in a ZW sex determination system. Journal of Molecular Evolution. 2006, 63 (4): 427-436. 10.1007/s00239-005-0263-4.PubMedView Article
- Batada NN, Hurst LD: Evolution of chromosome organization driven by selection for reduced gene expression noise. Nature Genetics. 2007, 39: 945-949. 10.1038/ng2071.PubMedView Article
- Fraser HB, Hirsh AE, Giaever G, Kumm J, Seisen MB: Noise minimization in eukaryotic gene expression. PLOS Biology. 2004, 2: 0834-
- Newman JRS, Ghaemmaghami S, Ihmels J, Breslow DK, Noble M, DeRisi JL, Weissman JS: Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise. Nature. 2006, 441: 840-846. 10.1038/nature04785.PubMedView Article
- Cohen BA, Mitra RD, Hughes JD, Church GM: A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression. Nature Genetics. 2000, 26 (2): 183-186. 10.1038/79896.PubMedView Article
- Hurst LD, Pal C, Lercher MJ: The evolutionary dynamics of eukaryotic gene order. Nature Reviews Genetics. 2004, 5 (4): 299-310. 10.1038/nrg1319.PubMedView Article
- Lercher MJ, Hurst LD: Co-expressed yeast genes cluster over a long range but are not regularly spaced. Journal of Molecular Biology. 2006, 359 (3): 825-831. 10.1016/j.jmb.2006.03.051.PubMedView Article
- Donofrio G, Mouchiroud D, Aissani B, Gautier C, Bernardi G: Correlations between the compositional properties of human Genes, codon usage, and amino acid composition of proteins. Journal of Molecular Evolution. 1991, 32 (6): 504-510. 10.1007/BF02102652.View Article
- Wright F: The effective number of codons used in a gene. Gene. 1990, 87: 23-29. 10.1016/0378-1119(90)90491-9.PubMedView Article
- Bruggeman V, van As P, Decuypere E: Developmental endocrinology of the reproductive axis in the chicken embryo. Comparative Biochemistry and Physiology, A. 2002, 131: 839-846. 10.1016/S1095-6433(02)00022-3.View Article
- Mank JE, Hultin-Rosenberg L, Zwahlen M, Ellegren H: Pleiotropic constraint hampers the resolution of sexual antagonism in vertebrate gene expression. The American Naturalist. 2008, 171: 35-43. 10.1086/523954.PubMedView Article
- Yang J, Su AI, Li WH: Gene expression evolves faster in narrowly than in broadly expressed mammalian genes. Molecular Biology and Evolution. 2005, 22 (10): 2113-2118. 10.1093/molbev/msi206.PubMedView Article
- Liao BY, Scott NM, Zhang JZ: Impacts of gene essentiality, expression pattern, and gene compactness on the evolutionary rate of mammalian proteins. Molecular Biology and Evolution. 2006, 23 (11): 2072-2080. 10.1093/molbev/msl076.PubMedView Article
- Liao BY, Zhang JZ: Low rates of expression profile divergence in highly expressed genes and tissue-specific genes during mammalian evolution. Molecular Biology and Evolution. 2006, 23 (6): 1119-1128. 10.1093/molbev/msj119.PubMedView Article
- Pal C, Papp B, Hurst LD: Highly expressed genes in yeast evolve slowly. Genetics. 2001, 158 (2): 927-931.PubMedPubMed Central
- Subramanian S, Kumar S: Gene expression intensity shapes evolutionary rates of the proteins encoded by the vertebrate genome. Genetics. 2004, 168 (1): 373-381. 10.1534/genetics.104.028944.PubMedPubMed CentralView Article
- Semon M, Duret L: Evolutionary origin and maintenance of coexpressed gene clusters in mammals. Molecular Biology and Evolution. 2006, 23 (9): 1715-1723. 10.1093/molbev/msl034.PubMedView Article
- Chamary JV, Parmley JL, Hurst LD: Hearing silence: non-neutral evolution at synonymous sites in mammals. Nature Reviews Genetics. 2006, 7 (2): 98-108. 10.1038/nrg1770.PubMedView Article
- Collins DW, Jukes TH: Relationship between G+C in silent sites of codons and amino acid composition of human proteins. Journal of Molecular Evolution. 1993, 36 (3): 201-213. 10.1007/BF00160475.PubMedView Article
- Robinson PN, Wollstein A, Böhme U, Beattie B: Ontologizing gene-expression microarray data: characterizing clusters with Gene Ontology. Bioinformatics. 2004, 20: 979-981. 10.1093/bioinformatics/bth040.PubMedView Article
- Gene Ontology Consortium: Gene ontology: a tool for the unification of biology. Nature Genetics. 2000, 25: 25-29. 10.1038/75556.View Article
- Ensembl. [http://www.ensembl.org]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.