N-terminus of Drosophila melanogaster MSL1 is critical for dosage compensation
eLife assessment
This is a potentially valuable contribution, reporting a deletion analysis of the MSL1 gene to assess how different parts of the protein product interact with the MSL2 protein and roX RNA to affect the association of the MSL complex with the male X chromosome of Drosophila. However, the framework that the MSL complex mediates dosage compensation is outdated and has flaws, and the evidence is currently considered inadequate to support the claims. Because there are many ways to alter viability, sex-specific viability is insufficient to make claims regarding dosage compensation.
https://doi.org/10.7554/eLife.93241.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Inadequate: Methods, data and analyses do not support the primary claims
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
The male-specific lethal complex (MSL), which consists of five proteins and two non-coding roX RNAs, is involved in the transcriptional enhancement of X-linked genes to compensate for the sex chromosome monosomy in Drosophila XY males compared with XX females. The MSL1 and MSL2 proteins form the heterotetrameric core of the MSL complex and are critical for the specific recruitment of the complex to the high-affinity ‘entry’ sites (HAS) on the X chromosome. In this study, we demonstrated that the N-terminal region of MSL1 is critical for stability and functions of MSL1. Amino acid deletions and substitutions in the N-terminal region of MSL1 strongly affect both the interaction with roX2 RNA and the MSL complex binding to HAS on the X chromosome. In particular, substitution of the conserved N-terminal amino-acids 3–7 in MSL1 (MSL1GS) affects male viability similar to the inactivation of genes encoding roX RNAs. In addition, MSL1GS binds to promoters such as MSL1WT but does not co-bind with MSL2 and MSL3 to X chromosomal HAS. However, overexpression of MSL2 partially restores the dosage compensation. Thus, the interaction of MSL1 with roX RNA is critical for the efficient assembly of the MSL complex on HAS of the male X chromosome.
Introduction
At present, it is unknown how transcription-regulating complexes bind only to certain chromatin regions that do not have a pronounced sequence specificity relative to other genome sites. Dosage compensation in Drosophila is a striking example of the specific recruitment of transcription complexes (Lucchesi, 2018). Dosage compensation equalizes the expression of different doses of X chromosomes in females (X/X) and males (X/Y), which is achieved by doubling the expression of genes on one X chromosome in males to match the level of gene expression on two chromosomes in females (Ferrari et al., 2014; Kuroda et al., 2016; Lucchesi, 2018; Samata and Akhtar, 2018). Currently, there is no consensus on the mechanisms that determine dosage compensation (Birchler and Veitia, 2021; Vensko and Stone, 2015). It is most likely that dosage compensation in males is the result of the simultaneous effect of activation of gene transcription on the X chromosome and the inverse dosage effect, which consists of a decrease in the expression of autosomal genes (Birchler, 2016; Ferrari et al., 2014; Kuroda et al., 2016; Lucchesi, 2018; Samata and Akhtar, 2018; Zhang et al., 2021).
A multisubunit complex that binds only to the X chromosome of males was discovered (Ferrari et al., 2014; Kuroda et al., 2016; Lucchesi, 2018; Samata and Akhtar, 2018), and it was proposed that this complex is involved in equalizing the expression of genes on the X chromosomes of males and females. The complex, named male-specific lethal, consists of five proteins, male-specific lethal (MSL)1, MSL2, MSL3, males-absent on the first protein (MOF), and maleless (MLE). Two non-coding RNAs are included, namely roX1 (3.7 kb) and roX2 (0.6 kb), which perform partially redundant functions (Kim et al., 2018; Kuroda et al., 2016; Meller and Rattner, 2002; Samata and Akhtar, 2018). The proteins MSL1, MSL3, MOF, and MLE are also expressed in females. These proteins, and are involved in the regulation of gene expression processes that are unrelated to X chromosome dosage compensation (Samata and Akhtar, 2018). The MSL2 protein is male-specific, and the expression of the msl-2 gene is controlled by the sex-lethal (SXL) protein, which binds to msl-2 mRNA in females and blocks protein translation (Beckmann et al., 2005). Therefore, the MSL2 protein is thought to play a key role in the specific binding of the MSL complex to the X chromosome in males.
The N-terminus of MSL2 contains a really interesting new gene (RING) domain, which functions as a ubiquitin E3 ligase (Villa et al., 2012; Wu et al., 2011). MSL2 has been shown to induce the ubiquitination of itself as well as other subunits of the MSL complex. The MSL1 protein functions as a scaffold for the assembly of the MSL complex (Gorman et al., 1993; Kadlec et al., 2011; Palmer et al., 1993; Scott et al., 2000), as the N-terminal coiled-coil domain (Hallacli et al., 2012) of MSL1 facilitates homodimerization and interacts with the N-terminal RING domain of MSL2. The PEHE domain located in the C-terminal region of MSL1 is responsible for interactions with MSL3 and MOF (Kadlec et al., 2011; Scott et al., 2000). MSL1 and MSL2 form the structural core of the complex, whereas MLE—an ATP-dependent RNA/DNA helicase belonging to the DEAD box family—interacts with the roX1 and roX2 RNAs to induce their unwinding (Ilik et al., 2017; Maenner et al., 2013; Meller et al., 2000). This interaction allows both RNAs to bind MSL2 and possibly MSL1 (Müller et al., 2020), combining MLE, MSL2, and other MSL proteins into a single complex.
In addition to the MSL complex, lysine acetyltransferase MOF forms the catalytic core of the nonspecific lethal complex (NSL), which is recruited to promoters located on all chromosomes of both sexes and is responsible for strong transcriptional activation (Sheikh et al., 2019). Unlike the MSL complex, in which MOF moderately stimulates transcription during RNA polymerase II elongation (Ferrari et al., 2014; Kuroda et al., 2016), the NSL complex strongly increases the transcription at the initiation step (Sheikh et al., 2019). In males, formation of the MSL complex titrates MOF from the NSL complex, resulting in a decrease in its activity at all promoters (Prestel et al., 2010; Sun et al., 2013). This may be one of the reasons for the inverse dose effect, which consists of a decrease in the expression of autosomal genes (Birchler and Veitia, 2021).
An incomplete MSL1-MSL2 subcomplex can bind to a reproducible set of several hundreds of sites on the X chromosome, referred to as the primary chromatin entry sites (CESs) (Alekseyenko et al., 2008) or HAS (Straub et al., 2008). A GA-rich motif within the HAS/CES, which serves as an MSL recognition element (MRE), is important for the MSL complex targeting (Alekseyenko et al., 2008; Straub et al., 2008). Structural analysis and an in vitro genome-wide DNA binding assay showed that MSL2, through its CXC and proline/basic-residue-rich domains, specifically recognizes and binds to a subclass of HAS-named pioneering sites on the X (PionX) (Fauth et al., 2010; Villa et al., 2016; Zheng et al., 2014).
The zinc-finger protein chromatin-linked adaptor for MSL proteins (CLAMP) was found to be significant for the specific recruitment of the MSL complex to the male X chromosome (Eggers et al., 2023; Larschan et al., 2012; Soruco et al., 2013). CLAMP is an essential protein in Drosophila that binds thousands of GA-rich sequences genome-wide in both sexes (Kuzu et al., 2016; Soruco et al., 2013; Urban et al., 2017). In particular, CLAMP binds to GA repeats in HAS and is involved in the organization of open chromatin in these genome regions (Albig et al., 2019; Urban et al., 2017). CLAMP has an N-terminal unstructured dimerization domain (Tikhonova et al., 2022a) and promotes long-range chromatin interactions that mediate the dosage compensation of the male X chromosome (Jordan and Larschan, 2021). CLAMP also directly interacts with the MSL2 (Albig et al., 2019; Eggers et al., 2023; Tikhonova et al., 2022b; Tikhonova et al., 2019) and MLE (Quinn et al., 2014; Tikhonova et al., 2024) proteins.
The CLAMP interaction and CXC domains of MSL2 have partially redundant functions for the recruitment of the MSL complex to the X chromosome (Tikhonova et al., 2019). A previous study (Li et al., 2005) showed that the N-terminal 265 amino acids region of MSL1 could bind to approximately 30 sites on polytene X chromosomes in females that exogenously express MSL2, and the first 26 aa are essential for such binding. The authors suggested that the N-terminal sequences of MSL1 can bind to DNA and facilitate the specific recruitment of MSL to the male X chromosome. However, further studies examining the mechanisms of the MSL complex recruitment to the male X chromosome have concentrated on MSL2 as this is the only dosage compensation protein that is exclusively expressed in males [see review Kuroda et al., 2016; Lucchesi, 2018; Samata and Akhtar, 2018]. Although MSL1 represents a core protein of the MSL complex (Gorman et al., 1993; Palmer et al., 1993; Scott et al., 2000), MSL1 can also independently interact with gene promoters (Chlamydas et al., 2016, p. 2; Straub et al., 2013).
Here, we tested the functional role of different regions in the N-terminal 1–85 aa region of MSL1. Deletion of the 1–15 aa results in the instability of MSL1. Unexpectedly, we found that the simultaneous replacement from 3 to 7 aa from the N-terminus or deletion from 41 to 65 aa of MSL1 strongly influences the binding of the MSL complex on the male X chromosome. These mutations also affect the interaction of MSL1 with roX2 RNA. Thus, the N-terminal region of MSL1 performs key functions in the recruitment of the MSL complex to the X chromosome.
Results
The N-terminal 85 amino acids of MSL1 are critical for the activity of the MSL complex
The principal goal of the work is to study the functional role of the N-terminal 85 aa MSL1 followed by the coiled-coil domain involved in homodimerization and interaction with MSL2. With the exception of the coiled-coil domain, the first 15 amino acids represent the most conserved region within the N-terminus of MSL1 among Drosophilidae (Figure 1A). However, among different families of Diptera (Figure 1—figure supplement 1), a high level of homology was only observed for the coiled-coil region, which is required for interactions with MSL2 (Hallacli et al., 2012). The N-terminal amino acids are also conserved among representatives within each of the Diptera families (Figure 1—figure supplement 1), but not across families. For example, the MSL1 of Anopheles mosquitoes, which are closely related to Drosophila melanogaster, features a completely different N-terminal domain.
-
Figure 1—source data 1
- https://cdn.elifesciences.org/articles/93241/elife-93241-fig1-data1-v1.zip
-
Figure 1—source data 2
- https://cdn.elifesciences.org/articles/93241/elife-93241-fig1-data2-v1.zip
-
Figure 1—source data 3
- https://cdn.elifesciences.org/articles/93241/elife-93241-fig1-data3-v1.xlsx
To evaluate the functional effects of N-terminal domain deletions in the MSL1 protein in vivo, we created transgenic flies expressing wild-type (MSL1WT) and mutant variants of MSL1 (MSL1Δ1-85, MSL1Δ1-15, MSL1Δ8-20, MSL1Δ24-39, MSL1Δ41-85, MSL1Δ41-65, and MSL1Δ66-85) (Figure 1A). The MSL1 variants were expressed under the control of a strong Ubi-p63E (Ubi) promoter (Figure 1—figure supplement 2). To avoid the influence of position effects on the expression of the MSL1 variants, all transgenes were inserted into the same genomic location (86F8) on the third chromosome using an φC31-based integration system (Bischof et al., 2007).
To determine the role played by the N-terminal sequence during the specific recruitment of the MSL complex in males, we crossed obtained transgenes into the null msl-1− background (the msl-1− background corresponds to the msl-1γ269/msl-1L60 heterozygote) (Palmer et al., 1993). As shown previously, the main band detected in immunoblot analysis corresponds to a 170-kD protein that is considerably larger than the predicted 1039 aa protein (Palmer et al., 1994). Anomalous electrophoretic migration was explained by specific clusters of proline, acidic, or basic residues as well as potential post-translational modifications in particular ubiquitination. MSL1 is expressed 2.5–3 times stronger under the control of the Ubi promoter in the transgenic line than in wild-type males (y1w1118 line) (Figure 1B, Figure 1—figure supplement 3A). Increased expression of MSL1 leads to the appearance of more MSL2, which is stabilized in a complex with MSL1. Similar results were obtained in females. In the absence of MSL2, the expression of MSL1 was lower in females compared to males (Figure 1B).
To compare the expression of the MSL1 variants in transgenic lines, we examined the levels of MSL1 expression in 2–3 day-old msl-1− females, which are homozygous for the transgenes expressing MSL1 variants. Immunoblot analysis showed that all MSL1 variants, with the exception of MSL1Δ1-85 and MSL1Δ1-15, were expressed in transgenic females at nearly equivalent levels (Figure 1C). Surprisingly, MSL1Δ1-85 and MSL1Δ1-15 could not be detected by immunoblotting, suggesting the instability of these MSL1 variants.
Homozygous msl-1− females remained viable, whereas males died during the late larva and early pupa stages (Gorman et al., 1993; Palmer et al., 1993). The transgenes, expressed MSL1WT and MSL1Δ66-85, compensated for the viability of msl-1− males, indicating that these MSL1 variants are functional (Figure 1D, Figure 1—source data 3). In MSL1Δ24-39 and MSL1Δ8-20 lines, male viability was reduced suggesting that dosage compensation is partially affected. In the end, all other tested mutant MSL1 proteins (MSL1Δ1-15, MSL1Δ1-85, MSL1Δ41-65, and MSL1Δ41-85) failed to compensate for the msl-1− mutation in males (Figure 1D). Taken together, these results suggested that at least two regions (1–15 aa and 41–65 aa) in the N-terminus of MSL1 are critical for the functional activity of the MSL complex.
We then asked how MSL2 expression might affect the stability of MSL1 variants. To answer this question, we measured MSL1, and MSL2 proteins in msl-1−; Ubi:msl-2WT-FLAG/ Ubi:msl-1* (Figure 1E). Again, immunoblotting showed that the amount of MSL1 was almost completely eliminated in females expressed either MSL1Δ1-15 and MSL2, or MSL1Δ1-85 and MSL2 (Figure 1E). All other MSL1 variants were present at the same level. Interestingly, there is no noticeable difference in the amount of MSL2 in all samples. Thus, the N-terminal 1–15 aa is essential for the stability of MSL1 independently of MSL2 presence.
The N-terminus of MSL1 regulates the binding of the MSL complex to the X chromosome in females expressing MSL2
We next examined the extent to which deletions in the N-terminal region of MSL1 would affect the recruitment efficiency of the MSL complex to the male X chromosome. Immunostaining of polytene chromosomes from the salivary glands of Drosophila larvae allows the visualization of proteins on interphase chromatin and has been used extensively to study dosage compensation (Dahlsveen et al., 2006; Demakova et al., 2003; Li et al., 2005; Lyman et al., 1997; Meller and Rattner, 2002; Palmer et al., 1994). In the MSL1Δ1-15, MSL1Δ41-85, MSL1Δ41-65, and MSL1Δ1-85 mutant males, rare larvae surviving to 2–3 stages had the polytene chromosomes with poor morphology. For this reason, we used a previously described sensitive model system to study the factors required for the recruitment of the MSL complex to the X chromosome in females (Chang and Kuroda, 1998; Li et al., 2008; Morra et al., 2011; Palmer et al., 1994; Zhou et al., 1995). It was shown (Tikhonova et al., 2019), that MSL1 and MSL2 were specifically recruited to the X chromosome in transgenic females expressing MSL2WT-FLAG under the control of the Ubi promoter inserted into attP-line 86Fb (Ubi:msl-2WT-FLAG).
When the MSL2WT-FLAG was expressed in the msl-1− background, we did not observe the binding of the MSL complex to the female X chromosome (data not shown). However, the binding of MSL1 and MSL2 to the X chromosome was detected in msl-1− females when MSL2WT-FLAG was expressed simultaneously with either MSL1WT or MSL1Δ66-85 (Figure 2, Figure 2—figure supplement 1). These results further supported that the 66–85-aa region is not essential for MSL1 function in dosage compensation.
The expression of all MSL1 variants with deletions of the first 15-aa region (MSL1Δ1-15 and MSL1Δ1-85) resulted in the near-complete loss of the MSL complex binding to the X chromosome in females (Figure 2, Figure 2—figure supplement 1). This result is expected as immunoblotting showed a strong decrease of MSL1 level in MSL1Δ1-15/MSL2WT and MSL1Δ1-85/MSL2WT females (Figure 1E). Previous studies (Dahlsveen et al., 2006; Demakova et al., 2003; Lyman et al., 1997; Meller and Rattner, 2002; Palmer et al., 1994) showed that the inactivation of MSL3, MLE, or roX RNAs resulted in the recruitment of the MSL complex to a small proportion of bands that correspond to the main HAS, including regions of the roX1 (3 F) and roX2 (10 C) genes. In polytene chromosomes in the MSL1Δ1-15 and MSL1Δ1-85 lines, we found only two intense and several faint MSL1/MSL2 bands (Figure 2—figure supplement 2). Interestingly, additional MSL1/MSL2 bands were identified on the polytene autosomes, suggesting that the MSL proteins began to bind to lower affinity sites on autosomes (Figure 2—figure supplement 1). Thus, the deletion of the N terminal 15 aa affects not only the stability of MSL1 but also the specificity of MSL binding to the X chromosome.
The deletion of the 41–85 aa or 41–65 aa regions did not affect the stability of MSL1, but at the same time significantly reduced the binding of MSL proteins to the X chromosome. MSL1Δ41-85 and MSL2 bind to several bands on the X chromosome, many of which can also be identified as weak bands in the MSL1Δ1-15 and MSL1Δ1-85 lines (Figure 2, Figure 2—figure supplements 1 and 2). In comparison with MSL1Δ1-15 and MSL1Δ1-85, MSL proteins bind more efficiently with the MSL1Δ41-85 and MSL1Δ41-65 polytene chromosomes. However, the binding of the MSL complex to the X chromosome in MSL1 mutants was much weaker than previously suggested for a minimal MSL complex lacking either MSL3 or MLE (Dahlsveen et al., 2006). As the 66–85-aa deletion did not affect the activity of MSL1, we concluded that the 41–65-aa region is critical for MSL1 function.
Deletion of 24–39 aa, which separates the 1–15 aa and 41–65 aa regions in MSL1, only moderately affects the binding of the MSL2 and MSL1Δ24-39 proteins on polytene chromosomes (Figure 2, Figure 2—figure supplement 1). This observation indicates that the 1–15 aa and 41–65 aa regions have at least partially independent functions in recruiting of the MSL complex to the X chromosome in males. In MSL1Δ8-20 female larvae patterns of MSL1 and MSL2 binding to polytene chromosome are the same as in larvae expressed MSL1Δ24-39 (Figure 2, Figure 2—figure supplements 1 and 2). Thus, a deletion from 8 aa up to 20 aa only partially affects the functions of MSL1, in contrast to the complete inactivation of this protein by deleting the terminal 15 aa (MSL1Δ1-15).
The N terminal 3-7 amino acids of MSL1 are critical for dosage compensation
Our results showed that MSL1Δ1-15 is unstable, while MSL1Δ8-20 only partially affects the specific recruitment of MSL to the X chromosome suggesting that the most terminal amino-acids are essential for functions. To test the role of the N-terminal amino acids (3-7) in the stability and functions of MSL1, we made several amino acids substitutions in this region. The N-terminal region contains well-conserved basic and aromatic amino acid residues (Figure 1A). We used data about the previously characterized MSL1 mutants (Li et al., 2005) in those three conserved basic amino acids, lysine 3, arginine 4, and lysine 6, or two conserved aromatic amino acids, phenylalanine 5 and tryptophan 7, were replaced by alanine. In our work, we replaced basic amino acids with serine, and aromatic amino acids – with alanine or glycine (Figure 3A). In result we obtained three different transgenes: MSL13S (K3S R4S K6S), MSL12A (F5A W7A), and MSL1GS (K3S R4S F5G K6S W7G).
-
Figure 3—source data 1
- https://cdn.elifesciences.org/articles/93241/elife-93241-fig3-data1-v1.zip
-
Figure 3—source data 2
- https://cdn.elifesciences.org/articles/93241/elife-93241-fig3-data2-v1.zip
-
Figure 3—source data 3
- https://cdn.elifesciences.org/articles/93241/elife-93241-fig3-data3-v1.xlsx
-
Figure 3—source data 4
- https://cdn.elifesciences.org/articles/93241/elife-93241-fig3-data4-v1.xlsx
To evaluate the functional effects of N-terminal amino acids of the MSL1 protein in vivo, we created transgenic flies expressing mutants MSL12A, MSL13S, and MSL1GS under the control of the Ubi promoter. All MSL1 variants were targeted by the HAx3 epitope at the C-termini (Figure 3A). MSL12A and MSL13S variants were expressed at levels comparable to MSL1WT in 2–3 day-old msl-1− females that are homozygous for the MSL1 transgenes (Figure 3B and C). At the same time, there was an approximately 1.5–2-fold decrease in the MSL1GS protein compared to MSL1WT.
The transgenes, expressed MSL12A and MSL13S, almost completely compensated for the viability of msl-1− males, indicating that these MSL1 variants are functional (Figure 3D, Figure 3—source data 3). In the msl-1− line expressing MSL1GS, the viability of males was strongly reduced. At the same time, increasing in MSL2 and a simultaneous decreasing in MSL1GS expression partially restore the viability of MSL1GS/MSL2 males (Figure 3E, Figure 3—source data 4). In a similar way, the viability of MSL1Δ8-20/MSL2 males also is close to wild-type, but MSL1Δ41-65/MSL2 and MSL1Δ1-15/MSL2 males are still lethal.
According to the results obtained on males, MSL12A and MSL13S showed binding patterns like MSL1WT on the larva polytene chromosome, while the MSL1GS protein failed to effectively bind on the X chromosome in complex with MSL2 (Figure 3F, Figure 3—figure supplements 1 and 2). There are also binding sites for only MSL2 or MSL1GS on the autosomes and the X chromosome. Taken the obtained results suggest that the N-terminal amino acids 3–7 additively contribute to specific recruitment of MSL on the X chromosome.
Next, we asked whether MSL1Δ1-15, MSL1Δ8-20, MSL1Δ24-39, MSL1Δ41-65, MSL1GS, MSL1Δ1-85, MSL1Δ41-85 variants exhibit similar stability when expressed in Drosophila male-like S2 cells. All MSL1 variants targeted by HAx3 were co-expressed with MSL2-FLAG in S2 cells (Figure 4A). MSL1WT was stable and effectively immunoprecipitated with MSL2 (Figure 4B). Similar results were obtained with MSL1Δ8-20, MSL1Δ24-39, MSL1Δ41-85, MSL12A, MSL13S, and MSL1GS (Figure 4B). Like in Drosophila, MSL1Δ1-15 and MSL1Δ1-85 were unstable and expressed at a very low level. These results confirm that the specific combinations of amino-acids at the N-terminus are required for the stability of the MSL1 protein.
The results in S2 cells showed that MSL1Δ41-85 and MSL1GS can co-immunoprecipitate with MSL2, as does MSL1WT (Figure 4B). These results suggest that mutant MSL1Δ41-65 and MSL1GS may be part of the MSL complex. Since S2 cells contain an endogenous MSL1 protein, it seems likely that the mutated MSL1 variants form heterodimers with wild-type ones. We then asked if these mutant protein variants could interact with endogenous MSL1 for specific recruitment to the X chromosome. To distinguish between WT and MSL1 mutants, we generated transgenic lines expressing MSL1Δ41-85-HA or MSL1GS-HA tagged with the HAx3 epitope at the msl-1+ background. In transgenic males and females expressed MSL2 and MSL1 variants in the heterozygous state, we observed complete co-localization between MSL2 and MSL1 variants marked by HAx3 epitope (+). Thus, in the presence of wild-type MSL1, MSL1GS and MSL1Δ41-85 variants enter the MSL complex, in which they specifically bind to the male X chromosome. These results suggest that MSL1GS and MSL1Δ41-85 effectively dimerize with the wild-type MSL1 via the coiled-coil domain.
The mutants MSL1GS and MSL1Δ41-85 disrupt the recruitment of the MSL complex to HAS
For a more complete understanding of the role played by the N-terminal region of MSL1 in the recruitment of the MSL complex to the X chromosome, we tested the binding of the MSL1, MSL2, and MSL3 proteins by chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) with a standard protocol for chromatin fragmentation based on sonication. We obtained MSL protein profiles from two biological replicates of 3-day-old adult males and females.
As the binding profiles of the MSL complex have been well-studied in cell lines, we compared MSL1, MSL2, and MSL3 binding efficiency at the HAS (Figure 5—figure supplement 1) with previously published data from S2 cells (Schauer et al., 2017; Straub et al., 2013). The ChIP-seq data from adult y1w1118 (WT) males showed MSL binding to approximately 90% of the previously identified HAS from three different studies (Alekseyenko et al., 2008; Ramírez et al., 2015; Straub et al., 2008; Figure 5—source data 1). As a control, the MSL1 and MSL3 proteins bind to autosomal sites but not to HAS in females (Figure 5—figure supplement 1A). These results indicate that most HAS are common between S2 cells and adult males.
Next, we compared the binding profiles of MSL1, MSL2, and MSL3 proteins obtained in msl-1−; Ubi:msl-2WT-FLAG/Ubi:msl-1WT females (M1[WT]+MSL2) and msl-1−; Ubi:msl-1WT/Ubi:msl-1WT males (M1[WT]) (Figure 5—figure supplement 2A). We found a good correlation between the obtained profiles, demonstrating that MSL2 expressed in females induces the binding of the MSL complex to the same HAS regions. We observed stronger binding of the MSL1, MSL2, and MSL3 proteins in M1[WT] (msl-1−; Ubi:msl-1WT/Ubi:msl-1WT) males compared with those in control y1w1118 males (Figure 5—figure supplement 1B, Figure 5—figure supplement 2A), which might be a consequence of the higher concentrations of MSL proteins in the transgenic line due to being expressed under the control of the Ubi promoter (Figure 1—figure supplement 3A). In M1[WT]+MSL2 females, the MSL proteins also bind to most HAS, but with lower efficiency (Figure 5—figure supplement 2A, B, C). These results are consistent with previous studies demonstrating that the ectopic expression of MSL2 in females is sufficient to trigger the correct targeting of MSL proteins to the X chromosome (Valsecchi et al., 2018).
To test the role played by the N-terminal region of MSL1, we compared the binding of MSL proteins in females expressing MSL2 combined with MSL1WT, MSL1Δ1-15, MSL1Δ41-85, and MSL1GS (Figure 5A). As shown above, MSL2 and MSL1WT effectively bind to most HAS in females. MSL3 was associated with extended regions surrounding the HAS, consistent with previously reported data (Larschan et al., 2007; Sural et al., 2008). In MSL1Δ1-15, MSL1Δ41-85, and MSL1GS flies binding of MSL1, MSL2, and MSL3 to X chromosomal HAS is strongly reduced (Figure 5A-E, Figure 5—figure supplement 3A). At the same time, we noticed that there are HAS in which MSL2 continues to bind effectively when MSL1 in MSL1Δ1-15, MSL1Δ41-85, and MSL1GS flies are strongly reduced (Figure 5—figure supplement 3B).
To investigate observed differences in localization of the MSL proteins at various genomic sites in female fly lines, we calculated the average logFC between normalized (RPKM) test signals and nonspecific IgG signals in regions where MSL1 and MSL2 colocalized in MSL1WT line on X chromosome (see Methods). Resulting list of peaks is mostly associated with HAS (217 out of 220 in total) and is highly enriched by the MSL1, MSL2, and MSL3 proteins in MSL1WT flies, suggesting the binding of the complete MSL complex. As a result of hierarchical clustering, these regions were divided into two groups based on the protein binding profiles of MSL1, MSL2, and MSL3. Cluster 1 (‘MSL2 Sensitive,’ N=133) contains sites to which proteins do not bind during expression of MSL1 mutants. Cluster 2 (‘MSL2 stable,’ N=87) includes sites to which MSL2 continues to bind even in the presence of MSL1 mutations (see Methods, Figure 6A, Figure 6—figure supplement 1, Figure 6—source data 1). In cluster 2 regions, in addition to MSL2, the mutant form of MSL1GS was also enriched. We decided to test whether the stability of MSL2 in the cluster 2 is associated with the presence of PionX sites, to which the MSL2 protein binds most effectively in vitro (Villa et al., 2016). Indeed, PionX sites made up a larger percentage of the peaks of the cluster 2 (17 out of 87), in contrast to cluster 1 (7 out of 133). Additionally, we tested the strength of MSL1 and MSL2 signals from female fly lines in peaks associated with PionX sites in comparison with all other regions from the cluster 2 (Figure 6B). As expected, we found that the MSL2 signal was significantly higher in the group of PionX sites in fly lines with the mutation (MSL1Δ1-15, MSL1Δ41-85, MSL1GS; Wilcoxon signed-rank test, p-value: 0.034, 7.3e-5, 2.2e-5, respectively), which was not noted for the MSL1WT (Wilcoxon signed-rank test, p-value: 0.71). Similar situation was observed for MSL1, whose signal at PionX sites was significantly greater for the MSL1GS females (Wilcoxon signed-rank test, p-value: 0.053) than for the MSL1WT females (Wilcoxon signed-rank test, p-value: 0.35). These results suggest that some of the PionX sites may be bound independently by the MSL2 protein.
Also, we analyzed regions with MSL1 alone peaks (no colocalization with MSL2) on all chromosomes (Figure 6—figure supplement 2). Previously, it was shown that MSL1 binds to promoters independently of the MSL complex (Chlamydas et al., 2016; Hallacli et al., 2012; Straub et al., 2013). Accordingly, MSL2 binding was weak or absent at all such sites. These sites are predominantly located in the promoter regions and MSL1GS binds to all sites with the same efficiency as MSL1WT, confirming that this mutant protein retains its functionality at the sites outside of HAS. The MSL1Δ1-15 and MSL1Δ41-85 are not associated with these sites. The absence of MSL1Δ1-15 can be explained by the low stability of this mutant protein. As MSL1Δ41-85 is expressed more strongly than MSL1GS, it seems likely that the 41–85 aa region facilitates the interaction of MSL1 with transcription factors associated with these promoters.
Two regions in the N-terminal domain are required for interaction with roX2 RNA
Our results suggest that MSL1GS and MSL1Δ41-85 can interact with wild-type MSL1 and form the complex. Previously it was suggested that the N-terminal part of MSL1 (191 aa) contributes to interaction with roX RNAs (Müller et al., 2020). It can be assumed that the ineffectiveness of MSL1 variants in dosage compensation is associated with their inability to interact with roX RNAs. To test this hypothesis, we obtained transgenic msl-1− (msl-1L60/msl-1γ269) females expressed either of MSL1 variants (MSL1WT-HA, MSL1GS-HA, and MSL1Δ41-85-HA) and roX2 RNA. To express roX2, we used the previously obtained roX2 gene under the control of the Ubi promoter in the 86Fb site (Tikhonova et al., 2022b). The roX2 transgene was expressed at the same level in all females expressed different MSL1 variants (Figure 7A).
The total extracts obtained from the 2–3 days old adult females (Figure 7) were used for immunoprecipitation with MSL1 antibodies. In all cases antibodies efficiently precipitate the tested MSL1 variants (Figure 7B). roX2 RNA was detected (Figure 7C) in precipitate obtained from MSL1WT extract (positive control) and was not found in precipitate obtained from msl-1− extract (negative control). In the case of MSL1GS and MSL1Δ41-85 precipitates, roX2 RNA was not detected as in the negative control. The obtained results suggest that MSL1WT interacts with roX2 in contrast to the MSL1GS and MSL1Δ41-85 mutants. These results are consistent with the hypothesis that failure to interact with roX RNAs is the primary cause of functional inactivation of MSL1GS and MSL1Δ41-85.
Discussion
This study focused on examining the role played by the MSL1 protein in recruiting the MSL complex in Drosophila to HAS. By UV-XL assay, it was shown that the N-terminal 191 aa of MSL1 and MSL2 are direct roX2 interactors (Müller et al., 2020). Here, we found that the 15 N-terminal aa are critical for the stability of the MSL1 protein. The substitution of the N-terminal amino-acids 3–7 (KRFKW) on the SSGSG in MSL1GS slightly reduces the stability of the mutant protein and strongly affects the interaction of MSL1GS with roX2. The 41–85 aa are also critical for MSL1 function and interaction with roX2 RNA in vivo. Both regions are separated by a linker that is less essential for MSL1 activity.
As functions of roX1 and roX2 are redundant (Meller and Rattner, 2002; Stuckenholz et al., 2003), it seems likely that roX1 interacts with the same N-terminal regions of MSL1 as roX2. Further study is required to understand the mechanisms for how N-terminal sequences of MSL1 determine the highly specific interaction with roX RNAs. In accordance with the key role played in the recruitment of the MSL complex to the male X chromosome, the N-terminal 15 amino acids of MSL1 are highly conserved among Drosophilidae. However, MSL1 orthologs in other families of Diptera have completely different N-terminal sequences despite the presence of high homology in the dimeric coiled-coil region. Even between Anopheles gambiae and Drosophila melanogaster, the X chromosome is thought to have formed independently from the same ancestral chromosome (Jiang et al., 2015; Pease and Hahn, 2012; Toups and Hahn, 2010).
The generally accepted model has suggested that most MSL complexes initially bind with numerous HAS on the male X chromosome (Kageyama et al., 2001; Kelley et al., 1997; Meller et al., 2000; Park et al., 2003; Park and Kuroda, 2001). Most HAS contain low-complexity GAGA motifs, referred to as the MRE (Alekseyenko et al., 2008; Straub et al., 2008). The CLAMP protein binds to MREs and is responsible for chromatin opening (Albig et al., 2019; Rieder et al., 2019). MSL2 contributes to the specific recognition of HAS through the direct binding of GAGA motifs (Villa et al., 2016) and interactions with CLAMP (Albig et al., 2019; Eggers et al., 2023; Tikhonova et al., 2019). Subsequently, the MSL complexes spread from the HAS to active X chromosomal genes through distance chromatin contacts, which are formed by preexisting chromosome architecture (Prayitno et al., 2019; Schauer et al., 2017). roX RNAs guide the MSL complex assembly and its spread from HAS into flanking chromatin (Ankush Jagtap et al., 2019; Ilik et al., 2017; Lv et al., 2019; Müller et al., 2020; Park et al., 2002; Valsecchi et al., 2021; Villa et al., 2021). It was shown that expression of roX RNA on autosome in the presence of MSL proteins excess can induce the local spreading of MSL complex from the site of the transgene localization (Kelley et al., 2008; Oh et al., 2003; Park et al., 2002).
Recent studies have shown that low-complexity C-terminal domains of MSL2 and roX RNAs form a stably condensed state that induces local ‘trapping’ of the MSL complex (Valsecchi et al., 2021). Experimental evidence shows that roX2 RNA plays an instructive role by modulating the MSL complex binding specificity (Villa et al., 2021). For example, Rieder et al., 2019 showed that at the initiation of embryogenesis, the MSL complex first binds to all chromosomes before becoming concentrated on the X chromosome. Weak interactions with chromatin sites on autosomes likely lead to the rapid dissociation of the MSL complex, resulting in the gradual concentration of the components of the complex on HAS. It was shown that complete MSL complexes preferentially interact with the sites on the X chromosome, while the MSL1 and MSL2 proteins can independently bind to many promoters located at all chromosomes (Chlamydas et al., 2016; Hallacli et al., 2012; Straub et al., 2013; Valsecchi et al., 2018).
The MSL1GS is the best for testing the functional role of the MSL1-roX interaction because MSL1Δ1-15 is unstable and MSL1Δ41-85 loses its ability to interact effectively with autosomal promoters. Thus, the N-terminal region between 41 and 85 aa is significant not only for interaction with RNA but also because plays an additional role in attracting MSL1 to active promoters independently from the MSL complex. Previous studies have shown that MSL1 alone or in combination with MSL3 or MOF also associated with many autosomal promoters independent of MSL2 (Chlamydas et al., 2016; Straub et al., 2013). Our ChIP seq data suggest that MSL1 and MSL3 are recruited to the promoters mainly independent of each other. In contrast with the recruitment of the MSL complex to HAS on the X chromosome, the binding of MSL2 to autosomal promoters does not require the MSL1 protein (Valsecchi et al., 2018). Thus, MSL proteins can function autonomously or as a part of alternative transcription complexes.
Our results with the MSL1GS mutant agree with the model that the interaction of MSL1 with roX RNAs is critical for the assembly of the MSL complex and its specific recruitment on the X chromosome. MSL1GS effectively binds to autosomal promoters like MSL1WT but disrupts the recruitment of the MSL complex on the X chromosomal HAS (Figure 8). As a result, males expressing MSL1GS have very low viability, as shown earlier for males that have both roX genes inactivated (roX-). In roX- males, MSL1 and MSL2 bind only to the strongest regions of the X chromosome and begin to bind to autosomal regions (Meller and Rattner, 2002; Park et al., 2003; Stuckenholz et al., 2003). However, overexpression of MSL1 and MSL2 can overcome a lack of roX RNAs (Oh et al., 2003). Consequently, the co-binding of MSL1 and MSL2 to the X chromosome is partially restored, and the viability of males is increased. We also found that the overexpression of MSL2 partially restores the viability of males expressed MSL1GS. Probably, a high concentration of MSL2 increases the likelihood of interaction with MSL1GS without the help of roX RNAs.
The MSL1GS and MSL1Δ41-85 variants can efficiently heterodimerize with MSL1WT, suggesting that efficient MSL1 dimerization does not require roX RNAs. Thus, the interaction of MSL2 with the MSL1 dimer is a critical step in the organization of the MSL complex, which is accelerated in the presence of roX RNA.
The structure of the complex has only been determined for human MSL1 and MSL2 (Hallacli et al., 2012), for which the interaction does not require an RNA mediator (Wu et al., 2011). In contrast, mutagenic analysis showed that in Drosophila, MSL2 preferentially interacts with the MSL1 coiled-coil domain via its RING finger (Copps et al., 1998). Thus, it is possible that interaction between MSL1 and MSL2 in Drosophila is weaker and requires roX RNAs for efficient association. Similar to the MSL1GS mutant, it has also been shown (Hallacli et al., 2012) that a mutant MSL1 dimer that has lost the ability to interact with MSL2 cannot effectively bind to the X chromosome.
Generally, our results are consistent with the model (Valsecchi et al., 2021; Villa et al., 2021) according to which roX RNAs provide efficient assembly of the MSL complex, which has a high affinity for HAS on the X chromosome. This process occurs through specific interactions of roX RNAs with MSL1, MSL2, MLE, and possibly MSL3 and MOF. The inability of MSL1GS to interact with roX RNA disrupts the formation of the MSL complex. Consequently, individual MSL proteins independently bind to regulatory elements (mainly promoters) or integrate into alternative transcription complexes. In general, the functions of roX RNA resemble the activity of the X inactive-specific transcript (Xist) in mammals, which specifically interacts with repression complexes involved in X chromosome inactivation (Brockdorff et al., 2020; Jacobson et al., 2022). In addition, Xist specifically targets repression complexes to the inactive X chromosome. However, for Drosophila, the significance of roX RNAs in specific MSL recruitment to the X chromosome remains an open question requiring further examination.
Materials and methods
Plasmid construction
Request a detailed protocolPCR-directed mutagenesis was used to make constructs with deletion variants of MSL1 (corresponding primers are given in Supplementary file 1). Different full-sized variants of MSL1 were cloned into an expression vector containing attB site for φC31-mediated recombination, Ubi-p63E promoter with its 5’UTR, 3’UTR with SV40 polyadenylation signal, intronless yellow gene as a reporter for detection of transformants (Figure 1—figure supplement 2).
Fly crosses and transgenic lines
Request a detailed protocolDrosophila strains were grown at 25 °C and 75% humidity under standard culture conditions. The transgenic constructs were injected into pre blastoderm embryos using the φC31-mediated site-specific integration system at locus 86Fb (Bischof et al., 2007). The emerging adults were crossed with the y ac w1118 flies, and the progeny carrying the transgene in the 86Fb region were identified by y+ pigmented cuticle.
The msl-1- background corresponds to msl-1γ269/msl-1L60 heterozygote. The msl-1L60/CyO (Chang and Kuroda, 1998) and msl-1γ269cn1bw1/CyO (Palmer et al., 1993) stocks with the null msl-1 mutations were kindly donated by M. Kuroda and J. Lucchesi. The w[1118]; sna[Sco]/CyO, P[ActGFP.w-]CC2 stock no. 9325 was obtained from the Bloomington stock center. To access the viability of males carrying MSL1 variants (MSL1*, MSL1Δ1-15, MSL1Δ41-85, and MSL1Δ1-85), 2–3 days old virgin msl-1L60/CyO, GFP; Ubi-msl-1*/ Ubi-msl-1* females were crossed with 2–3 days old msl-1γ269/CyO, GFP; Ubi-msl-1*/Ubi-msl-1* males. Ubi-msl-1* is any transgene expressing one of the tested MSL1 variants. msl-1L60/msl-1γ269; Ubi-msl-1*/Ubi-msl-1* males and msl-1L60/msl-1γ269; Ubi-msl-1*/Ubi-msl-1* females were counted after hatching. The percentage of male viability with different Ubi-msl-1* transgenes was estimated by taking the viability of females as 100%. The viability of males carrying MSL1 variants (MSL1*, MSL1WT, MSL1Δ66-85 and MSL1Δ24-39) was tested in msl-1L60/msl-1γ269; Ubi-msl-1*/Ubi-msl-1* lines.
To assess the viability of females expressed MSL2 and one of MSL1 variants (MSL*), 2–3 days old virgin Ubi-msl-1*/Ubi-msl-1* females were crossed with 2–3 days old Ubi:msl-2WT-FLAG/TM6,T viability was estimated by taking the viability of Ubi-msl-1*/TM6, Tb females as 100%.
To assess the role of MSL1 variants in the binding of MSL complex to polytene chromosomes, 2–3 days old virgin msl-1L60/msl-1γ269; Ubi-msl-1*/Ubi-msl-1* females were crossed with 2–3 days old msl-1L60/CyO, GFP; Ubi:msl-2WT-FLAG/ TM6, Tb males. msl-1L60/msl-1γ269; Ubi-msl-1*/ Ubi:msl-2WT-FLAG larvae were selected for further study.
b males. Ubi:msl-2WT-FLAG is transgene (86Fb) expressing the MSL2 cDNA under the control of the Ubi-p63E promoter at the 86F8 site (Tikhonova et al., 2019). The percentage of Ubi:msl-2WT-FLAG /Ubi-msl-1* female.
Protein extract preparation
Request a detailed protocolTwenty adult flies were cooled in liquid nitrogen, and homogenized for 30 s with a pestle in 200 μL of extraction buffer (20 mM HEPES, pH 7.5, 100 mM KCl, 5% glycerol, 10 mM EDTA, 1% NP-40, 1% sodium deoxycholate, 0.1% SDS, 1 mM DTT, 5 mM PMSF, and 1:100 Calbiochem Complete Protease Inhibitor Cocktails VII and V) and incubated on ice for 10 min. The suspension was sonicated in a Bioruptor (Diagenode, USA) for 3 min on setting H, 15 s ON/45 s OFF. Then 4 x SDS-PAGE sample buffer was added to the homogenate. Extracts were incubated for 10 min at 100 °C, centrifuged at 16,000 x g for 5 min, and loaded on a 6% SDS-PAGE gel.
Immunostaining of polytene chromosomes
Request a detailed protocolDrosophila hird instar larvae were cultured at 18 °C under standard conditions. Immunoblot analysis showed that the expression of MSL1 and MSL2 under the control of the Ubi-p63E promoter is significantly increased at 18 °C compared to 25 °C (Figure 1—figure supplement 3B). Polytene chromosome staining was performed as described (Murawska and Brehm, 2012). The following primary antibodies were used: rabbit anti-MSL1 at 1:500 dilution, rabbit anti-MSL2 at 1:500 dilution, mouse anti-HA at 1:50 dilution, and mouse anti-FLAG at 1:100 dilution. The secondary antibodies were Alexa Fluor 488 or 555 goat anti-rabbit and anti-mouse, at 1:2000 dilution (Invitrogen). The polytene chromosomes were co-stained with DAPI (AppliChem). Images were acquired on the Nikon Eclipse Ti fluorescent microscope using Nikon DS-Qi2 digital camera, processed with ImageJ 1.54 f and Fiji bundle 2.14.0. 3–4 independent staining and 4–5 samples of polytene chromosomes were performed with each MSL1-expressing transgenic line. Images of stretched chromosomes are obtained as follows: based on the DAPI channel, a chromosome trace is created using the ‘Segmented line’ tool, the thickness is selected (‘Edit-Selection-Properties’), saved as ‘Selection’ in.roi format. This file is opened on top of the desired channel image, and the ‘Edit-Selection-Straighten’ commands are applied.
Co-immunoprecipitation assay
Request a detailed protocolDrosophila S2 cells (obtained from DGRC, Stock 181; https://dgrc.bio.indiana.edu//stock/181) grown in SFX medium (HyClone) at 25 °C were co-transfected by plasmid expressing MSL1*–3xHA and plasmid expressing MSL2-3xFLAG with MACSFectin (Miltenyi Biotec), as recommended by the manufacturer. After transfection, the cells were incubated for 48 hr, washed once with cold 1x PBS, and resuspended in 20 packed cell volumes of hypotonic lysis buffer (20 mM Tris-HCl, pH 7.4, with 10 mM KCl, 10 mM MgCl2, 2 mM EDTA, 10% glycerol, 1% Triton X-100, 1 mM DTT, 0.5 mM PMSF and Halt Protease Inhibitor Cocktail). After incubation on ice for 10 min, the cells were sonicated in a Bioruptor (Diagenode, USA) for 2 min on setting L, 15 s ON/45 s OFF, NaCl was added to a final concentration of 420 mM, and incubation on ice continued for 60 min, with periodic mixing. Sonication was repeated as above to reduce viscosity, cell debris was pelleted by centrifugation at 10,000 g for 30 min at 4 °C, and the supernatant was collected for immunoprecipitation with anti-FLAG, anti-HA, and mouse IgG magnetic beads equilibrated in incubation buffer-150 (20 mM Tris-HCl, pH 7.4, with 150 mM NaCl, 10 mM MgCl2, 1 mM EDTA, 1 mM EGTA, 10% glycerol, and 0.1% NP-40). The protein extract (50 μg protein) was adjusted to a volume of 200 µL with buffer-150, mixed with magnetic beads (15 μL), and incubated on a rotary shaker overnight at 4 °C. The beads were then washed with five portions of buffer-150, resuspended in SDS-PAGE loading buffer, boiled, and analyzed by immunoblotting with anti-FLAG-HRP antibodies (Sigma) antibodies. Proteins were detected using the Crescendo Western Blotting substrate (Millipore).
RNA immunoprecipitation assay
Request a detailed protocolTo prepare the samples, 100 mg of adult flies were ground in a mortar in liquid nitrogen and resuspended in 0.5 mL of VNC-hypotonic lysis buffer (20 mM Tris-HCl, pH 7.4, with 10 mM KCl, 10 mM MgCl2, 10% glycerol, 1% Triton X-100, 1 mM DTT, 0.5 mM PMSF, 2 mM Ribonucleoside Vanadyl Complex, and 1xHalt Protease Inhibitor Cocktail). The suspension was homogenized in a Dounce (loose pestle, 20 strokes) homogenizer and filtered through 100 µm MACS SmartStrainers (Miltenyi Biotec, United States). After incubation on ice for 10 min, the suspension was sonicated in a Bioruptor (Diagenode, USA) for 3 min on setting L (15 s ON/45 sec OFF). NaCl was added to a final concentration of 420 mM, and the suspension was incubated on ice for an additional 60 min with periodic mixing. Sonication was repeated as above to reduce viscosity. Cell debris was pelleted by centrifugation at 10,000 g for 30 min at 4 °C, and the supernatant was collected for subsequent immunoprecipitation. The RNA-protein extracts were diluted to 150 mM of NaCl with VNC-incubation buffer-0 (20 mM Tris-HCl, pH 7.4, 10 mM MgCl2, 10% glycerol, 0.1% NP-40, 0.5 mM PMSF, 2 mM Ribonucleoside Vanadyl Complex, and 1xHalt Protease Inhibitor Cocktail). Rabbit anti-MSL1 (1:200) and non-specific IgG were incubated for 1 hr at room temperature with 20 µL aliquots of Protein A Dynabeads (Thermo Fisher, USA) mixed with 200 µL of PBST. The antibody–Dynabead complexes were then washed and equilibrated in VNC-incubation buffer-150 (20 mM Tris-HCl, pH 7.4, 150 mM NaCl, 10 mM MgCl2, 10% glycerol, 0.1% NP-40, 0.5 mM PMSF, 2 mM Ribonucleoside Vanadyl Complex, and 1xHalt Protease Inhibitor Cocktail). The RNA-protein extracts (10% aliquots were stored as input material for subsequent protein and RNA extractions) were mixed with magnetic beads and incubated on a rotary shaker overnight at 4 °C. The beads were then washed with four portions of VNC-incubation buffer-150 (10% aliquots of supernatants after immunoprecipitation were stored as output material for subsequent protein and RNA extractions). Half of the beads were resuspended in SDS-PAGE loading buffer, boiled, and analyzed by immunoblotting with anti-HA-HRP antibodies (Thermo Fisher, USA). Proteins were detected using the Crescendo Western Blotting substrate (Millipore). The other half of the beads were resuspended in the TRI reagent for subsequent RNA extraction and cDNA synthesis.
RT-PCR
Request a detailed protocolTotal RNA was isolated using the TRI reagent (Molecular Research Center, United States) according to the manufacturer’s instructions. RNA was treated with two units of Turbo DNase I (Ambion) for 30 min at 37 °C to eliminate genomic DNA. The synthesis of cDNA was performed according to the manufacturer’s instructions using 1 μg of RNA or entire sample volume after immunoprecipitation, 100 U of EpiScript reverse transcriptase (LGC Biosearch Technologies, UK), and 10 pM of random hexamers as a primer. The amounts of specific cDNA fragments were quantified by real-time PCR. At least three independent measurements were made for each RNA sample. Relative levels of mRNA expression were calculated in the linear amplification range by calibration to a standard genomic DNA curve to account for differences in primer efficiencies. Individual expression values were normalized with reference to RpL32 mRNA.
ChIP-seq
Request a detailed protocolChromatin was prepared from 2–3-day-old adult flies. Samples of 500 mg each of adult flies were ground in a mortar in liquid nitrogen and resuspended in 10 mL of buffer A (15 mM HEPES-KOH, pH 7.6, 60 mM KCl, 15 mM NaCl, 13 mM EDTA, 0.1 mM EGTA, 0.15 mM spermine, 0.5 mM spermidine, 0.5% NP-40, 0.5 mM DTT) supplemented with 0.5 mM PMSF and Calbiochem Complete Protease Inhibitor Cocktail V. The suspension was then homogenized in a Potter (10 strokes) and subsequently in a Dounce (tight pestle, 15 strokes) homogenizers and filtered through 70 µm MACS SmartStrainers (Miltenyi Biotec, United States). The homogenate was cross-linked with 1% formaldehyde for 15 min at room temperature. Cross-linking was stopped by adding glycine to a final concentration of 125 mM. The nuclei were washed with three 10 mL portions of wash buffer (15 mM HEPES-KOH, pH 7.6, 60 mM KCl, 15 mM NaCl, 1 mM EDTA, 0.1 mM EGTA, 0.1% NP-40, protease inhibitors) and one 5 mL portion of nuclear lysis basic buffer (15 mM HEPES, pH 7.6, 140 mM NaCl, 1 mM EDTA, 0.1 mM EGTA, 1% Triton X-100, 0.5 mM DTT, 0.1% sodium deoxycholate, and protease inhibitors), and resuspended in 1 mL of nuclear lysis buffer (15 mM HEPES, pH 7.6, 140 mM NaCl, 1 mM EDTA, 0.1 mM EGTA, 1% Triton X-100, 0.1% sodium deoxycholate, 0.5% SLS, 0.1% SDS, 0.5 mM PMSF, Calbiochem Complete Protease Inhibitor Cocktail V). The suspension was incubated for 30 min at 4 °C and then sonicated with Covaris ME220 focused-ultrasonicator (50 alternating 15 s ON and 45 s OFF intervals, at 6 °C, peak power 75.0, duty % factor 25, cycles/burst 1000), and 50 µL aliquots were used to test the extent of sonication and to measure DNA concentration. Chromatin was then transferred to 1.5 ml-Eppendorf tubes and centrifuged at 14,000 rpm for 10 min at 4 °C. Finally, supernatants were pooled, aliquoted, and frozen in liquid nitrogen.
For immunoprecipitation chromatin was pre-cleared with Protein A Dynabeads (Thermo Fisher, USA), blocked with 0.5% BSA. Corresponding antibodies were incubated for 1 hr at room temperature with 20 µL aliquots of Protein A (rabbit anti-MSL1, 1:200; anti-MSL2, 1:100; anti-MSL3, 1:200; rabbit non-specific IgG) mixed with 200 µL of PBST. Then antibody–Dynabead complexes were washed and equilibrated in nuclear lysis buffer. Chromatin samples containing 10–20 µg of DNA equivalent in 200 µL nuclear lysis buffer (2 µL aliquots of pre-cleared chromatin as input material) were incubated overnight at 4 °C with antibody–Dynabead complexes. After three rounds of washing with lysis buffer supplemented with 500 mM NaCl and TE buffer (10 mM Tris-HCl, pH 8; 1 mM EDTA), the DNA was eluted with elution buffer (50 mM Tris-HCl, pH 8.0; 1 mM EDTA, 1% SDS), treated with RNase A for 30 min and proteinase K overnight, incubated at 65 °C for 6 hr and extracted with ChIP DNA Clean&Concentrator kit (Zymo Research, USA).
The ChIP-seq libraries were prepared with NEBNext Ultra II DNA Library Prep kit, as described in the manufacturer’s instructions. In short, eluted DNA was end-repaired, and terminal adenosine residues were added using the NEBNext reagents. Indexed adapters were ligated, after which the material was size selected at ~200–600 bp with Ampure XP beads (Beckman Coulter). PCR amplification was performed using NEB primers for 15–16 cycles using the Q5 Hot Start HiFi PCR Master Mix (NEB). The PCR-amplified library was purified using Ampure XP beads and its quality was assessed on a Bioanalyzer 2100 system (Agilent). Diluted libraries were clustered on a pair-read flowcell and sequenced using a NovaSeq 6000 system (Illumina) (using «Genetiсo» facility).
ChIP-seq data processing and sequence analysis
Request a detailed protocolAll ChIP-seq raw data were presented as two biological replicates with paired-end reads (except MSL1 M1[wt]♂), for which only one biological replicate was obtained. Trimming and mapping were performed using cutadapt software (Martin, 2011) and Bowtie version 2 (Langmead and Salzberg, 2012), as described previously (Sabirov et al., 2021). The dm6 version of the Drosophila melanogaster genome was used as a reference genome. After merging replicates, coverage tracks (BedGraph) were obtained using the deepTools (Ramírez et al., 2014) bamCoverage function, with bin widths of 10 bp and extendReads option and normalized by reads per kilobase of the transcript, per million mapped reads (RPKM). Raw and processed data were deposited in the NCBI Gene Expression Omnibus (GEO) under accession number GSE243396.
Peak calling was performed for MSL proteins from flyes of ♂M1[wt] and ♀ M1[wt]+MSL2 lines, as described previously (Sabirov et al., 2021). Briefly:
duplicates were removed with Picard (http://broadinstitute.github.io/picard/);
blacklist filtration was performed;
peak calling was performed using MACS version 2 against nonspecific IgG (Zhang et al., 2008), with a soft p-value threshold of 1×10−2 for further irreproducible discovery rate (IDR) analysis.
Reproducibility was assessed using the IDR pipeline, with p-value thresholds of 0.05 for true replicates and 0.01 for pseudoreplicates. All samples showed ideal reproducibility both the rescue ratio (RR) and self-consistency ratio (SR) were less than 2; thus, an ‘optimal’ set of highly-reproduced peaks was chosen for each sample.
Further analysis was performed in R version 4.2.1 (http://www.r-project.org). Colocalization analysis was performed using the ChIPpeakAnno package (Zhu et al., 2010). Average signal calculation and heatmaps were constructed using the ChIPseeker package (Yu et al., 2015).
HAS regions for which ChIP-seq signal enrichment was calculated were obtained by combining HAS variants from three different studies (Alekseyenko et al., 2008; Ramírez et al., 2015; Straub et al., 2008). Only HAS regions located on the X chromosome were considered in this analysis, resulting in 301 HAS regions (Figure 5—source data 1). For Figure 5—figure supplements 1A and 2A visualization 10% autosomal regions from ♂M1[wt] with the highest log p-value were chosen for every MSL protein (MSL1 - 671, MSL2 - 288, MSL3 - 185 peaks).
To investigate MSL binding for female fly lines at different genomic locations, a list of unique MSL1, MSL2 peaks for ♀ M1[wt]+MSL2 fly line was generated using findOverlapsOfPeaks function (maxgap = 500). This list was subdivided into two groups, according to the mean logFC between protein signal and control at 500 bp region around the center of each peak: ‘MSL1 alone’ (MSL1 >threshold, MSL2
Data availability
Sequencing data have been deposited in GEO under accession code GSE243396.
-
NCBI Gene Expression OmnibusID GSE243396. N-terminus of Drosophila MSL1 is critical for dosage.
References
-
Factor cooperation for chromosome discrimination in DrosophilaNucleic Acids Research 47:1706–1724.https://doi.org/10.1093/nar/gky1238
-
Parallel universes for models of X chromosome dosage compensation in Drosophila: a reviewCytogenetic and Genome Research 148:52–67.https://doi.org/10.1159/000445924
-
One hundred years of gene balance: how stoichiometric issues affect gene expression, genome evolution, and quantitative traitsCytogenetic and Genome Research 161:529–550.https://doi.org/10.1159/000519592
-
Progress toward understanding chromosome silencing by Xist RNAGenes & Development 34:733–744.https://doi.org/10.1101/gad.337196.120
-
Functional interplay between MSL1 and CDK7 controls RNA polymerase II Ser5 phosphorylationNature Structural & Molecular Biology 23:580–589.https://doi.org/10.1038/nsmb.3233
-
Transcriptional control of a whole chromosome: emerging models for dosage compensationNature Structural & Molecular Biology 21:118–125.https://doi.org/10.1038/nsmb.2763
-
A lifelong duty: how Xist maintains the inactive X chromosomeCurrent Opinion in Genetics & Development 75:101927.https://doi.org/10.1016/j.gde.2022.101927
-
Complete dosage compensation in anopheles stephensi and the evolution of sex-biased genes in mosquitoesGenome Biology and Evolution 7:1914–1924.https://doi.org/10.1093/gbe/evv115
-
Structural basis for MOF and MSL3 recruitment into the dosage compensation complex by MSL1Nature Structural & Molecular Biology 18:142–149.https://doi.org/10.1038/nsmb.1960
-
Fast gapped-read alignment with Bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Incorporation of the noncoding roX RNAs alters the chromatin-binding specificity of the Drosophila MSL1/MSL2 complexMolecular and Cellular Biology 28:1252–1264.https://doi.org/10.1128/MCB.00910-07
-
Transcriptional modulation of entire chromosomes: dosage compensationJournal of Genetics 97:357–364.
-
Two-step mechanism for selective incorporation of lncRNA into a chromatin modifierNucleic Acids Research 48:7483–7501.https://doi.org/10.1093/nar/gkaa492
-
Immunostaining of Drosophila polytene chromosomes to investigate recruitment of chromatin-binding proteinsMethods in Molecular Biology 809:267–277.https://doi.org/10.1007/978-1-61779-376-9_18
-
Local spreading of MSL complexes from roX genes on the Drosophila X chromosomeGenes & Development 17:1334–1339.https://doi.org/10.1101/gad.1082003
-
Sex chromosomes evolved from independent ancestral linkage groups in winged insectsMolecular Biology and Evolution 29:1645–1653.https://doi.org/10.1093/molbev/mss010
-
deepTools: a flexible platform for exploring deep-sequencing dataNucleic Acids Research 42:W187–W191.https://doi.org/10.1093/nar/gku365
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
The MSL3 chromodomain directs a key targeting step for dosage compensation of the Drosophila melanogaster X chromosomeNature Structural & Molecular Biology 15:1318–1325.https://doi.org/10.1038/nsmb.1520
-
Dimerization activity of a disordered N-Terminal domain from Drosophila CLAMP proteinInternational Journal of Molecular Sciences 23:3862.https://doi.org/10.3390/ijms23073862
-
Divergent evolution toward sex chromosome-specific gene regulation in DrosophilaGenes & Development 35:1055–1070.https://doi.org/10.1101/gad.348411.121
-
Model-based analysis of ChIP-Seq (MACS)Genome Biology 9:R137.https://doi.org/10.1186/gb-2008-9-9-r137
Article and author information
Author details
Funding
Russian Science Foundation (21-14-00211)
- Oksana Maksimenko
Ministry of Science and Higher Education of the Russian Federation (075-15-2019-1661)
- Oksana Maksimenko
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Mitzi Kuroda (Harvard Medical School), and Erica Larschan (Brown University) for the kindly providing fly strains. We thank Farhod Hasanov and Aleksander Parshikov for fly injection. This study was performed using the equipment of the Center for Precision Genome Editing and Genetic Technologies for Biomedicine of IGB RAS supported by the Ministry of Science and Education of the Russian Federation. Some cartoons (Figures 7 and 8) were created with BioRender.com. This work was supported by the grant 21-14-00211 from Russian Science Foundation. Chip-Seq analysis was supported by grant 075-15-2019-1661 from the Ministry of Science and Higher Education of the Russian Federation.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.93241. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Babosha et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 392
- views
-
- 26
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Chromosomes and Gene Expression
- Evolutionary Biology
Gene regulation is essential for life and controlled by regulatory DNA. Mutations can modify the activity of regulatory DNA, and also create new regulatory DNA, a process called regulatory emergence. Non-regulatory and regulatory DNA contain motifs to which transcription factors may bind. In prokaryotes, gene expression requires a stretch of DNA called a promoter, which contains two motifs called –10 and –35 boxes. However, these motifs may occur in both promoters and non-promoter DNA in multiple copies. They have been implicated in some studies to improve promoter activity, and in others to repress it. Here, we ask whether the presence of such motifs in different genetic sequences influences promoter evolution and emergence. To understand whether and how promoter motifs influence promoter emergence and evolution, we start from 50 ‘promoter islands’, DNA sequences enriched with –10 and –35 boxes. We mutagenize these starting ‘parent’ sequences, and measure gene expression driven by 240,000 of the resulting mutants. We find that the probability that mutations create an active promoter varies more than 200-fold, and is not correlated with the number of promoter motifs. For parent sequences without promoter activity, mutations created over 1500 new –10 and –35 boxes at unique positions in the library, but only ~0.3% of these resulted in de-novo promoter activity. Only ~13% of all –10 and –35 boxes contribute to de-novo promoter activity. For parent sequences with promoter activity, mutations created new –10 and –35 boxes in 11 specific positions that partially overlap with preexisting ones to modulate expression. We also find that –10 and –35 boxes do not repress promoter activity. Overall, our work demonstrates how promoter motifs influence promoter emergence and evolution. It has implications for predicting and understanding regulatory evolution, de novo genes, and phenotypic evolution.
-
- Chromosomes and Gene Expression
- Genetics and Genomics
The enhancer-promoter looping model, in which enhancers activate their target genes via physical contact, has long dominated the field of gene regulation. However, the ubiquity of this model has been questioned due to evidence of alternative mechanisms and the lack of its systematic validation, primarily owing to the absence of suitable experimental techniques. In this study, we present a new MNase-based proximity ligation method called MChIP-C, allowing for the measurement of protein-mediated chromatin interactions at single-nucleosome resolution on a genome-wide scale. By applying MChIP-C to study H3K4me3 promoter-centered interactions in K562 cells, we found that it had greatly improved resolution and sensitivity compared to restriction endonuclease-based C-methods. This allowed us to identify EP300 histone acetyltransferase and the SWI/SNF remodeling complex as potential candidates for establishing and/or maintaining enhancer-promoter interactions. Finally, leveraging data from published CRISPRi screens, we found that most functionally verified enhancers do physically interact with their cognate promoters, supporting the enhancer-promoter looping model.