Reviewer #1 (Public Review):
Summary:
In this manuscript, the authors reveal a new role for SDG7 in the regulation of H3K36me2 and me3. SDG7 appear to be functionally redundant to SDG8 as the double mutant presents lower levels of H3K36me and stronger phenotypes than either single mutant, however, their mechanisms of action might differ as the proteins displayed different localization on their target genes, with SDG7 localizing preferentially to TSS and TES while SDG8 covers the gene body. SDG7 binds preferentially to PREs, which recruit PRC2 for H3K27me3 deposition. The authors therefore present an interesting model where SDG7 evicts PRC2 from silenced genes, leading to a loss of H3K27me3. This would allow the transcriptional activation of the genes and the deposition of H3K36me3.
Strengths:
Overall, the manuscript is well-written and organized, although some paragraphs need clarifications. The figures are clear and well designed and the proposed model is compelling. While the manuscript is already interesting as it is, I think addressing the following questions would elevate it even more and refine the proposed model:
Weaknesses/potential aspects to address:
(1) It is still unclear whether SDG7 directly catalyzes H3K36me or if it promotes its deposition simply by eviction of PRC2. The AlphaFold and structure analyses show a significant similarity between the catalytic domains which would support the first possibility, but some more experiments would be required to prove this more definitively.
(2) Does SDG7 directly recognize the PRE (as suggested by the model in Figure 5F) or is it recruited by some transcription factors? Is SDG7 known to interact with any of the PRC2 recruiters?
(3) Line154/Figure 2A: The metagene plot for H3K36me3 shows a lower level on the gene body but a higher peak in sdg7sdg8 double mutants compared to the Wild-type, which is a bit surprising, especially considering that the immunostaining in reference 19 showed a near complete loss of H3K36me3 signal in the same double mutant. Can this higher peak be an artifact from the normalization strategy, or due to the existence of different subpopulations of genes?
Indeed, on the genome tracks presented by the authors, the hypomethylated genes show a loss of signal on the entire gene body, and not a higher peak near the TSS. It might be interesting to generate metagene plots for H3K36me3 hypo and hyper-methylated genes, to see if the higher peak at the TSS is solely due to the hyper-methylated genes.
(4) Figure 2C: More than 40% of differentially methylated genes are actually hypermethylated, but the authors do not discuss this at all. What are those genes, are they targeted by SDG7 or 8? Could they be responsible for the higher peak at the TSS observed in the double mutant? (see previous comment).
(5) Figure 2C and D: The method section states that the ChIP-seq was performed on 5-day-old seedlings, while the legend of this figure mentions root and shoot samples but this does not appear in the figure itself. There is also mention of shoot and root samples in Supplementary Tables 1 and 2. The authors should clarify which tissue was used for the data presented in Figure 2 and correct the legends or the methods accordingly.
(6) Line 270/Fig 4K and L: The text mentions looking at the 838 genes "downregulated in clf sdg7 sdg8 relative to sdg7 sdg8" and in the overlap, the authors identified FLC. However, in Figure 5D, FLC is upregulated in clf sdg7-sdg8 compared to sdg7-sdg8, not downregulated as mentioned in line 270. The Venn diagram in Figure 4L mentions "sdg vs clf sdg up", which would fit the pattern seen in Figure 5D, but the number of genes (838) matches the number of downregulated genes in the sdg7sdg8 vs clf s dg7sdg8 volcano plot.
I would actually expect the phenotype rescue to be caused by genes that are up in Wt vs clf, down in Wt vs sdg7-sdg8, and back up in sdg7-sdg8 vs clf-sdg7-sdg8, not "up/down/down" as mentioned in the text: genes would be downregulated in sdg7-sdg8 because of a loss of H3K36me and therefore hypermethylation of H3K27, but in the absence of CLF, this hypermethylation is reversed and the genes are upregulated in the triple mutant compared to the sdg7-sdg8 mutant. This is also what the authors see and describe in their cluster analysis in Figure 4M and line 280, mentioning an upregulation in clf-sdg7-sdg8 vs sdg7-sdg8. Could the authors please clarify these discrepancies between the different subplots and within the text itself? Was there maybe some error plotting the volcano plot and/or Venn diagram?
In general, as this part is quite complicated, maybe it would benefit from a clearer explanation from the authors as to why they look at those particular overlaps, so that the reader can more easily follow their train of thought.
(7) Figure 4N/Line 286: How were these 828 genes identified? Is it stemming from a clf-sdg7-sdg8 vs sdg7-sdg8 comparison? The legend says "genes shown by white color in Fig. 4M", do the authors mean the two clusters previously described?
(8) Line 300: "suggesting that SDG8 primarily mediates target gene expression in conjunction with PAF1C". This statement is based on overlapping genes that are downregulated in sdg7-sdg8 double mutant and paf1c mutants but concludes only on the role of SDG8. I feel that to state that SDG8 regulates expression in conjunction with PAF1C, the authors should rather examine the genes downregulated in the sdg8 mutant, especially considering the reduced overlap between genes downregulated in sdg8 and sdg7-sdg8 (according to Figure 2C, only 30% of the genes downregulated in sdg8 are also downregulated in the double mutant), or this statement should be corrected to also include SDG7.
Maybe it would be easier to read the figure if the authors created a master list of genes downregulated in at least one of the paf1c mutants they examined (as they anyway do not examine in detail the contribution of each individual paf1c mutant), and overlap it with the genes downregulated in sdg7, sdg8 or sdg7-sdg8.
(9) Line 326: "We also discovered that SDG7 and SDG8 overcome PRC2-mediated silencing, leading to a switch from H3K27 methylation to H3K36 methylation during growth and development." While part of this statement is supported by the ChIP data presented in Figure 4E, I think a ChIP for H3K36me2 and/or me3 is necessary to prove the existence of a K27me to K36me switch.
(10) Line 347: The authors state that SDG8 is located at the TSS and 3' end of genes, but on line 187 they state that it occupies the gene body (which is supported by the plot in Figure 3A).
(11) Line 351: The authors suggest a role of RNApolII in the deposition of K36me, but their data are not sufficient to support this hypothesis. The transcriptome data show that both SDGs and PAF1C regulate a similar set of genes, but they do not show data demonstrating that RNApolII is necessary for the deposition of K36me. It might be interesting to examine H3K36me levels in a paf1c mutant to further consolidate their hypothesis.