Skip to main content

Male X-linked genes in Drosophila melanogaster are compensated independently of the Male-Specific Lethal complex

Abstract

Background

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.

Results

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.

Conclusions

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.

Background

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[13]. 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[48]. 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[1114]. 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[15]. 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[16]. 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[1619]. 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.[2022], 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[8]. 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).

Figure 1
figure 1

Male minus female expression in salivary glands (log 2 ) of all genes expressed in both sexes along chromosomes X and 3L. Dotted grey lines indicate ± one standard deviation from the mean of chromosome 3L.

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[1114]. 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[28]. 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).

Figure 2
figure 2

A significant part of the genes on the X chromosome are dosage compensated, unbound by MOF and not H4K16 acetylated in salivary glands. Male minus female expression (log2) of all genes expressed in both sexes plotted versus MOF gene binding values on (A) the X-chromosome and (B) chromosome 3L. Horizontal dotted grey lines indicate ± one standard deviation from the mean expression ratio for chromosome 3L and vertical lines the mean MOF binding on chromosome 3L plus one standard deviation. (C) Average H4K16ac in four non-overlapping groups with different levels of MOF binding, (Unbound n = 122, >3 n = 220, >4 n = 227, >5 n = 227).

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[31], 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[8] 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.[36]; 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[26], 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, 3741], 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.

Figure 3
figure 3

MOF unbound genes in salivary glands are less affected in MSL-complex mutants/knock-downs. (A) Expression in S2 MOF RNAi lines (log2 expression of MOF RNAi-GFP RNAi, Unbound n = 86, >3 n = 208, >4 n = 222, >5 n = 226). (B) Expression in roX1 roX2 mutant larvae (log2 expression of roX1 roX2-WT, Unbound n = 100, >3 n = 212, >4 n = 205, >5 n = 205).

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.[16], 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[42]. The general buffering of autosomal deficiencies reported by Lundberg et al. showed stronger buffering of long genes[18]; 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.

Figure 4
figure 4

Long genes are less dependent on the MSL-complex. (A) Average transcript length in four non-overlapping groups of expressed genes with different levels of MOF binding in salivary glands (Unbound n = 122, >3 n = 220, >4 n = 227, >5 n = 227). (B) Expression in msl2 mutant male larvae (relative to wild-type) plotted versus transcript length.

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.[29]. 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[12]. 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[25]) 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)[43]. 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.[7] and Straub et al.[44]) 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[45]. 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[8]. Recently, problems with the analysis pipeline applied to these data have been pointed out[4749]. However, even after correcting for the pre-processing issues, an approximately 20% increase of Pol2 initiation on MSL-bound genes was observed[4749]. 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.

Conclusions

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.

Methods

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[51] with the Bioconductor “affy” package[52]. 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[53], 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.[46], 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.[29].

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[54], 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.[7] and Straub et al.[44]. Before normalization, all data sets were tested for skewness as described by Landfors et al.[55].

Abbreviations

HAS:

High affinity sites

IP:

Immunoprecipitation

MOF:

Male absent on the first

MSL:

Male-specific Lethal

MSL1:

Male-specific Lethal 1

MSL2:

Male-specific Lethal 2

MSL3:

Male-specific Lethal 3.

References

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

    Article  PubMed Central  PubMed  Google Scholar 

  2. Oliver B: Sex, dose, and equality. PLoS Biol. 2007, 5: e340-10.1371/journal.pbio.0050340.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Zhang Y, Oliver B: Dosage compensation goes global. Curr Opin Genet Dev. 2007, 17: 113-120. 10.1016/j.gde.2007.02.002.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  9. Gelbart ME, Kuroda MI: Drosophila dosage compensation: a complex voyage to the X chromosome. Development. 2009, 136: 1399-1410. 10.1242/dev.029645.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Straub T, Becker PB: Dosage compensation: the beginning and end of generalization. Nat Rev Genet. 2007, 8: 47-57.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

  17. Stenberg P, Larsson J: Buffering and the evolution of chromosome-wide gene regulation. Chromosoma. 2011, 120: 213-225. 10.1007/s00412-011-0319-8.

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Georgiev P, Chlamydas S, Akhtar A: Drosophila dosage compensation: males are from Mars, females are from Venus. Fly. 2011, 5: 147-154.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    PubMed Central  CAS  PubMed  Google Scholar 

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

    PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Hiebert JC, Birchler JA: Effects of the maleless mutation on X and autosomal gene expression in Drosophila melanogaster. Genetics. 1994, 136: 913-926.

    PubMed Central  CAS  PubMed  Google Scholar 

  41. Bhadra U, Pal-Bhadra M, Birchler JA: Histone acetylation and gene expression analysis of sex lethal mutants in Drosophila. Genetics. 2000, 155: 753-763.

    PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Per Stenberg.

Additional information

Competing interests

The authors declare they have no competing interests.

Authors’ contributions

PP and PS designed the study, performed the analysis and wrote the paper. Both authors read and approved the final manuscript.

Electronic supplementary material

13072_2013_303_MOESM1_ESM.pdf

Additional file 1: Figure S1: Male minus female expression (log2) of all dosage compensated genes plotted versus MOF gene binding values on the X chromosome (salivary gland data). Coloring from red (strong) to dark blue (weak) based on MOF, MSL-1 and MSL-3 gene binding values in the S2 cell-line. All data are in log2. (PDF 510 KB)

13072_2013_303_MOESM2_ESM.pdf

Additional file 2: Figure S2: Expression of genes with variable MOF gene binding values in A) msl2 mutant versus wild type larvae (Unbound n = 93, >3 n = 208, >4 n = 211, >5 n = 202), B) msl2 RNAi versus control RNAi in S2 (Unbound n = 69, >3 n = 205, >4 n = 222, >5 n = 224), C) msl1 RNAi versus control RNAi in S2 (Unbound n = 86, >3 n = 208, >4 n = 218, >5 n = 223), and D) msl3 RNAi versus control RNAi in S2 (Unbound n = 86, >3 n = 208, >4 n = 222, >5 n = 226). All data are in log2. (PDF 280 KB)

13072_2013_303_MOESM3_ESM.pdf

Additional file 3: Figure S3: Male minus female expression in salivary glands (log2) of all genes expressed in both sexes along chromosomes 2L, 2R and 3R. (PDF 467 KB)

Authors’ original submitted files for images

Rights and permissions

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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Philip, P., Stenberg, P. Male X-linked genes in Drosophila melanogaster are compensated independently of the Male-Specific Lethal complex. Epigenetics & Chromatin 6, 35 (2013). https://doi.org/10.1186/1756-8935-6-35

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1756-8935-6-35

Keywords