- Jun 2023
-
www.biorxiv.org www.biorxiv.org
-
Author Response:
The following is the authors' response to the original reviews.
Reply to Public Reviews:
Reply to Reviewer #1:
This is a carefully performed and well-documented study to indicate that the FUS protein interacts with the GGGGCC repeat sequence in Drosophila fly models, and the mechanism appears to include modulating the repeat structure and mitigating RAN translation. They suggest FUS, as well as a number of other G-quadruplex binding RNA proteins, are RNA chaperones, meaning they can alter the structure of the expanded repeat sequence to modulate its biological activities.
Response: We would like to thank the reviewer for her/his time for evaluating our manuscript. We are very happy to see the reviewer for highly appreciating our manuscript.
1. Overall this is a nicely done study with nice quantitation. It remains somewhat unclear from the data and discussions in exactly what way the authors mean that FUS is an RNA chaperone: is FUS changing the structure of the repeat or does FUS binding prevent it from folding into alternative in vivo structure?
Response: We appreciate the reviewer’s constructive comments. Indeed, we showed that FUS changes the higher-order structures of GGGGCC [G4C2] repeat RNA in vitro, and that FUS suppresses G4C2 RNA foci formation in vivo. According to the established definition of RNA chaperone, RNA chaperones are proteins changing the structures of misfolded RNAs without ATP use, resulting in the maintenance of proper RNAs folding (Rajkowitsich et al., 2007). Thus, we consider that FUS is classified into RNA chaperone. To clarify these interpretations, we revised the manuscript as follows.
(1) On page 10, line 215-219, the sentence “These results were in good agreement with our previous study on SCA31 showing the suppressive effects of FUS and other RBPs on RNA foci formation of UGGAA repeat RNA as RNA chaperones …” was changed to “These results were in good agreement with … RNA foci formation of UGGAA repeat RNA through altering RNA structures and preventing aggregation of misfolded repeat RNA as RNA chaperones …”.
(2) On page 17, line 363-366, the sentence “FUS directly binds to G4C2 repeat RNA and modulates its G-quadruplex structure, as evident by CD and NMR analyses (Figure 5), suggesting its functional role as an RNA chaperone.” was changed to “FUS directly binds to G4C2 repeat RNA and modulates its G-quadruplex structure as evident by CD and NMR analyses (Figure 5, Figure 5—figure supplement 2), and suppresses RNA foci formation in vivo (Figures 3A and 3B), suggesting its functional role as an RNA chaperone.”
Reply to Reviewer #2:
Fuijino et al. provide interesting data describing the RNA-binding protein, FUS, for its ability to bind the RNA produced from the hexanucleotide repeat expansion of GGGGCC (G4C2). This binding correlates with reductions in the production of toxic dipeptides and reductions in toxic phenotypes seen in (G4C2)30+ expressing Drosophila. Both FUS and G4C2 repeats of >25 are associated with ALS/FTD spectrum disorders. Thus, these data are important for increasing our understanding of potential interactions between multiple disease genes. However, further validation of some aspects of the provided data is needed, especially the expression data.
Response: We would like to thank the reviewer for her/his time for evaluating our manuscript and also for her/his important comments that helped to strengthen our manuscript.
Some points to consider when reading the work:
1. The broadly expressed GMR-GAL4 driver leads to variable tissue loss in different genotypes, potentially confounding downstream analyses dependent on viable tissue/mRNA levels.
Response: We thank the reviewer for this constructive comment. In the RT-qPCR experiments (Figures 1E, 3C, 4G, 6D and Figure 1—figure supplement 1C), the amounts of G4C2 repeat transcripts were normalized to those of gal4 transcripts expressed in the same tissue, to avoid potential confounding derived from the difference in tissue viability between genotypes, as the reviewer pointed out. To clarify this process, we have made the following change to the revised manuscript.
(1) On page 30, line 548-550, the sentence “The amounts of G4C2 repeat transcripts were normalized to those of gal4 transcripts in the same sample” was changed to “The amounts of G4C2 repeat transcripts were normalized to those of gal4 transcripts expressed in the same tissue to avoid potential confounding derived from the difference in tissue viability between genotypes”.
2. The relationship between FUS and foci formation is unclear and should be interpreted carefully.
Response: We appreciate the reviewer’s important comment. We apologize for the lack of clarity. We showed the relationship between FUS and RNA foci formation in our C9-ALS/FTD fly, that is, FUS suppresses RNA foci formation (Figures 3A and 3B), and knockdown of endogenous caz, a Drosophila homologue of FUS, enhanced it conversely (Figures 4E and 4F). We consider that FUS suppresses RNA foci formation through altering RNA structures and preventing aggregation of misfolded G4C2 repeat RNA as an RNA chaperone. To clarify these interpretations, we revised the manuscript as follows.
(1) On page 10, line 215-219, the sentence “These results were in good agreement with our previous study on SCA31 showing the suppressive effects of FUS and other RBPs on RNA foci formation of UGGAA repeat RNA as RNA chaperones …” was changed to “These results were in good agreement with … RNA foci formation of UGGAA repeat RNA through altering RNA structures and preventing aggregation of misfolded repeat RNA as RNA chaperones …”.
(2) On page 17, line 363-366, the sentence “FUS directly binds to G4C2 repeat RNA and modulates its G-quadruplex structure, as evident by CD and NMR analyses (Figure 5), suggesting its functional role as an RNA chaperone.” was changed to “FUS directly binds to G4C2 repeat RNA and modulates its G-quadruplex structure as evident by CD and NMR analyses (Figure 5, Figure 5—figure supplement 2), and suppresses RNA foci formation in vivo (Figures 3A and 3B), suggesting its functional role as an RNA chaperone.”
Reply to Reviewer #3:
In this manuscript Fujino and colleagues used C9-ALS/FTD fly models to demonstrate that FUS modulates the structure of (G4C2) repeat RNA as an RNA chaperone, and regulates RAN translation, resulting in the suppression of neurodegeneration in C9-ALS/FTD. They also confirmed that FUS preferentially binds to and modulates the G-quadruplex structure of (G4C2) repeat RNA, followed by the suppression of RAN translation. The potential significance of these findings is high since C9ORF72 repeat expansion is the most common genetic cause of ALS/FTD, especially in Caucasian populations and the DPR proteins have been considered the major cause of the neurodegenerations.
Response: We would like to thank the reviewer for her/his time for evaluating our manuscript. We are grateful to the reviewer for the insightful comments, which were very helpful for us to improve the manuscript.
1. While the effect of RBP as an RNA chaperone on (G4C2) repeat expansion is supposed to be dose-dependent according to (G4C2)n RNA expression, the first experiment of the screening for RBPs in C9-ALS/FTD flies lacks this concept. It is uncertain if the RBPs of the groups "suppression (weak)" and "no effect" were less or no ability of RNA chaperone or if the expression of the RBP was not sufficient, and if the RBPs of the group "enhancement" exacerbated the toxicity derived from (G4C2)89 RNA or the expression of the RBP was excessive. The optimal dose of any RBPs that bind to (G4C2) repeats may be able to neutralize the toxicity without the reduction of (G4C2)n RNA.
Response: We appreciate the reviewer’s constructive comments. We employed the site-directed transgenesis for the establishment of RBP fly lines, to ensure the equivalent expression levels of the inserted transgenes. We also evaluated the toxic effects of overexpressed RBPs themselves by crossbreeding with control EGFP flies, showing in Figure 1A. To clarify them, we have made the following changes to the revised manuscript.
(1) On page 8, line 166-168, the sentence “The variation in the effects of these G4C2 repeat-binding RBPs on G4C2 repeat-induced toxicity may be due to their different binding affinities to G4C2 repeat RNA, and their different roles in RNA metabolism.” was changed to “The variation in the effects of these G4C2 repeat-binding RBPs on G4C2 repeat-induced toxicity may be due to their different binding affinities to G4C2 repeat RNA, and the different toxicity of overexpressed RBPs themselves.”.
(2) On page 29, line 519-522, the sentence “By employing site-specific transgenesis using the pUASTattB vector, each transgene was inserted into the same locus of the genome, and was expected to be expressed at the equivalent levels.” was added.
2. In relation to issue 1, the rescue effect of FUS on the fly expressing (G4C2)89 (FUS-4) in Figure 4-figure supplement 1 seems weaker than the other flies expressing both FUS and (G4C2)89 in Figure 1 and Figure 1-figure supplement 2. The expression level of both FUS protein and (G4C2)89 RNA in each line is important from the viewpoint of therapeutic strategy for C9-ALS/FTD.
Response: We appreciate the reviewer’s important comment. The FUS-4 transgene is expected to be expressed at the equivalent level to the FUS-3 transgene, since they are inserted into the same locus of the genome by the site-directed transgenesis. Thus, we suppose that the weaker suppressive effect of FUS-4 coexpression on G4C2 repeat-induced eye degeneration can be attributed to the C-terminal FLAG tag that is fused to FUS protein expressed in FUS-4 fly line. Since the caz fly expresses caz protein also fused to FLAG tag at the C-terminus, we used this FUS-4 fly line to directly compare the effect of caz on G4C2 repeat-induced toxicity to that of FUS.
3. While hallmarks of C9ORF72 are the presence of DPRs and the repeat-containing RNA foci, the loss of function of C9ORF72 is also considered to somehow contribute to neurodegeneration. It is unclear if FUS reduces not only the DPRs but also the protein expression of C9ORF72 itself.
Response: We thank the reviewer for this comment. We agree that not only DPRs, but also toxic repeat RNA and the loss-of-function of C9ORF72 jointly contribute to the pathomechanisms of C9-ALS/FTD. Since Drosophila has no homolog corresponding to the human C9orf72 gene, the effect of FUS on C9orf72 expression cannot be assessed. Our fly models are useful for evaluating gain-of-toxic pathomechanisms such as RNA foci formation and RAN translation, and the association between FUS and loss-of function of C9ORF72 is beyond the scope of this study.
4. In Figure 5E-F, it cannot be distinguished whether FUS binds to GGGGCC repeats or the 5' flanking region. The same experiment should be done by using FUS-RRMmut to elucidate whether FUS binding is the major mechanism for this translational control. Authors should show that FUS binding to long GGGGCC repeats is important for RAN translation.
Response: We would like to thank the reviewer for these insightful comments. Following the reviewer’s suggestion, we perform in vitro translation assay again using FUS-RRMmut, which loses the binding ability to G4C2 repeat RNA as evident by the filter binding assay (Figure 5A), instead of BSA. The results are shown in the figures of Western blot analysis below. The addition of FUS to the translation system suppressed the expression levels of GA-Myc efficiently, whereas that of FUS-RRMmut did not. FUS decreased the expression level of GA-Myc at as low as 10nM, and nearly eliminated RAN translation activity at 100nM. At 400nM, FUS-RRMmut weakly suppressed the GA-Myc expression levels probably because of the residual RNA-binding activity. These results suggest that FUS suppresses RAN translation in vitro through direct interactions with G4C2 repeat RNA.
Unfortunately, RAN translation from short G4C2 repeat RNA was not investigated in our translation system, although the previous study reported the low efficacy of RAN translation from short G4C2 repeat RNA (Green et al., 2017).
Author response image 1.
(A) Western blot analysis of the GA-Myc protein in the samples from in vitro translation.
(B) Quantification of the GA-Myc protein levels.
We have made the following changes to the revised manuscript.
(1) Figure 5F was replaced to new Figures 5F and 5G.
(2) On page 14-15, line 326-330, the sentence “Notably, the addition of FUS to this system decreased the expression level of GA-Myc in a dose-dependent manner, whereas the addition of the control bovine serum albumin (BSA) did not (Figure 5F).” was changed to “Notably, upon the addition to this translation system, FUS suppressed RAN translation efficiently, whereas FUS-RRMmut did not. FUS decreased the expression levels of GA-Myc at as low as 10nM, and nearly eliminated RAN translation activity at 100nM. At 400nM, FUS-RRMmut weakly suppressed the GA-Myc expression levels probably because of the residual RNA-binding activity (Figure 5F and 5G).”.
(3) On page 15, line 330-332, the sentence “Taken together, these results indicate that FUS suppresses RAN translation from G4C2 repeat RNA in vitro as an RNA chaperone.” was changed to “Taken together, these results indicate that FUS suppresses RAN translation in vitro through direct interactions with G4C2 repeat RNA as an RNA chaperone.”.
(4) On page 37, line 720-723, the sentence “For preparation of the FUS protein, the human FUS (WT) gene flanked at the 5¢ end with an Nde_I recognition site and at the 3¢ end with a _Xho_I recognition site was amplified by PCR from pUAST-_FUS.” was changed to “For preparation of the FUS proteins, the human FUS (WT) and FUS-RRMmut genes flanked at the 5¢ end with an Nde_I recognition site and at the 3¢ end with a _Xho_I recognition site was amplified by PCR from pUAST-_FUS and pUAST- FUS-RRMmut, respectively.”.
(5) On page 41, line 816-819, the sentence “FUS or BSA at each concentration (10, 100, and 1,000 nM) was added for translation in the lysate.” was changed to “FUS or FUS-RRMmut at each concentration (10, 100, 200, 400, and 1,000 nM) was preincubated with mRNA for 10 min to facilitate the interaction between FUS protein and G4C2 repeat RNA, and added for translation in the lysate.”.
5. It is not possible to conclude, as the authors have, that G-quadruplex-targeting RBPs are generally important for RAN translation (Figure 6), without showing whether RBPs that do not affect (G4C2)89 RNA levels lead to decreased DPR protein level or RNA foci.
Response: We appreciate the reviewer’s critical comment. Following the suggestion by the reviewer, we evaluate the effect of these G-quadruplex-targeting RBPs on RAN translation. We additionally performed immunohistochemistry of the eye imaginal discs of fly larvae expressing (G4C2)89 and these G-quadruplex-targeting RBPs. As shown in the figures of immunohistochemistry below, we found that coexpression of EWSR1, DDX3X, DDX5, and DDX17 significantly decreased the number of poly(GA) aggregates. The results suggest that these G-quadruplex-targeting RBPs regulate RAN translation as well as FUS.
Author response image 2.
(A) Immunohistochemistry of poly(GA) in the eye imaginal discs of fly larvae expressing (G4C2)89 and the indicated G-quadruplex-targeting RBPs.
(B) Quantification of the number of poly(GA) aggregates.
We have made the following changes to the revised manuscript.
(1) Figures 6E and 6F were added.
(2) On page 6-7, line 135-137, the sentence “In addition, other G-quadruplex-targeting RBPs also suppressed G4C2 repeat-induced toxicity in our C9-ALS/FTD flies.” was changed to “In addition, other G-quadruplex-targeting RBPs also suppressed RAN translation and G4C2 repeat-induced toxicity in our C9-ALS/FTD flies.”.
(3) On page 15, line 344-346, the sentence “As expected, these RBPs also decreased the number of poly(GA) aggregates in the eye imaginal discs (Figures 6E and 6F).” was added.
(4) On page 15, line 346-347, the sentence “Their effects on G4C2 repeat-induced toxicity and repeat RNA expression were consistent with those of FUS.” was changed to “Their effects on G4C2 repeat-induced toxicity, repeat RNA expression, and RAN translation were consistent with those of FUS.”
(5) On page 16, line 355-357, the sentence “Thus, some G-quadruplex-targeting RBPs regulate G4C2 repeat-induced toxicity by binding to and possibly by modulating the G-quadruplex structure of G4C2 repeat RNA.” was changed to “Thus, some G-quadruplex-targeting RBPs regulate RAN translation and G4C2 repeat-induced toxicity by binding to and possibly by modulating the G-quadruplex structure of G4C2 repeat RNA.”
(6) On page 19, line 417-421, the sentence “We further found that G-quadruplex-targeting RNA helicases, including DDX3X, DDX5, and DDX17, which are known to bind to G4C2 repeat RNA (Cooper-Knock et al., 2014; Haeusler et al., 2014; Mori et al., 2013a; Xu et al., 2013), also alleviate G4C2 repeat-induced toxicity without altering the expression levels of G4C2 repeat RNA in our Drosophila models.” was changed to “We further found that G-quadruplex-targeting RNA helicases, … ,also suppress RAN translation and G4C2 repeat-induced toxicity without altering the expression levels of G4C2 repeat RNA in our Drosophila models.”.
Reply to Recommendations For The Authors:
1) It is not clear from the start that the flies they generated with the repeat have an artificial vs human intronic sequence ahead of the repeat. It would be nice if they presented somewhere the entire sequence of the insert. The reason being that it seems they also tested flies with the human intronic sequence, and the effect may not be as strong (line 234). In any case, in the future, with a new understanding of RAN translation, it would be nice to compare different transgenes, and so as much transparency as possible would be helpful regarding sequences. Can they include these data?
Response: We thank the editors and reviewers for this comment. We apologize for the lack of clarity. We used artificially synthesized G4C2 repeat sequences when generating constructs for (G4C2)n transgenic flies, so these constructs do not contain human intronic sequence ahead of the G4C2 repeat in the C9orf72 gene, as explained in the Materials and Methods section. To clarify the difference between our C9-ALS/FTD fly models and LDS-(G4C2)44GR-GFP fly model (Goodman et al., 2019), we have made the following change to the revised manuscript.
(1) Schema of the LDS-(G4C2)44GR-GFP construct was presented in Figure 3—figure supplement 1.
Furthermore, to maintain transparency of the study, we have provided the entire sequence of the insert as the following source file.
(2) The artificial sequences inserted in the pUAST vector for generation of the (G4C2)n flies were presented in Figure 1—figure supplement 1—source data 1.
2) It is really nice how they quantitated everything and showed individual data points.
Response: We thank the editors and reviewers for appreciating our data analysis method. All individual data points and statistical analyses are summarized in source data files.
3) So when they call FUS an RNA chaperone, are they simply meaning it is changing the structure of the repeat, or could it just be interacting with the repeat to coat the repeat and prevent it from folding into whatever in vivo structures? Can they speculate on why some RNA chaperones lead to presumed decay of the repeat and others do not? Can they discuss these points in the discussion? Detailed mechanistic understanding of RNA chaperones that ultimately promote decay of the repeat might be of highly significant therapeutic benefit.
Response: We appreciate these critical comments. Indeed, we showed that FUS changes the higher-order structures of G4C2 repeat RNA in vitro, and that FUS suppresses G4C2 RNA foci formation. According to the established definition of RNA chaperone, RNA chaperones are proteins changing the structures of misfolded RNAs without ATP use, resulting in the maintenance of proper RNAs folding (Rajkowitsich et al., 2007). Thus, we consider that FUS is classified into RNA chaperone. To clarify these interpretations, we revised the manuscript as follows.
(1) On page 10, line 215-219, the sentence “These results were in good agreement with our previous study on SCA31 showing the suppressive effects of FUS and other RBPs on RNA foci formation of UGGAA repeat RNA as RNA chaperones …” was changed to “These results were in good agreement with … RNA foci formation of UGGAA repeat RNA through altering RNA structures and preventing aggregation of misfolded repeat RNA as RNA chaperones …”.
(2) On page 17, line 363-366, the sentence “FUS directly binds to G4C2 repeat RNA and modulates its G-quadruplex structure, as evident by CD and NMR analyses (Figure 5), suggesting its functional role as an RNA chaperone.” was changed to “FUS directly binds to G4C2 repeat RNA and modulates its G-quadruplex structure as evident by CD and NMR analyses (Figure 5, Figure 5—figure supplement 2), and suppresses RNA foci formation in vivo (Figures 3A and 3B), suggesting its functional role as an RNA chaperone.”
Besides these RNA chaperones, we observed the expression of IGF2BP1, hnRNPA2B1, DHX9, and DHX36 decreased G4C2 repeat RNA expression levels. In addition, we recently reported that hnRNPA3 reduces G4C2 repeat RNA expression levels, leading to the suppression of neurodegeneration in C9-ALS/FTD fly models (Taminato et al., 2023). We speculate these RBPs could be involved in RNA decay pathways as components of the P-body or interactors with the RNA deadenylation machinery (Tran et al., 2004; Katahira et al., 2008; Geissler et al., 2016; Hubstenberger et al., 2017), possibly contributing to the reduced expression levels of G4C2 repeat RNA. To clarify these interpretations, we revised the manuscript as follows.
(3) On page 18, line 392-398, the sentences “Similarly, we recently reported that hnRNPA3 reduces G4C2 repeat RNA expression levels, leading to the suppression of neurodegeneration in C9-ALS/FTD fly models (Taminato et al., 2023). Interestingly, these RBPs have been reported to be involved in RNA decay pathways as components of the P-body or interactors with the RNA deadenylation machinery (Tran et al., 2004; Katahira et al., 2008; Geissler et al., 2016; Hubstenberger et al., 2017), possibly contributing to the reduced expression levels of G4C2 repeat RNA.” was added.
4) What is the level of the G4C2 repeat when they knock down caz? Is it possible that knockdown impacts the expression level of the repeat? Can they show this (or did they and I miss it)?
Response: We thank the editors and reviewers for this comment. The expression levels of G4C2 repeat RNA in (G4C2)89 flies were not altered by the knockdown of caz, as shown in Figure 4G.
5) A puzzling point is that FUS is supposed to be nuclear, so where is FUS in the brain in their lines? They suggest it modulates RAN translation, and presumably, that is in the cytoplasm. Is FUS when overexpressed now in part in the cytoplasm? Is the repeat dragging it into the cytoplasm? Can they address this in the discussion? If FUS is never found in vivo in the cytoplasm, then it raises the point that the impact they find of FUS on RAN translation might not reflect an in vivo situation with normal levels of FUS.
Response: We appreciate these important comments. We agree with the editors and reviewers that FUS is mainly localized in the nucleus. However, FUS is known as a nucleocytoplasmic shuttling RBP that can transport RNA into the cytoplasm. Indeed, FUS is reported to facilitate transport of actin-stabilizing protein mRNAs to function in the cytoplasm (Fujii et al., 2005). Thus, we consider that FUS binds to G4C2 repeat RNA in the cytoplasm and suppresses RAN translation in this study.
6) When they are using 2 copies of the driver and repeat, are they also using 2 copies of FUS? These are quite high levels of transgenes.
Response: We thank the editors and reviewers for this comment. We used only 1 copy of FUS when using 2 copies of GMR-Gal4 driver. Full genotypes of the fly lines used in all experiments are described in Supplementary file 1.
7) In Figure5-S1, FUS colocalizing with (G4C2)RNA is not clear. High-magnification images are recommended.
Response: We appreciate this constructive comment on the figure. Following the suggestion, high-magnification images are added in Figure 5—figure supplement 1.
8) I also suggest that the last sentence of the Discussion be revised as follows: Thus, our findings contribute not only to the elucidation of C9-ALS/FTD, but also to the elucidation of the repeat-associated pathogenic mechanisms underlying a broader range of neurodegenerative and neuropsychiatric disorders than previously thought, and it will advance the development of potential therapies for these diseases.
Response: We appreciate this recommendation. We have made the following change based on the suggested sentence.
(1) On page 20-21, line 455-459, “Thus, our findings contribute not only towards the elucidation of repeat-associated pathogenic mechanisms underlying a wider range of neuropsychiatric diseases than previously thought, but also towards the development of potential therapies for these diseases.” was changed to “Thus, our findings contribute to the elucidation of the repeat-associated pathogenic mechanisms underlying not only C9-ALS/FTD, but also a broader range of neuromuscular and neuropsychiatric diseases than previously thought, and will advance the development of potential therapies for these diseases.”.
Authors’ comment on previous eLife assessment:
We thank the editors and reviewers for appreciating our study. We mainly evaluated the function of human FUS protein on RAN translation and G4C2 repeat-induced toxicity using Drosophila expressing human FUS in vivo, and the recombinant human FUS protein in vitro. To validate that FUS functions as an endogenous regulator of RAN translation, we additionally evaluated the function of Drosophila caz protein as well. We are afraid that the first sentence of the eLife assessment, that is, “This important study demonstrates that the Drosophila FUS protein, the human homolog of which is implicated in amyotrophic lateral sclerosis (ALS) and related conditions, …” is somewhat misleading. We would be happy if you modify this sentence like “This important study demonstrates that the human FUS protein, which is implicated in amyotrophic lateral sclerosis (ALS) and related conditions, …”.
-
- Feb 2023
-
www.biorxiv.org www.biorxiv.org
-
Note: This rebuttal was posted by the corresponding author to Review Commons. Content has not been altered except for formatting.
Learn more at Review Commons
Reply to the reviewers
We thank all the reviewers for having raised constructive criticism to fortify the main message and improve the clarity of the manuscript. We appreciate that all reviewers found that our work addresses an important topic and is of interest to a broad audience. We believe that we have thoroughly addressed the concerns of the reviewers, especially with regard to 1) performing another SMC3 chromatin immunoprecipitation and sequencing (ChIP-seq) replicate and control, 2) including a later time point for the transcriptional data, and 3) performing additional characterization of the growth phenotype of the SMC3 conditional knockdown.
Reviewer #1
(Evidence, reproducibility and clarity (Required)):*
Summary The present work by Rosa et al., provides convincing data about the presence and functional relevance of the cohesin complex in Plasmodium falciparum blood stages. In accordance with other organisms, the composition of the cohesin complex containing SMC1, SMC3 RAD21 and putatively STAG could be confirmed via pulldown and mass spectrometry. Basic characterization of endogenous tagged SMC3 demonstrated the expression and nuclear localization during IDC, as well as the relatively stable accumulation at centromeric regions, consistent with the known cohesin function in chromatid separation. Furthermore, dynamic and stage-dependent binding to intergenic regions observed in ChIPseq and major transcriptome aberrations upon knockdown of SMC3 (__Response: __As we regularly perform ChIP-seq experiments in the lab, we have generated multiple negative control datasets. In our opinion, the most stringent negative control for an HA-tagged protein is performing ChIP with an HA antibody in a WT strain. We have recently published an in-depth analysis of this (and other) negative ChIP-seq controls (Baumgarten & Bryant, 2022, https://doi.org/10.12688/openreseurope.14836.2). We show in this publication that non-specific ChIP-seq experiments (such as negative controls) result in an over-representation of HP1-heterochromatinized regions due to differences in sonication efficiency of heterochromatin and technical challenges with mapping regions with high levels of homology. In the anti-HA in WT ChIP negative control (performed at 12hpi), we do not see any enrichment at centromeric regions, but rather at heterochromatinized regions where clonally variant gene families are located. We performed peak calling analysis and found no significant overlap between the negative control ChIP-seq and the SMC3-3HA ChIP-seq data at 12hpi.
In addition, we have now performed a second biological replicate of the SMC3-3HA ChIP-seq with a different clone at all time points. We compared this data to that from the original clone and found significant overlap of the peaks called (see what is now Table 4 and Supp. Fig. 3A). We generated a stringent list of peaks that were shared between both clones at each time point and repeated all downstream analyses (see what are now Tables 5-8). We found that our conclusions were largely unchanged. Text describing these experiments and analyses have been added throughout the results section.
- Proposed mechanism of repressive effect of SMC3 early in IDC on genes, that get de-repressed in late stages: To claim this mode of function, it would be necessary to include a KD on late stage parasites. If there is an early repressive role of SMC3, upregulated genes should not be affected by late SMC3-KD. __Response: __To be clear, we are most interested in the transcriptional role of SMC3 during interphase, where results are not confounded by its potential role in mitosis. However, we did collect a 36hpi time point in the SMC3-3HA-glmS and WT strain, with and without glucosamine. We have added this last time point and the WT data from the other two time points to the manuscript (see Tables 11-13). Unfortunately, and for reasons unknown, the WT replicates treated with glucosamine showed a significantly advanced “transcriptional age” compared to the other replicates at 36hpi (see what is now Supp. Fig. 5B). Thus, we did not feel comfortable performing the RNA-seq analysis as we did with the other two time points (i.e. subtracting up- and down-regulated genes from the WT control from the SMC3-3HA-glmS data sets). We have added this information to the results section (Lines 256 and 261). As the WT parasites treated with glucosamine were approximately 8 hours in advance of the untreated WT parasites for the 36hpi time point, any up- and down-regulated genes might have been due to differences in the cell cycle rather than due to glucosamine treatment. The glmS system of inducible knockdown is widely used in P. falciparum; however, to our knowledge, no lab has investigated whether glucosamine treatment affects transcription in wildtype cells over the course of the IDC. Thus, for accurate phenotypic characterization of any protein with this system with regard to transcriptomics, we thought it was important to provide an RNA-seq dataset to define the cohort of genes affected by glucosamine treatment in WT parasites. We hope that our study will demonstrate the importance of using stringent controls when using inducible knockdown systems.
To address the question of whether genes that are upregulated upon depletion of SMC3 at early stages are affected at the 36hpi time point, we performed differential expression analysis of the SMC3-3HA-glmS parasites with and without glucosamine at 36hpi (we have added this data in Table 11). Again, significantly up- and down-regulated genes were not filtered using the WT dataset. With this analysis, we see only three genes from the list of invasion-related genes (Hu et al., 2010) that are up-regulated, but none of them have a significant q-value (Tab 5 of Table 18). Thus, depletion of SMC3 in late stage parasites does not lead to up-regulation of the same genes that are upregulated at 12 and 24hpi. We have added this information to the text (Line 273).
Furthermore, the hypothesized repressive effect of SMC3 does not explain the numerous genes downregulated in KD.
__Response: __As we state on line 350, we do not observe enrichment of SMC3 at downregulated genes, suggesting an indirect or secondary effect of SMC3 KD on these genes.
-
Due to the fact, that the KD was induced at the exact same timepoint and analysed 12h and 24h after induction it is possible that identified, differentially expressed genes at 24h are not directly regulated by SMC3, but rather due to a general deregulation of gene expression. Did the authors attempt to analyse gene expression upon induction at ring, trophozoite and schizont stage? Response: __As we state on line 230, in order to achieve SMC3 KD at the protein level, we had to treat the parasite with glucosamine for two cell cycles (approximately 96 hours). After two cell cycles of glucosamine treatment, the parasites were tightly synchronized and sampled 12 and 24 hours later. Thus, SMC3 KD takes place over the course of multiple days, but parasites are collected after stringent synchronization. Giemsa staining and bioinformatic analysis (line 250) of the RNA-seq data from parasites (with or without glucosamine) harvested at 12 and 24 hpi show that these parasites were synchronous and that there were no gross differences in genome-wide transcript levels. It is certainly possible that differentially expressed genes at 12 or 24hpi are not directly regulated by SMC3, and this is precisely why we perform ChIP-seq of SMC3: to provide evidence of direct involvement via binding, as stated on line 281. __
-
*Based on rapid parasite growth, the authors hypothesize a higher invasion rate due to upregulation of invasion genes. This hypothesis is not supported by quantitative invasion assays or quantification of invasion factors on the protein level. An alternative explanation could be a shorter cell cycle (__Response: __We have repeated the growth curve analysis with additional clones and no longer observe a growth phenotype in the SMC3 knockdown condition. We have added images of Giemsa-stained parasites from the knockdown time course we performed to what is now Supp. Fig. 5A. We see no obvious differences in cell morphology caused by glucosamine treatment in the WT or SMC3-3HA-glmS parasites.
-
Correlation of SMC3-occupancy/ATAC/expression profile of the exemplary genes rap2 and gap45 (Figure 4C,D,E): is this representative for all upregulated genes? __Response: __SMC3 occupancy shown at rap2 and gap45 is representative for all upregulated genes (see Fig. 4A and B). It is difficult to provide a general representation of the average expression profiles of all up-regulated genes over the course of the IDC, but Fig. 3E shows that the vast majority of up-regulated genes normally reach their peak expression in late stage parasites. With regard to ATAC-seq profiles, we have performed a metagene analysis of chromatin accessibility (data taken from (Toenhake et al., 2018)) at all up-regulated genes at time points that closely correspond to the time points used in our study: 15, 25, and 35, and 40 hpi (new Fig. 4C). This metagene analysis confirms what we observe at individual genes: increasing chromatin accessibility over the course of the IDC at these genes’ promoters. While metagene analyses offer important information, we always try to show the raw data (as in new Figs. 4D-F) from individual examples as proof of principle.
-
Given that SMC3 appears to be not essential for parasite growth, the authors could generate a null mutant for SMC3, which might allow for easier analysis of differences in gene regulation, cell cycle progression and/or invasion efficiency. __Response: __As we explain on line 327, very little cohesin is required for normal growth and/or mitosis in our study and two studies in S. cerevisiae and D. melanogaster. However, SMC3 is essential in S. cerevisiae. We were unable to knock out SMC3, and a recent mutagenesis study suggests that SMC3 and SMC1 are essential to the parasite during the intraerythrocytic developmental cycle (Zhang et al. Science, 2018). This is why we chose an inducible knockdown system.
*Reviewer #1 (Significance (Required)):
Own opinion The authors provide a basic characterization of the cohesin component SMC3 using NGS methods to investigate chromatin binding sites and its potential influence on gene expression. *
__Response: __We respectfully disagree that our study offers only a basic characterization of SMC3. We combine IFA, mass spectrometry, and both ChIP-seq and RNA-seq of SMC3 across the entire intraerythrocytic developmental cycle to provide the most detailed and comprehensive functional analysis of SMC3 in P. falciparum to date.
The localisation of SMC3 at centromers as described previously (Batugedara 2020) was confirmed. However, the dynamic binding to other regions in the genome, potentially mediated by other proteins, could not be resolved unequivocal with only one replicate of ChIPseq per time point.
__Response: __With regard to the replicates for ChIP-seq, please see our response to this same point above.
Similarly, the RNAseq data demonstrate the relevance of SMC3 for gene expression, but no clear picture of a regulatory mechanism can be drawn at his point. Lacking information about the mode of binding as well as the setup of transcriptome analysis (only two time-shifted sampling points after simultaneous glmS treatment for 96h resulting in incomplete knockdown) cannot definitely elucidate, if SMC3/cohesin is a chromatin factor that affects transcription of genes in general or a specific repressor of stage-specific genes. __Response: __We agree that we have not established a regulatory mechanism for how SMC3 achieves binding specificity. However, the combination of inducible knockdown (as SMC3 is essential to the cell cycle) and differential expression analysis with ChIP-seq from the same time points across the intraerythrocytic developmental cycle is the most stringent and standard approach in the field of epigenetics for determining the direct role of a chromatin-associated protein in gene expression. We provide a detailed explanation of how the transcriptome analysis was set up in the Results (lines 229-234) and Materials and Methods (lines 715-719) section. With regard to our sampling points being “time-shifted,” we provide bioinformatic analysis (line 246-251, what is now Supp. Fig. 5B) of the RNA-seq data from untreated and glucosamine-treated parasites showing highly similar “ages” with regard to progression through the intraerythrocytic developmental cycle. While we of course also monitor progression through the cell cycle with Giemsa staining (Supp. Fig. 5A), this bioinformatic analysis is the most stringent method of determining specific times in the cell cycle.
*The work will be interesting to a general audience, interested in gene regulation and chromatin remodelling
The reviewers are experts in Plasmodium cell biology and epigenetic regulation.*
Reviewer #2
(Evidence, reproducibility and clarity (Required)):
Rosa et al, Review Commons The manuscript by Rosa et al. addresses the function of the cohesion subunit Smc3 in gene regulation during the asexual life cycle of P. falciparum. Cohesin is a conserved protein complex involved in sister chromatin cohesion during mitosis and meiosis in eukaryotic cells. Cohesin also modulates transcription and DNA repair by mediating long range DNA interactions and regulating higher order chromatin structure in mammals and yeast. In P. falciparum, the Cohesin complex remains largely uncharacterized. In this manuscript, the authors present mass spectrometry data from co-IPs showing that Smc3 interacts with Smc1 and a putative Rad21 orthologue (Pf3D7_1440100, consistent with published data from Batugedara et al and Hilliers et al), as well as a putative STAG domain protein orthologue (PF3D7_1456500). Smc3 protein appears to be most abundant in schizonts, but ChIPseq indicates predominant enrichment of Smc3 in centromers in ring and trophozoite stages. In addition, Smc3 dynamically binds with low abundance to other loci across the genome; however, the enrichment is rather marginal and only a single replicate was conducted for each time point making the data interpretation difficult. Conditional knock-down using a GlmS ribozyme approach indicates that parasites with reduced levels of Smc3 have a mild growth advantage, which is only evident after five asexual replication cycles and which the authors attribute to the transcriptional upregulation of invasion-linked genes following Smc3 KD. Indeed, Smc3 seems to be more enriched upstream of genes that are upregulated after Smc3 KD in rings than in downregulated genes, indicating that Smc3/cohesin may have a function in supressing transcription of these schizont specific genes until they are needed. The manuscript is concise and very well written, however it suffers from the lack of experimental replicates for ChIP experiments and a better characterization of the phenotype of conditional KD parasites. * Major comments • In the mass spectrometry analysis, many seemingly irrelevant proteins are identified at similar abundance to the putative rad21 and ssc3 orthologues, and therefore the association with the cohesion complex seems to be based mostly on analogy to other species rather than statistical significance. Hence, it would be really nice to see a validation of the novel STAG domain and Rad21 proteins, for example by Co-IP using double transgenic parasites.*
__Response: __While our IP-MS data did not yield high numbers of peptides, the top most enriched proteins were SMC3 and SMC1. As we state on line 157, two previous studies have already shown a robust interaction between SMC1, SMC3, and RAD21 in Plasmodium, supporting the existence of a conserved cohesin complex. While the identification of the STAG domain-containing protein is interesting, the purpose of our IP-MS was less about redefining the cohesin complex in P. falciparum and more about confirming that the epitope-tagged SMC3 we generated was incorporated correctly into the cohesin complex and was specifically immunoprecipitated by the antibody we later use for western blot, immunofluorescence, and ChIP-seq analyses. However, to validate the results of ours and others’ mass spectrometry results, we generated two new parasite strains – SMC1-3HA-dd and STAG-3HA-dd – and an antibody against SMC3 (see what is now Supp. Fig. 1). We performed co-IP and western blot analysis with these strains and show an interaction between SMC1 and SMC3 and STAG and SMC3 (see what is now Supp. Fig. 2). This information has been added to the manuscript on lines 162-167.
- *The ChIPseq analysis presented here is based on single replicates for each of the three time points. The significance cutoffs for the peaks are rather high (q __Response: __In our experience, a significance cutoff of FDR As we regularly perform ChIP-seq experiments in the lab, we have generated multiple negative control datasets. In our opinion, the most stringent negative control for an HA-tagged protein is performing ChIP with an HA antibody in a WT strain. We have recently published an in-depth analysis of this (and other) negative ChIP-seq controls (Baumgarten & Bryant, 2022, https://doi.org/10.12688/openreseurope.14836.2). We show in this publication that non-specific ChIP-seq experiments (such as negative controls) result in an over-representation of HP1-heterochromatinized regions due to differences in sonication efficiency of heterochromatin and technical challenges with mapping regions with high levels of homology. In the anti-HA in WT ChIP negative control (performed at 12hpi), we do not see any enrichment at centromeric regions, but rather at heterochromatinized regions where clonally variant gene families are located. We performed peak calling analysis and found no significant overlap between the negative control ChIP-seq and the SMC3-3HA ChIP-seq data at 12hpi.
In addition, we have now performed a second biological replicate of the SMC3-3HA ChIP-seq with a different clone at all time points. We compared this data to that from the original clone and found significant overlap of the peaks called (see what is now Table 4 and Supp. Fig. 3A). We generated a stringent list of peaks that were shared between both clones at each time point and repeated all downstream analyses (see what are now Tables 5-8). We found that our conclusions were largely unchanged. Text describing these experiments and analyses have been added throughout the results section.
The SMC3 ChIP from Batugedara et al., 2020 was performed with an in-house generated antibody (not a commercially available, widely validated antibody as we use) at a single time point in the IDC: trophozoites. Batugedara et al. performed one replicate and did not have an input sample for normalization. Rather, it seems that they incubated beads, which were not bound by antibody or IgG, with their chromatin and used any sequenced reads from this beads sample to subtract from their SMC3 ChIP signal as means of normalization. According to ENCODE ChIP-seq standards, this is not a standard nor stringent way of performing ChIP-seq and the subsequent analysis. Because they did not generate a dataset for their ChIP input, it is not possible to call peaks as we do in our study and compare those peaks with ours.
-
The authors argue that during schizogony, cohesin may no longer be required at centromers, explaining the low ChIPsignal at this stage (Line 301). However, during schizogony parasites undergo repeated rounds of DNA replication (S-phase) and mitosis (M-phase) to generate multinucleated parasites; and concentrated spots of Smc3 are observed in each nucleus in schizonts by IFA. In turn, the strong presence of Smc3 at centromers in ring stage parasites is surprising, particularly since the Western Blot in Figure 1D shows most expression of Smc3 in schizonts and least in rings; and Smc3 is undetectable in rings by IFA. Yet, the ChIP signal shows very strong enrichment at centromers, long before S phase produces sister chromatids. What could be the reason for this discrepancy? Again, ChIP replicates and controls would be helpful in distinguishing technical problems with the ChIP from biologically relevant differences. __Response: __We discuss in lines 337-342 not that cohesin is no longer required at centromeres during schizogony, but that its removal from centromeres may be required specifically for separation of sister chromatids, as is seen in other eukaryotes. We also discuss that the unique asynchronous mitosis in Plasmodium may lead to a mixed population of parasites at the time point sampled where there may be some centromeres with SMC3 present and some where it is absent to promote sister chromatid separation. Even though SMC3 may be evicted from centromeres to promote sister chromatid separation, it is likely re-loaded onto centromeres once this process is complete. This is most likely why we see foci of SMC3 in each nucleus of mature schizonts by IFA. With regard to the discrepancy between SMC3 levels in rings seen in total nuclear extracts (by western blot) and at centromeres (by ChIP-seq): the total level of a protein in the nucleus does not necessarily dictate the genome-wide binding pattern or the level of enrichment of that protein at specific loci in the genome. Moreover, if one molecule of SMC3 binds to each centromere, 14 molecules would be needed in a ring stage parasite while over 500 would be needed in a schizont (assuming that there are ~36 merozoites present). SMC3 binds to centromeres in interphase cells in other eukaryotes as well, and we speculate that this binding may play a role in the nuclear organization of centromeres, as we discuss starting on line 333.
-
It is surprising that a conserved protein like Smc3 shows such a subtle phenotype, given that it is predicted to be essential and its orthologues have a function in mitosis. Generally, only limited data are presented to characterize the Smc3 KD parasites, and more detail should be included. For example validation of the parasite line using a PCR screen for integration and absence of wt, parasite morphology after KD, and/or analysis of the KD parasites for cell cycle status. __Response: __First, we have repeated our growth curve analysis several times and with more clones and have concluded that there is not a significant growth phenotype in SMC3 KD parasites (see what is now Supp. Fig. 4B). As we discuss on line 342, very little intact cohesin complex seems to be required for normal growth and mitosis in S. cerevisiae and D. melanogaster, which is probably why we do not see an obvious growth or morphological phenotype. Because we could not generate SMC3 knockout parasites, there may be just enough SMC3 left to perform its vital function in our KD strain. We have added PCR data to demonstrate integration of the 3HA tag- and glmS ribozyme-encoding sequence in the clonal strains we are using for all experiments (see what is now Supp. Fig. 1A). Sanger sequencing was performed on these PCR products to confirm correct sequences. We also added images of Giemsa-stained parasites in untreated and glucosamine-treated parasites at all time points to demonstrate a lack of an obvious morphological phenotype in SMC3 KD parasites (see what is now Supp. Fig. 5A).
-
Synchronization was performed at the beginning of the growth time course, which would be expected to result in a stepwise increase in parasitemia every 48 hours; however, the parasitemia according to Fig. 4F rises steadily, which would indicate that the parasites are actually not very synchronous. __Response: __We did indeed tightly synchronize these parasites and hope that the stepwise increase in parasitemia is seen better in our new growth curve analysis (see what is now Supp. Fig. 4B).
-
The question of whether Smc3 causes a shorter parasite life cycle (quicker progression) or more invasion is important and could be experimentally addressed by purifying synchronous schizont stage parasites and determining their invasion rates as well as morphological examination of the Giemsa smears over the time course. __Response: __We have repeated our growth curve analysis several times and with more clones and have concluded that there is not a significant growth phenotype in SMC3 KD parasites (see what is now Supp. Fig. 4B).
-
Please also compare Smc3 transcriptional levels in transgenic parasites to those in wt parasites to rule out that the genetic modification has lead to artificial upregulation of Smc3 transcription. __Response: __We have added this data to what is now Supp. Fig. 4C, showing that there is no significant difference in SMC3 transcript levels between WT and SMC3-3HA-glmS strains. We have added this information to the text of the manuscript (Line 243). As we also generated an SMC3 antibody, we could demonstrate that there is no appreciable difference in SMC3 protein levels between WT and SMC3-3HA-glmS strains (see what is now Supp. Fig. 1D).
-
According to Figure S2, even more genes were deregulated at the 12 hpi time point in the WT parasites than in Smc3 parasites, and even to a much higher extent. What "transcriptional age" did the WT control parasites have at each time point? __Response: __We have now included the transcriptional age of all strains, replicates, and treatments in what is now Supp. Fig. 5B. At the 12 hpi time point in particular, regardless of glucosamine treatment, the SMC3-3HA-glmS and WT parasites were highly synchronous. The only large discrepancy we see in transcriptional age is between untreated and glucosamine-treated WT parasites at 36 hpi, which is why we did not include this time point in our transcriptional analysis. We were also surprised by the number of genes that were de-regulated with simple glucosamine treatment. The glmS system of inducible knockdown is widely used in P. falciparum; however, to our knowledge, no lab has investigated whether glucosamine treatment affects transcription in wildtype cells over the course of the IDC. Thus, for accurate phenotypic characterization of any protein with this system with regard to transcriptomics, we thought it was important to provide an RNA-seq dataset to define the cohort of genes affected by glucosamine treatment in WT parasites. We hope that our study will demonstrate the importance of using stringent controls when using inducible knockdown systems.
-
A negative correlation with transcription is well established in S. cerevisiae, particularly at inducible genes. How does Smc3 enrichment generally look like for genes that show maximal expression at each of the time point? __Response: __We have performed a metagene analysis of SMC3 enrichment at all genes at each respective time point, which we divided into quartiles of expression based on their FPKM values in the RNA-seq data from the corresponding time point in untreated SMC3-3HA-glmS parasites. This quartile analysis considers all genes, including genes that are not transcribed at all and regardless of whether a gene has a significant SMC3 peak or is differentially expressed upon SMC3 knockdown. At the 12 hpi time point, we do see an inverse correlation between SMC3 enrichment and gene transcription level, but this enrichment is most pronounced across genes bodies. We see the highest SMC3 enrichment at genes in the 4th (lowest) quartile category. For the other two time points, we do not see any obvious pattern of SMC3 enrichment with regard to transcriptional status.
-
Line 590: according to the methods, a 36 hpi KD time point was also harvested. Why are the data not shown/analysed? __Response: __To be clear, we are most interested in the transcriptional role of SMC3 during interphase, where results are not confounded by its potential role in mitosis. However, we did collect a 36hpi time point in the SMC3-3HA-glmS and WT strain, with and without glucosamine. We have added this last time point and the WT data from the other two time points to the manuscript (see Tables 11-13). Unfortunately, and for reasons unknown, the WT replicates treated with glucosamine showed a significantly advanced “transcriptional age” compared to the other replicates at 36hpi (see what is now Supp. Fig. 5B). Thus, we did not feel comfortable performing the RNA-seq analysis as we did with the other two time points (i.e. subtracting up- and down-regulated genes from the WT control from the SMC3-3HA-glmS data sets). We have added this information to the results section (Lines 256 and 261). As the WT parasites treated with glucosamine were approximately 8 hours in advance of the untreated WT parasites for the 36hpi time point, any up- and down-regulated genes might have been due to differences in the cell cycle rather than due to glucosamine treatment. The glmS system of inducible knockdown is widely used in P. falciparum; however, to our knowledge, no lab has investigated whether glucosamine treatment affects transcription in wildtype cells over the course of the IDC. Thus, for accurate phenotypic characterization of any protein with this system with regard to transcriptomics, we thought it was important to provide an RNA-seq dataset to define the cohort of genes affected by glucosamine treatment in WT parasites. We hope that our study will demonstrate the importance of using stringent controls when using inducible knockdown systems.
Minor Comments • Line 103/104: the hinge domain and ATPase head domain are mentioned, please annotate these in Figure 1A.
__Response: __We have annotated the hinge and ATPase domains.
-
Figure 1D: the kDa scale is missing from the H3 WB. __Response: __We have added a kDa scale.
-
What is the scale indicated by different colors in Fig. 2A? __Response: __The different colors (blue, coral, and green) only represent the 12, 24, and 36hpi time points, respectively. This color scheme is used throughout the manuscript. If the reviewer is referring to the color gradation within each circos plot, this does not indicate a specific scale. The maximum y-axis value for all circos plots is 24, as indicated in the figure legend.
-
Line 189: it would also be interesting how many peaks are "conserved" between the different time points studied, so not only to compare the gene lists of closest genes but also the intersecting peaks and then the closest genes to the intersecting peaks. __Response: __We have added this information in Table 7 and in the manuscript starting on Line 203. Using the new dataset of consensus peaks between two replicates, there were 88 genes associated with an SMC3 peak across all three time points, most of which were close to a centromeric region.
-
What is the distribution of the peaks over diverse genetic elements, such as gene bodies, introns, convergent/ divergent/ tandem intergenic regions? In yeast, cohesion is particularly enriched in convergent intergenic regions, so it would be interesting to see how this behaves in P. falciparum. __Response: __We would have liked to define how many peaks were in intergenic versus genic regions of the genome, but the dataset of “genes” from PlasmoDB includes UTRs. Thus, we would need a better annotation of the genome to perform this analysis. Regardless, we calculated the average SMC3 peak enrichment (shared between both replicates) in intergenic regions between convergent and divergent genes (see what is now Supp. Fig. 3B and Table 6). As we now state in the manuscript on line 198, we see a slight enrichment in regions between convergent genes at all time points, but the differences were not significant.
-
Line 130 intra-chromosomal interactions (word missing) __Response: __Thank you for pointing this out. We have corrected this.
-
Contrary to Figure 1D, the WB in Figure 3A indicates strong expression of Smc3 in rings. Please comment on this discrepancy. __Response: __While extracts from all time points were run on the same western blot in Fig. 1D and thus developed for the same amount of time, this was not the case for Fig. 3A. In Fig. 3A, the samples were run on different blots and exposed for different times, so while we can compare SMC3-HA levels between – and + glucosamine for each time point, the levels at 12 hpi cannot be quantitatively compared to those at 24 or 36hpi.
-
What time point after glucosamine addition represents the WB in Fig. 3A? __Response: __The “12hpi” parasites were sampled approximately 108 hours post glucosamine addition and the “24hpi” parasites sampled approximately 120 hours post glucosamine addition. Basically, the parasites were treated with glucosamine for 96 hours, synchronized, and then harvested 12 and 24 hours later.
-
Line 233 / Suppl Figure 3: Isn't it a bit concerning that the untreated control parasites at 24 hpi statistically corresponded to 18-19 hpi? And to what timepoint did the wt parasites correspond? __Response: __We are not concerned by this, and we have included the WT parasites in what is now Supp. Fig. 5B for better comparison. In the analysis presented in Supp. Fig. 5B, regardless of glucosamine presence or absence, the differences among replicates and strains at 12 and 24hpi are, in our opinion, minimal, amounting to one or two hours of the 48-hour IDC. In our extensive experience with RNA-seq across the P. falciparum lDC, this synchronization is extremely tight. As we describe on line 430 of the Materials and Methods, there is a ±3 hour window in our synchronization method, meaning that parasites harvested at 24hpi could be anywhere from 21-27hpi. In addition, the dataset that was used for comparison (from Bozdech et al., 2003) was generated in 2003 in a different laboratory using different strains with microarray. While comparing more recent RNA-seq data to this classic study has become well-established practice and is useful for comparing transcriptional age between replicates and strains, it is inevitable that the calculated “hpi” from (Bozdech et al., 2003) will differ somewhat from our experimental “hpi”. We have indeed seen this small discrepancy in predicted transcriptional age in several of our RNA-seq datasets (unrelated to this study) from trophozoites harvested at 24hpi.
-
Line 264: "whether naturally or via knockdown" - the meaning of this sentence is not entirely clear __Response: __We are referring to depletion of SMC3 at promoters, either naturally (i.e. lack of binding at the promoter at 36hpi that is not the result of SMC3 knockdown, as we show in Fig. 4B) or via SMC3 knockdown, which is not natural but artificial.
-
Figure 4 Legend: A, B, C etc. are mixed up. Response: Thank you for pointing this out. We have corrected this.
-
Figure 4D: the differences seem to be marginally significant, even not significant at all (q=0.8) for gap45 at 12hpi. __Response: __If one defines a significance cutoff of q = 0.05 (as is common practice in differential expression analyses), then the differences are significant. For a small minority of invasion genes (such as gap45), we do observe significance at either 12 hpi or 24 hpi, but not both. Thus, we have removed the word “significant” from the descriptions of each dataset in Tab 1 of what is now Table 18. however, we do not believe that this rules out a role for SMC3 at such a gene during interphase. What is now Table 18 offers a longer list of invasion-related genes, most of which are more “significantly” affected than rap2 and gap45.
-
Figure 4F shows FACS data using SYBR green as a DNA stain. The authors could exploit this data to look at the relative DNA content per cell as a measure of parasite stage, since more mature parasites will have more DNA (mean fluorescence intensity). How did the corresponding parasite cultures look in Giemsa smears? Response: We have repeated our growth curve analysis several times and with more clones and have concluded that there is not a significant growth phenotype in SMC3 KD parasites (see what is now Supp. Fig. 4B). We have added images of Giemsa-stained parasites in untreated and glucosamine-treated parasites at all time points to demonstrate a lack of an obvious morphological phenotype in SMC3 KD parasites (see what is now Supp. Fig. 5A).
-
Are RNAseq replicates biological replicates from independent experiments or technical replicates? __Response: __RNA-seq replicates are technical replicates from the same parasite clone.
-
Why does the number of genes analysed for differential gene expression differ between the comparisons? __Response: __If the reviewer is referring to the discrepancy between the total number of genes for different time points [for example, between what is now Table 9 (12hpi) and Table 10 (24hpi)], this is because in the RNA-seq/differential expression analysis, there have to be reads mapping back to a gene in order for that gene to be included in the analysis. Thus, if a gene is not transcribed at a given time point in the treated or untreated samples, it will not be included in the analysis. Gene transcription fluctuates significantly over the course of the IDC, so different time points will have different total numbers of transcribed genes.
-
Line 372: Do you mean the proteins or the genes? AP2-I has a peak at 24 hpi and 36 hpi, and its interacting AP2 factor Pf3D7_0613800 at all time points. __Response: __We are referring to the genes. With the new ChIP-seq analysis including the second replicate, there are no consensus SMC3 peaks associated with ap2-I, bdp1, or Pf3D7_0613800 (see what is now Table 7).
-
Line 480: no aldolase was shown. __Response: __We have removed this sentence.
-
Line 838: include GO analysis in methods __Response: __We have added this.
Reviewer #2 (Significance (Required)): The paper addresses the function of the cohesin complex in gene regulation of malaria parasites for the first time. Due to the conserved nature of the complex, the data may be interesting for a broad audience of scientists interested in nuclear biology and cell division/ gene regulation.
Reviewer #3
(Evidence, reproducibility and clarity (Required)):
*Summary:
In the presented manuscript by Rosa et al. the authors investigate the longstanding question of how P. falciparum achieves the tight transcriptional regulation of its genome despite the apparent absence of many canonical sequence specific transcription factor families found in other eukaryotes. To do this the authors investigate the role of the spatial organization of the genome in this context, by performing a functional characterization of the conserved cohesion subunit SMC3 and its putative role in transcriptional regulation in P. falciparum. Using Cas9 mediated genome editing the authors generated a SMC3-3xHA-glmS parasite line, which they subsequently used to show expression of the protein over the asexual replication cycle by western blot and IFA analysis. In addition, using co-IP experiments coupled with mass spectrometry they identified the additional components of the cohesion complex also found in other eukaryotes as interaction partners of SMC3 in the parasite, thereby confirming the presence of the conserved cohesin complex in P. falciparum. By using a combination of ChIP-seq and RNA-seq experiments in SMC3 knockdown parasites the authors furthermore show that a reduction of SMC3 resulted in the up-regulation of a specific set of genes involved in invasion and egress in the early stages of the asexual replication cycle and that this up-regulation in transcription is correlated with a loss of SMC3 enrichment at these genes. From these observations the authors conclude, that SMC3 binds dynamically to a subset of genes and works as a transcriptional repressor, ensuring the timely expression of the bound genes. Overall, the presented data is intriguing, of high quality and very well presented. However, there are some points, which should be addressed to bolster the conclusions drawn by the authors.
Major points: I was not able to find the deposited datasets in the BioProject database under the given accession number. This should obviously be addressed and would have been nice to be able to have a look at these datasets also for the review process. *__Response: __We apologize for not giving the reviewers access. As the manuscript has been made available as a pre-print (which includes data accession numbers), but has not yet been published, we have not activated access to the data on the database.
*SMC3-ChIP-seq experiments:
"168 were bound by SMC3 across all three time points (Fig. 2D). However, most SMC3-bound genes showed a dynamic binding pattern, with a peak present at only one or two time points (Fig. 2B,D)."
Here it would be interesting to actually have more than one replicate of each of these ChIP-seq time points. This could provide a better idea of how "dynamic" these binding patterns actually are. Furthermore, I was missing a list of these 168 genes, which are constantly bound by SMC3. Anything special about those? What actually happens to this subset of genes in the SMC3 knockdown parasites? Do they show similar transcriptional changes?*
__Response: __We have now performed a second biological replicate of the SMC3-3HA ChIP-seq with a different clone at all time points. We compared this data to that from the original clone and found significant overlap of the peaks called (see what is now Table 4 and Supp. Fig. 3A). We generated a stringent list of peaks that were shared between both clones at each time point and repeated all downstream analyses (see what are now Tables 5-8). We found that our conclusions were largely unchanged. Text describing these experiments and analyses have been added throughout the results section. Using the new dataset of consensus peaks between two replicates, there were 88 genes associated with an SMC3 peak across all three time points (see what is now Table 7). The genes that are associated with an SMC3 peak at all time points are, in general, those closest to centromeric/pericentromeric regions and show no obvious functional relationship to each other. Out of these 88 genes, four are significantly up- or downregulated at 12 hpi and 26 are significantly up- or downregulated at 24 hpi. The most significantly downregulated of these genes in both datasets is smc3 itself.
*SMC3-knockdown experiments:
In Sup. Fig. 1 there is a double band in the HA-western blot in the 2nd cycle -GlcN. sample. This second band is absent in all other HA-western shown. Have the authors any idea where that second band comes from?*
__Response: __As the reviewer says, we do not see this second band in most of our western blots. It is possible that it is just a small amount of degradation in the lysate.
In Figure 3A, the WB data shown is slightly contrasting the RNA-seq quantification (3B). The knock-down on protein level seems to be stronger in the 12 hpi samples here than in the 24 hpi samples. Although the band for HA-SMC3 is stronger at the 12 hpi TP there's no band visible in the + GlcN. sample. There's however in the 24 hpi samples. Could the authors comment on this?
Response: __With regard to the discrepancy of the knockdown and protein versus RNA level, it is quite common for transcript levels to not agree with protein levels. This is why we always confirm a transcriptional knockdown with western blot analysis using appropriate loading controls. We are not sure why there is a more dramatic knockdown of SMC3 at 12hpi than at 24hpi, as these samples came from the same culture, but were simply harvested 12 hours apart. __
*"Comparison of our RNA-seq data to the time course transcriptomics data from (Painter et al., 2018) revealed that SMC3 depletion at 12 hpi caused downregulation of genes that normally reach their peak expression in the trophozoite stage (18-30 hpi), with the majority of upregulated genes normally reaching their peak expression in the schizont and very early ring stages (40-2 hpi) (Fig. 3E). At 24 hpi, a similar trend is observed, with most downregulated genes normally peaking in expression in trophozoite stage (24-32 hpi) and the majority of upregulated genes peaking in expression at very early ring stage (2 hpi) (Fig. 3F)."
I'm not fully convinced by these presented results/conclusions. This dataset would greatly benefit from the inclusion of additional later time points.*
__Response: __To be clear, we are most interested in the transcriptional role of SMC3 during interphase, where results are not confounded by its potential role in mitosis. However, we did collect a 36hpi time point in the SMC3-3HA-glmS and WT strain, with and without glucosamine. We have added this last time point and the WT data from the other two time points to the manuscript (see Tables 11-13). Unfortunately, and for reasons unknown, the WT replicates treated with glucosamine showed a significantly advanced “transcriptional age” compared to the other replicates at 36hpi (see what is now Supp. Fig. 5B). Thus, we did not feel comfortable performing the RNA-seq analysis as we did with the other two time points (i.e. subtracting up- and down-regulated genes from the WT control from the SMC3-3HA-glmS data sets). We have added this information to the results section (Lines 256 and 261). As the WT parasites treated with glucosamine were approximately 8 hours in advance of the untreated WT parasites for the 36hpi time point, any up- and down-regulated genes might have been due to differences in the cell cycle rather than due to glucosamine treatment. The glmS system of inducible knockdown is widely used in P. falciparum; however, to our knowledge, no lab has investigated whether glucosamine treatment affects transcription in wildtype cells over the course of the IDC. Thus, for accurate phenotypic characterization of any protein with this system with regard to transcriptomics, we thought it was important to provide an RNA-seq dataset to define the cohort of genes affected by glucosamine treatment in WT parasites. We hope that our study will demonstrate the importance of using stringent controls when using inducible knockdown systems.
We performed differential expression analysis of the SMC3-3HA-glmS parasites with and without glucosamine at 36hpi (we have added this data in Table 11). Again, significantly up- and down-regulated genes were not filtered using the WT dataset. With this analysis, we see only three genes from the list of invasion-related genes (Hu et al., 2010) that are up-regulated, but none of them have a significant q-value (Tab 5 of Table 18). Thus, depletion of SMC3 in late stage parasites does not lead to up-regulation of the same genes that are upregulated at 12 and 24hpi. We have added this information to the text (Line 277).
*The presented upregulation of the egress and invasion related genes is hard to pinpoint to be a direct effect of transcriptional changes due to the SMC3 knockdown. While there's a slight upregulation of these genes they still seem to be regulated in their normal overall transcriptional program as shown in Figure 4D/E. *
__Response: __We provide evidence of a direct effect of SMC3 binding by combining differential expression analysis performed upon SMC3 knockdown with SMC3 ChIP-seq at corresponding time points. As we show in what is now Fig. 4C and D, promoter accessibility of these egress/invasion genes correlates with their transcriptional activity. However, SMC3 binding to the promoters of these same genes shows inverse correlation with their transcriptional activity (what is now Fig. 4B and D). While we believe that SMC3 does contribute to the repression of these genes at specific time points during the cell cycle, it is highly likely that SMC3 is just one protein of many that regulates these genes. Moreover, and especially since we do not see a growth phenotype in the SMC3 KD, it is possible that another protein or even SMC1 could compensate for loss of SMC3 at these promoter regions. We now state these possibilities on lines 346 383 of the Discussion.
*So the changes could in theory also be explained by the differences in cell cycle progression which are present between +/- GlcN. cultures (Sup. Fig. 3). The presented normalization to the microarray data is a well-established practice to correct for this but, as presented seems to have its limitation with these parasite lines (line 233, glucosamine treated parasites harvested at 24 hpi correspond statistically to approximately 18-19 hpi (Supp. Fig. 3).) *
__Response: __In the analysis presented in what is now Supp. Fig. 5B, regardless of glucosamine presence or absence, the differences among replicates and strains at 12 and 24hpi are, in our opinion, minimal, amounting to one or two hours of the 48-hour IDC. In our extensive experience with RNA-seq across the P. falciparum lDC, this synchronization is extremely tight. As we describe on lines 416-421 of the Materials and Methods, there is a ±3 hour window in our synchronization method, meaning that parasites harvested at 24hpi could be anywhere from 21-27hpi. In addition, the dataset that was used for comparison (from Bozdech et al., 2003) was generated in 2003 in a different laboratory using different strains with microarray. While comparing more recent RNA-seq data to this classic study has become well-established practice and is useful for comparing transcriptional age between replicates and strains, it is inevitable that the calculated “hpi” from (Bozdech et al., 2003) will differ somewhat from our experimental “hpi”. We have indeed seen this small discrepancy in predicted transcriptional age in several of our RNA-seq datasets from trophozoites harvested at 24hpi.
By including additional later time points, one could actually follow the expression profiles over the whole cycle and elucidate if there's an actual transcriptional up-regulation of the genes, or if the + GlcN. parasites show a faster cell cycle progression, with a shifted peak expression timing compared to the - GlcN. parasites. __Response: __We did collect a 36hpi time point in the SMC3-3HA-glmS and WT strain, with and without glucosamine. We have added this last time point and the WT data from the other two time points to what is now Supp. Fig. 5. Unfortunately, and for reasons unknown, the WT replicates treated with glucosamine showed a significantly advanced “transcriptional age” compared to the other replicates at 36hpi. Thus, we did not feel comfortable performing the RNA-seq analysis as we did with the other two time points (i.e. subtracting up- and down-regulated genes from the WT control from the SMC3-3HA-glmS data sets). We have added this information to the results section (Lines 256 and 261). As the WT parasites treated with glucosamine were approximately 8 hours in advance of the untreated WT parasites for the 36hpi time point, any up- and down-regulated genes might have been due to differences in the cell cycle rather than due to glucosamine treatment. The glmS system of inducible knockdown is widely used in P. falciparum; however, to our knowledge, no lab has investigated whether glucosamine treatment affects transcription in wildtype cells over the course of the IDC. Thus, for accurate phenotypic characterization of any protein with this system with regard to transcriptomics, we thought it was important to provide an RNA-seq dataset to define the cohort of genes affected by glucosamine treatment in WT parasites. We hope that our study will demonstrate the importance of using stringent controls when using inducible knockdown systems.
*"These genes show SMC3 enrichment at their promoter regions at 12 and 24 hpi, but not at 36 hpi (Fig. 4C), and depletion of SMC3 resulted in upregulation at both 12 and 24 hpi (Fig. 4D). Comparison of the SMC3 ChIP-seq data with published Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data (Toenhake et al., 2018) and mRNA dynamics data (Painter et al., 2018) from similar time points in the IDC revealed that SMC3 binding at the promoter regions of these genes inversely correlates with chromatin accessibility (Fig. 4C) and their mRNA levels (Fig. 4E), which both peak in schizont stages. These data are consistent with a role of SMC3 in repressing this gene subset until their appropriate time of expression in the IDC."
The presented correlations certainly make an intriguing point towards the authors conclusion that SMC3/cohesin depletion from the promoter regions of the genes results in a de-repression of these genes and their transcriptional activation. However, the SMC3 knockdown is not complete and only up to 69% as presented on RNA level in these parasites. Therefore a control experiment which needs to be done is to actually show the loss of SMC3 from the presented activated example genes in the knockdown parasites. This could easily be done by ChIP-qPCR or even ChIP-seq, to get a global picture of the actual changes in SMC3 occupation in the knockdown parasites in correlation with changes in transcript levels. *__Response: __While SMC3-3HA-glmS knockdown is not complete at the RNA level, it is fairly robust at the protein level, especially at 12hpi (Fig. 3A).
*"These data suggest that SMC3 knockdown results in a faster progression through the cell cycle or a higher rate of egress/invasion."
The authors could greatly strengthen their conclusions by investigating this thoroughly. Pinpointing the observed phenotype to an actual increase in invasion or egress would add to the authors main conclusion that the loss of SMC3 de-regulates the timing of gene expression for these invasion related genes thereby increasing their transcript levels and thus leading to a higher rate of egress/invasion. To determine cell cycle progression simple comparisons between DNA content using a flow cytometer at timepoints together with visual inspection of Giemsa stained blood smears would give a ggod indication towards changes in cell cycle progression. In addition invasion/egress assays by counting newly invaded rings per schizont could reveal, if there are changes in the rate of egress/invasion upon SMC3 knockdown.*
Response: __We have repeated our growth curve analysis several times and with more clones and have concluded that there is not a significant growth phenotype in SMC3 KD parasites (see what is now Supp. Fig. 4B). We have added images of Giemsa-stained parasites from the knockdown time course we performed to what is now Supp. Fig. 5A. We see no obvious differences in cell morphology caused by glucosamine treatment in the WT or SMC3-3HA-glmS parasites. As we discuss on line 327, very little intact cohesin complex seems to be required for normal growth and mitosis in S. cerevisiae and D. melanogaster, which is probably why we do not see an obvious growth or morphological phenotype. We believe that SMC3 is probably only a part of a complex controlling transcription of these invasion or egress genes. Thus, the up-regulation of these genes upon SMC3 KD might not be enough to lead to a significant growth or invasion phenotype. __
*Minor points:
In the MM section on the Cas9 experiments it says dCas9 where it should be Cas9 (line 425)*
__Response: __Thank you for pointing this out. We have corrected this.
It would be great to add which HP1 antibody was used in which dilution in the IFAs to the MM section. __Response: __We have added this information to the Materials and Methods section.
In Figure 4C for the gap45 gene there's is some green peak floating around which should not be there. __Response: __Thank you for pointing this out, we have corrected it.
*Reviewer #3 (Significance (Required)):
Significance: The manuscript investigates a very timely topic by trying to uncover new molecular mechanisms of transcriptional regulation in P. falciparum. Investigating the role of the cohesin complex/SMC3 in this context provides valuable new insights to the field. While the first part with the description of the SMC3 cell line and the co-IP experiments largely confirms published data on the existence and composition of the cohesin complex in Plasmodium and its enrichment at the centromeres, the second part is especially intriguing since it investigates the molecular function of SMC3 in more detail. The results pointing to a role of SMC3/cohesin as a transcriptional repressor are of great interest to the field and will open up new concepts for future investigation.*
*Audience: The work is particularly interesting for people interested in gene regulatory processes in Plasmodium and Apicomplexan parasites in general. At the same time it also nicely points towards shared principles of gene regulation to other eukaryotes in relation to the spatial organization of the genome making the work also very interesting for a broader audience with interest in the general principles of gene regulatory processes in eukaryotic organisms.
Expertise: P. falciparum epignetics and chromatin biology / gene regulation / Cas9 gene editing*
CROSS-CONSULTATION COMMENTS
All reviewers agree that the paper addresses an important topic and provides convincing evidence for enrichment of the cohesin component Smc3 at P. falciparum centromers. In contrast, evidence for a function of Smc3 as a transcriptional repressor of genes in the first part of the parasite life cycle is less well supported. All reviewers agree that the statistical significance of the ChIP experiments needs to be impoved by including biological replicates. In addition, the phenotype of the conditional knock-down should be analysed in more detail by clarifying whether faster cell cycle progression or higher invasion rate are responsible for the observed growth adavantage. Inclusion of transcriptional data from a later time point in addition to the presented data for 12 hpi and 24 hpi was also requested by all reviewers. Finally, several inconsistencies require clarification.
-
- Dec 2022
-
www.biorxiv.org www.biorxiv.org
-
Note: This rebuttal was posted by the corresponding author to Review Commons. Content has not been altered except for formatting.
Learn more at Review Commons
Reply to the reviewers
Reviewer 1
Although this is an interesting, and generally well-performed study, it is primarily observational and there are few mechanistic insights provided into how MUC13 modulates barrier function. The authors propose a presumably direct interaction between MUC13 and PKC, which apparently sequesters PKC, preventing this kinase from triggering PKC-dependent increases in TJ barrier function; however, there is no evidence that a MUC13-PKC interaction occurs, that MUC13 is phosphorylated by PKC, or that phosphorylation of MUC13 has any impact on its function or overall barrier function. Thus, the hypothesis is not directly tested and all observations in this manuscript are generally correlative in nature.
While the MUC13 cytoplasmic tail contains a putative PKC-binding motif, we indeed do not show a direct interaction between MUC13 and a member of the PKC family in this manuscript. Unfortunately, we have so far not been able to successfully perform (co-)immunoprecipitation of MUC13 with our current anti-MUC13 antibodies.
To provide more insights into the possible MUC13-PKC interaction, we plan to perform several experiments.
- First, we will determine the expression levels of the different PKC isotypes (PKC alpha, beta, gamma, delta, epsilon, and zeta) in the HRT18 cell lines by western blot.
- Next, we will determine the localization of the relevant PKC isoforms and MUC13 by immunofluorescence microscopy. We are curious to see if we can find a colocalization between MUC13 and a PKC member on the lateral or apical membrane. If we can demonstrate a colocalization, we could follow up with a proximity ligation assay, but this would require the MUC13 antibody directed against the cytoplasmic tail (which only detects the lateral population) and might therefore be challenging.
- Furthermore, since PKC delta protein levels were upregulated in the total lysate of ∆MUC13 cells, we will test a PKC delta-specific inhibitor in the TEER assay.
Consider quantifying all blots (Fig. 5C, Fig. 6B).
As suggested, we will quantify both blots.
Consider using dot-plots for all quantified data.
The graphs will be altered to include individual measurement points.
Reviewer 2
Fig2E showed two bands with different size in the two MUC13 WT control cell lines. They hypothesized that this could be the consequences of glycosylation different patterns. A sample with untransfected HRT18 might be included in the western blot panel. Additionally, what is the 100kDa band?
Mucin blots are notoriously difficult and these MUC13 blots are the result of a lot of trial and error. We repeated the Western Blot with original HRT18 cells, HRT18 original cell line, as well as the two CRISPR control cells used in the study (WT 1 and WT 2) and one of the full-length MUC13 knockout cells. The higher band was absent from the MUC13 knockout cells, but a small shift in the MUC13 band size can be noted in the WT 1 cells compared to the original and the WT 2 cell lines, possibly indicating a change in the glycosylation pattern. The 100 kDa band remains detectable in all cell lines including the ∆MUC13 cell line, therefore we consider this to be an aspecific background band of the MUC13 antibody. We will add a more extensive Western Blot analysis to the manuscript.
Did the transfection of the inducible GFP-MUC13 plasmid induce any decrease of Claudin1/3/4 in HRT18 or Caco2 cells? Same question regarding PKCdelta.
These are indeed interesting questions. We will perform these experiments with our MUC13-overexpression HRT18 cells.
Reviewer 3
Moreover, the authors should determine if MUC13∆CT localize to TJs, as suggested by the working model in Figure 7C. The subcellular localization of MUC3∆CT could give critical clues for its function, but Figure 2G fails to provide any information and the authors do not present any additional data concerning the localization of MUC13∆CT. Detection of MUC13 in membrane fractions of WT, MUC13∆CT and cells lacking the mucin domain could be a feasible strategy forward.
We will perform additional immunofluorescence experiments to determine the subcellular localization of MUC13-∆CT more accurately. However, detection of the extracellular domain by western blot, as suggested, is not possible due to the incompatibility of the extracellular MUC13-directed hybridoma antibody with the western blot technique. We currently do not have a suitable antibody that recognizes the ED and can be used for western blot.
The authors introduce an inducible MUC13-GFP fusion protein into WT and ∆MUC13 cells and show that it reverses the enhanced TEER upon MUC13 deletion. Unfortunately, the "Materials and Methods" section lacks adequate information on how this fusion protein was designed. Critical questions are the position of the GFP tag within MUC13, whether the fusion protein is correctly processed in HRT18 cells, and if it localizes to the apical or apico-lateral membrane domains? Figure 2H is of low magnification and fails to provide information on the subcellular localization of the MUC13-GFP fusion protein.
The materials and methods section will be adjusted to describe all the design details of the fusion protein. The GFP tag was added to the MUC13 C-terminus with a GGGS linker sequence in between. Processing of the fusion protein seems correct as we observed MUC13-GFP localization to both lateral and apical membranes and no access intracellular build up. As suggested by the reviewer, we will add more detailed immunofluorescence pictures to the manuscript.
Figures 6B-C suggest that PKCdelta levels increase in ∆MUC13 cells, which correlates with higher enrichment of Claudins in membrane fractions. The authors then inhibited PKCdelta and observed reduced recruitment of Claudins to membrane fractions. Since the family of Claudins are differentially regulated by phosphorylation (PMID: 29186552), the authors should investigate the TEER phenotype of WT, ∆MUC13 and MUC13∆CT upon PKC inhibition.
We must clarify that figures 6C-D are done using the PKC inhibitor targeting all conventional PKCs (alpha, beta, gamma) as well as delta (https://www.tocris.com/products/gf-109203x_0741). We recently obtained a PKCdelta-specific inhibitor which we will test in the TEER build-up experiments.
Moreover, the authors predict phosphorylation sites in MUC13CT and suggest a link between PKC and MUC13 (Figure. 6A), however no evidence is presented to support this hypothesis. The authors should either determine if PKC phosphorylates MUC13 and if this modification has implication on MUC13 localization and TJ function, or remove statements regarding MUC13 phosphorylation. The data provided suggest that PKC regulates TJ proteins independent of MUC13.
We will adjust the manuscript to put less emphasis on the putative PKC motifs in the MUC13 cytoplasmic tail. For further details on how we will proceed regarding the possible MUC13-PKC interaction see question 1 from reviewer #1.
Figure 5C. Quantification of at least 3 independent experiments is required.
These data will be added to the manuscript.
Figure 6B. Quantification of at least 3 independent experiments is required.
These data will be added to the manuscript.
Reviewer 4
OPTIONAL: MUC13 is expressed both, in the basolateral membranes and in the apical membrane of intestinal epithelial cells (IECs). Does the authors check the relevance of MUC13 in the formation of microvilli in IECs? Are microvilli different (microvilli staining, number of positive cells to microvilli, length, width or distribution of microvilli) in ΔMUC13 and in MUC13-ΔCT? How the glycocalyx looks like in these cells genetically modified for MUC13?
HRT18 cells do not seem to develop microvilli. However, we plan to stain these cells with a microvilli-specific antibody (ACTUB). The HRT18 cells express mostly MUC13 and relatively low levels of the larger TM mucin MUC1. To study changes in the glycocalyx, we will stain using a MAL-II antibody which targets α-2,3 sialic acids, which are abundantly present in mucins. In this way, we will determine any big changes in the total glycocalyx that may occur in response to the removal of MUC13.
In the figure 1D would be nice to represent the co-localization of MUC13 together with occluding in a graph in each Z-stack so you can visualize in which part of the cell is maximum colocalization of these both components.
These data will be provided.
In the figure 1E, would be great to compare between the two different MUC13 antibodies the apical fraction stained in HRT18 and Caco-2. Specially in the HRT18 cell line since the first antibody did not label apical MUC13 expression meanwhile the second antibody detects the apical expression in these cells. How much lateral lateral stain the C terminal antibody compare with the extracellular antibody for MUC13 and how much stain apically the C terminal antibody compare with the extracellular antibody? Would be nice to see some comparative results using the intensity by Z-stack and plotting in a graph.
This is a good suggestion as it is quite intriguing that both MUC13 antibodies seem to target (partially) different MUC13 populations. We will perform co-staining with both MUC13 antibodies to provide information on which MUC13 populations are detected by each antibody (apical vs lateral membrane).
Manuscript would be improved if in the figure 2H to compare within the same cell line the number of MUC13 positive cells in the WT, number of MUC13 positive cells in WT+pMUC13 and the number of MUC13 positive cells in the ΔMUC13+pMUC13
We will quantify the percentage of MUC13-GFP positive cells in both the WT and ΔMUC13 backgrounds by either microscopy or flow cytometry.
In figure 5C would be helpful to plot in a graph the normalized expression of each TJ protein and compare between the different cells used (WT, ΔMUC13 and MUC13-ΔCT) as you did in figure 5A
We will provide the quantification data of three independent experiments.
Description of the revisions that have already been incorporated in the transferred manuscript
Reviewer 1
In addition, this model does not explain why all kinase inhibitors tested reverse the increase in TER observed in deltaMUC13 cell lines. Does this reflect the lack of inhibitor specificity or the likelihood that many kinases are involved?
As stated in the manuscript, we think that MLCK, ROCK, and PKC are all essential for TER buildup in the ∆MUC13 cells. Because the roles of MLCK and ROCK are well established, we choose to follow up on the PKC results. We adjusted the text to clarify this point.
The authors do observe that there is an increase in expression of several tight junction-associated proteins, including the claudins, in deltaMUC13 cells. Affected CLDNs include 1, 2, 3, 4, 7, 12. (1) While it appears the authors are arguing that this increased claudin expression results in increased barrier function, they do not sufficiently highlight the well-known role that CLDN2 has in cation transport, and both CLDN-4 and -7 have also been implicated in paracellular ion flux (although this is apparently cell-type specific). These observations would seem to argue against a simple correlation between claudin expression and tight junction barrier function.
The reviewer is right about the different functions of claudins. Claudin-2, -4 and -7 have (potentially) pore-forming properties, while the other claudins restrict paracellular passage. It has been previously demonstrated that the magnitude of paracellular ion and water flux is reflected by the specific repertoire of claudin family members (Shashikanth et al., 2022). In this paper, overexpression of claudin-4 was shown to mobilize and affect polymeric strands of claudin-2, thus blocking its channel activity. Our mass spectrometry data demonstrated a striking increase in claudin-1, -2, -3, -4, -7, and -12 in the MUC13 knockout membranes compared to WT. We hypothesize that the claudin repertoire in the MUC13 knockout cells leads to a more restricted paracellular route (as observed in the TEER and tracer experiments). The pore-forming claudins may be subject to “interclaudin interference” therefore leading to restriction of the total paracellular ion and water flux. We have adjusted the text of the manuscript to clarify this point.
We attempted to investigate claudin-2 expression levels in isolated membranes by Western Blot but were unsuccessful as the antibody did not detect any protein while claudins-1 and -4 could be detected with the same method.
Furthermore, the authors should note the disconnect between paracellular ion flux mediated by claudins and the flux of markers such as dextrans and lucifer yellow, which can be dissociated from claudin function.
We acknowledge that the flux of larger particles (the leak pathway) is not regulated by claudins (which regulates the pore pathway). We aimed to assess both the pore and the leak paracellular pathways, by using different techniques including TEER, small solutes (Lucifer Yellow CH), and larger molecules (4 and 70 kDa FITC-Dextrans). HRT18 wild type cells are already very restrictive to the pass of larger molecules (FITC-Dextrans) but are more permeable to smaller solutes such as Lucifer Yellow (400 Da). We observed that removal of the MUC13 cytoplasmic tail did not affect the TEER, but reduced the paracellular passage of Lucifer Yellow, demonstrating that manipulation of MUC13 can affect both the pore and leak pathways. We adjusted to text to include this point.
The increased expression of claudins in the nominally tail-minus MUC13 without a corresponding change in TER would again seem to argue against a simple correlation;
MUC13-dCT cells showed consistently increased levels of claudins-1 and -2, but not the other claudins. This claudin repertoire (with high claudins-1 and -2, but lower claudin-3, -4, -7, and -12) is apparently not enough to increase TEER. We think that this again reflects the importance of the total claudin composition for the control of the paracellular pathway.
Watch the use of decimal points instead of commas (lines 253 and 256).
Corrected.
Line 543: MilliQ is not a washing agent (or is it?). (Line 535) We use MilliQ as a final step before mounting the glass slides to remove any possible salt deposition that would affect the visualization by microscopy.
We have specified this in the text.
Line 553: TER is the product of total resistance times the area. The units are ohms times area.
Indeed, we have changed this mistake (line 545).
Line 630: Please provide the transfer conditions (voltage, amp, watts?) and transfer buffer when describing the Western blot protocol.
For immunoblotting of MUC13, protein lysates were transferred to 0.2 µm PVDF membranes using the Trans-Blot Turbo Transfer system (Biorad). The transfer was run using the protocol (High MW) which consisted in running for 10 min at 25 volts (V) and 1,3 amperes (A). These experimental data were added to the manuscript.
Reviewer 2
My main concern about this manuscript is that the authors analyzed MUC13 role in intestinal homeostasis and function using colorectal cancer cells. As helpful as cancer cells are, we should always be cautious about extrapolating roles in normal intestinal epithelium or IBD pathology. Obviously, these finding are also interesting in a cancer context. Using GEPIA (http://gepia.cancer-pku.cn/), I observed that MUC13 is overexpressed in colorectal cancer COAD-TCGA dataset (compared to normal colon from GTEX). Similar results were obtained previously by Gupta et al. (ref #10). I am aware that this would be difficult to confirm the main findings in a non-cancerous intestinal cell line but this limit (normal intestine using cancer cells) should be at least discussed in the manuscript.
We appreciate the reviewers’ comments and are aware of the downsides of using cancer-derived cell lines. We have performed the GEPIA analysis ourselves and have an ongoing project about the possible role of MUC13 in colorectal cancer progression. In a separate project, we are collaborating with the Gaultier Laboratory at the University of Virginia which has generated a MUC13 knockout mouse. This model will allow us to study the role of MUC13 in non-cancerous tissue. We recently received intestinal biopsies from these mice which will be stained with MUC13 and claudin antibodies to determine localization in healthy tissue. These experiments will reveal if MUC13 colocalizes with claudin on the lateral membrane in the healthy mouse intestinal tract. In future experiments, we will also address MUC13 localization and function in human intestinal organoids. We have adjusted the discussion to refer to the limitations of using cancer cell lines.
Massey et al (Micro 2021, PMC7014956) previously showed that MUC13 overexpression increased rigidity in PDAC cells and discussed involvement MUC13 link with EMT. MUC13-Her2 interaction was also associated with decrease of E-cadherin suggesting an EMT phenotype. This should be included in the discussion section.
The discussion has been adjusted to include the link with EMT.
The authors performed mass spectrometry analysis. Results are deposited on ProteomeXchange but are not yet publicly released. Among the 1189 membrane protein identified. Did the authors observed alteration of EMT proteins? (decrease of vimentin for example). In the discussion section (lane 347), the authors mentioned the relationship between other membrane bound mucins such as MUC1, MUC4, MUC16 or MUC17 and AJ/TJproteins. Did the authors observed any alteration of these mucin in the mass spectrometry data?
The mass spec analysis was performed on membrane fractions, therefore our dataset will not contain true cytosolic proteins. One of the key EMT proteins, Vimentin, is a cytosolic protein, and indeed it was not found in our dataset. Other EMT-related proteins are shown in the following table. TGF beta 1 was slightly decreased, while E-cadherin and Integrin beta 6 were slightly increased in the ∆MUC13 cells compared to WT cells.
Gene Name
Mean WT
Mean ∆MUC13
Mean MUC13-∆CT
TGFBI (TGB beta 1)
20,54
16,48
18,83
CDH1 (E-cadherin)
22,69
24,57
24,24
ITGB6 (Integrin beta 6)
18,86
21,74
19,19
Vimentin - Cytosolic
-
-
-
CDH2 (Cadherin-2, N-cadherin)
-
-
-
Mucins are large proteins comprised of densely O-glycosylated mucin domains, which makes them extremely challenging to study by mass spectrometry (MS) (Rangel-Angarita et al., 2021). We did not specifically employ mucin-directed technologies in this dataset, thus making the detection of mucins hard. No mucins other than MUC13 were detected. For MUC13, two peptides corresponding to the EGF-like domains in the extracellular domain, a region that is less densely glycosylated. We added a sentence to the description of the mass spec results to include the EMT proteins and other mucins.
Minor points:
Lane 126: HRT18 and Caco2 colon cancer cells instead of intestinal epithelial cells
Corrected.
Lane 181 and lane 514: add "full length" MUC13 DNA sequence
Corrected.
Lane 234: TEER was measured every 12h. How the authors did observed the largest increase at 42h? Was it 48h? Please clarify.
We aimed at measuring every 12 h, however the exact measurements were done at 18h, 24h, and 42 h post-infection. We have corrected this in the manuscript.
Reviewer 3
Line 43 and 46. "Enterocytes" should be replaced with "intestinal epithelial cells", since enterocytes are themselves a distinct subpopulation of IECs.
We have changed it in the manuscript.
Lines 58-60. References in support of the statements should be added.
We added a reference to this sentence.
Lines 188-190. Authors comment on "roundness" of different cell lines. If the parameter is critical for the manuscript, the authors should quantify this phenotype.
The parameter is not critical for the manuscript. We removed the sentence.
Figure 3A. Staining of cell lines should include panels showing localization of MUC13.
Co-staining of MUC13 with occludin in HRT18 cell lines can be found in figure 1D, and MUC13 with E-cadherin in supplementary figure 1.
Lines 323-327 and 390-392. Sentences on these lines contradict each other. The sentences should describe/discuss quantified data presented in Figure 6D.
The reviewer is right that we should be discussing the quantified data in 6D. We adjusted the sentence in line 323-327.
Proteomic data sets should be made publicly available on data depositories.
All proteomics raw data were deposited to the ProteomeXchange Consortium with the dataset identifier PXD029606.
Reviewer 4
OPTIONAL: In the figure 2E, is the extracellular antibody still detecting the MUC13-ΔCT?
No, unfortunately the antibody directed against the MUC13 ED is not compatible with western blot.
In the figure 2G, would be nice to comment possible reasons why the deletion in the first cell line of the MUC13-CT you can still detect with the extracellular antibody some lateral expression of MUC13 meanwhile in the second cell line, the same deletion (MUC13-CT) you cannot see any lateral MUC13 staining with the extracellular antibody.
Yes, this is indeed a puzzling finding, especially because the CRISPR deletion is the same in both cell lines. We will add a sentence about possible reduced stability of the MUC13 without CT domain that leads to a different outcome in both cell lines.
It would be nice that the results from Figure 3H are better explained since it is difficult to follow.
We adjusted the text to explain the experiment in more detail.
2. Description of analyses that authors prefer not to carry out
Reviewer 1
The authors may be overly reliant on TER measurements. Epithelial cells have two parallel resistive pathways: transcellular and paracellular. TER measure the contribution of both. Thus, an increase in TER could result from a decrease in transcellular ion transport. The authors need to measure transcellular ion flow or selectively measure the junctional resistance in a select set of experiments to rule this possibility out.
The reviewer is right that TEER is a sum of the resistance of the transcellular and paracellular pathways. However, due to the high resistance of cell membranes, the current predominantly travels via the paracellular route (Elbrecht et al., 2016). For this reason, TEER measurements are widely accepted techniques for the assessment of ions passage through the paracellular pathway (Shen et al., 2011).
Reviewer 3
Figure 1C. Caco2 and HRT18 cells exhibit distinct MUC13 expression patterns when probed with an antibody against the MUC13 CT; MUC13 localizes almost exclusively to lateral cell junction in HRT18 cells, while a higher portion of MUC13 is present on the apical surface of Caco2 cells. This observation has two possible explanations: 1) the two cell lines express distinct forms of MUC13, or 2) the two cell lines carry distinct machineries for anchoring MUC13 to apical versus apico-lateral membranes. Thus, The authors should take the opportunity to determine the impact of MUC13 deletion on TEER and TJ function in Caco2 cells. Proteomic analysis and functional assays in Caco2 cells may provide more a general mechanism for how MUC13 regulates TJ proteins.
Yes, this would be a great line of investigation. However, we aimed to knockout MUC13 in Caco-2 cell lines (with the same CRISPR/Cas9 protocol as the HRT18 cells) but were unable to obtain Caco-2 knockout clones. We think this might be a consequence of the poor capability of Caco-2 cells to grow as single colonies (a required step in the protocol). Another option is Caco-2 MUC13 knockout cells have reduced viability.
The authors generate cell lines that either lack MUC13 or express MUC13 lacking the cytoplasmic domain. Loss of MUC13 cells resulted in enhanced TEER and increased recruitment of TJ proteins to membrane fractions. MUC13∆CT cells show moderate recruitment of TJ proteins to membranes and no increase in TEER but inhibit paracellular diffusion of Luciferase Yellow across monolayers. Figure 3A suggests that Occludin redistributes to tricellular junctions in ∆MUC13 cells, whereas it is found more laterally in WT and MUC13∆CT cells. These finding suggest that full-length MUC13 interferes with TJ protein complexes. However the impact of the extracellular and intracellular (CT) domains is not fully elucidated. Does the O-glycosylated mucin domain interfere with the extracellular domains Occludin and Claudins? The authors should clarify the contribution of the mucin domain to the observed phenotype, for example by performing the described experiments in a cell line expressing MUC13 lacking the mucin domain.
Mucins are type I membrane proteins with the N-terminal part of the protein on the extracellular site. Therefore, a CRISPR method to specifically remove the glycosylated domain but leave the remainder of the protein in frame is challenging. An additional difficulty is that the ED contains a lot of repeats, complicating the design of specific guide RNAs. To specifically address the contribution of the glycosylated domain, we could complement the MUC13 knockout cell with a construct lacking the ED. However, this would not be comparable to the endogenous MUC13∆CT cell line presented in this manuscript. In future studies, we will strive to address the functions of the different MUC13 domains in more detail.
Figure 5A. Turnover of TJ proteins in membrane fractions occurs faster than over a period of 1-3 days (PMID: 18474622). The authors should determine TJ protein turnover over a period of minutes and hours.
We acknowledge the findings in this interesting paper concerning the continuous remodeling of tight junctions. However, the readout of our biotinylation assay is degradation and the timeframe of degradation turns out to be days and not hours. Within this timeframe remodeling is taking place but it cannot be captured in the total lysate.
Reviewer 4
OPTIONAL: The authors show that the probiotic Lactobacillus plantarum increase epithelial barrier independently of MUC13. Have the authors considered to use other probiotics as Lactobacillus paracasei (10.3389/fcimb.2015.00026), Akkermansia muciniphila (10.1038/emm.2017.282) or some metabolic products from intestinal microbiota as short-chain fatty acids (SCFAs) (10.3389/fphys.2021.650313) to check what is the role of MUC13 and if it is related with other microbe or microbiota metabolite?
Thank you for the suggestion. We have an ongoing project in which we investigate the impact of different probiotic bacteria and plan to investigate whether they have an impact on the epithelial barrier function in a MUC13-dependent manner. This study will lead to a separate publication.
OPTIONAL: The authors successfully delete MUC13 in IECs, both, full length and the cytosolic tail. Have the authors considered targeting the deletion of the PTS domain in MUC13? Could affect that something different from paracellular trafficking as the extracellular detection of microbes and microbial products?
Removal of a domain in the extracellular domain of MUC13 with CRISPR is challenging because mucins are type I membrane proteins, the repeats and possible frameshift, as described above.
-
- Aug 2022
-
www.biorxiv.org www.biorxiv.org
-
Note: This rebuttal was posted by the corresponding author to Review Commons. Content has not been altered except for formatting.
Learn more at Review Commons
Reply to the reviewers
Manuscript number: RC-2022-01481R
Corresponding author(s): Sebastian Voigt. Mirko Trilling, David Schwefel
1. General Statements [optional]
-
2. Description of the planned revisions
Reviewer #1: Evidence, reproducibility and clarity
Using proteome profiling of rat CMV infected cells, the authors of this study identify the E27 protein of rat cytomegalovirus as being crucial for proteasomal degradation of STAT2. Since E27 shares 56% sequence identity to the previously characterized STAT2 antagonist M27 of murine CMV the authors investigated association of E27 with the Cullin4-RING UbL CRL4. Using gel filtration chromatography they provide evidence that E27 forms a stable ternary complex with DDB1 and STAT2 suggesting that E27 bridges STAT2 to DDB1 which is further corroborated by data from cross-linking mass spectrometry. A cross-linked DDB1/DDA1/E27/STAT2 complex was then used for cryo-EM imaging experiments. The subsequent single particle analysis yielded a density map at 3.8 A resolution that was further used to generate an E27 molecular model. At this point it should be noted that resolution was not very high and data form AlphaFold2 prediction and CLMS experiments were necessary to build a model which was described as having "sufficient quality", however, no quality parameters are included for this model. In this model, a cryptic zinc-binding motif was identified that turned out to be well conserved in M27. At this point the study switches to a mutational analysis of M27: MCMV mutants either lacking M27 or bearing an AxAxxAA triple mutation were investigated both in cell culture and in animal models. Surprisingly, the M27-AxAxxA mutant while exhibiting attenuated IFN inhibition was still more active than an M27 deletion mutant. Later during the study it is postulated that this may be due to the fact that E27 binding to STAT2 abrogates the interaction with IRF9, however, this is only predicted from modeling and no experimental data are provided for this hypothesis. Furthermore, modeling approaches were used to predict how E27 replaces endogenous CRL4 substrate receptors and how E27 recruits STAT2 to mediate CRL4-catalysed ubiquitin transfer.
Reviewer #1: Significance
__Reviewer #1: __This is an interesting and well written paper describing for the first time in molecular detail how a cytomegalovirus-encoded interferon antagonist degrades STAT2 by mimicking the molecular surface properties of cellular CRL4 substrate receptors.
This study should be of broad interest for both virologists and structural biologists.
Authors Response: We thank the reviewer for the insightful and constructive evaluation. We are very grateful for highlighting the significance of our work.
Reviewer #1: Major points
__Reviewer #1: __To my opinion the authors should perform mutational analysis in the context of E27 and RCMV. I accept that switching to M27 may be easier due to established procedures for MCMV mutagenesis and analysis, however, since all structural work is primarily done on E27 it would be consequent to confirm these structural predictions in the context of E27 before switching to a related protein.
Authors Response: As the Reviewer appreciated, there were multiple reasons for the switch from RCMV-E E27 to MCMV M27. Most importantly, the MCMV in vivo infection model in mice is very well-established. Please also note that MCMV is applied far more often by virologists and immunologist as a standard model. Thus, the extension of our findings from RCMV to MCMV increases the relevance and outreach of the study. By performing the experiments in the MCMV context, we also aimed to emphasise that the function of the zinc-binding motif, which structurally organises the DDB1-binding domain, is functionally conserved among E27/M27-like proteins. Obviously, Reviewer #1 could ask why we do not solve the structure of M27 parallel to E27. With the sole exception of E27, none of the rodent M27 homologues could be produced recombinantly in a soluble form, preventing the purification and structure analysis of M27.
Since we agree with Reviewer #1 that the extension from E27 to M27 may read “a bit rough” without a mutational analysis in the E27 context, we will construct RCMV-E E27 mutants leading to Cys=>Ala exchanges in the Zn-binding motif. An analysis of the interaction between DDB1 and these E27 mutants will be included in the revised manuscript.
__Reviewer #1: __Moreover, data on the replication of the generated E27 deletion RCMV should be included in the manuscript (i.e. growth curves).
Authors Response: RCMV mutants lacking the E27 gene exhibit an impaired replication. According to the suggestion, the growth curves will be part of the revised manuscript.
Reviewer #1: The hypothesis that STAT2/E27 interaction is sterically incompatible with IRF9 binding is only based on structural prediction. It would help if the authors could present experimental evidence for such a mechanism.
Authors Response: The hypothesis is based on three lines of argumentation: (i) structural data regarding the binding interface between STAT2 and E27 covering the known STAT2-IRF9 interface (Fig. 7F) (Rengachari et al., 2018). (ii) The finding that M27 mutants incapable to bind DDB1 and induce STAT2 degradation along the ubiquitin proteasome pathway retain a residual capacity to inhibit ISRE signaling, suggesting that the binding of M27 to STAT2 suffices to elicit some signaling inhibitory functions (Fig. 7G). (iii) To elicit their function, CRL4 substrate receptors such as E27 interact with two partners. As we discussed elsewhere (Le-Trilling and Trilling, 2020), a simultaneous development of two independent traits violates evolutionary and probability theories. Thus, these receptors must acquire their binding interfaces sequentially, and the first interaction must provide an evolutionary advantage allowing the fixation of the allele in the population. Afterwards, the second binding interface evolves. Thus, a hypothesis in which E27/M27 precursors evolved the capacity to bind STAT2, preventing its association with IRF9 thereby establishing relevant but incomplete IFN inhibition (before the DDB1 interface was invented leading to STAT2 degradation by the proteasome), provides a parsimonious explanation for all these findings without violating evolutionary constraints. To corroborate our argumentation, we will analyse if E27 indeed displaces IRF9 from STAT2 by analytical gel filtration and/or co-immunoprecipitation experiments.
Reviewer #2: Evidence, reproducibility and clarity
__Reviewer #2: __The manuscript entitled "Structure and mechanism of a novel cytomegaloviral DCAF mediating interferon antagonism" by Dr. Schwefel and colleagues cleverly combines biochemistry, mass-spectrometry, Cryo-EM and cell biology to dissect how RCMV-E hijacks its hosts ubiquitylation machinery to mediate proteasomal degradation of STAT2, a key player driving the antiviral IFN response. They identify E27 as DDB1-binding element, which is able promote CRL4-dependent ubiquitylation of STAT2, and demonstrate its effect on STAT2 levels by knockout RCMV-E strains. These findings are supported by in vitro reconstitution of the DDB1/E27/STAT2 complex and analyses via XL-MS and Cryo-EM. The obtained data are then powerfully validated and analysed in mutational strains via infection of homologue in vivo models. The results collectively explain how E27 mimics endogenous CRL4 substrate receptors, thereby recruiting STAT2 to be targeted by CLR4 for ubiquitylation in a NEDD8-dependent manner.
Overall this is an important study that provides convincing insights on how rodent CMVs antagonize their host interferon response by exploiting its ubiquitin-proteasome system.
The manuscript is well written and its introduction is extraordinarily comprehensive. There are a few minor points for the authors to consider below.
Authors Response: We thank the reviewer for this very positive assessment.
Reviewer #2: Significance
Reviewer #2: The work of Schwefel and colleagues combines several powerful state-of-the art techniques to dissect the mechanism of the viral protein E27 and, for the first time, provides a rational for its ability to act as STAT2 antagonist. They performed outstanding structure-function analyses of the ubiquitin system, including the first global proteomic profiling of RCMV-infected cells, setting the standard for its human counterpart as rodent CMVs are commonly used as infection models. The manuscript is highly suitable for publication in any of the journals associated with the review commons platform.
Authors Response: Again, we thank the reviewer for these kind words and the appreciation of our work.
Reviewer #2: CROSS-CONSULTATION COMMENTS
Reviewer #2: This reviewer agrees that at least testing mutants in the E27 in some assays would be appropriate.
Authors Response: As detailed in the response to Reviewer #1, we will generate RCMV-E E27 mutants targeting the Zn-binding motif by site-directed mutagenesis. An analysis of the interaction between DDB1 and these E27 mutants will be included in the revised manuscript.
Reviewer #3: Evidence, reproducibility and clarity
__Reviewer #3: __Le-Trilling et al. present the first proteomic analysis of RCMV-infected cells, where they identified STAT2 as one of the most heavily downregulated (and degraded) proteins. This analysis showed that RCMV mediated degradation of STAT2 is conserved in closely related species used as animal models (rat and mouse) and human, despite the intra-host adaptation of each CMV. They also identify E27 as the RCMV factor that targets STAT2 for degradation, that exhibits ~50% homology with MCMV pM27. This study also identifies a Zinc binding motif in E27 using Cryo-EM which is conserved in other CMV species and is potentially involved in antagonising Type I and III responses.
Reviewer #3: Significance
__Reviewer #3: __The present work provides the first proteomics analysis of RCMV infection in rat cells, comparing infected vs non-infected rat fibroblasts to access potential RCMV targets. Then, it focuses on the characterisation of RCMV E27 and its role targeting and interacting with STAT2 (plus recruiting the Cul4 complex for STAT2 degradation). Finally, it provides the Cryo-EM structure of E27 and its CMV homologues, and the structure of the complex of E27 with elements of the CUL4 complex and STAT2. This is the first time that E27 function and structure are characterised. These are all novel findings - although the mouse homologue M27 has previously been found to interact with and degrade STAT2 (published by some of the same authors in Plos pathogens in 2011, (https://doi.org/10.1371/journal.ppat.1002069). Therefore the chief novel information is the structural studies.
The manuscript will be of interest to researchers working with human and animal herpesviruses.
My field of expertise is in Virology, Innate Immunity and host-virus interactions from an evolutionary perspective. I do not have expertise in Cryo-EM, so I could not evaluate the methods used in the section.
__Authors Response: __We thank the reviewer for the positive evaluation of our work and its significance.
Reviewer #3: Major points
__Reviewer #3: __1. The authors claim the identification of a Zinc-binding motif in the protein E27 (RCMV) using Cryo-EM, then validation of the phenotype with MCMV WT, delM27 and M27 AxAxxA. To justify the change to MCMV to perform the functional validation, they stated "MCMV M27, the closest E27 homologue, exhibits 56% and 76% amino acid sequence identity and similarity, respectively (Fig. S4B). E27 and M27 AlphaFold2 structure predictions are almost indistinguishable (RMSD of 1.195 Å, 6652 aligned atoms) (Figs. 3B, S4A), and structural alignment of these predictions demonstrated conservation of side chain positions involved in zinc-binding (Fig. 3C). Thus, M27 represents a valid model to study functional consequences of interference with the zinc coordination motif through site-directed mutagenesis, and to test the predictive power of our E27/M27 model". Although they rationalise the change to MCMV to validate the functional outcomes of the newly identified zinc binding motif with alignments and Cryo-EM data, it falls within the DDB1 binding region that is less conserved (Fig S4B). The addition of a mouse model here provides a solid result but given the aim of the paper is to provide a proper characterisation of RCMV and elucidate some inter-species adaptations, I strongly recommend the validation with E27 here given the potential impact of this motif. Rather than having to repeat this in a rat model (which would clearly be a large amount of work), this could simply be achieved by constructing the relevant deletion / mutant viruses and assessing in vitro in a relevant cell line (readout - either virus titre or luciferase assay as shown in Figure 3G/H).
__Authors Response: __Please also see our responses to the other reviewers. Briefly, we will apply side-directed mutagenesis to alter the CxCxxC motif in E27 that binds the zinc ion, and analyse the interaction of these E27 mutants with DDB1. In this context, we would like to add that almost two thirds of E27 residues in direct contact with DDB1 are at least type-conserved in M27, and the zinc-coordinating side chains are totally conserved (Fig. 3C). Together with a predicted similar structural organization of the respective binding regions (Fig. S11), and in light of our MCMV mutagenesis results (Fig. 7), it is highly likely that the DDB1-binding mode is conserved between E27 and M27. As mentioned above, we will put this assumption to the test in the revision process.
__Reviewer #3: __Furthermore, in Figure 2, the GF assay was performed using full-length DDB1, however CLMS was performed using DDB1 delBPB (interchange between these two proteins continues in the remainder of the paper). This should be at least justified, and preferably one or other of wt DDB1 and DDB1 delBPB used in the GF or CLMS assay where this has not yet been performed. Later on in the results section (Fig 5E), the authors use wt DDB1 while in fig 4 they used the delBPB to describe the interaction with E27 - would be relevant to have consistency across the paper and some supplementary data that could support using one or the other in each assay.
__Authors Response: __Protein complex preparations including full length DDB1 did not yield cryo-EM reconstructions at appropriate resolution for model building, almost certainly due to the known flexibility of the DDB1 BPB, impeding proper alignment of the cryo-EM particle images. This is why we switched to DDB1ΔBPB. Importantly, the structure model including full length DDB1 (Fig. S12B) clearly demonstrates that the BPB is located on the opposite side of the E27 binding interface on DDB1 (where it is situated to flexibly connect to the CUL4 scaffold to create the ubiquitination zone around immobilised substrates [Fig. 6]). This rules out an involvement of DDB1 BPB in E27- and/or STAT2-binding processes. Several previous studies have employed DDB1ΔBPB to facilitate structure determination, and have successfully applied the resulting structural models for functional follow-up experiments in the context of complete CRL4 assemblies (Bussiere et al., 2020; Petzold et al., 2016; Slabicki et al., 2020). Nevertheless, we will repeat GF experiments with DDB1ΔBPB for consistency and include these data in the revised manuscript.
Reviewer #3: Minor points
__Reviewer #3: __2. Although they present sufficient detail in the methods, further details in the text should be given as to the number of repeats performed in each case, and whether the data shown is representative or based on an average of repeats (preferably the latter; if representative, the data for other repeats should be shown in supplementary information).
Authors Response: We will add this information in the revised version of the manuscript.
3. Description of the revisions that have already been incorporated in the transferred manuscript
Reviewer #1: Major points
__Reviewer #1: __Resolution of the cryoEM structure is rather low and many predictions of the manuscript are based on modeling using AlphaFold2 prediction. The authors describe their model as of "sufficient quality", however, no quality measures are included in the manuscript. At least the discussion should address limitations of the used approach.
Authors Response: While we apologize for not sufficiently describing our quality measures, we respectfully disagree regarding the conclusion. Our resolution (3.8 Å, map 1) lies well within the 3–4 Å resolution range of the vast majority of structures deposited to the Electron Microscopy Data Bank during the last five years (https://www.emdataresource.org/statistics.html). Nevertheless, de novo modelling in this resolution regime is challenging. This is why we sought additional guidance through cross-linking mass spectrometry (XL-MS) restraints and AlphaFold2. Please also note that modelling of E27 was not based solely on the AlphaFold2 prediction. Instead, a partial model corresponding to the α-domain was manually built in map 1, guided by XL-MS information (see Methods - “Model building and refinement” and Fig. S5B, grey cartoon). This partial model proved to be in very good agreement with AlphaFold2 predictions (RMSD of 1.489 Å, 2764 aligned atoms). Only after this initial sanity check, the computational prediction was used for model completion, adjustment, and refinement.
We now added graphical overviews of model fits in Figs. S5 and S10. Furthermore, we included detailed views of the fit of relevant side chains involved in intermolecular interaction to the experimental density (Fig. S7, S9). We also calculated and listed quality indicators of the model-to-map fit in Table S1 (correlation coefficients and model resolution based upon model-map FSC). To ensure the validity of our atomic model using an alternative method besides cryo-EM and XL-MS, we have performed site-directed mutagenesis of critical binding regions in E27, followed by in vitro reconstitution and analytical GF (Fig. S7B, C, S9B, C). The text was revised accordingly (see p10 [ll22] and p14 [ll26]).
__Reviewer #1: __The authors identify a cryptic zinc-binding motif in E27 that is conserved in homologous proteins. For this reviewer it is not clear: is there experimental evidence for zinc binding of E27 or can the presence of zinc reliably be detected in their structural data? If not, it would be worth to confirm zinc binding.
Authors Response: Our structural data show a tetragonal metal coordination geometry, involving three cysteine side chains and one histidine side chain, with coordination bond lengths of 2.2 Å between the histidine nitrogen and the metal ion, and of 2.4 Å between the cysteine sulfurs and the metal ion. The density feature cannot be explained by another type of side chain interaction, e.g. a disulfide bond, because this would lead to a steric clash with the remaining adjacent side chains. Based on the knowledge on metal-binding sites in proteins and metal-coordination chemistry, these characteristics indicate the presence of a structural zinc-binding site for the following reasons: (i) after magnesium, zinc is the second most prevalent metal in the Protein Data Bank (https://metalpdb.cerm.unifi.it/getSummary), however, magnesium is coordinated octahedrally by oxygen ligands (Tang and Yang, 2013); (ii) the most abundant zinc ligands are cysteine and histidine; (iii) the most abundant zinc coordination number is four ligands; (iv) the average coordination bond lengths are 2.12±0.19 Å and 2.33±0.12Å for nitrogen-zinc and sulfur-zinc interactions, respectively (Ireland and Martin, 2019; Laitaoja et al., 2013), which is in very good agreement with our structural observations. We included this argumentation in the revised manuscript (see p9 [ll21]), and added Fig. S5C for visualization.
Reviewer #2: Minor points
Reviewer #2: Page 2, line 3. "Here," should be inserted before "Global proteome profiling..." to highlight the work of this manuscript.
Authors Response: We changed the text accordingly.
Reviewer #2: Page 3, line 21. "IFNs" instead of "IFN"
Authors Response: We changed the text accordingly.
Reviewer #2: Page 4, lines 9,15,27. "Ubiquitin Ligases (UbL)" is not a common abbreviation and could be mistaken for Ubl (Ubiquitin-like proteins). Possible abbreviation is "E3s" for Ubiquitin E3 ligases
Authors Response: We have amended the respective abbreviations accordingly.
Reviewer #2: Page 4 line 25. "RBX1" is the more common term for "ROC1"
Authors Response: This has been corrected throughout the manuscript.
Reviewer #2: Page 5 lines 1-9. Citing of the first structure of DDB1 in complex with a viral protein is recommended. (Ti Li et al. Cell 2006)
Authors Response: We thank the reviewer for this important suggestions and cited this landmark publication.
Reviewer #2: Figure 1 a) STAT2 dot is cut off in second panel. I recommend highlighting STAT2 in both panels.
We amended the figure accordingly. We furthermore additionally highlighted the “STAT2” text in both panels by increasing the font size and putting it in bold type.
Reviewer #2: Page 7 line 17. "Cross-linking MS (CLMS)" is commonly abbreviated as (XL-MS)
Authors Response: We changed the text accordingly.
Reviewer #2: Figure 2 a-c) These panels could benefit from thinner lines in order to increase visibility of chromatograms and cross-links.
Authors Response: The panels were changed accordingly.
Reviewer #2: Figure 2 a-b) Could the authors elaborate on why STAT2 is stoichiometrically
underrepresented in the SDS-PAGE of the E27/DDB1/STAT2 complex?
Authors Response: We applaud Reviewer #2 for their in-depth examination. Honestly, we were also puzzled by this. Based on the cryo-EM single particle analysis, we found an explanation: We separated a major contamination in silico during 2D classification (~12% of all particles). Out of curiosity, we reconstructed a density map from these particles (now shown in Fig. S3). The map was identical to a previous cryo-EM structure of the E. coli protein ArnA (Yang et al., 2019), a notorious contaminant in E. coli Ni-NTA protein purifications (Andersen et al., 2013). ArnA migrates similar to E27 on the SDS-PAGE, the band runs just a little bit faster (compare fraction 6 [ArnA] and fractions 8/9 [E27] from the SDS-PAGE of the analytical GF run of E27 in isolation, Fig. 2A, green trace). However, in analytical GF, ArnA elutes at higher molecular weight fractions, since it forms a hexamers (Ve~10.2 ml). Incidentally, this elution volume of the ArnA hexamer almost equals the one of DDB1 or DDB1ΔBPB/DDA1/E27/STAT2 complexes. This leads to a superposition of ArnA and E27 bands in the respective SDS-PAGE lanes corresponding to GF fraction 6. Accordingly, we conclude that it is actually not STAT2 that is underrepresented, but rather E27 seems overrepresented due to SDS-PAGE band overlap with the ArnA contaminant. We have now indicated the contaminant in Fig. 2A, amended the legend, and extended Fig. S3 to indicate at which point of the cryo-EM analysis the contaminating ArnA particles were separated, and to show the ArnA model to map fit.
In addition to this, it might be that potential STAT2 degradation products (marked by ** in Fig. 2), which seem to co-migrate with STAT2/E27 complexes, occupy FL STAT2 binding sites on E27.
Reviewer #2: Paragraph "The E27 structure.." page 9. Placing this paragraph after the overall
structure is recommended.
Authors Response: Accordingly, we have now moved this section to the end of the results section.
Reviewer #2: Figure 3 a) The grey mesh being laid over the ribbon structures is not contributing to the overall visibility. Adding a panel of the cryo-EM structure alone in cost of alphafold models is recommended.
Figure 4a) same issue with grey mesh
Authors Response: Thank you very much for the very good suggestions. We have removed the mesh representation, and included panels just showing the segmented cryo-EM map in the new Fig. 3A.
Reviewer #2: c) panels could benefit from fewer amino acids being labeled/shown
Authors Response: We understand the motives of the Reviewer. However, we would prefer to depict all relevant side chain interactions in these panels. The rearrangement of the figure, i.e. showing the overview of the interacting regions before the detailed panels, should make them more accessible (new Fig. 3B).
__Reviewer #2: __d) may want to avoid red-green coloring to improve for colorblindness
Authors Response: We are deeply sorry for our ignorance in this regard. We changed the colors accordingly (see new Fig. 3B, C).
__Reviewer #2: __Figure 6a) s.a grey mesh
Authors Response: We removed the mesh representations and included panels just showing the segmented cryo-EM density in the new Fig. 5C.
Reviewer #2: CROSS-CONSULTATION COMMENTS
__Reviewer #2: __A 3.8 A overall resolution map and the approach to fitting may be suitable, but it is unclear from the authors' figures whether the side-chains shown in the figures are clearly visible in the map or if they are modeled by some other approach. Side chains should ideally be visible in the maps if shown in figures, and if not, close-ups of the corresponding regions of the maps should be shown with sufficient depthcue to allow the reader to gauge how the map corresponds to the model.
Authors Response: This is a crucial point. As mentioned in the response to Reviewer #1, major point 2, we have now included very detailed views of the fit of relevant side chains involved in intermolecular interaction to the experimental density (Fig. S7, S9).
__Reviewer #2: __Along these lines, the figures with the mesh maps do not clearly show how well the model fits the map. This needs to be clearly visible in figures, and ideally maps and models provided to reviewers in order for the reviewers to gauge the level of accuracy of the fit.
Authors Response: Please see our response to Reviewer #1, major point 2. Briefly, we have now included graphical overviews of model fits in Figs. S5 and S10. We also calculated and listed quality indicators of the model-to-map fit in Table S1 (correlation coefficients and model resolution based upon model-map FSC). To ensure the validity of our atomic model using an alternative method besides cryo-EM and XL-MS, we have performed site-directed mutagenesis of critical binding regions in E27, followed by in vitro reconstitution and analytical GF (Fig. S7B, C, S9B, C). The text was extended accordingly (see p10 [ll22] and p14 [ll26]).
__Reviewer #2: __At minimum, the authors have nicely assembled proteomics and cell biological data indicating that E27 hijacks CRL4 to turn over Stat2 in rat cells in a manner paralagous to M27 hijacking in mouse cells, biophysical/structural data for a model of a CUL4-DDB1-E27-Stat2 complex, and mutagenesis of a putative zinc binding site in M27.
I feel most of the issues raised by all 3 reviewers could be addressed in the text, with more clarity about the structural models, and better explanation for why the construct with proteins from various organisms were used for structural studies (the authors had made human DDB1 before, and it expressed well, and perhaps didn't consider to make from rat? Or this mixture expressed, purified best? Gave best quality EM data?).
Authors Response: We thank Reviewer #2 for her/his overall assessment. As mentioned in the two cross-consultation comments before, and in the response to Reviewer #1, major point 2, we strived to provide adequate measures allowing to judge the quality of our structural models in the present updated version of the manuscript. In addition, as indicated in the response to reviewer #3, major point 2, we have now added Fig. S12 and extended the Discussion to explain and justify the use of different protein constructs.
__Reviewer #2: __Also, the presentation of the zinc binding site should come after the overall structure. As for the use of MCMV to assess the role of the zinc binding site, placing this last in the text might allow this to flow better.
Authors Response: Thank you very much for this suggestion. The manuscript has been restructured as recommended: details of the zinc-binding motif and the MCMV assays are now shown in Fig. 7 and described in the text just before the Discussion.
Reviewer #3: Major points
__Reviewer #3: __2. Given that previous data in mice showed that the E27 homologue pM27 binds a component of host Cullin4-RING UbLs (CRL4), to induce the poly-ubiquitination of STAT2, the current study also addressed if this mechanism was preserved in RCMV. Yet, they seemed to do this with E27, rnSTAT2 and hsDDB1 - Page 7 lines 1 to 3: "These results prompted us to explore the association of E27 with Rattus norvegicus (rn) STAT2 and Homo sapiens (hs) DDB1 in vitro. Importantly, 1128 of 1140 amino acids are identical between hsDDB1 and rnDDB1 (...)". They identify the residues and regions where the DDB1 is different between both species, but should provide a structure/alignment with this highlighted. In addition, DDB1 is a DNA damage protein that is annotated in the Rattus norvegicus genome. The authors should justify the assays between rnSTAT2-hsDDB1 instead of using the both proteins from rn, and present the equivalent data for rnDDB1 in the paper.
Authors Response: Among the 12 alterations between human and rat DDB1, 4 are type-conserved (Fig. S12A). Thus, >99% of amino acids are identical or similar. We mapped all exchanges on a model of full length human DDB1 bound to E27 and the rat STAT2 CCD. None are involved in intermolecular interactions (Fig. S12B, C). Please note that due to the high conservation of DDB1 across eukaryotes, this inter-species approach has been used by us and others to study DDB1-containing complexes (e.g., the SV5V, WHX, SIV Vpx and Vpr, zebrafish DDB2, and chicken CRBN proteins have been in vitro reconstituted with human DDB1 for structural characterisation) and valid biological conclusions have been drawn from these studies (Angers et al., 2006; Banchenko et al., 2021; Fischer et al., 2014; Fischer et al., 2011; Li et al., 2006; Li et al., 2010; Schwefel et al., 2015; Schwefel et al., 2014; Wu et al., 2015).
Reviewer #3: Minor points
__Reviewer #3: __1. In fig 5D, the authors present the H-box alignment, where it is clear that this motif is not conserved. The lack of H-box conservation should be discussed in the results and discussion, to provide an explanation for the competition/binding observed.
Authors Response: We respectfully disagree. There is conservation of amino acid side chains, regarding their physicochemical properties, observable in the H-box motif. Furthermore, the secondary structure is conserved. Please note, that the H-box is not our invention but rather represented a well-accepted motif known in the field, see e.g., (Li et al., 2010). We extended the discussion to cover this point (p21 [ll15]).
__Reviewer #3: __3. The authors commence their abstract justifying the study on the grounds of the usefulness of rodent HCMV counterparts as common infection models for HCMV. They should return to this theme in the discussion - what is the usefulness of their findings with regards to HCMV (particularly given the relatively low homology between E27 and HCMV pUL27, and the alternative mechanism for STAT2 antagonism encoded by HCMV UL145)?
Authors Response: We extended the discussion in this regard. Briefly, our data, to our knowledge for the first time, reveal that RCMV (like MCMV) exploits CRL4 to induce proteasomal degradation of STAT2. With pUL145, HCMV relies on an analogous protein. In clear contrast to HCMV, RMCV and MCMV are both amenable to in vivo experiments in small animal models. Over 40 years ago, HCMV has been called the troll of transplantation due to its grim impact on immunosuppressed individuals after transplantation surgery (Balfour, 1979). Despite tremendous efforts, HCMV still harms and kills graft recipients. While MCMV allows various experiments regarding general principles of cytomegaloviral pathogenesis and antiviral immunity, one shortcoming is that the mouse obviously is a rather small animal, preventing various chirurgical and solid organ transplantation (SOT) procedures. In clear contrast, SOT procedures that are indispensable for human medicine can be recapitulated in rat models. Thus, according to our opinion, our work lays the molecular foundation for future studies addressing the relevance of STAT2 and CMV-induced STAT2 degradation in rat SOT models.
4. Description of analyses that authors prefer not to carry out
-
- *
References
Andersen, K.R., Leksa, N.C., and Schwartz, T.U. (2013). Optimized E. coli expression strain LOBSTR eliminates common contaminants from His-tag purification. Proteins 81, 1857-1861.
Angers, S., Li, T., Yi, X., MacCoss, M.J., Moon, R.T., and Zheng, N. (2006). Molecular architecture and assembly of the DDB1-CUL4A ubiquitin ligase machinery. Nature 443, 590-593.
Balfour, H.H., Jr. (1979). Cytomegalovirus: the troll of transplantation. Arch Intern Med 139, 279-280.
Banchenko, S., Krupp, F., Gotthold, C., Burger, J., Graziadei, A., O'Reilly, F.J., Sinn, L., Ruda, O., Rappsilber, J., Spahn, C.M.T., et al. (2021). Structural insights into Cullin4-RING ubiquitin ligase remodelling by Vpr from simian immunodeficiency viruses. PLoS pathogens 17, e1009775.
Bussiere, D.E., Xie, L., Srinivas, H., Shu, W., Burke, A., Be, C., Zhao, J., Godbole, A., King, D., Karki, R.G., et al. (2020). Structural basis of indisulam-mediated RBM39 recruitment to DCAF15 E3 ligase complex. Nat Chem Biol 16, 15-23.
Fischer, E.S., Bohm, K., Lydeard, J.R., Yang, H., Stadler, M.B., Cavadini, S., Nagel, J., Serluca, F., Acker, V., Lingaraju, G.M., et al. (2014). Structure of the DDB1-CRBN E3 ubiquitin ligase in complex with thalidomide. Nature 512, 49-53.
Fischer, E.S., Scrima, A., Bohm, K., Matsumoto, S., Lingaraju, G.M., Faty, M., Yasuda, T., Cavadini, S., Wakasugi, M., Hanaoka, F., et al. (2011). The molecular basis of CRL4DDB2/CSA ubiquitin ligase architecture, targeting, and activation. Cell 147, 1024-1039.
Ireland, S.M., and Martin, A.C.R. (2019). ZincBind-the database of zinc binding sites. Database (Oxford) 2019.
Laitaoja, M., Valjakka, J., and Janis, J. (2013). Zinc coordination spheres in protein structures. Inorg Chem 52, 10983-10991.
Le-Trilling, V.T.K., and Trilling, M. (2020). Ub to no good: How cytomegaloviruses exploit the ubiquitin proteasome system. Virus Res 281, 197938.
Li, T., Chen, X., Garbutt, K.C., Zhou, P., and Zheng, N. (2006). Structure of DDB1 in complex with a paramyxovirus V protein: viral hijack of a propeller cluster in ubiquitin ligase. Cell 124, 105-117.
Li, T., Robert, E.I., van Breugel, P.C., Strubin, M., and Zheng, N. (2010). A promiscuous alpha-helical motif anchors viral hijackers and substrate receptors to the CUL4-DDB1 ubiquitin ligase machinery. Nature structural & molecular biology 17, 105-111.
Petzold, G., Fischer, E.S., and Thoma, N.H. (2016). Structural basis of lenalidomide-induced CK1alpha degradation by the CRL4(CRBN) ubiquitin ligase. Nature 532, 127-130.
Rengachari, S., Groiss, S., Devos, J.M., Caron, E., Grandvaux, N., and Panne, D. (2018). Structural basis of STAT2 recognition by IRF9 reveals molecular insights into ISGF3 function. Proceedings of the National Academy of Sciences of the United States of America 115, E601-E609.
Schwefel, D., Boucherit, V.C., Christodoulou, E., Walker, P.A., Stoye, J.P., Bishop, K.N., and Taylor, I.A. (2015). Molecular Determinants for Recognition of Divergent SAMHD1 Proteins by the Lentiviral Accessory Protein Vpx. Cell host & microbe 17, 489-499.
Schwefel, D., Groom, H.C., Boucherit, V.C., Christodoulou, E., Walker, P.A., Stoye, J.P., Bishop, K.N., and Taylor, I.A. (2014). Structural basis of lentiviral subversion of a cellular protein degradation pathway. Nature 505, 234-238.
Slabicki, M., Kozicka, Z., Petzold, G., Li, Y.D., Manojkumar, M., Bunker, R.D., Donovan, K.A., Sievers, Q.L., Koeppel, J., Suchyta, D., et al. (2020). The CDK inhibitor CR8 acts as a molecular glue degrader that depletes cyclin K. Nature 585, 293-297.
Tang, S., and Yang, J.J. (2013). Magnesium Binding Sites in Proteins. In Encyclopedia of Metalloproteins, R.H. Kretsinger, V.N. Uversky, and E.A. Permyakov, eds. (New York, NY: Springer New York), pp. 1243-1250.
Wu, Y., Koharudin, L.M., Mehrens, J., DeLucia, M., Byeon, C.H., Byeon, I.J., Calero, G., Ahn, J., and Gronenborn, A.M. (2015). Structural Basis of Clade-specific Engagement of SAMHD1 (Sterile alpha Motif and Histidine/Aspartate-containing Protein 1) Restriction Factors by Lentiviral Viral Protein X (Vpx) Virulence Factors. The Journal of biological chemistry 290, 17935-17945.
Yang, M., Chen, Y.S., Ichikawa, M., Calles-Garcia, D., Basu, K., Fakih, R., Bui, K.H., and Gehring, K. (2019). Cryo-electron microscopy structures of ArnA, a key enzyme for polymyxin resistance, revealed unexpected oligomerizations and domain movements. J Struct Biol 208, 43-50.
-
- Jan 2021
-
dimartinobooth.com dimartinobooth.com
-
In the meantime, I invite you to partake of this week’s tribute to the current generation of traders who’ve withstood the destruction of price discovery at the hands of overly-intrusive central banking policy. Please enjoy, The Unbearable Lightness of Trading.
destruction of price discovery
-