- Open Access
Male X-linked genes in Drosophila melanogaster are compensated independently of the Male-Specific Lethal complex
Epigenetics & Chromatin volume 6, Article number: 35 (2013)
In organisms where the two sexes have unequal numbers of X-chromosomes, the expression of X-linked genes needs to be balanced not only between the two sexes, but also between X and the autosomes. In Drosophila melanogaster, the Male-Specific Lethal (MSL) complex is believed to produce a 2-fold increase in expression of genes on the male X, thus restoring this balance.
Here we show that almost all the genes on the male X are effectively compensated. However, many genes are compensated without any significant recruitment of the MSL-complex. These genes are very weakly, if at all, affected by mutations or RNAi against MSL-complex components. In addition, even the genes that are strongly bound by MSL rely on mechanisms other than the MSL-complex for proper compensation. We find that long, non-ubiquitously expressed genes tend to rely less on the MSL-complex for their compensation and genes that in addition are far from High Affinity Sites tend to not bind the complex at all or very weakly.
We conclude that most of the compensation of X-linked genes is produced by an MSL-independent mechanism. Similar to the case of the MSL-mediated compensation we do not yet know the mechanism behind the MSL-independent compensation that appears to act preferentially on long genes. Even if we observe similarities, it remains to be seen if the mechanism is related to the buffering that is observed in autosomal aneuploidies.
In organisms where females have two X-chromosomes and males have only one, some type of dosage compensation is needed. The expression of X-linked genes needs to be balanced not only between the two sexes, but also between X and the autosomes, in order to maintain the balance of metabolic networks[1–3]. In Drosophila melanogaster, this is achieved by an approximately 2-fold increase in expression of the genes on the male X-chromosome. It is commonly stated in the literature that the Male-Specific Lethal (MSL) complex is responsible for this 2-fold increase by binding to expressed genes on the male X[4–8]. The MSL-complex consists of five protein components, namley the H4K16 histone acetyltransferase Male absent on the first (MOF), Maleless, Male-Specific Lethal 1 (MSL1), Male-Specific Lethal 2 (MSL2), and Male-Specific Lethal 3 (MSL3), and two non-coding RNAs (roX1 and roX2)[9, 10]. The components of the MSL-complex clearly have important functions since mutations are male-specific lethal. However, many (10–25%) of the expressed genes on the male X are not bound by the complex[11–14]. There have been different explanations for this; the most common explanation has been that these unbound genes are only transiently bound by the complex and thus binding escapes detection. Although it has generally been assumed that the genes that are unbound but still expressed are dosage compensated, this assumption has not been investigated thoroughly on a genome-wide level. Zhang et al. found that genes on the X-chromosome were, on average, dosage compensated approximately 1.35-fold by an MSL-dependent mechanism, irrespective of the gene dose, in a male cell-line. It has been speculated that the remaining compensation, up to two-fold, could be mediated by a more general buffering mechanism acting on all genomic regions that are present in single copy[16–19]. The idea that the MSL-complex is not responsible for the bulk of chromosome X dosage compensation has been previously proposed by Birchler et al.[20–22], suggesting that an “inverse dosage effect” causes an up-regulation of the single male X-chromosome mainly through the loss of negative regulators. However, to date, no study has attempted a genome-wide scale correlation between compensation and the binding of the MSL-complex at the individual gene level.
If the MSL-complex is the mediator of the full, 2-fold, dosage compensation, we would expect to observe an approximately 2-fold reduction in the expression of genes on the X-chromosome if components of the complex are removed. Arguing against this prediction, several studies have shown that expression of male X-linked genes is only reduced to about 80% of the normal value after RNAi knock-down or mutation of MSL-complex components[5, 6, 23, 24].
In this paper, we use publicly available data which we combine and present in novel ways in order to investigate the proportion of expressed genes on the male X that are dosage compensated. We also analyze how much the dosage compensated genes are dependent on the MSL-complex.
Results and discussion
Most expressed genes on the male X are dosage compensated
We first set out to investigate how many of the expressed genes on the male X-chromosome are dosage compensated. To do this, we used the (to our knowledge) only available set of genome wide expression data from a single tissue (salivary glands) where both males and females were studied separately. It is important to perform this analysis on a sample including as few different cell types as possible since we would otherwise be unable to distinguish gene expression differences between the sexes from, for example, variation in cell type proportions. Similarly, we cannot use cell-lines as male and female cell-lines could have very different origins. Since unexpressed genes or genes expressed below the detection limit of the technique (in this case array) will appear to be dosage compensated, we excluded all genes not expressed in males and/or females, using the same approach as in previous studies[18, 25, 26]. We next plotted the male/female expression ratio across the major chromosome arms (chromosome X and 3L are shown in Figure 1, the other major chromosome arms are shown in Additional file1: Figure S1) and calculated the ratio mean and standard deviation (the male/female ratio appears to be normally distributed on all chromosomes). Interestingly the chromosomes show a similar mean and standard deviation (mean -0.044 to -0.003, SD 0.33 to 0.51). The X-chromosome seems to be well compensated (mean -0.020) and shows low variation in the expression ratio (SD 0.38).
We then defined dosage compensated genes as genes on chromosome X that are expressed in both males and females and where the log2 ratio was within one standard deviation of the mean for an autosomal chromosome; for this study we used chromosome 3L (> - 0.55, <0.47). According to this definition, 90.4% of all expressed genes (n = 882) on the X-chromosome are dosage compensated (Figure 1), a value which is comparable to those of the autosomes where 88.2% to 91.1% of the genes are expressed at the same levels in males and females.
Based on these results, we conclude that most genes on the male X-chromosome are effectively compensated and that there are similar numbers of differentially expressed genes when comparing males and females across all major chromosome arms. Although this effect has not been formally demonstrated previously, it is in line with current opinion in the field. Expression levels in males and females have been compared in other studies[24, 27]. However, these comparisons were done using whole flies, making it impossible to separate dosage compensatory effects acting on individual genes from differences in expression between tissues. Moreover, such studies (or studies of whole embryos) do not allow direct comparison of expression effects at an individual gene level with binding levels of the MSL-complex.
About 15% of expressed genes on male X are dosage compensated but unbound by MSL and have significantly reduced levels of H4K16 acetylation
Various figures have been reported for the number of genes that are expressed but not bound by MSL[11–14]. To determine how many of the compensated genes are significantly bound by the MSL-complex we used available MOF binding data from male salivary glands. MOF binds to both the promoters and the gene bodies of MSL-targeted genes on the X-chromosome[26, 29, 30]. After determining the extent of binding to all gene bodies excluding the first 200 bp of exon sequence (see Methods), we plotted the male/female expression ratio versus MOF binding for all chromosomes (chromosome X and 3L are shown in Figure 2A,B).
To determine which genes undergo binding by MOF to the gene body, we calculated the mean and standard deviation of MOF binding to expressed genes on chromosome 3L (MOF binding values were found to be normally distributed on 3L). Genes on X with a binding value higher than the 3L mean plus one standard deviation (mean = 1.00, SD = 0.48) were considered to be bound. This cut-off probably identifies some unbound or very weakly bound genes as bound, since 147 genes on chromosome 3L would be considered bound by this definition. It is known that MOF is also bound to autosomes as part of the Non-Specific Lethal complex, but since this binding is to the promoters of genes it will not influence our calculated gene binding values for MOF (see Methods). Using this definition, we find that 122 genes on the male X-chromosome are dosage compensated although not bound at any detectable level by the MSL-complex. The fact that many expressed genes are not bound by the MSL-complex has been reported several times. However, the level to which these genes are compensated has been unknown[29, 32, 33]. The usual explanation for the lack of MSL binding to these expressed genes has been that they only transiently interact with the complex and they do show H4K16 acetylation marks. However, when we look at the H4K16ac data from male salivary glands in the gene bodies of all dosage compensated X-linked genes, we see drastically reduced levels of acetylation of the genes that we classified as unbound but still compensated (Figure 2C). Although we cannot exclude the possibility of the MSL-complex influencing these genes on the basis of these results, we argue that it is very unlikely.
The data we are using is from salivary glands which have multiple copies of each chromosome. It is known that parts of the genome are under-replicated in salivary glands and the female X-chromosome have more such regions when compared to the male X-chromosome[34, 35]. It is therefore possible that MSL unbound genes in our data set will appear to be compensated, when in fact they may have similar copy numbers in males and females. To investigate this, we first counted how many expressed and compensated genes map to the under-replicated regions on the X-chromosome reported by Makunin et al.; 27% of the MSL unbound, expressed and compensated genes map to under-replicated regions, which is very similar to the 28% of strongly bound (MOF binding value >5), expressed and compensated genes that map to under-replicated regions. In addition we calculated the copy number ratio between males and females of all compensated genes using the input reads of the MOF ChIP-seq experiments. The male and female input samples varied in sequencing depth but we observed that on average the X-linked genes clearly have a lower male to female ratio when compared to the autosomes (Man-Whitney U Test, P <0.05). The unbound compensated genes did not have more similar copy number in males versus females. In this analysis we saw a weak trend that the more MOF bound the higher the male/female copy number ratio. We therefore conclude that the compensation of the MSL unbound genes is not caused by similar copy numbers in males and females.
Genes unbound by MSL but still compensated are not affected by MSL RNAi
In order to determine whether the genes identified as unbound (or bound below the detection limit) by MSL are in fact independent of the MSL-complex, we turned to MSL-complex component RNAi knock-down and mutant data available from S2 cells and male larvae, respectively. Although these data sets are not from salivary glands, we and others have shown that MSL-complex binding is extremely similar between different cell types[13, 26, 32] (Additional file2: Figure S2). Although the MSL binding levels of some genes are likely to differ between, for example, S2 cells and salivary glands, we believe that comparison of data from these sources will be meaningful (we exclude genes that are not expressed in the wild-type or mock RNAi).
When we correlate the effect of RNAi knock-down or mutations in MSL-complex components with MOF binding to compensated genes, we see that the expression of genes that are unbound but still compensated is very weakly affected by MSL (Figure 3; Additional file3: Figure S3). In the mutant data sets from male larvae, the unbound genes are virtually unaffected, and in the RNAi experiments these genes show a weak reduction in expression (significantly different from autosomes, four RNAi data sets, Man-Whitney U Test, all P <0.05). The small effect on the unbound genes in the RNAi experiments could indicate that they are weakly bound (below what we can detect with ChIP-seq) or that these effects are secondary gene network effects. We note that the most strongly bound genes show less reduction in expression in some experiments (Additional file3: Figure S3). In the RNAi knock-down experiments this could be due to the fact that knock-down is not complete and that weakly bound genes will lose the MSL-complex before those that are more strongly bound. If so, that would make it unlikely that the genes we defined as unbound have weak interactions with the MSL-complex as those should be lost first. Unbound genes also score lower than bound genes (Man-Whitney U Test, all P <0.05) on sequence signatures previously found to be enriched in MSL bound genes in S2 cells, further supporting that unbound expressed genes do not recruit MSL. We note that the genes which show the strongest reduction in expression are still expressed at about 80% of wild-type levels. We conclude that a significant proportion of dosage compensated genes on the male X-chromosome is compensated by a mechanism independent of MSL. A few such genes have been identified in earlier studies[21, 22, 37–41], but in salivary glands, these genes (such as runt, giant and buttonhead) are either not expressed according to our definitions, or, if expressed, bound by MOF.
MSL independent mechanisms probably play an important role in compensating the male X-chromosome
An MSL-independent mechanism probably also plays a significant role for MSL-bound genes, as the latter are only mildly affected by the loss of MSL-complex components. This was previously observed by (for example) Zhang et al., and could be explained either by the presence of a general buffering mechanism, which recognizes genomic regions present in fewer copies and which acts on all chromosomes[18, 25], or by gene network effects. The general buffering of autosomal deficiencies reported by Lundberg et al. showed stronger buffering of long genes; if the same mechanism also acts on the X-chromosome, we expect that long genes would be less MSL dependent. We found that unbound genes are on average much longer than bound genes (Figure 4A) suggesting that long genes are less dependent on the MSL-complex.
If this mechanism applies, we would expect long bound genes to be less affected by reduction or loss of the MSL-complex. When long genes (>10 kb, the approximate median length of unbound genes) are compared to short genes (<10 kb), MSL-bound long genes are less affected by RNAi/mutation of MSL components (five data sets, Man-Whitney U Test, all P <0.05), with the exception of the MSL-3 RNAi which was investigated by Kind et al.. The differences between long and short genes are most prominent in the msl2 and roX1 roX2 mutant datasets. For bound genes, there is a correlation between transcript length and expression in msl2 mutants (R = 0.48, P <0.05, Figure 4B). We conclude that long X-linked genes are less, if at all, affected by loss of the MSL-complex.
The finding that long genes are less dependent on MSL for dosage compensation is intriguing and is in line with the fact that it has previously been shown that MSL-bound genes are mainly housekeeping genes. We therefore looked at the genes that are unbound but still compensated and found that only 29% of them are classified as housekeeping (according to the definition in) compared to 70% of the bound genes. Kharchenko et al. found that long genes are preferentially found in regions enriched for the enhancer type of chromatin (chromatin state 3). We speculate that long genes are more dynamically regulated and are less dependent on the MSL-complex for compensation of their expression in males. The binding of MSL to some long genes is probably due to their proximity to high affinity sites. Unbound long genes (>10 kb) are on average positioned 48 kb from a high affinity site (defined by combining the sites reported by Alekseyenko et al. and Straub et al.) while bound long genes are on average 19 kb away.
Recently, two contradictory studies tried to identify the mechanism behind MSL-dependent compensation[8, 45]. The results of Larschan et al. indicated enhanced elongation efficiency, but it was later questioned whether enhanced elongation efficiency on the X-chromosome could be demonstrated using the same data. When chromosome X is compared to individual autosomes, no enhanced elongation efficiency can be seen[8, 46]. Conrad et al. provided evidence for increased RNA Pol2 initiation on MSL-bound genes. Recently, problems with the analysis pipeline applied to these data have been pointed out[47–49]. However, even after correcting for the pre-processing issues, an approximately 20% increase of Pol2 initiation on MSL-bound genes was observed[47–49]. This is compatible with an MSL-dependent mechanism that recruits Pol2 and with the fact that the loss of MSL results in only an approximately 20% drop in expression. Although additional evidence is needed, it has been proposed that the MSL-complex sequesters MOF from autosomes[21, 22] and constrains its activity on the X-chromosome[20, 50]. Clearly, the MSL-complex has an essential function, as mutations affecting its components are male-specific lethal, and we speculate that this function is to enhance the expression of genes that lack effective feed-back regulation; these seem to be mostly short housekeeping genes. Although the mechanisms of the MSL independent compensation mechanism remain elusive, they clearly make a significant contribution to dosage compensation of the X-chromosome and merit more focused investigation.
The Drosophila male X-chromosome is extremely well compensated and a significant number of genes are compensated without being bound by any detectable levels of MSL-complex. Taken together with the fact that even the most strongly affected genes (those with robust MSL binding) in MSL mutants/knock-downs are still expressed at about 80% of wild-type (50% is expected if the MSL-complex is the only player in dosage compensation), the majority of the dosage compensation of all male X-linked genes is likely mediated by an MSL-independent mechanism. Similar to the MSL-mediated part of dosage compensation, we do not yet know the mechanism behind the MSL-independent compensation other than that it appears to act preferentially on long genes. We hope that future studies not exclusively focused on the MSL-complex will shed light on these important MSL-independent mechanisms.
Micro-array expression data analysis
Gene expression levels in wild-type, RNAi, and mutant samples were computed from Affymetrix .CEL files of raw gene expression data using Robust Multi-Array Average with the Bioconductor “affy” package. Each probeset was then mapped to the genomic Release 5 coordinates using the latest library files from Affymetrix. Only probeset mapping to a single gene were included in this study (in the main text these are referred to as genes) and the medians of biological replicates were used for downstream analysis. We defined expressed genes as those with expression values of 6 or more (log2) as previously described[18, 25, 26]. The following datasets were used for the analysis: (male, female wild-type – E-MEXP-3506, mof, msl1, msl3 RNAi – E-MEXP-1505, msl2 RNAi – http://compbio.med.harvard.edu/Supplements/GD05.html, msl2 mutant – GSE12054, roX1 roX2 mutant – GSE3990).
ChIP data analysis
Reads from MOF and H4K16ac immunoprecipitation (IP) and input samples (E-MTAB-911) were aligned against the Drosophila melanogaster reference sequence (release 5) using the Bowtie software package, mapping reads only to unique locations in the genome (parameters: -n 2 –k 1 --best). Read counts per nucleotide were calculated and median smoothed using a window size of 10 bp, and smoothed read counts were extracted at intervals of 10 bp across the genome. For each position in the IP sample with at least one mapped read, a ratio of IP-input (on a log2 scale) was calculated. If the read count in the input was below the read-count mean (in the input sample) it was set to the mean. If the input mean was below four the minimum input value was set to four (to avoid division by near-zero values). All ratio values were then adjusted by reducing each value by the median of the ratios. This linear adjustment was carried out in order to compensate for differences in IP and input sequencing depth. Ratio values were further median smoothed using a window size of 200 bp. Windows with fewer than 10 data points were discarded. Gene binding values for MOF and H4K16ac from salivary glands were calculated by first extracting all data points (at 10 bp resolution) within annotated exons of each gene. If the number of data points for a given gene exceeded 10, they were sorted and the gene binding value calculated as the mean of the 50% highest values as described by Johansson et al., with the exception that here we excluded the first 200 bp exonic sequence of each gene. Introns were excluded since both MOF and H4K16ac are enriched in exons. Using the top 50% of the exon values makes the mean less sensitive to alternatively spliced exons. We exclude genes with fewer than 10 data points to avoid uncertain binding values from short genes. We also calculated gene binding values for MOF, MSL1, and MSL3 in S2 cells from data (E-MEXP-1508) in the study by Kind et al..
To calculate copy number ratio for each gene between males and females, read counts per nucleotide in male and female input samples were first median smoothed using a window size of 1 kb and smoothed read counts were extracted at intervals of 10 bp across the genome. All such 10-step positions within compensated genes were averaged and a male/female ratio (on a log2 scale) was calculated.
Statistical analysis and visualization of data
All statistical analyses and visualizations were performed using Statistica (Statsoft, USA) and Evince (UmBio, Sweden). Transcript length was calculated, using Flybase annotation Release 5.43, as the difference between Gene Start and End positions. To calculate the distances between genes and High Affinity Sites (HAS), we used 188 HAS defined by the union of the HAS reported by Alekseyenko et al. and Straub et al.. Before normalization, all data sets were tested for skewness as described by Landfors et al..
High affinity sites
Male absent on the first
Male-specific Lethal 1
Male-specific Lethal 2
Male-specific Lethal 3.
Deng X, Disteche CM: Genomic responses to abnormal gene dosage: the X chromosome improved on a common strategy. PLoS Biol. 2010, 8: e1000318-10.1371/journal.pbio.1000318.
Oliver B: Sex, dose, and equality. PLoS Biol. 2007, 5: e340-10.1371/journal.pbio.0050340.
Zhang Y, Oliver B: Dosage compensation goes global. Curr Opin Genet Dev. 2007, 17: 113-120. 10.1016/j.gde.2007.02.002.
Bai X, Alekseyenko AA, Kuroda MI: Sequence-specific targeting of MSL complex regulates transcription of the roX RNA genes. EMBO J. 2004, 23: 2853-2861. 10.1038/sj.emboj.7600299.
Hamada FN, Park PJ, Gordadze PR, Kuroda MI: Global regulation of X chromosomal genes by the MSL complex in Drosophila melanogaster. Genes Dev. 2005, 19: 2289-2294. 10.1101/gad.1343705.
Straub T, Gilfillan GD, Maier VK, Becker PB: The Drosophila MSL complex activates the transcription of target genes. Genes Dev. 2005, 19: 2284-2288. 10.1101/gad.1343105.
Alekseyenko AA, Peng S, Larschan E, Gorchakov AA, Lee OK, Kharchenko P, McGrath SD, Wang CI, Mardis ER, Park PJ, Kuroda MI: A sequence motif within chromatin entry sites directs MSL establishment on the Drosophila X chromosome. Cell. 2008, 134: 599-609. 10.1016/j.cell.2008.06.033.
Conrad T, Cavalli FM, Vaquerizas JM, Luscombe NM, Akhtar A: Drosophila dosage compensation involves enhanced Pol II recruitment to male X-linked promoters. Science. 2012, 337: 742-746. 10.1126/science.1221428.
Gelbart ME, Kuroda MI: Drosophila dosage compensation: a complex voyage to the X chromosome. Development. 2009, 136: 1399-1410. 10.1242/dev.029645.
Straub T, Becker PB: Dosage compensation: the beginning and end of generalization. Nat Rev Genet. 2007, 8: 47-57.
Alekseyenko AA, Larschan E, Lai WR, Park PJ, Kuroda MI: High-resolution ChIP-chip analysis reveals that the Drosophila MSL complex selectively identifies active genes on the male X chromosome. Genes Dev. 2006, 20: 848-857. 10.1101/gad.1400206.
Gilfillan GD, Straub T, De Wit E, Greil F, Lamm R, Van Steensel B, Becker PB: Chromosome-wide gene-specific targeting of the Drosophila dosage compensation complex. Genes Dev. 2006, 20: 858-870. 10.1101/gad.1399406.
Legube G, McWeeney SK, Lercher MJ, Akhtar A: X-chromosome-wide profiling of MSL-1 distribution and dosage compensation in Drosophila. Genes Dev. 2006, 20: 871-883. 10.1101/gad.377506.
Larschan E, Alekseyenko AA, Gortchakov AA, Peng S, Li B, Yang P, Workman JL, Park PJ, Kuroda MI: MSL complex is attracted to genes marked by H3K36 trimethylation using a sequence-independent mechanism. Mol Cell. 2007, 28: 121-133. 10.1016/j.molcel.2007.08.011.
Gelbart ME, Larschan E, Peng S, Park PJ, Kuroda MI: Drosophila MSL complex globally acetylates H4K16 on the male X chromosome for dosage compensation. Nat Struct Mol Biol. 2009, 16: 825-832. 10.1038/nsmb.1644.
Zhang Y, Malone JH, Powell SK, Periwal V, Spana E, Macalpine DM, Oliver B: Expression in aneuploid Drosophila S2 cells. PLoS Biol. 2010, 8: e1000320-10.1371/journal.pbio.1000320.
Stenberg P, Larsson J: Buffering and the evolution of chromosome-wide gene regulation. Chromosoma. 2011, 120: 213-225. 10.1007/s00412-011-0319-8.
Lundberg LE, Figueiredo ML, Stenberg P, Larsson J: Buffering and proteolysis are induced by segmental monosomy in Drosophila melanogaster. Nucleic Acids Res. 2012, 40: 5926-5937. 10.1093/nar/gks245.
Prestel M, Feller C, Becker PB: Dosage compensation and the global re-balancing of aneuploid genomes. Genome Biol. 2010, 11: 216-10.1186/gb-2010-11-8-216.
Sun L, Fernandez HR, Donohue RC, Li J, Cheng J, Birchler JA: Male-specific lethal complex in Drosophila counteracts histone acetylation and does not mediate dosage compensation. Proc Natl Acad Sci USA. 2013, 110: E808-E817. 10.1073/pnas.1222542110.
Bhadra U, Pal-Bhadra M, Birchler JA: Role of the male specific lethal (msl) genes in modifying the effects of sex chromosomal dosage in Drosophila. Genetics. 1999, 152: 249-268.
Bhadra MP, Bhadra U, Kundu J, Birchler JA: Gene expression analysis of the function of the male-specific lethal complex in Drosophila. Genetics. 2005, 169: 2061-2074. 10.1534/genetics.104.036020.
Deng X, Meller VH: roX RNAs are required for increased expression of X-linked genes in Drosophila melanogaster males. Genetics. 2006, 174: 1859-1866. 10.1534/genetics.106.064568.
Deng X, Koya SK, Kong Y, Meller VH: Coordinated regulation of heterochromatic genes in Drosophila melanogaster males. Genetics. 2009, 182: 481-491. 10.1534/genetics.109.102087.
Stenberg P, Lundberg LE, Johansson AM, Rydén P, Svensson MJ, Larsson J: Buffering of segmental and chromosomal aneuploidies in Drosophila melanogaster. PLoS Genet. 2009, 5: e1000465-10.1371/journal.pgen.1000465.
Philip P, Pettersson F, Stenberg P: Sequence signatures involved in targeting the Male-Specific Lethal complex to X-chromosomal genes in Drosophila melanogaster. BMC Genomics. 2012, 13: 97-10.1186/1471-2164-13-97.
Gupta V, Parisi M, Sturgill D, Nuttall R, Doctolero M, Dudko OK, Malley JD, Eastman PS, Oliver B: Global analysis of X-chromosome dosage compensation. J Biol. 2006, 5: 3-10.1186/jbiol30.
Conrad T, Cavalli FM, Holz H, Hallacli E, Kind J, Ilik I, Vaquerizas JM, Luscombe NM, Akhtar A: The MOF chromobarrel domain controls genome-wide H4K16 acetylation and spreading of the MSL complex. Dev Cell. 2012, 22: 610-624. 10.1016/j.devcel.2011.12.016.
Kind J, Vaquerizas JM, Gebhardt P, Gentzel M, Luscombe NM, Bertone P, Akhtar A: Genome-wide analysis reveals MOF as a key regulator of dosage compensation and gene expression in Drosophila. Cell. 2008, 133: 813-828. 10.1016/j.cell.2008.04.036.
Straub T, Zabel A, Gilfillan GD, Feller C, Becker PB: Different chromatin interfaces of the Drosophila dosage compensation complex revealed by high-shear ChIP-seq. Genome Res. 2013, 23: 473-485. 10.1101/gr.146407.112.
Lam KC, Mühlpfordt F, Vaquerizas JM, Raja SJ, Holz H, Luscombe NM, Manke T, Akhtar A: The NSL complex regulates housekeeping genes in Drosophila. PLoS Genet. 2012, 8: e1002736-10.1371/journal.pgen.1002736.
Georgiev P, Chlamydas S, Akhtar A: Drosophila dosage compensation: males are from Mars, females are from Venus. Fly. 2011, 5: 147-154.
Alekseyenko AA, Ho JW, Peng S, Gelbart M, Tolstorukov MY, Plachetka A, Kharchenko PV, Jung YL, Gorchakov AA, Larschan E, Gu T, Minoda A, Riddle NC, Schwartz YB, Elgin SC, Karpen GH, Pirrotta V, Kuroda MI, Park PJ: Sequence-specific targeting of dosage compensation in Drosophila favors an active chromatin context. PLoS Genet. 2012, 8: e1002646-10.1371/journal.pgen.1002646.
Belyaeva ES, Goncharov FP, Demakova OV, Kolesnikova TD, Boldyreva LV, Semeshin VF, Zhimulev IF: Late replication domains in polytene and non-polytene cells of Drosophila melanogaster. PLoS One. 2012, 7: e30035-10.1371/journal.pone.0030035.
Alekseyenko AA, Demakova OV, Belyaeva ES, Makarevich GF, Kotlikova IV, Nöthiger R, Zhimulev IF: Dosage compensation and intercalary heterochromatin in X chromosomes of Drosophila melanogaster. Chromosoma. 2002, 111: 106-113. 10.1007/s00412-002-0191-7.
Makunin IV, Volkova EI, Belyaeva ES, Nabirochkina EN, Pirrotta V, Zhimulev IF: The Drosophila suppressor of underreplication protein binds to late-replicating regions of polytene chromosomes. Genetics. 2002, 160: 1023-1034.
Gergen JP: Dosage compensation in Drosophila: evidence that daughterless and sex-lethal control X chromosome activity at the blastoderm stage of embryogenesis. Genetics. 1987, 117: 477-485.
Kelley RL, Solovyeva I, Lyman LM, Richman R, Solovyev V, Kuroda MI: Expression of msl-2 causes assembly of dosage compensation regulators on the X chromosomes and female lethality in Drosophila. Cell. 1995, 81: 867-877. 10.1016/0092-8674(95)90007-1.
Lott SE, Villalta JE, Schroth GP, Luo S, Tonkin LA, Eisen MB: Noncanonical compensation of zygotic X transcription in early Drosophila melanogaster development revealed through single-embryo RNA-seq. PLoS Biol. 2011, 9: e1000590-10.1371/journal.pbio.1000590.
Hiebert JC, Birchler JA: Effects of the maleless mutation on X and autosomal gene expression in Drosophila melanogaster. Genetics. 1994, 136: 913-926.
Bhadra U, Pal-Bhadra M, Birchler JA: Histone acetylation and gene expression analysis of sex lethal mutants in Drosophila. Genetics. 2000, 155: 753-763.
Malone JH, Cho DY, Mattiuzzo NR, Artieri CG, Jiang L, Dale RK, Smith HE, McDaniel J, Munro S, Salit M, Andrews J, Przytycka TM, Oliver B: Mediation of Drosophila autosomal dosage effects and compensation by network interactions. Genome Biol. 2012, 13: r28-10.1186/gb-2012-13-4-r28.
Kharchenko PV, Alekseyenko AA, Schwartz YB, Minoda A, Riddle NC, Ernst J, Sabo PJ, Larschan E, Gorchakov AA, Gu T, Linder-Basso D, Plachetka A, Shanower G, Tolstorukov MY, Luquette LJ, Xi R, Jung YL, Park RW, Bishop EP, Canfield TK, Sandstrom R, Thurman RE, MacAlpine DM, Stamatoyannopoulos JA, Kellis M, Elgin SC, Kuroda MI, Pirrotta V, Karpen GH, Park PJ: Comprehensive analysis of the chromatin landscape in Drosophila melanogaster. Nature. 2011, 471: 480-485. 10.1038/nature09725.
Straub T, Grimaud C, Gilfillan GD, Mitterweger A, Becker PB: The chromosomal high-affinity binding sites for the Drosophila dosage compensation complex. PLoS Genet. 2008, 4: e1000302-10.1371/journal.pgen.1000302.
Larschan E, Bishop EP, Kharchenko PV, Core LJ, Lis JT, Park PJ, Kuroda MI: X chromosome dosage compensation via enhanced transcriptional elongation in Drosophila. Nature. 2011, 471: 115-118. 10.1038/nature09757.
Johansson AM, Stenberg P, Allgardsson A, Larsson J: POF regulates the expression of genes on the fourth chromosome in Drosophila melanogaster by binding to nascent RNA. Mol Cell Biol. 2012, 32: 2121-2134. 10.1128/MCB.06622-11.
Straub T, Becker PB: Comment on “Drosophila dosage compensation involves enhanced Pol II recruitment to male X-linked promoters”. Science. 2013, 340: 273-Discussion 273
Ferrari F, Jung YL, Kharchenko PV, Plachetka A, Alekseyenko AA, Kuroda MI, Park PJ: Comment on “Drosophila dosage compensation involves enhanced Pol II recruitment to male X-linked promoters”. Science. 2013, 340: 273-Discussion 273
Vaquerizas JM, Cavalli FM, Conrad T, Akhtar A, Luscombe NM: Response to Comments on “Drosophila Dosage Compensation Involves Enhanced Pol II Recruitment to Male X-Linked Promoters”. Science. 2013, 340: 273.
Prestel M, Feller C, Straub T, Mitlöhner H, Becker PB: The activation potential of MOF is constrained for dosage compensation. Mol Cell. 2010, 38: 815-826. 10.1016/j.molcel.2010.05.022.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.
Gautier L, Cope L, Bolstad BM, Irizarry RA: affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Tweedie S, Ashburner M, Falls K, Leyland P, McQuilton P, Marygold S, Millburn G, Osumi-Sutherland D, Schroeder A, Seal R, Zhang H, FlyBase Consortium: FlyBase: enhancing Drosophila Gene Ontology annotations. Nucleic Acids Res. 2009, 37: D555-D559. 10.1093/nar/gkn788.
Landfors M, Philip P, Rydén P, Stenberg P: Normalization of high dimensional genomics data where the distribution of the altered variables is skewed. PLoS One. 2011, 6: e27942-10.1371/journal.pone.0027942.
This work was supported by the Magnus Bergvall and Åke Wiberg foundations to PS and the JC Kempe, Nilsson-Ehle, and Lars Hierta foundations to PP.
The authors declare they have no competing interests.
PP and PS designed the study, performed the analysis and wrote the paper. Both authors read and approved the final manuscript.