A SUMO E3 ligase promotes long non-coding RNA transcription to regulate small RNA-directed DNA elimination
Abstract
Small RNAs target their complementary chromatin regions for gene silencing through nascent long non-coding RNAs (lncRNAs). In the ciliated protozoan Tetrahymena, the interaction between Piwi-associated small RNAs (scnRNAs) and the nascent lncRNA transcripts from the somatic genome has been proposed to induce target-directed small RNA degradation (TDSD), and scnRNAs not targeted for TDSD later target the germline-limited sequences for programmed DNA elimination. In this study, we show that the SUMO E3 ligase Ema2 is required for the accumulation of lncRNAs from the somatic genome and thus for TDSD and completing DNA elimination to make viable sexual progeny. Ema2 interacts with the SUMO E2 conjugating enzyme Ubc9 and enhances SUMOylation of the transcription regulator Spt6. We further show that Ema2 promotes the association of Spt6 and RNA polymerase II with chromatin. These results suggest that Ema2-directed SUMOylation actively promotes lncRNA transcription, which is a prerequisite for communication between the genome and small RNAs.
Editor's evaluation
This important study demonstrates that protein SUMOylation is essential for programmed DNA elimination guided by small RNAs during conjugation in Tetrahymena ciliates. The authors present convincing evidence that the E3 SUMO ligase Ema2 is necessary for the production of long non-coding RNAs from the somatic nucleus, targeted small RNA degradation, and DNA elimination. The authors also show that the transcription regulator Spt6 is a SUMOylation target of Ema2, though the relevance of this is not completely established. This paper is of broad significance and will appeal to those interested in non-coding RNA biology, the control of programmed genome rearrangements, or ciliate biology.
https://doi.org/10.7554/eLife.95337.sa0Introduction
Small RNAs of approximately 20–30 nucleotides that are complexed with Argonaute family proteins target either mRNAs for post-transcriptional silencing or chromatin regions for transcriptional gene silencing (Holoch and Moazed, 2015; Wilson and Doudna, 2013). For the latter process, small RNAs are generally considered to recognize their genomic targets via nascent lncRNAs to induce heterochromatin formation. Therefore, chromatin regions that are targeted for silencing by small RNAs must be paradoxically transcribed to provide nascent lncRNAs.
In fission yeast, small interfering RNAs (siRNAs) mediate the deposition of histone 3 lysine 9 di- and trimethylation (H3K9me2/3) for heterochromatin assembly at centromeric repeats (Hall et al., 2002; Volpe et al., 2002). While the HP1 protein Swi6 binds to H3K9me2/3 for transcriptional silencing, phosphorylation of histone H3 serine 10 at the M phase of the cell cycle evicts Swi6, passively allowing lncRNA transcription from the centromeric repeats at G1 and S phase and thus promoting H3K9me2/3 deposition to newly assembled nucleosomes (Chen et al., 2008; Kloc et al., 2008). In contrast, the HP1 paralog Rhino in fruit flies, which is specifically enriched at PIWI-interacting RNA (piRNA) clusters by binding to H3K9me2/3 (Le Thomas et al., 2014), actively recruits dedicated variants of basal transcription factors, allowing lncRNA transcription from heterochromatin while preventing mRNA transcription from transposons in the same loci (Andersen et al., 2017). A similar active heterochromatin-dependent lncRNA transcription has been reported in plants, where SHH1, a reader of H3K9me2/3 as well as mono-methylated H3K9, recruits the plant-specific RNA polymerase IV (Law et al., 2013; Law et al., 2011).
Because small RNA-producing loci are also small RNA targets in most of the studied small RNA-directed heterochromatin formation processes, it poses a challenge to separately investigate lncRNA transcription for small RNA biogenesis and that for small RNA-dependent recruitment of downstream effectors in these processes. In contrast, the source and target loci of small RNAs reside in different nuclei during programmed DNA elimination in some ciliated protozoans such as Tetrahymena thermophila and Paramecium tetraurelia, which provides a unique system to study the mechanism and roles of lncRNA transcription in small RNA-directed chromatin regulation.
In most ciliates, each cell contains two types of nuclei, the diploid germline micronucleus (MIC) and the polyploid somatic macronucleus (MAC). During conjugation, a sexual reproduction process in ciliates, the MIC undergoes meiosis and fertilization, followed by the formation of the new MIC and MAC, while the parental MAC is degraded (Chalker et al., 2013). Programmed DNA elimination occurs in the new MAC in most ciliates. In Tetrahymena, this process downsizes the 200 Mb MIC genome to the 103 Mb MAC genome by removing ~12,000 internal elimination sequences (IESs), many of which include transposons, followed by re-ligation of the remaining macronuclear-destined sequences (MDSs) forming the MAC chromosomes (Noto and Mochizuki, 2018).
DNA elimination in Tetrahymena and Paramecium is regulated by three types of lncRNAs, which occur in different nuclei at distinct times (Aronica et al., 2008; Schoeberl and Mochizuki, 2011). During meiotic prophase in Tetrahymena, lncRNAs are transcribed bidirectionally in the MIC in a genome-wide manner by RNA polymerase II and dedicated conjugation-specific Mediator-associated proteins (Chalker and Yao, 2001; Garg et al., 2019; Mochizuki and Gorovsky, 2004a; Schoeberl et al., 2012; Tian et al., 2019). The sexual reproduction-specific Spt4 and Spt5 paralogs are specifically involved in the MIC transcription in Paramecium (Gruchota et al., 2017; Owsian et al., 2022), while Tetrahymena genome does not encode such specialized Spt4/Spt5 paralogs. These micronuclear long non-coding RNAs (MIC-lncRNAs) are processed into small RNAs (~29 nt in Tetrahymena and 25-nt in Paramecium), called scnRNAs, by Dicer homologs (Lepère et al., 2009; Malone et al., 2005; Mochizuki and Gorovsky, 2005) and loaded into Piwi-clade Argonaute proteins (Bouhouche et al., 2011; Mochizuki et al., 2002; Noto et al., 2010).
The Piwi-scnRNA complex then moves into the parental MAC, where parental macronuclear long non-coding RNAs (pMAC-lncRNAs) are transcribed bidirectionally during the mid-conjugation stages (Woo et al., 2016). The RNA helicase Ema1 in Tetrahymena promotes the interaction between pMAC-lncRNAs and scnRNAs (Aronica et al., 2008). Reminiscently of target-directed micro RNA degradation (TDMD) (Han and Mendell, 2023), this interaction induces TDSD, leading to the selective retention of IES-derived scnRNAs (Aronica et al., 2008; Mochizuki and Gorovsky, 2004b; Noto and Mochizuki, 2018; Noto et al., 2015; Schoeberl et al., 2012). A similar TDSD has also been suggested in Paramecium, and the importance of pMAC-lncRNAs for DNA elimination was demonstrated by disrupting certain pMAC-lncRNAs by RNAi in this ciliate (Lepère et al., 2008). Although mRNAs are transcribed in the parental MAC, it remains unclear whether they also can induce TDSD and how mRNAs and pMAC-lncRNAs can be transcribed from overlapping locations. Also, although trimethylation of histone 3 lysine 27 (H3K27me3) occurs in the parental MAC in a scnRNA- and Polycomb repressive complex 2 (PRC2)-dependent manner in both Tetrahymena and Paramecium (Lhuillier-Akakpo et al., 2014; Liu et al., 2007), the role of H3K27me3 in the parental MAC, if any, is unclear.
When the new MAC develops, the retained Piwi-scnRNA complexes translocate to the new MAC, where yet another type of lncRNA, new macronuclear non-coding RNAs (nMAC-lncRNAs), is transcribed. The interaction between nMAC-lncRNA and scnRNA, which is also dependent on Ema1 (Aronica et al., 2008), is believed to recruit PRC2, which catalyzes both H3K9me2/3 and H3K27me3 for IES-specific heterochromatin assembly (Frapporti et al., 2019; Liu et al., 2007; Miró-Pina et al., 2022; Wang et al., 2022; Xu et al., 2021), and facilitate the secondary production of scnRNAs from nMAC-lncRNAs, which further promotes heterochromatin assembly (Allen et al., 2017; Noto et al., 2015). Then, IESs are excised by domesticated PiggyBac transposases (Baudry et al., 2009; Bischerour et al., 2018; Cheng et al., 2010; Vogt and Mochizuki, 2013). While nMAC-lncRNA production occurs prior to the excision of IESs in Tetrahymena (Mutazono et al., 2019), it also occurs from excised IESs that are concatenated and circularized in Paramecium (Allen et al., 2017).
Among the above three lncRNAs, pMAC-lncRNA does not produce small RNAs and is believed to be specialized for receptor function. Therefore, pMAC-lncRNA transcription and the following TDSD in the parental MAC provide a unique paradigm to investigate how the genome is transcribed to communicate with small RNAs. In this study, we show that Ema2, the conjugation-specific E3 ligase for a small ubiquitin-like modifier (SUMO), is required for pMAC-lncRNA transcription in Tetrahymena and thus provides a tool to dissect the role of and the molecular mechanism for the transcription of this lncRNA.
Results
Ema2 is exclusively expressed during conjugation and localized in the MAC
As part of our systematic investigation into genes highly upregulated during conjugation (Loidl, 2021), we explored the function of EMA2 (TTHERM_00113330). EMA2 mRNA is exclusively expressed during conjugation (Figure 1A). The encoded Ema2 protein tagged with HA at the endogenous EMA2 locus was not detectable in the vegetative cells (Figure 1B, Vg) and first appeared in the MAC during conjugation (Figure 1B, 3 hr post-induction of mating [hpm] and 6 hpm). At the onset of new MAC development, Ema2 disappeared from the parental MAC and appeared in the new MAC (Figure 1B, 8 hpm), which later faded away (Figures 1B, 12, and 14 hpm). The conjugation-specific expression and the localization switch from the parental to the new MAC are reminiscent of the factors involved in DNA elimination such as the Piwi protein Twi1, which is loaded by scnRNAs, and PRC2 (Liu et al., 2007; Mochizuki et al., 2002; Noto et al., 2010).
Ema2 is required for completing DNA elimination
DNA elimination of exconjugants (sexual progeny) at 36 hpm was analyzed by DNA fluorescent in situ hybridization (FISH) using probes complementary to the transposable element Tlr1 (Wuitschick et al., 2002). DNA elimination is completed by ~14–18 hpm in wild-type cells (Austerberry et al., 1984; Mutazono et al., 2019), and the Tlr1 element was detected only in the MICs in the exconjugants from wild-type cells (Figure 2A, WT). In contrast, the Tlr1 element was detected in both the MICs and the MACs in the exconjugants from the EMA2 somatic KO strains, in which all EMA2 copies in the MAC were disrupted (Shehzada and Mochizuki, 2022; Figure 2A, KO). The intensity of the FISH signal in the new MACs was lower in the exconjugants from the EMA2 KO cells than in those from the TWI1 KO cells (Figure 2B), in the latter of which DNA elimination is known to be completely blocked (Noto et al., 2015). Therefore, DNA elimination was partially blocked in the exconjugants of EMA2 KO cells. Consistent with the requirement of DNA elimination in the viability of sexual progeny (Cheng et al., 2010; Vogt and Mochizuki, 2013), EMA2 KO cells did not produce viable progeny (Figure 2C). Altogether, we conclude that maternally expressed EMA2 is required for completing DNA elimination.
Ema2 is required for target-directed small RNA degradation (TDSD)
We next asked whether EMA2 is involved in TDSD. The production of scnRNAs occurs from IESs and their surrounding MDS regions in the MIC at the early conjugation stages (~2–3 hpm), and scnRNAs complementary to the MAC genome (=MDSs) are subjected to TDSD in the parental MAC in the mid-conjugation stages (~3–7 hpm), resulting in the selective retention of IES-derived scnRNAs that later target IESs for DNA elimination in the new MAC (Aronica et al., 2008; Mochizuki and Gorovsky, 2004b; Noto et al., 2015; Schoeberl et al., 2012). We first compared scnRNAs at different conjugation stages by northern blot analysis using the 50-nt Mi-9 probe that is complementary to a repetitive sequence found in both IESs and MDSs (Aronica et al., 2008). Because scnRNAs detected by this probe are complementary to the MAC genome, in wild-type cells, the Mi-9-complementary scnRNAs were detected at 3 hpm, reduced at 4.5 hpm, and became undetectable at 6 hpm due to TDSD (Figure 3A, WT). In contrast, the Mi-9-complementary scnRNAs remained detectable at 6 hpm in the EMA2 KO cells (Figure 3A, KO), suggesting that Ema2 is required for TDSD.
We also analyzed TDSD by small RNA sequencing at 3 and 8 hpm. Sequenced 26- to 32-nt small RNAs, which correspond to scnRNAs (Schoeberl et al., 2012), were mapped to the genomic tiles of MDSs and IESs (see Materials and methods). Although scnRNAs that mapped many of the MDS tiles were greatly reduced from 3 to 8 hpm due to TDSD in the wild-type cells (Figure 3B, left, WT), those in the EMA2 KO cells were only slightly reduced by 8 hpm (Figure 3B, left, EMA2 KO). Because TDSD takes place concurrently with the scnRNA production (Schoeberl et al., 2012), the increased abundance of MDS-complementary scnRNAs at 3 hpm in the EMA2 KO cells compared to the wild-type cells can also be attributed to the necessity of Ema2 in TDSD. This TDSD defect in the EMA2 KO cells seems milder than that in the EMA1 KO cells, in which the amount of MDS-complementary scnRNAs remained constant until 8 hpm (Figure 3B, left, EMA1 KO), indicating that Ema1 and Ema2 act differently in TDSD. In contrast, the amounts of scnRNAs complementary to IES tiles at 3 hpm and 8 hpm were comparable in all strains (Figure 3B, right), which is consistent with previous observations that only scnRNAs complementary to MDSs are targeted for TDSD. These results suggest that Ema2 is required for TDSD at the genome-wide level. Because it is known that loss of TDSD results in a partial block of DNA elimination (Aronica et al., 2008), the mild DNA elimination defect in EMA2 KO cells (Figure 2B) can be explained by the requirement of Ema2 in TDSD.
Ema2 is required for the accumulation of SUMOylated proteins during conjugation
Ema2 possesses an SP-RING domain that has been found in many SUMO E3 ligases (Hochstrasser, 2001). The SP-RING domain of Ema2 is atypical (Figure 4A) in that the first cysteine of the zinc ion-binding residues in typical counterparts is replaced by histidine and some of the stabilizer residues (Duan et al., 2009; Yunus and Lima, 2009) are likely missing. SUMOylation is catalyzed by the sequential actions of E1, E2, and in most cases, E3 enzymes and thus Ema2 should interact with the E2 enzyme if it acts in SUMOylation. We indeed found that Ema2 can directly interact with the Tetrahymena SUMO E2 enzyme Ubc9 in vitro (Figure 4B), suggesting that Ema2 is a bona fide SUMO E3 ligase.
To examine the role of Ema2 in SUMOylation, we expressed HA-tagged Smt3 (HA-Smt3) in an EMA2 KO strain. SUMO is solely encoded by SMT3 in Tetrahymena (Nasir et al., 2015), and HA-Smt3 could replace the essential function of endogenous Smt3 (Figure 4—figure supplement 1). We then crossed this strain with either a wild-type strain (called the WT cross) or another EMA2 KO strain (called the EMA2 KO cross). Because proteins and mRNAs are exchanged between two mating pairs through the conjugation junction, EMA2 mRNA/Ema2 protein expressed in the wild-type partner of the WT cross is expected to move into the EMA2 KO partner and restore the EMA2 KO phenotypes.
Total proteins were harvested at 4.5 and 6 hpm, and SUMOylated proteins were detected by western blotting using an anti-HA antibody (Figure 4C). At both time points, high molecular weight proteins (mainly >200 kDa) were detected in the WT cross (Figure 4C, WT), and they were reduced to ~50% in the EMA2 KO cross (Figure 4C, KO). These results indicate that Ema2 is the major SUMO E3 ligase during the mid-conjugation stages. The remaining Ema2-independent SUMOylation is likely mediated by other SUMO E3 ligases (including the SP-RING containing proteins TTHERM_00227730, TTHERM_00442270 and TTHERM_00348490), and/or E3-independent SUMOylation (Sampson et al., 2001). The requirement of protein SUMOylation in DNA elimination was previously demonstrated in Paramecium by RNAi knockdown of UBA2, the gene encoding the SUMO E1 enzyme, and SUMO (Matsuda and Forney, 2006). Therefore, the involvement of a SUMO pathway in DNA elimination is likely conserved among ciliates.
Ema2 is required for SUMOylation of Spt6
Next, to identify the SUMOylation target(s) of Ema2, we introduced a construct expressing His-tagged Smt3, which can also replace the essential function of Smt3 (Figure 4—figure supplement 1 ), into an EMA2 KO strain and crossed it with a wild-type strain (WT cross) or another EMA2 KO strain (KO cross). Then, SUMOylated proteins at 6 hpm were concentrated using nickel-NTA beads and identified by mass spectrometry. We additionally examined the mating of wild-type strains without His-Smt3 expression and excluded any proteins identified with a log2 label-free quantification (LFQ) score above 25 or those possessing more than six consecutive histidine residues in this analysis, considering them as proteins binding to the nickel-NTA beads without His-Smt3 conjugation. Although most of the proteins were detected similarly between the WT cross and EMA2 KO cross, Spt6, the most abundantly detected protein in the WT cross, was detected at a much lower level in the EMA2 KO cross (Figure 5A). This result suggests that Spt6 is the major Ema2 target for SUMOylation.
To confirm the Ema2-dependent SUMOylation of Spt6, we introduced a construct expressing HA-tagged Spt6 (Spt6-HA) from the endogenous SPT6 locus into an EMA2 KO strain and mated it with a wild-type strain (WT cross) or another EMA2 KO strain (EMA2 KO cross). Then, total proteins harvested at 4.5 and 6 hpm were analyzed by western blotting using an anti-HA antibody (Figure 5B). In the WT cross, a slower migrating population of Spt6-HA was detected in addition to a band corresponding to unmodified Spt6-HA at both time points (Figure 5B, WT). In contrast, slower migrating Spt6-HA was barely detectable in the EMA2 KO cross (Figure 5B, KO). Then, to examine the timing of the appearance of the slower migrating Spt6 species, we introduced the same Spt6-HA-expressing construct into a wild-type strain and Spt6-HA was analyzed by western blotting (Figure 5—figure supplement 1). Consistent with the Ema2-dependent appearance of the slower migrating Spt6-HA, they were not detected in growing and starved vegetative wild-type cells (Figure 5—figure supplement 1, Veg and 0 hpm, respectively) when Ema2 was not expressed (Figure 1). The slower migrating Spt6-HA was also detected at 8 hpm when the new MAC was already formed (Figure 5—figure supplement 1, 8 hpm) suggesting that Spt6 is possibly SUMOylated also in the new MAC.
The nature of slower migrating species of Spt6-HA was further examined by immunoprecipitating Spt6-HA using an anti-HA antibody. Among the total purified Spt6-HA (Figure 5C, left, IP-WT), the slower migrating species were detected by an anti-Smt3 antibody in the WT cross (Figure 5C, right, IP-WT), and such SUMOylated Spt6-HA species were greatly reduced in the EMA2 KO cross (Figure 5C, KO). These results indicate that the slower migrating species of Spt6 are SUMOylated Spt6 and Ema2 is required for the majority of SUMOylation of Spt6 during the mid-conjugation stages. The remaining SUMOylation observed on Spt6 in the absence of Ema2 is likely facilitated by other SUMO E3 ligases and/or E3-independent SUMOylation, as discussed earlier for the other instances of Ema2-independent SUMOylations.
Ema2 is required for the accumulation of lncRNA in the parental MAC
Spt6 is a conserved regulator of several steps of transcription in various eukaryotes. Because TDSD was proposed to be triggered by the base-pairing interaction between scnRNAs and nascent pMAC-lncRNAs in the parental MAC (Aronica et al., 2008; Noto and Mochizuki, 2018), we hypothesized that Ema2-dependent Spt6 SUMOylation promotes pMAC-lncRNA transcription.
To examine pMAC-lncRNAs, we amplified transcripts spanning IES-MDS borders by RT‒PCR in which MAC-lncRNAs can be distinguished from MIC-lncRNAs by their lengths (Figure 6A). For the three loci examined, pMAC-lncRNAs were detected in wild-type cells but not in EMA2 KO cells at 6 hpm, although the control RPL21 mRNA was detected in both conditions (Figure 6A). The lack of detection of pMAC-lncRNAs in EMA2 KO cells was not due to loss of the primer binding sites by alternative DNA elimination in their prior sexual reproductions, as genomic PCR with the same primer sets detected the corresponding MAC loci in the EMA2 KO cells (Figure 6—figure supplement 1). We, therefore, conclude that Ema2 is required for the accumulation of pMAC-lncRNAs at least for the three tested loci.
The accumulation of lncRNAs was also examined cytologically. As the consequence of the bidirectional lncRNA transcription in the three different nuclei, double-stranded RNAs are accumulated in these nuclei (Aronica et al., 2008; Malone et al., 2005; Mochizuki and Gorovsky, 2005; Woo et al., 2016). As previously reported, immunostaining of wild-type cells using the long double-stranded (ds)RNA-specific J2 antibody (Schönborn et al., 1991; Figure 6B, WT) detected long dsRNAs first in the MIC at its premeiotic stage (3 hpm), then in the parental MAC at the mid stages of conjugation (6 hpm), and finally in the new MAC at late conjugation (9 hpm). In the EMA2 KO cells, long dsRNAs were detected in the MIC at 3 hpm and in the new MAC at 9 hpm but were undetectable in the parental MAC at 6 hpm (Figure 6B, EMA2 KO), indicating that Ema2 is specifically required for the accumulation of pMAC-lncRNAs.
In contrast, EMA1 was dispensable for the accumulation of dsRNAs in the parental MAC (Figure 6B, EMA1 KO), which is consistent with the previous notions that Ema1 is required for pMAC-lncRNAs to interact with scnRNAs but not their transcription (Aronica et al., 2008). Although it is unclear whether lncRNAs are single or double-stranded when Ema1 promotes the lncRNA-scnRNAs interaction, the less severe TDSD defect observed in the EMA2 KO cells compared to the EMA1 KO cells (Figure 3B) indicates that certain Ema1-dependent TDSD may be initiated by single-stranded lncRNAs or mRNAs that are transcribed independently of Ema2.
In parallel to TDSD, scnRNA-dependent accumulation of H3K27me3, which is catalyzed by the histone methyltransferase Ezl1 within PRC2, occurs in the parental MAC (Liu et al., 2007). We found that H3K27me3 was undetectable in the parental MAC in the EMA2 KO cells at 6 hpm or in the EMA1, TWI1, and EZL1 KO cells (Figure 6C). Therefore, Ema2 is required for the scnRNA-directed deposition of H3K27me3 in the parental MAC, which can be explained by the loss of the chromatin recruitment of PRC2 by the scnRNA-lncRNA interaction. Altogether, these results indicate that Ema2 is required for the genome-wide accumulation of pMAC-lncRNAs.
Ema2 facilitates the chromatin association of Spt6 and RNA polymerase II
Next, we analyzed the localization of Spt6 by immunofluorescence staining. In both WT and EMA2 KO cross described above, Spt6-HA was similarly detected in the MAC at 4.5 hpm (Figure 7A, left) and in the new MAC at 9 hpm (Figure 7—figure supplement 1). We also found that the MAC localization of Rpb3, the third largest subunit of RNA polymerase II (RNAPII) was not affected in the absence of Ema2 (Figure 7A, right). Therefore, at least at the cytological level, the absence of Ema2 does not affect the localization of Spt6 and RNAPII.
We then examined their subnuclear localization by cell fractionation (Figure 7B, left), in which cells at 4.5 hpm were first incubated with lysis buffer to extract soluble cyto- and nucleoplasmic proteins, and then, the insoluble fraction was sonicated to extract chromatin-bound proteins (Ali et al., 2018). In the WT cross (Figure 7B, WT), Spt6-HA was detected in both the soluble (S1) and chromatin-bound (S2) fractions, but slower-migrating SUMOylated Spt6-HA was detected only in the latter fraction, indicating that SUMOylated Spt6 is associated with chromatin. A similar distribution was detected for Rpb3. This chromatin association of Spt6 and RNAPII does not require RNA, as they were retained in the insoluble fraction after RNase A treatment without sonication (Figure 7C), although Twi1, which interacts with chromatin through nascent pMAC-lncRNA (Aronica et al., 2008), was mostly detected in the soluble fraction after the treatment (Figure 7C). In the EMA2 KO cross, Spt6-HA and Rpb3 in the chromatin-bound fraction were greatly reduced (Figure 7B, KO, S2). These findings were further validated in two additional replicated experiments (Figure 7—figure supplements 2 and 3). We, therefore, conclude that Ema2 promotes the association of Spt6 and RNAPII with chromatin in the parental MAC at the mid-conjugation stage.
Spt6 is not necessary to be SUMOylated in Ema2-directed lncRNA transcription
Next, to directly examine the importance of SUMOylation of Spt6, we intended to produce a SUMOylation-defective Spt6 mutant. Tetrahymena Spt6 consists of an N-terminal domain of low complexity (DOL) and a downstream region that is conserved among eukaryotes (Figure 8A, WT). We produced constructs expressing HA-tagged Spt6 with lysine to arginine substitutions either for all lysine residues in the DOL (HA-SPT6-DOL-KR) or for lysine residues that were on the surface of an AlphaFold predicted Tetrahymena Spt6 structure in one of the three non-DOL parts (N, M, and C in Figure 8A) in addition to those that were also predicted to be SUMOylation targets by the JASSA algorithm (Beauclair et al., 2015) in all three parts (HA-SPT6-N-KR, HA-SPT6-M-KR, and HA-SPT6-C-KR, Figure 8A).
We introduced these mutant constructs as well as the construct expressing wild-type Spt6 (HA-SPT6-WT) into SPT6 KO cells and found that all of them were able to restore the lethality of SPT6 KO cells (Figure 8—figure supplement 1). However, the cells rescued by HA-SPT6-N-KR and HA-SPT6-M-KR showed severe defects in meiotic progression and mating initiation, respectively, making their SUMOylation status during conjugation uninvestigable. In contrast, the cells rescued by HA-SPT6-DOL-KR and HA-SPT6-C-KR showed a seemingly normal progression of conjugation, although the HA-SPT6-C-KR cells showed lower mating efficiency. We found that while HA-Spt6-DOL-KR was SUMOylated like HA-Spt6-WT, SUMOylation was not detected on HA-Spt6-C-KR (Figure 8B). Therefore, Spt6-C-KR represents a SUMOylation-defective Spt6 mutant, exhibiting at least a reduced level of SUMOylation compared to Spt6 in the absence of Ema2 (compare Figure 8B and Figure 5B).
Cell fractionation experiments showed that HA-Spt6 and RNAPII were detected in the chromatin-bound fraction in HA-SPT6-C-KR cells as in the HA-SPT6-WT cells consistently in three independent experiments (Figure 8C, Figure 8—figure supplement 2), and long dsRNAs were detected in the parental MAC (Figure 8D) as well as in the new MAC (Figure 8—figure supplement 3) in the HA-SPT6-C-KR cells. These results suggest that, contrary to our expectation, Spt6 SUMOylation per se is not required for the lncRNA transcription in the parental MAC. Nonetheless, exconjugants from the HA-SPT6-C-KR cells showed a mild DNA elimination defect (Figure 8E), indicating that Ema2-directed Spt6-SUMOylation plays a role in DNA elimination other than promoting pMAC-lncRNA transcription.
Discussion
In this study, we showed that the conjugation-specific SUMO E3 ligase Ema2 is required for the accumulation of lncRNAs (Figure 6), TDSD (Figure 3A and B), and heterochromatin formation (Figure 3C) in the parental MAC and eventually for completing DNA elimination (Figure 1) in Tetrahymena. We found that Ema2 is responsible for SUMOylation of the transcriptional regulator Spt6 (Figure 5) and promotes the interaction of Spt6 and RNAPII with chromatin (Figure 7). We, therefore, conclude that Ema2 facilitates genome-wide lncRNA transcription in the parental MAC, which is a prerequisite for scnRNA-chromatin communication and thus for the downstream TDSD that regulates DNA elimination.
Our observation that Ema2 enhances the chromatin interaction of Spt6 and RNAPII (Figure 7) seems to contradict the fact that Spt6 and RNAPII are essential for vegetative cell viability (Figure 8—figure supplement 1; Mochizuki and Gorovsky, 2004a), but Ema2 is expressed exclusively during conjugation (Figure 2). Moreover, as EMA2 KO cells did not significantly impede the progression of conjugation processes, any essential mRNA transcriptions for these processes must take place in the parental MAC during conjugation even in the absence of Ema2. Therefore, the observed loss of the majority of Spt6 and RNAPII from chromatin in the absence of Ema2 (Figure 7B) must be a temporal event during the mid-conjugation stage. This suggests that Spt6 and RNAPII might be specifically engaged in pMAC-lncRNA transcription at this particular time window in wild-type cells. It is likely that some radical change in the chromatin environment in the parental MAC is required for the global transcription of pMAC-lncRNAs, and Ema2-dependent SUMOylation might maintain Spt6 and RNAPII on chromatin for transcription in such an environment. Alternatively, because the expression of many of the transcriptional machineries, including Spt6 and RNAPII subunits, is highly upregulated during conjugation (Miao et al., 2009; Mochizuki and Gorovsky, 2004a; Xiong et al., 2012), the transcription of pMAC-lncRNAs might require a greater amount of transcriptional machinery than other types of transcription, and Ema2-dependent SUMOylation might be required for committing these excess machineries to chromatin for transcription.
Our investigation of the SUMOylation-defective Spt6 mutant suggested that although Spt6 is SUMOylated in an Ema2-dependent manner (Figure 5), SUMOylation of Spt6 by itself is not required for pMAC-lncRNA transcription (Figure 8). This observation suggests that Ema2 may facilitate pMAC-lncRNA transcription through SUMOylation of some protein(s) other than Spt6. Although, we failed to detect any other protein that was SUMOylated in an Ema2-dependent manner (Figure 5A), some Ema2-dependent SUMOylation events might not be detected in our current approach if Ema2-dependent and Ema2-independent SUMOylation occur in the same proteins and/or if Ema2-dependent SUMOylation results in protein degradation, such as through SUMO-targeted ubiquitination (Praefcke et al., 2012; Staudinger, 2017). Alternatively, SUMOylation of Spt6 may promote pMAC-lncRNA transcription by competing with some other modification that occurs at the same lysine residue(s) and our SUMOylation-defective Spt6 mutant might actually mimic the SUMOylated state of Spt6 by preventing such modification. To further investigate the role of Ema2, future studies will need to identify individual Ema2-dependent SUMOylated lysine residues under conditions of in vivo proteasome inhibition and to compare the status of other modifications of Spt6 in the presence and absence of Ema2. These are also important to understand the role of Spt6 SUMOylation in DNA elimination (Figure 8E).
It was reported that SUMOylation of the heterochromatin components Swi6, Chp2, and Clr4 is important for RNAi-directed heterochromatin silencing in fission yeast (Shin et al., 2005). Additionally, SUMOylation of the histone deacetylase HDAC1 promotes small RNA-directed transcriptional silencing in C. elegans (Kim et al., 2021). Furthermore, in the piRNA-directed transposon silencing of Drosophila, the SUMO E3 ligase Su(var)2–10 plays an essential role in the recruitment of the histone methyltransferase complex SetDB1, which deposits H3K9me3 (Ninova et al., 2020a; Ninova et al., 2020b), and Su(var)2–10-independent SUMOylation of Panoramix (Panx) also promotes transcriptional silencing by recruiting the general heterochromatin effector Sov (Andreev et al., 2022). All these SUMOylation events seem to require pre-existing lncRNA transcription. In contrast, this study showed that Ema2-dependent SUMOylation during the process of DNA elimination in Tetrahymena occurs upstream of lncRNA transcription (Figure 6A and B). In addition, the fact that the loss of the Piwi protein Twi1 blocks the accumulation of H3K27me3 in the parental MAC (Liu et al., 2007; Xu et al., 2021) without affecting lncRNA accumulation (Woo et al., 2016) indicates that heterochromatin formation is not a prerequisite for lncRNA transcription in Tetrahymena, in contrast to heterochromatin-dependent lncRNA transcription of piRNA clusters in Drosophila (Andersen et al., 2017). Altogether, Tetrahymena likely uses a unique evolutionary solution for transcribing lncRNAs to promote TDSD and small RNA-directed heterochromatin formation. We, therefore, believe that further investigations of the role of Ema2-dependent SUMOylation in lncRNA transcription will clarify an undescribed layer of regulatory mechanisms in small RNA-directed chromatin regulation.
Materials and methods
Strains and cell culture condition
Request a detailed protocolThe Tetrahymena thermophila wild-type strains B2086, CU427, and CU428 and the MIC defective ‘star’ strains B*VI and B*VII were obtained from the Tetrahymena stock center. The production of the EMA2 KO strains was described previously (Shehzada and Mochizuki, 2022). The other transgenic Tetrahymena strains are described below. The Tetrahymena strains reported in this manuscripts are available from the corresponding author upon reasonable request. Cells were grown in SPP medium (Gorovsky et al., 1975) containing 2% proteose peptone at 30 °C. For conjugation, growing cells (~5–7 × 105 /mL) of two different mating types were washed, pre-starved (~12–24 hr), and mixed in 10 mM Tris (pH 7.5) at 30 °C.
Oligo DNA, synthetic DNA, and antibodies
Request a detailed protocolThe primers and synthetic DNAs used in this study are listed in Supplementary file 1. The following antibodies were obtained from manufacturers: anti-H3K27me3 (Millipore, 07–449, RRID:AB_310624); anti-HA HA11 (Covance, clone 6B12, RRID:AB_291231); anti-H3K9me3 (Active Motif, 39162); anti-GST (BD Biosciences, 554805, RRID:AB_395536); anti-His (Proteintech, 66005–1-Ig, RRID:AB_11232599); anti-long dsRNA J2 (Jena Bioscience RNT-SCI-10010200, RRID:AB_2651015); anti-α-tubulin 12G10 (Developmental Studies Hybridoma Bank of the University of Iowa). The anti-Rpb3 antibody was described previously (Kataoka and Mochizuki, 2017). The anti-Smt3 antibody was a gift from Dr. James Forney (Purdue University) and described previously (Nasir et al., 2015).
Establishment of EMA2-HA strains
Request a detailed protocolThe C-terminal part of the EMA2 coding sequence was amplified from the genomic DNA of B2086 using EMA2-C-FW and EMA2-HA-C-RV. Similarly, the 3' flanking sequence of EMA2 was amplified using EMA2-F-FW and EMA2-F-RV. The two PCR products were combined by overlapping PCR using EMA2-C-FW and EMA2-F-RV and cloned into XbaI and XhoI sites of pBlueScriptSK (+) to obtain pEMA2. Next, the HA-neo3 construct was excised out from pHA-neo3 (Kataoka et al., 2010) using SalI and BamHI and cloned into the same restriction sites of pEMA2 to obtain pEMA2-HA-neo3. Finally, pEMA2-HA-neo3 were digested with XbaI and XhoI and introduced into the EMA2 MAC locus of B2086 and CU428 by homologous recombination using biolistic transformation (Bruns and Cassidy-Hanley, 2000) followed by the selection with 100 μg/mL paromomycin and 1 μg/mL CdCl2. Then, the endogenous EMA2 MAC copies were replaced by the EMA2-HA-neo3 locus through phenotypic assortment (Hamilton and Orias, 2000).
Immunofluorescent staining and western blotting
Request a detailed protocolImmunofluorescent staining of long dsRNA was performed as described (Woo et al., 2016). The other immunofluorescent staining experiments were performed as described (Loidl and Scherthan, 2004). For immunofluorescent staining, all the primary and respective secondary antibodies (Alexa-488-coupled anti-rabbit or mouse IgG, Invitrogen) were diluted 1:1000. For western blotting, all the primary and respective secondary antibodies (HRP-coupled anti-rabbit or mouse IgG, Jackson ImmunoResearch Lab) were diluted 1:10,000.
DNA elimination assay by FISH
Request a detailed protocolDNA elimination assay by FISH was performed as described (Loidl and Scherthan, 2004). The plasmid DNAs pMBR 4C1, pMBR 2, and Tlr IntB, which contain different parts of the Tlr1element were used to produce Cy5-labeled FISH probes by nick translation (Wuitschick et al., 2002). Images of stained cells were analyzed by Fiji software. Whole areas of the new MAC (A) and MIC (I) and a part of the cytoplasm (C) were manually encircled in the images; and the intensities of DAPI and Cy5 (FISH) signals in each area were measured. Then, the IES retention index was calculated as follows: 0.35 × {[(mean Cy3 intensity A - mean Cy3 intensity C)]/[(mean Cy3 intensity I - mean Cy3 intensity C)]}/{[(mean DAPI intensity A - mean DAPI intensity C) ]/[(mean DAPI intensity I - mean DAPI intensity C)]}. Based on our previous observation, DNA elimination is completely blocked in the absence of TWI1 KO, therefore, the correction factor of 0.35 was obtained experimentally to make the average IES retention index of TWI1 KO cells to 1.
Viability test for progeny
Request a detailed protocolTo test the progeny viability, mating pairs were isolated into drops of 1 x SPP at 6 hpm and incubated for ~3 days. Then, the completion of the conjugation of grown-up cells from wild-type and EMA2 KO pairs were examined for their 6-methylpurine (6-mp) resistance and paromomycin sensitivity, respectively as described previously (Mochizuki et al., 2002).
Small RNA analyses
Request a detailed protocolFor northern hybridization analysis of scnRNAs, 10 µg of total RNAs isolated from conjugating wild-type and EMA2 KO cells using TRIzol reagent were separated in a 15% acrylamide-urea gel and analyzed by northern blot using the radio-labeled Mi-9 probe as previously described (Aronica et al., 2008). The hybridized probe was detected by Typhoon IP Phosphorimager (GE healthcare). High-throughput sequencing and analyses of small RNAs were performed as described (Mutazono et al., 2019). The 2020 version of the MAC genome assembly (Sheng et al., 2020) was fragmented into 10 kb pieces. Each 10 kb fragment that contains longer than 3 kb mappable sequence was used as an MDS tile (total 10,235 tiles). Each Type-A IES (Noto et al., 2015) was used as an IES tile (total of 4691 tiles). Normalized numbers (RPKM [reads per kilobase of unique sequences per million]) of 26- to 32-nt small RNAs that uniquely matched one of the MDS and IES tiles were counted.
GST pull down assay
Request a detailed protocolE. coli codon-optimized synthetic EMA2 and UBC9 genes (EMA2-Ec and UBC9-Ec) were cloned into the BamHI and XhoI sites of pGEX-4T-1 and NdeI and BamHI sites of pET28a (+), respectively to obtain pGEX-4T-1-EMA2-Ec and pET28-UBC9-Ec. pET28-UBC9-Ec was introduced into BL21 E. coli strain, His-Ubc9 expression was induced ~16 hr at 25 °C with 0.05 mM isopropylthio-β-galactoside (IPTG), the cells were lysed in HEPES buffer (50 mM HEPES pH 7.5, 800 mM KCL, 10% glycerol, 0.2 mM PMSF, 1 x cOmplete ETDA-free [Roche]) including 10 mM imidazole, incubated with Ni-NTA beads for 2 hr at 4 °C, washed four times with HEPES buffer including 20 mM imidazole, eluted with HEPES buffer including 250 mM imidazole, dialyzed with 50 mM Tris pH 7.5, 150 mM NaCl, 2 mM MgCl2, 1 mM DTT, 50% glycerol and finally passed through Microcon 10 kDa filter (Merck). For GST pull-down assay, BL21 E. coli cells were transformed with pGEX-4T-1 or pGEX-4T-1-EMA2-Ec, and expression of GST-alone or GST-Ema2 was induced ~16 hr at 25 °C, lysed with 50 mM Tris–HCl, pH 7.5, 100 mM NaCl, 1 mM DTT, 0.2 mM PMSF, 1x cOmplete EDTA-free and incubated with Glutathione Sepharose 4B for 90 min at 4 °C. Then, the beads were incubated with 10 µg of His-Ubc9 in GST-pull down buffer (20 mM Tris–HCl, pH 7.5, 100 mM NaCl, 1% NP-40, 5 mM MgCl2, 1 mM DTT, 0.2 mM PMSF, 1 x cOmplete EDTA-free) for 90 min at 4 °C, washed four times with GST-pull down buffer and eluted with 2 x SDS sample buffer. The eluted proteins were analyzed by western blotting.
Functionality test of HA-SMT3 and His-SMT3 using SMT3 germline KO strains
Request a detailed protocolThe SMT3 homozygous heterokaryon KO strain GC6 (Nasir et al., 2015) was obtained from Dr. James Forney (Purdue University). This strain was crossed with B2086, and sexual progeny were selected for their paromomycin resistance. After sexual maturation, clones showing paromomycin sensitivity (i.e. the SMT3 KO locus was removed from the MAC by phenotypic assortment) were established as SMT3 heterozygous heterokaryon KO strains. Then they were crossed with B*VII to obtain their round I exconjugants. Finally, their genotypes were examined by PCR to select two SMT3 homozygous heterokaryon KO strains, B1-1 and B5-17. They were then mated and pBCMB1-His-SMT3 (or pBCMB1-His for negative control) or pBCS3B1-HA-SMT3 (or pBCS3B1-HA for negative control) were introduced into the new MAC at 8 hpm by particle gun. The cells were incubated in 10 mM Tris pH 7.5 at 30 °C until 24 hpm and fed by adding the equal volume of 2x SPP. After 3 hr incubation at 30 °C, the neo expression (and His-Smt3 expression for pBCMB1-His-SMT3) was induced by adding 1 µg/mL CdCl2 and cells were incubated an additional 1 hr. Then 100 µg/mL paromomycin was added and aliquoted to 96-well plates (150 µL/well, seven plates). The plates were incubated at 30 °C for 4–5 days and paromomysin-resistant cells were further examined by genomic PCR to confirm genetic rescue by the introduced constructs (but not by the endogenous SMT3 copy from the parental MAC).
Establishment of HA-Smt3 and His-Smt3 expressing strains
Request a detailed protocolSTM3 cDNA was amplified by RT-PCR from the total RNA of vegetatively growing CU428 using SMT3_HA_FW and SMT3_HisExt_RV and cloned into the BamHI and SpeI sites of pBNMB1-HA using NEBuilder HiFi DNA assembly kit to obtain pBNMB1-HA-SMT3. Next, neo5 was excised with SalI and XmaI and replaced with pur6, which was amplified from pBP6MB1 with BTU1-LCFW and MTT1_5UTR-SeqRV, using NEBuilder HiFi DNA assembly kit to obtain pBP6MB1-HA-SMT3. For His-SMT3 expression, the HA-tag encoding sequence of pBP6MB1-HA-SMT3 was replaced by a 6x His tag-encoding sequence. Then, pBP6MB1-HA-SMT3 and pBP6MB1-His-SMT3 were digested by XhoI and introduced into the MAC BTU1 locus of the EMA2 KO cells. Transformation, selection, and phenotypic assortment were performed as described for the SPT6-HA strains.
Affinity purification of SUMOylated proteins
Request a detailed protocolFor purification of His-Smt3-conjugated proteins, EMA2 KO cells containing the BP6MB1-His-SMT3 construct were mated with either a wild-type or another EMA2 KO strain and His-SMT3 was induced by adding 0.05 μg/mL CdCl2. The purification was performed as descried (Hendriks and Vertegaal, 2016; Liebelt et al., 2019). Briefly, the cells were lysed in Guanidine Lysis Buffer (6 M Guanidine-HCl, 93.2 mM Na2HPO4, 6.8 mM NaH2PO4, 10 mM Tris-HCl pH 8.0) by sonication (2 × 25 times, 20% power, 60% pulse) using a probe sonicator (Omni-Ruptor 250 Ultrasonic Homogenizer). After centrifugation, the cleared lysate was supplemented with 50 mM imidazole and 5 mM β-mercaptoethanol and incubated with Ni-NTA agarose beads. The beads were washed, and the bound proteins were eluted with Elution buffer (6 M Urea, 58 mM Na2HPO4, 42 mM NaH2PO4, 10 mM Tris-HCl pH 8.0, 500 mM imidazole). The eluted proteins were precipitated in 10% TCA, resuspended in 1 x SDS buffer, and analyzed by western blot and mass-spectrometry. Wild-type cells that do not express His-Smt3 was also analyzed in parallel to identify proteins that intrinsically bind to the Ni-NTA agarose beads. For purification of HA-Smt3-conjugated proteins, total proteins were precipitated by incubating cells in 6% trichloroacetic (TCA) on ice for 10 min and then spinned down at 13,000 rpm for 5 min at 4 °C. The precipitate was dissolved in 1 x SDS sample buffer, neutralized with 2 M Tris, and then incubated for 5 min at 95 °C. The dissolved proteins (600 µl) were then suspended in 7.4 mL of 20 mM Tris pH 7.5, 100 mM NaCl, 2 mM MgCl2, 2 mM CaCl2 and incubated with Ezview Red anti-HA agarose beads (Sigma) for overnight at 4 °C. Then the beads were washed five times with 20 mM Tris pH 7.5, 100 mM NaCl, 2 mM MgCl2, 2 mM CaCl2, 0.1% tween 20 and the bound proteins were eluted by incubating with 1 x SDS buffer for 5 min at 95 °C. The eluted proteins were analyzed by western blotting.
Establishment of Spt6-HA expressing strains
Request a detailed protocolTo obtain the SPT6-HA construct, a part of the coding sequence of SPT6 was amplified from the genomic DNA of B2086 using SPT6-FW1-SacII and SPT6-RV2-overlap. Similarly, the 3' flanking sequence was amplified using SPT6-3F-FW2-overlap and SPT6-3F-RV1-KpnI. The two PCR products were combined by overlapping PCR using SPT6-FW1-SacII and SPT6-3F-RV1-KpnI and cloned into the SacII and KpnI sites of pBlueScriptSK (+) to obtain pSPT6. Next, the NheI and XhoI fragment of pHA-neo3 was cloned into the same restriction sites of pSPT6 to obtain pSPT6-HA-neo3. Then, neo3 was excised from pSPT6-HA-neo3 with XmaI and SalI and replaced by pur6 from pBP6MB1-TR-TUBE to obtain pSPT6-HA-pur6. Finally, the pSPT6-HA-pur6 plasmid were digested by SacII and KpnI and introduced into the MAC SPT6 locus of an EMA2 KO and a wild-type strain by biolistic transformation. The transformed cells were selected in 300 μg/mL puromycin (Invitrogen) and 1 μg/mL CdCl2 and cells were phenotypically assorted as previously described (Hamilton and Orias, 2000), until cells grew in 1.2 mg/mL of puromycin.
RT-PCR analyses of lncRNA transcripts
Request a detailed protocolTotal RNA from mating wild-type and EMA2 KO cells were isolated at 6 hpm using TRIzol reagent (Invitrogen). The RNA was treated with Turbo DNase (Ambion) and cDNA was synthesized from 1 µg of RNA using SuperScript II (Invitrogen) with random 6-mer as primers. The lncRNA transcripts from the MAC loci containing the M-, L8-, and R2 IES boundaries (see also Figure 6A) were amplified by nested PCR (95 °C for 20 s/50 °C for 30 s/68 °C for 1 min; 30 cycles each) using the primers M5'–3+M3'–3 (M-1st), M5'–4+M5'–4 (M-2nd), L8_5'–1+L8_3'–1 (L8-1st), L8_5'–2+L8_3'–2 (L8-2nd), R2_5'–1+R2_3'–1 (R2-1st), and R2_5'–2+R2_3'–2 (R2-2nd) and analyzed in 1% agarose gel. As a positive control, the constitutively expressed RPL21 mRNA was detected with RPL21-FW and RPL21-RV.
Cell fractionation
Request a detailed protocolCell fractionation was performed as described (Ali et al., 2018). 1.4 x 106 cells were incubated with 1 mL of Lysis buffer (50 mM Tris pH 8.0, 100 mM NaCl, 5 mM MgCl2, 1 mM EDTA, 0.5% Triton X-100, 1 x Complete Protease Inhibitor) on ice for 30 min and spun down at 15,000 g for 10 min at 4 °C. The proteins in 500 µL of supernatant were precipitated by adding 10% (final) TCA, washed three times with ice-cold acetone, and dissolved in 70 µL 1x SDS sample buffer, while the precipitate was dissolved in 2 mL of fresh Lysis buffer, sonicated (five times of 20% power-50% pulse-5 pulses with >20 s interval) using a probe sonicator (Omni-Ruptor 250 Ultrasonic Homogenizer), and then centrifuged at 20,000 g for 10 min at 4 °C. Proteins in 1 mL of the supernatant (solubilized chromatin) were precipitated by adding 10% (final) TCA, washed three times with ice-cold acetone and dissolved in 70 µL 1x SDS sample buffer. For cell Fractionation with RNase treatment, 1.4 × 106 cells were incubated with 1 mL of Lysis buffer containing 20 µg/mL RNase A on ice for 30 min and spun own at 15,000 g for 10 min at 4 °C. The proteins in 1 mL of the supernatant were precipitated by adding 10% (final) TCA, washed three times with ice-cold acetone, and dissolved in 70 µL 1x SDS sample buffer, while the precipitate was directly dissolved in 140 µL of 1 x SDS sample buffer.
Establishment of SPT6 germline KO strains
Request a detailed protocolTo create the knockout construct for SPT6, the 5’ and 3’ flanking regions of the SPT6 gene were amplified by PCR using the primer sets SPT6-KO-5F-FW/SPT6-KO-5F-RV and SPT6-KO-3F-FW/SPT6-KO-3F-RV, respectively and concatenated by overlap PCR using SPT6-KO-5F-FW and SPT6-KO-3F-RV. The PCR product was cloned into the XhoI and BamHI digested pBlueScriptSK(+). Then the obtained plasmid was digested with BamHI and the pm-resistant cassette neo4, which was isolated from the pNeo4 plasmid (Aronica et al., 2008) by SmaI digestion (Bruns and Cassidy-Hanley, 2000), was inserted. All assembly reactions were done using NEBuilder HiFi DNA assembly mix (NEB). The resulting plasmid, pSPT6-KO-neo4, was digested with SacII and XhoI to release the targeting construct from the vector backbone and introduced into the MIC of mating UMPS214 and UMPS811 cells (Vogt and Mochizuki, 2013) by biolistic transformation. Heterozygous progenies were selected with 100 µg/mL pm and 500 µg/mL 5-Fluoroorotic Acid (5-FOA, Fermentas) in the presence of 1 μg/ml CdCl2. The resulting heterozygous strains were cultured for 10 passages without pm and 5-FOA, and pm-sensitive (pm-s) heterozygous strains were isolated. The pm-s heterozygous cells were then mated with B*VI or B*VII to obtain their round I exconjugants. Finally, their genotypes were examined by PCR to select SPT6 homozygous heterokaryon KO strains.
Establishment of wild-type and mutant SPT6 strains
Request a detailed protocolSPT6 cDNA was amplified from total RNA of mating B2086 and CU428 at 6 hpm with SPT6-FW-Turbo and SPT6-RV-Turbo, cloned into SpeI and BamHI sites of pBCMB1-HA-Turbo to obtain pBCMB1-HA-Turbo-SPT6. Then, the MTT1 promoter and the HA-Turbo-tag were removed using AgeI and XmaI and replaced by the SPT6 promoter and the N-terminal HA tag, using the overlapping product of two genomic PCR fragments amplified with SPT6-5F-FW, SPT6-5F-RV and with SPT6-N-FW-HA and SPT6-N-RV to obtain pBCSB-HA-SPT6. To produce pBCSB-HA-SPT6-DOL-KR, the SPT6 cDNA lacking the first 248 amino acids were amplified using SPT6_M249_FW and SPT6_M249_RV, replaced with the NdeI-BamHI fragment of pBCSB-HA-SPT6 and then the synthetic DNA SPT6-DOL-KtoR was inserted into the BamHI site. To make pBCSB-HA-SPT6-N-KR, pBCSB-HA-SPT6-M-KR and pBCSB-HA-SPT6-C-KR, the conserved region of SPT6 was excised from pBCSB-HA-SPT6 using NdeI and AgeI and replaced respectively with Synthetic DNA Set 1 (SPT6-gBlock-118KR-1, SPT6-gBlock-42KR-2, and SPT6-gBlock-42KR-3), Synthetic DNA Set 2 (SPT6-gBlock-42KR-1, SPT6-gBlock-118KR-2, and SPT6-gBlock-42KR-3), or Synthetic DNA Set 3 (SPT6-gBlock-42KR-1, SPT6-gBlock-42KR-2, and SPT6-gBlock-118KR-3). All assembly reactions were done using NEBuilder HiFi DNA assembly mix (NEB). The constructs were digested with XhoI and introduced into the BTU1 locus of the new MAC of conjugating cells of two SPT6 KO homozygous heterokaryon strains at 10 hpm by electroporation. Progeny rescued by the introduced constructs were selected in 100 μg/mL paromomycin and 1 µg/mL CdCl2, and cells were phenotypically assorted until cells grew in ~100–120 µg/mL cycloheximide.
Data availability
The small RNA sequencing data have been deposited in the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo) with the accession no. GSE243435. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the MassIVE partner repository (https://massive.ucsd.edu/) with the database identifier MSV000092977.
-
NCBI Gene Expression OmnibusID GSE243435. A SUMO E3 ligase promotes long non-coding RNA transcription to regulate small RNA-directed DNA elimination.
-
MassIVE partner repositoryID MSV000092977. Comparison of SUMOylated proteins between wild-type and EMA2 KO Tetrahymena cells.
References
-
Panoramix SUMOylation on chromatin connects the piRNA pathway to the cellular heterochromatin machineryNature Structural & Molecular Biology 29:130–142.https://doi.org/10.1038/s41594-022-00721-x
-
Biolistic transformation of macro- and micronucleiMethods in Cell Biology 62:501–512.https://doi.org/10.1016/s0091-679x(08)61553-8
-
Epigenetics of ciliatesCold Spring Harbor Perspectives in Biology 5:a017764.https://doi.org/10.1101/cshperspect.a017764
-
Isolation of micro- and macronuclei of Tetrahymena pyriformisMethods in Cell Biology 9:311–327.https://doi.org/10.1016/s0091-679x(08)60080-1
-
A meiosis-specific Spt5 homolog involved in non-coding transcriptionNucleic Acids Research 45:4722–4732.https://doi.org/10.1093/nar/gkw1318
-
Genetic crosses: setting up crosses, testing progeny, and isolating phenotypic assortantsMethods in Cell Biology 62:219–228.https://doi.org/10.1016/s0091-679x(08)61532-0
-
MicroRNA turnover: a tale of tailing, trimming, and targetsTrends in Biochemical Sciences 48:26–39.https://doi.org/10.1016/j.tibs.2022.06.005
-
Label-free identification and quantification of SUMO target proteinsMethods in Molecular Biology 1475:171–193.https://doi.org/10.1007/978-1-4939-6358-4_13
-
RNA-mediated epigenetic regulation of gene expressionNature Reviews. Genetics 16:71–84.https://doi.org/10.1038/nrg3863
-
Modules for C-terminal epitope tagging of Tetrahymena genesJournal of Microbiological Methods 82:342–346.https://doi.org/10.1016/j.mimet.2010.07.009
-
Heterochromatin aggregation during DNA elimination in Tetrahymena is facilitated by a prion-like proteinJournal of Cell Science 130:480–489.https://doi.org/10.1242/jcs.195503
-
Silencing-associated and meiosis-specific small RNA pathways in Paramecium tetraureliaNucleic Acids Research 37:903–915.https://doi.org/10.1093/nar/gkn1018
-
Organization and pairing of meiotic chromosomes in the ciliate Tetrahymena thermophilaJournal of Cell Science 117:5791–5801.https://doi.org/10.1242/jcs.01504
-
Tetrahymena meiosis: Simple yet ingeniousPLOS Genetics 17:e1009627.https://doi.org/10.1371/journal.pgen.1009627
-
The transient Spt4-Spt5 complex as an upstream regulator of non-coding RNAs during developmentNucleic Acids Research 50:2603–2620.https://doi.org/10.1093/nar/gkac106
-
SUMO playing tag with ubiquitinTrends in Biochemical Sciences 37:23–31.https://doi.org/10.1016/j.tibs.2011.09.002
-
The small ubiquitin-like modifier-1 (SUMO-1) consensus sequence mediates Ubc9 binding and is essential for SUMO-1 modificationThe Journal of Biological Chemistry 276:21664–21669.https://doi.org/10.1074/jbc.M100006200
-
Keeping the soma free of transposons: programmed DNA elimination in ciliatesThe Journal of Biological Chemistry 286:37045–37052.https://doi.org/10.1074/jbc.R111.276964
-
Monoclonal antibodies to double-stranded RNA as probes of RNA structure in crude nucleic acid extractsNucleic Acids Research 19:2993–3000.https://doi.org/10.1093/nar/19.11.2993
-
PIWI-Directed DNA elimination for tetrahymena geneticsMethods in Molecular Biology 2509:53–68.https://doi.org/10.1007/978-1-0716-2380-0_3
-
The molecular interface between the SUMO and ubiquitin systemsAdvances in Experimental Medicine and Biology 963:99–110.https://doi.org/10.1007/978-3-319-50044-7_6
-
Molecular mechanisms of RNA interferenceAnnual Review of Biophysics 42:217–239.https://doi.org/10.1146/annurev-biophys-083012-130404
-
A novel family of mobile genetic elements is limited to the germline genome in Tetrahymena thermophilaNucleic Acids Research 30:2524–2537.https://doi.org/10.1093/nar/30.11.2524
Article and author information
Author details
Funding
Agence Nationale de la Recherche (ANR-10-INBS-04)
- Kazufumi Mochizuki
Agence Nationale de la Recherche (ANR-10-LABX-12-01)
- Kazufumi Mochizuki
Agence Nationale de la Recherche (ANR-16-ACHN-0017)
- Kazufumi Mochizuki
Fondation pour la Recherche Médicale (FDT20210601285)
- Salman Shehzada
Fondation pour la Recherche Médicale (EQU202203014651)
- Kazufumi Mochizuki
Fondation ARC (ARCPJA2021060003830)
- Kazufumi Mochizuki
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We acknowledge the MRI facility, member of the national infrastructure France-BioImaging supported by the French National Research Agency (ANR-10-INBS-04), the Proteomics Platform (PPM) of BioCampus Montpellier, and the NGS unit of Vienna BioCenter Core Facilities. This work was supported by an Advanced Grant from the 'Investissements d’avenir' program Labex EpiGenMed (ANR-10-LABX-12–01) and an 'Accueil de Chercheurs de Haut Niveau' grant (ANR-16-ACHN-0017) from the French National Research Agency, Equipes a FRM 2022 grant from the Fondation Recherche pour Médicale (FRM, EQU202203014651), an ARC 2021 PJA3 grant from the ARC Foundation (ARCPJA2021060003830) to KM, and a 'Fin de these' program fellowship from Fondation pour la Recherche Médicale (FRM, FDT20210601285) to SS.
Copyright
© 2024, Shehzada 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
-
- 1,016
- views
-
- 281
- 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
- Developmental Biology
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.